The Milankovitch theory and climate sensitivity: 2. Interaction between the northern hemisphere ice sheets and the climate system

1988 ◽  
Vol 93 (D9) ◽  
pp. 11175 ◽  
Author(s):  
Binyamin U. Neeman ◽  
George Ohring ◽  
Joachim H. Joseph
2021 ◽  
Author(s):  
Denis-Didier Rousseau ◽  
Witold Bagniewski ◽  
Michael Ghil

Abstract. Abrupt climate changes constitute a relatively new field of research, which addresses variations occurring in a relatively short time interval of tens to a hundred years. Such time scales do not correspond to the tens or hundreds of thousands of years that the astronomical theory of climate addresses. The latter theory involves parameters that are external to the climate system and whose multi-periodic variations are reliably known and almost constant for a large extent of Earth history. Abrupt changes, conversely, appear to involve fast processes that are internal to the climate system; these processes varied considerably during the past 2.6 Myr, and yielded more irregular fluctuations. In this paper, we re-examine the main climate variations determined from the U1308 North Atlantic marine record, which yields a detailed calving history of the Northern Hemisphere ice sheets over the past 3.2 Myr. The magnitude and periodicity of the ice-rafted debris (IRD) events observed in the U1308 record allow one to determine the timing of several abrupt climate changes, the larger ones corresponding to the massive iceberg discharges labeled Heinrich events (HEs). In parallel, abrupt warmings, called Dansgaard-Oeschger (DO) events, have been identified in the Greenland records of the last glaciation cycle. Combining the HE and DO observations, we study a complex mechanism that may lead to the observed millennial-scale variability corresponding to the abrupt climate changes of last 0.9 Myr. This mechanism relies on amended Bond cycles, which group DO events and the associated Greenland stadials into a trend of increased cooling, with IRD events embedded into every stadial, the latest of these being an HE. These Bond cycles may have occurred during the last 0.9 Ma when Northern Hemisphere ice sheets reached their maximum extent and volume, thus becoming a major player in this time interval’s climate dynamics. Since the waxing and waning of ice sheets during the Quaternary period are orbitally paced, we conclude that the abrupt climate changes observed during the Mid and Upper Pleistocene are therewith indirectly linked to the astronomical theory of climate.


1993 ◽  
Vol 39 (131) ◽  
pp. 45-49 ◽  
Author(s):  
André Berger ◽  
Hubert Gallée ◽  
Christian Tricot

Abstract A two-dimensional model which links the atmosphere, the mixed layer of the ocean, the sea ice, the continents, the ice sheets and their underlying bedrock has been used to test the Milankovitch theory over the last glacial—interglacial cycle. It was found that the orbital variations alone can induce, in such a system, feed-backs sufficient to generate the low-frequency part of the climatic variations over the last 122 kyear. These simulated variations at the astronomical time-scale are broadly in agreement with ice volume and sea-level reconstructions independently obtained from geological data. Imperfections in the simulated climate were the insufficient southward extent of the ice sheets and the too small hemispheric cooling during the last glacial maximum. These deficiencies were partly remedied in a further experiment (Gallée and others, in press) by using the time-dependent CO2 atmospheric concentration given by the Vostok ice core in addition to the astronomical forcing. For this second experiment, the main mechanisms and feedbacks responsible for the glaciation and the deglaciation in the model are discussed here.


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.


2005 ◽  
Vol 48 (4) ◽  
pp. 253-273 ◽  
Author(s):  
Sylvie Charbit ◽  
Masa Kageyama ◽  
Didier Roche ◽  
Catherine Ritz ◽  
Gilles Ramstein

2021 ◽  
Vol 17 (1) ◽  
pp. 253-267
Author(s):  
Jiang Zhu ◽  
Christopher J. Poulsen

Abstract. Equilibrium climate sensitivity (ECS) has been directly estimated using reconstructions of past climates that are different than today's. A challenge to this approach is that temperature proxies integrate over the timescales of the fast feedback processes (e.g., changes in water vapor, snow, and clouds) that are captured in ECS as well as the slower feedback processes (e.g., changes in ice sheets and ocean circulation) that are not. A way around this issue is to treat the slow feedbacks as climate forcings and independently account for their impact on global temperature. Here we conduct a suite of Last Glacial Maximum (LGM) simulations using the Community Earth System Model version 1.2 (CESM1.2) to quantify the forcing and efficacy of land ice sheets (LISs) and greenhouse gases (GHGs) in order to estimate ECS. Our forcing and efficacy quantification adopts the effective radiative forcing (ERF) and adjustment framework and provides a complete accounting for the radiative, topographic, and dynamical impacts of LIS on surface temperatures. ERF and efficacy of LGM LIS are −3.2 W m−2 and 1.1, respectively. The larger-than-unity efficacy is caused by the temperature changes over land and the Northern Hemisphere subtropical oceans which are relatively larger than those in response to a doubling of atmospheric CO2. The subtropical sea-surface temperature (SST) response is linked to LIS-induced wind changes and feedbacks in ocean–atmosphere coupling and clouds. ERF and efficacy of LGM GHG are −2.8 W m−2 and 0.9, respectively. The lower efficacy is primarily attributed to a smaller cloud feedback at colder temperatures. Our simulations further demonstrate that the direct ECS calculation using the forcing, efficacy, and temperature response in CESM1.2 overestimates the true value in the model by approximately 25 % due to the neglect of slow ocean dynamical feedback. This is supported by the greater cooling (6.8 ∘C) in a fully coupled LGM simulation than that (5.3 ∘C) in a slab ocean model simulation with ocean dynamics disabled. The majority (67 %) of the ocean dynamical feedback is attributed to dynamical changes in the Southern Ocean, where interactions between upper-ocean stratification, heat transport, and sea-ice cover are found to amplify the LGM cooling. Our study demonstrates the value of climate models in the quantification of climate forcings and the ocean dynamical feedback, which is necessary for an accurate direct ECS estimation.


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.


2016 ◽  
Author(s):  
David A. Hodell ◽  
James E.T. Channell

Abstract. We present a 3.2-Myr record of stable isotopes and physical properties at IODP Site U1308 (re-occupation of DSDP Site 609) located within the ice-rafted detritus (IRD) belt of the North Atlantic. We compare the isotope and lithological proxies at Site U1308 with other North Atlantic records (e.g., Sites 982, 607/U1313 and U1304) to reconstruct the history of orbital and millennial-scale climate variability during the Quaternary. The Site U1308 record documents a progressive increase in the intensity of Northern Hemisphere glacial-interglacial cycles during the late Pliocene and Quaternary with mode transitions at ~ 2.7, 1.5, 0.9 and 0.65 Ma. These transitions mark times of change in the growth and stability of Northern Hemisphere ice sheets. They also coincide with increases in vertical carbon isotope gradients between the intermediate and deep ocean, suggesting changes in deep carbon storage and atmospheric CO2. Orbital and millennial climate variability co-evolved during the Quaternary such that the trend towards larger ice sheets was accompanied by changes in the style, frequency and intensity of millennial-scale variability. This co-evolution may be important for explaining the observed patterns of Quaternary climate change.


2018 ◽  
Vol 4 (1/2) ◽  
pp. 19-36 ◽  
Author(s):  
Alex G. Libardoni ◽  
Chris E. Forest ◽  
Andrei P. Sokolov ◽  
Erwan Monier

Abstract. Historical time series of surface temperature and ocean heat content changes are commonly used metrics to diagnose climate change and estimate properties of the climate system. We show that recent trends, namely the slowing of surface temperature rise at the beginning of the 21st century and the acceleration of heat stored in the deep ocean, have a substantial impact on these estimates. Using the Massachusetts Institute of Technology Earth System Model (MESM), we vary three model parameters that influence the behavior of the climate system: effective climate sensitivity (ECS), the effective ocean diffusivity of heat anomalies by all mixing processes (Kv), and the net anthropogenic aerosol forcing scaling factor. Each model run is compared to observed changes in decadal mean surface temperature anomalies and the trend in global mean ocean heat content change to derive a joint probability distribution function for the model parameters. Marginal distributions for individual parameters are found by integrating over the other two parameters. To investigate how the inclusion of recent temperature changes affects our estimates, we systematically include additional data by choosing periods that end in 1990, 2000, and 2010. We find that estimates of ECS increase in response to rising global surface temperatures when data beyond 1990 are included, but due to the slowdown of surface temperature rise in the early 21st century, estimates when using data up to 2000 are greater than when data up to 2010 are used. We also show that estimates of Kv increase in response to the acceleration of heat stored in the ocean as data beyond 1990 are included. Further, we highlight how including spatial patterns of surface temperature change modifies the estimates. We show that including latitudinal structure in the climate change signal impacts properties with spatial dependence, namely the aerosol forcing pattern, more than properties defined for the global mean, climate sensitivity, and ocean diffusivity.


Sign in / Sign up

Export Citation Format

Share Document