Technical Comments

North American Carbon Sink

Science  19 Mar 1999:
Vol. 283, Issue 5409, pp. 1815
DOI: 10.1126/science.283.5409.1815a

S. Fan et al. (1) suggest that 1.4 ± 0.4 petagrams of carbon per year (Pg C year−1) is taken up by the forests in North America, in contrast to 0.1 Pg C year−1 taken up in Eurasia. Fan et al. invoke reforestation and regrowth, fertilization by anthropogenic N deposition, global warming, and CO2 fertilization to substantiate their inverse-model calculation. However, mechanistic models and measurements in terrestrial ecosystems do not agree with either the magnitude or spatial distribution of the CO2 sink proposed by Fan et al.

Direct estimates of forest C uptake—which reflect the interacting effects of rising CO2, N fertilization, climatic changes, as well as reforestation and regrowth—indicate that forests in 28 eastern U.S. states during the late 1980s to early 1990s had an estimated net C uptake of only 0.17 Pg year−1 aboveground (2). This estimate includes the southeastern forests, which are among the most productive forest in the United States. European and Russian forest inventories (3) suggest comparable C uptake, which contradicts the North American distribution of the sink proposed in the report (1).

Anthropogenic N deposition stimulates C uptake in North America by 0.29 to 0.35 Pg C year−1 with significantly more uptake in Eurasia (0.67 to 0.86 Pg C year−1) according to model calculations (4). Moreover, recent field experiments by Nadelhoffer et al. (5) suggest that even this estimate of fertilized CO2 uptake may be too high. Modeled CO2 and climate effects generate a Northern Hemisphere sink with no significant east-west bias and a magnitude of ∼0.58 Pg C year−1 (6).

A robust understanding of the global carbon budget requires reconciliation of ecological mechanisms with inverse estimates of the spatial distribution of the so-called missing sink. Fan et al. concede that their analysis approaches the limits of uncertainty, but the lack of corroboration from independent observations suggest that they may have overextended these limits.

REFERENCES AND NOTES

Table 21

Estimated terrestrial C uptake (Pg C year−1) in the mid-latitude Northern Hemisphere, based on various data and models.

View this table:

The report by S. Fan et al. (1) has generated an important debate about their inference of an annual C sink as large as 1.7 ± 0.5 Pg in terrestrial ecosystems of North America. A terrestrial sink flux of this magnitude could completely offset a continental emission source from fossil fuel of 1.6 Pg C year−1. Fan et al. applied atmospheric constraints in their inverse model approach and used monthly C fluxes from the “equilibrium” version (without year-to-year variability) of the Carnegie-Ames-Stanford Approach (CASA) biosphere model mainly to estimate the magnitude of the seasonal “rectifier effect” in the atmosphere. We took a fundamentally different approach to study these questions (2). We used newly derived terrestrial C fluxes, predicted directly from forward-modeling CASA simulations in a nonequilibrium mode (using observed interannual variability for surface climate and satellite imagery). Our results (2) imply a much different conclusion from that of the report (1) with respect to the North America land sink during the late 1980s.

In accordance with the findings by Fan et al. at the global scale, we calculate that the worldwide C sink from net ecosystem production (NEP) can vary between 0.4 and 2.6 Pg C year−1in the terrestrial biosphere. However, by making NEP estimates at 1° resolution directly from the CASA biosphere model, we found that C sink fluxes in terrestrial ecosystems of the United States and Canada totaled only 0.12 and 0.10 Pg C in 1987, respectively, and 0.05 and 0.17 Pg C in 1988. The land area covered by states of the former Soviet Union shows a larger NEP sink of 0.4 to 0.6 Pg C year-1during 1987 to 1988. These years may in fact represent two of most favorable of the 1980s for North American sink fluxes, owing to temperature warming trends.

We also included in the CASA model routines to determine aboveground biomass with the use of typical residence times in pools of wood and leaves to estimate changes in forest regrowth and deforestation C fluxes. This analysis, which relies on country-by-country changes in global forest cover for the years 1990 to 1995 as compiled by the United Nations (3), suggests that net forest regrowth can add about 0.09 Pg C annually to the North America C sink. Thus, it appears from forward CASA ecosystem modeling that total annual accumulation of atmospheric C in terrestrial ecosystems of North America could offset only about 20% of continental C source from fossil fuel burning in the late 1980s.

While analysis of CASA's satellite-driven net primary production (NPP) and soil heterotrophic CO2 fluxes suggests that regional warming in Northern Hemisphere forests can enhance ecosystem production notably, our prediction of the regional distribution for this NEP sink over North America and Eurasia is not in agreement with the geographic patterns reported by Fan et al. for the late 1980s. It seems unlikely that the intercontinental balance of a large terrestrial C sink could shift so rapidly over just a few years, as suggested by results in the report. In any case, the potential for high interannual variability implies merely a transient sink pattern in North America. We infer that (in relation to results from the CASA model) climate variability, recent forest regrowth, and increased greenness as recorded in the satellite data capture the major processes of a terrestrial ecosystem sink.

REFERENCES AND NOTES

Response: Our report (1) confirmed previous estimates of the size of a terrestrial C sink in the mid-latitude Northern Hemisphere and suggested that most of the sink occurred in North America from 1988 to 1992. Partitioning of the terrestrial sink among the continents is controversial because the CO2 observations are sparse and the models are imperfect. In our report we mentioned possible causes of the sink, but we did not identify it with any specific mechanism because our method does not allow us to do so. Potter et al. and Holland and Brown [in these comments and elsewhere (2) with colleagues] estimated terrestrial C uptake by invoking specific causes. The difference between these estimates and ours is the subject of intense research and may be resolved by a growing network of atmospheric and ecological observations.

Current research on the sources and sinks of atmospheric CO2 consists of four major activities: (i) long-term continuous measurement of C fluxes between the atmosphere and the biosphere by eddy correlation, (ii) repeated measurements over time of C inventories in terrestrial ecosystems and in the ocean, (iii) long-term monitoring of atmospheric CO2 in a global air sampling network, and (iv) development of mechanistical ecosystem models that predict past, present, and future C cycles on land and in the sea. These activities have been conducted separately in different laboratories and have produced independent estimates of the global carbon budget.

The estimates of oceanic uptake of CO2 by different approaches appear to converge. For instance, a global uptake of about 2 Pg C year−1 is estimated for the 1980s separately by ocean circulation and biogeochemistry models (from changes in dissolved inorganic C over time) and by measurements of the air-sea difference of CO2 partial pressure. The estimated air-sea exchange fluxes show similar spatial distributions. More recent measurements of O2/N2 and δ13CO2 in the atmosphere are consistent with the ocean sink estimates.

However, estimates of terrestrial C sequestration by different methods are disparate in both their magnitudes and spatial distributions. Terrestrial ecosystem models (such as the CASA model discussed by Potter) relate the NEP to environmental parameters (light; air; soil temperature, moisture, and N; and the ambient CO2 mixing ratio) and to ecological properties (land use history, stand age, species of vegetation, and leaf area index), and predict C fluxes resulting from succession after disturbance, interannual changes in climate pattern, increasing CO2 in the atmosphere, and input of N from atmospheric deposition. Some of these models predict large global NEP or C uptake rates as a result of CO2fertilization. Most of this uptake is predicted to occur in the tropical forests, with a small contribution from the mid-latitude Northern Hemisphere (for example, 3, 4). But the response of plant photosynthesis to ambient CO2concentration is uncertain in the natural environment. The treatment of land-use changes and natural disturbances is difficult in a time-dependent terrestrial ecosystem model because of a lack of historical data in most of the world. As a result, we cannot assess the reliability of NEP estimates from any of the global ecosystem models including that of Potter and Klooster (4).

A fertilization of the ecosystems by atmospheric N deposition may also cause a terrestrial C sink. Nitrogen deposition is heaviest in Europe and eastern North America. Earlier estimates of this sink are as large as 1.5 Pg C year−1 on a global scale, if one assumes a 100% utilization of the N for plant CO2 assimilation (5). However, recent studies indicate that only one-third of the deposited N may be used by plants, while the rest may be tied up in soil organic matter (6). Plant tissues have a much higher ratio of C:N than does soil organic matter.

Measurements of forest C inventory can provide a direct estimate of forest C uptake. The estimate of Brown and Schroeder (as cited in the comment by Holland and Brown) is similar to previous analyses of forest inventory data: the forest ecosystems in North America sequester carbon at a rate < 0.2 Pg C year−1. The global C sequestration is estimated to be 0.8 Pg C year−1 by forest ecosystems, more than half of the estimated sink is located in Europe and Russia.

In contrast, atmospheric observations of CO2 and the ratios of O2/N2 and13CO2/12CO2, combined with numerical models of atmospheric transport, have consistently implied the presence of a terrestrial C sink in the mid-latitude Northern Hemisphere of 1 to 3 Pg C year−1 (Table 1). Our report supports the many previous estimates of this sink (1.4 ± 0.2 Pg C year−1 by SKYHI and 2.2 ± 0.2 Pg C year−1 by GCTM). Table 1 compares various estimates of the terrestrial C uptake in the mid-latitude Northern Hemisphere. The uptake rate implied by atmospheric and oceanic CO2 data and models tends to be much larger than estimates based on forest inventory data and mechanistic terrestrial ecosystem models.

We stated in our report that the spatial distribution of the mid-latitude terrestrial C sink is not as well constrained by the atmospheric CO2 data. As shown in figure 2 in our report, we cannot reject the presence of a substantial Eurasian sink. One-half of the estimated Northern Hemisphere terrestrial sink is attributable to Eurasia at the limit of 1 standard deviation (SD), and two-thirds at the limit of 2 SD. However, we estimate a very low probability for the North American sink to be as small as 0.2 Pg C year−1 for 1988 to 1992, as suggested by Potter and Klooster (4) and by Brown and Schroeder (2).

The intercomparison between our inverse estimates (based on atmospheric and oceanic data and models) and estimates based on mechanistic models and ground data must take into consideration the large temporal variability of the terrestrial C cycling (7). The growth of atmospheric CO2 decreased from 4.9 to 1.4 Pg C year−1 from 1988 to 1992, while the global fossil CO2 emission remained nearly constant (5.9 Pg C year−1 in 1988 to 6.1 Pg C year−1 in 1992). This implies an increase of the global ocean and terrestrial C sink from 1.0 Pg C year−1 in 1988 to 4.7 Pg C year−1 in 1992. According to model calculations, the oceanic uptake could not have changed by more than 1 Pg C year−1. The large interannual CO2 variability must have been caused mainly by terrestrial ecosystem productivity and respiration.

The estimate of Potter and Klooster covers a different period (1987 and 1988). Given the large interannual variability of terrestrial NEP, their estimate should not agree with ours even if both are correct. Also, it is possible that our estimate, if correct, might apply to the restricted years studied and might disagree with correct long-term estimates from inventory data.

Last, but not least, direct CO2 flux measurements by eddy correlation have shown evidence of significant C uptake in a deciduous forest in New England that had its last major disturbance by a hurricane some 60 years ago (8). A network of the tower-based measurements has been implemented in North America (AmeriFlux) and Europe (EuroFlux), and is being extended to South America (and remains to be extended to Africa and Aia). Long-term C flux and ecological data collected from these tower sites will provide critical information on the size and causes of the terrestrial C sink from landscape to global scales. This is an exciting development in C cycle research that we hope will soon resolve the discrepancies between the various estimates of terrestrial C sequestration.

REFERENCES

Navigate This Article