ALBERT

Description: In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. These fields are used to force a high resolution ice sheet model through the annual mean surface temperature and mass balance. The melt component of the mass balance is computed here using both a conventional positive degree day approach and a more physically-based alternative. As a validation of the model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography of Greenland and compared them with empirical data. For the present-day climate, simulated surface boundary conditions for the GIS do not differ significantly from those of a simple parameterization used in many previous simulations. However, for a prescribed, ice-free state, the differences in simulated climatology and surface mass balance between the regional energy-moisture balance model and the conventional approach become significant, with our model showing much stronger summer warming. When coupled to a high resolution, three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow sufficiently realistic simulations of the present-day GIS. However, when starting from ice-free conditions, two different equilibrium states are achieved, depending on the choice of melt scheme.

Description: In this paper, we propose a new sub-grid scale parameterization for the ice discharge into the ocean through outlet glaciers and inspect the role of different observational and palaeo constraints for the choice of an optimal set of model parameters. This parameterization was introduced into the polythermal ice-sheet model SICOPOLIS, which is coupled to the regional climate model of intermediate complexity REMBO. Using the coupled model, we performed large ensemble simulations over the last two glacial cycles. We exploit two major parameters: a melt parameter in the surface melt scheme of REMBO and an ice discharge parameter in our parameterization of ice discharge. Our constraints are the present-day Greenland ice sheet surface elevation, surface mass balance partition (ratio between ice discharge and total precipitation) and the Eemian interglacial elevation drop relative to present-day in the vicinity of the NEEM ice core. We show that the ice discharge parameterization enables us to simulate both the correct ice-sheet shape and mass balance partition at the same time without explicitly resolving the Greenland outlet glaciers. For model verification, we compare simulated total and sectoral ice discharge with those from other findings, including observations. For the model versions, which are inside the range of observational and palaeo constraints, our simulated Greenland ice sheet contribution to Eemian sea level rise relative to present-day amounts to 1.4 m on average (in the range of 0.6 and 2.5 m).

Description: In this paper, we propose a new sub-grid scale parameterization for the ice discharge into the ocean through outlet glaciers and inspect the role of different observational and palaeo constraints for the choice of an optimal set of model parameters. This parameterization was introduced into the polythermal ice-sheet model SICOPOLIS, which is coupled to the regional climate model of intermediate complexity REMBO. Using the coupled model, we performed large ensemble simulations over the last two glacial cycles by varying two major parameters: a melt parameter in the surface melt scheme of REMBO and a discharge scaling parameter in our parameterization of ice discharge. Our empirical constraints are the present-day Greenland ice sheet surface elevation, the surface mass balance partition (ratio between total ice discharge and total precipitation) and the Eemian interglacial elevation drop relative to present day in the vicinity of the NEEM ice core. We show that the ice discharge parameterization enables us to simulate both the correct ice-sheet shape and mass balance partition at the same time without explicitly resolving the Greenland outlet glaciers. For model verification, we compare the simulated total and sectoral ice discharge with other estimates. For the model versions that are consistent with the range of observational and palaeo constraints, our simulated Greenland ice sheet contribution to Eemian sea-level rise relative to present-day amounts to 1.4 m on average (in the range of 0.6 and 2.5 m).

Description: Here we use a very simple conceptual model in an attempt to reduce essential parts of the complex nonlinearity of abrupt glacial climate changes (the so-called Dansgaard-Oeschger events) to a few simple principles, namely (i) the existence of two different climate states, (ii) a threshold process and (iii) an overshooting in the stability of the system at the start and the end of the events, which is followed by a millennial-scale relaxation. By comparison with a so-called Earth system model of intermediate complexity (CLIMBER-2), in which the events represent oscillations between two climate states corresponding to two fundamentally different modes of deep-water formation in the North Atlantic, we demonstrate that the conceptual model captures fundamental aspects of the nonlinearity of the events in that model. We use the conceptual model in order to reproduce and reanalyse nonlinear resonance mechanisms that were already suggested in order to explain the characteristic time scale of Dansgaard-Oeschger events. In doing so we identify a new form of stochastic resonance (i.e. an overshooting stochastic resonance) and provide the first explicitly reported manifestation of ghost resonance in a geosystem, i.e. of a mechanism which could be relevant for other systems with thresholds and with multiple states of operation. Our work enables us to explicitly simulate realistic probability measures of Dansgaard-Oeschger events (e.g. waiting time distributions, which are a prerequisite for statistical analyses on the regularity of the events by means of Monte-Carlo simulations). We thus think that our study is an important advance in order to develop more adequate methods to test the statistical significance and the origin of the proposed glacial 1470-year climate cycle.

Description: During glacial-interglacial cycles, atmospheric CO2 concentration varied by about 100 ppmv in amplitude. While testing mechanisms that have led to the low glacial CO2 level could be done in equilibrium model experiments, an ultimate goal is to explain CO2 changes in transient simulations through the complete glacial-interglacial cycle. The computationally efficient Earth System model of intermediate complexity CLIMBER-2 is used to simulate global biogeochemistry over the last glacial cycle (126 kyr). The physical core of the model (atmosphere, ocean, land and ice sheets) is driven by orbital changes and reconstructed radiative forcing from greenhouses gases, ice, and aeolian dust. The carbon cycle model is able to reproduce the main features of the CO2 changes: a 50 ppmv CO2 drop during glacial inception, a minimum concentration at the last glacial maximum 80 ppmv lower than the Holocene value, and an abrupt 60 ppmv CO2 rise during the deglaciation. The model deep ocean δ13C also resembles reconstructions from deep-sea cores. The main drivers of atmospheric CO2 evolve in time: changes in sea surface temperatures and in the volume of bottom water of southern origin control atmospheric CO2 during the glacial inception and deglaciation; changes in carbonate chemistry and marine biology are dominant during the first and second parts of the glacial cycle, respectively. These feedback mechanisms could also significantly impact the ultimate climate response to the anthropogenic perturbation.

Description: We present results from modelling of the mid-Pliocene warm period (3.3–3 million years ago) using the Earth system model of intermediate complexity CLIMBER-2 analysing the effect of changes in boundary conditions as well as of orbital forcing on climate. First we performed equilibrium experiments following the PlioMIP (Pliocene Model Intercomparison Project) protocol with a CO2 concentration of 405 ppm, reconstructed mid-Pliocene orography and vegetation and a present-day orbital configuration. Simulated global Pliocene warming is about 2.5 °C, fully consistent with results of atmosphere–ocean general circulation model simulations performed for the same modelling setup. A factor separation analysis attributes 1.5 °C warming to CO2, 0.3 °C to orography, 0.2 °C to ice sheets and 0.4 °C to vegetation. Transient simulations for the entire mid-Pliocene warm period with time-dependent orbital forcing as well as interactive ice sheets and vegetation give a global warming varying within the range 1.9–2.8 °C. Ice sheet and vegetation feedbacks in synergy act as amplifiers of the orbital forcing, transforming seasonal insolation variations into an annual mean temperature signal. The effect of orbital forcing is more significant at high latitudes, especially during boreal summer, when the warming over land varies in the wide range from 0 to 10 °C. The modelled ice-sheet extent and vegetation distribution also show significant temporal variations. Modelled and reconstructed data for Northern Hemisphere sea-surface temperatures and vegetation distribution show the best agreement if the reconstructions are assumed to be representative for the warmest periods during the orbital cycles. This suggests that low-resolution Pliocene palaeoclimate reconstructions can reflect not only the impact of increased CO2 concentrations and topography changes but also the effect of orbital forcing. Therefore, the climate (Earth system) sensitivity estimates from Pliocene reconstructions which do not account for the effect of orbital forcing can be biased toward high values.

Description: A new version of the Earth system model of intermediate complexity, CLIMBER-2, which includes the three-dimensional polythermal ice-sheet model SICOPOLIS, is used to simulate the last glacial cycle forced by variations of the Earth's orbital parameters and atmospheric concentration of major greenhouse gases. The climate and ice-sheet components of the model are coupled bi-directionally through a physically-based surface energy and mass balance interface. The model accounts for the time-dependent effect of aeolian dust on planetary and snow albedo. The model successfully simulates the temporal and spatial dynamics of the major Northern Hemisphere (NH) ice sheets, including rapid glacial inception and strong asymmetry between the ice-sheet growth phase and glacial termination. Spatial extent and elevation of the ice sheets during the last glacial maximum agree reasonably well with palaeoclimate reconstructions. A suite of sensitivity experiments demonstrates that simulated ice-sheet evolution during the last glacial cycle is very sensitive to some parameters of the surface energy and mass-balance interface and dust module. The possibility of a considerable acceleration of the climate ice-sheet model is discussed.

Description: Permafrost influences a number of processes which are relevant for local and global climate. For example, it is well known that permafrost plays an important role in global carbon and methane cycles. Less is known about the interaction between permafrost and ice sheets. In this study a permafrost module is included in the Earth system model CLIMBER-2, and the coupled Northern Hemisphere (NH) permafrost–ice-sheet evolution over the last glacial cycle is explored. The model performs generally well at reproducing present-day permafrost extent and thickness. Modeled permafrost thickness is sensitive to the values of ground porosity, thermal conductivity and geothermal heat flux. Permafrost extent at the Last Glacial Maximum (LGM) agrees well with reconstructions and previous modeling estimates. Present-day permafrost thickness is far from equilibrium over deep permafrost regions. Over central Siberia and the Arctic Archipelago permafrost is presently up to 200–500 m thicker than it would be at equilibrium. In these areas, present-day permafrost depth strongly depends on the past climate history and simulations indicate that deep permafrost has a memory of surface temperature variations going back to at least 800 ka. Over the last glacial cycle permafrost has a relatively modest impact on simulated NH ice sheet volume except at LGM, when including permafrost increases ice volume by about 15 m sea level equivalent in our model. This is explained by a delayed melting of the ice base from below by the geothermal heat flux when the ice sheet sits on a porous sediment layer and permafrost has to be melted first. Permafrost affects ice sheet dynamics only when ice extends over areas covered by thick sediments, which is the case at LGM.

Description: The impacts of various scenarios of snow and glaciers developing over the Tibetan Plateau on climate change in Afro-Asian monsoon region and other regions during the Holocene (9 kyr BP–0 kyr BP) are studied by using the coupled climate model of intermediate complexity, CLIMBER-2. The simulations show that the imposed snow and glaciers over the Tibetan Plateau in the mid-Holocene induce global summer temperature decreases, especially in the northern parts of Europe, Asia, and North America. At the same time, with the imposed snow and glaciers, summer precipitation decreases strongly in North Africa and South Asia as well as northeastern China, while it increases in Southeast Asia and the Mediterranean. For the whole period of Holocene (9 kyr BP–0 kyr BP), the response of vegetation cover to the imposed snow and glaciers cover over the Tibetan Plateau is not synchronous in South Asia and in North Africa, showing an earlier and a more rapid decrease in vegetation cover in North Africa from 9 to 6 kyr BP while it has only minor influence on that in South Asia until 5 kyr BP. Imposed gradually increased snow and glacier cover over the Tibetan Plateau causes temperature increases in South Asia and it decreases in North Africa and Southeast Asia during 6 kyr BP to 0 kyr BP. The precipitation decreases rapidly in North Africa and South Asia while it decreases slowly or unchanged during 6 kyr BP to 0 kyr BP with imposed snow and glacier cover over the Tibetan Plateau. The different scenarios of snow and glacier developing over the Tibetan Plateau would result in differences in variation of temperature, precipitation and vegetation cover in North Africa, South Asia and Southeast Asia. The model results show that the response of climate change in African-Asian monsoon region to snow and glacier cover over the Tibetan Plateau is in the way that the snow and glaciers amplify the effect of vegetation feedback and, hence, further amplify orbital forcing.