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.

Present-day stress field pattern in the Vrancea seismic zone (Romania) deduced from earthquake focal mechanism inversion

Vrancea seismogenic zone in the South-Eastern Carpathians is characterized by localized intermediate-depth seismicity. Due to its complex geodynamics and large strain release, Vrancea represents a key element in the Carpatho-Pannonian system. Data from a recently compiled catalogue of fault plane solutions (REFMC) are inverted to evaluate stress regime in Vrancea on depth. A single predominant downdip extensive regime is obtained in all considered clusters, including the crustal layers located above the Vrancea slab. The prevalent stress regime confirms previous investigations and requires some mantle-crust coupling. The S3 principal stress is close to vertical, while S1 and S2 are horizontal, oriented perpendicularly and respectively tangentially to the Carpathians Arc bend. This configuration is present at any depth level. According to seismicity patterns, there are two main active segments in the Vrancea intermediate-depth domain, at 55 – 105 km and 105 – 180 km, both able to generate major events. The configuration of the tectonic stresses as resulted from inversion is similar in both segments. Also, high fault instability (I > 0.95) is characterizing the segments. The only notable difference is given by the friction and stress ratio parameters which drop down in the bottom segment from μ = 0.95 to μ = 0.55 and from R = 0.51 to R = 0.29. This variation is attributed to possible weakening processes activated below 100 km depth and can explain the intensification of seismicity production as earthquake rate and average energy release in the lower segment versus the upper segment.

CHINGIZ EARTHQUAKE on January 20, 2015, KR=12.2, mb=5.6, I_0^P=5–6 NEAR the SEMIPALATINSK POLYGON (Eastern Kazakhstan)

The work presents the results of data analysis of the earthquake on January 20, 2015, mb=5.6 occurred at the low-active region of Eastern Kazakhstan near the border of the former Semipalatinsk Test Site. The work briefly describes the seismicity of the Semipalatinsk Test Site related to aseismic regions according to the acting map of general seismic zoning. The earthquake was recorded by permanent seismic stations of the Institute of Geophysical Research of the Ministry of Energy of the Republic of Kazakhstan and by two strong-motion accelerometers. The earthquake was felt at distances up to 300 km from the epicentre, its maximum intensity was noted at Medeu settlement where the shaking intensity was 56 by the MSK-64 scale. The earthquake focal mechanism is shown, the parameters of seismic effects and response spectra are calculated.

Weak, Seismogenic Faults Inherited From Mesozoic Rifts Control Mountain Building in the Andean Foreland

New earthquake focal mechanism and centroid depth estimates show that the deformation style in the forelands of the Andes is spatially correlated with rift systems that stretched the South American lithosphere in the Mesozoic. Where the rifts trend sub-parallel to the Andean range front, normal faults inherited from the rifts are being reactivated as reverse faults, causing the 30--45 km thick seismogenic layer to break up. Where the rift systems are absent from beneath the range front, the seismogenic layer is bending and being thrust beneath the Andes like a rigid plate. Force-balance calculations show that the faults inerhited from former rift zones have an effective coefficient of static friction < 0.2. In order for these frictionally-weak faults to remain seismogenic in the lower crust, their wall rocks are likely to be formed of dry granulite. Xenolith data support this view, and suggest that parts of the lower crust are now mostly metastable, having experienced temperatures at least 75--250 degrees hotter than present. The conditions in the lower crust make it unlikely that highly-pressurised free water, or networks of intrinsically-weak phyllosilicate minerals, are the cause of their low effective friction, as, at such high temperatures, both mechanisms would cause the faults to deform through viscous creep and not frictional slip. Therefore pre-existing faults in the Andean forelands have remained weak and seismogenic after reactivation, and have influenced the style of mountain building in South America. However, the controls on their mechanical properties in the lower crust remain unclear.

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>

A functional tool to explore the reliability of micro-earthquake focal mechanism solution 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 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 operated in the Campania-Lucania Apennines (Southern Italy) to monitor the complex normal fault system activated during the Ms 6.9, 1980 earthquake. We demonstrate that the method we propose can have 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.