skip to main content

This content will become publicly available on November 27, 2024

Title: Constructing a 3‐D Radially Anisotropic Crustal Velocity Model for Oklahoma Using Full Waveform Inversion

Over the past decade, the seismicity rate in the state of Oklahoma has increased significantly, which has been linked to industrial operations, such as saltwater injection and hydraulic fracturing. Taking advantage of induced earthquakes and recently deployed seismometers, we construct a 3‐D radially anisotropic seismic velocity model for the crust of Oklahoma by using full waveform inversion. To mitigate the well‐known cycle‐skipping problem, we use misfit functions based on phase and waveform differences in several frequency bands. Relative velocity perturbations in the inverted model allow us to delineate major geological provinces in Oklahoma, such as the Anadarko Basin and the Cherokee Platform/Shelf. In addition, radial anisotropy in the inverted model reflects deformation within the crust of Oklahoma, which might correlate with sedimentary layering, microcracks/fractures, as well as dominant orientation of anisotropic minerals. The crystalline basement beneath Oklahoma can be inferred from the new velocity model, which enables us to better classify induced seismicity in current earthquake catalogs. Furthermore, synthetic experiments suggest that the new velocity model enables us to better constrain earthquake locations in Oklahoma, especially for determining their depths, which are important for investigating induced seismicity.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    To complement velocity distributions, seismic attenuation provides additional important information on fluid properties of hydrocarbon reservoirs in exploration seismology, as well as temperature distributions, partial melting, and water content within the crust and mantle in earthquake seismology. Full waveform inversion (FWI), as one of the state‐of‐the‐art seismic imaging techniques, can produce high‐resolution constraints for subsurface (an)elastic parameters by minimizing the difference between observed and predicted seismograms. Traditional waveform inversion for attenuation is commonly based on the standard‐linear‐solid (SLS) wave equation, in which case the quality factor (Q) has to be converted to stress and strain relaxation times. When using multiple attenuation mechanisms in the SLS method, it is difficult to directly estimate these relaxation time parameters. Based on a time domain complex‐valued viscoacoustic wave equation, we present an FWI framework for simultaneously estimating subsurfacePwave velocity and attenuation distributions. BecauseQis explicitly incorporated into the viscoacoustic wave equation, we directly derivePwave velocity andQsensitivity kernels using the adjoint‐state method and simultaneously estimate their subsurface distributions. By analyzing the Gauss‐Newton Hessian, we observe strong interparameter crosstalk, especially the leakage from velocity toQ. We approximate the Hessian inverse using a preconditioned L‐BFGS method in viscoacoustic FWI, which enables us to successfully reduce interparameter crosstalk and produce accurate velocity and attenuation models. Numerical examples demonstrate the feasibility and robustness of the proposed method for simultaneously mapping complex velocity andQdistributions in the subsurface.

    more » « less
  2. Abstract

    The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

    more » « less
  3. Abstract

    The magnitude‐frequency distribution (MFD) describes the relative proportion of earthquake magnitudes and provides vital information for seismic hazard assessment. Theb‐value, derived from the MFD, is commonly used to estimate the probability that a future earthquake will exceed a specified magnitude threshold. Improved MFD andb‐value estimates are of great importance in the central and eastern United States where high volumes of fluid injection have contributed to a significant rise in seismicity over the last decade. In this study, we recalculate the magnitudes of 8,775 events for the 2011 Prague, Oklahoma sequence using a relative magnitude approach that depends only on waveform data to calculate magnitudes. We also compare the distribution of successive magnitude differences to the MFD and show that a combination of the magnitude difference distribution (MDFD) and relative magnitudes yields a reliable estimate ofb‐value. Using the MDFD and relative magnitudes, we examine the temporal and spatial variations in theb‐value and show thatb‐value ranges between ∼0.6 and 0.85 during the aftershock sequence for at least 5 months after theM5.7 mainshock, though areas surrounding the northeast part of the sequence experience higherb‐values (0.7–0.85) than the southwestern part of the Meeker‐Prague fault whereb‐value is the lowest (0.6–0.7). We also identify a cluster of off‐fault events with the highestb‐values in the catalog (0.85). These new estimates of MFD andb‐value will contribute to understanding of the relations between induced and tectonic earthquake sequences and promote discussion regarding the use ofb‐value in induced seismic hazard estimation.

    more » « less

    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
  5. Abstract

    Oklahoma is a multihazard environment where both natural (e.g., tornadoes) and technological hazards (e.g., induced seismicity) are significant, making Oklahoma a unique setting to better understand how to manage and prepare for multiple hazards. While studies have attempted to understand drivers of hazard adjustments, few have focused on the overall number of adjustments undertaken instead of individual adjustments or adjustments in a multihazard environment. To address these gaps, we employ a survey sample of 866 households in Oklahoma to understand households’ danger control responses (protective hazard adjustments) for tornado and earthquake risks in Oklahoma. We apply the extended parallel processing model (EPPM) to categorize respondents according to their relative level of perceived threat and efficacy of protective actions in predicting the number of hazard adjustments they intend to or have adopted in response to tornadoes and induced earthquakes. In line with the EPPM, we found that households have the highest number of danger control responses when their perceived threat and efficacy are both high. Counter to the EPPM literature, we found low threat coupled with high efficacy moved some individuals toward the adoption of danger control responses in response to both tornadoes and earthquakes. When households have high efficacy, threat appraisals matter in tornado danger control responses but not in earthquake danger control responses. This EPPM categorization opens new research approaches for studies of natural and technological hazards. This study also provides information for local officials and emergency managers making mitigation and preparedness investments and policies.

    more » « less