scholarly journals Coupling the U.K. Earth System Model to dynamic models of the Greenland and Antarctic ice sheets

Author(s):  
Robin Smith ◽  
Pierre Mathiot ◽  
Antony Siahaan ◽  
Victoria Lee ◽  
Stephen Cornford ◽  
...  

<p>In this presentation we describe how models of the Greenland and Antarctic ice sheets have been incorporated in the global U.K. Earth System model (UKESM1) with a two-way coupling that passes fluxes of energy, water and the locations of ice surfaces between the component models. Offline, file-based coupling is used throughout to pass information between the components, which is both physically appropriate and convenient within the UKESM1 structure. Ice sheet surface mass balance is computed in the land surface model using sub-gridscale multi-layer snowpacks. Icebergs calved from the ice sheets are fed into a Langrangian iceberg drift scheme in the ocean. Ice shelf basal melt is explicitly calculated in cavities resolved by the ocean model, and ice sheet and shelf geometries are kept consistent in all components. We demonstrate that our coupled model remains stable when simulating changes in ice sheet height, extent and grounding-line position of hundreds of kilometres.</p>


Author(s):  
Robin S. Smith ◽  
Pierre Mathiot ◽  
Antony Siahaan ◽  
Victoria Lee ◽  
Stephen L. Cornford ◽  
...  


2014 ◽  
Vol 28 (1) ◽  
pp. 272-291 ◽  
Author(s):  
Daniela Dalmonech ◽  
Sönke Zaehle ◽  
Gregor J. Schürmann ◽  
Victor Brovkin ◽  
Christian Reick ◽  
...  

Abstract The capacity of earth system models (ESMs) to make reliable projections of future atmospheric CO2 and climate is strongly dependent on the ability of the land surface model to adequately simulate the land carbon (C) cycle. Defining “adequate” performance of the land model requires an understanding of the contributions of climate model and land model errors to the land C cycle. Here, a benchmarking framework is applied based on significant, observed characteristics of the land C cycle for the contemporary period, for which sufficient evaluation data are available, to test the ability of the JSBACH land surface component of the Max Planck Institute Earth System Model (MPI-ESM) to simulate land C trends. Particular attention is given to the role of potential effects caused by climate biases, and therefore investigation is made of the results of model configurations in which JSBACH is interactively “coupled” to atmosphere and ocean components and of an “uncoupled” configuration, where JSBACH is driven by reconstructed meteorology. The ability of JSBACH to simulate the observed phase of phenology and seasonal C fluxes is not strongly affected by climate biases. Contrarily, noticeable differences in the simulated gross primary productivity and land C stocks emerge between coupled and uncoupled configurations, leading to significant differences in the decadal terrestrial C balance and its sensitivity to climate. These differences are strongly controlled by climate biases of the MPI-ESM, in particular those affecting soil moisture. To effectively characterize model performance, the potential effects of climate biases on the land C dynamics need to be considered during the development and calibration of land surface models.



2021 ◽  
Author(s):  
Sarah L Bradley ◽  
Michele Petrini ◽  
Raymond Sellevold ◽  
Miren Vizcaino ◽  
William H. Lipscomb ◽  
...  

<p>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.</p><p> </p><p>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).</p><p> </p><p>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.</p>



2014 ◽  
Vol 7 (5) ◽  
pp. 2003-2013 ◽  
Author(s):  
D. Barbi ◽  
G. Lohmann ◽  
K. Grosfeld ◽  
M. Thoma

Abstract. We present first results from a coupled model setup, consisting of the state-of-the-art ice sheet model RIMBAY (Revised Ice Model Based on frAnk pattYn), and the community earth system model COSMOS. We show that special care has to be provided in order to ensure physical distributions of the forcings as well as numeric stability of the involved models. We demonstrate that a suitable statistical downscaling is crucial for ice sheet stability, especially for southern Greenland where surface temperatures are close to the melting point. The downscaling of net snow accumulation is based on an empirical relationship between surface slope and rainfall. The simulated ice sheet does not show dramatic loss of ice volume for pre-industrial conditions and is comparable with present-day ice orography. A sensitivity study with high CO2 level is used to demonstrate the effects of dynamic ice sheets onto climate compared to the standard setup with prescribed ice sheets.



2020 ◽  
Author(s):  
Adam M. Schneider ◽  
Charles S. Zender ◽  
Stephen F. Price

Abstract. Earth's largest island, Greenland, and the Antarctic continent are both covered by massive ice sheets. A large fraction of their surfaces consist of multi-year snow, known as firn, which has undergone a process of densification since falling from the atmosphere. Until now this firn densification has not been fully accounted for in the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). Here, we expand the E3SM Land Model (ELM) snowpack from 1 m to up to 60 m to enable more accurate simulation of snowpack evolution. We test four densification models in a series of century-scale land surface simulations forced by atmospheric re-analyses, and evaluate these parameterizations against empirical density-versus-depth data. To tailor candidate densification models for use across the ice sheets' dry-snow zones, we optimize parameters using a regularized least squares algorithm applied to two distinct stages of densification. We find that a dynamic implementation of a semi-empirical compaction model, originally calibrated to measurements from the Antarctic peninsula, gives results more consistent with ice core measurements from the cold, dry snow zones of Greenland and Antarctica, compared to when using the original ELM snow compaction physics. In its latest release, the Community Land Model (CLM) (version 5) provides updated snow compaction physics that we test in ELM, resulting in top 10 m firn densities that are in better agreement with observations than densities simulated with the semi-empirical model. Below 10 m, however, the semi-empirical model gives results that more closely match observations, while the current CLM(v5) compaction physics predict firn densities that increase too slowly with depth and are thus unable to simulate pore close off (a phenomenon of particular interest to paleoclimate studies). Because snow and firn density play roles in snowpack albedo, liquid water storage, and ice sheet surface mass balance, these improvements will contribute to broader E3SM efforts to simulate the response of land ice to atmospheric forcing and the resulting impacts on global sea level.



2010 ◽  
Vol 3 (2) ◽  
pp. 603-633 ◽  
Author(s):  
H. Goosse ◽  
V. Brovkin ◽  
T. Fichefet ◽  
R. Haarsma ◽  
P. Huybrechts ◽  
...  

Abstract. The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The ocean component is CLIO3, which consists of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is of 3° by 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the ocean carbon cycle is represented by LOCH, a comprehensive model that takes into account both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the ice-atmosphere and ice-ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, and an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.



2021 ◽  
Author(s):  
Darin Comeau ◽  
Xylar Asay-Davis ◽  
Carolyn Begeman ◽  
Matthew Hoffman ◽  
Wuyin Lin ◽  
...  

<p>The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS) -- ice-shelf basal melting and iceberg calving -- are generally poorly represented in current Earth System Models (ESMs). Here, we document the first effort to date at simulating the ocean circulation and exchanges of heat and freshwater within ice-shelf cavities in a coupled ESM, the Department of Energy's Energy Exascale Earth System Model (E3SM). As a step towards full ice-sheet coupling, we implemented static Antarctic ice-shelf cavities and the ability to calculate ice-shelf basal melt rates from the heat and freshwater fluxes computed by the ocean model. In addition, we added the capability to prescribe forcing from iceberg melt, allowing us to realistically represent the other dominant mass loss process from the AIS. In global, low resolution (i.e., non-eddying ocean) simulations, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a tipping point to high melt regimes under certain ice shelves, presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy-induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic and stable ice-shelf basal melt rates across the continent under pre-industrial climate forcing. We also show preliminary results using an ocean/sea-ice grid that makes use of E3SM’s regional-refinement capability, where increased resolution (down to 12km) is placed in the Southern Ocean around Antarctica, bypassing the need for parameterization of eddy-induced transport in this region. The accurate representation of these processes within a coupled ESM is an important step towards reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea-level rise.</p>



2013 ◽  
Vol 10 (3) ◽  
pp. 1659-1674 ◽  
Author(s):  
R. J. Schuldt ◽  
V. Brovkin ◽  
T. Kleinen ◽  
J. Winderlich

Abstract. Since the Last Glacial Maximum, boreal wetlands have accumulated substantial amounts of peat, estimated at 180–621 Pg of carbon. Wetlands have significantly affected the atmospheric greenhouse gas composition in the past and will play a significant role in future changes of atmospheric CO2 and CH4 concentrations. In order to investigate those changes with an Earth system model, biogeochemical processes in boreal wetlands need to be accounted for. Thus, a model of peat accumulation and decay was developed and included in the land surface model JSBACH of the Max Planck Institute Earth System Model (MPI-ESM). Here we present the evaluation of model results from 6000 yr BP to the pre-industrial period. Over this period of time, 240 Pg of peat carbon accumulated in the model in the areas north of 40° N. Simulated peat accumulation rates agree well with those reported for boreal wetlands. The model simulates CH4 emissions of 49.3 Tg CH4 yr−1 for 6000 yr BP and 51.5 Tg CH4 yr−1 for pre-industrial times. This is within the range of estimates in the literature, which range from 32 to 112 Tg CH4 yr−1 for boreal wetlands. The modelled methane emission for the West Siberian Lowlands and Hudson Bay Lowlands agree well with observations. The rising trend of methane emissions over the last 6000 yr is in agreement with measurements of Antarctic and Greenland ice cores.



2021 ◽  
Author(s):  
Alessandro Anav ◽  
Adriana Carillo ◽  
Massimiliano Palma ◽  
Maria Vittoria Struglia ◽  
Ufuk Utku Turuncoglu ◽  
...  

Abstract. In this study, a new regional Earth system model is developed and applied to the Med-CORDEX region. The ENEA-REG system is made up of two interchangeable regional climate models as atmospheric components (RegCM and WRF), a river model (HD), and an ocean model (MITgcm); processes taking place at the land surface are represented within the atmospheric models with the possibility to use several land surface schemes of different complexity. The coupling between these components is performed through the RegESM driver. Here, we present and describe our regional Earth system model and evaluate its components using a multidecadal hindcast simulation over the period 1980–2013 driven by ERA-INTERIM reanalysis. We show how the atmospheric components are able to correctly reproduce both large-scale and local features of the Euro-Mediterranean climate, although some remarkable biases are relevant for some variables. In particular, WRF has a significant cold bias during winter over North-Eastern bound of the domain, while RegCM systematically overestimates the wind speed over the Mediterranean Sea. This latter bias has severe consequences on the ocean component: we show that when WRF is used as the atmospheric component of the Earth system, the performances of the ocean model are remarkably better compared with the RegCM version. Our regional Earth system model allows studying the Euro-Mediterranean climate system and can be applied to both hindcast and scenario simulations.



Sign in / Sign up

Export Citation Format

Share Document