Contemporary ice sheet thinning drives subglacial groundwater exfiltration with potential feedbacks on glacier flow.

Observations indicate that groundwater-laden sedimentary aquifers are extensive beneath large portions of the Greenland and Antarctic ice sheets. A reduction in the mechanical loading of aquifers is known to lead to groundwater exfiltration, a discharge of groundwater from the aquifer. Here, we provide a simple expression predicting exfiltration rates under a thinning ice sheet. Using contemporary satellite altimetry observations, we predict that exfiltration rates may reach tens to hundreds of millimeters per year under the fastest thinning parts of the Antarctic Ice Sheet. In parts of West Antarctica, predicted rates of exfiltration would cause the total subglacial water discharge rate to be nearly double what is currently predicted from subglacial basal melting alone. Continued Antarctic Ice Sheet thinning into the future guarantees that the rate and potential importance of exfiltration will only continue to grow. Such an increase in warm, nutrient-laden subglacial water discharge would cause changes in ice sliding, melt of basal ice and marine biological communities.

Biogeochemical and historical drivers of microbial community composition and structure in sediments from Mercer Subglacial Lake, West Antarctica.

Ice streams that flow into Ross Ice Shelf are underlain by water-saturated sediments, a dynamic hydrological system, and subglacial lakes that intermittently discharge water downstream across grounding zones of West Antarctic Ice Sheet (WAIS). A 2.06 m composite sediment profile was recently recovered from Mercer Subglacial Lake, a 15 m deep water cavity beneath a 1087 m thick portion of the Mercer Ice Stream. We examined microbial abundances, used 16S rRNA gene amplicon sequencing to assess community structures, and characterized extracellular polymeric substances (EPS) associated with distinct lithologic units in the sediments. Bacterial and archaeal communities in the surficial sediments are more abundant and diverse, with significantly different compositions from those found deeper in the sediment column. The most abundant taxa are related to chemolithoautotrophs capable of oxidizing reduced nitrogen, sulfur, and iron compounds with oxygen, nitrate, or iron. Concentrations of dissolved methane and total organic carbon together with water content in the sediments are the strongest predictors of taxon and community composition. δ¹³C values for EPS (-25 to -30‰) are consistent with the primary source of carbon for biosynthesis originating from legacy marine organic matter. Comparison of communities to those in lake sediments under an adjacent ice stream (Whillans Subglacial Lake) and near its grounding zone provide seminal evidence for a subglacial metacommunity that is biogeochemically and evolutionarily linked through ice sheet dynamics and the transport of microbes, water, and sediments beneath WAIS.

A dynamic saline groundwater system mapped beneath an Antarctic ice stream.

Antarctica's fast-flowing ice streams drain the ice sheet, with their velocity modulated by subglacial water systems. Current knowledge of these water systems is limited to the shallow portions near the ice-bed interface, but hypothesized deeper groundwater could also influence ice streaming. Here, we use magnetotelluric and passive seismic data from Whillans Ice Stream, West Antarctica, to provide the first observations of deep sub-ice stream groundwater. Our data reveal a volume of groundwater within a >1-kilometer-thick sedimentary basin that is more than an order of magnitude larger than the known subglacial system. A vertical salinity gradient indicates exchange between paleo seawater at depth and contemporary basal meltwater above. Our results provide new constraints for subglacial water systems that affect ice streaming and subglacial biogeochemical processes.

Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves.

Ocean-driven basal melting of Antarctica's floating ice shelves accounts for about half of their mass loss in steady-state, where gains in ice shelf mass are balanced by losses. Ice shelf thickness changes driven by varying basal melt rates modulate mass loss from the grounded ice sheet and its contribution to sea level, and the changing meltwater fluxes influence climate processes in the Southern Ocean. Existing continent-wide melt rate datasets have no temporal variability, introducing uncertainties in sea level and climate projections. Here, we combine surface height data from satellite radar altimeters with satellite-derived ice velocities and a new model of firn-layer evolution to generate a high-resolution map of time-averaged (2010-2018) basal melt rates, and time series (1994-2018) of meltwater fluxes for most ice shelves. Total basal meltwater flux in 1994 (1090±150 Gt/yr) was not significantly different from the steady-state value (1100±60 Gt/yr), but increased to 1570±140 Gt/yr in 2009, followed by a decline to 1160±150 Gt/yr in 2018. For the four largest "cold-water" ice shelves we partition meltwater fluxes into deep and shallow sources to reveal distinct signatures of temporal variability, providing insights into climate forcing of basal melting and the impact of this melting on the Southern Ocean.

Pervasive ice sheet mass loss reflects competing ocean and atmosphere processes.

Quantifying changes in Earth's ice sheets and identifying the climate drivers are central to improving sea level projections. We provide unified estimates of grounded and floating ice mass change from 2003 to 2019 using NASA's Ice, Cloud and land Elevation Satellite (ICESat) and ICESat-2 satellite laser altimetry. Our data reveal patterns likely linked to competing climate processes: Ice loss from coastal Greenland (increased surface melt), Antarctic ice shelves (increased ocean melting), and Greenland and Antarctic outlet glaciers (dynamic response to ocean melting) was partially compensated by mass gains over ice sheet interiors (increased snow accumulation). Losses outpaced gains, with grounded-ice loss from Greenland (200 billion tonnes per year) and Antarctica (118 billion tonnes per year) contributing 14 millimeters to sea level. Mass lost from West Antarctica's ice shelves accounted for more than 30% of that region's total.

Multidecadal observations of the Antarctic ice sheet from restored analog radar records.

Airborne radar sounding can measure conditions within and beneath polar ice sheets. In Antarctica, most digital radar-sounding data have been collected in the last 2 decades, limiting our ability to understand processes that govern longer-term ice-sheet behavior. Here, we demonstrate how analog radar data collected over 40 y ago in Antarctica can be combined with modern records to quantify multidecadal changes. Specifically, we digitize over 400,000 line kilometers of exploratory Antarctic radar data originally recorded on 35-mm optical film between 1971 and 1979. We leverage the increased geometric and radiometric resolution of our digitization process to show how these data can be used to identify and investigate hydrologic, geologic, and topographic features beneath and within the ice sheet. To highlight their scientific potential, we compare the digitized data with contemporary radar measurements to reveal that the remnant eastern ice shelf of Thwaites Glacier in West Antarctica had thinned between 10 and 33% between 1978 and 2009. We also release the collection of scanned radargrams in their entirety in a persistent public archive along with updated geolocation data for a subset of the data that reduces the mean positioning error from 5 to 2.5 km. Together, these data represent a unique and renewed extensive, multidecadal historical baseline, critical for observing and modeling ice-sheet change on societally relevant timescales.

A decade of progress in observing and modelling Antarctic subglacial water systems.

In the decade since the discovery of active Antarctic subglacial water systems by detection of subtle surface displacements, much progress has been made in our understanding of these dynamic systems. Here, we present some of the key results of observations derived from ICESat laser altimetry, CryoSat-2 radar altimetry, Operation IceBridge airborne laser altimetry, satellite image differencing and ground-based continuous Global Positioning System (GPS) experiments deployed in hydrologically active regions. These observations provide us with an increased understanding of various lake systems in Antarctica: Whillans/Mercer Ice Streams, Crane Glacier, Recovery Ice Stream, Byrd Glacier and eastern Wilkes Land. In several cases, subglacial water systems are shown to control ice flux through the glacier system. For some lake systems, we have been able to construct more than a decade of continuous lake activity, revealing internal variability on time scales ranging from days to years. This variability indicates that continuous, accurate time series of altimetry data are critical to understanding these systems. On Whillans Ice Stream, our results from a 5-year continuous GPS record demonstrate that subglacial lake flood events significantly change the regional ice dynamics. We also show how models for subglacial water flow have evolved since the availability of observations of lake volume change, from regional-scale models of water routeing to process models of channels carved into the subglacial sediment instead of the overlying ice. We show that progress in understanding the processes governing lake drainage now allows us to create simulated lake volume time series that reproduce time series from satellite observations. This transformational decade in Antarctic subglacial water research has moved us significantly closer to understanding the processes of water transfer sufficiently for inclusion in continental-scale ice-sheet models.

Arsenic speciation in food chains from mid-Atlantic hydrothermal vents.

Arsenic concentration and speciation were determined in benthic fauna collected from the Mid-Atlantic Ridge hydrothermal vents. The shrimp species, Rimicaris exoculata, the vent chimney-dwelling mussel, Bathymodiolus azoricus, Branchipolynoe seepensis, a commensal worm of B. azoricus, and the gastropod Peltospira smaragdina showed variations in As concentration and in stable isotope (δ13C and δ15N) signature between species, suggesting different sources of As uptake. Arsenic speciation showed arsenobetaine to be the dominant species in R. exoculata, whereas in B. azoricus and B. seepensis arsenosugars were most abundant, although arsenobetaine, dimethylarsinate, and inorganic arsenic were also observed, along with several unidentified species. Scrape samples from outside the vent chimneys, covered with microbial mat, which is a presumed food source for many vent organisms, contained high levels of total As, but organic species were not detectable. The formation of arsenosugars in pelagic environments is typically attributed to marine algae, and the pathway to arsenobetaine is still unknown. The occurrence of arsenosugars and arsenobetaine in these deep sea organisms, where primary production is chemolithoautotrophic and stable isotope analyses indicate food sources are of vent origin, suggests that organic arsenicals can occur in a food web without algae or other photosynthetic life.

Elucidating the effect of additives on the growth and stability of Cu2O surfaces via shape transformation of pre-grown crystals.

A new strategy of using pre-grown crystals to study preferential adsorption of various additives is demonstrated for the electrocrystallization of Cu2O. In this method, micron-size Cu2O crystals with well-defined cubic and octahedral shapes were first electrochemically grown, and their crystallization was resumed in a medium containing the additive to be investigated (e.g., Na+, NH4+, SO42-, Cl-, dodecyl sulfate). This method makes it possible to systematically study the interaction of additives with specific planes (e.g., {100} of a cube and {111} of an octahedron) already present. By observing shape transformation over time, the relative stabilities of {100}, {111}, and {110} planes of Cu2O in various growth media could be determined. During this study, a general scheme of forming new crystal shapes containing crystallographic planes that cannot be directly stabilized by preferential adsorption alone was also established (i.e., rhombicuboctahedral shape of Cu2O containing {110} planes). This method can be extended to other crystal systems, which will enable us to classify common features of additives (e.g., charges, type of atoms) and crystallographic planes (e.g., atomic arrangement, surface termination, surface charge) required to allow for strong preferential adsorption.