Influence of fast ice on future ice shelf melting in the Totten Glacier area, East Antarctica

The Totten Glacier in East Antarctica is of major climatic interest because of the large fluctuations in its grounding line and potential vulnerability to climate change. Here, we use a series of high-resolution, regional NEMO-LIM-based (Nucleus for European Modelling of the Ocean coupled with the Louvain-la-Neuve sea ice model) experiments, which include an explicit treatment of ocean–ice shelf interactions, as well as a representation of grounded icebergs and fast ice, to investigate the changes in ocean–ice interactions in the Totten Glacier area between the recent past (1995–2014) and the end of the 21st century (2081–2100) under SSP4–4.5 climate change conditions. By the end of the 21st century, the wide areas of multiyear fast ice simulated in the recent past are replaced by small patches of first year fast ice along the coast, which decreases the total summer sea ice extent. The Antarctic Slope Current is accelerated by about 116 %, which decreases the heat exchange across the shelf and tends to reduce the ice shelf basal melt rate, but this effect is counterbalanced by the effect of the oceanic warming. As a consequence, despite the accelerated Antarctic Slope Current, the Totten ice shelf melt rate is increased by 91 % due to the intrusion of warmer water into its cavity. The representation of fast ice dampens the ice shelf melt rate increase throughout the 21st century, as the Totten ice shelf melt rate increase reaches 136 % when fast ice is not taken into account. The Moscow University ice shelf melt rate increase is even more impacted by the representation of fast ice, with a 36 % melt rate increase with fast ice, compared to a 75 % increase without a fast ice representation. This influence of the representation of fast ice in our simulations on the basal melting rate trend over the 21st century is explained by the large impact of the fast ice for present-day conditions (∼25 % difference in m yr−1), while the impact decreases significantly at the end of the 21st century (∼4 % difference in m yr−1). As a consequence, the reduction in the fast ice extent in the future induces a decrease in the fast ice effect on the ice shelf melt rate that partly compensates for the increase due to warming of the ocean. This highlights the importance of including a representation of fast ice to simulate realistic ice shelf melt rate increase in East Antarctica under warming conditions.</p

Influence of fast ice on future ice shelf melting in the Totten Glacier area, East Antarctica

The Totten Glacier in East Antarctica is of major climatic interest because of the large fluctuations in its grounding line and potential vulnerability to climate change. Here, we use a series of high-resolution, regional NEMO-LIM-based (Nucleus for European Modelling of the Ocean coupled with the Louvain-la-Neuve sea ice model) experiments, which include an explicit treatment of ocean–ice shelf interactions, as well as a representation of grounded icebergs and fast ice, to investigate the changes in ocean–ice interactions in the Totten Glacier area between the recent past (1995–2014) and the end of the 21st century (2081–2100) under SSP4–4.5 climate change conditions. By the end of the 21st century, the wide areas of multiyear fast ice simulated in the recent past are replaced by small patches of first year fast ice along the coast, which decreases the total summer sea ice extent. The Antarctic Slope Current is accelerated by about 116 %, which decreases the heat exchange across the shelf and tends to reduce the ice shelf basal melt rate, but this effect is counterbalanced by the effect of the oceanic warming. As a consequence, despite the accelerated Antarctic Slope Current, the Totten ice shelf melt rate is increased by 91 % due to the intrusion of warmer water into its cavity. The representation of fast ice dampens the ice shelf melt rate increase throughout the 21st century, as the Totten ice shelf melt rate increase reaches 136 % when fast ice is not taken into account. The Moscow University ice shelf melt rate increase is even more impacted by the representation of fast ice, with a 36 % melt rate increase with fast ice, compared to a 75 % increase without a fast ice representation. This influence of the representation of fast ice in our simulations on the basal melting rate trend over the 21st century is explained by the large impact of the fast ice for present-day conditions (∼25 % difference in m yr−1), while the impact decreases significantly at the end of the 21st century (∼4 % difference in m yr−1). As a consequence, the reduction in the fast ice extent in the future induces a decrease in the fast ice effect on the ice shelf melt rate that partly compensates for the increase due to warming of the ocean. This highlights the importance of including a representation of fast ice to simulate realistic ice shelf melt rate increase in East Antarctica under warming conditions.</p

Interactions between wind-blown snow redistribution and melt ponds in a coupled ocean–sea ice model

Introducing a parameterization of the interactions between wind-driven snow depth changes and melt pond evolution allows us to improve large scale models. In this paper we have implemented an explicit melt pond scheme and, for the first time, a wind dependant snow redistribution model and new snow thermophysics into a coupled ocean–sea ice model. The comparison of long-term mean statistics of melt pond fractions against observations demonstrates realistic melt pond cover on average over Arctic sea ice, but a clear underestimation of the pond coverage on the multi-year ice (MYI) of the western Arctic Ocean. The latter shortcoming originates from the concealing effect of persistent snow on forming ponds, impeding their growth. Analyzing a second simulation with intensified snow drift enables the identification of two distinct modes of sensitivity in the melt pond formation process. First, the larger proportion of wind-transported snow that is lost in leads directly curtails the late spring snow volume on sea ice and facilitates the early development of melt ponds on MYI. In contrast, a combination of higher air temperatures and thinner snow prior to the onset of melting sometimes make the snow cover switch to a regime where it melts entirely and rapidly. In the latter situation, seemingly more frequent on first-year ice (FYI), a smaller snow volume directly relates to a reduced melt pond cover. Notwithstanding, changes in snow and water accumulation on seasonal sea ice is naturally limited, which lessens the impacts of wind-blown snow redistribution on FYI, as compared to those on MYI. At the basin scale, the overall increased melt pond cover results in decreased ice volume via the ice-albedo feedback in summer, which is experienced almost exclusively by MYI

Impact of ice sheet meltwater fluxes on the climate evolution at the onset of the Last Interglacial

Large climate perturbations occurred during Termination II when the ice sheets retreated from their glacial configuration. Here we investigate the impact of ice sheet changes and associated freshwater fluxes on the climate evolution at the onset of the Last Interglacial. The period from 135 to 120 kyr BP is simulated with the Earth system model of intermediate complexity LOVECLIM v.1.3 with prescribed evolution of the Antarctic ice sheet, the Greenland ice sheet and the other Northern Hemisphere ice sheets. Variations in meltwater fluxes from the Northern Hemisphere ice sheets lead to North Atlantic temperature changes and modifications of the strength of the Atlantic meridional overturning circulation. By means of the interhemispheric see-saw effect, variations in the Atlantic meridional overturning circulation also give rise to temperature changes in the Southern Hemisphere, which are modulated by the direct impact of Antarctic meltwater fluxes into the Southern Ocean. Freshwater fluxes from the melting Antarctic ice sheet lead to a millennial time scale oceanic cold event in the Southern Ocean with expanded sea ice as evidenced in some ocean sediment cores, which may be used to constrain the timing of ice sheet retreat

PARASO, a circum-Antarctic fully coupled ice-sheet–ocean–sea-ice–atmosphere–land model involving f.ETISh1.7, NEMO3.6, LIM3.6, COSMO5.0 and CLM4.5

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.This research has been supported by the Fonds De La Recherche Scientifique – FNRS (grant no. O0100718F).Peer ReviewedArticle signat per 23 autors/es: Charles Pelletier (1), Thierry Fichefet (1), Hugues Goosse (1), Konstanze Haubner (2), Samuel Helsen (3), Pierre-Vincent Huot (1), Christoph Kittel (4), François Klein (1), Sébastien Le clec'h (5), Nicole P. M. van Lipzig (3), Sylvain Marchi (3), François Massonnet (1), Pierre Mathiot (6,7), Ehsan Moravveji (3,8), Eduardo Moreno-Chamarro (9), Pablo Ortega (9), Frank Pattyn (2), Niels Souverijns (3,10), Guillian Van Achter (1), Sam Vanden Broucke (3), Alexander Vanhulle (5), Deborah Verfaillie (1), and Lars Zipf (2) // (1) Earth and Life Institute (ELI), UCLouvain, Louvain-la-Neuve, Belgium / (2) Laboratoire de Glaciologie, Université Libre de Bruxelles, Brussels, Belgium / (3) Department of Earth and Environmental Sciences, KU Leuven, Leuven, Belgium / (4) Laboratory of Climatology, Department of Geography, SPHERES, University of Liège, Liège, Belgium / (5) Earth System Science and Departement Geografie, Vrije Universiteit Brussel, Brussels, Belgium, (6) Met Office, Exeter, United Kingdom / (7) Université Grenoble Alpes/CNRS/IRD/G-INP, IGE, Grenoble, France / (8) ICTS, KU Leuven, Leuven, Belgium / (9) Barcelona Supercomputing Center (BSC), Barcelona, Spain / (10) Environmental Modelling Unit, Flemish Institute for Technological Research (VITO), Mol, BelgiumPostprint (published version

Impact of model resolution on Arctic sea ice and North Atlantic Ocean heat transport

Arctic sea-ice area and volume have substantially decreased since the beginning of the satellite era. Concurrently, the pole-ward heat transport from the North Atlantic Ocean into the Arctic has increased, partly contributing to the loss of sea ice. Increasing the horizontal resolution of general circulation models (GCMs) improves their ability to represent the complex interplay of processes at high latitudes. Here, we investigate the impact of model resolution on Arctic sea ice and Atlantic Ocean heat transport (OHT) by using five different state-of-the-art coupled GCMs (12 model configurations in total) that include dynamic representations of the ocean, atmosphere and sea ice. The models participate in the High Resolution Model Intercomparison Project (HighResMIP) of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). Model results over the period 1950–2014 are compared to different observational datasets. In the models studied, a finer ocean resolution drives lower Arctic sea-ice area and volume and generally enhances Atlantic OHT. The representation of ocean surface characteristics, such as sea-surface temperature (SST) and velocity, is greatly improved by using a finer ocean reso-lution. This study highlights a clear anticorrelation at interannual time scales between Arctic sea ice (area and volume) and Atlantic OHT north of 60 ◦N in the models studied. However, the strength of this relationship is not systematically impacted by model resolution. The higher the latitude to compute OHT, the stronger the relationship between sea-ice area/volume and OHT. Sea ice in the Barents/Kara and Greenland–Iceland–Norwegian (GIN) Seas is more strongly connected to Atlantic OHT than other Arctic seas

Long-term climate change commitment and reversibility: an EMIC intercomparison

This paper summarizes the results of an intercomparison project with Earth System Models of Intermediate Complexity (EMICs) undertaken in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). The focus is on long-term climate projections designed to: (i) quantify the climate change commitment of different radiative forcing trajectories, and (ii) explore the extent to which climate change is reversible on human timescales. All commitment simulations follow the four Representative Concentration Pathways (RCPs) and their extensions to 2300. Most EMICs simulate substantial surface air temperature and thermosteric sea level rise commitment following stabilization of the atmospheric composition at year-2300 levels. The meridional overturning circulation (MOC) is weakened temporarily and recovers to near pre-industrial values in most models for RCPs 2.6–6.0. The MOC weakening is more persistent for RCP 8.5. Elimination of anthropogenic CO2 emissions after 2300 results in slowly decreasing atmospheric CO2 concentrations. At year 3000 atmospheric CO2 is still at more than half its year-2300 level in all EMICs for RCPs 4.5–8.5. Surface air temperature remains constant or decreases slightly and thermosteric sea level rise continues for centuries after elimination of CO2 emissions in all EMICs. Restoration of atmospheric CO2 from RCP to pre-industrial levels over 100–1000 years requires large artificial removal of CO2 from the atmosphere and does not result in the simultaneous return to pre-industrial climate conditions, as surface air temperature and sea level response exhibit a substantial time lag relative to atmospheric CO2

Solar radiation and global climate change: Some experiments with a two-dimensional climate model

n/