P-wave traveltime and polarization tomography of VSP data

The presence of anisotropy requires that tomographic methods be generalized to account for anisotropy. This generalization allows geological structure to be correctly imaged and allows the anisotropic parameters to be estimated. Use of isotropic inversion for imaging anisotropic structures gives systematic trends in the traveltime and polarization residuals. However, due to the limited directional coverage, the traveltimes along may not be sufficient to study the anisotropic properties of the structure. Polarizations can provide independent information on the structure. Traveltime and polarization inversion are applied to synthetic examples simulating VSP experiments. Transverse isotropy and 1-D structure are assumed. Plots of traveltime and polarization residuals are an important tool to detect the anomalies due to the presence of anisotropy. For receivers located in anisotropic layers, polarization residuals display consistent anomalies of several degrees. The synthetic examples show that even the simple 1-D problem is difficult, when using direct arrivals only. Large a posteriori errors in anisotropic parameters are obtained by traveltime inversion in layers where available incidence angles are less than 45°. Resolution of the tomographic image of VSP data is greatly improved by a combination of traveltime and polarization information. In order to obtain accurate inversion results, the measurement error of polarization data should be kept to within a few degrees.     

Palaeoseismicity of the Oquirrh fault, Utah from shallow seismic tomography

The magnitude and frequency of normal-fault palaeoearthquakes are usually determined by trenching studies that ascertain the size and number of colluvial wedges along the fault. Such information can be invaluable in predicting the seismic hazard and potential for a future earthquake in that region. Digging trenches across normal faults, however, is environmentally intrusive, expensive and limited in the penetration depth. To overcome these problems we propose the use of 3-D seismic tomography as a means to identify the shapes and sizes of colluvial wedges along normal faults. As an example,2-D and 3-D seismic surveys were conducted across the Oquirrh fault, Utah with the purpose of imaging the normal-fault structure to a depth of about 10  m. Results show that the 3-D tomogram clearly delineates the fault zone and a colluvial wedge, both of which correlate extremely well with the geological cross-section interpreted from an adjacent trench. The thickness of the colluvial wedge image is used in conjunction with a seismic section to compute an estimate of a 6.8 moment magnitude earthquake for the most recent event on this fault, which is in close agreement with the 7.0 estimate based on a nearby trenching study. These tomographic results demonstrate, for the first time, that seismic imaging methods can be used in some cases to estimate unambiguously the shapes of colluvial wedges and the sizes of prehistoric earthquakes. Thus, seismic tomography has the possibility of providing cheaper, deeper and wider, but less resolved, images of fault systems than the intrusive excavation of trenches across faults.     

Simultaneous velocity and interface tomography of normal-incidence and wide-aperture seismic traveltime data

A tomographic inversion technique that inverts traveltimes to obtain a model of the subsurface in terms of velocities and interfaces is presented. It uses a combination of refraction, wide-angle reflection and normal-incidence data, it simultaneously inverts for velocities and interface depths, and it is able to quantify the errors and trade-offs in the final model. The technique uses an iterative linearized approach to the non-linear traveltime inversion problem. The subsurface is represented as a set of layers separated by interfaces, across which the velocity may be discontinuous. Within each layer the velocity varies in two dimensions and has a continuous first derivative. Rays are traced in this medium using a technique based on ray perturbation theory, and two-point ray tracing is avoided by interpolating the traveltimes to the receivers from a roughly equidistant fan of rays. The calculated traveltimes are inverted by simultaneously minimizing the misfit between the data and calculated traveltimes, and the roughness of the model. This 'smoothing regularization' stabilizes the solution of the inverse problem. In practice, the first iterations are performed with a high level of smoothing. As the inversion proceeds, the level of smoothing is gradually reduced until the traveltime residual is at the estimated level of noise in the data. At this point, a minimum-feature solution is obtained, which should contain only those features discernible over the noise.
The technique is tested on a synthetic data set, demonstrating its accuracy and stability and also illustrating the desirability of including a large number of different ray types in an inversion.     

Crustal and uppermost mantle structure in southern Africa revealed from ambient noise and teleseismic tomography

Rayleigh wave phase velocity maps in southern Africa are obtained at periods from 6 to 40 s using seismic ambient noise tomography applied to data from the Southern Africa Seismic Experiment (SASE) deployed between 1997 and 1999. These phase velocity maps are combined with those from 45 to 143 s period which were determined previously using a two-plane-wave method by Li & Burke. In the period range of overlap (25–40 s), the ambient noise and two-plane-wave methods yield similar phase velocity maps. Dispersion curves from 6 to 143 s period were used to estimate the 3-D shear wave structure of the crust and uppermost mantle on an 1°× 1° grid beneath southern Africa to a depth of about 100 km. Average shear wave velocity in the crust is found to vary from 3.6 km s^–1 at 0–10 km depths to 3.86 km s^–1 from 20 to 40 km, and velocity anomalies in these layers correlate with known tectonic features. Shear wave velocity in the lower crust is on average low in the Kaapvaal and Zimbabwe cratons and higher in the surrounding Proterozoic terranes, such as the Limpopo and the Namaqua-Natal belts, which suggests that the lower crust underlying the Archean cratons is probably less mafic than beneath the Proterozoic terranes. Crustal thickness estimates agree well with a previous receiver function study of Nair et al. . Archean crust is relatively thin and light and underlain by a fast uppermost mantle, whereas the Proterozoic crust is thick and dense with a slower underlying mantle. These observations are consistent with the southern African Archean cratons having been formed by the accretion of island arcs with the convective removal of the dense lower crust, if the foundering process became less vigorous in arc environments during the Proterozoic.     

Shallow velocity structure along the Hirapur–Mandla profile using traveltime inversion of wide-angle seismic data, and its tectonic implications

In order to investigate the velocity structure, and hence shed light on the related tectonics, across the Narmada–Son lineament, traveltimes of wide-angle seismic data along the 240 km long Hirapur–Mandla profile in central India have been inverted. A blocky, laterally heterogeneous, three-layer velocity model down to a depth of 10 km has been derived. The first layer shows a maximum thickness of the upper Vindhyans (4.5 km s⁻¹ ) of about 1.35 km and rests on top of normal crystalline basement, represented by the 5.9 km s⁻¹ velocity layer. The anomalous feature of the study is the absence of normal granitic basement in the great Vindhyan Graben, where lower Vindhyan sediments (5.3 km s⁻¹ ) were deposited during the Precambrian on high-velocity (6.3 km s⁻¹ ) metamorphic rock. The block beneath the Narmada–Son lineament represents a horst feature in which high-velocity (6.5 km s⁻¹ ) lower crustal material has risen to a depth of less than 2 km. South of the lineament, the Deccan Traps were deposited on normal basement during the upper Cretaceous period and attained a maximum thickness of about 800 m.