An official website of the United States government
SUMMARY The uneven distribution of earthquakes and stations in seismic tomography leads to slower convergence of nonlinear inversions and spatial bias in inversion results. Including dense regional arrays, such as USArray or Hi-Net, in global tomography causes severe convergence and spatial bias problems, against which conventional pre-conditioning schemes are ineffective. To save computational cost and reduce model bias, we propose a new strategy based on a geographical weighting of sources and receivers. Unlike approaches based on ray density or the Voronoi tessellation, this method scales to large full-waveform inversion problems and avoids instabilities at the edges of dense receiver or source clusters. We validate our strategy using a 2-D global waveform inversion test and show that the new weighting scheme leads to a nearly twofold reduction in model error and much faster convergence relative to a conventionally pre-conditioned inversion. We implement this geographical weighting strategy for global adjoint tomography.
more »
« less
Accelerating full waveform inversion via source stacking and cross-correlations
SUMMARY Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography. This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility of designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography. Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function. The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.
more »
« less
- Award ID(s):
- 1758198
- PAR ID:
- 10123100
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 220
- Issue:
- 1
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. 308-322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The cross‐correlation of a diffuse or random wavefield at two points has been demonstrated to recover an empirical estimate of the Green's function under a wide variety of source conditions. Over the past two decades, the practical development of this principle, termed ambient noise interferometry, has revolutionized the fields of seismology and acoustics. Yet, because of the spatial sparsity of conventional water column and seafloor instrumentation, such array‐based processing approaches have not been widely utilized in oceanography. Ocean‐bottom distributed acoustic sensing (OBDAS) repurposes pre‐existing optical fibers laid in seafloor cables as dense arrays of broadband strain sensors, which observe both seismic waves and ocean waves. The thousands of sensors in an OBDAS array make ambient noise interferometry of ocean waves straightforward for the first time. Here, we demonstrate the application of ambient noise interferometry to surface gravity waves observed on an OBDAS array near the Strait of Gibraltar. We focus particularly on a 3‐km segment of the array on the continental shelf, containing 300 channels at 10‐m spacing. By cross‐correlating the raw strain records, we compute empirical ocean surface gravity wave Green's functions for each pair of stations. We first apply beamforming to measure the time‐averaged dispersion relation along the cable. Then, we exploit the non‐reciprocity of waves propagating in a flow to recover the depth‐averaged current velocity as a function of time using a waveform stretching method. The result is a spatially continuous matrix of current velocity measurements with resolution <100 m and <1 hr.more » « less
-
SUMMARY Ultra-low velocity zones (ULVZs) have been identified as regions of extremely low velocity anomalies in the Earth's lowermost mantle using seismic observations from reflected, refracted and diffracted arrivals along the mantle side of the core–mantle boundary. Estimation of ULVZ geometrical (i.e. shape and size) and elastic (i.e. velocity and density) parameters with uncertainties is crucial in understanding the role of ULVZs in the ongoing dynamic processes within the Earth's mantle; however, these parameters are still poorly known due to uncertainties and tradeoffs of the seismically resolved ULVZ geometries and elastic parameters. Computation of synthetic waveforms for 2-D and 3-D ULVZs shapes is currently computationally feasible, but past studies utilize higher dimensional waveform modelling of mostly only low-frequency diffracted waves. Most studies focusing on high-frequency core-reflected waveforms (e.g. ScP) still use 1-D modelling approaches to determine ULVZ properties. This approach might lead to wrong results if the imaged structures have inherently 3-D geometries. This study investigates high-frequency synthetic ScP waveforms for various 2.5-D ULVZ geometries showing that additional seismic arrivals are generated even when the ScP geometrical ray path does not directly strike the location of the ULVZ. The largest amplitude additional phases in the 2.5-D models are post-cursor arrivals that are generated at the edges of the finite-length ULVZs. These newly identified ScP post-cursors can arrive within the ScsP post-cursor time window traditionally analysed in 1-D ULVZ studies. These post-cursors might then be misidentified or constructively/destructively interfere with the ScsP post-cursor, leading to incorrect estimation of ULVZ parameters. In this study we investigate the bias introduced by the 2.5-D morphologies on the 1-D estimated ULVZ elastic properties in a Bayesian waveform inversion scheme. We further expand the Bayesian method by including the data noise covariance matrix in the inversion and compare it to an autoregressive noise model that was utilized in previous studies. From the application to the observed ScP data, we find that the new approach converges faster, particularly for the inversion of data from multiple events, and the new algorithm retrieves ULVZ parameters with more realistic uncertainties. The inversion of 2.5-D synthetic ScP waveforms suggests that the retrieved ULVZ parameters can be misleading with unrealistically high confidence if we do not consider the data noise covariance matrix in the inversion. Our new approach can also retrieve the shape of a multidimensional Gaussian ULVZ if its length is 12° or longer in the great circle arc direction. However, 2.5-D synthetic waveforms show additional waveform complexity which can constructively interfere with the ScP wavefield. Hence, in many cases the estimation of ULVZ properties using 1-D forward modelling can provide incorrect ULVZ parameters. Hence previous ULVZ modelling efforts using 1-D parameter estimation methods may be incorrect.more » « less
-
SUMMARY We present a new 3-D radially anisotropic seismic velocity model EARA2024 of the crust and mantle beneath East Asia and the northwestern Pacific using adjoint full-waveform inversion tomography. We construct the EARA2024 model by iteratively minimizing the waveform similarity misfit between the synthetic and observed waveforms from 142 earthquakes recorded by about 2000 broad-band stations in East Asia. Compared to previous studies, this new model renders significantly improved images of the subducted oceanic plate in the upper mantle, mantle transition zone, and uppermost lower mantle along the Kuril, Japan, Izu-Bonin and Ryukyu Trenches. Complex slab deformation and break-offs are observed at different depths. Moreover, our model provides new insights into the origins of intraplate volcanoes in East Asia, including the Changbaishan, Datong-Fengzhen, Tengchong and Hainan volcanic fields.more » « less
-
SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.more » « less
