skip to main content


Title: Long Short-Term Memory Networks for Pattern Recognition of Synthetical Complete Earthquake Catalog
Exploring the spatiotemporal distribution of earthquake activity, especially earthquake migration of fault systems, can greatly to understand the basic mechanics of earthquakes and the assessment of earthquake risk. By establishing a three-dimensional strike-slip fault model, to derive the stress response and fault slip along the fault under regional stress conditions. Our study helps to create a long-term, complete earthquake catalog. We modelled Long-Short Term Memory (LSTM) networks for pattern recognition of the synthetical earthquake catalog. The performance of the models was compared using the mean-square error (MSE). Our results showed clearly the application of LSTM showed a meaningful result of 0.08% in the MSE values. Our best model can predict the time and magnitude of the earthquakes with a magnitude greater than Mw = 6.5 with a similar clustering period. These results showed conclusively that applying LSTM in a spatiotemporal series prediction provides a potential application in the study of earthquake mechanics and forecasting of major earthquake events.  more » « less
Award ID(s):
1918126
PAR ID:
10253866
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Sustainability
Volume:
13
Issue:
9
ISSN:
2071-1050
Page Range / eLocation ID:
4905
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Faults are usually surrounded by damage zones associated with localized deformation. Here we use fully dynamic earthquake cycle simulations to quantify the behaviors of earthquakes in fault damage zones. We show that fault damage zones can make a significant contribution to the spatial and temporal seismicity distribution. Fault stress heterogeneities generated by fault zone waves persist over multiple earthquake cycles that, in turn, produce small earthquakes that are absent in homogeneous simulations with the same friction conditions. Shallow fault zones can produce a bimodal depth distribution of earthquakes with clustering of seismicity at both shallower and deeper depths. Fault zone healing during the interseismic period also promotes the penetration of aseismic slip into the locked region and reduces the sizes of fault asperities that host earthquakes. Hence, small and moderate subsurface earthquakes with irregular recurrence intervals are commonly observed in immature fault zone simulations with interseismic healing. To link our simulation results to geological observations, we will use simulated fault slip at different depths to infer the timing and recurrence intervals of earthquakes and discuss how such measurements can affect our understanding of earthquake behaviors. We will also show that the maturity and material properties of fault damage zones have strong influence on whether long-term earthquake characteristics are represented by single events. 
    more » « less
  2. ABSTRACT Although the Brune source model describes earthquake moment release as a single pulse, it is widely used in studies of complex earthquakes with multiple episodes of high moment release (i.e., multiple subevents). In this study, we investigate how corner frequency estimates of earthquakes with multiple subevents are biased if they are based on the Brune source model. By assuming complex sources as a sum of multiple Brune sources, we analyze 1640 source time functions of Mw 5.5–8.0 earthquakes in the seismic source characteristic retrieved from deconvolving teleseismic body waves catalog to estimate the corner frequencies, onset times, and seismic moments of subevents. We identify more subevents for strike-slip earthquakes than dip-slip earthquakes, and the number of resolvable subevents increases with magnitude. We find that earthquake corner frequency correlates best with the corner frequency of the subevent with the highest moment release (i.e., the largest subsevent). This suggests that, when the Brune model is used, the estimated corner frequency and, therefore, the stress drop of a complex earthquake is determined primarily by the largest subevent rather than the total rupture area. Our results imply that, in addition to the simplified assumption of a radial rupture area with a constant rupture velocity, the stress variation of asperities, rather than the average stress change of the whole fault, contributes to the large variance of stress-drop estimates. 
    more » « less
  3. Abstract

    Earthquake stress drop is an important source parameter that directly links to strong ground motion and fundamental questions in earthquake physics. Stress drop estimations may contain significant uncertainties due to such factors as variations in material properties and data limitations, which limit the applications of stress drop interpretations. Using a high‐resolution borehole network, we estimate stress drop for 4551 (M0‐4) earthquakes on the San Andreas Fault at Parkfield, California, between 2001 and 2016 using spectral decomposition and an improved stacking method. To evaluate the influence of spatiotemporal variations of material properties on stress drop estimations, we apply different strategies to account for spatial variations of velocity and attenuation changes, and divide earthquakes into three separate time periods to correct temporal variations of attenuation. These results show that appropriate corrections can significantly reduce the scatter in stress drop estimates, and decrease apparent depth and magnitude dependence. We find that insufficient bandwidth can cause systematic underestimation of stress drop estimates and increased scatter. The stress drop measurements from the high‐frequency borehole recordings exhibit complex stable spatial patterns with no clear correlation with the nature of fault slip, or the slip distribution of the 2004 M6 earthquake. Temporal variations are significantly smaller, less well resolved and varying spatially. They do not affect the long‐term stress drop spatial variations, suggesting local material properties may control the spatial heterogeneity of stress drop.

     
    more » « less
  4. null (Ed.)
    SUMMARY Earthquake ruptures are generally considered to be cracks that propagate as fracture or frictional slip on pre-existing faults. Crack models have been used to describe the spatial distribution of fault offset and the associated static stress changes along a fault, and have implications for friction evolution and the underlying physics of rupture processes. However, field measurements that could help refine idealized crack models are rare. Here, we describe large-scale laboratory earthquake experiments, where all rupture processes were contained within a 3-m long saw-cut granite fault, and we propose an analytical crack model that fits our measurements. Similar to natural earthquakes, laboratory measurements show coseismic slip that gradually tapers near the rupture tips. Measured stress changes show roughly constant stress drop in the centre of the ruptured region, a maximum stress increase near the rupture tips and a smooth transition in between, in a region we describe as the earthquake arrest zone. The proposed model generalizes the widely used elliptical crack model by adding gradually tapered slip at the ends of the rupture. Different from the cohesive zone described by fracture mechanics, we propose that the transition in stress changes and the corresponding linear taper observed in the earthquake arrest zone are the result of rupture termination conditions primarily controlled by the initial stress distribution. It is the heterogeneous initial stress distribution that controls the arrest of laboratory earthquakes, and the features of static stress changes. We also performed dynamic rupture simulations that confirm how arrest conditions can affect slip taper and static stress changes. If applicable to larger natural earthquakes, this distinction between an earthquake arrest zone (that depends on stress conditions) and a cohesive zone (that depends primarily on strength evolution) has important implications for how seismic observations of earthquake fracture energy should be interpreted. 
    more » « less
  5. The 2019 Ridgecrest, CA earthquake sequence has provided a unique opportunity and a rich dataset to understand earthquake source properties and near-fault structure. Using the high-quality seismic data provided by the SCEC Stress Drop Validation group, we first estimate the corner frequency of M2.0-4.5 earthquakes by applying the spectral ratio method based on empirical Green’s function (Liu et al., 2020). We relate corner frequency estimates to stress drops assuming the Brune source model and circular cracks. Our preliminary results show increasing median stress drops with magnitude for both P and S waves, from 1 MPa for M2.0 events to 10 MPa for M4.0 events, though the limited frequency bandwidth may cause underestimation for small events. The estimated moment magnitude is proportional to the catalog magnitude by a factor of 0.72, which is close to 0.74 estimated by Trugman (2020) for the Ridgecrest earthquake sequence. In the second part of the study, we examine the impact of fault zone structure on the azimuthal variation of the source spectra. Using kinematic simulations and observations of the 2003 Big Bear earthquake sequence, Huang et al. (2016) showed that fault damage zones can act as an effective wave guide and cause high-frequency wave amplification along directions close to fault strike. We use clusters of M1.5-3 earthquakes in the Ridgecrest region to further examine the azimuthal variation of the stacked source spectra and investigate if the near-source structure can affect our corner frequency estimates. We aim to develop robust methods that utilize high-quality seismic data to illuminate earthquake source processes and fault zone properties. 
    more » « less