Downscaling CESM2 in CLM5 to Hindcast Pre-Industrial Equilibrium Line Altitude for Tropical Mountains

Author(s):  
Nicholas Heavens

<p>Highland environments are rarely preserved in the geological record, particularly from as early as the Paleozoic Era. However, several stratigraphic locations are now known which definitely or potentially preserve such environments near the paleoequator during the Late Carboniferous and Early Permian Periods, during which the Earth was in the depths of an icehouse climate like that of the Pliocene and Pleistocene Epochs, the Late Paleozoic Ice Age (LPIA). Several of these locations contain evidence of mountain glaciation at altitudes below 2000 m, leading to questions about the significance of tropical mountain glaciation for global climate during this interval of geologic time. However, climate model simulations for the LPIA have not been able to simulate mountain glaciation like that inferred from the geological record, possibly because of low resolution, incorrect boundary conditions, or climate model bias resulting from incomplete representation of moist convective processes impacting tropical lapse rates. </p><p>The overarching purpose of this study is to develop a climate modeling framework that enables the significance of mountain glaciation for global paleoclimate to be evaluated. Ideally, such a framework would allow low-resolution global model output to be downscaled to the scale of a mountain range to calculate the equilibrium line altitude and similar parameters, enabling evidence of mountain glaciation in the deep past to be used to constrain/tune the low-resolution global models. While this study was designed to inform a specific problem in deep time paleoclimate, its results are likely broadly applicable to assessing how well mountain glaciation is captured by global climate modeling of the past, present, and future.   </p><p>Here, I present a framework in which the CMIP6 pre-industrial control simulation for the Community Earth System Model version 2 (CESM2) at 0.9°x1.25° resolution is used to generate a data atmosphere for the Community Land Model version 5 (CLM5) run at 0.01° resolution in 10 tropical and 1 mid-latitude domain to study the surface mass balance over the domain. For computational reasons, glaciation is assumed to cover a small portion of each grid cell, but surface mass balance still can be evaluated. Topographic boundary conditions come from GMTED2010, but most other information is directly interpolated from the CESM2 simulation. CLM5 simulations require a fixed lapse rate to be assumed, which is varied in each CLM5 simulation across six different values. The CLM5 simulation output along with the mean tropical lapse rate in the CESM2 simulation is then used to evaluate the various biases of this framework in comparison with estimated pre-industrial equilibrium line altitudes for the studied domains.</p><p>This work is supported by the National Science Foundation (USA) under grant EAR-1849754. </p>

2012 ◽  
Vol 6 (2) ◽  
pp. 255-272 ◽  
Author(s):  
M. M. Helsen ◽  
R. S. W. van de Wal ◽  
M. R. van den Broeke ◽  
W. J. van de Berg ◽  
J. Oerlemans

Abstract. It is notoriously difficult to couple surface mass balance (SMB) results from climate models to the changing geometry of an ice sheet model. This problem is traditionally avoided by using only accumulation from a climate model, and parameterizing the meltwater run-off as a function of temperature, which is often related to surface elevation (Hs). In this study, we propose a new strategy to calculate SMB, to allow a direct adjustment of SMB to a change in ice sheet topography and/or a change in climate forcing. This method is based on elevational gradients in the SMB field as computed by a regional climate model. Separate linear relations are derived for ablation and accumulation, using pairs of Hs and SMB within a minimum search radius. The continuously adjusting SMB forcing is consistent with climate model forcing fields, also for initially non-glaciated areas in the peripheral areas of an ice sheet. When applied to an asynchronous coupled ice sheet – climate model setup, this method circumvents traditional temperature lapse rate assumptions. Here we apply it to the Greenland Ice Sheet (GrIS). Experiments using both steady-state forcing and glacial-interglacial forcing result in realistic ice sheet reconstructions.


2020 ◽  
Author(s):  
Nicolas Ghilain ◽  
Stéphane Vannitsem ◽  
Quentin Dalaiden ◽  
Hugues Goosse

<p>Over recent decades, the Antarctic Ice Sheet has witnessed large spatial variations at its surface through the surface mass balance (SMB). Since the complex Antarctic topography, working at high resolution is crucial to represent accurately the dynamics of SMB. While ice cores provide a mean to infer the SMB over centuries, the view is very spatially constrained. Global Climate models estimate the spatial distribution of SMB over centuries, but with a too coarse resolution with regards to the large variations due to local orographic effects. We have therefore explored a methodology to statistically downscale the SMB components from the climate model historical simulations (1850-present day). An analogue method is set up over a period of 30 years with the ERA-Interim reanalysis (1979-2010 AD) and associated with SMB components from the Regional Atmospheric Climate Model (RACMO) at 5 km spatial resolution over Dronning Maud in East Antarctica. The same method is then applied to the period from 1850 to present days using an ensemble of 10 simulations from the CESM2 model. This method enables to derive a spatial distribution of SMB. In addition, the changes in precipitation delivery mechanisms can be unveiled.</p>


2016 ◽  
Author(s):  
Michiel Helsen ◽  
Roderik Van de Wal ◽  
Thomas Reerink ◽  
Richard Bintanja ◽  
Marianne Sloth Madsen ◽  
...  

Abstract. The albedo of the surface of ice sheets changes as a function of time, due to the effects of deposition of new snow, ageing of dry snow, melting and runoff. Currently, the calculation of the albedo of ice sheets is highly parameterized within the Earth System Model EC-Earth, by taking a constant value for areas with thick perennial snow cover. This is one of the reasons that the surface mass balance (SMB) of the Greenland ice sheet (GrIS) is poorly resolved in the model. To improve this, eight snow albedo schemes are evaluated here. The resulting SMB is downscaled from the lower resolution global climate model topography to the higher resolution ice sheet topography of the GrIS, such that the influence of these different SMB climatologies on the long-term evolution of the GrIS is tested by ice sheet model simulations. This results in an optimised albedo parameterization that can be used in future EC-Earth simulations with an interactive ice sheet component.


2015 ◽  
Vol 9 (3) ◽  
pp. 3165-3204 ◽  
Author(s):  
R. J. Braithwaite

Abstract. Translated into modern terminology, Kurowski suggested in 1891 that the equilibrium line altitude (ELA) of a glacier is equal to the mean altitude of the glacier when the whole glacier is in balance between accumulation and ablation. Kurowski's method has been widely misunderstood, partly due to inappropriate use of statistical terminology by later workers, and has been little tested except by Braithwaite and Müller in a 1980 paper (for 32 glaciers). I now compare Kurowski's mean altitude with balanced-budget ELA calculated for 103 modern glaciers with measured surface mass balance data. Kurowski's mean altitude is significantly higher (at 95% level) than balanced-budget ELA for 19 outlet and 42 valley glaciers, but not significantly higher for 34 mountain glaciers. The error in Kurowski mean altitude as a predictor of balanced-budget ELA might be due to generally lower balance gradients in accumulation area compared with ablation areas for many glaciers, as suggested by several workers, but some glaciers have higher gradients, presumably due to precipitation increase with altitude. The relatively close agreement between balanced-budget ELA and mean altitude for mountain glaciers (mean error −8 m with standard deviation 59 m) may reflect smaller altitude ranges for these glaciers such that there is less room for effects of different balance gradients to manifest themselves.


2019 ◽  
Vol 9 (1) ◽  
pp. 1-12
Author(s):  
Zbyněk Engel ◽  
Filip Hrbáček ◽  
Kamil Láska ◽  
Daniel Nývlt ◽  
Zdeněk Stachoň

This study presents surface mass balance of two small glaciers on James Ross Island calculated using constant and zonally-variable conversion factors. The density of 500 and 900 kg·m–3 adopted for snow in the accumulation area and ice in the ablation area, respectively, provides lower mass balance values that better fit to the glaciological records from glaciers on Vega Island and South Shetland Islands. The difference between the cumulative surface mass balance values based on constant (1.23 ± 0.44 m w.e.) and zonally-variable density (0.57 ± 0.67 m w.e.) is higher for Whisky Glacier where a total mass gain was observed over the period 2009–2015. The cumulative surface mass balance values are 0.46 ± 0.36 and 0.11 ± 0.37 m w.e. for Davies Dome, which experienced lower mass gain over the same period. The conversion approach does not affect much the spatial distribution of surface mass balance on glaciers, equilibrium line altitude and accumulation-area ratio. The pattern of the surface mass balance is almost identical in the ablation zone and very similar in the accumulation zone, where the constant conversion factor yields higher surface mass balance values. The equilibrium line altitude and accumulation-area ratio determined for the investigated glaciers differ by less than 2m and 0.01, respectively. The annual changes of equilibrium line altitude and the mean values determined over the period 2009–2015 for Whisky Glacier (311 ± 16 m a.s.l.) and Davies Dome (393 ± 18 m a.s.l.) coincide with the values reported from Bahía del Diablo Glacier on Vega Island but differ from the glaciological records on South Shetland Islands.


2021 ◽  
Author(s):  
Fredrik Boberg ◽  
Ruth Mottram ◽  
Nicolaj Hansen ◽  
Shuting Yang ◽  
Peter L. Langen

Abstract. The future rates of ice sheet melt in Greenland and Antarctica are an important factor when making estimates of the likely rate of sea level rise. Global climate models that took part in the fifth Coupled Model Intercomparison Project (CMIP5) have generally been unable to replicate observed rates of ice sheet melt. With the advent of the sixth Coupled Model Intercomparison Project (CMIP6), with a general increase in the equilibrium climate sensitivity, we here compare two versions of the global climate model EC-Earth using the regional climate model HIRHAM5 downscaling EC-Earth for Greenland and Antarctica. One version (v2) of EC-Earth is taken from CMIP5 for the high-emissions Representative Concentration Pathways (RCP8.5) scenario and the other (v3) from CMIP6 for the comparable high-emissions Shared Socioeconomic Pathways (SSP5-8.5) scenario). For Greenland, we downscale the two versions of EC-Earth for the historical period 1991–2010 and for the scenario period 2081–2100. For Antarctica, the periods are 1971–2000 and 2071–2100, respectively. For the Greenland Ice Sheet, we find that the mean change in temperature is 5.9 °C when downscaling EC-Earth v2 and 6.8 °C when downscaling EC-Earth v3. Corresponding values for Antarctica are 4.1 °C for v2 and 4.8 °C for v3. The mean change in surface mass balance at the end of the century under these high emissions scenarios is found to be −210 Gt yr−1 (v2) and −1150 Gt yr−1 (v3) for Greenland and 420 Gt yr−1 (v2) and 80 Gt yr−1 (v3) for Antarctica. These distinct differences in temperature change and particularly surface mass balance change are a result of the higher equilibrium climate sensitivity in EC-Earth v3 (4.3 K) compared with 3.3 K in EC-Earth v2 and the differences in greenhouse gas concentrations between the RCP8.5 and the SSP5-8.5 scenarios.


2015 ◽  
Vol 9 (6) ◽  
pp. 2135-2148 ◽  
Author(s):  
R. J. Braithwaite

Abstract. Translated into modern terminology, Kurowski suggested in 1891 that the equilibrium line altitude (ELA) of a glacier is equal to the mean altitude of the glacier when the whole glacier is in balance between accumulation and ablation. Kurowski's method has been widely misunderstood, partly due to inappropriate use of statistical terminology by later workers, and has only been tested by Braithwaite and Müller in a 1980 paper (for 32 glaciers). I now compare Kurowski's mean altitude with balanced-budget ELA calculated for 103 present-day glaciers with measured surface mass-balance data. Kurowski's mean altitude is significantly higher (at 95 % level) than balanced-budget ELA for 19 outlet and 42 valley glaciers, but not significantly higher for 34 mountain glaciers. The error in Kurowski mean altitude as a predictor of balanced-budget ELA might be due to generally lower balance gradients in accumulation areas compared with ablation areas for many glaciers, as suggested by several workers, but some glaciers have higher gradients, presumably due to precipitation increase with altitude. The relatively close agreement between balanced-budget ELA and mean altitude for mountain glaciers (mean error – 8 m with standard deviation 59 m) may reflect smaller altitude ranges for these glaciers such that there is less room for effects of different balance gradients to manifest themselves.


2020 ◽  
Author(s):  
Fredrik Boberg ◽  
Ruth Mottram ◽  
Nicolaj Hansen ◽  
Shuting Yang ◽  
Peter L. Langen

Abstract. The future rates of ice sheet melt in Greenland and Antarctica are an important factor when making estimates of the likely rate of sea level rise. Global climate models that took part in the fifth Coupled Model Intercomparison Project (CMIP5) have generally been unable to replicate observed rates of ice sheet melt. With the advent of the sixth Coupled Model Intercomparison Project (CMIP6), with a general increase in the equilibrium climate sensitivity, we here compare two versions of the global climate model EC-Earth using the regional climate model HIRHAM5 downscaling EC-Earth for Greenland and Antarctica. One version (v2) of EC-Earth is taken from CMIP5 for the high-emissions Representative Concentration Pathways (RCP8.5) scenario and the other (v3) from CMIP6 for the comparable high-emissions Shared Socioeconomic Pathways (SSP5-8.5) scenario). For Greenland, we downscale the two versions of EC-Earth for the historical period 1991–2010 and for the scenario period 2081–2100. For Antarctica, the periods are 1971–2000 and 2071–2100, respectively. For the Greenland Ice Sheet, we find that the mean change in temperature is 5.9 °C when downscaling EC-Earth v2 and 6.8 °C when downscaling EC-Earth v3. Corresponding values for Antarctica are 4.1 °C for v2 and 4.9 °C for v3. The mean change in surface mass balance at the end of the century under these high emissions scenarios is found to be −210 Gt yr−1 (v2) and −1150 Gt yr−1 (v3) for Greenland and +150 Gt yr−1 (v2) and −710 Gt yr−1 (v3) for Antarctica. These distinct differences in temperature change and particularly surface mass balance change are a result of the higher equilibrium climate sensitivity in EC-Earth v3 (4.3 K) compared with 3.3 K in EC-Earth v2 and the differences in greenhouse gas concentrations between the RCP8.5 and the SSP5-8.5 scenarios.


2008 ◽  
Vol 54 (185) ◽  
pp. 307-314 ◽  
Author(s):  
Antoine Rabatel ◽  
Jean-Pierre Dedieu ◽  
Emmanuel Thibert ◽  
Anne Letréguilly ◽  
Christian Vincent

AbstractAnnual equilibrium-line altitude (ELA) and surface mass balance of Glacier Blanc, Ecrins region, French Alps, were reconstructed from a 25 year time series of satellite images (1981–2005). The remote-sensing method used was based on identification of the snowline, which is easy to discern on optical satellite images taken at the end of the ablation season. In addition, surface mass balances at the ELA were reconstructed for the same period using meteorological data from three nearby weather stations. A comparison of the two types of series reveals a correlation of r > 0.67 at the 0.01 level of significance. Furthermore, the surface mass balances obtained from remote-sensing data are consistent with those obtained from field measurements on five other French glaciers (r = 0.76, p < 0.01). Also consistent for Glacier Blanc is the total mass loss (10.8 m w.e.) over the studied period. However, the surface mass balances obtained with the remote-sensing method show lower interannual variability. Given that the remote-sensing method is based on changes in the ELA, this difference probably results from the lower sensitivity of the surface mass balance to climate parameters at the ELA.


Sign in / Sign up

Export Citation Format

Share Document