Greenland melt drives continuous export of methane from the ice-sheet bed.

Ice sheets are currently ignored in global methane budgets1,2. Although ice sheets have been proposed to contain large reserves of methane that may contribute to a rise in atmospheric methane concentration if released during periods of rapid ice retreat3,4, no data exist on the current methane footprint of ice sheets. Here we find that subglacially produced methane is rapidly driven to the ice margin by the efficient drainage system of a subglacial catchment of the Greenland ice sheet. We report the continuous export of methane-supersaturated waters (CH4(aq)) from the ice-sheet bed during the melt season. Pulses of high CH4(aq) concentration coincide with supraglacially forced subglacial flushing events, confirming a subglacial source and highlighting the influence of melt on methane export. Sustained methane fluxes over the melt season are indicative of subglacial methane reserves that exceed methane export, with an estimated 6.3 tonnes (discharge-weighted mean; range from 2.4 to 11 tonnes) of CH4(aq) transported laterally from the ice-sheet bed. Stable-isotope analyses reveal a microbial origin for methane, probably from a mixture of inorganic and ancient organic carbon buried beneath the ice. We show that subglacial hydrology is crucial for controlling methane fluxes from the ice sheet, with efficient drainage limiting the extent of methane oxidation5 to about 17 per cent of methane exported. Atmospheric evasion is the main methane sink once runoff reaches the ice margin, with estimated diffusive fluxes (4.4 to 28 millimoles of CH4 per square metre per day) rivalling that of major world rivers6. Overall, our results indicate that ice sheets overlie extensive, biologically active methanogenic wetlands and that high rates of methane export to the atmosphere can occur via efficient subglacial drainage pathways. Our findings suggest that such environments have been previously underappreciated and should be considered in Earth's methane budget.

Algal photophysiology drives darkening and melt of the Greenland Ice Sheet.

Blooms of Zygnematophycean "glacier algae" lower the bare ice albedo of the Greenland Ice Sheet (GrIS), amplifying summer energy absorption at the ice surface and enhancing meltwater runoff from the largest cryospheric contributor to contemporary sea-level rise. Here, we provide a step change in current understanding of algal-driven ice sheet darkening through quantification of the photophysiological mechanisms that allow glacier algae to thrive on and darken the bare ice surface. Significant secondary phenolic pigmentation (11 times the cellular content of chlorophyll a) enables glacier algae to tolerate extreme irradiance (up to ∼4,000 µmol photons⋅m-2⋅s-1) while simultaneously repurposing captured ultraviolet and short-wave radiation for melt generation. Total cellular energy absorption is increased 50-fold by phenolic pigmentation, while glacier algal chloroplasts positioned beneath shading pigments remain low-light-adapted (Ek ∼46 µmol photons⋅m-2⋅s-1) and dependent upon typical nonphotochemical quenching mechanisms for photoregulation. On the GrIS, glacier algae direct only ∼1 to 2.4% of incident energy to photochemistry versus 48 to 65% to ice surface melting, contributing an additional ∼1.86 cm water equivalent surface melt per day in patches of high algal abundance (∼104 cells⋅mL-1). At the regional scale, surface darkening is driven by the direct and indirect impacts of glacier algae on ice albedo, with a significant negative relationship between broadband albedo (Moderate Resolution Imaging Spectroradiometer [MODIS]) and glacier algal biomass (R2 = 0.75, n = 149), indicating that up to 75% of the variability in albedo across the southwestern GrIS may be attributable to the presence of glacier algae.

Decadal slowdown of a land-terminating sector of the Greenland Ice Sheet despite warming.

Ice flow along land-terminating margins of the Greenland Ice Sheet (GIS) varies considerably in response to fluctuating inputs of surface meltwater to the bed of the ice sheet. Such inputs lubricate the ice-bed interface, transiently speeding up the flow of ice. Greater melting results in faster ice motion during summer, but slower motion over the subsequent winter, owing to the evolution of an efficient drainage system that enables water to drain from regions of the ice-sheet bed that have a high basal water pressure. However, the impact of hydrodynamic coupling on ice motion over decadal timescales remains poorly constrained. Here we show that annual ice motion across an 8,000-km(2) land-terminating region of the west GIS margin, extending to 1,100 m above sea level, was 12% slower in 2007-14 compared with 1985-94, despite a 50% increase in surface meltwater production. Our findings suggest that, over these three decades, hydrodynamic coupling in this section of the ablation zone resulted in a net slowdown of ice motion (not a speed-up, as previously postulated). Increases in meltwater production from projected climate warming may therefore further reduce the motion of land-terminating margins of the GIS. Our findings suggest that these sectors of the ice sheet are more resilient to the dynamic impacts of enhanced meltwater production than previously thought.

Greenland ice sheet motion insensitive to exceptional meltwater forcing.

Changes to the dynamics of the Greenland ice sheet can be forced by various mechanisms including surface-melt-induced ice acceleration and oceanic forcing of marine-terminating glaciers. We use observations of ice motion to examine the surface melt-induced dynamic response of a land-terminating outlet glacier in southwest Greenland to the exceptional melting observed in 2012. During summer, meltwater generated on the Greenland ice sheet surface accesses the ice sheet bed, lubricating basal motion and resulting in periods of faster ice flow. However, the net impact of varying meltwater volumes upon seasonal and annual ice flow, and thus sea level rise, remains unclear. We show that two extreme melt events (98.6% of the Greenland ice sheet surface experienced melting on July 12, the most significant melt event since 1889, and 79.2% on July 29) and summer ice sheet runoff ~3.9 σ above the 1958-2011 mean resulted in enhanced summer ice motion relative to the average melt year of 2009. However, despite record summer melting, subsequent reduced winter ice motion resulted in 6% less net annual ice motion in 2012 than in 2009. Our findings suggest that surface melt-induced acceleration of land-terminating regions of the ice sheet will remain insignificant even under extreme melting scenarios.

Alpine glacier algal bloom during a record melt year.

Glacier algal blooms dominate the surfaces of glaciers and ice sheets during summer melt seasons, with larger blooms anticipated in years that experience the greatest melt. Here, we characterize the glacier algal bloom proliferating on Morteratsch glacier, Switzerland, during the record 2022 melt season, when the Swiss Alps lost three times more ice than the decadal average. Glacier algal cellular abundance (cells ml-1), biovolume (µm3 cell-1), photophysiology (Fv/Fm, rETRmax), and stoichiometry (C:N ratios) were constrained across three elevations on Morteratsch glacier during late August 2022 and compared with measurements of aqueous geochemistry and outputs of nutrient spiking experiments. While a substantial glacier algal bloom was apparent during summer 2022, abundances ranged from 1.78 × 104 to 8.95 × 105 cells ml-1 of meltwater and did not scale linearly with the magnitude of the 2022 melt season. Instead, spatiotemporal heterogeneity in algal distribution across Morteratsch glacier leads us to propose melt-water-redistribution of (larger) glacier algal cells down-glacier and presumptive export of cells from the system as an important mechanism to set overall bloom carrying capacity on steep valley glaciers during high melt years. Despite the paradox of abundant glacier algae within seemingly oligotrophic surface ice, we found no evidence for inorganic nutrient limitation as an important bottom-up control within our study site, supporting our hypothesis above. Fundamental physical constraints may thus cap bloom carrying-capacities on valley glaciers as 21st century melting continues.

Increasing surface runoff from Greenland's firn areas.

At high elevations of ice sheets, melting snow generally percolates and refreezes, so does not contribute to the shrinking of the ice sheet. Here, we systematically map the runoff area of the Greenland ice sheet, using surface rivers visible on satellite imagery. Between 1985 and 2020, the maximum runoff elevation rose by 58-329 metres, expanding the runoff area by 29% (-8%/+6%). Excess melt beyond the refreezing capacity of pores in snowfall has created near-impermeable ice slabs that sustain surface runoff even in cooler summers. We show that two surface mass balance models over-estimate the runoff area by 16-30%. Once restricted to our observed areas they indicate that 5-10% of recent runoff likely comes from the expanded runoff area. Runoff from higher elevations is sensitive to projected warming as further increases in the runoff limit will increase the runoff area disproportionately.

Mineral phosphorus drives glacier algal blooms on the Greenland Ice Sheet.

Melting of the Greenland Ice Sheet is a leading cause of land-ice mass loss and cryosphere-attributed sea level rise. Blooms of pigmented glacier ice algae lower ice albedo and accelerate surface melting in the ice sheet's southwest sector. Although glacier ice algae cause up to 13% of the surface melting in this region, the controls on bloom development remain poorly understood. Here we show a direct link between mineral phosphorus in surface ice and glacier ice algae biomass through the quantification of solid and fluid phase phosphorus reservoirs in surface habitats across the southwest ablation zone of the ice sheet. We demonstrate that nutrients from mineral dust likely drive glacier ice algal growth, and thereby identify mineral dust as a secondary control on ice sheet melting.

Greater Greenland Ice Sheet contribution to global sea level rise in CMIP6.

Future climate projections show a marked increase in Greenland Ice Sheet (GrIS) runoff during the 21st century, a direct consequence of the Polar Amplification signal. Regional climate models (RCMs) are a widely used tool to downscale ensembles of projections from global climate models (GCMs) to assess the impact of global warming on GrIS melt and sea level rise contribution. Initial results of the CMIP6 GCM model intercomparison project have revealed a greater 21st century temperature rise than in CMIP5 models. However, so far very little is known about the subsequent impacts on the future GrIS surface melt and therefore sea level rise contribution. Here, we show that the total GrIS sea level rise contribution from surface mass loss in our high-resolution (15 km) regional climate projections is 17.8  ±  7.8 cm in SSP585, 7.9 cm more than in our RCP8.5 simulations using CMIP5 input. We identify a +1.3 °C greater Arctic Amplification and associated cloud and sea ice feedbacks in the CMIP6 SSP585 scenario as the main drivers. Additionally, an assessment of the GrIS sea level contribution across all emission scenarios highlights, that the GrIS mass loss in CMIP6 is equivalent to a CMIP5 scenario with twice the global radiative forcing.

The silicon cycle impacted by past ice sheets.

Globally averaged riverine silicon (Si) concentrations and isotope composition (δ30Si) may be affected by the expansion and retreat of large ice sheets during glacial-interglacial cycles. Here we provide evidence of this based on the δ30Si composition of meltwater runoff from a Greenland Ice Sheet catchment. Glacier runoff has the lightest δ30Si measured in running waters (-0.25 ± 0.12‰), significantly lower than nonglacial rivers (1.25 ± 0.68‰), such that the overall decline in glacial runoff since the Last Glacial Maximum (LGM) may explain 0.06-0.17‰ of the observed ocean δ30Si rise (0.5-1.0‰). A marine sediment core proximal to Iceland provides further evidence for transient, low-δ30Si meltwater pulses during glacial termination. Diatom Si uptake during the LGM was likely similar to present day due to an expanded Si inventory, which raises the possibility of a feedback between ice sheet expansion, enhanced Si export to the ocean and reduced CO2 concentration in the atmosphere, because of the importance of diatoms in the biological carbon pump.

Decreasing cloud cover drives the recent mass loss on the Greenland Ice Sheet.

The Greenland Ice Sheet (GrIS) has been losing mass at an accelerating rate since the mid-1990s. This has been due to both increased ice discharge into the ocean and melting at the surface, with the latter being the dominant contribution. This change in state has been attributed to rising temperatures and a decrease in surface albedo. We show, using satellite data and climate model output, that the abrupt reduction in surface mass balance since about 1995 can be attributed largely to a coincident trend of decreasing summer cloud cover enhancing the melt-albedo feedback. Satellite observations show that, from 1995 to 2009, summer cloud cover decreased by 0.9 ± 0.3% per year. Model output indicates that the GrIS summer melt increases by 27 ± 13 gigatons (Gt) per percent reduction in summer cloud cover, principally because of the impact of increased shortwave radiation over the low albedo ablation zone. The observed reduction in cloud cover is strongly correlated with a state shift in the North Atlantic Oscillation promoting anticyclonic conditions in summer and suggests that the enhanced surface mass loss from the GrIS is driven by synoptic-scale changes in Arctic-wide atmospheric circulation.

Ice sheets as a missing source of silica to the polar oceans.

Ice sheets play a more important role in the global silicon cycle than previously appreciated. Input of dissolved and amorphous particulate silica into natural waters stimulates the growth of diatoms. Here we measure dissolved and amorphous silica in Greenland Ice Sheet meltwaters and icebergs, demonstrating the potential for high ice sheet export. Our dissolved and amorphous silica flux is 0.20 (0.06-0.79) Tmol year-1, ∼50% of the input from Arctic rivers. Amorphous silica comprises >95% of this flux and is highly soluble in sea water, as indicated by a significant increase in dissolved silica across a fjord salinity gradient. Retreating palaeo ice sheets were therefore likely responsible for high dissolved and amorphous silica fluxes into the ocean during the last deglaciation, reaching values of ∼5.5 Tmol year-1, similar to the estimated export from palaeo rivers. These elevated silica fluxes may explain high diatom productivity observed during the last glacial-interglacial period.

Ice sheets as a significant source of highly reactive nanoparticulate iron to the oceans.

The Greenland and Antarctic Ice Sheets cover ~ 10% of global land surface, but are rarely considered as active components of the global iron cycle. The ocean waters around both ice sheets harbour highly productive coastal ecosystems, many of which are iron limited. Measurements of iron concentrations in subglacial runoff from a large Greenland Ice Sheet catchment reveal the potential for globally significant export of labile iron fractions to the near-coastal euphotic zone. We estimate that the flux of bioavailable iron associated with glacial runoff is 0.40-2.54 Tg per year in Greenland and 0.06-0.17 Tg per year in Antarctica. Iron fluxes are dominated by a highly reactive and potentially bioavailable nanoparticulate suspended sediment fraction, similar to that identified in Antarctic icebergs. Estimates of labile iron fluxes in meltwater are comparable with aeolian dust fluxes to the oceans surrounding Greenland and Antarctica, and are similarly expected to increase in a warming climate with enhanced melting.