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.

Clay hydroxyl isotopes show an enhanced hydrologic cycle during the Paleocene-Eocene Thermal Maximum

The Paleocene-Eocene Thermal Maximum (PETM) was an abrupt global warming event associated with a large injection of carbon into the ocean-atmosphere system, as evidenced by a diagnostic carbon isotope excursion (CIE). Evidence also suggests substantial hydrologic perturbations, but details have been hampered by a lack of appropriate proxies. To address this shortcoming, here we isolate and measure the isotopic composition of hydroxyl groups (OH^-) in clay minerals from a highly expanded PETM section in the North Sea Basin, together with their bulk oxygen isotope composition. At this location, we show that hydroxyl O- and H-isotopes are less influenced than bulk values by clay compositional changes due to mixing and/or inherited signals and thus better track hydrologic variability. We find that clay OH^- hydrogen-isotope values (δ²H_OH) decrease slowly prior to the PETM and then abruptly by ∼8‰ at the CIE onset. Coincident with an increase in relative kaolinite content, this indicates increased rainfall and weathering and implies an enhanced hydrologic cycle response to global warming, particularly during the early stages of the PETM. Subsequently, δ²H_OH returns to pre-PETM values well before the end of the CIE, suggesting hydrologic changes in the North Sea were short-lived relative to carbon-cycle perturbations.

Changes in North Atlantic deep-water oxygenation across the Middle Pleistocene Transition

The oxygen concentrations of oceanic deep-water and atmospheric carbon dioxide (pCO₂) are intrinsically linked through organic carbon remineralization and storage as dissolved inorganic carbon in the deep sea. We present a high-resolution reconstruction of relative changes in oxygen concentration in the deep North Atlantic for the past 1.5 million years using the carbon isotope gradient between epifaunal and infaunal benthic foraminifera species as a proxy for paleo-oxygen. We report a significant (>40 micromole per kilogram) reduction in glacial Atlantic deep-water oxygenation at ~960 thousand to 900 thousand years ago that coincided with increased continental ice volume and a major change in ocean thermohaline circulation. Paleo-oxygen results support a scenario of decreasing deep-water oxygen concentrations, increased respired carbon storage, and a reduction in glacial pCO₂ across the Middle Pleistocene Transition.

Continuous and simultaneous measurement of triple-oxygen and hydrogen isotopes of liquid and vapor during evaporation experiments

RATIONALE

Oxygen and hydrogen isotopes are important tools for studying the modern and past hydrological cycle. Previous evaporation experiments used episodic measurement of liquid and/or vapor or did not measure all isotopologues of water. Here, we describe an evaporation experimental system that allows all isotopologues of liquid and water vapor to be measured simultaneously and near-continuously at high precision using cavity ring-down laser spectroscopy (CRDS).

METHODS

Evaporating liquid is periodically sampled from a closed recirculating loop by a syringe pump that delivers a constant supply of water to the vaporizer, achieving a water vapor concentration of 20,000 ppmV H₂ O (±132, 1σ). Vapor is sampled directly from the evaporation chamber. Isotope ratios are measured simultaneously with a Picarro L2140-i CRDS instrument.

RESULTS

For liquid measurements, Allan variance analysis indicates an optimum data collection window of 34 min for oxygen isotopes and 27 min for hydrogen isotopes. During these periods, the mean standard error is ±0.0081‰ for δ¹⁷ O values, ±0.0081‰ for δ¹⁸ O values, and ±0.019‰ for δ² H values. For the derived parameters ¹⁷ O-excess and d-excess, the standard error of the mean is 5.8 per meg and 0.07‰, respectively. For the vapor phase a 12.5 min data window for all isotopologues results in a mean standard error of ±0.012‰ for δ¹⁷ O values, ±0.011‰ for δ¹⁸ O values, and ±0.023‰ for δ² H values. For the derived parameters, the standard error of the mean is 9.2 per meg for ¹⁷ O-excess and 0.099‰ for d-excess. These measurements result in consistently narrow 95% confidence limits for the slopes of ln(δ¹⁷ O + 1) vs ln(δ¹⁸ O + 1) and ln(δ² H + 1) vs ln(δ¹⁸ O + 1).

CONCLUSIONS

The experimental method permits measurement of fractionation of triple-oxygen and hydrogen isotopes of evaporating water under varying controlled conditions at high precision. Application of this method will be useful for testing theoretical models of evaporation and conducting experiments to simulate evaporation and isotopic equilibration in natural systems.

An astronomically dated record of Earth’s climate and its predictability over the last 66 million years

Much of our understanding of Earth’s past climate comes from the measurement of oxygen and carbon isotope variations in deep-sea benthic foraminifera. Yet, long intervals in existing records lack the temporal resolution and age control needed to thoroughly categorize climate states of the Cenozoic era and to study their dynamics. Here, we present a new, highly resolved, astronomically dated, continuous composite of benthic foraminifer isotope records developed in our laboratories. Four climate states—Hothouse, Warmhouse, Coolhouse, Icehouse—are identified on the basis of their distinctive response to astronomical forcing depending on greenhouse gas concentrations and polar ice sheet volume. Statistical analysis of the nonlinear behavior encoded in our record reveals the key role that polar ice volume plays in the predictability of Cenozoic climate dynamics.

Interglacial instability of North Atlantic Deep Water ventilation

Disrupting North Atlantic Deep Water (NADW) ventilation is a key concern in climate projections. We use (sub)centennially resolved bottom water δ¹³C records that span the interglacials of the last 0.5 million years to assess the frequency of and the climatic backgrounds capable of triggering large NADW reductions. Episodes of reduced NADW in the deep Atlantic, similar in magnitude to glacial events, have been relatively common and occasionally long-lasting features of interglacials. NADW reductions were triggered across the range of recent interglacial climate backgrounds, which demonstrates that catastrophic freshwater outburst floods were not a prerequisite for large perturbations. Our results argue that large NADW disruptions are more easily achieved than previously appreciated and that they occurred in past climate conditions similar to those we may soon face.

The residence time of Southern Ocean surface waters and the 100,000-year ice age cycle

From 1.25 million to 700,000 years ago, the ice age cycle deepened and lengthened from 41,000- to 100,000-year periodicity, a transition that remains unexplained. Using surface- and bottom-dwelling foraminifera from the Antarctic Zone of the Southern Ocean to reconstruct the deep-to-surface supply of water during the ice ages of the past 1.5 million years, we found that a reduction in deep water supply and a concomitant freshening of the surface ocean coincided with the emergence of the high-amplitude 100,000-year glacial cycle. We propose that this slowing of deep-to-surface circulation (i.e., a longer residence time for Antarctic surface waters) prolonged ice ages by allowing the Antarctic halocline to strengthen, which increased the resistance of the Antarctic upper water column to orbitally paced drivers of carbon dioxide release.

Transient hydrodynamic effects influence organic carbon signatures in marine sediments

Ocean dynamics served an important role during past dramatic climate changes via impacts on deep-ocean carbon storage. Such changes are recorded in sedimentary proxies of hydrographic change on continental margins, which lie at the ocean-atmosphere-earth interface. However, interpretations of these records are challenging, given complex interplays among processes delivering particulate material to and from ocean margins. Here we report radiocarbon (14C) signatures measured for organic carbon in differing grain-size sediment fractions and foraminifera in a sediment core retrieved from the southwest Iberian margin, spanning the last ~25,000 yr. Variable differences of 0-5000 yr in radiocarbon age are apparent between organic carbon in differing grain-sizes and foraminifera of the same sediment layer. The magnitude of 14C differences co-varies with key paleoceanographic indices (e.g., proximal bottom-current density gradients), which we interpret as evidence of Atlantic-Mediterranean seawater exchange influencing grain-size specific carbon accumulation and translocation. These findings underscore an important link between regional hydrodynamics and interpretations of down-core sedimentary proxies.

Quantification of drought during the collapse of the classic Maya civilization

The demise of Lowland Classic Maya civilization during the Terminal Classic Period (~800 to 1000 CE) is a well-cited example of how past climate may have affected ancient societies. Attempts to estimate the magnitude of hydrologic change, however, have met with equivocal success because of the qualitative and indirect nature of available climate proxy data. We reconstructed the past isotopic composition (δ¹⁸O, δD, ¹⁷O-excess, and d-excess) of water in Lake Chichancanab, Mexico, using a technique that involves isotopic analysis of the structurally bound water in sedimentary gypsum, which was deposited under drought conditions. The triple oxygen and hydrogen isotope data provide a direct measure of past changes in lake hydrology. We modeled the data and conclude that annual precipitation decreased between 41 and 54% (with intervals of up to 70% rainfall reduction during peak drought conditions) and that relative humidity declined by 2 to 7% compared to present-day conditions.

Reducing Uncertainties in Carbonate Clumped Isotope Analysis Through Consistent Carbonate-Based Standardization

About a decade after its introduction, the field of carbonate clumped isotope thermometry is rapidly expanding because of the large number of possible applications and its potential to solve long-standing questions in Earth Sciences. Major factors limiting the application of this method are the very high analytical precision required for meaningful interpretations, the relatively complex sample preparation procedures, and the mass spectrometric corrections needed. In this paper we first briefly review the evolution of the analytical and standardization procedures and discuss the major remaining sources of uncertainty. We propose that the use of carbonate standards to project the results to the carbon dioxide equilibrium scale can improve interlaboratory data comparability and help to solve long-standing discrepancies between laboratories and temperature calibrations. The use of carbonates reduces uncertainties related to gas preparation and cleaning procedures and ensures equal treatment of samples and standards. We present a set of carbonate standards of diverse composition, discuss how they can be used to correct for mass spectrometric biases, and demonstrate that their use significantly improves the comparability among four laboratories. We propose that the use of these standards or of a similar set of carbonate standards will improve the comparability of data across laboratories.

Online Differential Thermal Isotope Analysis of Hydration Water in Minerals by Cavity Ringdown Laser Spectroscopy

We have developed a new method for measuring the isotopic composition (δ¹⁸O and δD) of different types of bonded water (e.g., molecular water, hydroxyl) contained in hydrated minerals by coupling a thermal gravimeter (TG) and a cavity ringdown laser spectrometer (CRDS). The method involves precisely step-heating a mineral sample, allowing the separation of the different types of waters that are released at different temperatures. Simultaneously, the water vapor evolved from the mineral sample is analyzed for oxygen and hydrogen isotopes by CRDS. Isotopic values for the separate peaks are calculated by integrating the product of the water amounts and its isotopic values, after correcting for background. We provide examples of the application of the differential thermal isotope analysis (DTIA) method to a variety of hydrous minerals and mineraloids including gypsum, clays, and amorphous silica (opal). The isotopic compositions of the total water evolved from a set of natural gypsum samples by DTIA are compared with the results of a conventional offline water extraction method followed by CRDS analysis. The results from both methods are in excellent agreement, and precisions (1σ) for δ¹⁸O (±0.12‰) and δD (±0.8‰) of the total gypsum hydration water from the DTIA method are comparable to that obtained by the offline method. A range of analytical challenges and solutions (e.g., spectroscopic interferences produced by VOCs in natural samples, isotopic exchange with structural oxygen, etc.) are discussed. The DTIA method has wide ranging applications for addressing fundamental problems across many disciplines in earth and planetary sciences, including paleoclimatology, sedimentology, volcanology, water exchange between the solid earth and hydrosphere, and water on Mars and other planetary bodies.

Lead (Pb) Isotope Baselines for Studies of Ancient Human Migration and Trade in the Maya Region

We examined the potential use of lead (Pb) isotopes to source archaeological materials from the Maya region of Mesoamerica. The main objectives were to determine if: 1) geologic terrains throughout the Maya area exhibit distinct lead isotope ratios (²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb), and 2) a combination of lead and strontium ratios can enhance sourcing procedures in the Mesoamerica region. We analyzed 60 rock samples for lead isotope ratios and a representative subset of samples for lead, uranium, and thorium concentrations across the Maya region, including the Northern Lowlands of the Mexican Yucatan Peninsula, the Southern Lowlands of Guatemala and Belize, the Volcanic Highlands, the Belizean Maya Mountains, and the Metamorphic Province/Motagua Valley. Although there is some overlap within certain sub-regions, particularly the geologically diverse Metamorphic Province, lead isotopes can be used to distinguish between the Northern Lowlands, the Southern Lowlands, and the Volcanic Highlands. The distinct lead isotope ratios in the sub-regions are related to the geology of the Maya area, exhibiting a general trend in the lowlands of geologically younger rocks in the north to older rocks in the south, and Cenozoic volcanic rocks in the southern highlands. Combined with other sourcing techniques such as strontium (⁸⁷Sr/⁸⁶Sr) and oxygen (δ¹⁸O), a regional baseline for lead isotope ratios can contribute to the development of lead isoscapes in the Maya area, and may help to distinguish among geographic sub-regions at a finer scale than has been previously possible. These isotope baselines will provide archaeologists with an additional tool to track the origin and movement of ancient humans and artifacts across this important region.

Evolution of South Atlantic density and chemical stratification across the last deglaciation

Explanations of the glacial-interglacial variations in atmospheric pCO2 invoke a significant role for the deep ocean in the storage of CO2. Deep-ocean density stratification has been proposed as a mechanism to promote the storage of CO2 in the deep ocean during glacial times. A wealth of proxy data supports the presence of a "chemical divide" between intermediate and deep water in the glacial Atlantic Ocean, which indirectly points to an increase in deep-ocean density stratification. However, direct observational evidence of changes in the primary controls of ocean density stratification, i.e., temperature and salinity, remain scarce. Here, we use Mg/Ca-derived seawater temperature and salinity estimates determined from temperature-corrected δ(18)O measurements on the benthic foraminifer Uvigerina spp. from deep and intermediate water-depth marine sediment cores to reconstruct the changes in density of sub-Antarctic South Atlantic water masses over the last deglaciation (i.e., 22-2 ka before present). We find that a major breakdown in the physical density stratification significantly lags the breakdown of the deep-intermediate chemical divide, as indicated by the chemical tracers of benthic foraminifer δ(13)C and foraminifer/coral (14)C. Our results indicate that chemical destratification likely resulted in the first rise in atmospheric pCO2, whereas the density destratification of the deep South Atlantic lags the second rise in atmospheric pCO2 during the late deglacial period. Our findings emphasize that the physical and chemical destratification of the ocean are not as tightly coupled as generally assumed.

Simultaneous analysis of (17) O/(16) O, (18) O/(16) O and (2) H/(1) H of gypsum hydration water by cavity ring-down laser spectroscopy

RATIONALE

The recent development of cavity ring-down laser spectroscopy (CRDS) instruments capable of measuring (17) O-excess in water has created new opportunities for studying the hydrologic cycle. Here we apply this new method to studying the triple oxygen ((17) O/(16) O, (18) O/(16) O) and hydrogen ((2) H/(1) H) isotope ratios of gypsum hydration water (GHW), which can provide information about the conditions under which the mineral formed and subsequent post-depositional interaction with other fluids.

METHODS

We developed a semi-automated procedure for extracting GHW by slowly heating the sample to 400°C in vacuo and cryogenically trapping the evolved water. The isotopic composition (δ(17) O, δ(18) O and δ(2) H values) of the GHW is subsequently measured by CRDS. The extraction apparatus allows the dehydration of five samples and one standard simultaneously, thereby increasing the long-term precision and sample throughput compared with previous methods. The apparatus is also useful for distilling brines prior to isotopic analysis. A direct comparison is made between results of (17) O-excess in GHW obtained by CRDS and fluorination followed by isotope ratio mass spectrometry (IRMS) of O2 .

RESULTS

The long-term analytical precision of our method of extraction and isotopic analysis of GHW by CRDS is ±0.07‰ for δ(17) O values, ±0.13‰ for δ(18) O values and ±0.49‰ for δ(2) H values (all ±1SD), and ±1.1‰ and ±8 per meg for the deuterium-excess and (17) O-excess, respectively. Accurate measurement of the (17) O-excess values of GHW, of both synthetic and natural samples, requires the use of a micro-combustion module (MCM). This accessory removes contaminants (VOCs, H2 S, etc.) from the water vapour stream that interfere with the wavelengths used for spectroscopic measurement of water isotopologues. CRDS/MCM and IRMS methods yield similar isotopic results for the analysis of both synthetic and natural gypsum samples within analytical error of the two methods.

CONCLUSIONS

We demonstrate that precise and simultaneous isotopic measurements of δ(17) O, δ(18) O and δ(2) H values, and the derived deuterium-excess and (17) O-excess, can be obtained from GHW and brines using a new extraction apparatus and subsequent measurement by CRDS. This method provides new opportunities for the application of water isotope tracers in hydrologic and paleoclimatologic research.

A stagnation event in the deep South Atlantic during the last interglacial period

During the last interglacial period, global temperatures were ~2°C warmer than at present and sea level was 6 to 8 meters higher. Southern Ocean sediments reveal a spike in authigenic uranium 127,000 years ago, within the last interglacial, reflecting decreased oxygenation of deep water by Antarctic Bottom Water (AABW). Unlike ice age reductions in AABW, the interglacial stagnation event appears decoupled from open ocean conditions and may have resulted from coastal freshening due to mass loss from the Antarctic ice sheet. AABW reduction coincided with increased North Atlantic Deep Water (NADW) formation, and the subsequent reinvigoration in AABW coincided with reduced NADW formation. Thus, alternation of deep water formation between the Antarctic and the North Atlantic, believed to characterize ice ages, apparently also occurs in warm climates.

Iron fertilization of the Subantarctic ocean during the last ice age

John H. Martin, who discovered widespread iron limitation of ocean productivity, proposed that dust-borne iron fertilization of Southern Ocean phytoplankton caused the ice age reduction in atmospheric carbon dioxide (CO2). In a sediment core from the Subantarctic Atlantic, we measured foraminifera-bound nitrogen isotopes to reconstruct ice age nitrate consumption, burial fluxes of iron, and proxies for productivity. Peak glacial times and millennial cold events are characterized by increases in dust flux, productivity, and the degree of nitrate consumption; this combination is uniquely consistent with Subantarctic iron fertilization. The associated strengthening of the Southern Ocean's biological pump can explain the lowering of CO2 at the transition from mid-climate states to full ice age conditions as well as the millennial-scale CO2 oscillations.

Rapid reductions in North Atlantic Deep Water during the peak of the last interglacial period

Deep ocean circulation has been considered relatively stable during interglacial periods, yet little is known about its behavior on submillennial time scales. Using a subcentennially resolved epibenthic foraminiferal δ(13)C record, we show that the influence of North Atlantic Deep Water (NADW) was strong at the onset of the last interglacial period and was then interrupted by several prominent centennial-scale reductions. These NADW transients occurred during periods of increased ice rafting and southward expansions of polar water influence, suggesting that a buoyancy threshold for convective instability was triggered by freshwater and circum-Arctic cryosphere changes. The deep Atlantic chemical changes were similar in magnitude to those associated with glaciations, implying that the canonical view of a relatively stable interglacial circulation may not hold for conditions warmer and fresher than at present.

Correlating the ancient Maya and modern European calendars with high-precision AMS 14C dating

The reasons for the development and collapse of Maya civilization remain controversial and historical events carved on stone monuments throughout this region provide a remarkable source of data about the rise and fall of these complex polities. Use of these records depends on correlating the Maya and European calendars so that they can be compared with climate and environmental datasets. Correlation constants can vary up to 1000 years and remain controversial. We report a series of high-resolution AMS (14)C dates on a wooden lintel collected from the Classic Period city of Tikal bearing Maya calendar dates. The radiocarbon dates were calibrated using a Bayesian statistical model and indicate that the dates were carved on the lintel between AD 658-696. This strongly supports the Goodman-Martínez-Thompson (GMT) correlation and the hypothesis that climate change played an important role in the development and demise of this complex civilization.

Two modes of change in Southern Ocean productivity over the past million years

Export of organic carbon from surface waters of the Antarctic Zone of the Southern Ocean decreased during the last ice age, coinciding with declining atmospheric carbon dioxide (CO(2)) concentrations, signaling reduced exchange of CO(2) between the ocean interior and the atmosphere. In contrast, in the Subantarctic Zone, export production increased into ice ages coinciding with rising dust fluxes, thus suggesting iron fertilization of subantarctic phytoplankton. Here, a new high-resolution productivity record from the Antarctic Zone is compiled with parallel subantarctic data over the past million years. Together, they fit the view that the combination of these two modes of Southern Ocean change determines the temporal structure of the glacial-interglacial atmospheric CO(2) record, including during the interval of "lukewarm" interglacials between 450 and 800 thousand years ago.

Rapid Acidification of the Ocean During the Paleocene-Eocene Thermal Maximum

The Paleocene-Eocene thermal maximum (PETM) has been attributed to the rapid release of approximately 2000 x 10(9) metric tons of carbon in the form of methane. In theory, oxidation and ocean absorption of this carbon should have lowered deep-sea pH, thereby triggering a rapid (<10,000-year) shoaling of the calcite compensation depth (CCD), followed by gradual recovery. Here we present geochemical data from five new South Atlantic deep-sea sections that constrain the timing and extent of massive sea-floor carbonate dissolution coincident with the PETM. The sections, from between 2.7 and 4.8 kilometers water depth, are marked by a prominent clay layer, the character of which indicates that the CCD shoaled rapidly (<10,000 years) by more than 2 kilometers and recovered gradually (>100,000 years). These findings indicate that a large mass of carbon (>>2000 x 10(9) metric tons of carbon) dissolved in the ocean at the Paleocene-Eocene boundary and that permanent sequestration of this carbon occurred through silicate weathering feedback.

Abrupt Temperature Changes in the Western Mediterranean over the Past 250,000 Years

A continuous high-resolution Western Mediterranean sea surface temperature (SST) alkenone record spanning the past 250,000 years shows that abrupt changes were more common at warming than at cooling. During marine isotope stage (MIS) 6, SST oscillated following a stadial-interstadial pattern but at lower intensities and rates of change than in the Dansgaard/Oeschger events of MIS 3. Some of the most prominent events occurred over MISs 5 and 7, after prolonged warm periods of high stability. Climate during the whole period was predominantly maintained in interglacial-interstadial conditions, whereas the duration of stadials was much shorter.

Solar forcing of drought frequency in the Maya lowlands

We analyzed lake-sediment cores from the Yucatan Peninsula, Mexico, to reconstruct the climate history of the region over the past 2600 years. Time series analysis of sediment proxies, which are sensitive to the changing ratio of evaporation to precipitation (oxygen isotopes and gypsum precipitation), reveal a recurrent pattern of drought with a dominant periodicity of 208 years. This cycle is similar to the documented 206-year period in records of cosmogenic nuclide production (carbon-14 and beryllium-10) that is thought to reflect variations in solar activity. We conclude that a significant component of century-scale variability in Yucatan droughts is explained by solar forcing. Furthermore, some of the maxima in the 208-year drought cycle correspond with discontinuities in Maya cultural evolution, suggesting that the Maya were affected by these bicentennial oscillations in precipitation.

Millennial-scale instability of the antarctic ice sheet during the last glaciation

Records of ice-rafted detritus (IRD) concentration in deep-sea cores from the southeast Atlantic Ocean reveal millennial-scale pulses of IRD delivery between 20,000 and 74,000 years ago. Prominent IRD layers correlate across the Polar Frontal Zone, suggesting episodes of Antarctic Ice Sheet instability. Carbon isotopes (delta(13)C) of benthic foraminifers, a proxy of deepwater circulation, reveal that South Atlantic IRD events coincided with strong increases in North Atlantic Deep Water (NADW) production and inferred warming (interstadials) in the high-latitude North Atlantic. Sea level rise or increased NADW production associated with strong interstadials may have resulted in destabilization of grounded ice shelves and possible surging in the Weddell Sea region of West Antarctica.