Imaging global mantle discontinuities: a test using full-waveforms and adjoint kernels

Summary We present a novel approach for imaging global mantle discontinuities based on full-waveform inversion (FWI). Over the past decades, extensive research has been done on imaging mantle discontinuities at approximately 400 km and 670 km depth. Accurate knowledge of their topography can put strong constraints on thermal and compositional variations and hence geodynamic modelling. So far, however, there is little consensus on their topography. We present an approach based on adjoint tomography, which has the advantage that Fréchet derivatives for discontinuities and measurements, to be inverted for, are fully consistent. Rather than working with real data, we focus on synthetic tests, where the answer is known in order to be able to evaluate the performance of the developed method. All calculations are based on the community code SPECFEM3D_GLOBE. We generate data in fixed 1-D or 3-D elastic background models of mantle velocity. Our ‘data’ to be inverted contain topography along the 400 km and 670 km mantle discontinuities. To investigate the approach, we perform several tests: (i) In a situation where we know the elastic background model 1-D or 3-D, we recover the target topography fast and accurately, (ii) The exact misfit is not of great importance here, except in terms of convergence speed, similar to a different inverse algorithm, (iii) In a situation where the background model is not known, the convergence is markedly slower, but there is reasonable convergence towards the correct target model of discontinuity topography. It has to be noted that our synthetic test is idealised and in a real data situation, the convergence to and uncertainty of the inferred model is bound to be larger. However, the use of data consistent with Fréchet kernels seems to pay off and might improve our consensus on the nature of mantle discontinuities. Our workflow could be incorporated in future FWI mantle models to adequately infer boundary interface topography.

Estimating bend-faulting and mantle hydration at the Marianas trench from seismic anisotropy

<p><span>Bend faults formed in oceanic lithosphere approaching deep ocean trenches promote water circulation and the formation of hydrous minerals. As the plate subducts, these minerals can dehydrate into the mantle wedge, generating the melts that feed arc volcanoes, or subduct fully into the deeper mantle. Balancing the global water budget requires an estimate of the amount of water recycled to the mantle by subduction, but current estimates for water fluxes at subduction zones span several orders of magnitude, mainly because of large uncertainties in the amount of water carried in the lithospheric mantle of the incoming plate. </span></p><p><span>We use active source seismic refraction data collected on the incoming plate at the Marianas trench to measure azimuthal seismic anisotropy in the uppermost mantle, and assess the degree of faulting and associated serpentinization of the uppermost mantle based on spatial variations in the observed anisotropy. We find that the fast direction of anisotropy varies with distance from the trench, rotating from APM-parallel at the eastern side of the study area to approximately fault-parallel near the trench. The fast direction orientations suggest that a coherent set of bend-faults are beginning to form at least 200 km out from the trench, although the extrinsic anisotropy signal from the faults does not substantially overprint the signal from preexisting mineral fabrics until the plate is ~100 km from the trench. The average (isotropic) mantle velocity decreases slightly as the plate nears the trench. Preliminary interpretation suggests that the observed spatial variations in anisotropy can be explained by serpentinization localized along pervasive, trench-parallel faults or joints.</span></p>

Uppermost Mantle Velocity beneath the Mid-Atlantic Ridge and Transform Faults in the Equatorial Atlantic Ocean

ABSTRACT Seismic rays traveling just below the Moho provide insights into the thermal and compositional properties of the upper mantle and can be detected as Pn phases from regional earthquakes. Such phases are routinely identified in the continents, but in the oceans, detection of Pn phases is limited by a lack of long-term instrument deployments. We present estimates of upper-mantle velocity in the equatorial Atlantic Ocean from Pn arrivals beneath, and flanking, the Mid-Atlantic Ridge and across several transform faults. We analyzed waveforms from 50 earthquakes with magnitude Mw>3.5, recorded over 12 months in 2012–2013 by five autonomous hydrophones and a broadband seismograph located on the St. Peter and St. Paul archipelago. The resulting catalog of 152 ray paths allows us to resolve spatial variations in upper-mantle velocities, which are consistent with estimates from nearby wide-angle seismic experiments. We find relatively high velocities near the St. Paul transform system (∼8.4 km s−1), compared with lower ridge-parallel velocities (∼7.7 km s−1). Hence, this method is able to resolve ridge-transform scale velocity variations. Ray paths in the lithosphere younger than 10 Ma have mean velocities of 7.9±0.5 km s−1, which is slightly lower than those sampled in the lithosphere older than 20 Ma (8.1 km±0.3 s−1). There is no apparent systematic relationship between velocity and ray azimuth, which could be due to a thickened lithosphere or complex mantle upwelling, although uncertainties in our velocity estimates may obscure such patterns. We also do not find any correlation between Pn velocity and shear-wave speeds from the global SL2013sv model at depths <150 km. Our results demonstrate that data from long-term deployments of autonomous hydrophones can be used to obtain rare and insightful estimates of uppermost mantle velocities over hundreds of kilometers in otherwise inaccessible parts of the deep oceans.

2D crust-upper mantle velocity structure along a seismic section in Nanling, South China

<p>    There exists an important polymetallic ore belt in Nanling of the southeastern China. Previous studies suggest that the mineralization of Nanling is probably related to the bottom intrusion of magmatic rocks in the late Mesozoic. In this study, a natural seismic section was installed by using 81 portable stations with an interval of 5 km from July 2017 to August 2019, which runs across the Nanling belt in the south of Fujian and Jiangxi provinces. As a result, we have picked up 3,818 relative residual data from 215 teleseismic events with magnitude greater than 5.5. And we have applied the teleseismic full-waveform tomography and the teleseismic travel-time tomography to study the crust and the mantle velocity structure beneath the Nanling metallogenic belt, respectively. Our preliminary results show that: (1) a clear low-velocity anomaly exists in the crust beneath the Zhenghe-Dapu fault and its east side, which might be related to the rich ore deposits in Nanling; (2) some high-velocity anomalies in the uppermost mantle beneath the Wuyi metallogenic belt may be relevant to the igneous rock cooling and the lithospheric thickening; (3) there are obvious low-velocity anomalies at the upper mantle beneath the Wuyi and Nanling metallogenic belts, which are speculated to be hot materials from asthenosphere upwelling into the bottom of the lithosphere. Our results provide a new insight for investigating the deep structures and deep dynamic processes of Nanling tectonic belt.</p>