Towards a Middle Pleistocene terrestrial climate reconstruction based on herpetofaunal assemblages from the Iberian Peninsula : state of the art and perspectives

This paper has been mainly funded by the University of Reading and also by projects CGL2016-80000-P (Spanish Ministry of Economy and Competitiveness) and 2017SGR-859 (Generalitat de Catalunya).The pattern of the varying climatic conditions in southern Europe over the last million years is well known from isotope studies on deep-ocean sediment cores and the long pollen records that have been produced for lacustrine and marine sedimentary sequences from Greece, Italy and the Iberian margin. However, although relative glacial and interglacial intensities are well studied, there are still few proxies that permit quantitative terrestrial temperature and precipitation reconstruction. In this context, fauna-based climate reconstructions based on evidence preserved in archaeological or palaeontological sites are of great interest, even if they only document short windows of that climate variability, because (a) they provide a range of temperature and precipitation estimates that are understandable in comparison with present climate; (b) they may allow the testing of predicted temperature changes under scenarios of future climate change; and (c) quantitative temperature and precipitation estimates for past glacials and interglacials for specific regions/latitudes can help to understand their effects on flora, fauna and hominids, as they are directly associated with those cultural and/or biological events. Moreover such reconstructions can bring further arguments to the discussion about important climatic events like the Mid-Bruhnes Event, a climatic transition between moderate warmths and greater warmths during interglacials. In this paper we review a decade of amphibian- and reptile-based climate reconstructions carried out for the Iberian Peninsula using the Mutual Ecogeographic Range method in order to present a regional synthesis from MIS 22 to MIS 6, discuss the climate pattern in relation to the Mid-Bruhnes Event and the thermal amplitude suggested by these estimates and finally to identify the chronological gaps that have still to be investigated.PostprintPeer reviewe

Land-ocean changes on orbital and millennial time scales and the penultimate glaciation

Past glacials can be thought of as natural experiments in which variations in boundary conditions influenced the character of climate change. However, beyond the last glacial, an integrated view of orbital- and millennial-scale changes and their relation to the record of glaciation has been lacking. Here, we present a detailed record of variations in the land-ocean system from the Portuguese margin during the penultimate glacial and place it within the framework of ice-volume changes, with particular reference to European ice-sheet dynamics. The interaction of orbital- and millennial-scale variability divides the glacial into an early part with warmer and wetter overall conditions and prominent climate oscillations, a transitional mid-part, and a late part with more subdued changes as the system entered a maximum glacial state. The most extreme event occurred in the mid-part and was associated with melting of the extensive European ice sheet and maximum discharge from the Fleuve Manche river. This led to disruption of the meridional overturning circulation, but not a major activation of the bipolar seesaw. In addition to stadial duration, magnitude of freshwater forcing, and background climate, the evidence also points to the influence of the location of freshwater discharges on the extent of interhemispheric heat transport

Extreme glacial cooling likely led to hominin depopulation of Europe in the Early Pleistocene

The oldest known hominin remains in Europe [~1.5 to ~1.1 million years ago (Ma)] have been recovered from Iberia, where paleoenvironmental reconstructions have indicated warm and wet interglacials and mild glacials, supporting the view that once established, hominin populations persisted continuously. We report analyses of marine and terrestrial proxies from a deep-sea core on the Portugese margin that show the presence of pronounced millennial-scale climate variability during a glacial period ~1.154 to ~1.123 Ma, culminating in a terminal stadial cooling comparable to the most extreme events of the last 400,000 years. Climate envelope–model simulations reveal a drastic decrease in early hominin habitat suitability around the Mediterranean during the terminal stadial. We suggest that these extreme conditions led to the depopulation of Europe, perhaps lasting for several successive glacial-interglacial cycles

Fast and slow components of interstadial warming in the North Atlantic during the last glacial

The abrupt nature of warming events recorded in Greenland ice-cores during the last glacial has generated much debate over their underlying mechanisms. Here, we present joint marine and terrestrial analyses from the Portuguese Margin, showing a succession of cold stadials and warm interstadials over the interval 35–57 ka. Heinrich stadials 4 and 5 contain considerable structure, with a short transitional phase leading to an interval of maximum cooling and aridity, followed by slowly increasing sea-surface temperatures and moisture availability. A climate model experiment reproduces the changes in western Iberia during the final part of Heinrich stadial 4 as a result of the gradual recovery of the Atlantic meridional overturning circulation. What emerges is that Greenland ice-core records do not provide a unique template for warming events, which involved the operation of both fast and slow components of the coupled atmosphere–ocean–sea-ice system, producing adjustments over a range of timescales

Fast and slow components of interstadial warming in the North Atlantic during the last glacial

Abstract: The abrupt nature of warming events recorded in Greenland ice-cores during the last glacial has generated much debate over their underlying mechanisms. Here, we present joint marine and terrestrial analyses from the Portuguese Margin, showing a succession of cold stadials and warm interstadials over the interval 35–57 ka. Heinrich stadials 4 and 5 contain considerable structure, with a short transitional phase leading to an interval of maximum cooling and aridity, followed by slowly increasing sea-surface temperatures and moisture availability. A climate model experiment reproduces the changes in western Iberia during the final part of Heinrich stadial 4 as a result of the gradual recovery of the Atlantic meridional overturning circulation. What emerges is that Greenland ice-core records do not provide a unique template for warming events, which involved the operation of both fast and slow components of the coupled atmosphere–ocean–sea-ice system, producing adjustments over a range of timescales

A 1.5-million-year record of orbital and millennial climate variability in the North Atlantic

Climate during the last glacial period was marked by abrupt instability on millennial timescales that included large swings of temperature in and around Greenland (Daansgard-Oeschger events) and smaller, more gradual changes in Antarctica (AIM events). Less is known about the existence and nature of similar variability during older glacial periods, especially during the early Pleistocene when glacial cycles were dominantly occurring at 41 kyr intervals compared to the much longer and deeper glaciations of the more recent period. Here, we report a continuous millennially resolved record of stable isotopes of planktic and benthic foraminifera at IODP Site U1385 (the "Shackleton Site") from the southwestern Iberian margin for the last 1.5 million years, which includes the Middle Pleistocene Transition (MPT). Our results demonstrate that millennial climate variability (MCV) was a persistent feature of glacial climate, both before and after the MPT. Prior to 1.2 Ma in the early Pleistocene, the amplitude of MCV was modulated by the 41 kyr obliquity cycle and increased when axial tilt dropped below 23.5° and benthic δ18O exceeded ∼3.8 ‰ (corrected to Uvigerina), indicating a threshold response to orbital forcing. Afterwards, MCV became focused mainly on the transitions into and out of glacial states (i.e. inceptions and terminations) and during times of intermediate ice volume. After 1.2 Ma, obliquity continued to play a role in modulating the amplitude of MCV, especially during times of glacial inceptions, which are always associated with declining obliquity. A non-linear role for obliquity is also indicated by the appearance of multiples (82, 123 kyr) and combination tones (28 kyr) of the 41 kyr cycle. Near the end of the MPT (∼0.65 Ma), obliquity modulation of MCV amplitude wanes as quasi-periodic 100 kyr and precession power increase, coinciding with the growth of oversized ice sheets on North America and the appearance of Heinrich layers in North Atlantic sediments. Whereas the planktic δ18O of Site U1385 shows a strong resemblance to Greenland temperature and atmospheric methane (i.e. Northern Hemisphere climate), millennial changes in benthic δ18O closely follow the temperature history of Antarctica for the past 800 kyr. The phasing of millennial planktic and benthic δ18O variation is similar to that observed for MIS 3 throughout much of the record, which has been suggested to mimic the signature of the bipolar seesaw - i.e. an interhemispheric asymmetry between the timing of cooling in Antarctica and warming in Greenland. The Iberian margin isotopic record suggests that bipolar asymmetry was a robust feature of interhemispheric glacial climate variations for at least the past 1.5 Ma despite changing glacial boundary conditions. A strong correlation exists between millennial increases in planktic δ18O (cooling) and decreases in benthic δ13C, indicating that millennial variations in North Atlantic surface temperature are mirrored by changes in deep-water circulation and remineralization of carbon in the abyssal ocean. We find strong evidence that climate variability on millennial and orbital scales is coupled across different timescales and interacts in both directions, which may be important for linking internal climate dynamics and external astronomical forcing

Extreme glacial cooling likely led to hominin depopulation of Europe in the Early Pleistocene

The oldest known hominin remains in Europe [~1.5 to ~1.1 million years ago (Ma)] have been recovered from Iberia, where paleoenvironmental reconstructions have indicated warm and wet interglacials and mild glacials, supporting the view that once established, hominin populations persisted continuously. We report analyses of marine and terrestrial proxies from a deep-sea core on the Portugese margin that show the presence of pronounced millennial-scale climate variability during a glacial period ~1.154 to ~1.123 Ma, culminating in a terminal stadial cooling comparable to the most extreme events of the last 400,000 years. Climate envelope-model simulations reveal a drastic decrease in early hominin habitat suitability around the Mediterranean during the terminal stadial. We suggest that these extreme conditions led to the depopulation of Europe, perhaps lasting for several successive glacial-interglacial cycles.We gratefully acknowledge financial support from The Leverhulme Trust grant RPG-2014-417 (P.C.T., V.M., and D.A.H.); the Catalan Government, Research Group 2021SGR00986 (J.O.G.); IBS, South Korea grant IBS-R028-D1 (A.T., K.-S.Y., and H.K.); the Human Origins Research Fund (C.B.S.); and the Calleva Foundation (C.B.S., S.A.P., and N.M.A.).Peer reviewe

A reference time scale for Site U1385 (Shackleton Site) on the SW Iberian Margin

We produced a composite depth scale and chronology for Site U1385 on the SW Iberian Margin. Using log(Ca/Ti) measured by core scanning XRF at 1-cm resolution in all holes, a composite section was constructed to 166.5 meters composite depth (mcd) that corrects for stretching and squeezing in each core. Oxygen isotopes of benthic foraminifera were correlated to a stacked δ^18O reference signal (LR04) to produce an oxygen isotope stratigraphy and age model. Variations in sediment color contain very strong precession signals at Site U1385, and the amplitude modulation of these cycles provides a powerful tool for developing an orbitally-tuned age model. We tuned the U1385 record by correlating peaks in L* to the local summer insolation maxima at 37^oN. The benthic δ^18O record of Site U1385, when placed on the tuned age model, generally agrees with other time scales within their respective chronologic uncertainties. The age model is transferred to down-core data to produce a continuous time series of log(Ca/Ti) that reflect relative changes of biogenic carbonate and detrital sediment. Biogenic carbonate increases during interglacial and interstadial climate states and decreases during glacial and stadial periods. Much of the variance in the log(Ca/Ti) is explained by a linear combination of orbital frequencies (precession, tilt and eccentricity), whereas the residual signal reflects suborbital climate variability. The strong correlation between suborbital log(Ca/Ti) variability and Greenland temperature over the last glacial cycle at Site U1385 suggests that this signal can be used as a proxy for millennial-scale climate variability over the past 1.5 Ma. Millennial climate variability, as expressed by log(Ca/Ti) at Site U1385, was a persistent feature of glacial climates over the past 1.5 Ma, including glacial periods of the early Pleistocene (‘41-kyr world’) when boundary conditions differed significantly from those of the late Pleistocene (‘100-kyr world’). Suborbital variability was suppressed during interglacial stages and enhanced during glacial periods, especially when benthic δ^18O surpassed ~ 3.3-3.5‰. Each glacial inception was marked by appearance of strong millennial variability and each deglaciation was preceded by a terminal stadial event. Suborbital variability may be a symptomatic feature of glacial climate or, alternatively, may play a more active role in the inception and/or termination of glacial cycles.This research was supported by the Natural Environmental Research Council Grant NE/K005804/1 to DH and LS and NE/J017922/1 to DH.This is the final version of the article. It first appeared from Elsevier via http://dx.doi.org/10.1016/j.gloplacha.2015.07.00

A reference time scale for Site U1385 (Shackleton Site) on the SW Iberian Margin

Variations in sediment color contain very strong precession signals at Site U1385, and the amplitude modulation of these cycles provides a powerful tool for developing an orbitally-tuned age model. We tuned the U1385 record by correlating peaks in L* to the local summer insolation maxima at 37°N. The benthic δ18O record of Site U1385, when placed on the tuned age model, generally agrees with other time scales within their respective chronologic uncertainties. The age model is transferred to down-core data to produce a continuous time series of log(Ca/Ti) that reflect relative changes of biogenic carbonate and detrital sediment. Biogenic carbonate increases during interglacial and interstadial climate states and decreases during glacial and stadial periods. Much of the variance in the log(Ca/Ti) is explained by a linear combination of orbital frequencies (precession, tilt and eccentricity), whereas the residual signal reflects suborbital climate variability. The strong correlation between suborbital log(Ca/Ti) variability and Greenland temperature over the last glacial cycle at Site U1385 suggests that this signal can be used as a proxy for millennial-scale climate variability over the past 1.5 Ma. Millennial climate variability, as expressed by log(Ca/Ti) at Site U1385, was a persistent feature of glacial climates over the past 1.5 Ma, including glacial periods of the early Pleistocene (‘41-kyr world’) when boundary conditions differed significantly from those of the late Pleistocene (‘100-kyr world’). Suborbital variability was suppressed during interglacial stages and enhanced during glacial periods, especially when benthic δ18O surpassed ~ 3.3–3.5‰. Each glacial inception was marked by appearance of strong millennial variability and each deglaciation was preceded by a terminal stadial event. Suborbital variability may be a symptomatic feature of glacial climate or, alternatively, may play a more active role in the inception and/or termination of glacial cycles