Performance of the Italian strong motion network during the 2009, L��Aquila seismic sequence (central Italy)

On April 6, 2009, the town of L’Aquila in the Abruzzo region (central Italy) was struck by a seismic event at 01:32 (UTC), of magnitude M_W = 6.3. The mainshock was followed by a long period of intense seismic activity and within seven days after the mainshock there were seven events of magnitude M_W ≥ 5 that occurred from April 6 to April 13. This long seismic sequence was characterized by a complex rupture mechanism that involved two major normal faults of the central Apennines: the Paganica and the Gorzano faults. The strong-motions of the mainshock were recorded by 64 stations of the Italian Strong-motion Network (RAN) operated by the National Civil Protection Department (DPC). Six stations of a local strong-motion array were working in NW L’Aquila suburb area. One of them, located at about 6 km from the Paganica fault surface tip-line, set up in trigger mode, recorded continuously for more than 20 min the mainshock and the aftershocks. Besides the mainshock, the RAN stations recorded in total 78 foreshocks and aftershocks of M_L ≥ 3.5, during the period from January to December 2009. The corresponding waveforms provide the most extensive digital strong ground motion data set ever recorded in Italy. Moreover, the 48 three-component observations of events of magnitude M_W ≥ 5, recorded at a distance less than 15 km from each of the major involved faults, provide a significant increasing of near-field records available for the Italian territory. Six days after the mainshock, the strong-motion dataset, referred to preliminary locations of the events with M_L ≥ 4.0, was made available on the DPC web site () and at the same time it was delivered to the ITACA database (). This dataset has been used by many authors in scientific papers and by engineers, geophysicists and geologists for professional technical works. In this paper, the present-day available strong-motion signals from the L’Aquila sequence and the performance of the Italian strong-motion network in terms of the number and quality of recorded data, the geometry and data transmission system are described. In addition the role of the temporary network that represents an extension of the permanent Italian strong-motion network, supporting the emergency response by civil protection authorities and improving the network coverage has been evaluated.     

Peak and integral seismic parameters of L��Aquila 2009 ground motions: observed versus code provision values

The M _w 6.3 L’Aquila earthquake of April 6, 2009 hit a wide area of the Abruzzo region (Central Italy). The epicentre of the main shock was very close to the urban centre of L’Aquila, the regional capital, with an epicentral distance less than 10 km. It was the strongest earthquake ever recorded in Italy which has provided ground motion recordings from accelerometric stations located in close proximity to the epicentre. Because of this, several remarkable results can be achieved by analysing the strong motion recorded signals in terms of peak (PGA, PGV and PGD) and integral (Housner Intensity, I _H) seismic parameters. Additionally, an alternative time-domain representation of recorded signals has been used to furnish a rapid comparison of traces recorded at different stations and along different directions. Some comparisons between the response spectra derived from the recordings and the elastic demand spectra provided in the new seismic Italian code have also been performed. PGA recorded values are very high and generally higher than code values for seismic actions with return period T _R = 475 years. In some cases, this also happens for seismic actions with T _R up to 2,475 years. With regard to I _H, recorded values are generally higher for T _R = 475 years, whilst they are remarkably lower for T _R = 2,475 years. Accurate analyses have been carried out in the article to better understand the above differences and their significance and implications.     

Characteristics of the strong ground motions from the 6 April 2009 L’Aquila earthquake,Italy

An M_w 6.25 earthquake occurred on April 6, 2009 at 03:33 a.m. local time, in the Abruzzo region (Central Italy), close to the city of L’Aquila. The earthquake ruptured a North-West (NW)–South-East (SE) oriented normal fault dipping toward the South-West (SW), with the city of L’Aquila lying a few kilometers away on the hanging wall.The main shock has been recorded by fifty-eight accelerometric stations: the highest number of digital recordings ever obtained in Italy for a single earthquake, one of the best-recorded earthquakes with a normal fault mechanism. Very high values of peak ground acceleration (0.3–0.65 g) were observed close to the center of L’Aquila (6 stations at zero JB distance from the fault). The earthquake caused severe loss of lives (299 victims and 1500 injured) and damage (about 18000 unusable buildings) in the epicentral area.In this study we analyze the ground motion characteristics of both the main shock in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and pseudo-acceleration response spectra (5% of damping ratio). In particular, we compare the pseudo-acceleration response spectra for horizontal directions with the EC8 design spectrum and the new Italian building code (NTC08). In order to understand the characteristics of the ground motions induced by L’Aquila earthquake, we also study the source-related effects and site response of the strong motion stations that recorded the seismic sequence. A novel method is used for the analysis of inter-station and site-specific H/V spectral ratios for the main event and for 12 aftershocks.     

Frequency variation in site response as observed from strong motion data of the L��Aquila (2009) seismic sequence

Previous works based mainly on strong-motion recordings of large Japanese earthquakes showed that site amplification and soil fundamental frequency could vary over long and short time scales. These phenomena were attributed to non-linear soil behaviour: the starting fundamental frequency and amplification were both instantaneously decreasing and then recovering for a time varying from few seconds to several months. The recent April 6, 2009 earthquake (M _W 6.3), occurred in the L’Aquila district (central Italy), gave us the possibility to test hypotheses on time variation of amplification function and soil fundamental frequency, thanks to the recordings provided by a pre-existing strong-motion array and by a large number of temporary stations. We investigated the intra- and inter-event soil frequency variations through different spectral analyses, including time-frequency spectral ratios and S-Transform (Stockwell et al. in IEEE Trans Signal Process 44:998–1001, 1996). Finally, analyses on noise recordings were performed, in order to study the soil behaviour in linear conditions. The results provided puzzling evidences. Concerning the long time scale, little variation was observed at the permanent stations of the Aterno Valley array. As for the short time-scale variation, the evidence was often contrasting, with some station showing a time-varying behavior, while others did not change their frequency with respect to the one evaluated from noise measurements. Even when a time-varying fundamental frequency was observed, it was difficult to attribute it to a classical, softening non-linear behaviour. Even for the strongest recorded shocks, with peak ground acceleration reaching 0.7 g, variations in frequency and amplitude seems not relevant from building design standpoint. The only exception seems to be the site named AQV, where the analyses evidence a fundamental frequency of the soil shifting from 3 Hz to about 1.5 Hz during the mainshock.     

Preliminary report of the September 5, 2022 MS 6.8 Luding earthquake,Sichuan, China

The 2022 M_S 6.8 Luding earthquake is the strongest earthquake in Sichuan Province, Western China, since the 2017 M_S 7.0 Jiuzhaigou earthquake. It occurred on the Moxi fault in the southeastern segment of the Xianshuihe fault, a tectonically active and mountainous region with severe secondary earthquake disasters. To better understand the seismogenic mechanism and provide scientific support for future hazard mitigation, we summarize the preliminary results of the Luding earthquake, including seismotectonic background, seismicity and mainshock source characteristics and aftershock properties, and direct and secondary damage associated with the mainshock. The peak ground displacements in the NS and EW directions observed by the nearest GNSS station SCCM are ～35 mm and ～55 mm, respectively, resulting in the maximum coseismic dislocation of 20 mm along the NWW direction, which is consistent with the sinistral slip on the Xianshuihe fault. Back-projection of teleseismic P waves suggest that the mainshock rupture propagated toward south-southeast. The seismic intensity of the mainshock estimated from the back-projection results indicates a Mercalli scale of VIII or above near the ruptured area, consistent with the results from instrumental measurements and field surveys. Numerous aftershocks were reported, with the largest being M_S 4.5. Aftershock locations (up to September 18, 2022) exhibit 3 clusters spanning an area of 100 km long and 30 km wide. The magnitude and rate of aftershocks decreased as expected, and the depths became shallower with time. The mainshock and two aftershocks show left-lateral strike-slip focal mechanisms. For the aftershock sequence, the b-value from the Gutenberg-Richter frequency-magnitude relationship, h-value, and p-value for Omori’s law for aftershock decay are 0.81, 1.4, and 1.21, respectively, indicating that this is a typical mainshock-aftershock sequence. The low b-value implies high background stress in the hypocenter region. Analysis from remote sensing satellite images and UAV data shows that the distribution of earthquake-triggered landslides was consistent with the aftershock area. Numerous small-size landslides with limited volumes were revealed, which damaged or buried the roads and severely hindered the rescue process.     

Separation of source and site effects by generalized inversion technique using the aftershock recordings of the 2009 L��Aquila earthquake

We exploit S-wave spectral amplitudes from 112 aftershocks (3.0 ≤ M_L ≤ 5.3) of the L’Aquila 2009 seismic sequence recorded at 23 temporary stations in the epicentral area to estimate the source parameters of these events, the seismic attenuation characteristics and the site amplification effects at the recording sites. The spectral attenuation curves exhibit a very fast decay in the first few kilometers that could be attributed to the large attenuation of waves traveling trough the highly heterogeneous and fractured crust in the fault zone of the L’Aquila mainshock. The S-waves total attenuation in the first 30 km can be parameterized by a quality factor Q_S(f) = 23f ^0.58 obtained by fixing the geometrical spreading to 1/R. The source spectra can be satisfactorily modeled using the omega-square model that provides stress drops between 0.3 and 60 MPa with a mean value of 3.3±2.8 MPa. The site responses show a large variability over the study area and significant amplification peaks are visible in the frequency range from 1 to more than 10 Hz. Finally, the vertical component of the motion is amplified at a number of sites where, as a consequence, the horizontal-to-vertical spectral ratios (HVSR) method fails in detecting the amplitude levels and in few cases the resonance frequencies.     

Predicted ground motion after the L��Aquila 2009 earthquake (Italy, M w 6.3): input spectra for seismic microzoning

After the April 6th 2009 L’Aquila earthquake (M _w 6.3), where 306 people died and a further 60,000 were displaced, seismic microzoning investigations have been carried out for towns affected by a macroseismic intensity equal to or greater than 7 MCS. Based upon seismotectonic data, historical seismicity and strong motion records, we defined input spectra to be used in the numerical simulations of seismic microzoning in four key municipalities, including the town of L’Aquila. We adopted two main approaches: uniform hazard response spectra are obtained by a probabilistic seismic hazard assessment introducing some time-dependency for individual faults on the study area; a deterministic design spectrum is computed from magnitude/distance pairs extracted by a stationary probabilistic analysis of historical intensities. The uniform hazard spectrum of the present Italian building code represents the third, less restrictive, response spectrum to be used for the numerical simulations in seismic microzoning. Strong motions recordings of the main shock of the L’Aquila sequence enlighten the critical role played by both the local response and distances metric for sites located above a seismogenic fault; however, these time-histories are compatible with the uncertainties of a deterministic utilization of ground motion predictive equations. As recordings at very near field are rare, they cannot be neglected while defining the seismic input. Disaggregation on the non-Possonian seismotectonic analysis and on the stationary site-intensity estimates reach very similar results in magnitude-distance pairs identification; we interpret this convergence as a validation of the geology-based model by historical observations.     

Ionospheric precursor of a destructive earthquake that occurred on April 6, 2009 at L’Aquila (Italy)

On the basis of the 15-min data from a series of ground-based vertical ionospheric sounding stations, a study of variations of the foF2 critical frequency before the strong earthquake (M = 6.3) that occurred on April 6, 2009 at L’Aquila (Italy) was carried out. The earthquake epicenter was located 85 km north-eastward from Rome. Approximately 20 h prior to the earthquake, a well-pronounced statistically significant effect of foF2 increase relative to the average background for magnetically quiet days was observed for almost 1.5 h at the Rome ionospheric station. In this case, at control stations distanced from the earthquake epicenter, no statistically significant deviations of foF2 from the background values were detected during the same observations period. This fact provides grounds for consideration of the foF2 increase observed at Rome station as a possible ionospheric precursor of this earthquake.     

Analysis of the liquefaction phenomena in the village of Vittorito (L��Aquila)

This paper describes a case-history of liquefaction occurred near the village of Vittorito after the April 6, 2009 L’Aquila earthquake (moment magnitude M_w = 6.3), approximately 45 km far from the epicentre. In the document, first, an estimation of the seismic motion in the area has been made. Thereafter, the performed geotechnical investigation is described, followed by the application of some fast assessment criteria for the occurrence of liquefaction, recently proposed by the new Italian Building Code. A careful assessment of all the parameters involved in conventional Seed and Idriss (1971) liquefaction analyses is considered. The cyclic resistance ratio CRR is evaluated by cone penetration tests CPT and by in situ seismic dilatometer tests SMDT; in the latter case CRR is evaluated by different empirical correlations with shear wave velocity V_s and horizontal stress index K_D. Analytical data confirmed the observed occurrence of the liquefaction in Vittorito, even if the acceleration field in the area, produced by the L’Aquila earthquake, was very low.     

Scaling characteristics of SEGMA magnetic field data around the Mw 6.3 Aquila earthquake

We apply detrended fluctuation analysis (DFA) on fluxgate and search-coil data in ULF range (scales 10–90 s or 0.1–0.011 Hz) for the months January–April 2009 available from the South European GeoMagnetic Array stations: Castello Tesino (CST), Ranchio (RNC), and L’Aquila (AQU) in Italy; Nagycenk (NCK) in Hungary; and Panagyuriste (PAG) in Bulgaria. DFA is a data processing method that allows for the detection of scaling behaviors in observational time series even in the presence of non-stationarities. The H and Z magnetic field components at night hours (00-03 UT, 01–04 LT) and their variations at the stations CST, AQU, NCK, and PAG have been examined and their scaling characteristics are analyzed depending on geomagnetic and local conditions. As expected, the scaling exponents are found to increase when the K _p index increases, indicating a good correlation with geomagnetic activity. The scaling exponent reveals also local changes (at L’Aquila), which include an increase for the Z (vertical) component, followed by a considerable decrease for the X (horizontal) component in the midst of February 2009. Attempts are made to explain this unique feature with artificial and/or natural sources including the enhanced earthquake activity in the months January–April 2009 at the L’Aquila district.     

Recorded ground motion and site effects evaluation for the April 6, 2009 L��Aquila earthquake

On April 6, 2009 a M_L = 5.8 earthquake hit the city of L’Aquila on the Apennine chain in central Italy. Notwithstanding the moderate-size event the L’Aquila city and several small villages along the Aterno river valley suffered severe damage, because of the unusual strong motions, mainly due to proximity to the fault (estimated hypocentral depth of about 10 km). In this paper the main features of the recorded motion are discussed. Four accelerometric stations were located within the surface projection of the fault and recorded peak values ranging from 0.4 to 0.6 g. The recorded motions were characterised by short durations and high peak accelerations both in the horizontal and vertical directions. The strong portions of vertical and horizontal motions occurred almost simultaneously due to the short travel paths of P and S waves from the fault to the ground surface near the fault area. Hence site response analyses were performed for the sites where recording stations were located. The geotechnical subsoil model was derived by boreholes, in situ dynamic tests (D-H and SDMT) and by laboratory tests (RCT). One-dimensional numerical analyses were carried out employing the well known computer code EERA. The numerical model was calibrated, in the linear equivalent range, by comparing numerical results with the horizontal acceleration recorded components.     

Using the match-and-locate method to characterize foreshocks of the July 2019 MW6.4 Ridgecrest,California earthquake

The July 2019 M_W6.4 Ridgecrest, California earthquake and its distinct foreshocks were well recorded by local and regional stations, providing a great opportunity to characterize its foreshocks and investigate the nucleation mechanisms of the mainshock. In this study, we utilized the match-and-locate (M&L) method to build a high-precision foreshock catalog for this M_W6.4 earthquake. Compared with the sequential location methods (matched-filter + cross-correlation-based hypoDD), our new catalog contains more events with higher location accuracy. The M_W6.4 mainshock was preceded by 40 foreshocks within ～2 h (on July 4, 2019 from 15:35:29 to 17:32:52, UTC). Their spatiotemporal distribution revealed a complex seismogenic structure consisting of multiple fault strands, which were connected as a throughgoing fault by later foreshocks and eventually accommodated the 2019 M_W6.4 mainshock. To better understand the nucleation mechanism, we determined the rupture dimension of the largest M_L4.0 foreshock by calculating its initial rupture and centroid points using the M&L method. By estimating Coulomb stress change we suggested that the majority of foreshocks following the M_L4.0 event and M_W6.4 mainshock occurred within regions of increasing Coulomb stress, indicating that they were triggered by stress transfer. The nucleation process before the M_L4.0 event remains unclear due to the insufficient sampling rate of waveforms and small magnitude of events. Thus, our study demonstrates that the M&L method has superior detection and location ability, showing potential for studies that require high-precision location (e.g., earthquake nucleation).     

Quick moment tensor solution for 6 April 2009, L'Aquila,Italy, earthquake

A M_W6.4 earthquake occurred in L'Aquila, central Italy at 1:32:42 (UTC), April 6, 2009. We quickly obtained the moment tensor solution of the earthquake by inverting the P waveforms of broadband recordings from the global seismographic network (GSN) stations using the quick technique of moment tensor inversion, and further inferred that the nodal plane of strike 132°, dip 53° and rake ?103° is the seismogenic fault.     

Seismic moment tensor inversion using 3D velocity model and its application to the 2013 Lushan earthquake sequence

Source inversion of small-magnitude events such as aftershocks or mine collapses requires use of relatively high frequency seismic waveforms which are strongly affected by small-scale heterogeneities in the crust. In this study, we developed a new inversion method called gCAP3D for determining general moment tensor of a seismic source using Green's functions of 3D models. It inherits the advantageous features of the “Cut-and-Paste” (CAP) method to break a full seismogram into the Pnl and surface-wave segments and to allow time shift between observed and predicted waveforms. It uses grid search for 5 source parameters (relative strengths of the isotropic and compensated-linear-vector-dipole components and the strike, dip, and rake of the double-couple component) that minimize the waveform misfit. The scalar moment is estimated using the ratio of L₂ norms of the data and synthetics. Focal depth can also be determined by repeating the inversion at different depths. We applied gCAP3D to the 2013 M_s 7.0 Lushan earthquake and its aftershocks using a 3D crustal-upper mantle velocity model derived from ambient noise tomography in the region. We first relocated the events using the double-difference method. We then used the finite-differences method and reciprocity principle to calculate Green's functions of the 3D model for 20 permanent broadband seismic stations within 200 km from the source region. We obtained moment tensors of the mainshock and 74 aftershocks ranging from M_w 5.2 to 3.4. The results show that the Lushan earthquake is a reverse faulting at a depth of 13–15 km on a plane dipping 40–47° to N46° W. Most of the aftershocks occurred off the main rupture plane and have similar focal mechanisms to the mainshock's, except in the proximity of the mainshock where the aftershocks' focal mechanisms display some variations.     

High-precision relocation of the aftershock sequence of the January 8, 2022, MS6.9 Menyuan earthquake

The 2022 Menyuan M_S6.9 earthquake, which occurred on January 8, is the most destructive earthquake to occur near the Lenglongling (LLL) fault since the 2016 Menyuan M_S6.4 earthquake. We relocated the mainshock and aftershocks with phase arrival time observations for three days after the mainshock from the Qinghai Seismic Network using the double-difference method. The total length and width of the aftershock sequence are approximately 32 km and 5 km, respectively, and the aftershocks are mainly concentrated at a depth of 7–12 km. The relocated sequence can be divided into 18 km west and 13 km east segments with a boundary approximately 5 km east of the mainshock, where aftershocks are sparse. The east and west fault structures revealed by aftershock locations differ significantly. The west fault strikes EW and inclines to the south at a 71º–90º angle, whereas the east fault strikes 133º and has a smaller dip angle. Elastic strain accumulates at conjunctions of faults with different slip rates where it is prone to large earthquakes. Based on surface traces of faults, the distribution of relocated earthquake sequence and surface ruptures, the mainshock was determined to have occurred at the conjunction of the Tuolaishan (TLS) fault and LLL fault, and the west and east segments of the aftershock sequence were on the TLS fault and LLL fault, respectively. Aftershocks migrate in the early and late stages of the earthquake sequence. In the first 1.5 h after the mainshock, aftershocks expand westward from the mainshock. In the late stage, seismicity on the northeast side of the east fault is higher than that in other regions. The migration rate of the west segment of the aftershock sequence is approximately 4.5 km/decade and the afterslip may exist in the source region.     

Seismic relocation,focal mechanism and crustal seismic anisotropy associated with the 2010 Yushu MS7.1 earthquake and its aftershocks

The 2010 Yushu M_S7.1 earthquake occurred in Ganzi-Yushu fault, which is the south boundary of Bayan Har block. In this study, by using double difference algorithm, the locations of mainshock (33.13°N, 96.59°E, focal depth 10.22 km) and more than 600 aftershocks were obtained. The focal mechanisms of the mainshock and some aftershocks with M_S>3.5 were estimated by jointly using broadband velocity waveforms from Global Seismic Network (GSN) and Qinghai Seismic Network as well. The focal mechanisms and relocation show that the strike of the fault plane is about 125° (WNW-ESE), and the mainshock is left-laterally strikeslip. The parameters of shear-wave splitting were obtained at seismic stations of YUS and L6304 by systematic analysis method of shear-wave splitting (SAM) method. Based on the parameters of shear-wave splitting and focal mechanism, the characteristics of stress field in seismic source zone were analyzed. The directions of polarization at stations YUS and L6304 are different. It is concluded that after the mainshock and the M_S6.3 aftershock on April 14, the stress-field was changed.     

A summary of seismic activities in and around China in 2021

In this article, we review the general characteristics of seismicity in and around China and the overall statistics of earthquake damage in 2021, focusing on several significant events and related scientific topics. Among them, the largest event is the M_S 7.4 Madoi earthquake in Qinghai Province, northwest China. The event marks another M_S ?≥ ?7 earthquake occurring near the boundary of the Bayan Har Block that has ended a remarkable quiescence of the M_S ?≥ ?7 earthquakes within the Chinese mainland. In addition, the M_S 6.4 Yangbi earthquake in Yunnan Province, southwest China draws the most attention because of its abundant foreshocks, which are well recorded by the densely distributed seismic stations in the surrounding regions. Regarding this event, we review several recent publications focusing on the Gutenberg-Richter b-value change and the physical mechanism of foreshocks associated with this sequence. The M_S 6.0 Luxian earthquake in Sichuan Province, southwest China has caused serious damage with a relatively low magnitude, partly because the focal depth of the mainshock is relatively shallow (3.5 ?km). It is another strong earthquake occurring within the southeast Sichuan basin with low historical seismicity yet has increased significantly since 2015, probably due to shale gas development and associated hydraulic fracturing.     

Segmented seismicity of the M w 6.2 Baladeh earthquake sequence (Alborz Mountains, Iran) revealed from regional moment tensors

The M _w 6.2 Baladeh earthquake occurred on 28 May 2004 in the Alborz Mountains, northern Iran. This earthquake was the first strong shock in this intracontinental orogen for which digital regional broadband data are available. The Baladeh event provides a rare opportunity to study fault geometry and ongoing deformation processes using modern seismological methods. A joint inversion for hypocentres and a velocity model plus a surface-wave group dispersion curve analysis were used to obtain an adapted velocity model, customised for mid- and long-period waveform modelling. Based on the new velocity model, regional waveform data of the mainshock and larger aftershocks (M _w ?≥3.3) were inverted for moment tensors. For the Baladeh mainshock, this included inversion for kinematic parameters. All analysed earthquakes show dominant thrust mechanisms at depths between 14 and 26 km, with NW–SE striking fault planes. The mainshock ruptured a 28° south-dipping area of 24 × 21 km along a north-easterly direction. The rupture plane of the mainshock does not coincide with the aftershock distribution, neither in map view nor with respect to depth. The considered aftershocks form two main clusters. The eastern cluster is associated with the mainshock. The western cluster does not appear to be connected with the rupture plane of the mainshock but, instead, indicates a second activated fault plane dipping at 85° towards the north.     

Comparison of two earthquake early warning location methods

According to earthquake catalog records of Fujian Seismic Network, the T _now method and the four-station continuous location method put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations in each earthquake. It shows that the four-station continuous location method can locate more seismic events than the T _now method. By analyzing the results, it is concluded that the reason for this is that the T _now method makes use of information from stations without being triggered, while some stations failed to be reflected in earthquake catalog because of discontinuous records or unclear records of seismic phases. For seismic events whose location results can be given, there is no obvious difference in location results of the two methods and positioning deviation of most seismic events is also not significant. For earthquakes outside the network, the positioning deviation may amplify as the epicentral distance enlarges, which may relate to the situation that the seismic stations are centered on one side of epicenter and the opening angle between seismic stations used for location and epicenter is small.     

Source parameter and rupture process of the MW6.3 early strong aftershock immediately following the 2016 MW7.8 Kaikoura earthquake (New Zealand)

The 2016 M_W7.8 Kaikoura (New Zealand) earthquake was the most complex event ever instrumentally recorded and geologically investigated, as it ruptured on more than 12 fault segments of various geometries. To study the mainshock rupture characteristics, geodetic methods like InSAR and GPS play an essential role in providing satisfactory spatial resolution. However, early strong aftershocks may cause extra ground deformation which bias the mainshock rupture inversion result. In this paper, we will focus on studying the M_W 6.3 aftershock, which is the only M6+ thrust slip aftershock that occurred only 30 minutes after the Kaikoura mainshock. We will relocate the hypocenter of this event using the hypo2000 method, make the finite fault model (FFM) inversion for the detailed rupture processes and calculate the synthetic surface displacement to compare with the observed GPS data and figure out its influence on the mainshock study. Although we are not able to resolve the real ruptured fault of this event because of limited observation data, we infer that it is a west-ward dipping event of oblique slip mechanism, consistent with the subfault geometries of the Kaikoura mainshock. According to the inverted FFM, this event can generate 10–20 cm ground surface displacement and affect the ground displacement observation at nearby GPS stations.