 Award ID(s):
 2012046
 Publication Date:
 NSFPAR ID:
 10278867
 Journal Name:
 GEOPHYSICS
 Volume:
 86
 Issue:
 4
 Page Range or eLocationID:
 T179 to T192
 ISSN:
 00168033
 Sponsoring Org:
 National Science Foundation
More Like this

We have developed a Liouville partialdifferentialequation (PDE)based method for computing complexvalued eikonals in real phase space in the multivalued sense in attenuating media with frequencyindependent qualify factors, where the new method computes the real and imaginary parts of the complexvalued eikonal in two steps by solving Liouville equations in real phase space. Because the earth is composed of attenuating materials, seismic waves usually attenuate so that seismic data processing calls for properly treating the resulting energy losses and phase distortions of wave propagation. In the regime of highfrequency asymptotics, the complexvalued eikonal is one essential ingredient for describing wave propagationmore »

Firstarrival traveltime tomography is an essential method for obtaining nearsurface velocity models. The adjointstate firstarrival traveltime tomography is appealing due to its straightforward implementation, low computational cost, and low memory consumption. Because solving the pointsource isotropic eikonal equation by either ray tracers or eikonal solvers intrinsically corresponds to emanating discrete rays from the source point, the resulting traveltime gradient is singular at the source point, and we denote such a singular pattern the imprint of rayillumination. Because the adjointstate equation propagates traveltime residuals back to the source point according to the negative traveltime gradient, the resulting adjoint state will inheritmore »

SUMMARY Accurate synthetic seismic wavefields can now be computed in 3D 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 highamplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography.
This can bemore »
Here we present the results of proof of concept testing of such an approach for a synthetic 3component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3D 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 timeseries, starting from a 1D model, using source stacked waveforms and stationpair crosscorrelations 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 crosstalk 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 tradeoffs 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 singleevent wavefield computations.

Embedding properties of network realizations of dissipative reduced order models Jörn Zimmerling, Mikhail Zaslavsky,Rob Remis, Shasri Moskow, Alexander Mamonov, Murthy Guddati, Vladimir Druskin, and Liliana Borcea Mathematical Sciences Department, Worcester Polytechnic Institute https://www.wpi.edu/people/vdruskin Abstract Realizations of reduced order models of passive SISO or MIMO LTI problems can be transformed to tridiagonal and blocktridiagonal forms, respectively, via dierent modications of the Lanczos algorithm. Generally, such realizations can be interpreted as ladder resistorcapacitorinductor (RCL) networks. They gave rise to network syntheses in the rst half of the 20th century that was at the base of modern electronics design and consecutively to MORmore »

SUMMARY Knowledge of attenuation structure is important for understanding subsurface material properties. We have developed a doubledifference seismic attenuation (DDQ) tomography method for highresolution imaging of 3D attenuation structure. Our method includes two main elements, the inversion of eventpair differential ${t^*}$ ($d{t^*}$) data and 3D attenuation tomography with the $d{t^*}$ data. We developed a new spectral ratio method that jointly inverts spectral ratio data from pairs of events observed at a common set of stations to determine the $d{t^*}$ data. The spectral ratio method cancels out instrument and site response terms, resulting in more accurate $d{t^*}$ data compared to absolutemore »