Emergence of Low-Frequency Aftershocks of the 2019 Ridgecrest Earthquake Sequence

ABSTRACT Low-frequency earthquakes (LFEs) generally have relatively stronger spectral components in the lower frequency range compared with what is expected for regular earthquakes based on their magnitude. LFEs generally occur in volcanic systems or deep (>∼15 km) in plate boundary fault zones; however, LFEs have also been observed in nonvolcanic, upper crustal settings. Because there are few studies that explore the spatiotemporal behaviors of LFEs in the shallow crust, it remains unclear whether the shallow-crustal LFEs reflect local attenuation in their immediate vicinity or differences in their source mechanism. Therefore, it is important to identify shallow-crustal LFEs and to characterize their spatiotemporal activity, which may also improve our understanding of LFEs. In this study, we focus on detecting shallow-crustal LFEs and explore the possible generation mechanisms. We analyze 29,646 aftershocks in the 2019 Ridgecrest, California, earthquake sequence, by measuring the frequency index (FI) to identify candidate low-frequency aftershocks (LFAs), while accounting for the magnitude dependency of the FI. Using small earthquakes (ML 1–3) recorded in the borehole stations to minimize the attenuation effects in near-surface layers, we identify 68 clear LFAs in total. Based on their distribution and comparisons with other seismic parameters measured by Trugman (2020), the LFAs possess distinct features from regular events in the same depths range, including low corner frequencies and low stress drops. Events in the close vicinity of LFAs exhibit lower average FI values than regular aftershocks, particularly if the hypocentral distance between an LFA and its neighbors is less than 1 km. Our results suggest that LFAs are related to local heterogeneity or a highly fractured fault zone correlated with an abundance of cross faults induced by the aftershock sequence at shallow depths. Zones of high pore-fluid pressure in intensely fractured fault zones could cause the bandlimited nature of LFAs and LFEs in general.

A study of aftershocks of 20 October 1991 Uttarkashi earthquake

The Uttarkashi earthquake of 20 October 1991, which caused widespread damage in the Galhwal Himalayan region, was followed by a prominent aftershock. activity extending over a period of about two months. The aftershock activity was monitored using temporary networks established after the mainshock and the permanent stations in operation in the region. About 142 aftershocks could be located accurately using the data of these stations. The b-value of the Gutenberg-Richter's relationship for the aftershock sequence works out to be 0.6. The temporal distribution of the aftershocks suggests a hyperbolic decay with a decay constant (p) of 1.17. Macroseismic observations derived from field surveys show good agreement with the instrumentally determined source parameters.

Detection of Building Response Changes Using Deconvolution Interferometry: A Case Study in Bogota, Colombia

Seismic structural health monitoring allows for the continuous evaluation of engineering structures by monitoring changes in the structural response that can potentially localize associated damage that has occurred. For the first time in Colombia, a permanent and continuous monitoring network has been deployed in a 14-story ecofriendly steel-frame building combined with a reinforced concrete structure in downtown Bogota. The six three-component ETNA-2 accelerometers recorded continuously for 225 days between July 2019 and February 2020. We use deconvolution-based seismic interferometry to calculate the impulse response function (IRF) using earthquake and ambient-vibration data and a stretching technique to estimate velocity variations before and after the Ml 6.0 Mesetas earthquake and its aftershock sequence. A consistent and probably permanent velocity variation (2% reduction) is detected for the building using ambient-vibration data. In contrast, a 10% velocity reduction is observed just after the mainshock using earthquake-based IRFs showing a quick recovery to about 2%. A combination of both earthquake-based and ambient-vibration-based deconvolution interferometry provides a more complete picture of the state of health of engineering structures.

Combining a Locally Restricted Anisotropic Spatial Kernel with an ETAS Incomplete Model for Better Forecasts of the 2019 Ridgecrest Sequence

Strong earthquakes cause aftershock sequences that are clustered in time according to a power decay law, and in space along their extended rupture, shaping a typically elongate pattern of aftershock locations. A widely used approach to model seismic clustering is the Epidemic Type Aftershock Sequence (ETAS) model, that shows three major biases: First, the conventional ETAS approach assumes isotropic spatial triggering, which stands in conflict with observations and geophysical arguments for strong earthquakes. Second, the spatial kernel has unlimited extent, allowing smaller events to exert disproportionate trigger potential over an unrealistically large area. Third, the ETAS model assumes complete event records and neglects inevitable short-term aftershock incompleteness as a consequence of overlapping coda waves. These three effects can substantially bias the parameter estimation and particularly lead to underestimated cluster sizes. In this article, we combine the approach of Grimm (2021), which introduced a generalized anisotropic and locally restricted spatial kernel, with the ETAS-Incomplete (ETASI) time model of Hainzl (2021), to define an ETASI space-time model with flexible spatial kernel that solves the abovementioned shortcomings. We apply different model versions to a triad of forecasting experiments of the 2019 Ridgecrest sequence, and evaluate the prediction quality with respect to cluster size, largest aftershock magnitude and spatial distribution. The new model provides the potential of more realistic simulations of on-going aftershock activity, e.g.~allowing better predictions of the probability and location of a strong, damaging aftershock, which might be beneficial for short term risk assessment and desaster response.

THE HINDU KUSH EARTHQUAKE on October 26, 2015 with Mw=7.5, I0~7: PRELIMINARY SEISMICITY and AFTERSHOCK SEQUENCE

On October 26, 2015, a strong Hindu Kush earthquake with KR=17.0, Mw=7.5 occurred in the Afghan Pamir-Hindu Kush subzone at a depth of hpP=217 km. Shakes of varying intensity caused by this earthquake were recorded in settlements of 14 states: Afghanistan, Tajikistan, Pakistan, Turkmenistan, India, Kyrgyzstan, Uzbekistan, Kazakhstan, China, Iran, Nepal, United Arab Emirates, Russia, Qatar and Bangladesh with a total area of S=14106 km2. The earthquake was preceded by three large (KR=12.5, 12.1, 14.0) foreshocks and was accompanied by a series of more than 1400 aftershocks unprecedented for aftershocks of deep earthquakes with KR=9–13. The energy step between the mainshock and the maximum foreshock is Kfor=3.0, between the mainshock and the maximum (KR=12.8) aftershock – Kaft=4.2. The aftershock recurrence graph has a slope =–0.67, which in absolute value is higher than the average value in the region =0.50. The attenuation para-meter  of the Omori law in the initial phase of attenuation, =–1.26, in absolute value is also higher than the average =1.0 for strong earthquakes in the World. Based on the results of a joint analysis of the focal me-chanism solutions of different agencies and vertical sections along and across the aftershock cloud, it was con-cluded that an upthrust movement occurred in the source along a steep east-south-east nodal plane, dipping to the south. The reason for the activity at the site of the earthquake is the movement of the Indian continent to the north and its collision with Eurasia, as a result of which the separation and subduction of the Hindu Kush plate continue. The Hindu Kush earthquake on October 26, 2015, and its aftershocks are just one of the events of successive deformation and stress relief in the latitudinal zone, marked in 2015 by the migration of earthquake epicenters with KR=13–17 from east to west.

Seismic Anisotropy and Stress of the Canterbury Plains

<p>The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.</p>

Seismic Anisotropy and Stress of the Canterbury Plains

<p>The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.</p>

Short-Term Foreshock and Aftershock Patterns of the 2021 Ms 6.4 Yangbi Earthquake Sequence

An Ms 6.4 earthquake struck Yangbi County in western Yunnan province, China, on 21 May 2021, causing damage in the nearby region. Intensive foreshock activity started three days before the mainshock, and numerous aftershocks followed along a northwest–southeast-trending right-lateral main rupture fault. Double-difference relocation of the foreshock and aftershock sequence shortly before and after the Ms 6.4 mainshock is conducted using the phase picks from the local seismic network. The focal mechanisms of relatively large foreshocks and aftershocks are also derived. The results not only delineate the ruptured fault geometry during the mainshock but also indicate the mechanism of static stress transfer according to the spatiotemporal evolution of foreshocks. The low background b-values around the mainshock are also consistent with the occurrence of the Yangbi earthquake sequence.

Seismic Assessment and Rehabilitation of a RC Structure Under Mainshock-Aftershock Seismic Sequences Using Beam-Column Bonded CFRP Strategy

A mainshock is usually accompanied by a group of ground motions. In many design codes, the effects of the seismic sequences have been neglected or underestimated. Aftershocks can increase structural damage or even cause failure. The current study evaluated the seismic behavior of a rehabilitated and as-built RC structure under real scaled mainshock-aftershocks using nonlinear analysis. Verification was done in two modes. The inter-story drift ratio, maximum residual, and relative displacements were studied. The seismic study of the as-built structure showed that the residual displacement grew, on average, more than 90% under the mainshock-aftershock sequence compared to the mainshock-only record. A beam-column bonded CFRP rehabilitation strategy using six layers of T-700 CFRP was chosen based on the specific performance level. The strategy showed the ability to transfer the plastic strain from the columns to the beams, which could be considered as a change from a weak column-strong beam concept to a strong column-weak beam concept. Compared to the as-built structure under the seismic sequences, the rehabilitated structure showed an average growth of 78% in the first-story drift ratio, which was significant. Unlike the as-built structure, seismic sequences caused no growth drift in the rehabilitated structure. It also was observed that the ratio of aftershock PGA to mainshock PGA could have an intensive effect on the seismic behavior of both rehabilitated and as-built structures.

Crustal Structure of the Seismogenic Volume of the 2010–2014 Pollino (Italy) Seismic Sequence From 3D P- and S-Wave Tomographic Images

A tomographic analysis of Mt. Pollino area (Italy) has been performed using earthquakes recorded in the area during an intense seismic sequence that occurred between 2010 and 2014. 870 local earthquakes with magnitude ranging from 1.8 to 5.0 were selected considering the number of recording stations, the signal quality, and the hypocenter distribution. P- and S-wave arrival times were manually picked and used to compute 3D velocity models through tomographic seismic inversion. The resulting 3D distributions of VP and VS are characterized by high resolution in the central part of the investigated area and from surface to about 10 km below sea level. The aim of the work is to obtain high-quality tomographic images to correlate with the main lithological units that characterize the study area. The results will be important to enhance the seismic hazard assessment of this complex tectonic region. These images show the ductile Apennine platform (VP = 5.3 km/s) overlaying the brittle Apulian platform (VP = 6.0 km/s) at depth of around 5 km. The central sector of the area shows a clear fold and thrust interface. Along this structure, most of the seismicity occurred, including the strongest event of the sequence (MW 5.0). High VP (&gt;6.8 km/s) and high VP/VS (&gt;1.9) patterns, intersecting the southern edge of this western seismogenic volume, have been interpreted as water saturated rocks, in agreement with similar geological context in the Apennines. These fluids could have played a role in nucleation and development of the seismic sequence. A recent study revealed the occurrence of clusters of earthquakes with similar waveforms along the same seismogenic volume. The hypocenters of these cluster events have been compared with the events re-located in this work. Jointly, they depict a 10 km × 4 km fault plane, NW-SE oriented, deepening towards SW with a dip angle of 40–45°. Instead, the volume of seismicity responsible for the ML 4.3 earthquake developed as a mainshock-aftershock sequence, occurring entirely within the average-to-low VP/VS Apennine platform. Our results agree with other independent geophysical analyses carried out in this area, and they could significantly improve the actual knowledge of the main lithologic units of this complex tectonic area.