Precisely constraining the source parameters of large earthquakes is one of the primary objectives of seismology. However, the quality of the results relies on the quality of synthetic earth response. Although earth structure is laterally heterogeneous, particularly at shallow depth, most earthquake source studies at the global scale rely on the Green's functions calculated with radially symmetric (1-D) earth structure. To avoid the impact of inaccurate Green's functions, these conventional source studies use a limited set of seismic phases, such as long-period seismic waves, broad-band P and S waves in teleseismic distances (30° < ∆ < 90°), and strong ground motion records at close-fault stations. The enriched information embedded in the broad-band seismograms recorded by global and regional networks is largely ignored, limiting the spatiotemporal resolution. Here we calculate 3-D strain Green's functions at 30 GSN stations for source regions of 9 selected global earthquakes and one earthquake-prone area (California), with frequency up to 67 mHz (15 s), using SPECFEM3D_GLOBE and the reciprocity theorem. The 3-D SEM mesh model is composed of mantle model S40RTS, crustal model CRUST2.0 and surface topography ETOPO2. We surround each target event with grids in horizontal spacing of 5 km and vertical spacing of 2.0–3.0 km, allowing us to investigate not only the main shock but also the background seismicity. In total, the response at over 210 000 source points is calculated in simulation. The number of earthquakes, including different focal mechanisms, centroid depth range and tectonic background, could further increase without additional computational cost if they were properly selected to avoid overloading individual CPUs. The storage requirement can be reduced by two orders of magnitude if the output strain Green's functions are stored for periods over 15 s. We quantitatively evaluate the quality of these 3-D synthetic seismograms, which are frequency and phase dependent, for each source region using nearby aftershocks, before using them to constrain the focal mechanisms and slip distribution. Case studies show that using a 3-D earth model significantly improves the waveform similarity, agreement in amplitude and arrival time of seismic phases with the observations. The limitations of current 3-D models are still notable, dependent on seismic phases and frequency range. The 3-D synthetic seismograms cannot well match the high frequency (>40 mHz) S wave and (>20 mHz) Rayleigh wave yet. Though the mean time-shifts are close to zero, the standard deviations are notable. Careful calibration using the records of nearby better located earthquakes is still recommended to take full advantage of better waveform similarity due to the use of 3-D models. Our results indicate that it is now feasible to systematically study global large earthquakes using full 3-D earth response in a global scale.
Recent work published by Lindsey et al. find evidence for a deep and compact seismic source for the sunquake associated with the 2011 July 30 M9.3 flare, as well as seismic emission extending up to 10 mHz. We examine the sunquake independently, and a possible wave front is found in the 8 mHz band, though no wave front is easily discernible in the 10 mHz band. Additionally, we perform numerical simulations of seismic excitation modeled with the reported parameters and changes in the power spectra with increasing depth of the excitation source are examined. It is found that the peak frequency decreases for increasing depths, but a shallow minimum is indicated between
- Award ID(s):
- 1916509
- PAR ID:
- 10372352
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 937
- Issue:
- 2
- ISSN:
- 2041-8205
- Format(s):
- Medium: X Size: Article No. L26
- Size(s):
- Article No. L26
- Sponsoring Org:
- National Science Foundation
More Like this
-
SUMMARY -
A dual‐band scandium‐aluminum nitride (Sc0.28Al0.72N) Lamb‐wave bandpass filter is implemented with intrinsically configurable center frequencies at 628 MHz and 1.07 GHz. The filter is created from electrical coupling of dual‐mode resonators with complementary‐switchable operation states defined by spatial polarization tuning in ferroelectric Sc0.28Al0.72N film. The resonators are created with patterning interdigitated transducers (IDT) on the top surface and exploit the floating bottom electrode for application of low‐frequency polarization switching pulses. Switching polarization direction under consecutive electrode fingers in unified and alternating directions results in excitation of Lamb modes with different wavelengths equal and twice the IDT pitch size, with a frequency ratio of ≈1.70. The intrinsically configurable dual‐band filter shows bandwidths of 11.9 MHz and 10.8 MHz, insertion losses of 2.6 and 2.8 dB, and on‐off switching isolations of 12 and 9 dB, when operating at 628 MHz and 1.07 GHz, respectively. The intrinsic configuration of center frequency, enabled by ferroelectric polarization switching, highlights the potential of dual‐band scandium‐aluminum nitride Lamb‐wave filters for creation of compact multifrequency spectral processors.
-
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
-
null (Ed.)Abstract Ice shelves play an important role in buttressing land ice from reaching the sea, thus restraining the rate of grounded ice loss. Long-period gravity-wave impacts excite vibrations in ice shelves that can expand pre-existing fractures and trigger iceberg calving. To investigate the spatial amplitude variability and propagation characteristics of these vibrations, a 34-station broadband seismic array was deployed on the Ross Ice Shelf (RIS) from November 2014 to November 2016. Two types of ice-shelf plate waves were identified with beamforming: flexural-gravity waves and extensional Lamb waves. Below 20 mHz, flexural-gravity waves dominate coherent signals across the array and propagate landward from the ice front at close to shallow-water gravity-wave speeds (~70 m s −1 ). In the 20–100 mHz band, extensional Lamb waves dominate and propagate at phase speeds ~3 km s −1 . Flexural-gravity and extensional Lamb waves were also observed by a 5-station broadband seismic array deployed on the Pine Island Glacier (PIG) ice shelf from January 2012 to December 2013, with flexural wave energy, also detected at the PIG in the 20–100 mHz band. Considering the ubiquitous presence of storm activity in the Southern Ocean and the similar observations at both the RIS and the PIG ice shelves, it is likely that most, if not all, West Antarctic ice shelves are subjected to similar gravity-wave excitation.more » « less
-
ABSTRACT. Ice shelves play an important role in buttressing land ice from reaching the sea, thus restrain- ing the rate of grounded ice loss. Long-period gravity-wave impacts excite vibrations in ice shelves that can expand pre-existing fractures and trigger iceberg calving. To investigate the spatial amplitude vari- ability and propagation characteristics of these vibrations, a 34-station broadband seismic array was deployed on the Ross Ice Shelf (RIS) from November 2014 to November 2016. Two types of ice-shelf plate waves were identified with beamforming: flexural-gravity waves and extensional Lamb waves. Below 20 mHz, flexural-gravity waves dominate coherent signals across the array and propagate land- ward from the ice front at close to shallow-water gravity-wave speeds (∼70 m s−1). In the 20–100 mHz band, extensional Lamb waves dominate and propagate at phase speeds ∼3 km s−1. Flexural- gravity and extensional Lamb waves were also observed by a 5-station broadband seismic array deployed on the Pine Island Glacier (PIG) ice shelf from January 2012 to December 2013, with flexural wave energy, also detected at the PIG in the 20–100 mHz band. Considering the ubiquitous presence of storm activity in the Southern Ocean and the similar observations at both the RIS and the PIG ice shelves, it is likely that most, if not all, West Antarctic ice shelves are subjected to similar gravity- wave excitation.more » « less