skip to main content


Title: Tracking Crustal Permeability and Hydrothermal Response During Seafloor Eruptions at the East Pacific Rise, 9°50’N
Abstract

Permeability controls energy and matter fluxes in deep‐sea hydrothermal systems fueling a 'deep biosphere' of microorganisms. Here, we indirectly measure changes in sub‐seafloor crustal permeability, based on the tidal response of high‐temperature hydrothermal vents at the East Pacific Rise 9°50’N preceding the last phase of volcanic eruptions during 2005–2006. Ten months before the last phase of the eruptions, permeability decreased, first rapidly, and then steadily as the stress built up, until hydrothermal flow stopped altogether ∼2 weeks prior to the January 2006 eruption phase. This trend was interrupted by abrupt permeability increases, attributable to dike injection during last phase of the eruptions, which released crustal stress, allowing hydrothermal flow to resume. These observations and models suggest that abrupt changes in crustal permeability caused by magmatic intrusion and volcanic eruption can control first‐order hydrothermal circulation processes. This methodology has the potential to aid eruption forecasting along the global mid‐ocean ridge network.

 
more » « less
Award ID(s):
1948936 1949485
NSF-PAR ID:
10367745
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
3
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The loss of volatiles from a magma reservoir affects the magmatic overpressure responsible for driving ground deformation and eruptions. Although the high‐temperature metamorphic aureole around a magma chamber is typically considered to have low permeability, recent theoretical, experimental, and field studies have highlighted the role of transient permeability in magmatic systems. Also, direct measurements suggest that passive degassing is a significant component of total volatile loss in both basaltic and silicic volcanoes. Consequently, the effective permeability of the crust when magma is present in the system can be many orders of magnitude larger than that of exhumed rock samples. We develop a fully coupled porothermoelastic framework to account for both the flow of volatiles as well as associated effects on the stress state of the crust and calculate an analytical solution for spherical geometry. We then combine a magma chamber box model with these solutions to analyze eruption dynamics in magmatic systems. We find that in addition to viscous relaxation, magma recharge, and cooling timescales, the pore pressure diffusion timescale exerts a first‐order control on volcanic eruptions with moderately high crustal permeabilities of order 10−17to 10−19m2. We describe a parameter space to identify which components dominate in different regimes for volcanic eruptions according to these different timescales.

     
    more » « less
  2. Abstract

    Seafloor volcanic eruptions are difficult to directly observe due to lengthy eruption cycles and the remote location of mid‐ocean ridges. Volcanic eruptions in 2005–2006 at 9°50′N on the East Pacific Rise have been well documented, but the lava volume and flow extent remain uncertain because of the limited near‐bottom bathymetric data. We present near‐bottom data collected during 19 autonomous underwater vehicle (AUV)Sentrydives at 9°50′N in 2018, 2019, and 2021. The resulting 1 m‐resolution bathymetric grid and 20 cm‐resolution sidescan sonar images cover 115 km2, and span the entire area of the 2005–2006 eruptions, including an 8 km2pre‐eruption survey collected with AUVABEin 2001. Pre‐ and post‐eruption surveys, combined with sidescan sonar images and seismo‐acoustic impulsive events recorded during the eruptions, are used to quantify the lava flow extent and to estimate changes in seafloor depth caused by lava emplacement. During the 2005–2006 eruptions, lava flowed up to ∼3 km away from the axial summit trough, covering an area of ∼20.8 km2; ∼50% larger than previously thought. Where pre‐ and post‐eruption surveys overlap, individual flow lobes can be resolved, confirming that lava thickness varies from ∼1 to 10 m, and increases with distance from eruptive fissures. The resulting lava volume estimate indicates that ∼57% of the melt extracted from the axial melt lens probably remained in the subsurface as dikes. These observations provide insights into recharge cycles in the subsurface magma system, and are a baseline for studying future eruptions at the 9°50′N area.

     
    more » « less
  3. Abstract

    Conduit models of volcanic eruptions simulate magma evolution through phase transitions and material changes during ascent. We present a time‐dependent one‐dimensional model of a chamber‐conduit system to examine the temporal evolution of dome‐forming eruptions. As magma ascends, volatiles exsolve and may escape vertically through the column or laterally through the conduit walls. Magma solidifies which increases viscosity, leading to a natural transition from viscous flow at depth to frictional sliding along the conduit walls near the surface, resulting in the extrusion of a semisolid plug. The model evaluates time‐ and depth‐dependent pressure, velocity, porosity, and relative amounts of exsolved water to carbon dioxide. Transient effects arise when magma outflux from the chamber appreciably decreases pressure over the magma ascent timescale. For low magma permeability, transient effects increase porosity and velocity relative to steady‐state solutions. For high magma permeability, efficient vertical and lateral gas escape depresses porosity and velocity at later times. We use the model to predict three time series data sets from the 2004–2008 eruption of Mount St. Helens: extruded volume, ground deformation, and carbon dioxide emissions. We quantify sensitivity of model predictions to input parameters using the distance‐based generalized sensitivity analysis. Chamber volatile content, volume, and excess pressure influence the amplitude of observables, while conduit radius, frictional rate dependence and magma permeability influence temporal evolution. High magma permeability can cause marked departures from exponentially decaying flux and may explain the unique temporal evolution of deformation observed at the only nearby continuous GPS station in operation at the eruption onset.

     
    more » « less
  4. Abstract

    Although episodes of surface uplift and elevated seismicity precede many volcanic eruptions, their temporal evolution is often complex, and apparently in contradiction to simple trends predicted by mechanical deformation models. Here, we use continuous global positioning system and seismic data recorded at Sierra Negra volcano, Galápagos Islands, to show how the edifice responded to stress changes driven by magma accumulation in a shallow sill. The rate of uplift varied during the 13 years and 6.5 m of inflation before the 2018 eruption. The number of earthquakes per unit of uplift increased exponentially with total uplift as the differential stress increased. Accordingly, the temporal seismicity rate varied in time as a function of both the total uplift and the uplift rate. The Gutenberg‐Richterb‐value decreased as a function of total uplift. In the final six months before the eruption, a sequence of large (M > 4) earthquakes regulated the state of stress on the fault, each being followed by 2–3 days of postseismic quiescence, and retarding the increase in seismicity rate. These earthquakes did not affect the overall uplift rate. Subsidence of 8.5 m accompanied the 2‐month eruption. On resumption of uplift, the number of earthquakes per unit of uplift was very low, and theb‐value high, reflecting the relaxed stress state of the fault system. These observations show that crustal deformation becomes increasingly brittle at higher stress states, and supports theoretical models based on elastic‐brittle mechanics. They suggest that joint interpretation of deformation and seismicity is key for forecasting future eruptions in similar volcanic settings.

     
    more » « less
  5. Abstract

    We present deformation measurements of the Kirishima volcanic complex from ALOS and ALOS‐2 Interferometric Synthetic Aperture Radar (InSAR) time series during 2006–2019. Shinmoe‐dake deflated ∼6 cm during the 2008–2010 phreatic eruptions and inflated ∼5 cm prior to the 2017 magmatic eruption. Iwo‐yama inflated ∼19 cm within the crater since January 2015 and ∼7 cm around the southern and western vents since 4 months before the 2018 eruption. These deformations can be modeled as ellipsoids at ∼700 m depth beneath Shinmoe‐dake and as a sphere on top of an ellipsoid at ∼130 and ∼340 m depths beneath Iwo‐yama. Combining geodetic, geoelectric, geochemical, and petrological analysis, we interpret the hydrothermal origin of the deflation at Shinmoe‐dake and inflation at Iwo‐yama; the hydrothermal–magmatic transition during the 2011 Shinmoe‐dake eruption; and water‐boiling and bottom‐up pressurization as driving mechanisms of the inflation at Iwo‐yama. The study highlights the imaging potential of InSAR time series on complex hydrothermal systems.

     
    more » « less