Sensitivity of the Terrestrial Ecosystem to Precipitation and Temperature Variability over China
SUN Guo-Dong
The State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
Corresponding author: SUN Guo-Dong,sungd@mail.iap.ac.cn
Abstract

In this study, the sensitivities of net primary production (NPP), soil carbon, and vegetation carbon to precipitation and temperature variability over China are discussed using the state-of-the-art Lund-Potsdam-Jena dynamic global vegetation model (LPJ DGVM). The impacts of the sensitivities to precipitation variability and temperature variability on NPP, soil carbon, and vegetation carbon are discussed. It is shown that increasing precipitation variability, representing the frequency of extreme precipitation events, leads to losses in NPP, soil carbon, and vegetation carbon over most of China, especially in North and Northeast China where the dominant plant functional types (i.e., those with the largest simulated areal cover) are grass and boreal needle-leaved forest. The responses of NPP, soil carbon, and vegetation carbon to decreasing precipitation variability are opposite to the responses to increasing precipitation variability. The variations in NPP, soil carbon, and vegetation carbon in response to increasing and decreasing precipitation variability show a nonlinear asymmetry. Increasing precipitation variability results in notable interannual variation of NPP. The sensitivities of NPP, soil carbon, and vegetation carbon to temperature variability, whether negative or positive, meaning frequent hot and cold days, are slight. The present study suggests, based on the LPJ model, that precipitation variability has a more severe impact than temperature variability on NPP, soil carbon, and vegetation carbon.

Keyword: climate variability; net primary production (NPP); soil carbon; vegetation carbon; sensitivity
1 Introduction

Climate variability exhibits distinct differences on various timescales and in different regions according to observational data ( Easterling et al., 2000; Schneider et al., 2004; Zveryaev, 2004; Gu and Adler, 2006; Parker et al., 2007; Efthymiadis et al., 2007). Changes in climate extreme events, such as floods, droughts, and heat waves, may be one of the important causes of variations in climate variability ( Liu et al., 2006; Moberg et al., 2006; Brown et al., 2008). There is also increasing evidence from climate models suggesting further variation in climate variability in the future ( Mishra and Singh, 2009).

As one of the most interesting topics in this field of research, characterizing the response of the terrestrial ecosystem to climate variability has recently begun to receive much attention. Some investigations have indicated that climate variability not only has a severe impact on ecosystem structure ( Claussen et al., 1999; Zeng and Neelin, 2000; Liu et al., 2006), but also leads to changes in soil and vegetation carbon stocks ( Tian et al., 1998; Botta and Foley, 2002). In addition, the level of uncertainty in estimating carbon flux increases as a result of the variation in precipitation variability, as demonstrated, for example, in northern mixed grass prairie ( Mitchell and Csillag, 2001). There are different views about the relationship between climate variability and terrestrial ecosystems ( Fang et al., 2001). For example, Botta and Foley (2002) demonstrated that climate variability plays an important role in the variations of ecosystem structure, soil carbon, and vegetation carbon. Mitchell and Csillag (2001) also pointed out that grassland is sensitive to climate variability.

The response of terrestrial carbon to climate variability over China has been examined previously. The interannual variation of carbon fluxes in China has been demonstrated using different vegetation models ( Cao et al., 2003; Gao and Liu, 2008). Recent estimates have shown that there are notable interannual and decadal variations in China's carbon exchange as a consequence of climate variability ( Mao et al., 2009). However, there has been little discussion about the sensitivity of the terrestrial ecosystem to changing temperature and precipitation variability over China.

In the present study, the sensitivities of net primary production (NPP), soil carbon, and vegetation carbon to temperature and precipitation variability over China are discussed by analyzing results from the state-of-the-art Lund-Potsdam-Jena dynamic global vegetation model ( LPJ DGVM; Sitch et al., 2003). We also show whether the variability in temperature or the variability in precipitation has a more important impact on the terrestrial ecosystem.

2 Model and method
2.1 The model

Some previous studies have demonstrated that terrestrial vegetation models are effective tools to assess the responses of terrestrial ecosystems to climate change and land use ( Gao and Yu, 1998; Friend and White, 2000; Botta and Foley, 2002). In our research, the LPJ model is used to explore the sensitivity of the terrestrial ecosystem over China to climate variability.

The LPJ model, which originates from the biome model family ( Prentice et al., 1992), is process-based model that can represent terrestrial vegetation dynamics and land- atmosphere carbon and water exchanges. The LPJ DGVM explicitly considers photosynthesis, mortality, fire disturbance, and soil heterotrophic respiration, among other processes. The distribution of plant functional types (PFTs) simulated by the LPJ model is represented reasonably well ( Sun, 2009); the carbon flux falls within an acceptable range ( Sun, 2009); and the model has been used previously to study terrestrial ecosystem dynamics ( Beer et al., 2007; Jung et al., 2007; Sitch et al., 2005). For a detailed description and evaluation of the model, see Sitch et al. (2003).

2.2 Input

The simulation of any grid cell with the LPJ model is driven by the input of monthly climatology, including temperature, precipitation, wet day frequency, and cloud cover. The data provided by the Climate Research Unit ( CRU; Mitchell and Jones, 2005) covering the period 1901-1998 (resolution: 0.5°) are used in this study. In addition, eight soil texture classes (coarse, medium, fine, medium-coarse, fine-coarse, fine-medium, fine-medium- coarse, and organic) are used, which are also available globally on a 0.5° grid.

2.3 Modeling protocol

The LPJ model is run from 1901 to 1998 as the reference experiment using the CRU data, preceded by a 1000-year spin-up period, repeatedly using 30-year transient climate data. This achieves a state of equilibrium in terms of a preindustrial stable vegetation structure and carbon pools, which are then employed as the initial conditions for all the sensitivity experiments ( Sitch et al., 2003; Zaehle et al., 2005). To solely examine the sensitivities of NPP, soil carbon, and vegetation carbon to temperature or precipitation variability, we change the variability presented by the standard deviation without changing the mean state.

Several experiments are independently conducted to measure the impact of precipitation or temperature variability on the terrestrial ecosystem. According to the change in climate variability during 1961-2000 based on the monthly precipitation and mean surface air temperature data of 160 terrestrial stations provided by the China Meteorological Administration (Fig. 1), we perform 13 simulations with the same initial state'the reference experiment and 12 disturbance experiments. Table 1 shows the experimental design, in which the precipitation variability and temperature variability during the period of climate variability (1901-1998) are increased or decreased by 20%, 40%, and 60%.

To distinguish the difference between the sensitivities of NPP, soil carbon, and vegetation to precipitation and temperature variability in different regions over China, the whole region is divided into five approximate regions: North China (42-48°N, 107-118°E); Northeast China (42-48°N, 120-128°E); South China (22-33°N, 105- 120°E); Northwest China (35-42°N, 80-97°E); and Southwest China (27-34°N, 80-97°E).

Figure 1 The relative change of variability denoted by standard deviation between 1961-1980 and 1980-2000 using the Chinese Meteorology Administration (a) precipitation and (b) temperature.

Table 1 The variational scope of precipitation and temperature variability denoted by standard deviation.
3 Sensitivities of NPP, soil carbon, and vegeta-tion carbon to climate variability
3.1 The impact of precipitation variability

Figure 2a shows the precipitation variations due to the different experimental perturbations. The precipitation perturbation compared with the reference precipitation is positive for increasing standard deviation and is negative for decreasing standard deviation when the precipitation is larger than the mean state. However, when the precipitation is less than the mean state, the precipitation perturbation is opposite to that when the precipitation is larger than the mean state. An increase of the standard deviation in terms of precipitation means increasing wet and drought events, such as heavy precipitation events. On the contrary, decreasing precipitation variability leads to low interannual precipitation variation.

Figure 2 (a) Annual precipitation for the changing precipitation variability during the period 1901-1998 as the reference period; (b) annual temperature for the changing temperature variability during the period 1901-1998 as the reference period; (c) annual NPP for the changing precipitation variability during the period 1901-1998 as reference period; (d) annual NPP for the changing temperature variability during the period 1901-1998 as the reference period.

To discuss the sensitivity of NPP to precipitation variability, we change the precipitation variability only, without changing the mean state and other input factors. The NPP in all of China shows slight variation. However, the variations in different regions are clearly distinct (Table 2). The average NPP during the experimental period decreases by 10.32-39.55 gC m-2 yr-1 (approxi-mately 2.55%-9.77%) in North China, where grassland is the dominant PFT, and by 4.68-15.31 gC m-2 yr-1(approximately 0.97%-3.16%) in Northeast China, where boreal forest is the dominant PFT, when the increase in precipitation variability rises from 20% to 60% compared to the reference state (Table 2). When the precipitation variability is reduced, the average NPP has an adverse impact compared to when it is increased. The NPP shows a large gain of 11.78-25.09 gC m-2 yr-1 in North China, and of 5.12-9.78 gC m-2 yr-1 in Northeast China, when the reduction in precipitation variability rises from 20% to 60%. The amount that NPP varies as a result of increasing the precipitation variability is greater than that induced by strengthening the decrease in precipitation variability. Therefore, these numerical results demonstrate a nonlinear asymmetry in the variation of NPP with changing precipitation variability. In the other three regions, the variation of NPP is negligible.

To analyze the cause of the variation in NPP, we also show the variation in gross primary production (GPP) and autotrophic respiration ( Ra). In the LPJ model, the difference between GPP and Ra is defined as the NPP. For the five regions, the average GPP and Ra during the exp-erimental period reduce with greater increased precipitation variability, and increase with stronger decreased precipitation variability. These numerical results suggest that increasing precipitation variability restricts photo-synthesis and autotrophic respiration, while decreasing precipitation variability accelerates these processes. The amounts of variation in GPP and Ra in North China, Northeast China, and South China are greater than those in Northwest China and Southwest China. In South China, GPP decreases by 12.20 gC m-2 yr-1 and Ra decreases by 14.90 gC m-2 yr-1, meaning NPP has a slight loss when precipitation variability increases by 60%. In North and Northeast China, the fact that the decrease in the amount of GPP is greater than that of Ra causes the serious loss in NPP. This result is analogous to that reported by Cao et al. (2003).

The sensitivity of NPP variability to precipitation variability is also analyzed. Figure 3b shows that the NPP variability changes little in Northwest China, Southwest China, and South China with increasing precipitation variability, while there is a marked change in North and Northeast China. These results demonstrate that NPP in North and Northeast China is sensitive to precipitation variability. The annual variation of NPP increases with increasing precipitation variability, and decreases with decreasing precipitation variability (Figs. 2c and 3b). The simulation indicates that the growth of grass and boreal forest are highly dependent on water supply in North and Northeast China ( Gao et al., 2000).

Table 2 The variation in net primary production (NPP), gross primary production (GPP), Ra, soil carbon (SoilC), and vegetation carbon (VegeC) for the different regions in China due to different precipitation variability. The reference experiment is for no impact of variation of precipitation variability on the terrestrial ecosystem. PV-60% means that the precipitation variability decreases by 60%. PV+60% means that the precipitation variability increases by 60%.

There is marked variation in soil carbon not only in North China and Northeast China, but also in South China with changing precipitation variability (Table 2). Soil carbon consists of three components: below-ground litter, the short-term carbon pool, and the long-term carbon pool. The below-ground litter and long-term carbon pool components are insensitive to changing precipitation variability; however, the decomposition of litter responds rapidly to changing soil moisture caused by precipitation variability. In addition, the variation in soil carbon shows a nonlinear asymmetry. In North China, the soil carbon increases by 1.78% to 3.04% when the decrease in precipitation variability rises from 20% to 60%, and decreases by 0.72% to 3.17% when the increase in precipitation variability rises from 20% to 60%. The situation is similar in Northeast China. However, in South China, soil carbon decreases by 9.55-21.53 gC m-2 yr-1 when the decrease in precipitation variability rises from 20% to 60%, and increases by 19.32-54.98 gC m-2 yr-1 when the increase in precipitation variability rises from 20% to 60%.

The vegetation carbon is restrained with increasing precipitation variability in the whole region (Table 2). When the increase in precipitation variability rises from 20% to 60%, the vegetation carbon decreases by 8.08% to 25.25%, by 2.68% to 10.07%, and by 0.44% to 0.94% in North China, Northeast China, and South China, respectively. When the decrease in precipitation varia-bility rises from 20% to 60%, the vegetation carbon gradually increases. The variation of vegetation carbon originates from leaf, sapwood, heartwood, and root variations. Because the variation of soil moisture is responsive to changing precipitation variability, the variations of sapwood, heartwood, and root are noticeable. It is clear that nonlinear asymmetries exist in the soil and vegetation carbon.

3.2 Temperature variability

There is a similar variation in temperature perturbation with that in precipitation perturbation when the temperature variability changes (Fig. 2b). The perturbation causes high temperatures when the temperature is greater than the mean state, and causes low temperatures when the temperature is less than the mean state when the temperature variability rises. This means that hot and cold extreme events become more frequent.

Figure 3 (a) Simulated NPP variability during the period 1901-1998 (units: mm); (b) the difference between NPP variability for decreasing precipitation variability (PV-60%) and the NPP variability for the reference state; (c) the difference between NPP variability for increasing precipitation variability (PV+60%) and the NPP variability for the reference state.

In terms of the variation in soil carbon, vegetation carbon, and NPP, there is little sensitivity to the hot and cold events in the whole region; for example, the different impacts of temperature variability and precipitation variability on soil carbon in Northeast China. The relative change of soil carbon is 3.17% when the precipitation variability increases by 60%. However, the relative change of soil carbon is 0.09% when the temperature variability increases by 60%. The soil carbon consists of three components: below-ground litter, the short-term carbon pool, and the long-term carbon pool. The effects of precipitation and temperature on the below-ground litter and long-term carbon pool components are similar. However, the change of precipitation variability accelerates the decomposition of litter into soil carbon for the short-term carbon pool component due to the rapid response of soil moisture to precipitation. The soil temperature responds slowly to temperature variability. The decomposition of litter into soil carbon is also slow. From the results, we find that NPP, soil carbon, and vegetation carbon are more sensitive to precipitation variability than temperature variability.

4 Discussion and summary

In this study, the sensitivities of NPP, soil carbon, and vegetation carbon to precipitation and temperature vari-ability have been discussed. Within the LPJ DGVM, the numerical results show that precipitation variability plays an active role in vegetation dynamics over China. In-creasing precipitation variability restrains NPP, GPP, Ra,soil carbon, and vegetation carbon, especially in North China where the dominant PFT is grassland and in Northeast China where the dominant PFT is boreal forest. Decreasing precipitation variability has the opposite effect on the terrestrial ecosystem. The results also show that there is a nonlinear asymmetry with the changing pre-cipita-tion variability. The variational amounts in carbon flux induced by increasing precipitation variability are greater than those caused by decreasing precipitation variability. However, the response of the terrestrial eco-system to temperature variability is slight. These results indicate that the interannual precipitation variability may be a crucial factor influencing the terrestrial ecosystem in North and Northeast China, where precipitation is scarce (Cao et al., 2003). These results also confirm that terres-trial ecosystems may be quite sensitive to precipitation perturbations ( Lotsch et al., 2003).

The numerical results show that the interannual NPP rises with an increase in precipitation variability. This means that when there are more flood or drought events, the NPP has a severe response. However, despite the precipitation and temperature variability, the dominant PFTs remain the same in China.

In this study, the variation in climate variability is linear. However, climate change is nonlinear, so it is neces- sary to discuss the sensitivity of the terrestrial ecosystem to nonlinear climate change by employing a nonlinear approach. The LPJ model, in which the feedback to changing climate is not considered, is an offline model. It is therefore reasonable to discuss the terrestrial ecosystem through coupling the LPJ DGVM into general circulation models. However, this research is only the initial step for discussing the variation in the terrestrial ecosystem in China. One of the next steps will be to introduce nonlinear variation into the terrestrial ecosystem using a nonlinear method, such as the conditional nonlinear optimal perturbation (CNOP) approach ( Mu et al., 2003).

Acknowledgments. The authors thank the reviewers for their valuable suggestions. Funding was provided by grants from the National Basic Research Program of China (Grant No. 2012CB955202), the National Natural Science Foundation of China (Grant No. 41375111), the LASG Free Exploration Fund, and the LASG State Key Laboratory Special Fund.

Reference
1 Beer C. , W. Lucht, D. Gerten, et al. , 2007: Effects of soil freezing and thawing on vegetation carbon density in Siberia: A modeling analysis with the Lund-Potsdam-Jena Dynamic Global Vegetation Model (LPJ-DGVM), Glob. Biogeochem. Cycles, 21, GB1012, doi: DOI:10.1029/2006GB002760.
2 Botta A. , J. A. Foley, 2002: Effects of climate variability and disturbances on the Amazonian terrestrial ecosystems dynamics, Glob. Biogeochem. Cycles, 16(4), 1070, doi: DOI:10.1029/2000GB-001338.
3 Brown S. J. , J. Caesar, C. A. T. Ferro, 2008: Global changes in extreme daily temperature since 1950, J. Geophys. Res. , 113, D05115, doi: DOI:10.1029/2006JD008091.
4 Cao M. K. , S. D. Prince, K. Li, et al. , 2003: Response of terrestrial carbon uptake to climate interannual variability in China, Glob. Change Biol. , 9, 536-546.
5 Claussen M. , C. Kubatzki, V. Brovkin, et al. , 1999: Simulation of an abrupt change in Saharan vegetation in the Mid-Holocene, Geophys. Res. Lett. , 26(14), 2037-2040.
6 Easterling D. R. , T. Karl, K. Gallo, et al. , 2000: Observed climate variability and change of relevance to the biosphere, J. Geophys. Res. , 105, D15, 20101-20114.
7 Efthymiadis D. , P. D. Jones, K. R. Briffa, et al. , 2007: Influence of large-scale atmospheric circulation on climate variability in the Greater Alpine Region of Europe, J. Geophys. Res. , 112, D12104, doi: DOI:10.1029/2006JD008021.
8 Fang J. Y. , S. L. Piao, Z. Y. Tang, et al. , 2001: Interannual variability in net primary productivity and precipitation, Science, 293, 1723, doi: DOI:10.1126/science.293.5536.1723a.
9 Friend A. D. , A. White, 2000: Evaluation and analysis of a dynamic terrestrial ecosystem model under preindustrial conditions at the global scale, Glob. Biogeochem. Cycles, 14(4), 1173-1190.
10 Gao Z. Q. , J. Y. Liu, 2008: Simulation study of China's net primary production, Chinese Sci. Bull, 53(3), 434-443.
11 Gao Q. , M. Yu, 1998: A model of regional vegetation dynamics and its application to the study of northeast China transect (NECT) responses to global change, Glob. Biogeochem. Cycles, 12(2), 329-344.
12 Gao Q. , M. Yu, X. Yang, 2000: An analysis of sensitivity of terrestrial ecosystems of China to climatic change using spatial simulation, Climatic Change, 47, 373-400.
13 Gu G. , R. F. Adler, 2006: Interannual rainfall variability in the tropical Atlantic region, J. Geophys. Res. , 111, D02106, doi: DOI:10.1029/2005JD005944.
14 Jung M. , M. Vetter, M. Herold, et al. , 2007: Uncertainties of modeling gross primary productivity over Europe: A systematic study on the effects of using different drivers and terrestrial biosphere models, Glob. Biogeochem. Cycles, 21, GB4021, doi: DOI:10.1029/2006GB002915.
15 Liu X. , Z. Y. Yin, X. Shao, et al. , 2006: Temporal trends and variability of daily maximum and minimum, extreme temperature events, and growing season length over the eastern and central Tibetan Plateau during 1961-2003, J. Geophys. Res. , 111, D19109, doi: DOI:10.1029/2005JD006915.
16 Liu Z. , Y. Wang, R. Gallimore, et al. , 2006: On the cause of abrupt vegetation collapse in North Africa during the Holocene: Climate
17 variability vs. vegetation feedback, Geophys. Res. Lett. , 33, L22709, doi: DOI:10.1029/2006GL028062.
18 Lotsch A. , M. A. Friedl, B. T. Anderson, et al. , 2003: Coupled vegetation-precipitation variability observed from satellite and climate records, Geophys. Res. Lett. , 30(14), 1774, doi: DOI:10.1029/2003-GL017506.
19 Mao J. F. , B. Wang, Y. J. Dai, 2009: Sensitivity of the carbon storage of potential vegetation to historical climate variability and CO2 in continental China, Adv. Atmos. Sci. , 26(1), 87-100, doi: DOI:10.1007/s00376-009-0087-z.
20 Mishra A. K. , V. P. Singh, 2009: Analysis of drought severity-area-frequency curves using a general circulation model and scenario uncertainty, J. Geophys. Res. , 114, D06120, doi: DOI:10.1029/2008JD010986.
21 Mitchell T. D. , P. D. Jones, 2005: An improved method of constructing a database of monthly climate observations and associated high-resolution grids, Int. J. Climatol. , 25(6), 693-712.
22 Mitchell S. W. , F. Csillag, 2001: Assessing the stability and uncertainty of predicted vegetation growth under climatic variability: Northern mixed grass prairie, Ecol. Model. , 139(2-3), 101-121.
23 Moberg A. , P. D. Jones, D. Lister, et al. , 2006: Indices for daily temperature and precipitation extremes in Europe analyzed for the period 1901-2000, J. Geophys. Res. , 111, D22106, doi: DOI:10.1029/2006JD007103.
24 Mu M. , W. S. Duan, B. Wang, 2003: Conditional nonlinear optimal perturbation and its applications, Nonlinear Processes Geophys. , 10, 493-501.
25 Parker D. , C. Folland , A. Scaife, et al. , 2007: Decadal to multidecadal variability and the climate change background, J. Geophys. Res. , 112, D18115, doi: DOI:10.1029/2007JD008411.
26 Prentice I. C. , W. Cramer, S. P. Harrisonet, et al. , 1992: A global biome model based on plant physiology and dominance, soil properties and climate, J. Biogeogr. , 19(2), 117-134.
27 Schneider D. P. , E. J. Steig, J. C. Comiso, 2004: Recent climate variability in Antarctica from satellite-derived temperature data, J. Climate, 17(7), 1569-1583.
28 Sitch S. , B. Smith, I. C. Prentice, et al. , 2003: Evaluation of ecosystem dynamics, plant geography and terrestrial carbon cycling in the LPJ Dynamic Vegetation Model, Glob. Change Biol. , 9, 161-185.
29 Sitch S. , V. Brovkin, W. Von Bloh, et al. , 2005: Impacts of future land cover changes on atmspheric CO2 and climate, Glob. Biogeochem. Cycles, 19, GB2013, doi: DOI:10.1029/2004GB002311.
30 Sun G. D. , 2009: Simulation of potential vegetation distribution and estimation of carbon flux in China from 1981 to 1998 with LPJ dynamic global vegetation model, Climatic Environ. Res. (in Chinese), 14(4), 341-351.
31 Tian H. , J. M. Melillo, D. W. Kicklighter, et al. , 1998: Effect of interannual climate variability on carbon storage in Amazonian ecosystems, Nature, 396, 664-667.
32 Zaehle S. , S. Sitch, B. Smith, et al. , 2005: Effects of parameter uncertainties on the modeling of terrestrialbiosphere dynamics, Glob. Biogeochem. Cycles, 19, GB3020, doi: DOI:10.1029/2004-GB002395.
33 Zeng N. , J. D. Neelin, 2000: The role of vegetation-climate interaction and interannual variability in shaping the African Savanna, J. Climate, 13, 2665-2670.
34 Zveryaev I. I. , 2004: Seasonality in precipitation variability over Europe, J. Geophys. Res. , 109, D05103, doi: DOI:10.1029/2003-JD003668.