Declining Temporal Effectiveness of Carbon Sequestration: Implications for Compliance with the United National Framework Convention on Climate Change

Carbon sequestration is increasingly being promoted as a potential response to the risks of unrestrained emissions of CO₂, either in place of or as a complement to reductions in the use of fossil fuels. However, the potential role of carbon sequestration as an (at-least partial) substitute for reductions in fossil fuel use can be properly evaluated only in the context of a long-term acceptable limit (or range of limits) to the increase in atmospheric CO₂ concentration, taking into account the response of the entire carbon cycle to artificial sequestration. Under highly stringent emission-reduction scenarios for non-CO₂ greenhouse gases, 450 ppmv CO₂ is the equivalent, in terms of radiative forcing of climate,to a doubling of the pre-industrial concentration of CO₂. It is argued in this paper that compliance with the United Nations Framework Convention on Climate Change (henceforth, the UNFCCC) implies that atmospheric CO₂ concentration should be limited, or quickly returned to, a concentration somewhere below 450 ppmv. A quasi-one-dimensional coupled climate-carbon cycle model is used to assess the response of the carbon cycle to idealized carbon sequestration scenarios. The impact on atmospheric CO₂ concentration of sequestering a given amount of CO₂ that would otherwise be emitted to the atmosphere, either in deep geological formations or in the deep ocean, rapidly decreases over time. This occurs as a result of a reduction in the rate of absorption of atmospheric CO₂ by the natural carbon sinks (the terrestrial biosphere and oceans) in response to the slower buildup of atmospheric CO₂ resulting from carbon sequestration. For 100 years of continuous carbon sequestration, the sequestration fraction (defined as the reduction in atmospheric CO₂ divided by the cumulative sequestration) decreases to 14% 1000 years after the beginning of sequestration in geological formations with no leakage, and to 6% 1000 years after the beginning of sequestration in the deep oceans. The difference (8% of cumulative sequestration) is due to an eflux from the ocean to the atmosphere of some of the carbon injected into the deep ocean.The coupled climate-carbon cycle model is also used to assess the amount of sequestration needed to limit or return the atmospheric CO₂ concentration to 350–400 ppmv after phasing out all use of fossil fuels by no later than 2100. Under such circumstances, sequestration of 1–2 Gt C/yr by the latter part of this century could limit the peak CO₂ concentration to 420–460 ppmv, depending on how rapidly use of fossilfuels is terminated and the strength of positive climate-carbon cycle feedbacks. To draw down the atmospheric CO₂ concentration requires creating negative emissions through sequestration of CO₂ released as a byproduct of the production of gaseous fuels from biomass primary energy. Even if fossil fuel emissions fall to zero by 2100, it will be difficult to create a large enough negative emission using biomass energy to return atmospheric CO₂ to 350 ppmv within 100 years of its peak. However, building up soil carbon could help in returning CO₂ to 350 ppmv within 100 years of its peak. In any case, a 100-year period of climate corresponding to the equivalent of a doubled-CO₂ concentration would occur before temperatures decreased. Nevertheless, returning the atmospheric CO₂concentration to 350 ppmv would reduce longterm sea level rise due to thermal expansion and might be sufficient to prevent the irreversible total melting of the Greenland ice sheet, collapse of the West Antarctic ice sheet, and abrupt changes in ocean circulation that might otherwise occur given a prolonged doubled-CO₂ climate. Recovery of coral reef ecosystems, if not already driven to extinction, could begin.     

Commentary on “The Anthropogenic Greenhouse Era Began Thousands of Years Ago”

Bill Ruddiman (Climatic Change, 61, 261–293, 2003) recently suggested that early civilisations could have saved us from an ice age because land management over substantial areas caused an increase in atmospheric CO₂ concentration. Ruddiman suggests a decreasing “natural course” of the Holocene greenhouse gases concentrations and sea-level by referring to analogous situations in the past, namely the last three interglacials. An examination of marine isotopic stage 11 would perhaps make Ruddiman’s argument even more thought-challenging. Yet, the hypothesis of a natural lowering of CO₂ during the Holocene contradicts recent numerical simulations of the Earth carbon cycle during this period. We think that the only way to resolve this conflict is to properly assimilate the palæoclimate information in numerical climate models. As a general rule, models are insufficiently tested with respect to the wide range of climate situations that succeeded during the Pleistocene. In this comment, we present three definitions of palæoclimate information assimilation with relevant examples. We also present original results with the Louvain-la-Neuve climate-ice sheet model suggesting that if, indeed, the Holocene atmospheric CO₂ increase is anthropogenic, a late Holocene glacial inception is plausible, but not certain, depending on the exact time evolution of the atmospheric CO₂ concentration during this period.     

The danger of overvaluing methane’s influence on future climate change

Minimizing the future impacts of climate change requires reducing the greenhouse gas (GHG) load in the atmosphere. Anthropogenic emissions include many types of GHG’s as well as particulates such as black carbon and sulfate aerosols, each of which has a different effect on the atmosphere, and a different atmospheric lifetime. Several recent studies have advocated for the importance of short timescales when comparing the climate impact of different climate pollutants, placing a high relative value on short-lived pollutants, such as methane (CH₄) and black carbon (BC) versus carbon dioxide (CO₂). These studies have generated confusion over how to value changes in temperature that occur over short versus long timescales. We show the temperature changes that result from exchanging CO₂ for CH₄ using a variety of commonly suggested metrics to illustrate the trade-offs involved in potential carbon trading mechanisms that place a high value on CH₄ emissions. Reducing CH₄ emissions today would lead to a climate cooling of approximately ~0.5 °C, but this value will not change greatly if we delay reducing CH₄ emissions by years or decades. This is not true for CO₂, for which the climate is influenced by cumulative emissions. Any delay in reducing CO₂ emissions is likely to lead to higher cumulative emissions, and more warming. The exact warming resulting from this delay depends on the trajectory of future CO₂ emissions but using one business-as usual-projection we estimate an increase of 3/4 °C for every 15-year delay in CO₂ mitigation. Overvaluing the influence of CH₄ emissions on climate could easily result in our “locking” the earth into a warmer temperature trajectory, one that is temporarily masked by the short-term cooling effects of the CH₄ reductions, but then persists for many generations.     

The effects of carbon cycle model error in calculating future atmospheric carbon dioxide levels

Empirical investigations have indicated that projections of future atmospheric carbon dioxide concentrations of a quality quite adequate for practical questions regarding the environmental threat of anthropogenic carbon dioxide emissions and its relationship to energy use policy could be made with the simple assumption that a constant fraction of these emissions would be retained by the atmosphere. By analysis of the structural behavior of equations describing the transfer of carbon and carbon dioxide between their several reservoirs we have been able to demonstrate that this characteristic can be explained to result from approximately linear behavior and exponentially growing carbon dioxide release rates, combined with fitting of carbon cycle model parameters to the last twenty years of observed atmospheric carbon dioxide growth. These conclusions are independent of the details of carbon cycle model structure for projections up to 100 years into the future as long as the growth in atmospheric carbon dioxide release rates is sufficiently high, of the order of 1.5% per annum or more, as referenced to p re-industrial (steady state) conditions. At low rates of growth, when the longer response times of the carbon cycling system become important, for most energy use projections the resultant CO₂ induced climate changes are small and the uncertainties in predicted atmospheric carbon dioxide level are thus not important. A possible exception to this condition occurs for scenarios of future fossil fuel use rates designed to avoid atmospheric CO₂ levels exceeding a chosen threshold. In this instance details of carbon cycle model structure could significantly affect conclusions that might be drawn concerning future energy use policies; however, it is possible that such a result stems from inappropriate specification of a criterion for an environmental threat, rather than from inherent inadequacy of current carbon cycle models. Recent carbon cycle model developments postulate transfer processes of carbon into the deep ocean, large carbon storage reservoir at rates much higher than in the models we have analysed. If the existence of such mechanisms is confirmed, and they are found to be sufficiently rapid and large, some of our conclusions regarding the use of the constant fractional retention assumption may have to be modified. Currently at the Gas Research Institute, 8600 West Bryn, Mawr Ave., Chicago, IL 60631, U.S.A.     

Accounting for the missing carbon-sink with the CO2-fertilization effect

A terrestrial-biosphere carbon-sink has been included in global carbon-cycle models in order to reproduce past atmospheric CO₂, ¹³C and ¹⁴C concentrations. The sink is of large enough magnitude that its effect on projections of future CO₂ levels should not be ignored. However, the cause and mechanism of this sink are not well understood, contributing to uncertainty of projections. The estimated magnitude of the biospheric sink is examined with the aid of a global carbon-cycle model. For CO₂ emissions scenarios, model estimates are made of the resulting atmospheric CO₂ concentration. Next, the response of this model to CO₂-emission impulses is broken down to give the fractions of the impulse which reside in the atmosphere, oceans, and terrestrial biosphere - all as a perturbation to background atmospheric CO₂ concentration time-profiles that correspond to different emission scenarios. For a biospheric sink driven by the CO₂-fertilization effect, we find that the biospheric fraction reaches a maximum of roughly 30% about 50 years after the impulse, which is of the same size as the oceanic fraction at that time. The dependence of these results on emission scenario and the year of the impulse are reported.     

Interactions between perturbations to different Earth system components simulated by a fully-coupled climate model

We examine the global mean surface temperature and carbon cycle responses to the A1B emissions scenario for a new 57 member perturbed-parameter ensemble of simulations generated using the fully coupled atmosphere-ocean-carbon cycle climate model HadCM3C. The model variants feature simultaneous perturbation to parameters that control atmosphere, ocean, land carbon cycle and sulphur cycle processes in this Earth system model, and is the first experiment of its kind. The experimental design, based on four earlier ensembles with parameters varied within each individual Earth system component, allows the effects of interactions between uncertainties in the different components to be explored. A large spread in response is obtained, with atmospheric CO₂ at the end of the twenty-first century ranging from 615 to 1,100 ppm. On average though, the mean effect of the parameter perturbations is to significantly reduce the amount of atmospheric CO₂ compared to that seen in the standard HadCM3C model. Global temperature change for 2090–2099 relative to the pre-industrial period ranges from 2.2 to 7.5 °C, with large temperature responses occurring when atmospheric model versions with high climate sensitivities are combined with carbon cycle components that emit large amounts of CO₂ to the atmosphere under warming. A simple climate model, tuned to reproduce the responses of the separate Earth system component ensembles, is used to demonstrate that interactions between uncertainties in the different components play a significant role in determining the spread of responses in global mean surface temperature. This ensemble explores a wide range of interactions and response, and therefore provides a useful resource for the provision of regional climate projections and associated uncertainties.     

Simulated changes in vegetation distribution, land carbon storage, and atmospheric CO2 in response to a collapse of the North Atlantic thermohaline circulation

It is investigated how abrupt changes in the North Atlantic (NA) thermohaline circulation (THC) affect the terrestrial carbon cycle. The Lund–Potsdam–Jena Dynamic Global Vegetation Model is forced with climate perturbations from glacial freshwater experiments with the ECBILT-CLIO ocean–atmosphere–sea ice model. A reorganisation of the marine carbon cycle is not addressed. Modelled NA THC collapses and recovers after about a millennium in response to prescribed freshwater forcing. The initial cooling of several Kelvin over Eurasia causes a reduction of extant boreal and temperate forests and a decrease in carbon storage in high northern latitudes, whereas improved growing conditions and slower soil decomposition rates lead to enhanced storage in mid-latitudes. The magnitude and evolution of global terrestrial carbon storage in response to abrupt THC changes depends sensitively on the initial climate conditions. These were varied using results from time slice simulations with the Hadley Centre model HadSM3 for different periods over the past 21 kyr. Changes in terrestrial storage vary between −67 and +50 PgC for the range of experiments with different initial conditions. Simulated peak-to-peak differences in atmospheric CO₂ are 6 and 13 ppmv for glacial and late Holocene conditions. Simulated changes in δ¹³C are between 0.15 and 0.25‰. These simulated carbon storage anomalies during a NA THC collapse depend on their magnitude on the CO₂ fertilisation feedback mechanism. The CO₂ changes simulated for glacial conditions are compatible with available evidence from marine studies and the ice core CO₂ record. The latter shows multi-millennial CO₂ variations of up to 20 ppmv broadly in parallel with the Antarctic warm events A1 to A4 in the South and cooling in the North.     

Constraints on the temperature sensitivity of global soil respiration from the observed interannual variability in atmospheric CO2

Carbon cycle feedbacks have been shown to be very important in predicting climate change over the next century. The response of the terrestrial carbon cycle to climate change depends on the competition between increased respiration due to warmer temperatures and increased uptake due to elevated CO₂levels. Whether the terrestrial carbon cycle remains a sink for anthropogenic carbon, or switches to become a source, depends particularly on the response of soil respiration to temperature. Here we use observed global atmospheric CO₂concentration to constrain the behaviour of soil respiration in a coupled climate–carbon cycle GCM.     

Dependence of climate forcing and response on the altitude of black carbon aerosols

Black carbon aerosols absorb solar radiation and decrease planetary albedo, and thus can contribute to climate warming. In this paper, the dependence of equilibrium climate response on the altitude of black carbon is explored using an atmospheric general circulation model coupled to a mixed layer ocean model. The simulations model aerosol direct and semi-direct effects, but not indirect effects. Aerosol concentrations are prescribed and not interactive. It is shown that climate response of black carbon is highly dependent on the altitude of the aerosol. As the altitude of black carbon increases, surface temperatures decrease; black carbon near the surface causes surface warming, whereas black carbon near the tropopause and in the stratosphere causes surface cooling. This cooling occurs despite increasing planetary absorption of sunlight (i.e. decreasing planetary albedo). We find that the trend in surface air temperature response versus the altitude of black carbon is consistent with our calculations of radiative forcing after the troposphere, stratosphere, and land surface have undergone rapid adjustment, calculated as “regressed” radiative forcing. The variation in climate response from black carbon at different altitudes occurs largely from different fast climate responses; temperature dependent feedbacks are not statistically distinguishable. Impacts of black carbon at various altitudes on the hydrological cycle are also discussed; black carbon in the lowest atmospheric layer increases precipitation despite reductions in solar radiation reaching the surface, whereas black carbon at higher altitudes decreases precipitation.     

The role of ecosystem-atmosphere interactions in simulated Amazonian precipitation decrease and forest dieback under global climate warming

Summary A suite of simulations with the HadCM3LC coupled climate-carbon cycle model is used to examine the various forcings and feedbacks involved in the simulated precipitation decrease and forest dieback. Rising atmospheric CO₂ is found to contribute 20% to the precipitation reduction through the physiological forcing of stomatal closure, with 80% of the reduction being seen when stomatal closure was excluded and only radiative forcing by CO₂ was included. The forest dieback exerts two positive feedbacks on the precipitation reduction; a biogeophysical feedback through reduced forest cover suppressing local evaporative water recycling, and a biogeochemical feedback through the release of CO₂ contributing to an accelerated global warming. The precipitation reduction is enhanced by 20% by the biogeophysical feedback, and 5% by the carbon cycle feedback from the forest dieback. This analysis helps to explain why the Amazonian precipitation reduction simulated by HadCM3LC is more extreme than that simulated in other GCMs; in the fully-coupled, climate-carbon cycle simulation, approximately half of the precipitation reduction in Amazonia is attributable to a combination of physiological forcing and biogeophysical and global carbon cycle feedbacks, which are generally not included in other GCM simulations of future climate change. The analysis also demonstrates the potential contribution of regional-scale climate and ecosystem change to uncertainties in global CO₂ and climate change projections. Moreover, the importance of feedbacks suggests that a human-induced increase in forest vulnerability to climate change may have implications for regional and global scale climate sensitivity.     

An estimate of equilibrium sensitivity of global terrestrial carbon cycle using NCAR CCSM4

Increasing concentrations of atmospheric CO₂ influence climate, terrestrial biosphere productivity and ecosystem carbon storage through its radiative, physiological and fertilization effects. In this paper, we quantify these effects for a doubling of CO₂ using a low resolution configuration of the coupled model NCAR CCSM4. In contrast to previous coupled climate-carbon modeling studies, we focus on the near-equilibrium response of the terrestrial carbon cycle. For a doubling of CO₂, the radiative effect on the physical climate system causes global mean surface air temperature to increase by 2.14 K, whereas the physiological and fertilization on the land biosphere effects cause a warming of 0.22 K, suggesting that these later effects increase global warming by about 10 % as found in many recent studies. The CO₂-fertilization leads to total ecosystem carbon gain of 371 Gt-C (28 %) while the radiative effect causes a loss of 131 Gt-C (~10 %) indicating that climate warming damps the fertilization-induced carbon uptake over land. Our model-based estimate for the maximum potential terrestrial carbon uptake resulting from a doubling of atmospheric CO₂ concentration (285–570 ppm) is only 242 Gt-C. This highlights the limited storage capacity of the terrestrial carbon reservoir. We also find that the terrestrial carbon storage sensitivity to changes in CO₂ and temperature have been estimated to be lower in previous transient simulations because of lags in the climate-carbon system. Our model simulations indicate that the time scale of terrestrial carbon cycle response is greater than 500 years for CO₂-fertilization and about 200 years for temperature perturbations. We also find that dynamic changes in vegetation amplify the terrestrial carbon storage sensitivity relative to a static vegetation case: because of changes in tree cover, changes in total ecosystem carbon for CO₂-direct and climate effects are amplified by 88 and 72 %, respectively, in simulations with dynamic vegetation when compared to static vegetation simulations.     

A robust sequential CO2 emissions strategy based on optimal control of atmospheric CO2 concentrations

This paper formally introduces the concept of mitigation as a stochastic control problem. This is illustrated by applying a digital state variable feedback control approach known as Non-Minimum State Space (NMSS) control to the problem of specifying carbon emissions to control atmospheric CO₂ concentrations in the presence of uncertainty. It is shown that the control approach naturally lends itself to integrating both anticipatory and reflexive mitigation strategies within a single unified framework. The framework explicitly considers the closed-loop nature of climate mitigation, and employs a policy orientated optimisation procedure to specify the properties of this closed-loop system. The product of this exercise is a control law that is suitably conditioned to regulate atmospheric CO₂ concentrations through assimilating online information within a 25-year review cycle framework. It is shown that the optimal control law is also robust when faced with significant levels of uncertainty about the functioning of the global carbon cycle.     

Global transport and inter-reservoir exchange of carbon dioxide with particular reference to stable isotopic distributions

A two-dimensional model of global atmospheric transport is used to relate estimated air-to-surface exchanges of carbon dioxide (CO₂) to spatial and temporal variations of atmospheric CO₂ concentrations and isotopic composition. The atmospheric model coupled with models of the biosphere and mixed layer of the ocean describes the gross features of the global carbon cycle. In particular this paper considers the change in isotopic composition due to interreservoir exchanges and thus the potential application and measurement requirements of new isotopic observational programs.A comparison is made between the model-generated CO₂ concentration variation and those observed on secular, interannual and seasonal time scales and spatially through the depth of the troposphere and meridionally from pole-to-pole.The relationship between isotopic and concentration variation on a seasonal time-scale is discussed and it is shown how this can be used to quantitatively estimate relative contributions of biospheric and oceanic CO₂ exchange. Further, it is shown that the interhemispheric gradient of concentration and isotopic ratio results primarily from the redistribution of fossil fuel CO₂. Both isotopic and concentration data indicate that tropical deforestation contributes less than 2 Gt yr^-1 of carbon to the atmosphere.The study suggests that changes in the rate of change of the ratio of ¹³C to ¹²C in the atmosphere of less than 0.03 yr^-1 might be expected if net exchanges with the biosphere are the cause of interannual variations of CO₂ concentrations.     

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.     

A nonlinear impulse response model of the coupled carbon cycle-climate system (NICCS)

 Impulse-response-function (IRF) models are designed for applications requiring a large number of climate change simulations, such as multi-scenario climate impact studies or cost-benefit integrated-assessment studies. The models apply linear response theory to reproduce the characteristics of the climate response to external forcing computed with sophisticated state-of-the-art climate models like general circulation models of the physical ocean-atmosphere system and three-dimensional oceanic-plus-terrestrial carbon cycle models. Although highly computer efficient, IRF models are nonetheless capable of reproducing the full set of climate-change information generated by the complex models against which they are calibrated. While limited in principle to the linear response regime (less than about 3 ^∘C global-mean temperature change), the applicability of the IRF model presented has been extended into the nonlinear domain through explicit treatment of the climate system's dominant nonlinearities: CO₂ chemistry in ocean water, CO₂ fertilization of land biota, and sublinear radiative forcing. The resultant nonlinear impulse-response model of the coupled carbon cycle-climate system (NICCS) computes the temporal evolution of spatial patterns of climate change for four climate variables of particular relevance for climate impact studies: near-surface temperature, cloud cover, precipitation, and sea level. The space-time response characteristics of the model are derived from an EOF analysis of a transient 850-year greenhouse warming simulation with the Hamburg atmosphere-ocean general circulation model ECHAM3-LSG and a similar response experiment with the Hamburg carbon cycle model HAMOCC. The model is applied to two long-term CO₂ emission scenarios, demonstrating that the use of all currently estimated fossil fuel resources would carry the Earth's climate far beyond the range of climate change for which reliable quantitative predictions are possible today, and that even a freezing of emissions to present-day levels would cause a major global warming in the long term. Received: 28 January 2000 / Accepted: 9 March 2001     

不同陆地生态系统碳通量对GEOS-Chem模型模拟全球CO2浓度的影响

大气中CO₂含量的增加速率已经超过了自然界所能吸收的速度,并逐步影响到全球气候变暖。利用模型模拟分析已经成为一个重要的工具用以深入对碳循环的理解。本文使用2008～2010年的生物模型SiB3（Simple Biosphere version 3）与优化后的CT2016（Carbon Tracker 2016）陆地生态系统碳通量驱动GEOS-Chem大气化学传输模型模拟全球CO₂浓度。通过分析模拟CO₂浓度的空间分布与季节变化,加深对全球碳源汇分布特点的理解,探究陆地生态系统碳通量不确定性对模拟结果的影响,进而认识陆地生态系统碳通量反演精度提升的重要性。SiB3与优化后的CT2016陆地生态系统碳通量都具有明显的季节变化,但在欧洲地区碳源汇的表现相反,其全球总量与空间分布也存在极大的不确定性。模拟CO₂浓度结果表明:在人为活动较少地区,陆地生态系统碳通量对近地面CO₂浓度空间分布起主导作用,尤其在南半球和欧洲地区模拟浓度有明显差异,且两种模拟结果的季节差异依赖于陆地生态系统碳通量的季节变化。将模拟结果与9个观测站点资料进行对比,以期选用合适的陆地生态系统碳通量来提升GEOS-Chem模拟CO₂浓度的精度。实验结果表明:两种模拟结果均能较好的模拟CO₂浓度的季节变化及其峰谷值,但CT2016模拟的CO₂浓度在多数站点处更接近观测资料,模拟准确性更高。     

The role of closed ecological systems in carbon cycle modelling

Acquiring a mechanistic understanding of the role of biotic feedbacks for the links between atmospheric CO₂ concentrations and temperature is essential for trustworthy climate predictions. Currently, computer-based simulations are the only available tool to estimate the global impact of biotic feedbacks on future atmospheric CO₂ and temperatures. Here we propose an alternative and complementary approach by using materially closed, energetically open analogue/physical models of the carbon cycle. We argue that there is unexplored potential in using a materially closed approach to improve our understanding of the magnitude and direction of many biotic carbon feedbacks and that recent technological advances make this feasible. We also suggest how such systems could be designed and discuss the advantages and limitations of establishing physical models of the global carbon cycle.     

Carbon cycle–climate feedback sensitivity to parameter changes of a zero-dimensional terrestrial carbon cycle scheme in a climate model of intermediate complexity

Summary A series of sensitivity runs have been performed with a coupled climate–carbon cycle model. The climatic component consists of the climate model of intermediate complexity IAP RAS CM. The carbon cycle component is formulated as a simple zero-dimensional model. Its terrestrial part includes gross photosynthesis, and plant and soil respirations, depending on temperature via Q ₁₀-relationships (Lenton, 2000). Oceanic uptake of anthropogenic carbon is formulated is a bi-linear function of tendencies of atmospheric concentration of CO₂ and globally averaged annual mean sea surface temperature. The model is forced by the historical industrial and land use emissions of carbon dioxide for the second half of the 19th and the whole of the 20th centuries, and by the emission scenario SRES A2 for the 21st century. For the standard set of the governing parameters, the model realistically captures the main features of the Earth’s observed carbon cycle. A large number of simulations have been performed, perturbing the governing parameters of the terrestrial carbon cycle model. In addition, the climate part is perturbed, either by zeroing or artificially increasing the climate model sensitivity to the doubling of the atmospheric CO₂ concentration. Performing the above mentioned perturbations, it is possible to mimic most of the range found in the C4MIP simulations. In this way, a wide range of the climate–carbon cycle feedback strengths is obtained, differing even in the sign of the feedback. If the performed simulations are subjected to the constraints of a maximum allowed deviation of the simulated atmospheric CO₂ concentration (pCO_2(a)) from the observed values and correspondence between simulated and observed terrestrial uptakes, it is possible to narrow the corresponding uncertainty range. Among these constraints, considering pCO_2(a) and uptakes are both important. However, the terrestrial uptakes constrain the simulations more effectively than the oceanic ones. These constraints, while useful, are still unable to rule out both extremely strong positive and modest negative climate–carbon cycle feedback.     

The effect of atmospheric CO2 and ice sheet topography on LGM climate

The role of reduced atmospheric CO₂ concentration and ice sheet topography plus its associated land albedo on the LGM climate is investigated using a coupled atmosphere-ocean-sea ice climate system model. The surface cooling induced by the reduced CO₂ concentration is larger than that by the ice sheet topography plus other factors by about 30% for the surface air temperature and by about 100% for the sea surface temperature. A large inter-hemispheric asymmetry in surface cooling with a larger cooling in the Northern Hemisphere is found for both cases. This asymmetric inter-hemispheric temperature response is consistent in the ice sheet topography case with earlier studies using an atmospheric model coupled with a mixed-layer ocean representation, but contrasts with these results in the reduced CO₂ case. The incorporation of ocean dynamics presumably leads to a larger snow and sea ice feedback as a result of the reduction in northward ocean heat transport, mainly as a consequence of the decrease in the North Atlantic overturning circulation by the substantial freshening of the North Atlantic convection regions. A reversed case is found in the Southern Ocean. Overall, the reduction in atmospheric CO₂ concentration accounts for about 60% of the total LGM climate change.