Dynamic Hydrological Discharge and Lake Modelling for Coupled Climate Model Simulations of the Last Glacial Cycle

Author(s):  
Thomas Riddick ◽  
Victor Brovkin ◽  
Stefan Hagemann ◽  
Uwe Mikolajewicz

<p>The continually evolving large ice sheets present in the Northern Hemisphere during the last glacial cycle caused significant changes to river pathways both through directly blocking rivers and through glacial isostatic adjustment. Associated with these changing river pathways was the formation and evolution of large glacial lakes such as Lake Agassiz. Studies have shown this changing hydrology had a significant impact on the ocean circulation through changing the pattern of freshwater discharge into the oceans. A coupled Earth system model (ESM) simulation of the last glacial cycle thus requires a hydrological discharge and lake model that uses a set of river pathways and lakes that evolve with Earth's changing orography while being able to reproduce the known present-day river network given the present-day orography. Here, we present a method for dynamically modelling rivers and lakes by applying predefined corrections to an evolving fine-scale orography (accounting for the changing ice sheets and isostatic rebound) each time the river directions and lakes basins are recalculated. The corrected orography thus produced is then used to create a set of fine-scale river pathways and these are then upscaled to a coarser scale on which an existing present-day hydrological discharge model within the JSBACH land surface model simulates the river flow. The associated glacial lakes are delineated from the same corrected fine scale orography; lake inflow and outflow being linked to the river flow model.</p>

2018 ◽  
Author(s):  
Thomas Riddick ◽  
Victor Brovkin ◽  
Stefan Hagemann ◽  
Uwe Mikolajewicz

Abstract. The continually evolving large ice sheets present in the Northern Hemisphere during the last glacial cycle caused significant changes to river pathways both through directly blocking rivers and through glacial isostatic adjustment. Studies have shown these river pathway changes had a significant impact on the ocean circulation through changing the pattern of freshwater discharge into the oceans. A coupled Earth System Model (ESM) simulation of the last glacial cycle thus requires a hydrological discharge model that uses a set of river pathways that evolve with Earth's changing orography while being able to reproduce the known present-day river network given the present-day orography. Here we present a method for dynamically modelling river pathways that meets such requirements by applying pre-defined corrections to an evolving fine scale orography (accounting for the changing ice sheets and isostatic rebound) each time the river directions are recalculated. The corrected orography thus produced is then used to create a set of fine scale river pathways and these are then upscaled to a coarser scale on which an existing present-day hydrological discharge model within the JSBACH land surface model simulates the river flow. Tests show that this procedure reproduces the known present-day river network to a sufficient degree of accuracy and is able to simulate plausible paleo river networks. It has also been shown this procedure can be run successfully multiple times as part of a transient coupled climate model simulation.


2018 ◽  
Vol 11 (10) ◽  
pp. 4291-4316 ◽  
Author(s):  
Thomas Riddick ◽  
Victor Brovkin ◽  
Stefan Hagemann ◽  
Uwe Mikolajewicz

Abstract. The continually evolving large ice sheets present in the Northern Hemisphere during the last glacial cycle caused significant changes to river pathways both through directly blocking rivers and through glacial isostatic adjustment. Studies have shown these river pathway changes had a significant impact on the ocean circulation through changing the pattern of freshwater discharge into the oceans. A coupled Earth system model (ESM) simulation of the last glacial cycle thus requires a hydrological discharge model that uses a set of river pathways that evolve with Earth's changing orography while being able to reproduce the known present-day river network given the present-day orography. Here, we present a method for dynamically modelling river pathways that meets such requirements by applying predefined corrections to an evolving fine-scale orography (accounting for the changing ice sheets and isostatic rebound) each time the river directions are recalculated. The corrected orography thus produced is then used to create a set of fine-scale river pathways and these are then upscaled to a coarser scale on which an existing present-day hydrological discharge model within the JSBACH land surface model simulates the river flow. Tests show that this procedure reproduces the known present-day river network to a sufficient degree of accuracy and is able to simulate plausible paleo-river networks. It has also been shown this procedure can be run successfully multiple times as part of a transient coupled climate model simulation.


2009 ◽  
Vol 5 (3) ◽  
pp. 329-345 ◽  
Author(s):  
S. Bonelli ◽  
S. Charbit ◽  
M. Kageyama ◽  
M.-N. Woillez ◽  
G. Ramstein ◽  
...  

Abstract. A 2.5-dimensional climate model of intermediate complexity, CLIMBER-2, fully coupled with the GREMLINS 3-D thermo-mechanical ice sheet model is used to simulate the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle and to investigate the ice sheets responses to both insolation and atmospheric CO2 concentration. This model reproduces the main phases of advance and retreat of Northern Hemisphere ice sheets during the last glacial cycle, although the amplitude of these variations is less pronounced than those based on sea level reconstructions. At the last glacial maximum, the simulated ice volume is 52.5×1015 m3 and the spatial distribution of both the American and Eurasian ice complexes is in reasonable agreement with observations, with the exception of the marine parts of these former ice sheets. A set of sensitivity studies has also been performed to assess the sensitivity of the Northern Hemisphere ice sheets to both insolation and atmospheric CO2. Our results suggest that the decrease of summer insolation is the main factor responsible for the early build up of the North American ice sheet around 120 kyr BP, in agreement with benthic foraminifera δ18O signals. In contrast, low insolation and low atmospheric CO2 concentration are both necessary to trigger a long-lasting glaciation over Eurasia.


2012 ◽  
Vol 6 (6) ◽  
pp. 4897-4938 ◽  
Author(s):  
S. Charbit ◽  
C. Dumas ◽  
M. Kageyama ◽  
D. M. Roche ◽  
C. Ritz

Abstract. Since the original formulation of the positive-degree-day (PDD) method, different PDD calibrations have been proposed in the literature in response to the increasing number of observations. Although these formulations provide a satisfactory description of the present-day Greenland geometry, they have not all been tested for paleo ice sheets. Using the climate-ice sheet model CLIMBER-GRISLI coupled with different PDD models, we evaluate how the parameterization of the ablation may affect the evolution of Northern Hemisphere ice sheets in the transient simulations of the last glacial cycle. Results from fully coupled simulations are compared to time-slice experiments carried out at different key periods of the last glacial period. We find large differences in the simulated ice sheets according to the chosen PDD model. These differences occur as soon as the onset of glaciation, therefore affecting the subsequent evolution of the ice system. To further investigate how the PDD method controls this evolution, special attention is given to the role of each PDD parameter. We show that glacial inception is critically dependent on the representation of the impact of the temperature variability from the daily to the inter-annual time scale, whose effect is modulated by the refreezing scheme. Finally, an additional set of sensitivity experiments has been carried out to assess the relative importance of melt processes with respect to initial ice sheet configuration in the construction and the evolution of past Northern Hemisphere ice sheets. Our analysis reveals that the impacts of the initial ice sheet condition may range from quite negligible to explaining about half of the LGM ice volume depending on the representation of stochastic temperature variations which remain the main driver of the evolution of the ice system.


1995 ◽  
Vol 21 ◽  
pp. 103-110 ◽  
Author(s):  
G. S. Boulton ◽  
N. Hulton ◽  
M. Vautravers

A numerical model is used to simulate ice-sheet behaviour in Europe through the last glacial cycle. It is used in two modes: a forward mode, in which the model is driven by a proxy palaeoclimate record and the output compared with a geological reconstruction of ice-sheet fluctuation; and an inverse mode, in which we determine the climate function that would be required to simulate geologically reconstructed ice-sheet fluctuations. From these simulations it is concluded that extra-glacial climates may be poor predictors of ice-sheet surface climates, and that climatic transitions during the glacial period may have been much more rapid and the intensity of warming during the early Holocene much greater than hitherto supposed. Stronger climate forcing is required to drive ice-sheet expansion when sliding occurs at the bed compared with a non-sliding bed. Sliding ice sheets grow more slowly and decay more rapidly than non-sliding ice sheets with the same climate forcing.


2019 ◽  
Vol 65 (252) ◽  
pp. 645-661 ◽  
Author(s):  
LU NIU ◽  
GERRIT LOHMANN ◽  
SEBASTIAN HINCK ◽  
EVAN J. GOWAN ◽  
UTA KREBS-KANZOW

ABSTRACTThe evolution of Northern Hemisphere ice sheets through the last glacial cycle is simulated with the glacial index method by using the climate forcing from one General Circulation Model, COSMOS. By comparing the simulated results to geological reconstructions, we first show that the modelled climate is capable of capturing the main features of the ice-sheet evolution. However, large deviations exist, likely due to the absence of nonlinear interactions between ice sheet and other climate components. The model uncertainties of the climate forcing are examined using the output from nine climate models from the Paleoclimate Modelling Intercomparison Project Phase III. The results show a large variability in simulated ice sheets between the different models. We find that the ice-sheet extent pattern resembles summer surface air temperature pattern at the Last Glacial Maximum, confirming the dominant role of surface ablation process for high-latitude Northern Hemisphere ice sheets. This study shows the importance of the upper boundary condition for ice-sheet modelling, and implies that careful constraints on climate output is essential for simulating realistic glacial Northern Hemisphere ice sheets.


2014 ◽  
Vol 10 (4) ◽  
pp. 1453-1471 ◽  
Author(s):  
M. Löfverström ◽  
R. Caballero ◽  
J. Nilsson ◽  
J. Kleman

Abstract. We present modelling results of the atmospheric circulation at the cold periods of marine isotope stage 5b (MIS 5b), MIS 4 and the Last Glacial Maximum (LGM), as well as the interglacial. The palaeosimulations are forced by ice-sheet reconstructions consistent with geological evidence and by appropriate insolation and greenhouse gas concentrations. The results suggest that the large-scale atmospheric winter circulation remained largely similar to the interglacial for a significant part of the glacial cycle. The proposed explanation is that the ice sheets were located in areas where their interaction with the mean flow is limited. However, the LGM Laurentide Ice Sheet induces a much larger planetary wave that leads to a zonalisation of the Atlantic jet. In summer, the ice-sheet topography dynamically induces warm temperatures in Alaska and central Asia that inhibits the expansion of the ice sheets into these regions. The warm temperatures may also serve as an explanation for westward propagation of the Eurasian Ice Sheet from MIS 4 to the LGM.


2010 ◽  
Vol 6 (2) ◽  
pp. 229-244 ◽  
Author(s):  
A. Ganopolski ◽  
R. Calov ◽  
M. Claussen

Abstract. 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.


Sign in / Sign up

Export Citation Format

Share Document