Drought-Induced Reduction in Global Terrestrial Net Primary Production from 2000 Through 2009

See allHide authors and affiliations

Science  20 Aug 2010:
Vol. 329, Issue 5994, pp. 940-943
DOI: 10.1126/science.1192666


Terrestrial net primary production (NPP) quantifies the amount of atmospheric carbon fixed by plants and accumulated as biomass. Previous studies have shown that climate constraints were relaxing with increasing temperature and solar radiation, allowing an upward trend in NPP from 1982 through 1999. The past decade (2000 to 2009) has been the warmest since instrumental measurements began, which could imply continued increases in NPP; however, our estimates suggest a reduction in the global NPP of 0.55 petagrams of carbon. Large-scale droughts have reduced regional NPP, and a drying trend in the Southern Hemisphere has decreased NPP in that area, counteracting the increased NPP over the Northern Hemisphere. A continued decline in NPP would not only weaken the terrestrial carbon sink, but it would also intensify future competition between food demand and proposed biofuel production.

Terrestrial ecosystems are a major sink in the global carbon cycle, sequestering carbon and slowing the increasing CO2 concentration in the atmosphere (1). Terrestrial net primary production (NPP), the initial step of the carbon cycle in which carbon is fixed as biomass, increased from 1982 through 1999, in part due to eased climatic constraints on plant growth (2). The World Meteorological Organization (WMO), National Oceanic and Atmospheric Administration (NOAA), and NASA all reported that 2000 to 2009 was the warmest decade since instrumental measurements of temperatures began in the 1880s (3). We questioned whether the warming climate of the past decade continued to increase NPP, or if different climate constraints were more important.

Between 2000 and 2008, CO2 emissions from fossil fuel combustion continued to increase at a rate consistent with the average of the highest-emissions family of scenarios, A1FI, used by the Intergovernmental Panel on Climate Change in the Fourth Assessment (1). Carbon-budget methods show that the land is becoming a stronger carbon sink, whereas large uncertainties exist in the partitioning of ocean and land carbon-sink components (1, 4). Satellite data can generally provide realistic information on vegetation dynamics, including land cover change (5, 6), disturbances, and recovery (7), which may help to reduce uncertainties in carbon-budget estimates. In this study, we investigate terrestrial NPP and climate variability over the past decade (2000 to 2009) by analyzing satellite data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board NASA’s Terra satellite and global climate data.

We used the global MODIS NPP algorithm (8) [see supporting online material (SOM) text S1] to examine spatially explicit NPP changes from 2000 through 2009. We used collection 5 (C5) 8-day composite 1-km fraction of photosynthetically active radiation (FPAR) and leaf area index (LAI) data from the MODIS sensor (9) as remotely sensed vegetation property dynamic inputs to the algorithm. Data gaps in the 8-day temporal MODIS FPAR/LAI caused by cloudiness were filled with information from accompanying quality-assessment fields (SOM text S2) (10). For daily meteorological data required to drive the algorithm, we used a reanalysis data set from National Center for Environmental Prediction (NCEP) (SOM text S3) (11). A Palmer Drought Severity Index (PDSI) (12) at 0.5° resolution was used as a surrogate of soil moisture (13) to measure environmental water stress by combining information from both evaporation and precipitation (SOM text S4). A lower PDSI generally implies a drier climate.

Global NPP slightly decreased for the past decade by –0.55 Pg C (Fig. 1). Interannual variations of the global NPP were negatively correlated with the global atmospheric CO2 growth rates (correlation coefficient r = –0.89, p < 0.0006) (Fig. 1) (14), suggesting that global terrestrial NPP is a major driver of the interannual CO2 growth rate. Carbon isotopic measurements have indicated that the exchange of CO2 with terrestrial ecosystems is the dominant cause of the CO2 interannual growth rate (15). Though NPP is a part of carbon exchange between the land and atmosphere, the strong correlation may imply that the process of heterotrophic respiration depends ultimately on the substrate supply from NPP (16), as shown by the strong correlation between gross primary production (GPP) and ecosystem respiration derived from eddy-flux towers (17). Although also a major contributor, global fire emissions do not correlate as strongly as NPP with interannual CO2 growth rate (SOM text S5) (18).

Fig. 1

Interannual variations from 2000 through 2009 in anomalies of annual total global terrestrial NPP (green circles) and inverted global atmospheric CO2 annual growth rate [red squares and (14)]. Global average annual total NPP is 53.5 Pg C/yr.

Regionally, negative annual NPP anomalies were mainly caused by large-scale droughts (fig. S7 and table S2). In 2000, droughts reduced NPP in North America and China; in 2002, droughts reduced NPP in North America and Australia; in 2003, drought caused by a major heat wave reduced NPP in Europe (19); in 2005, severe droughts in the Amazon (20), Africa, and Australia greatly reduced both regional and global NPP (Fig. 1 and fig. S7); and from 2007 through 2009, over large parts of Australia, continuous droughts reduced continental NPP.

NPP in the tropics (23.5°S to 23.5°N) explains 93% of variations in the global NPP, of which tropical rainforests explain 61% of global NPP variations (SOM text S8 and table S4). For the three major tropical rainforests (Amazon, Africa and Asia, see fig. S11), Amazon NPP alone explains 66% of the global NPP variations, though it accounts for only 14% of the global total, whereas the other two rainforests have no significant relation with global NPP variations (table S4). The growth of rainforests is mainly constrained by solar radiation due to severe cloudiness (fig. S10) (2, 21), yet our study reveals that, except for rainforests in the islands of southeast Asia, the growth of rainforests in both the Amazon and Africa is also water limited (fig. S10 and table S5). Of the three major rainforests, only Africa had an increasing trend in NPP (0.189 Pg per decade), which was mostly due to decreased vapor pressure deficit (VPD) (figs. S9 and S10). The other two had decreasing trends, with the Amazon at –0.424 Pg C per decade and Asia at –0.562 Pg C per decade (table S4). The decreased NPP in Asian rainforests was mainly caused by reductions in solar radiation, whereas reduction of NPP in the Amazon was mainly caused by increasing air temperature, which greatly increased autotrophic respiration (0.38 Pg C per decade), and also by a slight drying trend for the past decade, especially a severe drought in 2005 (figs. S12 and S14 and table S5).

There has been an ongoing debate about whether there was “green-up” of Amazon rainforests in the intense drought of 2005 (22, 23). Hence, we examined MODIS C5 FPAR/LAI, related climate variables, and the resulting GPP and NPP over Amazon rainforests (SOM text S9). FPAR/LAI were higher in 2005, implying green-up of Amazon during drought (figs. S13 and S14). However, annual calculated NPP greatly decreased because of the combined effects of severe water stress and high-temperature–induced high-autotrophic respiration (figs. S13 and S14). The reduced NPP agrees with reported large decreases in plant growth (20). Though there was no green-up of C5 MODIS enhanced vegetation index (EVI) from the drought, there was no clear “browning” either (23). It is unclear what caused green-up of FPAR/LAI or the lack of browning of EVI; however, the reduction in NPP due to severe drought reveals that radiometric greenness is not necessarily a proxy for NPP in some cases. Therefore, remotely sensed NPP models should account for environment stresses and plant autotrophic respiration.

Figure 2 is the spatial pattern of NPP trends over the past decade: NPP increased over large areas in the Northern Hemisphere (NH) while decreasing in the Southern Hemisphere (SH). Over the NH, 65% of vegetated land area had increased NPP, including large areas of North America, Western Europe, India, China, and the Sahel. Regions with decreased NPP include Eastern Europe, central Asia, and high latitudes of west Asia. In the SH, 70% of vegetated land areas had decreased NPP, including large parts of South America, Africa, and Australia.

Fig. 2

Spatial pattern of terrestrial NPP linear trends from 2000 through 2009 (SOM text S1) (8, 10).

For northern high latitudes (>47.5°N), warming climate lengthens growing seasons, promoting plant growth (24). However, continuing warming may offset these benefits of an earlier spring and decrease carbon sequestration in a dryer summer and warmer autumn (25, 26). Over the past decade, warming temperatures continued increasing NPP, albeit with a concurrent drying trend (figs. S8 and S9 and table S3). For northern mid-latitudes (22.5°N to 47.5°N), where there are large areas of short-rooted grasslands and croplands (fig. S2), water availability is a dominant control on plant growth. NPP had significant correlations with both growing-season total precipitation (r = 0.96, p < 0.0001) and average PDSI (r = 0.76, p < 0.05) (table S3), and a wet trend over the past decade increased NPP (fig. S8 and S9 and table S3). For the Tibet plateau with an average elevation of more than 4500 m, similar to the high latitudes, temperature is the dominant control on vegetation growth, and the recent warming climate increased NPP (Fig. 2 and figs. S9 and S10). For northern low latitudes (0°N to 22.5°N), reduced atmospheric VPD in Africa and India caused increased NPP in the zone (Fig. 2, fig. S9, and table S3). Overall, the NH had an increasing trend in NPP (Figs. 2 and 3, fig. S8, and table S3). Nitrogen deposition and fertilizer use also contributed to increased NPP in the NH (SOM text S5) (27).

Fig. 3

Interannual variations in total NPP (solid lines) and growing-season average PDSI (dotted lines) over the NH (green) and SH (red) (SOM text S1 and S4).

For vegetated land areas in the SH, temperature is not a dominant factor controlling the length of growing season, as indicated by an average of only 7.5 days of snow cover annually, in contrast to 125 days in the NH, as observed by MODIS (fig. S3). Instead, high temperature causes high VPD and evaporation, generally reducing water availability for vegetation growth. Therefore, there is a significant correlation (r = 0.87, p < 0.001) between NPP and PDSI in the SH, though only a weak correlation (r = 0.39, p < 0.27) in the NH (Fig. 3). Despite increased precipitation over large parts of the SH (fig. S9 and table S3), the warming trend induced a much higher evaporative demand, leading to a drying trend (fig. S9 and table S3). As a result, NPP decreased because of warming-associated drying trends (Figs. 2 and 3, and figs. S9 and S10). Though NPP in the SH accounts for 41% of the global NPP, the stronger decreasing NPP trend (–1.83 Pg C per decade) counteracted a weaker increasing trend in the NH (1.28 Pg C per decade), resulting in a slight reduction in the global NPP (Figs. 1 and 3 and table S3).

To understand why there are different responses of NPP to a warming climate for two hemispheres, we examine two basic variables governing water and carbon fluxes: (i) saturation vapor pressure (SVP) and (ii) a maintenance respiration index (MRI) (SOM text S10). VPD is the difference between SVP and air vapor pressure, and SVP is a nonlinear function of air temperature (T). VPD is one major factor controlling surface evapotranspiration (ET) (28); thus, a high T value can increase both VPD and subsequent ET and lead to a dryer environment, eventually reducing NPP. MRI is a power function of T, and a warming climate will increase autotrophic respiration and potentially decrease NPP. Because of the amplification nature of these nonlinear functions to increasing temperature, for different temperature bases, the same increase in air temperature will have much higher increases of SVP and MRI at higher temperature bases than at lower temperatures. Figure 4 shows how SVP and MRI and their derivations respond to air temperatures of the two hemispheres. For vegetated land area, the annual mean T of the SH is more than 12°C higher than that of the NH. An equivalent increase of temperature will consequently have changes in SVP and MRI that are twice as strong in the SH than in the NH. Only during the NH summer (June, July, and August) will there be similar levels of response to warming compared to those of the SH. For the NH, summer is the major growing season; thus, NPP is also susceptible to warming and drought.

Fig. 4

SVP of the air as a function of T and derivative (A). MRI of plants as a function of T and derivative (B). Twelve-month average air temperature and derivatives of SVP and MRI with respect to T for both the NH and SH are also shown (table S6). For MRI, there are two equations, one with Q10 (the difference in respiration rates over a 10°C interval) as a constant value of 2.0, and the other with a temperature-acclimated Q10. Annual mean air temperature is 9.51°C for vegetated NH and 21.56°C for vegetated SH (table S6). See SOM text S10 for details of the equations.

Recent simulations from coupled climate–carbon cycle models have shown that there is a positive feedback between the carbon cycle and the climate system (2931), and also that future biological carbon sinks could eventually level off and subsequently decline to zero (32). Though Fig. 4 only shows MRI responses to temperature, soil respiration largely follows a similar nonlinear temperature response in these models (31). These nonlinear amplification responses of water and carbon processes to warming are the major mechanisms responsible for the positive feedback between the carbon and climate systems.

Though the warming climate during this period continuously increased NPP over areas of high latitude and high elevations, these warming-benefited areas only account for 16% of the global total NPP and 24% of global vegetated land area (SOM text S11). NPP over large land areas of lower latitude and altitude is negatively correlated with temperature (fig. S10A), mostly due to the warming-related increases in water stress and autotrophic respiration, especially for the SH (r = –0.94, p < 0.0001). Globally, interannual NPP is negatively correlated with air temperature over vegetated land (r = –0.64, p < 0.05) (table S3).

Over the past 10 years, large-scale periodic regional droughts and a general drying trend over the SH reduced global terrestrial NPP. Under a changing climate, severe regional droughts have become more frequent, a trend expected to continue for the foreseeable future (13, 33, 34). The warming-associated heat and drought not only decrease NPP, but also may trigger many more ecosystem disturbances (6, 35, 36), releasing carbon to the atmosphere (18, 37). Reduced NPP potentially threatens global food security and future biofuel production and weakens the terrestrial carbon sink. Continuous global monitoring of NPP will be essential to determining whether the reduced NPP over the past 10 years is a decadal variation or a turning point to a declining terrestrial carbon sequestration under changing climate.

Supporting Online Material

SOM Text

Figs. S1 to S14

Tables S1 to S6

References and Notes

References and Notes

  1. We thank the reviewers and L. Jones, C. Cleveland, and Q. Mu for helpful comments. This work would not have been possible without continuous financial support from the NASA Earth Observing System MODIS project (grant NNX08AG87A). Moreover, we gratefully acknowledge all MODIS land product science team members for their efforts in generating high-quality data sets. NCEP_Reanalysis 2 data were provided by the NOAA/OAR/Earth System Research Laboratory Physical Sciences Division, Boulder, Colorado, USA, from their Web site at
View Abstract

Navigate This Article