Search

Three-dimensional coseismic surface deformation fields are important for quantifying the geometric and kinematic characteristics of earthquake rupture faults. However, traditional geodetic techniques are constrained by intrinsic limitations: Interferometric synthetic aperture radar (InSAR) can only extract far-field deformation fields owing to incoherence; global navigation satellite systems (GNSSs) can only acquire displacement at discrete points. The recently developed optical pixel correlation technique, which is based on high-resolution remote sensing images, can acquire near-field coseismic horizontal deformation. In this study, InSAR line-of-sight (LOS) and azimuth direction far-field deformation, horizontal near-field deformation determined using optical pixel correlation based on pre- and post-earthquake GaoFen (GF)-2/7 images, and vertical deformation determined by differencing pre- and post-earthquake GF-7 digital elevation models (DEMs) were combined to comprehensively provide the three-dimensional deformation field of the 2022 Mw 6.6 Menyuan earthquake. The results show that the near-field deformation field calculated by optical pixel correlation quantified displacements distributed over the rupture fault zone, which were not available from the InSAR deformation maps. We identified significant vertical displacements of ~1–1.5 m at a bend region, which were induced by local compressive stress. The maximum uplift (>2.0 m) occurred near the epicenter, on the southern sides of the main and secondary faults along the middle segment of the ruptured Lenglongling fault. In addition, surface two-dimensional strain derived from the displacement maps calculated by optical pixel correlation revealed high strain concentration on the rupture fault zone. The method described herein provides a new tool for a better understanding of the characteristics of coseismic surface deformation and rupture patterns of faults. [ABSTRACT FROM AUTHOR]

Recurrence of large-scale ancient landslides often manifests as catastrophic events, so it is important to assess their stability by monitoring the deformation process. In this study, the Da'ao ancient landslide along the Jinsha River in China was investigated and analyzed. Field investigation showed that the landslide body is being experienced deforming with many new failures. Three sliding events can be identified according to different main scarps. The Permanent Scatterer Synthetic Aperture Radar Interferometry (PS-InSAR) data further revealed that the most obvious deforming area is in the middle and toe part of the sliding body with the largest cumulative displacement of 55 mm during 2017 to 2021. Additionally, the relatively stable zone and the accelerating deforming zone are divided by displacement rate. This indicates that the Da'ao landslide is becoming unstable. At the same time, the Particle Flow Code (PFC) simulation verified that the risk of landslide recurrence is likely when subjected to strong earthquakes or heavy rainfall, and the resulting sidings could dam the Jinsha River. In addition, the analysis results suggested that the displacement rates were approximately 5 mm/year and 13.3 mm/year before and after the Ms 5.2 earthquake, respectively, which occurred on 04/21/2019. The deformation changed with seasonal variation during 2017–2021 and lagged behind precipitation, while the earthquake can increase the amount and rate of displacement. The results presented here can help us to understand how landslides recur over time and their deformation and failure mechanisms, which is useful for disaster prevention and mitigation. [ABSTRACT FROM AUTHOR]

The coseismic deformation field on both sides of the fault, especially the distribution and change characteristics of near-field deformation, not only provides important constraints for the fine inversion of the slip distribution model but also serves as an important basis for the anti-disruption defense of the cross-fault linear engineering facilities. In this paper, we used Sentinel-1 satellite data to obtain the coseismic deformation field of the Maduo earthquake by using InSAR and offset techniques. We quantitatively compared the coseismic displacement of the three types of data: InSAR, offset, and optical images. The results show that optical images and offset provided more robust near-fault (<2 km) deformation insights than InSAR, which exhibited irregular deformation patterns due to incoherence near the fault. The maximum relative displacements for InSAR and offset observations are ~2.8 m and 4 m, respectively. Then we tested various fault slip models with different data constraints, revealing that a combined inversion of GPS, InSAR, and offset data offers superior constraints on slip distribution. This integrative approach effectively captured both shallow and deep fault slip, particularly near the fault zone. The eastern branch fault model, jointly constrained by GPS, InSAR, and offset data, is the optimal coseismic slip distribution model for the Maduo earthquake, and the maximum slip is 5.55 m. [ABSTRACT FROM AUTHOR]

Understanding the kinematic characteristics of the Longitudinal Valley Fault Zone (LVFZ) can help us to better understand the evolution of orogens. The 2022 Mw7.0 Taitung earthquake that occurred in Taiwan provides us with a good opportunity to understand the motion characteristics of the Central Range Fault (CRF) and the strain partitioning pattern within the Longitudinal Valley Fault (LVF). We obtained the coseismic displacement and slip distribution of the 2022 Taiwan earthquake based on the strong-motion and GPS data available. The causative fault of this earthquake is the west-dipping Central Range Fault, which is buried beneath the western boundary of the LVF. The coseismic displacement field exhibits a quadrant distribution pattern, indicating a left-lateral strike-slip mechanism with a maximum displacement exceeding 1.25 m. The joint inversion results show that the size of the main asperity is 40 km × 20 km, and the maximum slip amount of 2.6 m is located at a depth of 10 km, equivalent to an earthquake of Mw7.04. The LVFZ is composed of LVF and CRF, which accommodates nearly half of the oblique convergence rate between the Philippine Sea Plate and the Eurasian Plate. There is a phenomenon of strain partitioning in the southern segment of the Longitudinal Valley Fault Zone. The Central Mountain Range Fault is primarily responsible for accommodating strike-slip motion, while the Longitudinal Valley Fault is mainly responsible for accommodating thrust motion. [ABSTRACT FROM AUTHOR]

Strain accumulation in foreland zones is important for understanding mountain‐building processes and seismic hazards. Material heterogeneity in this border zone, particularly the contrast between sedimentary cover and basement, affects strain accommodation and fault behavior. One manifestation of this effect is depth separation of Interferometric Synthetic Aperture Radar (InSAR)‐derived slip models and corresponding aftershock clusters; hypotheses for this vertical separation remain controversial. In this study, we investigated strain accumulation in the Kalpin fold‐and‐thrust belt (KFTB) of southwestern Tian Shan, China, using the integration of InSAR measurements, teleseismic body‐waves, near‐field strong motion data (SM), and relocated aftershocks of the 2020 Mw 6.0 Kalpin earthquake. The SM modeling and analysis of post‐seismic InSAR time series are performed for the first time for this earthquake. Our results confirm a vertical separation of the mainshock and aftershock cluster, with the former on a weak décollement and the latter on faults within the basement. InSAR time series analysis shows that post‐seismic deformation was dominated by afterslip on a splay fault directly above the ruptured décollement. Static stress transfer of co‐seismic rupture cannot explain these observations. We speculate that fluid flow and high pore pressure along pre‐existing fault planes may have reduced the fault strength and been involved in the evolution of this aftershock cluster. We conclude that compressive strain in the KFTB is accommodated by a mixture of thin‐ and thick‐skin faulting and seismic deformation across the entire crustal thickness. Specifically, we suggest that the observed shortening is partly accommodated by infrequent large earthquakes on the weak décollement. Plain Language Summary: In foreland areas, strain can be accommodated by earthquakes, aseismic slip, or both. The contrasting material properties of the sediment and basement certainly affects this process, for example, the depth separation of different fault behavior. However, it is unclear how the fault behavior varies in depth from basement to weak sediment throughout different stages of the seismic cycle. 2020 Mw 6.0 Kalpin earthquake occurred in Kalpin fold‐and‐thrust belt, southwestern Tian Shan provides some insight on this matter. Here, we combined Interferometric Synthetic Aperture Radar, aftershock relocation, teleseismic body‐waves, and strong motion data to investigate the strain release behavior during and immediately after the earthquake. We found that strain release behavior differs between sediment and basement during and immediately after the earthquake. The mainshock ruptured the weak décollement and the aftershocks were primarily rooted in the basement, while the transient afterslip was mainly distributed in the sediments. Two findings are worth noting: (a) the afterslip occurred on a fault, whose stress was partially unloaded during the earthquake, suggesting the fault strength was reduced, and (b) the earthquake rupture of the décollement highlights how destructive earthquakes can initiate on the weak décollement. Key Points: The moderate 2020 Mw 6.0 Kalpin earthquake was initiated on a weak décollement in southwest Tian ShanSeismological data and Interferometric Synthetic Aperture Radar postseismic time series analysis confirm a vertical separation of the mainshock and aftershocksCombination of thin‐skin and thick‐skin deformation processes throughout the different stages of the seismic cycle [ABSTRACT FROM AUTHOR]

We obtained the rupture process and slip distribution of the 2022 Mw6.6 Menyuan earthquake by jointly inverting accelerogram data and InSAR measurements. The near-field InSAR measurements provide good constraints on the shallow slip distributions (<6 km). The accelerogram data enable us to better resolve the deeper coseismic slip (>6 km). The combination of two types of data provided improved constrains on slip distribution of the 2022 Menyuan earthquake. The results from joint inversion of InSAR and accelerogram data reveal a 26-km-long rupture length, which roughly agrees with the mapped length from the optically identified surface rupture trace and the InSAR deformation field. We imaged a major asperity with a dimension of 14 × 6 km at 4 km depth updip of the hypocenter. The maximum slip is estimated to be 3.8 m at 4 km depth. The duration of the 2022 Menyuan earthquake is ~14 s, and 90% of the seismic moment is released in the first 10 s. The total seismic moment is estimated to be 1.31 × 1 × 1019 N·m, equivalent to a moment magnitude of Mw6.7. Our results highlight that the moderate but shallow rupture during the 2022 Menyuan earthquake could intensify the seismic damage on the surface, confirmed by field investigations. [ABSTRACT FROM AUTHOR]

We derived the joint slip models of the three major events in the 2016 Central Italy earthquake sequence by inverting strong-motion and InSAR datasets. b-values and the historic earthquake scarp offset were also investigated after processing the earthquake catalog and near-field digital elevation model data. The three major events gradually released seismic moments of 1.6 × 1018 Nm (Mw 6.1), 1.5 × 1018 Nm (Mw 6.1), and 1.1 × 1019 Nm (Mw 6.7), respectively. All the ruptures exhibit both updip and along-strike directivity, but differ in the along-strike propagation direction. The high b-value found beneath three mainshock hypocenters suggests possible fluid intrusions, explaining the cascading earthquake behavior. The cumulative surface scarp from past earthquakes shows rupturing features that are consistent with the 2016 earthquake sequence, suggesting a characteristic fault behavior. Under the assumption of the Gutenberg–Richter law, the slip budget closure test gives a maximum magnitude of Mw 6.7 and implies the seismic hazard from the largest event has been released in this sequence. [ABSTRACT FROM AUTHOR]

The challenge of ruling out potential rupture nodal planes with opposite dip orientations during interferometric synthetic aperture radar (InSAR)-based kinematic inversions has been widely reported. Typically, slip on two or more different fault planes can match the surface deformation measurements equally well. The ambiguous choice of the nodal plane for the InSAR-based models was thought to be caused by InSAR's 1-D measurement and polar orbiting direction, leading to its poor sensitivity to north–south crustal motion. Through synthetic experiments and simulations, this paper quantitatively demonstrates the main reason of the ambiguous InSAR-based models, which confuse researchers in the small-to-moderate thrust earthquake cases investigation. We propose the inherent 1-D measurement is not the principle cause of the fault plane ambiguity, since models derived from the same InSAR data predict similar, but not identical, 3-D deformation patterns. They key to differentiating between these different models is to be able to resolve the small asymmetry in the surface deformation pattern, which may be smaller in amplitude than the typical noise levels in InSAR measurements. We investigate the fault geometry resolvability when using InSAR data with different noise levels through 'R' value. We find that the resolvability does not only rely on the InSAR noise, but also on the fault geometry itself (i.e. depth, dips angle and strike). Our result shows that it is impossible to uniquely determine the dip orientation of thrust earthquakes with M w < 6.0 and depth > 5.0 km with InSAR data at a noise level that is typical for mountain belts. This inference is independent from the specific data set (i.e. interferogram or time-series) and allows one to assess if one can expect to be able to resolve the correct fault plane at all. [ABSTRACT FROM AUTHOR]

We investigate the coseismic and post-seismic deformation due to the 6 February 2018 M w 6.4 Hualien earthquake to gain improved insights into the fault geometries and complex regional tectonics in this structural transition zone. We generate coseismic deformation fields using ascending and descending Sentinel-1A/B InSAR data and GPS data. Analysis of the aftershocks and InSAR measurements reveal complex multifault rupture during this event. We compare two fault model joint inversions of SAR, GPS and teleseismic body waves data to illuminate the involved seismogenic faults, coseismic slip distributions and rupture processes. Our preferred fault model suggests that both well-known active faults, the dominantly left-lateral Milun and Lingding faults, and previously unrecognized oblique-reverse west-dipping and north-dipping detachment faults, ruptured during this event. The maximum slip of ∼1.6 m occurred on the Milun fault at a depth of ∼2–5 km. We compute post-seismic displacement time series using the persistent scatterer method. The post-seismic range-change fields reveal large surface displacements mainly in the near-field of the Milun fault. Kinematic inversions constrained by cumulative InSAR displacements along two tracks indicate that the afterslip occurred on the Milun and Lingding faults and the west-dipping fault just to the east. The maximum cumulative afterslip of 0.4–0.6 m occurred along the Milun fault within ∼7 months of the main shock. The main shock-induced static Coulomb stress changes may have played an important role in driving the afterslip adjacent to coseismic high-slip zones on the Milun, Lingding and west-dipping faults. [ABSTRACT FROM AUTHOR]

Pre-seismic thermal anomalies have been widely reported, however, there is still no unified understanding of their causes or physical mechanisms. To clarify whether subsurface fluid convection and rock deformation directly related to seismic fault stress release can form thermal anomalies before an earthquake, we carried out thermo-hydro-mechanical (THM) coupled modeling of the Wenchuan earthquake, and considered the influence of maximum horizontal principal stress direction (MPSD) and fault dip. The results show that convective thermal flux (CT) anomalies of subsurface fluid and the deformation induced thermal flux (DT) of the rock mass caused by stress release of the seismic fault can explain the observed thermal anomalies along the fault. The calculated and observed thermal anomalies are comparable in spatial distribution, evolution, and magnitude. Specifically, for the Wenchuan earthquake, seismic fault stress release-induced CT and DT anomalies occur in the fault zone and in the adjacent hanging wall area. In the fault zone, CT and DT induced thermal anomalies were comparable, and superposition of the two was characterized by strong initial warming, followed by weaker warming and then short-term cooling. However, outside of the fault zone, CT induced thermal anomalies had much higher intensities than those induced by DT, and were characterized by band-like warming along the hanging wall of the fault. In addition, variations in the MPSD and fault dip cause the segmentation of thermal anomalies along the Wenchuan earthquake fault. From the southern to the northern segments of the fault, the intensity of thermal anomalies is significantly weakened with rotation of MPSD and an increase in fault dip. Our study provides evidence that seismic fault stress release is the cause of pre-seismic thermal anomalies; furthermore, the results indicate the CT and DT could be important physical mechanisms of pre-seismic thermal anomalies. [ABSTRACT FROM AUTHOR]

On January 21st, 2016, a Ms 6.4 earthquake hit Menyuan County, Qinghai province, China. The nearest known fault is the Leng Long Ling (LLL) fault which is located approximately 7 km north of the epicenter. This fault has mainly shown sinistral strike-slip movement since the late Quaternary Period. However, the focal mechanism indicates that it is a thrust earthquake, which is different from the well-known strike-slip feature of the LLL fault. In this study, we determined the focal mechanism and primary nodal plane through multi-step inversions in the frequency and time domain by using the broadband regional seismic waveforms recorded by the China Digital Seismic Network (CDSN). Our results show that the rupture duration was short, within 0–2 s after the earthquake, and the rupture expanded upwards along the fault plane. Based on these fault parameters, we then solve for variable slip distribution on the fault plane using the InSAR data. We applied a three-segment fault model to simulate the arc-shaped structure of the northern LLL fault, and obtained a detailed slip distribution on the fault plane. The inversion results show that the maximum slip is 0.43 m, and the average slip angle is 78.8°, with a magnitude of Mw 6.0 and a focal depth of 9.38 km. With the geological structure and the inversion results taken into consideration, it can be suggested that this earthquake was caused by the arc-shaped secondary fault located at the north side of the LLL fault. The secondary fault, together with the LLL fault, forms a normal flower structure. The main LLL fault extends almost vertically into the base rock and the rocks between the two faults form a bulging fault block. Therefore, we infer that this earthquake is the manifestation of a neotectonics movement, in which the bulging fault block is lifted further up under the compresso-shear action caused by the Tibetan Plateau pushing towards the northwest direction. [ABSTRACT FROM AUTHOR]

The July 3rd, 2015 Pishan, China Mw 6.4 earthquake occurred within the Hetian fold belt, a frontal thrust region between the Northwestern Tibetan Plateau and the Tarim basin. We investigate the fault geometry and the rupture process of the Pishan earthquake based on joint inversion of teleseismic body waves and InSAR measurements. Our results show that the top of the fault that ruptured in the Pishan earthquake is buried 5 ± 2 km beneath the surface and that the earthquake ruptured only a segment of a deeply-seated detachment lying under the Hetian fold belt. Our inferred coseismic slip model shows two slip asperities, located in the depth range of about 7–9 km and 10–11 km. Additionally, we resolve a slip deficit region between the slip asperities, which might indicate friction heterogeneity on the fault plane. The earthquake ruptured for ∼15 s, releasing a total seismic moment of 4.7 × 10 18 N m. We also show that the latest Pishan earthquake caused positive Coulomb stress changes on back thrust faults, potentially decreasing the time to the next earthquake on those faults. Based on our Pishan coseismic slip model, along with both the trend of the anticline system and the pattern of coseismic deformation, we envision an ongoing east-west extensional growth of the anticlines in the Western Kunlun frontal thrust region in addition to the northern underthrusting of the Northwestern Tibetan Plateau below the Tarim basin. [ABSTRACT FROM AUTHOR]

On September 16, 2015, a magnitude 8.3 earthquake struck west of Illapel, Chile. We analyzed Sentinel-1A/IW InSAR data on the descending track acquired before and after the Chile Mw8.3 earthquake of 16 September 2015. We found that the coseismic deformation field of this event consists of many semi circular fringes protruding to east in an approximately 300km long and 190km wide region. The maximum coseismic displacement is about 1.33m in LOS direction corresponding to subsidence or westward shift of the ground. We inverted the coseismic fault slip based on a small-dip single plane fault model in a homogeneous elastic half space. The inverted coseismic slip mainly concentrates at shallow depth above the hypocenter with a symmetry shape. The rupture length along strike is about 340 km with maximum slip of about 8.16m near the trench. The estimated moment is 3.126×1021 N.m(Mw8.27), the maximum depth of coseismic slip near zero appears to 50km. We also analyzed the postseismic deformation fields using four interferograms with different time intervals. The results show that postseismic deformation occurred in a narrow area of approximately 65km wide with maximum slip 11cm, and its predominant motion changes from uplift to subsidence with time. that is to say, at first, the postseismic deformation direction is opposite to that of coseismic deformation, then it tends to be consistent with coseismic deformation. It maybe indicates the differences and changes in the velocity between the Nazca oceanic plate and the South American continental plate. [ABSTRACT FROM AUTHOR]

The complete surface deformation of 2015 Mw 8.3 Illapel, Chile earthquake is obtained using SAR interferograms obtained for descending and ascending Sentinel-1 orbits. We find that the Illapel event is predominantly thrust, as expected for an earthquake on the interface between the Nazca and South America plates, with a slight right-lateral strike slip component. The maximum thrust-slip and right-lateral strike slip reach 8.3 and 1.5 m, respectively, both located at a depth of 8 km, northwest to the epicenter. The total estimated seismic moment is 3.28 × 10 N.m, corresponding to a moment magnitude Mw 8.27. In our model, the rupture breaks all the way up to the sea-floor at the trench, which is consistent with the destructive tsunami following the earthquake. We also find the slip distribution correlates closely with previous estimates of interseismic locking distribution. We argue that positive coulomb stress changes caused by the Illapel earthquake may favor earthquakes on the extensional faults in this area. Finally, based on our inferred coseismic slip model and coulomb stress calculation, we envision that the subduction interface that last slipped in the 1922 Mw 8.4 Vallenar earthquake might be near the upper end of its seismic quiescence, and the earthquake potential in this region is urgent. [ABSTRACT FROM AUTHOR]

The 2013 Mw6.8 Lushan, China earthquake occurred in the southwestern end of the Longmenshan fault zone. We jointly invert local strong motion data and geodetic measurements of coseismic surface deformation, including GPS and InSAR, to obtain a robust model of the rupture process of the 2013 Lushan earthquake. Our joint inversion best model involves the rupture of two opposing faults during the Lushan earthquake, a main fault and a secondary fault. It is only when the secondary fault is included that both the GPS and InSAR measurements are fit along with the near-field strong motion. Over 75% of the computed moment was released in slip on the main fault segment, a northwest dipping, listric thrust fault, with buried thrust and dextral strike-slip at hypocenter depths, and with only minor slip closer to the surface. The secondary fault mainly involved oblique thrust slip or pure dextral strike-slip at shallower depths, and accounts for just under 24% of the moment released in the Lushan earthquake. Coulomb stress changes of about 0.5 MPa on the secondary fault segment at the time coseismic slip initiated on that fault indicate that slip was likely triggered by the coseismic slip on the main blind thrust fault. Our coseismic slip model is consistent with a sub-horizontal and east–west to southeast–northwest trending most compressive stress. Our inferred coseismic slip model is also consistent with previous GPS derived models of strain accumulation on the Longmenshan fault system. [ABSTRACT FROM AUTHOR]

The 2019 Ridgecrest, California earthquake sequence ruptured along a complex fault system and triggered seismic and aseismic slips on intersecting faults. To characterize the surface rupture kinematics and fault slip distribution, we used optical images and Interferometric Synthetic Aperture Radar (InSAR) observations to reconstruct the displacement caused by the earthquake sequence. We further calculated curl and divergence from the north-south and east-west components, to effectively identify the surface rupture traces. The results show that the major seismogenic fault had a length of ~55 km and strike of 320 ° and consisted of five secondary faults. On the basis of the determined multiple-fault geometries, we inverted the coseismic slip distributions by InSAR measurements, which indicates that the Mw7.1 mainshock was dominated by the right-lateral strike-slip (maximum strike-slip of ~5.8 m at the depth of ~7.5 km), with a small dip-slip component (peaking at ~1.8 m) on an east-dipping fault. The Mw6.4 foreshock was dominated by the left-lateral strike-slip on a north-dipping fault. These earthquakes triggered obvious aseismic creep along the Garlock fault (117.3° W–117.5° W). These results are consistent with the rupture process of the earthquake sequence, which featured a complicated cascading rupture rather than a single continuous rupture front propagating along multiple faults. [ABSTRACT FROM AUTHOR]

We obtained high-resolution (10 m) horizontal displacement fields from pre- and post-seismic Sentinel-2 optical images of the 2018 Mw7.5 Palu earthquake using subpixel image correlation. From these, we calculated the curl, divergence, and shear strain fields from the north-south (NS) and east-west (EW) displacement fields. Our results show that the surface rupture produced by the event was distributed within the Sulawesi neck (0.0974–0.6632°S) and Palu basin (0.8835–1.4206°S), and had a variable strike of 313.0–355.2° and strike slip of 2.00–6.62 m. The NS and EW displacement fields within the Palu basin included fine-scale displacements in both the near- and far-fault, the deformation patterns included a small restraining bend (localized shortening), a distributed rupture zone, and a major releasing bend (net extension) from the curl, divergence, and shear strain. Surface rupture was dominated by left-lateral strike-slip from initiation to termination, with a localized normal slip component peaking at ~3.75 m. The characteristics and geometric variation of the ruptured fault controlled both the formation of these surface deformation patterns and sustained supershear rupture. [ABSTRACT FROM AUTHOR]

The 12 November 2017 Darbandikhan earthquake (Mw 7.3) occurred along the converence zone. Despite the extensive research on this earthquake, none of this work explained whether this earthquake rupture was limited to the thick sedimentary cover or it extends to the underlying crystalline basement rock (or both). Besides, whether this region will generate devastating earthquakes again and whether there is a one-to-one correlation between these anticlines and blind-reverse faults need further investigation. In this study, we derived the co-seismic interferograms from the Sentinel-1A/B data and successfully described the surface deformation of the main seismic zone. The fringe patterns of both the ascending and descending interferograms show that the co-seismic deformation is dominated by horizontal movements. Then, using the along- and across-track deformation fields of different orbits, we retrieved the three-dimensional deformation field, which suggests that the Darbandikhan earthquake may be a blind thrust fault close to the north–south direction. Finally, we inverted the geometrical parameters of the seismogenic fault and the slip distribution of the fault plane. The results show that the source fault has an average strike of 355.5° and a northeast dip angle of −17.5°. In addition, the Darbandikhan earthquake has an average rake of 135.5°, with the maximum slip of 4.5 m at 14.5 km depth. On the basis of the derived depth and the aftershock information provided by the Iranian Seismological Center, we inferred that this event primarily ruptured within the crystalline basement and the seismogenic fault is the Zagros Mountain Front Fault (MFF). The seismogenic region has both relatively low historical seismicity and convergent strain rate, which suggests that the vicinity of the epicenter may have absorbed the majority of the energy released by the convergence between the Arabian and the Eurasian plates and may generate Mw > 7 earthquakes again. Moreover, the Zagros front fold between the Lurestan Arc and the Kirkuk Embayment may be generated by the long-distance slippage of the uppermost sedimentary cover in response to the sudden shortening of the MFF basement. We thus conclude that the master blind thrust may control the generation of the Zagros front folding. [ABSTRACT FROM AUTHOR]

During an earthquake, crustal deformation, fluid flow, and temperature variation are coupled; however, earthquake-related land surface temperature (LST) variations remain unclear. To determine whether post-seismic fluid migration can cause changes in LST, and taking the Mw 7.3 2017 Iran earthquake as an example, we modeled surface cooling (CA) and warming (WA) areas induced by co-seismic slip and fluid migration using a thermo-hydro-mechanical (THM) coupled numerical simulation. Moreover, using nighttime LST data with 15-min resolution, the daily attenuation coefficient k of nighttime LST was extracted by attenuation function fitting, and the trend of the k time series was analyzed using the Mann–Kendall and Sen's methods. Based on the comparison of k trends between the post-seismic and 2010–2016 periods, we obtained cooling and warming trends for the modeled CA and WA. The numerical simulation and observational data show good consistency, and both indicate that fluid migration caused by crustal deformation can lead to changes in LST. The numerical simulations show that after the Iran earthquake, the surface projection area of co-seismic slip correlated with a cooling area (CA), while the surrounding area correlated with a warming area (WA). For the LST observational data, the post-seismic k trends of the calculated CA and WA are positive and negative, indicating sustained cooling and warming processes, respectively. This study provides evidence that LST variation is caused by co-seismic crustal deformation and fluid migration and reveals the coupled evolution of deformation, fluid, and temperature fields. The results provide new insights into the mechanisms of seismic thermal anomalies. [ABSTRACT FROM AUTHOR]

The 2001 Mw7.8 Kokoxili earthquake, which occurred in the north Tibetan Plateau, ruptured ~400 km of the westernmost portion of the Kunlun fault and produced significant time-dependent postseismic deformation over a large area around the rupture zone and nearby regions. To analyze the postseismic deformation features along different sections of the coseismic surface rupture, we describe the total cumulative postseismic deformation near the center of the rupture and produce velocity maps for the whole observation period and six sub-periods, using InSAR observations (ENVISAT/ASAR, 2003–2010) on five descending tracks. The results indicate that the postseismic deformation is asymmetrically distributed across the fault over a very broad area of ~300 km × 500 km. The south side of the fault exhibits larger displacements and a wider area of deformation that is steadily decaying from near-field to far-field, while the north side displays a narrow, rapidly diminishing deformation field. The maximum cumulative displacement in 2003–2010 reaches up to ~45–60 mm and the LOS peak-to-trough average velocity offset in 2003–2010 reaches ~13–16 mm/yr at ~92.5°E. The short-term postseismic velocity estimates in the six sub-periods reflect significant spatial variation and temporal differences on different sections. Motions to the south of the two ends of the rupture zone show more rapid velocity decay compared to near the main central rupture zone. The time- and distance-dependent timeseries of postseismic surface displacement reveal exponential decay in the near-field and a nearly linear trend in the far-field of the fault. [ABSTRACT FROM AUTHOR]

[Display omitted] [ABSTRACT FROM AUTHOR]

Investigating relationships between temporally- and spatially-related continental earthquakes is important for a better understanding of the crustal deformation, the mechanism of earthquake nucleation and occurrence, and the triggering effect between earthquakes. Here we utilize Global Positioning System (GPS) velocities before and after the 2008 Mw 7.9 Wenchuan earthquake to invert the fault coupling of the Longmenshan Fault (LMSF) and investigate the impact of the 2008 Mw 7.9 Wenchuan earthquake on the 2013 Mw 6.6 Lushan earthquake. The results indicate that, before the 2008 Mw 7.9 Wenchuan earthquake, fault segments were strongly coupled and locked at a depth of ~18 km along the central and northern LMSF. The seismic gap between the two earthquake rupture zones was only locked at a depth < 5 km. The southern LMSF was coupled at a depth of ~10 km. However, regions around the hypocenter of the 2013 Mw 6.6 Lushan earthquake were not coupled, with an average coupling coefficient ~0.3. After the 2008 Mw 7.9 Wenchuan earthquake, the central and northern LMSF, including part of the seismic gap, were decoupled, with an average coupling coefficient smaller than 0.2. The southern LMSF, however, was coupled to ~20 km depth. Regions around the hypocenter of the 2013 Mw 6.6 Lushan earthquake were also coupled. Moreover, by interpreting changes of the GPS velocities before and after the 2008 Mw 7.9 Wenchuan earthquake, we find that the upper crust of the eastern Tibet (i.e., the Bayan Har block), which was driven by the postseismic relaxation of the 2008 Mw 7.9 Wenchuan earthquake, thrust at an accelerating pace to the Sichuan block and result in enhanced compression and shear stress on the LMSF. Consequently, downdip coupling of the fault, together with the rapid accumulation of the elastic strain, lead to the occurrence of the 2013 Mw 6.6 Lushan earthquake. Finally, the quantity analysis on the seismic moment accumulated and released along the southern LMSF show that the 2013 Mw 6.6 Lushan earthquake should be defined as a “delayed” aftershock of the 2008 Mw 7.9 Wenchuan earthquake. The seismic risk is low along the seismic gap, but high on the unruptured southwesternmost area of the 2013 Mw 6.6 Lushan earthquake. [ABSTRACT FROM AUTHOR]

We have generated a more than 500 km long postseismic deformation rate map and cumulative displacement time series in the central Bayankala Block of the Tibetan Plateau using ENVISAT/ASAR data from 2003 to 2010 by the π-RATE stacking algorithm. This rate map spans a period of ∼7.2 years and reveals that postseismic motion of 2001 Kokoxili earthquake exhibits a striking signal, dominating crustal deformation of the central Bayankala Block with a cross-fault magnitude ∼9–11 mm/yr in line of sight (LOS) (∼93.1°E). The southern and northern parts of the postseismic deformation field exhibit different patterns and variable magnitudes, reflecting asymmetry of the displacement distribution. Postseismic motion affects eastward extrusion of the central Bayankala Block, which reaches ∼15 km north of the Ganzi-Yushu fault. To further investigate viscoelastic relaxation of the lower crust in this area after approximately two years, E-M, E-M-M and E-S models are constructed. The result shows that the best fit viscosity for the lower crust is about 1 × 10 19  Pa·s. Comparison between cumulative displacements resolved by these three models shows that viscoelastic relaxation of the lower crust makes the most significant contribution to postseismic stress relaxation after 2001 event. [ABSTRACT FROM AUTHOR]

On 8 August 2017, a Ms7.0 earthquake stroke the city of Jiuzhaigou, Sichuan, China. The Jiuzhaigou earthquake occurred on a buried fault in the vicinity of three well-known active faults and this event offers a unique opportunity to study tectonic structures in the epicentral region and stress transferring. Here we present coseismic displacement field maps for this earthquake using descending and ascending Sentinel-1A Interferometric Synthetic Aperture Radar (InSAR) data. Deformation covered an area of approximately 50 × 50 km, with a maximum line-of-sight (LOS) displacement of ~22 cm in ascending and ~14 cm in descending observations on the west side of the source fault. Based on InSAR and Global Positioning System (GPS) measurements, both separately and jointly, we constructed a one-segment model to invert the coseismic slip distribution and dip angle of this event. Our final fault slip model suggests that slip was concentrated at an upper depth of 15 km; there was a maximum slip of ~1.3 m and the rupture was dominated by a left-lateral strike-slip motion. The inverted geodetic moment was approximately 6.75 × 10 18  Nm, corresponding to a moment magnitude of Mw6.5, consistent with seismological results. The calculated static Coulomb stress changes indicate that most aftershocks occurred in stress increasing zones caused by the mainshock rupture; the Jiuzhaigou earthquake has brought the western part of the Tazang fault 0.1–0.4 MPa closer to failure, indicating an increasing seismic hazard in this region. The Coulomb stress changes caused by the 2008 Mw7.8 Wenchuan earthquake suggest that stress loading from this event acted as a trigger for the Jiuzhaigou earthquake. [ABSTRACT FROM AUTHOR]

We invert InSAR from ascending and descending passes of ENVISAT/ASAR data to estimate the distribution of coseismic slip for the 2008 Mw 7.2 Yutian earthquake by separately using uniform finite element models (FEMs), non-uniform FEMs, and analytical elastic half-space (Okada) models. We use the Steepest Descent method and Laplacian smoothing to regularize and estimate the slip distribution for each of these different models. Fault surface ruptures interpolated from Quickbird optical images constrain the strike of the fault. The uniform FEM and Okada models assume a fault embedded in a homogeneous and elastic finite domain or half-space, respectively, while the non-uniform FEM describes a domain with spatially variable elastic properties according to the geological structure of Tibet plateau and Tarim. We find that the estimated slip distribution of the non-uniform FEM is similar to that of the uniform FEM but has a slightly larger moment (~2 %). However, more noticeable discrepancies occur between the slip distributions of the uniform FEM and Okada model, where the slip pattern estimated with the FEM looks more scattered and locates at a lower depth. We conclude that, for the case of the Yutian earthquake, the improvement brought by simulating the distribution of geologic material properties is insignificant. This suggests that the uniform geologic structure of Okada models may be sufficient for simulating large intra-continental earthquakes. [ABSTRACT FROM AUTHOR]

Three-dimensional (3-D) fault geometry and surface slip are key factors for understanding fault behaviour and earthquake mechanics. Previous studies of the 2008 Yutian earthquake occurred on the northwestern Tibetan Plateau, and are limited by a lack of detailed measurements of fault structures and near-fault motion due to the inaccessibility and extreme weather of the area and severe decorrelation of InSAR. In this paper, we use a very-high resolution (VHR) DEM of the 2008 rupture zone extracted from the Pléiades satellite stereo imagery to measure the 3-D fault geometry and on-fault slip. The maximum vertical and horizontal coseismic displacement is 3.6 m and 2.8 m respectively. Despite the overall normal faulting, near-fault reverse vertical motion was observed, ranging from 0.5 m to 1.9 m, caused by the transformation of strike-slip motion into localised shortening due to local fault geometric variations. Measurements of the near-surface dip for the three fault segments (55°, 62°, and 58°) show good consistency, different from previous studies (43°, 60°, and 52°). Joint inversion of InSAR and the near-fault measurements of surface slip also show that the fault geometry does not change much along strike. The local ratio of the cumulative to the coseismic vertical offset recorded on an old alluvial fan and an ice tongue ranges from 1.86 to 2.79, suggesting at least 1–2 historical events occurred here before the 2008 earthquake. • 3D fault geometry and slip of Yutian earthquake from very-high resolution DEM • Near-fault reverse vertical motion caused by local fault geometric variations • Measurements and inversion results show good consistency for dips along strike. [ABSTRACT FROM AUTHOR]