RESUMEN
The study of palaeoclimates relies on information sampled in natural archives such as deep sea cores. Scientific investigations often use such information in multi-stage analyses, typically with an age model being fitted to a core to convert depths into ages at stage one. These age estimates are then used as inputs to develop, calibrate or select climate models in a second stage of analysis. Here, we show that such multi-stage approaches can lead to misleading conclusions, and develop a joint inferential approach for climate reconstruction, model calibration and age estimation. As an illustration, we investigate the glacial-interglacial cycle, fitting both an age model and dynamical climate model to two benthic sediment cores spanning the past 780 kyr. To show the danger of a multi-stage analysis, we sample ages from the posterior distribution, then perform model selection conditional on the sampled age estimates, mimicking standard practice. Doing so repeatedly for different samples leads to model selection conclusions that are substantially different from each other, and from the joint inferential analysis. We conclude that multi-stage analyses are insufficient when dealing with uncertainty, and that to draw sound conclusions the full joint inferential analysis should be performed.
RESUMEN
The pacing of glacial-interglacial cycles during the Quaternary period (the past 2.6 million years) is attributed to astronomically driven changes in high-latitude insolation. However, it has not been clear how astronomical forcing translates into the observed sequence of interglacials. Here we show that before one million years ago interglacials occurred when the energy related to summer insolation exceeded a simple threshold, about every 41,000 years. Over the past one million years, fewer of these insolation peaks resulted in deglaciation (that is, more insolation peaks were 'skipped'), implying that the energy threshold for deglaciation had risen, which led to longer glacials. However, as a glacial lengthens, the energy needed for deglaciation decreases. A statistical model that combines these observations correctly predicts every complete deglaciation of the past million years and shows that the sequence of interglacials that has occurred is one of a small set of possibilities. The model accounts for the dominance of obliquity-paced glacial-interglacial cycles early in the Quaternary and for the change in their frequency about one million years ago. We propose that the appearance of larger ice sheets over the past million years was a consequence of an increase in the deglaciation threshold and in the number of skipped insolation peaks.
RESUMEN
Large proglacial lakes cool regional summer climate because of their large heat capacity, and have been shown to modify precipitation through mesoscale atmospheric feedbacks, as in the case of Lake Agassiz. Several large ice-dammed lakes, with a combined area twice that of the Caspian Sea, were formed in northern Eurasia about 90,000 years ago, during the last glacial period when an ice sheet centred over the Barents and Kara seas blocked the large northbound Russian rivers. Here we present high-resolution simulations with an atmospheric general circulation model that explicitly simulates the surface mass balance of the ice sheet. We show that the main influence of the Eurasian proglacial lakes was a significant reduction of ice sheet melting at the southern margin of the Barents-Kara ice sheet through strong regional summer cooling over large parts of Russia. In our simulations, the summer melt reduction clearly outweighs lake-induced decreases in moisture and hence snowfall, such as has been reported earlier for Lake Agassiz. We conclude that the summer cooling mechanism from proglacial lakes accelerated ice sheet growth and delayed ice sheet decay in Eurasia and probably also in North America.
