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 Response to Ice Unloading: Applications to Volcanism in the West Antarctic Rift System
Abstract Volcanic activity has been shown to affect Earth's climate in a myriad of ways. One such example is that eruptions proximate to surface ice will promote ice melting. In turn, the crustal unloading associated with melting an ice sheet affects the internal dynamics of the underlying magma plumbing system. Geochronologic data from the Andes over the last two glacial cycles suggest that glaciation and volcanism may interact via a positive feedback loop. At present, accurate sea‐level predictions hinge on our ability to forecast the stability of the West Antarctic Ice Sheet, and thus require consideration of two‐way subglacial volcano‐deglaciation processes. The West Antarctic Ice Sheet is particularly vulnerable to collapse, yet its position atop an active volcanic rift is seldom considered. Ice unloading deepens the zone of melting and alters the crustal stress field, impacting conditions for dike initiation, propagation, and arrest. However, the consequences for internal magma chamber dynamics and long‐term eruption behavior remain elusive. Given that unloading‐triggered volcanism in West Antarctica may contribute to the uncertainty of ice loss projections, we adapt a previously published thermomechanical magma chamber model and simulate a shrinking ice load through a prescribed lithostatic pressure decrease. We investigate the impacts of varying unloading scenarios on magma volatile partitioning and eruptive trajectory. Considering the removal of km‐thick ice sheets, we demonstrate that the rate of unloading influences the cumulative mass erupted and consequently the heat released into the ice. These findings provide fundamental insights into the complex volcano‐ice interactions in West Antarctica and other subglacial volcanic settings.  more » « less
Award ID(s):
2121570 2444709
PAR ID:
10557234
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
25
Issue:
12
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. While the effects of volcanism on Earth’s climate are well understood, the volcano-ice sheet system hosts a two-way feedback. Volcanic activity promotes ice melting, which in turn affects the internal dynamics of the magma chamber below. At present, accurate forecasts of sea-level rise hinge on the stability of the West Antarctic Ice Sheet, and thus require consideration of subglacial volcano-deglaciation feedbacks. The West Antarctic Ice Sheet, grounded below sea-level, is particularly vulnerable to collapse, yet its position atop an active volcanic rift is seldom considered. Ice unloading raises the geotherm and alters the crustal stress field, impacting dike propagation. However, the consequences on internal magma chamber dynamics and thus long-term eruption behavior remain elusive. Given potential for unloading-triggered volcanism in West Antarctica to accelerate ice retreat, we adapt the thermomechanical magma chamber model of Scholz et al. (2023) for West Antarctic Rift basalts, simulating a shrinking ice load through a prescribed decrease of lithostatic pressure. Examining different unloading scenarios, we investigate the impacts on volatile partitioning within the magma and eruptive trajectory across a wide range of initial magma chamber conditions. Pressurization of a magma chamber beyond a critical threshold results in eruption, delivering enthalpy to the ice. Considering the removal of km-thick ice sheets, we demonstrate the rate of unloading is dominant in influencing the cumulative mass erupted and consequently, heat released to the ice. These findings provide fundamental insights into the complex volcano-ice interactions in West Antarctica and other subglacial volcanic settings. 
    more » « less
  2. Troll, Valentin (Ed.)
    Continental alkaline magmatism produces a wide variety of igneous rock types because of varying degrees of partial melting of heterogenous mantle sources, fractional crystallization, and magma contamination during transit through the continental crust. The Mount Overlord Volcanic Field (MOVF) is a continental alkaline volcanic province in northern Victoria Land, Antarctica. Mount Overlord and the associated vents that make up the volcanic field are some of the least-explored volcanic rocks in the western Ross Sea. The MOVF sits within the Transantarctic Mountains, which form the rift shoulder of the extensive West Antarctic Rift System. The compositions of volcanic rocks in the MOVF range widely from basanite to evolved trachyte and comendite with a suite of intermediate rock types. Here we present 40Ar/39Ar ages, petrography, and whole-rock and mineral geochemistry to establish the temporal and magmatic evolution of the magmatic system. Volcanic activity occurred from 21.2 to 6.9 Ma, making it one of the longest records of volcanism in the western Ross Sea area. Mount Rittmann, an active volcano that is part of the MOVF, is not discussed here but extends the timing of volcanism of the MOVF into the Holocene. At Mount Overlord and surrounding areas, there were eruptions of lava flows, domes, and pyroclastic rocks. Localized deposits of hyaloclastites formed by magma-ice interactions provide an insight into former ice levels. Geochemically the MOVF shows a single magma differentiation trend except for Navigator Nunatak lavas which have a potassic affinity rarely seen in northern Victoria Land. Partial melting of an amphibole-bearing mantle lithology at or near the base of the continental lithospheric mantle (CLM) was the main source of the parental basaltic magmas. Polybaric crystal fractionation of the primary basaltic magmas mainly occurred at lower crustal depths and involved fractionation of clinopyroxene, olivine, kaersutite, feldspars, biotite, Fe–Ti oxides, apatite, and sodalite. Crustal assimilation of c. 10% granite harbor igneous complex granitoids was important in the evolution of intermediate composition magmas. Trachyte, phonolite, and comendite magmas stagnated and evolved at shallow crustal depths (c. <8 km). Over 95% crystal fractionation was required to generate the comendites. Extraction of the comendite melt from a felsic crystal mush was an important process. The potassic Navigator Nunatak magma required partial melting of phlogopite-bearing metasomatized CLM. The metasomes had ‘HIMU-like’ or FOZO isotopic compositions that ultimately originated from recycling of materials in the mantle. The MOVF displays a stronger affinity toward FOZO than other northern Victoria Land basaltic rocks. This suggests that the interaction between parental melt and juvenile CLM was limited, which is similar to volcanic rocks from the oceanic Adare Basin seamounts. Our result emphasizes the critical importance of a thick CLM for the genesis of diverse alkaline magma compositions in a continental rift system. 
    more » « less
  3. Ice shelf basal melting is the primary mechanism driving mass loss from the Antarctic Ice Sheet, yet it is unknown how the localized melt enhancement from subglacial discharge will affect future Antarctic glacial retreat. We develop a parameterization of ice shelf basal melt that accounts for both ocean and subglacial discharge forcing and apply it in future projections of Denman and Scott Glaciers, East Antarctica, through 2300. In forward simulations, subglacial discharge accelerates the onset of retreat of these systems into the deepest continental trench on Earth by 25 years. During this retreat, Denman Glacier alone contributes 0.33 millimeters per year to global sea level rise, comparable to half of the contemporary sea level contribution of the entire Antarctic Ice Sheet. Our results stress the importance of resolving complex interactions between the ice, ocean, and subglacial environments in future Antarctic Ice Sheet projections. 
    more » « less
  4. Mocho-Choshuenco volcano (39.9°S, 72.0°W) produced ∼75 explosive eruptions following retreat of the >1.5-km-thick Patagonian Ice Sheet associated with the local Last Glacial Maximum (LGM, from 35 to 18 ka). Here, we extend this record of volcanic evolution to include pre- and syn-LGM lavas that erupted during the Pleistocene. We establish a long-term chronology of magmatic and volcanic evolution and evaluate the relationship between volcanism and loading/unloading of the Patagonian Ice Sheet via twenty-four 40Ar/39Ar and two 3He age determinations integrated with stratigraphy and whole-rock compositions of lava flows and glass compositions of tephra. Our findings reveal that the edifice is much younger than previously thought and preserves 106 km3 of eruptive products, of which 50% were emplaced immediately following the end of the penultimate glaciation and 20% after the end of the LGM. A period of volcanic inactivity between 37 and 26 ka, when glaciers expanded, was followed by the eruption of incompatible element-rich basaltic andesites. Several of these syn-LGM lavas dated between 26 and 16 ka, which crop out at 1500−1700 m above sea level, show ice contact features that are consistent with emplacement against a 1400- to 1600-m-thick Patagonian Ice Sheet. Small volume dacitic eruptions and two explosive rhyolitic eruptions dominate the volcanic output from 18 to 8 ka, when the Patagonian Ice Sheet began to retreat rapidly. We hypothesize that increased lithostatic loading as the Patagonian Ice Sheet grew prohibited dike propagation, thus stalling the ascent of magma, promoting growth of at least three discrete magma reservoirs, and enhancing minor crustal assimilation to generate incompatible element-rich basaltic andesitic to dacitic magmas that erupted between 26 and 17 ka. From an adjacent reservoir, incompatible element-poor dacites erupted from 17 to 12 ka. These lava flows were followed by the caldera-forming eruption at 11.5 ka of 5.3 km3 of rhyolite from a deeper reservoir atop which a silicic melt lens had formed and expanded. Subsequent eruptions of oxidized dacitic magmas from the Choshuenco cone from 11.5 to 8 ka were followed by andesitic to dacitic eruptions at the more southerly Mocho cone, as well as small flank vent eruptions of basaltic andesite at 2.5 and 0.5 ka. This complex history reflects a multi-reservoir plumbing system beneath Mocho-Choshuenco, which is characterized by depths of magma storage, oxidation states, and trace element compositions that vary over short periods of time (<2 k.y.). 
    more » « less
  5. Antarctic subglacial lakes can play an important role in ice sheet dynamics, biology, geology, and oceanography, but it is difficult to definitively constrain their character and locations. Subglacial lake locations are related to factors including heat flux, ice surface slope, ice thickness, and bed topography, though these relationships are not fully quantified. Bed topography is particularly important for determining where water flows and accumulates, but digital elevation models of the ice sheet bed rely on interpolation and are unrealistically smooth, biasing estimates of subglacial lake location and surface area. To address this issue, we use geostatistical methods to simulate realistically rough bed topography. We use our simulated topography to predict subglacial lake distribution across the continent using a binomial logistic regression, which uses physical parameters and known lake locations to calculate the probabilities of lake occurrences. Our results suggest that topography models interpolated without appropriate geostatistics overestimate subglacial lake surface area and that total lake surface area is lower than previously predicted. We find that radar‐detected lakes are more likely to occur in the interior of East Antarctica, while altimetry‐detected (active) lakes are expected to be found in West Antarctica and near the grounding line. We observe that radar‐detected lakes have a high correlation with heat flux and ice thickness, while active lakes are associated with higher ice velocity. 
    more » « less