A Retrospective Analysis of b-Value Changes Preceding Strong Earthquakes

Author(s):  
Nicolas D. DeSalvio ◽  
Maxwell L. Rudolph

Abstract Earthquake precursors have long been sought as a means to predict earthquakes with very limited success. Recently, it has been suggested that a decrease in the Gutenberg–Richter b-value after a magnitude 6 earthquake is predictive of an imminent mainshock of larger magnitude, and a three-level traffic-light system has been proposed. However, this method is dependent on parameters that must be chosen by an expert. We systematically explore the parameter space to find an optimal set of parameters based on the Matthews correlation coefficient. For each parameter combination, we analyze the temporal changes in the frequency–magnitude distribution for every M ≥ 6 earthquake sequence in the U.S. Geological Survey Comprehensive Earthquake Catalog for western North America. We then consider smaller events, those with a foreshock magnitude as small as 5, and repeat the analysis to assess its performance for events that modify stresses over smaller spatial regions. We analyze 25 M ≥ 6 events and 88 M 5–6 events. We find that no perfect parameter combination exists. Although the method generates correct retrodictions for some M 5 events, the predictions are dependent on the retrospectively selected parameters. About 80%–95% of magnitude 5–6 events have too little data to generate a result. Predictions are time dependent and have large uncertainties. Without a precise definition of precursory b-value changes, this and similar prediction schemes are incompatible with the IASPEI criteria for evaluating earthquake precursors. If limitations on measuring precursory changes in seismicity and relating them to the state of stress in the crust can be overcome, real-time forecasting of mainshocks could reduce the loss of lives.

2014 ◽  
Vol 08 (04) ◽  
pp. 1450010 ◽  
Author(s):  
Santi Pailoplee

In this study, the geospatial frequency–magnitude distribution (FMD) b-value images of the prospect sources of upcoming earthquakes were investigated along the Indonesian Sunda Margin (ISM) that strikes parallel to and near the Indonesian Island chain. After enhancing the completeness and stability of the earthquake catalogue, the seismicity data were separated according to their seismotectonic setting into shallow crustal and Intraslab earthquakes. In order to verify the spatial relationship between the b-values and the occurrence of subsequent major earthquakes, the complete shallow crustal seismicity dataset (1980–2005) was truncated into the 1980–2000 sub-dataset. Utilizing the suitable assumption of fixed-number of earthquakes, retrospective tests of both the complete and truncated datasets supported that areas of comparatively low b-values could reasonably be expected to predict likely hypocenters of future earthquakes. As a result, the present-day distributions of b-values derived from the complete (1980–2005) shallow crustal and Intraslab seismicity datasets revealed eight and six earthquake-prone areas, respectively, along the ISM. Since most of these high risk areas proposed here are quite close to the major cities of Indonesia, attention should be paid and mitigation plans should be developed for both seismic and tsunami hazards.


2022 ◽  
Author(s):  
Marcus Herrmann ◽  
Ester Piegari ◽  
Warner Marzocchi

Abstract The Magnitude–Frequency-Distribution (MFD) of earthquakes is typically modeled with the (tapered) Gutenberg–Richter relation. The main parameter of this relation, the b-value, controls the relative rate of small and large earthquakes. Resolving spatiotemporal variations of the b-value is critical to understanding the earthquake occurrence process and improving earthquake forecasting. However, this variation is not well understood. Here we present unexpected MFD variability using a high-resolution earthquake catalog of the 2016–2017 central Italy sequence. Isolation of seismicity clusters reveals that the MFD differs in nearby clusters, varies or remains constant in time depending on the cluster, and features an unexpected b-value increase in the cluster where the largest event will occur. These findings suggest a strong influence of the heterogeneity and complexity of tectonic structures on the MFD. Our findings raise the question of the appropriate spatiotemporal scale for resolving the b-value, which poses a serious obstacle to interpreting and using the MFD in earthquake forecasting.


2021 ◽  
Author(s):  
Rodrigo Estay ◽  
Claudia Pavez

<p>The Gutenberg – Richter’s b-value is commonly used to analyze the frequency-magnitude distribution of earthquakes, describing the proportion of small and large seismic events as the first estimation of seismic hazard. Additionally, the b-value has been used as a stress meter, giving some insights into the stress regime in different regions around the world. In this research, a grid-based spatial distribution for the b – value was estimated in three different areas of Norway: northern (74°-81° N/ 12°-26° E), southern (57°-64°N/3°-12° E), and the ridge zones of Mohns and Knipovich. For this, we used a complete catalog from the years 2000 to 2019, which was obtained from the Norwegian National Seismic Network online database. The magnitude of completeness was estimated separately for each zone both in time and space, covering a total area of ~425,000 km<sup>2</sup>. Our results show a regional variation of the mean b-value for northern (b<sub>north</sub> = 0.79) and southern (b<sub>south</sub> = 1.03) Norway, and the Ridge (b<sub>ridge</sub> = 0.73), which can be interpreted in terms of the predominant stress regime in the different zones. So far, a few calculations regarding the b-value were previously done in Norway to analyze local intraplate sequences. Then, according to our knowledge, this research corresponds to the first estimation of a regional spatial variation of the b – value in the country.</p>


1980 ◽  
Vol 70 (5) ◽  
pp. 1607-1633
Author(s):  
Christopher D. Stephens ◽  
John C. Lahr ◽  
Kent A. Fogleman ◽  
Robert B. Horner

abstract The St. Elias, Alaska, earthquake (Ms 7.1) of February 28, 1979 occurred beneath the Chugach and St. Elias Mountains of southeastern Alaska and southwestern Yukon Territory. The main shock and aftershocks were recorded at regional high-gain, high-frequency seismographs operated by the U.S. and Canada. Hypocenters and magnitudes are presented for 308 aftershocks that occurred between February 28 and March 31, 1979. These data contain a nearly complete record of events of magnitude 3.5 and larger starting about 20 min after the main shock. The largest aftershock has a poorly determined magnitude slightly above 5, and the frequency-magnitude distribution has a b value of 1.36. A t−p inverse power law with an unusually low value of 0.93 for p adequately describes the decay with time in the frequency of occurrence of large aftershocks. The aftershocks occurred in a broad zone that extends about 115 km southeast from the epicenter of the main shock. Events tend to form clusters within this zone. One of the most remarkable features in the distribution of epicenters is that relatively few aftershocks were located near the epicenter of the main shock, and the highest rate of activity was centered about 50 km southeast of the epicenter of the main shock. Within the accuracy of the data, the depths of the aftershocks are all less than about 20 km. In the few areas where good depth control is available, the seismicity appears not to extend to the Earth's surface. Additional data from temporary stations operated in the aftershock zone during July and August 1979 indicate that the seismicity in some areas may be confined to a zone less than 6 km in vertical thickness. Focal mechanisms determined from p-wave first motions for some of the larger aftershocks all indicate northward-directed compression, which is consistent with the focal mechanism of the main shock. A review of the regional seismicity during the 6-month period preceding the St. Elias earthquake indicates that, relative to a comparable 6-month period 1 yr earlier, there was a 45 per cent increase in the rate of activity for events of magnitude 1.8 and larger, and possibly a decrease in the b value during the same period. Also, a prominent cluster of events with magnitudes less than 3.3 occurred at the southeast corner of the aftershock zone during the 6 months prior to the earthquake. The seismic record from the USGS network since 1974 is not yet complete in time, so it is not possible to determine how unusual the seismic activity preceding this earthquake has been.


Author(s):  
Alireza Babaie Mahani

Critical analysis of induced earthquake occurrences requires comprehensive datasets obtained by dense seismographic networks. In this study, using such datasets, I take a detailed investigation into induced seismicity that occurred in the Montney play of northeast British Columbia, mostly caused by hydraulic fracturing. The frequency-magnitude distribution (FMD) of earthquakes in several temporal and spatial clusters, show fundamental discrepancies between seismicity in the southern Montney play (2014-2018) and the northern area (2014-2016). In both regions, FMDs follow the linear Gutenberg-Richter (G-R) relationship for magnitudes up to 2-3. While in the southern Montney, within the Fort St. John graben complex, the number of earthquakes at larger magnitudes falls off rapidly below the G-R line, within the northern area with a dominant compressional regime, the number of events increases above the G-R line. This systematic difference may have important implications with regard to seismic hazard assessments from induced seismicity in the two regions, although caution in the interpretation is warranted due to local variabilities. While for most clusters within the southern Montney area, the linear or truncated G-R relationship provide reliable seismicity rates for events below magnitude 4, the G-R relationship underestimates the seismicity rate for magnitudes above 3 in northern Montney. Using a well-located dataset of earthquakes in southern Montney, one can observe generally that 1) seismic productivity correlates well with the injected volume during hydraulic fracturing and 2) there is a clear depth dependence for the G-R b-value; clusters with deeper median depths show lower b-values than those with shallower depths.


Author(s):  
Leila Mizrahi ◽  
Shyam Nandan ◽  
Stefan Wiemer

Abstract Declustering aims to divide earthquake catalogs into independent events (mainshocks), and dependent (clustered) events, and is an integral component of many seismicity studies, including seismic hazard assessment. We assess the effect of declustering on the frequency–magnitude distribution of mainshocks. In particular, we examine the dependence of the b-value of declustered catalogs on the choice of declustering approach and algorithm-specific parameters. Using the catalog of earthquakes in California since 1980, we show that the b-value decreases by up to 30% due to declustering with respect to the undeclustered catalog. The extent of the reduction is highly dependent on the declustering method and parameters applied. We then reproduce a similar effect by declustering synthetic earthquake catalogs with known b-value, which have been generated using an epidemic-type aftershock sequence model. Our analysis suggests that the observed decrease in b-value must, at least partially, arise from the application of the declustering algorithm on the catalog, rather than from differences in the nature of mainshocks versus fore- or aftershocks. We conclude that declustering should be considered as a potential source of bias in seismicity and hazard studies.


Sign in / Sign up

Export Citation Format

Share Document