Nonlinear Conjugate Gradients Algorithm For 2-D Magnetotelluric Inversion

We investigate a new algorithm for computing regularized solutions of the two-dimensional magnetotelluric inverse problem. The algorithm employs a nonlinear conjugate gradients (NLCG) scheme to minimize an objective function that penalizes data residuals and second spatial derivatives of resistivity. We compare this algorithm theoretically and numerically to two previous algorithms for constructing such 'minimum-structure' models: the Gauss-Newton method, which solves a sequence of linearized inverse problems and has been the standard approach to nonlinear inversion in geophysics, and an algorithm due to Mackie and Madden, which solves a sequence of linearized inverse problems incompletely using a (linear) conjugate gradients technique. Numerical experiments involving synthetic and field data indicate that the two algorithms based on conjugate gradients (NLCG and Mackie-Madden) are more efficient than the GaussNewton algorithm in terms of both computer memory requirements and CPU time needed to find accurate solutions to problems of realistic size. This owes largely to the fact that the conjugate gradients-based algorithms avoid two computationally intensive tasks that are performed at each step of a Gauss-Newton iteration: calculation of the full Jacobian matrix of the forward modeling operator, and complete solution of a linear system on the model space. The numerical tests also show that the Mackie-Madden algorithm reduces the objective function more quickly than our new NLCG algorithm in the early stages of minimization, but NLCG is more effective in the later computations. To help understand these results, we describe the Mackie-Madden and new NLCG algorithms in detail and couch each as a special case of a more general conjugate gradients scheme for nonlinear inversion

Three-dimensional magnetotelluric modeling and inversion with applications to the California Basin and Range Province

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, 1991.Includes bibliographical references.by Randall Lee Mackie.Ph.D

3-D resistivity forward modeling and inversion using conjugate gradients

We have developed rapid 3-D dc resistivity forward modeling and inversion algorithms that use conjugate gradient relaxation techniques. In the forward network modeling calculation, an incomplete Cholesky decomposition for preconditioning and sparse matrix routines combine to produce a fast and efficient algorithm (approximately 2 minutes CPU time on a Sun SPARC‐station 2 for 50 × 50 × 20 blocks). The side and bottom boundary conditions are scaled impedance conditions that take into account the local current flow at the boundaries as a result of any configuration of current sources. For the inversion, conjugate gradient relaxation is used to solve the maximum likelihood inverse equations. Since conjugate gradient techniques only require the results of the sensitivity matrix [tilde under A] or its transpose [tilde under A][superscript T] multiplying a vector, we are able to bypass the actual computation of the sensitivity matrix and the inversion of [tilde under A][superscript T] [tilde under A], thus greatly decreasing the time needed to do 3-D inversions. We demonstrate 3-D resistivity tomographic imaging using pole‐pole resistivity data collected during an experiment for a leakage monitoring system near evaporation ponds at the Mojave Generating Station in Laughlin, Nevada.United States. Environmental Protection Agency (grant #CR-821516

Electrical conductivity of the Pampean Shallow Subduction Region of Argentina near 33 S: evidence for a slab window

We present a three-dimensional (3-D) interpretation of 117 long period (20–4096 s) magnetotelluric (MT) sites between 31°S and 35°S in western Argentina. They cover the most horizontal part of the Pampean shallow angle subduction of the Nazca Plate and extend south into the more steeply dipping region. Sixty-two 3-D inversions using various smoothing parameters and data misfit goals were done with a nonlinear conjugate gradient (NLCG) algorithm. A dominant feature of the mantle structure east of the horizontal slab is a conductive plume rising from near the top of the mantle transition zone at 410 km to the probable base of the lithosphere at 100 km depth. The subducted slab is known to descend to 190 km just west of the plume, but the Wadati-Benioff zone cannot be traced deeper. If the slab is extrapolated downdip it slices through the plume at 250 km depth. Removal of portions of the plume or blocking vertical current flow at 250 km depth significantly changes the predicted responses. This argues that the plume is not an artifact and that it is continuous. The simplest explanation is that there is a “wedge”-shaped slab window that has torn laterally and opens down to the east with its apex at the plume location. Stress within the slab and seismic tomography support this shape. Its northern edge likely explains why there is no deep seismicity south of 29°S.Fil: Burd, Aurora I.. University of Washington; Estados UnidosFil: Booker, John R.. University of Washington; Estados UnidosFil: Mackie, Randall. Land General Geophysics; ItaliaFil: Pomposiello, Maria Cristina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Geocronología y Geología Isotopica. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Geocronología y Geología Isotopica; ArgentinaFil: Favetto, Alicia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Geocronología y Geología Isotopica. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Geocronología y Geología Isotopica; Argentin

Fe-oxide grain coatings support bacterial Fe-reducing metabolisms in 1.7−2.0 km-deep subsurface quartz arenite sandstone reservoirs of the Illinois Basin (USA)

The Cambrian-age Mt. Simon Sandstone, deeply buried within the Illinois Basin of the midcontinent of North America, contains quartz sand grains ubiquitously encrusted with iron-oxide cements and dissolved ferrous iron in pore-water. Although microbial iron reduction has previously been documented in the deep terrestrial subsurface, the potential for diagenetic mineral cementation to drive microbial activity has not been well studied. In this study, two subsurface formation water samples were collected at 1.72 and 2.02 km, respectively, from the Mt. Simon Sandstone in Decatur, Illinois. Low-diversity microbial communities were detected from both horizons and were dominated by Halanaerobiales of Phylum Firmicutes. Iron-reducing enrichment cultures fed with ferric citrate were successfully established using the formation water. Phylogenetic classification identified the enriched species to be related to Vulcanibacillus from the 1.72 km depth sample, while Orenia dominated the communities at 2.02 km of burial depth. Species-specific quantitative analyses of the enriched organisms in the microbial communities suggest that they are indigenous to the Mt. Simon Sandstone. Optimal iron reduction by the 1.72 km enrichment culture occurred at a temperature of 40oC (range 20 to 60oC) and a salinity of 25 parts per thousand (range 25-75 ppt). This culture also mediated fermentation and nitrate reduction. In contrast, the 2.02 km enrichment culture exclusively utilized hydrogen and pyruvate as the electron donors for iron reduction, tolerated a wider range of salinities (25-200 ppt), and exhibited only minimal nitrate- and sulfate-reduction. In addition, the 2.02 km depth community actively reduces the more crystalline ferric iron minerals goethite and hematite. The results suggest evolutionary adaptation of the autochthonous microbial communities to the Mt. Simon Sandstone and carries potentially important implications for future utilization of this reservoir for CO2 injection

Electrical structure beneath the northern MELT line on the East Pacific Rise at 15°45′S

Author Posting. © American Geophysical Union, 2006. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 33 (2006): L22301, doi:10.1029/2006GL027528.The electrical structure of the upper mantle beneath the East Pacific Rise (EPR) at 15°45′S is imaged by inverting seafloor magnetotelluric data obtained during the Mantle ELectromagnetic and Tomography (MELT) experiment. The electrical conductivity model shows no evidence for a conductive region immediately beneath the ridge, in contrast to the model previously obtained beneath the EPR at 17°S. This observation can be explained by differences in current melt production along the ridge, consistent with other observations. The mantle to the east of the ridge at 60 –100 km depth is anisotropic, with higher conductivity in the spreading direction compared to the along-strike direction, similar to the 17°S region. The high conductivity in the spreading direction can be explained by a hydrated mantle with strain-induced lattice preferred orientation of olivine or by partial melt preferentially connected in the spreading direction.This work was supported by NSF grant OCE0118254

Compositional controls on oceanic plates : geophysical evidence from the MELT area

Author Posting. © The Authors, 2005. This is the author's version of the work. It is posted here by permission of Nature Publishing Group for personal use, not for redistribution. The definitive version was published in Nature 437 (2005): 249-252, doi:10.1038/nature04014.Magnetotelluric (MT) and seismic data, collected during the MELT experiment at the Southern East Pacific Rise (SEPR) constrain the distribution of melt beneath this mid-ocean-ridge spreading center and also the evolution of the oceanic lithosphere during its early cooling history. In this paper, we focus on structure imaged at distances ~100 to 350 km east of the ridge crest, corresponding to seafloor ages of ~1.3 to 4.5 Ma, where the seismic and electrical conductivity structure is nearly constant, independent of age. Beginning at a depth of about 60 km, there is a large increase in electrical conductivity and a change from isotropic to transversely anisotropic electrical structure with higher conductivity in the direction of fast propagation for seismic waves. Because conductive cooling models predict structure that increases in depth with age, extending to about 30 km at 4.5 Ma, we infer that the structure of young oceanic plates is instead controlled by a decrease in water content above 60 km induced by the melting process beneath the spreading center.US participation in the MELT experiment and subsequent analysis was funded by NSF grants through the Marine Geology and Geophysics Program, Ocean Sciences Division

The impact of breathing motion versus heterogeneity effects in lung cancer treatment planning

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134883/1/mp3427.pd