Earthquake Source Research Articles

Reactivated Seismogenic Faults and Earthquake Source Mechanisms in the Snyder Area of Texas

Abstract Stretching across New Mexico and Texas of the United States, the greater Permian basin is composed of two subunits—the Delaware and the Midland basins. Induced seismicity in the greater Permian basin has significantly increased since 2008, which has revealed previously unmapped seismogenic structures in several geographic regions. Among them, the Snyder area of northwest Texas has a long history of oil and gas activities, resulting in a higher rate of induced seismicity. In this study, we investigated these previously unknown seismogenic structures using three main approaches: (1) relocated and delineated seismicity, (2) performed waveform moment tensor inversion to determine earthquake source mechanisms, as well as (3) conducted stress inversion to assess the stress state. The results show that the overall depth range of seismicity is 0–5.5 km and concentrated in a range of 2–3 km below mean sea level, in the top portion of the crystalline basement. As we have determined 297 source mechanisms, their collective pattern presents a mix of strike-slip and normal faulting, suggesting an extensional strain field at the edge of the Midland basin. We have identified nine significant seismogenic episodes by distinctive increases of seismic moment release in 2017–March 2024. The results also demonstrate a temporal variation of b-value spanning across the seismogenic episodes, associated with the progression of fault reactivation initiated by fluid injection.

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

AbstractWe develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Use of decision tree ensembles for crustal structure imaging from receiver functions

SUMMARY Mode conversion of P waves at the boundary between Earth's crust and upper mantle, when analysed using receiver functions (RFs), allows characterization of Earth structure where seismic station density is high and earthquake sources are favourably distributed. We applied two ensemble decision tree algorithms—Random Forest (RanFor) and eXtreme Gradient Boost (XGBoost)—to synthetic and real RF data to assess these machine learning techniques' potential for crustal imaging when available data are sparse. The synthetic RFs, entailing both sharp increases in seismic velocity across the Moho and gradational Moho structures, calculated with and without added random noise, correspond to idealized crustal structures: a dipping Moho, Moho offset by crustal-scale faults, anti- and synform Moho structures and combinations of these. The RanFor/XGBoost algorithm recovers input structures well regardless of event-station distributions. Useful crustal and upper mantle seismic velocities can also be determined using RanFor and XGBoost, making it possible to image crustal thickness and P- and S-wave velocities simultaneously from RFs alone. We applied the trained RanFor/XGBoost to RFs determined from real seismic data recorded in the contiguous United States, producing a map of the Moho and P- and S-wave velocities of the lowermost crust and uppermost mantle. Use of XGBoost, which evaluates residuals between input RFs and ground-truth to update the decision tree using the gradient of a penalty function, improves the crustal thickness estimates.

Massively parallel Bayesian estimation with Sequential Monte Carlo sampling for simultaneous estimation of earthquake fault geometry and slip distribution

In inverse analysis, Bayesian estimation is useful for understanding the reliability of the estimation result or prediction based on it because it can estimate not only optimal parameters but also their uncertainties. The estimation of an earthquake source fault based on observed crustal deformation is a typical inverse problem in the field of earthquake research. In this study, a method for simultaneous Bayesian estimation of earthquake fault plane geometry and spatially variable slip distribution on the plane has been developed. The developed method can be applied to stochastic models with arbitrary probability distribution settings, and it enables to incorporate appropriate constraints for slip distribution in the estimation process, which can lead to enhanced robustness and stability of estimation. Since this method is computationally more expensive than conventional methods, large-scale parallel computing was introduced to cope with the increased computational cost and a supercomputer was used for the analysis. To validate the proposed method, simultaneous Bayesian estimation of fault geometry and slip distribution with slip constraints was performed using the crustal deformation observed in the 2018 Hokkaido Eastern Iburi earthquake. Hierarchical parameterization and massively parallelized Bayesian inference used in this study have broad applicability not only in earthquake research but also in other scientific and engineering fields.

Uncovering a Seismogenic Fault in Southern Iran through Co-Seismic Deformation of the Mw 6.1 Doublet Earthquake of 14 November 2021

On 14 November 2021, a doublet earthquake, each event of which had an Mw of 6.1, struck near Fin in the Simply Folded Belt (SFB) in southern Iran. The first quake occurred at 12:07:04 UTC, followed by a second one just a minute and a half later. The SFB is known for its blind thrust faults, typically not associated with surface ruptures. These earthquakes are usually linked to the middle and lower layers of the sedimentary cover. Identifying the faults that trigger earthquakes in the region remains a significant challenge and is subject to high uncertainty. This study aims to identify and determine the fault(s) that may have caused the doublet earthquake. To achieve this goal, we utilized the DInSAR method using Sentinel-1 to detect deformation, followed by finite-fault inversion and magnetic interpretation to determine the location, geometry, and slip distribution of the fault(s). Bayesian probabilistic joint inversion was used to model the earthquake sources and derive the geometric parameters of potential fault planes. The study presents two potential fault solutions—one dipping to the north and the other to the south. Both solutions showed no significant difference in strike and fault location, suggesting a single fault. Based on the results of the seismic inversion, it appears that a north-dipping fault with a strike, dip, and rake of 257°, 74°, and 77°, respectively, is more consistent with the geological setting of the area. The fault plane has a width of roughly 3.6 km, a length of 13.4 km, and a depth of 5.6 km. Our results revealed maximum displacements along the radar line of sight reaching values of up to −360 mm in the ascending orbit, indicating an unknown fault with horizontal displacements at the surface ranging from −144 to 170 mm and maximum vertical displacements between −204 and 415 mm. Aeromagnetic data for Iran were utilized with an average flight-line spacing of 7.5 km. The middle of the data observation period was considered to apply the RTP filter, and the DRTP method was used. We calculated the gradient of the residual anomaly in the N-S direction due to the direction of the existing faults and folds. The gradient map identified the fault and potential extension of the observed anomalies related to a fault with an ENE-WSW strike, which could extend to the ~ E-W. We suggest that earthquakes occur in the sedimentary cover of the SFB where subsurface faulting is involved, with Hormuz salt acting as an important barrier to rupture. The multidisciplinary approach used in this study, including InSAR and magnetic data, underscores the importance of accurate fault characterization. These findings provide valuable insights into the seismic hazard of the area.

Earthquake on January 24, 2024 near Krasnodar city, Mw=4.1, I0=5

The article presents instrumental and macroseismic data on the earthquake on January 24, 2024 at 11:19 (UTC), Mw=4.1, h=8 km. The earthquake parameters were determined using instrumental data from the network of regional seismic stations in the North Caucasus of the GS RAS. The earthquake parameters were determined from instrumental data from the network of regional seismic stations in the North Caucasus of the GS RAS. The source spectra of three regional seismic stations (GOYR, MRNR and GUZR) were calculated, from which the seismic moment М0, the stress drop Δσ and the source radius R were determined. The solution of the focal mechanism was obtained from the polarization in P-waves at 71 seismic stations and the type of source fault with a right-lateral strike-slip component. The epicentral zone of the earthquake on 24.01.2024 is located within the high-magnitude Akhtyrka zone of potential earthquake source zones (ESZ) with Mmax=6.5-6.8 and represents the boundary between the Northwestern part of the structures of the Greater Caucasus and the West Kuban trough. According to macroseismic survey in the regions, 300 respondents were interviewed in 20 settlements. The maximum observed intensity was I=V MSK-64.

Prediction of Regional Broadband Strong Ground Motions Using a Teleseismic Source Model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, Earthquake

ABSTRACT To estimate predicted ground motion from a teleseismic slip model, we use a low- and high-frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude (Mw) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite-fault, kinematic model is used to define the earthquake source, taking into account slip-model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low-frequency ground motions using a finite-element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full-wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low- and high-frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation-independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong-motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground-motion prediction relations utilized in the assessment of seismic hazard.

Seismotectonics-Driven Estimation of b-Value: Implications for Seismic Hazard Assessment

Abstract Taking full advantage of good-quality historical earthquake and seismogenic source data available for Italy (Catalogo Parametrico dei Terremoti Italiani [CPTI15] and Database of Individual Seismogenic Sources [DISS] databases), we tested the regional variability of the b-value of the Gutenberg–Richter law, focusing on the dominantly extensional, nearly 1200-km-long seismogenic corridor straddling the Apennines chain; an 18,200 km2 area generating about 45% of the seismic moment released countrywide. We carefully chose and tested the most appropriate completeness interval and completeness magnitude (Mc), and used the Lilliefors method, the Utsu test, and a bootstrap technique to verify the quality of our results and inferences. The 0.65 b-value we obtained is substantially lower than b-values used in Italian seismic hazard models, often close to 1.0; yet it predicts more accurately the observed Apennines earthquake record, suggesting that most currently adopted magnitude-frequency distributions underpredict events in the Mw range 6.0–7.1 and overpredict those in the range 5.2–5.5. We also highlight that the time, space, and magnitude distribution of the largest Italian earthquakes hardly follows the Gutenberg–Richter law: it rather favors a characteristic behavior. We contend that the b-value is too critical of a parameter for being calculated using statistically weak data sets, such as those resulting from overly detailed area-source models: one may end up characterizing low-hazard areas more accurately than high-hazard areas. Conversely, we advocate the use of fewer, larger area sources that fully exploit the current understanding of regional seismotectonics, as a way of obtaining more realistic regional hazard models.

Simulation of large plausible tsunami scenarios associated with the 2019 Durres (Albania) earthquake source and adjacent seismogenic zones

We present an analysis of the hazards of potential earthquake-generated tsunamis along the Albanian–Adriatic coast. The study adopts a case study approach to model plausible tsunamigenic events associated with the 2019 Mw 6.4 Durres (Albania) earthquake source zone. The approach combines current findings on regional tectonics and scenario-based calculations of potential tsunami impacts. The study’s goal is to analyse the propagation of tsunami waves generated by identified seismogenic sources (namely ALCS002 [Lushnje] and ALCS018 [Shijak]) and determine the tsunami risk assessment for Durres City on the Albanian–Adriatic coast. The sources can generate earthquakes with maximum moment magnitudes of Mw 7.5 and Mw 6.8, which are likely to trigger tsunamis that could cause significant impacts in the region. The modelling is performed deterministically with the NAMI DANCE numerical code, including scenarios associated with the largest plausible earthquake. The model integrates bathymetry and topography datasets of large and medium resolutions. Each tsunami scenario simulation is based on the solution of the non-linear shallow water equations used to generate maximum positive wave amplitudes (water elevation), travel time, and tsunami inundation maps. In Durres City, modelling indicates that medium-sized waves could reach up to 2.5 m inland, posing a significant danger to the city’s low-lying areas. The most substantial tsunami waves are expected to impact the area within the first 10 to 20 min. Combining inundation maps and information on exposed assets allows for identifying areas where damages can be expected. In terms of human impact, a preliminary analysis shows that the study area is prone to tsunami threat, with more than 138,000 inhabitants living in vulnerable urban areas of Durres City by 2036. The model’s capacity to capture details related to the presence of buildings is limited due to constraints posed by the resolution of bathymetry and topography datasets available during this study. If refined with high-resolution bathymetry and topography datasets, our results can be considered a backbone for exposure and resilience assessment features to be integrated into preparedness or new urban development plans.

High‐Frequency Ground Motions of Earthquakes Correlate With Fault Network Complexity

AbstractUnderstanding the generation of damaging, high‐frequency ground motions during earthquakes is essential both for fundamental science and for effective hazard preparation. Various theories exist regarding the origin of high‐frequency ground motions, including the standard paradigm linked to slip heterogeneity on the rupture plane, and alternative perspectives associated with fault complexity. To assess these competing hypotheses, we measure ground motion amplitudes in different frequency bands for 3 ≤ M ≤ 5.8 earthquakes in Southern California and compare them to empirical ground motion models. We utilize a Bayesian inference technique called the Integrated Nested Laplace Approximation (INLA) to identify earthquake source regions that produce higher or lower ground motions than expected. Our analysis reveals a strong correlation between fault complexity measurements and the high‐frequency ground motion event terms identified by INLA. These findings suggest that earthquakes on complex faults (or fault networks) lead to stronger‐than‐expected ground motions at high frequencies.

Evidence for active faults along the Sinop trough close to the Mid Black Sea Ridge Using Multi-Channel Seismic and Multi-Beam Bathymetric Data

Abstract Focal mechanism solutions of the earthquakes that occurred in the Black Sea in the last 100 years show that thrust and strike-slip faults are predominantly active in the region. One of these clusters is developing off the coast of Samsun where the submarine Sinop basin is located. In order to investigate the seismic sources of earthquakes in this area, 14 high-resolution multi-channel seismic sections and multi-beam bathymetric data were processed and evaluated. The oldest seismic stratigrafic units (U4) is referred to as coustic basement with wavy reflectors, while its top is marked by high ampitude reflection indicating an unconformity. This uniconformity is overlain by parallel or less parallel Unit 3 deposits. Upper surface of Unit 2 is marked by an erosional surface which is overlain by parallel reflector of Unit 1. Unit 1 is the youngest sediments and truncated by Yeşilırmak canyon. Three fault types of different ages were determined. The oldest fault could reach only top of unit 4 has been interpreted as inactive fault. Faults that border the Sinop basin and reach the sea floor by cutting all seismic units are considered in the active fault class. However, the faults that remain within Unit 1 and do not reach the seafloor are considered as faults that are coeval with sedimentation. According to these data, faulting in the Central/Eastern Pontide structural block should be reconsidered in the context of the seismic activity of the region and re-evaluated as an earthquake hazard that will pose a risk in the future.

Cold Diffusion Model for Seismic Denoising

AbstractSeismic waves contain information about the earthquake source, the geologic structure they traverse, and many forms of noise. Separating the noise from the earthquake is a difficult task because optimal parameters for filtering noise typically vary with time and, if chosen inappropriately, may strongly alter the original seismic waveform. Diffusion models based on Deep Learning have demonstrated remarkable capabilities in restoring images and audio signals. However, those models assume a Gaussian distribution of noise, which is not the case for typical seismic noise. Motivated by the effectiveness of “cold” diffusion models in speech enhancement, medical anomaly detection, and image restoration, we present a cold variant for seismic data restoration. We describe the first Cold Diffusion Model for Seismic Denoising (CDiffSD), including key design aspects, model architecture, and noise handling. Using metrics to quantify the performance of CDiffSD models compared to previous works, we demonstrate that it provides a new standard in performance. CDiffSD significantly improved the Signal to Noise Ratio by about 18% compared to previous models. It also enhanced Cross‐correlation by 6%, showing a better match between denoised and original signals. Moreover, testing revealed a 50% increase in the recall of P‐wave picks for seismic picking. Our work show that CDiffSD outperforms existing benchmarks, further underscoring its effectiveness in seismic data denoising and analysis. Additionally, the versatility of this model suggests its potential applicability across a range of tasks and domains, such as GNSS, Lab Acoustic Emission, and Distributed Acoustic Sensing data, offering promising avenues for further utilization.

Seismotectonics of the Philippine and Taiwan Subduction Systems and Implications for Seismic Hazards

AbstractThe seismicity of the Philippines and Taiwan provides insight into the tectonics and seismic hazards of a region characterized by subduction and collision. We summarize the seismotectonics of the Philippines and Taiwan by documenting the distribution of hypocenters for earthquakes of magnitude M ≥ 4.6 and focal mechanisms for earthquakes of magnitude M ≥ 5.0 over ∼21 years. We quantify seismicity rates (earthquake frequency) and compare seismicity distributions with proposed tectonic and faulting models. 6,187 earthquakes of magnitude M ≥ 4.6 occurred between 1 January 2000 and 31 March 2021, 79% at depths <70 km and 70% having magnitudes M < 5.0. Approximately 88 earthquakes of magnitude M ≥ 5.0 occur per year, with 12 events of magnitude M ≥ 7.0 occurring since January 2000. Seismic activity decreases exponentially between 50 and 210 km depth at a rate ∼10% faster than the global average. Intermediate and deep earthquakes at depths >70 km trace the Wadati‐Benioff zones of subducting slabs, most of which are only seismically active to depths of ∼250 km. The distribution of earthquakes at depths >70 km is likely influenced by the subduction of young lithosphere, slab tearing, and phase boundary interactions between depths of 410 and 660 km. Shallow earthquakes at depths ≤70 km are generated by megathrust, crustal, and intraslab faulting. Crustal thrust and strike‐slip faulting are the most abundant and prevalent sources of damaging earthquakes. The Philippines and Taiwan are subject to high seismic risk, similar to nearby Indonesia.

Patterns of Causative Faults of Normal Earthquakes in the Fluid‐Rich Outer Rise of Northeastern Japan, Constrained With 3D Teleseismic Waveform Modeling

AbstractAccurate earthquake source parameters are crucial for understanding plate tectonics, yet, it is difficult to determine these parameters precisely for offshore events, especially for outer‐rise earthquakes, as the limited availability of direct P or S wave data sets from land‐based seismic networks and the unsuitability of simplified 1D methods for the complex 3D structures of subducting systems. To overcome these challenges, we employ an efficient hybrid numerical simulation method to model these 3D structural effects on teleseismic P/SH and P‐coda waves and determine the reliable centroid locations and focal mechanisms of outer‐rise normal‐faulting earthquakes in northeastern Japan. Two M6+ events with reliable locations from ocean bottom seismic observations are utilized to calibrate the 3D velocity structure. Our findings indicate that 3D synthetic waveforms are sensitive to both event location, thanks to bathymetry and water reverberation effects, and the shallow portion of the lithospheric structure. With our preferred velocity model, which has Vs ∼16% lower than the global average, event locations are determined with uncertainties of <5 km for horizontal position and <1 km for depth. The refined event locations in a good match between one of the nodal strikes and the high‐resolution bathymetry, enabling the determination of the causative fault plane. Our results reveal that trench‐ward dipping normal faults are more active, with three parallel to the trench as expected, while five are associated with the abyssal hills. The significant velocity reduction in the uppermost lithosphere suggests abundant water migrating through active normal faults, enhancing both mineral alteration and pore density.

Analysis of the Stress Field Around Concealed Active Fault From Minor Faults‐Slip Data Collected by Geological Survey: An Example in the 1984 Western Nagano Earthquake Region

AbstractEarthquakes with magnitudes of 6–7 have been reported even in various active tectonic settings where fault deformation topography have not been detected. Therefore, delineating concealed active faults generating such earthquakes is necessary to reduce earthquake damage; however, few studies exist to provide its clues regarding such faults. The 1984 Western Nagano Earthquake in Japan was a main shock with a magnitude of Mj 6.8 and a depth of 2 km at the source. Solid bedrocks are well‐exposed in the earthquake source region; however, no surface rupture have been identified, and the active fault is known to be concealed. In this study, we collected data on striations observed in fractures by geological survey around the source area of the 1984 Western Nagano Earthquake. Using the collected data, the multiple inverse method was used to estimate the stresses that affected the striation formation. Consequently, stresses similar to acting faults in this area were detected in minor faults around the known concealed active fault. This suggests that the minor faults might be part of the damage zone that has been developed around the concealed active fault. Some minor faults were recognized in Quaternary volcanic rocks, confirming that they experienced displacements recently. This study indicates the possibility of detecting concealed active faults in the bedrock by geological survey.

Site-specific probabilistic seismic hazard analysis in Nuweiba, Gulf of Aqaba, Egypt: Combining area source model and active faults

Gulf of Aqaba in Eastern Egypt is characterized by significant seismic activity and high deformation rate, making it a major source of destructive earthquakes in the region. While previous studies have produced probabilistic seismic hazard analysis (PSHA) maps for the region, they have not adequately accounted for active faults or site response in their calculations. This study introduces a newly developed seismic source model for the Gulf of Aqaba, considering area source, active fault source, and a combination of both models to improve seismic hazard assessment for the region. The study incorporates the previously measured average shear-wave velocity for the upper 30 m (Vs30), fundamental frequency (F0), and amplification factor (A0) at some selected sites in Nuweiba city. Accordingly, site-specific PSHA is carried out for the bedrock condition (Vs30 = 800 m/s, according to the Egyptian Building Code) and for the surface of soft sediments using F0, A0, and Vs30. Uncertainties in slip rate, maximum magnitude, and recurrence model are accounted during PSHA computation. Ground motion parameters such as peak ground acceleration (PGA) and spectral acceleration (SA) for various probabilities of exceedance in 50-year lifetime are calculated for both rock site conditions and soft sediments. The maximum calculated PGA and SA are 0.65 g and 1.25 g, respectively, for a 10% probability of exceedance over 50 years at bedrock conditions.

Earthquake source impacts on the generation and propagation of seismic infrasound to the upper atmosphere

SUMMARY Earthquakes with moment magnitude (Mw) ranging from 6.5 to 7.0 have been observed to generate sufficiently strong acoustic waves (AWs) in the upper atmosphere. These AWs are detectable in Global Navigation Satellite System satellite signals-based total electron content (TEC) observations in the ionosphere at altitudes ∼250–300 km. However, the specific earthquake source parameters that influence the detectability and characteristics of AWs are not comprehensively understood. Here, we extend our approach of coupled earthquake-atmosphere dynamics modelling by combing dynamic rupture and seismic wave propagation simulations with 2-D and 3-D atmospheric numerical models, to investigate how the characteristics of earthquakes impact the generation and propagation of AWs. We developed a set of idealized dynamic rupture models varying faulting types and fault sizes, hypocentral depths and stress drops. We focus on earthquakes of Mw 6.0–6.5, which are considered the smallest detectable with TEC, and find that the resulting AWs undergo non-linear evolution and form acoustic shock N waves reaching thermosphere at ∼90–140 km. The results reveal that the magnitude of the earthquakes is not the sole or primary factor determining the amplitudes of AWs in the upper atmosphere. Instead, various earthquake source characteristics, including the direction of rupture propagation, the polarity of seismic wave imprints on the surface, earthquake mechanism, stress drop and radiated energy, significantly influence the amplitudes and periods of AWs. The simulation results are also compared with observed TEC fluctuations from AWs generated by the 2023 Mw 6.2 Suzu (Japan) earthquake, finding preliminary agreement in terms of model-predicted signal periods and amplitudes. Understanding these nuanced relationships between earthquake source parameters and AW characteristics is essential for refining our ability to detect and interpret AW signals in the ionosphere.

The Impact of the Three-Dimensional Structure of a Subduction Zone on Time-dependent Crustal Deformation Measured by HR-GNSS

Accurately modeling time-dependent coseismic crustal deformation as observed on high-rate Global Navigation Satellite System (HR-GNSS) lends insight into earthquake source processes and improves local earthquake and tsunami early warning algorithms. Currently, time-dependent crustal deformation modeling relies most frequently on simplified 1D radially symmetric Earth models. However, for shallow subduction zone earthquakes, even low-frequency shaking is likely affected by the many strongly heterogeneous structures such as the subducting slab, mantle wedge, and the overlying crustal structure. We demonstrate that including 3D structure improves the estimation of key features of coseismic HR-GNSS time series, such as the peak ground displacement (PGD), the time to PGD (tPGD), static displacements (SD), and waveform cross-correlation values. We computed synthetic 1D and 3D, 0.25 Hz and 0.5 Hz waveforms at HR-GNSS stations for four M7.3+ earthquakes in Japan using MudPy and SW4, respectively. From these synthetics, we computed intensity-measure residuals between the synthetic and observed GNSS waveforms. Comparing 1D and 3D residuals, we observed that the 3D simulations show better fits to the PGD and SD in the observed waveforms than the 1D simulations for both 0.25 Hz and 0.5 Hz simulations. We find that the reduction in PGD residuals in the 3D simulations is a combined effect of both shallow and deep 3D structures; hence incorporating only the upper 30 km of 3D structure will still improve the fit to the observed PGD values. Our results demonstrate that 3D simulations significantly improve models of GNSS waveform characteristics and will not only help understand the underlying processes, but also improve local tsunami warning.

Numerical Simulation of Electromagnetic Responses to an Earthquake Source Due To the Piezoelectric Effect of ‘∞m’ Symmetry

AbstractIt is reported that earthquakes are accompanied by electromagnetic (EM) anomalies. These anomalies are thought to be caused by earthquakes but their generation mechanism is still unclear. The piezoelectric effect has been proposed as a possible mechanism, particularly in quartz‐rich rocks. However, the EM responses to earthquakes due to such an effect have not been well understood. In this article, we study the EM signals generated by an earthquake source due to the piezoelectric effect. We develop a semi‐analytical method to simulate the EM responses to earthquake sources in a 3‐D layered model and conduct numerical simulations to investigate the characteristics of the EM fields. The results show that the piezoelectric effect in quartz‐rich rocks could be the responsible mechanism for the EM signals recorded during the earthquake. Two kinds of EM signals can be generated by the earthquake, namely, the early EM wave arriving before the seismic waves and the coseismic EM field accompanying the seismic waves. For an Mw 6.1 earthquake, the coseismic electric field is ∼0.1 μV/m and the magnetic field can only reach ∼10−5 nT. We also study the sensitivity of the coseismic EM fields to the rock conductivity. The results indicate that the coseismic EM fields are mainly affected by the conductivity of the shallow layer where the receiver is located, and if the top layer is thin, they are also affected by the conductivity of the deeper layer due to the contribution of the evanescent EM wave created at the interface.

Low‐Frequency Earthquakes Downdip of Deep Slow Slip Beneath the North Island of New Zealand

AbstractWe report the first catalog of low‐frequency earthquakes in the Hikurangi subduction zone, located beneath the Kaimanawa Range of the North Island at 50 km depth, downdip of regularly recurring (every 4–5 years) deep M7 slow slip events. To systematically detect low‐frequency earthquakes within the regional continuous seismic data, we utilized a matched‐filter approach with template waveforms derived from previous observations of tectonic tremor. We built our catalog of 36 low‐frequency earthquake sources, that produced almost 21,000 events over more than a decade, with two matched‐filter search iterations. In each iteration, the detections were gathered into families and their coherent waveforms processed and stacked to extract high‐quality waveforms, allowing us to pick seismic phase arrivals to locate the low‐frequency earthquakes. We highlight three characteristic features to validate that our detected events are indeed low‐frequency earthquakes: the eponymous deficit of high frequencies in their seismic waveforms, the episodic swarms of activity that define their activity through time, and their location at the plate boundary with a double‐couple source mechanism and geometry consistent with the subduction interface. Considering the observed low‐frequency earthquakes' relationship to neighboring slow slip, we observe the event swarms to occur much more frequently than the M7 slow slip events located just updip. Similar to other deep low‐frequency earthquakes in other subduction zones, we suggest that this characteristic clustering in time is driven by more frequent, smaller slow slip events that are not clearly observable at the surface.