Changes in the Terrestrial Carbon Cycle of China during the 2010 Drought
Received Date: Apr 30, 2013 / Accepted Date: Jun 13, 2013 / Published Date: Jun 18, 2013
Satellite remote sensing was combined with the NASA-CASA (Carnegie Ames Stanford Approach) carbon cycle simulation model to evaluate the impact of the 2010 drought throughout China. Results indicated that net primary production (NPP) for 2010 declined most notably in the provinces of Xinjiang, Hebei, and Zhejiang, predominantly in cropland and mixed forest areas. Annual NPP in the most drought-impacted areas of southwestern China declined by an average of 13% compared to the previous 10-year average. The greatest decline in total ecosystem carbon gain during 2010 (compared to 2000 to 2009) was detected in the provinces of Yunnan (-28.1 Tg C), Tibet (-18.8 Tg C), and Guizhou, (-6.3 Tg C) predominantly in cropland, grassland, and mixed forest areas. Despite these drought impacts over the southwest and northern regions, the CASA model estimated that China gained total ecosystem carbon in 2010 at double the mean annual rate (computed for the period 2000 to 2009). Cropland and forest areas of Heilongjiang and Inner Mongolia remained productive in 2010 and these regions were estimated, along with several other provinces, to offset and exceed declines in plant production and total ecosystem carbon losses in the most drought-impacted areas of Tibet, Guizhou, Yunnan, Xinjiang, and Zhejiang provinces.
Keywords: China; Drought; carbon cycle; MODIS
China experienced a severe drought during 2009 and 2010, particularly in southwestern and northern provinces, where the accumulated precipitation from May 2009 to April 2010 was 25% below normal . This included the longest recorded period of zero precipitation during a winter season (October–February) in the past 50 years in some provinces, and the lowest percentage rainfall anomaly over the same season since 1880. The most extreme drought conditions were experienced by the end of March 2010 in the southern provincial areas of Tibet, Qinghai, Yunnan, Sichuan, Guizhou, and Guangxi, and moderate drought conditions extended into much of them northern provincial areas of Shanxi, Hebei, and Inner Mongolia by the end of July 2010 (Figure 1) .
One of the most important impacts of severe drought on a region is the loss of plant production due to extreme heat and soil moisture stress . Such losses can lead to crop failures, reduction in livestock herds, hardening of soils, and wildfires, all in association with depletion of river flow and groundwater aquifers. Consequently, changes in plant production and carbon cycling between baseline years (with more normal weather conditions) and years of extreme drought can be a valuable assessment indicator.
Net primary production (NPP) by vegetation provides the chemical energy that drives most biotic processes on Earth. NPP represents much of the carbon that is consumed by microbes and animals. Climate controls on plant production are an issue of central relevance to society, mainly because of concerns about the extent to which NPP in managed ecosystems can provide human populations with adequate forage and fiber products. In areas where NPP is not harvested for human use, the carbon in dead plant biomass returns to the atmosphere mainly as CO2 from decomposition by microbial respiration (Rh). The difference between NPP and Rh on the same time scale is called net ecosystem production (NEP), which represents the net sink or source of atmospheric CO2 from terrestrial ecosystems and the most complete indicator of carbon exchange in land areas.
In this study, we summarized results from the CASA (Carnegie- Ames-Stanford Approach) model for estimated NPP and NEP across all of China during 2010, and compared these fluxes to 10 previous baseline years of NPP and NEP estimates. Direct input of satellite vegetation index “greenness” data from the NASA moderate resolution imaging spectroradiometer (MODIS) satellite sensor into the simulation model with regional climate records was used to estimate spatial variability in monthly NPP and NEP at a ground resolution of 8-km [4,5]. These CASA model results have been designed to be spatially detailed enough to support drought impact assessments in different vegetation and land use types, e.g., grasslands, woodlands, forests, and wetland areas .
Monthly NPP of vegetation from the CASA model was predicted using the relationship between greenness reflectance properties and the fraction of absorption of photosynthetically active radiation (fPAR), assuming that net conversion efficiencies of PAR to plant carbon can be approximated for different ecosystems or are nearly constant across all ecosystems [7-9]. For this study, we used MODIS collection 5 of the Enhanced Vegetation Index  as model inputs for PAR interception, aggregated for regional assessments to an 8-km spatial resolution. Monthly 0.05° data (MOD13C2 version 005) from the USGS Land Processes data center (LP DAAC) was multiplied by 0.0001 to scale the data to the standard MODIS EVI value range. The MOD13C2 product provided the greatest assurance of high-quality, cloud-free EVI inputs to our carbon cycle model. Gap-filling and temporal smoothing procedures are commonly accepted methods for unifying multi-year satellite EVI data sets in cloudy regions . In order to calculate annual CASA NPP, these procedures must be applied to fill cloud- (or aerosol-) contaminated observations, since plant photosynthesis continues during daylight hours regardless of cloud cover .
As documented in Potter , monthly production of plant biomass is estimated as a product of time-varying surface solar irradiance, Sr, and EVI (for fPAR) from the MODIS sensor, plus a constant light utilization efficiency term (emax) that is modified by time-varying stress scalar terms for temperature (T) and moisture (W) effects (Equation 1, NPP in units of g C m-2 mo-1).
NPP = Sr EVI emax T W (1)
The CASA emax term was set uniformly at 0.55 g C MJ-1 PAR, an approach that derives from calibration of predicted annual NPP to previous field estimates . This model setting has been successfully validated globally (with a high level of statistical significance, i.e., R2>0.9) by comparing predicted annual NPP to more than 1900 field measurements of NPP , and against numerous Fluxnet eddy covariance tower site measurements of NPP from 2000-2007 . Gridded monthly climate inputs for these CASA runs were from National Center for Environmental Prediction (NCEP) reanalysis products . Land cover settings were aggregated from the MODIS global 1-km product .
The T stress scalar in Equation 1 was computed with reference to derivation of optimal temperatures (Topt) for plant production. The Topt setting will vary by latitude and longitude, ranging from just above 0°C in the Arctic to the middle thirties in low latitude deserts. The W stress scalar was estimated from monthly water deficits, based on a comparison of moisture supply (precipitation and stored soil water) to potential evapotranspiration (PET) demand using the methods described by Potter . The soil model design includes three-layer heat and moisture content computations: surface organic matter, topsoil (0.3 m), and subsoil to rooting depth (1 to 10 m). Maximum rooting depth for cropland and grassland cover types was set at 1 m, whereas non-tropical forest was set at 2 m and tropical forest was set at 10 m . These layers can differ in soil texture, moisture holding capacity, and carbon-nitrogen dynamics. Water balance in the soil was modeled as the difference between precipitation or volumetric percolation inputs, monthly estimates of PET, and the drainage output for each layer.
Inputs from rainfall can recharge the soil layers to field capacity. Excess water percolates through to lower layers and may eventually leave the system as seepage and runoff. The CASA model is designed to couple daily and seasonal patterns in soil nutrient mineralization and soil heterotropic respiration (Rh) of CO2 from soils. The CASA soil carbon model uses a set of compartmental difference equations . First-order decay equations simulate exchanges of decomposing plant residue (metabolic and structural fractions) at the soil surface. The model also simulates surface soil organic matter (SOM) fractions that presumably vary in age and chemical composition. Turnover of active (microbial biomass and labile substrates), slow (chemically protected), and passive (physically protected) fractions of the SOM are represented. Monthly and annual NEP was computed as NPP minus Rh fluxes, excluding the effects of small-scale fires and other localized disturbances or vegetation regrowth patterns on carbon fluxes. Global CASA NEP predictions have been validated previously against numerous Fluxnet eddy covariance tower site measurements of NEP from 2000-2007 .
Annual NPP on a unit area basis (g C m-2) was estimated to be highest over the period 2000-2010 within evergreen forest, savanna and cropland vegetation classes of China (Figure 2). The highest annual NPP on a unit area basis was detected in these vegetation classes within the provinces of Guangxi, Guizhou, Guangdong, Tibet, Hainan, and Yunnan (Table 1). The annual NPP ranking of provinces changes somewhat when total land area is factored into the assessment, with the leading provinces ordered (with highest total NPP first) as Yunnan, Sichuan, Heilongjiang, Inner Mongolia, Tibet, and Guangxi (Figure 3a). At the national level, the vegetation classes ranked with highest total 2010 NPP were (in order of highest first) were cropland, evergreen forest, grassland, and savanna. Total NPP for 2010 was estimated at 2.9 Pg [1015 g] C for the entire country.
|Evergreen Broadleaf forest||Evergreen Needleleaf forest||Deciduous Needleleaf forest||Mixed forest||savana||Grassland||Barren||Cropland||Deciduous Broadleaf forest||Open Shrubland|
Table 1: CASA average NPP flux (grams Cm-2 yr-1) in 2012 for chinese provinces by vegetation cover class. Shaded cells indicate the highest mean NPP flux values for the country.
The greatest decline in total annual NPP during 2010 (compared to the mean total NPP from 2000 to 2009) was detected in the provinces of Xinjiang (-4.6 Tg [1012 g] C) and Zhejiang (-3.1 Tg C), predominantly in cropland and mixed forest areas (Figure 3b). Cropland NPP declined similarly in 2010 within the provinces of Sichuan, Hebei, Shandong, Tianjin, Jiangxi, Jiangsu, and Shanghai (between -0.1 and -1.9 Tg C). These declines in total annual NPP were offset however by increases in total 2010 NPP (compared to the mean total NPP from 2000 to 2009) in cropland, grassland, and mixed forest vegetation classes in the other provinces of the country. Notable increases in total 2010 NPP were detected in cropland areas of Heilongjiang, Inner Mongolia, and Jilin, and in grassland cover areas of Qinghai and Tibet. The net increase in total NPP for 2010 was estimated at 132 Tg C for the entire country, relative to the mean total NPP from 2000 to 2009.
The greatest decline in total annual NEP during 2010 (compared to the mean total NEP from 2000 to 2009) was detected in the provinces of Yunnan (-28.1 Tg C), Tibet (-18.8 Tg C), and Guizhou, (-6.3 Tg C) predominantly in cropland, grassland, and mixed forest areas (Figure 3c). These declines in total annual NEP were offset however by increases in total 2010 NEP (compared to the mean total NEP from 2000 to 2009) of cropland areas in the other provinces of the country (Figure 2) (Zhao and Running, 2010). Modest increases in total 2010 NEP of between +2 and + 9 Tg C were detected in cropland areas of Heilongjiang, Inner Mongolia, Jilin, Henan, Liaoning, Shaanxi, and Shanxi. Mixed forest areas of Heilongjiang and Inner Mongolia were detected with similar increases of total NEP for 2010. China’s total NEP flux for 2010 was estimated at +12.6 Tg C for the entire country (Table 2), a net increase of +6.3 Tg C, compared to the mean total NEP from 2000 to 2009 represented in Figure 2.
|Evergreen Broadleaf Forest||Evergreen Needle leaf forest||Decidious Needle leaf forest||Mixed Forest||savanna||Grasslands||Barren||Cropland||Decidious Broadleaf forest||Open Shrubland||SUM|
Table 2: CASA total NEp flyx (grams C yr-1)in 2012 for chinise proviences by vegetation class. Net source fluxes of C to the atmosphere are represented as negative values, where net sink fluxes of C in to ecosystems from the atmisphere are represented in positive values.
Examples of CASA monthly NEP flux predictions for selected provinces with among the highest areas of forest, grassland, and cropland cover showed that drought-impacted ecosystems generally had lower estimated net carbon gains from January to June 2010 than did lesser drought-impacted ecosystems (Figure 4). For instance, average forest NEP fluxes in eilongjiang far exceeded (in the positive direction) those in Yunnan from April to July.
Average cropland NEP fluxes in Sichuan exceeded those in Guizhou from March to May.
The drought of 2010 was the most severe during the first half of the year and would have impacted vegetation growth mainly during the months of March to June, as temperatures were becoming warm enough to promote green-up of forest and crop cover types. Despite severe drought over several provinces of the southwest and northern regions, our CASA model estimated that China gained ecosystem carbon at double the mean annual rate in 2010 compared to the past 10 years. Cropland and forest areas of Heilongjiang and Inner Mongolia remained productive in 2010 and these ecosystems were estimated (along with several other provinces listed in Table 2) to offset and exceed declines in NPP and total ecosystem carbon losses in the most drought-impacted areas of Tibet, Guizhou, Yunnan, Xinjiang, and Zhejiang provinces.
Changes in terrestrial NPP from the 2010 drought in the southern provinces of Yunnan, Guizhou, Guangxi, and Sichuan were reviewed by Zhang et al.  based on the global MODIS NPP algorithm (hereafter abbreviated MNA; Zhao and Running, 2010). In contrast to CASA model results, this MNA analysis suggested that annual NPP decreased by over more than 60% of this southern region at between 100 to 200 g C m-2yr-1 lower than the mean annual NPP over the period 2000 to 2009. Such MNA summaries of NPP change are unreliable, however, due to limitations of the MODIS processing algorithms applied by Zhao and Running (2010), which were evaluated by Medlyn  and Samanta et al.  to include generally weak correlations with field observations and extreme sensitivity to increasing air temperatures that can unrealistically increase autotrophic respiration costs and artificially drive down the MNA NPP predictions during drought periods. Our CASA model results suggested instead that 2010 annual NPP in the most drought-impacted areas of southwestern China only declined by an average of between 20 to 50 g C m-2yr-1 (equivalent to 13% CASA mean annual flux, compared to the period 2000 to 2009) .
For national and regional carbon cycle assessments, NEP (rather than NPP) should be considered the most complete measure to apply [13,16]. Because it takes into account the losses of CO2 from decaying organic matter and soils, annual NEP uniquely represents the net sink or source of atmospheric CO2 from terrestrial ecosystems .
Based on the results from this study of China in 2010, we hypothesize that unusually arid periods experienced during severe droughts affect (diminish) plant production to a significantly greater degree than soil respiration fluxes of CO2 from ecosystems.
This work was supported by funding from NASA Ames Research Center.
- Barriopedro D, Gouveia C, Trigo R, Wang l (2012) The 2009/10 drought in China: Possible causes and impacts on vegetation. J Hydrometeorology 13: 1251-1267.
- National Climatic Data Center (2010) State of the Climate, Global Analysis, Annual 2010, National Oceanic and Atmospheric Administration.
- National Research Council (2010) Monitoring Climate Change Impacts: Metrics at the Intersection of the Human and Earth Systems. Committee on Indicators for Understanding Global Climate Change. Washington DC, USA.
- Potter CS, Randerson JT, Field CB, Matson PA, Vitousek PM, et al. (1993) Terrestrial ecosystem production: A process model based on global satellite and surface data, Global Biogeochem, Cycles 7: 811-841.
- Potter C, Klooster S, Myneni R, Genovese V, Tan P, Kumar V (2003) Continental scale comparisons of terrestrial carbon sinks estimated from satellite data and ecosystem modeling 1982-98. Global Planet Change 39: 201-213.
- Zhao M, Running S W (2010) Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. 329: 940-943.
- Goetz SJ, Prince SD, (1998) Variability in carbon exchange and light utilization and net primary production in boreal forest stands. Can J Forest Res 28: 375-389.
- Samanta A, Costa MH, Nunes EL, Vieira SA, Xu L, et al. (2011) Comment on “Drought- Induced reduction in Global Terrestrial Net Primary Production from 2000 Through 2009”. 333:1093.
- Zhang LJ, Xiao J, Li K, Wang L, Lei, et al. (2012) The 2010 spring drought reduced primary productivity in southwestern China, Environ Res Lett 7.
- Huete AR, Didan K, Shimabukuro YE, Ratana P, Saleska SP, et al. (2006) Amazon rainforests green-up with sunlight in dry season. Geophys Res Lett 33.
- Jönsson P, Eklundh LL (2004) TIMESAT: A program for analyzing time-series of satellite sensor data. Com and Geo 30: 833-845.
- Potter CS (1999) Terrestrial biomass and the effects of deforestation on the global carbon cycle. Bio Science. 49: 769-778.
- Potter C, Klooster S, Genovese V (2012) Net primary production of terrestrial ecosystems from 2000 to 2009. Clim Cha.
- Kistler RR, Collins W, Saha S, White G, Woollen J (2001) The NCEP-NCAR 50-year reanalysis: monthly means CD-ROM and documentation, Bull. Am Meteorol Soc 82: 247-68.
- Potter CS, Klooster SA, Brooks V (1999) Interannual variability in terrestrial net primary production: Exploration of trends and controls on regional to global scales. Ecos 11: 36-48.
- Potter C, Klooster S, Genovese V, Hiatt C, Boriah S, et al. (2012) Terrestrial ecosystem carbon fluxes predicted from MODIS satellite data and large-scale disturbance modeling. Int J Geosci.
- Medlyn BE (2011) Comment on “Drought-Induced reduction in Global Terrestrial Net Primary Production from 2000 Through 2009. Sci 333: 1093.
- Sellers PJ, Tucker CJ, Collatz GJ, Los SO, Justice CO, et al. (1994) A global 1x1 NDVI data set for climate studies. Part 2: the generation of global fields of terrestrial biophysical parameters from the NDVI. Int J Remote Sen 15: 3519-3545.
- Running SW, Nemani RR, Heinsch FA, Zhao M, Reeves M, et al. (2004) A continuous satellite-derived measure of global terrestrial primary production. BioScience 54: 547-560.
- Nemani RR, Keeling CD, Hashimoto H, Jolly WM, Piper SC, et al. (2003) Climate driven increases in global terrestrial net primary production from 1982 to 1999. Science 300: 1560-1563.
Citation: Potter C, Li S, Klooster S, Genovese V (2013) Changes in the Terrestrial Carbon Cycle of China during the 2010 Drought. J Earth Sci Clim Change 4: 141. Doi: 10.4172/2157-7617.1000141
Copyright: ©2013 Potter C, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Select your language of interest to view the total content in your interested language
Share This Article
- Total views: 12140
- [From(publication date): 7-2013 - Jul 17, 2019]
- Breakdown by view type
- HTML page views: 8333
- PDF downloads: 3807