Next Article in Journal
Evaluation of Agricultural Bare Soil Properties Retrieval from Landsat 8, Sentinel-2 and PRISMA Satellite Data
Next Article in Special Issue
The Forest Change Footprint of the Upper Indus Valley, from 1990 to 2020
Previous Article in Journal
The Influence of Climate Change on Forest Fires in Yunnan Province, Southwest China Detected by GRACE Satellites
Previous Article in Special Issue
Assessing the Effects of Time Interpolation of NDVI Composites on Phenology Trend Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Environmental Factors on the Changes in MODIS NPP along DEM in Global Terrestrial Ecosystems over the Last Two Decades

1
State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University, Xining 810016, China
2
Key Laboratory of Restoration Ecology for Cold Regions Laboratory in Qinghai, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining 810008, China
3
Key Laboratory of Adaptation and Evolution of Plateau Biota, Chinese Academy of Sciences, Xining 810008, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(3), 713; https://doi.org/10.3390/rs14030713
Submission received: 16 December 2021 / Revised: 26 January 2022 / Accepted: 31 January 2022 / Published: 2 February 2022

Abstract

:
Global warming has exerted widespread impacts on the terrestrial ecosystem in the past three decades. Vegetation is an important part of the terrestrial ecosystem, and its net primary productivity (NPP) is an important variable in the exchange of materials and energy in the terrestrial ecosystem. However, the effect of climate variation on the spatial pattern of zonal distribution of NPP has remained unclear over the past two decades. Therefore, we analyzed the spatiotemporal patterns and trends of MODIS NPP and environmental factors (temperature, radiation, and soil moisture) derived from three sets of reanalysis data. The moving window method and digital elevation model (DEM) were used to explore their changes along elevation gradients. Finally, we explored the effect of environmental factors on the changes in NPP and its elevation distribution patterns. Results showed that nearly 60% of the global area exhibited an increase in NPP with increasing elevation. Soil moisture has the largest uncertainty either in the spatial pattern or inter-annual variation, while temperature has the smallest uncertainty among the three environmental factors. The uncertainty of environmental factors is also reflected in its impact on the elevation distribution of NPP, and temperature is still the main dominating environmental factor. Our research results imply that the carbon sequestration capability of vegetation is becoming increasingly prominent in high-elevation regions. However, the quantitative evaluation of its carbon sink (source) functions needs further research under global warming.

1. Introduction

Since the nineteenth century, the global near-surface temperature has continued to increase according to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC). The temperature in the last 10 years after 2000 resulted in the hottest 10 years in history, and 1983–2012 was the hottest 30 years in the past 800 years. The widespread impact of global warming has caused a series of negative ecological consequences, such as drought [1,2], melting [3], rising sea levels [4], and frequent extreme climates [5,6]. The accelerating global warming has become a major challenge that is restricting the sustainable development of human society [7,8,9].
Vegetation is an important part of the terrestrial ecosystem and plays a crucial role in sequestering carbon and mitigating climate change [10]. The vegetation ecosystem is found to be more vulnerable and sensitive to climate change than the other ecosystems. As a variable that reflects the efficiency of vegetation fixation and conversion of light energy, net primary productivity (NPP) is widely used in the monitoring of vegetation dynamics [11,12]. It is related to life activities such as vegetation growth, development, and reproduction, and it also provides an indispensable material basis for the life activities of other biological members in the entire ecosystem. Most of the studies focus on the spatiotemporal changes in NPP [13,14] and the impact of climate and human activities on NPP [15,16], modeling organic carbon storage with NPP as input data, and changes in carbon sources and sinks of the terrestrial ecosystem [12,17,18,19]. Studies on the changes in NPP with increasing elevation gradients (EG) are mainly conducted in local regions, and all come from instantaneous surveys [20,21,22,23,24]. Therefore, we still lack the knowledge and a full picture of the changes in NPP along EG at the global scale.
Typically, NPP declines with increasing elevation, which can largely be explained by the decreasing temperature as elevation increases [20,22]. However, the impact of temperature on the NPP of vegetation may be more significant in high-elevation areas with the intensification of global warming [25,26], because the warming rate in high-elevation areas is often greater than that in low-elevation areas [25,27]. Thus far, the warming rate has been intensifying at the global scale for decades [26,28,29,30], which is likely to cause a completely different spatial pattern of NPP on EG [31]. A few studies have focused on determining the spatial pattern of changes in vegetation greenness along EG in recent years [29,32,33], and their results indicated that the signal of the vegetation greenness increases with increasing elevation is found at the global scale and regional scale. However, vegetation greenness is not completely equal to vegetation productivity, and it does not directly participate in the process of the carbon cycle. This situation raises a scientific question as to whether the elevation pattern of NPP has changed under the influence of global warming.
Considering the importance of the effect of environmental variation on vegetation NPP, and the defects and gaps in current scientific research, we analyzed the inter-annual variation of MODIS NPP product (MOD17A2HGF, version 6.1) and environmental factors (air temperature, solar radiation, and soil moisture that derived from the reanalysis data of ERA5, MERRA2, and NCEP2). We calculated the changes in NPP, temperature, radiation, and soil moisture along EG by using the digital elevation model (DEM), which was also used to verify whether the elevation pattern of NPP has been altered under environmental variation. Furthermore, we analyzed the effect of environmental factors on the spatiotemporal changes and elevation distribution of the NPP from 2001 to 2020. We hope that the results of this study can provide references for the evaluation of terrestrial ecosystem carbon source and sink functions and the improvement and development of carbon cycle models. This study will contribute to our understanding of the impact of environmental factors on the elevational distribution of NPP in recent decades and will help us formulate strategies for mitigating climate change. We also expect that the outcomes of this study can provide references for the evaluation of terrestrial ecosystem carbon sink (source) functions and the improvement of current carbon cycle models.

2. Materials and Methods

2.1. Datasets

2.1.1. NPP Data

The NPP data used in this study come from the MODIS (Moderate Resolution Imaging Spectroradiometer) remote sensing product of MOD17A3HGF Version 6.1 [34], which will be generated at the end of each year when the entire yearly 8-day MOD15A2H is available. MOD17A3HGF has two data fields named Npp_500m and Npp_QC_500m, which represent NPP data and their quality control (QC). The poor-quality inputs were cleaned from the 8-day leaf area index and the fraction of photosynthetically active radiation based on the QC label for every pixel. The data type of Npp_QC_500m is an unsigned 1-byte integer (uint8) with a valid range from 0 to 100 (Units = 100%). Finally, we use the QC file to select pixels with good quality to participate in the analysis. The annual NPP, with a spatial resolution of 500 m, is derived from the sum of all 8 days of NPP in a certain year from 2001 to 2020. We further converted the coordinate system to World Geodetic System 1984 and resampled the spatial resolution to 0.008° (≈1 km) by using the bicubic method to match the resolution of the DEM data.

2.1.2. DEM Data

The Shuttle Radar Topography Mission (SRTM) is a joint project between the National Geospatial-Intelligence Agency and the National Aeronautics and Space Administration. The SRTM30_PLUS DEM dataset with a spatial resolution of 30-arc second (≈1 km) was developed by the Scripps Institute of Oceanography, University of California San Diego [35]. The coordinate system is the World Geodetic System 1984. The land data are based on the 1 km averages of topography derived from the United States Geological Survey (USGS) SRTM30 (30 arc-sec) grided DEM data product created with data from the NASA SRTM [36]. Global 30 Arc-Second Elevation (GTOPO30) data are used for high latitudes where SRTM data are not available. GTOPO30 is a global DEM with a horizontal grid spacing of 30 arc seconds (approximately 1 km) [37]. SRTM has been extensively used and provides a good representation of the topography [38], which performs better than Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM) (the spatial resolution is 500 m) in terms of micro-topography, hydrologic-network and structural information characterization [39,40]. Considering the high quality of the SRTM data and the computational efficiency of the data, we chose SRTM30_PLUS as the DEM data in this study.

2.1.3. Reanalysis Data

We selected three sets of reanalysis data (Table 1) to study the spatiotemporal changes of environmental factors (temperature, radiation, and soil moisture), their distribution along EG, and the uncertainties of the correlation between environmental factors and NPP. The three sets of reanalysis data are the fifth generation of European ReAnalysis (ERA5), the recent Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA2), and National Centers for Environmental Prediction Reanalysis 2 (NCEP2). These reanalysis data can completely cover the data period of MODIS NPP and contain three kinds of environmental factor data at the same time. In data processing, we convert the unit of temperature from K to °C by subtracting 273.15, and the unit of radiation was converted from W/m2 to MJ/m2 because the unit of absorbed photosynthetically active radiation is MJ/m2 in the NPP simulation. Therefore, we can intuitively observe the magnitude of the spatial distribution of radiation. The soil moisture has different names in each reanalysis dataset, and we selected the variables with similar meanings and consistent units. Furthermore, we supplemented the monthly data with mean annual temperature (MAT), radiation (MAR), and soil moisture (MASM) raster data. Then, we extracted these environmental data from 2001 to 2020 for analysis. To keep the coordinate system of the reanalysis datasets consistent with the DEM, we set the resolution of the MERRA2 and NCEP2 data to 0.5° and 1.875°, respectively. After that, we upscaled the DEM data to match them, and then performed the moving window operation. We explain the reason for this in Section 2.2. However, trend and correlation analysis were performed at the resolution of 0.008° to keep the data resolution consistent throughout the main text, and also to be consistent with the NPP data.

2.2. Calculation of the Changes in NPP and Environmental Factors along EG

The changes in NPP along EG (NPPEG) were calculated in three steps that are the same as those for the change of environmental factors along EG [29,32,33]. Step 1: We selected a 9 × 9 moving window to traverse the DEM and NPP (or environmental factors) raster data (Figure 1). The size of the moving window determines the amount of data used for analysis for each moving step. A previous study indicated that the difference caused by the window size is extremely small and recommended choosing a window size of 9 × 9 for calculation [29,33]. We also tested the window sizes of 5, 7, 9, and 11, and then expanded the size of the moving window to 19, 29, 39, and 49. We found that the changes in the areas of positive NPPEG were 0.0039% from the window size of 5 to 11, while the changes in the areas were reached 1% from the window size of 19 to 49. Therefore, we followed the recommendation of previous studies and chose a window size of 9 × 9 for the study. Step 2: The data extracted by the moving window have two dimensions. Thus, we need to reduce the dimensions and arrange them in pairs for linear regression analysis. A 9 × 9 moving window contains 81 data after dimensionality reduction. Step 3: We assign the regression coefficient obtained by linear regression analysis to the center pixel of the moving window and traverse the entire raster data accordingly to obtain the distribution map of the change of the variable along EG. The statistical significance of the slope is tested by a t-test. It should be noted that we set the NPP value of the ocean to no data. Therefore, when the moving window contains pixels of the ocean, the center pixel of the moving window will be filled with no data. In addition, we upscaled the DEM data to match the original resolution of environmental factors derived from three reanalysis datasets, such that we can avoid the regression analysis was established between the same environmental factors (coarser grid) and different DEM data (fine grid) in a moving window.

2.3. Trend Analysis

We use the ordinary least squares regression method [15,44,45,46] to determine the variation of NPP and environmental factors in time series:
Slope   =   n × i = 1 n i × ( VAR ) i - i = 1 n i i = 1 n ( VAR ) i n × i = 1 n i 2 - i = 1 n i 2 ,
where VAR can be NPP, MAT, MAR, MASM, and their changes along EG; i is the sequence number of the year (from 2001 to 2020); n represents the total number of years. The significance of the trend was determined by the T-test.

2.4. Correlation Analysis

We use Pearson correlation analysis to calculate the correlation between NPP and environmental factors (MAT, MAR, and MASM). The correlation coefficient (r) of the two variables can be calculated by Equation (2).
r xy = i = 1 n x i   x ¯ y i   y ¯ i = 1 n x i   x ¯ 2 i = 1 n y i   y ¯ 2 ,
The correlation coefficient between NPP and MAT was taken as an example. Variables xi and yi denote the NPP and MAT in year i, respectively; x ¯ and y ¯ are the mean values of NPP and MAT from 2001 to 2020, respectively.

3. Results

3.1. Trends of NPP and Environmental Factors

NPP exhibited an increasing trend (1.75 gC/m2/yr) from 2001 to 2020 (Figure 2a), which was almost below the mean value for the first 10 years, and above the mean value for the next 10 years. The trends of the three kinds of reanalysis data showed that MAT was significantly increased. However, ERA5 was significantly decreased, while both MERRA2 and NCEP2 showed significant increasing trends. By contrast, the trends of MAR exhibited a decreasing trend, but only MERRA2 was statistically significant (Figure 2b–d). The spatial pattern of mean annual NPP is consistent with our perception that tropical forests have the highest NPP, while NPP is relatively low in alpine and arid regions (Figure 3a). The spatial pattern of the NPP trend showed that 70.63% of the global areas (16.21% were statistically significant) presented an increasing trend from 2001 to 2020, which were mainly found in central and western Canada, parts of central and northern China, central and southern South America, and central Africa (Figure 3b). By contrast, the regions with a decreasing trend of NPP occupied 29.37% of the global areas, and 3.02% of them were statistically significant, which were mainly distributed in northern South America.
The uncertainties of the trend of environmental factors were reflected in the spatial distribution (Figure 4). MAT showed a significant increasing trend in high−latitude regions of Asia and Europe, but opposite trends were found in central Africa (Figure 4a,d,g). The consistency of the positive trend of MAT accounted for 56.80% of the global area, whereas the areas with a negative trend accounted for 7.28%, mainly in northeastern North America, the Iranian plateau, and parts of the region across central India (Figure 4 j). Strong spatial heterogeneities are observed in the spatial pattern of MAR trends of the three sets of reanalysis data (Figure 4b,e,h). The inconsistent trends of MAR accounted for 61% of the global area. By contrast, the consistency of negative trends (27.55%) was greater than that of positive trends (11.65%). It is difficult to find areas where the MASM trend is consistent (Figure 4c,f,i), and MASM trends are inconsistent across most regions of the globe (73.84%) (Figure 4l). We found similar spatial patterns of MAT, MAR, and MASM (Figure 5). MAT is lower in high−elevation areas (Figure 5a,d,g). However, the cold climate of the Qinghai–Tibet Plateau is not reflected in NCEP2 MAT data (Figure 5g). Therefore, we speculated that the NCEP2 MAT probably fail to capture the temperature distribution. The higher uncertainty of MAT includes western North America, the Andes in South America, Greenland, the Mongolian Plateau, and the Qinghai–Tibetan Plateau (Figure 5j). MAR in high latitudes is lower than that in other regions (Figure 5b,e,h), which is highly uncertain in Greenland and the Sahara Desert (Figure 5k). MASM has relatively high values in the high latitudes of the northern hemisphere, the eastern United States, southeastern China, and tropical forests (Figure 5c,f,i). By contrast, it is lower in the western United States, the Mongolian plateau, the Sahara Desert, southern Africa, and Australia. The higher uncertainty of MASM is mainly found in the high latitudes of the northern hemisphere, eastern China, and tropical forests (Figure 5l).

3.2. Elevational Distribution of NPP and Environmental Factors

Changes in NPP and its slope along EG (NPPEG and NPP EG slope ) are shown in Figure 6. The positive NPPEG accounted for 59.98% of the global area, and 31.32% of them reached a significance level of 0.05. The positive NPPEG was scattered around the Great Lakes, the eastern foothills of the Andes, the eastern Brazilian plateau, the sub−Saharan African continent, the Indochina Peninsula, and eastern Australia. By contrast, the region with negative NPPEG accounted for 40.02%, and the significantly negative region occupied 11.18%. These regions are mainly found in northeastern North America, southern South America, and Central Asia (Figure 6a). The NPP EG slope showed strong spatial heterogeneity during the study period (Figure 6b), and there is no obvious spatial distribution feature exists. The region with positive NPP EG slope accounted for 53.88%. By contrast, 46.12% of the global area had negative NPP EG slope .
The spatial pattern of MAT, MAR, and MASM along EG (MATEG, MAREG, MASMEG) is shown in Figure 7. MAT EG ERA 5 and MAT EG ERA 2 follow the natural law that temperature decreases with the increase in elevation, and the area of the two reached 84.36% and 83.67% (Figure 7a,d), with the statistically significant area even reaching up to 78.57% and 74.35%, respectively. By contrast, only 32.90% of the area of MAT EG NCEP 2 conforms to the natural law, which indicates that the temperature data of NECP2 are likely to be wrong in elevational distribution (Figure 7g). Therefore, we excluded NCEP2 when calculating the spatial uncertainty of MAT. The results showed that the spatial consistency of MAT EG ERA 5 and MATMERRA2 reached 78.84%. MAREG does not show obvious spatial distribution characteristics (Figure 7b,e,h), and the areas with positive MAREG (59.73%) were larger than those with negative MAREG (40.27%). The spatial consistency of the MAREG is 39.06%, which is mainly distributed in central and eastern North America, most of Europe, and central and western Russia. However, the inconsistent regions still accounted for most of the global area (Figure 7k). The areas with positive MASMEG have strong spatial heterogeneity, although it occupied (64.09 ± 7.02)% of the global area (Figure 7c,f,i), and the spatial inconsistency of MASMEG is up to 62.73% (Figure 7l).
We further investigated the spatial pattern of changes in the slope of MAT, MAR, and MASM along EG ( MAT EG slope , MAR EG slope and MASM EG slope ) (Figure 8). The positive and negative MAT EG slope accounted for 52.03% and 47.97% of the global area, respectively (Figure 8a,d,g). The area with consistent changes of MAT EG slope accounts for 56.28% of the global area, and the remaining regions are highly uncertain (Figure 8j). The areas with positive MAR EG slope occupied (45.50 ± 3.00)% of the global area, which is less than that of negative MAR EG slope (54.50 ± 3.00)% (Figure 8b,e,h). Anyway, MAR EG slope is highly uncertain in most of the global area (68.11%) (Figure 8k). The areas of positive MASM EG slope (51.49%) are slightly more than that of negative MASM EG slope (48.51%) (Figure 8c,f,i). The inconsistent region covered up to 72.09% of the global area (Figure 8l).

3.3. Effect of Environmental Factors on the Elevational Distribution of NPP

Figure 9 illustrates spatial patterns of the dominating environmental factor on the changes in NPP from 2001 to 2020. ERA5 data showed that MAT and MAR seem to be the dominating factors on the changes in NPP in the north of 30°N, central and northern parts of South America, central Africa, and Southeast Asia, whereas MASM mainly affects the changes in NPP in the central and eastern United States, the eastern part of the Brazilian plateau, southern Africa, and Australia (Figure 9a). By contrast, the effects of environmental factors of MERRA2 on the changes in NPP have a clear spatial distribution pattern. MAT mainly affects the changes in NPP in the Qinghai–Tibet Plateau and parts of central South America. MAR is the dominant factor in the high latitudes of the northern hemisphere, southern South America, and Southeast Asia. MASM has a wider range of influence, including southern North America, Eurasia from 30°N to 60°N, the Indian peninsula, most of Africa, and Australia (Figure 9b). The MAT and MAR of NCEP2 are the dominant factors in the changes in NPP in the north of 30°N. MAR is also the main environmental factor that affects NPP changes in northwestern South America, central Africa, and Southeast Asia. The areas where MASM presented a dominant environmental factor include most of the United States, northern South America and most of Brazil, the southern edge of the Sahara Desert and southern Africa, most of Eurasia from 30°N to 60°N, the Indian Peninsula, and Australia (Figure 9c). MASM has the highest spatial consistency (25.15%). These regions are mainly distributed in the United States, eastern and southern South America, southern Africa, Eastern Europe, Central Asia, parts of East Asia, India, and Australia, followed by MAR (9.74%) and MAT (7.70%), which can be found in high latitudes in the northern hemisphere (Figure 9d).
We found that NPP has an obvious turning point at the elevation of 3050 m and divided the elevation into high and low gradients based on the turning point (Figure A1 in Appendix A). Then, we explored the impact of environmental factors on NPP in high− and low−elevation areas (Figure 10a). The three kinds of reanalysis datasets showed that MAT still maintained the highest R2 in the three environmental factors. We further fitted the regression equation curve and uncertainty range of NPP and environmental factors (Figure 10). The effect of MAT on NPP is linear and non−linear at low and high elevations, respectively (Figure 10a). The R2 between MAT and NPP exceeds 0.9 and is statistically significant at the 0.01 level. In general, the uncertainty of the effect of MAT on NPP gradually decreases as the temperature increases, and it decreases more at low elevations. The effect of MAR on NPP is nonlinear at both low and high elevations (Figure 10b). NPP increases as MAR increases at low elevations, whereas opposite trends were observed at high elevations. The uncertainty of the impact of MAR on NPP gradually decreases with the increase in MAR at low elevations. By contrast, the uncertainty does not decrease significantly as the MAR decreases at high elevations (Figure 10c). We found a linear and nonlinear effect of MASM on NPP at low and high elevations, respectively. NPP increases as MASM increases at low elevations. However, the upper and lower limits of uncertainty trend toward the opposite direction. The same situation occurs at high elevations that larger uncertainty of the effect of MASM on NPP, and MASM remains stable with the increase in NPP (Figure 10d). The R2 of all fitted equations is statistically significant at both high and low elevations.

4. Discussion

4.1. Elevational Patterns of Environmental Factors and NPP

In this article, we use the moving window method to explore the changes in NPP and environmental factors along EG at the global scale. Unlike field surveys, this method is based on remote sensing and DEM data, which can quickly and economically monitor the changes of variables along EG. In addition, the moving window method used in this article is intended to investigate the change of NPP and environmental factors at the relative EG within the moving window, which is different from the studies on the absolute EG conducted through a field survey. The accuracy of this method can be indirectly proved by the natural law that the temperature decreases with the increase in elevation. Nearly 80% of the regions showed a decrease in temperature with increasing elevation, and more than 70% of the regions are statistically significant. This result means that the moving window method can effectively monitor the change of the variable in the EG and has high accuracy. Undeniably some regions (although very few) have a local micro-climate, the temperature of which will not decrease with the increase in elevation. We can also evaluate the data quality of temperature based on the natural law of temperature–elevation and the moving window method. For example, NCEP2 fails to capture the elevation distribution of MAT because it is contrary to the natural law that the temperature decreases with the increase in elevation. Unfortunately, we are still unable to assess the data quality of radiation and soil moisture because the spatial distributions of the reanalysis data and their trends are highly uncertain, and there is no corresponding natural law to follow. Typically, NPP decreased with increasing elevation (negative NPPEG) because of the limitation of low temperature. However, we found that the areas with positive NPPEG accounted for 59.98% of the global area, which means the elevation pattern of NPP in our perception has changed. The spatial pattern of NPP slope demonstrated that, compared with NPP in low-elevation areas, NPP in high-elevation areas has a higher increase rate, and this phenomenon has become more common worldwide.

4.2. Uncertainty of Environmental Factors and Their Effect on NPP

We found that MAR and MASM had larger uncertainties than MAT in the inter-annual variations and the spatial distribution along EG. Such a large uncertainty makes it difficult to assess their effect on the changes in NPP and further forms a lower spatial correlation consistency. Factors affecting MAR mainly include cloud height, thickness, shape, and aerosols. [47]. The factors that affect soil moisture include soil texture [48], the data assimilation process [49], the sample size of observation data [50], and land use type [51].
We found that the effect of radiation on NPP is nonlinear regardless of whether it is in high- or low-elevation regions. The positive correlation between NPP and radiation in low-elevation is likely to be that, as the elevation increases, temperature and water inhibit the photosynthesis of plants, thereby reducing the use of photosynthetically active radiation. Moreover, the radiation available for plants decreases, thereby forming a non-linear positive correlation between NPP and radiation. However, the non-linear negative correlation between radiation and NPP in high-elevation regions suggests that radiation is not the main driving force for changes in NPP. We suppose that radiation cannot be fully utilized by plants in high-elevation areas, and temperature is still the main environmental factor that determines changes in NPP along EG. The effect of soil moisture on NPP shows a non-linear positive correlation in low elevation areas, but it is highly uncertain because of the large differences in the trend of soil moisture among ERA5, NCEP2, and MERRA2. Anyway, the positive correlation between NPP and soil moisture has also been confirmed in the arid region [52], whereas soil moisture in high-elevation areas is higher than that in low-elevation areas, it has not caused an increase in NPP. Therefore, we hold the view that soil moisture determines the lower limit of NPP, and the upper limit of NPP is determined by both soil moisture and temperature.
Overall, the temperature is still the dominating climate factor that determines the spatial patterns of NPP along EG, and it is generally positively correlated with NPP whether in high- or low-elevation areas. However, NPP is more sensitive to temperature in low-elevation areas, which indicates that temperatures are more conducive to vegetation photosynthesis in this region. By contrast, the temperature gradually decreases with the increase in high-elevation, and the photosynthesis of vegetation is also limited by low temperature. Thus, NPP appears to rise first and gradually stabilizes. The effect of temperature on NPP along EG has been confirmed in regional-scale studies [53,54].

4.3. Implications for Carbon Cycle

NPP is the source of energy in the carbon cycle of the terrestrial ecosystem. Our research shows that nearly 60% of the global area exhibits an increase in NPP with increasing elevation, which means that vegetation in high-elevation areas plays an increasingly prominent role in absorbing atmospheric CO2 and mitigating climate change. However, this increase in NPP is caused by global climate change, especially the increase in temperature in high mountain areas [33], which improves the photosynthesis capacity of plants. However, when the temperature exceeds the optimal temperature of the plants, it will inhibit the photosynthesis of plants, and even cause the death of local species because they cannot adapt to the rapid warming. This kind of research on vegetation degradation caused by warming has been widely reported [12]. What we want to emphasize is that this phenomenon of increasing NPP with elevation may be beneficial to CO2 fixation in the short term, and an uncontrolled increase in temperature will inevitably lead to vegetation degradation and even ecosystem collapse. The increase in temperature will also cause the increase in plant autotrophic respiration and soil heterotrophic respiration, and the CO2 produced by the respiration process is directly discharged out of the vegetation–soil system. With the high degree of uncertainty in soil respiration, much uncertainty exists in the quantitative evaluation of the carbon source and sink functions of the ecosystem.
Environmental factors have a strong influence on the terrestrial carbon cycle. Temperature is the basic input data to establish its impact on the soil carbon cycle. In this study, the difference in the spatial distribution of ERA5 and MERRA2 temperature is very small. However, the algorithm for the effect of temperature on the decomposition of soil organic carbon (shown as a nonlinear positive correlation effect) has slight differences. The differences in the algorithms are reflected in space, mainly located in high latitude areas and the Qinghai–Tibet Plateau with great uncertainty. This condition means that the uncertainty of the temperature algorithm in the alpine region will greatly affect the decomposition of soil organic carbon in the carbon cycle, which in turn affects the spatial distribution of soil organic carbon storage in the alpine region, leading to distortions in the evaluation of the carbon source and sink functions of these regions.
Numerous studies have shown that soil moisture is crucial to the carbon cycle process of terrestrial ecosystems [55,56,57]. However, the three sets of reanalysis data show that the spatial distribution of soil moisture has great uncertainty. Similarly, there are large discrepancies in the algorithm for the effect of soil moisture on the decomposition of soil organic carbon. The spatial uncertainty of soil moisture involved tropical rain forests and the Sahara Desert [17]. The fundamental reason is that the understanding of the effect of soil moisture on the decomposition of soil organic carbon varies, especially as to whether high-humidity environments will inhibit the decomposition of soil organic matter. An anaerobic environment created by high humidity will inhibit the decomposition of soil organic matter in the carbon cycle model, which will lead to the accumulation of soil organic carbon (high soil organic carbon), and vice versa, leading to the decomposition of soil organic carbon (low soil organic carbon). The uncertainty of the original input data of soil moisture and the algorithm makes it difficult for us to evaluate its effect on the carbon cycle. Considering the importance of soil moisture in the carbon cycle of terrestrial ecosystems, we propose that future studies need to further strengthen the observation of soil moisture and improve its simulation accuracy to more accurately assess the carbon source and sink functions of terrestrial ecosystems.

5. Conclusions

This study analyzed the spatiotemporal changes of the MODIS NPP product and environmental factors (temperature, radiation, and soil moisture data derived from the reanalysis data ERA5, MERRA2, and NCEP2) and their distribution along EG. We also identified the spatial uncertainty of environmental factors and their effects on the elevation distribution of NPP. We found that nearly 60% of the global area presented an increase in NPP with increasing elevation, which implied that the elevation pattern of NPP has changed, and the carbon sequestration capacity of vegetation is increasing elevation. However, soil respiration is likely to increase as well. Quantitatively evaluating the carbon sink (source) function of vegetation remains to be further studied in high-elevation regions. The temperature of NCEP2 failed to capture the alpine environment of the Qinghai–Tibet Plateau, and it does not clearly show the natural law that the temperature decreases with the increase in elevation. Soil moisture has the largest areas of spatial consistency in affecting the spatiotemporal changes in NPP among the three environmental factors. However, its spatial pattern and variation are the most uncertain among the three environmental factors, even though it is essential to the carbon cycle of terrestrial ecosystems. NPP has obvious elevation differentiation with an elevation of 3060 m as the demarcation point, which divides the elevation into low and high. MAT is the main driving force that affects the elevation distribution of NPP, with its effect on NPP exhibiting a significant linear and nonlinear positive correlation at low and high elevations, respectively. The results of this study are expected to contribute to our understanding of the changes in NPP along EG and provide references for the development of terrestrial ecosystem carbon cycle models.

Author Contributions

Conceptualization, Z.W. and X.H.; data curation, Z.W., K.Z. and X.L.; Formal analysis, Z.W.; investigation, Z.W., H.W. and T.W.; methodology, Z.W. and L.W.; writing—original draft, Z.W.; writing—review and editing, Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 32160278); the Open Project of State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University (Grant No. 2020-ZZ-14, 2021-KF-08); the National Natural Science Foundation of China (Grant No. NSFC32001188); the Natural Science Foundation of Qinghai Province of China (Grant No. 2021-ZJ-973Q); the National Natural Science Foundation of China (Grant No. 42067070).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in (NASA EOSDIS Land Processes DAAC) at (doi.org/10.5067/MODIS/MOD17A3HGF.061), reference number [33].

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Changes in NPP and environmental factors along EG from 2001 to 2020. The elevation of 3060 m is used to distinguish high and low elevation (the vertical black solid line in Figure (aj)). The dashed line represents the average value of the variables at high and low elevations. The shaded areas with different colors in each figure represent ±1 SD.
Figure A1. Changes in NPP and environmental factors along EG from 2001 to 2020. The elevation of 3060 m is used to distinguish high and low elevation (the vertical black solid line in Figure (aj)). The dashed line represents the average value of the variables at high and low elevations. The shaded areas with different colors in each figure represent ±1 SD.
Remotesensing 14 00713 g0a1

References

  1. Xu, C.; McDowell, N.G.; Fisher, R.A.; Wei, L.; Sevanto, S.; Christoffersen, B.O.; Weng, E.; Middleton, R.S. Increasing impacts of extreme droughts on vegetation productivity under climate change. Nat. Clim. Chang. 2019, 9, 948–953. [Google Scholar] [CrossRef] [Green Version]
  2. Vicente-Serrano, S.M.; Quiring, S.M.; Peña-Gallardo, M.; Yuan, S.; Domínguez-Castro, F. A review of environmental droughts: Increased risk under global warming? Earth Sci. Rev. 2019, 201, 102953. [Google Scholar] [CrossRef]
  3. Kraaijenbrink, P.D.A.; Bierkens, M.F.P.; Lutz, A.F.; Immerzeel, W.W. Impact of a global temperature rise of 1.5 degrees Celsius on Asia’s glaciers. Nature 2017, 549, 257–260. [Google Scholar] [CrossRef] [PubMed]
  4. Mengel, M.; Levermann, A.; Frieler, K.; Robinson, A.; Marzeion, B.; Winkelmann, R. Future sea level rise constrained by observations and long-term commitment. Proc. Natl. Acad. Sci. USA 2016, 113, 2597–2602. [Google Scholar] [CrossRef] [Green Version]
  5. Cai, W.; Wang, G.; Dewitte, B.; Wu, L.; Santoso, A.; Takahashi, K.; Yang, Y.; Carréric, A.; McPhaden, M.J. Increased variability of eastern Pacific El Niño under greenhouse warming. Nature 2018, 564, 201–206. [Google Scholar] [CrossRef]
  6. Wang, G.; Cai, W.; Gan, B.; Wu, L.; Santoso, A.; Lin, X.; Chen, Z.; McPhaden, M.J. Continued increase of extreme El Niño frequency long after 1.5 °C warming stabilization. Nat. Clim. Chang. 2017, 7, 568–572. [Google Scholar] [CrossRef]
  7. Zhu, Z.; Piao, S.; Yan, T.; Ciais, P.; Bastos, A.; Zhang, X.; Wang, Z. The accelerating land carbon sink of the 2000s may not be driven predominantly by the warming hiatus. Geophys. Res. Lett. 2018, 45, 1402–1409. [Google Scholar] [CrossRef]
  8. Liu, Y.; Piao, S.; Gasser, T.; Ciais, P.; Yang, H.; Wang, H.; Keenan, T.F.; Huang, M.; Wan, S.; Song, J.; et al. Field-experiment constraints on the enhancement of the terrestrial carbon sink by CO2 fertilization. Nat. Geosci. 2019, 12, 809–814. [Google Scholar] [CrossRef] [Green Version]
  9. IPCC. Climate Change 2013—The Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  10. Alward, R.D.; Detling, J.K.; Milchunas, D.G. Grassland vegetation changes and nocturnal global warming. Science 1999, 283, 229. [Google Scholar] [CrossRef] [Green Version]
  11. Cramer, W.; Kicklighter, D.W.; Bondeau, A.; Iii, B.M.; Churkina, G.; Nemry, B.; Ruimy, A.; Schloss, A.L. Comparing global models of terrestrial net primary productivity (NPP): Overview and key results. Glob. Chang. Biol. 1999, 5, 1–15. [Google Scholar] [CrossRef]
  12. Wang, Z.; Zhang, Y.; Yang, Y.; Zhou, W.; Gang, C.; Zhang, Y.; Li, J.; An, R.; Wang, K.; Odeh, I.; et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai–Tibet Plateau, in China. Ecol. Inf. 2016, 33, 32–44. [Google Scholar] [CrossRef]
  13. Cuo, L.; Zhang, Y.; Xu, R.; Zhou, B. Decadal change and inter-annual variability of net primary productivity on the Tibetan Plateau. Clim. Dyn. 2021, 56, 1837–1857. [Google Scholar] [CrossRef]
  14. Mao, F.; Du, H.; Li, X.; Ge, H.; Cui, L.; Zhou, G. Spatiotemporal dynamics of bamboo forest net primary productivity with climate variations in Southeast China. Ecol. Indic. 2020, 116, 106505. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Wang, Q.; Wang, Z.; Yang, Y.; Li, J. Impact of human activities and climate change on the grassland dynamics under different regime policies in the Mongolian Plateau. Sci. Total Environ. 2020, 698, 134304. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Yue, H.; Peng, Q.; He, C.; Hong, S.; Bryan, B.A. Recent responses of grassland net primary productivity to climatic and anthropogenic factors in Kyrgyzstan. Land Degrad. Dev. 2020, 31, 2490–2506. [Google Scholar] [CrossRef]
  17. Wang, Z.; Yang, Y.; Li, J.; Zhang, C.; Chen, Y.; Wang, K.; Odeh, I.; Qi, J. Simulation of terrestrial carbon equilibrium state by using a detachable carbon cycle scheme. Ecol. Indic. 2017, 75, 82–94. [Google Scholar] [CrossRef]
  18. Wang, Z. Estimating of terrestrial carbon storage and its internal carbon exchange under equilibrium state. Ecol. Model. 2019, 401, 94–110. [Google Scholar] [CrossRef]
  19. Wang, Z.; Chang, J.; Peng, S.; Piao, S.; Ciais, P.; Betts, R. Changes in productivity and carbon storage of grasslands in China under future global warming scenarios of 1.5 °C and 2 °C. J. Plant Ecol. 2019, 12, 804–814. [Google Scholar] [CrossRef]
  20. Raich, J.W.; Russell, A.E.; Vitousek, P.M. Primary productivity and ecosystem development along an elevational gradient on Mauna Loa, Hawaii. Ecology 1997, 78, 707–721. [Google Scholar] [CrossRef] [Green Version]
  21. Sierra Cornejo, N.; Leuschner, C.; Becker, J.N.; Hemp, A.; Schellenberger Costa, D.; Hertel, D. Climate implications on forest above- and belowground carbon allocation patterns along a tropical elevation gradient on Mt. Kilimanjaro (Tanzania). Oecologia 2021, 195, 797–812. [Google Scholar] [CrossRef]
  22. Luo, T.; Pan, Y.; Ouyang, H.; Shi, P.; Luo, J.; Yu, Z.; Lu, Q. Leaf area index and net primary productivity along subtropical to alpine gradients in the Tibetan Plateau. Glob. Ecol. Biogeogr. 2004, 13, 345–358. [Google Scholar] [CrossRef]
  23. Girardin, C.A.J.; Malhi, Y.; Aragão, L.E.O.C.; Mamani, M.; Huaraca Huasco, W.; Durand, L.; Feeley, K.J.; Rapp, J.; Silva-Espejo, J.E.; Silman, M.; et al. Net primary productivity allocation and cycling of carbon along a tropical forest elevational transect in the Peruvian Andes. Glob. Chang. Biol. 2010, 16, 3176–3192. [Google Scholar] [CrossRef] [Green Version]
  24. Lin, H.; Zhang, Y. Evaluation of six methods to predict grassland net primary productivity along an altitudinal gradient in the Alxa Rangeland, Western Inner Mongolia, China. Grassl. Sci. 2013, 59, 100–110. [Google Scholar] [CrossRef]
  25. Rangwala, I.; Miller, J.R. Climate change in mountains: A review of elevation-dependent warming and its possible causes. Clim. Chang. 2012, 114, 527–547. [Google Scholar] [CrossRef]
  26. Bertrand, R.; Lenoir, J.; Piedallu, C.; Riofrío-Dillon, G.; de Ruffray, P.; Vidal, C.; Pierrat, J.-C.; Gégout, J.-C. Changes in plant community composition lag behind climate warming in lowland forests. Nature 2011, 479, 517–520. [Google Scholar] [CrossRef] [PubMed]
  27. Qin, J.; Yang, K.; Liang, S.; Guo, X. The altitudinal dependence of recent rapid warming over the Tibetan Plateau. Clim. Change 2009, 97, 321–327. [Google Scholar] [CrossRef]
  28. Pepin, N.; Bradley, R.; Diaz, H.; Baraër, M.; Caceres, E.; Forsythe, N.; Fowler, H.; Greenwood, G.; Hashmi, M.; Liu, X. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 2015, 5, 424–430. [Google Scholar] [CrossRef] [Green Version]
  29. Gao, M.; Piao, S.; Chen, A.; Yang, H.; Liu, Q.; Fu, Y.H.; Janssens, I.A. Divergent changes in the elevational gradient of vegetation activities over the last 30 years. Nat. Commun. 2019, 10, 2970. [Google Scholar] [CrossRef]
  30. Pearson, R.G.; Phillips, S.J.; Loranty, M.M.; Beck, P.S.; Damoulas, T.; Knight, S.J.; Goetz, S.J. Shifts in Arctic vegetation and associated feedbacks under climate change. Nat. Clim. Chang. 2013, 3, 673–677. [Google Scholar] [CrossRef]
  31. Shen, M.; Zhang, G.; Cong, N.; Wang, S.; Kong, W.; Piao, S. Increasing altitudinal gradient of spring vegetation phenology during the last decade on the Qinghai–Tibetan Plateau. Agric. For. Meteorol. 2014, 189, 71–80. [Google Scholar] [CrossRef]
  32. Wang, Z.; Liu, X.; Wang, H.; Zheng, K.; Li, H.; Wang, G.; An, Z. Monitoring vegetation greenness in response to climate variation along the elevation gradient in the three-river source region of China. ISPRS Int. J. Geo. Inf. 2021, 10, 193. [Google Scholar] [CrossRef]
  33. Wang, Z.; Cui, G.; Liu, X.; Zheng, K.; Lu, Z.; Li, H.; Wang, G.; An, Z. Greening of the Qinghai–Tibet plateau and its response to climate variations along elevation gradients. Remote Sens. 2021, 13, 3712. [Google Scholar] [CrossRef]
  34. Running, S.W.; Zhao, M. MODIS/Terra Net Primary Production Gap-Filled Yearly L4 Global 500 m Sin Grid V061. 2021. Available online: https://lpdaac.usgs.gov/products/mod17a3hgfv006/ (accessed on 14 December 2021). [CrossRef]
  35. Becker, J.J.; Sandwell, D.T.; Smith, W.H.F.; Braud, J.; Binder, B.; Depner, J.; Fabre, D.; Factor, J.; Ingalls, S.; Kim, S.H.; et al. Global bathymetry and elevation data at 30 arc seconds resolution: SRTM30_PLUS. Mar. Geod. 2009, 32, 355–371. [Google Scholar] [CrossRef]
  36. Farr, T.G.; Rosen, P.A.; Caro, E.; Crippen, R.; Duren, R.; Hensley, S.; Kobrick, M.; Paller, M.; Rodriguez, E.; Roth, L.; et al. The Shuttle Radar Topography Mission. Rev. Geophys. 2007, 45, RG2004. [Google Scholar] [CrossRef] [Green Version]
  37. Miliaresis, G.C.; Argialas, D.P. Segmentation of physiographic features from the global digital elevation model/GTOPO30. Comput. Geosci. 1999, 25, 715–728. [Google Scholar] [CrossRef]
  38. Grohmann, C.H. Evaluation of TanDEM-X DEMs on selected Brazilian sites: Comparison with SRTM, ASTER GDEM and ALOS AW3D30. Remote Sens. Environ. 2018, 212, 121–133. [Google Scholar] [CrossRef] [Green Version]
  39. Bannari, A.; Kadhem, G.; El-Battay, A.; Hameid, N. Comparison of SRTM-V4. 1 and ASTER-V2. 1 for accurate topographic attributes and hydrologic indices extraction in flooded areas. J. Earth Sci. Eng. 2018, 8, 8–30. [Google Scholar] [CrossRef] [Green Version]
  40. Han, H.; Zeng, Q.; Jiao, J. Quality assessment of TanDEM-X DEMs, SRTM and ASTER GDEM on selected Chinese sites. Remote Sens. 2021, 13, 1304. [Google Scholar] [CrossRef]
  41. Muñoz-Sabater, J.; Dutra, E.; Agustí-Panareda, A.; Albergel, C.; Arduini, G.; Balsamo, G.; Boussetta, S.; Choulga, M.; Harrigan, S.; Hersbach, H.; et al. ERA5-Land: A state-of-the-art global reanalysis dataset for land applications. Earth Syst. Sci. Data 2021, 13, 4349–4383. [Google Scholar] [CrossRef]
  42. Gelaro, R.; McCarty, W.; Suárez, M.J.; Todling, R.; Molod, A.; Takacs, L.; Randles, C.A.; Darmenov, A.; Bosilovich, M.G.; Reichle, R.; et al. The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). J. Clim. 2017, 30, 5419–5454. [Google Scholar] [CrossRef]
  43. Kanamitsu, M.; Ebisuzaki, W.; Woollen, J.; Yang, S.-K.; Hnilo, J.; Fiorino, M.; Potter, G. NCEP-DOE AMIP-II reanalysis (R-2). Bull. Am. Meteorol. Soc. 2002, 83, 1631–1644. [Google Scholar] [CrossRef]
  44. Gang, C.; Wang, Z.; Chen, Y.; Yang, Y.; Li, J.; Cheng, J.; Qi, J.; Odeh, I. Drought-induced dynamics of carbon and water use efficiency of global grasslands from 2000 to 2011. Ecol. Indic. 2016, 67, 788–797. [Google Scholar] [CrossRef]
  45. Zhang, Y.; Zhang, C.; Wang, Z.; Chen, Y.; Gang, C.; An, R.; Li, J. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012. Sci. Total Environ. 2016, 563, 210–220. [Google Scholar] [CrossRef] [PubMed]
  46. Chen, J.; Yan, F.; Lu, Q. Spatiotemporal variation of vegetation on the Qinghai–Tibet plateau and the influence of climatic factors and human activities on vegetation trend (2000–2019). Remote Sens. 2020, 12, 3150. [Google Scholar] [CrossRef]
  47. Wang, K.; Dickinson, R.E. Global atmospheric downward longwave radiation at the surface from ground-based observations, satellite retrievals, and reanalyses. Rev. Geophys. 2013, 51, 150–185. [Google Scholar] [CrossRef]
  48. Yang, S.; Li, R.; Wu, T.; Hu, G.; Xiao, Y.; Du, Y.; Zhu, X.; Ni, J.; Ma, J.; Zhang, Y.; et al. Evaluation of reanalysis soil temperature and soil moisture products in permafrost regions on the Qinghai-Tibetan Plateau. Geoderma 2020, 377, 114583. [Google Scholar] [CrossRef]
  49. Shangguan, W.; Dai, Y.; Liu, B.; Zhu, A.; Duan, Q.; Wu, L.; Ji, D.; Ye, A.; Yuan, H.; Zhang, Q.; et al. A China data set of soil properties for land surface modeling. J. Adv. Modeling Earth Syst. 2013, 5, 212–224. [Google Scholar] [CrossRef]
  50. Du, Y.; Li, R.; Zhao, L.; Yang, C.; Wu, T.; Hu, G.; Xiao, Y.; Zhu, X.; Yang, S.; Ni, J.; et al. Evaluation of 11 soil thermal conductivity schemes for the permafrost region of the central Qinghai-Tibet Plateau. Catena 2020, 193, 104608. [Google Scholar] [CrossRef]
  51. Zucco, G.; Brocca, L.; Moramarco, T.; Morbidelli, R. Influence of land use on soil moisture spatial–temporal variability and monitoring. J. Hyd. 2014, 516, 193–199. [Google Scholar] [CrossRef]
  52. Yue, D.; Zhou, Y.; Guo, J.; Chao, Z.; Guo, X. Relationship between net primary productivity and soil water content in the Shule River Basin. Catena 2022, 208, 105770. [Google Scholar] [CrossRef]
  53. Guan, X.; Shen, H.; Li, X.; Gan, W.; Zhang, L. Climate control on net primary productivity in the complicated mountainous area: A case study of Yunnan, China. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 4637–4648. [Google Scholar] [CrossRef]
  54. Xu, H.-J.; Zhao, C.-Y.; Wang, X.-P. Elevational differences in the net primary productivity response to climate constraints in a dryland mountain ecosystem of northwestern China. Land Degrad. Dev. 2020, 31, 2087–2103. [Google Scholar] [CrossRef]
  55. Humphrey, V.; Berg, A.; Ciais, P.; Gentine, P.; Jung, M.; Reichstein, M.; Seneviratne, S.I.; Frankenberg, C. Soil moisture–atmosphere feedback dominates land carbon uptake variability. Nature 2021, 592, 65–69. [Google Scholar] [CrossRef] [PubMed]
  56. Green, J.K.; Seneviratne, S.I.; Berg, A.M.; Findell, K.L.; Hagemann, S.; Lawrence, D.M.; Gentine, P. Large influence of soil moisture on long-term terrestrial carbon uptake. Nature 2019, 565, 476–479. [Google Scholar] [CrossRef] [PubMed]
  57. Stocker, B.D.; Zscheischler, J.; Keenan, T.F.; Prentice, I.C.; Peñuelas, J.; Seneviratne, S.I. Quantifying soil moisture impacts on light use efficiency across biomes. N. Phytol. 2018, 218, 1430–1449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic diagram of the calculation of the NPP or environmental factors changes along EG. We take a 3 × 3 window as an example.
Figure 1. Schematic diagram of the calculation of the NPP or environmental factors changes along EG. We take a 3 × 3 window as an example.
Remotesensing 14 00713 g001
Figure 2. Trend of the NPP (a), MAT, MAR, and MASM (bd) anomaly from 2001 to 2020. The dashed straight lines denote the trendlines in subfigure (a). Asterisks (**) denote that the slope is statistically significant at the 0.01 level.
Figure 2. Trend of the NPP (a), MAT, MAR, and MASM (bd) anomaly from 2001 to 2020. The dashed straight lines denote the trendlines in subfigure (a). Asterisks (**) denote that the slope is statistically significant at the 0.01 level.
Remotesensing 14 00713 g002
Figure 3. Spatial pattern of global NPP (a) and its trend (b) from 2001 to 2020. The frequency of the uncertainty value is in the left of each subfigure. The regions with black dots in (b) indicate that the trend is statistically significant at the 0.05 level.
Figure 3. Spatial pattern of global NPP (a) and its trend (b) from 2001 to 2020. The frequency of the uncertainty value is in the left of each subfigure. The regions with black dots in (b) indicate that the trend is statistically significant at the 0.05 level.
Remotesensing 14 00713 g003
Figure 4. Spatial patterns of the trends of MAT (a, d, and g), MAR (b,e,h), MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2), and their spatial consistencies (j,k,l) from 2001 to 2020. The frequency of the uncertainty value is in the left of subfigures (j,k,l). The regions with black dots indicate that the trend is statistically significant at the 0.05 level. “Pos” and “Neg” denote the regions with the positive and negative agreement, respectively. “Non” denotes that the regions have not reached an agreement (the same meaning thereafter).
Figure 4. Spatial patterns of the trends of MAT (a, d, and g), MAR (b,e,h), MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2), and their spatial consistencies (j,k,l) from 2001 to 2020. The frequency of the uncertainty value is in the left of subfigures (j,k,l). The regions with black dots indicate that the trend is statistically significant at the 0.05 level. “Pos” and “Neg” denote the regions with the positive and negative agreement, respectively. “Non” denotes that the regions have not reached an agreement (the same meaning thereafter).
Remotesensing 14 00713 g004
Figure 5. Spatial patterns of the MAT (a,d,g), MAR (b,e,h), MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2), and their spatial consistencies (j,k,l) from 2001 to 2020.
Figure 5. Spatial patterns of the MAT (a,d,g), MAR (b,e,h), MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2), and their spatial consistencies (j,k,l) from 2001 to 2020.
Remotesensing 14 00713 g005
Figure 6. Spatial pattern of the changes in NPP (a) and its slope (b) along EG from 2001 to 2020. The value frequency is in the left of each subfigure. The regions with black dots indicate that the trend is statistically significant at the 0.05 level.
Figure 6. Spatial pattern of the changes in NPP (a) and its slope (b) along EG from 2001 to 2020. The value frequency is in the left of each subfigure. The regions with black dots indicate that the trend is statistically significant at the 0.05 level.
Remotesensing 14 00713 g006
Figure 7. Spatial patterns of MAT (a,d,g), MAR (b,e,h), MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2) along EG and their spatial consistencies (j,k,l) from 2001 to 2020. The frequency of the uncertainty value is in the left of subfigures (j,k,l).
Figure 7. Spatial patterns of MAT (a,d,g), MAR (b,e,h), MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2) along EG and their spatial consistencies (j,k,l) from 2001 to 2020. The frequency of the uncertainty value is in the left of subfigures (j,k,l).
Remotesensing 14 00713 g007
Figure 8. Spatial patterns of the slope of MAT (a,d,g), MAR (b,e,h), and MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2) along EG and their spatial consistencies (j,k,l) from 2001 to 2020. The frequency of the uncertainty value is in the left of subfigures (j,k,l).
Figure 8. Spatial patterns of the slope of MAT (a,d,g), MAR (b,e,h), and MASM (c,f,i) (derived from ERA5, MERRA2, and NCEP2) along EG and their spatial consistencies (j,k,l) from 2001 to 2020. The frequency of the uncertainty value is in the left of subfigures (j,k,l).
Remotesensing 14 00713 g008
Figure 9. Spatial patterns of the effect of MAT, MAR, and MASM derived from ERA5 (a), MERRA2 (b), and NCEP2 (c) on the changes in NPP and their spatial consistencies (d) from 2001 to 2020.
Figure 9. Spatial patterns of the effect of MAT, MAR, and MASM derived from ERA5 (a), MERRA2 (b), and NCEP2 (c) on the changes in NPP and their spatial consistencies (d) from 2001 to 2020.
Remotesensing 14 00713 g009
Figure 10. Regression curves, functions, and their R2 between NPP and environmental factors. (a) represent R2 between NPP and environmental factors. The regression curves and functions can be found in (bd). Subscript L of R2 means low elevations, and H means high elevations, and ** denotes the R2 is significant at the 0.01 level.
Figure 10. Regression curves, functions, and their R2 between NPP and environmental factors. (a) represent R2 between NPP and environmental factors. The regression curves and functions can be found in (bd). Subscript L of R2 means low elevations, and H means high elevations, and ** denotes the R2 is significant at the 0.01 level.
Remotesensing 14 00713 g010
Table 1. Summary of reanalysis data products used in this study.
Table 1. Summary of reanalysis data products used in this study.
NameCategoryTimespanSpatial ResolutionTemporal ResolutionReferencesData Acquisition
ERA5Temperature (K)1981–present0.1° × 0.1°Monthly[41]https://cds.climate.copernicus.eu (accessed on 10 December 2021)
Radiation (W/m2)
Volumetric soil water (m3 m−3)
MERRA2Temperature (K)1980–present0.5° × 0.625°Monthly[42]https://disc.gsfc.nasa.gov/ (accessed on 10 December 2021)
Radiation (W/m2)
Water root zone (m3 m−3)
NCEP2Temperature (K)1979–present1.875° × 1.904°Monthly[43]https://psl.noaa.gov/ (accessed on 10 December 2021)
Radiation (W/m2)
Volumetric Soil Moisture (m3 m−3)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Wang, H.; Wang, T.; Wang, L.; Huang, X.; Zheng, K.; Liu, X. Effects of Environmental Factors on the Changes in MODIS NPP along DEM in Global Terrestrial Ecosystems over the Last Two Decades. Remote Sens. 2022, 14, 713. https://doi.org/10.3390/rs14030713

AMA Style

Wang Z, Wang H, Wang T, Wang L, Huang X, Zheng K, Liu X. Effects of Environmental Factors on the Changes in MODIS NPP along DEM in Global Terrestrial Ecosystems over the Last Two Decades. Remote Sensing. 2022; 14(3):713. https://doi.org/10.3390/rs14030713

Chicago/Turabian Style

Wang, Zhaoqi, Hong Wang, Tongfang Wang, Lina Wang, Xiaotao Huang, Kai Zheng, and Xiang Liu. 2022. "Effects of Environmental Factors on the Changes in MODIS NPP along DEM in Global Terrestrial Ecosystems over the Last Two Decades" Remote Sensing 14, no. 3: 713. https://doi.org/10.3390/rs14030713

APA Style

Wang, Z., Wang, H., Wang, T., Wang, L., Huang, X., Zheng, K., & Liu, X. (2022). Effects of Environmental Factors on the Changes in MODIS NPP along DEM in Global Terrestrial Ecosystems over the Last Two Decades. Remote Sensing, 14(3), 713. https://doi.org/10.3390/rs14030713

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