skip to main content


Title: Joint inversion for 1-D crustal seismic S - and P -wave velocity structures with interfaces and its application to the Wabash Valley Seismic Zone
SUMMARY Interfaces are important part of Earth’s layering structure. Here, we developed a new model parametrization and iterative linearized inversion method that determines 1-D crustal velocity structure using surface wave dispersion, teleseismic P-wave receiver functions and Ps and PmP traveltimes. Unlike previous joint inversion methods, the new model parametrization includes interface depths and layer Vp/Vs ratios so that smoothness constraint can be conveniently applied to velocities of individual layers without affecting the velocity discontinuity across the interfaces. It also allows adding interface-related observation such as traveltimes of Ps and PmP in the joint inversion to eliminate the trade-off between interface depth and Vp/Vs ratio and therefore to reduce the uncertainties of results. Numerical tests show that the method is computationally efficient and the inversion results are robust and independent of the initial model. Application of the method to a dense linear array across the Wabash Valley Seismic Zone (WVSZ) produced a high-resolution crustal image in this seismically active region. The results show a 51–55-km-thick crust with a mid-crustal interface at 14–17 km. The crustal Vp/Vs ratio varies from 1.69 to 1.90. There are three pillow-like, ∼100 km apart high-velocity bodies sitting at the base of the crust and directly above each of them are a low-velocity anomaly in the middle crust and a high-velocity anomaly in the upper crust. They are interpreted to be produced by mantle magmatic intrusions and remelting during rifting events in the end of the Precambrian. The current diffuse seismicity in the WVSZ might be rooted in this ancient distributed rifting structure.  more » « less
Award ID(s):
1661519
NSF-PAR ID:
10296010
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
226
Issue:
1
ISSN:
0956-540X
Page Range / eLocation ID:
47 to 55
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present two new seismic velocity models for Alaska from joint inversions of body-wave and ambient-noise-derived surface-wave data, using two different methods. Our work takes advantage of data from many recent temporary seismic networks, including the Incorporated Research Institutions for Seismology Alaska Transportable Array, Southern Alaska Lithosphere and Mantle Observation Network, and onshore stations of the Alaska Amphibious Community Seismic Experiment. The first model primarily covers south-central Alaska and uses body-wave arrival times with Rayleigh-wave group-velocity maps accounting for their period-dependent lateral sensitivity. The second model results from direct inversion of body-wave arrival times and surface-wave phase travel times, and covers the entire state of Alaska. The two models provide 3D compressional- (VP) and shear-wave velocity (VS) information at depths ∼0–100  km. There are many similarities as well as differences between the two models. The first model provides a clear image of the high-velocity subducting plate and the low-velocity mantle wedge, in terms of the seismic velocities and the VP/VS ratio. The statewide model provides clearer images of many features such as sedimentary basins, a high-velocity anomaly in the mantle wedge under the Denali volcanic gap, low VP in the lower crust under Brooks Range, and low velocities at the eastern edge of Yakutat terrane under the Wrangell volcanic field. From simultaneously relocated earthquakes, we also find that the depth to the subducting Pacific plate beneath southern Alaska appears to be deeper than previous models. 
    more » « less
  2. 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
  3. We constrain the viscosity of the lower crust through a joint inversion of seismic P-wave (Vp) and S-wave (Vs) velocities. Previous research has demonstrated robust relationships between seismic velocity and crustal composition, as well as between the composition and viscosity of the lower crust. Here we extend these analyses, showing seismic velocity can be used as a robust indicator of crustal viscosity. First, we use the Gibbs free energy minimization routine Perple_X to calculate equilibrium mineral assemblages for a global compilation of crustal rocks at various pressures and temperatures. Second, we use a rheological mixing model that combines single-phase flow laws for major crust-forming minerals to calculate bulk viscosity from the predicted mineral assemblages incorporating the effects of strain rate, temperature, pressure, and water activity. We apply our method to regional seismic and heat flow data across East Asia in order to separate the relative variations in mid-crustal viscosity associated with composition and temperature. In some regions, temperature variations are the dominant influence on viscosity; e.g., we predict a 3 order of magnitude increase in viscosity between the low heat flow Sichuan Basin and higher heat flow surrounding regions. These viscosity variations are consistent with those previously inferred to produce the different topographic gradients in these areas [1]. However in constant heat flow regions, compositional variations exert the primary influence on viscosity; e.g., the North China Craton and the Yangtze Craton are predicted to have compositionally-controlled viscosities ranging from 1022–1023 Pa×s. Finally, the regional Vp/Vs ratios in the Tibetan Plateau cannot be explained by thermal and/or compositional variations alone, possibly indicating the presence of melt, which would lead to additional viscosity reductions. [1] Clark & Royden, Geology, 2000. 
    more » « less
  4. Abstract

    We invertPg,PmP, andPntraveltimes from an active‐source, multiscale tomography experiment to constrain the three‐dimensional isotropic and anisotropicPwave velocity structure of the topmost oceanic mantle and crust and crustal thickness variations beneath the entire Endeavour segment of the Juan de Fuca Ridge. The isotropic velocity structure is characterized by a semicontinuous, narrow (5‐km‐wide) crustal low‐velocity volume that tracks the sinuous ridge axis. Across the Moho, the low‐velocity volume abruptly broadens to approximately 20 km in width and displays a north‐south linear trend that connects the two overlapping spreading centers bounding the segment. From the seismic results, we estimate the thermal structure and melt distribution beneath the Endeavour segment. The thermal structure indicates that the observed skew, or lateral offset, between the crustal and mantle magmatic systems is a consequence of differences in mechanisms of heat transfer at crustal and mantle depths, with the crust and mantle dominated by advection and conduction, respectively. Melt volume estimates exhibit significant along‐axis variations that coincide with the observed skew between the mantle and crustal magmatic systems, with sites of enhanced crustal melt volumes and vigorous hydrothermal activity corresponding to regions where the mantle and crustal magmatic systems are vertically aligned. These results contradict models of ridge segmentation that predict enhanced and reduced melt supply beneath the segment center and ends, respectively. Our results instead support a model in which segment‐scale skew between the crustal and mantle magmatic systems governs magmatic and hydrothermal processes at mid‐ocean ridges.

     
    more » « less
  5. Abstract

    The Hawaiian‐Emperor seamount chain in the Pacific Ocean has provided fundamental insights into hotspot generated intraplate volcanism and the long‐term strength of oceanic lithosphere. However, only a few seismic experiments to determine crustal and upper mantle structure have been carried out on the Hawaiian Ridge, and no deep imaging has ever been carried out along the Emperor seamounts. Here, we present the results of an active source seismic experiment using 29 Ocean‐Bottom Seismometers (OBS) carried out along a strike profile of the seamounts in the region of Jimmu and Suiko guyots. Joint reflection and refraction tomographic inversion of the OBS data show the upper crust is highly heterogeneous withPwave velocities <4–5 km s−1, which are attributed to extrusive lavas and clastics. In contrast, the lower crust is remarkably homogeneous with velocities of 6.5–7.2 km s−1, which we attribute to oceanic crust and mafic intrusions. Moho is identified by a strongPmParrival at offsets of 20–80 km, yielding depths of 13–16 km. The underlying mantle is generally homogeneous with velocities in the range 7.9–8.0 km s−1. The crust and mantle velocity structure has been verified by gravity modeling. While top of oceanic crust prior to volcano loading is not recognized as a seismic or gravity discontinuity, flexural modeling reveals a ∼5.0–5.5 km thick preexisting oceanic crust that is overlain by a ∼8 km thick volcanic edifice. Unlike at the Hawaiian Ridge, we find no evidence of magmatic underplating.

     
    more » « less