Both crustal velocity and density models provide important constraints for understanding crustal tectonics, composition and magmatic system. However, conventional gravity inversion methods for density models are troubled by poor depth-resolution and non-uniqueness. We apply a joint inversion method to full-waveform ambient noise and gravity data for simultaneously deriving high-resolution 3-D crustal S-wave velocity (Vs) and density models. By constraints of seismic and gravity data, our joint inversion could significantly improve the resolution of density models and reduce the uncertainty on the inversion results. Our method is suitable for any areas especially with low seismicity and can be extended to basin-scale. We apply the method on the seismic and gravity data around the Jingpohu (JPH) volcanic area in Northeast China to obtain high-resolution 3-D crustal Vs and density models, which are subsequently used to build the crustal lithological model. Our models demonstrate that the JPH volcanic group is located near the junction of the Mudanjiang–Yilan and Dunhua–Mishan fault zones and was tectonically controlled by them. A steep crustal-scale mafic intrusion with high Vs and density is present beneath the JPH volcanic group, and we interpret it as the same product as the JPH volcanic basaltic rocks outcropped on the surface by the mantle-sourced magmatic upwelling. Furthermore, none significant molten magma reservoirs with low Vs and density are observed currently within the crust of the JPH volcanic group. Thus, our results support that the JPH volcanic group is a single volcanic genesis with a mantle-only magma source.
more » « less- PAR ID:
- 10563457
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 240
- Issue:
- 2
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 1168-1180
- Size(s):
- p. 1168-1180
- Sponsoring Org:
- National Science Foundation
More Like this
-
null (Ed.)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
-
Santorini volcano in the South Aegean Volcanic Arc has a detailed history of ongoing volcanic and seismic activity, making it a prime location for studying magma storage and transport at arc volcanoes. The shallow magmatic system (<5 km depth) is well constrained by geophysical studies, but the deeper crustal structure is not. Located 15 km NE of Santorini, the Kolumbo seamount is also an active edifice, with consistently more seismicity and hydrothermal venting than Santorini. Geochemical studies indicate that Santorini and Kolumbo are fed by separate mantle and crustal magma sources, but prior seismic studies suggest otherwise (Dimitriadis et al, 2010; McVey et al, 2020). This study addresses the nature of lower-crustal magma structure beneath arc volcanoes and whether evolved volcanoes and nearby vents are connected through their plumbing. Tomographic inversion of P-wave Moho reflection (PmP) and turning P-wave (Pg) traveltimes is used to create 3-D models of Moho depth and P-wave velocity (Vp) down to depths of ~25 km. The PROTEUS experiment provides an exceptionally dense and large aperture traveltime dataset from an amphibious array of ~150 seismometers and ~14,000 active marine sources. The data are ~33,000 manually picked PmP arrivals and ~256,000 Pg arrivals from existing studies. Results show a low Vp anomaly extending from the Moho to the surface. This anomaly starts at the base of the crust under the NW Santorini caldera and extends up to the east. It is most pronounced at 10-15 km depth, where it is offset from both Santorini and Kolumbo. Limited resolution prevents imaging of a connection between this mid-crustal anomaly and the known shallow magma storage region under the Santorini caldera. A high-velocity core beneath Santorini is not found, a feature interpreted at other volcanoes as a cooled intrusive complex. Because no additional low Vp anomalies are found in the lower crust, we infer that a common mantle source and mid-crustal plumbing system is actively feeding both Santorini and Kolumbo. The spatial offset and elongated nature of magma storage implies a complex relationship between evolving magmatic structures and tectonics.more » « less
-
Abstract Seismic tomography of the crust is an essential tool for studying the three‐dimensional structure of magmatic plumbing systems feeding active volcanoes, but it is often limited in resolution by the absence of deep local seismicity. Teleseismic receiver functions can be used to illuminate local structural variations, but typically do not account for the effects of three‐dimensional velocity heterogeneities. Here we harness the complementary strengths of both techniques by processing
Ps‐P delay times derived from teleseismic receiver functions in a tomographicS wave inversion. Using our inversion technique, we produce the first tomographic crustal velocity model beneath Cleveland Volcano, identifying a vertically extensive highV P /V S anomaly beneath the volcano that likely signifies a middle‐to‐lower crustal magma reservoir. The observation is the first of its kind in the central Aleutians, illustrating the potential of our technique to advance our understanding of crustal magmatic systems without broad seismic networks or distributed local seismicity. -
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
-
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