Evaluation of six geothermal heat flux maps for the Antarctic Lambert-Amery glacial system

Basal thermal conditions play an important role in ice sheet dynamics, and they are sensitive to geothermal heat flux (GHF). Here we estimate the basal thermal conditions, including basal temperature, basal melt rate, and friction heat underneath the Lambert-Amery glacier system in east Antarctica, using a combination of a forward model and an inversion from a 3D ice flow model. We assess the sensitivity and uncertainty of basal thermal conditions using six different GHFs. We evaluate the modelled results using all available observed subglacial lakes. There are very large differences in modelled spatial pattern of temperate basal conditions using the different GHFs. The two most-recent GHF fields inverted from aerial geomagnetic observations have higher values of GHF in the region, produce a larger warm-based area, and match the observed subglacial lakes better than the other GHFs. The fast flowing glacier region has a lower modelled basal friction coefficient, faster basal velocity, with higher basal frictional heating in the range of 50–2000 mW m−2 than the base under slower flowing glaciated areas. The modelled basal melt rate reaches ten to hundreds of mm per year locally in Lambert, Lepekhin and Kronshtadtskiy glaciers feeding the Amery ice shelf, and ranges from 0–5 mm yr−1 on the temperate base of the vast inland region.

Did Holocene climate changes drive West Antarctic grounding line retreat and readvance?

Knowledge of past ice sheet configurations is useful for informing projections of future ice sheet dynamics and for calibrating ice sheet models. The topology of grounding line retreat in the Ross Sea sector of Antarctica has been much debated, but it has generally been assumed that the modern ice sheet is as small as it has been for more than 100 000 years (Conway et al., 1999; Lee et al., 2017; Lowry et al., 2019; McKay et al., 2016; Scherer et al., 1998). Recent findings suggest that the West Antarctic Ice Sheet (WAIS) grounding line retreated beyond its current location earlier in the Holocene and subsequently readvanced to reach its modern position (Bradley et al., 2015; Kingslake et al., 2018). Here, we further constrain the post-LGM (Last Glacial Maximum) grounding line retreat and readvance in the Ross Sea sector using a two-phase model of radiocarbon input and decay in subglacial sediments from six sub-ice sampling locations. In addition, we reinterpret high basal temperature gradients, measured previously at three sites in this region (Engelhardt, 2004), which we explain as resulting from recent ice shelf re-grounding accompanying grounding line readvance. At one location – Whillans Subglacial Lake (SLW) – for which a sediment porewater chemistry profile is known, we estimate the grounding line readvance by simulating ionic diffusion. Collectively, our analyses indicate that the grounding line retreated over SLW 4300-2500+1500 years ago, and over sites on Whillans Ice Stream (WIS), Kamb Ice Stream (KIS), and Bindschadler Ice Stream (BIS) 4700-2300+1500, 1800-700+2700, and 1700-600+2800 years ago, respectively. The grounding line only recently readvanced back over those sites 1100-100+200, 1500-200+500, 1000-300+200, and 800±100 years ago for SLW, WIS, KIS, and BIS, respectively. The timing of grounding line retreat coincided with a warm period in the mid-Holocene to late Holocene. Conversely, grounding line readvance is coincident with cooling climate in the last 1000–2000 years. Our estimates for the timing of grounding line retreat and readvance are also consistent with relatively low carbon-to-nitrogen ratios measured in our subglacial sediment samples (suggesting a marine source of organic matter) and with the lack of grounding zone wedges in front of modern grounding lines. Based on these results, we propose that the Siple Coast grounding line motions in the mid-Holocene to late Holocene were primarily driven by relatively modest changes in regional climate, rather than by ice sheet dynamics and glacioisostatic rebound, as hypothesized previously (Kingslake et al., 2018).

Fjord network in Namibia: A snapshot into the dynamics of the late Paleozoic glaciation

Fjords are glacially carved estuaries that profoundly influence ice-sheet stability by draining and ablating ice. Although abundant on modern high-latitude continental shelves, fjord-network morphologies have never been identified in Earth’s pre-Cenozoic glacial epochs, hindering our ability to constrain ancient ice-sheet dynamics. We show that U-shaped valleys in northwestern Namibia cut during the late Paleozoic ice age (LPIA, ca. 300 Ma), Earth’s penultimate icehouse, represent intact fjord-network morphologies. This preserved glacial morphology and its sedimentary fill permit a reconstruction of paleo-ice thicknesses, glacial dynamics, and resulting glacio-isostatic adjustment. Glaciation in this region was initially characterized by an acme phase, which saw an extensive ice sheet (1.7 km thick) covering the region, followed by a waning phase characterized by 100-m-thick, topographically constrained outlet glaciers that shrank, leading to glacial demise. Our findings demonstrate that both a large ice sheet and highland glaciers existed over northwestern Namibia at different times during the LPIA. The fjords likely played a pivotal role in glacier dynamics and climate regulation, serving as hotspots for organic carbon sequestration. Aside from the present-day arid climate, northwestern Namibia exhibits a geomorphology virtually unchanged since the LPIA, permitting unique insight into this icehouse.

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ∼2 km, as opposed to ∼0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ∼2 km (arguably 4 km) when the connectivity is low or moderate and ∼1 km (arguably 2 km) when the connectivity is strong.

Inversion of Geothermal Heat Flux under the Ice Sheet of Princess Elizabeth Land, East Antarctica

Antarctic geothermal heat flux is a basic input variable for ice sheet dynamics simulation. It greatly affects the temperature and mechanical properties at the bottom of the ice sheet, influencing sliding, melting, and internal deformation. Due to the fact that the Antarctica is covered by a thick ice sheet, direct measurements of heat flux are very limited. This study was carried out to estimate the regional heat flux in the Antarctic continent through geophysical inversion. Princess Elizabeth Land, East Antarctica is one of the areas in which we have a weak understanding of geothermal heat flux. Through the latest airborne geomagnetic data, we inverted the Curie depth, obtaining the heat flux of bedrock based on the one-dimensional steady-state heat conduction equation. The results indicated that the Curie depth of the Princess Elizabeth Land is shallower than previously estimated, and the heat flux is consequently higher. Thus, the contribution of subglacial heat flux to the melting at the bottom of the ice sheet is likely greater than previously expected in this region. It further provides research clues for the formation of the developed subglacial water system in Princess Elizabeth Land.

Red Sea tectonics unveil the largest known terrestrial ice stream: New constraints on Late Ordovician ice sheet dynamics

Mega-streamlined landforms on Earth and Mars have been attributed to aeolian, glaciogenic, fluvial, and tectonic processes. Identifying the forces that shaped these landforms is paramount for understanding landscape evolution and constraining paleo-climate models and ice sheet reconstructions. In Arabia, east-northeast, kilometer-scale streamlined landforms were interpreted to have been formed by Quaternary aeolian erosion. We provide field and satellite-based evidence for a Late Ordovician glacial origin for these streamlined landforms, which were exhumed during the Red Sea–related uplift. Then we use Late Ordovician paleo-topographic data to reconstruct the Late Ordovician ice sheet using identified and previously reported glacial deposits and landforms. Our reconstruction suggests these glacial features are part of a major, topographically controlled, marine-terminating ice stream, twice the length of the largest known terrestrial ice streams. Our results support models that advocate for a single, major, and highly dynamic ice sheet and provide new morphological-based constraints for Late Ordovician climate models.