The millennial atmospheric lifetime of anthropogenic CO<Subscript>2</Subscript>

The notion is pervasive in the climate science community and in the public at large that the climate impacts of fossil fuel CO₂ release will only persist for a few centuries. This conclusion has no basis in theory or models of the atmosphere/ocean carbon cycle, which we review here. The largest fraction of the CO₂ recovery will take place on time scales of centuries, as CO₂ invades the ocean, but a significant fraction of the fossil fuel CO₂, ranging in published models in the literature from 20–60%, remains airborne for a thousand years or longer. Ultimate recovery takes place on time scales of hundreds of thousands of years, a geologic longevity typically associated in public perceptions with nuclear waste. The glacial/interglacial climate cycles demonstrate that ice sheets and sea level respond dramatically to millennial-timescale changes in climate forcing. There are also potential positive feedbacks in the carbon cycle, including methane hydrates in the ocean, and peat frozen in permafrost, that are most sensitive to the long tail of the fossil fuel CO₂ in the atmosphere.     

Transient simulation of the last glacial inception. Part II: sensitivity and feedback analysis

The sensitivity of the last glacial-inception (around 115 kyr BP, 115,000 years before present) to different feedback mechanisms has been analysed by using the Earth system model of intermediate complexity CLIMBER-2. CLIMBER-2 includes dynamic modules of the atmosphere, ocean, terrestrial biosphere and inland ice, the last of which was added recently by utilising the three-dimensonal polythermal ice-sheet model SICOPOLIS. We performed a set of transient experiments starting at the middle of the Eemiam interglacial and ran the model for 26,000 years with time-dependent orbital forcing and observed changes in atmospheric CO₂ concentration (CO₂ forcing). The role of vegetation and ocean feedback, CO₂ forcing, mineral dust, thermohaline circulation and orbital insolation were closely investigated. In our model, glacial inception, as a bifurcation in the climate system, appears in nearly all sensitivity runs including a run with constant atmospheric CO₂ concentration of 280 ppmv, a typical interglacial value, and simulations with prescribed present-day sea-surface temperatures or vegetation cover—although the rate of the growth of ice-sheets growth is smaller than in the case of the fully interactive model. Only if we run the fully interactive model with constant present-day insolation and apply present-day CO₂ forcing does no glacial inception appear at all. This implies that, within our model, the orbital forcing alone is sufficient to trigger the interglacial–glacial transition, while vegetation, ocean and atmospheric CO₂ concentration only provide additional, although important, positive feedbacks. In addition, we found that possible reorganisations of the thermohaline circulation influence the distribution of inland ice.     

Impact of anthropogenic CO 2 on the next glacial cycle

The model of Paillard and Parrenin (Earth Planet Sci Lett 227(3–4):263–271, 2004) has been recently optimized for the last eight glacial cycles, leading to two different relaxation models with model-data correlations between 0.8 and 0.9 (García-Olivares and Herrero (Clim Dyn 1–25, 2012b)). These two models are here used to predict the effect of an anthropogenic CO ₂ pulse on the evolution of atmospheric CO ₂, global ice volume and Antarctic ice cover during the next 300 kyr. The initial atmospheric CO ₂ condition is obtained after a critical data analysis that sets 1300 Gt as the most realistic carbon Ultimate Recoverable Resources (URR), with the help of a global compartmental model to determine the carbon transfer function to the atmosphere. The next 20 kyr will have an abnormally high greenhouse effect which, according to the CO ₂ values, will lengthen the present interglacial by some 25 to 33 kyr. This is because the perturbation of the current interglacial will lead to a delay in the future advance of the ice sheet on the Antarctic shelf, causing that the relative maximum of boreal insolation found 65 kyr after present (AP) will not affect the developing glaciation. Instead, it will be the following insolation peak, about 110 kyr AP, which will find an appropriate climatic state to trigger the next deglaciation.     

Role of the ocean in controlling atmospheric CO2 concentration in the course of global glaciations

Responses of ocean circulation and ocean carbon cycle in the course of a global glaciation from the present Earth conditions are investigated by using a coupled climate-biogeochemical model. We investigate steady states of the climate system under colder conditions induced by a reduction of solar constant from the present condition. A globally ice-covered solution is obtained under the solar constant of 92.2% of the present value. We found that because almost all of sea water reaches the frozen point, the ocean stratification is maintained not by temperature but by salinity just before the global glaciation (at the solar constant of 92.3%). It is demonstrated that the ocean circulation is driven not by the surface cooling but by the surface freshwater forcing associated with formation and melting of sea ice. As a result, the deep ocean is ventilated exclusively by deep water formation in southern high latitudes where sea ice production takes place much more massively than northern high latitudes. We also found that atmospheric CO₂ concentration decreases through the ocean carbon cycle. This reduction is explained primarily by an increase of solubility of CO₂ due to a decrease of sea surface temperature, whereas the export production weakens by 30% just before the global glaciation. In order to investigate the conditions for the atmospheric CO₂ reduction to cause global glaciations, we also conduct a series of simulations in which the total amount of carbon in the atmosphere?Cocean system is reduced from the present condition. Under the present solar constant, the results show that the global glaciation takes place when the total carbon decreases to be 70% of the present-day value. Just before the glaciation, weathering rate becomes very small (almost 10% of the present value) and the organic carbon burial declines due to weakened biological productivity. Therefore, outgoing carbon flux from the atmosphere?Cocean system significantly decreases. This suggests the atmosphere?Cocean system has strong negative feedback loops against decline of the total carbon content. The results obtained here imply that some processes outside the atmosphere?Cocean feedback loops may be required to cause global glaciations.     

Glacial-Interglacial Atmospheric CO2 Change—The Glacial Burial Hypothesis

Organic carbon buried under the great ice sheets of the Northern Hemisphere is suggested to be the missing link in the atmospheric CO2 change over the glacial-interglacial cycles. At glaciation, the advancement of continental ice sheets buries vegetation and soil carbon accumulated during warmer pe-riods. At deglaciation, this burial carbon is released back into the atmosphere. In a simulation over two glacial-interglacial cycles using a synchronously coupled atmosphere-land-ocean carbon model forced by reconstructed climate change, it is found that there is a 547-Gt terrestrial carbon release from glacial maximum to interglacial, resulting in a 60-Gt (about 30-ppmv) increase in the atmospheric CO2, with the remainder absorbed by the ocean in a scenario in which ocean acts as a passive buffer. This is in contrast to previous estimates of a land uptake at deglaciation. This carbon source originates from glacial burial,continental shelf, and other land areas in response to changes in ice cover, sea level, and climate. The input of light isotope enriched terrestrial carbon causes atmospheric δ^13C to drop by about 0.3‰ at deglaciation,followed by a rapid rise towards a high interglacial value in response to oceanic warming and regrowth on land. Together with other ocean based mechanisms such as change in ocean temperature, the glacial burial hypothesis may offer a full explanation of the observed 80-100-ppmv atmospheric CO2 change.     

Ice and climate modeling: An editorial essay

The growth and decay of ice sheets are driven by forces affecting the seasonal cycles of snowfall and snowmelt. The external forces are likely to be variations in the earth's orbit which cause differences in the solar radiation received. Radiational control of snowmelt is modulated by the seasonal cycles of snow albedo and cloud cover. The effects of orbital changes can be magnified by feedbacks involving atmospheric CO₂ content, ocean temperatures and desert areas. Climate modeling of the causes of the Pleistocene ice ages involves modeling the interactions of all components of the climate system; snow, sea ice, glacier ice, the ocean, the atmosphere, and the solid earth. Such modeling is also necessary for interpreting oxygen isotope records from ice and ocean as paleoclimatic evidence.     

World ocean temperature lag time: an analysis based on glaciation data for the last two million years

Summary ?It is postulated that before the influence of glaciation, it was the amount of cloud cover and the thermal inertia of the ocean that controlled the Earth’s temperature. The control system went into oscillation 37 myr BP when Antarctica started moving into its present position, the temperature of the ocean and that of the rest of the environment opposing each other in antisymmetric mode. Support for this theory is provided by the observation of periods of enhanced glaciation at regular intervals. The enhancement, being attributed to harmonics with the Earth’s 22,000 yr-precession and 41,000 yr-nutation cycles, allows the calculation of 23,500 yr for the period of the ocean/atmosphere-temperature cycle. The corresponding lag time between atmosphere and ocean is 11,750 yr. Received February 17, 2002; revised March 19, 2002; accepted April 9, 2002     

CO2 threshold for millennial-scale oscillations in the climate system: implications for global warming scenarios

We present several equilibrium runs under varying atmospheric CO₂ concentrations using the University of Victoria Earth System Climate Model (UVic ESCM). The model shows two very different responses: for CO₂ concentrations of 400 ppm or lower, the system evolves into an equilibrium state. For CO₂ concentrations of 440 ppm or higher, the system starts oscillating between a state with vigorous deep water formation in the Southern Ocean and a state with no deep water formation in the Southern Ocean. The flushing events result in a rapid increase in atmospheric temperatures, degassing of CO₂ and therefore an increase in atmospheric CO₂ concentrations, and a reduction of sea ice cover in the Southern Ocean. They also cool the deep ocean worldwide. After the flush, the deep ocean warms slowly again and CO₂ is taken up by the ocean until the stratification becomes unstable again at high latitudes thousands of years later. The existence of a threshold in CO₂ concentration which places the UVic ESCM in either an oscillating or non-oscillating state makes our results intriguing. If the UVic ESCM captures a mechanism that is present and important in the real climate system, the consequences would comprise a rapid increase in atmospheric carbon dioxide concentrations of several tens of ppm, an increase in global surface temperature of the order of 1–2°C, local temperature changes of the order of 6°C and a profound change in ocean stratification, deep water temperature and sea ice cover.     

Transient simulation of the last glacial inception. Part I: glacial inception as a bifurcation in the climate system

We study the mechanisms of glacial inception by using the Earth system model of intermediate complexity, CLIMBER-2, which encompasses dynamic modules of the atmosphere, ocean, biosphere and ice sheets. Ice-sheet dynamics are described by the three-dimensional polythermal ice-sheet model SICOPOLIS. We have performed transient experiments starting at the Eemiam interglacial, at 126 ky BP (126,000 years before present). The model runs for 26 kyr with time-dependent orbital and CO₂ forcings. The model simulates a rapid expansion of the area covered by inland ice in the Northern Hemisphere, predominantly over Northern America, starting at about 117 kyr BP. During the next 7 kyr, the ice volume grows gradually in the model at a rate which corresponds to a change in sea level of 10 m per millennium. We have shown that the simulated glacial inception represents a bifurcation transition in the climate system from an interglacial to a glacial state caused by the strong snow-albedo feedback. This transition occurs when summer insolation at high latitudes of the Northern Hemisphere drops below a threshold value, which is only slightly lower than modern summer insolation. By performing long-term equilibrium runs, we find that for the present-day orbital parameters at least two different equilibrium states of the climate system exist—the glacial and the interglacial; however, for the low summer insolation corresponding to 115 kyr BP, we find only one, glacial, equilibrium state, while for the high summer insolation corresponding to 126 kyr BP only an interglacial state exists in the model.

---

Sea‐ice dynamics and CO2 sensitivity in a global climate model

Abstract

Present‐day results and CO₂ sensitivity are described for two versions of a global climate model (genesis) with and without sea‐ice dynamics. Sea‐ice dynamics is modelled using the cavitating‐fluid method of Flato and Hibler (1990, 1992). The atmospheric general circulation model originated from the NCAR Community Climate Model version 1, but is heavily modified to include new treatments of clouds, penetrative convection, planetary boundary‐layer mixing, solar radiation, the diurnal cycle and the semi‐Lagrangian transport of water vapour. The surface models include an explicit model of vegetation (similar to BATS and SiB), multilayer models of soil, snow and sea ice, and a slab ocean mixed layer.

When sea‐ice dynamics is turned off, the CO₂‐induced warming increases drastically around ～60–80°S in winter and spring. This is due to the much greater (and unrealistic) compactness of the Antarctic ice cover without dynamics, which is reduced considerably when CO₂ is doubled and exposes more open ocean to the atmosphere. With dynamics, the winter ice is already quite dispersed for 1 × CO₂ so that its compactness does not decrease as much when CO₂ is doubled.     

The Pangean ice age: studies with a coupled climate-ice sheet model

 Application of an ice sheet model developed for the Pleistocene to the extensive Carboniferous glaciation on Gondwana yields an ice sheet which has several features consistent with observations. While complete deglaciation is not achieved without CO₂ changes, the Milankovich-induced fluctuations in ice sheet volume are comparable to Pleistocene glacial/ interglacial signals. This result is shown to hold for a large fraction of physically reasonable parameter space. The model also exhibits multiple equilibria and sharp bifurcations, as infinitesimal changes in the solar constant or precipitation can lead to a qualitatively different climate. The success of the model in predicting ice location in an environment quite different from the Pleistocene provides additional support for the robustness of the basic model physics and suggests that the model can be applied with some confidence to other pre-Pleistocene glaciations. Received: 30 June 1998 / Accepted: 5 January 1999     

A 30,000 yr precession-related cycle affecting climate

Summary ?Evidence is presented for a previously unknown climate cycle of 30,000 yr period. The cycle is deemed to be related to the gyroscopic precession, or wobble of the Earth axis. Since it inhibits glaciation, the 30,000 yr cycle is called the summer cycle while its counterpart, the 22,000 yr “precession” cycle, is the winter cycle. Because of the aspect presented to the Sun, summer is effectively longer than winter. This is used to explain the difference between summer and winter cycle periods and that of the wobble. Some of the problems encountered in interpreting oxygen ratio glaciation data are resolved by knowledge of the existence of the summer cycle and a means is devised for determining the relative volume of glacial ice. An argument is made that essentially eliminates the 100,000 yr-orbit cycle, by itself or in combination with other attitude cycles, as a possible cause of the glacial/interglacial cycles. Received August 17, 2002; accepted September 3, 2002     

Climate model simulation of anthropogenic influence on greenhouse-induced climate change (early agriculture to modern): the role of ocean feedbacks

We present further steps in our analysis of the early anthropogenic hypothesis (Ruddiman, Clim Change 61:261–293, 2003) that increased levels of greenhouse gases in the current interglacial, compared to lower levels in previous interglacials, were initiated by early agricultural activities, and that these increases caused a warming of climate long before the industrial era (～1750). These steps include updating observations of greenhouse gas and climate trends from earlier interglacials, reviewing recent estimates of greenhouse gas emissions from early agriculture, and describing a simulation by a climate model with a dynamic ocean forced by the low levels of greenhouse gases typical of previous interglacials in order to gauge the magnitude of the climate change for an inferred (natural) low greenhouse gas level relative to a high present day level. We conduct two time slice (equilibrium) simulations using present day orbital forcing and two levels of greenhouse gas forcing: the estimated low (natural) levels of previous interglacials, and the high levels of the present (control). By comparing the former to the latter, we estimate how much colder the climate would be without the combined greenhouse gas forcing of the early agriculture era (inferred from differences between this interglacial and previous interglacials) and the industrial era (the period since ～1750). With the low greenhouse gas levels, the global average surface temperature is 2.7 K lower than present day—ranging from ～2 K lower in the tropics to 4–8 K lower in polar regions. These changes are large, and larger than those reported in a pre-industrial (～1750) simulation with this model, because the imposed low greenhouse gas levels (CH₄ = 450 ppb, CO₂ = 240 ppm) are lower than both pre-industrial (CH₄ = 760 ppb, CO₂ = 280 ppm) and modern control (CH₄ = 1,714 ppb, CO₂ = 355 ppm) values. The area of year-round snowcover is larger, as found in our previous simulations and some other modeling studies, indicating that a state of incipient glaciation would exist given the current configuration of earth’s orbit (reduced insolation in northern hemisphere summer) and the imposed low levels of greenhouse gases. We include comparisons of these snowcover maps with known locations of earlier glacial inception and with locations of twentieth century glaciers and ice caps. In two earlier studies, we used climate models consisting of atmosphere, land surface, and a shallow mixed-layer ocean (Ruddiman et al., Quat Sci Rev 25:1–10, 2005; Vavrus et al., Quat Sci Rev 27:1410–1425, 2008). Here, we replaced the mixed-layer ocean with a complete dynamic ocean. While the simulated climate of the atmosphere and the surface with this improved model configuration is similar to our earlier results (Vavrus et al., Quat Sci Rev 27:1410–1425, 2008), the added information from the full dynamical ocean is of particular interest. The global and vertically-averaged ocean temperature is 1.25 K lower, the area of sea ice is larger, and there is less upwelling in the Southern Ocean. From these results, we infer that natural ocean feedbacks could have amplified the greenhouse gas changes initiated by early agriculture and possibly account for an additional increment of CO₂ increase beyond that attributed directly to early agricultural, as proposed by Ruddiman (Rev Geophys 45:RG4001, 2007). However, a full test of the early anthropogenic hypothesis will require additional observations and simulations with models that include ocean and land carbon cycles and other refinements elaborated herein.     

Snow cover validation and sensitivity to CO2 in the UVic ESCM

Abstract

The University of Victoria's (UVic) Earth System Climate model is used to conduct equilibrium atmospheric CO₂ sensitivity experiments over the range 200–1600 ppm in order to explore changes in northern hemisphere snow cover and feedbacks on terrestrial surface air temperature (SAT). Simulations of warmer climates predict a retreat of snow cover over northern continents, in a northeasterly direction. The decline in northern hemisphere global snow mass is estimated to reach 33% at 600 ppm and 54% at 1200 ppm. In the most northerly regions, annual mean snow depth increases for simulations with CO₂ levels higher than present day. The shift in the latitude of maximum snowfall is estimated to be inversely proportional to the CO₂ concentration. The northern hemisphere net shortwave radiation changes are found to be greater over land than over the ocean, suggesting a stronger albedo feedback from changes in terrestrial snow cover than from changes in sea ice. Results also reveal high sensitivity of the snow mass balance under low CO₂ conditions. The amplification feedback (defined as the zonal SAT anomaly caused by doubling CO₂ divided by the equatorial anomaly) is greatest for scenarios with less than 300 ppm, reaching 1.9 at the pole for 250 ppm. The stronger feedback is attributed to the significant albedo changes over land areas. The simulation with 200 ppm triggers continuous accumulation of snow ('glaciation') in regions which, according to paleo‐reconstructions, were covered by ice during the last glacial cycle (the Canadian Arctic, Scandinavia, and the Taymir Peninsula).     

The influence of continental ice,atmospheric CO2, and land albedo on the climate of the last glacial maximum

The contributions of expanded continental ice, reduced atmospheric CO₂, and changes in land albedo to the maintenance of the climate of the last glacial maximum (LGM) are examined. A series of experiments is performed using an atmosphere-mixed layer ocean model in which these changes in boundary conditions are incorporated either singly or in combination. The model used has been shown to produce a reasonably realistic simulation of the reduced temperature of the LGM (Manabe and Broccoli 1985b). By comparing the results from pairs of experiments, the effects of each of these environmental changes can be determined.Expanded continental ice and reduced atmospheric CO₂ are found to have a substantial impact on global mean temperature. The ice sheet effect is confined almost exclusively to the Northern Hemisphere, while lowered CO₂ cools both hemispheres. Changes in land albedo over ice-free areas have only a minor thermal effect on a global basis. The reduction of CO₂ content in the atmosphere is the primary contributor to the cooling of the Southern Hemisphere. The model sensitivity to both the ice sheet and CO₂ effects is characterized by a high latitude amplification and a late autumn and early winter maximum.Substantial changes in Northern Hemisphere tropospheric circulation are found in response to LGM boundary conditions during winter. An amplified flow pattern and enhanced westerlies occur in the vicinity of the North American and Eurasian ice sheets. These alterations of the tropospheric circulation are primarily the result of the ice sheet effect, with reduced CO₂ contributing only a slight amplification of the ice sheet-induced pattern.     

Analysis of measurement data on carbon dioxide concentration in the near-ice surface atmosphere at North Pole-35 drifting ice station (2007–2008)

Significant CO₂ concentration disturbances with an amplitude of 150 ppm are observed in October–November in the region of North Pole-35 drifting ice station over the continental slope of the Arctic Ocean. A local maximum of the CO₂ concentration disturbances is observed in April; it is connected with a change in the regime of the sea ice deformation at that time of the year in the western Arctic basin. Changes in the CO₂ concentration have oscillations at the frequency of oceanic tides. Several episodes of the initial phase of CO₂ release into the near-ice surface atmosphere in October 2007 are recorded. The CO₂ amount released into the atmosphere as a result of ice formation in the Arctic Ocean is a significant quantity in the total balance of the CO₂ release on the global scale and accounts for approximately 30% of the anthropogenic release of CO₂ per year.     

Sensitivity of the last glacial inception to initial and surface conditions

We investigate the sensitivity of simulations of the last glacial inception (LGI) with respect to initial (size of the Greenland ice sheet) and surface (state of ocean/vegetation) conditions and two different CO₂ reconstructions. Utilizing the CLIMBER-2 Earth system model, we obtain the following results: (a) ice-sheet expansion in North America at the end of the Eemian can be reduced or even completely suppressed when pre-industrial or Eemian ocean/vegetation is prescribed. (b) A warmer surrounding ocean and, in particular, a large Laurentide ice sheet reduce the size of the Greenland ice sheet before and during the LGI. (c) A changing ocean contributes much stronger to the expansion of the Laurentide ice sheet when we apply the CO₂ reconstruction according to Barnola et al. (Nature 329:408–414, 1987) instead of Petit et al. (Nature 399:429–436, 1999). (d) In the fully coupled model, the CO₂ reconstruction used has only a small impact on the simulated ice sheets but it does impact the course of the climatic variables. (e) For the Greenland ice sheet, two equilibrium states exist under the insolation and CO₂ forcing at 128,000 years before present (128 kyear BP); the one with an ice sheet reduced by about one quarter as compared to its simulated pre-industrial size and the other with nearly no inland ice in Greenland. (f) Even the extreme assumption of no ice sheet in Greenland at the beginning of our transient simulations does not alter the simulated expansion of northern hemispheric ice sheets at the LGI.     

Implications of CO2 global warming on great lakes ice cover

Statistical ice cover models were used to project daily mean basin ice cover and annual ice cover duration for Lakes Superior and Erie. Models were applied to a 1951–80 base period and to three 30-year steady double carbon dioxide (2 × CO₂) scenarios produced by the Geophysical Fluid Dynamics Laboratory (GFDL), the Goddard Institute of Space Studies (GISS), and the Oregon State University (OSU) general circulation models. Ice cover estimates were made for the West, Central, and East Basins of Lake Erie and for the West, East, and Whitefish Bay Basins of Lake Superior. Average ice cover duration for the 1951– 80 base period ranged from 13 to 16 weeks for individual lake basins. Reductions in average ice cover duration under the three 2 × CO₂ scenarios for individual lake basins ranged from 5 to 12 weeks for the OSU scenario, 8 to 13 weeks for the GISS scenario, and 11 to 13 weeks for GFDL scenario. Winters without ice formation become common for Lake Superior under the GFDL scenario and under all three 2 × CO₂ scenarios for the Central and East Basins of Lake Erie. During an average 2 × CO₂ winter, ice cover would be limited to the shallow areas of Lakes Erie and Superior. Because of uncertainties in the ice cover models, the results given here represent only a first approximation and are likely to represent an upper limit of the extent and duration of ice cover under the climate change projected by the three 2 × CO₂scenarios. Notwithstanding these limitations, ice cover projected by the 2 × CO₂ scenarios provides a preliminary assessment of the potential sensitivity of the Great Lakes ice cover to global warming. Potential environmental and socioeconomic impacts of a 2 × CO₂ warming include year-round navigation, change in abundance of some fish species in the Great Lakes, discontinuation or reduction of winter recreational activities, and an increase in winter lake evaporation.     

CO2 feedbacks and the 100 K year cycle in climate

Summary In an earlier paper (Lindzen, 1986), it was shown that allowing CO₂ to vary with snow/sea ice position could lead to a greatly enhanced response in glaciation to 100 K year orbital forcing—even when 20 K year forcing was much stronger. In that model, snow/sea ice position (SSIP) and glaciation were different: the former was the forcing for the latter. However, SSIP and glaciation were not decorrelated. Observations (Berner et al., 1979; Lorius et al., 1985; Neftel et al., 1982) suggest that CO₂ may be independently related to both SSIP and glaciation. In the present paper, we allow (in a highly simplified manner) such independent dependence, and show how it alters the earlier results. Briefly, the dependence of CO₂ on glaciation can contribute to and even cause a highly enhanced response to the 100 K year component of the forcing. However, the CO₂ dependence on SSIP is, on the whole, more effective in this regard. Thus, we expect time series of CO₂ to show variation on the faster time scales than does glaciation.With 5 Figures