Regional 3-D heterogeneities by waveform inversion – application to the Atlantic area

The waveform inversion method described in Woodhouse & Dziewonski (1984) was modified to retrieve regional scale 3-D heterogeneities by using the minor arc part of seismograms. The lateral heterogeneities are expanded horizontally into blocks (10°× 10°) and radially into Legendre polynomials up to order 3 (0–670 km), and thus the results show much fine details of lateral variation than previous global scale studies. We assumed that the heterogeneities produce the perturbation of eigenfrequencies which are the minor arc average of local eigenfrequency shift. We applied the method to the upper mantle beneath the Atlantic Ocean and its environments. Care was taken about the weighting of the data set. We found that the fit of each seismogram became better when the weighting of each seismogram is proportional to the inverse of initial data residuals. Resolution is good in the triangular region surrounded by South America, Europe, and North America. Resolution is not good in the South Atlantic because of the poor path coverage. Depth resolution is not clear, because of the use of Legendre polynomials, though the results suggest a broad half-width of the order of 200 km or more. We found some similarities between previous global studies and our results. For example, low velocities beneath the East Pacific Rise, Chile Rise and Azores triple junction and a high velocity Canadian shield are obtained. However, there are also differences; the high-velocity zone beneath the Brazilian shield at shallow depth is not a prominent feature in this study. Instead, we found a somewhat unexpected feature near the Romanche and Vema fracture zones where shallow positive anomalies exist. Smoothed results calculated by the spherical harmonic expansion are also shown for the purpose of comparison with global studies.     

Three-dimensional magnetotelluric inversion using non-linear conjugate gradients

We have formulated a 3-D inverse solution for the magnetotelluric (MT) problem using the non-linear conjugate gradient method. Finite difference methods are used to compute predicted data efficiently and objective functional gradients. Only six forward modelling applications per frequency are typically required to produce the model update at each iteration. This efficiency is achieved by incorporating a simple line search procedure that calls for a sufficient reduction in the objective functional, instead of an exact determination of its minimum along a given descent direction. Additional efficiencies in the scheme are sought by incorporating preconditioning to accelerate solution convergence. Even with these efficiencies, the solution's realism and complexity are still limited by the speed and memory of serial processors. To overcome this barrier, the scheme has been implemented on a parallel computing platform where tens to thousands of processors operate on the problem simultaneously. The inversion scheme is tested by inverting data produced with a forward modelling code algorithmically different from that employed in the inversion algorithm. This check provides independent verification of the scheme since the two forward modelling algorithms are prone to different types of numerical error.     

3-D magnetotelluric inversion and model validation with gravity data for the investigation of flood basalts and associated volcanic rifted margins

20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.     

Two efficient algorithms for iterative linearized inversion of seismic waveform data

We present two equivalent algorithms for iterative linearized waveform inversion for 3-D Earth structure with respect to an arbitrary 3-D starting model; one is a matrix formulation, and the second is a wavefield formulation. Both algorithms require the computation of accurate synthetic seismograms, but neither requires that any particular method be used to compute the synthetics. The matrix formulation is equivalent to our previously published algorithm (Hara, Tsuboi & Geller 1991), but requires less than 10 per cent of the CPU time of the previous algorithm. The wavefield algorithm is equivalent to that of Tarantola (1986) and Mora (1987), but appears to be substantially more efficient.