Subglacial Environments Research Articles

Hematite U-Pb dating of Snowball Earth meltwater events

The Snowball Earth hypothesis predicts global ice cover; however, previous descriptions of Cryogenian (720-635 Ma) glacial deposits are limited to continental margins and shallow marine basins. The Tavakaiv (Tava) sandstone injectites and ridges in Colorado, USA, preserve a rare terrestrial record of Cryogenian low-latitude glaciation. Injectites, ridges, and chemically weathered crystalline rock display features characteristic of fluidization and pervasive deformation in a subglacial environment due to glacial loading, fluid overpressure, and repeated sand injection during meltwater events. In situ hematite U-Pb geochronology on hematite-quartz veins, which crosscut and are cut by Tava dikes, constrain sand injection at ~690-660 Ma. We attribute early Tava sand injection episodes to basal melting associated with rifting and geothermal heating, and later injections to meltwater generation during ~661 Ma Sturtian deglaciation. A modern analog is provided by the Ross Embayment of Antarctica, where rift-related faults border sediment-filled basins, overpressurized fluids circulate in confined aquifers below ice, and extensive preglacial topography is preserved. Field evidence and geochronology in Colorado further highlight that deep chemical weathering of Proterozoic bedrock and denudation associated with the Great Unconformity predate Cryogenian injection of fluidized sand, consistent with limited glacial erosion.

Profiling The Growth Conditions and Persistent Organic Pollutants (POPS) Tolerance of Phenoliferia glacialis USM-PSY62

Antarctica is characterized by extreme cold, isolated, and unique ecosystems. Nevertheless, Antarctica harbors diverse species of microorganisms, particularly in its ice-covered lakes and subglacial environments. These microorganisms have special adaptations to extreme cold and low-nutrient conditions. Some extremophiles, like psychrophiles can thrive in these harsh environments. Phenoliferia glacialis USM-PSY62, previously identified as Rhodotorula sp. USM-PSY62 is a psychrophilic yeast isolated from the ice brine of Antarctica. However, there is very little information on this psychrophilic yeast. This study aims to characterize the P. glacialis USM-PSY62 through the identification of the optimum growth parameters in different media (Yeast Peptone Dextrose, YPD & Yeast Malt, YM), temperature (4°C, 15°C, 20°C) and pH (6, 7, 8, 9) as well as their ability in carbon assimilation and extracellular enzyme production. It has an optimal growth in YPD compared to YM broth media. P. glacialis USM-PSY62 grows optimally at 15°C and pH 7.0. This Antarctic yeast enters the stationary phase on day six of incubation under optimum conditions. It appeared mainly as elongated-shape and oval-shaped with budding formation and was found to produce extracellular enzymes such as protease and amylase in the presence of 2% glucose concentration in YM media. P. glacialis USM-PSY62 also can assimilate various types of carbon sources including raffinose, arabinose, and maltose. Interestingly, the psychrophilic yeast presented growth in media supplemented with Persistent Organic Pollutants (POPs) such as dichlorophenyldichloroethylene (DDE) and polychlorinated biphenyl (PCB). These preliminary findings suggest that P. glacialis USM-PSY62 has tremendous potential for bioremediation application in polluted cold regions, as well as deepening our knowledge of its optimal growth conditions.

Characteristics of methane and carbon dioxide in ice caves at a high-mountain glacier of China

The exploration of the spatiotemporal distribution of greenhouse gas (GHG) exchange in the cryosphere (including ice sheet, glaciers, and permafrost) is important for understanding its future feedback to the atmosphere. Mountain glaciers and ice sheets may be potential sources of GHG emissions, but the magnitude and distribution of GHG emissions from glaciers and ice sheets remain unclear because observation data are lacking. In this study, in situ CH4 and CO2 and the mixing ratios of their carbon isotope signatures in the air inside an ice cave were measured, and CH4 and CO2 exchange in the meltwater of Laohugou glacier No. 12, a high-mountain glacier in an arid region of western China, was also analyzed and compared with the exchange in downstream rivers and a reservoir. The results indicated elevated CH4 mixing ratios (up to 5.7 ppm) and depleted CO2 (down to 168 ppm) in the ice cave, compared to ambient levels during field observations. The CH4 and CO2 fluxes in surface meltwater of the glacier were extremely low compared with their fluxes in rivers from the Tibetan Plateau (TP). CH4 and CO2 mixing ratios in the air inside the ice cave were mainly controlled by local meteorological conditions (air temperature, wind speed and direction) and meltwater runoff. The carbon isotopic compositions of CH4 and CO2 in the ice cave and terminus meltwater indicated δ13C-CH4 depletion compared to ambient air, suggesting an acetate fermentation pathway. The abundances of key genes for methanogenic archaea/genes encoding methyl coenzyme M reductase further indicated the production of CH4 by methanogenic archaea from the subglacial meltwater of high-mountain glaciers. The discovery of CH4 emissions from even small high-mountain glaciers indicates a more prevalent characteristic of glaciers to produce and release CH4 from the subglacial environment than previously believed. Nevertheless, further research is required to understand the relationship between this phenomenon and glacial dynamics in the third pole.

Insights into glacial processes from micromorphology of silt-sized sediment

Abstract. Silt-rich meltwater plume deposits (MPDs) analyzed from marine sediment cores have elucidated relationships that are clearly connected, yet difficult to constrain, between subglacial hydrology, ice-marginal landforms, and grounding-zone retreat patterns for several glacial catchments. Few attempts have been made to infer details of subglacial hydrology, such as flow regime, geometry of drainage pathways, and mode(s) of sediment transport through time, from grain-scale characteristics of MPDs. Using sediment samples from MPD, till, and grounding-zone proximal diamicton collected offshore of six modern and relict glacial catchments in both hemispheres, we examine grain shape distributions and microtextures (collectively, grain micromorphology) of the silt fraction to explore whether grains are measurably altered from their subglacial sources via meltwater action. We find that 75 % of all imaged grains (n = 9400) can be described by 25 % of the full range of measured shape morphometrics, indicating grain shape homogenization through widespread and efficient abrasive processes in subglacial environments. Although silt grains from MPDs exhibit edge rounding more often than silt grains from tills, grain surface textures indicative of fluvial transport (e.g., v-shaped percussions) occur in only a modest number of grains. Furthermore, MPD grain surfaces retain several textures consistent with transport beneath glacial ice (e.g., straight or arcuate steps, (sub)linear fractures) in comparable abundances to till grains. Significant grain shape alteration in MPDs compared to their till sources is observed in sediments from glacial regions where (1) high-magnitude, potentially catastrophic meltwater drainage events are inferred from marine sediment records and (2) submarine landforms suggest supraglacial melt contributed to the subglacial hydrological budget. This implies that quantifiable grain shape alteration in MPDs could reflect a combination of high-energy flow of subglacial meltwater, persistent sediment entrainment, and/or long sediment transport distances through subglacial drainage pathways. Integrating grain micromorphology into analysis of MPDs in site-specific studies could therefore aid in distinguishing periods of persistent, well-connected subglacial discharge from periods of sluggish or disorganized drainage. In the wider context of deglacial marine sedimentary and bathymetric records, a grain micromorphological approach may bolster our ability to characterize ice response to subglacial meltwater transmission through time. This work additionally demonstrates that glacial and fluvial surface textures are retained on silt-sized quartz grains in adequate amounts for microtexture analysis, which has heretofore been conducted exclusively on the sand fraction. Therefore, grain microtextures can be examined on silt-rich glaciogenic deposits that contain little to no sand as a means to evaluate sediment transport processes.

Subglacial bedform and moat initiation beneath Rutford Ice Stream, West Antarctica

The accuracy of sea-level rise predictions is highly dependent on reliably understanding the subglacial environment beneath ice streams. Bedforms result from the interaction between ice and its substrate, and therefore have the potential to shed light on poorly understood basal conditions influencing ice dynamics. However, theoretical models of bedform formation are often based on observations from deglaciated areas or from sparsely-sampled geophysical surveys over glaciated regions. Here, we use high-resolution three- and two-dimensional radar and seismic data to reveal details of the initiation and evolution of a subglacial bedform beneath Rutford Ice Stream, West Antarctica. Radar surveys at 20 m and 50 m line spacing allow detailed imaging of bed topography, including a moat up to 55 m deep, surrounding the upstream end of a 50-m high and >18-km long bedform. Many models rely on either a topographical or a locally resistant seed point to initiate bedform formation. The bedform described here is mostly composed of soft sediment (porosity >0.3) and lacks a boulder or outcrop, suggesting the bedform initiated without a topographic seed point. Sediment at the upstream end of bedforms appears stiffer over a distance of 2.3 km. We suggest sediment inhomogeneities in the initially flat bed cause the deposition of sediment, which, assuming resistant enough, acts as a seed point for bedform extension and moat erosion. The moat's geometry and its truncation of other bedforms suggest that it was eroded after the deposition of surrounding bedforms. These observations from a modern ice stream deliver information of subglacial processes involved in the initiation as well as in situ high-resolution topography and properties of bedforms and moats. Using these observations numerical models can be tested and developed accordingly.

Marine redox fluctuations during the Marinoan glaciation

The proliferation of eukaryotes preceding the Cryogenian Marinoan glaciation (650–635 Ma) and the subsequent radiation of early animals imply a conceivable linkage between global glaciation and complex life evolution. However, the marine redox condition remains enigmatic during the Cryogenian, impeding our understanding of the role of O2 on this biological innovation. To fill this gap, we comprehensively reported mineralogy, iron isotope (δ56Fe), and iron speciation in the Cryogenian Nantuo Formation at Liulongshan, South China. This section encompasses glacial marine and non-glacial deposits, where the redbeds situated in the middle of the section, in conjunction with the diamictites adjacent to them, divide the Nantuo Formation into three distinct glacial sedimentary cycles. Micro-particle hematite within the Nantuo Formation indicates that the primary Fe-oxide was authigenic and trapped in water column. Furthermore, the Fe isotopic compositions of these authigenic components (δ56FeHR) calculated by mass balance equation reveal the redox state of the ocean coinciding with dynamic glaciation. The lower diamictite shows an increasing trend of δ56FeHR, Fe/Al, and FeHR/FeT, indicative of gradually more anoxic seawater conditions coupled with ice advance. Although redbeds intervals exhibit higher FeHR/FeT ratios (∼ 0.51), traditionally indicative of anoxic environments, our integration of mineralogical evidence suggested that the increased Feox content in redbeds might represent increased oxidation of Fe(II) at the redox chemocline, rather than the redox conditions of the bottom waters. This interpretation is consistent with the observed decreasing trend in our Fe isotope data, as well as the sedimentological evolution. Subsequently, the positive and increasing δ56FeHR values suggest that the oceanic redox conditions shifted to anoxic as the second glacial episode re-emerged. These findings suggest dynamic fluctuations in the redox condition of the ocean during the Marinoan Ice Age. Specifically, the formation of the redbeds may reflect a locally oxic environment. A modern analogy resembling this situation is blood falls, where the oxidation of iron-rich brine plumes under the Antarctic ice sheets. Such subglacial environments could have potentially provided a habitable “refuge” for early life to survive during the Neoproterozoic glaciation.

Seasonal to decadal dynamics of supraglacial lakes on debris-covered glaciers in the Khumbu region, Nepal

Abstract. Supraglacial lakes (SGLs) play an important role in debris-covered glacier (DCG) systems by enabling efficient interactions between the supraglacial, englacial, and subglacial environments. Developing a better understanding of the short-term and long-term development of these features is needed to constrain DCG evolution and the hazards posed to downstream communities, ecosystems, and infrastructure from rapid drainage. In this study, we present an analysis of supraglacial lakes on eight DCGs in the Khumbu region of Nepal by automating SGL identification in PlanetScope, Sentinel-2, and Landsat 5–9 images. We identify a regular annual cycle in SGL area, with lakes covering approximately twice as much area during their maximum annual extent (in the pre-monsoon season) than their minimum annual extent (in the post-monsoon season). The high spatiotemporal resolution of PlanetScope imagery (∼ daily, 3 m) shows that this cycle is driven by the appearance and expansion of small lakes in the upper debris-covered regions of these glaciers throughout the winter. Decadal-scale expansion of large, near-terminus lakes was identified on four of the glaciers (Khumbu, Lhotse, Nuptse, and Ambulapcha), while the remaining four showed no significant increases over the study period. The seasonal variation in SGL area is of comparable or greater magnitude as decadal-scale changes, highlighting the importance of accounting for this seasonality when interpreting long-term records of SGL changes from sparse observations. The complex spatiotemporal patterns revealed in our analysis are not captured in existing regional-scale glacial lake databases, suggesting that more targeted efforts are needed to capture the true variability of SGLs on large scales.

Accelerating Subglacial Hydrology for Ice Sheet Models With Deep Learning Methods

AbstractSubglacial drainage networks regulate the response of ice sheet flow to surface meltwater input to the subglacial environment. Simulating subglacial hydrology evolution is critical to projecting ice sheet sensitivity to climate, and contribution to sea‐level change. However, current numerical subglacial hydrology models are computationally expensive, and, consequently, evolving subglacial hydrology is neglected in large‐scale ice sheet simulations. We present a deep learning emulator of a state‐of‐the‐art subglacial hydrology model, trained at multiple Greenland glaciers. Our emulator performs strongly in both temporal (R2 > 0.99) and spatial (R2 > 0.95) generalization, offers high computational savings, and can be used to force numerical ice sheet models. This will enable century‐ and large‐scale ice sheet model simulations, including interactions between ice flow and increased meltwater input to the subglacial environment. Generally, our work demonstrates that machine learning can further improve ice sheet models, reduce computational bottlenecks, and exploit information from high‐fidelity models and novel observational platforms.

Tidally Modulated Glacial Slip and Tremor at Helheim Glacier, Greenland

AbstractNumerical modeling of ice sheet motion and hence projections of global sea level rise require information about the evolving subglacial environment, which unfortunately remains largely unknown due to its difficulty of access. Here we advance such subglacial observations by reporting multi‐year observations of seismic tremor likely associated with glacier sliding at Helheim Glacier. This association is confirmed by correlation analysis between tremor power and multiple environmental forcings on different timescales. Variations of the observed tremor power indicate that different factors affect glacial sliding on different timescales. Effective pressure may control glacial sliding on long (seasonal/annual) timescales, while tidal forcing modulates the sliding rate and tremor power on short (hourly/daily) timescales. Polarization results suggest that the tremor source comes from an upstream subglacial ridge. This observation provides insights on how different factors should be included in ice sheet modeling and how their timescales of variability play an essential role.

Methane emissions from subglacial meltwater of three alpine glaciers in Yukon, Canada

ABSTRACT Subglacial meltwater of land-terminating glaciers in Greenland and Iceland are sources of methane (CH4) to the atmosphere, but sparse empirical data exist about the spatial distribution of subglacial CH4 production and emission from glaciers in other regions of the world. This study presents the first measurements of CH4 emissions from the subglacial meltwater of three outlet glaciers of the St. Elias Mountains in Yukon, Canada. Dissolved CH4 concentrations were highly elevated at 45, 135, and 250 times compared to the atmospheric equilibrium concentration in the meltwater of Dusty, Kluane, and Donjek glaciers, respectively. Dissolved CO2 concentrations were depleted relative to the atmospheric equilibrium. This points to the meltwater being a source of CH4 and a sink of CO2. Stable carbon (13C) and hydrogen (2H) isotopic signatures of the subglacial CH4 were depleted compared to atmospheric CH4 at all sites, indicating both biotic and abiotic sources and possible alteration from bacterial CH4 oxidation in the meltwater. No relation was found between CH4 concentrations in the meltwater and the meltwater chemistry or the size of the glaciers in this study. These findings suggest that CH4 emissions from subglacial environments under alpine glaciers may be a more common phenomenon than previously thought.

Outbursts from an ice-marginal lake in Antarctica in 1969–1971 and 2017, revealed by aerial photographs and satellite data

The liquid water around the Antarctic Ice Sheet plays a key role in modulating both the vulnerability of ice shelves to hydrofracturing and ice discharge from outlet glaciers. Therefore, it needs to be adequately constrained for precise future projections of ice-mass loss and global sea-level rise. Although glacial lake outburst floods (GLOFs) pose one of the greatest risks in glacierized mountainous regions, any long-term monitoring of Antarctic ice-marginal lakes and their associated potential for GLOFs has been neglected until recently owing to the limited number of such events reported in Antarctica. Here we present direct evidence of repeated GLOFs from Lake Kaminotani-Ike, an ice-sheet-dammed lake in East Antarctica, via an analysis of historical aerial photographs and recent satellite data. Two GLOFs occurred in 1969–1971 and 2017, with discharge volumes of (8.6 ± 1.5) × 107 and (7.1 ± 0.4) × 107 m3, respectively, making them two of the largest GLOFs in Antarctica. A southerly oceanward pathway beneath the ice sheet is the most likely drainage route of these GLOF events based on the available surface- and bed-elevation datasets. Furthermore, the 2017 event occurred during the austral winter, thereby implying the possibility of year-round active subglacial networks in Antarctica. Our results highlight that studies on Antarctic ice-marginal lakes provide an opportunity to better understand Antarctic hydrological processes and emphasize the need for both detailed monitoring of ice-marginal lakes and detailed surveying of the subglacial environments of the Antarctic Ice Sheet.

Subglacial discharge accelerates future retreat of Denman and Scott Glaciers, East Antarctica.

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.

Glacial secrets uncovered: Revealing the modes of survival of metabolically active microbial communities entrapped in polar glacial ice

Glacial secrets uncovered: Revealing the modes of survival of metabolically active microbial communities entrapped in polar glacial ice

Methylotrophic Communities Associated with a Greenland Ice Sheet Methane Release Hotspot

Subglacial environments provide conditions suitable for the microbial production of methane, an important greenhouse gas, which can be released from beneath the ice as a result of glacial melting. High gaseous methane emissions have recently been discovered at Russell Glacier, an outlet of the southwestern margin of the Greenland Ice Sheet, acting not only as a potential climate amplifier but also as a substrate for methane consuming microorganisms. Here, we describe the composition of the microbial assemblage exported in meltwater from the methane release hotspot at Russell Glacier and its changes over the melt season and as it travels downstream. We found that a substantial part (relative abundance 27.2% across the whole dataset) of the exported assemblage was made up of methylotrophs and that the relative abundance of methylotrophs increased as the melt season progressed, likely due to the seasonal development of the glacial drainage system. The methylotrophs were dominated by representatives of type I methanotrophs from the Gammaproteobacteria; however, their relative abundance decreased with increasing distance from the ice margin at the expense of type II methanotrophs and/or methylotrophs from the Alphaproteobacteria and Betaproteobacteria. Our results show that subglacial methane release hotspot sites can be colonized by microorganisms that can potentially reduce methane emissions.

Constraints on subglacial melt fluxes from observations of active subglacial lake recharge

AbstractActive subglacial lakes provide a rare glimpse of the subglacial environment and hydrological processes at play. Several studies contributed to establishing active subglacial lake inventories and document lake drainage and connection, but few focused on the period between lake drainage when the melt production and transport contribute to the refilling of these lakes. In this study, we employ high-resolution CryoSat-2 altimetry data from 2010 to 2021 to compile an inventory of recharging lakes across Antarctica. We extract recharge rates from these lakes, which serve as a lower limit on subglacial melt production. These recharge rates are compared against predictions obtained by routing modelled subglacial meltwater at the ice-sheet's base. Our findings indicate that modelled recharge rates are consistent with observations in all but one of the investigated lakes, providing a lower bound on geothermal heat fluxes. Lake Cook E2 displays recharge rates far exceeding predictions, indicating that processes are taking place that are currently unaccounted for. Considering recharge in hydrologically connected lake networks instead of individually provides a stricter constraint on melt production. Recharge rates extracted from the Thwaites Lake system suggest that subglacial melt production has been underestimated.

Antarctic Sedimentary Basins and Their Influence on Ice‐Sheet Dynamics

AbstractKnowledge of Antarctica's sedimentary basins builds our understanding of the coupled evolution of tectonics, ice, ocean, and climate. Sedimentary basins have properties distinct from basement‐dominated regions that impact ice‐sheet dynamics, potentially influencing future ice‐sheet change. Despite their importance, our knowledge of Antarctic sedimentary basins is restricted. Remoteness, the harsh environment, the overlying ice sheet, ice shelves, and sea ice all make fieldwork challenging. Nonetheless, in the past decade the geophysics community has made great progress in internationally coordinated data collection and compilation with parallel advances in data processing and analysis supporting a new insight into Antarctica's subglacial environment. Here, we summarize recent progress in understanding Antarctica's sedimentary basins. We review advances in the technical capability of radar, potential fields, seismic, and electromagnetic techniques to detect and characterize basins beneath ice and advances in integrated multi‐data interpretation including machine‐learning approaches. These new capabilities permit a continent‐wide mapping of Antarctica's sedimentary basins and their characteristics, aiding definition of the tectonic development of the continent. Crucially, Antarctica's sedimentary basins interact with the overlying ice sheet through dynamic feedbacks that have the potential to contribute to rapid ice‐sheet change. Looking ahead, future research directions include techniques to increase data coverage within logistical constraints, and resolving major knowledge gaps, including insufficient sampling of the ice‐sheet bed and poor definition of subglacial basin structure and stratigraphy. Translating the knowledge of sedimentary basin processes into ice‐sheet modeling studies is critical to underpin better capacity to predict future change.

Episodic Subglacial Drainage Cascades Below the Northeast Greenland Ice Stream

AbstractSubglacial hydrology can exert an important control on ice flow by affecting friction at the ice‐bedrock interface. Here, we report on a series of subglacial drainage events along the Northeast Greenland Ice Stream (NEGIS), initiating as far inland as 500 km from the margin of Zachariae Isstrøm. The drainage events exhibit local transient uplift, followed by prolonged subsidence, measured by differential satellite synthetic aperture radar interferometry (DInSAR). In downstream regions, drainage events are associated with temporary acceleration in ice flow. The high spatiotemporal resolution of the DInSAR measurements allows for a detailed mapping of the drainage propagation pathway. We show that multiple drainage cascades have occurred along the same pathway over the years 2020–2022. Finally, the propagation speed of subglacial water flow is found to vary greatly along NEGIS, suggesting that fundamental differences could exist in the subglacial environment.

A thermomechanical model for frost heave and subglacial frozen fringe

Ice-infiltrated sediment, or frozen fringe, is responsible for phenomena such as frost heave, ice lenses and metres of debris-rich ice under glaciers. Understanding frozen fringes is important as frost heave is responsible for damaging infrastructure at high latitudes and sediment freeze-on at the base of glaciers can modulate subglacial friction, influencing the rate of global sea level rise. Here we describe the thermomechanics of liquid water flow and freezing in ice-saturated sediments, focusing on the conditions relevant for subglacial environments. The force balance that governs the frozen fringe thickness depends on the weight of the overlying material, the thermomolecular force between ice and sediments across liquid premelted films and the water pressure required by Darcy flow. We combine this mechanical model with an enthalpy method that conserves energy across phase change interfaces on a fixed computational grid. The force balance and enthalpy model together determine the evolution of the frozen fringe thickness and our simulations predict frost heave rates and ice lens spacing. Our model accounts for premelting at ice–sediment contacts, partial ice saturation of the pore space, water flow through the fringe, the thermodynamics of the ice–water–sediment interface and vertical force balance. We explicitly account for the formation of ice lenses, regions of pure ice that cleave the fringe at the depth where the interparticle force vanishes. Our model results allow us to predict the thickness of a frozen fringe and the spacing of ice lenses in subaerial and subglacial sediments.

Trace element, rare earth element and trace carbon compounds in Subglacial Lake Whillans, West Antarctica

Whillans Subglacial Lake (SLW) lies beneath 801 m of ice in the lower portion of the Whillans Ice Stream (WIS) in West Antarctica and is part of an extensive and active subglacial drainage network. Here, the geochemical characterization of SLW rare earth elements (REE), trace elements (TE), free amino acids (FAA), and phenolic compounds (PC) measured in lakewater and sediment porewater are reported. The results show, on average, higher values of REEs in the lakewater than in the porewater, and clear changes in all REE concentrations and select redox sensitive trace element concentrations in porewaters at a depth of ~15 cm in the 38 cm lake sediment core. This is consistent with prior results on the lake sediment redox conditions based on gas chemistry and microbiological data. Low concentrations of vanillyl phenols were measured in the SLW water column with higher concentrations in porewater samples and their concentration profiles in the sediments may also reflect changing redox conditions in the sediments. Vanillin concentrations increased with depth in the sediments as oxygenation decreases, while the concentrations of vanillic acid, the more oxidized component, were higher in the more oxygenated surface sediments. Collectively these results indicate redox changes occurring with the upper 38 cm of sediment in SLW and provide support for the existence of a seawater source, already hypothesized, in the sediments below the lowest measured depth, and of a complex and dynamic geochemical system beneath the West Antarctic Ice Sheet.Our results are the first to detail geochemical properties from an Antarctic subglacial environment using direct sampling technology. Due to their isolation from the wider environment, subglacial lakes represent one of our planets last pristine environments that provide habitats for microbial life and natural biogeochemical cycles but also impact the basal hydrology and can cause ice flow variations.

Flash heating boosts the potential for mechanochemical energy sources for subglacial ecosystems

Subglacial environments harbour a diversity of microbial ecosystems capable of influencing biogeochemical cycles. However, the darkness and isolation of subglacial environments limit the energy sources available for microbial metabolism. A recently recognised energy source for these microbes in wet-based regions is the rock-water reactions that occur after the mechanical fracturing of glacial bedrock. These mechanochemical reactions produce H2 and H2O2 at 0°C from reactions with mineral surface defects (Si• and SiO•) and release Fe from within the mineral structures, providing electron donors and acceptors for microbial metabolism. However, the production of H2O2 and H2 may be underestimated as temperatures at rock abrasion sites can increase substantially above 0°C as glaciers “slip and grind” rocks, potentially accelerating the rates of mechanochemical reactions. Despite this, the effect of rapid heating on subsequent low-temperature mechanochemical reactions has yet to be examined. Here, we investigate H2, H2O2, and Fe production during low-temperature (0 °C) incubations of water with a range of ground rocks and minerals following “flash heating” to 30, 60, or 121 °C. We show that transient increases (as little as 5–10 min of heating) to moderate temperatures (30 or 60 °C) can significantly increase the rate of H2 production, while short-term heating to 121 °C generates larger bursts of H2. In addition, pyrite is easily crushed, potentially releasing large quantities of Fe2+ into subglacial systems and promoting mechanochemical reactions due to the resulting large surface area (10× larger than other materials). We provide the first evidence for H2 production from water reactions with crushed pyrite and suggest that crushed pyrite has a greater influence on subglacial H2O2 production than silicates. We conclude that electron donors in the form of Fe2+ and H2 bursts can be produced in subglacial ecosystems, which may be coupled to substantial concentrations of H2O2 produced from crushed pyrite. This suggests that rock–water mechanochemical reactions may be a greater source of energy for subglacial environments than previously recognised.