Elsevier

Ecological Indicators

Volume 154, October 2023, 110610
Ecological Indicators

Quantitative analysis of spatiotemporal changes and driving forces of vegetation net primary productivity (NPP) in the Qimeng region of Inner Mongolia

环境科学与生态学TOPEI检索SCI升级版 环境科学与生态学2区
https://doi.org/10.1016/j.ecolind.2023.110610 Get rights and content
Under a Creative Commons license
open access

Highlights

  • This paper aims at quantifying the factors affecting vegetation growth in the western of Inner Mongolia, and judging the relative contribution of natural factors and man-made factors.
  • It is obtained through geodetectors (GD) that Rainfall is the most important driving factor, with the q-value of 0.77 and GDP also having a huge impact on vegetation growth, q-value of 0.64.
  • Different from most studies, human activity also has a positive effect on vegetation growth.
  • The relative contribution of human activities reached 51.75%, slightly exceeding the influence of natural factors.

Abstract

Vegetation is an essential component of terrestrial ecosystems, and understanding the drivers of vegetation change is of great importance for ecological management. In recent years, vegetation growth has increased under the combined effect of global warming and human activities in Inner Mongolia. The net primary productivity (NPP) was used as an indicator to study the spatial and temporal changes in vegetation in the Qimeng Region (QR). The residual trend analysis method was used to analyze the relative contributions of climate variations (CV) and human activities (HA) to NPP changes across the QR, while their drivers were explored using a geographical detector approach to quantify the driving forces of NPP. The results show that (1) NPP exhibited a fluctuating growth trend from 2003 to 2020, with an overall growth rate of 2.91%/year. (2) Precipitation, GDP and population density were the dominant driving factors for the spatial distribution of NPP. The combined explanatory power of any two dominant factors exceeded the power of any dominant individual factor, and the interaction between climate and human factors had a significant effect on NPP. (3) The change in NPP was influenced by the combined effect of HA and CV, accounting for 37.69% of the total area, with the relative contribution of HA being 51.75%. Finally, the relative contribution of human activities was slightly higher than that of climate change, confirming the initial success of the Grain to Green Program as well as ecological conservation projects. This paper provides a scientific basis for the local government to carry out the conversion of cropland to forest.

Keywords

NPP
Driving force
Geographical detectors
Inner Mongolia

1. Introduction

As the main component of the terrestrial ecosystem, vegetation plays a pivotal role in protecting the living environment of humans and stabilizing the global carbon cycle (Frank et al., 2015, Negi et al., 2019, Wang et al., 2021b). To research the impact of global warming and industrialization on vegetation, the net primary productivity (NPP) of vegetation is used as an indicator to investigate variation in plant growth (Chen et al., 2021a). The NPP of vegetation refers to the total organic matter produced by vegetation through photosynthesis in a unit of time obtained by subtracting the organic matter consumed by plant respiration from the accumulated organic matter (Lin et al., 2012, Turner et al., 2006). Vegetation is influenced by climate change and human activities, which can be directly reflected in the vegetation NPP in a short time period. Therefore, NPP can be used to analyze the health status and sustainable development of regional ecosystems.
Vegetation changes are affected by numerous factors, most of which can be attributed to climate variations and human activities (Di Pasquale et al., 2020, Pan et al., 2011, Shi et al., 2021). There is no doubt that climate change, associated with temperature and precipitation changes, will have a very large impact on vegetation NPP (Alkama et al., 2022, Guo et al., 2020, Qu et al., 2020, Sun et al., 2015a). For example, Liang He et al. analyzed the most significant factors affecting the NDVI spatial distribution on the Loess Plateau. Studies have shown that precipitation plays a key role in plant growth and the spatial distribution of vegetation, especially in arid and semiarid regions (He et al., 2022). Chenli Liu et al. employed a geographical detector to identify the primary factors affecting vegetation changes, and temperature had the highest q value and explanatory power for the NDVI (Liu et al., 2021b). Vegetation is also affected by human activities, especially government policies. Zenghui Sun's research on deserts proves that the project of returning farmland to forests has been a major and vital influence in restoring deserts (Sun et al., 2021). He Liang and Guo Jianbin's studies on the Loess Plateau also demonstrate that human activities have a key influence on vegetation distribution (Asefa et al., 2020). Finally, other factors also affect vegetation distribution, such as slope, elevation, aspect, and sunshine hours, and the importance of different factors will also be significantly variable among different regions (Liu et al., 2021b, Yang et al., 2020).
According to the above description, we know that vegetation changes will be affected by many factors, divided mainly into two categories: climate change and human activities. Therefore, there are two questions to address. The first question is how to determine the impact of various factors on NPP and select the most influential factors. In previous studies, linear correlation analysis methods were usually used to qualitatively analyze the driving factors of vegetation cover (Liu et al., 2018, Whetton et al., 2017, Zhang et al., 2018b). However, there are many factors that affect NPP, and they are not completely independent. Therefore, the results obtained by simple linear or correlation analysis methods cannot fully reflect what is actually driving the relationship. The geographic detector proposed by Wang et al. can overcome the shortcomings of traditional correlation analysis methods [38]. It can quantify the driving effect of each element and the interaction between each factor and geographic elements, which is based on the relationship of driving factors and target variables in space. It is already used to detect the driving factors affecting vegetation change, such as precipitation, temperature and land use type. The second question is how to evaluate the roles of natural changes and human activities in NPP changes. Some scholars use regression models to establish the relationships between NPP and its drivers and then evaluate the relative importance of each factor (Xin et al., 2008). Others use residual trend analysis methods, such as Yanling Sun's research on the relative role of climate change and human activities in vegetation change in North China from 1998 to 2001 (Sun et al., 2015b). The relative contribution of human activities is 69.18%, and the relative contribution of natural factors is 30.82% (Wu et al., 2015). The direct relationship between natural change and human activities is complex. Both natural factors and human activities represent numerous driving factors. It is inaccurate to use only linear relationships to determine the overall impact of these factors. In this paper, nonlinear residual trend analysis is used to more accurately assess the relative contributions of human activities and climate variations.
In the past few decades, the Qimeng region has undergone a transition from ecological degradation to ecological stability, making it an appropriate study area for investigating the driving factors of ecological changes. Before 2000, an increase in population and the modernization of lifestyles resulted in more energy being consumed, along with more ecological pressure (Wang et al., 2017, Wu et al., 2015). Land degradation due to overgrazing is an increasingly serious problem. After 2000, to solve the problem of land degradation, a series of policies were introduced. Typical examples are the Three-North Shelter Forest and programs to return farmland to forest (Du et al., 2016, Yu et al., 2021). To date, the current ecosystem has been significantly restored. Previous studies on vegetation in Inner Mongolia often covered whole provinces or even the Mongolian Plateau (Yang et al., 2019). The purposes of this paper are to (1) illustrate the spatiotemporal variation in NPP from 2003 to 2020, (2) assess the influence of natural and anthropogenic factors and their interactions on NPP spatial distribution, and (3) quantitatively measure the relative contributions of human activities and climate variations.

2. Materials and methods

2.1. Study area

The Qimeng Region (QR, 37.40 N-43.38 N 97.17 E-114.80E) (Fig. 1.) is in southwestern Inner Mongolia. Inner Mongolia is a province with a considerable longitudinal span. The administrative region includes the seven cities of Hohhot, Baotou, Ordos, Ulanqab, Bayannaoer, Wuhai and Alxa League, with a total area of 524,933 square kilometers. It is composed of 59.06% desert and 40.94 grasslands, belonging to arid and semiarid regions (Lu, 2018). It is a crucial area of the ecological protection buffer zone in North China. The topography is higher in the northeastern and middle parts and that of the southwestern part is relatively gentle. The overall altitude ranges from 721 m to 3531 m. The QR is located at a high latitude far from the ocean and is blocked by mountains along the edge. The climate is mainly a temperate continental monsoon climate. The average annual precipitation in the QR is less than 400 mm in most areas, with precipitation gradually decreasing from the southeast to the northwest. Approximately 60% to 80% of the precipitation occurs in the summer, and the distribution of precipitation is also uneven over time. The average yearly temperature of the QR in Inner Mongolia is 6.3–7.6 °C, and the temperature difference between day and night is considerable. This contrasts with the spatial precipitation distribution in that the temperature gradually decreases from northwest to southeast (Mu et al., 2013). Therefore, the QR region is characterized by two main types of belt-shaped grasslands: typical grasslands and desert grasslands. Typical grasslands usually grow under semiarid conditions with an annual precipitation of 200 to 400 mm and an average annual temperature of 0 to 8 °C. They are primarily in the northeastern part of the QR. These grasslands are rich in herbivorous animals, such as Mongolian gazelles, wild horses, sheep, and cattle. Desert grasslands, on the other hand, receive an annual precipitation ranging from 150 to 200 mm and have an average annual temperature of 5 to 10 °C. The biomass in these areas is minimal. In even drier desert regions, the characteristic features include sand dunes, sparse vegetation, and extreme temperature fluctuations. The biodiversity is low, with only a few plant and animal species adapted to desert conditions, such as saxaul trees, jerboas, and Bactrian camels. In recent years, a series of ecological protection projects, such as returning farmland to forest and returning grazing land to grassland, has been carried out in the QR, which has significantly improved the ecological environment.
  1. Download: Download high-res image (583KB)
  2. Download: Download full-size image

Fig. 1. The geographical location and elevation of the Qimeng Region (QR).

2.2. Data sources and processing

The net primary productivity (NPP) dataset used comes from MOD17A3HGF (https://lpdaac.usgs.gov) and covers the period 2003 to 2020. It has a spatial resolution of 500 m and a temporal resolution of 1 year. The preprocessing of the NPP dataset was implemented through Google Earth Engine by the JavaScript application programming interface. Google's platform provides a wealth of free satellite data and products to help researchers perform better research. The elevation data were obtained from the SRTM dataset with a resolution of 30 m, and the slope and aspect were obtained using the functions of Google Earth Engine (https://cmr.earthdata.nasa.gov). The landcover data were derived from the MC- D12Q106 dataset, with a 500 m spatial resolution (Rabus et al., 2003).
The annual average precipitation and temperature data of the QR with a resolution of 1 km were selected from the Resource and Environmental Sciences Data Center of the Chinese Academy of Sciences (https://www.resdc.cn/). Night light data were obtained from the global NPP-VIIRS-like nighttime light dataset from the Harvard Dataverse platform (10.7910/ DVN/YGIVCD). This dataset is based on the DMSP-OLS and NPP-VIIRS nighttime light data correction scheme and is the first “NPP-VIIRS-like” nighttime light dataset with a resolution of 500 m in the world (Chen et al., 2021c). Soil type data, economic data, landform data and other data were selected from the National Earth System Science Data Center (https://www.geodata.cn). A total of 38,000 sample points were extracted from the same grid interval, such as extracting one sample point every 20 grids by using MATLAB. More sample points can be obtained faster by using this approach than by using fishing net tools. The raw precipitation and temperature data had a temporal resolution of one month. Projection transformations, synthesis with the average, resampling, and other processing steps were applied to the raw precipitation and temperature data using ArcGIS to ensure that the output results had the same spatiotemporal resolution. NPP is affected by the main climate and human activity factors listed in Table 1.

Table 1. Indicators of climate and human factors. (The symbol/indicates the absence of temporal resolution).

FactorFactor typeAbbreviationOriginal
temporal resolution
Original
spatial resolution
Source
TemperatureClimate FactorTempMonthly1 kmRESDC
(https://www.resdc.cn/)
PrecipitationClimate FactorPreMonthly1 km
Landform typeClimate FactorLandf/1 km
AspectClimate FactorAspect/30 mNASA/CGIAR
(https://srtm.csi.cgiar.org/)
ElevationClimate FactorDem/30 m
SlopeClimate FactorSlope/30 m
Solar radiationClimate FactorSolarMonthly27.83 kmECMWF (https://cds.climate.copernicus.eu/)
Soil temperatureClimate FactorS_tempMonthly27.83 km
Soil moistureClimate FactorS_moiMonthly27.83 km
Soil typeClimate FactorSoil_t/1 kmNational Earth System Science Data Center (https://www.geodata.cn)
Population densityHuman FactorQmpYearly1 km
GDPHuman FactorGDPYearly1 km
Nighttime lightHuman FactorNltYearly500 mHarvard Dataverse
(https://doi.org/10.7910/DVN/YGIVCD)
Climate factors refer to the various elements or variables that influence the overall climate patterns and conditions in a specific region or on a global scale. In this article, this term mainly refers to the factors that have a significant impact on vegetation growth in the Qimeng city area. These factors operate without direct human intervention or influence. In addition to temperature and precipitation, aspect, elevation and slope were also classified as climate factors. Landform type, solar radiation, soil temperature, soil moisture, and soil type also meet the definition of climate factors. Therefore, this article categorizes these factors as climate factors. Human activity factors refer to the elements directly influenced by human production and daily life. Population density, GDP, and nighttime lights are undisputedly classified as human activity factors.

2.3. Methods

2.3.1. Theil–Sen median trend analysis

The Theil–Sen median method, also known as Sen’s slope estimation, is a robust nonparametric statistical trend method. The linear regression method is susceptible to data noise interference and therefore requires a normal data distribution (Gocic and Trajkovic, 2013, Sharma and Saha, 2017). The Theil–Sen median method is less affected by data noise interference than the other methods. At the same time, it is not clear what the distribution and fluctuations in the annual average NPP are, so it is very suitable to use this method to calculate the trend of NPP (Andaryani et al., 2023, Huang and Kong, 2016). The Sen slope is determined by calculating the permutation and combination for each time point in the time series, where the resulting median of these slopes is the Sen slope. This paper adopts the pixel-by-pixel method to apply the Sen slope, and the calculation formula is as follows:(1)SNPP=Median(NPPj-NPPij-i)2003ij2020
In the formula, NPPi and NPPj represent the i and j years of the annual average NPP value corresponding to the pixel, respectively. When SNPP is greater than 0, the NPP trends upward during the study period, and the environment is more suitable for vegetation. When it is less than 0, the NPP trends downward during the study period.

2.3.2. Mann–Kendall method

Mann–Kendall's test method is a nonparametric method. The significance of the time series trend is determined based on statistical theory, and the data do not need to conform to a normal distribution (Rusciano et al., 2022, Tosunoglu and Kisi, 2017, Zhang and Zwiers, 2004). Mann–Kendall's test statistic was calculated with the following equations:(2)Z=S-1Var(S),S>00,S=0S+1Var(S),S<0(3)S=i=1n-1j=i+1nsgn(NPPj-NPPi)(4)sgn(NPPj-NPPi)=1·ifNPPj-NPPi>00·ifNPPj-NPPi=0-1·ifNPPj-NPPi>0(5)VarS=nn-12n+5-i=1mtiti-12ti+518
In formulas (2) and (3), NPPi and NPPj represent the i and j years for the annual average NPP value corresponding to the pixel, respectively. n represents the length of the time series, n = 18 in this paper, m is the number of tied groups, sgn is a symbolic function, and the value range of the statistic Z is (-∞, +∞). The studied time series can be considered to have changed significantly when we select α=0.05 and Z>u1-α2 to reject the null hypothesis. The results are divided into two categories: significant changes (z>1.96) and slight changes (z<1.96).

2.3.3. Geographical detector

Geographical detectors (Geodetectors) are a set of statistical methods used to detect spatial variability and reveal the driving forces behind it (Liu et al., 2017a, Su et al., 2011, Zhang et al., 2019). In essence, they compare the local variance and the overall variance. If the study area is divided according to a certain driving factor, the local variance of the divided area is calculated. If the value of the dependent variable NPP is close to or similar to the divided area variance, the local variance is much smaller than the global variance. This means that the driving factor has a strong influence on, or correlation with, the NPP. Accordingly, the q value is measured by the relationship between the factor and the analysis variable. This paper focuses on the use of factor detectors and interaction detectors. The factor detector is used to find the most critical driving factors that change the spatial differentiation of NPP, and the calculation formula is as follows:(6)q=1-1nσH2i=1mniσi2
In the formula, the q value indicates the explanatory power of each driving factor for the spatial distribution of NPP, and a larger q value indicates a larger impact of the driving factor on NPP. n refers to all samples in the study area, σH2 refers to the overall variance of the dependent variable NPP in the entire study area, ni is the number of samples in the zone (category) of driving factor X, m is the number of small segmented regions according to a certain driving factor, and σi2 is the local variance in the small segmented region. The interaction detector is used to identify the impact on the dependent variable NPP of the interaction between different driving factors and whether the combined effect of the driving factors will enhance or weaken the explanatory power of the driving factors. The ecological detector is used to assess the significant impact of factors X1 and X2 on the spatial distribution of NPP, and the F test is employed to measure this impact. The corresponding formula is as follows:(7)F=NX1(NX2-1)SSWX1NX2(NX1-1)SSWX2(8)SSWX1=h=1L1Nhσh2SSWX2=h=1L2Nhσh2
NX1 and NX2 represent the sample sizes of two factors, X1 and X2, respectively. The numbers of categories formed by factors X1 and X2 are labelled L1 and L2. SSWX1 and SSWX2 represent the sum of within-group variances formed by X1 and X2, respectively. The null hypothesis assumes that the two factors have the same sum of within-group variances, indicating similar effects on the spatial distribution of NPP. If the null hypothesis is rejected at the significance level α, a significant difference exists between the effects of the two factors on NPP distribution.

2.3.4. Relative contributions of driving forces

The residual trend analysis method was first proposed by Evans and Geerken to reflect the relative effects of climate change and human activities on vegetation (Wang et al., 2021a). In previous studies, rainfall, temperature, and solar radiation were used as key natural factors, and multiple linear regression was used to predict these factors pixel by pixel (He et al., 2022). That is, with precipitation, temperature, and solar radiation as independent variables and the NDVI as the dependent variable, the NDVI under the influence of precipitation, temperature, and solar radiation NDVIcli is calculated through multiple linear regression. The predicted NDVIcli is regarded as the NDVI value affected only by natural factors. In this paper, geographic detectors are used to select the natural factors that have the greatest impact on NPP, including precipitation, soil temperature, and soil type. These factors and the dependent variable NPP do not have simple linear relationships, so traditional linear fitting is replaced by nonlinear fitting to obtain better prediction results. To more accurately describe the influence of natural factors on NPP, which is expressed by NPPcli, the corresponding Equation (9). The corresponding formula for linear fitting is given by Equation (10).(9)NPPcli=a1×pre+a2×pre2+a3×pre3+b1×temps+b2×temps2+b3×temps3+c1×ST+c2×ST2+c3×ST3+d1×pre×temps+d2×temps×ST+d3×pre×ST+d4×pre×temps×ST+e(10)NPPcli=a×pre+b×temps+c×ST+d
NPPcli is the predicted NPP value, Pre refers to the average annual precipitation, temps is the soil temperature, and ST refers to soil type. The terms a, b, c, d, and e represent the multivariate nonlinear regression coefficients. Then, the observed value of NPP is subtracted from the predicted value to obtain the residual.(11)NPPres=NPPobs-NPPcli
NPPres indicates the impact of human activities on NPP, NPPobs is the observed value of NPP, and NPPcli is the predicted NPP value, representing the response of NPP to climate change. This article studies NPP (net primary productivity) changes from 2003 to 2020. Therefore, the format of each grid corresponding to NPPcli is 18*1, where 18 represents the 18 years from 2003 to 2020, and 1 refers to the fitted value for each year. The corresponding sloperes and slopecli are calculated using a simple linear regression method for the yearly NPPres and NPPcli. If the calculated slopes are all greater than 0, both human activities and natural factors have promoted the growth of plants. If the slopes are less than 0, these factors have hindered the growth of plants. This paper uses the p test to detect the significance of sloperes and slopecli. The relative contribution of the driving forces is calculated according to the method proposed by Sun et al. Table 2 provides specific methods. In Table 2, CV is used to represent climate variations, and HA is used to represent human activities. Based on the table, an objective assessment can be made to determine whether the main driving force is climate variations or human activities.

Table 2. Residual trend analysis formula.

ScenarioslopeclisloperesRelative role
of CV/%
Relative role
of HA/%
Driving force
slopeNPP>0>0>0slopecli/slopeNPPsloperes/slopeNPPCV&HA
>0<01000CV
>0>00100HA
slopeNPP<0<0<0slopecli/slopeNPPsloperes/slopeNPPCV&HA
<0>01000CA
>0<00100HA

2.3.5. Evaluation indicators

The coefficient of determination (R2) is commonly used to evaluate the goodness of fit of a model. Its formula is shown below:(12)R2=1-yi-y^i2yi-y¯2
In this context, yi refers to the observed value at time i, y^i represents the simulated value at time i, and y¯ denotes the average of the observed values. This metric is used to evaluate the effectiveness of nonlinear fitting. In the formula for residual trend analysis, NPPcli represents the results obtained by simulating the fitting of NPP values from 2003 to 2020 at each grid using a nonlinear approach. In this study, we utilize R2 to explain the results of the fitting process. The R2 value ranges from 0 to 1, where a value closer to 1 indicates a higher level of goodness of fit, implying that the variability of the dependent variable is better explained by the independent variables.

3. Results

3.1. Temporal-spatial variation characteristics of plant NPP

The average value of vegetation NPP showed an increasing trend from 2003 to 2020 in the QR (Fig. 2). The multiyear average NPP in the study area was 167.31 gC·m-2·a-1. Among them, the annual mean value was the lowest in 2005, only 131.17 gC·m-2·a-1. The NPP annual average sharply rose to 195.55 gC·m-2·a-1 in 2012 and then fluctuated at approximately 175.48 gC·m-2·a-1. The annual average NPP was generally lower than the multiyear average NPP (167.31) in the study area until 2011. In the QC, the annual average NPP was higher than the multiyear average (167.31), except in 2015.
  1. Download: Download high-res image (142KB)
  2. Download: Download full-size image

Fig. 2. Temporal variation in the annual average plant NPP.

Fig. 3 shows that the average value of NPP for many years varies greatly in its spatial distribution. Snpp represents the Sen slope of NPP, and |z| is an indicator for statistical tests. Most of the Alxa League area is desert, and its corresponding NPP is extremely low. Overall, the average NPP gradually increases from west to east. The annual average precipitation also increases gradually from west to east, which is consistent with the spatial distribution of NPP. High-NPP value areas with NPP over 400 gC·m−2·a−1 were distributed in the south and east of the QC area. The areas that passed the MK significance test were mainly in Ordos city in the southwestern portion of the study area, with an average altitude of 3500 m and a climate characterized by low temperature and rainfall.
  1. Download: Download high-res image (359KB)
  2. Download: Download full-size image

Fig. 3. (a) The spatial distribution of the NPP variation trend in the Qimeng Region from 2003 to 2020. (b) Spatial distribution of NPP in 2020 in the Qimeng region.

The vegetation in that area is mainly typical warm-season grassland vegetation and afforestation vegetation in Ordos city, so its corresponding NPP range is 100–200 gC·m-2·a-1. Most regions of the study area are arid and semiarid, with large temperature differences between day and night, and the rainfall distribution pattern is uneven in time and space, which leads to the vegetation being mainly grassland, shrubs, and the artificial vegetation created by returning farmland to forest.

3.2. Dominant factors in the spatial distribution of NPP for each period

The factor detector in the geographic detector can quantitatively evaluate the influence of each driving factor on the spatial distribution of NPP through the q value. Fig. 5 shows the effect of various driving factors on the spatial differentiation of NPP in 2005, 2010, 2015, and 2019. Due to the lack of GDP data in 2020, 2019 was selected instead of 2020 for the last year. The final results illustrate that the dominant factors are precipitation, population density, and GDP in the years mentioned above. In 2005, factors such as precipitation, population density, and GDP contributed the most to the spatial differentiation of NPP. The influences of these factors on the spatial differentiation of NPP were of the same order of magnitude, and there was no noticeable difference, but in 2019, soil type replaced rainfall as the main driving factor.
The secondary driving factors of NPP are soil temperature, soil moisture, and soil type, which are related to soil properties. The influence of soil type and soil temperature is prominently higher than that of soil moisture on NPP. The effects of soil type and soil temperature on NPP are basically the same. The q values corresponding to slope, aspect, elevation, landform, nighttime light, and solar radiation range from 0.0034 to 0.1732, which have very little influence on the spatial distribution of NPP.
In 2010 and 2015, precipitation, population density, and GDP were the main driving forces for changes in the spatial distribution of NPP, and the q values of these factors were significantly higher than those of other factors. This impact is obviously different from that of soil temperature and soil type on NPP. The highest q value among the driving factors is observed for precipitation, with a q value of 0.7705. Precipitation is the determining factor for NPP distribution, reflecting that water was insufficient for vegetation in QC in 2010. Water is the most critical factor restricting vegetation growth.

3.3. Detection of striking differences in the main driving factors

Ecological factors were considered in the geographical detectors to ascertain whether there were striking differences in the influence of various driving factors on the spatial distribution of NPP. The only nonsignificant relationship is found between GDP and population (Fig. 4). The effect of annual average precipitation on the spatial distribution of NPP is markedly different from that of other factors in the study area. Factor detection also illustrated that average yearly precipitation was the crucial factor in NPP variation in QC. At the same time, ecological driver detection also proved that the impact of precipitation on NPP is much greater than that of other driving factors. Similarly, the q-values of population density and GDP are very close in factor detection. The similarity of the two effects on NPP is fully demonstrated in ecological detection.
  1. Download: Download high-res image (123KB)
  2. Download: Download full-size image

Fig. 4. Striking differences in the main driving factors.

3.4. Impacts of climate change and human activities

The impact of climate factors on NPP (Fig. 6). The overall trend of change in the global climate is an increase in temperature, a gradual upward shift of the temperate zone, and an increase in rainfall in North China. In the past 18 years, climate change has been very small, and precipitation has increased slightly to 11.53 mm/10 yr-1 in QC. The annual precipitation peaked of 18.972 mm during 2018–2020, and the average annual rainfall was 10.277 mm in 2007, which was the lowest value in the past 18 years. The temperatures were generally mild, with the average temperature dropping rapidly from 2009 to 2011. In 2009–2011, the annual average temperature reached the lowest value of 8.17 °C, and outside the period of 2009–2011, the temperature showed a slight upward trend. In 2016–2018, the annual average temperature reached a maximum of 15.56℃. Soil temperature, a factor that has a pronounced effect on the spatial heterogeneity of NPP, was identified by the geodetector and increased steadily over 18 years with a slope of 0.28 yr-1. The situation where the annual temperature suddenly drops sharply in some years with a difference between the maximum and minimum values of the annual average temperature reaching more than 7℃ occurred once. The soil temperature peaked of 10.77 °C in the past 18 years. In 2010–2012, the annual average soil temperature reached a minimum value of 9.13 °C, and the difference between the maximum and minimum values was 1.64 °C.
  1. Download: Download high-res image (584KB)
  2. Download: Download full-size image

Fig. 5. Changes in explanatory power (q values) in 2005, 2010, 2015, and 2019.

  1. Download: Download high-res image (530KB)
  2. Download: Download full-size image

Fig. 6. The impact of climate variations and human activities on NPP during 2003–2020.

The impact of human activities on NPP. In this paper, the residual trend analysis method is used to study the impact of human activities on NPP driving factors. The increase in population density is more stable, with a slope of 0.26, which corresponds to an increase in the population of the entire study area of 190,000 people per year. In the population density data, the pixel value corresponding to each pixel is the population density of the corresponding area, and the change rate is multiplied by the total number of pixels to obtain the population change in the study area. According to the ecological detector included in the geographic detector, there is no significant difference observed in the impact of population density and GDP on NPP. Fig. 7 shows that from 2003 to 2020, the population continued to increase in the study area, the rate of increase continued to accelerate, and the slope became increasingly steep. At the same time, the environment has not deteriorated due to population growth; in contrast, it has continued to improve. This phenomenon indicates that the government is actively promoting ecological protection and returning farmland to forests, and the construction of the Three-North Shelter Forest Project has achieved remarkable results. The area where NPPF increased due to human activities accounted for 39.377% of the entire study area, as shown in Fig. 6. Among them, the areas with obviously increased NPP are mainly distributed in Baotou city, the central and western parts of Ulanqab and most of Ordos, accounting for 17.362% of the entire study area. The regions with a slight increase in NPP are concentrated mostly in Hohhot and the southeastern part of Ulanqab, accounting for 22.015% of the entire QC. The area with declining NPP accounted for 1.563% of Qimeng city, mainly in the northwestern part of Bayannaoer and the eastern part of Ulanqab.
  1. Download: Download high-res image (343KB)
  2. Download: Download full-size image

Fig. 7. Time series variation in climatic factors and residual values over the study area from 2003 to 2020: (a) precipitation, (b) temperature, (c) soil temperature, and (d) population density.

3.5. Analysis of NPP change drivers

Fig. 8 shows the spatial distribution of the main drivers of NPP, human activities and climate variations in the 2003–2020 period.
  1. Download: Download high-res image (421KB)
  2. Download: Download full-size image

Fig. 8. Spatial distribution of drivers affecting changes in NPP.

Under the combined effect of human activities and the climate variations, 37.695% of the area shows an upwardly trending NPP in the study area, which is mainly distributed in almost all areas of Ordos, Baotou, Hohhot, and Ulanqab. In the abovementioned area, the NPP was distributed in subregions where climate change was the main driving force, and the increasing trend was in the eastern and northernmost parts of Ulanqab and the southern part of Bayannaoer, accounting for 0.927% of the study area. Correspondingly, with human activities as the main driving force, the areas with an increasing trend of NPP are mainly concentrated in the northern part of Ulanqa and the southeastern part of Bayannaoer. Under the combined effect of human activities and the climate variations, 0.422% of the study area shows downwardly trending NPP, and the point-like division is in the southeastern part of Bayannaoer. The areas where NPP has declined mainly due to climate variations are in the southeastern part of Bayannaoer and appear as discrete point-like divisions in the Fig. 8. In particular, the decline in NPP due to human activities is very rare in QC. The corresponding area cannot even be seen on the map, accounting for only 0.214% of the entire study area. In general, the changes in NPP are mainly driven by the combination of human activities and climate variations in Qimeng city. The individual driving effects of human activities and climate change on NPP are almost negligible.

3.6. Relative contribution of drivers to changes in NPP

This study quantitatively analyzes the relative contributions of climate variations and human activities to NPP based on multivariate nonlinear regression and residual analysis. Table 2 presents a detailed description of the calculation method and related parameters. The relative contribution of the different driving factors to NPP is used to analyze the relationships between the original NPP trend, the NPP trend predicted by the environment, and the NPP trend affected by human activities. Then, calculated by the formula, the result is shown in Fig. 9.
  1. Download: Download high-res image (814KB)
  2. Download: Download full-size image

Fig. 9. The relative contributions of climate variations and human activities to NPP.

The relative contribution of climate change to NPP is 48.25%, while the relative contribution of human activities is 51.75%. In Fig. 9. above, the original trend of NPP is specifically divided into two categories, increase and decrease, and the relative contributions of climate variations and human activities to NPP change are calculated. In the case of an increasing NPP trend, climate variations and human activities accounted for 53.452% and 46.548% of the NPP increase, respectively. Similarly, in the case of a declining NPP trend, the relative contributions of climate change and human activities were 48.310% and 51.689%, respectively. In areas where NPP trends upward, the areas with climate change as the dominant driving factor (relative contribution greater than 50%) are mainly distributed in the eastern part of Ordos, most of Hohhot, and the southeastern part of Ulanqab. Correspondingly, the areas dominated by human activities are distributed in the northwestern part of Ulanqab, the middle part of Baotou and the southwestern part of Bayannaoer.

3.7. Comparison between nonlinear fitting and linear fitting

Equation (9) represents the nonlinear fitting formula for NPPcli. We calculated the R2 values between the fitted values of NPP (NPPcli) and the observed values of NPP (NPPobs) on a grid-by-grid basis. The corresponding R2 distribution is shown in Fig. 10. Overall, the R2 values for nonlinear regression are higher than those for linear regression, and the distribution of R2 values is relatively uniform. In terms of spatial distribution, the R2 values for linear regression are unevenly distributed. In regions such as Ordos, Baotou, and the southern part of Ulanqab, the R2 values are all less than 0.2. In the northern part of Ulanqab and the southern part of Bayannur, the corresponding R2 values are represented by green color in Fig. 10(b), ranging from 0.4 to 0.8. The remaining areas with R2 values greater than 0.4 are sporadically distributed in the southeastern part of Wuhai. In contrast, in nonlinear regression, most regions have R2 values greater than 0.4, with only Baotou, the northern part of Ulanqab, and the border area between Bayannur and Ordos having R2 values less than 0.4. In the eastern part of Wuhai and the northernmost part of Bayannur, there are sporadic areas with R2 values exceeding 0.8.
  1. Download: Download high-res image (1MB)
  2. Download: Download full-size image

Fig. 10. (a) The spatial distribution of R2 with a nonlinear regression. (b) The spatial distribution of R2 with a linear regression.

Within the entire study area, three representative regions were selected: Gander Mountain Nature Reserve, Jiufeng Mountain Nature Reserve, and the urban area of Baotou city. Both linear regression and nonlinear regression yielded R2 values greater than 0.6 for the Gander Mountain Nature Reserve, indicating a good fit. However, in the urban area of Baotou city and the Jiufeng Mountain Nature Reserve, the R2 values corresponding to linear regression were significantly lower than those of nonlinear regression, suggesting that the nonlinear regression outperformed the linear regression in these regions.

4. Discussion

Global warming has led to climate variations in various regions in the past few decades. Under the combined effect of climate variations and human activities, the degree of greening in North China has generally increased. The NDVI in North China has shown an upward trend in the past four decades, and the degree of vegetation cover has increased. This study shows that the rising trend of NPP in Qimeng city is consistent with the increasing trend of the NDVI on the.
Loess Plateau during the same period (He et al., 2022). In this paper, geographic detectors and residual trend analysis were used to quantify the importance of each driving factor and its impact on NPP and the relative impact of human activities and climate change on NPP in the Qimeng city area. Below, the impacts and relative contributions of human activities and climate change are discussed.

4.1. Influence and relative contribution of climate variations

Many studies have proven that vegetation cover has increased significantly at the middle and high latitudes of the Northern Hemisphere, including most of the northern part of China (Pei et al., 2018, Xiao and Moody, 2004). The temperature of the Mongolian Plateau is rising faster than the rate of global warming (Guo et al., 2021, Hao et al., 2021). Similarly, the vegetation cover of the Mongolian Plateau is also increasing, and this greening trend is being experienced over almost the entire Mongolian Plateau (Zhao et al., 2015). The climate variations of Qimeng City is harsh with the Alxa Desert and severe land desertification in Inner Mongolia. At the same time, approximately 67% or 195 million acres (78 million hectares) of the total land in Inner Mongolia is classified as pasture, which greatly limits the ecological carrying capacity (Lee & Sohn, 2011). Some scholars have studied the influence of precipitation and temperature on NPP in Inner Mongolia, finding that precipitation is more influential (Gerten et al., 2008). This is very different from the drivers of vegetation change in Guangxi and Tibet, where temperature is the determining driver (Liu et al., 2021b, Liu et al., 2017b). The QC area is a typical arid and semiarid area, and the impact of precipitation on vegetation growth is far greater than that of temperature. Water resources are the key factor limiting vegetation growth. The rainfall in this area cannot meet the needs of vegetation growth, so the spatial distribution of NPP is consistent with the spatial distribution of rainfall and increases with increasing rainfall. When the maximum precipitation was 30.6–36.7 mm, the average value of NPP also reached a maximum of 179.39. Higher temperatures will reduce soil moisture, increase evaporation, and make water resources scarcer, so NPP has a significant negative correlation with temperature. The Qimeng city area mainly has two types of land cover: grassland and desert. According to Haiyan Zhang’s research, soil stability is significantly weakened under the influence of climate variations, which indirectly alters the soil environment and affects vegetation growth in the Qimeng city area of Inner Mongolia (Zhang et al., 2018a). The interaction of precipitation and temperature significantly enhanced the effect of precipitation alone on vegetation and considerably improved the hydrothermal conditions for vegetation growth. Yuling Zhao focused on the change in grassland NPP in Inner Mongolia from 2000 to 2014 and concluded that the NPP of grassland showed an increasing trend, which is consistent with the findings of this research (Zhao et al., 2019). However, the research conducted by Jay Angerer concluded that the NPP in Inner Mongolia showed a slow downward trend, which is different from the conclusion of this study (Jay et al., 2008), and this discrepancy is likely caused by the different time periods selected for the research. In general, the Qimeng city area of Inner Mongolia is arid and water deficient, and the vegetation lacks sufficient water sources, so rainfall has always been the dominant influencing factor. Soil temperature will affect soil water content and soil erosion, which will indirectly affect the growth of vegetation. Different soil types have different physical and chemical properties and water-holding capacities, which play an important role in determining vegetation growth. The relative contribution of climate factors was slightly higher than that of human activities. The rainfall in Qimeng city showed a slight increasing trend, but it had a significant impact on the growth of vegetation, so the effect of the climate variations on NPP was meaningful. In some areas with more precipitation, such as the southeastern part of Ordos, climate variations are the dominant factor.

4.2. The impact and relative contribution of human activities

The rapid economic growth of Qimeng city has had an important impact on vegetation, and the significant development of coal, natural gas and oil resources has also brought great challenges to the ecological environment in Inner Mongolia (Liu et al., 2021a). Relevant studies have shown that with economic development, population growth, and urbanization, the vegetation in the urban center and its adjacent areas will be rapidly degraded (Du et al., 2020, Liu et al., 2021a). In contrast, the vegetation coverage in Qimeng city did not decline rapidly due to the population influx but increased with increasing GDP. The economy's improvement allowed the government to invest more funds toward improving the ecological environment (Chen et al., 2021b). This fact is consistent with Hu's view that economic development and ecological protection can form a virtuous circle (Hu & Xia, 2019). First, because a large amount of wasteland has been developed into cultivated land, the NPP of cultivated land is higher than that of unused land. Compared with grasslands, crops tend to grow well with human irrigation and pest control. Second, economic development has played an essential role in the government's ecological protection projects (Kang et al., 2021, Matsushita et al., 2007). According to Liang He's research, the improvement of the NDVI on the Loess Plateau and the cumulative area of farmland converted to forest in various provinces in the Loess Plateau have a strong correlation with the NDVI, and the main factor affecting the change in the NDVI is human activities (He et al., 2022). There are some differences in this study. Specifically, human activities were found to be an important factor driving NPP changes, slightly more important than environmental factors, with a relative contribution of 51.75%. This is because compared with other provinces on the Loess Plateau, Inner Mongolia has a smaller, less dense population, so the relative contribution of human activities to NPP is lower than that on the Loess Plateau. A similar example is observed in the Sanjiangyuan area. The population density of the Sanjiangyuan area is lower, and the conclusion is that environmental factors are much stronger drivers than human activities (Gao et al., 2021). This result shows that the greater the population density is, the greater the relative impact of human activities on vegetation growth. The spatial distribution of population density and economic development is consistent with NPP, not only in terms of the role of humans in improving ecological health but also because the most suitable growth environments for humans and plants overlap, and the environment where humans live is often accompanied by more time for plant growth.

4.3. The evaluation of nonlinear fitting and linear fitting

Nonlinear fitting requires determining the degree of the polynomial to obtain the best possible results. In this study, R2 was used to evaluate the fitting. In general simulations, a higher R2 value is considered better. However, in this study, this is not the case. NPPcli attempted to calculate NPP values based solely on climate factors. In regions with minimal or no human activities, a higher simulated R2 value is preferred. However, in regions with significant human activities, excessively high R2 values are not desirable. Two representative regions, the Gander Mountain and Jiufeng Mountain Nature Reserves, were selected and are shown in Fig. 10. In these two regions, human activities are negligible. In the Gander Mountain Nature Reserve, the R2 value for nonlinear regression is slightly higher than that for linear regression, while in the Jiufeng Mountain Nature Reserve, the R2 value for linear regression is significantly lower than that for nonlinear regression, even below 0.2 overall. This indicates that linear regression fails to accurately describe the relationship between climate factors and NPP values in these regions. On the other hand, nonlinear regression performs well and better aligns with the actual situation. However, importantly, overfitting can easily occur with nonlinear fitting, especially considering the 18-year duration of this study. In other words, if the degree of the polynomial is too high with too many terms, the R2 values will be high throughout the entire study area. Therefore, we further compare the R2 values between linear and nonlinear regression in the urban area of Baotou city. In urban areas, an excessively high R2 value in nonlinear regression indicates overfitting. Although the R2 value for nonlinear regression is slightly higher than that for linear regression in the urban area, the overall R2 values are below 0.4, and the R2 values are generally low near the urban area. This indicates that nonlinear regression does not suffer from overfitting. Through the comparison of these three representative regions, it can be mathematically demonstrated that nonlinear fitting performs better than linear fitting.

4.4. Effectiveness, limitations, and future directions

There are some limitations in this study. First, the temporal resolution of the NPP dataset used is annual, while the driving factors for NPP can vary in different months within the same region. Second, it does not encompass all the factors that influence vegetation growth changes. We only considered drivers of climate variations and human activities obtained through publicly available datasets, with less consideration given to factors such as built-up land changes and changes in biological populations. Third, in this paper, nonlinear multivariate regression was used in residual trend analysis instead of linear regression, which showed innovation and provided an example of the superiority of nonlinear multivariate regression over linear regression. However, it still needs further validation and recognition by more researchers. In the geographic detector, this study only used parameter optimization algorithms for the discretization of continuous variables without employing natural classification methods. Different discretization methods can yield different detection results; hence, these differences should be explored in future research. There should be greater utilization of the online calculation and analysis capabilities of the Google Earth Engine (GEE) platform, which allows for studying larger areas using data with high spatial and temporal resolutions without considering local device performance issues (Gorelick et al., 2017). Regarding the consideration of human factors, population and GDP were selected, but the driving factors of animal husbandry were not included. The study area of this paper primarily focuses on the southwestern part of Inner Mongolia, which is predominantly characterized by typical grasslands and desert grasslands. The impact of animal husbandry on this study area is limited, but if the study area is expanded to encompass all of Inner Mongolia or even China, then animal husbandry should be included as a driving factor.
The idea of nonlinear residual regression is very interesting. Simulating the value of NPP should consider the influence of only climate variables using polynomial models, and the results should be compared with the actual NPP to assess the impact of human activities on NPP in the area. If we take nonlinear regression further and directly use machine learning or deep learning to simulate NPP based on climate variables while incorporating more driving factors, we can obtain more accurate predictive results (Ahmad et al., 2023). By including all driving factors in deep learning models, we can establish models to simulate and predict NPP, enabling a precise and scientific assessment of the impact of environmental changes on NPP. This represents a transition from qualitative to quantitative analysis and from simulation to prediction.

5. Conclusion

The Qimeng city area is a key area for sand control and is the site of the Three-North Shelter Forest. From 2003 to 2020, when the climate tended to be humid and rainfall increased, combined with the local government's active efforts to convert farmland to forests and invest in ecological protection projects, the NPP in the Qimeng city area increased significantly. The increase in NPP in Qimeng city is mainly due to the combined effect of human activities and the natural environment, accounting for 37.695% of the study area, and the impact of human activities is slightly higher than that of the environment at 51.75%. The influence of each factor on NPP in order of importance is annual precipitation > population density > GDP > soil type > soil temperature > soil moisture > temperature > elevation > solar radiation. Among natural factors, precipitation and soil type played an important role in the spatial distribution of NPP, and rainfall was the dominant factor in QC. The interaction of driving factors enhances its explanatory power for the spatial distribution of NPP through linear or nonlinear effects. The combined effect of soil types and precipitation has the highest explanatory power for the spatial distribution of NPP. Annual precipitation and temperature had a combined explanatory power of more than 70%, and the interaction of soil types and precipitation had the highest explanatory power of more than 80%. Soil type and soil temperature had an explanatory power of more than 40% and explained NPP changes well. Among the human activity factors, the population and economy had similar effects on NPP. Population density and GDP had an explanatory power of more than 60%. With the increases in population density and GDP, the average value of NPP is also rising because locals actively implement afforestation and ecological protection projects.
MOD17A3HGF, SRTM dataset, and MCD12Q106 used in this paper can be downloaded from Google Earth Engine. Precipitation and temperature data are selected from the Resource and Environmental Sciences Data Center of the Chinese Academy of Sciences. The detailed URLs of the data sources are described in Section 2.2. (Data sources and processing). All data are from public datasets.

CRediT authorship contribution statement

Huazhu Xue: Conceptualization, Methodology, Validation, Software. Yunpeng Chen: Data curation, Writing – original draft, Visualization, Software. Guotao Dong: Investigation, Writing – review & editing. Jinyu Li: Visualization.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This research was funded by the National Natural Science Foundation of China (grant number 51779099), the Key Scientific and Technological Project of Henan Province (232102320247), and the Fundamental Research Funds for the Universities of Henan Province (NSFRF230631).

Data availability

All data are from public datasets.

References

Cited by (4)