Kirchhoff synthetics for seismic reflections and diving waves in two dimensions and application to the D" layer

The Kirchhoff-Helmholtz (KH) integration has been used to model the reflected and the diving waves from an interface with a positive velocity gradient. The modelling is carried out for a spherical boundary and for a sinusoidal topography with a long-scale wavelength.
An artefact, which is a major problem in modelling the seismic response using the KH integration, has been reduced by introducing a Hilbert transform sign manipulation. Cleaner synthetic seismograms with correct amplitudes have been produced by this method. A discretization in larger surface elements has been made possible by introducing a smoothing factor that suppresses the noise that normally follows the constructed signal if a large element size is taken.     

Modelling of scattering of seismic waves from a corrugated rough sea surface: a comparison of three methods

We compare three numerical methods to model the sea surface interaction in a marine seismic reflection experiment (the frequencies considered are in the band 10–100 Hz): the finite-difference method (FDM), the spectral element method (SEM) and the Kirchhoff method (KM). A plane wave is incident at angles of 0° and 30° with respect to the vertical on a rough Pierson–Moskowitz surface with 2 m significant wave height and the response is synthesized at 6, 10 and 50 m below the average height of the sea surface. All three methods display an excellent agreement for the main reflected arrival. The FDM and SEM also agree very well all through the scattered coda. The KM shows some discrepancies, particularly in terms of amplitudes.     

Seismic imaging of the laterally varying D" region beneath the Cocos Plate

We use an axisymmetric, spherical Earth finite difference algorithm to model SH -wave propagation through cross-sections of laterally varying lower mantle models beneath the Cocos Plate derived from recent data analyses. Synthetic seismograms with dominant periods as short as 4 s are computed for several models: (1) a D" reflector 264 km above the core–mantle boundary with laterally varying S -wave velocity increases of 0.9–2.6 per cent, based on localized structures from a 1-D double-array stacking method; (2) an undulating D" reflector with large topography and uniform velocity increase obtained using a 3-D migration method and (3) cross-sections through the 3-D mantle S -wave velocity tomography model TXBW. We apply double-array stacking to assess model predictions of data. Of the models explored, the S -wave tomography model TXBW displays the best overall agreement with data. The undulating reflector produces a double Scd arrival that may be useful in future studies for distinguishing between D" volumetric heterogeneity and D" discontinuity topography. Synthetics for the laterally varying models show waveform variability not observed in 1-D model predictions. It is challenging to predict 3-D structure based on localized 1-D models when lateral structural variations are on the order of a few wavelengths of the energy used, particularly for the grazing geometry of our data. Iterative approaches of computing synthetic seismograms and adjusting model characteristics by considering path integral effects are necessary to accurately model fine-scale D" structure.     

The accuracy of models derived by WKBJ waveform inversion

Summary . In this paper the accuracy of velocity-depth profiles derived by matching WKBJ seismograms to observations is quantitatively evaluated. Seismograms computed with the WKBJ method are generally quite reliable but possess predictable, systematic inaccuracies in the presence of strong velocity gradients. The effects of these inaccuracies on models derived through WKBJ waveform inversion are studied, using reflectivity seismograms as 'data'. The velocity structure used is an oceanic lithosphere model that contains several transition regions separated by relatively homogeneous layers, producing partially-reflected reverberations in the reflectivity synthetics that are absent from the WKBJ seismograms. The inversion incorporates the 'jumping' strategy to solve for the smoothest models consistent with the data. We find these solutions to be independent of the starting model and to have a stable basic structure that agrees well with the correct model. The differences, everywhere less than a seismic wavelength, depend on the frequency content of the seismograms. Reverberations in the reflectivity seismograms that are well separated from WKBJ arrivals are treated as 'noise' in the inversion.     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.