scholarly journals FAMOUS version xotzt (FAMOUS-ice): a general circulation model (GCM) capable of energy- and water-conserving coupling to an ice sheet model

2021 ◽  
Vol 14 (9) ◽  
pp. 5769-5787
Author(s):  
Robin S. Smith ◽  
Steve George ◽  
Jonathan M. Gregory

Abstract. The physical interactions between ice sheets and their surroundings are major factors in determining the state of the climate system, yet many current Earth system models omit them entirely or approximate them in a heavily parameterised manner. In this work we have improved the snow and ice sheet surface physics in the FAMOUS climate model, with the aim of improving the representation of polar climate and implementing a bidirectional coupling to the Glimmer dynamic ice sheet model using the water and energy fluxes calculated by FAMOUS. FAMOUS and Glimmer are both low-resolution, computationally affordable models used for multi-millennial simulations. Glaciated surfaces in the new FAMOUS-ice are modelled using a multi-layer snow scheme capable of simulating compaction of firn and the percolation and refreezing of surface melt. The low horizontal resolution of FAMOUS compared to Glimmer is mitigated by implementing this snow model on sub-grid-scale tiles that represent different elevations on the ice sheet within each FAMOUS grid box. We show that with this approach FAMOUS-ice can simulate relevant physical processes on the surface of the modern Greenland ice sheet well compared to higher-resolution climate models and that the ice sheet state in the coupled FAMOUS-ice–Glimmer system does not drift unacceptably. FAMOUS-ice coupled to Glimmer is thus a useful tool for modelling the physics and co-evolution of climate and grounded ice sheets on centennial and millennial timescales, with applications to scientific questions relevant to both paleoclimate and future sea level rise.

2020 ◽  
Author(s):  
Robin S. Smith ◽  
Steve George ◽  
Jonathan M. Gregory

Abstract. The physical interactions between ice sheets and their surroundings are major factors in determining the state of the climate system, yet many current Earth System models omit them entirely or approximate them in a heavily parameterised manner. In this work we have improved the snow and ice sheet surface physics in the FAMOUS climate model, with the aim of improving the representation of polar climate and implementing a bidirectional coupling to the Glimmer dynamic ice sheet model using the water and energy fluxes calculated by FAMOUS. FAMOUS and Glimmer are both low resolution, computationally affordable models used for multi-millennial simulations. Glaciated surfaces in the new FAMOUS-ice are modelled using a multi-layer snow scheme capable of simulating compaction of firn and the percolation and refreezing of surface melt. The low horizontal resolution of FAMOUS compared to Glimmer is mitigated by implementing this snow model on sub-gridscale tiles that represent different elevations on the ice sheet within each FAMOUS grid-box. We show that with this approach FAMOUS-ice can simulate relevant physical processes on the surface of the modern Greenland ice sheet well compared to higher resolution climate models, and that the ice sheet state in the coupled FAMOUS-ice-Glimmer system does not drift unacceptably. FAMOUS-ice coupled to Glimmer is thus a useful tool for modelling the physics and co-evolution of climate and grounded ice sheets on centennial and millennial timescales, with applications to scientific questions relevant to both paleoclimate and future sea level rise.


2012 ◽  
Vol 6 (2) ◽  
pp. 1037-1083 ◽  
Author(s):  
A. Quiquet ◽  
H. J. Punge ◽  
C. Ritz ◽  
X. Fettweis ◽  
M. Kageyama ◽  
...  

Abstract. The prediction of future climate and ice sheet evolution requires coupling of ice sheet and climate models. Before proceeding to a coupled setup, we propose to analyze the impact of model simulated climate on an ice sheet. Here, we undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary condition to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyr of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed ones, there are considerable deviations among the ice sheets on regional scales. These can be explained by difficulties in modelling local temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations of each climate model are also due to the differences in the atmospheric general circulation. In the context of coupling ice sheet and climate models, we conclude that appropriate downscaling methods will be needed and systematic corrections of the climatic variables at the interface may be required in some cases to obtain realistic results for the Greenland ice sheet (GIS).


2015 ◽  
Vol 11 (3) ◽  
pp. 403-424 ◽  
Author(s):  
A. M. Dolan ◽  
S. J. Hunter ◽  
D. J. Hill ◽  
A. M. Haywood ◽  
S. J. Koenig ◽  
...  

Abstract. During an interval of the Late Pliocene, referred to here as the mid-Pliocene Warm Period (mPWP; 3.264 to 3.025 million years ago), global mean temperature was similar to that predicted for the end of this century, and atmospheric carbon dioxide concentrations were higher than pre-industrial levels. Sea level was also higher than today, implying a significant reduction in the extent of the ice sheets. Thus, the mPWP provides a natural laboratory in which to investigate the long-term response of the Earth's ice sheets and sea level in a warmer-than-present-day world. At present, our understanding of the Greenland ice sheet during the mPWP is generally based upon predictions using single climate and ice sheet models. Therefore, it is essential that the model dependency of these results is assessed. The Pliocene Model Intercomparison Project (PlioMIP) has brought together nine international modelling groups to simulate the warm climate of the Pliocene. Here we use the climatological fields derived from the results of the 15 PlioMIP climate models to force an offline ice sheet model. We show that mPWP ice sheet reconstructions are highly dependent upon the forcing climatology used, with Greenland reconstructions ranging from an ice-free state to a near-modern ice sheet. An analysis of the surface albedo variability between the climate models over Greenland offers insights into the drivers of inter-model differences. As we demonstrate that the climate model dependency of our results is high, we highlight the necessity of data-based constraints of ice extent in developing our understanding of the mPWP Greenland ice sheet.


2018 ◽  
Author(s):  
Constantijn J. Berends ◽  
Bas de Boer ◽  
Roderik S. W. van de Wal

Abstract. Fully coupled ice-sheet-climate modelling over 10,000–100,000-year time scales on high spatial and temporal resolution remains beyond the capability of current computational systems. Hybrid GCM-ice-sheet modelling offers a middle ground, balancing the need to accurately capture both long-term processes, in particular circulation driven changes in precipitation, and processes requiring a high spatial resolution like ablation. Here, we present and evaluate a model set-up that forces the ANICE 3D thermodynamic ice-sheet-shelf model calculating all ice on Earth, with pre-calculated output from several steady-state simulations with the HadCM3 general circulation model (GCM), using a so-called matrix method of coupling both components, where simulations with various levels of pCO2 and ice-sheet configuration are combined to form a time-continuous transient climate forcing consistent with the modelled ice-sheets. We address the difficulties in downscaling low-resolution GCM output to the higher-resolution grid of an ice-sheet model, and account for differences between GCM and ice-sheet model surface topography ranging from interglacial to glacial conditions. As a benchmark experiment to assess the validity of this model set-up, we perform a simulation of the entire last glacial cycle, from 120 kyr ago to present-day. The simulated eustatic sea-level drop at the Last Glacial maximum (LGM) for the combined Antarctic, Greenland, Eurasian and North-American ice-sheets amounts to 100 m, in line with many other studies. The simulated ice-sheets at LGM agree well with the ICE-5G reconstruction and the more recent DATED-1 reconstruction in terms of total volume and geographical location of the ice sheets. Moreover, modelled benthic oxygen isotope abundance and the relative contributions from global ice volume and deep-water temperature agree well with available data, as do surface temperature histories for the Greenland and Antarctic ice-sheets. This model strategy can be used to create time-continuous ice-sheet distribution and sea-level reconstructions for geological periods up to several millions of years in duration, capturing climate model driven variations in the mass balance of the ice sheet.


2014 ◽  
Vol 27 (13) ◽  
pp. 4835-4856 ◽  
Author(s):  
Richard I. Cullather ◽  
Sophie M. J. Nowicki ◽  
Bin Zhao ◽  
Max J. Suarez

Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS-5), atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980–2008 GrIS SMB average is 24.7 ± 4.5 cm yr−1 water-equivalent (w.e.) at ½° model grid spacing, and 18.2 ± 3.3 cm yr−1 w.e. for 2° grid spacing. The spatial variability and seasonal cycle of the ½° simulation compare favorably to recent studies using regional climate models, while results from 2° integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser-resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.


2018 ◽  
Vol 12 (4) ◽  
pp. 1499-1510 ◽  
Author(s):  
Marcus Lofverstrom ◽  
Johan Liakka

Abstract. Coupled climate–ice sheet simulations have been growing in popularity in recent years. Experiments of this type are however challenging as ice sheets evolve over multi-millennial timescales, which is beyond the practical integration limit of most Earth system models. A common method to increase model throughput is to trade resolution for computational efficiency (compromise accuracy for speed). Here we analyze how the resolution of an atmospheric general circulation model (AGCM) influences the simulation quality in a stand-alone ice sheet model. Four identical AGCM simulations of the Last Glacial Maximum (LGM) were run at different horizontal resolutions: T85 (1.4∘), T42 (2.8∘), T31 (3.8∘), and T21 (5.6∘). These simulations were subsequently used as forcing of an ice sheet model. While the T85 climate forcing reproduces the LGM ice sheets to a high accuracy, the intermediate resolution cases (T42 and T31) fail to build the Eurasian ice sheet. The T21 case fails in both Eurasia and North America. Sensitivity experiments using different surface mass balance parameterizations improve the simulations of the Eurasian ice sheet in the T42 case, but the compromise is a substantial ice buildup in Siberia. The T31 and T21 cases do not improve in the same way in Eurasia, though the latter simulates the continent-wide Laurentide ice sheet in North America. The difficulty to reproduce the LGM ice sheets in the T21 case is in broad agreement with previous studies using low-resolution atmospheric models, and is caused by a substantial deterioration of the model climate between the T31 and T21 resolutions. It is speculated that this deficiency may demonstrate a fundamental problem with using low-resolution atmospheric models in these types of experiments.


2018 ◽  
Vol 11 (11) ◽  
pp. 4657-4675 ◽  
Author(s):  
Constantijn J. Berends ◽  
Bas de Boer ◽  
Roderik S. W. van de Wal

Abstract. Fully coupled ice-sheet–climate modelling over 10 000–100 000-year timescales at high spatial and temporal resolution remains beyond the capability of current computational systems. Forcing an ice-sheet model with precalculated output from a general circulation model (GCM) offers a middle ground, balancing the need to accurately capture both long-term processes, in particular circulation-driven changes in precipitation, and processes requiring a high spatial resolution like ablation. Here, we present and evaluate a model set-up that forces the ANICE 3-D thermodynamic ice-sheet–shelf model calculating the four large continental ice sheets (Antarctica, Greenland, North America, and Eurasia) with precalculated output from two steady-state simulations with the HadCM3 (GCM) using a so-called matrix method of coupling both components, whereby simulations with various levels of pCO2 and ice-sheet configuration are combined to form a time-continuous transient climate forcing consistent with the modelled ice sheets. We address the difficulties in downscaling low-resolution GCM output to the higher-resolution grid of an ice-sheet model and account for differences between GCM and ice-sheet model surface topography ranging from interglacial to glacial conditions. Although the approach presented here can be applied to a matrix with any number of GCM snapshots, we limited our experiments to a matrix of only two snapshots. As a benchmark experiment to assess the validity of this model set-up, we perform a simulation of the entire last glacial cycle from 120 kyr ago to present day. The simulated eustatic sea-level drop at the Last Glacial Maximum (LGM) for the combined Antarctic, Greenland, Eurasian, and North American ice sheets amounts to 100 m, in line with many other studies. The simulated ice sheets at the LGM agree well with the ICE-5G reconstruction and the more recent DATED-1 reconstruction in terms of total volume and geographical location of the ice sheets. Moreover, modelled benthic oxygen isotope abundance and the relative contributions from global ice volume and deep-water temperature agree well with available data, as do surface temperature histories for the Greenland and Antarctic ice sheets. This model strategy can be used to create time-continuous ice-sheet distribution and sea-level reconstructions for geological periods up to several million years in duration, capturing climate-model-driven variations in the mass balance of the ice sheet.


2012 ◽  
Vol 6 (4) ◽  
pp. 743-762 ◽  
Author(s):  
C. H. Reijmer ◽  
M. R. van den Broeke ◽  
X. Fettweis ◽  
J. Ettema ◽  
L. B. Stap

Abstract. Retention and refreezing of meltwater are acknowledged to be important processes for the mass budget of polar glaciers and ice sheets. Several parameterizations of these processes exist for use in energy and mass balance models. Due to a lack of direct observations, validation of these parameterizations is difficult. In this study we compare a set of 6 refreezing parameterizations against output of two Regional Climate Models (RCMs) coupled to an energy balance snow model, the Regional Atmospheric Climate Model (RACMO2) and the Modèle Atmosphérique Régional (MAR), applied to the Greenland ice sheet. In both RCMs, refreezing is explicitly calculated in a snow model that calculates vertical profiles of temperature, density and liquid water content. Between RACMO2 and MAR, the ice sheet-integrated amount of refreezing differs by only 4.9 mm w.e yr−1 (4.5 %), and the temporal and spatial variability are very similar. For consistency, the parameterizations are forced with output (surface temperature, precipitation and melt) of the RCMs. For the ice sheet-integrated amount of refreezing and its inter-annual variations, all parameterizations give similar results, especially after some tuning. However, the spatial distributions differ significantly and the spatial correspondence between the RCMs is better than with any of the parameterizations. Results are especially sensitive to the choice of the depth of the thermally active layer, which determines the cold content of the snow in most parameterizations. These results are independent of which RCM is used to force the parameterizations.


2014 ◽  
Vol 10 (4) ◽  
pp. 3483-3535 ◽  
Author(s):  
A. M. Dolan ◽  
S. J. Hunter ◽  
D. J. Hill ◽  
A. M. Haywood ◽  
S. J. Koenig ◽  
...  

Abstract. During the mid-Pliocene Warm Period (3.264 to 3.025 million years ago), global mean temperature was similar to that predicted for the end of this century, and atmospheric carbon dioxide concentrations were higher than pre-industrial levels. Sea level was also higher than today, implying a significant reduction in the extent of the ice sheets. Thus, the mid-Pliocene Warm Period provides a natural laboratory in which to investigate the long-term response of the Earth's ice sheets and sea level in a warmer-than-modern world. At present, our understanding of the Greenland ice sheet during the warmest intervals of the mid-Pliocene is generally based upon predictions using single climate and ice sheet models. Therefore, it is essential that the model dependency of these results is assessed. The Pliocene Model Intercomparison Project (PlioMIP) has brought together nine international modelling groups to simulate the warm climate of the Pliocene. Here we use the climatological fields derived from the results of the PlioMIP climate models to force an offline ice sheet model. We show Pliocene ice sheet reconstructions are highly dependent upon the forcing climatology used, with Greenland reconstructions ranging from an ice-free state to a near modern ice sheet. An analysis of surface albedo differences between the climate models over Greenland offers insights into the drivers of inter-model differences. As we demonstrate that the climate model dependency of our results is high, we highlight the necessity of data-based constraints in developing our understanding of the Pliocene Greenland ice sheet.


2012 ◽  
Vol 6 (5) ◽  
pp. 999-1018 ◽  
Author(s):  
A. Quiquet ◽  
H. J. Punge ◽  
C. Ritz ◽  
X. Fettweis ◽  
H. Gallée ◽  
...  

Abstract. Predicting the climate for the future and how it will impact ice sheet evolution requires coupling ice sheet models with climate models. However, before we attempt to develop a realistic coupled setup, we propose, in this study, to first analyse the impact of a model simulated climate on an ice sheet. We undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary conditions to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyrs of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed one, there are considerable deviations among the ice sheets on regional scales. These deviations can be explained by biases in temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations between the climate models are also due to the differences in the atmospheric general circulation. To account for these differences in the context of coupling ice sheet models with climate models, we conclude that appropriate downscaling methods will be needed. In some cases, systematic corrections of the climatic variables at the interface may be required to obtain realistic results for the Greenland ice sheet (GIS).


Sign in / Sign up

Export Citation Format

Share Document