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).@en

The Last Glacial Maximum (LGM, from ∼26 to 20 ka BP) was the most recent period with large ice sheets in Eurasia and North America. At that time, global temperatures were 5–7 ∘C lower than today, and sea level ∼125 m lower. LGM simulations are useful to understand earth system dynamics, including climate–ice sheet interactions, and to evaluate and improve the models representing those dynamics. Here, we present two simulations of the Northern Hemisphere ice sheet climate and surface mass balance (SMB) with the Community Earth System Model v2.1 (CESM2.1) using the Community Atmosphere Model v5 (CAM5) with prescribed ice sheets for two time periods that bracket the LGM period: 26 and 21 ka BP. CESM2.1 includes an explicit simulation of snow/firn compaction, albedo, refreezing, and direct coupling of the ice sheet surface energy fluxes with the atmosphere. The simulated mean snow accumulation is lowest for the Greenland and Barents–Kara Sea ice sheets (GrIS, BKIS) and highest for British and Irish (BIIS) and Icelandic (IcIS) ice sheets. Melt rates are negligible for the dry BKIS and GrIS, and relatively large for the BIIS, North American ice sheet complex (NAISC; i.e. Laurentide, Cordilleran, and Innuitian), Scandinavian ice sheet (SIS), and IcIS, and are reduced by almost a third in the colder (lower temperature) 26 ka BP climate compared with 21 ka BP. The SMB is positive for the GrIS, BKIS, SIS, and IcIS during the LGM (26 and 21 ka BP) and negative for the NAISC and BIIS. Relatively wide ablation areas are simulated along the southern (terrestrial), Pacific and Atlantic margins of the NAISC, across the majority of the BIIS, and along the terrestrial southern margin of the SIS. The integrated SMB substantially increases for the NAISC and BIIS in the 26 ka BP climate, but it does not reverse the negative sign. Summer incoming surface solar radiation is largest over the high interior of the NAISC and GrIS, and minimum over the BIIS and southern margin of NAISC. Summer net radiation is maximum over the ablation areas and minimum where the albedo is highest, namely in the interior of the GrIS, northern NAISC, and all of the BKIS. Summer sensible and latent heat fluxes are highest over the ablation areas, positively contributing to melt energy. Refreezing is largest along the equilibrium line altitude for all ice sheets and prevents 40 %–50 % of meltwater entering the ocean. The large simulated melt for the NAISC suggests potential biases in the climate simulation, ice sheet reconstruction, and/or highly non-equilibrated climate and ice sheet at the LGM time.@en

Earth system/ice-sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea-level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy-balance-based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice-sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.@en

The Greenland ice sheet (GrIS) is now losing mass at a rate of 0.7 mm of sea level rise (SLR) per year. Here we explore future GrIS evolution and interactions with global and regional climate under high greenhouse gas forcing with the Community Earth System Model version 2.1 (CESM2.1), which includes an interactive ice sheet component (the Community Ice Sheet Model v2.1 [CISM2.1]) and an advanced energy balance-based calculation of surface melt. We run an idealized 350-year scenario in which atmospheric CO2 concentration increases by 1% annually until reaching four times pre-industrial values at year 140, after which it is held fixed. The global mean temperature increases by 5.2 and 8.5 K by years 131–150 and 331–350, respectively. The projected GrIS contribution to global mean SLR is 107 mm by year 150 and 1,140 mm by year 350. The rate of SLR increases from 2 mm yr−1 at year 150 to almost 7 mm yr−1 by year 350. The accelerated mass loss is caused by rapidly increasing surface melt as the ablation area expands, with associated albedo feedback and increased sensible and latent heat fluxes. This acceleration occurs for a global warming of approximately 4.2 K with respect to pre-industrial and is in part explained by the quasi-parabolic shape of the ice sheet, which favors rapid expansion of the ablation area as it approaches the interior “plateau.”.@en