The inversion of electromagnetic sounding data does not yield a unique solution, but inevitably a single model to interpret the observations is sought. We recommend that this model be as simple, or smooth, as possible, in order to reduce the temptation to overinterpret the data and to eliminate arbitrary discontinuities in simple layered models.To obtain smooth models, the nonlinear forward problem is linearized about a starting model in the usual way, but it is then solved explicitly for the desired model rather than for a model correction. By parameterizing the model in terms of its first or second derivative with depth, the minimum norm solution yields the smoothest possible model.Rather than fitting the experimental data as well as possible (which maximizes the roughness of the model), the smoothest model which fits the data to within an expected tolerance is sought. A practical scheme is developed which optimizes the step size at each iteration and retains the computational efficiency of layered models, resulting in a stable and rapidly convergent algorithm. The inversion of both magnetotelluric and Schlumberger sounding field data, and a joint magnetotelluric-resistivity inversion, demonstrate the method and show it to have practical application.

Magnetotelluric (MT) data are inverted for smooth 2-D models using an extension of the existing 1-D algorithm, Occam’s inversion. Since an MT data set consists of a finite number of imprecise data, an infinity of solutions to the inverse problem exists. Fitting field or synthetic electromagnetic data as closely as possible results in theoretical models with a maximum amount of roughness, or structure. However, by relaxing the misfit criterion only a small amount, models which are maximally smooth may be generated. Smooth models are less likely to result in overinterpretation of the data and reflect the true resolving power of the MT method. The models are composed of a large number of rectangular prisms, each having a constant conductivity. [Formula: see text] information, in the form of boundary locations only or both boundary locations and conductivity, may be included, providing a powerful tool for improving the resolving power of the data. Joint inversion of TE and TM synthetic data generated from known models allows comparison of smooth models with the true structure. In most cases, smoothed versions of the true structure may be recovered in 12–16 iterations. However, resistive features with a size comparable to depth of burial are poorly resolved. Real MT data present problems of non‐Gaussian data errors, the breakdown of the two‐dimensionality assumption and the large number of data in broadband soundings; nevertheless, real data can be inverted using the algorithm.

The lithosphere-asthenosphere boundary (LAB) separates rigid oceanic plates from the underlying warm ductile asthenosphere. Although a viscosity decrease beneath this boundary is essential for plate tectonics, a consensus on its origin remains elusive. Seismic studies identify a prominent velocity discontinuity at depths thought to coincide with the LAB but disagree on its cause, generally invoking either partial melting or a mantle dehydration boundary as explanations. Here we use sea-floor magnetotelluric data to image the electrical conductivity of the LAB beneath the edge of the Cocos plate at the Middle America trench offshore of Nicaragua. Underneath the resistive oceanic lithosphere, the magnetotelluric data reveal a high-conductivity layer confined to depths of 45 to 70 kilometres. Because partial melts are stable at these depths in a warm damp mantle, we interpret the conductor to be a partially molten layer capped by an impermeable frozen lid that is the base of the lithosphere. A conductivity anisotropy parallel to plate motion indicates that this melt has been sheared into flow-aligned tube-like structures. We infer that the LAB beneath young plates consists of a thin, partially molten, channel of low viscosity that acts to decouple the overlying brittle lithosphere from the deeper convecting mantle. Because this boundary layer has the potential to behave as a lubricant to plate motion, its proximity to the trench may have implications for subduction dynamics.

Marine controlled-source electromagnetic ͑CSEM͒ surveying has been in commercial use for predrill reservoir appraisal and hydrocarbon exploration for 10 years. Although a recent decrease has occurred in the number of surveys and publications associated with this technique, the method has become firmly established as an important geophysical tool in the offshore environment. This is a consequence of two important aspects associated with the physics of the method: First, it is sensitive to high electrical resistivity, which, although not an unambiguous indicator of hydrocarbons, is an important property of economically viable reservoirs. Second, although the method lacks the resolution of seismic wave propagation, it has a much better intrinsic resolution than potential-field methods such as gravity and magnetic surveying, which until now have been the primary nonseismic data sets used in offshore exploration. Although by many measures marine CSEM is still in its infancy, the reliability and noise floors of the instrument systems have improved significantly over the last decade, and interpretation methodology has progressed from simple anomaly detection to 3D anisotropic inversion of multicomponent data using some of the world's fastest supercomputers. Research directions presently include tackling the airwave problem in shallow water by applying time-domain methodology, continuous profiling tools, and the use of CSEM for reservoir monitoring during production.