Initial Upper Palaeolithic Homo sapiens from Bacho Kiro Cave, Bulgaria.

The Middle to Upper Palaeolithic transition in Europe witnessed the replacement and partial absorption of local Neanderthal populations by Homo sapiens populations of African origin1. However, this process probably varied across regions and its details remain largely unknown. In particular, the duration of chronological overlap between the two groups is much debated, as are the implications of this overlap for the nature of the biological and cultural interactions between Neanderthals and H. sapiens. Here we report the discovery and direct dating of human remains found in association with Initial Upper Palaeolithic artefacts2, from excavations at Bacho Kiro Cave (Bulgaria). Morphological analysis of a tooth and mitochondrial DNA from several hominin bone fragments, identified through proteomic screening, assign these finds to H. sapiens and link the expansion of Initial Upper Palaeolithic technologies with the spread of H. sapiens into the mid-latitudes of Eurasia before 45 thousand years ago3. The excavations yielded a wealth of bone artefacts, including pendants manufactured from cave bear teeth that are reminiscent of those later produced by the last Neanderthals of western Europe4-6. These finds are consistent with models based on the arrival of multiple waves of H. sapiens into Europe coming into contact with declining Neanderthal populations7,8.

A new view on abrupt climate changes and the bipolar seesaw based on paleotemperatures from Iberian Margin sediments.

The last glacial cycle provides the opportunity to investigate large changes in the Atlantic Meridional Overturning Circulation (AMOC) beyond the small fluctuations evidenced from direct measurements. Paleotemperature records from Greenland and the North Atlantic show an abrupt variability, called Dansgaard-Oeschger (DO) events, which is associated with abrupt changes of the AMOC. These DO events also have Southern Hemisphere counterparts via the thermal bipolar seesaw, a concept describing the meridional heat transport leading to asynchronous temperature changes between both hemispheres. However, temperature records from the North Atlantic show more pronounced DO cooling events during massive releases of icebergs known as Heinrich (H) events, contrary to ice-core-based temperature records from Greenland. Here, we present high-resolution temperature records from the Iberian Margin and a Bipolar Seesaw Index to discriminate DO cooling events with and without H events. We show that the thermal bipolar seesaw model generates synthetic Southern Hemisphere temperature records that best resemble Antarctic temperature records when using temperature records from the Iberian Margin as inputs. Our data-model comparison emphasizes the role of the thermal bipolar seesaw in the abrupt temperature variability of both hemispheres with a clear enhancement during DO cooling events with H events, implying a relationship that is more complex than a simple flip-flop between two climate states linked to a tipping point threshold.

A radiocarbon spike at 14 300 cal yr BP in subfossil trees provides the impulse response function of the global carbon cycle during the Late Glacial.

We present new 14C results measured on subfossil Scots Pines recovered in the eroded banks of the Drouzet watercourse in the Southern French Alps. About 400 new 14C ages have been analysed on 15 trees sampled at annual resolution. The resulting Δ14C record exhibits an abrupt spike occurring in a single year at 14 300-14 299 cal yr BP and a century-long event between 14 and 13.9 cal kyr BP. In order to identify the causes of these events, we compare the Drouzet Δ14C record with simulations of Δ14C based on the 10Be record in Greenland ice used as an input of a carbon cycle model. The correspondence with 10Be anomalies allows us to propose the 14.3 cal kyr BP event as a solar energetic particle event. By contrast, the 14 cal kyr BP event lasted about a century and is most probably a common Maunder-type solar minimum linked to the modulation of galactic cosmic particles by the heliomagnetic field. We also discuss and speculate about the synchroneity and the possible causes of the 14 cal kyr BP event with the brief cold phase called Older Dryas, which separates the Bølling and Allerød millennium-long warm phases of the Late Glacial period. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Making the case for an International Decade of Radiocarbon.

Radiocarbon (14C) is a critical tool for understanding the global carbon cycle. During the Anthropocene, two new processes influenced 14C in atmospheric, land and ocean carbon reservoirs. First, 14C-free carbon derived from fossil fuel burning has diluted 14C, at rates that have accelerated with time. Second, 'bomb' 14C produced by atmospheric nuclear weapon tests in the mid-twentieth century provided a global isotope tracer that is used to constrain rates of air-sea gas exchange, carbon turnover, large-scale atmospheric and ocean transport, and other key C cycle processes. As we write, the 14C/12C ratio of atmospheric CO2 is dropping below pre-industrial levels, and the rate of decline in the future will depend on global fossil fuel use and net exchange of bomb 14C between the atmosphere, ocean and land. This milestone coincides with a rapid increase in 14C measurement capacity worldwide. Leveraging future 14C measurements to understand processes and test models requires coordinated international effort-a 'decade of radiocarbon' with multiple goals: (i) filling observational gaps using archives, (ii) building and sustaining observation networks to increase measurement density across carbon reservoirs, (iii) developing databases, synthesis and modelling tools and (iv) establishing metrics for identifying and verifying changes in carbon sources and sinks. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Extended dilation of the radiocarbon time scale between 40,000 and 48,000 y BP and the overlap between Neanderthals and Homo sapiens.

The new radiocarbon calibration curve (IntCal20) allows us to calculate the gradient of the relationship between 14C age and calendar age over the past 55 millennia before the present (55 ka BP). The new gradient curve exhibits a prolonged and prominent maximum between 48 and 40 ka BP during which the radiocarbon clock runs almost twice as fast as it should. This radiocarbon time dilation is due to the increase in the atmospheric 14C/12C ratio caused by the 14C production rise linked to the transition into the Laschamp geomagnetic excursion centered around 41 ka BP. The major maximum in the gradient from 48 to 40 ka BP is a new feature of the IntCal20 calibration curve, with far-reaching impacts for scientific communities, such as prehistory and paleoclimatology, relying on accurate ages in this time range. To illustrate, we consider the duration of the overlap between Neanderthals and Homo sapiens in Eurasia.

Multiradionuclide evidence for an extreme solar proton event around 2,610 B.P. (∼660 BC).

Recently, it has been confirmed that extreme solar proton events can lead to significantly increased atmospheric production rates of cosmogenic radionuclides. Evidence of such events is recorded in annually resolved natural archives, such as tree rings [carbon-14 (14C)] and ice cores [beryllium-10 (10Be), chlorine-36 (36Cl)]. Here, we show evidence for an extreme solar event around 2,610 years B.P. (∼660 BC) based on high-resolution 10Be data from two Greenland ice cores. Our conclusions are supported by modeled 14C production rates for the same period. Using existing 36Cl ice core data in conjunction with 10Be, we further show that this solar event was characterized by a very hard energy spectrum. These results indicate that the 2,610-years B.P. event was an order of magnitude stronger than any solar event recorded during the instrumental period and comparable with the solar proton event of AD 774/775, the largest solar event known to date. The results illustrate the importance of multiple ice core radionuclide measurements for the reliable identification of short-term production rate increases and the assessment of their origins.

Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation.

The covariation of carbon dioxide (CO(2)) concentration and temperature in Antarctic ice-core records suggests a close link between CO(2) and climate during the Pleistocene ice ages. The role and relative importance of CO(2) in producing these climate changes remains unclear, however, in part because the ice-core deuterium record reflects local rather than global temperature. Here we construct a record of global surface temperature from 80 proxy records and show that temperature is correlated with and generally lags CO(2) during the last (that is, the most recent) deglaciation. Differences between the respective temperature changes of the Northern Hemisphere and Southern Hemisphere parallel variations in the strength of the Atlantic meridional overturning circulation recorded in marine sediments. These observations, together with transient global climate model simulations, support the conclusion that an antiphased hemispheric temperature response to ocean circulation changes superimposed on globally in-phase warming driven by increasing CO(2) concentrations is an explanation for much of the temperature change at the end of the most recent ice age.

Ice-sheet collapse and sea-level rise at the Bølling warming 14,600 years ago.

Past sea-level records provide invaluable information about the response of ice sheets to climate forcing. Some such records suggest that the last deglaciation was punctuated by a dramatic period of sea-level rise, of about 20 metres, in less than 500 years. Controversy about the amplitude and timing of this meltwater pulse (MWP-1A) has, however, led to uncertainty about the source of the melt water and its temporal and causal relationships with the abrupt climate changes of the deglaciation. Here we show that MWP-1A started no earlier than 14,650 years ago and ended before 14,310 years ago, making it coeval with the Bølling warming. Our results, based on corals drilled offshore from Tahiti during Integrated Ocean Drilling Project Expedition 310, reveal that the increase in sea level at Tahiti was between 12 and 22 metres, with a most probable value between 14 and 18 metres, establishing a significant meltwater contribution from the Southern Hemisphere. This implies that the rate of eustatic sea-level rise exceeded 40 millimetres per year during MWP-1A.

Abrupt drainage cycles of the Fennoscandian Ice Sheet.

Continental ice sheets are a key component of the Earth's climate system, but their internal dynamics need to be further studied. Since the last deglaciation, the northern Eurasian Fennoscandian Ice Sheet (FIS) has been connected to the Black Sea (BS) watershed, making this basin a suitable location to investigate former ice-sheet dynamics. Here, from a core retrieved in the BS, we combine the use of neodymium isotopes, high-resolution elemental analysis, and biomarkers to trace changes in sediment provenance and river runoff. We reveal cyclic releases of meltwater originating from Lake Disna, a proglacial lake linked to the FIS during Heinrich Stadial 1. Regional interactions within the climate-lake-FIS system, linked to changes in the availability of subglacial water, led to abrupt drainage cycles of the FIS into the BS watershed. This phenomenon raised the BS water level by ∼100 m until the sill of the Bosphorus Strait was reached, flooding the vast northwestern BS shelf and deeply affecting the hydrology and circulation of the BS and, probably, of the Marmara and Aegean Seas.

Migration of the subtropical front as a modulator of glacial climate.

Ice cores extracted from the Antarctic ice sheet suggest that glacial conditions, and the relationship between isotopically derived temperatures and atmospheric PCO(2) have been constant over the last 800,000 years of the Late Pleistocene epoch. But independent lines of evidence, such as the extent of Northern Hemisphere ice sheets, sea level and other temperature records, point towards a fluctuating severity of glacial periods, particularly during the more extreme glacial stadials centred around 340,000 and 420,000 years ago (marine isotope stages 10 and 12). Previously unidentified mechanisms therefore appear to have mediated the relationship between insolation, CO(2) and climate. Here we test whether northward migration of the subtropical front (STF) off the southeastern coast of South Africa acts as a gatekeeper for the Agulhas current, which controls the transport of heat and salt from the Indo-Pacific Ocean to the Atlantic Ocean. Using a new 800,000-year record of sea surface temperature and ocean productivity from ocean sediment core MD962077, we demonstrate that during cold stadials (particularly marine isotope stages 10 and 12), productivity peaked and sea surface temperature was up to 6 degrees C cooler than modern temperatures. This suggests that during these cooler stadials, the STF moved northward by up to 7 degrees latitude, nearly shutting off the Agulhas current. Our results, combined with faunal assemblages from the south Atlantic show that variable northwards migration of the Southern Hemisphere STF can modulate the severity of each glacial period by altering the strength of the Agulhas current carrying heat and salt to the Atlantic meridional overturning circulation. We show hence that the degree of northwards migration of the STF can partially decouple global climate from atmospheric partial pressure of carbon dioxide, P CO(2), and help to resolve the long-standing puzzle of differing glacial amplitudes within a consistent range of atmospheric PCO(2).

Moisture transport across Central America as a positive feedback on abrupt climatic changes.

Moisture transport from the Atlantic to the Pacific ocean across Central America leads to relatively high salinities in the North Atlantic Ocean and contributes to the formation of North Atlantic Deep Water. This deep water formation varied strongly between Dansgaard/Oeschger interstadials and Heinrich events-millennial-scale abrupt warm and cold events, respectively, during the last glacial period. Increases in the moisture transport across Central America have been proposed to coincide with northerly shifts of the Intertropical Convergence Zone and with Dansgaard/Oeschger interstadials, with opposite changes for Heinrich events. Here we reconstruct sea surface salinities in the eastern equatorial Pacific Ocean over the past 90,000 years by comparing palaeotemperature estimates from alkenones and Mg/Ca ratios with foraminiferal oxygen isotope ratios that vary with both temperature and salinity. We detect millennial-scale fluctuations of sea surface salinities in the eastern equatorial Pacific Ocean of up to two to four practical salinity units. High salinities are associated with the southward migration of the tropical Atlantic Intertropical Convergence Zone, coinciding with Heinrich events and with Greenland stadials. The amplitudes of these salinity variations are significantly larger on the Pacific side of the Panama isthmus, as inferred from a comparison of our data with a palaeoclimate record from the Caribbean basin. We conclude that millennial-scale fluctuations of moisture transport constitute an important feedback mechanism for abrupt climate changes, modulating the North Atlantic freshwater budget and hence North Atlantic Deep Water formation.

Chronological and genetic analysis of an Upper Palaeolithic female infant burial from Borsuka Cave, Poland.

Six infant human teeth and 112 animal tooth pendants from Borsuka Cave were identified as the oldest burial in Poland. However, uncertainties around the dating and the association of the teeth to the pendants have precluded their association with an Upper Palaeolithic archaeological industry. Using <67 mg per tooth, we combined dating and genetic analyses of two human teeth and six herbivore tooth pendants to address these questions. Our interdisciplinary approach yielded informative results despite limited sampling material, and high levels of degradation and contamination. Our results confirm the Palaeolithic origin of the human remains and herbivore pendants, and permit us to identify the infant as female and discuss the association of the assemblage with different Palaeolithic industries. This study exemplifies the progress that has been made toward minimally destructive methods and the benefits of integrating methods to maximize data retrieval from precious but highly degraded and contaminated prehistoric material.

Tracing the mobility of a Late Epigravettian (~ 13 ka) male infant from Grotte di Pradis (Northeastern Italian Prealps) at high-temporal resolution.

We present the results of a multi-disciplinary investigation on a deciduous human tooth (Pradis 1), recently recovered from the Epigravettian layers of the Grotte di Pradis archaeological site (Northeastern Italian Prealps). Pradis 1 is an exfoliated deciduous molar (Rdm2), lost during life by an 11-12-year-old child. A direct radiocarbon date provided an age of 13,088-12,897 cal BP (95% probability, IntCal20). Amelogenin peptides extracted from tooth enamel and analysed through LC-MS/MS indicate that Pradis 1 likely belonged to a male. Time-resolved 87Sr/86Sr analyses by laser ablation mass spectrometry (LA-MC-ICPMS), combined with dental histology, were able to resolve his movements during the first year of life (i.e. the enamel mineralization interval). Specifically, the Sr isotope ratio of the tooth enamel differs from the local baseline value, suggesting that the child likely spent his first year of life far from Grotte di Pradis. Sr isotopes are also suggestive of a cyclical/seasonal mobility pattern exploited by the Epigravettian human group. The exploitation of Grotte di Pradis on a seasonal, i.e. summer, basis is also indicated by the faunal spectra. Indeed, the nearly 100% occurrence of marmot remains in the entire archaeozoological collection indicates the use of Pradis as a specialized marmot hunting or butchering site. This work represents the first direct assessment of sub-annual movements observed in an Epigravettian hunter-gatherer group from Northern Italy.

A predominantly tropical influence on late Holocene hydroclimate variation in the hyperarid central Sahara.

The climate history of the Sahara desert during recent millennia is obscured by the near absence of natural climate archives, hampering insight in the relative importance of southerly (tropical) and northerly (midlatitude) weather systems at submillennial time scales. A new lake sediment record from Ounianga Serir oasis in northern Chad, spanning the Late Holocene without interruption, confirms that immediately before ca 4200 years ago, the Sahara experienced an episode of hyperaridity even more extreme than today's desert climate. The hypersaline terminal lake which formed afterwards never desiccated during the late Holocene due to continuous inflow of fossil groundwater, yet its water balance was sensitive to temporal variation in local rainfall and lake surface evaporation. Our in-lake geochemical proxies show that, during the last 3000 years, century-scale hydroclimate variation in the central Sahara primarily tracked the intensity of the tropical West African monsoon, modulated at shorter time scales by weather patterns linked to shifts in midlatitude Atlantic Ocean circulation.

Shut down of the South American summer monsoon during the penultimate glacial.

We analysed changes in mean annual air temperature (MAAT), vegetation and biomass burning on a long and continuous lake-peat sediment record from the Colônia basin, southeastern Brazil, examining the responses of a wet tropical rainforest over the last 180 ka. Stronger southern atmospheric circulation up to the latitude of Colônia was found for the penultimate glacial with lower temperatures than during the last glacial, while strengthening of the South American summer monsoon (SASM) circulation started during the last interglacial and progressively enhanced a longer wet summer season from 95 ka until the present. Past MAAT variations and fire history were possibly modulated by eccentricity, although with signatures which differ in average and in amplitude between the last 180 ka. Vegetation responses were driven by the interplay between the SASM and southern circulation linked to Antarctic ice volume, inferred by the presence of a cool mixed evergreen forest from 180 to 45 ka progressively replaced by a rainforest. We report cooler temperatures during the marine isotope stage 3 (MIS 3: 57-29 ka) than during the Last Glacial Maximum (LGM: 23-19 ka). Our findings show that tropical forest dynamics display different patterns than mid-latitude during the last 180 ka.

A 14C chronology for the Middle to Upper Palaeolithic transition at Bacho Kiro Cave, Bulgaria.

The stratigraphy at Bacho Kiro Cave, Bulgaria, spans the Middle to Upper Palaeolithic transition, including an Initial Upper Palaeolithic (IUP) assemblage argued to represent the earliest arrival of Upper Palaeolithic Homo sapiens in Europe. We applied the latest techniques in 14C dating to an extensive dataset of newly excavated animal and human bones to produce a robust, high-precision radiocarbon chronology for the site. At the base of the stratigraphy, the Middle Palaeolithic (MP) occupation dates to >51,000 yr BP. A chronological gap of over 3,000 years separates the MP occupation from the occupation of the cave by H. sapiens, which extends to 34,000 cal BP. The extensive IUP assemblage, now associated with directly dated H. sapiens fossils at this site, securely dates to 45,820-43,650 cal BP (95.4% probability), probably beginning from 46,940 cal BP (95.4% probability). The results provide chronological context for the early occupation of Europe by Upper Palaeolithic H. sapiens.

Liquid chromatographic isolation of individual carbohydrates from environmental matrices for stable carbon analysis and radiocarbon dating.

Carbohydrates are among the most abundant organic molecules in both aquatic and terrestrial ecosystems; however, very few studies have addressed their isotopic signature using compound-specific isotope analysis, which provides additional information on their origin (δ13C) and fate (Δ14C). In this study, semi-preparative liquid chromatography with refractive index detection (HPLC-RI) was employed to produce pure carbohydrate targets for subsequent offline δ13C and Δ14C isotopic analysis. δ13C analysis was performed by elemental analyzer-isotope ratio mass spectrometer (EA-IRMS) whereas Δ14C analysis was performed by an innovative measurement procedure based on the direct combustion of the isolated fractions using an elemental analyzer coupled to the gas source of a mini carbon dating system (AixMICADAS). In general, four successive purifications with Na+, Ca2+, Pb2+, and Ca2+ cation-exchange columns were sufficient to produce pure carbohydrates. These carbohydrates were subsequently identified using mass spectrometry by comparing their mass spectra with those of authentic standards. The applicability of the proposed method was tested on two different environmental samples comprising marine particulate organic matter (POM) and total suspended atmospheric particles (TSP). The obtained results revealed that for the marine POM sample, the δ13C values of the individual carbohydrates ranged from -18.5 to -16.8‰, except for levoglucosan and mannosan, which presented values of -27.2 and -26.2‰, respectively. For the TSP sample, the δ13C values ranged from -26.4 to -25.0‰. The galactose and glucose Δ14C values were 19 and 43‰, respectively, for the POM sample. On the other hand, the levoglucosan radiocarbon value was 33‰ for the TSP sample. These results suggest that these carbohydrates exhibit a modern age in both of these samples. Radiocarbon HPLC collection window blanks, measured after the addition of phthalic acid (14C free blank), ranged from -988 to -986‰ for the abovementioned compounds, indicating a very small background isotopic influence from the whole purification procedure. Overall, the proposed method does not require derivatization steps, produces extremely low blanks, and may be applied to different types of environmental samples.

Permafrost thawing as a possible source of abrupt carbon release at the onset of the Bølling/Allerød.

One of the most abrupt and yet unexplained past rises in atmospheric CO2 (>10 p.p.m.v. in two centuries) occurred in quasi-synchrony with abrupt northern hemispheric warming into the Bølling/Allerød, ~14,600 years ago. Here we use a U/Th-dated record of atmospheric Δ(14)C from Tahiti corals to provide an independent and precise age control for this CO2 rise. We also use model simulations to show that the release of old (nearly (14)C-free) carbon can explain these changes in CO2 and Δ(14)C. The Δ(14)C record provides an independent constraint on the amount of carbon released (~125 Pg C). We suggest, in line with observations of atmospheric CH4 and terrigenous biomarkers, that thawing permafrost in high northern latitudes could have been the source of carbon, possibly with contribution from flooding of the Siberian continental shelf during meltwater pulse 1A. Our findings highlight the potential of the permafrost carbon reservoir to modulate abrupt climate changes via greenhouse-gas feedbacks.

Pronounced interannual variability in tropical South Pacific temperatures during Heinrich Stadial 1.

The early last glacial termination was characterized by intense North Atlantic cooling and weak overturning circulation. This interval between ~18,000 and 14,600 years ago, known as Heinrich Stadial 1, was accompanied by a disruption of global climate and has been suggested as a key factor for the termination. However, the response of interannual climate variability in the tropical Pacific (El Niño-Southern Oscillation) to Heinrich Stadial 1 is poorly understood. Here we use Sr/Ca in a fossil Tahiti coral to reconstruct tropical South Pacific sea surface temperature around 15,000 years ago at monthly resolution. Unlike today, interannual South Pacific sea surface temperature variability at typical El Niño-Southern Oscillation periods was pronounced at Tahiti. Our results indicate that the El Niño-Southern Oscillation was active during Heinrich Stadial 1, consistent with climate model simulations of enhanced El Niño-Southern Oscillation variability at that time. Furthermore, a greater El Niño-Southern Oscillation influence in the South Pacific during Heinrich Stadial 1 is suggested, resulting from a southward expansion or shift of El Niño-Southern Oscillation sea surface temperature anomalies.

Deglacial meltwater pulse 1B and Younger Dryas sea levels revisited with boreholes at Tahiti.

Reconstructing sea-level changes during the last deglaciation provides a way of understanding the ice dynamics that can perturb large continental ice sheets. The resolution of the few sea-level records covering the critical time interval between 14,000 and 9,000 calendar years before the present is still insufficient to draw conclusions about sea-level changes associated with the Younger Dryas cold event and the meltwater pulse 1B (MWP-1B). We used the uranium-thorium method to date shallow-living corals from three new cores drilled onshore in the Tahiti barrier reef. No significant discontinuity can be detected in the sea-level rise during the MWP-1B period. The new Tahiti sea-level record shows that the sea-level rise slowed down during the Younger Dryas before accelerating again during the Holocene.