Timing of abrupt climate change at the end of the Younger Dryas interval from thermally fractionated gases in polar ice

Into the Holocene, anatomy of the Younger Dryas cold reversal and preboreal oscillation

During the most recent deglaciation, the upwards trend of warmer Northern Hemisphere (NH) temperatures was punctuated by a rapid and intense return to glacial conditions: the Younger Dryas (YD). The end of this event marks the beginning of the Holocene. Using the University of Toronto version of CCSM4, a model of the climate prior to the YD was created with correct boundary conditions. Various amounts of freshwater forcing were then applied to the Beaufort Gyre for forcing intervals ranging from 1 to 125 years. In several cases, this was sufficient to collapse the Atlantic Meridional Overturning Circulation (AMOC) and cause significant cooling over the NH. Crucially, after the forcing was ceased, the AMOC stayed in an off state for approximately a millennium before mounting a rapid recover to pre-YD levels. This recovery, which permanently reduced the extent of NH sea ice, occurred through the mechanism of a Polynya opening in the Irminger Sea during winter and led to a pronounced "overshoot" of the AMOC, during which NH temperatures were higher than before the YD.

Bipolar impact and phasing of Heinrich-type climate variability

During the last ice age, the Laurentide Ice Sheet exhibited extreme iceberg discharge events that are recorded in North Atlantic sediments1. These Heinrich events have far-reaching climate impacts, including widespread disruptions to hydrological and biogeochemical cycles2-4. They occurred during Heinrich stadials-cold periods with strongly weakened Atlantic overturning circulation5-7. Heinrich-type variability is not distinctive in Greenland water isotope ratios, a well-dated site temperature proxy8, complicating efforts to assess their regional climate impact and phasing against Antarctic climate change. Here we show that Heinrich events have no detectable temperature impact on Greenland and cooling occurs at the onset of several Heinrich stadials, and that both types of Heinrich variability have a distinct imprint on Antarctic climate. Antarctic ice cores show accelerated warming that is synchronous with increases in methane during Heinrich events, suggesting an atmospheric teleconnection9, despite the absence of a Greenland climate signal. Greenland ice-core nitrogen stable isotope ratios, a sensitive temperature proxy, indicate an abrupt cooling of about three degrees Celsius at the onset of Heinrich Stadial 1 (17.8 thousand years before present, where present is defined as 1950). Antarctic warming lags this cooling by 133 ± 93 years, consistent with an oceanic teleconnection. Paradoxically, proximal sites are less affected by Heinrich events than remote sites, suggesting spatially complex event dynamics.

Interpreting spatially explicit variation in dietary proxies through species distribution modeling reveals foraging preferences of mammoth (Mammuthus) and American mastodon (Mammut americanum)

Introduction

The end Pleistocene was a time of considerable ecological upheaval. Recent work has explored the megafauna extinction’s role in altering ecosystem processes. Analyses of functional traits withing communities reveal hidden consequences of the megafauna extinction beyond declines in taxonomic diversity. Functional diversity analyses offer new insight into our understanding of past ecosystems and may even inform future rewilding efforts. However, the utility of functional diversity may be hampered by the use of discrete, taxon-level functional traits, such as dietary categories, that mask variation in functional diversity over space and time.

Methods

We present an approach in which species distribution modeling, in Maxent, provides context for interpreting variation in two widely used proxies for diet among fossil taxa: stable isotope analysis and dental microwear texture analysis. We apply this approach to two ecologically distinct taxa, the American mastodon (Mammut americanum) and mammoths (Mammuthus) and investigate their resource use over space and time from the last glacial maximum to the end Pleistocene (25–11.7 thousand years before present).

Results

Mammoth dietary behavior varies by context across their geographic distribution, despite possessing evolutionary adaptations that facilitate grazing. Mammoths exhibit a preference for grazing where species distribution modeling predicts the highest likelihood of occurrence but engage in more mixed-feeding outside of core likelihood areas. In contrast, dietary preferences for mastodon are less resolved and our analyses were unable to identify significant differences in diet across their distribution.

Discussion

The ecological roles of some species are context specific and need to be critically evaluated when planning for management of reintroductions or introducing novel species to restore lost ecological function.

A Tibetan ice core covering the past 1,300 years radiometrically dated with 39Ar

Ice cores from alpine glaciers are unique archives of past global and regional climate conditions. However, recovering climate records from these ice cores is often hindered by the lack of a reliable chronology, especially in the age range of 100 to 500 anni (a) for which radiometric dating has not been available so far. We report on radiometric ³⁹Ar dating of an ice core from the Tibetan Plateau and the construction of a chronology covering the past 1,300 a using the obtained ³⁹Ar ages. This is made possible by advances in the analysis of ³⁹Ar using the laser-based detection method atom trap trace analysis, resulting in a twofold increase in the upper age limit of ³⁹Ar dating. By measuring the anthropogenic ⁸⁵Kr along with ³⁹Ar we quantify and correct modern air contamination, thus removing a major systematic uncertainty of ³⁹Ar dating. Moreover, the ⁸⁵Kr data for the top part of the ice core provide information on firn processes, including the age difference between the ice and its enclosed gas. This first application of ³⁹Ar and ⁸⁵Kr to an ice core facilitates further ice cores from nonpolar glaciers to be used for recovering climate records of the Common Era, a period including pronounced anomalies such as the Little Ice Age and the Medieval Warm Period.

Multiple carbon cycle mechanisms associated with the glaciation of Marine Isotope Stage 4

Here we use high-precision carbon isotope data (δ¹³C-CO₂) to show atmospheric CO₂ during Marine Isotope Stage 4 (MIS 4, ~70.5-59 ka) was controlled by a succession of millennial-scale processes. Enriched δ¹³C-CO₂ during peak glaciation suggests increased ocean carbon storage. Variations in δ¹³C-CO₂ in early MIS 4 suggest multiple processes were active during CO₂ drawdown, potentially including decreased land carbon and decreased Southern Ocean air-sea gas exchange superposed on increased ocean carbon storage. CO₂ remained low during MIS 4 while δ¹³C-CO₂ fluctuations suggest changes in Southern Ocean and North Atlantic air-sea gas exchange. A 7 ppm increase in CO₂ at the onset of Dansgaard-Oeschger event 19 (72.1 ka) and 27 ppm increase in CO₂ during late MIS 4 (Heinrich Stadial 6, ~63.5-60 ka) involved additions of isotopically light carbon to the atmosphere. The terrestrial biosphere and Southern Ocean air-sea gas exchange are possible sources, with the latter event also involving decreased ocean carbon storage.

Summary for general audience: We used carbon stable isotope data from an Antarctic ice core to evaluate which mechanisms caused changes in atmospheric carbon dioxide 74-59 thousand years ago, including a ~40 ppm decrease at the beginning of the last ice age.

Genomic insights into the recent chromosome reduction of autopolyploid sugarcane Saccharum spontaneum

Saccharum spontaneum is a founding Saccharum species and exhibits wide variation in ploidy levels. We have assembled a high-quality autopolyploid genome of S. spontaneum Np-X (2n = 4x = 40) into 40 pseudochromosomes across 10 homologous groups, that better elucidates recent chromosome reduction and polyploidization that occurred circa 1.5 million years ago (Mya). One paleo-duplicated chromosomal pair in Saccharum, NpChr5 and NpChr8, underwent fission followed by fusion accompanied by centromeric split around 0.80 Mya. We inferred that Np-X, with x = 10, most likely represents the ancestral karyotype, from which x = 9 and x = 8 evolved. Resequencing of 102 S. spontaneum accessions revealed that S. spontaneum originated in northern India from an x = 10 ancestor, which then radiated into four major groups across the Indian subcontinent, China, and Southeast Asia. Our study suggests new directions for accelerating sugarcane improvement and expands our knowledge of the evolution of autopolyploids.

Automatic detection of abrupt transitions in paleoclimate records

Bifurcations and tipping points (TPs) are an important part of the Earth system's behavior. These critical points represent thresholds at which small changes in the system's parameters or in the forcing abruptly switch it from one state or type of behavior to another. Current concern with TPs is largely due to the potential of slow anthropogenic forcing to bring about abrupt, and possibly irreversible, change to the physical climate system and impacted ecosystems. Paleoclimate proxy records have been shown to contain abrupt transitions, or "jumps," which may represent former instances of such dramatic climate change events. These transitions can provide valuable information for identifying critical TPs in current and future climate evolution. Here, we present a robust methodology for detecting abrupt transitions in proxy records that is applied to ice core and speleothem records of the last climate cycle. This methodology is based on the nonparametric Kolmogorov-Smirnov (KS) test for the equality, or not, of the probability distributions associated with two samples drawn from a time series, before and after any potential jump. To improve the detection of abrupt transitions in proxy records, the KS test is augmented by several other criteria and it is compared with recurrence analysis. The augmented KS test results show substantial skill when compared with more subjective criteria for jump detection. This test can also usefully complement recurrence analysis and improve upon certain aspects of its results.

Antarctic surface temperature and elevation during the Last Glacial Maximum

Water-stable isotopes in polar ice cores are a widely used temperature proxy in paleoclimate reconstruction, yet calibration remains challenging in East Antarctica. Here, we reconstruct the magnitude and spatial pattern of Last Glacial Maximum surface cooling in Antarctica using borehole thermometry and firn properties in seven ice cores. West Antarctic sites cooled ~10°C relative to the preindustrial period. East Antarctic sites show a range from ~4° to ~7°C cooling, which is consistent with the results of global climate models when the effects of topographic changes indicated with ice core air-content data are included, but less than those indicated with the use of water-stable isotopes calibrated against modern spatial gradients. An altered Antarctic temperature inversion during the glacial reconciles our estimates with water-isotope observations.

Abrupt Heinrich Stadial 1 cooling missing in Greenland oxygen isotopes

Abrupt climate changes during the last deglaciation have been well preserved in proxy records across the globe. However, one long-standing puzzle is the apparent absence of the onset of the Heinrich Stadial 1 (HS1) cold event around 18 ka in Greenland ice core oxygen isotope δ¹⁸ O records, inconsistent with other proxies. Here, combining proxy records with an isotope-enabled transient deglacial simulation, we propose that a substantial HS1 cooling onset did indeed occur over the Arctic in winter. However, this cooling signal in the depleted oxygen isotopic composition is completely compensated by the enrichment because of the loss of winter precipitation in response to sea ice expansion associated with AMOC slowdown during extreme glacial climate. In contrast, the Arctic summer warmed during HS1 and YD because of increased insolation and greenhouse gases, consistent with snowline reconstructions. Our work suggests that Greenland δ¹⁸ O may substantially underestimate temperature variability during cold glacial conditions.

Volcanic origin for Younger Dryas geochemical anomalies ca. 12,900 cal B.P

The Younger Dryas (YD) abrupt cooling event ca. 12.9 ± 0.1 ka is associated with substantial meltwater input into the North Atlantic Ocean, reversing deglacial warming. One controversial and prevailing hypothesis is that a bolide impact or airburst is responsible for these environmental changes. Here, highly siderophile element (HSE; Os, Ir, Ru, Pt, Pd, and Re) abundances and 187Os/188Os ratios were obtained in a well-dated sediment section at Hall's Cave, TX, USA to test this hypothesis. In Hall's Cave, layers below, above, and in the YD have 187Os/188Os ratios consistent with incorporation of extraterrestrial or mantle-derived material. The HSE abundances indicate that these layers contain volcanic gas aerosols and not extraterrestrial materials. The most likely explanation is that episodic, distant volcanic emissions were deposited in Hall's Cave sediments. Coupled 187Os/188Os ratios and HSE concentration data at close stratigraphic intervals are required to effectively differentiate between bolide and volcanic origins.

Differing climatic mechanisms control transient and accumulated vegetation novelty in Europe and eastern North America

Understanding the mechanisms of climate that produce novel ecosystems is of joint interest to conservation biologists and palaeoecologists. Here, we define and differentiate transient from accumulated novelty and evaluate four climatic mechanisms proposed to cause species to reshuffle into novel assemblages: high climatic novelty, high spatial rates of change (displacement), high variance among displacement rates for individual climate variables, and divergence among displacement vector bearings. We use climate simulations to quantify climate novelty, displacement and divergence across Europe and eastern North America from the last glacial maximum to the present, and fossil pollen records to quantify vegetation novelty. Transient climate novelty is consistently the strongest predictor of transient vegetation novelty, while displacement rates (mean and variance) are equally important in Europe. However, transient vegetation novelty is lower in Europe and its relationship to climatic predictors is the opposite of expectation. For both continents, accumulated novelty is greater than transient novelty, and climate novelty is the strongest predictor of accumulated ecological novelty. These results suggest that controls on novel ecosystems vary with timescale and among continents, and that the twenty-first century emergence of novelty will be driven by both rapid rates of climate change and the emergence of novel climate states. This article is part of a discussion meeting issue 'The past is a foreign country: how much can the fossil record actually inform conservation?'

Earth's radiative imbalance from the Last Glacial Maximum to the present

The energy imbalance at the top of the atmosphere determines the temporal evolution of the global climate, and vice versa changes in the climate system can alter the planetary energy fluxes. This interplay is fundamental to our understanding of Earth's heat budget and the climate system. However, even today, the direct measurement of global radiative fluxes is difficult, such that most assessments are based on changes in the total energy content of the climate system. We apply the same approach to estimate the long-term evolution of Earth's radiative imbalance in the past. New measurements of noble gas-derived mean ocean temperature from the European Project for Ice Coring in Antarctica Dome C ice core covering the last 40,000 y, combined with recent results from the West Antarctic Ice Sheet Divide ice core and the sea-level record, allow us to quantitatively reconstruct the history of the climate system energy budget. The temporal derivative of this quantity must be equal to the planetary radiative imbalance. During the deglaciation, a positive imbalance of typically +0.2 W⋅m^-2 is maintained for ∼10,000 y, however, with two distinct peaks that reach up to 0.4 W⋅m^-2 during times of substantially reduced Atlantic Meridional Overturning Circulation. We conclude that these peaks are related to net changes in ocean heat uptake, likely due to rapid changes in North Atlantic deep-water formation and their impact on the global radiative balance, while changes in cloud coverage, albeit uncertain, may also factor into the picture.

Millennial-Scale Climate Variability and Dinoflagellate-Cyst-Based Seasonality Changes Over the Last ~150 kyrs at "Shackleton Site" U1385

During the last glacial period, climate conditions in the North Atlantic region were determined by the alternation of relatively warm interstadials and relatively cool stadials, with superimposed rapid warming (Dansgaard-Oeschger) and cooling (Heinrich) events. So far little is known about the impact of these rapid climate shifts on the seasonal variations in sea surface temperature (SST) within the North Atlantic region. Here, we present a high-resolution seasonal SST record for the past 152 kyrs derived from Integrated Ocean Drilling Program "Shackleton" Site U1385, offshore Portugal. Assemblage counts of dinoflagellates cysts (dinocysts) in combination with a modern analog technique (MAT), and regression analyses were used for the reconstructions. We compare our records with previously published SST records from the same location obtained from the application of MAT on planktonic foraminifera. Our dinocyst-based reconstructions confirm the impression of the Greenland stadials and interstadials offshore the Portuguese margin and indicate increased seasonal contrast of temperature during the cold periods of the glacial cycle (average 9.0 °C, maximum 12.2 °C) with respect to present day (5.1 °C), due to strong winter cooling by up to 8.3 °C. Our seasonal temperature reconstructions are in line with previously published data, which showed increased seasonality due to strong winter cooling during the Younger Dryas and the Last Glacial Maximum over the European continent and North Atlantic region. In addition, we show that over longer time scales, increased seasonal contrasts of temperature remained characteristic of the colder phases of the glacial cycle.

Mean global ocean temperatures during the last glacial transition

Little is known about the ocean temperature's long-term response to climate perturbations owing to limited observations and a lack of robust reconstructions. Although most of the anthropogenic heat added to the climate system has been taken up by the ocean up until now, its role in a century and beyond is uncertain. Here, using noble gases trapped in ice cores, we show that the mean global ocean temperature increased by 2.57 ± 0.24 degrees Celsius over the last glacial transition (20,000 to 10,000 years ago). Our reconstruction provides unprecedented precision and temporal resolution for the integrated global ocean, in contrast to the depth-, region-, organism- and season-specific estimates provided by other methods. We find that the mean global ocean temperature is closely correlated with Antarctic temperature and has no lead or lag with atmospheric CO₂, thereby confirming the important role of Southern Hemisphere climate in global climate trends. We also reveal an enigmatic 700-year warming during the early Younger Dryas period (about 12,000 years ago) that surpasses estimates of modern ocean heat uptake.

Volcanic influence on centennial to millennial Holocene Greenland temperature change

Solar variability has been hypothesized to be a major driver of North Atlantic millennial-scale climate variations through the Holocene along with orbitally induced insolation change. However, another important climate driver, volcanic forcing has generally been underestimated prior to the past 2,500 years partly owing to the lack of proper proxy temperature records. Here, we reconstruct seasonally unbiased and physically constrained Greenland Summit temperatures over the Holocene using argon and nitrogen isotopes within trapped air in a Greenland ice core (GISP2). We show that a series of volcanic eruptions through the Holocene played an important role in driving centennial to millennial-scale temperature changes in Greenland. The reconstructed Greenland temperature exhibits significant millennial correlations with K⁺ and Na⁺ ions in the GISP2 ice core (proxies for atmospheric circulation patterns), and δ¹⁸O of Oman and Chinese Dongge cave stalagmites (proxies for monsoon activity), indicating that the reconstructed temperature contains hemispheric signals. Climate model simulations forced with the volcanic forcing further suggest that a series of large volcanic eruptions induced hemispheric-wide centennial to millennial-scale variability through ocean/sea-ice feedbacks. Therefore, we conclude that volcanic activity played a critical role in driving centennial to millennial-scale Holocene temperature variability in Greenland and likely beyond.

A Pleistocene ice core record of atmospheric O2 concentrations

The history of atmospheric O₂ partial pressures (Po₂) is inextricably linked to the coevolution of life and Earth's biogeochemical cycles. Reconstructions of past Po₂ rely on models and proxies but often markedly disagree. We present a record of Po₂ reconstructed using O₂/N₂ ratios from ancient air trapped in ice. This record indicates that Po₂ declined by 7 per mil (0.7%) over the past 800,000 years, requiring that O₂ sinks were ~2% larger than sources. This decline is consistent with changes in burial and weathering fluxes of organic carbon and pyrite driven by either Neogene cooling or increasing Pleistocene erosion rates. The 800,000-year record of steady average carbon dioxide partial pressures (Pco₂) but declining Po₂ provides distinctive evidence that a silicate weathering feedback stabilizes Pco₂ on million-year time scales.

Greenland temperature response to climate forcing during the last deglaciation

Greenland ice core water isotopic composition (δ(18)O) provides detailed evidence for abrupt climate changes but is by itself insufficient for quantitative reconstruction of past temperatures and their spatial patterns. We investigate Greenland temperature evolution during the last deglaciation using independent reconstructions from three ice cores and simulations with a coupled ocean-atmosphere climate model. Contrary to the traditional δ(18)O interpretation, the Younger Dryas period was 4.5° ± 2°C warmer than the Oldest Dryas, due to increased carbon dioxide forcing and summer insolation. The magnitude of abrupt temperature changes is larger in central Greenland (9° to 14°C) than in the northwest (5° to 9°C), fingerprinting a North Atlantic origin. Simulated changes in temperature seasonality closely track changes in the Atlantic overturning strength and support the hypothesis that abrupt climate change is mostly a winter phenomenon.

Younger Dryas deglaciation of Scotland driven by warming summers

The Younger Dryas Stadial (YDS; ∼ 12,900-11,600 y ago) in the Northern Hemisphere is classically defined by abrupt cooling and renewed glaciation during the last glacial-interglacial transition. Although this event involved a global reorganization of atmospheric and oceanic circulation [Denton GH, Alley RB, Comer GC, Broecker WS (2005) Quat Sci Rev 24:1159-1182], the magnitude, seasonality, and geographical footprint of YDS cooling remain unresolved and pose a challenge to our understanding of abrupt climate change. Here, we present a deglacial chronology from Scotland, immediately downwind of the North Atlantic Ocean, indicating that the Scottish ice cap disintegrated during the first half of the YDS. We suggest that stratification of the North Atlantic Ocean resulted in amplified seasonality that, paradoxically, stimulated a severe wintertime climate while promoting warming summers through solar heating of the mixed layer. This latter process drove deglaciation of downwind landmasses to completion well before the end of the YDS.

Role of megafauna and frozen soil in the atmospheric CH4 dynamics

Modern wetlands are the world’s strongest methane source. But what was the role of this source in the past? An analysis of global ¹⁴C data for basal peat combined with modelling of wetland succession allowed us to reconstruct the dynamics of global wetland methane emission through time. These data show that the rise of atmospheric methane concentrations during the Pleistocene-Holocene transition was not connected with wetland expansion, but rather started substantially later, only 9 thousand years ago. Additionally, wetland expansion took place against the background of a decline in atmospheric methane concentration. The isotopic composition of methane varies according to source. Owing to ice sheet drilling programs past dynamics of atmospheric methane isotopic composition is now known. For example over the course of Pleistocene-Holocene transition atmospheric methane became depleted in the deuterium isotope, which indicated that the rise in methane concentrations was not connected with activation of the deuterium-rich gas clathrates. Modelling of the budget of the atmospheric methane and its isotopic composition allowed us to reconstruct the dynamics of all main methane sources. For the late Pleistocene, the largest methane source was megaherbivores, whose total biomass is estimated to have exceeded that of present-day humans and domestic animals. This corresponds with our independent estimates of herbivore density on the pastures of the late Pleistocene based on herbivore skeleton density in the permafrost. During deglaciation, the largest methane emissions originated from degrading frozen soils of the mammoth steppe biome. Methane from this source is unique, as it is depleted of all isotopes. We estimated that over the entire course of deglaciation (15,000 to 6,000 year before present), soils of the mammoth steppe released 300–550 Pg (10¹⁵ g) of methane. From current study we conclude that the Late Quaternary Extinction significantly affected the global methane cycle.

Integrating multiple lines of evidence into historical biogeography hypothesis testing: a Bison bison case study

One of the grand goals of historical biogeography is to understand how and why species' population sizes and distributions change over time. Multiple types of data drawn from disparate fields, combined into a single modelling framework, are necessary to document changes in a species's demography and distribution, and to determine the drivers responsible for change. Yet truly integrated approaches are challenging and rarely performed. Here, we discuss a modelling framework that integrates spatio-temporal fossil data, ancient DNA, palaeoclimatological reconstructions, bioclimatic envelope modelling and coalescence models in order to statistically test alternative hypotheses of demographic and potential distributional changes for the iconic American bison (Bison bison). Using different assumptions about the evolution of the bioclimatic niche, we generate hypothetical distributional and demographic histories of the species. We then test these demographic models by comparing the genetic signature predicted by serial coalescence against sequence data derived from subfossils and modern populations. Our results supported demographic models that include both climate and human-associated drivers of population declines. This synthetic approach, integrating palaeoclimatology, bioclimatic envelopes, serial coalescence, spatio-temporal fossil data and heterochronous DNA sequences, improves understanding of species' historical biogeography by allowing consideration of both abiotic and biotic interactions at the population level.

Diet of upper paleolithic modern humans: evidence from microwear texture analysis

This article presents the results of the occlusal molar microwear texture analysis of 32 adult Upper Paleolithic modern humans from a total of 21 European sites dating to marine isotope stages 3 and 2. The occlusal molar microwear textures of these specimens were analyzed with the aim of examining the effects of the climatic, as well as the cultural, changes on the diets of the Upper Paleolithic modern humans. The results of this analysis do not reveal any environmentally driven dietary shifts for the Upper Paleolithic hominins indicating that the climatic and their associated paleoecological changes did not force these humans to significantly alter their diets in order to survive. However, the microwear texture analysis does detect culturally related changes in the Upper Paleolithic humans' diets. Specifically, significant differences in diet were found between the earlier Upper Paleolithic individuals, i.e., those belonging to the Aurignacian and Gravettian contexts, and the later Magdalenian ones, such that the diet of the latter group was more varied and included more abrasive foods compared with those of the former.

Onset of deglacial warming in West Antarctica driven by local orbital forcing

The cause of warming in the Southern Hemisphere during the most recent deglaciation remains a matter of debate. Hypotheses for a Northern Hemisphere trigger, through oceanic redistributions of heat, are based in part on the abrupt onset of warming seen in East Antarctic ice cores and dated to 18,000 years ago, which is several thousand years after high-latitude Northern Hemisphere summer insolation intensity began increasing from its minimum, approximately 24,000 years ago. An alternative explanation is that local solar insolation changes cause the Southern Hemisphere to warm independently. Here we present results from a new, annually resolved ice-core record from West Antarctica that reconciles these two views. The records show that 18,000 years ago snow accumulation in West Antarctica began increasing, coincident with increasing carbon dioxide concentrations, warming in East Antarctica and cooling in the Northern Hemisphere associated with an abrupt decrease in Atlantic meridional overturning circulation. However, significant warming in West Antarctica began at least 2,000 years earlier. Circum-Antarctic sea-ice decline, driven by increasing local insolation, is the likely cause of this warming. The marine-influenced West Antarctic records suggest a more active role for the Southern Ocean in the onset of deglaciation than is inferred from ice cores in the East Antarctic interior, which are largely isolated from sea-ice changes.

Younger Dryas cooling and the Greenland climate response to CO2

Greenland ice-core δ(18)O-temperature reconstructions suggest a dramatic cooling during the Younger Dryas (YD; 12.9-11.7 ka), with temperatures being as cold as the earlier Oldest Dryas (OD; 18.0-14.6 ka) despite an approximately 50 ppm rise in atmospheric CO(2). Such YD cooling implies a muted Greenland climate response to atmospheric CO(2), contrary to physical predictions of an enhanced high-latitude response to future increases in CO(2). Here we show that North Atlantic sea surface temperature reconstructions as well as transient climate model simulations suggest that the YD over Greenland should be substantially warmer than the OD by approximately 5 °C in response to increased atmospheric CO(2). Additional experiments with an isotope-enabled model suggest that the apparent YD temperature reconstruction derived from the ice-core δ(18)O record is likely an artifact of an altered temperature-δ(18)O relationship due to changing deglacial atmospheric circulation. Our results thus suggest that Greenland climate was warmer during the YD relative to the OD in response to rising atmospheric CO(2), consistent with sea surface temperature reconstructions and physical predictions, and has a sensitivity approximately twice that found in climate models for current climate due to an enhanced albedo feedback during the last deglaciation.

Application of the N(2)/Ar technique to measuring soil-atmosphere N(2) fluxes

RATIONALE

The emission of dinitrogen (N(2) ) gas from soil is the most poorly constrained flux in terrestrial nitrogen (N) budgets because the high background atmospheric N(2) concentration makes soil N(2) emissions difficult to measure. In this study, we tested the theoretical and analytical feasibility of using the N(2) /Ar technique to measure soil-atmosphere N(2) fluxes.

METHODS

Dual inlet isotope ratio mass spectrometry was used to measure δAr/N(2) values of gas sampled from surface flux chambers. In laboratory experiments using dry sand in a diffusion box, we induced a known steady-state flux of N(2) , and then measured the change in the N(2) /Ar ratio of chamber headspace air samples to test our ability to reconstruct this flux. We m\odeled solubility, thermal, and water vapor flux fractionation effects on the N(2) /Ar ratio to constrain physical effects on the measured N(2) flux.

RESULTS

In dry sand, an actual N(2) flux of 108 mg N m(-2)  day(-1) was measured as 111 ± 19 mg N m(-2)  day(-1) (± standard error (SE)). In wet sand, an actual N(2) flux of 160 mg N m(-2)  day(-1) was measured as 146 ± 20 mg N m(-2)  day(-1) when solubility and water vapor flux fractionation were taken into account. Corrections for thermal fractionation did not improve estimates of N(2) fluxes.

CONCLUSIONS

We conclude that our application of the N(2) /Ar technique to soil surface fluxes is valid only above a detection limit of approximately 108 mg N m(-2)  day(-1) . The N(2) /Ar method is currently best used as a validation tool for other methods in ecosystems with high soil N(2) fluxes, but, with future improvements, it holds promise to provide high-resolution measurements in systems with low soil N(2) fluxes.

Faunal histories from Holocene ancient DNA

Recent studies using ancient DNA have been instrumental in advancing understanding of the impact of Holocene climate change on biodiversity. Ancient DNA has been used to track demography, migration and diversity, and is providing new insights into the long-term dynamics of species and population distributions. The Holocene is key to understanding how the past has impacted on the present, as it bridges the gap between contemporary phylogeographic studies and those with inference on Pleistocene patterns, based on ancient DNA studies. Here, we examine the major patterns of Holocene faunal population dynamics and connectivity; highlighting the dynamic nature of species and population responses to Holocene climatic change, thereby providing an 'analogue' for understanding potential impacts of future change.

Greenhouse gases in the Earth system: a palaeoclimate perspective

While the trends in greenhouse gas concentrations in recent decades are clear, their significance is only revealed when viewed in the context of a longer time period. Fortunately, the air bubbles in polar ice cores provide an unusually direct method of determining the concentrations of stable gases over a period of (so far) 800,000 years. Measurements on different cores with varying characteristics, as well as an overlap of ice-core and atmospheric measurements covering the same time period, show that the ice-core record provides a faithful record of changing atmospheric composition. The mixing ratio of CO(2) is now 30 per cent higher than any value observed in the ice-core record, while methane is more than double any observed value; the rate of change also appears extraordinary compared with natural changes. Before the period when anthropogenic changes have dominated, there are very interesting natural changes in concentration, particularly across glacial/interglacial cycles, and these can be used to understand feedbacks in the Earth system. The phasing of changes in temperature and CO(2) across glacial/interglacial transitions is consistent with the idea that CO(2) acts as an important amplifier of climate changes in the natural system. Even larger changes are inferred to have occurred in periods earlier than the ice cores cover, and these events might be used to constrain assessments of the way the Earth could respond to higher than present concentrations of CO(2), and to a large release of carbon: however, more certainty about CO(2) concentrations beyond the time period covered by ice cores is needed before such constraints can be fully realized.

Climate change and evolutionary adaptations at species' range margins

During recent climate warming, many insect species have shifted their ranges to higher latitudes and altitudes. These expansions mirror those that occurred after the Last Glacial Maximum when species expanded from their ice age refugia. Postglacial range expansions have resulted in clines in genetic diversity across present-day distributions, with a reduction in genetic diversity observed in a wide range of insect taxa as one moves from the historical distribution core to the current range margin. Evolutionary increases in dispersal at expanding range boundaries are commonly observed in virtually all insects that have been studied, suggesting a positive feedback between range expansion and the evolution of traits that accelerate range expansion. The ubiquity of this phenomenon suggests that it is likely to be an important determinant of range changes. A better understanding of the extent and speed of adaptation will be crucial to the responses of biodiversity and ecosystems to climate change.

No evidence of nanodiamonds in Younger-Dryas sediments to support an impact event

The causes of the late Pleistocene megafaunal extinctions in North America, disappearance of Clovis paleoindian lithic technology, and abrupt Younger-Dryas (YD) climate reversal of the last deglacial warming in the Northern Hemisphere remain an enigma. A controversial hypothesis proposes that one or more cometary airbursts/impacts barraged North America ≈12,900 cal yr B.P. and caused these events. Most evidence supporting this hypothesis has been discredited except for reports of nanodiamonds (including the rare hexagonal polytype) in Bølling-Ållerod-YD-boundary sediments. The hexagonal polytype of diamond, lonsdaleite, is of particular interest because it is often associated with shock pressures related to impacts where it has been found to occur naturally. Unfortunately, previous reports of YD-boundary nanodiamonds have left many unanswered questions regarding the nature and occurrence of the nanodiamonds. Therefore, we examined carbon-rich materials isolated from sediments dated 15,818 cal yr B.P. to present (including the Bølling-Allerod-YD boundary). No nanodiamonds were found in our study. Instead, graphene- and graphene/graphane-oxide aggregates are ubiquitous in all specimens examined. We demonstrate that previous studies misidentified graphene/graphane-oxide aggregates as hexagonal diamond and likely misidentified graphene as cubic diamond. Our results cast doubt upon one of the last widely discussed pieces of evidence supporting the YD impact hypothesis.

The last glacial termination

A major puzzle of paleoclimatology is why, after a long interval of cooling climate, each late Quaternary ice age ended with a relatively short warming leg called a termination. We here offer a comprehensive hypothesis of how Earth emerged from the last global ice age. A prerequisite was the growth of very large Northern Hemisphere ice sheets, whose subsequent collapse created stadial conditions that disrupted global patterns of ocean and atmospheric circulation. The Southern Hemisphere westerlies shifted poleward during each northern stadial, producing pulses of ocean upwelling and warming that together accounted for much of the termination in the Southern Ocean and Antarctica. Rising atmospheric CO2 during southern upwelling pulses augmented warming during the last termination in both polar hemispheres.

Isotope fractionation in silicate melts by thermal diffusion

The phenomenon of thermal diffusion (mass diffusion driven by a temperature gradient, known as the Ludwig-Soret effect) has been investigated for over 150 years, but an understanding of its underlying physical basis remains elusive. A significant hurdle in studying thermal diffusion has been the difficulty of characterizing it. Extensive experiments over the past century have established that the Soret coefficient, S(T) (a single parameter that describes the steady-state result of thermal diffusion), is highly sensitive to many factors. This sensitivity makes it very difficult to obtain a robust characterization of thermal diffusion, even for a single material. Here we show that for thermal diffusion experiments that span a wide range in composition and temperature, the difference in S(T) between isotopes of diffusing elements that are network modifiers (iron, calcium and magnesium) is independent of the composition and temperature. On the basis of this finding, we propose an additive decomposition for the functional form of S(T) and argue that a theoretical approach based on local thermodynamic equilibrium holds promise for describing thermal diffusion in silicate melts and other complex solutions. Our results lead to a simple and robust framework for characterizing isotope fractionation by thermal diffusion in natural and synthetic systems.

Clathrate hydrates in nature

Scientific knowledge of natural clathrate hydrates has grown enormously over the past decade, with spectacular new findings of large exposures of complex hydrates on the sea floor, the development of new tools for examining the solid phase in situ, significant progress in modeling natural hydrate systems, and the discovery of exotic hydrates associated with sea floor venting of liquid CO2. Major unresolved questions remain about the role of hydrates in response to climate change today, and correlations between the hydrate reservoir of Earth and the stable isotopic evidence of massive hydrate dissociation in the geologic past. The examination of hydrates as a possible energy resource is proceeding apace for the subpermafrost accumulations in the Arctic, but serious questions remain about the viability of marine hydrates as an economic resource. New and energetic explorations by nations such as India and China are quickly uncovering large hydrate findings on their continental shelves.

The contemporary degassing rate of 40Ar from the solid Earth

Knowledge of the outgassing history of radiogenic (40)Ar, derived over geologic time from the radioactive decay of (40)K, contributes to our understanding of the geodynamic history of the planet and the origin of volatiles on Earth's surface. The (40)Ar inventory of the atmosphere equals total (40)Ar outgassing during Earth history. Here, we report the current rate of (40)Ar outgassing, accessed by measuring the Ar isotope composition of trapped gases in samples of the Vostok and Dome C deep ice cores dating back to almost 800 ka. The modern outgassing rate (1.1 +/- 0.1 x 10(8) mol/yr) is in the range of values expected by summing outgassing from the continental crust and the upper mantle, as estimated from simple calculations and models. The measured outgassing rate is also of interest because it allows dating of air trapped in ancient ice core samples of unknown age, although uncertainties are large (+/-180 kyr for a single sample or +/-11% of the calculated age, whichever is greater).