An official website of the United States government
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
Constructing a 3‐D Radially Anisotropic Crustal Velocity Model for Oklahoma Using Full Waveform Inversion
Abstract 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):
- 2042098
- PAR ID:
- 10475991
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 128
- Issue:
- 11
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract The interpretation of precursory seismicity can depend on a critical nucleation length scale h*, yet h* is largely unconstrained in the seismogenic crust. To estimate h* and associated earthquake nucleation processes at 2–7 km depths in Oklahoma, we studied seismic activity occurring prior to nine M 2.5–3.0 earthquakes that are aftershocks of the 3 September 2016 M 5.8 Pawnee, Oklahoma, earthquake. Four of the nine M 2.5–3.0 aftershocks studied did not have detectable seismicity within a 2 km radius of their hypocenters in the preceding 16 hr time windows. For the other five events, which did exhibit foreshock sequences, we estimated the static stress changes associated with each event of each sequence based on precise earthquake relocations and magnitude estimates. By carefully examining the spatiotemporal characteristics, we found all five of these M 2.5–3.0 aftershocks, and 70% of our studied events were plausibly triggered via static stress transfer from nearby earthquakes occurring hours to seconds earlier, consistent with the cascade nucleation model and a small h* in this region. The smallest earthquakes we could quantitatively study were M −1.5 events, which likely have 1–2 m rupture dimensions. The existence of these small events also supports a small nucleation length scale h*≤1 m, consistent with laboratory estimates. However, our observations cannot rule out more complicated earthquake initiation processes involving interactions between foreshocks and slow slip. Questions also remain as to whether aftershocks initiate differently from more isolated earthquakes.more » « less
-
Abstract The lithospheric structure of the contiguous US and surrounding regions offers clues into the tectonic history, including interactions between subducting slabs and cratons. In this paper, we present a new radially anisotropic shear wave speed model of the upper mantle (70–410 km) of the contiguous US and surrounding regions, constrained by seismic full‐waveform inversion. The new model (named CUSRA2021) utilizes frequency‐dependent travel time measurements, from 160 earthquake events recorded by 5,280 stations. The data coverage in eastern US is improved by incorporating more intraplate earthquakes. The final model exhibits clear and detailed shear wave speed anomalies correlating well with tectonic units such as North America Craton (high‐Vs), Cascadia subduction zones (high‐Vs), Columbia Plateau (low‐Vs), Basin and Range (low‐Vs), etc. In particular, the detailed structure of the North America Craton beneath Illinois basin is revealed. The depth of high‐Vs anomaly beneath the North America Craton correlates well with S‐to‐P receiver function and SH reflection results. Besides, the radial anisotropy in the Craton lithosphere shows a layering structure, which may relate to the process of lithospheric accretion and the origin of mid‐lithosphere discontinuities.more » « less
-
J.P. Hambleton, R. Makhnenko (Ed.)Deep wastewater injection-induced seismicity has led to over a thousand Mw>3 earthquakes and four Mw>5 earthquakes in Oklahoma over the last ten years. The 3 September 2016, Mw5.8 Pawnee, Oklahoma, earthquake was the first induced seismic event worldwide, that the authors are aware of, where liquefaction was observed and documented, raising concerns regarding the liquefaction risk posed by future induced earthquakes. Determining the suitability of current variants of the simplified procedure for evaluating the regional liquefaction hazard caused by induced earthquakes is part of an ongoing study. A detailed site characterization campaign examining profiles where liquefaction was and was not observed following the 2016 Mw5.8 Pawnee, Oklahoma, earthquake is part of this study. The purpose of this paper is to present an overview of the sites targeted as part of this testing, a summary of preliminary results from the site characterization campaign, and a description of planned future testing.more » « less
