Next Article in Journal
Regional Climate Impacts of Irrigation in Northern Italy Using a High Resolution Model
Previous Article in Journal
Remote Sensing of Atmospheric Methane with IR OPO Lidar System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Projection of Net Primary Productivity under Global Warming Scenarios of 1.5 °C and 2.0 °C in Northern China Sandy Areas

1
State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
School of Geographical Sciences, Nanjing University of Information Science and Technology, Nanjing 210044, China
4
School of Geographic Sciences, Xinyang Normal University, Xinyang 464000, China
5
School of Computer and Information Engineering, Beijing Technology and Business University, Beijing 100048, China
*
Author to whom correspondence should be addressed.
Atmosphere 2020, 11(1), 71; https://doi.org/10.3390/atmos11010071
Submission received: 11 December 2019 / Revised: 27 December 2019 / Accepted: 4 January 2020 / Published: 6 January 2020
(This article belongs to the Section Biosphere/Hydrosphere/Land–Atmosphere Interactions)

Abstract

:
Empirical evidence suggests that variations in climate affect the net primary productivity (NPP) across sandy areas over time. However, little is known about the relative impacts of climate change on NPP with global warming of 1.5 and 2.0 °C (GW_1.5 °C_2.0 °C) relative to pre-industrial levels. Here, we used a new set of climate simulations from four Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP 2b) datasets, modified the Carnegie-Ames-Stanford approach (CASA) model and assessed the spatio-temporal variation in NPP in sandy areas of northern China (SAONC). Compared with the reference period (RP, 1986–2005), the NPP variation under four emission scenarios showed clear rising trends and increased most significantly under RCP8.5 with an annual average increase of 2.34 g C/m2. The estimated annual NPP under global warming of 1.5 °C (GW_1.5 °C) increased by 14.17, 10.72, 8.57, and 26.68% in different emission scenarios, and under global warming of 2.0 °C (GW_2.0 °C) it increased by 20.87, 24.01, 29.31, and 39.94%, respectively. In terms of seasonal change, the NPP value under the four emission scenarios changed most significantly in the summer relative to RP, exhibiting a growth of 16.48%. Temperature changes (p > 0.614) had a greater impact on NPP growth than precipitation (p > 0.017), but solar radiation showed a certain negative impact in the middle- and low-latitude regions. NPP showed an increasing trend that changed from the southeast to the central and western regions at GW_1.5 to GW_2.0 °C. NPP was consistent with the spatial change in climate factors and had a promoting role in high latitudes in SAONC, but it was characterized by a certain inhibitory effect at middle and low latitudes in SAONC. The uncertainty of NPP under the four models ranged from 16.29 to 26.52%. Our findings suggest that the impact of GW_1.5 °C is relatively high compared with the current conditions, whereas GW_2.0 °C implies significantly lower projected NPP growth in all areas.

1. Introduction

Net primary productivity (NPP), which indicates the degree of accumulation of atmospheric CO2 in terrestrial ecosystems, has an important role in global climate change [1,2]. In addition, NPP is an important indicator of terrestrial carbon (C) and modulator of ecological processes [3]. Recently, many studies have focused on terrestrial NPP at regional to global scales as well as the driving factors of terrestrial NPP [4,5,6]. Previous works have estimated NPP based on empirical evidence to further study the terrestrial ecosystem evolution and responses to global change [7,8]. Sandy areas are important within the global C cycle [9,10,11], and they are believed to represent a large C pool containing the vast majority of C that is missing from terrestrial ecosystems [12,13,14]. Thus, it is important to understand the amount and variation of NPP in sandy areas [15,16,17].
However, because of the sparse vegetation and relatively low biological productivity in sandy areas, NPP estimates in these areas over a large scale are difficult to perform using the method of direct harvest [18]. Currently, three types of models are predominantly used to simulate NPP: models based on climate [19], light-use efficiency models [3], and mechanistic ecological process models [20]. Although the parameters of climate-based models are few and easily obtained, simulation results can reflect the zonal distribution of vegetation NPP; however, because climate-based models do not involve biological factors that affect changes in NPP, the estimation errors are relatively large [21,22,23,24]. Advances in remote sensing have resulted in the increased development of models that use satellite data to simulate the NPP of terrestrial ecosystems [25,26,27]. The Carnegie-Ames-Stanford approach (CASA) is a light-use efficiency model that simulates NPP based on its relationships with the characteristics of vegetation and environmental parameters or indicators, including solar radiation, temperature, and precipitation [3,28]. Several researchers have successfully implemented the CASA model to simulate NPP over North and South America, Australia, Eurasia, and Africa at various temporal and spatial scales [22,23,28]. Many studies based on the CASA model have estimated the distribution of China’s terrestrial NPP and the responses of NPP to global climate change [23,27,29,30,31]. However, fewer studies have applied this model to sandy areas with sparse vegetation distributions and low precipitation, especially to estimate the potential variation in NPP in sandy ecosystems when global warming is stabilized at 1.5 or 2.0 °C.
Climate change will have extremely adverse impacts on ecosystems and humans [32]. Recently, signatories to the Paris Agreement agreed to hold “the increase in the global average temperature to well below 2.0 °C above preindustrial levels and to adopt measures to limit further increases in temperature to no more than 1.5 °C above preindustrial levels” [33]. The lowered median Mediterranean region water availability resulting from 2.0 °C warming was predicted to be almost twice that resulting from 1.5 °C warming [34,35]. However, regional impact assessments targeting global climate policy (1.5 and 2.0 °C) remain incomplete [36,37,38]. In China, sandy land dominates the changes occurring to land use/land use cover in arid and semiarid regions, which account for 18.12% of China’s land area [39,40]. Climate change directly affects NPP variations [41,42], which in turn affect the extent of desertification in northern China. Therefore, to predict and evaluate the impact of future global C cycles on human social development and formulate scientific climate policies, investigations of the total and spatial changes of NPP in sandy areas under future climatic changes are of great importance.
The Coupled Model Inter-comparison Project Phase 5 (CMIP5), which was facilitated and put into practice by the World Climate Research Programme (WCRP), contributed a key framework for better understanding the mechanism of a changing climate and predicting future scenarios in sandy areas [43]. Previous studies have shown that CMIP5 is to a certain extent able to simulate changes in climate over regional and global scales [44,45,46,47] and is also able to predict the spatial-temporal distribution of arid regions worldwide. However, certain uncertainties in drought simulation on some regional scales remain, such as the central region of North America and the western regions of the Amazon and Africa [35,48,49,50,51]. The Inter-Sectoral Impact Model Inter-comparison Project (ISIMIP) was initiated in 2010 and is currently in its third phase. In this study, we selected ISIMIP2b datasets with better simulation effects in the typical global semi-arid and arid regions [33], simulated the variations of NPP in the northern sandy areas of China well and analyzed the future NPP variations under the different typical concentration paths (RCPS) of ISIMIP2b. First, we focused on comparing the differences of NPP in sandy ecosystems between global increases in temperature of 1.5 and 2.0 °C and that of the reference year (1986–2005). Additionally, we analyzed the effects of the main meteorological factors (precipitation, air temperature, and solar radiation) on NPP variations in future warming scenarios. In addition, the response mechanism of the desert ecosystem NPP to future climate change, which was proposed as “the target of temperature control” in the Paris Agreement, is further discussed to provide a theoretical basis for global C cycle research and desertification controlling.
Consequently, our objectives are to (1) document the temporal and spatial variation of NPP in sandy areas of northern China by modified model, (2) determine the climate change impacts on NPP formation when stabilizing global temperatures increases at 1.5 and 2.0 °C, and (3) understand the large gaps remaining in the study of the regional impacts of temperature increases of 1.5 and 2.0 °C (hereafter referred to as the warming scenarios) by analyzing quantifiable sources of uncertainty.

2. Data and Methods

2.1. Study Area

Northern China’s sandy areas are located from are located from 75°–125° E, 35°–50° N (Figure 1) and form a belt of arcuate desert stretching from the Tarim Basin’s western edge to the eastern Songnen Plain and extending across northwestern China’s warm temperate semiarid, arid, and semi-humid continental climatic zones [40,43]. This zone extends 4500 km from the west to the east and 600 km from the south to the north [44,45], effectively covering 28% of the land area of China while crossing several climate zones, namely, extreme arid, arid, semiarid, and semi-humid climatic areas. The desert zone is one of the eight largest areas of desert globally and a major global source of sandstorms [46,47,48]. The range of precipitation is 30 mm in the west to 450 mm in the east.
Northern China’s sandy areas consist of several regions, including a transitional zone of desert wilderness; the deserts of Guerbantunggurt, Taklimakan, Qaidam Kumtag, Badain Jaran, Tengger, Ulanbuh, and Kubuqi; and the sandy lands of Mu Us, Otindag, Horqin, Hulun Buir, and Qinghai Gonghe. The SAONC are among the most densely populated areas within the semiarid and arid climatic zones, which are characterized as highly sensitive and vulnerable [49]. Despite substantial reductions in desertification land over the past decade, the net increase has remained high compared with that in the 1950s because of the two regions of the Tengger Desert [50].

2.2. Observed Data and Vegetation Data

Observed meteorological and flux data were collected from five flux sites (Figure 1) administered by the Chinese FLUX Observation and Research Network (ChinaFLUX) (http://www.Chinaflux.org) and the China Ecosystem Research Network (CERN) (http://www.cern.ac.cn), including the Atmospheric Environment Monitoring Station in the Taklimakan Desert (http://www.idm.cn/main.asp). We used the open-path eddy-covariance system to measure the flux [52] at a 10-Hz frequency. Fluxes of CO2 and H2O were calculated and recorded every 30 min, and then traditional data quality control measures were applied [53], including a canopy storage calculation, three-dimensional (3D) rotation [54], the Webb-Pearman-Leuning (WPL) correction [55], and the removal of spurious data.
The soil mechanical composition in northern China corresponds well with the bioclimatic zone distribution [56]. Large differences in soil, vegetation, and climate occur across the sandy regions depending on the bioclimatic zone in which they occur [40,57]. Therefore, we chose five sand flux sites representative of different bioclimatic zones and indicate the distribution of NPP across the sandy regions of northern China (Table 1).
We used the normalized difference vegetation index (NDVI) to trace the vegetation distribution and dynamics. The GIMMS (Global Inventory Modeling and Mapping Studies) 15-day (1982–2015) composite NDVI3g dataset, which was applied in the present study, has been shown to have higher accuracy in relation to the GIMMS NDVI for phonological change and vegetation activity monitoring [61], and it has a spatial resolution of 8 km (https://glam1.gsfc.nasa.gov). We also adjusted the NDVI data for nonvegetative factors to represent the interannual variations in SAONC. The future NDVI data (2016–2100) are from http://data.ess.tsinghua.edu.cn; these data are mainly used to evaluate the global land cover under different emission scenarios, and they have been applied in many fields. For vegetation type data, we adopted the 1:100,000 China vegetation type maps, which include statistical information from the Resources and Environment Science Data Center of the Chinese Academy of Sciences of China (http://www.resdc.cn). The type code for major vegetation was regarded as the benchmark for merging and processing [62].

2.3. Climate Projections under the Warming Scenarios

Four models were selected to generate the daily data considering a single run per scenario from ISI-MIP 2b. Table 2 provides the information on each of the four models, and additional information can be obtained at https://www.isimip.org. The simulations included several variables, such as the minimum and maximum air temperatures, wind speed, relative humidity, atmospheric pressure, precipitation, shortwave and longwave radiation, snowfall flux, and specific humidity. The reference period (RP) of 1986 to 2005 was used to calculate the historical NPP. In addition, the projected changes in NPP from 2006 to 2100 were investigated under the RCP2.6, RCP4.5, RCP6.0, and RCP8.5 scenarios.
The UNFCCC defines one year of rising temperature as an increase in the global average temperature to values exceeding certain preindustrial temperatures [63], however, disagreement over the specific definition of rising temperature remains [64], and each definition considers multiple variables, including post-stabilization, peak and transition temperatures. Compared with the temperature increase to above preindustrial levels calculated by the UNFCCC, the IPCC Fifth Assessment Report (AR5) (i.e., the first and third working group reports) assessed future temperature changes relative to RP (1986–2005). These temperature changes can therefore be superimposed onto a temperature increase of 0.61 °C relative to the pre-industrialization level based on observations by HadCRUT4 [65]. This computing method has the advantage of consistent changes in the modern climate relative to the preindustrial. Therefore, a large number of studies have adopted this method, which has resulted in high comparability among studies. By adopting the ISI-MIP 2b approach, Warszawski et al. [66] found that the global temperature increases by 0.89 °C in the RCP2.6 scenario relative to RP with global warming of 1.5 °C (GW_1.5 °C) above the preindustrial level but increases by 1.39 °C with global warming of 2.0 °C (GW_2.0 °C) above the preindustrial level in the RCP4.5 scenario. Thus, GW_1.5 °C is likely to occur during the period 2020–2039 under the 2.6 scenario, and GW_2.0 °C is likely to occur during the period 2040–2059 under the RCP4.5 scenario.

2.4. CASA Model Overview

Potter et al. [3] and Field et al. [28] first developed the CASA model based on the light-use efficiency (LUE) model (Figure 2). The version modified by Zhu et al. [21] was used in the present study. The CASA model incorporates the product of photosynthetic available radiation absorbed by green vegetation (APAR) and the actual LUE (ε):
NPP x , t = A P A R x , t × ε x , t
where NPP x , t is the vegetation NPP in the geographic coordinate system of a given location x and time t.
APAR can be described as follows:
A P A R x , t = S O L x , t × F P A R x , t × 0.5
where S O L x , t represents the total solar radiation at pixel x in month t, F P A R x , t represents the fraction of the incident photosynthetic active radiation absorbed by vegetation, and 0.5 represents 50% of the incoming solar radiation within photosynthetic active radiation (with a wavelength range of 0.38–0.71 μm).
F P A R x , t is based on the relationship between the FPAR and NDVI as well as the simple ratio (SR).
  F P A R x , t = N D V I x , t N D V I m i n N D V I m a x N D V I m i n × F P A R m a x F P A R m i n + F P A R m i n
where N D V I m a x and N D V I m i n are the maximum and minimum values of the NDVI, respectively, and F P A R m a x and F P A R m i n are defined as constants of 0.001 and 0.95, respectively.
Meanwhile, a linear function between FPAR and SR also exists:
  F P A R x , t = S R x , t S R m i n S R m a x S R m i n × F P A R m a x F P A R m i n + F P A R m i n
where S R m a x and S R m i n represent the 95 and 5% values of NDVI, respectively, for different vegetation types. S R x , t is defined as a linear function of the NDVI:
  S R x , t = 1 + N D V I x , t 1 N D V I x , t
The light-use efficiency algorithm can be expressed as follows:
ε x , t = T ε 1 x , t × T ε 2 x , t × W ε x , t × ε m a x
where T ε 1 x , t denotes the limitations of extreme low and high temperatures on LUE, and T ε 2 x , t denotes the decreasing trend of LUE when the environmental temperature deviates from the optimum temperature T o p t x to extreme low and high temperatures.
T ε 1 x , t = 0.8 + 0.2 × T o p t x 0.0005 × T o p t x 2
T ε 2 x , t = 1.184 1 + e 0.2 × T o p t x 10 T x , t × 1 1 + e 0.3 × T o p t x 10 + T x , t
W ε x , t represents the monthly water deficit calculated by the comparing actual (EET) to potential (PET) evapotranspiration data.
W ε x , t = 0.5 + 0.5 × E E T x , t P E T x , t
The EET and PET values are calculated using the method of Yu et al. [67]. ε m a x is the LUE maximum of vegetation under optimal environmental conditions [68], and in the traditional CASA model the ε m a x value was set to 0.389 g C·MJ−1. We refer to the results of Piao et al. [23] and Zhu et al. [21] to define the ε m a x values of typical vegetation types in China.
A number of factors, including the observed data and simulation technology, introduce a certain amount of uncertainty into the simulation result [69,70]. To further quantify the uncertainty of the simulation results, a local sensitivity analysis method was used to determine the uncertainty of the simulation results [66]. The method is to calculate the rate of variation (VR + 10% and VR − 10%) of the simulation results with the chosen parameter ( R U N P ) of the CASA model ±10% (±2 °C) without changing the reference ( R U N R ). Finally, the maximum value of the two represents the influence of this parameter on the simulation results.
VR ± 10 % = R U N R R U N P ± 10 % R U N R × 100   %
RU = 90   % P u 90   % P l P m
where the maximum variation (%) between the VR values was used as a representation of the effect of the initial condition or parameter on the annual NPP. The initial conditions and parameters were then sorted from largest to smallest based on their effect on NPP. All parameters and initial conditions showing effects of >10% (2 °C) on annual NPP were retained for the analysis of uncertainty. In addition, the ratio of 90% of the predicted variation amplitude (predicted upper limit ( P u ) predicted lower limit ( P l )) to the predicted mean ( P m ) was used as the relative uncertainty ( RU ) to evaluate the amplitude of the variation in NPP.

2.5. Trend Analysis

The absolute interannual rate of change of the NPP residual was calculated based on the pixels using the method of unitary linear regression analysis:
θ s l o p e = n × i = 1 n i × Δ N P P i i = 1 n i i = 1 n Δ N P P i n × i = 1 n i 2 i = 1 n i 2
where n represents the time scale of the study; Δ NPPi represents the residual NPP of one pixel in year i; θslope > 0 and θslope < 0 represent increasing and decreasing trends, respectively; and |θslope| ≈ 0 indicates no change in the regional NPP residual.

3. Results

3.1. Verification of Net Primary Production Estimates

We used the CASA model to compare the NPP simulation (Figure 3a) and the MODIS data product (Figure 3b) for the sandy regions of northern China from 1983 to 2005. According to the spatial distribution, the simulation results can reflect the fluctuation of NPP in SAONC well, and the simulation results are relatively good for certain areas of sparse vegetation. The simulated data fit the MODIS data well under four climate models (R2 > 0.77, slope > 0.85). The field monitoring data based on 5 flux observation stations were compared with the simulation results (Figure 3d), and the correlation between the flux observation data and the simulation results was better (R2 > 0.8), indicating that the prediction of NPP based on climate model data was more suitable in SAONC.

3.2. Temporal Variation in Net Primary Productivity under Different Warming Scenarios

The future NPP values of the four emission scenarios under the temporal variation in NPP showed a rising trend compared with those in RP, and the increase in RCP8.5 was the greatest (Figure 4). Compared with the value in RP, the simulated NPP value of RCP8.5 increased by 26.68%, and the other three emission scenarios (RCP2.6, RCP4.5, and RCP6.0) showed little differences, with increases of 14.17, 10.72, and 8.57%, respectively; however, under GW_2.0 °C, the simulated NPP values of these RCPs increased by 20.87, 24.01, 29.31, and 39.94% relative to RP, respectively. The NPP changes in different emission scenarios are proportional to the level of the emission path. In a high-emission scenario, the content of NPP is relatively high, whereas in a low-emission scenario, the corresponding NPP content is relatively low. A difference in the corresponding NPP changes was observed between the two warming scenarios. Under GW_2.0 °C, the NPP contents in the high- and low-emission scenarios increased by 11.33 and 10.47% compared with GW_1.5 °C, respectively. From 1980 to 2100, the NPP contents in different emission scenarios showed fluctuating and increasing trends, and the values increased by 138.07, 165.75, 192.22, and 252.42 g C/m2·a in the RCP2.6, RCP4.5, RCP6.0, and RCP8.5 scenarios compared with that in 1981, respectively. In 2006–2029 and 2037–2042, marked increases in NPP were observed, whereas after 2050, this trend became flat.
On a seasonal scale, the NPP contents in the four emission scenarios changed to the greatest extent in summer, with an increase of 16.48% compared with that in RP under GW_1.5 °C. Compared with the NPP content in RP, the range of variation in winter was the smallest and tended to represent a straight line. Under GW_1.5 °C, the NPP contents in spring and autumn of the four emission scenarios were similar to those of RP and increased by 13.16 and 14.72%, respectively, compared with those in RP. Under GW_2.0 °C, the seasonal variation in NPP in the four emission scenarios was similar to the annual change, and the increase tended to be gradual. The most typical performances in summer were 15.23 and 33.27% higher than that in RP.
Overall, the variation ranges of NPP in SAONC will occur mainly between 2017 and 2036 (465–518 g C/m2, except RCP8.5) and 2042 and 2100 (527–614 g C/m2, except RCP8.5). The fluctuations in NPP under the RCP8.5 scenario increased significantly; however, after 2052 the increase tended to be gradual (634–668 g C/m2·a). In the future, seasonal variations in NPP in SAONC will occur in spring and autumn (287–418 g C/m2·a), summer (69–104 g C/m2·a), and winter (11–19 g C/m2·a).

3.3. Spatial Variation in Net Primary Productivity under Warming Scenarios

The spatial distribution of the annual NPP estimated by the four climate models was compared with that in RP at GW_1.5 °C (Figure 5). The areas showing relative decreasing trends in NPP contents accounted for <20% of the study area and were predominantly located in the northwest and northeast, whereas the areas showing increasing trends accounted for 45% of the study area and were mainly observed in the southeast; they presented relative increases between 0 and 50 g C/m2. The relative NPP contents simulated by the four climate models in spring mainly showed increasing trends, with increases ranging between 0 and 25 g C/m2, and these areas were mainly distributed in the central and eastern regions of the study area. The regions showing relative decreasing trends were predominantly located in the northwest and southeast. The IPSL-CM5A-LR model showed the greatest relative decrease, with the area experiencing a decreasing trend accounting for 37.52% of the study area. In summer, the range of relative NPP contents simulated by the four climate models was mainly concentrated between 0 and 50 g C/m2 and distributed throughout the study area, and this region accounted for >68% of the study area. The ranges of relative NPP contents simulated in autumn and winter by the four climate models were 0–25 g C/m2 and 0–10 g C/m2, respectively; however, the area showing a relative decreasing NPP content in autumn had a wider distribution than that in winter.
The spatial distribution of the relative NPP contents was simulated by the four climate models under GW_2.0 °C (Figure 6). The interannual variation shows that the primary range of the relative increase in NPP was 50 to 100 g C/m2, and these values were mainly located in the central and western regions of the study area. The region showing a relative decrease was mainly located in the northeastern region of the study area, especially in Heilongjiang province. The relative spatial variation in the NPP contents in spring, summer, and autumn decreased in the east and increased in the west, with relative increases mainly ranging from 25–50, 50–100, and 25–50 g C/m2, respectively. During winter, the MIROC5 model simulated relative decreases in NPP contents covering a proportional area of 11.90%, whereas the other three models mainly showed relative increases and only observed relative decreases over <1% of the study area.
Relative to GW_1.5 °C, the spatial distribution of estimated annual NPP under GW_2.0 °C showed a regular change from the east to the west in the total area. This pattern may occur because the spatial heterogeneity of the NPP content becomes more obvious with the further increase in temperature, thus exhibiting a certain longitudinal zonality.

4. Discussion

4.1. Impacts of Climate Variations on Net Primary Productivity in SAONC

The correlations between the changes in NPP and temperature, precipitation, and radiation under different future emission scenarios (2006–2100) under the warming scenarios are shown in Table 3. Temperature had the strongest correlation with the changes in NPP: in both the low (RCP2.6) and high (RCP8.5) emission scenarios, significant positive correlations were observed between NPP and temperature under the four climate models. Precipitation was found to be the second-most important factor (R > 0.017). These results are supported by those of Li et al. [71] and Li et al. [72], who used the Climate Model Intercomparison Project (CMIP5) multipattern to explore the response of global terrestrial NPP to climate change during the historical period from 1850–2005. However, our study also found that as emissions increased, the correlation between solar radiation and alterations in NPP content changed from positive under the low-emission scenario (RCP2.6) to negative under the high-emission scenario (RCP8.5).
The spatial variation analyses of these climatic factors under different warming scenarios can help illuminate the mechanisms driving future NPP changes in terrestrial ecosystems [18,73]. Under the warming scenarios (Figure 7), precipitation mainly increases in SAONC, although certain regions show decreasing trends, predominantly in SAONC and the Mu Us and Tengger Deserts in the central arid zone. Under the warming scenarios, the relative spatial variations in precipitation under the four climate models present a regularly decreasing trend from east to west. The MIROC5 climate model shows a relative trend of decreasing precipitation across the greatest proportion of the study area, reaching values of up to 23.56 and 24.88% under GW_1.5 and GW_2.0 °C, respectively. Under the warming scenarios, the relative spatial variation in solar radiation also shows a decreasing trend from east to west. However, the area that shows a decrease in solar radiation under the four climate models accounts for >43.52% of the total study area, and it is mainly located in the central and western regions of the study area.
Differences are observed among the models in terms of the areas that exhibit temperature increases under the warming scenarios (Figure 8). The GFDL-ESM2M model projections show that only small parts of the sandy regions demonstrated lower temperatures under global warming, and they were predominantly concentrated in inner Mongolia in northeastern China. The model results of HadGEM2-ES show temperature increase greater than 1.5 and 2.0 °C for nearly all of SAONC, including the northwestern sandy region of Xinjiang, which will experience the greatest warming, with elevated temperatures predicted to reach 3.0 and 4.0 °C. The results of the IPSL-CM5A-LR model show that areas are predominantly distributed in the southern Tarim Basin and southeast Qinghai province. The results from the MIROC5 model show small areas exhibiting warming below the global aggregated values of 1.5 and 2.0 °C located in northeastern China. Considered collectively, the four model-specific climate patterns under the warming scenarios indicate that the greatest increase in arid areas will occur across northwestern China. In some regions, the increase from the preindustrial temperature is projected to exceed 4.0 °C. These areas are mainly distributed in the northwestern sandy area of the Guerbantunggurt Desert and in the northwestern and central Taklimakan Desert.
Under warming ranging from the current level to the target set by the Paris Agreement, the increase in atmospheric CO2 concentration is the dominant factor driving the increase in the total volume of global NPP, whereas the contributions of temperature, precipitation, and radiation are relatively weak [74,75]. However, a clear spatial differentiation of changes in NPP with solar radiation occurs under the different warming targets. In the high-latitude regions of SAONC, the spatial changes in solar radiation are consistent with those of NPP, and solar radiation plays a certain role in promoting the change in NPP, whereas in the middle- and low-latitude regions solar radiation shows a certain negative effect. This effect may be explained by the fact that changes in NPP in high-latitude ecosystems in the north are mainly temperature dependent [26]. In these areas, the relatively rapid increase in temperature [32] has eased the temperature limit of vegetation to some extent, thus promoting the growth of NPP in these regions [18,76]. In contrast, the temperature increase at low and medium latitudes may lead to environmental temperatures that exceed the optimum temperature for vegetation growth, thus reducing the resistance of vegetation to drought and limiting the growth of NPP in the ecosystem [41]. With increases in emissions, the correlation between solar radiation and NPP changes gradually weakens, whereas that between precipitation and NPP changes increases, thus reflecting the importance of precipitation. Therefore, precipitation is the key factor dominating vegetation growth and succession in sandy land, which in turn affects the NPP changes [77,78].

4.2. Progress of Net Primary Productivity Simulations in SAONC

We have presented a study of the spatial and temporal patterns of and climate controls on NPP obtained from ISI-MIP2b climate models in SAONC from 1980 to 2100 through the use of a carbon model based on remote-sensing data, i.e., the CASA model [3]. This model has been successfully applied worldwide to map NPP patterns, including in mainland China [18,22]. Previous studies have reported similar analyses performed within shorter time periods [24]. Table 4 provides a comparison of the NPP prediction results of the present study with those of previous studies. The previous time scales simulated for SAONC were relatively short, and the results varied widely among simulations and presented an approximate range of 483–843 g C/m2·a. The present study simulated a value of 604.52 g C/m2·a, which extends beyond the range predicted by other studies that used the CASA model but is lower than the NPP change estimated by Gang et al. [79] using the CMIP5 dataset. This difference may be due to the relatively short time scales employed in previous studies and the absence of future changes in NPP in the context of global warming. However, future climate warming can directly affect future changes in NPP [23], which in turn can affect the average annual NPP.
Many studies have focused on how climate change will affect the future terrestrial ecosystem productivity in China; however, these studies have yielded widely different results. Existing studies have indicated that under the B2 scenario the total NPP within China’s terrestrial ecosystems will increase from 2.94 to 3.99 Pg C/a in the next 100 years [92]. In contrast, Zhao and Wu [87], who used the Lund-Potsdam-Jena (LPJ) model, observed a declining but fluctuating trend in natural vegetation NPP in China from 1961 to 1980. Wu et al. [93], who used the AVIMI model, found a similarly declining trend in the terrestrial ecosystem NPP in China. Ju et al. [94] used the Integrated Terrestrial Ecosystem Carbon (INTEC) model to predict changes in the forest ecosystem NPP in China under the A2 and B2 scenarios and observed an increasing trend from 2001 to 2100, with NPP increasing by 0.19 Pg C/a for the first 50 years and increasing by 0.15 Pg C/a for the second 50 years. Wen et al. [95] applied the Crop-C model to predict the NPP of China’s farmland ecosystem from 2000 to 2050 and predicted that NPP will increase at a rate of 0.0006 Pg C/a under the A1B scenario. Tao and Zhang [96] simulated the changes in China’s terrestrial ecosystem NPP until the end of this century under eight climatic scenarios; they observed an initial increase to a maximum in 2090, followed by a decrease, which they attributed to “drought stress”.
In conclusion, China’s future terrestrial ecosystem NPP may show an initial increasing trend followed by a decrease; however, the results of different regional studies vary widely and even provide contradictory results. The predicted responses of different vegetation types to future climate change also vary. The present study focused on desert vegetation and explored the changes in NPP of SAONC under future scenarios. The results indicate that the NPP in SAONC will initially increase, and then the rate of increase will be slow.

4.3. Quantifiable Sources of Uncertainty

Against the background of global warming, accurate estimations of NPP are relevant to global ecological systems [97]. However, considerable differences are observed among different studies. For example, the estimated value of global NPP is 56.2 ± 14.3 Gt C/a, with an uncertainty of approximately 25% [98]. The main sources of uncertainty in the NPP evaluation process of terrestrial ecosystems are the model structure, model parameters, and data sources. In terms of the data source, an uncertainty analysis of five flux observation sites was conducted in the current study (Figure 9a). The Tazhong station had the largest RU, possibly because it was located in an extremely arid area and was greatly influenced by the climate factors such as rainfall and temperature [60], which is important given that the meteorological conditions are the main factors affecting the change in NPP [99,100,101]. When using multiple models, the differences in the structures of the models used in different studies have a certain impact on the estimation of NPP [102,103]. Cramer, Kicklighter, Bondeau, Moore, Churkina, Nemry, Ruimy, Schloss, and Participants Potsdam [4] estimated the global NPP; although there was a 20% difference between their estimates, none of the models could explain this difference. Adams et al. [104] studied the relationships between different climate variables such as light, temperature, and moisture in different models with the function of NPP, but they failed to directly connect these differences with differences in NPP. The present study simulated a value of 604.52 g C/m2·a, which was close to the mean NPP value obtained by other studies. However, it is noteworthy that the NPP values of global and regional scales exhibit obvious interannual differences [43,105]. Therefore, the annual differences among different studies also introduce additional uncertainties [106,107]. The model parameters exert an important influence on the model simulation results, and a change in one parameter can lead to an obvious increase or decrease in the simulation results [108,109]. With other parameters fixed, the parameters requiring sensitivity analysis increased or decreased. According to Figure 9c, εmax, FPAR, NDVI, and SOL are direct linear variables of the model, and they are consistent with the amplitude of NPP. Among them, the value of εmax has the greatest influence on the estimated results of NPP. This result is consistent with those of [110], who considered that it was the most important parameter in the LUE model and an important factor in the low accuracy of the vegetation productivity model. The values of Tε1, Tε2, Wε, and Topt have a certain influence on the NPP output results. Among them, the value of Topt has the greatest influence on the estimated NPP results and leads to higher NPP amplitudes than the direct linear variable of the model.
In the present study, the changes in NPP compared with the changes in RP in the future global warming scenarios were compared and analyzed by synthesizing the simulated results of four climate models. Although the present study improved the accuracy of the research results to some extent, it should be noted that large differences persist among the different model parameters. Based on the above results, an uncertainty analysis of NPP under the four models was implemented (Figure 10), and the values ranged from 16.29 to 26.52%. The uncertainty of the HadGEM2-ES model was the largest, and the uncertainty of the MIROC5 model was the smallest. These results are consistent with those obtained by Li, Lu, Zhang, Liu, Gao, and Ao [71] and Fu et al. [111]. In fact, multiple studies have found large uncertainty associated with the NPP of the global terrestrial ecosystem. However, because of the lack of long-term global NPP data [112,113], the validation of modeled NPP values cannot be as rigorously implemented at regional and global scales as observed for validations of climatic simulations. In the present study, we utilized ChinaFlux data and MODIS NPP data to compare the NPP values of the ISI-MIP2b models, and the results highlighted the need to improve our knowledge of the factors governing variability within model estimations of NPP in sandy regions.

5. Conclusions

We investigated the spatiotemporal features of NPP in SAONC using four ISI-MIP2b datasets and the CASA model. The SAONC exceeded the global warming level at 0–7.1 and 0–12.5%, respectively, with some areas experiencing 4.0 °C increases under GW_1.5 °C or under GW_2.0 °C. The spatiotemporal variation in NPP was documented with a modified model, and an assessment of the impacts of climate change on NPP formation when global temperatures were stabilized under the two warming scenarios was carried out using four ISI-MIP2b datasets. The results indicated greater spatiotemporal variation in NPP in SAONC when global temperatures were stabilized under the warming scenarios. The temporal variations in annual NPP showed an increasing trend relative to the reference period. Under GW_1.5 °C and GW_2.0 °C, the estimated annual NPP increased by 14.17, 10.72, 8.57, and 26.68% and by 20.87, 24.01, 29.31, and 39.94% in the RCP2.6, RCP4.5, RCP6.0, and RCP8.5 scenarios, respectively. The spatial distribution of the estimated annual NPP under GW_1.5 °C showed clear increasing trends in the southeast (accounting for 45% of the total area) and decreasing trends in the northwest and northeast (accounting for 20%). Relative to GW_1.5 °C, the spatial distribution of the estimated annual NPP under GW_2.0 °C showed a regular change from the east to the west. However, our results showed a robust nonlinear relationship between temperature and NPP growth by accounting for intrayear monthly variability in temperature and precipitation as well as monthly temperature and precipitation maximum and minimum values and outliers. The changes in temperature (R > 0.614) had a greater impact on the NPP growth than did precipitation (R > 0.017), and solar radiation showed a negative impact in the middle- and low-latitude regions.
The projected impacts on NPP under GW_1.5 °C relative to no additional warming (and relative to 2.0 °C) are uncertain because of the large range of likely outcomes within the model when accounting for both estimation uncertainties regarding climate and NPP. The uncertainty of NPP under the four models ranged from 16.29 to 26.52%. The results imply a large impact of GW_1.5 °C compared with the current conditions, whereas GW_2.0 °C leads to a significantly lower growth in projected NPP over all study areas. It must be emphasized that the results solely represent the projected net effect of temperature under the assumption of stable estimated relationships, and large overall uncertainty in NPP growth is observed under the hypothetical temperature scenarios. However, the results suggest that considerable benefits may arise because of the lower warming levels. This finding further motivates efforts to limit GMST warming to 1.5 °C through the use of stringent emission adjustments.

Author Contributions

Conceptualization, X.M., C.Z. and T.H.; formal analysis, C.Z.; funding acquisition, C.Z.; methodology, X.M., T.H. and X.Z.; resources, C.Z.; supervision, C.Z.; visualization, C.Z. and X.M.; writing—original draft, X.M., W.Y.; writing—review and editing, C.Z., W.Y. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Project of National Natural Science Foundation (41671030) and the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA2006030201).

Acknowledgments

The authors would like to express their appreciation to China FLUX Observed and Research Network (ChinaFLUX) and China Ecosystem Research Network (CERN) for the observed data from Naiman, Yanchi, Shapotou, Fukang, Tazhong and Aksu stations. We also wish to thank Tong Jiang and Buda Su, who helped improve this manuscript through their thorough review.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keeling, C.D.; Chin, J.F.S.; Whorf, T.P. Increased activity of northern vegetation inferred from atmospheric CO2 measurements. Nature 1996, 382, 146–149. [Google Scholar] [CrossRef]
  2. Roxburgh, S.H.; Berry, S.L.; Buckley, T.N.; Barnes, B.; Roderick, M.L. What is NPP? Inconsistent accounting of respiratory fluxes in the definition of net primary production. Funct. Ecol. 2005, 19, 378–382. [Google Scholar] [CrossRef]
  3. Potter, C.S.; Randerson, J.T.; Field, C.B.; Matson, P.A.; Vitousek, P.M.; Mooney, H.A.; Klooster, S.A. Terrestrial ecosystem production: A process model based on global satellite and surface data. Glob. Biogeochem. Cycle. 1993, 7, 811–841. [Google Scholar] [CrossRef]
  4. Cramer, W.; Kicklighter, D.W.; Bondeau, A.; Moore, B.; Churkina, G.; Nemry, B.; Ruimy, A.; Schloss, A.L.; The Participants of The Potsdam NpP Model Intercomparison. Comparing global models of terrestrial net primary productivity (NPP): Overview and key results. Glob. Change Biol. 1999, 5, 1–15. [Google Scholar] [CrossRef]
  5. Lehuger, S.; Gabrielle, B.; Cellier, P.; Loubet, B.; Roche, R.; Beziat, P.; Ceschia, E.; Wattenbach, M. Predicting the net carbon exchanges of crop rotations in Europe with an agro-ecosystem model. Agric. Ecosyst. Environ. 2010, 139, 384–395. [Google Scholar] [CrossRef] [Green Version]
  6. Xu, X.; Niu, S.L.; Sherry, R.A.; Zhou, X.H.; Zhou, J.Z.; Luo, Y.Q. Interannual variability in responses of belowground net primary productivity (NPP) and NPP partitioning to long-term warming and clipping in a tallgrass prairie. Glob. Change Biol. 2012, 18, 1648–1656. [Google Scholar] [CrossRef]
  7. Eisfelder, C.; Klein, I.; Niklaus, M.; Kuenzer, C. Net primary productivity in Kazakhstan, its spatio-temporal patterns and relation to meteorological variables. J. Arid. Environ. 2014, 103, 17–30. [Google Scholar] [CrossRef] [Green Version]
  8. Chang, J.F.; Ciais, P.; Wang, X.H.; Piao, S.L.; Asrar, G.; Betts, R.; Chevallier, F.; Dury, M.; Francois, L.; Frieler, K.; et al. Benchmarking carbon fluxes of the ISIMIP2a biome models. Environ. Res. Lett. 2017, 12, 16. [Google Scholar] [CrossRef] [Green Version]
  9. Li, Y.; Wang, Y.G.; Houghton, R.A.; Tang, L.S. Hidden carbon sink beneath desert. Geophys. Res. Lett. 2015, 42, 5880–5887. [Google Scholar] [CrossRef]
  10. Li, J.; Zhao, C.; Zhu, H.; Li, Y.; Wang, F. Effect of plant species on shrub fertile island at an oasis-desert ecotone in the South Junggar Basin, China. J. Arid. Environ. 2007, 71, 350–361. [Google Scholar] [CrossRef]
  11. Liu, Z.; Yao, Z.; Huang, H.; Wu, S.; Liu, G. Land use and climate changes and their impacts on runoff in the Yarlung Zangbo river basin, China. Land Degrad. Dev. 2014, 25, 203–215. [Google Scholar] [CrossRef]
  12. Stone, R. Ecosystems-Have desert researchers discovered a hidden loop in the carbon cycle? Science 2008, 320, 1409–1410. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, X.N.; Zhao, C.Y.; Wang, J.Y.; Stahr, K.; Kuzyakov, Y. CaCO3 recrystallization in saline and alkaline soils. Geoderma 2016, 282, 1–8. [Google Scholar] [CrossRef]
  14. Zhang, C.; Lu, D.S.; Chen, X.; Zhang, Y.M.; Maisupova, B.; Tao, Y. The spatiotemporal patterns of vegetation coverage and biomass of the temperate deserts in Central Asia and their relationships with climate controls. Remote Sens. Environ. 2016, 175, 271–281. [Google Scholar] [CrossRef]
  15. McKeon, G.M.; Stone, G.S.; Syktus, J.I.; Carter, J.O.; Flood, N.R.; Ahrens, D.G.; Bruget, D.N.; Chilcott, C.R.; Cobon, D.H.; Cowley, R.A.; et al. Climate change impacts on northern Australian rangeland livestock carrying capacity: A review of issues. Rangel. J. 2009, 31, 1–29. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, P.J.; Xie, D.H.; Zhou, Y.Y.; Youhao, E.; Zhu, Q.J. Estimation of net primary productivity using a process-based model in Gansu Province, Northwest China. Environ. Earth Sci. 2014, 71, 647–658. [Google Scholar] [CrossRef]
  17. Li, J.; Zhao, C.Y.; Song, Y.J.; Sheng, Y.; Zhu, H. Spatial patterns of desert annuals in relation to shrub effects on soil moisture. J. Veg. Sci. 2010, 21, 221–232. [Google Scholar] [CrossRef]
  18. Piao, S.; Fang, J.; Ciais, P.; Peylin, P.; Huang, Y.; Sitch, S.; Wang, T. The carbon balance of terrestrial ecosystems in China. Nature 2009, 458, 1009–1013. [Google Scholar] [CrossRef]
  19. Lieth, H. Modeling the Primary Productivity of the World. In Primary Productivity of the Biosphere; Lieth, H., Whittaker, R.H., Eds.; Springer: Berlin/Heidelberg, Germany, 1975; pp. 237–263. [Google Scholar]
  20. Parton, W.J.; Scurlock, J.M.O.; Ojima, D.S.; Gilmanov, T.G.; Scholes, R.J.; Schimel, D.S.; Kirchner, T.; Menaut, J.C.; Seastedt, T.; Moya, E.G.; et al. Observations and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide. Glob. Biogeochem. Cycle. 1993, 7, 785–809. [Google Scholar] [CrossRef] [Green Version]
  21. Zhu, W.Q.; Pan, Y.Z.; Liu, X.; Wang, A.L. Spatio-temporal distribution of net primary productivity along the Northeast China Transect and its response to climatic change. J. For. Res. 2006, 17, 93–98. [Google Scholar] [CrossRef]
  22. Zhang, M.L.; Lal, R.; Zhao, Y.Y.; Jiang, W.L.; Chen, Q.G. Estimating net primary production of natural grassland and its spatio-temporal distribution in China. Sci. Total Environ. 2016, 553, 184–195. [Google Scholar] [CrossRef] [PubMed]
  23. Piao, S.L.; Fang, J.Y.; Zhou, L.M.; Zhu, B.; Tan, K.; Tao, S. Changes in vegetation net primary productivity from 1982 to 1999 in China. Glob. Biogeochem. Cycle. 2005, 19, 19. [Google Scholar] [CrossRef] [Green Version]
  24. Liang, W.; Yang, Y.T.; Fan, D.M.; Guan, H.D.; Zhang, T.; Long, D.; Zhou, Y.; Bai, D. Analysis of spatial and temporal patterns of net primary production and their climate controls in China from 1982 to 2010. Agric. For. Meteorol. 2015, 204, 22–36. [Google Scholar] [CrossRef]
  25. Potter, C. Microclimate influences on vegetation water availability and net primary production in coastal ecosystems of Central California. Landsc. Ecol. 2014, 29, 677–687. [Google Scholar] [CrossRef]
  26. Nemani, R.R.; Keeling, C.D.; Hashimoto, H.; Jolly, W.M.; Piper, S.C.; Tucker, C.J.; Myneni, R.B.; Running, S.W. Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science 2003, 300, 1560–1563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Piao, S.L.; Ciais, P.; Friedlingstein, P.; Peylin, P.; Reichstein, M.; Luyssaert, S.; Margolis, H.; Fang, J.Y.; Barr, A.; Chen, A.P.; et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 2008, 451, 49–52. [Google Scholar] [CrossRef]
  28. Field, C.B.; Randerson, J.T.; Malmstrom, C.M. Global net primary production: Combining ecology and remote sensing. Remote Sens. Environ. 1995, 51, 74–88. [Google Scholar] [CrossRef] [Green Version]
  29. Li, A.N.; Bian, J.H.; Lei, G.B.; Huang, C.Q. Estimating the Maximal Light Use Efficiency for Different Vegetation through the CASA Model Combined with Time-Series Remote Sensing Data and Ground Measurements. Remote Sens. 2012, 4, 3857–3876. [Google Scholar] [CrossRef] [Green Version]
  30. Gao, Z.Q.; Liu, J.Y. Simulation study of China’s net primary production. Chin. Sci. Bull. 2008, 53, 434–443. [Google Scholar] [CrossRef]
  31. Yu, D.; Shao, H.B.; Shi, P.J.; Zhu, W.Q.; Pan, Y.Z. How does the conversion of land cover to urban use affect net primary productivity? A case study in Shenzhen city, China. Agric. For. Meteorol. 2009, 149, 2054–2060. [Google Scholar]
  32. Stocker, T.F.; Plattner, G.K. Rethink IPCC reports. Nature 2014, 513, 163–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Frieler, K.; Lange, S.; Piontek, F.; Reyer, C.P.O.; Schewe, J.; Warszawski, L.; Zhao, F.; Chini, L.; Denvil, S.; Emanuel, K.; et al. Assessing the impacts of 1.5 degrees C global warming-simulation protocol of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2b). Geosci. Model Dev. 2017, 10, 4321–4345. [Google Scholar] [CrossRef] [Green Version]
  34. Schleussner, C.F.; Lissner, T.K.; Fischer, E.M.; Wohland, J.; Perrette, M.; Golly, A.; Rogelj, J.; Childers, K.; Schewe, J.; Frieler, K.; et al. Differential climate impacts for policy-relevant limits to global warming: The case of 1.5 degrees C and 2 degrees C. Earth Syst. Dynam. 2016, 7, 327–351. [Google Scholar] [CrossRef] [Green Version]
  35. Huang, J.P.; Yu, H.P.; Dai, A.G.; Wei, Y.; Kang, L.T. Drylands face potential threat under 2 degrees C global warming target. Nat. Clim. Chang. 2017, 7, 417–422. [Google Scholar] [CrossRef]
  36. Su, B.D.; Jian, D.N.; Li, X.C.; Wang, Y.J.; Wang, A.Q.; Wen, S.S.; Tao, H.; Hartmann, H. Projection of actual evapotranspiration using the COSMO-CLM regional climate model under global warming scenarios of 1.5 degrees C and 2.0 degrees C in the Tarim River basin, China. Atmos. Res. 2017, 196, 119–128. [Google Scholar] [CrossRef]
  37. Schaeffer, M.; Hare, W.; Rahmstorf, S.; Vermeer, M. Long-term sea-level rise implied by 1.5 degrees C and 2 degrees C warming levels. Nat. Clim. Chang. 2012, 2, 867–870. [Google Scholar] [CrossRef]
  38. Karmalkar, A.V.; Bradley, R.S. Consequences of Global Warming of 1.5 degrees C and 2 degrees C for Regional Temperature and Precipitation Changes in the Contiguous United States. PLoS ONE 2017, 12, e0168697. [Google Scholar] [CrossRef] [Green Version]
  39. Wang, X.P.; Young, M.H.; Yu, Z.; Li, X.R.; Zhang, Z.S. Long-term effects of restoration on soil hydraulic properties in revegetation-stabilized desert ecosystems. Geophys. Res. Lett. 2007, 34, 1061–1064. [Google Scholar] [CrossRef]
  40. Li, X.R.; Ma, F.Y.; Xiao, H.L.; Wang, X.P.; Kim, K.C. Long-term effects of revegetation on soil water content of sand dunes in arid region of Northern China. J. Arid. Environ. 2004, 57, 1–16. [Google Scholar] [CrossRef]
  41. Sitch, S.; Smith, B.; Prentice, I.C.; Arneth, A.; Bondeau, A.; Cramer, W.; Kaplan, J.O.; Levis, S.; Lucht, W.; Sykes, M.T.; et al. Evaluation of ecosystem dynamics, plant geography and terrestrial carbon cycling in the LPJ dynamic global vegetation model. Glob. Chang. Biol. 2003, 9, 161–185. [Google Scholar] [CrossRef]
  42. Ajaj, Q.M.; Pradhan, B.; Noori, A.M.; Jebur, M.N. Spatial Monitoring of Desertification Extent in Western Iraq using Landsat Images and GIS. Land Degrad. Dev. 2017, 28, 2418–2431. [Google Scholar] [CrossRef]
  43. Taylor, K.E.; Stouffer, R.J.; Meehl, G.A. An overview of CMIP5 and the experiment design. Bull. Amer. Meteorol. Soc. 2012, 93, 485–498. [Google Scholar] [CrossRef] [Green Version]
  44. Elguindi, N.; Grundstein, A.; Bernardes, S.; Turuncoglu, U.; Feddema, J. Assessment of CMIP5 global model simulations and climate change projections for the 21 (st) century using a modified Thornthwaite climate classification. Clim. Chang. 2014, 122, 523–538. [Google Scholar] [CrossRef]
  45. Chadwick, R.; Boutle, I.; Martin, G. Spatial Patterns of Precipitation Change in CMIP5: Why the Rich Do Not Get Richer in the Tropics. J. Clim. 2013, 26, 3803–3822. [Google Scholar] [CrossRef]
  46. Fischer, E.M.; Knutti, R. Observed heavy precipitation increase confirms theory and early models. Nat. Clim. Chang. 2016, 6, 986–991. [Google Scholar] [CrossRef]
  47. Miao, C.Y.; Duan, Q.Y.; Sun, Q.H.; Huang, Y.; Kong, D.X.; Yang, T.T.; Ye, A.Z.; Di, Z.H.; Gong, W. Assessment of CMIP5 climate models and projected temperature changes over Northern Eurasia. Environ. Res. Lett. 2014, 9, 055007. [Google Scholar] [CrossRef]
  48. Hoerling, M.P.; Eischeid, J.K.; Quan, X.W.; Diaz, H.F.; Webb, R.S.; Dole, R.M.; Easterling, D.R. Is a Transition to Semipermanent Drought Conditions Imminent in the US Great Plains? J. Clim. 2012, 25, 8380–8386. [Google Scholar] [CrossRef]
  49. Wuebbles, D.; Meehl, G.; Hayhoe, K.; Karl, T.R.; Kunkel, K.; Santer, B.; Wehner, M.; Colle, B.; Fischer, E.M.; Fu, R.; et al. CMIP5 climate model analyses climate extremes in the United States. Bull. Amer. Meteorol. Soc. 2014, 95, 571–583. [Google Scholar] [CrossRef] [Green Version]
  50. Marx, A.; Kumar, R.; Thober, S.; Rakovec, O.; Wanders, N.; Zink, M.; Wood, E.F.; Pan, M.; Sheffield, J.; Samaniego, L. Climate change alters low flows in Europe under global warming of 1.5, 2, and 3 degrees C. Hydrol. Earth Syst. Sci. 2018, 22, 1017–1032. [Google Scholar] [CrossRef] [Green Version]
  51. Liu, M.; Dries, L.; Heijman, W.; Huang, J.K.; Zhu, X.Q.; Hu, Y.N.; Chen, H.B. The impact of ecological construction programs on grassland conservation in Inner Mongolia, China. Land Degrad. Dev. 2018, 29, 326–336. [Google Scholar] [CrossRef] [Green Version]
  52. Yu, G.R.; Wen, X.F.; Sun, X.M.; Tanner, B.D.; Lee, X.H.; Chen, J.Y. Overview of ChinaFLUX and evaluation of its eddy covariance measurement. Agric. For. Meteorol. 2006, 137, 125–137. [Google Scholar] [CrossRef]
  53. Wang, K.C.; Dickinson, R.E. A review of global terrestrial evapotranspiration: Observation, modeling, climatology, and climatic variability. Rev. Geophys. 2012, 50, 54. [Google Scholar] [CrossRef]
  54. Aubinet, M.; Grelle, A.; Ibrom, A.; Rannik, U.; Moncrieff, J.; Foken, T.; Kowalski, A.S.; Martin, P.H.; Berbigier, P.; Bernhofer, C.; et al. Estimates of the annual net carbon and water exchange of forests: The EUROFLUX methodology. Adv. Ecol. Res. 2000, 30, 113–175. [Google Scholar]
  55. Webb, E.K.; Pearman, G.I.; Leuning, R. Correction of flux measurements for density effects due to heat and water-vapor transfer. Q. J. R. Meteorol. Soc. 1980, 106, 85–100. [Google Scholar] [CrossRef]
  56. Zomer, R.J.; Trabucco, A.; Metzger, M.J.; Wang, M.C.; Oli, K.P.; Xu, J.C. Projected climate change impacts on spatial distribution of bioclimatic zones and ecoregions within the Kailash Sacred Landscape of China, India, Nepal. Clim. Chang. 2014, 125, 445–460. [Google Scholar] [CrossRef]
  57. Feng, Q.; Cheng, G.D.; Endo, K.N. Water content variations and respective ecosystems of sandy land in China. Environ. Geol. 2001, 40, 1075–1083. [Google Scholar]
  58. Zheng, H.; Yu, G.R.; Wang, Q.F.; Zhu, X.J.; Yan, J.H.; Wang, H.M.; Shi, P.L.; Zhao, F.H.; Li, Y.N.; Zhao, L.; et al. Assessing the ability of potential evapotranspiration models in capturing dynamics of evaporative demand across various biomes and climatic regimes with ChinaFLUX measurements. J. Hydrol. 2017, 551, 70–80. [Google Scholar] [CrossRef]
  59. Fu, Y.L.; Yu, G.R.; Sun, X.M.; Li, Y.N.; Wen, X.F.; Zhang, L.M.; Li, Z.Q.; Zhao, L.; Hao, Y.B. Depression of net ecosystem CO2 exchange in semi-arid Leymus chinensis steppe and alpine shrub. Agric. For. Meteorol. 2006, 137, 234–244. [Google Scholar] [CrossRef]
  60. Yang, X.H.; He, Q.; Mamtimin, A.; Huo, W.; Liu, X.C. Diurnal variations of saltation activity at Tazhong: The hinterland of Taklimakan Desert. Meteorol. Atmos. Phys. 2013, 119, 177–185. [Google Scholar] [CrossRef] [Green Version]
  61. Tucker, C.J.; Pinzon, J.E.; Brown, M.E.; Slayback, D.A.; Pak, E.W.; Mahoney, R.; Vermote, E.F.; El Saleous, N. An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. Int. J. Remote Sens. 2005, 26, 4485–4498. [Google Scholar] [CrossRef]
  62. Zhang, Y.; Zhang, X.L. Estimation of net primary productivity of different forest types based on improved CASA model in Jing-Jin-Ji region, China. J. Sustain. For. 2017, 36, 568–582. [Google Scholar] [CrossRef]
  63. UNFCCC. UNFCCC Conference of the Parties: Adoption of the Paris Agreement; UNFCCC: Paris, France, 2015; Available online: https://unfccc.int/resource/docs/2015/cop21/eng/l09r01.pdf (accessed on 11 December 2019).
  64. Huang, J.P.; Ji, M.X.; Xie, Y.K.; Wang, S.S.; He, Y.L.; Ran, J.J. Global semi-arid climate change over last 60 years. Clim. Dyn. 2016, 46, 1131–1150. [Google Scholar] [CrossRef] [Green Version]
  65. Hartmann, D.L.; Tank, A.M.G.K.; Rusticucci, M.; Alexander, L.V.; Brönnimann, S.; Charabi, Y.; Dentener, F.J.; Dlugokencky, E.J.; Easterling, D.R.; Kaplan, A.; et al. Observations: Atmosphere and Surface. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013; pp. 159–254. [Google Scholar]
  66. Warszawski, L.; Frieler, K.; Huber, V.; Piontek, F.; Serdeczny, O.; Schewe, J. The Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP): Project framework. Proc. Natl. Acad. Sci. USA 2014, 111, 3228–3232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Yu, Z.; Wang, J.X.; Liu, S.R.; Rentch, J.S.; Sun, P.S.; Lu, C.Q. Global gross primary productivity and water use efficiency changes under drought stress. Environ. Res. Lett. 2017, 12, 10. [Google Scholar] [CrossRef] [Green Version]
  68. Bradford, J.B.; Hicke, J.A.; Lauenroth, W.K. The relative importance of light-use efficiency modifications from environmental conditions and cultivation for estimation of large-scale net primary productivity. Remote Sens. Environ. 2005, 96, 246–255. [Google Scholar] [CrossRef]
  69. Mahamah, D.S. Simplified sensitivity analysis applied to a nutrient-biomass model. Ecol. Model. 1988, 42, 103–109. [Google Scholar] [CrossRef]
  70. Yan, Y.A.; Wang, S.Q.; Wang, Y.D.; Wu, W.X.; Wang, J.Y.; Chen, B.; Yang, F.T. Assessing productivity and carbon sequestration capacity of subtropical coniferous plantations using the process model PnET-CN. J. Geogr. Sci. 2011, 21, 458–474. [Google Scholar] [CrossRef]
  71. Li, S.S.; Lu, S.H.; Zhang, Y.J.; Liu, Y.P.; Gao, Y.H.; Ao, Y.H. The change of global terrestrial ecosystem net primary productivity (NPP) and its response to climate change in CMIP5. Theor. Appl. Climatol. 2015, 121, 319–335. [Google Scholar] [CrossRef]
  72. Li, S.S.; Lu, S.H.; Liu, Y.P.; Gao, Y.H.; Ao, Y.H. Variations and trends of terrestrial NPP and its relation to climate change in the 10 CMIP5 models. J. Earth Syst. Sci. 2015, 124, 395–403. [Google Scholar] [CrossRef] [Green Version]
  73. Xu, X.; Sherry, R.A.; Niu, S.L.; Li, D.J.; Luo, Y.Q. Net primary productivity and rain-use efficiency as affected by warming, altered precipitation, and clipping in a mixed-grass prairie. Glob. Change Biol. 2013, 19, 2753–2764. [Google Scholar] [CrossRef]
  74. Beer, C.; Reichstein, M.; Tomelleri, E.; Ciais, P.; Jung, M.; Carvalhais, N.; Rodenbeck, C.; Arain, M.A.; Baldocchi, D.; Bonan, G.B.; et al. Terrestrial Gross Carbon Dioxide Uptake: Global Distribution and Covariation with Climate. Science 2010, 329, 834–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Zhao, C.Y.; Wang, Y.C.; Hu, S.J.; Li, Y. Effects of spatial variability on estimation of evapotranspiration in the continental river basin. J. Arid. Environ. 2004, 56, 373–382. [Google Scholar] [CrossRef]
  76. Fang, J.Y.; Piao, S.L.; Field, C.B.; Pan, Y.D.; Guo, Q.H.; Zhou, L.M.; Peng, C.H.; Tao, S. Increasing net primary production in China from 1982 to 1999. Front. Ecol. Environ. 2003, 1, 293–297. [Google Scholar] [CrossRef]
  77. Li, F.; Zhao, W.Z.; Liu, H. Productivity responses of desert vegetation to precipitation patterns across a rainfall gradient. J. Plant Res. 2015, 128, 283–294. [Google Scholar] [CrossRef]
  78. Zhao, C.Y.; Wang, Y.C.; Chen, X.; Li, B.G. Simulation of the effects of groundwater level on vegetation change by combining FEFLOW software. Ecol. Model. 2005, 187, 341–351. [Google Scholar] [CrossRef]
  79. Gang, C.C.; Wang, Z.Q.; Zhou, W.; Chen, Y.Z.; Li, J.L.; Cheng, J.M.; Guo, L.; Odeh, I.; Chen, C. Projecting the dynamics of terrestrial net primary productivity in response to future climate change under the RCP2.6 scenario. Environ. Earth Sci. 2015, 74, 5949–5959. [Google Scholar] [CrossRef]
  80. Sun, G.D.; Mu, M. Understanding variations and seasonal characteristics of net primary production under two types of climate change scenarios in China using the LPJ model. Clim. Chang. 2013, 120, 755–769. [Google Scholar] [CrossRef]
  81. Li, Z.; Pan, J.H. Spatiotemporal changes in vegetation net primary productivity in the arid region of Northwest China, 2001 to 2012. Front. Earth Sci. 2018, 12, 108–124. [Google Scholar] [CrossRef]
  82. Piao, S.; Fang, J. Terrestrial net primary production and its spatio-temporal patterns in Qinghai-Xizang Plateau, China during 1982-1999. J. Nat. Resour. 2002, 17, 373–380. [Google Scholar]
  83. Chen, L.; Liu, G.; Li, H. Estimating Net Primary Productivity of Terrestrial Vegetation in China Using Remote Sensing. J. Remote Sens. 2002, 6, 129–135. [Google Scholar]
  84. Cao, M.K.; Tao, B.; Li, K.R.; Shao, X.M.; Prience, S.D. Interannual variation in terrestrial ecosystem carbon fluxes in China from 1981 to 1998. Acta Bot. Sin. 2003, 45, 552–560. [Google Scholar]
  85. Feng, X.; Liu, G.; Chen, J.M.; Chen, M.; Liu, J.; Ju, W.M.; Sun, R.; Zhou, W. Net primary productivity of China’s terrestrial ecosystems from a process model driven by remote sensing. J. Environ. Manag. 2007, 85, 563–573. [Google Scholar] [CrossRef] [PubMed]
  86. Chen, B.; Shaoqiang, W.; Ronggao, L.; Ting, S. Study on Modeling and Spatial Pattern of Net Primary Production in China’s Terrestrial Ecosystem. Resour. Sci. 2007, 29, 45–53. [Google Scholar]
  87. Zhao, D.S.; Wu, S.H. Vulnerability of natural ecosystem in China under regional climate scenarios: An analysis based on eco-geographical regions. J. Geogr. Sci. 2014, 24, 237–248. [Google Scholar] [CrossRef]
  88. Mao, J.F.; Dan, L.; Wang, B.; Dai, Y.J. Simulation and evaluation of terrestrial ecosystem NPP with M-SDGVM over continental China. Adv. Atmos. Sci. 2010, 27, 427–442. [Google Scholar] [CrossRef]
  89. Chen, F.; Shen, Y.; Li, Q.; Guo, Y.; Xu, L. Spatio-temporal Variation Analysis of Ecological Systems NPP in China in Past 30 years. Sci. Geogr. Sin. 2011, 31, 1409–1414. [Google Scholar]
  90. Liu, Y.B.; Ju, W.M.; He, H.L.; Wang, S.Q.; Sun, R.; Zhang, Y.D. Changes of net primary productivity in China during recent 11 years detected using an ecological model driven by MODIS data. Front. Earth Sci. 2013, 7, 112–127. [Google Scholar] [CrossRef]
  91. Pei, F.S.; Li, X.; Liu, X.P.; Lao, C.H. Assessing the impacts of droughts on net primary productivity in China. J. Environ. Manag. 2013, 114, 362–371. [Google Scholar] [CrossRef]
  92. Ji, J.J.; Huang, M.; Li, K.R. Prediction of carbon exchanges between China terrestrial ecosystem and atmosphere in 21st century. Sci. China Ser. D-Earth Sci. 2008, 51, 885–898. [Google Scholar] [CrossRef]
  93. Wu, S.H.; Yin, Y.H.; Zhao, D.S.; Huang, M.; Shao, X.M.; Dai, E.F. Impact of future climate change on terrestrial ecosystems in China. Int. J. Climatol. 2010, 30, 866–873. [Google Scholar] [CrossRef]
  94. Ju, W.M.; Chen, J.M.; Harvey, D.; Wang, S. Future carbon balance of China’s forests under climate change and increasing CO2. J. Environ. Manag. 2007, 85, 538–562. [Google Scholar] [CrossRef]
  95. Wen, Z.; Yao, H.; Wenjuan, S.J.; Yongqiang, Y.Q. Simulating crop net primary production in China from 2000 to 2050 by linking the crop-C model with a FGOALS’s model climate change scenario. Adv. Atmos. Sci. 2007, 24, 845–854. [Google Scholar]
  96. Tao, F.L.; Zhang, Z. Dynamic responses of terrestrial ecosystems structure and function to climate change in China. J. Geophys. Res.-Biogeosci. 2010, 115, 58–72. [Google Scholar] [CrossRef]
  97. Pretis, F.; Schwarz, M.; Tang, K.; Haustein, K.; Allen, M.R. Uncertain impacts on economic growth when stabilizing global temperatures at 1.5°C or 2°C warming. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2018, 376, 20160460. [Google Scholar] [CrossRef] [PubMed]
  98. Ito, A. Decadal Variability in the Terrestrial Carbon Budget Caused by the Pacific Decadal Oscillation and Atlantic Multidecadal Oscillation. J. Meteorol. Soc. Jpn. 2011, 89, 441–454. [Google Scholar] [CrossRef] [Green Version]
  99. Shafran Nathan, R.; Svoray, T.; Perevolotsky, A. The resilience of annual vegetation primary production subjected to different climate change scenarios. Clim. Chang. 2013, 118, 227–243. [Google Scholar] [CrossRef]
  100. Shao, P.; Zeng, X.B.; Sakaguchi, K.; Monson, R.K.; Zeng, X.D. Terrestrial Carbon Cycle: Climate Relations in Eight CMIP5 Earth System Models. J. Clim. 2013, 26, 8744–8764. [Google Scholar] [CrossRef]
  101. Finstad, A.G.; Hein, C.L. Migrate or stay: Terrestrial primary productivity and climate drive anadromy in Arctic char. Glob. Change Biol. 2012, 18, 2487–2497. [Google Scholar] [CrossRef]
  102. Gurney, K.R.; Law, R.M.; Denning, A.S.; Rayner, P.J.; Baker, D.; Bousquet, P.; Bruhwiler, L.; Chen, Y.H.; Ciais, P.; Fan, S.M.; et al. TransCom 3 CO2 inversion intercomparison: 1. Annual mean control results and sensitivity to transport and prior flux information. Tellus Ser. B-Chem. Phys. Meteorol. 2003, 55, 555–579. [Google Scholar] [CrossRef] [Green Version]
  103. Bondeau, A.; Kicklighter, D.W.; Kaduk, J.; The Participants of The Potsdam NpP Model Intercomparison. Comparing global models of terrestrial net primary productivity (NPP): Importance of vegetation structure on seasonal NPP estimates. Glob. Change Biol. 1999, 5, 35–45. [Google Scholar] [CrossRef]
  104. Adams, B.; White, A.; Lenton, T.M. An analysis of some diverse approaches to modelling terrestrial net primary productivity. Ecol. Model. 2004, 177, 353–391. [Google Scholar] [CrossRef]
  105. Baldocchi, D. Breathing of the terrestrial biosphere: Lessons learned from a global network of carbon dioxide flux measurement systems. Aust. J. Bot. 2008, 56, 1–26. [Google Scholar] [CrossRef]
  106. Patra, P.K.; Ishizawa, M.; Maksyutov, S.; Nakazawa, T.; Inoue, G. Role of biomass burning and climate anomalies for land-atmosphere carbon fluxes based on inverse modeling of atmospheric CO2. Glob. Biogeochem. Cycle. 2005, 19, 15. [Google Scholar] [CrossRef] [Green Version]
  107. Sarmiento, J.L.; Lequere, C.; Pacala, S.W. Limiting future atmospheric carbon dioxide. Glob. Biogeochem. Cycle. 1995, 9, 121–137. [Google Scholar] [CrossRef]
  108. Zhu, Q.; Zhuang, Q.L. Parameterization and sensitivity analysis of a process-based terrestrial ecosystem model using adjoint method. J. Adv. Model. Earth Syst. 2014, 6, 315–331. [Google Scholar] [CrossRef]
  109. Huang, N.; Niu, Z.; Wu, C.Y.; Tappert, M.C. Modeling net primary production of a fast-growing forest using a light use efficiency model. Ecol. Model. 2010, 221, 2938–2948. [Google Scholar] [CrossRef]
  110. Doughty, C.E.; Metcalfe, D.B.; Girardin, C.A.J.; Amezquita, F.F.; Cabrera, D.G.; Huasco, W.H.; Silva-Espejo, J.E.; Araujo-Murakami, A.; da Costa, M.C.; Rocha, W.; et al. Drought impact on forest carbon dynamics and fluxes in Amazonia. Nature 2015, 519, 78–82. [Google Scholar] [CrossRef] [Green Version]
  111. Fu, W.W.; Randerson, J.T.; Moore, J.K. Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in the CMIP5 models. Biogeosciences 2016, 13, 5151–5170. [Google Scholar] [CrossRef] [Green Version]
  112. Zheng, D.L.; Prince, S.; Wright, R. Terrestrial net primary production estimates for 0.5 degrees grid cells from field observations-a contribution to global biogeochemical modeling. Glob. Chang. Biol. 2003, 9, 46–64. [Google Scholar] [CrossRef]
  113. Zhao, M.S.; Heinsch, F.A.; Nemani, R.R.; Running, S.W. Improvements of the MODIS terrestrial gross and net primary production global data set. Remote Sens. Environ. 2005, 95, 164–176. [Google Scholar] [CrossRef]
Figure 1. Map showing the location of the sandy areas in northern China.
Figure 1. Map showing the location of the sandy areas in northern China.
Atmosphere 11 00071 g001
Figure 2. Framework of the Carnegie-Ames-Stanford approach (CASA) estimation model.
Figure 2. Framework of the Carnegie-Ames-Stanford approach (CASA) estimation model.
Atmosphere 11 00071 g002
Figure 3. Simulation and validation of the CASA model; (a) represents the CASA simulation of the multi-annual change in mean net primary productivity (NPP); (b) represents a multi-annual NPP variation chart based on MODIS data (http://luna.ntsg.umt.edu/data); (c) is based on the CASA model simulation of different climate models of the correlation between the model data and MODIS data; and (d) is the Taylor diagram between the field monitoring data and the simulated data based on the flux observation site from CASA model.
Figure 3. Simulation and validation of the CASA model; (a) represents the CASA simulation of the multi-annual change in mean net primary productivity (NPP); (b) represents a multi-annual NPP variation chart based on MODIS data (http://luna.ntsg.umt.edu/data); (c) is based on the CASA model simulation of different climate models of the correlation between the model data and MODIS data; and (d) is the Taylor diagram between the field monitoring data and the simulated data based on the flux observation site from CASA model.
Atmosphere 11 00071 g003
Figure 4. Temporal variation in NPP changes under GW_1.5 °C_2.0 °C levels from 1980 to 2100 and the difference in NPP between warming levels for the whole year, spring (March, April, and May), summer (June, July, and August), autumn (September, October and November), and winter (December, January, and February). The horizonal lines in each boxplot represent the minimum, 25th percentile, median, 75th percentile and maximum NPP, respectively, and the stars represent the average NPP. All results were derived from CASA model simulations.
Figure 4. Temporal variation in NPP changes under GW_1.5 °C_2.0 °C levels from 1980 to 2100 and the difference in NPP between warming levels for the whole year, spring (March, April, and May), summer (June, July, and August), autumn (September, October and November), and winter (December, January, and February). The horizonal lines in each boxplot represent the minimum, 25th percentile, median, 75th percentile and maximum NPP, respectively, and the stars represent the average NPP. All results were derived from CASA model simulations.
Atmosphere 11 00071 g004
Figure 5. Spatial distribution of NPP changes under GW_1.5 °C for the whole year, spring (March, April, May), summer (June, July, August), autumn (September, October, November), and winter (December, January, February). Each letter represents a different climate model; G represents the GFDL-ESM2M model, H represents the HadGEM2-ES model, I represents the IPSL-CM5A-LR model, and M represents the MIROC5 model.
Figure 5. Spatial distribution of NPP changes under GW_1.5 °C for the whole year, spring (March, April, May), summer (June, July, August), autumn (September, October, November), and winter (December, January, February). Each letter represents a different climate model; G represents the GFDL-ESM2M model, H represents the HadGEM2-ES model, I represents the IPSL-CM5A-LR model, and M represents the MIROC5 model.
Atmosphere 11 00071 g005
Figure 6. Spatial distribution of the changes in NPP under GW_2.0 °C for annual and seasonal projections. Notation used for climate models and definition of seasons—as in caption to Figure 5.
Figure 6. Spatial distribution of the changes in NPP under GW_2.0 °C for annual and seasonal projections. Notation used for climate models and definition of seasons—as in caption to Figure 5.
Atmosphere 11 00071 g006
Figure 7. Spatial variation of precipitation, solar radiation and warming level in SAONC. Each letter represents a different climate model; G represents the GFDL-ESM2M model, H represents the HadGEM2-ES model, I represents the IPSL-CM5A-LR model, and M represents the MIROC5 model.
Figure 7. Spatial variation of precipitation, solar radiation and warming level in SAONC. Each letter represents a different climate model; G represents the GFDL-ESM2M model, H represents the HadGEM2-ES model, I represents the IPSL-CM5A-LR model, and M represents the MIROC5 model.
Atmosphere 11 00071 g007
Figure 8. The percentages of areas where the temperature increase is below the global values of 1.5 and 2.0 °C (B1.5_2.0_area, Units: %).
Figure 8. The percentages of areas where the temperature increase is below the global values of 1.5 and 2.0 °C (B1.5_2.0_area, Units: %).
Atmosphere 11 00071 g008
Figure 9. Analysis of uncertainty based on the CASA model. (a) the uncertainty based on observation site; (b) the results of this study were compared with other studies; (c) uncertainty of model parameters.
Figure 9. Analysis of uncertainty based on the CASA model. (a) the uncertainty based on observation site; (b) the results of this study were compared with other studies; (c) uncertainty of model parameters.
Atmosphere 11 00071 g009
Figure 10. Uncertainty in the estimation of NPP based on CASA model in four climate models.
Figure 10. Uncertainty in the estimation of NPP based on CASA model in four climate models.
Atmosphere 11 00071 g010
Table 1. Information on the five flux sites in the sandy areas of northern China (SAONC).
Table 1. Information on the five flux sites in the sandy areas of northern China (SAONC).
Site NameLat (° N)Lon (° E)Precipitation (mm)Height (m)Bioclimate RegionReference
Naiman42.93120.7366.4361Semi-humidZheng et al. [58]
Yanchi37.4107.12295442Semi-aridFu et al. [59]
Shapotou37.53105.8186.61227AridZheng et al. [58]
Fukang42.2887.92160482AridYu et al. [52]
Tazhong38.9783.6522.81082Extremely aridYang et al. [60]
Table 2. List of the four ISI-MIP 2b models used in the analysis.
Table 2. List of the four ISI-MIP 2b models used in the analysis.
ModelInstitutionCountryResolution (Lon × Lat)
GFDL-ESM2MGeophysical Fluid Dynamics LaboratoryUSA144 × 90
HadGEM2-ESMet Office Hadley CenterUK145 × 192
IPSL-CM5A-LRL’ Institute Pierre-Simon LaplaceFrance96 × 96
MIROC5Model for Interdisciplinary Research on ClimateJapan256 × 128
Table 3. Correlation of precipitation, temperature, solar radiation, and NPP.
Table 3. Correlation of precipitation, temperature, solar radiation, and NPP.
Emission ScenariosISI-MIP 2bPrecipitationTemperatureSolar Radiation
RCP2.6GFDL-ESM2M0.198 *0.672 **0.481 **
HadGEM2-ES0.505 **0.812 **0.493 **
IPSL-CM5A-LR0.494 **0.787 **0.425 **
MIROC50.555 **0.847 **0.647 **
RCP4.5GFDL-ESM2M0.355 **0.853 **0.215 *
HadGEM2-ES0.500 **0.882 **0.356 **
IPSL-CM5A-LR0.453 **0.883 **0.512 **
MIROC50.537 **0.88 3**0.352 **
RCP6.0GFDL-ESM2M0.0170.642 **−0.239
HadGEM2-ES0.2020.702 **−0.192
IPSL-CM5A-LR0.0720.742 **−0.325 *
MIROC50.0170.614 **−0.490 **
RCP8.5GFDL-ESM2M0.390 **0.863 **0.062
HadGEM2-ES0.674 **0.850 **−0.223 *
IPSL-CM5A-LR0.598 **0.850 **0.213 *
MIROC50.752 **0.870 **−0.087
* at the 0.05 level; ** at the 0.01 level.
Table 4. Summary of studies of NPP in SAONC.
Table 4. Summary of studies of NPP in SAONC.
MethodsStudy PeriodsNPP Ranges (g C/m2)Reference
LPJ model1961–1970686.75Sun and Mu [80]
CASA model2000–2012556.29Li and Pan [81]
CASA model1982–1999510.72Piao and Fang [82]
LUE model1990584.75Chen et al. [83]
CEVSE model1981–1998420.49Cao et al. [84]
CASA model1982–1999523.59Fang et al. [76]
BEPS model2001558.29Feng et al. [85]
C-Fix model2003654.17Chen et al. [86]
CASA model1989–1993483.19Zhu et al. [21]
GLO-PEM model1981–2000710.49Gao and Liu [30]
CEVSA model1980–2000667.45Gao and Liu [30]
GEOPRO model2000683.29Gao and Liu [30]
GEO-LUE model2000–2004724.54Gao and Liu [30]
LPJ model1961–2080564.42Zhao and Wu [87]
M-SDGVM1981–2000537.24Mao et al. [88]
CASA model1981–2008487.69Chen et al. [89]
BEPS model2000–2010684.29Liu et al. [90]
CASA model2000–2010526.47Pei et al. [91]
CASA model1982–2010545.29Liang et al. [24]
CSCS model2030 2050 2070843.27Gang et al. [79]
CASA model1980–2100604.52This study

Share and Cite

MDPI and ACS Style

Ma, X.; Huo, T.; Zhao, C.; Yan, W.; Zhang, X. Projection of Net Primary Productivity under Global Warming Scenarios of 1.5 °C and 2.0 °C in Northern China Sandy Areas. Atmosphere 2020, 11, 71. https://doi.org/10.3390/atmos11010071

AMA Style

Ma X, Huo T, Zhao C, Yan W, Zhang X. Projection of Net Primary Productivity under Global Warming Scenarios of 1.5 °C and 2.0 °C in Northern China Sandy Areas. Atmosphere. 2020; 11(1):71. https://doi.org/10.3390/atmos11010071

Chicago/Turabian Style

Ma, Xiaofei, Tianci Huo, Chengyi Zhao, Wei Yan, and Xun Zhang. 2020. "Projection of Net Primary Productivity under Global Warming Scenarios of 1.5 °C and 2.0 °C in Northern China Sandy Areas" Atmosphere 11, no. 1: 71. https://doi.org/10.3390/atmos11010071

APA Style

Ma, X., Huo, T., Zhao, C., Yan, W., & Zhang, X. (2020). Projection of Net Primary Productivity under Global Warming Scenarios of 1.5 °C and 2.0 °C in Northern China Sandy Areas. Atmosphere, 11(1), 71. https://doi.org/10.3390/atmos11010071

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop