skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Probabilistic inversion of seafloor compliance for oceanic crustal shear velocity structure using mixture density neural networks
Title: Probabilistic inversion of seafloor compliance for oceanic crustal shear velocity structure using mixture density neural networks
SUMMARY Measurements of various physical properties of oceanic sediment and crustal structures provide insight into a number of geological and geophysical processes. In particular, knowledge of the shear wave velocity (VS) structure of marine sediments and oceanic crust has wide ranging implications from geotechnical engineering projects to seismic mantle tomography studies. In this study, we propose a novel approach to nonlinearly invert compliance signals recorded by colocated ocean-bottom seismometers and high-sample-rate pressure gauges for shallow oceanic shear wave velocity structure. The inversion method is based on a type of machine learning neural network known as a mixture density neural network (MDN). We demonstrate the effectiveness of the MDN method on synthetic models with a fixed deployment depth of 2015 m and show that among 30 000 test models, the inverted shear wave velocity profiles achieve an average error of 0.025 km s−1. We then apply the method to observed data recorded by a broad-band ocean-bottom station in the Lau basin, for which a VS profile was estimated using Monte Carlo sampling methods. Using the mixture density network approach, we validate the method by showing that our VS profile is in excellent agreement with the previous result. Finally, we argue that the mixture density network approach to compliance inversion is advantageous over other compliance inversion methods because it is faster and allows for standardized measurements.  more » « less
Award ID(s):
1658214
PAR ID:
10333930
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
227
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
1879 to 1892
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Journal of Geophysical Research: Solid Earth)
    Abstract We present models of compressional and shear velocity structure of the oceanic sediments and upper crust surrounding the Hawaiian islands. The models were derived from analysis of seafloor compliance data and measurements of Ps converted phases originating at the sediment‐bedrock interface. These data were estimated from continuous broadband ocean bottom seismometer acceleration and pressure records collected during the Plume‐Lithosphere Undersea Mantle Experiment, an amphibious array of wideband and broadband instruments with an aperture of over 1,000 km. Our images result from a joint inversion of compliance and Ps delay data using a nonlinear inversion scheme whereby deviation from a priori constraints is minimized. In our final model, sediment thickness increases from 50 m at distal sites to over 1.5 km immediately adjacent to the islands. The sedimentary shear velocity profiles exhibit large regional variations. While sedimentary structure accounts for the majority of the compliance signal, we infer variations in shear velocity in the uppermost bedrock on the order of ±5%. We also require relatively high values of Poisson's ratio in the uppermost crust. Lower crustal velocities are generally seen to the north and west of the islands but do not appear well correlated with the Hawaiian Swell bathymetry. A region of strong low velocity anomalies to the northeast of Hawaii may be associated with the Molokai fracture zone. 
    more » « less
  2. ; ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract A shallow sub‐seafloor seismic model that includes well‐determined seismic velocities and clarifies sediment‐crust discontinuities is needed to characterize the physical properties of marine sediments and the oceanic crust and to serve as a reference for deeper seismic modeling endeavors. This study estimates the seismic structure of marine sediments and the shallow oceanic crust of the Alaska‐Aleutian subduction zone at the Alaska Peninsula, using data from the Alaska Amphibious Community Seismic Experiment (AACSE). We measure seafloor compliance and Ps converted wave delays from AACSE ocean‐bottom seismometers (OBS) and seafloor pressure data and interpret these measurements using a joint Bayesian Monte Carlo inversion to produce a sub‐seafloor S‐wave velocity model beneath each available OBS station. The sediment thickness across the array varies considerably, ranging from about 50 m to 2.80 km, with the thickest sediment located in the continental slope. Lithological composition plays an important role in shaping the seismic properties of seafloor sediment. Deep‐sea deposits on the incoming plate, which contain biogenic materials, tend to have reduced S‐wave velocities, contrasting with the clay‐rich sediments in the shallow continental shelf and continental slope. A difference in S‐wave velocities is observed for upper oceanic crust formed at fast‐rate (Shumagin) and intermediate‐rate (Semidi) spreading centers. The reduced S‐wave velocities in the Semidi crust may be caused by increased faulting and possible lithological variations, related to a previous period of intermediate‐rate spreading. 
    more » « less
  3. ; ; (, The Seismic Record)
    Abstract We construct a 3D shear velocity model of the Salt Lake Valley using Rayleigh waves excited by the 31 March 2020 Mw 6.5 central Idaho earthquake recorded on a 168-station temporary nodal geophone network and the 49-station permanent regional network. The temporary array—deployed in response to the March 18 Mw 5.7 Magna earthquake—serendipitously recorded clear surface waves between 10 and 20 s period from the Idaho event at ∼500 km epicentral distance, from which we measure both Rayleigh wave phase velocity and ellipticity (H/V ratio). In addition, we employ multicomponent earthquake coda cross correlation to extend the measurements down to 5 s period. Because Rayleigh wave ellipticity features outstanding shallow sensitivity, we invert for a 3D upper crust VS model of the Salt Lake Valley. Our model shows basin structure in general agreement with and complements the current Community Velocity Model, which is mostly constrained by borehole and gravity measurements. Our model thus provides critical information for future earthquake hazard assessment studies, which require detailed shallow velocity structure. 
    more » « less
  4. ; ; ; ; ; (, Geophysical Research Letters)
    Abstract Near‐surface seismic velocity structure plays a critical role in ground motion amplification during large earthquakes. In particular, the local Vp/Vs ratio strongly influences the amplitude of Rayleigh waves. Previous studies have separately imaged 3D seismic velocity and Vp/Vs ratio at seismogenic depth, but lack regional coverage and/or fail to constrain the shallowest structure. Here, we combine three datasets with complementary sensitivity in a Bayesian joint inversion for shallow crustal shear velocity and near‐surface Vp/Vs ratio across Southern California. Receiver functions–including with an apparent delayed initial peak in sedimentary basins, and long considered a nuisance in receiver function imaging studies–highly correlate with short‐period Rayleigh wave ellipticity measurements and require the inclusion of a Vp/Vs parameter. The updated model includes near‐surface low shear velocity more in line with geotechnical layer estimates, and generally lower than expected Vp/Vs outside the basins suggesting widespread shallow fracturing and/or groundwater undersaturation. 
    more » « less
  5. ; ; ; (, The Leading Edge)
    Elastic full-waveform inversion (EFWI) is a state-of-the-art seismic tomographic method. Recent advances in technology and instrumentation, combining crosstalk-free source-encoded FWI (SE-FWI) with multicomponent marine data acquisition using ocean-bottom nodes (OBNs), enable full-physics wave propagation and parameter inversion without the computational burden of traditional FWI. With OBN acquisition, P waves, S waves, and P-to-S conversions are recorded. It is not well understood to what extent adding horizontal components to SE-FWI improves the resolution of subsurface modeling. We assess their potential for the reconstruction of shear and compressional wave speeds (VPand VS) by using a synthetic data set modeled after a recently acquired OBN survey in the North Sea. We perform synthetic inversion tests to design suitable strategies that leverage the information recorded in the horizontal components of the data to improve the reconstructed model resolution laterally and in depth. We advocate for a hierarchical inversion approach to recover the elastic parameters. We exploit the P and P-to-S converted waves recorded on the horizontal components to robustly reconstruct both VPand VS. Adding horizontal components to the SE-FWI modeling workflow results in improved spatial resolution, enhanced depth coverage, and more accurate elastic wave speed estimates. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email