No abstract

S U M M A R YWe determine finite-frequency sensitivity kernels for seismic interferometry based upon noise cross-correlation measurements. Under the assumptions that noise is spatially uncorrelated but non-uniform, we determine ensemble-averaged cross correlations between synthetic seismograms at two geographically distinct locations. By minimizing a measure of the difference between observed and simulated ensemble cross correlations-subject to the constraint that the simulated wavefield satisfies the seismic wave equation-we obtain ensemble sensitivity kernels. These ensemble kernels reflect the sensitivity of ensemble cross-correlation measurements to variations in model parameters, for example, mass density, shear and compressional wave speeds and the spatial distribution of noise. Ensemble kernels are calculated based upon the interaction between two wavefields: an ensemble forward wavefield and an ensemble adjoint wavefield. To obtain the ensemble forward wavefield, one first calculates a generating wavefield obtained by inserting a signal determined by the characteristics of the noise at the location of the first receiver, saving the results of this calculation at locations where noise is generated, that is, typically on (a portion of) the Earth's surface. Next, one uses this generating wavefield as the source of the ensemble forward wavefield associated with the first receiver. The ensemble adjoint wavefield is obtained by using measurements between simulated and observed ensemble cross correlations as a seismic source located at the second receiver. The interaction between ensemble forward and adjoint wavefields 'paints' ensemble sensitivity kernels. We illustrate the construction of ensemble kernels and their nature in two and three dimensions using a spectral-element method. In addition to a 'banana-doughnut' feature connecting the two receivers, as in traditional finite-frequency earthquake tomography, some noise cross-correlation sensitivity kernels exhibit hyperbolic 'jets' protruding from each receiver in a direction away from the other receiver. Ensemble sensitivity kernels for long-period (T > ∼50 s) non-uniform noise in global models exhibit sensitivity along the minor and major arcs. These kernels reflect the fact that measurements typically involve long time-series that include multi-orbit surface waves. Like free oscillations, such measurements are sensitive to structure along the great circle through the two receivers. From the perspective of noise cross-correlation tomography, we discuss strategies for inversions in terrestrial and helioseismology.

International audienceWe present forward and adjoint spectral-element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite-frequency Fr'echet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite-frequency kernels generated by Version 2.0 'Sesame' of our widely used open source spectral-element package SPECFEM3D

We present the first-generation global tomographic model constructed based on adjoint tomography, an iterative full-waveform inversion technique. Synthetic seismograms were calculated using GPU-accelerated spectral-element simulations of global seismic wave propagation, accommodating effects due to 3-D anelastic crust & mantle structure, topography & bathymetry, the ocean load, ellipticity, rotation, and self-gravitation. Fréchet derivatives were calculated in 3-D anelastic models based on an adjoint-state method. The simulations were performed on the Cray XK7 named 'Titan', a computer with 18 688 GPU accelerators housed at Oak Ridge National Laboratory. The transversely isotropic global model is the result of 15 tomographic iterations, which systematically reduced differences between observed and simulated three-component seismograms. Our starting model combined 3-D mantle model S362ANI with 3-D crustal model Crust2.0. We simultaneously inverted for structure in the crust and mantle, thereby eliminating the need for widely used 'crustal corrections'. We used data from 253 earthquakes in the magnitude range 5.8 ≤ M w ≤ 7.0. We started inversions by combining ∼30 s body-wave data with ∼60 s surface-wave data. The shortest period of the surface waves was gradually decreased, and in the last three iterations we combined ∼17 s body waves with ∼45 s surface waves. We started using 180 min long seismograms after the 12th iteration and assimilated minor-and major-arc body and surface waves. The 15th iteration model features enhancements of well-known slabs, an enhanced image of the Samoa/Tahiti plume, as well as various other plumes and hotspots, such as Caroline, Galapagos, Yellowstone and Erebus. Furthermore, we see clear improvements in slab resolution along the Hellenic and Japan Arcs, as well as subduction along the East of Scotia Plate, which does not exist in the starting model. Point-spread function tests demonstrate that we are approaching the resolution of continentalscale studies in some areas, for example, underneath Yellowstone. This is a consequence of our multiscale smoothing strategy in which we define our smoothing operator as a function of the approximate Hessian kernel, thereby smoothing gradients less wherever we have good ray coverage, such as underneath North America.

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.