Search

The Greenland ice sheet (GrIS) is currently losing mass at an accelerated rate, due to atmospheric and ocean warming causing respectively enhanced melt and ice discharge to the ocean. A large part of the uncertainty on future GrIS contribution to sea level rise relates to unknown atmospheric and ocean circulation change. For the later, AR6 models project a weakening of the North Atlantic Meridional Overturning Circulation (NAMOC) during the 21st century. The magnitude of this weakening depends on the greenhouse gas scenario and model, but none of the models project a complete collapse.Projections of future GrIS evolution in the last IPCC report AR6 are mostly based on simulations with ice sheet models forced with the output of climate models (e.g., Goelzer et al. (2020)). This method permits large ensembles of simulations, however the coupling between climate and GrIS is not represented. Here, we use a coupled Earth System and Ice Sheet Model (ESM-ISM), the CESM2-CISM2 (Muntjewerf et al. 2021) to examine the multi-millennial evolution of the GrIS surface mass balance for a middle-of-the-road CO2 scenario. The model couples realistic simulation of global climate (Danabasoglu et al. 2020), surface processes (van Kampenhout et al. 2020) and ice dynamics (Lipscomb et al. 2019). We use an idealized scenario of 1% CO2 increase until stabilization at two times pre-industrial values. compare our results with pre-industrial and 1% to 4xCO2 simulations (Muntjewerf et al. 2020).We find small increases and even reduction of annual temperatures in the GrIS area in connection with strong NAMOC weakening in the first two centuries of simulation. Summer temperatures and surface melt increase moderately with respect to pre-industrial. From simulation year 500, the NAMOC recovers, resulting in strong increases in GrIS melt rates and contribution to sea level rise. We compare the deglaciation pattern over a period of 3,000 years with deglaciation simulations with the same model for the last interglacial (Sommers et al. 2021).

Land ice discharge from the Antarctic continent into the ocean is restrained by ice shelves, floating extensions of grounded ice that buttress the glacier outflow. The ongoing thinning of these ice shelves – largely due to enhanced melting at their base in response to global warming – is known to accelerate the release of glacier meltwater into the world oceans, augmenting global sea level. Mechanisms of ocean heat intrusion under the ice base are therefore crucial to project the future of Antarctic ice shelves. Furthermore, ice shelves are weakened by the presence of kilometer-wide full-thickness ice rifts, which are observed all around Antarctica. However, their impact on ocean circulation around and below ice shelves has been largely unexplored as ocean models are commonly characterized by resolutions that are too coarse to resolve their presence. Here, we apply the Massachusetts Institute of Technology general circulation model at high resolution to investigate the sensitivity of sub-shelf ocean dynamics and ice-shelf melting to the presence of a kilometer-wide rift in proximity of the ice front. We find that (a) the rift curtails water and heat intrusion beneath the ice-shelf base and (b) the basal melting of a rifted ice shelf is on average 20 % lower than for an intact ice shelf under identical forcing. Notably, we calculate a significant reduction in melting rates of up to 30 % near the grounding line of a rifted ice shelf. We therefore posit that rifts and their impact on the sub-shelf dynamics are important to consider in order to accurately reproduce and project pathways of heat intrusion into the ice-shelf cavity.

The Arctic is the region on Earth that is warming the fastest. At the same time, Arctic sea ice is reducing while the Greenland ice sheet (GrIS) is losing mass at an accelerated pace. Here, we study the seasonal impact of reduced Arctic sea ice on GrIS surface mass balance (SMB), using the Community Earth System Model version 2.1 (CESM2), which features an advanced, interactive calculation of SMB. Addressing the impact of sea-ice reductions on the GrIS SMB from observations is difficult due to the short observational records. Also, signals detected using transient climate simulations may be aliases of other forcings. Here, we analyze dedicated simulations from the Polar Amplification Model Intercomparison Project with reduced Arctic sea ice and compare them with preindustrial sea ice simulations while keeping all other forcings constant. In response to reduced sea ice, the GrIS SMB increases in winter due to increased precipitation, driven by the more humid atmosphere and increasing cyclones. In summer, surface melt increases due to a warmer, more humid atmosphere providing increased energy transfer to the surface through the sensible and latent heat fluxes, which triggers the melt-albedo feedback. Further, warming occurs throughout the entire troposphere over Baffin Bay. This deep warming results in regional enhancement of the 500 hPa geopotential heights over the Baffin Bay and Greenland, increasing blocking and heat advection over the GrIS’ surface. This anomalous circulation pattern has been linked to recent increases in the surface melt of the GrIS. Supplementary Information The online version contains supplementary material available at 10.1007/s00382-021-05897-4.

The elevation feedback on melt has been identified as a key process to explain (Gregoire et al, Nature, 2012) and project (Aschwanden et al., Sci Advances, 2019; Ridley et al, J Clim, 2005; Vizcaino et al, Clim Dyn, 2008) long-term deglaciation. It is also central to the theory of ice sheet evolution hysteresis, deglaciation thresholds/tipping points, and the problem of reversibility (Garbe et al, Nature, 2020; Gregory et al, TC, 2020; Gregory et al, Nature, 2004, Robinson et al, Nature Clim. Change, 2012). Ice-sheet-model-only estimates of this feedback rely on a largely unexplored parameter, the so-called “lapse rate”. This parameter is defined as the rate at which the near-surface atmospheric temperature changes strictly due to surface elevation change.In this work, we use coupled an uncoupled ice sheet and climate simulations with two different General Circulation/Earth System Models (Vizcaino et al, GRL, 2015; Muntjewerf et al, JAMES, 2019) to estimate the temperature lapse rate over the Greenland ice sheet as it deglaciates. We find that this lapse rate is highly variable over seasons, with much reduced lapse rates during summer over melting surfaces. We propose that uncoupled state-of-the-art projections are likely overestimating deglaciation rates due to too high summer lapse rates over the ablation area.

Previous coupled climate-ice sheet modeling studies indicate that the warming threshold leading to multi-millennial, large-scale deglaciation of the Greenland Ice Sheet (GrIS) is in the range of 1.6-3.0 K above the pre-industrial climate. These studies either used an intermediate complexity RCM (Robinson et al. 2012) or a low resolution GCM (Gregory et al., 2020) coupled to a zero-order ISM. Here, we investigate the warming threshold and long-term response time of the GrIS using the higher-order Community Ice Sheet Model version 2 (CISM2, Lipscomb et al. 2019), forced with surface mass balance (SMB) calculated with the Community Earth System Model version 2 (CESM2, Danabasoglu et al. 2020). We use different forcing climatologies from a coupled CESM2/CISM2 simulation under high greenhouse gas forcing (Muntjewerf et al. 2020), where each climatology corresponds to a different global warming level in the range of 1-8.5 K above the pre-industrial climate. The SMB, which is calculated in CESM2 using an advanced energy balance scheme at multiple elevation classes (Muntjewerf et al. 2020), is downscaled during runtime to CISM2, thus allowing to account for the surface elevation feedback. In all the simulations the forcing is cycled until the ice sheet is fully deglaciated or has reached a new equilibrium. In a first set of simulations, we find that for a warming level higher than 5.2 K above pre-industrial the ice sheet will disappear, with the timing ranging between 2000 (+8.5 K) and 6000 years (+5.2 K). At a warming level of 2.8 K above pre-industrial, the ice loss does not exceed 2 m SLE, and most of the retreat occurs in the first 10,000 years in the south-west and central-west basins. In contrast, with a higher warming level of 3.6 K above pre-industrial as much as 7 m SLE of ice are loss in 20,000 years, with primary contributions from the western, northern and north-eastern basins. We will conclude by showing preliminary results from a second set of simulations focusing on the 2.8-3.6 K warming above pre-industrial interval.

Accelerated surface melt of the Greenland ice sheet (GrIS) is currently a large contributor to sea level rise, and the primary process of GrIS mass loss. Projections of future GrIS melt are limited by the lack of explicit melt calculations within most global climate models and the high computational cost of dynamical downscaling with regional models. To translate global climate evolution to GrIS surface melt, we train artificial neural networks (ANNs) with the output of the explicit melt calculation of the Community Earth System Model 2.1 (CESM2). ANNs are well suited for this task, as they are capable of learning complex, non-linear relationships, and they are fast to run.Our results show that the ANNs accurately project GrIS surface melt when evaluated against regional climate simulations. Further, the ANNs recognize patterns already established in litterature as important for surface melt, and use bases the projections on these patterns. Using the global climate simulations from the CMIP6 archive, the ANNs project a GrIS surface melt increase ranging from 414 Gt yr-1 to 1,378 Gt yr-1 by the end of the 21st century. The main source of projection uncertainty throughout the 21st century is due to the spread in the models’ climate sensitivity.

The Greenland ice sheet is currently losing mass at a rate of 0.8 mm of global mean sea level rise (SLR) per year. Here, we simulate its future evolution under an idealized scenario of high greenhouse gas forcing (1% increase per year until four times pre-industrial CO2). To this end, we use the newly, bi-directionally coupled Community Earth System Model version 2 and Community Ice Sheet Model version 2 (CESM2-CISM2, Muntjewerf et al, GRL, 2020), that includes an advanced calculation of the surface mass balance in the land component with elevation classes downscaling to CISM. Deglaciation rates increase from 2 mm SLR per year by simulation year 140 (or time of CO2 stabilization) to 7 mm SLR per year two centuries later as the ablation areas expand and net solar radiation and turbulent (latent, sensible) heat fluxes become the dominant energy sources for melt. The ice sheet retreats to an ice cap in the interior of the northern half of Greenland, that melts completely by simulation year 1,700. We compare the Greenland climate evolution with a CESM2 simulation with fixed topography, and evaluate the role of vegetation, clouds, precipitation, and surface energy fluxes on the relatively fast decay of the ice sheet. In addition, we use a set of CISM2 simulations forced with CESM2 SMB to estimate the global warming/forcing threshold for complete deglaciation.

The last deglaciation provides as unique a framework to investigate the processes of ice sheet and climate interaction during periods of mass loss as in the current climate. Here we simulate the Last Glacial Maximum (LGM) northern hemisphere ice sheets climate, surface mass balance (SMB), and dynamics with the Community Earth System Model version 2 (CESM2, Danabasoglu et al., 2020)) and the Community Ice Sheet Model version 2 (CISM2, Lipscomb et al., 2019). This LGM simulation will be later used as starting point for coupled CESM2-CISM2 simulations of the last deglaciation. CESM2 is run at the nominal resolution used for IPCC-type projections (approx. 1 degree for all components). The model includes an advanced snow/firn and SMB calculation (van Kampenhout et al, 2019; Sellevold et al, 2019) the land component (CLM, cite) that has been evaluated and applied to the simulation of the future Greenland melt (van Kampenhout et al, 2020, Muntjewerf et al., 2020a,b, Sellevold & Vizcaino, 2020). Our analysis examines how the global, Arctic, and North Atlantic climate result in the simulated radiative and turbulent heat fluxes over the ice sheets, and the mass fluxes from precipitation, refreezing, runoff, and sublimation. We also examine the simulated ice streams in CISM2, which is run at 8 km under a higher-order approximation for ice flow.

Hydrological regimes of alpine catchments are expected to be strongly affected by climate change mostly due to their dependence on snow and ice dynamics. While seasonal changes have been studied extensively, studies on changes in the timing and magnitude of annual extremes remain rare. This study investigates the effects of climate change on runoff patterns in six contrasting alpine catchments in Austria using a process-based semi-distributed hydrological model and projections from 14 regional climate and global climate model combinations for RCP 4.5 and RCP 8.5. The study catchments represent a spectrum of different hydrological regimes, from pluvial-nival to nivo-glacial, as well as distinct topographies and land forms, characterizing different elevation zones across the Eastern Alps to provide a comprehensive picture of future runoff changes. The climate projections are used to model river runoff in 2071–2100, which are then compared to the 1981–2010 reference period for all study catchments. Changes in timing and magnitude of annual maximum and minimum flows as well as in monthly runoff and snow melt are quantified and analyzed. Our results indicate a substantial shift to earlier occurrences in annual maximum flows by 9 to 31 days and an extension of the potential flood season by one to three months for high-elevation catchments. For low-elevation catchments, changes in timing of annual maximum flows are less pronounced. Magnitudes of annual maximum flows are likely to increase by 2–18 % under RCP 4.5, while no clear changes are projected for four catchments under RCP 8.5. The latter is caused by a pronounced increase in evaporation and decrease in snow melt contributions which offset increases in precipitation. Minimum annual runoff occur 13–31 days earlier in the winter months for high-elevation catchments, whereas for low-elevation catchments a shift from winter to autumn by about 15–100 days is projected. While all catchments show an increase in mean magnitude of minimum flows by 7–30 % under RCP 4.5, this is only the case for four catchments under RCP 8.5. Our results suggest a relationship between the elevation of catchments and changes in timing of annual maximum and minimum flows. For the magnitude of the extreme flows, a relationship is found between catchment elevation and annual minimum flows, whereas this relationship is lacking between elevation and annual maximum flow.

Future Greenland ice sheet (GrIS) melt projections are limited by the lack of explicit melt calculations within most global climate models and the high computational cost of dynamical downscaling with regional climate models (RCMs). Here, we train artificial neural networks (ANNs) to obtain relationships between quantities consistently available from global climate model simulations and annually integrated GrIS surface melt. To this end, we train the ANNs with model output from the Community Earth System Model 2.1 (CESM2), which features an interactive surface melt calculation based on a downscaled surface energy balance. We find that ANNs compare well with an independent CESM2 simulation and RCM simulations forced by a CMIP6 subset. The ANNs estimate a melt increase for 2,081–2,100 ranging from 414 (Formula presented.) 275 Gt (Formula presented.) (SSP1-2.6) to 1,378 (Formula presented.) 555 Gt (Formula presented.) (SSP5-8.5) for the full CMIP6 suite. The primary source of uncertainty throughout the 21st century is the spread of climate model sensitivity.

Hydrological regimes of alpine catchments are expected to be strongly affected by climate change, mostly due to their dependence on snow and ice dynamics. While seasonal changes have been studied extensively, studies on changes in the timing and magnitude of annual extremes remain rare. This study investigates the effects of climate change on runoff patterns in six contrasting Alpine catchments in Austria using a process-based, semi-distributed hydrological model and projections from 14 regional and global climate model combinations for two representative concentration pathways, namely RCP4.5 and RCP8.5. The study catchments represent a spectrum of different hydrological regimes, from pluvial–nival to nivo-glacial, as well as distinct topographies and land forms, characterizing different elevation zones across the eastern Alps to provide a comprehensive picture of future runoff changes. The climate projections are used to model river runoff in 2071–2100, which are then compared to the 1981–2010 reference period for all study catchments. Changes in the timing and magnitude of annual maximum and minimum flows, as well as in monthly runoff and snowmelt, are quantified and analyzed. Our results indicate a substantial shift to earlier occurrences in annual maximum flows by 9 to 31 d and an extension of the potential flood season by 1 to 3 months for high-elevation catchments. For low-elevation catchments, changes in the timing of annual maximum flows are less pronounced. Magnitudes of annual maximum flows are likely to increase by 2 %–18 % under RCP4.5, while no clear changes are projected for four catchments under RCP8.5. The latter is caused by a pronounced increase in evaporation and decrease in snowmelt contributions, which offset increases in precipitation. In the future, minimum annual runoff will occur 13–31 d earlier in the winter months for high-elevation catchments, whereas for low-elevation catchments a shift from winter to autumn by about 15–100 d is projected, with generally larger changes for RCP8.5. While all catchments show an increase in mean magnitude of minimum flows by 7–30% under RCP4.5, this is only the case for four catchments under RCP8.5. Our results suggest a relationship between the elevation of catchments and changes in the timing of annual maximum and minimum flows. For the magnitude of the extreme flows, a relationship is found between catchment elevation and annual minimum flows, whereas this relationship is lacking between elevation and annual maximum flow., info:eu-repo/semantics/published

The Greenland ice sheet (GrIS) has been losing mass in the last several decades, with a current contributing of around 0.7 mm per year to global mean sea level rise (SLR). Projections of future melt rates are often derived from standalone ice sheet models, forced by data from global or regional climate models. In many cases, the surface mass balance parameterization relies on simplified schemes that relate melt with surface temperature.In this study, we present a mass and energy conserving, 350-year simulation with the Community Earth System Model version 2.1 (CESM2.1) bidirectionally coupled to the Community Ice Sheet Model version 2.1 (CISM2.1). In this simulation, the carbon dioxide concentration is initially increasing by 1% per year from pre-industrial levels (287 ppmv), to a quadrupling (1140 ppmv) and stabilization after year 140. The model simulates a global warming of 5.3 K and 8.5 K with respect to preindustrial by years 131-150 and 331-150, respectively, and a strong decline in the North Atlantic Meridional Overturning Circulation that is initiated before GrIS runoff substantially increases. 91% of the total GrIS contribution to global mean sea level rise (SLR, 1140 mm) is simulated in the two centuries following CO2 stabilization, as the mass loss increases from 2.2 mm SLR per year in 131-150 to 6.6 mm SLR per year in 331-351. This increase is caused by melt acceleration as the ablation areas expand, and Greenland summer surface temperatures predominantly approach melt conditions when the global warming exceeds a certain threshold (around 4.2 K). This enhances the albedo and turbulent heat fluxes contribution to total melt energy.

With the Community Earth System Model version 2.1 (CESM2.1) interactively coupled to the evolving Greenland Ice Sheet as simulated by the Community Ice Sheet Model version 2.1 (CISM2.1), we examine the Greenland Ice Sheet (GrIS) mass balance. The model has been run for the period 1850-2100 with historical and SSP5-8.5 scenario forcing, contributing to the coupled experiments within the framework of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) (Nowicki et al., 2016).CESM2.1-CISM2.1 simulates a relatively strong global warming signal and strong weakening of meridional overturning circulation by 2100 compared to CMIP5 models. In our projection, the GrIS contributes 23 mm sea-level equivalent by 2050, and 109 mm by 2100, to global mean sea level rise. The southern GrIS drainage basins contribute 73% of the mass loss by mid-century, but the contribution decreases to 55% by 2100, as surface runoff in the northern basins progressively increases.

We use the Community Earth System model 2.1 to investigate the response of the Greenland Ice sheet (GrIS) surface mass balance (SMB) to an idealized high CO2 forcing scenario (1% per year increase to four-times-preindustrial). The SMB calculation is coupled with the atmospheric model, using a physically-based surface energy balance scheme for melt, explicit calculation of snow albedo, and a realistic treatment of polar snow and firn compaction. The SMB becomes negative for a global mean temperature increase of 2.7 K compared to pre-industrial temperature, and the surface mass loss accelerates. Longwave radiation is the primary contributor to melt energy before acceleration. A decrease of the albedo due to ablation area expansion together with turbulent heat flux increase due to the surface of the ice sheet nearing melting point, are the main contributors at/after acceleration. Further, trends towards more positive North Atlantic Oscillation and more negative Greenland Blocking Index partially reduces future melt increase.

Coupled climate models predict a weakening of the Atlantic Meridional Overturning (AMOC) circulation in the future. However, it is not clear what is the cause of the AMOC weakening. Studies have suggested that the freshwater (FW) is an important factor in the AMOC reduction. There are different sources of FW that may play a role, such as, river discharge, sea ice melt, and precipitation. Currently, due to global warming, the Greenland Ice Sheet (GrIS) melt rate is rising, which increases the amount of freshwater (ice discharge) into the ocean. Thus, it is possible that this input of freshwater would affect the ocean circulation on a regional and global scale. Hence, the GrIS freshwater cannot be neglected. The goal of this study is to understand the impact of the freshwater from GrIS on the North AMOC (NAMOC) strength in the future. We used the Community Earth System Model (CESM) version 2.1, which contains a fully coupled and an active ice sheet, to simulate an idealized greenhouse gas scenario (1% CO2). The CO2 concentration is 1140 ppm at the end of the simulation. The results show that GrIS delivers, on average, about 0.062 Sv/yr of FW to the Subpolar North Atlantic Ocean. The bulk of the total freshwater input comes from the southeastern and southwestern parts of the ice sheet: the regions where some fast-flowing marine-terminating glaciers are located (e.g. Helheim and Kangerlussuaq). The NAMOC index (maximum barotropic stream function from above 28°N and from 500 m to 5500 m depth) was calculated. It displays a fast weakening, approximately 16.7 Sv (0.11 Sv/yr), during the first 150 yrs. After that, the NAMOC reaches a stable state where the index is around 5.7 Sv (year 350). When the NAMOC index was compared to the FW from GrIS time series, we observed that change in AMOC occurs before the FW starts to increase (from year 200). Our results thus suggest that the FW input from GrIS does not cause significant changes in the AMOC strength. It is necessary to further investigate other possible causes for the strong NAMOC decline in this model.

Recent research shows increasing decadal ice mass losses from the Greenland and Antarctic Ice Sheets and more generally from glaciers worldwide in the light of continued global warming. Here, in an update of our previous ISMASS paper (Hanna et al. 2013), we review recent observational estimates of ice sheet and glacier mass balance, and their related uncertainties, first briefly considering relevant monitoring methods. Focusing on the response to climate change during 1992–2018, and especially the post-IPCC AR5 period, we discuss recent changes in the relative contributions of ice sheets and glaciers to sea-level change. We assess recent advances in understanding of the relative importance of surface mass balance and ice dynamics in overall ice-sheet mass change. We also consider recent improvements in ice-sheet modelling, highlighting data-model linkages and the use of updated observational datasets in ice-sheet models. Finally, by identifying key deficiencies in the observations and models that hamper current understanding and limit reliability of future ice-sheet projections, we make recommendations to the research community for reducing these knowledge gaps. Our synthesis aims to provide a critical and timely review of the current state of the science in advance of the next Intergovernmental Panel on Climate Change Assessment Report that is due in 2021., SCOPUS: re.j, info:eu-repo/semantics/published

The Greenland Ice Sheet (GrIS) mass balance is examined with an Earth system/ice sheet model that interactively couples the GrIS to the broader Earth system. The simulation runs from 1850 to 2100, with historical and SSP5-8.5 forcing. By the mid-21st century, the cumulative GrIS contribution to global mean sea level rise (SLR) is 23 mm. During the second half of the 21st century, the surface mass balance becomes negative in all drainage basins, with an additional SLR contribution of 86 mm. The annual mean GrIS mass loss in the last two decades is 2.7-mm sea level equivalent (SLE) year−1. The increased SLR contribution from the surface mass balance (3.1 mm SLE year−1) is partly offset by reduced ice discharge from thinning and retreat of outlet glaciers. The southern GrIS drainage basins contribute 73% of the mass loss in mid-century but 55% by 2100, as surface runoff increases in the northern basins.

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin‐up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet‐enabled simulations with CESM2‐CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6)., Key Points We describe a computationally tractable, iterative procedure for spinning up a coupled Earth system‐ice sheet modelEquilibrium state from the iterative procedure is similar to a more expensive traditional model spin‐up with prescribed ice sheetsThe procedure is used for developing initial conditions for transient, fully coupled simulations in the Coupled Model Intercomparison Project phase 6

The Greenland Ice Sheet (GrIS) is currently losing ice mass as the result of changes in the complex ice-climate interactions that have been driven by global climate change. In order to accurately predict future sea level rise, the mechanisms driving the observed mass loss must be better understood. Here, we combine data from the satellite gravimetry mission GRACE, surface mass balance (SMB) output of RACMO 2.3, and ice discharge estimates to analyze the mass budget of Greenland at various temporal and spatial scales. Firstly, in agreement with previous estimates, we find that the rate of mass loss from Greenland observed by GRACE was between −277 and −269 Gt/yr in 2003–2012. This estimate is consistent with the sum of individual contributions: surface mass balance (SMB, around 216 ± 122 Gt/yr) and ice discharge (520 ± 31 Gt/yr), indicating a good performance of the regional climate model. Secondly, we examine the average accelerations of mass anomalies in Greenland over 2003–2012, suggesting that the SMB (−23.3 ± 2.7 Gt/yr2) contributes 75 % to the total acceleration observed by GRACE. The remaining contributions to the mass loss acceleration for entire Greenland are statistically insignificant. Finally and most importantly, this study suggests the existence of a substantial meltwater storage during summer, with a peak value of 80–120 Gt in July. The robustness of this estimate is demonstrated by using both different GRACE-based solutions and different meltwater runoff estimates (namely, RACMO 2.3 and SNOWPACK). Meltwater storage in the ice sheet occurs primarily due to storage in the high-accumulation regions of the southeast (SE) and northwest (NW) parts of Greenland. Analysis of seasonal variations in outlet glacier discharge shows that the contribution of ice discharge to the observed signal is minor (at the level of only a few Gt) and does not explain the intra-annual differences between the total mass and SMB signals.

Here we present a long-term (1850–2200) best estimate of Greenland ice sheet (GrIS) freshwater runoff that improves spatial detail of runoff locations and temporal resolution. Ice discharge is taken from observations since 2000 and assumed constant in time. Surface meltwater runoff is retrieved from regional climate model output for the recent past and parameterized for the future based on significant correlations between runoff and midtropospheric (500 hPa) summer temperature changes over the GrIS. The simplicity of this approach enables assimilation of meltwater runoff into coupled climate models, which is demonstrated here in a case study with the medium-resolution (1?) Community Earth System Model. The model results suggest that the decrease in Atlantic Meridional Overturning Circulation (AMOC) is dominated by warming of the surface ocean and enhanced GrIS freshwater forcing leads to a slightly enhanced (?1.2 sverdrup in the 21st century) weakening of the AMOC.

One of the major impacts of anthropogenic climate change is sea level rise. Reliable estimates of the contribution of ice sheets to future sea level rise are important to policy makers and the civil society. In addition to sea level rise, ice sheet changes can affect the global climate through modified freshwater fluxes in the areas of deep-water convection. Also, ice sheets modify local and large-scale climate through changes in surface albedo and in their own topography. In the past, ice sheets have played a fundamental role in shaping climate and climate transitions. Despite their strong interactions with the climate system, they are not yet standard components of climate models. First attempts have been made in this direction, and it is foreseeable that in several years ice sheets will be included as interactive components of most models. The main challenges for this coupling are related to spatial and temporal resolution, ice sheet initialization, model climate biases, the need for explicit representation of snow/ice surface physics (e.g., albedo evolution, surface melt, refreezing, compaction), and coupling to the ocean component. This article reviews the main processes contributing to the ice sheet mass budget, the suite of ice sheet–climate interactions, and the requirements for modelling them in a coupled system. Focus is given to four major subjects: surface mass balance, ice sheet flow, ocean–ice sheet interaction, and challenges in coupling ice sheet models to climate models. For further resources related to this article, please visit the WIREs website. Conflict of interest: The author has declared no conflicts of interest for this article.

The Glimmer Community Ice Sheet Model (Glimmer-CISM) has been implemented in the Community Earth System Model (CESM). Glimmer-CISM is forced by a surface mass balance (SMB) computed in multiple elevation classes in the CESM land model and downscaled to the ice sheet grid. Ice sheet evolution is governed by the shallow-ice approximation with thermomechanical coupling and basal sliding. This paper describes and evaluates the initial model implementation for the Greenland Ice Sheet (GIS). The ice sheet model was spun up using the SMB from a coupled CESM simulation with preindustrial forcing. The model's sensitivity to three key ice sheet parameters was explored by running an ensemble of 100 GIS simulations to quasi equilibrium and ranking each simulation based on multiple diagnostics. With reasonable parameter choices, the steady-state GIS geometry is broadly consistent with observations. The simulated ice sheet is too thick and extensive, however, in some marginal regions where the SMB is anomalously positive. The top-ranking simulations were continued using surface forcing from CESM simulations of the twentieth century (1850–2005) and twenty-first century (2005–2100, with RCP8.5 forcing). In these simulations the GIS loses mass, with a resulting global-mean sea level rise of 16 mm during 1850–2005 and 60 mm during 2005–2100. This mass loss is caused mainly by increased ablation near the ice sheet margin, offset by reduced ice discharge to the ocean. Projected sea level rise is sensitive to the initial geometry, showing the importance of realistic geometry in the spun-up ice sheet.

A complex earth system model, simulating atmosphere and ocean dynamics, marine biogeochemistry, terrestrial vegetation and ice sheets, was used to study feedbacks between the terrestrial biosphere and climate with a set of long-term climate change ensemble experiments. CO2 emissions were assigned according to historical data and the IPCC SRES scenarios B1, A1B and A2, followed by an exponential decay of the emissions for the period 2100-3000. The experiments give a reasonable reconstruction of the measured CO2 concentrations between 1750 and 2000. Maximum atmospheric CO2 concentrations of 520 ppm (B1), 860 ppm (A1B) and 1680 ppm (A2) were reached between 2200 and 2500. Additional experiments were performed with CO2 emissions and suppressed climate change, as well as an experiment with a prescribed land surface. The experiments were repeated with the vegetation model driven offline, to investigate the effects of climate and CO2 changes separately. The biogeochemical and biogeophysical interactions between terrestrial biosphere and atmosphere were quantified and compared. A decrease of albedo at high latitudes was the most important biogeophysical change. For the A2 scenario experiment, it causes an additional temperature increase of 1 to 2 K for some high latitude regions by the year 3000, but the changes are minor compared to the heating due to CO2 increase. The terrestrial biosphere takes up between 15 and 30% of the CO2 emissions, depending on the scenario and the period considered. The carbon is stored in the tropics and subtropics, where carbon is stored fast, and in the high latitudes, where carbon storage, partly due to forest expansion, is Much slower. By the year 3000, the storage of terrestrial carbon results in a decrease of atmospheric CO2 concentration of 400 ppm, which in turn decreases the global temperature increase by 0.4 K. (C) 2008 Elsevier B.V. All rights reserved. [References: 46]

Transient experiments for the Eemian (128–113 ky BP) were performed with a complex, coupled earth system model, including atmosphere, ocean, terrestrial biosphere and marine biogeochemistry. In order to investigate the effect of land surface parameters (background albedo, vegetation and tree fraction and roughness length) on the simulated changes during the Eemian, simulations with interactive coupling between climate and vegetation were compared with additional experiments in which these feedbacks were suppressed. The experiments show that the influence of land surface on climate is mainly caused by changes in the albedo. For the northern hemisphere high latitudes, land surface albedo is changed partially due to the direct albedo effect of the conversion of grasses into forest, but the indirect effect of forests on snow albedo appears to be the major factor influencing the total absorption of solar radiation. The Western Sahara region experiences large changes in land surface albedo due to the appearance of vegetation between 128 and 120 ky BP. These local land surface albedo changes can be as much as 20%, thereby affecting the local as well as the global energy balance. On a global scale, latent heat loss over land increases more than 10% for 126 ky BP compared to present-day.

Recent observations indicate a high sensitivity of the Greenland Ice Sheet (GrIS) to climate change. We examine the coupling between the GrIS surface mass balance, elevation, and dynamical flow with one of the few coupled GrIS and atmosphere-ocean general circulation models. Bidirectional coupling from the early Holocene reveals a growing present-day GrIS in the absence of anthropogenic forcing. We identify atmospheric sources of biases in the simulated present-day GrIS and assess the GrIS sensitivity to future greenhouse gas forcing through three Representative Concentration Pathways and their extensions and to climate variability. The elevation-surface mass balance feedback contributes to future GrIS mass loss with 8–11% (by 2100), depending on the forcing scenario, and 24–31% (by 2300). Climate variability causes a 2.5 times spread in the magnitude of the simulated present-day GrIS mass trends in a three-member ensemble. Our results represent a first step toward more advanced higher resolution coupled modeling of GrIS and climate evolution.

A new complex earth system model consisting of an atmospheric general circulation model, an ocean general circulation model, a three-dimensional ice sheet model, a marine biogeochemistry model, and a dynamic vegetation model was used to study the long-term response to anthropogenic carbon emissions. The prescribed emissions follow estimates of past emissions for the period 1751–2000 and standard IPCC emission scenarios up to the year 2100. After 2100, an exponential decrease of the emissions was assumed. For each of the scenarios, a small ensemble of simulations was carried out. The North Atlantic overturning collapsed in the high emission scenario (A2) simulations. In the low emission scenario (B1), only a temporary weakening of the deep water formation in the North Atlantic is predicted. The moderate emission scenario (A1B) brings the system close to its bifurcation point, with three out of five runs leading to a collapsed North Atlantic overturning circulation. The atmospheric moisture transport predominantly contributes to the collapse of the deep water formation. In the simulations with collapsed deep water formation in the North Atlantic a substantial cooling over parts of the North Atlantic is simulated. Anthropogenic climate change substantially reduces the ability of land and ocean to sequester anthropogenic carbon. The simulated effect of a collapse of the deep water formation in the North Atlantic on the atmospheric CO2 concentration turned out to be relatively small. The volume of the Greenland ice sheet is reduced, but its contribution to global mean sea level is almost counterbalanced by the growth of the Antarctic ice sheet due to enhanced snowfall. The modifications of the high latitude freshwater input due to the simulated changes in mass balance of the ice sheet are one order of magnitude smaller than the changes due to atmospheric moisture transport. After the year 3000, the global mean surface temperature is predicted to be almost constant due to the compensating effects of decreasing atmospheric CO2 concentrations due to oceanic uptake and delayed response to increasing atmospheric CO2 concentrations before.

A complex earth system model (atmosphere and ocean general circulation models, ocean biogeochemistry and terrestrial biosphere) was used to perform transient simulations of two interglacial sections (Eemian, 128–113 ky B.P., and Holocene, 9 ky B.P.–present). The changes in terrestrial carbon storage during these interglacials were studied with respect to changes in the earth's orbit. The effects of different climate factors on changes in carbon storage were studied in offline experiments in which the vegetation model was forced only with temperature, hydrological parameters, radiation, or CO2 concentration from the transient runs. The largest anomalies in terrestrial carbon storage were caused by temperature changes. However, the increase in storage due to forest expansion and increased photosynthesis in the high latitudes was nearly balanced by the decrease due to increased respiration. Large positive effects on carbon storage were caused by an enhanced monsoon circulation in the subtropics between 128 and 121 ky B.P. and between 9 and 6 ky B.P., and by increases in incoming radiation during summer for 45° to 70° N compared to a control simulation with present-day insolation. Compared to this control simulation, the net effect of these changes was a positive carbon storage anomaly in the terrestrial biosphere of about 200 Pg C for 125 ky B.P. and 7 ky B.P., and a negative anomaly around 150 Pg C for 116 ky B.P. Although the net increases for Eemian and Holocene were rather similar, the magnitudes of the processes causing these effects were different. The decrease in terrestrial carbon storage during the experiments was the main driver of an increase in atmospheric CO2 concentration during both the Eemian and the Holocene.

Surface mass balance (SMB) trends influence observed Greenland Ice Sheet (GrIS) mass loss, but the component of these trends related to anthropogenic forcing is unclear. Here we study the simulated spatial pattern of emergence of an anthropogenically derived GrIS SMB signal between 1850 and 2100 using the Community Earth System Model. We find emergence timing heterogeneity, with a bimodal structure reflecting interior snowfall increases against a background of low SMB variability, and peripheral surface melting increases against a backdrop of high SMB variability. We also find a nonemerging intermediate region. We conclude that (1) a bimodal pattern of GrIS SMB change will unambiguously reflect the impact of anthropogenic forcing; (2) present-day peripheral and interior SMB trends likely have an underlying anthropogenically forced component; (3) local emergence occurs well before emergence of a spatially integrated signal; and (4) the GrIS summit region may be an ideal location for monitoring regional/global climate change.

The integrated surface mass balance (SMB) of the Greenland ice sheet (GrIS) has large interannual variability. Long-term future changes to this variability will affect GrIS dynamics, freshwater fluxes, regional oceanography, and detection of changes in ice volume trends. Here we analyze a simulated 1850–2100 GrIS SMB time series from the Community Earth System Model, currently the only global climate model that realistically simulates GrIS SMB. We find a significant increase in interannual integrated SMB variability over time, which we attribute primarily to a shift to a high-variability melt-dominated SMB regime due to GrIS ablation area growth. We find temporal increases to characteristic ablation and accumulation area-specific SMB variabilities to be of secondary importance. Since ablation area SMB variability is driven largely by variability in summer surface melt, variability in the climate processes regulating the energy fluxes that control melting will likely increasingly determine future GrIS SMB variability.

The future evolution of global ice sheets under anthropogenic greenhouse forcing and its impact on the climate system, including the regional climate of the ice sheets, are investigated with a comprehensive earth system model consisting of a coupled Atmosphere–Ocean General Circulation Model, a dynamic vegetation model and an ice sheet model. The simulated control climate is realistic enough to permit a direct coupling of the atmosphere and ice sheet components, avoiding the use of anomaly coupling, which represents a strong improvement with respect to previous modelling studies. Glacier ablation is calculated with an energy-balance scheme, a more physical approach than the commonly used degree-day method. Modifications of glacier mask, topographic height and freshwater fluxes by the ice sheets influence the atmosphere and ocean via dynamical and thermodynamical processes. Several simulations under idealized scenarios of greenhouse forcing have been performed, where the atmospheric carbon dioxide stabilizes at two and four times pre-industrial levels. The evolution of the climate system and the ice sheets in the simulations with interactive ice sheets is compared with the simulations with passively coupled ice sheets. For a four-times CO2 scenario forcing, a faster decay rate of the Greenland ice sheet is found in the non-interactive case, where melting rates are higher. This is caused by overestimation of the increase in near-surface temperature that follows the reduction in topographic height. In areas close to retreating margins, melting rates are stronger in the interactive case, due to changes in local albedo. Our results call for careful consideration of the feedbacks operating between ice sheets and climate after substantial decay of the ice sheets.

[1] A complex Earth system model has been forced by insolation changes for the last interglacial (LIG) for which geological records evidence a climate cooling that culminated in the last glacial inception. During an early warm period (127–125 ka B.P.) the Arctic/North Atlantic drainage basin north of 30°N receives ∼27,000 m3/s more fresh water mainly imported across the watershed of North America compared to a cooler period centered around 115 ka B.P. Together with a strong surface warming by up to 2 K, this results in a weaker North Atlantic overturning (reduced by ∼20%). Sea ice cover is reduced by 3 × 106 km2 in the Northern Hemisphere, stimulating a stronger productivity in the Barents and Kara seas. In the North Pacific a lowered salinity and a surface warming of more than 3 K increase the stratification. The eastern Pacific becomes nutrient poorer, which leads to a reduction in the export production in areas with vigorous upwelling. The climate cooling in the course of the LIG induces a carbon release by the terrestrial biosphere of 308 Gt, which is counteracted by an oceanic/sedimentary uptake of 290 Gt. The rest remains in the atmosphere increasing the pCO2 by 9 ppm. The large oceanic CO2 uptake corresponds to a rise of the calcite lysocline of about 300 m in the midlatitude Atlantic. The chemically interactive carbonate sediment pool in the Atlantic is reduced by 10%. For nearly the entire North Atlantic, our results are contradictory to the common interpretation of carbonate preservation proxies as a reliable tool for monitoring changes in deep water circulation.

[1] We investigate the effect of ice sheets on climate change under elevated atmospheric CO2 concentrations with an atmosphere ocean general circulation model (AOGCM) coupled to a thermomechanical ice sheet model and a vegetation model. The effect of increased meltwater fluxes from ice sheets turned out to be negligible in the phase of initial weakening of the North Atlantic meridional overturning circulation (AMOC), and more important during the recovery in subsequent centuries. Lower surface height of the Greenland ice sheet (GRIS) leads locally to a warming, especially in winter, and remotely to a cooling over northern Eurasia due to modified atmospheric circulation. With quadrupling of the atmospheric CO2 concentration the entire GRIS is exposed to surface melt in summer. On formerly ice-covered grid points climate locally warms strongly via increased albedos, with positive feedbacks due to boreal forest expansion.

The presently developed MPI/UW 3D — Earth system model for long-term integrations is applied to simulate the climate of the last interglacial. The model consists of an atmospheric and oceanic general circulation model, a dynamical terrestrial vegetation model and a marine carbon cycle model. The model was forced with time-varying insolation from 129000 to 113000 years before present (yr BP), revealing substantial feedbacks from the land biosphere on climate. These turned out to be important both for the simulated temperature in high northern latitudes and for the precipitation in the northern hemisphere monsoon belt. During the Eemian warm period, the simulated boreal forest extends in many places till the Arctic Ocean, and during the following cold period the transition between tundra and taiga migrated further south. Furthermore, associated albedo changes strongly amplify the simulated temperature changes. The intensified summer insolation during the Eemian leads to a higher precipitation over continents in the northern hemisphere. The strongest response is seen in the tropics and the African-Asian monsoon belt due to increased land-sea temperature contrasts. Vegetation is established in the Western Sahara desert. Compared to simulations with land vegetation prescribed at presentday pattern, the amount of precipitation in the Sahara is more than twice as large. The simulations show strong impacts on the time-transient climate response by triggering nonlinear delays and accelerations seen in various atmospheric and oceanic temperature time series. The simulated storage of carbon in the terrestrial biosphere is relatively large. The carbon storage in land vegetation is increased by more than 10% during the Eemian compared to the following cold period. The associated changes in storage in soil and litter account for more than 100 GtC.

[1] A complex Earth system model including atmosphere, ocean, ice sheets, marine carbon cycle and terrestrial vegetation was used to study the long-term response (100–1000 yrs) of the climate to different increased atmospheric CO2 concentrations. A 3.2 K global mean surface temperature increase is simulated for a 3xCO2 experiment. The freshwater input by melting of the Greenland Ice Sheet due to global warming is of minor importance compared to hydrological changes in the atmosphere. Increased equatorial upwelling enhances the tropical outgassing of CO2 from the oceans, lowering the total marine carbon uptake by 16–22%. On land, carbon release due to increase in soil temperature reduces the anthropogenic carbon uptake from CO2 fertilization up to 43%. Thus, we show that both marine and terrestrial carbon cycle have a positive feedback on climate, which has to be considered for future carbon emission scenarios.

Several multi-century and multi-millennia simulations have been performed with a complex Earth System Model (ESM) for different anthropogenic climate change scenarios in order to study the long-term evolution of sea level and the impact of ice sheet changes on the climate system. The core of the ESM is a coupled coarse-resolution Atmosphere-Ocean General Circulation Model (AOGCM). Ocean biogeochemistry, land vegetation and ice sheets are included as components of the ESM. The Greenland Ice Sheet (GrIS) decays in all simulations, while the Antarctic ice sheet contributes negatively to sea level rise, due to enhanced storage of water caused by larger snowfall rates. Freshwater flux increases from Greenland are one order of magnitude smaller than total freshwater flux increases into the North Atlantic basin (the sum of the contribution from changes in precipitation, evaporation, run-off and Greenland meltwater) and do not play an important role in changes in the strength of the North Atlantic Meridional Overturning Circulation (NAMOC). The regional climate change associated with weakening/collapse of the NAMOC drastically reduces the decay rate of the GrIS. The dynamical changes due to GrIS topography modification driven by mass balance changes act first as a negative feedback for the decay of the ice sheet, but accelerate the decay at a later stage. The increase of surface temperature due to reduced topographic heights causes a strong acceleration of the decay of the ice sheet in the long term. Other feedbacks between ice sheet and atmosphere are not important for the mass balance of the GrIS until it is reduced to 3/4 of the original size. From then, the reduction in the albedo of Greenland strongly accelerates the decay of the ice sheet.