Normal-faulting Earthquakes Research Articles

Contemporary crustal strain rates derived from GPS measurements in southwestern China

Surface strain rates in southwestern China can provide constraints on tectonic evolution of the Tibetan Plateau. Unfortunately, different strain rate fields were obtained from different authors although almost the same GPS data were exploited. Based on updated GPS data, in this paper, we calculate the strain rates in southwestern China using the method proposed by Zhu et al. (2005, 2006). The results confirm that the strain rates in the Sichuan Basin and the South China Block are very small, and the high values are largely concentrated along the Xianshuihe-Anninghe-Xiaojiang fault system and the Sagaing Fault. Furthermore, the highest principal strain rates are located around the eastern Himalayan syntaxis (EHS), with compressive orientations perpendicular to the strike of the Main Thrust Belt. The basic characteristic of the strain rate distribution is in accordance with the tectonic structures from the geological and geophysical investigations. In particular, the two calculated extensive deformation regions basically match the locations of normal-faulting earthquakes in the Sichuan-Yunnan Rhombic Block (SYRB), with one in N-S extension trending in the northern SYRB, and the other nearly E-W orientation in the southern SYRB. We suggest that the extensive deformation and the normal faulting earthquakes in southwestern China is mainly controlled by the lower crustal channel flow with a bend direction, rather than by gravitational spreading alone, as previous authors proposed.

Intensity Prediction Equations Based on the Environmental Seismic Intensity (ESI-07) Scale: Application to Normal Fault Earthquakes

Earthquake environmental effects may significantly contribute to the damage caused by seismic events; similar to ground motion, the environmental effects are globally stronger in the vicinity and decrease moving away from the epicenter or seismogenic source. To date, a single intensity prediction equation (IPE) has been proposed in the Italian Apennines for intensity scale dealings with environmental effects: the Environmental Seismic Intensity (ESI-07). Here, we evaluate the sensitivity of the IPE with respect to input data and methodological choices and we propose IPEs with global validity for crustal normal faults. We show the strong influence of input data on the obtained attenuation investigating the 1980 Irpinia–Basilicata (Southern Italy) earthquake. We exploit a dataset of 26 earthquakes to build an IPE considering the epicentral distance. We also propose an IPE considering the distance from the fault rupture, which is derived from a dataset of 10 earthquakes. The proposed equations are valid for normal faults up to 40 km from the epicenter/fault and may flank other models predicting ground motion or damage to the built environment. Our work thus contributes to the use of the ESI-07 scale for hazard purposes.

Constraints on the 2013 Saravan intraslab earthquake using permanent GNSS, InSAR and seismic data

SUMMARY On 16 April 2013, an Mw = 7.7 earthquake struck the border of Iran and Pakistan in the central part of the Makran subduction zone with a reported depth of 80 km by USGS. This rare event in this poorly instrumented region helps to shed light on the kinematics of the subducting slab. We investigate source parameters of the Saravan intraslab normal earthquake using RADARSAT-2 SAR images in three ascending tracks, nine permanent GNSS sites and teleseismic data. The maximum coseismic displacement occurred at the SRVN GNSS station with 54.1 mm southeast horizontal and 42.7 mm upward vertical displacements. The coseismic ascending InSAR displacement maps illustrate a continuous and smooth NE-trending elliptical shape deformation pattern with a maximum of ∼29 cm of displacement away from the satellite. We use 25 broad-band teleseismic P-waveforms to estimate the focal mechanism of the main shock. A joint uniform inversion of InSAR, GNSS and teleseismic data reveals a NW-dipping SW-striking fault and a primarily normal-faulting earthquake with a minor right-lateral strike-slip component. The static slip distribution of the InSAR coseismic maps localizes variable slip at depths between 50 and 81 km with a maximum amplitude >3 m at 60–75.5 km depth, rupturing the oceanic crust of the subducted slab. The kinematic slip distribution exhibits a well-constrained slip pattern with a nucleation depth of 65 km. The source time function indicates that the earthquake reaches its maximum moment tensor release at ∼8 and ∼16 s. The NE-trend of the Saravan earthquake slip pattern, the orientation of the volcanic arc, and the distribution of the intraslab intermediate-depth normal earthquakes provide new insights into slab geometry in the central Makran subduction zone. We suggest that the slab bending at the hinge of subducting Arabian Plate is oblique along a NE–SW direction parallel to the volcanic arc rather than the shoreline or deformation front, and it is likely to be the reason for an oblique volcanic arc in the Makran subduction zone. These new constraints on the Makran slab geometry will help further studies in establishing realistic coupling maps for seismic hazard assessment.

Resonance modeling of the tsunami caused by the Aegean Sea Earthquake (Mw7.0) of October 30, 2020

The resonance of tsunami waves in semi-enclosed bays is paramount in understanding and mitigating the impact of seismic events on coastal communities. Semi-enclosed bays, characterized by their partial enclosure, can amplify the effects of incoming tsunami waves due to resonance behavior, where the natural frequencies of the bay correspond to those of the incoming waves. This resonance phenomenon can significantly increase wave height and inundation levels, posing an increased risk to nearby settlements and infrastructure. Understanding the resonance patterns in these bays is crucial for accurate hazard assessment, early warning systems, and effective disaster preparedness and response strategies. On October 30, 2020, an earthquake occurred between the Turkish Bay of Seferihisar Bay and the Greek island of Samos in the Aegean Sea. Long waves generated by the normal-faulting earthquake caused notable damage to settlements within Seferihisar Bay and the north coast of Samos Island. According to the measurements of the Syros mareograph stations, the wave heights were between 2 and 20 cm and wave periods between 9 and 20 seconds. Based on on-site survey reports conducted after the earthquake, inundation was reported in six settlements within Seferihisar Bay. However, inundation was notably higher in Sığacık and Akarca, reaching 2–3 times higher than in other locations, and the water level reached 2 m high. Given that the variance in inundation levels is attributed to resonance phenomena in Sığacık and Akarca rather than the propagation of tsunami waves, this study focused on conducting wave resonance modeling in Seferihisar Bay. The resonance modeling was performed using the RIDE wave model. Furthermore, the research has been expanded to assess the resonance patterns that might emerge in the event of an alternative earthquake or underwater landslide along the fault line responsible for the seismic event, encompassing wave periods ranging from T = 1–9 minutes and T = 20–30 minutes. Modeling results revealed that on the day of the earthquake, wave heights in Sığacık Marina and Akarca surged by 8.5 times in comparison to the wave height at the epicenter. This increase is notably higher, ranging from 2 to 2.5 times, compared to calculations made for other locations (Demircili, Altınköy, and Tepecik). Consequently, it was concluded that one of the reasons for the heightened effectiveness of inundation in Sığacık and Akarca was attributable to resonance. Moreover, supplementary investigations have indicated that waves with a period of T<9 minutes will pose higher risks for Demircili, Altınköy, Sığacık Marina, and Tepecik compared to the day of the earthquake. By comprehensively studying wave resonance in semi-enclosed bays, researchers and policymakers can better anticipate the potential impact of tsunami events and take measures to protect coastal communities, ultimately increasing resilience and reducing the loss of life and property in vulnerable regions.

Global seismic energy scaling relationships based on the type of faulting

Abstract. We derived scaling relationships for different seismic energy metrics for earthquakes around the globe with MW &gt; 6.0 from 1990 to 2022. The seismic energy estimations were derived with two methodologies, the first based on the velocity flux integration and the second based on finite-fault models. In the first case, we analyzed 3331 reported seismic energies derived by integrating far-field waveforms. In the latter methodology, we used the total moment rate functions and the approximation of the overdamped dynamics to quantify seismic energy from 231 finite-fault models (Emrt and EO, EU, respectively). Both methodologies provide compatible energy estimates. The radiated seismic energies estimated from the slip models and integration of velocity records are also compared for different types of focal mechanisms (R, reverse; R-SS, reverse–strike-slip; SS, strike-slip; SS-R, strike-slip–reverse; SS-N, strike-slip–normal; N, normal; and N-SS, normal–strike-slip), and then used to derive converting scaling relations among the different energy types. Additionally, the behavior of radiated seismic energy (ER), energy-to-moment ratio (ER/M0), and apparent stress (τα) for different rupture types at a global scale is examined by considering depth variations in mechanical properties, such as seismic velocities, rock densities, and rigidities. For this purpose, we used a 1-D global velocity model. The ER/M0 ratio is, based on statistical t tests, largest for strike-slip earthquakes, followed by normal-faulting events, with the lowest values for reverse earthquakes for hypocentral depths &lt; 90 km. Not enough data are available for statistical tests at deeper intervals except for the 90 to 120 km range, where we can satisfactorily conclude that ER/M0 for R-SS and SS-R types is larger than for N-type faulting, which also conforms to the previous assumption. In agreement with previous studies, our results exhibit a robust variation in τα with the focal mechanism. Regarding the behavior of τα with depth, our results agree with the existence of a bimodal distribution with two depth intervals where the apparent stress is maximum for normal-faulting earthquakes. At depths in the range of 180–240 km, τα for reverse earthquakes is higher than for normal-faulting events. We find the trend EU &gt; Emrt &gt; EO for all mechanism types based on statistical t tests. Finite-fault energy estimations also support focal mechanism dependence of apparent stress but only for shallow earthquakes (Z &lt; 30 km). The slip distribution population used was too small to conclude that finite-fault energy estimations support the dependence of average apparent stress on rupture type at different depth intervals.

Insight into the 1 December 2016 Mw 6.2 Juliaca Earthquake, Southern Peru, by InSAR Observations and Field Investigation

On 1 December 2016, an Mw 6.2 earthquake characterized by normal faulting occurred in the highlands of the central Andes in southern Peru, marking the region’s largest shallow event. The occurrence of the earthquake provides a significant chance to gain insight into the regional tectonic deformation and the seismogenic mechanism of the shallow normal-faulting earthquake, as well as the regional potential seismic risk. Here, we first utilize Sentinel-1A interferometric synthetic aperture radar (InSAR) data to extract the coseismic and postseismic deformation associated with this earthquake and then determine the detailed coseismic slip and postseismic afterslip distribution of this event. Coseismic modeling results indicate that the coseismic rupture is mainly characterized by normal faulting with some dextral strike-slip components. Most coseismic slip is confined to a depth range of 2–12 km, indicating an obvious slip deficit area in the shallow fault part. Further postseismic modeling reveals that the majority of afterslip is concentrated at depths of 0 to 5.4 km. The relatively shallow postseismic afterslip makes up for the coseismic slip deficit area to some extent. Through a joint analysis of the inversions, seismic data, and regional geology and geomorphology, we infer that the occurrence of this 2016 normal-faulting event is a result of regional gravitational collapse. In addition, we investigate the relationship between the 2016 earthquake and great historical earthquakes near the subduction zone of the central Andes and find that the 2016 event is likely promoted in advance by these events through our calculations of the coseismic and postseismic Coulomb stress changes. Finally, we should pay more attention to the nearby Falla Huaytacucho-Condoroma fault and the western segment of the Vilcanota Fault because of their relatively high stress loading.

Geologic and Geophysical Evidence of Repeated Normal Faulting on the Itozawa Fault—the Source of the Mw 6.6 Iwaki Earthquake Triggered by the Mw 9.0 Tohoku-Oki Megathrust Earthquake, Northeast Japan

ABSTRACT The 2011 Mw 9.0 Tohoku-Oki earthquake ruptured the 500 km long and 200 km wide convergent plate margin between the North American and Pacific plates, and changed the crustal stress field and triggered widespread seismic activity in northeast Japan. In particular, many crustal earthquakes struck the southern Fukushima area. The largest inland normal-faulting earthquake with Mw 6.6 occurred on the northwest-trending Yunodake fault and the north-northwest-trending Itozawa fault on 11 April 2011. The coseismic surface ruptures appeared along the previously identified active and presumed active faults. To investigate the near-surface structure of the Itozawa fault, we conducted ground-penetrating radar (GPR) profiling across the fault, and two drilling surveys on its hanging-wall and footwall blocks. The GPR survey line about 50 m long was located near the previous trenching site. The GPR two-way travel-time sections were collected by common-offset modes by 50 and 100 MHz GPR systems. We acquired common midpoint ensembles on both sides of the surface rupture for depth conversion of the time sections. The resulting GPR sections show detailed subsurface structures to a depth of about 7 m. We identified event horizons for the 2011, penultimate, antepenultimate, and preantepenultimate earthquakes, indicating that the Itozawa fault has ruptured repeatedly as a normal fault. The vertical displacements during these earthquake events, obtained from the GPR sections, were about 0.8, 0.3, 0.4, and 0.7 m, respectively. These observations suggest that the coseismic surface offset amount on the Itozawa fault varies greatly. The recurrence intervals of surface-rupturing earthquakes on the Itozawa fault, calculated from the previous research and this study, are much longer than those of giant interplate earthquakes along the Japan trench, such as the 2011 Tohoku-Oki earthquake. We conclude that not every giant earthquake along the Japan trench triggers a rupture of the Itozawa fault.

The mechanism of the present-day crustal deformation in southeast Tibet: from numerical modelling and geodetic observations

SUMMARY The mechanism of present-day crustal deformation in southeast Tibet remains controversial. 3-D high-precision geodetic data can provide significant clues to analyse the key driving forces. Here, we conduct a series of 3-D finite-element modelling to investigate the influences of gravitational collapse, tectonic extrusion and mid-to-lower crustal flow on crustal deformation in southeast Tibet. The numerical results show that the gravitational collapse leads to predominant N-S extension and surface subsidence in the northern region, and predominant NW-SE compression and uplift in the southern region, which can explain the normal-faulting earthquakes in the interior. The gravity-driven horizontal velocity depends on the upper crustal viscosity, while the vertical velocity is determined by mid-to-lower crustal viscosity. The eastward tectonic extrusion causes slight southeastward rotation and predominant E-W compression in the northern region but has a little effect on the deformation in the southern region. By considering the joint effects of gravitational collapse and tectonic extrusion, we simulate the crustal deformation that reconciles with present-day geodetic observations. Both the two driving forces lead to positive shear strain rates along the major fault zones, with more contributions from the tectonic extrusion of the Tibetan Plateau. Constrained by the 3-D geodetic observations, the numerical results argue against the presence of massive fast mid-to-lower crustal flow from the Tibetan Plateau. Overall, the present-day crustal deformation in southeast Tibet is jointly driven by gravitational collapse and tectonic extrusion, which play distinct roles in shaping the faulting kinematics and regional strain partitioning.

Fault geometry of M6-class normal-faulting earthquakes in the outer trench slope of Japan Trench from ocean bottom seismograph observations

Since the 2011 Mw 9.0 Tohoku-oki earthquake, intra-plate normal-faulting earthquakes, including several M7-class earthquakes, have occurred in the outer trench slope area from the trench to the outer rise along the Japan Trench. Concerns regarding large earthquakes and associated tsunamis have also arisen. Based on aftershock distributions, several outer trench slope normal-faulting earthquakes (hereinafter referred to as outer-rise earthquakes) are likely related to the rupture of multiple faults. However, few observations have clearly shown how multiple faults act during outer-rise earthquakes. During the ocean bottom seismograph (OBS) observations in the outer trench slope area of the central Japan Trench from September 2017 to July 2018, three M6-class normal-faulting earthquakes (Mw 6.2 on September 20, Mw 6.2 on October 06, and Mw 6.0 on November 12) occurred around the OBS network. The near-field OBS observations provided detailed information on hypocenter locations and focal mechanisms of the mainshocks and aftershocks, including immediately after the mainshocks. We investigated the fault configurations of normal-faulting earthquakes based on OBS observations. During the September 2017 earthquake, the mainshock ruptured high-angle normal faults with a dip angle of 65°. Off-fault aftershock activities that were not directly related to the mainshock rupture and could be explained by the stress changes caused by the mainshock were confirmed. However, hypocenter distributions and focal mechanisms of the main and aftershocks of the October and November 2017 earthquakes suggest that the mainshock ruptured multiple faults with various dipping directions, angles, and strike orientations. The complicated fault geometry should be considered a possible fault model for large outer-rise earthquakes and related tsunamis.

Coseismic and early postseismic deformation of the 2020 Nima Mw 6.4 earthquake, central Tibet, from InSAR and GNSS observations

The 2020 Mw6.4 Nima earthquake is one of the largest normal-faulting earthquakes recently occurring north of the Banggong suture zone in remote central Tibet, where geologic investigation of active faults is extremely limited. We analyze jointly InSAR and GNSS observations over 9 months after the Nima earthquake and calculate the coseismic and postseismic displacement. The optimal coseismic slip model suggests this event is the result of moderate-angle down-dip slip on a complex reversed “S-shape” three-segmented structure at fault junctions of the West Yibu-Chaka fault, the Heishi fault, and an unmapped blind fault, with a small component of left-lateral slip. The superposition of seismic waves from faults with different strikes and dips accounted for a large non-double-couple component in the long-period point-source solutions. The geodetic moment released by the mainshock is 6.4 × 1018 N⋅m, equivalent to Mw 6.42. Coseismic rupture concentrated at a depth of 4–15 km, with a peak slip of 1.36 m at 8.5 km depth. The cumulative afterslip moment within 9 months after the mainshock is 1 × 1018 N m, about 15.6% of that released by the mainshock coseismic slip. The afterslips contributed largely to the release of additional strain energy. In addition, shallow creep on the northern part of the blind fault, and deep uplift on the east normal fault system are promoted by stress perturbations. A significant proportion of down-dip coseismic slip spreading to more than 20 km beneath the surface, and deep up-slip afterslip have implications for the rheology of down-dip extension of the dipping faults in northern Tibet. Two obvious stress loading zones of more than 1 bar highlight seismic hazards in the region, especially in the junction between normal faults and ends of the large-size sinistral Riganpei-Co and Jiangai-Zangbu faults. It is necessary to forecast accurately by longer-term afterslip observation over timescales of years for the faults. Compared with previous studies, our results suggest a more complex subsurface fault geometry linking the normal and strike-slip faults and dynamic stress adjustment in this poorly-known region of Tibet.

Active Faults Revealed and New Constraints on Their Seismogenic Depth from a High-Resolution Regional Focal Mechanism Catalog in Myanmar (2016–2021)

ABSTRACTWe derive a new earthquake focal mechanism catalog for 86 Mw&amp;gt;4.0 earthquakes that occurred in the Myanmar region from 2016 to 2021. We apply the generalized Cut-and-Paste inversion method to a new set of regional broadband waveform data to obtain the earthquake focal mechanism and centroid depth with uncertainties estimated in a bootstrapping manner. Compared with global earthquake catalogs, our results are better aligned with mapped, active faults and reveal seismic activity along unmapped, blind faults. Our new catalog shows that the Sagaing Fault is more active in its northern segment with deeper seismogenic zone (∼27 km) compared to its southern segment that has a shallower seismogenic zone (∼10 km), sandwiching a seismic gap in its central segment. Earthquakes that occurred on the unmapped, blind faults beneath the Central Myanmar Basin at shallow depths (3–12 km) suggest a dominating northeast–southwest compressional stress field. Shallow earthquakes beneath the Indo-Myanmar Range (IMR) are rare, instead, north–south-oriented strike-slip faults are active within the deep accretionary wedge or lower crust of the Myanmar plate between depths of 20 and 40 km. At the eastern edge of the IMR, earthquakes with high-angle thrust mechanisms occurred between depths of 30 and 48 km, likely along steep faults separating the accretionary wedge from the Myanmar forearc crust. High-resolution intraslab focal mechanisms show that to the north of 22° N, slab deformation is dominated by strike-slip earthquakes with subvertical fault planes down to a depth of ∼25 km beneath the slab, suggesting lateral shear within the slab due to the northward motion of the Indian plate. To the south, more normal-faulting earthquakes suggest a stronger role of plate-bending processes in the slab deformation.

Modeling of Continental Normal Fault Earthquakes

AbstractThe magnitude of earthquakes on continental normal faults rarely exceeds 7.0 Mw. However, because of their vicinity to large population centers they can be highly destructive. Long recurrence time, relatively small deformations, and limited observations hinder our understanding of the deformation patterns and mechanisms controlling the magnitude of events. Here, this problem is addressed with 2D thermomechanical modeling of normal fault seismic cycles. The 2020 Samos, Greece Mw7.0 earthquake is used as an example as it is one of the largest and most studied continental normal fault earthquakes. The modeling approach employs visco‐elasto‐plastic rheology, compressibility, free surface, and a rate‐and‐state friction law for the fault. Modeling of the Samos earthquake suggests the pore fluid pressure ratio on the fault ranges from 0 to 0.7. The model demonstrates that most of the deformation during interseismic and coseismic periods, besides on the fault, occurs in the hanging wall and footwall below the seismogenic part of the fault. The largest vertical surface displacement during the earthquake is the subsidence of the hanging wall in the vicinity of the fault, while the uplift of the footwall and remote part of the hanging wall is significantly smaller. Modeling of the seismic cycles on normal faults with different setups shows the dependency of the magnitude on the thermal profile and dipping angle of the fault; low heat flow and low dipping angle are favorable conditions for the largest events, while steep normal faults in the areas of high heat flow tend to have the smallest magnitudes.

A Deeper Look Into the 2021 Tyrnavos Earthquake Sequence (TES) Reveals Coseismic Breaching of an Unrecognized Large‐Scale Fault Relay Zone in Continental Greece

AbstractLarge magnitude (Mw ∼ ≥6) earthquakes in extensional settings are often associated with simultaneous rupture of multiple normal faults as a result of static and/or dynamic stress transfer. Here, we report details of the coseismic breaching of a previously unrecognized large‐scale fault relay zone in central Greece, through three successive normal fault earthquakes of moderate magnitude (Mw 5.7–6.3) that occurred over a period of ∼10 days in March 2021. Specifically, joint analysis of InSAR, GNSS and seismological data, coupled with detailed field and digital fault mapping, reveals that the Tyrnavos Earthquake Sequence (TES) was accommodated at the northern end of a ∼100 km wide transfer structure, by faults largely unbroken during the Holocene. By contrast, the southern section of this relay zone appears to have accrued significant slip during Holocene. InSAR‐derived displacements agree with the loci of eight subtle, previously undetected, faults that accommodated coseismic and/or syn‐seismic normal fault slip during the TES. Kinematic modeling coupled with fault mapping suggests that all involved faults are interconnected at depth, with their conjugate fault‐intersections acting largely as barriers to coseismic rupture propagation. We also find that the TES mainshocks were characterized by unusually high (&gt;6 MPa) stress‐drop values that scale inversely with rupture length and earthquake magnitude. These findings, collectively suggest that the TES propagated north‐westward to rupture increasingly stronger asperities at fault intersections, transferring slip between the tips of a well‐established, but previously unrecognized, relay structure. Fault relay zones may be prone to high stress‐drop earthquakes and associated elevated seismic hazard.

The impact of faulting complexity and type on earthquake rupture dynamics

The statistical properties of seismicity are known to be affected by several factors such as the rheological parameters of rocks. We analysed the earthquake double-couple as a function of the faulting type. Here we show that it impacts the moment tensors of earthquakes: thrust-faulting events are characterized by higher double-couple components with respect to strike-slip- and normal-faulting earthquakes. Our results are coherent with the stress dependence of the scaling exponent of the Gutenberg-Richter law, which is anticorrelated to the double-couple. We suggest that the structural and tectonic control of seismicity may have its origin in the complexity of the seismogenic source marked by the width of the cataclastic damage zone and by the slip of different fault planes during the same seismic event; the sharper and concentrated the slip as along faults, the higher the double-couple. This phenomenon may introduce bias in magnitude estimation, with possible impact on seismic forecasting.

FOCAL MECHANISMS OF EARTHQUAKES IN THE SUBDUCTION ZONE OF THE WESTERN PACIFIC PLATE

Deformation features of the subducting Pacific lithospheric plate are considered according to the data on earthquake focal mechanisms. The territory includes the convergent boundaries between the Pacific Plate and the North American (in the Aleutian arc region), the Okhotsk, the Eurasian and the Philippine plates.It has been shown that the angle of subducting Pacific Plate in the Aleutian subduction zone affects the focal mechanisms of earthquakes that occurred in the upper, 35 km part of the oceanic plate in the zone of its bending. There occur normal-fault earthquakes at a steep-angle subduction and rare thrust earthquakes at a shallow-angle subduction. The azimuthal orientation of P-axes of the focal mechanism solutions in the upper (1–70 km) contact zone corresponds to the Pacific Plate displacement vector when the plate fragments are subducting west-northwestwards. There occurs a change in azimuthal orientation of the compression axes in the subducting plate at a depth of more than 70 km: the axes occupy different azimuthal sectors showing difference in the orientation of their slope, with the orientations of the T-axes become multidirectional.The calculation of seismotectonic deformations was carried out based on the data on focal mechanisms of 7768 earthquakes. It was revealed that the Exx and Ezz deformation fields are the most homogeneous at depths of 1–70 km. The pattern of seismotectonic deformations changes abruptly for deep parts of the subducting plate (105–200, 200–400, and 400–700 km), there are observed heterogeneous deformation fields Exx, Eyy and Еzz with alternating episodes of extension and shortening.There has been proposed the author’s scheme of the influence of the upper mantle convection structure on the geometry of the subducting plate (slab) as a potential catalyst for the processes responsible for the separation of seismic activity zones and the change of earthquake types with depth and in different parts of the extended subduction zone.

The cryptic seismic potential of the Pichilemu blind fault in Chile revealed by off-fault geomorphology

The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.

Empirical Evidence of Frequency‐Dependent Directivity Effects From Small‐To‐Moderate Normal Fault Earthquakes in Central Italy

AbstractRupture directivity and its potential frequency dependence is an open issue within the seismological community, especially for small‐to‐moderate events. Here, we provide a statistical overview based on empirical evidence of seismological observations, thanks to the large amount of high‐quality seismic recordings (more than 30,000 waveforms) from Central Italy, which represents an excellent and almost unique natural laboratory of normal faulting earthquakes in the magnitude range between 3.4 and 6.5 within the time frame 2008–2018. In order to detect an anisotropic distribution of ground motion amplitudes due to the rupture directivity, we fit the smoothed Fourier Amplitude Spectra (FAS) cleared of source‐, site‐ and path‐effects. According to our criteria, about 36% of the analyzed events (162 out of 456) are directive and the distribution of rupture direction is aligned with the strikes of the major faults of the Central Apennines. We find that the directivity is a band‐limited phenomenon whose width may extend up to five times the corner frequency. The results of this research provide useful insights to parameterize directivity, to be explicitly implemented in future ground motion modeling and scenario predictions.

A geomorphic-process-based cellular automata model of colluvial wedge morphology and stratigraphy

Abstract. The development of colluvial wedges at the base of fault scarps following normal-faulting earthquakes serves as a sedimentary record of paleoearthquakes and is thus crucial in assessing seismic hazard. Although there is a large body of observations of colluvial wedge development, connecting this knowledge to the physics of sediment transport can open new frontiers in our understanding. To explore theoretical colluvial wedge evolution, we develop a cellular automata model driven by the production and disturbance (e.g., bioturbative reworking) of mobile regolith and fault-scarp collapse. We consider both 90 and 60∘ dipping faults and allow the colluvial wedges to develop over 2000 model years. By tracking sediment transport time, velocity, and provenance, we classify cells into analogs for the debris and wash sedimentary facies commonly described in paleoseismic studies. High values of mobile regolith production and disturbance rates produce relatively larger and more wash-facies-dominated wedges, whereas lower values produced relatively smaller, debris-facies-dominated wedges. Higher lateral collapse rates lead to more debris facies relative to wash facies. Many of the modeled colluvial wedges fully developed within 2000 model years after the earthquake, with many being much faster when process rates are high. Finally, for scenarios with the same amount of vertical displacement, differently sized colluvial wedges developed depending on the rates of geomorphic processes and fault dip. A change in these variables, say by environmental change such as precipitation rates, could theoretically result in different colluvial wedge facies assemblages for the same characteristic earthquake rupture scenario. Finally, the stochastic nature of collapse events, when coupled with high disturbance, illustrates that multiple phases of colluvial deposition are theoretically possible for a single earthquake event.

Preliminary Analysis of Ionospheric Anomalies before Strong Earthquakes in and around Mainland China

The aims of the present study were to use Langmuir Probe payload electron density data and Plasma Analyzer Package payload O+ density data from the Zhangheng-1 electromagnetic satellite to statistically analyze anomalies in electron and oxygen ion densities before strong earthquakes (Ms ≥ 6.0) in western China and its neighboring areas. The goal was to investigate the physical mechanisms underlying electron and oxygen ion generation by evaluating the correlations between such anomalies and the seismic activity before the 6.4-magnitude earthquake in Yangbi, Yunnan, China, on 21 May 2021. Nine (75%) of the twelve earthquakes that occurred during the study period and were not affected by magnetic storms were preceded by anomalous electron or oxygen ion densities of 1.1–4.5 × 1010/m3 and 2.8–6.0 × 1010/m3, respectively. The anomalies were generally observed within the two weeks preceding the earthquakes and were associated with most strike-slip and thrust earthquakes, which were mainly located on the southeastern and northwestern edges of the Tibetan Plateau—but not normal fault earthquakes. The anomalies were likely the result of acoustic-gravity waves generated by slow vibrations of the Earth’s surface reaching the ionosphere, where they cause oscillations in ionospheric electron and ion densities. In addition, the association between ionospheric anomalies and strong earthquakes was confirmed by the observation of other atmospheric anomalies before the Yangbi earthquake.

The July 2020 Mw 6.3 Nima Earthquake, Central Tibet: A Shallow Normal-Faulting Event Rupturing in a Stepover Zone

Abstract On 22 July 2020, an Mw 6.3 earthquake with a predominantly normal-faulting mechanism struck the Yibug Caka fault zone, central Tibet, where the overall tectonic environment is characterized by left-lateral strike-slip motion. This event offers a chance to gain insight into the tectonic deformation and the cause of shallow normal-faulting earthquakes in this little studied region. Here, we use Sentinel-1A/B Interferometric Synthetic Aperture Radar data to investigate the coseismic and postseismic deformation related to this earthquake. The earthquake ruptured a previously mapped West Yibug Caka fault and is dominated by normal slip with a peak value of 1.9 m at depth of 6.9 km. Postseismic deformation analysis indicates that the observed subsidence signals of up to ∼4.7 cm are a consequence of afterslip. Most of the afterslip is confined at depths between 0.8 and 8.4 km, peaking at 0.27 m at depth of 6.1 km. The significant coseismic slip and afterslip involved in the earthquake highlights a complex interaction between the major normal fault and the secondary synthetic fault. By an integrated analysis of the inversions, regional geology geomorphology, fault kinematics, and seismicity background, we propose a tectonic model that attributes the occurrence of this normal-faulting event to the release of extensional stress in a stepover zone controlled by the northeast-striking sinistral strike-slip Riganpei Co fault and Bu Zang Ai fault. Compared with that the structural stepover often acts as a barrier to affect the propagation of earthquake rupture, our study demonstrates that the failure of a stepover may potentially induce the occurrence of earthquake along the bounding strike-slip faults.