skip to main content

Title: Radiated seismic energy and source damage evolution from the analysis of simulated dynamic rupture and far-field seismograms

We present results on radiated seismic energy during simulations of dynamic ruptures in a continuum damage-breakage rheological model incorporating evolution of damage within the seismic source region. The simulations vary in their initial damage zone width and rate of damage diffusion with parameter values constrained by observational data. The radiated energy recorded at various positions around the source is used to calculate seismic potency and moment. We also calculate the normalized radiated energy from the source, in a way that allows comparing between results of different simulations and highlighting aspects related to the dilatational motion during rupture. The results show that at high-frequencies, beyond the dominant frequency of the source ($( {f > 3{f}_d} )$, the damage process produces an additional burst of energy mainly in the Pwaves. This eccess of high-frequency energy is observed by comparing the radiated energy to a standard Brune's model with a decay slope of the radiated energy of n = 2. While the Swaves show good agreement with the n = 2 slope, the Pwaves have a milder slope of n = 1.75 or less depending on the damage evolution at the source. In the used damage-breakage rheology, the rate of damage diffusivity governs the damage evolution perpendicular to the rupture direction and dynamic changes of the damage zone width. For increasing values of damage diffusivity, dilatational energy becomes more prominent during rupture, producing a high-frequency dilatational signature within the radiation pattern. The high-frequency radiation pattern of the Pwaves includes two main lobes perpendicular to the rupture direction, reflecting high-rate local tensile cracking during the overall shear rupture process. Analysing the possible existence and properties of such high-frequency radiation pattern in observed Pwaves could provide important information on earthquake source processes.

more » « less
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Medium: X Size: p. 1705-1726
p. 1705-1726
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Frictional heating during earthquake rupture raises the fault-zone fluid pressure, which affects dynamic rupture and seismic radiation. Here, we investigate two key parameters governing thermal pressurization of pore fluids – hydraulic diffusivity and shear-zone half-width – and their effects on earthquake rupture dynamics, kinematic source properties, and ground motions. We conduct 3D strike-slip dynamic rupture simulations assuming a rate-and-state dependent friction law with strong velocity weakening coupled to thermal-pressurization of pore fluids. Dynamic rupture evolution and ground shaking are densely evaluated across the fault and Earth’s surface to analyze the variations of rupture parameters (slip, peak slip rate, rupture speed, and rise time), correlations among rupture parameters, and variability of peak ground velocity. Our simulations reveal how variations in thermal-pressurization affect earthquake rupture properties. We find that the mean slip and rise time decrease with increasing hydraulic diffusivity, whereas mean rupture speed and peak slip-rate remain almost constant. Mean slip, peak slip-rate, and rupture speed decrease with increasing shear-zone half-width, whereas mean rise time increases. Shear-zone half-width distinctly affects the correlation between rupture parameters, especially for parameter pairs (slip, rupture speed), (peak slip-rate, rupture speed), and (rupture speed, rise time). Hydraulic diffusivity has negligible effects on these correlations. Variations in shear-zone half-width primarily impact rupture speed, which then may affect other rupture parameters. We find a negative correlation between slip and peak slip-rate, unlike simpler dynamic rupture models. Mean peak ground velocities decrease faster with increasing shear-zone half-width than with increasing hydraulic diffusivity, whereas ground-motion variability is similarly affected by both the parameters. Our results show that shear-zone half-width affects rupture dynamics, kinematic rupture properties, and ground shaking more strongly than hydraulic diffusivity. We interpret the importance of shear-zone half-width based on the characteristic time of diffusion. Our findings may inform pseudodynamic rupture generators and guide future studies on how to account for thermal-pressurization effects. 
    more » « less
  2. Abstract

    Earthquake radiation is broadband, and its temporal evolution carries the seismic signature of the rupture process. I use established and develop new observational metrics to quantify the temporal variation in radiated energy (seismic power) and in a time‐dependent scaled radiated energy. Universal functional forms of moment rate, seismic power, and scaled energy arise from a global database of source time functions. Earthquakes radiate mostly and most efficiently in the early stage of development of the rupture that is in the first 10–30% of the overall rupture duration. This result is independent of earthquake magnitude, depth, and focal mechanism. Beyond depth dependence in elastic properties with depth that explain variations in source duration and radiated energy, subtle differences in functional forms exist between shallow and deep earthquakes. Deep events have a slower development than shallow earthquakes, but they carry a higher peak in seismic power.

    more » « less
  3. Abstract

    Megathrust earthquakes exhibit a ubiquitous seismic radiation style: low‐frequency (LF) seismic energy is efficiently emitted from the shallowest portion of the fault, whereas high‐frequency (HF) seismic energy is efficiently emitted from the deepest part of the fault. Although this is observed in many case‐specific studies, we show that it is ubiquitous in global megathrust earthquakes between 1995 and 2021. Previous studies have interpreted this as an effect of systematic depth variation in either the plate interface frictional properties (Lay et al., 2012) or the P wavespeeds (Sallarès & Ranero, 2019). This work suggests an alternative hypothesis: the interaction between waves and ruptures due to the Earth's free surface is the leading mechanism that generates this behavior. Two‐dimensional dynamic rupture simulations of subduction zone earthquakes support this hypothesis. Our simulations show that the interaction between the seismic waves reflected at the Earth's free surface and the updip propagating rupture results in LF radiation at the source. In contrast, the downdip propagation of rupture is less affected by the free surface and is thus dominated by HF radiation typical of buried faults. To a second degree, the presence of a realistic Earth structure derived from P‐wave velocity (VP) tomographic images and realistic VP/VSratio estimated in boreholes further enhances the contrast in source radiation. We conclude that the Earth's free surface is necessary to explain the observed megathrust earthquake radiation style, and the realistic structure of subduction zone is necessary to better predict earthquake ground motion and tsunami potential.

    more » « less
  4. Fault damage zones can influence various aspects of the earthquake cycle, such as the recurrence intervals and magnitudes of large earthquakes. The properties and structure of fault damage zones are often characterized using dense arrays of seismic stations located directly above the faults. However, such arrays may not always be available. Hence, our research aims to develop a novel method to image fault damage zones using broadband stations at relatively larger distances. Previous kinematic simulations and a case study of the 2003 Big Bear earthquake sequence demonstrated that fault damage zones can act as effective waveguides, amplifying high-frequency waves along directions close to fault strike via multiple reflections within the fault damage zone. The amplified high-frequency energy can be observed by stacking P-wave spectra of earthquake clusters with highly-similar waveforms (Huang et al., 2016), and the frequency band which is amplified may be used to estimate the width and velocity contrast of the fault damage zone. We attempt to identify the high-frequency peak associated with fault zone waves in stacked spectra by conducting a large-scale study of small earthquakes (M1.5–3). We use high quality broadband data from seismic stations at hypocentral distances of 20-80 km in the 2019 Ridgecrest earthquake regions. First, we group the Ridgecrest earthquakes in clusters by their locations and their waveform similarity, and then stack their velocity spectra to average the source effects of individual earthquakes. Our results show that the stations close to the fault strike record more high-frequency energies around the characteristic frequency of fault zone reflections. We find that the increase in the amount of high-frequencies is consistent across clusters with average magnitudes ranging from 1.6-2.4, which suggests that the azimuthal variation in spectra is caused by fault zone amplification rather than rupture directivity. We will apply our method to other fault zones in California, in order to search for fault damage zone structures and estimate their material properties. 
    more » « less
  5. ABSTRACT A depletion of high-frequency ground motions on soil sites has been observed in recent large earthquakes and is often attributed to a nonlinear soil response. Here, I show that the reduced amplitudes of high-frequency horizontal-to-vertical spectral ratios (HVSRs) on soil can also be caused by a smooth crustal velocity model with low shear-wave velocities underneath soil sites. I calculate near-fault ground motions using both 2D dynamic rupture simulations and point-source models for both rock and soil sites. The 1D velocity models used in the simulations are derived from empirical relationships between seismic wave velocities and depths in northern California. The simulations for soil sites feature lower shear-wave velocities and thus larger Poisson’s ratios at shallow depths than those for rock sites. The lower shear-wave velocities cause slower shallow rupture and smaller shallow slip, but both soil and rock simulations have similar rupture speeds and slip for the rest of the fault. However, the simulated near-fault ground motions on soil and rock sites have distinct features. Compared to ground motions on rock, horizontal ground acceleration on soil is only amplified at low frequencies, whereas vertical ground acceleration is deamplified for the whole frequency range. Thus, the HVSRs on soil exhibit a depletion of high-frequency energy. The comparison between smooth and layered velocity models demonstrates that the smoothness of the velocity model plays a critical role in the contrasting behaviors of HVSRs on soil and rock for different rupture styles and velocity profiles. The results reveal the significant role of shallow crustal velocity structure in the generation of high-frequency ground motions on soil sites. 
    more » « less