scholarly journals The East Aegean Sea strong earthquake sequence of October–November 2005: lessons learned for earthquake prediction from foreshocks

2006 ◽  
Vol 6 (6) ◽  
pp. 895-901 ◽  
Author(s):  
G. A. Papadopoulos ◽  
I. Latoussakis ◽  
E. Daskalaki ◽  
G. Diakogianni ◽  
A. Fokaefs ◽  
...  
Keyword(s):  
Seismic Analysis ◽  
Strong Shock ◽  
Aegean Sea ◽  
B Value ◽  
Aftershock Sequence ◽  
Lessons Learned ◽  
Time Interval ◽  
Aftershock Sequences ◽  
Duration Magnitude ◽  
Foreshock Activity

Abstract. The seismic sequence of October–November 2005 in the Samos area, East Aegean Sea, was studied with the aim to show how it is possible to establish criteria for (a) the rapid recognition of both the ongoing foreshock activity and the mainshock, and (b) the rapid discrimination between the foreshock and aftershock phases of activity. It has been shown that before the mainshock of 20 October 2005, foreshock activity is not recognizable in the standard earthquake catalogue. However, a detailed examination of the records in the SMG station, which is the closest to the activated area, revealed that hundreds of small shocks not listed in the standard catalogue were recorded in the time interval from 12 October 2005 up to 21 November 2005. The production of reliable relations between seismic signal duration and duration magnitude for earthquakes included in the standard catalogue, made it possible to use signal durations in SMG records and to determine duration magnitudes for 2054 small shocks not included in the standard catalogue. In this way a new catalogue with magnitude determination for 3027 events was obtained while the standard catalogue contains 1025 events. At least 55 of them occurred from 12 October 2005 up to the occurrence of the two strong foreshocks of 17 October 2005. This implies that foreshock activity developed a few days before the strong shocks of 17 October 2005 but it escaped recognition by the routine procedure of seismic analysis. The onset of the foreshock phase of activity is recognizable by the significant increase of the mean seismicity rate which increased exponentially with time. According to the least-squares approach the b-value of the magnitude-frequency relation dropped significantly during the foreshock activity with respect to the b-value prevailing in the declustered background seismicity. However, the maximum likelihood approach does not indicate such a drop of b. The b-value found for the aftershocks that followed the strong shock of 20 October 2005 is significantly higher than in foreshocks. The significant aftershock-foreshock difference in b-value is valid not only if the entire aftershock sequence is considered but also if only the segment of aftershocks that occurred within the first 24-h or the first 48-h after the mainshock of 20 October 2005 are taken into account. This difference in b-value should be examined further in other foreshock-aftershock sequences because it could be used as a diagnostic of the mainshock occurrence within a few hours after its generation.

Operational Earthquake Forecasting during the 2019 Ridgecrest, California, Earthquake Sequence with the UCERF3-ETAS Model

2020 ◽  
Vol 91 (3) ◽  
pp. 1567-1578 ◽  
Author(s):  
Kevin R. Milner ◽  
Edward H. Field ◽  
William H. Savran ◽  
Morgan T. Page ◽  
Thomas H. Jordan
Keyword(s):  
High Performance ◽  
Aftershock Sequence ◽  
Lessons Learned ◽  
Ease Of Use ◽  
Earthquake Forecasting ◽  
Earthquake Rupture ◽  
Etas Model ◽  
Epidemic Type Aftershock Sequence ◽  
Aftershock Sequences ◽  
High Performance Computing Cluster

Abstract The first Uniform California Earthquake Rupture Forecast, Version 3–epidemic-type aftershock sequence (UCERF3-ETAS) aftershock simulations were running on a high-performance computing cluster within 33 min of the 4 July 2019 M 6.4 Searles Valley earthquake. UCERF3-ETAS, an extension of the third Uniform California Earthquake Rupture Forecast (UCERF3), is the first comprehensive, fault-based, epidemic-type aftershock sequence (ETAS) model. It produces ensembles of synthetic aftershock sequences both on and off explicitly modeled UCERF3 faults to answer a key question repeatedly asked during the Ridgecrest sequence: What are the chances that the earthquake that just occurred will turn out to be the foreshock of an even bigger event? As the sequence unfolded—including one such larger event, the 5 July 2019 M 7.1 Ridgecrest earthquake almost 34 hr later—we updated the model with observed aftershocks, finite-rupture estimates, sequence-specific parameters, and alternative UCERF3-ETAS variants. Although configuring and running UCERF3-ETAS at the time of the earthquake was not fully automated, considerable effort had been focused in 2018 on improving model documentation and ease of use with a public GitHub repository, command line tools, and flexible configuration files. These efforts allowed us to quickly respond and efficiently configure new simulations as the sequence evolved. Here, we discuss lessons learned during the Ridgecrest sequence, including sensitivities of fault triggering probabilities to poorly constrained finite-rupture estimates and model assumptions, as well as implications for UCERF3-ETAS operationalization.

Two Foreshock Sequences Post Gulia and Wiemer (2019)

2020 ◽  
Vol 91 (5) ◽  
pp. 2843-2850 ◽  
Author(s):  
Kelian Dascher-Cousineau ◽  
Thorne Lay ◽  
Emily E. Brodsky
Keyword(s):  
Puerto Rico ◽  
Real Time ◽  
Exact Result ◽  
Monte Carlo Sampling ◽  
B Value ◽  
Aftershock Sequence ◽  
Traffic Light ◽  
B Values ◽  
Aftershock Sequences ◽  
Abrupt Changes

Abstract Recognizing earthquakes as foreshocks in real time would provide a valuable forecasting capability. In a recent study, Gulia and Wiemer (2019) proposed a traffic-light system that relies on abrupt changes in b-values relative to background values. The approach utilizes high-resolution earthquake catalogs to monitor localized regions around the largest events and distinguish foreshock sequences (reduced b-values) from aftershock sequences (increased b-values). The recent well-recorded earthquake foreshock sequences in Ridgecrest, California, and Maria Antonia, Puerto Rico, provide an opportunity to test the procedure. For Ridgecrest, our b-value time series indicates an elevated risk of a larger impending earthquake during the Mw 6.4 foreshock sequence and provides an ambiguous identification of the onset of the Mw 7.1 aftershock sequence. However, the exact result depends strongly on expert judgment. Monte Carlo sampling across a range of reasonable decisions most often results in ambiguous warning levels. In the case of the Puerto Rico sequence, we record significant drops in b-value prior to and following the largest event (Mw 6.4) in the sequence. The b-value has still not returned to background levels (12 February 2020). The Ridgecrest sequence roughly conforms to expectations; the Puerto Rico sequence will only do so if a larger event occurs in the future with an ensuing b-value increase. Any real-time implementation of this approach will require dense instrumentation, consistent (versioned) low completeness catalogs, well-calibrated maps of regionalized background b-values, systematic real-time catalog production, and robust decision making about the event source volumes to analyze.

The foreshock-aftershock sequence of the March 20 1966 earthquake in the Republic of Congo

1970 ◽  
Vol 60 (4) ◽  
pp. 1245-1258 ◽  
Author(s):  
John Lahr ◽  
Paul W. Pomeroy
Keyword(s):  
Dynamic Range ◽  
B Value ◽  
High Gain ◽  
Aftershock Sequence ◽  
Confidence Limits ◽  
B Values ◽  
Aftershock Sequences ◽  
Successful Prediction ◽  
Seismograph Station ◽  
The Republic

abstract The activity associated with the Congo earthquake of March 20 1966 (mb = 6.5 to 7) was studied with emphasis on the time and magnitude distributions. The data were recorded at the Abéché, Chad, seismograph station operated by Lamont-Doherty Geological Observatory. Over a period of about 70 days, 815 earthquakes with magnitude (mb) greater than or equal to 3.3 were recorded, and they form the basis for this study. The aftershocks are distributed with magnitude (mb) according to the formula long n = a - bm with b = 1.05 ± 0.07 at the 95 per cent confidence limits. The foreshocks have b = 1.06 ± 0.35 at the 95 per cent confidence limits. These b values are in general agreement with b values derived from other aftershock sequences throughout the world. Some authors have suggested that foreshocks may have a lower b value than background activity and that this difference might be used in earthquake prediction. In this paper, an evaluation is made of the limitations of this method of prediction. Assuming that such a difference in b values does exist, it is found that a closely spaced network of high-gain seismographs with wide dynamic range would be required to assure successful prediction.

scholarly journals Statistical analysis of aftershock sequences related with two major Nepal earthquakes: April 25, 2015, MW 7.8, and May 12, 2015, MW 7.2

2016 ◽  
Vol 59 (5) ◽  
Author(s):  
Prasanta Chingtham ◽  
Babita Sharma ◽  
Sumer Chopra ◽  
Pareshnath SinghaRoy
Keyword(s):  
Scaling Laws ◽  
Main Shock ◽  
B Value ◽  
Temporal Variations ◽  
Aftershock Sequence ◽  
Earthquake Catalog ◽  
Study Region ◽  
Nepal Himalaya ◽  
B Values ◽  
Aftershock Sequences

Present study describes the statistical properties of aftershock sequences related with two major Nepal earthquakes (April 25, 2015, MW 7.8, and May 12, 2015, MW 7.2) and their correlations with the tectonics of Nepal Himalaya. The established empirical scaling laws such as the Gutenberg–Richter (GR) relation, the modified Omori law, and the fractal dimension for both the aftershock sequences of Nepal earthquakes have been investigated to assess the spatio-temporal characteristics of these sequences. For this purpose, the homogenized earthquake catalog in moment magnitude, MW is compiled from International Seismological Center (ISC) and Global Centroid Moment Tensor (GCMT) databases during the period from April 25 to October 31, 2015. The magnitude of completeness, MC, a and b-values of Gutenberg–Richter relationship for the first aftershock sequence are found to be 3.0, 4.74, 0.75 (±0.03) respectively whereas the MC, a and b-values of the same relationship for the second aftershock sequence are calculated to be 3.3, 5.46, 0.90 (±0.04) respectively. The observed low b-values for both the sequences, as compared to the global mean of 1.0 indicate the presence of high differential stress accumulations within the fractured rock mass of Nepal Himalaya. The calculated p-values of 1.01 ± 0.05 and 0.95 ± 0.04 respectively for both the aftershock sequences also imply that the aftershock sequence of first main-shock exhibits relatively faster temporal decay pattern than the aftershock sequence of second main-shock. The fractal dimensions, DC values of 1.84 ± 0.05 and 1.91 ± 0.05 respectively for both the aftershock sequences of Nepal earthquakes also reveal the clustering pattern of earthquakes and signifies that the aftershocks are scattered all around the two dimensional space of fractured fault systems of the Nepal region. The low b-value and low DC observed in the temporal variations of b-value and DC before the investigated earthquake (MW 7.2) suggest the presence of high-stress concentrations in the thrusting regimes of the Nepal region before the failure of faults. Moreover, the decrease of b-value with the corresponding decrease of DC observed in their temporal variations can primarily act as an indicator for possible prediction of major earthquakes in the study region.

scholarly journals Temporal characteristics of some aftershock sequences in Bulgaria

1999 ◽  
Vol 42 (5) ◽  
Author(s):  
S. Simeonova ◽  
D. Solakov
Keyword(s):  
Maximum Likelihood ◽  
Temporal Distribution ◽  
Aftershock Sequence ◽  
Maximum Likelihood Estimates ◽  
Likelihood Method ◽  
Cumulative Number ◽  
Aftershock Activity ◽  
Aftershock Sequences ◽  
Secondary Aftershock ◽  
Foreshock Activity

We apply statistical analysis to study the temporal distribution of aftershocks in aftershock sequences of five earthquakes which occurred in Bulgaria. We use the maximum likelihood method to estimate the parameters of the modified Omori formula for aftershock sequences which is directly based on a time series. We find that: the maximum likelihood estimates of the parameter p show a regional variation, with lower values of the decay rate in North Bulgaria; the modified Omori formula provides an appropriate representation of temporal variation of the aftershock activity in North Bulgaria; the aftershock sequences in South Bulgaria are best modeled by the combination of an ordinary aftershock sequence with secondary aftershock activity. A plot of the cumulative number of events versus the frequency-linearized time t clearly demonstrates a transition from aftershock to foreshock activity prior to the second 1986 Strazhitsa (North Bulgaria) earthquake.

Earthquake migration in the Fairbanks, Alaska seismic zone

1980 ◽  
Vol 70 (1) ◽  
pp. 223-241
Author(s):  
Larry Gedney ◽  
Steve Estes ◽  
Nirendra Biswas
Keyword(s):  
B Value ◽  
Aftershock Sequence ◽  
Stress Axis ◽  
Seismic Zone ◽  
Epicentral Area ◽  
Focal Mechanism Solutions ◽  
Relative Stress ◽  
Earthquake Migration ◽  
Dislocation Process ◽  
High Level

abstract Since a series of moderate earthquakes near Fairbanks, Alaska in 1967, the “Fairbanks seismic zone” has maintained a consistently high level of seismicity interspersed with sporadic earthquake swarms. Five swarms occurring since 1970 demonstrate that tightly compacted centers of activity have tended to migrate away from the epicentral area of the 1967 earthquakes. Comparative b-coefficients of the first four swarms indicate that they occurred under different relative stress conditions than the last episode, which exhibited a higher b-value and was, in fact, a main shock of magnitude 4.6 with a rapidly decaying aftershock sequence. This last recorded sequence in February 1979 was an extension to greater depths along a lineal seismic zone whose first recorded activation occurred during a swarm two years earlier. Focal mechanism solutions indicate a north-south orientation of the greatest principal stress axis, σ1, in the area. A dislocation process related to crustal spreading between strands of a right-lateral fault, similar to that which has been inferred for southern California, is suggested.

Rampart seismic zone of central Alaska

1983 ◽  
Vol 73 (3) ◽  
pp. 813-829
Author(s):  
P. Yi-Fa Huang ◽  
N. N. Biswas
Keyword(s):  
Main Shock ◽  
Fault Plane ◽  
B Value ◽  
Aftershock Sequence ◽  
Strike Slip ◽  
Lithospheric Plate ◽  
Normal Faults ◽  
Seismic Zone ◽  
The West ◽  
Fault Plane Solutions

abstract This paper describes the characteristics of the Rampart seismic zone by means of the aftershock sequence of the Rampart earthquake (ML = 6.8) which occurred in central Alaska on 29 October 1968. The magnitudes of the aftershocks ranged from about 1.6 to 4.4 which yielded a b value of 0.96 ± 0.09. The locations of the aftershocks outline a NNE-SSW trending aftershock zone about 50 km long which coincides with the offset of the Kaltag fault from the Victoria Creek fault. The rupture zone dips steeply (≈80°) to the west and extends from the surface to a depth of about 10 km. Fault plane solutions for a group of selected aftershocks, which occurred over a period of 22 days after the main shock, show simultaneous occurrences of strike-slip and normal faults. A comparison of the trends in seismicity between the neighboring areas shows that the Rampart seismic zone lies outside the area of underthrusting of the lithospheric plate in southcentral and central Alaska. The seismic zone outlined by the aftershock sequence appears to represent the formation of an intraplate fracture caused by regional northwest compression.

Observations of the Coyote Lake, California earthquake sequence of August 6, 1979

1980 ◽  
Vol 70 (2) ◽  
pp. 559-570 ◽  
Author(s):  
R. A. Uhrhammer
Keyword(s):  
San Francisco ◽  
Fault Zone ◽  
Strong Earthquake ◽  
Main Shock ◽  
B Value ◽  
Aftershock Sequence ◽  
The South ◽  
The North ◽  
Temporal And Spatial Characteristics ◽  
Largest Aftershock

abstract At 1705 UTC on August 6, 1979, a strong earthquake (ML = 5.9) occurred along the Calaveras fault zone south of Coyote Lake about 110 km southeast of San Francisco. This strong earthquake had an aftershock sequence of 31 events (2.4 ≦ ML ≦ 4.4) during August 1979. No foreshocks (ML ≧ 1.5) were observed in the 3 months prior to the main shock. The local magnitude (ML = 5.9) and the seismic moment (Mo = 6 × 1024 dyne-cm from the SH pulse) for the main shock were determined from the 100 × torsion and 3-component ultra-long period seismographs located at Berkeley. Local magnitudes are determined for the aftershocks from the maximum trace amplitudes on the Wood-Anderson torsion seismograms recorded at Berkeley (Δ ≊ 110 km). Temporal and spatial characteristics of the aftershock sequence are presented and discussed. Some key observations are: (1) the first six aftershocks (ML ≧ 2.4) proceed along the fault zone progressively to the south of the main shock; (2) all of the aftershocks (ML ≧ 2.4) to the south of the largest aftershock (ML = 4.4) have a different focal mechanism than the aftershocks to the north; (3) no aftershocks (ML ≧ 2.4) were observed significantly to the north of the main shock for the first 5 days of the sequence; and (4) the b-value (0.70 ± 0.17) for the aftershock sequence is not significantly different from the average b-value (0.88 ± 0.08) calculated for the Calaveras fault zone from 16 yr of data.

scholarly journals Dependences of the Omori and Gutenberg–Richter parameters

2019 ◽  
pp. 149-165 ◽  
Author(s):  
V. B. Smirnov ◽  
A. V. Ponomarev ◽  
S. A. Stanchits ◽  
M. G. Potanina ◽  
A. V. Patonin ◽  
...  
Keyword(s):  
Time Delay ◽  
Stress State ◽  
Axial Loading ◽  
B Value ◽  
Exponential Dependence ◽  
Parameter Estimates ◽  
Aftershock Activity ◽  
Uniform Compression ◽  
Effective Strength ◽  
Aftershock Sequences

Laboratory experiments on studying the aftershock regime are carried out with sandstone specimens under different axial loading and uniform compression and constant pore pressure. The aftershock sequences are modeled by the scenario of stepwise increasing axial loading of a specimen with strain control ensuring regular generation of aftershock sequences. The experiments are conducted on intact specimens and on the specimens with preliminarily formed shear macrofractures simulating natural faults. The experiments were conducted with multichannel recording of the acoustic emission (AE) signals which made it possible to locate the AE sources. Several types of the dependence of the acoustic activity relaxation parameters (parameters p and c of the modified Omori law and the Gutenberg–Richter b-value) on the level of acting stresses are revealed. The b-value decreases with the growth of axial stresses at all levels of uniform compression. In the case of fracture on the preexisting fault, the Omori relaxation parameter p increases with the growth of axial stresses whereas parameter c (the time delay before the onset of relaxation) decreases with the growth of axial stresses and increases with the rise of the level of uniform compression. In the case of a fracture of an undamaged specimen, parameter p remains unchanged as the axial stresses grow, whereas parameter c increases slightly. Parameter variations in the case of a complex stress state with both varying deviatoric (differential stresses) and spherical parts (effective pressure) of the stress tensor take on a unified form when expressed in terms of Coulomb stresses. It is hypothesized that the time delay of the aftershock activity relaxation is determined by the kinetics of fracture in accordance with the kinetic concept of strength in solids. This hypothesis is supported by exponential dependence of parameter c on stresses and on the effective strength of the medium revealed in the experiments. Under this hypothesis, the dependences of parameter c on the Coulomb stresses can be unified for different effective strength values with the use of Zhurkov’s formula for durability of materials. The obtained parameter estimates for the dependence of c on strength and stresses suggest that the c value is determined by the difference of the strength and the acting stresses, indicating how far the stress state of the medium is from the critical state corresponding to the ultimate strength.

scholarly journals Modeling active fault systems and seismic events by using a fiber bundle model – example case: the Northridge aftershock sequence

Solid Earth ◽  
2019 ◽  
Vol 10 (5) ◽  
pp. 1519-1540
Author(s):  
Marisol Monterrubio-Velasco ◽  
F. Ramón Zúñiga ◽  
José Carlos Carrasco-Jiménez ◽  
Víctor Márquez-Ramírez ◽  
Josep de la Puente
Keyword(s):  
Fiber Bundle ◽  
Active Faults ◽  
Aftershock Sequence ◽  
Statistical Characteristics ◽  
Real Sequence ◽  
Random Load ◽  
Fault Systems ◽  
Aftershock Sequences ◽  
Fiber Bundle Model ◽  
Initial Load

Abstract. Earthquake aftershocks display spatiotemporal correlations arising from their self-organized critical behavior. Dynamic deterministic modeling of aftershock series is challenging to carry out due to both the physical complexity and uncertainties related to the different parameters which govern the system. Nevertheless, numerical simulations with the help of stochastic models such as the fiber bundle model (FBM) allow the use of an analog of the physical model that produces a statistical behavior with many similarities to real series. FBMs are simple discrete element models that can be characterized by using few parameters. In this work, the aim is to present a new model based on FBM that includes geometrical characteristics of fault systems. In our model, the faults are not described with typical geometric measures such as dip, strike, and slip, but they are incorporated as weak regions in the model domain that could increase the likelihood to generate earthquakes. In order to analyze the sensitivity of the model to input parameters, a parametric study is carried out. Our analysis focuses on aftershock statistics in space, time, and magnitude domains. Moreover, we analyzed the synthetic aftershock sequences properties assuming initial load configurations and suitable conditions to propagate the rupture. As an example case, we have modeled a set of real active faults related to the Northridge, California, earthquake sequence. We compare the simulation results to statistical characteristics from the Northridge sequence determining which range of parameters in our FBM version reproduces the main features observed in real aftershock series. From the results obtained, we observe that two parameters related to the initial load configuration are determinant in obtaining realistic seismicity characteristics: (1) parameter P, which represents the initial probability order, and (2) parameter π, which is the percentage of load distributed to the neighboring cells. The results show that in order to reproduce statistical characteristics of the real sequence, larger πfrac values (0.85<πfrac<0.95) and very low values of P (0.0<P≤0.08) are needed. This implies the important corollary that a very small departure from an initial random load configuration (computed by P), and also a large difference between the load transfer from on-fault segments than by off-faults (computed by πfrac), is required to initiate a rupture sequence which conforms to observed statistical properties such as the Gutenberg–Richter law, Omori law, and fractal dimension.

Sign in / Sign up

Export Citation Format

Share Document