Search

Magnetotelluric (MT) data image the bulk resistivity of the subsurface which can be used to infer magma reservoir conditions beneath volcanoes. The bulk resistivity of magma depends primarily on the melt volume fraction, temperature, and water content. These variables are coupled thermodynamically, yet mixing relations for bulk resistivity implicitly treat them as independent. Here, we use a parameterization of the rhyolite‐MELTS thermodynamic model to constrain relationships between melt fraction, temperature, dissolved water content and bulk resistivity for rhyolitic magmas. This method results in MT interpretations which are (a) thermodynamically consistent at near‐equilibrium conditions, (b) independent of temperature and water content estimates derived from erupted products, and (c) able to consider saturated melts containing a volatile (i.e., aqueous fluid) phase. The utility of the method is demonstrated with three case studies of silicic systems: Mono Basin, Newberry volcano and the Laguna del Maule Volcanic Field (LdMVF). The moderately conductive feature at Mono Basin can be explained by under‐saturated partial melt (6–15 vol%) at <775°C, indicating relatively stable magma storage conditions since the last eruption. However, a relatively resistive feature at Newberry Volcano requires lower temperatures (<750°C) than previous estimates, suggesting that the system has cooled since the last eruption. A conductive feature at the LdMVF cannot be explained by saturated or under‐saturated rhyolitic melt and requires additional conductive phases. These results demonstrate the potential of this new method to reduce uncertainty in MT interpretations and highlight the need for additional coupling strategies between petrology, geophysics, and thermo‐mechanical models to better understand magmatic systems. Plain Language Summary: Determining the magma storage conditions beneath a volcano can give insights into the volcanic hazards. The magnetotelluric (MT) method measures the electrical resistivity of the subsurface, which is valuable because molten rock (i.e., melt) is less resistive than solid rock (or crystals). A mixture of melt and crystals will have a particular resistivity value which means that MT data can infer the relative proportion of melt (i.e., melt fraction). Magmas with higher melt fractions are more mobile and more likely to erupt. In addition, hotter melt with more dissolved water is more conductive than cooler, drier melt, so the electrical resistivity can also be used to infer the temperature and water content. However, there are more unknown variables (temperature, melt fraction, and water content) than known variables (observed resistivity), so more constraints are needed. Thermodynamic models indicate that hotter, wetter conditions have higher melt fractions than cooler, drier conditions, so the unknown variables are not independent. We use thermodynamic models to provide additional constraints on the interpretation of magma reservoir conditions from bulk resistivity, which enables us to provide realistic bounds on the range of plausible conditions and better explain the electrical resistivity observed at three different volcanoes. Key Points: A parameterization of the thermodynamic rhyolite‐MELTS model is coupled to empirical relations for the bulk resistivity of rhyolitic magmaThe coupling can significantly alter previous interpretations of magma reservoir conditions from magnetotelluric dataAllows for a comparison between past magma reservoir conditions inferred from erupted products and current magma reservoir conditions [ABSTRACT FROM AUTHOR]

Krafla is an active volcanic field and a high-temperature geothermal system in northeast Iceland. As part of a program to produce more energy from higher temperature wells, the IDDP-1 well was drilled in 2009 to reach supercritical fluid conditions below the Krafla geothermal field. However, drilling ended prematurely when the well unexpectedly encountered rhyolite magma at a depth of 2.1 km. In this paper we re-examine the magnetotelluric (MT) data that were used to model the electrical resistivity structure at Krafla. We present a new 3-D resistivity model that differs from previous inversions due to (1) using the full impedance tensor data and (2) a finely discretized mesh with horizontal cell dimensions of 100 m by 100 m. We obtained similar resistivity models from using two different prior models: a uniform half-space, and a previously published 1-D resistivity model. Our model contains a near-surface resistive layer of unaltered basalt and a low resistivity layer of hydrothermal alteration (C1). A resistive region (R1) at 1 to 2 km depth corresponds to chlorite-epidote alteration minerals that are stable at temperatures of about 220 to 500 °C. A low resistivity feature (C2) coincides with the Hveragil fault system, a zone of increased permeability allowing interaction of aquifer fluids with magmatic fluids and gases. Our model contains a large, low resistivity zone (C3) below the northern half of the Krafla volcanic field that domes upward to a depth of about 1.6 km b.s.l. C3 is partially coincident with reported low S -wave velocity zones which could be due to partial melt or aqueous fluids. The low resistivity could also be attributed to dehydration and decomposition of chlorite and epidote that occurs above 500 °C. As opposed to previously published resistivity models, our resistivity model shows that IDDP-1 encountered rhyolite magma near the upper edge of C3, where it intersects C2. In order to assess the sensitivity of the MT data to melt at the bottom of IDDP-1, we added hypothetical magma bodies with resistivities of 0.1 to 30 Ωm to our resistivity model and compared the synthetic MT data to the original inversion response. We used two methods to compare the MT data fit: (1) the change in r.m.s. misfit and (2) an asymptotic p -value obtained from the Kolmogorov–Smirnov (K–S) statistical test on the two sets of data residuals. We determined that the MT data can only detect sills that are unrealistically large (2.25 km3) with very low resistivities (0.1 or 0.3 Ωm). Smaller magma bodies (0.125 and 1 km3) were not detected; thus the MT data are not sensitive to small rhyolite magma bodies near the bottom of IDDP-1. Our tests gave similar results when evaluating the changes in r.m.s. misfit and the K–S test p -values, but the K–S test is a more objective method than appraising a relative change in r.m.s. misfit. Our resistivity model and resolution tests are consistent with the idea of rhyolite melt forming by re-melting of hydrothermally altered basalt on the edges of a deeper magma body. [ABSTRACT FROM AUTHOR]

The Laguna del Maule (LdM) volcanic field comprises the greatest concentration of postglacial rhyolite in the Andes and includes the products of ~40 km3 of explosive and effusive eruptions. Recent observations at LdM by interferometric synthetic aperture radar and global navigation satellite system geodesy have revealed inflation at rates exceeding 20 cm/year since 2007, capturing an ongoing period of growth of a potentially large upper crustal magma reservoir. Moreover, magnetotelluric and gravity studies indicate the presence of fluids and/or partial melt in the upper crust near the center of inflation. Petrologic observations imply repeated, rapid extraction of rhyolitic melt from crystal mush stored at depths of 4–6 km during at least the past 26 ka. We utilize multiple types of surface‐wave observations to constrain the location and geometry of low‐velocity domains beneath LdM. We present a three‐dimensional shear‐wave velocity model that delineates a ~450‐km3 shallow magma reservoir ~2 to 8 km below surface with an average melt fraction of ~5%. Interpretation of the seismic tomography in light of existing gravity, magnetotelluric, and geodetic observations supports this model and reveals variations in melt content and a deeper magma system feeding the shallow reservoir in greater detail than any of the geophysical methods alone. Geophysical imaging of the LdM magma system today is consistent with the petrologic inferences of the reservoir structure and growth during the past 20–60 kyr. Taken together with the ongoing unrest, a future rhyolite eruption of at least the scale of those common during the Holocene is a reasonable possibility. Key Points: Results from a major deployment of seismic instruments at Laguna del Maule volcanic center are presentedWe use a novel application of multiple types of surface‐wave dispersion data to infer 3D shear velocity structureWe present a unified interpretation of results from seismology, geodesy, gravity, magnetotellurics, and petrochronology [ABSTRACT FROM AUTHOR]