Explore the vast range of titles available.

Finding suitable potential sites for an undisturbed record of million-year old ice in Antarctica requires slow-moving ice (preferably an ice divide) and basal conditions that are not disturbed by large topographic variations. Furthermore, ice should be thick and cold basal conditions should prevail, since basal melting would destroy the bottom layers. However, thick ice (needed to resolve the signal at sufficient high resolution) increases basal temperatures, which is a conflicting condition for finding a suitable drill site. In addition, slow moving areas in the center of ice sheets are also low-accumulation areas, and low accumulation reduces potential cooling of the ice through vertical advection. While boundary conditions such as ice thickness and accumulation rates are relatively well constrained, the major uncertainty in determining basal thermal conditions resides in the geothermal heat flow (GHF) underneath the ice sheet. We explore uncertainties in existing GHF data sets and their effect on basal temperatures of the Antarctic Ice Sheet, and propose an updated method based on Pattyn (2010) to improve existing GHF data sets in agreement with known basal temperatures and their gradients to reduce this uncertainty. Both complementary methods lead to a better comprehension of basal temperature sensitivity and a characterization of potential ice coring sites within these uncertainties. The combination of both modeling approaches show that the most likely oldest ice sites are situated near the divide areas (close to existing deep drilling sites, but in areas of smaller ice thickness) and across the Gamburtsev Subglacial Mountains.

Ice-shelf channels are long curvilinear tracts of thin ice found on Antarctic ice shelves. Many of them originate near the grounding line, but their formation mechanisms remain poorly understood. Here we use ice-penetrating radar data from Roi Baudouin Ice Shelf, East Antarctica, to infer that the morphology of several ice-shelf channels is seeded upstream of the grounding line by large basal obstacles indenting the ice from below. We interpret each obstacle as an esker ridge formed from sediments deposited by subglacial water conduits, and calculate that the eskers’ size grows towards the grounding line where deposition rates are maximum. Relict features on the shelf indicate that these linked systems of subglacial conduits and ice-shelf channels have been changing over the past few centuries. Because ice-shelf channels are loci where intense melting occurs to thin an ice shelf, these findings expose a novel link between subglacial drainage, sedimentation and ice-shelf stability. The formation mechanisms of ice-shelf channels remain poorly understood. Here, using ice-penetrating radar data, the authors propose that ice-shelf channel morphology in the Roi Baudouin Ice Shelf, East Antarctica, is seeded by esker ridges, indenting the ice from below.

Mangrove forests prominently occupy an intertidal boundary position where the effects of sea level rise will be fast and well visible. This study in East Africa (Gazi Bay, Kenya) addresses the question of whether mangroves can be resilient to a rise in sea level by focusing on their potential to migrate towards landward areas. The combinatory analysis between remote sensing, DGPS-based ground truth and digital terrain models (DTM) unveils how real vegetation assemblages can shift under different projected (minimum (+9 cm), relative (+20 cm), average (+48 cm) and maximum (+88 cm)) scenarios of sea level rise (SLR). Under SLR scenarios up to 48 cm by the year 2100, the landward extension remarkably implies an area increase for each of the dominant mangrove assemblages except for Avicennia marina and Ceriops tagal, both on the landward side. On the one hand, the increase in most species in the first three scenarios, including the socio-economically most important species in this area, Rhizophora mucronata and C. tagal on the seaward side, strongly depends on the colonisation rate of these species. On the other hand, a SLR scenario of +88 cm by the year 2100 indicates that the area flooded only by equinoctial tides strongly decreases due to the topographical settings at the edge of the inhabited area. Consequently, the landward Avicennia-dominated assemblages will further decrease as a formation if they fail to adapt to a more frequent inundation. The topography is site-specific; however non-invadable areas can be typical for many mangrove settings.

Surface melt has been tied to the collapse of Antarctic Peninsula ice shelves. This study illustrates that warmer temperatures associated with katabatic winds drive similar processes in an East Antarctic ice shelf, highlighting vulnerability to disintegration. Surface melt and subsequent firn air depletion can ultimately lead to disintegration of Antarctic ice shelves 1 , 2 causing grounded glaciers to accelerate 3 and sea level to rise. In the Antarctic Peninsula, foehn winds enhance melting near the grounding line 4 , which in the recent past has led to the disintegration of the most northerly ice shelves 5 , 6 . Here, we provide observational and model evidence that this process also occurs over an East Antarctic ice shelf, where meltwater-induced firn air depletion is found in the grounding zone. Unlike the Antarctic Peninsula, where foehn events originate from episodic interaction of the circumpolar westerlies with the topography, in coastal East Antarctica high temperatures are caused by persistent katabatic winds originating from the ice sheet’s interior. Katabatic winds warm and mix the air as it flows downward and cause widespread snow erosion, explaining > 3 K higher near-surface temperatures in summer and surface melt doubling in the grounding zone compared with its surroundings. Additionally, these winds expose blue ice and firn with lower surface albedo, further enhancing melt. The in situ observation of supraglacial flow and englacial storage of meltwater suggests that ice-shelf grounding zones in East Antarctica, like their Antarctic Peninsula counterparts, are vulnerable to hydrofracturing 7 .

Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving \"shelfy stream\" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients.

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

Sea level rise (SLR) is a long-lasting consequence of climate change because global anthropogenic warming takes centuries to millennia to equilibrate for the deep ocean and ice sheets. SLR projections based on climate models support policy analysis, risk assessment and adaptation planning today, despite their large uncertainties. The central range of the SLR distribution is estimated by process-based models. However, risk-averse practitioners often require information about plausible future conditions that lie in the tails of the SLR distribution, which are poorly defined by existing models. Here, a community effort combining scientists and practitioners builds on a framework of discussing physical evidence to quantify high-end global SLR for practitioners. The approach is complementary to the IPCC AR6 report and provides further physically plausible high-end scenarios. High-end estimates for the different SLR components are developed for two climate scenarios at two timescales. For global warming of +2°C in 2100 (RCP2.6/SSP1-2.6) relative to pre-industrial values our high-end global SLR estimates are up to 0.9 m in 2100 and 2.5 m in 2300. Similarly, for a (RCP8.5/SSP5-8.5), we estimate up to 1.6 m in 2100 and up to 10.4 m in 2300. The large and growing differences between the scenarios beyond 2100 emphasize the long-term benefits of mitigation. However, even a modest 2°C warming may cause multi-meter SLR on centennial time scales with profound consequences for coastal areas. Earlier high-end assessments focused on instability mechanisms in Antarctica, while here we emphasize the importance of the timing of ice shelf collapse around Antarctica. This is highly uncertain due to low understanding of the driving processes. Hence both process understanding and emission scenario control high-end SLR.

The magnitude of the Antarctic ice sheet's contribution to global sea-level rise is dominated by the potential of its marine sectors to become unstable and collapse as a response to ocean (and atmospheric) forcing. This paper presents Antarctic sea-level response to sudden atmospheric and oceanic forcings on multi-centennial timescales with the newly developed fast Elementary Thermomechanical Ice Sheet (f.ETISh) model. The f.ETISh model is a vertically integrated hybrid ice sheet–ice shelf model with vertically integrated thermomechanical coupling, making the model two-dimensional. Its marine boundary is represented by two different flux conditions, coherent with power-law basal sliding and Coulomb basal friction. The model has been compared to existing benchmarks. Modelled Antarctic ice sheet response to forcing is dominated by sub-ice shelf melt and the sensitivity is highly dependent on basal conditions at the grounding line. Coulomb friction in the grounding-line transition zone leads to significantly higher mass loss in both West and East Antarctica on centennial timescales, leading to 1.5 m sea-level rise after 500 years for a limited melt scenario of 10 m a−1 under freely floating ice shelves, up to 6 m for a 50 m a−1 scenario. The higher sensitivity is attributed to higher ice fluxes at the grounding line due to vanishing effective pressure. Removing the ice shelves altogether results in a disintegration of the West Antarctic ice sheet and (partially) marine basins in East Antarctica. After 500 years, this leads to a 5 m and a 16 m sea-level rise for the power-law basal sliding and Coulomb friction conditions at the grounding line, respectively. The latter value agrees with simulations by DeConto and Pollard (2016) over a similar period (but with different forcing and including processes of hydrofracturing and cliff failure). The chosen parametrizations make model results largely independent of spatial resolution so that f.ETISh can potentially be integrated in large-scale Earth system models.

Climate model projections are often aggregated into multi-model averages of all models participating in an intercomparison project, such as the Coupled Model Intercomparison Project (CMIP). The \"multi-model\" approach provides a sensitivity test to the models' structural choices and implicitly assumes that multiple models provide additional and more reliable information than a single model, with higher confidence being placed on results that are common to an ensemble. A first initiative of the ice sheet modeling community, SeaRISE, provided such multi-model average projections of polar ice sheets' contribution to sea-level rise. The SeaRISE Antarctic numerical experiments aggregated results from all models devoid of a priori selection, based on the capacity of such models to represent key ice-dynamical processes. Here, using the experimental setup proposed in SeaRISE, we demonstrate that correctly representing grounding line dynamics is essential to infer future Antarctic mass change. We further illustrate the significant impact on the ensemble mean and deviation of adding one model with a known bias in its ability of modeling grounding line dynamics. We show that this biased model can hardly be identified from the ensemble only based on its estimation of volume change, as ad hoc and untrustworthy parametrizations can force any modeled grounding line to retreat. However, tools are available to test parts of the response of marine ice sheet models to perturbations of climatic and/or oceanic origin (MISMIP, MISMIP3d). Based on recent projections of Pine Island Glacier mass loss, we further show that excluding ice sheet models that do not pass the MISMIP benchmarks decreases the mean contribution and standard deviation of the multi-model ensemble projection by an order of magnitude for that particular drainage basin.

We introduce PARASO, a novel five-component fully coupled regional climate model over an Antarctic circumpolar domain covering the full Southern Ocean. The state-of-the-art models used are the fast Elementary Thermomechanical Ice Sheet model (f.ETISh) v1.7 (ice sheet), the Nucleus for European Modelling of the Ocean (NEMO) v3.6 (ocean), the Louvain-la-Neuve sea-ice model (LIM) v3.6 (sea ice), the COnsortium for Small-scale MOdeling (COSMO) model v5.0 (atmosphere) and its CLimate Mode (CLM) v4.5 (land), which are here run at a horizontal resolution close to 1/4∘. One key feature of this tool resides in a novel two-way coupling interface for representing ocean–ice-sheet interactions, through explicitly resolved ice-shelf cavities. The impact of atmospheric processes on the Antarctic ice sheet is also conveyed through computed COSMO-CLM–f.ETISh surface mass exchange. In this technical paper, we briefly introduce each model's configuration and document the developments that were carried out in order to establish PARASO. The new offline-based NEMO–f.ETISh coupling interface is thoroughly described. Our developments also include a new surface tiling approach to combine open-ocean and sea-ice-covered cells within COSMO, which was required to make this model relevant in the context of coupled simulations in polar regions. We present results from a 2000–2001 coupled 2-year experiment. PARASO is numerically stable and fully operational. The 2-year simulation conducted without fine tuning of the model reproduced the main expected features, although remaining systematic biases provide perspectives for further adjustment and development.