To: “Linear inversion of body wave data — Part I: Velocity structure from traveltimes and ranges,” LeRoy M. Dorman and R. S. Jacobson (GEOPHYSICS, v. 46, p. 138–151).

Geophysics ◽  
1981 ◽  
Vol 46 (6) ◽  
pp. 958-958

The authors of “Linear inversion of body wave data — Part I: Velocity structure from traveltimes and ranges,” LeRoy M. Dorman and R. S. Jacobson, have forwarded a change in their paper which appeared in the February 1981 issue of (Geophysics, v. 46, p. 138–151).

2020 ◽  
Vol 222 (3) ◽  
pp. 1639-1655
Author(s):  
Xin Zhang ◽  
Corinna Roy ◽  
Andrew Curtis ◽  
Andy Nowacki ◽  
Brian Baptie

SUMMARY Seismic body wave traveltime tomography and surface wave dispersion tomography have been used widely to characterize earthquakes and to study the subsurface structure of the Earth. Since these types of problem are often significantly non-linear and have non-unique solutions, Markov chain Monte Carlo methods have been used to find probabilistic solutions. Body and surface wave data are usually inverted separately to produce independent velocity models. However, body wave tomography is generally sensitive to structure around the subvolume in which earthquakes occur and produces limited resolution in the shallower Earth, whereas surface wave tomography is often sensitive to shallower structure. To better estimate subsurface properties, we therefore jointly invert for the seismic velocity structure and earthquake locations using body and surface wave data simultaneously. We apply the new joint inversion method to a mining site in the United Kingdom at which induced seismicity occurred and was recorded on a small local network of stations, and where ambient noise recordings are available from the same stations. The ambient noise is processed to obtain inter-receiver surface wave dispersion measurements which are inverted jointly with body wave arrival times from local earthquakes. The results show that by using both types of data, the earthquake source parameters and the velocity structure can be better constrained than in independent inversions. To further understand and interpret the results, we conduct synthetic tests to compare the results from body wave inversion and joint inversion. The results show that trade-offs between source parameters and velocities appear to bias results if only body wave data are used, but this issue is largely resolved by using the joint inversion method. Thus the use of ambient seismic noise and our fully non-linear inversion provides a valuable, improved method to image the subsurface velocity and seismicity.


Geophysics ◽  
1981 ◽  
Vol 46 (2) ◽  
pp. 138-151 ◽  
Author(s):  
LeRoy M. Dorman ◽  
R. S. Jacobson

In a laterally homogeneous medium, the traveltime (T) and distance (X) for a ray with horizontal slowness p are linearly related to the depth Z(v) at which the velocity v = 1/p occurs. In order to exploit this linearity, we must infer the inverse velocity p from the observations of X, T pairs. Uncertainty in the determination of p causes correlation between the X and T observations. This correlation can be eliminated by rotation of the data into a coordinate system in which the covariance matrix is diagonal. These independent coordinates are, except for a scaling factor, the well‐known intercept time [Formula: see text] and a new variable [Formula: see text] The derivatives of T and X with respect to a depth‐velocity model contain singularities and so do those for ζ. These singularities can be quelled by representing the model as a stack of layers, each of which has a constant velocity gradient. Depth is then obtained by integration of the gradients. The sharpness of the partial derivatives of ζ w.r.t. the layer gradients indicates that ζ contains information in a more concentrated form than does τ. This manifests itself in smaller error bounds on the solution when ζ observations are used to supplement τ data. In the determination of ζ(p) from X,T data, an uncertainty principle or tradeoff applies. The delta‐like nature of the zeta partial derivatives means that the uncertainty in ζ will be closely related to the solution uncertainty and that we should choose in the parameterization the ζ, p pair which minimizes the uncertainty in ζ. This will avoid degrading the ultimate depth resolution achievable while still in the parameterization stage. We have applied these methods to sea floor hydrophone and surface buoy data from the Bengal Fan, and, we derive a model whose gradient is [Formula: see text] at the surface reaching [Formula: see text] at 500 m and remaining constant to at least 5.5 km.


1971 ◽  
Vol 61 (5) ◽  
pp. 1369-1379 ◽  
Author(s):  
Nezihi Canitez ◽  
M. Nafi Toksöz

abstract The determination of focal depth and other source parameters by the use of first-motion data and surface-wave spectra is investigated. It is shown that the spectral ratio of Love to Rayleigh waves (L/R) is sensitive to all source parameters. The azimuthal variation of the L/R spectral ratios can be used to check the fault-plane solution as well as for focal depth determinations. Medium response, attenuation, and source finiteness seriously affect the absolute spectra and introduce uncertainty into the focal depth determinations. These effects are nearly canceled out when L/R amplitude ratios are used. Thus, the preferred procedure for source mechanism studies of shallow earthquakes is to use jointly the body-wave data, absolute spectra of surface waves, and the Love/Rayleigh spectral ratios. With this procedure, focal depths can be determined to an accuracy of a few kilometers.


1969 ◽  
Vol 59 (5) ◽  
pp. 2071-2078
Author(s):  
Tom Landers ◽  
Jon F. Claerbout

abstract The inability of simple layered models to fit both Rayleigh wave and Love wave data has led to the proposal of an upper mantle interleaved with thin soft horizontal layers. Since surface-wave dispersion is not sensitive to the distribution of soft material but only to the fraction of soft material a variety of models is possible. The solution to this indeterminancy is found through body-wave analysis. It is shown that body waves are dispersed according to the thinness and softness of the layers. Three models, each of which satisfy all surface-wave data, are examined. Transmission seismograms calculated for these models show one to be impossible, one improbable and the other possible. Synthesis of the seismograms is accomplished through the use of time domain theory as the complicated frequency response of the models makes a frequency oriented Haskell-Thompson approach impractical.


2021 ◽  
Author(s):  
Nevra Bulut ◽  
Valerie Maupin ◽  
Hans Thybo

<p><span>We present a seismic tomographic image of Fennoscandia based on data from the ScanArray project in Norway, Sweden, and Finland, which operated during 2012-2017, together with data from earlier projects and stationary stations. We use relative traveltime residuals of P- and S- waves in high- and low-frequency bands and apply the frequency-dependent crustal correction. We use seismic signals from earthquakes at epicentral distances between 30° and 104° and magnitudes larger than 5.5. The general purpose of this study is to understand the possible causes of the high topography in Scandinavia along the passive continental margins in the North Atlantic as well as the interrelation between structure at the surface and in the lithospheric mantle.</span></p><p><span>We present an upper-mantle velocity structure for most Fennoscandia derived for the depth range 50-800 km with a 3D multiscale parameterization for an inversion mesh-grid with dimensions </span><em><span>dx</span></em><span>=</span><em><span>dy</span></em><span>=17.38 km and </span><em><span>dz</span></em><span>=23.44 km. In all body-wave tomography methods, smearing of anomalies is expected. Therefore resolution tests are critical for assessing the resolution of the parameters determined in the velocity models. The resolution of the models depends on several factors, including the noise level and general quality of data, the density of observations, the distance and back-azimuthal distribution of sources, the damping applied, and the model parameterization. We use checkerboard and model-driven (block and cylindrical) tests for assessing the resolution of our models.</span></p><p><span>Seismic models derived in this study are compared to existing and past topography to contribute to understanding mechanisms responsible for the topographic changes in the Fennoscandian region. The models also provide a basis for deriving high-resolution models of temperature and compositional anomalies that may contribute to understanding the observed, enigmatic topography.</span></p>


Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1917-1929 ◽  
Author(s):  
Joseph P. Stefani

Turning‐ray tomography is useful for estimating near‐surface velocity structure in areas where conventional refraction statics techniques fail because of poor data or lack of smooth refractor/velocity structure. This paper explores the accuracy and inherent smoothing of turning‐ray tomography in its capacity to estimate absolute near‐surface velocity and the statics times derived from these velocities, and the fidelity with which wavefields collapse to point diffractors when migrated through these estimated velocities. The method comprises nonlinear iterations of forward ray tracing through triangular cells linear in slowness squared, coupled with the LSQR linear inversion algorithm. It is applied to two synthetic finite‐ difference data sets of types that usually foil conventional refraction statics techniques. These models represent a complex hard‐rock overthrust structure with a low‐velocity zone and pinchouts, and a contemporaneous near‐shore marine trench filled with low‐ velocity unconsolidated deposits exhibiting no seismically apparent internal structure. In both cases velocities are estimated accurately to a depth of one‐ fifth the maximum offset, as are the associated statics times. Of equal importance, the velocities are sufficiently accurate to correctly focus synthetic wavefields back to their initial point sources, so migration/datuming applications can also use these velocities. The method is applied to a real data example from the Timbalier Trench in the Gulf of Mexico, which exhibits the same essential features as the marine trench synthetic model. The Timbalier velocity inversion is geologically reasonable and yields long and short wavelength statics that improve the CMP gathers and stack and that correctly align reflections to known well markers. Turning‐ray tomography estimates near‐surface velocities accurately enough for the three purposes of lithology interpretation, statics calculations, and wavefield focusing for shallow migration and datuming.


Geophysics ◽  
2008 ◽  
Vol 73 (4) ◽  
pp. D17-D33 ◽  
Author(s):  
Bing Zhou ◽  
Stewart Greenhalgh ◽  
Alan Green

Crosshole seismic tomography often is applied to image the velocity structure of an interwell medium. If the rocks are anisotropic, the tomographic technique must be adapted to the complex situation; otherwise, it leads to a false interpretation. We propose a nonlinear kinematic inversion method for crosshole seismic tomography in composite transversely isotropic media with known dipping symmetry axes. This method is based on a new version of the first-order traveltime perturbation equation. It directly uses the derivative of the phase velocity rather than the eigenvectors of the body-wave modes to overcome the singularity problem for application to the two quasi-shear waves. We applied an iterative nonlinear solver incorporating our kinematic ray-tracing scheme and directly compute the Jacobian matrix in an arbitrary reference medium. This reconstructs the five elastic moduli or Thomsen parameters from the first-arrival traveltimes of the three seismic body waves (qP, qSV, qSH) in strongly and weakly anisotropic media. We conducted three synthetic experiments that involve determining anisotropic parameters for a homogeneous rock, reconstructing a fault embedded in a strongly anisotropic background, and imaging a complicated four-layer model containing a small channel and a buried dipping interface. We compared results of our nonlinear inversion method with isotropic tomography and the traditional linear anisotropic inversion scheme, which showed the capability and superiority of the new scheme for crosshole tomographic imaging.


Sign in / Sign up

Export Citation Format

Share Document