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.     

Neoproterozoic 'snowball Earth' simulations with a coupled climate/ice-sheet model

Ice sheets may have reached the Equator in the late Proterozoic era (600-800 Myr ago), according to geological and palaeomagnetic studies, possibly resulting in a 'snowball Earth'. But this period was a critical time in the evolution of multicellular animals, posing the question of how early life survived under such environmental stress. Here we present computer simulations of this unusual climate stage with a coupled climate/ice-sheet model. To simulate a snowball Earth, we use only a reduction in the solar constant compared to present-day conditions and we keep atmospheric CO2 concentrations near present levels. We find rapid transitions into and out of full glaciation that are consistent with the geological evidence. When we combine these results with a general circulation model, some of the simulations result in an equatorial belt of open water that may have provided a refugium for multicellular animals.     

Snowball Earth prevention by dissolved organic carbon remineralization

The 'snowball Earth' hypothesis posits the occurrence of a sequence of glaciations in the Earth's history sufficiently deep that photosynthetic activity was essentially arrested. Because the time interval during which these events are believed to have occurred immediately preceded the Cambrian explosion of life, the issue as to whether such snowball states actually developed has important implications for our understanding of evolutionary biology. Here we couple an explicit model of the Neoproterozoic carbon cycle to a model of the physical climate system. We show that the drawdown of atmospheric oxygen into the ocean, as surface temperatures decline, operates so as to increase the rate of remineralization of a massive pool of dissolved organic carbon. This leads directly to an increase of atmospheric carbon dioxide, enhanced greenhouse warming of the surface of the Earth, and the prevention of a snowball state.     

Stable isotopic evidence for methane seeps in Neoproterozoic postglacial cap carbonates

The Earth's most severe glaciations are thought to have occurred about 600 million years ago, in the late Neoproterozoic era. A puzzling feature of glacial deposits from this interval is that they are overlain by 1-5-m-thick 'cap carbonates' (particulate deep-water marine carbonate rocks) associated with a prominent negative carbon isotope excursion. Cap carbonates have been controversially ascribed to the aftermath of almost complete shutdown of the ocean ecosystems for millions of years during such ice ages--the 'snowball Earth' hypothesis. Conversely, it has also been suggested that these carbonate rocks were the result of destabilization of methane hydrates during deglaciation and concomitant flooding of continental shelves and interior basins. The most compelling criticism of the latter 'methane hydrate' hypothesis has been the apparent lack of extreme isotopic variation in cap carbonates inferred locally to be associated with methane seeps. Here we report carbon isotopic and petrographic data from a Neoproterozoic postglacial cap carbonate in south China that provide direct evidence for methane-influenced processes during deglaciation. This evidence lends strong support to the hypothesis that methane hydrate destabilization contributed to the enigmatic cap carbonate deposition and strongly negative carbon isotopic anomalies following Neoproterozoic ice ages. This explanation requires less extreme environmental disturbance than that implied by the snowball Earth hypothesis.     

Carbon dioxide release from the North Pacific abyss during the last deglaciation

Atmospheric carbon dioxide concentrations were significantly lower during glacial periods than during intervening interglacial periods, but the mechanisms responsible for this difference remain uncertain. Many recent explanations call on greater carbon storage in a poorly ventilated deep ocean during glacial periods, but direct evidence regarding the ventilation and respired carbon content of the glacial deep ocean is sparse and often equivocal. Here we present sedimentary geochemical records from sites spanning the deep subarctic Pacific that--together with previously published results--show that a poorly ventilated water mass containing a high concentration of respired carbon dioxide occupied the North Pacific abyss during the Last Glacial Maximum. Despite an inferred increase in deep Southern Ocean ventilation during the first step of the deglaciation (18,000-15,000 years ago), we find no evidence for improved ventilation in the abyssal subarctic Pacific until a rapid transition approximately 14,600 years ago: this change was accompanied by an acceleration of export production from the surface waters above but only a small increase in atmospheric carbon dioxide concentration. We speculate that these changes were mechanistically linked to a roughly coeval increase in deep water formation in the North Atlantic, which flushed respired carbon dioxide from northern abyssal waters, but also increased the supply of nutrients to the upper ocean, leading to greater carbon dioxide sequestration at mid-depths and stalling the rise of atmospheric carbon dioxide concentrations. Our findings are qualitatively consistent with hypotheses invoking a deglacial flushing of respired carbon dioxide from an isolated, deep ocean reservoir, but suggest that the reservoir may have been released in stages, as vigorous deep water ventilation switched between North Atlantic and Southern Ocean source regions.     

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.     

High CO2 levels in the Proterozoic atmosphere estimated from analyses of individual microfossils

Solar luminosity on the early Earth was significantly lower than today. Therefore, solar luminosity models suggest that, in the atmosphere of the early Earth, the concentration of greenhouse gases such as carbon dioxide and methane must have been much higher. However, empirical estimates of Proterozoic levels of atmospheric carbon dioxide concentrations have not hitherto been available. Here we present ion microprobe analyses of the carbon isotopes in individual organic-walled microfossils extracted from a Proterozoic ( approximately 1.4-gigayear-old) shale in North China. Calculated magnitudes of the carbon isotope fractionation in these large, morphologically complex microfossils suggest elevated levels of carbon dioxide in the ancient atmosphere--between 10 and 200 times the present atmospheric level. Our results indicate that carbon dioxide was an important greenhouse gas during periods of lower solar luminosity, probably dominating over methane after the atmosphere and hydrosphere became pervasively oxygenated between 2 and 2.2 gigayears ago.     

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.     

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.     

Evidence for decoupling of atmospheric CO2 and global climate during the Phanerozoic eon

Atmospheric carbon dioxide concentrations are believed to drive climate changes from glacial to interglacial modes, although geological and astronomical mechanisms have been invoked as ultimate causes. Additionally, it is unclear whether the changes between cold and warm modes should be regarded as a global phenomenon, affecting tropical and high-latitude temperatures alike, or if they are better described as an expansion and contraction of the latitudinal climate zones, keeping equatorial temperatures approximately constant. Here we present a reconstruction of tropical sea surface temperatures throughout the Phanerozoic eon (the past approximately 550 Myr) from our database of oxygen isotopes in calcite and aragonite shells. The data indicate large oscillations of tropical sea surface temperatures in phase with the cold-warm cycles, thus favouring the idea of climate variability as a global phenomenon. But our data conflict with a temperature reconstruction using an energy balance model that is forced by reconstructed atmospheric carbon dioxide concentrations. The results can be reconciled if atmospheric carbon dioxide concentrations were not the principal driver of climate variability on geological timescales for at least one-third of the Phanerozoic eon, or if the reconstructed carbon dioxide concentrations are not reliable.     

Quantifying the risk of extreme seasonal precipitation events in a changing climate

Increasing concentrations of atmospheric carbon dioxide will almost certainly lead to changes in global mean climate. But because--by definition--extreme events are rare, it is significantly more difficult to quantify the risk of extremes. Ensemble-based probabilistic predictions, as used in short- and medium-term forecasts of weather and climate, are more useful than deterministic forecasts using a 'best guess' scenario to address this sort of problem. Here we present a probabilistic analysis of 19 global climate model simulations with a generic binary decision model. We estimate that the probability of total boreal winter precipitation exceeding two standard deviations above normal will increase by a factor of five over parts of the UK over the next 100 years. We find similar increases in probability for the Asian monsoon region in boreal summer, with implications for flooding in Bangladesh. Further practical applications of our techniques would be helped by the use of larger ensembles (for a more complete sampling of model uncertainty) and a wider range of scenarios at a resolution adequate to analyse average-size river basins.     

Indirect radiative forcing of climate change through ozone effects on the land-carbon sink

The evolution of the Earth's climate over the twenty-first century depends on the rate at which anthropogenic carbon dioxide emissions are removed from the atmosphere by the ocean and land carbon cycles. Coupled climate-carbon cycle models suggest that global warming will act to limit the land-carbon sink, but these first generation models neglected the impacts of changing atmospheric chemistry. Emissions associated with fossil fuel and biomass burning have acted to approximately double the global mean tropospheric ozone concentration, and further increases are expected over the twenty-first century. Tropospheric ozone is known to damage plants, reducing plant primary productivity and crop yields, yet increasing atmospheric carbon dioxide concentrations are thought to stimulate plant primary productivity. Increased carbon dioxide and ozone levels can both lead to stomatal closure, which reduces the uptake of either gas, and in turn limits the damaging effect of ozone and the carbon dioxide fertilization of photosynthesis. Here we estimate the impact of projected changes in ozone levels on the land-carbon sink, using a global land carbon cycle model modified to include the effect of ozone deposition on photosynthesis and to account for interactions between ozone and carbon dioxide through stomatal closure. For a range of sensitivity parameters based on manipulative field experiments, we find a significant suppression of the global land-carbon sink as increases in ozone concentrations affect plant productivity. In consequence, more carbon dioxide accumulates in the atmosphere. We suggest that the resulting indirect radiative forcing by ozone effects on plants could contribute more to global warming than the direct radiative forcing due to tropospheric ozone increases.     

A carbon isotope challenge to the snowball Earth

The snowball Earth hypothesis postulates that the planet was entirely covered by ice for millions of years in the Neoproterozoic era, in a self-enhanced glaciation caused by the high albedo of the ice-covered planet. In a hard-snowball picture, the subsequent rapid unfreezing resulted from an ultra-greenhouse event attributed to the buildup of volcanic carbon dioxide (CO(2)) during glaciation. High partial pressures of atmospheric CO(2) (pCO2; from 20,000 to 90,000?p.p.m.v.) in the aftermath of the Marinoan glaciation (～635?Myr ago) have been inferred from both boron and triple oxygen isotopes. These pCO2 values are 50 to 225 times higher than present-day levels. Here, we re-evaluate these estimates using paired carbon isotopic data for carbonate layers that cap Neoproterozoic glacial deposits and are considered to record post-glacial sea level rise. The new data reported here for Brazilian cap carbonates, together with previous ones for time-equivalent units, provide estimates lower than 3,200?p.p.m.v.--and possibly as low as the current value of ～400?p.p.m.v. Our new constraint, and our re-interpretation of the boron and triple oxygen isotope data, provide a completely different picture of the late Neoproterozoic environment, with low atmospheric concentrations of carbon dioxide and oxygen that are inconsistent with a hard-snowball Earth.     

A lower limit for atmospheric carbon dioxide levels 3.2 billion years ago

The quantification of greenhouse gases present in the Archaean atmosphere is critical for understanding the evolution of atmospheric oxygen, surface temperatures and the conditions for life on early Earth. For instance, it has been argued that small changes in the balance between two potential greenhouse gases, carbon dioxide and methane, may have dictated the feedback cycle involving organic haze production and global cooling. Climate models have focused on carbon dioxide as the greenhouse gas responsible for maintaining above-freezing surface temperatures during a time of low solar luminosity. However, the analysis of 2.75-billion-year (Gyr)-old palaeosols--soil samples preserved in the geologic record--have recently provided an upper constraint on atmospheric carbon dioxide levels well below that required in most climate models to prevent the Earth's surface from freezing. This finding prompted many to look towards methane as an additional greenhouse gas to satisfy climate models. Here we use model equilibrium reactions for weathering rinds on 3.2-Gyr-old river gravels to show that the presence of iron-rich carbonate relative to common clay minerals requires a minimum partial pressure of carbon dioxide several times higher than present-day values. Unless actual carbon dioxide levels were considerably greater than this, climate models predict that additional greenhouse gases would still need to have a role in maintaining above-freezing surface temperatures.     

Proterozoic low orbital obliquity and axial-dipolar geomagnetic field from evaporite palaeolatitudes

Palaeomagnetism of climatically sensitive sedimentary rock types, such as glacial deposits and evaporites, can test the uniformitarianism of ancient geomagnetic fields and palaeoclimate zones. Proterozoic glacial deposits laid down in near-equatorial palaeomagnetic latitudes can be explained by 'snowball Earth' episodes, high orbital obliquity or markedly non-uniformitarian geomagnetic fields. Here I present a global palaeomagnetic compilation of the Earth's entire basin-scale evaporite record. Magnetic inclinations are consistent with low orbital obliquity and a geocentric-axial-dipole magnetic field for most of the past two billion years, and the snowball Earth hypothesis accordingly remains the most viable model for low-latitude Proterozoic ice ages. Efforts to reconstruct Proterozoic supercontinents are strengthened by this demonstration of a consistently axial and dipolar geomagnetic reference frame, which itself implies stability of geodynamo processes on billion-year timescales.     

Reduction of soil carbon formation by tropospheric ozone under increased carbon dioxide levels

In the Northern Hemisphere, ozone levels in the troposphere have increased by 35 per cent over the past century, with detrimental impacts on forest and agricultural productivity, even when forest productivity has been stimulated by increased carbon dioxide levels. In addition to reducing productivity, increased tropospheric ozone levels could alter terrestrial carbon cycling by lowering the quantity and quality of carbon inputs to soils. However, the influence of elevated ozone levels on soil carbon formation and decomposition are unknown. Here we examine the effects of elevated ozone levels on the formation rates of total and decay-resistant acid-insoluble soil carbon under conditions of elevated carbon dioxide levels in experimental aspen (Populus tremuloides) stands and mixed aspen-birch (Betula papyrifera) stands. With ambient concentrations of ozone and carbon dioxide both raised by 50 per cent, we find that the formation rates of total and acid-insoluble soil carbon are reduced by 50 per cent relative to the amounts entering the soil when the forests were exposed to increased carbon dioxide alone. Our results suggest that, in a world with elevated atmospheric carbon dioxide concentrations, global-scale reductions in plant productivity due to elevated ozone levels will also lower soil carbon formation rates significantly.     

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.     

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.     

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.