skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, February 13 until 2:00 AM ET on Friday, February 14 due to maintenance. We apologize for the inconvenience.


Title: Analysis and Modeling of High-frequency Emission and Deep Seismic Sources of Sunquakes
Abstract

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 betweenz= 0 andz= −840 km. Analysis of the suspected wave front of the M9.3 sunquake finds that the power spectrum of the reported seismic emission is close to that of background oscillations, though with a peak frequency noticeably lower than the background peak. Additionally, it is found that the amplitude of the source estimated by Lindsey et al. is too low to produce the observed wave front.

 
more » « less
Award ID(s):
1916509
PAR ID:
10372352
Author(s) / Creator(s):
;
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
  1. SUMMARY

    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.

     
    more » « less
  2.  
    more » « less
  3. 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
  4. 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
  5. 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