skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Magma Chamber Growth During Intercaldera Periods: Insights From Thermo‐Mechanical Modeling With Applications to Laguna del Maule, Campi Flegrei, Santorini, and Aso
Abstract Crustal magma chambers can grow to be hundreds to thousands of cubic kilometers, potentially feeding catastrophic caldera‐forming eruptions. Smaller volume chambers are expected to erupt frequently and freeze quickly; a major outstanding question is how magma chambers ever grow to the sizes required to sustain the largest eruptions on Earth. We use a thermo‐mechanical model to investigate the primary factors that govern the extrusive:intrusive ratio in a chamber, and how this relates to eruption frequency, eruption size, and long‐term chamber growth. The model consists of three fundamental timescales: the magma injection timescaleτin, the cooling timescaleτcool, and the timescale for viscous relaxation of the crustτrelax. We estimate these timescales using geologic and geophysical data from four volcanoes (Laguna del Maule, Campi Flegrei, Santorini, and Aso) to compare them with the model. In each of these systems,τinis much shorter thanτcooland slightly shorter thanτrelax, conditions that in the model are associated with efficient chamber growth and simultaneous eruption. In addition, the model suggests that the magma chambers underlying these volcanoes are growing at rates between ~10−4and 10−2 km3/year, speeding up over time as the chamber volume increases. We find scaling relationships for eruption frequency and size that suggest that as chambers grow and volatiles exsolve, eruption frequency decreases but eruption size increases. These scaling relationships provide a good match to the eruptive history from the natural systems, suggesting that the relationships can be used to constrain chamber growth rates and volatile saturation state from the eruptive history alone.  more » « less
Award ID(s):
1760004
PAR ID:
10457136
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
20
Issue:
3
ISSN:
1525-2027
Page Range / eLocation ID:
p. 1574-1591
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Magmatic volatiles drive pressure, temperature, and compositional changes in upper crustal magma chambers and alter the physical properties of stored magmas. Previous studies suggest that magmatic H2O content influences the growth and longevity of silicic chambers through regulating the size and frequency of eruptions and impacting the crystallinity‐temperature curve. However, there has been comparatively little exploration of how CO2impacts the evolution of magma chambers despite the strong influence of CO2on H2O solubility and the high concentrations of CO2often present in mafic systems. In this study, we integrate the thermodynamic effects of dissolved and exsolved H2O and CO2with the mechanics of open‐system magma chambers that interact thermally and mechanically with the crust. We applied this model to investigate how intrinsic variations in magmatic H2O‐CO2content influence the growth and longevity of silicic and mafic magma chambers. Our findings indicate that even with a tenfold increase in CO2content (up to 10,000 ppm), CO2plays a minimal role in long‐term chamber growth and longevity. While CO2content affects the magma compressibility, the resulting changes in eruption mass are balanced out by a commensurate change in eruption frequency so that the time‐averaged eruptive flux and long‐term chamber behavior remain similar. In contrast, H2O content strongly influences chamber growth and longevity. In silicic systems, high H2O contents hinder magma chamber growth by increasing the total eruptive flux and steepening the slope of the crystallinity‐temperature curve. In mafic systems, high H2O contents promote magma chamber growth by flattening the slope of the crystallinity‐temperature curve. 
    more » « less
  2. Abstract Water and carbon dioxide are the most abundant volatile components in terrestrial magmas. As they exsolve into magmatic vapour, they promote magma buoyancy, accelerating ascent and modulating eruptive dynamics. It is commonly thought that an increase in pre-eruptive volatile content produces an increase in eruption intensity. Using a conduit model for basaltic eruptions, covering the upper 6 km of conduit, we show that for the same chamber conditions mass eruption rate is not affected by CO2content, whereas an increase in H2O up to 10 wt.% produces an increase in eruption rate of an order of magnitude. It is only when CO2is injected in the magma reservoir from an external source that the resulting pressurisation will generate a strong increase in eruption rate. Results also show that ascent velocity and fragmentation depth are strongly affected by pre-eruptive volatile contents demonstrating a link between volatile content and eruptive style. 
    more » « less
  3. Abstract We present a model for a coupled magma chamber–dike system to investigate the conditions required to initiate volcanic eruptions and to determine what controls the size of eruptions. The model combines the mechanics of dike propagation with internal chamber dynamics including crystallization, volatile exsolution, and the elastic response of the magma and surrounding crust to pressure changes within the chamber. We find three regimes for dike growth and eruptions: (1) below a critical magma chamber size, eruptions are suppressed because chamber pressure drops to lithostatic before a dike reaches the surface; (2) at an intermediate chamber size, the erupted volume is less than the dike volume (“dike-limited” eruption regime); and (3) above a certain chamber size, dikes can easily reach the surface and the erupted volume follows a classic scaling law, which depends on the attributes of the magma chamber (“chamber-limited” eruption regime). The critical chamber volume for an eruption ranges from ∼0.01 km3 to 10 km3 depending on the water content in the magma, depth of the chamber, and initial overpressure. This implies that the first eruptions at a volcano likely are preceded by a protracted history of magma chamber growth at depth, and that the crust above the magma chamber may have trapped several intrusions or “failed eruptions.” Model results can be combined with field observations of erupted volume, pressure, and crystal and volatile content to provide tighter constraints on parameters such as the eruptible chamber size. 
    more » « less
  4. Abstract Explosive volcanic eruptions radiate seismic waves as a consequence of pressure and shear traction changes within the conduit/chamber system. Kinematic source inversions utilize these waves to determine equivalent seismic force and moment tensor sources, but relation to eruptive processes is often ambiguous and nonunique. In this work, we provide an alternative, forward modeling approach to calculate moment tensor and force equivalents of a model of eruptive conduit flow and chamber depressurization. We explain the equivalence of two seismic force descriptions, the first in terms of traction changes on conduit/chamber walls, and the second in terms of changes in magma momentum, weight, and momentum transfer to the atmosphere. Eruption onset is marked by a downward seismic force, associated with loss of restraining shear tractions from fragmentation. This is followed by a much larger upward seismic force from upward drag of ascending magma and reduction of magma weight remaining in the conduit/chamber system. The static force is upward, arising from weight reduction. We calculate synthetic seismograms to examine the expression of eruptive processes at different receiver distances. Filtering these synthetics to the frequency band typically resolved by broadband seismometers produces waveforms similar to very long period seismic events observed in strombolian and vulcanian eruptions. However, filtering heavily distorts waveforms, accentuating processes in early, unsteady parts of eruptions and eliminating information about longer (ultra long period time scale depressurization and weight changes that dominate unfiltered seismograms. Our workflow can be utilized to directly and quantitatively connect eruption models with seismic observations. 
    more » « less
  5. Abstract Recent observations have reported that magnetosonic waves can exhibit rising‐tone structures in the frequency‐time spectrogram. However, the generation mechanism has not been identified yet. In this paper, we investigate the generation of rising‐tone magnetosonic waves in the terrestrial magnetosphere using 1‐D particle‐in‐cell (PIC) simulations, in which the plasma consists of three components: cool electrons, cool protons and ring distribution protons. We find that the magnetosonic waves excited by the ring distribution protons can form a rising‐tone structure with frequency of the structure ranging from about0.5Ωlhto nearlyΩlh, whereΩlhis the lower hybrid frequency. It is further demonstrated that the rising frequency of magnetosonic waves can be accounted for by the scattering of ring distribution protons. Moreover, the rising‐tone timescale obtained by PIC simulation is compared with the satellite observation. Our findings provide some new insights to understand the nonlinear evolution of plasma waves in the Earth's magnetosphere. 
    more » « less