A functional tool to explore the reliability of micro-earthquake focal mechanism solutions for seismotectonic purposes

Improving the knowledge of seismogenic faults requires the integration of geological, seismological, and geophysical information. Among several analyses, the definition of earthquake focal mechanisms plays an essential role in providing information about the geometry of individual faults and the stress regime acting in a region. Fault plane solutions can be retrieved by several techniques operating in specific magnitude ranges, both in the time and frequency domain and using different data. For earthquakes of low magnitude, the limited number of available data and their uncertainties can compromise the stability of fault plane solutions. In this work, we propose a useful methodology to evaluate how well a seismic network, used to monitor natural and/or induced micro-seismicity, estimates focal mechanisms as a function of magnitude, location, and kinematics of seismic source and consequently their reliability in defining seismotectonic models. To study the consistency of focal mechanism solutions, we use a Bayesian approach that jointly inverts the P/S long-period spectral-level ratios and the P polarities to infer the fault plane solutions. We applied this methodology, by computing synthetic data, to the local seismic network operating in the Campania–Lucania Apennines (southern Italy) aimed to monitor the complex normal fault system activated during the Ms 6.9, 1980 earthquake. We demonstrate that the method we propose is effective and can be adapted for other case studies with a double purpose. It can be a valid tool to design or to test the performance of local seismic networks, and more generally it can be used to assign an absolute uncertainty to focal mechanism solutions fundamental for seismotectonic studies.

Spatiotemporal Variations of Focal Mechanism Solutions and Stress Field of the 2019 Changning Ms 6.0 Earthquake Sequence

The June 17, 2019, MS 6.0 Changning earthquake is the largest recorded event in the Sichuan basin, spatiotemporal variations of stress field may shed light on the seismogenic mechanism of the earthquake. We determined the focal mechanism solutions (FMSs) of 124 earthquakes with MS ≥ 3.0 occurring in the Changning area from April 1, 2007, to February 29, 2020, and analyzed changes of FMSs and stress field before and after Changning earthquake. The Changning aftershocks were predominantly thrust fault earthquakes, followed by strike slip. The P-axis azimuths of the aftershock FMSs were oriented predominantly in the NEE direction, notably differing from the NWW-oriented P-axis azimuths of pre-earthquake FMSs; it shows the rotation of local stress field before and after the Changning earthquake, it is speculated that the change of stress field in Changning area may be caused by long-term water injection and salt mining activities. From the southeast to the northwest of the aftershock zone, the azimuths of principal compressive stress (S1) change from NEE to near-EW in both horizontal and vertical planes. Significant changes occurred in the FMS types and stress field of the aftershock zone following the Changning earthquake, the FMSs became diverse, the S1 azimuth of the Changning area changed from NWW to NEE, and then EW, the plunge and stress tensor variances increased, it reflects that the stress field of the Changning area adjusts continually with time.

Fault Geometry and Mechanism of the Mw 5.7 Nakchu Earthquake in Tibet Inferred from InSAR Observations and Stress Measurements

Different types of focal mechanism solutions for the 19 March 2021 Mw 5.7 Nakchu earthquake, Tibet, limit our understanding of this earthquake’s seismogenic mechanism and geodynamic process. In this study, the coseismic deformation field was determined and the geometric parameters of the seismogenic fault were inverted via Interferometric Synthetic Aperture Radar (InSAR) processing of Sentinel-1 data. The inversion results show that the focal mechanism solutions of the Nakchu earthquake are 237°/69°/−70° (strike/dip/rake), indicating that the seismogenic fault is a NEE-trending, NW-dipping fault dominated by the normal faulting with minor sinistral strike-slip components. The regional tectonic stress field derived from the in-situ stress measurements shows that the orientation of maximum principal compressive stress around the epicenter of the Nakchu earthquake is NNE, subparallel to the fault strike, which controlled the dominant normal faulting. The occurrence of seven M ≥ 7.0 historical earthquakes since the M 7.0 Shenza earthquake in 1934 caused a stress increase of 1.16 × 105 Pa at the hypocenter, which significantly advanced the occurrence of the Nakchu earthquake. Based on a comprehensive analysis of stress fields and focal mechanisms of the Nakchu earthquake, we propose that the dominated normal faulting occurs to accommodate the NE-trending compression of the Indian Plate to the Eurasian Plate and the strong historical earthquakes hastened the process. These results provide a theoretical basis for understanding the geometry and mechanics of the seismogenic fault that produced the Nakchu earthquake.

Earthquake monitoring of the Baribis Fault near Jakarta, Indonesia, using borehole seismometers

The geological setting of Jakarta and its immediate surroundings are poorly understood, yet it is one of the few places in Indonesia that is impacted by earthquakes from both the Java subduction zone and active faults on land. In this study, a borehole seismic experiment with low noise characteristics was deployed to record seismic activity on the ~ E-W oriented Baribis Fault, which is ~ 130 km long, passes to the south of Jakarta, and is only ~ 20 km away at its nearest point. A primary objective of this study is to determine whether this fault is seismically active, and therefore, whether it might pose a threat to nearby population centers, including Jakarta in particular. A total of seven broadband instruments that spanned Jakarta and the surrounding region were installed between the end of July 2019 and August 2020, during which time we were able to detect and locate 91 earthquakes. Two earthquakes were located close to the Baribis Fault line, one of which was felt in Bekasi (southeast of Jakarta) where it registered II-III on the Modified Mercalli Intensity (MMI) scale. The focal mechanism solutions of these events indicate the presence of a thrust fault, which is in good agreement with previous studies, and suggest that the Baribis Fault is active.

Multisource stress data constraints on Cretaceous—present regional tectonic stress field evolution in the southern Jinzhou area, North China Craton

Regional tectonic stress fields are key crustal stress elements that drive tectonic movements and are associated with regional tectonics and geological resources. Regional tectonic stress field evolution of the Jinzhou area, located in the eastern block of the North China Craton (NCC), may provide a deeper understanding of tectonics of western Liaoning and the NCC. This work conducted borehole television, hydraulic fracturing and focal mechanism solutions to invert the paleo and present regional tectonic stress fields. Four groups of tensile fracture in the southern Jinzhou area were identified via borehole television, and their azimuths were NNW–SSE, NWW–SEE, nearly W–E and NE–SW in temporal order representing four stages of extensional tectonic events. Hydraulic fracturing and focal mechanism solutions showed that the stress status was normal fault and strike-slip, revealing that the southern Jinzhou area is undergoing NEE–SWW-oriented compression and nearly N–S-oriented extension in accordance with the strike-slip mechanism. From the Early Cretaceous to the present, the direction of the regional extensional stress in the southern Jinzhou area has evolved counterclockwise and sequentially from NNW–SSE to NWW–SEE, W–E, NE–SW and nearly N–S, and the regional tectonic mechanism has transited from extension to extension-strike-slip to strike-slip, leading to the current tectonic framework.

Directional Statistics, Bayesian Methods of Earthquake Focal Mechanism Estimation, and Their Application to New Zealand Seismicity Data

<p>A focal mechanism is a geometrical representation of fault slip during an earthquake. Reliable earthquake focal mechanism solutions are used to assess the tectonic characteristics of a region, and are required as inputs to the problem of estimating tectonic stress. We develop a new probabilistic (Bayesian) method for estimating the distribution of focal mechanism parameters based on seismic wave polarity data. Our approach has the advantage of enabling us to incorporate observational errors, particularly those arising from imperfectly known earthquake locations, allowing exploration of the entire parameter space, and leads to natural point estimates of focal mechanism parameters. We investigate the use of generalised Matrix Fisher distributions for parameterising focal mechanism uncertainties by minimising the Kullback-Leibler divergence. We present here the results of our method in two situations. We first consider the case in which the seismic velocity of the region of interest (described by a velocity model) is presumed to be precisely known, with application to seismic data from the Raukumara Peninsula, New Zealand. We then consider the case in which the velocity model is imperfectly known, with application to data from the Kawerau region, New Zealand. We find that our estimated focal mechanism solutions for the most part are consistent with all available polarity data, and correspond closely to solutions obtained using established methods. Further, the generalised Matrix Fisher distributions we examine provide a good fit to our Bayesian posterior PDF of the focal mechanism parameters, enabling the posterior PDF to be succinctly summarised by reporting the estimated parameters of the fitted distribution.</p>

Directional Statistics, Bayesian Methods of Earthquake Focal Mechanism Estimation, and Their Application to New Zealand Seismicity Data

<p>A focal mechanism is a geometrical representation of fault slip during an earthquake. Reliable earthquake focal mechanism solutions are used to assess the tectonic characteristics of a region, and are required as inputs to the problem of estimating tectonic stress. We develop a new probabilistic (Bayesian) method for estimating the distribution of focal mechanism parameters based on seismic wave polarity data. Our approach has the advantage of enabling us to incorporate observational errors, particularly those arising from imperfectly known earthquake locations, allowing exploration of the entire parameter space, and leads to natural point estimates of focal mechanism parameters. We investigate the use of generalised Matrix Fisher distributions for parameterising focal mechanism uncertainties by minimising the Kullback-Leibler divergence. We present here the results of our method in two situations. We first consider the case in which the seismic velocity of the region of interest (described by a velocity model) is presumed to be precisely known, with application to seismic data from the Raukumara Peninsula, New Zealand. We then consider the case in which the velocity model is imperfectly known, with application to data from the Kawerau region, New Zealand. We find that our estimated focal mechanism solutions for the most part are consistent with all available polarity data, and correspond closely to solutions obtained using established methods. Further, the generalised Matrix Fisher distributions we examine provide a good fit to our Bayesian posterior PDF of the focal mechanism parameters, enabling the posterior PDF to be succinctly summarised by reporting the estimated parameters of the fitted distribution.</p>

The Seismogenic Structure of the 2017 Mw 6.9 Milin, Tibet, Earthquake: A Possible Newly Active Fault at the Eastern Himalayan Syntaxis

On 18 November 2017, an Mw 6.9 earthquake occurred in Milin, Tibet, with the epicenter at the top of the Namche Barwa syntaxis. This event did not produce surface ruptures, and its seismogenic structure remains unclear or controversial. Using the locations of the Milin mainshock and aftershocks, locations of regional small earthquakes and focal mechanism solutions from 2007 to 2009, this work analyzed the causative fault and tectonic setting of the Milin earthquake and assessed the regional seismic risk. The results suggest that the seismogenic structure of the Milin earthquake was a secondary fault, the southern branch of the XiXingla fault (XXLF). Within 28 hr after the mainshock, the aftershocks of the Milin event spread northeastward to the secondary north branch fault of the XXLF and the secondary south branch fault of the Palong–Pangxin fault. Across the top of the Namche Barwa syntaxis (Namche Barwa block) and the Chayu block in the southeast, an earthquake dense belt (EDB) has developed. This EDB has similar deep structures beneath the two blocks, in which several northeast-dipping structural planes exit, and different portions of the EDB imply a unified tectonic stress field. Combining these data with the foreshock–mainshock–aftershock data for the 1950 Mw 8.6 Chayu, Tibet, earthquake, we speculate that the structural planes produced by the EDB at depth in the two blocks have already been connected or tended to connect, resulting in a new fault system trending northwest and approximately 280 km long. The 2017 Mw 6.9 Milin earthquake occurred at the northwestern end of this fault system. At present, the development stage, maturity, and fine structure of this new fault system remain unclear but should receive additional attention. Based on its maximum rupture area, this new fault system is capable of generating an Mw 7.7 earthquake in the future.

Focal Mechanism Analysis of the September 25th, 2019 Mw 6.5 Ambon Earthquake and Its Implication for Seismotectonics

On September 25th, 2019, an Mw 6.5 earthquake occurred in Ambon, Maluku Province, Indonesia, and caused casualties and infrastructures damages. The epicenter located in a tectonically active region with the potential strike-slip and thrust faulting earthquake sources, yet the responsible fault is still not well understood. Based on focal mechanism solutions from available seismological agencies, i.e. USGS, GFZ, GCMT, and BMKG, the earthquake has a similar strike-slip focal mechanism, although there are discrepancies on detailed source parameters. To provide a better understanding of the earthquake mechanism and seismotectonic, we apply the Cut-and-Paste (CAP) focal mechanism inversion method to broadband seismic waveforms from regional and teleseismic distances. The CAP inversion results on the regional data grouped in different distance ranges show a robust strike-slip solution. We then refine the earthquake focal depth by performing the CAPtele inversion and resulted in a depth of 12 km with similar fault plane solution as the regionals. The ruptured fault plane is resolved by a directivity analysis using azimuthal pattern of the apparent source durations, which indicates an obvious unilateral rupture propagation toward SSE direction. Our result suggests the NNW-SSE orientated fault is the ruptured fault plane, which is also consistent with the near N-S distributed aftershocks. This fault is located in a narrow sea between Seram, Ambon and Haruku island and was not reported yet in previous studies. The Coulomb failure stress (CFS) changes analysis of the mainshock shows that the Ambon earthquake has promoted the off-fault aftershocks which occurred to the west of the ruptured fault.

The key role of conjugate fault system in importing earthquakes into the eastern flank of the Red Sea

This study aims to synthesize seismic observations with gravity and magnetic data and to suggest a new scenario on the development of the Harrat Lunayyir (HL) tectonic system on the eastern Red Sea coastline, Saudi Arabia. Gravity and aeromagnetic anomalies distinctly mapped the NE and NW trends, while the InSAR data depict a small NW–SE graben and an NW–SE dyke. High-resolution relocations, which are well-consistent with the focal mechanism solutions for events with magnitudes greater than 3.0, admit two distinctly fault styles of different orientations. Thus, leading to the NE and NW fault planes’ reactivation related to the Precambrian basement faults and the Red Sea rift system, respectively. The spatiotemporal distributions of epicenters and focal mechanism solutions suggest a new seismic deformation scenario of the 2009 earthquake seismic activity. The low static frictions of 0.2–0.35 obtained from the stress inversion indicates reactivation of preexisting faults in the respective seismogenic zones. The obtained results give rise to a swarm-like sequence of tectonic implications, two activated fault styles differently oriented, and an NE conjugate fault system inherited in the region, which plays a vital role in transferring the ambient stress regime into the Red Sea’s eastern flank.