Entity

Time filter

Source Type

Wellington, New Zealand

Kordy M.,University of Utah | Wannamaker P.,University of Utah | Maris V.,University of Utah | Cherkaev E.,University of Utah | And 2 more authors.
Geophysical Journal International | Year: 2016

Following the creation described in Part I of a deformable edge finite-element simulator for 3-D magnetotelluric (MT) responses using direct solvers, in Part II we develop an algorithm named HexMT for 3-D regularized inversion of MT data including topography. Direct solvers parallelized on large-RAM, symmetric multiprocessor (SMP) workstations are used also for the Gauss-Newton model update. By exploiting the data-space approach, the computational cost of the model update becomes much less in both time and computer memory than the cost of the forward simulation. In order to regularize using the second norm of the gradient, we factor thematrix related to the regularization termand apply its inverse to the Jacobian, which is done using the MKL PARDISO library. For dense matrix multiplication and factorization related to the model update, we use the PLASMA library which shows very good scalability across processor cores. A synthetic test inversion using a simple hill model shows that including topography can be important; in this case depression of the electric field by the hill can cause false conductors at depth ormask the presence of resistive structure.With a simplemodel of two buried bricks, a uniform spatial weighting for the norm of model smoothing recovered more accurate locations for the tomographic images compared to weightings which were a function of parameter Jacobians.We implement joint inversion for static distortionmatrices tested using the Dublin secret model 2, for which we are able to reduce nRMS to ~1.1 while avoiding oscillatory convergence. Finally we test the code on field data by inverting full impedance and tipper MT responses collected around Mount St Helens in the Cascade volcanic chain. Among several prominent structures, the north-south trending, eruption-controlling shear zone is clearly imaged in the inversion. © The Authors 2015. Source


Hill G.J.,Antarctica Scientific Ltd | Hill G.J.,University of Canterbury | Bibby H.M.,Institute of Geological & Nuclear Sciences | Ogawa Y.,Tokyo Institute of Technology | And 5 more authors.
Earth and Planetary Science Letters | Year: 2015

The dynamics of magma reservoirs (the main repositories for eruptible magma) play a fundamental role in the style and behaviour of volcanic systems. A key first step in understanding these systems is to identify their location and size accurately. We present results from a broadband magnetotelluric study of the Tongariro Volcanic system and discuss how the results fit within current petrological models. The Tongariro Volcanic system is a composite andesitic cone complex, located at the southern end of the Taupo Volcanic Zone in the central North Island of New Zealand. We use data from 136 broadband magnetotelluric soundings within a 25. ×. 35 km area covering the volcanic system to construct a 3D image of the magmatic system of the Tongariro Volcanic Complex including Mount Ngauruhoe. The structure of the Tongariro magmatic system has been determined from 3D forward and inverse modelling of the magnetotelluric data and allowed for an estimation of the melt fraction present within the system. 3D inverse modelling of the magnetotelluric data shows: a well-developed shallow low resistivity zone outlining the geothermal system; a zone of even lower resistivity representing a shallow crustal magma accumulation zone located at a depth of ~4-12 km offset to the east of the Tongariro vent system; and a zone with a slightly higher resistivity connecting these two components of the magmatic system providing the path for magmatic fluids from the deeper source region to reach the surface during eruptive events. © 2015. Source


Kordy M.,University of Utah | Wannamaker P.,University of Utah | Maris V.,University of Utah | Cherkaev E.,University of Utah | And 2 more authors.
Geophysical Journal International | Year: 2016

We have developed an algorithm, which we call HexMT, for 3-D simulation and inversion of magnetotelluric (MT) responses using deformable hexahedral finite elements that permit incorporation of topography. Direct solvers parallelized on symmetric multiprocessor (SMP), single-chassis workstations with large RAM are used throughout, including the forward solution, parameter Jacobians and model parameter update. In Part I, the forward simulator and Jacobian calculations are presented. We use first-order edge elements to represent the secondary electric field (E), yielding accuracy O(h) for E and its curl (magnetic field). For very low frequencies or small material admittivities, the E-field requires divergence correction. With the help of Hodge decomposition, the correction may be applied in one step after the forward solution is calculated. This allows accurate E-field solutions in dielectric air. The system matrix factorization and source vector solutions are computed using the MKL PARDISO library, which shows good scalability through 24 processor cores. The factorized matrix is used to calculate the forward response as well as the Jacobians of electromagnetic (EM) field and MT responses using the reciprocity theorem. Comparison with other codes demonstrates accuracy of our forward calculations.We consider a popular conductive/resistive double brick structure, several synthetic topographic models and the natural topography of Mount Erebus in Antarctica. In particular, the ability of finite elements to represent smooth topographic slopes permits accurate simulation of refraction of EM waves normal to the slopes at high frequencies. Run-time tests of the parallelized algorithm indicate that for meshes as large as 176 × 176 × 70 elements, MT forward responses and Jacobians can be calculated in ~1.5 hr per frequency. Together with an efficient inversion parameter step described in Part II, MT inversion problems of 200-300 stations are computable with total run times of several days on such workstations. © The Authors 2015. Source

Discover hidden collaborations