Organic carbon buried in sediments as coal, natural gas, and oil over literally hundreds of millions of years is being consumed as a result of human activities and returned to the atmosphere as carbon dioxide (CO2) on a time scale of a few centuries. The energy harvested from this transformation of fossil fuels supplies us with electricity, heat, transportation, and industrial power. The clearing of forests for agricultural lands and the harvesting of wood, both of which remove carbon-bearing vegetation, have also added CO2 to the atmosphere, in amounts equivalent to more than half of the fossil fuel source. The CO2 added to the atmosphere because of man’s activities, and the way it is currently distributed within the land, air, and sea, is depicted in the carbon cycle diagram shown in figure 1.
Because of anthropogenic emissions, atmospheric CO2 has climbed to levels that are presently more than 30% higher than before the industrial revolution, 1,2 as seen in figure 2. Indeed, geochemical measurements made on ancient ocean sediments suggest that atmospheric CO2 levels over the past 20 million years were never as high as they are today. 3
The increase in atmospheric CO2 has drawn a great deal of attention because of the impact it has on the trapping of long wavelength radiation emitted from Earth’s surface. More than half of the increase in the direct trapping that has occurred since preindustrial times is attributed to CO2, with the rest coming from other gases such as methane, nitrous oxide, and chlorofluorocarbons. The effect of this increased trapping on Earth’s climate depends on a number of complex feedbacks. Nevertheless, the strong consensus of the scientific community is that the increased trapping will lead to global warming; it probably accounts for most of the 0.6 ± 0.2 °C warming that occurred during the last century. 3 Humankind thus appears to be playing a significant role in altering Earth’s climate.
Because CO2 is nonreactive in the atmosphere, it has a relatively long residence time there. However, its growth rate is presently less than half of what would be expected if all the CO2 released by fossil-fuel burning and land-use change remained in the atmosphere (figure 3). The growth rate is lower because the terrestrial biosphere (plants and soils) and the ocean are taking up a significant amount of anthropogenic CO2, that is, acting as “sinks.” The scientific community has made much progress in establishing the relative role of these two major natural sinks on a global scale, and it appears that the missing carbon is about equally divided between them. However, scientists continue to debate aspects of the spatial distribution and mechanisms of these sinks. The future behavior of the sinks turns out to be highly sensitive to whatever mechanisms we assume. Thus, better understanding of their behaviors is key to predicting, and hopefully mitigating, the future impact of anthropogenic CO2. An important starting point for forecasting the future behavior is to understand its past.
Global carbon balance in the recent past
The most recent decade for which we have an estimate of all the carbon sources and sinks 3 is the 1980s, during which the average fossil-fuel emissions were estimated at 5.4 ± 0.3 petagrams of carbon per year (1 Pg = 1015 g) and the atmospheric growth rate at 3.3 ± 0.1 Pg C/yr. The difference of 2.1 ± 0.3 Pg C/yr must be taken up by the ocean and terrestrial biosphere. New developments using observations of oxygen and carbon isotopes in the atmosphere allow us to partition the sink between these two reservoirs, as described in box 1. The partitioning gives an average net oceanic uptake of 1.9 ± 0.6 Pg C/yr and a net land carbon uptake of only 0.2 ± 0.7 Pg C/yr. However, there is an additional terrestrial source of CO2 to the atmosphere, mainly due to deforestation in the tropics, estimated to be about 1.7 Pg C/yr. One can infer that the land must be taking up 1.9 Pg C/yr outside the areas affected by deforestation and hence appears to be as large a carbon sink as the oceans, as shown in figure 1. Note, however, that the uncertainty in the estimated terrestrial source is very large, with a reported range of 0.6–2.5 Pg C/yr. The uncertainty in the land sink, which is estimated by differencing all the other components of the carbon budget, is correspondingly large.
The carbon balance is not static in time. As seen in figure 3, the atmospheric growth rate varies by a large amount from year to year. Most of the interannual variability is correlated with the El Niño southern oscillation climate mode, with higher growth rates generally being associated with El Niño (warm climate) episodes. The climate cooling caused by the Mt. Pinatubo eruption in the early 1990s appears to have contributed to reduced atmospheric growth rates. The primary cause of the variability remains controversial, but is probably due mostly to the response of terrestrial vegetation to climate variability, with a smaller contribution due to the oceanic response. 4
The total capacity of the oceans for taking up anthropogenic CO2 is a function primarily of the solubility of CO2 and the chemical buffering capacity of seawater. The chemical buffering stems from the reaction of CO2 molecules with carbonate ions to produce bicarbonate ions. If a moderately sized pulse of CO2 were added to the atmosphere today, about 85% of it would dissolve in the ocean, but the process would take more than 1000 years because of the sluggishness of vertical exchange between the surface and interior of the ocean. For example, the oldest ocean water in the world, which is found in the deep North Pacific, has been out of contact with the surface of the ocean for about 1000 years. Estimates of the oceanic concentration of anthropogenic carbon, shown in figure 4, demonstrate that, except in a few locations, carbon from human activities generally has not yet penetrated in significant amounts below a depth of about 1000 m. Surface waters carrying anthropogenic carbon descend into the abyssal ocean primarily in the North Atlantic and Southern Oceans, as part of the ocean’s thermohaline circulation, the so-called “conveyor belt.” From there, those waters spread, over the next centuries, throughout the entire deep ocean. Currently one can see a modest signal of the anthropogenic carbon in deep waters only in the Atlantic north of the equator.
Model simulations suggest that the Southern Ocean around the Antarctic (south of 35 °S) accounts for nearly half of the net air-sea flux of anthropogenic carbon. Another sixth of the total uptake occurs in the tropics (13 °S to 13 °N). Those are both regions with large surface areas, where upwelling (in the tropics and Southern Ocean) and deep vertical exchange (in the Southern Ocean) bring to the surface older, relatively uncontaminated deep water that has a high capacity for uptake of CO2. From those regions, a large part of the absorbed CO2 is transported laterally toward the subtropics, producing a relatively uniform distribution of anthropogenic carbon, as shown in figure 4. Different methods for estimating the uptake of anthropogenic CO2 agree reasonably well with each other and with observational constraints in most regions except the Southern Ocean, where the relatively small-scale processes that destabilize the water column and lead to vertical overturning and formation of deep water are poorly understood and particularly difficult to model.
Understanding the land sink
The land sink for carbon is the subject of considerable controversy at present, concerning not only its magnitude but also its cause. For many years, researchers have believed that the dominant sink mechanism is the fertilizing effects of increased CO2 concentrations in the atmosphere and the addition to soils of fixed nitrogen from fossil-fuel burning and agricultural fertilizers. This fertilization mechanism has been incorporated into most existing models of the terrestrial biosphere that are used to predict future concentrations of atmospheric CO2. However, a recent analysis of long-term observations of the change in biomass and growth rates, made by the US Forest Service, suggests that such fertilization effects are much too small to explain more than a small fraction of the observed sink in the US. 5 In addition, long-term experiments in which small forest patches and other land ecosystems have been exposed to elevated CO2 levels for extended periods show a rapid decrease of the fertilization effect after an initial enhancement. 6
What other mechanisms might then explain the existence of a large land sink? One of the few land regions of the world where the observational coverage and model analysis are sufficient for putting together a carbon budget is the coterminous US. A recent estimate of carbon uptake in this region 7 shows that the largest sink is due to regrowth in abandoned farmland and areas that had previously been logged. 5 As shown in figure 5, a large fraction of the forests in the eastern US were cut down before 1920. As the use of agricultural land shifted westward over time, land in the east that had been used for agriculture was abandoned and started to recover its forests. Moreover, as the logging pressure has decreased in the east over time, there has been a large increase in secondary vegetation.
After regrowth in farmlands, the second most important sink mechanism identified by the study of carbon uptake in the US is the encroachment of woody growth into areas where fires have been suppressed, especially in the southwest, as seen in figure 5. The third largest mechanism is the growth of food and wood that is exported from the US before being consumed. But that is not a true sink because the CO2 is returned to the atmosphere elsewhere in the world. Other sinks include storage of carbon in wood products and cropland soils, as well as in freshwater reservoirs, alluvium, and colluvium. Some carbon is taken up by weathering, such as the reaction of rock with CO2 and water to form dissolved bicarbonate and calcium, but it is exported from the US by rivers.
What do we know about the land sink in the rest of the world? For the past decade, “inverse modeling” 8 has been one of the primary tools to improve our understanding of the large-scale behavior of both the terrestrial and the oceanic sinks. Inverse models use information about how CO2 is currently distributed in the atmosphere and knowledge of how the atmosphere transports CO2 to infer the spatial distribution of carbon sources and sinks at Earth’s surface. From such models, as discussed in box 2, we have learned that more than half of the sink for anthropogenic carbon must be in the Northern Hemisphere. Furthermore, a majority of the Northern Hemisphere sink appears to be in the terrestrial biosphere.
While many attempts have been made to determine the location of this Northern Hemisphere land sink more precisely, the results show considerable variation among the various inversion models. Unfortunately, in situ verification of those estimates has proven to be extremely challenging. Despite the development of innovative new methods to measure in situ fluxes (as seen on the cover and as described at http://www.as.harvard.edu/chemistry/hf) and detailed distributions of atmospheric CO2 from aircraft, it continues to be difficult to close the gap between lower in situ estimates and the large sinks inferred from the models. The estimate of the land sink in the coterminous US discussed previously gives a total carbon sink that is on the low side but within the uncertainty of what is implied by atmospheric inversions. Although the US data suggest that the land sink there is driven primarily by land-use history and not fertilization by atmospheric CO2 or atmospheric nitrogen deposition, it is not clear if the same findings hold true for other land regions. If they do, however, they have major implications about the future land sink: The global CO2 uptake capacity from changes in land-use such as forest regrowth is much more modest than that predicted for fertilization. 9
To establish the redistribution of the anthropogenic carbon dioxide within the global carbon cycle, one ideally would measure the changes in the carbon inventory of all relevant reservoirs—atmosphere, ocean, and terrestrial biosphere—and then compare these changes with the total anthropogenic emissions. Unfortunately, the atmosphere is the only reservoir that is readily amenable to such an approach: It is relatively well mixed so that observations at a few key locations give a very accurate measure of its changes over time, as seen in figure 2. By contrast, mixing in the ocean is very slow, and the invasion of anthropogenic CO2 there is masked by a substantial amount of natural background variability in the ocean’s inorganic carbon system. Furthermore, very accurate measurements of dissolved inorganic carbon in the ocean have become available only over the past two decades, and are still insufficient to address changes over time. Measuring the carbon inventory changes in the terrestrial biosphere is even more difficult, given its inherent large spatial and temporal inhomogeneity. Therefore, up to now, researchers have had to rely primarily on indirect methods to determine the partitioning of anthropogenic CO2 between the ocean and terrestrial biosphere.
Until the early 1990s, ocean models were the primary tool for investigating how the carbon was partitioned. The many indirect approaches that have been developed since then have led to a substantial improvement in the constraints available. Among the most powerful of these is the atmospheric oxygen method developed by Ralph Keeling of the Scripps Institution of Oceanography. 13 Oxygen and CO2 are exchanged in relatively fixed stoichiometric ratios during the burning of fossil fuels as well as during photosynthesis and respiration by plants, animals, and bacteria. By contrast, the exchange of CO2 between the ocean and the atmosphere is independent of the exchange of oxygen between the two fluids. Therefore, only uptake of CO2 by the terrestrial biosphere will leave an imprint on atmospheric oxygen. By measuring the changes in atmospheric oxygen and knowing the fossil-fuel emissions and the exact values of the stoichiometric ratios, one can separate the ocean uptake from that of the terrestrial biosphere.
Although the basic principle of the oxygen method is very elegant, it could not be implemented until the development of an instrument capable of measuring the resulting tiny changes in atmospheric oxygen. Ralph Keeling, while working as a graduate student at Harvard University, came up with an interferometric technique that gave sufficient precision. More recently, the development of other methods based on mass spectrometry or the paramagnetic properties of oxygen has substantially decreased the complexities involved in measuring the subtle changes in atmospheric oxygen. The oxygen-based analysis of the carbon budget has been adopted by the Intergovernmental Panel on Climate Change (IPCC) as the primary means to determine the redistribution of anthropogenic CO2. Unfortunately, the method may have somewhat larger uncertainties than previously thought because it assumes no net exchange of oxygen between the ocean and atmosphere. By contrast, recent modeling and data-based studies are finding that such air–sea oxygen fluxes are fairly sensitive to temperature changes.
A second important method for determining the partitioning of anthropogenic CO2 is based on measurements of the ratio of two isotopes, carbon-13 to carbon-12, in the atmosphere and ocean. Fossil-fuel CO2 is slightly depleted in 13C relative to 12C, leading to a long-term decrease of the 13C/12C isotopic ratio of atmospheric CO2. As anthropogenic CO2 is redistributed within the global carbon cycle, it brings with it this isotopic “color,” which permits us to distinguish anthropogenic CO2 from the natural CO2. Several variations of this basic concept have been developed and successfully applied, but they are currently limited by the scarcity of data, the exquisite precision and calibration required of the 13C/12C measurements, and the limited quantitative knowledge of the fractionation of 13C relative to 12C by land plants.
Despite the uncertainties associated with the different indirect methods for detecting the magnitude of the terrestrial and oceanic carbon sinks, it is encouraging to note that almost all of these methods have converged toward the values adopted by IPCC. 3 We believe that the uncertainties will further decrease in the future because of the ongoing development of improved methods for estimating fluxes and inventories over land, and because a recently completed large-scale CO2 survey of the ocean has given scientists, for the first time, an oceanwide high-quality data set of the oceanic CO2 distribution that can serve as the basis for direct documentation of the cumulative increase in anthropogenic CO2 in the ocean. Given enough time, we will be able to account for anthropogenic CO2 in both the atmosphere and the ocean directly, and hence reduce the need to rely on indirect methods and likely reduce the uncertainties associated with the anthropogenic CO2 budget.
Predicting the future
The usual approach for evaluating the future trajectory of atmospheric CO2 is to predict potential concentrations under each of several possible emissions scenarios. The emissions scenarios are beset with uncertainties: How will population changes, economic growth, technology developments, and so forth affect the amount of CO2 emitted? And how likely is it that nations will agree to curb those emissions? Models for determining the resulting atmospheric concentrations are further beset by uncertainties about the land and ocean sinks.
A common starting point for climate predictions is a scenario for anthropogenic emissions of CO2 and other greenhouse gases known as IS92a; it was developed by the Intergovernmental Panel on Climate Change (IPCC) to represent “business-as-usual,” with no intentional efforts to mitigate the greenhouse gas. This scenario estimates total fossil-fuel emissions between 2000 and 2100 of 1450 Pg C—a substantial fraction of the estimated remaining resource base of 3500 Pg C. 10
To estimate how much of this CO2 stays in the atmosphere, one must use models of the global carbon cycle. Until about the mid-1990s, most such models assumed that the climate remains stationary and that the land sink is primarily driven by CO2 fertilization. Many such models project that about half of the CO2 emissions will have accumulated in the atmosphere by the year 2100, and will lead to atmospheric CO2 concentrations of around 700 parts per million (ppm). How the rest of the carbon is partitioned between land and ocean varies widely with the model applied: 3 In one particular simulation by the British Met Office’s Hadley Centre, 11 60% of the remaining CO2, or 450 Pg C, ends up on land and the remaining 300 Pg C in the ocean.
Such predictions can change dramatically if one considers feedbacks due to global warming. In another run of their climate model, Hadley Centre researchers calculated that an increase of atmospheric CO2 concentrations to 713 ppm by 2100 would cause global temperatures to rise by 4.0 °C, (5.5 °C over land) with an associated intensification of the hydrological cycle. In that simulation, the climate change is allowed to influence the terrestrial and oceanic carbon sinks but not the atmospheric concentrations, which remain fixed. As the climate warms, the land no longer absorbs carbon, but emits it; in other words, the land sink turns into a source, adding 60 Pg C to the atmosphere by 2100. At the same time, the oceanic sink becomes less effective, soaking up 250 Pg C, rather than 300 Pg C. The changed nature of the land sink in this model results first from the disappearance of the Amazonian rainforest due to increased dryness of soils and second from an increased rate of bacterial breakdown of organic carbon in soils due to increased soil temperatures. The oceans lose some solubility due to higher water temperatures, slowed vertical exchange, and changes in the way that biological processes redistribute carbon within the ocean. The combined terrestrial and oceanic carbon sinks predicted by the Hadley model for 2100 are thus 190 Pg C rather than 750 Pg C, and this simulation is out of balance with the IS92a emissions by about 560 Pg C.
In a third simulation, Hadley Centre researchers put the excess 560 Pg C back into the atmosphere using a new fully coupled simulation of the carbon cycle and climate. Concentrations of atmospheric CO2 climb to 980 ppm instead of 713 ppm, and the surface air temperature warming is 5.5 °C (8.0 °C on land). The greater warming induces the land to release even more carbon, but that release is counterbalanced by an increased oceanic absorption caused by the faster growth rate of atmospheric CO2, as illustrated in the top panel of figure 6. All the additional CO2 remains in the atmosphere.
The results are as uncertain as they are disconcerting. The bottom panel of Figure 6 shows the results of another simulation, done by the Institut Pierre Simon Laplace (IPSL), in France, which also predicts the impact of changing climate on the atmospheric carbon concentrations. 12 The IPSL model uses a newer emissions scenario, similar to the IS92a. The model’s oceanic uptake is almost two times more effective than that of the Hadley model ocean, with most of the difference being in the difficult-to-simulate Southern Ocean. More important, however, the land uptake in the IPSL model is less than half as sensitive to temperature as the Hadley model, and does not allow the distribution of terrestrial vegetation to change. The IPSL model with global warming thus predicts an atmospheric CO2 level in 2100 that only reaches 780 ppm, as contrasted with 980 ppm for the Hadley model, as shown in the lower panel of figure 6.
We scientists cannot hope to narrow the range between such simulations until we have a better understanding of the fundamental mechanisms that control the carbon sinks. The uncertainties that must be addressed include the magnitude of CO2 uptake in the Southern Ocean and the worldwide impact of fertilization and land-use changes.
The spatial and temporal distribution of carbon dioxide sources and sinks produces small gradients in the observed distribution of its atmospheric concentration. Much of what we currently know about the spatiotemporal pattern of fluxes of CO2 between the atmosphere and Earth’s surface, in particular for land regions, comes from analysis of the observed distributions. In the inversion method, the land and ocean surfaces are first divided into a modest number of regions. Next, atmospheric transport models are used to estimate the concentration footprint that would result if a given region emitted a single unit flux of CO2. The transport models are tested using tracers with known source-sink distributions. Various regression methods are then used to determine the weighted sum of footprints that best reproduces the observations. The weights represent the flux contributions from each region.
By far the largest signal in the annual mean atmospheric CO2 distribution is the north-south gradient, which comes about because almost 90% of the global fossil-fuel emissions occur north of the tropics. The measured concentration gradient is about 3 parts per million (ppm) between the two poles. That’s substantially smaller than the difference of about 5 ppm one would expect if the carbon sinks were evenly distributed around the world. This discrepancy implies the existence of strong sinks for atmospheric CO2 in the Northern Hemisphere. Two contrary interpretations have been put forward to explain the nature of these Northern Hemisphere sinks.
One interpretation, put forward in 1989 by Charles D. Keeling of the Scripps Institution of Oceanography, working with Martin Heimann and colleagues, 14 is that this apparent sink is not driven by anthropogenic processes, but to a large extent reflects part of the preindustrial (natural) carbon cycle in the ocean. They postulated that an amount of the order of 1 Pg C/yr of CO2 is taken up in the northern North Atlantic and transported toward the Southern Ocean as part of a large-scale ocean circulation. As the waters upwell in the Southern Ocean, they tend to give this CO2 back to the atmosphere. To force a return transport pathway back to the Northern Hemisphere, the CO2 must have built up to higher levels in the atmosphere in the Southern Hemisphere than the Northern Hemisphere in preindustrial times. This hypothesis seemed to be confirmed by an extrapolation of the interhemispheric CO2 gradient back to a time when fossil-fuel emissions were zero.
Around the same time that Keeling and his colleagues made their proposal, an independent explanation for the reduced gradient was proposed by Pieter Tans of the National Oceanic and Atmospheric Administration in Boulder, Colorado, Inez Fung, now at the University of California, Berkeley, and Taro Takahashi of Columbia University’s Lamont Doherty Earth Observatory. 15 Using observations of the difference in the partial pressure of CO2 between the atmosphere and the surface ocean in the North Atlantic and in a few other oceanic regions as additional constraints in their atmospheric inversions, they concluded that the apparent Northern Hemispheric sink cannot be attributed to the ocean and must be almost exclusively driven by the terrestrial biosphere.
These two strongly diverging views regarding the relative role of the oceans and terrestrial biosphere in the global carbon cycle dominated most of the carbon-cycle research discussions during the 1990s. Recently, a combined analysis of new oceanic and atmospheric observations has lent support to the view that this sink is primarily driven by the terrestrial biosphere, with a smaller contribution by the ocean. So far, there is little evidence for a strong interhemispheric transport of CO2 by the oceans.
Atmospheric inversion studies have raised additional puzzles, which remain to be solved. One puzzle is posed by evidence indicating that the Southern Ocean is relatively neutral with respect to the atmospheric carbon budget, even though model simulations suggest it is a strong sink. As another puzzle, tropical lands also appear to be neutral, despite the large deforestation occurring there. Hopefully, extended inversion models now under development, when combined with improved observations, will provide insight into these puzzles.
No magic bullets
Although we don’t yet fully understand the global carbon cycle, it’s safe to say there are no magic bullets in the carbon sinks to rescue the world from high atmospheric CO2 levels any time in the next few centuries. Quite to the contrary, most of the feedbacks between the global carbon cycle and global warming seem to be positive—that is, global warming reduces the sink strengths. Furthermore, the land carbon sink, if it is due primarily to land-use changes, will likely saturate much earlier than if it were due largely to fertilization. Thus, the emission reduction and sequestration measures required to stabilize future atmospheric CO2 at levels that will minimize impacts are likely to be more stringent than estimated with models that did not include carbon-climate feedbacks and that did include fertilization as the main mechanism for the terrestrial carbon sink.
Some researchers have proposed enhancing the natural carbon sinks by, for example, reforestation, fertilizing with missing nutrients such as iron to stimulate the oceanic biological pump, and injecting CO2 directly into the deep ocean. But the projected business-as-usual emissions of order 1450 Pg C over this century are far too large for the first two measures to alter significantly the final fate of CO2. Deep ocean injection might reduce the future peak value of CO2 in the atmosphere, but would not affect the final atmosphere–ocean equilibration. Moreover, it would require extremely large injections, with unknown technological challenges and environmental consequences. Geological sequestration is also being explored. Currently, the only method that can be guaranteed to mitigate the expected increases in atmospheric CO2 is reduced fossil-energy utilization.
The past two decades have seen an improvement in our understanding of the global carbon cycle. We now know the history of atmospheric CO2 for many hundreds of thousands of years. We have a good grasp of the partitioning of the present carbon sink between land and ocean, and strong constraints on the latitudinal distribution of those sinks. However, there is much need for improvement. On land, new findings call into question the fundamental mechanisms to which the historic terrestrial carbon sink had been ascribed. If the sink is primarily due to land-use changes, it will likely saturate much earlier than if due to CO2 fertilization.
Regarding the oceans, we understand their chemistry and most of the fundamental mechanisms driving the ocean circulation and can develop sophisticated models of them. However, we have to settle for only a simplified representation of such important related processes as convective overturning and the formation and melting of sea ice, although progress is being made in these areas. We also face limitations in computer capabilities. Consequently, simulations of both the present uptake and response of the ocean to future warming using different models disagree substantially. The problem is particularly acute in the Southern Ocean around Antarctica, which we believe may account for 40–50% of present uptake of anthropogenic CO2 and will continue to dominate the uptake in the future. In addition, we have only begun to understand how the cycling of carbon by organisms in the ocean will respond to climate and chemical changes.
The international scientific community is well aware of the problems we have identified and has developed long-term observational and modeling strategies for tackling them. Those strategies include proposals for improved observational capabilities on land and in the ocean, with a particular focus on monitoring the changes that the carbon system will undergo in response to warming. Improved models of the coupled carbon–climate system are being developed at a number of locations around the world. Considerable additional progress can be expected over the next decade.
REFERENCES
Jorge L. Sarmiento is a professor in the department of geosciences at Princeton University. Nicolas Gruber is an assistant professor in the department of atmospheric sciences and in the Institute of Geophysics and Planetary Physics at the University of California, Los Angeles.