Abstract
Northwest China has experienced dramatic changes in vegetation cover over the past few decades with the Yellow River Basin (YRB) being the most representative area. As the major climate-sensitive area in China, vegetation cover change is one reason for its impact on surface air temperature (SAT). This study uses the observation minus reanalysis (OMR) method to reveal the spatio-temporal variations of vegetation cover and its impact on SAT change over non-urban areas of the YRB from 1982 to 2015. The Global Inventory Modeling and Mapping Studies dataset, SAT derived from meteorological stations, and European Reanalysis (ERA)-interim reanalysis temperature data were used to analyze the relationship between normalized difference vegetation index (NDVI) and temperature variation caused by vegetation change. The NDVI trend (SlopeNDVI) of the entire YRB reached 1.11×10−2 decade−1, which indicated the recovery of vegetation in general. The impact of variation in vegetation conditions on SAT change during 1982–2015 was estimated to be 0.037 °C decade−1, which contributed 7.62% to the temperature change. The mean annual NDVI (MNDVI) and SlopeNDVI in the YRB were significantly negatively correlated (P<0.001) with OMR temperature variation. A negative correlation was exhibited in semi-arid and semi-humid regions, whereas a positive correlation was found in the arid region. The observed changes in vegetation and SAT in the YRB support the theory of the impact of vegetation variation on SAT in China.
HIGHLIGHTS
The relationship between vegetation and surface air temperature is different across the climate regions. The spatio-temporal changes of vegetation in non-urban areas of the Yellow River Basin were studied.
The contribution rate of the surface air temperature changes caused by vegetation changes during 1982–2015 is 7.62%. The vegetation changes warmed the surface temperature by 0.037 °C decade−1 during 1982–2015.
Graphical Abstract
INTRODUCTION
The global surface temperature in the first two decades of the 21st century (2001–2020) was 0.99 °C higher than in 1850–1900, on average (IPCC 2021). The impact of human activities since industrialization on the global climate system has been further affirmed in the The Intergovernmental Panel on Climate Change (IPCC) AR6. As one of these human activities, land use/land cover change (LULC) influences the biophysical properties of the surface; the surface controls the exchange of moisture, kinetic energy, and energy in the atmosphere, which is finally expressed as local or even regional surface temperature change (Bonan 2008). Therefore, when a large-scale change in vegetation cover occurs in a certain area, it directly impacts the surface energy balance, which may cause warming or cooling effects on the area. Over the past 50 years, the Chinese government has promoted the implementation of a series of environmental protection projects, such as the Grain for Green Project, the construction of the Three-North Shelterbelt Program, and the protection of natural forest resources, which have led to a significant change in vegetation cover in the YRB (Hu et al. 2020; Lu et al. 2021). The spatio-temporal distribution of vegetation and its response to the climate in the Yellow River Basin (YRB) is an essential part of promoting the water and carbon cycle in the ecosystem and has consistently been a topic in recent scientific research (Guo et al. 2021; Ye et al. 2021). Such research provides reference and guidance for the implementation of China's afforestation strategy. However, few studies have reported the impact of vegetation cover changes on climate systems.
Surface air temperature (SAT) is a parameter that can characterize the physical process of a surface, and it is the driving factor in surface–atmosphere energy exchange. It was used to quantitatively evaluate the influence of vegetation on climate in this study. The effects of variation in vegetation on the change of SAT varied substantially among regions, with complex factors, playing warming or cooling roles. If regional vegetation coverage increases, vegetation reduces the albedo of the area. The absorption of more radiation is accompanied by a heating effect on the surface. An increase in vegetation coverage means a larger leaf area index and more severe transpiration, which is the main cause of the cooling effect of the surface (Chi et al. 2021). The dynamic balance of the two effects leads to overall heating or cooling effects of vegetation on SAT. Related studies have concluded that the effect of vegetation on SAT is influenced by complex factors, such as vegetation cover types, seasonal variations, and geographical location; consequently, different conclusions may be drawn from different studies. As early as 2010, Yang et al. (2010) classified the study area with the normalized difference vegetation index (NDVI) as the index and found that the heating effect on the surface was obvious with a low NDVI, but with the increase in the NDVI, the change in SAT gradually decreased and even exhibited a negative value. Liu et al. (2013) pointed out that vegetation variation is negatively correlated with SAT in eastern China, and the impact of vegetation on SAT differs at different latitudes (Wickham et al. 2012; Wang et al. 2014; Shen et al. 2015). To evaluate the effects of variation of vegetation on change in SAT in the middle latitudes, this study considered the YRB as the study area, and the dataset selected was Global Inventory Modeling and Mapping Studies (GIMMS) NDVI 3rd generation (NDVI3 g) V1.0 (Soomro et al. 2021). The aim was to provide scientific evidence and support for the effect of vegetation variation on SAT in the middle latitudes, which may contribute to the development of vegetation construction policies.
One of the most commonly used methods for studying the influence of vegetation on climate is to build a climate model (Deborah Lawrence 2014). However, because of the coarse spatial resolution of a global climate model and the uncertainty of the physical processes, parameterization, and input data, a global climate model cannot reliably reproduce local climate effects. When data continuity and reliability are assured, the observation minus reanalysis (OMR) method can separate vegetation from complex factors and analyze local surface temperature changes. The OMR method can eliminate the changes in SAT caused by large-scale climate changes to allow the investigation of variation in vegetation (Kalnay & Cai 2003). The reanalysis data used in the OMR method can represent large-scale climate changes but are insensitive to regional surface processes because little or no surface data or information (surface observations of temperature, moisture, and wind overland) about land surface changes have been utilized in the data assimilation process (Wang & Yan 2015; Ao et al. 2020). Assuming that the LULC could affect the observed temperature, the difference between the reanalysis temperature and the observed temperature could indicate the influence of LULC on SAT change. After the basic idea of OMR was proposed in 2003, it was initially used to describe the impact of urbanization and land-use types on the SAT (Zhou et al. 2004; Lim et al. 2005). The application experiences have sped up the widespread application of OMR methods in research on the impact of LULC on climate in multiple countries and regions (Sun et al. 2011; Jin et al. 2020), which has verified the reliability of the OMR method. The current study uses the OMR method as the primary research method based on the observed SAT from the China Meteorological Data Service Center and ERA-interim reanalysis temperature.
The study aimed to analyze the temporal and spatial variation of vegetation in the YRB and the impact of vegetation change on SAT using the OMR method to provide theoretical support for understanding vegetation cover dynamics in the YRB and the bidirectional coupling between the atmosphere and vegetation.
MATERIALS AND METHODS
Study area
The geographical location of the YRB and its three climatic sub-regions.
Methods
Considering that the impact of urbanization on SAT in the process of the global warming cannot be ignored, we extracted unurbanized areas for this study (IPCC 2007). For different buffer radii of meteorological observation stations, the correlation between annual temperature trends and urban land expansion varied with high correlation values at 10–20 km (Ge et al. 2013). This rule applied to grid data after the temperature was interpolated. The YRB was divided into 20 km2×20 km2 using 1980 and 2015 land-use data; squares with changes in the construction land area fraction larger than 1% were defined as urban areas. The remaining areas were non-urban, allowing the sole consideration of the impacts of vegetation changes and background climate change on SAT.
This study used variations in NDVI, SAT, OMR, and ERA-interim reanalysis data as evaluation indexes when analyzing trends in the YRB and its three climate regions during 1982–2015 (He et al. 2017). This study analyzed the grid data of the YRB using the ordinary least-squares regression method to obtain the spatial distributions of SlopeNDVI, SlopeRT, and SlopeSAT. In addition, SlopeOMR was obtained by ‘SlopeSAT minus SlopeRT’. The contribution of vegetation changes was quantified by the ‘contribution rate of NDVI influence’, which is the proportion of vegetation change with respect to temperature change (Ao et al. 2020).
NDVI data
To quantify vegetation conditions, the NDVI data used in this study were the GIMMS NDVI3 g V1.0 dataset (Sun et al. 2021), which was obtained from the National Aeronautics and Space Administration (NASA) Ames Ecological Forecasting Lab (www.ecocast.arc.nasa.gov). The dataset was acquired by the advanced very high-resolution radiometer (AVHRR) of the National Oceanic and Atmospheric Administration (NOAA). It is one of the most extended NDVI time series globally, with a spatial resolution of 8 km and a temporal resolution of 15 days. In this study, data from 1982 to 2015 were selected for analysis (Ni et al. 2020; Zhai et al. 2020).
The dataset was processed in three steps: (1) the maximum value composite (MVC) method was used to merge the 15-day raster files into monthly NDVI data (Sun et al. 2021); (2) the mean annual NDVIs from 1982 to 2105 were obtained by averaging monthly NDVI values (Wu et al. 2020); (3) the actual value of the data ranged from −1 to 1. To eliminate the impact of other factors, the study focused on variation in vegetation. Grid cells with negative values were assigned the value of zero to eliminate the effects of water.
Surface air temperature data
The spatial interpolation dataset of the annual mean SAT in China was taken from the China Meteorological Data Service Center (http://data.cma.cn/en). The spatial interpolation dataset is based on the daily observation data of over 2,400 meteorological stations in China. It interpolated the data to a grid with a 1 km spatial resolution using the ANUSPLIN software (Zhou 2019). The annual mean SAT data from 1982 to 2015 in the YRB was applied for the study.
Reanalysis temperature data
The spatial distribution of correlation coefficients between observed temperatures and ERA-interim reanalysis during 1982–2015 in the YRB.
The spatial distribution of correlation coefficients between observed temperatures and ERA-interim reanalysis during 1982–2015 in the YRB.
The annual mean reanalysis of air temperature with a spatial resolution of 7.5 arc min was calculated from the monthly means of daily means 2 m temperature data. Because of the different spatial resolutions of the datasets, the bilinear interpolation method was used to resample and unify the data resolution for better subsequent data analysis.
Land-use data
To extract the areas affected by urbanization in the YRB from 1982 to 2015, we selected land-use data with resolving 1 km provided by the Data Center for Resources and Environmental Sciences (https://www.resdc.cn). Land-use data for 1980 and 2015 were generated by manual visual interpretation using the Landsat thematic mapper (Landsat TM)/enhanced thematic mapper (ETM) and Landsat8 remote sensing images (Deng et al. 2015). These data serve as the basis for determining the urbanization areas and non-urbanization areas.
RESULTS
Excluding the SAT impact of urbanization
Percentage of land-use types in 1980 and 2015
. | Farmland . | Forest . | Grassland . | Water . | Constructed land . | Unused land . |
---|---|---|---|---|---|---|
1980 (%) | 26.2 | 12.8 | 48.1 | 1.8 | 2 | 9.1 |
2015 (%) | 25.8 | 13.1 | 47.5 | 1.7 | 2.9 | 9 |
. | Farmland . | Forest . | Grassland . | Water . | Constructed land . | Unused land . |
---|---|---|---|---|---|---|
1980 (%) | 26.2 | 12.8 | 48.1 | 1.8 | 2 | 9.1 |
2015 (%) | 25.8 | 13.1 | 47.5 | 1.7 | 2.9 | 9 |
The spatial distribution of land-use types in: (a) 1980 and (b) 2015 in the YRB. The regions in the white squares are areas that experienced urbanization between 1980 and 2015.
The spatial distribution of land-use types in: (a) 1980 and (b) 2015 in the YRB. The regions in the white squares are areas that experienced urbanization between 1980 and 2015.
Vegetation dynamics
The spatial distribution of: (a) NDVI trend (SlopeNDVI) and (b) OMR temperature trend (SlopeOMR) in the YRB during 1982–2015.
The spatial distribution of: (a) NDVI trend (SlopeNDVI) and (b) OMR temperature trend (SlopeOMR) in the YRB during 1982–2015.
Figure 6(b) shows the trend of OMR temperature change in the YRB, with SlopeOMR values between −0.3 and 0.5 °C decade−1. The OMR temperature in the southeastern triangular area declined, and the north of the arid region also showed a clear downward trend. The OMR temperature west of the YRB exhibited an upward trend.
Temporal changes in the NDVI and temperature
The inter-annual variation of the annual mean NDVI from 1982 to 2015 on: (a) the Yellow River Basin (YRB) and its (b) arid; (c) semi-arid; and (d) semi-humid regions. The inter-annual variation of annual mean SAT and reanalysis temperatures during 1982–2015 in (e) the YRB and its (f) arid, (g) semi-arid; and (h) semi-humid regions. The blue and red lines refer to the linear trends of SAT and reanalysis temperature, respectively. (i)–(l) refer to the inter-annual variation in annual mean OMR temperature from 1982 to 2015 in the YRB, arid, semi-arid, and semi-humid regions. Please refer to the online version of this paper to see this figure in colour: http://dx.doi.org/10.2166/wcc.2022.037.
The inter-annual variation of the annual mean NDVI from 1982 to 2015 on: (a) the Yellow River Basin (YRB) and its (b) arid; (c) semi-arid; and (d) semi-humid regions. The inter-annual variation of annual mean SAT and reanalysis temperatures during 1982–2015 in (e) the YRB and its (f) arid, (g) semi-arid; and (h) semi-humid regions. The blue and red lines refer to the linear trends of SAT and reanalysis temperature, respectively. (i)–(l) refer to the inter-annual variation in annual mean OMR temperature from 1982 to 2015 in the YRB, arid, semi-arid, and semi-humid regions. Please refer to the online version of this paper to see this figure in colour: http://dx.doi.org/10.2166/wcc.2022.037.
Similarly, Figure 7(e)–7(h) describes the changing trends of observed and reanalyzed temperatures. The changing trends in the two types of datasets were consistent, and the annual mean reanalysis temperature was higher than the annual mean surface observation temperature. The temperature of the YRB had clear upward trends of between 0.40 and 0.50 °C decade−1 in different climate regions.
Figure 7(i)–7(l) shows the change in the OMR temperature, which is caused by changes in vegetation. The trend of increasing OMR temperature in the YRB was 0.04 °C decade−1 (R2 = 0.11716; Figure 7(i)). The OMR trends (SlopeOMR) in the three climatic sub-regions were estimated to be 0.023, 0.024, and 0.047 °C decade−1 in the arid, semi-humid, and semi-arid regions, respectively. The R2 for the arid, semi-arid, and semi-humid regions were 0.02427, 0.15167, and 0.04142, respectively.
The impact of variations in the NDVI on SAT was consistent across the YRB (Table 2). The YRB SlopeOMR was 0.037 °C decade−1, contributing 7.62% to the SAT. Among the three climate regions, the contribution of NDVI influence was highest in semi-arid areas, accounting for 9.36%, followed by arid and semi-humid areas, with contributions of 5.35 and 5.05%, respectively. This means that the contribution of vegetation change to SAT in the YRB was no greater than 10%.
NDVI trend (SlopeNDVI), OMR trend (SlopeOMR), SAT trend (SlopeSAT), and contribution rate of NDVI influence of non-urban region during 1982–2015
1982–2015 . | The Yellow River Basin . | Arid region . | Semi-arid region . | Semi-humid region . |
---|---|---|---|---|
SlopeNDVI (10−2 decade−1) | 1.11 | 0.59 | 1.15 | 1.32 |
SlopeOMR (°C decade−1) | 0.037 | 0.023 | 0.047 | 0.024 |
SlopeSAT (°C decade−1) | 0.48 | 0.43 | 0.50 | 0.48 |
Contribution rate of NDVI influence (%) | 7.62 | 5.35 | 9.36 | 5.05 |
1982–2015 . | The Yellow River Basin . | Arid region . | Semi-arid region . | Semi-humid region . |
---|---|---|---|---|
SlopeNDVI (10−2 decade−1) | 1.11 | 0.59 | 1.15 | 1.32 |
SlopeOMR (°C decade−1) | 0.037 | 0.023 | 0.047 | 0.024 |
SlopeSAT (°C decade−1) | 0.48 | 0.43 | 0.50 | 0.48 |
Contribution rate of NDVI influence (%) | 7.62 | 5.35 | 9.36 | 5.05 |
Note: Contribution rate of NDVI influence = SlopeOMR/SlopeSAT*100%.
Relationship between the NDVI and OMR temperature
(a)–(d) Scatter diagrams between the SlopeOMR and mean annual NDVI (MNDVI) during 1982–2015 in (a) the YRB; (b) arid region; (c) semi-arid region; and (d) semi-humid region. The blue lines represent the mean annual NDVI. (e) and (f) Scatter diagrams between the SlopeOMR and the SlopeNDVI during 1982–2015 in (e) the YRB; (f) arid region; (g) semi-arid region; and (h) semi-humid region. The figures in blue and red colors represent the percent of the grids, with SlopeOMR being larger and less than 0 decade−1, respectively. Please refer to the online version of this paper to see this figure in colour: http://dx.doi.org/10.2166/wcc.2022.037.
(a)–(d) Scatter diagrams between the SlopeOMR and mean annual NDVI (MNDVI) during 1982–2015 in (a) the YRB; (b) arid region; (c) semi-arid region; and (d) semi-humid region. The blue lines represent the mean annual NDVI. (e) and (f) Scatter diagrams between the SlopeOMR and the SlopeNDVI during 1982–2015 in (e) the YRB; (f) arid region; (g) semi-arid region; and (h) semi-humid region. The figures in blue and red colors represent the percent of the grids, with SlopeOMR being larger and less than 0 decade−1, respectively. Please refer to the online version of this paper to see this figure in colour: http://dx.doi.org/10.2166/wcc.2022.037.
In the non-urban areas of the YRB from 1982 to 2015, SlopeOMR decreased (Figure 8(e)), with a SlopeNDVI negative correlation coefficient of −0.582 (P < 0.001). In addition, the percentage of grids with SlopeOMR less than and larger than 0 decade−1 was approximately 38 and 62%, respectively. Figure 8(g) and 8(h) shows that a negative correlation still occurred between SlopeNDVI and SlopeOMR in the semi-arid and semi-humid regions, with correlation coefficients of −0.702 (P < 0.001) and −0.589 (P < 0.001), respectively. However, SlopeNDVI and SlopeOMR in arid regions were positively correlated, with a correlation coefficient of 0.142 (P < 0.01). The ratio of a grid with SlopeOMR was less or larger than 0 decade−1. The greatest imbalances in the semi-arid regions were 68 and 32%, respectively. Semi-humid regions were more uniform with 54 and 46% ratios, and the ratio in arid regions was in the middle with 60 and 40%, respectively.
DISCUSSION
Vegetation cover and SAT change over the YRB
Based on NDVI data from 1982 to 2015, this study analyzed the spatial and temporal variation of vegetation in non-urban areas of the YRB. The results show that the vegetation change trend in the YRB was 1.11 × 10−2 decade−1 from 1982 to 2015, which is similar to the vegetation change trend of 1.58 × 10−2 decade−1 in the same period in the Loess Plateau (Jin et al. 2020). The difference between the results was caused by large spatial differences in vegetation cover in the YRB. The source of the Yellow River is on the Qinghai-Tibet Plateau, where the vegetation trend is much lower than that in the central and eastern regions of the YRB (Guan et al. 2015; Ma et al. 2018). Compared with the global NDVI change trend, which was 0.46 × 10−2 decade−1 in 1982–2012, the vegetation change trend in the YRB was much higher. It demonstrates the effectiveness of a series of ecological environmental protection policies implemented by the Chinese government (Liu et al. 2015). Recapping vegetation in the arid region of the YRB was poor and showed the lowest changing trend (SlopeNDVI = 0.59 × 10−2 decade−1). However, the trend was still higher than the result for NDVI during the 1982–2013 growth season in Inner Mongolia, which was SlopeNDVI = 0.3 × 10−2 decade−1 (Tong et al. 2017). The differences between these studies and Tong's results are related to the climate sub-region division pattern, whether urbanized areas are excluded and how the annual mean NDVI is calculated.
The trend of temperature change in the YRB affected by vegetation changes during 1982–2015 was 0.037 °C decade−1, which showed that vegetation restoration in the YRB had a heating effect on SAT. The heating impact of vegetation change on the surface is mainly due to the radiation effect. An increase in vegetation leads to a decrease in albedo, which manifests as a heating effect. Vegetation restoration in the YRB is mainly manifested in the expansion of the forest, and the albedo of the forest is lower than those of grassland and farmland, which may be a cause of local warming.
In the YRB, although the NDVI changes varied in different climate regions from 1982 to 2015, the contributions of NDVI influences were similar and concentrated between 5 and 10%. This means that among the factors affecting the SAT, the effects of vegetation cover changes were limited, and their impact was not dominant among the many complicated factors. The change in SAT may be mainly derived from the change in the climate background and the increase in carbon dioxide emissions caused by human activities (IPCC 2021).
Analysis of the relationship between vegetation change and SAT change
The difference between the observed and reanalyzed temperatures was used to determine the temperature affected by vegetation change using the OMR method. According to the results, the NDVI (MNDVI) in the YRB was negatively correlated with the SlopeOMR (Ning et al. 2017). However, the correlation between MNDVI and SlopeOMR was not always negative in different climate regions. The negative correlation (r = −0.067) was already very weak when the mean MNDVI was 0.29 in a semi-arid region. When the average value of MNDVI reached only 0.16 in the arid region, MNDVI and SlopeOMR showed a weak positive correlation (r = 0.154). This study showed that when the vegetation coverage reached a certain degree, the better the vegetation coverage was, and the more restrained the warming effect on SAT. When the vegetation coverage level was low, the influence of vegetation on SAT was less consistent and showed a certain regionality.
In Liu's research, the OMR temperature showed a cooling trend in eastern China when the annual mean NDVI was greater than 0.55 (Liu et al. 2013). However, the value of MNDVI in semi-humid areas with the best coverage in the YRB reached only 0.45. Therefore, the cooling effect of vegetation restoration on the SAT has not yet been confirmed. Although the results still showed a warming effect (SlopeOMR = 0.024 °C decade−1) on the surface, the negative correlation between SlopeNDVI and SlopeOMR was consistent with the findings of Liu.
With respect to exploring the correlation between SlopeNDVI and SlopeOMR, we found a significant negative correlation in semi-arid and semi-humid regions with high latent heat exchange levels and a positive correlation in arid regions with significant sensible heat exchange. In addition, rapid vegetation restoration in areas with suitable vegetation cover may have a cooling effect on SAT.
The heating and cooling effects of vegetation change on SAT are dynamic. Areas with sparse vegetation have higher levels of sensible heat exchange. When the heat cannot be absorbed by sufficient water evaporation, it accumulates on the earth's surface, and the SAT is higher. With the improvement in vegetation cover, sensible heat exchange is gradually replaced by latent heat exchange, and the transpiration of vegetation increases to achieve the effect of cooling the SAT (Gou et al. 2018). However, the improvement in vegetation cover will lead to a decrease in surface reflectivity because the albedo of bare land is higher than that of vegetation. Areas with increased vegetation cover absorbed more solar radiation, which had a warming effect on the surface. The two effects interact in the dynamic response of SAT to vegetation changes. The YRB is a geographically diverse river basin with specific spatial differences in vegetation types, and cover levels are one reason for the observed differences in relationships among the different study regions.
The analyses in this study were based on different climatic conditions. However, the YRB has a wide range of complex terrains. The relationship between vegetation and temperature cannot be completely described based on the division of climate characteristics. In subsequent studies, taking vegetation type or land-use type as the basis for division will be a direction for future work.
CONCLUSION
Based on the OMR method, we analyzed the spatio-temporal variation of vegetation cover and its impact on SAT change in non-urban areas of the YRB from 1982 to 2015. This provides a basis for understanding the characteristics of vegetation cover in the YRB. It also provides a reference for deepening the understanding of the coupling relationship between vegetation and climate.
The vegetation cover in the YRB showed clear spatial differences. Areas with suitable vegetation cover were mainly concentrated in the southeastern triangle of the basin. The source of the Yellow River and the northwestern region of the basin are poor. In the time series of this study, the change trends of the NDVI increased from north to south. The change trends in the arid, semi-arid, and semi-humid regions were estimated to be 0.59, 1.15, and 1.32 × 10−1 decade−1, respectively. Overall, the recovery rate of the YRB was 1.11 × 10−1 decade−1. Thus, the vegetation recovery results were more significant in wetter regions.
During the 34 years from 1982 to 2015, SAT in the YRB showed a warming trend of 0.48 °C decade−1. The temperature change trend affected by vegetation was 0.037 °C decade−1, with a contribution of 7.62% by the NDVI. Vegetation changes had different warming effects on the surface temperature over the three climatic regions, with the results of SlopeOMR in arid, semi-arid, and semi-humid areas being 0.023, 0.047, and 0.024 °C decade−1, respectively. These results show that vegetation changes in the YRB from 1982 to 2015 had warming effects on the SAT, but the contribution of NDVI was not high (less than 10%). Therefore, the restoration of vegetation guided by a series of vegetation restoration policies in the YRB has neither had a significant impact on the local climate nor is it the main reason for the warming of the YRB.
After exploring the correlation between NDVI and SAT, the following conclusions were drawn: the strong negative correlations among MNDVI, SlopeNDVI, and SlopeOMR can be used as a reference only in areas with good vegetation coverage. In areas with general vegetation coverage and even poor conditions, such as the arid region of the YRB divided in this study, this negative correlation is not obvious with only a weak correlation; this correlation law has strong regional characteristics. Therefore, local vegetation coverage and climate conditions should be fully considered when formulating policies such as regional vegetation protection.
AUTHORS CONTRIBUTIONS
For this research paper with several authors, a short paragraph specifying their individual contributions was provided. Z.H. developed the original idea and contributed to the research design for the study. Q.W. was responsible for data collection. C.H. and R.L. provided guidance and contributed to the research design. D.H. and S.-S. provided guidance for the writing of the article. All authors have read and approved the final manuscript.
FUNDING
This work was funded by the National Natural Science Foundation of China (51979250), the National key research priorities program of China (2019YFC1510703), and Key projects of the National Natural Science Foundation of China (51739009).
DATA AVAILABILITY STATEMENT
Data cannot be made publicly available; readers should contact the corresponding author for details.
CONFLICT OF INTEREST
The authors declare there is no conflict.