High Stress Drop Research Articles

Dual-initiation ruptures in the 2024 Noto earthquake encircling a fault asperity at a swarm edge.

To reveal the connections between the 2024 moment magnitude (Mw) 7.5 Noto earthquake in Japan and the seismicity swarms that preceded it, we investigated its rupture process through near-source waveform analysis and source imaging techniques, combining seismic and geodetic datasets. We found notable complexity in the initial rupture stages. A strong fault asperity, which remained unbroken in preceding seismic swarms, slowed down the rupture. Then, a second rupture initiated at the opposite edge of the asperity, and the asperity succumbed to double-pincer rupture fronts. The failure of this high-stress drop asperity drove the earthquake into a large-scale event. Our observations help unravel the crucial role of fault asperities in controlling swarm migration and rupture propagation and underscore the need for detailed seismological and interdisciplinary studies to assess seismic risk in swarm-prone regions.

Insights into the Dead Sea Transform Activity through the study of fracture-induced electromagnetic radiation (FEMR) signals before the Syrian-Turkey earthquake (Mw-6.3) on 20.2.2023

Observations of fracture-induced electromagnetic radiation (FEMR) were conducted along the Dead Sea Transform (DST) from Sodom to Jericho, coinciding with a magnitude (Mw) 6.3 aftershock earthquake (EQ) in the Turkey-Syrian region on February 20, 2023. The FEMR parameters (“hits,” Benioff strain release, frequency, rise-time, energy) and associated crack dimensions were analyzed, focusing on trends leading up to the EQ. This study investigated the Benioff Strain plot and other parameters in three consecutive earthquake nucleation stages leading to the catastrophe. The first stage showed increased FEMR hits and frequency, decreased rise time (T′), and crack dimensions. In the second stage, FEMR hits and crack width decreased while other parameters continued to rise, accumulating the second-highest energy, likely due to high-stress drop. The third stage exhibited steadily increasing FEMR hits and energy and a notable increase in crack dimensions, suggesting an imminent macro failure event. The cyclic trend in FEMR hits indicates alternating periods of high activity and silence, potentially linked to stress changes during crack propagation. Taken shortly before the earthquake, these measurements offer valuable insights into how FEMR parameters vary before seismic events, bridging the gap between lab-scale studies of rock collapses under stress and large-scale failure phenomena.

Portevin-Le Chatelier characterization of annealed Al-Mg alloy sheets with different Mg contents

The stress-strain flow curves of solid solution treated and annealed Al-(2.86-9.41)Mg aluminum alloy sheets were investigated during tensile process at room temperature at a strain rate of 1 × 10−3 s−1, in order to reveal the effect of solute Mg atoms and β-phase nanoparticles in the α-Al matrix on their Portevin-Le Chatelier (PLC) characteristics. The results show that with the increase of Mg contents of alloys, the concentration of solute Mg atoms and the volume fraction of β-phase nanoparticles in the α-Al matrix of the annealed alloy sheets increase. The β-phase nanoparticles promote the PLC effect during the tensile process, causing the increase of stress drop amplitude Δσ of the tensile curves of the alloy sheets. However, the β-phase nanoparticles have little influence on the strain period ΔεTmax of rectangular waves. The flow strain difference Δεc−σ between the yield point and the beginning of PLC effect decreases monotonically from 5.4 × 10−4 to 5 × 10−5 as the Mg contents of alloy sheets increase from 2.86 wt.% to 9.41 wt.%. The strain period ΔεTmax of the rectangular waves corresponding to the high flow stress and stress drop amplitude Δσ obtained for the serrated stress-strain curves of the alloy sheets gradually increase with the increasing Mg contents of alloy sheets.

Broadband source model of the 2021 MW 7.1 Fukushima‐ken Oki earthquake in Japan based on seafloor and onshore strong‐motion records

AbstractThe 2021 MW 7.1 Fukushima‐ken Oki earthquake was the largest down‐dip compressional intraslab earthquake since August 2016, which is when S‐net records from the seafloor started to be provided. In this study, the empirical Green's function method was used with strong‐motion records from S‐net on the seafloor and KiK‐net onshore to estimate the broadband (0.2–10 Hz) source model of this earthquake. Three strong‐motion generation areas (SMGAs) with very high stress drops were estimated around the edge of the fault. One SMGA with a stress drop of 125 MPa was located in the oceanic crust in the shallow part of the Pacific Plate. The other two SMGAs with stress drops of 313 and 188 MPa were located in the deep part where previous source models without S‐net estimated a small slip. The two deep SMGAs increased the contribution of S‐net to more than that of KiK‐net. The short‐period spectral level was as high as that of the 2011 MW 7.1 Miyagi‐ken Oki earthquake, which was also a down‐dip compressional intraslab earthquake. The estimated broadband source model simulated the horizontal and vertical motions well at stations within 200 km of the hypocenter, except for surface waves at a few S‐net stations.

A Comprehensive Stress Drop Map From Trench to Depth in the Northern Chilean Subduction Zone

AbstractWe compute stress drops for earthquakes in Northern Chile recorded between 2007 and 2021. By applying two analysis techniques, (a) the spectral ratio (SR) method and (b) the spectral decomposition (SDC) method, a stress drop map for the subduction zone consisting of 51,510 stress drop values is produced. We build an extended set of empirical Green’s functions (EGF) for the SR method by systematic template matching. Outputs are used to compare with results from the SDC approach, where we apply cell‐wise obtained global EGF's to compensate for the structural heterogeneity of the subduction zone. We find a good consistency of results of the two methods. The increased spatial coverage and quantity of stress drop estimates from the SDC method facilitate a consistent stress drop mapping of the different seismotectonic domains. Albeit only small differences of median stress drop, strike‐perpendicular depth sections clearly reveal systematic variations, with earthquakes at different seismotectonic locations exhibiting distinct values. In particular, interface seismicity is characterized by the lowest observed median value, whereas upper plate earthquakes show noticeably higher stress drop values. Intermediate depth earthquakes show comparatively high average stress drop and a rather strong depth‐dependent increase of median stress drop. Additionally, we observe spatio‐temporal variability of stress drops related to the occurrence of the two megathrust earthquakes in the study region. The presented study is the first coherent large scale 3D stress drop mapping of the Northern Chilean subduction zone. It provides an important component for further detailed analysis of the physics of earthquake ruptures.

Depth-dependent seismic anomalies and potential asperity linked to fluid-driven crustal structure in Garhwal region, NW Himalaya

Characterizing the stress heterogeneity and crustal complexities are essential to understand the important aspects of rupture development and its propagation in a tectonically active region. We analyze a continuous seismic data (1999–2020) to estimate stress drops, b-value, fractal dimension (Dc) and focal mechanisms using a local network of 14 broadband stations. Seismic periodicity on a Himalayan major thrust above locked zone, involves long period of stress accumulation and interseismic strain. There exist a good spatio-temporal-depth correlation between the estimated results which indicates the positive anomaly at 12–14 km (stress transition) beneath the Chamoli region. Interestingly, this region shows maximum geodetic strain rate with comparatively high-stress drop, low coefficient friction value and very low b- (0.583 ± 0.02) and Dc- (1.20 ± 0.01) values. Majority of moderate to strong earthquakes and swarms activity have occurred in this region with maximum crustal accommodation. It is also observed that the significant temporal decrease in b-value always followed by moderate earthquakes. Here shallow crustal potential seismogenic asperity controls the earthquake's stress drop and rupture propagation. The observed correlations strongly support the evidence of entrapped fluid-faults interaction which alters pore pressure and generates numerous ruptures. Also, the role of pressurized entrapped fluids and roughness of fault-plane responsible for low friction zone and microseismicity. Continuous stress accumulation along the fault-plane within the asperity could release significant amount of energy in future.

Re-examination and comparison of activation volume, strain rate sensitivity, and dislocation density for coarse grain and CGP-processed UFG materials

The frequently used metallurgical parameters such as activation volume, internal stress, strain rate sensitivity, and dislocation density are identified in ultrafine-grained (UFG) materials through erudite analytical relations. It is difficult to analyze such parameters optically since the material consists of highly distorted and inhomogeneous microstructure due to the large strain imposed during the constrained groove pressing (CGP) process. Hence, transient tests like stress relaxation tests are advantageous over the optical to use. Apart from the identification of parameters in low-carbon steel materials using uniaxial loading and stress relaxation tests, the relationship between the parameters is nobly investigated in as-received and UFG structure samples produced through the CGP process. In the relationship, the exponential-power law constitutive to describe the accurate stress–strain curves, identification of mobile and forest dislocation density during strain hardening, and exponential variation between grain size and activation volume are proposed. The replenishment or exhaustion of mobile and forest dislocation density significantly affects the metallurgical parameters during uniaxial tensile loading followed by stress relaxation testing. In this experimental and analytical work, the CGP process produces a UFG material that has lower activation volume and higher stress drop, mobile dislocation density, forest dislocations density, strain rate sensitivity, and internal stress.

Spectral Characteristics of Hydraulic Fracturing-Induced Seismicity Can Distinguish between Activation of Faults and Fractures

Abstract Analysis of earthquake spectra can aid in understanding source characteristics like stress drop and rupture complexity. There is growing interest in probing the similarities and differences of fault rupture for natural and human-induced seismic events. Here, we analyze waveform data from a shallow, buried geophone array that recorded seismicity during a hydraulic fracturing operation near Fox Creek, Alberta. Starting from a quality-controlled catalog of 4000 events between magnitude 0 and 3.2, we estimate source-spectral corner frequencies using methods that account for the band-limited nature of the sensor response. The stress-drop values are found to be approximately self-similar, but with a slight magnitude dependence in which larger events have higher stress drop (∼10 MPa). Careful analysis of the relative corner frequencies shows that individual fault and fracture segments experienced systematic variations in relative corner frequency over time, indicating a possible change in the stress state. Clustering analysis of source spectra based on the relative proportion of high- and low-frequency content relative to the Brune model further shows that event complexity evolves over time. In addition, the faults produce earthquakes with systematically larger stress-drop values than the fractures. Combined, these results indicate that the features activated by hydraulic fracturing experience observable changes in source behavior over time and exhibit different properties depending on the orientation, scale, and fabric of the structural feature on which they occur.

The Rocks That Did Not Fall: A Multidisciplinary Analysis of Near‐Source Ground Motions From an Active Normal Fault

AbstractOn 8 July 2021 a M6.0 normal faulting earthquake rocked the community of Walker and the surrounding region near the California‐Nevada border. In the 1990s, field surveys of nearby Meadowcliff Canyon identified numerous precarious rocks deemed likely to topple in the event of strong shaking. Despite their proximity (∼6 km) to the 2021 earthquake, the precarious rocks still remain standing. In this work, we combine advanced source and ground motion characterization techniques to help unravel this mystery. High‐precision hypocentral locations reveal a clear north/south‐striking, east‐dipping rupture plane along the southern extension of the Slinkard Valley fault. The mainshock nucleated near the base of the fault, triggering thousands of aftershocks. Bayesian source spectral analyses indicate that the mainshock had a moderately‐high stress drop (∼17 MPa), and that aftershocks with deeper hypocenters have higher stress drops. Peak Ground Acceleration (PGA) recordings at regional stations agree well with existing ground motion models, predicting PGA of ∼0.3 g in Meadowcliff Canyon, a level sufficient to topple precarious rocks based on PGA‐derived stability criteria. We demonstrate that despite these large ground accelerations, the pulse duration in Meadowcliff Canyon is too short to supply the impulse necessary to damage these features, observations which support the application of dynamic toppling models that account for the joint effects of pulse amplitude and duration when assessing rock fragility. This study provides a unique vantage point from which to interpret rarely‐observed strong‐motion recordings from close to an active normal fault, one of many that dominate hazard along the eastern Sierra.

Eastward expansion of the Tibetan plateau: Insights from stress drops of the 2021 Ms 6.4 Yangbi, Yunnan and Ms 7.4 Maduo, Qinghai earthquake sequences in China

The eastward expansion of the Tibetan Plateau has resulted in extensive seismic activities in the eastern Tibetan Plateau, along with various types of faulting processes across this region. However, the crustal stress status related to these phenomena remains unclear. On 21 May 2021, the Ms 6.4 Yangbi and Ms 7.4 Maduo earthquakes with diverse focal mechanisms occurred consecutively in the Chuandian and Bayan Har blocks in the southern and northern parts of the eastern Tibetan Plateau, providing a distinct opportunity to directly compare the stresses in both flanks of the laterally expanding plateau. In this study, we estimate the stress drops of the Yangbi and Maduo earthquake sequences for all M≥3.0 events from the Lg-wave spectra. Both the path attenuation and site effect are properly removed based on an established broadband Lg-wave attenuation model. Then, the distribution of the stress drops is refined through high-precision earthquake relocation. Quick decays in stress drops were observed after both mainshocks. The Yangbi sequence has a relatively high median stress drop value, with significantly high stress drops for some foreshocks that possibly indicate a cascade triggering mechanism for the nucleation of the mainshock. In comparison, the Maduo sequence is characterized by a higher stress drop for the mainshock and sustained large fluctuation in stress drops for aftershocks. The Maduo aftershocks extend both the eastern and western directions from the mainshock, with relatively low stress drops in the east, suggesting a large amount of energy has been spent to drive the supershear rupture during the mainshock. We also find that the stress drops are associated with fault junctions for the Yangbi sequence, and are associated with variations in strength along the fault for the Maduo sequence, which are very likely linked with patterns of crustal motion and deformation in the eastern Tibetan Plateau.

Rupture Scenarios for the 3 June 1770 Haiti Earthquake

ABSTRACT The 2010 M 7.0 Haiti earthquake provided the impetus to reconsider historical earthquakes in Hispaniola (Bakun et al., 2012). That earthquake also shed new light on complex fault systems along Haiti’s southern peninsula (Douilly et al., 2013; Saint Fleur et al., 2015). Recently, the 2021 M 7.2 Nippes earthquake (Calais et al., 2022; Douilly et al., 2022), and a recent study reconsidering the 1860 sequence (Martin et al., 2022) further underscored the complexity of fault systems and large earthquake ruptures along the peninsula. Motivated by these studies and recent geological investigations (Prentice et al., 2010; Saint Fleur et al., 2020), we reconsider the 3 June 1770 Haiti earthquake to explore the conventional assumption that it was the last major (M ≥7.5) earthquake along the Enriquillo–Plantain Garden fault (EPGF). Accounts provide compelling evidence for substantial liquefaction in the Cul-de-Sac plain, one or more likely landslide-driven tsunami in Gonaïves Bay, and extensive landsliding that created at least three documented landslide dams. We consider three end-member rupture scenarios that are consistent with available constraints: two scenarios with M 7.7 and rupture lengths of 150–170 km, and one scenario with a ∼90 km rupture and M 7.5. Absent future work to identify and date paleoevents along the southern peninsula, none of these scenarios can be ruled out. Our preferred rupture model extends from the Miragoâne pull-apart to near la Selle mountain, with a rupture length of 127 km, M 7.6, and a high stress drop. Rupture could have been on the EPGF or on an oblique thrust fault associated with overthrusting of the Massif de la Selle. The results do support the conclusion that the 1770 earthquake was the last major earthquake in southern Haiti, with a magnitude upward of M 7.5 and significantly more severe shaking in southern Haiti than during the 2010 earthquake.

Characteristics of stress field and crustal kinematics of the southern region of Sichuan-Yunnan block, southern Tibet Plateau

The Red River fault is a massive tectonic shear zone of the Sichuan-Yunnan block, southern Tibet Plateau. The target area of this study is the Red River fault's middle-northern section and adjacent areas. Considering the spatial difference of tectonic structures and the inhomogeneity of focal mechanism solutions, we divided the target area into three structural units: the northern extension area, the middle strike-slip movement area, and the west complex rupture area. The study is carried out on the basis of zoning. Using the CAP method to invert the focal mechanism of moderate earthquakes and the principal compressive and tensile stress. Restoring the seismic source spectrum of small and medium earthquakes, then inversion the stress drops in three sub-units based on Brune model using digital waveform records. Our study suggests that: (1) Many high-stress drop earthquakes occurred in the transitional areas from the Weixi-Qiaohou fault to the Red River fault and from the Chuxiong fault to the Qujiang fault, also occurred in the Changning-Baoshan-Yongping area, revealing a large accumulation of local tectonic stress. (2) Some events with Δσ ≥ 20 MPa have occurred in the prominent low-speed anomaly area, indicating that a passive deformation crossed the Red River fault and extended to the west. Yangbi MS 6.4 earthquake (2021) occurred in a strongly locked area with high potential earthquake risk. (3) In the north of the study region, the movement of the Sichuan-Yunnan block to SE-SSE orientation was governed by the dextral strike-slip movement of the NW strike fault; in the middle region, the SE-orientated movement of the block relative to the South China block was governed by the sinistral slip movement of the SN-strike faults; and in the west region, the tectonic stress was governed by clockwise rotation of the block and the sinistral slip movement of the NE strike faults. We suggest that strike slip and rotation of the Red River fault's middle section may drive and aggravate passive movement in the north section, resulting in locally geophysical field responses. Therefore, the potential risk of a large earthquake in the middle-northern section of the Red River fault should be monitored.

Stress dissipation and seismic potential in the central seismic gap of the north-west Himalaya

• One-third of the earthquakes in CSG show stress drop <1 bar. • High stress drop events are absent in the upper 10 km of the crust. • Source parameter values are lower in Kumaon than CSG. • Kangra and Kumaon segments can be tipped to host a future great earthquake. The central seismic gap (CSG) of the Himalaya continues to pose the ever more intriguing question of whether it has the potential to host a devastating earthquake. To address the issue, we examine the distribution and dissipation of post-seismic stress in the northwest segment of the CSG. This study computes the source parameters of 47 local earthquake events (3.5 ≤ M L ≤ 5.5) in the Garhwal-Kumaon Himalaya during the period 2016–2018 using P-wave spectral analysis. In addition, we compile the stress drop results for 326 earthquake events from numerous studies carried out across the seismic gap, to estimate the source parameter variation. The estimated seismic moments (10 12 ≤ M o ≤ 10 16 N-m), source radii (0.3 ≤ r ≤ 1.3 km) and stress drops (0.1 ≤ σ ≤ 40.6 bar) for Kumaon are observed to be significantly lower than the overall values for the CSG (10 10 ≤ Mo ≤ 10 20 N-m; 0.1 ≤ r ≤ 13.2 km; 0.01 ≤ σ ≤ 77 bar), showing incomplete dissipation of the accumulated stress. About one-third of the earthquakes in the CSG show very low stress drop (σ < 1 bar). The high stress drop events are predominant along the Main Himalayan Thrust at mid-crustal depths, with the upper ∼ 10 km of the brittle crust rarely hosting any big event. Based on the analysis, the scaling relation M o f c 3.677 = 6E + 15 N-m/s 3 has been proposed for the Kumaon region. Other scaling relations have been developed between important source parameters (M L , M w , M o , and σ) that may serve as valuable inputs for assessing the earthquake hazard in the region. The estimated b-value for the Kumaon Himalaya estimated from the local earthquake catalog is 0.64 ± 0.08, which is low compared to the Garhwal and rest of the NW Himalaya. The overall b-value of the CSG region computed from the ISC catalog is 0.62 ± 0.06. From the stress drop and b-value variation in the CSG, the Kangra and the Kumaon segments are identified to be potentially hazardous zones for a future great earthquake.

Source Parameters of the Mw 5.7 Pica Crustal Earthquake in Northern Chile

Abstract On 10 September, 2008, an Mw 5.7 earthquake occurred under the Central Valley of northern Chile near the town of Pica at a depth of ∼33 km within the continental crust of the South America plate. We find this earthquake to be a high stress-drop, reverse-oblique event that generated unusually high ground accelerations of up to 0.67g. Overall, its observed ground motion intensities are considerably larger than those predicted by ground motion models, particularly at short periods. The source properties inferred through waveform modeling indicate reverse-oblique fault motion on a ∼75 km2 plane dipping to the northeast, which is corroborated by the located aftershock distribution. Stress-drop values of the mainshock and larger aftershocks were estimated through S-wave spectrum modeling, with values up to ∼250 MPa for the mainshock. The event occurred in a cold section of the continental crust under the Central Valley, and its fault kinematics and orientation are consistent with the dominant style of faulting and stress field under the neighboring Coastal Cordillera. Although our recurrence analysis shows that crustal events in the region occur at a lower rate than interplate and inslab events, crustal events of similar or higher magnitude than the Pica earthquake have occurred, on average, approximately once every three years in northern Chile, which could pose an important hazard to nearby populations or critical infrastructure.

Portevin-Le Chatelier Characterization of Quenched Al-Mg Alloy Sheet with Different Mg Concentrations.

In the present study, the PLC characteristic parameters and DSA mechanism of Al-(2.86~9.41) Mg alloy sheets were investigated during tensile testing at room temperature with a tensile rate of 1 × 10−3 s−1. On the basis of the solution Mg concentrations in the α-Al matrix, the initial vacancy concentration, the second-phase particle configuration and the recrystallized grain configuration are almost the same by quenching treatment. The results show that the type of room-temperature tensile stress–strain curves of quenched Al-(2.86~9.41) Mg alloy sheets varied according to the Mg content. The type of stress–strain curve of the Al-2.86 Mg alloy sheet was B + C, while the type of stress–strain curve of the Al-(4.23~9.41) Mg alloy sheets was C. When the quenched Al-(2.86~9.41) Mg alloy sheets were stretched at room temperature, the strain cycle of the rectangular waves corresponding to the high stress flow and stress drop amplitude on the zigzag stress–strain curve of alloy sheets increased with increasing the Mg content. Moreover, the strain cycle of and on the stress–strain curve of alloy sheets increased gradually with increasing tensile deformation. The yield stress of quenched Al-(2.86~9.41) Mg alloy sheets increased gradually with increasing the Mg content. Moreover, the critical strain corresponding to yield stress and the critical strain corresponding to the occurrence of the PLC shearing band of alloy sheets both increased with increasing the Mg content. However, the difference in flow strain value between and of alloy sheets decreased gradually with increasing the Mg content.

Radiated seismic energy from the 2021 ML 5.8 and ML 6.2 Shoufeng (Hualien), Taiwan, earthquakes and their aftershocks

We calculated radiated seismic energy (ES), seismic moment (M0), and moment magnitude (MW) and then determined the ES–ML, ES–M0, M0–ML, and MW–ML relations for the 2021 Shoufeng earthquake sequence (2.5 < ML < 6.3), where ML is the local magnitude. Notably, a crossover magnitude was detected as ML = 4.0 for the earthquake sequence. For ML < 4.0, we obtained logM0 ∝ ML, MW ∝ 0.67ML, and a low ES/M0, indicating a low average stress drop; for ML > 4.0, logM0 ∝ 0.67ML, MW ∝ ML, and a high ES/M0 were present, and then there was a high average stress drop. These derived relations implied that source duration (T) is independent of M0 for ML < 4.0. Moreover, the M0 ∝ T3 relation seemed able to interpret those relations for ML > 4.0. Nevertheless, the ES–ML relation remains logES ∝ 2.0ML for 2.5 < ML < 6.3. From this study, the derived relations could not predict the source parameters for ML > 6.3 events. This might indicate that ML saturates beyond ML 6.3. Through such analyses, we not only established the relations among source parameters but also elucidated the basic physics of the earthquake sequence.

Analysis of the April 10, 2021 (Mw 6.1) destructive intra-slab earthquake, East Java, Indonesia

A destructive earthquake (Mw 6.1) struck East Java, Indonesia, on 10 April 2021, resulting in at least 10 deaths and damage to over 3000 buildings. Preliminary hypocenter determinations indicate that it was located ~20 km south of Java, at a depth of ~60 km, suggesting that it likely occurred within the subducting oceanic slab beneath the Sunda Arc. In order to better understand the origins of this important event, we have undertaken non-linear hypocenter determination of the mainshock and its aftershocks through picking of P- and S-wave arrival times recorded by nearby BMKG seismic stations, followed by precise relative relocation using the double-difference method. We have also conducted source mechanism analysis to determine the type of fault slip, and obtained a rupture model of the mainshock by performing a finite fault inversion. Our results show that the mainshock was intra-slab, with an estimated depth of 58 km and was followed by, on average, deeper aftershocks that ranged in depth from 57 to 69 km. The source mechanism solution of selected aftershocks reveal that they rupture along a thrust fault oriented almost parallel to the trench, which is consistent with the mainshock source solution. The rupture modelling results indicate that the earthquake is characterised by a high stress-drop and compact rupture area, which together with its depth of nearly 60 km, may explain the relatively small number of aftershocks that were recorded.

Integrated Analysis of the 2020 Mw 7.4 La Crucecita, Oaxaca, Mexico, Earthquake from Joint Inversion of Geodetic and Seismic Observations

ABSTRACT On 23 June 2020, a major thrust faulting earthquake ruptured the megathrust near a sharp transition of the plate boundary beneath the Mexico subduction zone, overlapping the rupture area of the 1965 Mw 7.5 (Ms 7.6) event. A joint inversion of geodetic and seismic observations and prediction of coastal uplifts indicate a compact rupture with a peak slip of ∼3.6 m. The cumulative seismic moment is 1.41×1020 N·m within 20 s of coseismic rupture, giving a moment magnitude of Mw 7.4. The rupture propagated with a relatively high rupture speed (∼3.5 km/s) downdip at depths from 20 to 30 km, terminating at the edge of a slow slip region. The energy-related static stress drop is ∼8.1 MPa. The aftershock distribution and the region where slow slip events occurred show a clear complementary relationship with the coseismic rupture zone of the 2020 Oaxaca event. The high-stress drop and associated deep compact rupture pattern suggest a homogeneous asperity with high strength beneath the Oaxaca region. Comparisons of rupture areas, mechanisms, and seismic waveforms show that the 2020 Mw 7.4 Oaxaca earthquake is a quasi-repeat of the 1965 earthquake. The inferred recurrence interval is ∼55 yr. Static stress transfer from mainshock may trigger near trench aftershocks but has small loading on the megathrust interface in the Tehuantepec seismic gap.

Source properties of the 2019ML6.3 Hualien, Taiwan, earthquake, determined by the local strong motion networks

SUMMARYThe 2019 ML 6.3 Hualien earthquake struck the northern Longitudinal Valley (LV) and generated not only large strong motions (intensity of 7, as defined by the Central Weather Bureau, Taiwan) locally but also widespread strong shaking in metropolises in northern Taiwan. In this study, we analyse strong motion records from local seismic networks to understand the source properties of the 2019 event. We first obtain the centroid location of the 2019 event using the source-scanning algorithm (SSA) technique by applying the unfiltered records. The determined centroid location is 121.55°E, 24.10°N, with a depth of 22.5 km. This location is 5.5 km north–northwest of and 3.8 km deeper than the Central Weather Bureau hypocentre, suggesting that the 2019 event occurred on the high-angle west-dipping plane of the focal solution. The centroid time delay is 3.35 s. Then, we obtain strong motion generation areas (SMGAs) of the 2019 event using the empirical Green's function method by considering the broad-band waveforms (0.4–10 Hz). Unlike other moderate-sized earthquakes in Taiwan, which have one SMGA, we determine that there were two SMGAs in the 2019 event. SMGA1 initiated at the CWB hypocentre with a size of 4.00 km2, and SMGA2 initiated at the centroid location determined by the SSA approach with a size of 3.63 km2. Such small areas cause high stress drops of 13.7 and 27.4 MPa for SMGA1 and SMGA2, respectively. We infer that the localized high stress drop of SMGAs is one of the important factors responsible for high peak-ground accelerations (PGAs) in Taiwan in addition to a strong directivity effect coupled with the radiation pattern reported by the previous study. Furthermore, previous moderate-sized earthquakes on an active structure called the Xiulin segment revealed similar source properties with a high stress drop and generated large PGA locally as well as in the metropolises of northern Taiwan. Considering the stored moment deficit, the probability of a future large earthquake in the northern LV region remains high. It is essential to consider seismic hazard assessment and mitigation for this not-well-known but high-seismic-potential region.

The 2020Mw 6.4 Petrinja earthquake: a dextral event with large coseismic slip highlights a complex fault system in northwestern Croatia

SUMMARYOn 2020 December 29, the Mw 6.4 Petrinja earthquake hit the Kupa Valley region and set a record for the largest earthquake in northwestern (NW) Croatia. The coseismic surface displacements are well obtained on three pairs of interferometric synthetic aperture radar images from Sentinel-1 satellites. The interferograms exhibit coseismic ground deformation with a maximum line-of-sight displacement of 0.4 m. Based on the coseismic deformation field, we investigate both the fault geometry and the coseismic slip distribution. The results show a dextral event with a peak slip of 3.50 m at a depth of 3.47 km. The shallow depth and unusually large coseismic slip correspond to obvious ground deformation and serious damage in the epicentral zone. The 2020 earthquake highlights an unmapped, steeply dipping strike-slip fault, which possibly enabled a potential ‘curve cut-off’ process on the bending segment of the Pokupsko fault in the context of ∼N–S compression in NW Croatia. The large coseismic slip and high stress drop associated with the Mw 6.4 Petrinja earthquake are likely products of the geometrically complex fault zones and immature seismotectonic environment in NW Croatia.