Climate sensitivity constrained by temperature reconstructions over the past seven centuries

The magnitude and impact of future global warming depends on the sensitivity of the climate system to changes in greenhouse gas concentrations. The commonly accepted range for the equilibrium global mean temperature change in response to a doubling of the atmospheric carbon dioxide concentration, termed climate sensitivity, is 1.5-4.5 K (ref. 2). A number of observational studies, however, find a substantial probability of significantly higher sensitivities, yielding upper limits on climate sensitivity of 7.7 K to above 9 K (refs 3-8). Here we demonstrate that such observational estimates of climate sensitivity can be tightened if reconstructions of Northern Hemisphere temperature over the past several centuries are considered. We use large-ensemble energy balance modelling and simulate the temperature response to past solar, volcanic and greenhouse gas forcing to determine which climate sensitivities yield simulations that are in agreement with proxy reconstructions. After accounting for the uncertainty in reconstructions and estimates of past external forcing, we find an independent estimate of climate sensitivity that is very similar to those from instrumental data. If the latter are combined with the result from all proxy reconstructions, then the 5-95 per cent range shrinks to 1.5-6.2 K, thus substantially reducing the probability of very high climate sensitivity.     

The long-term carbon cycle, fossil fuels and atmospheric composition

The long-term carbon cycle operates over millions of years and involves the exchange of carbon between rocks and the Earth's surface. There are many complex feedback pathways between carbon burial, nutrient cycling, atmospheric carbon dioxide and oxygen, and climate. New calculations of carbon fluxes during the Phanerozoic eon (the past 550 million years) illustrate how the long-term carbon cycle has affected the burial of organic matter and fossil-fuel formation, as well as the evolution of atmospheric composition.     

Eocene bipolar glaciation associated with global carbon cycle changes

The transition from the extreme global warmth of the early Eocene 'greenhouse' climate approximately 55 million years ago to the present glaciated state is one of the most prominent changes in Earth's climatic evolution. It is widely accepted that large ice sheets first appeared on Antarctica approximately 34 million years ago, coincident with decreasing atmospheric carbon dioxide concentrations and a deepening of the calcite compensation depth in the world's oceans, and that glaciation in the Northern Hemisphere began much later, between 10 and 6 million years ago. Here we present records of sediment and foraminiferal geochemistry covering the greenhouse-icehouse climate transition. We report evidence for synchronous deepening and subsequent oscillations in the calcite compensation depth in the tropical Pacific and South Atlantic oceans from approximately 42 million years ago, with a permanent deepening 34 million years ago. The most prominent variations in the calcite compensation depth coincide with changes in seawater oxygen isotope ratios of up to 1.5 per mil, suggesting a lowering of global sea level through significant storage of ice in both hemispheres by at least 100 to 125 metres. Variations in benthic carbon isotope ratios of up to approximately 1.4 per mil occurred at the same time, indicating large changes in carbon cycling. We suggest that the greenhouse-icehouse transition was closely coupled to the evolution of atmospheric carbon dioxide, and that negative carbon cycle feedbacks may have prevented the permanent establishment of large ice sheets earlier than 34 million years ago.     

Net carbon dioxide losses of northern ecosystems in response to autumn warming

The carbon balance of terrestrial ecosystems is particularly sensitive to climatic changes in autumn and spring, with spring and autumn temperatures over northern latitudes having risen by about 1.1 degrees C and 0.8 degrees C, respectively, over the past two decades. A simultaneous greening trend has also been observed, characterized by a longer growing season and greater photosynthetic activity. These observations have led to speculation that spring and autumn warming could enhance carbon sequestration and extend the period of net carbon uptake in the future. Here we analyse interannual variations in atmospheric carbon dioxide concentration data and ecosystem carbon dioxide fluxes. We find that atmospheric records from the past 20 years show a trend towards an earlier autumn-to-winter carbon dioxide build-up, suggesting a shorter net carbon uptake period. This trend cannot be explained by changes in atmospheric transport alone and, together with the ecosystem flux data, suggest increasing carbon losses in autumn. We use a process-based terrestrial biosphere model and satellite vegetation greenness index observations to investigate further the observed seasonal response of northern ecosystems to autumnal warming. We find that both photosynthesis and respiration increase during autumn warming, but the increase in respiration is greater. In contrast, warming increases photosynthesis more than respiration in spring. Our simulations and observations indicate that northern terrestrial ecosystems may currently lose carbon dioxide in response to autumn warming, with a sensitivity of about 0.2 PgC degrees C(-1), offsetting 90% of the increased carbon dioxide uptake during spring. If future autumn warming occurs at a faster rate than in spring, the ability of northern ecosystems to sequester carbon may be diminished earlier than previously suggested.     

Stable sea surface temperatures in the western Pacific warm pool over the past 1.75 million years

About 850,000 years ago, the period of the glacial cycles changed from 41,000 to 100,000 years. This mid-Pleistocene climate transition has been attributed to global cooling, possibly caused by a decrease in atmospheric carbon dioxide concentrations. However, evidence for such cooling is currently restricted to the cool upwelling regions in the eastern equatorial oceans, although the tropical warm pools on the western side of the ocean basins are particularly sensitive to changes in radiative forcing. Here we present high-resolution records of sea surface temperatures spanning the past 1.75 million years, obtained from oxygen isotopes and Mg/Ca ratios in planktonic foraminifera from the western Pacific warm pool. In contrast with the eastern equatorial regions, sea surface temperatures in the western Pacific warm pool are relatively stable throughout the Pleistocene epoch, implying little long-term change in the tropical net radiation budget. Our results challenge the hypothesis of a gradual decrease in atmospheric carbon dioxide concentrations as a dominant trigger of the longer glacial cycles since 850,000 years ago. Instead, we infer that the temperature contrast across the equatorial Pacific Ocean increased, which might have had a significant influence on the mid-Pleistocene climate transition.     

Negligible glacial-interglacial variation in continental chemical weathering rates

Chemical weathering of the continents is central to the regulation of atmospheric carbon dioxide concentrations, and hence global climate. On million-year timescales silicate weathering leads to the draw-down of carbon dioxide, and on millennial timescales chemical weathering affects the calcium carbonate saturation state of the oceans and hence their uptake of carbon dioxide. However, variations in chemical weathering rates over glacial-interglacial cycles remain uncertain. During glacial periods, cold and dry conditions reduce the rate of chemical weathering, but intense physical weathering and the exposure of carbonates on continental shelves due to low sea levels may increase this rate. Here we present high-resolution records of the lead isotope composition of ferromanganese crusts from the North Atlantic Ocean that cover the past 550,000 years. Combining these records with a simple quantitative model of changes in the lead isotope composition of the deep North Atlantic Ocean in response to chemical weathering, we find that chemical weathering rates were two to three times lower in the glaciated interior of the North Atlantic Region during glacial periods than during the intervening interglacial periods. This decrease roughly balances the increase in chemical weathering caused by the exposure of continental shelves, indicating that chemical weathering rates remained relatively constant on glacial-interglacial timescales. On timescales of more than a million years, however, we suggest that enhanced weathering of silicate glacial sediments during interglacial periods results in a net draw-down of atmospheric carbon dioxide, creating a positive feedback on global climate that, once initiated, promotes cooling and further glaciation.     

Triple oxygen isotope evidence for elevated CO2 levels after a Neoproterozoic glaciation

Understanding the composition of the atmosphere over geological time is critical to understanding the history of the Earth system, as the atmosphere is closely linked to the lithosphere, hydrosphere and biosphere. Although much of the history of the lithosphere and hydrosphere is contained in rock and mineral records, corresponding information about the atmosphere is scarce and elusive owing to the lack of direct records. Geologists have used sedimentary minerals, fossils and geochemical models to place constraints on the concentrations of carbon dioxide, oxygen or methane in the past. Here we show that the triple oxygen isotope composition of sulphate from ancient evaporites and barites shows variable negative oxygen-17 isotope anomalies over the past 750 million years. We propose that these anomalies track those of atmospheric oxygen and in turn reflect the partial pressure of carbon dioxide (P(CO2)) in the past through a photochemical reaction network linking stratospheric ozone to carbon dioxide and to oxygen. Our results suggest that P(CO2) was much higher in the early Cambrian than in younger eras, agreeing with previous modelling results. We also find that the (17)O isotope anomalies of barites from Marinoan (approximately 635 million years ago) cap carbonates display a distinct negative spike (around -0.70 per thousand), suggesting that by the time barite was precipitating in the immediate aftermath of a Neoproterozoic global glaciation, the P(CO2) was at its highest level in the past 750 million years. Our finding is consistent with the 'snowball Earth' hypothesis and/or a massive methane release after the Marinoan glaciation.     

Impacts of orbital forcing and atmospheric carbon dioxide on Miocene ice-sheet expansion

The processes causing the middle Miocene global cooling, which marked the Earth's final transition into an 'icehouse' climate about 13.9 million years ago (Myr ago), remain enigmatic. Tectonically driven circulation changes and variations in atmospheric carbon dioxide levels have been suggested as driving mechanisms, but the lack of adequately preserved sedimentary successions has made rigorous testing of these hypotheses difficult. Here we present high-resolution climate proxy records, covering the period from 14.7 to 12.7 million years ago, from two complete sediment cores from the northwest and southeast subtropical Pacific Ocean. Using new chronologies through the correlation to the latest orbital model, we find relatively constant, low summer insolation over Antarctica coincident with declining atmospheric carbon dioxide levels at the time of Antarctic ice-sheet expansion and global cooling, suggesting a causal link. We surmise that the thermal isolation of Antarctica played a role in providing sustained long-term climatic boundary conditions propitious for ice-sheet formation. Our data document that Antarctic glaciation was rapid, taking place within two obliquity cycles, and coincided with a striking transition from obliquity to eccentricity as the drivers of climatic change.     

High-resolution carbon dioxide concentration record 650,000-800,000 years before present

Changes in past atmospheric carbon dioxide concentrations can be determined by measuring the composition of air trapped in ice cores from Antarctica. So far, the Antarctic Vostok and EPICA Dome C ice cores have provided a composite record of atmospheric carbon dioxide levels over the past 650,000 years. Here we present results of the lowest 200 m of the Dome C ice core, extending the record of atmospheric carbon dioxide concentration by two complete glacial cycles to 800,000 yr before present. From previously published data and the present work, we find that atmospheric carbon dioxide is strongly correlated with Antarctic temperature throughout eight glacial cycles but with significantly lower concentrations between 650,000 and 750,000 yr before present. Carbon dioxide levels are below 180 parts per million by volume (p.p.m.v.) for a period of 3,000 yr during Marine Isotope Stage 16, possibly reflecting more pronounced oceanic carbon storage. We report the lowest carbon dioxide concentration measured in an ice core, which extends the pre-industrial range of carbon dioxide concentrations during the late Quaternary by about 10 p.p.m.v. to 172-300 p.p.m.v.     

Persistent near-tropical warmth on the Antarctic continent during the early Eocene epoch

The warmest global climates of the past 65 million years occurred during the early Eocene epoch (about 55 to 48 million years ago), when the Equator-to-pole temperature gradients were much smaller than today and atmospheric carbon dioxide levels were in excess of one thousand parts per million by volume. Recently the early Eocene has received considerable interest because it may provide insight into the response of Earth's climate and biosphere to the high atmospheric carbon dioxide levels that are expected in the near future as a consequence of unabated anthropogenic carbon emissions. Climatic conditions of the early Eocene 'greenhouse world', however, are poorly constrained in critical regions, particularly Antarctica. Here we present a well-dated record of early Eocene climate on Antarctica from an ocean sediment core recovered off the Wilkes Land coast of East Antarctica. The information from biotic climate proxies (pollen and spores) and independent organic geochemical climate proxies (indices based on branched tetraether lipids) yields quantitative, seasonal temperature reconstructions for the early Eocene greenhouse world on Antarctica. We show that the climate in lowland settings along the Wilkes Land coast (at a palaeolatitude of about 70° south) supported the growth of highly diverse, near-tropical forests characterized by mesothermal to megathermal floral elements including palms and Bombacoideae. Notably, winters were extremely mild (warmer than 10?°C) and essentially frost-free despite polar darkness, which provides a critical new constraint for the validation of climate models and for understanding the response of high-latitude terrestrial ecosystems to increased carbon dioxide forcing.     

Atmospheric carbon dioxide concentrations over the past 60 million years

Knowledge of the evolution of atmospheric carbon dioxide concentrations throughout the Earth's history is important for a reconstruction of the links between climate and radiative forcing of the Earth's surface temperatures. Although atmospheric carbon dioxide concentrations in the early Cenozoic era (about 60 Myr ago) are widely believed to have been higher than at present, there is disagreement regarding the exact carbon dioxide levels, the timing of the decline and the mechanisms that are most important for the control of CO2 concentrations over geological timescales. Here we use the boron-isotope ratios of ancient planktonic foraminifer shells to estimate the pH of surface-layer sea water throughout the past 60 million years, which can be used to reconstruct atmospheric CO2 concentrations. We estimate CO2 concentrations of more than 2,000 p.p.m. for the late Palaeocene and earliest Eocene periods (from about 60 to 52 Myr ago), and find an erratic decline between 55 and 40 Myr ago that may have been caused by reduced CO2 outgassing from ocean ridges, volcanoes and metamorphic belts and increased carbon burial. Since the early Miocene (about 24 Myr ago), atmospheric CO2 concentrations appear to have remained below 500 p.p.m. and were more stable than before, although transient intervals of CO2 reduction may have occurred during periods of rapid cooling approximately 15 and 3 Myr ago.     

Coupling of surface temperatures and atmospheric CO2 concentrations during the Palaeozoic era

Atmospheric carbon dioxide concentrations seem to have been several times modern levels during much of the Palaeozoic era (543-248 million years ago), but decreased during the Carboniferous period to concentrations similar to that of today. Given that carbon dioxide is a greenhouse gas, it has been proposed that surface temperatures were significantly higher during the earlier portions of the Palaeozoic era. A reconstruction of tropical sea surface temperatures based on the delta18O of carbonate fossils indicates, however, that the magnitude of temperature variability throughout this period was small, suggesting that global climate may be independent of variations in atmospheric carbon dioxide concentration. Here we present estimates of sea surface temperatures that were obtained from fossil brachiopod and mollusc shells using the 'carbonate clumped isotope' method-an approach that, unlike the delta18O method, does not require independent estimates of the isotopic composition of the Palaeozoic ocean. Our results indicate that tropical sea surface temperatures were significantly higher than today during the Early Silurian period (443-423 Myr ago), when carbon dioxide concentrations are thought to have been relatively high, and were broadly similar to today during the Late Carboniferous period (314-300 Myr ago), when carbon dioxide concentrations are thought to have been similar to the present-day value. Our results are consistent with the proposal that increased atmospheric carbon dioxide concentrations drive or amplify increased global temperatures.     

Large temperature drop across the Eocene-Oligocene transition in central North America

The Eocene-Oligocene transition towards a cool climate (approximately 33.5 million years ago) was one of the most pronounced climate events during the Cenozoic era. The marine record of this transition has been extensively studied. However, significantly less research has focused on continental climate change at the time, yielding partly inconsistent results on the magnitude and timing of the changes. Here we use a combination of in vivo stable isotope compositions of fossil tooth enamel with diagenetic stable isotope compositions of fossil bone to derive a high-resolution (about 40,000 years) continental temperature record for the Eocene-Oligocene transition. We find a large drop in mean annual temperature of 8.2 +/- 3.1 degrees C over about 400,000 years, the possibility of a small increase in temperature seasonality, and no resolvable change in aridity across the transition. The large change in mean annual temperature, exceeding changes in sea surface temperatures at comparable latitudes and possibly delayed in time with respect to marine changes by up to 400,000 years, explains the faunal turnover for gastropods, amphibians and reptiles, whereas most mammals in the region were unaffected. Our results are in agreement with modelling studies that attribute the climate cooling at the Eocene-Oligocene transition to a significant drop in atmospheric carbon dioxide concentrations.     

A 'snowball Earth' climate triggered by continental break-up through changes in runoff

Geological and palaeomagnetic studies indicate that ice sheets may have reached the Equator at the end of the Proterozoic eon, 800 to 550 million years ago, leading to the suggestion of a fully ice-covered 'snowball Earth'. Climate model simulations indicate that such a snowball state for the Earth depends on anomalously low atmospheric carbon dioxide concentrations, in addition to the Sun being 6 per cent fainter than it is today. However, the mechanisms producing such low carbon dioxide concentrations remain controversial. Here we assess the effect of the palaeogeographic changes preceding the Sturtian glacial period, 750 million years ago, on the long-term evolution of atmospheric carbon dioxide levels using the coupled climate-geochemical model GEOCLIM. In our simulation, the continental break-up of Rodinia leads to an increase in runoff and hence consumption of carbon dioxide through continental weathering that decreases atmospheric carbon dioxide concentrations by 1,320 p.p.m. This indicates that tectonic changes could have triggered a progressive transition from a 'greenhouse' to an 'icehouse' climate during the Neoproterozoic era. When we combine these results with the concomitant weathering effect of the voluminous basaltic traps erupted throughout the break-up of Rodinia, our simulation results in a snowball glaciation.     

The influence of rivers on marine boron isotopes and implications for reconstructing past ocean pH

Ocean pH is particularly sensitive to atmospheric carbon dioxide content. Records of ocean pH can therefore be used to estimate past atmospheric carbon dioxide concentrations. The isotopic composition of boron (delta11B) contained in the carbonate shells of marine organisms varies according to pH, from which ocean pH can be reconstructed. This requires independent estimates of the delta11B of dissolved boron in sea water through time. The marine delta11B budget, however, is still largely unconstrained. Here we show that, by incorporating the global flux of riverine boron (as estimated from delta11B measurements in 22 of the world's main rivers), the marine boron isotope budget can be balanced. We also derive ocean delta11B budgets for the past 120 Myr. Estimated isotope compositions of boron in sea water show a remarkable consistency with records of delta11B in foraminiferal carbonates, suggesting that foraminifera delta11B records may in part reflect changes in the marine boron isotope budget rather than changes in ocean pH over the Cenozoic era.     

Effect of natural iron fertilization on carbon sequestration in the Southern Ocean

The availability of iron limits primary productivity and the associated uptake of carbon over large areas of the ocean. Iron thus plays an important role in the carbon cycle, and changes in its supply to the surface ocean may have had a significant effect on atmospheric carbon dioxide concentrations over glacial-interglacial cycles. To date, the role of iron in carbon cycling has largely been assessed using short-term iron-addition experiments. It is difficult, however, to reliably assess the magnitude of carbon export to the ocean interior using such methods, and the short observational periods preclude extrapolation of the results to longer timescales. Here we report observations of a phytoplankton bloom induced by natural iron fertilization--an approach that offers the opportunity to overcome some of the limitations of short-term experiments. We found that a large phytoplankton bloom over the Kerguelen plateau in the Southern Ocean was sustained by the supply of iron and major nutrients to surface waters from iron-rich deep water below. The efficiency of fertilization, defined as the ratio of the carbon export to the amount of iron supplied, was at least ten times higher than previous estimates from short-term blooms induced by iron-addition experiments. This result sheds new light on the effect of long-term fertilization by iron and macronutrients on carbon sequestration, suggesting that changes in iron supply from below--as invoked in some palaeoclimatic and future climate change scenarios--may have a more significant effect on atmospheric carbon dioxide concentrations than previously thought.     

Projected increase in continental runoff due to plant responses to increasing carbon dioxide

In addition to influencing climatic conditions directly through radiative forcing, increasing carbon dioxide concentration influences the climate system through its effects on plant physiology. Plant stomata generally open less widely under increased carbon dioxide concentration, which reduces transpiration and thus leaves more water at the land surface. This driver of change in the climate system, which we term 'physiological forcing', has been detected in observational records of increasing average continental runoff over the twentieth century. Here we use an ensemble of experiments with a global climate model that includes a vegetation component to assess the contribution of physiological forcing to future changes in continental runoff, in the context of uncertainties in future precipitation. We find that the physiological effect of doubled carbon dioxide concentrations on plant transpiration increases simulated global mean runoff by 6 per cent relative to pre-industrial levels; an increase that is comparable to that simulated in response to radiatively forced climate change (11 +/- 6 per cent). Assessments of the effect of increasing carbon dioxide concentrations on the hydrological cycle that only consider radiative forcing will therefore tend to underestimate future increases in runoff and overestimate decreases. This suggests that freshwater resources may be less limited than previously assumed under scenarios of future global warming, although there is still an increased risk of drought. Moreover, our results highlight that the practice of assessing the climate-forcing potential of all greenhouse gases in terms of their radiative forcing potential relative to carbon dioxide does not accurately reflect the relative effects of different greenhouse gases on freshwater resources.     

Southern Ocean dust-climate coupling over the past four million years

Dust has the potential to modify global climate by influencing the radiative balance of the atmosphere and by supplying iron and other essential limiting micronutrients to the ocean. Indeed, dust supply to the Southern Ocean increases during ice ages, and 'iron fertilization' of the subantarctic zone may have contributed up to 40?parts per million by volume (p.p.m.v.) of the decrease (80-100 p.p.m.v.) in atmospheric carbon dioxide observed during late Pleistocene glacial cycles. So far, however, the magnitude of Southern Ocean dust deposition in earlier times and its role in the development and evolution of Pleistocene glacial cycles have remained unclear. Here we report a high-resolution record of dust and iron supply to the Southern Ocean over the past four million years, derived from the analysis of marine sediments from ODP Site 1090, located in the Atlantic sector of the subantarctic zone. The close correspondence of our dust and iron deposition records with Antarctic ice core reconstructions of dust flux covering the past 800,000 years (refs 8, 9) indicates that both of these archives record large-scale deposition changes that should apply to most of the Southern Ocean, validating previous interpretations of the ice core data. The extension of the record beyond the interval covered by the Antarctic ice cores reveals that, in contrast to the relatively gradual intensification of glacial cycles over the past three million years, Southern Ocean dust and iron flux rose sharply at the Mid-Pleistocene climatic transition around 1.25 million years ago. This finding complements previous observations over late Pleistocene glacial cycles, providing new evidence of a tight connection between high dust input to the Southern Ocean and the emergence of the deep glaciations that characterize the past one million years of Earth history.     

Quantification of modelling uncertainties in a large ensemble of climate change simulations

Comprehensive global climate models are the only tools that account for the complex set of processes which will determine future climate change at both a global and regional level. Planners are typically faced with a wide range of predicted changes from different models of unknown relative quality, owing to large but unquantified uncertainties in the modelling process. Here we report a systematic attempt to determine the range of climate changes consistent with these uncertainties, based on a 53-member ensemble of model versions constructed by varying model parameters. We estimate a probability density function for the sensitivity of climate to a doubling of atmospheric carbon dioxide levels, and obtain a 5-95 per cent probability range of 2.4-5.4 degrees C. Our probability density function is constrained by objective estimates of the relative reliability of different model versions, the choice of model parameters that are varied and their uncertainty ranges, specified on the basis of expert advice. Our ensemble produces a range of regional changes much wider than indicated by traditional methods based on scaling the response patterns of an individual simulation.     

Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model

The continued increase in the atmospheric concentration of carbon dioxide due to anthropogenic emissions is predicted to lead to significant changes in climate. About half of the current emissions are being absorbed by the ocean and by land ecosystems, but this absorption is sensitive to climate as well as to atmospheric carbon dioxide concentrations, creating a feedback loop. General circulation models have generally excluded the feedback between climate and the biosphere, using static vegetation distributions and CO2 concentrations from simple carbon-cycle models that do not include climate change. Here we present results from a fully coupled, three-dimensional carbon-climate model, indicating that carbon-cycle feedbacks could significantly accelerate climate change over the twenty-first century. We find that under a 'business as usual' scenario, the terrestrial biosphere acts as an overall carbon sink until about 2050, but turns into a source thereafter. By 2100, the ocean uptake rate of 5 Gt C yr(-1) is balanced by the terrestrial carbon source, and atmospheric CO2 concentrations are 250 p.p.m.v. higher in our fully coupled simulation than in uncoupled carbon models, resulting in a global-mean warming of 5.5 K, as compared to 4 K without the carbon-cycle feedback.