Next Article in Journal
A New Framework for Generating Indoor 3D Digital Models from Point Clouds
Next Article in Special Issue
Desertification Mitigation in Northern China Was Promoted by Climate Drivers after 2000
Previous Article in Journal
Changes in Vegetation Cover and the Relationship with Surface Temperature in the Cananéia–Iguape Coastal System, São Paulo, Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Evaluation of Ecosystem Quality and Its Response to Aridity on the Qinghai–Tibet Plateau

1
School of Ecology and Nature Conservation, Beijing Forestry University, Beijing 100083, China
2
School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798, Singapore
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(18), 3461; https://doi.org/10.3390/rs16183461
Submission received: 7 July 2024 / Revised: 24 August 2024 / Accepted: 12 September 2024 / Published: 18 September 2024

Abstract

:
Exploring the response of spatial and temporal characteristics of ecological quality change to aridity on the Qinghai–Tibet Plateau (QTP) can provide valuable information for regional ecological protection, water resource management, and climate change adaptation. In this study, we constructed the Remote Sensing Ecological Index (RSEI) and Standardized Precipitation Evapotranspiration Index (SPEI) based on the Google Earth Engine (GEE) platform with regional characteristics and completely analyzed the spatial and temporal variations of aridity and ecological quality on the QTP in the years 2000, 2005, 2010, 2015, and 2020. Additionally, we explored the responses of ecological quality to aridity indices at six different time scales. The Mann–Kendall test, correlation analysis, and significance test were used to study the spatial and temporal distribution characteristics of meteorological aridity at different time scales on the QTP and their impacts on the quality of the ecological environment. The results show that the ecological environmental quality of the QTP has a clear spatial distribution pattern. The ecological environment quality is significantly better in the south-east, while the Qaidam Basin and the west have lower ecological environment quality indices, but the overall trend of environmental quality is getting better. The Aridity Index of the QTP shows a differentiated spatial and temporal distribution pattern, with higher Aridity Indexes in the north-eastern and south-western parts of the plateau and lower Aridity Indexes in the central part of the plateau at shorter time scales. Monthly, seasonal, and annual-scale SPEI values showed an increasing trend. There is a correlation between aridity conditions and ecological quality on the QTP. The areas with significant positive correlation between the RSEI and SPEI in the study area were mainly concentrated in the south-eastern, south-western, and northern parts of the QTP, where the ecological quality of the environment is more seriously affected by meteorological aridity.

1. Introduction

According to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC), as global warming continues, more regions are expected to experience an increased frequency and severity of aridity. Aridity, as one of the most devastating and widespread extreme climatic conditions, has severe impacts on water resources, ecosystems, agricultural production, and economic development [1]. The main types of aridity include meteorological, agricultural, hydrological, and socio-economic aridity [2,3]. The degree of aridity significantly influences vegetation cover types and structures, thereby affecting the quality of local ecosystems [4]. With ongoing global warming, both the intensity and frequency of aridity are expected to increase, further impacting vegetation productivity. Prolonged and intense aridity can even lead to substantial changes in ecosystem quality [5]. Therefore, understanding the response of ecological environment quality to meteorological aridity is crucial for regional ecological protection and sustainable development.
The ecological environment forms the foundation of human survival, economic development, and social progress. However, with advances in science and technology and the intensification of human activities, ecosystems globally have suffered considerable damage. This has led to an increasing frequency of ecological problems and poses serious threats to regional ecological security. Conducting ecological quality analyses and evaluations can provide targeted recommendations for environmental protection and is vital for promoting harmony between humanity and nature, as well as sustainable development. The Qinghai–Tibet Plateau (QTP) is characterized by a complex landscape, harsh natural environment, and extremely fragile ecosystem. In recent years, significant climate change and intensified human activities have resulted in varying degrees of degradation of the plateau’s ecosystem. Consequently, the ecological quality of the QTP has garnered considerable attention from both domestic and international scholars, leading to numerous studies on vegetation changes, grassland degradation, and related issues. For example, recent research has analyzed the spatio-temporal patterns of Net Primary Productivity (NPP) changes in the QTP’s vegetation over the past two decades, clarifying its relationship with various influencing factors [6]. Other studies have examined the spatio-temporal variations in water use efficiency and the effects of climatic factors on the QTP [7], as well as the ecological vulnerability of the QTP, identifying vegetation as a major driver [8]. Additionally, research has explored the combined effects of climate change and human activities on the spatio-temporal dynamics of vegetation in the QTP by conducting correlation analyses using data from meteorological stations and regional economic statistics [9]. Studies have demonstrated the effectiveness of the Remote Sensing Ecological Index (RSEI) in monitoring and assessing the ecological conditions of the QTP. For instance, study has highlighted significant spatio-temporal variations in ecological quality across the region over the past two decades, with periods of ecological degradation primarily driven by rising temperatures and overgrazing. Notably, between 2015 and 2020, overgrazing in the south-western regions of the QTP was identified as a major factor contributing to poor ecological conditions, as indicated by lower RSEI values [10]. Additionally, studies focused on ecological assessments within the QTP, such as the research conducted in the Changtang Nature Reserve, have used the RSEI to evaluate changes in ecological quality over extended periods. These studies have revealed both periods of ecological deterioration and subsequent recovery, underscoring the value of the RSEI in tracking long-term ecological trends [11]. Furthermore, study on ecological vulnerability across the QTP has utilized the RSEI to map and analyze patterns of vulnerability, demonstrating its utility in guiding environmental protection and policy-making [12]. These studies underscore the RSEI’s effectiveness as a reliable tool for assessing environmental quality and informing ecological management strategies in the Qinghai–Tibet Plateau. However, focusing on changes in specific ecosystem features alone makes it difficult to synthesize a comprehensive assessment of the plateau’s overall ecological state.
Aridity is a multi-scale phenomenon. The variability in water availability, including soil moisture, groundwater, reservoir storage, river discharge, and snowpack, exemplifies the complexity of aridity [13]. The time lag between water input and the availability of usable resources varies significantly. Therefore, understanding the accumulation time scales of water scarcity, which functionally distinguishes between hydrological, environmental, agricultural, and other forms of aridity, is critical. For instance, the hydrological system’s response to precipitation can differ widely over time [14], determined by the varying frequencies of hydrological and climatic variables [15]. Consequently, the Aridity Index must be associated with a specific time scale to facilitate the monitoring and management of different water resources. The Aridity Index is an effective method for assessing aridity, simplifying the complex phenomenon and measuring the severity of aridity events [16]. Numerous aridity indices have been proposed internationally, with commonly used indices for meteorological aridity monitoring including the Relative Humidity Index (M), Standardized Precipitation Index (SPI), and Standardized Precipitation Evapotranspiration Index (SPEI) [17]. Among these, the SPEI is widely used due to its integration of multiple climate factors and its applicability across different timescales. The Standardized Precipitation Evapotranspiration Index (SPEI) evolves from the Normalized Precipitation Index (NPI) [18], incorporating factors such as precipitation, temperature, and potential evapotranspiration to provide a comprehensive measure of water balance and aridity conditions [19].
In the study, the SPEI at several time scales of 1, 3, 6, 9, 12, and 24 are selected to reflect different characteristics of meteorological drought. By analyzing different time scales, the severity and impact range of drought can be assessed more comprehensively. The one-month SPEI (SPEI-1) serves as a critical indicator for monitoring short-term meteorological anomalies, such as seasonal fluctuations or abrupt arid events. As the time scale expands, the three-month SPEI (SPEI-3) delineates the drought conditions over a complete season, providing a more integrated view of seasonal aridity patterns. The six-month SPEI (SPEI-6) extends this assessment to a semi-annual basis, offering insights into medium-term drought dynamics. The annual SPEI (SPEI-12) synthesizes the year-long precipitation and evaporation trends, reflecting the annual drought status and its implications for agricultural planning and water resource management. Furthermore, the two-year SPEI (SPEI-24) encapsulates long-term drought conditions, revealing persistent drought patterns and the potential long-term impacts of climate change on hydrological cycles and ecosystem health.
In recent years, numerous approaches have been developed for monitoring ecosystem status and detecting changes. Remote sensing techniques, celebrated for their effectiveness and precision, have become a staple in ecological assessments. These include indices such as the Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Normalized Difference Water Index (NDWI), Normalized Difference Moisture Index (NDMI), and Land Surface Temperature (LST). The Remote Sensing-based Ecological Index (RSEI) [20] was introduced as a comprehensive ecological assessment tool, integrating these remote sensing indices with socio-economic dimensions. The RSEI combines four key remote sensing indicators: NDVI for vegetation greenness, Wetness Index (WET) for moisture, NDBSI for aridity, and LST for thermal conditions. The specific contributions of each factor to ecological quality are determined by covariance-based Principal Component Analysis (PCA). The synthesis of these indicators into the RSEI offers an objective and meaningful approach to assessing ecological quality across different regions and ecosystems [21,22,23]. Therefore, the RSEI is a feasible and reliable tool for assessing the changing environmental quality in the QTP region. Additionally, the acquisition of long time-series and high-resolution remote sensing datasets, coupled with the complicated data preprocessing and RSEI calculation, is often time-consuming and costly. However, the cloud platform Google Earth Engine (GEE) offers high-performance resources for processing large geospatial datasets, providing a significant advantage for large-scale and long time-series remote sensing research.
In this study, the SPEI was selected as the indicator for monitoring meteorological aridity, while the RSEI data from the GEE platform was used to monitor ecological environment quality. The analysis focused on the spatial and temporal variations of aridity and ecological environment quality in the QTP from 2000 to 2020. By combining multi-scale meteorological aridity data, this study aims to elucidate the impact mechanism of long-term aridity on ecological environment quality in the QTP. This study will provide a theoretical foundation and scientific basis for the restoration and construction of vegetation cover, as well as the management and utilization of water resources in the QTP. Meteorological aridity affects regional ecosystems in various ways, significantly impacting vegetation growth and recovery, land degradation, and forest fire incidence.

2. Materials and Methods

2.1. Study Area

The QTP is located in the southern part of the Asian continent and covers a total area of about 2.62 million km2, encompassing Tibet, Qinghai, Xinjiang, Gansu, Sichuan, and Yunnan Provinces. The average elevation of the entire plateau is greater than 4000 m above sea level, showing a trend of terrain tilting from north-west to south-east. The complexity and diversity of the terrain form the basic geological structure of the QTP (Figure 1). The natural environment of the QTP is very complex, with a rich and diverse vegetation. The main vegetation types include alpine grasslands and meadows, which form the unique ecological landscape of the plateau [24].
As a result of the high average altitude, the temperature on the plateau is relatively low, with an overall distribution characteristic of “low in the middle and high on both sides” that decreases from north-east to south-west. At the same time, precipitation shows large spatial differences, gradually decreasing from over 2000 mm in the south-eastern Shannan region to even less than 25 mm in the north-western Qaidam Basin. The QTP region is densely populated with lakes and crisscrossed with glaciers, and natural landscapes such as wetlands and swamps are widely distributed. These lakes, glaciers, and wetlands are important carriers of the country’s water resources storage. The water systems of the QTP play a vital role in supplying water resources in the region, maintaining ecological balance and supporting the livelihoods of the local population. The QTP is not only a geographically important region in China but also plays an important role in the global ecosystem [25]. Therefore, an in-depth study of the geographic, climatic, and ecological characteristics of the QTP, with particular attention to the effects of meteorological aridity on vegetation and ecosystems, is essential to the scientific understanding and effective conservation of the ecosystems of this region [26].

2.2. Data Sources

The dataset for evapotranspiration in the Chinese region is derived from the Hargreaves formula, which calculates potential evapotranspiration using monthly average, minimum, and maximum temperatures [27]. Spanning the years 1901 to 2022, this dataset boasts a spatial resolution of about 1 km by 1 km.
In contrast, the precipitation dataset for China is crafted through the Delta spatial downscaling method, utilizing the 0.5-degree global climate data from the Climate Research Unit (CRU) and the high-resolution climate data supplied by WorldClim. This dataset mirrors the temporal extent of 1901 to 2022 and maintains a similar spatial resolution of roughly 1 km by 1 km [28]. The temperature dataset was generated using the Delta spatial downscaling program. This dataset is based on the 0.5 global climate dataset published by CRU (Climatic Research Unit) and the global high-resolution climate dataset published by WorldClim. In addition, this dataset has been validated by 496 independent meteorological observation sites, and the validation results are credible [29]. The total length of the dataset is 1901~2021, and the spatial resolution is about 1 km × 1 km.
The RSEI assessment system encompasses four critical indicators: vegetation cover (greenness), soil moisture (humidity), surface temperature (heat), and building cover (desiccation). Xu et al. [20] reported that vegetation cover is measured by the Normalized Difference Vegetation Index (NDVI) [30], and soil moisture is indicated by the Wet Component of the Coma Cap Transformation (WET) [31,32]. Building cover is evaluated using the Normalized Differential Impervious Surface Index (NDBSI) [33], incorporating the Index-based Built-up area Index (IBI) and Soil Index (SI) [33]. Surface temperature is indicated by the Land Surface Temperature (LST). These indices are extensively used to evaluate the ecological status across various regions [34,35,36].
The construction of the RSEI indices was facilitated by utilizing a suite of MODIS (Moderate Resolution Imaging Spectroradiometer) products, encompassing MOD13A1, MOD11A2, and MOD09A1 datasets. These datasets were sourced from the United States Geological Survey (USGS) and have been seamlessly integrated into the Google Earth Engine (GEE) cloud computing platform, covering a period from 1 January 2001 to 31 December 2020. The MODIS satellite imagery selected for this study was specifically captured during the peak growing season, which spans from July to September. The spatial and temporal resolutions of these images are depicted in Figure 1.

2.3. Methods

2.3.1. Standardized Precipitation Evapotranspiration Index

The Standardized Precipitation Evapotranspiration Index (SPEI) is an Aridity Index based on the Standardized Precipitation Index (SPI). The SPEI can be calculated on different time scales (n = 1, 3, 6, 12), where the 1-month time scale (SPEI-1) is calculated using the precipitation and potential evapotranspiration data of the current month, and the n-month time scale is calculated using the sum of the cumulative precipitation and potential evapotranspiration data of the current month and the forward continuation of the n-1 months. The SPEI values for spring (March–May), summer (June–August), autumn (September–November), and winter (December–February) in the study area reflect the seasonal aridity conditions. The SPEI values for the 6-month time scale (SPEI-6) reflect the bi-annual aridity conditions in the study area, and the SPEI values for the 12-month time scale (SPEI-12) reflect the inter-annual variations in aridity.
SPEI has been extensively applied across various fields, including climate change analysis, environmental quality assessment, and aridity evaluation [37]. The core principle of SPEI is to gauge the level of aridity in a region by determining the difference between precipitation and potential evapotranspiration. The deviation of this difference from the mean value indicates the degree of aridity in the area. The formulation of SPEI typically involves three main steps: the calculation of potential evapotranspiration, the calculation of the difference between precipitation and actual evapotranspiration, and the normalization of the data, which are calculated using the following Formulas (1)–(12).
P E T = 16.0 × 10 T i H A
H = i = 1 12   T i 5 1.514
D i = P i P E T i
where PET is the potential evapotranspiration; T i is the mean monthly air temperature; T i is the Annual Heat Index; and A is a function of H , A = 0.49 + 0.179 H − 0.0000771 H 2 + 0.000000675 H 3 . Potential evapotranspiration in this study uses the data sets that have been obtained (Table 1). D i is the difference between precipitation and potential evapotranspiration. The normalization of the sequence D i can be calculated by Equations (4)–(7):
F x = 1 + α x γ β 1
α = ω 0 2 ω 1 β Γ 1 + 1 β Γ 1 1 β
β = 2 ω 1 ω 0 6 ω 1 ω 0 6 ω 2
γ = ω 0 α Γ 1 + 1 β Γ 1 1 β
where F x is the probability distribution function (log-logistic); G is the factorial function; and ω 0 , ω 1 , and ω 2 are the probability-weighted moments of D i , computed as Equations (8) and (9):
ω s = 1 N i = 1 N   ( 1 F i ) s D i
F i = i 0.35 N
where N denotes the count of months included in the calculations. Following this, the subsequent stage is the normalization of the accumulated probability density, a crucial step for ensuring the data’s standardization:
P = 1 F x
when the cumulative probability P ≤ 0.5:
ω = 2 ln P
  S P E I = ω c 0 + c 1 ω + c 2 ω 2 1 + d 1 ω + d 2 ω 2 + d 3 ω 3
where the SPEI is the Standard Precipitation Evapotranspiration Index. c 0 , c 1 , c 2 , d 1 , d 2 , and   d 3   have the values of 2.515517, 0.802853, 0.010328, 1.432788, 0.189269, and 0.001308, respectively. In this study the SPEI was constructed using Arcgis 10.8 and Matlab 2021b. The criteria for classifying the degree of aridity according to the SPEI values are shown in Table 2.

2.3.2. Remote Sensing Ecological Index

Over the last 20 years of the study period, data from 2000, 2005, 2010, 2015, and 2020 were selected for this study and analyzed to produce the results. The RSEI was calculated in two main steps [35].
Firstly, the NDVI, WET, LST, and NDBSI are calculated, which represent easily interpretable greenness, humidity, heat, and dryness, respectively. Finally, the RSEI is calculated as Equation (13):
R S E I = f N D V I , W E T , L S T , N D B S I
A higher RSEI indicates higher environmental quality. The RSEI methodology improves the practicality of environmental assessment indicators by reducing the impact of subjective human judgment in assigning factor values. It also utilizes a significant amount of medium-resolution satellite data, effectively broadening its scope.
Since 2002, NDVI has been widely used for vegetation characterization, such as regional vegetation cover, plant biomass, vegetation cover, etc. It is also an important ecological indicator [38,39,40]. The calculation formula is as follows, where ρ R e d and ρ N I R   are the spectral reflectance in the infrared and near-infrared bands, respectively.
N D V I = ρ N I R ρ R e d ρ N I R + ρ R e d
Temperature plays a central role in influencing plant growth and is also one of the key indicators of environmental change. In this study, grayscale values from remote sensing data were converted to degrees Celsius to determine the distribution of Land Surface Temperature across the study area. To determine this distribution, the Land Surface Temperature (LST) was calculated using the following formula:
L S T = 0.02 × D N S 273.15
where DNS is the gray value of the surface temperature image. Studies have indicated that the humidity component derived from MODIS’s Tasseled Cap Transformation (TCT) effectively mirrors the combined moisture of soil and vegetation in the QTP. Consequently, this study employs a method that multiplies the Tasseled Cap Transformation’s Humidity Index with surface reflectance to calculate WET, with the formula for this calculation being as follows:
W E T = 0.1147 ρ 1 + 0.2489 ρ 2 + 0.2408 ρ 3 + 0.3132 ρ 4 0.3122 ρ 5 0.6416 ρ 6 0.5087 ρ 7
where ρ i (i = 1, 2, 3, …, n) is the surface reflectance product reflectance for each band.
As urbanization advances, the expanse of impervious surfaces grows annually, exacerbating the environmental damage inflicted by economic activities. The increasing area of exposed natural soils also plays a role in intensifying surface aridity. For this study, we determine the Normalized Difference Built-up Soil Index (NDBSI) by applying a weighted mean of the Bare Soil Index (SI) and the Impervious Surface Index (IBI), as illustrated in the formula presented in Equations (17)–(19).
S I = ρ S W I R 1 + ρ R e d ρ B l u e + ρ N I R ρ S W I R 1 + ρ R e d + ρ B l u e + ρ N I R
I B I & = 2 ρ S W I R 1 ρ S W I R 1 + ρ N I R ρ N I R ρ N I R + ρ R e d + ρ G r e e n ρ G r e e n + ρ S W I R 1 2 ρ S W I R 1 ρ S W I R 1 + ρ N I R + ρ N I R ρ N I R + ρ R e d + ρ G r e e n ρ G r e e n + ρ S W I R 1
N D B S I = S I + I B I 2
where ρ S W I R 1 , ρ R e d , ρ B l u e , ρ G r e e n , and ρ N I R are the spectral reflectance in the short-wave infrared, infrared, extra-blue, extra-green, and near-infrared bands, respectively. In ArcGIS, the RSEI can be classified into five classes by equidistance, which are best, good, fair, poor, and worst, respectively [41].
In this step, all four indices, i.e., the NDVI, LST, WET, and NDBSI, are exemplified by the year 2000. The data sources for these indices are MOD13A1 for the NDVI, MOD11A2 for LST, and MOD09A1 for both the WET and NDBSI. The acquisition intervals for these images are 16 days for MOD13A1 and 8 days for MOD11A2 and MOD09A1. All MODIS satellite imagery is captured during the growing season of the QTP, which spans from July to September. After undergoing noise reduction, geometric correction, image enhancement, and grayscale processing, we obtained 4 NDVI, 7 LST, and 7 each of WET and NDBSI images for the year 2000. The average values of these indices represent the composite measure for that period. Following the acquisition of these averages, Principal Component Analysis (PCA) was employed to compute the RSEI for the year in question. After calculating the mean values for the four indices, the RSEI for each year was determined through Principal Component Analysis (PCA). This process was repeated for the years 2000, 2005, 2010, 2015, and 2020, amounting to a total of five RSEI computations over the study period. During this time, each image was cloud-masked, and the median value was taken for analysis. The cloud cover rates for the respective years were as follows: 2.37% in 2000, 3.47% in 2005, 2.00% in 2010, 0.17% in 2015, and 0.30% in 2020. All steps were re-implemented within the Google Earth Engine (GEE) platform, where the operator is required to input the vector boundaries of the study area. Upon completion of the process, the spatial characteristics and distribution of the RSEI are obtained, along with the proportional contributions of each index to the annual calculation.

2.4. Mann–Kendall Test

The nonparametric Mann–Kendall test, developed by Mann and Kendall, first saw publication in the year 1945, as a technique for detecting and forecasting time-series data [42]. It is widely used in environmental science, climatology, and other fields. The Mann–Kendall test can be used to determine whether a mutation has occurred in a climatic sequence and, if so, when such a mutation occurred [43]. The Mann–Kendall test is a widely recognized tool for detecting trends in time series data attributed to climate change. As a non-parametric method, it has the advantage of not assuming any specific distribution for the sample data, thereby being less susceptible to the influence of outliers. Consequently, it is particularly well suited for analyzing both categorical and sequential variables [44].
The formulation of the Mann–Kendall test is based on a set of Equations (20)–(23):
s g n A i A j = 1 , A i A j > 0 0 , A i A j = 0 1 , A i A j < 0
S = i = 2 n   j = 1 i 1   s g n A i A j
V a r S = n n 1 2 n + 5 18
The Mann–Kendall test incorporates a step function, represented by “sgn”, which is applied to the time series data. In this study, the forest carbon stock represents the time series variable, denoted as A i . The test statistic follows a normal distribution, indicated by “ S ”. The “ V a r ” term refers to the variance calculation, which is an essential part of the test’s methodology. Lastly, “ n ” signifies the sample size of the time series data under examination.
Once the variance calculation and the normal distribution test have been conducted, the standardized test statistic, denoted as Z:
Z = S 1 V a r S , S > 0 0 , S = 0 S + 1 V a r S , S < 0
At a given confidence level α, the original hypothesis is invalid if | Z | ≥ Z 1 − α/2. If Z > 0, the extent of fluctuation within the time series expands; if Z < 0, the range of variation of the time series decreases, then it indicates a decrease in the range of variation. In general, if | Z | ≥ 2.58, the trend of change is deemed to be of extreme significance. If | Z | ≤ 2.58, the trend is deemed significantly pronounced. If 1.96 ≤ | Z | < 2.58, the trend is regarded as statistically significant. If 0 ≤ | Z | < 1.96, the trend is considered insignificant. In this study, the ecological quality of the QTP region was assessed using the Mann–Kendall method, with calculations performed in ArcGIS Pro version 3.0.

2.5. Pearson Correlation Analysis

Using the Pearson correlation coefficient method, the correlation analysis between the Remote Sensing Ecological Index (RSEI), the Aridity Index, and the SPEI of each pixel was carried out based on the pixel scale [45]. Pearson’s correlation coefficient is employed to quantify the degree of linear association between a pair of independent variables [46]. The formula for its calculation is as follows:
r = i = 1 n   X i X ¯ Y i Y ¯ i = 1 n   X i X ¯ 2 i = 1 n   Y i Y ¯ 2
In this context, r signifies the correlation coefficient, while X i and Y i denote the two variables under consideration, and X ¯ and Y ¯ represent the average values of the respective variables. r is distributed in the interval [–1, 1]. In this study, we employed Pearson’s correlation analysis to investigate the relationship between forest carbon stocks and various environmental elements. The strength of a positive correlation is indicated by values approaching 1, while a strong negative correlation is signified by values near −1. Values within the range of −1 to 1 reflect the degree of correlation between the variables [47]. The computational phase of this analysis was conducted using Matlab 2021b, ensuring a rigorous examination of the correlation coefficients. When conducting correlation analysis on long-term time-series raster data with Matlab 2021b, it is crucial to be mindful of the varying resolutions present in the datasets. Data must undergo preprocessing to align the dimensions, ensuring consistency in the number of rows and columns across all frames before the analysis can be initiated.

2.6. Hurst Index Analysis

Self-similarity and long-range dependence are common phenomena in nature and are widely used in hydrology, climate, geology, and other fields. British scholar Hurst proposed the Hurst Index by studying the changes of the Nile flow time series, which is used to quantitatively characterize the continuity or long-range dependence of the time series. The Hurst Index is estimated by a variety of methods, generally using the R/S analysis method, which is more reliable. In recent years, R/S analysis has been widely used in the fields of environmental change and geoscience.
For the time-series RSEI, i = 1, 2, 3, …, n, and the time series is defined as follows:
R S E I τ = 1 τ 1 τ   R S E I τ τ = 1 ,   2 , , n
X τ = ι = 1 τ ( NDVI ι N D V ¯ I τ ) 1 ι τ
R τ = m a x 1 t τ   X τ m i n 1 t τ   X τ τ = 1 ,   2 , , n
S τ = 1 τ t = 1 τ   R S E I τ R S E I τ 2 τ = 1 ,   2 , , n
For the ratio R ( τ ) / S ( τ ) R / S , a Hurst phenomenon exists in the time series if the following relationship exists R / S τ H . H is the Hurst exponent, which can be obtained by using least-squares fitting according to l o g ( R / S ) n = a + H × l o g ( n ) . The size of the H value can judge the persistence of the RSEI sequence, including three forms: 0.5 < H < 1, suggesting that the time series exhibits characteristics of long-term correlation within a persistent sequence, with higher values of H approaching 1 indicating greater persistence; H = 0.5, then the RSEI sequence is a stochastic sequence; 0 < H < 0.5, then the RSEI sequence has an anti-persistence, and the closer H is to 0, the stronger the anti-persistence.

3. Results

3.1. Characterization of the Spatial and Temporal Distribution of SPEI

The annual mean Standardized Precipitation Evapotranspiration Index (SPEI) values exhibit varying patterns across different time scales as depicted in Figure 2. For the SPEI-1 (as shown in Figure 2a), the values oscillate between −0.17 and 0.12, averaging at −0.07. Moving on to the SPEI-3 (Figure 2b), the range of fluctuation is from −0.12 to 0.30, with an average of 0.03. The SPEI-6 (Figure 2c) shows a broader fluctuation, spanning from −0.09 to 0.43, and has an average value of 0.20. Further extending the time scale, the SPEI-9 (Figure 2d) demonstrates a wider range, oscillating between −0.16 and 0.65, with the mean settling at 0.18. The SPEI-12 (Figure 2e) follows suit with a fluctuation from −0.25 to 0.75, averaging at 0.16. Finally, the SPEI-24 (Figure 2f), significantly influenced by long-term precipitation and potential evapotranspiration, displays the most extensive variation, ranging from −0.45 to 1.03, with an average value of 0.43.
Overall, as illustrated in Figure 3, the SPEI values of the six different time scales show more differentiated characteristics. One- and three-month-scale SPEI values show a more consistent trend, showing that the regional value of the Qaidam Basin is significantly higher than that of other regions, and the SPEI value of some regions in the south-east Hengduan Mountains is higher. Nine- and twelve-month-scale SPEI values show a more consistent trend, showing that the regional value of the Qaidam Basin is significantly higher, and some regions in the south-west are higher than other regions. The SPEI values at the 6-month scale show a more uniform distribution, with lower values in the eastern region, and the SPEI values at the 24-month scale are higher in the whole, with slightly lower values in the Qaidam Basin than in other regions.

3.2. Characteristics and Trends of Spatial and Temporal Distribution of RSEI

3.2.1. RSEI Components

The RSEI consists of four evaluation indicators, NDVI, LST, WET, and NDBSI, which represent the ease of greenness, heat, humidity, and dryness, respectively. Figure 4 presents the four indicators in time and space for the years 2000, 2005, 2010, and 2020, respectively. Over the past two decades, various environmental indices have shown distinct spatial and temporal trends. The NDVI has demonstrated a stable trend, with higher values consistently observed in the south-east, which decrease as one moves towards the north-west [48]. Similarly, the WET values have maintained a general stability, with a decline from the south-east to the north-west [49]. On the other hand, the LST has shown higher values near the northern and south-western regions, with a decrease towards the south-east. The NDBSI has exhibited a general decrease over the 20-year period, with lower values in the south-east and higher values in the north-west.

3.2.2. RSEI Spatial and Temporal Distribution

The RSEI, categorized into five classes of best, good, normal, bad, and worst according to the equidistance values in ArcGIS (Figure 5a–e), shows the spatial and temporal distribution characteristics of the RSEI in the QTP region. Figure 5f–j is the total area of best, good, normal, bad, and worst levels in each study period.
Overall, the trend analysis, along with the spatial and temporal attributes of the RSEI, revealed that the ecological status of the QTP was predominantly characterized by normal and best. The normal and best land areas were 2,214,760 km2 in 2015 and reached a maximum of 2,365,378 km2 in 2000 (Figure 5f). In the rest of the years, the proportion of land characterized by normal and best land area was stable at about 70%. The analysis implies that the ecological state of the QTP is generally favorable.
Over the past two decades, the spatial distribution of the RSEI across the QTP has exhibited a relatively consistent pattern, demonstrating a gradient of decline from the south-eastern to the north-western regions (Figure 5). The areas with “bad” and “worst” RSEI ratings are mainly located in the south-western and northern parts of the QTP, and the areas with “good” and “normal” RSEI ratings are mainly located in south-eastern Tibet, Central Yushu Prefecture, and northern Golog, while the RSEI “best” regions are mainly located in the eastern and southern parts of the Plateau.
The northern region has the largest and longest duration of worst in 20 years. The region is the Qaidam Basin, with high surface temperatures and scarce surface water resources. The RSEI indicates that the area is ecologically vulnerable and necessitates the implementation of suitable conservation strategies to safeguard against ecological decline. In terms of temporal characteristics, the proportion of bad area decreased from 160,422 k m 2 in 2000 (Figure 5f) to 146,297 km2 in 2005 (Figure 5g) and has remained below 5% since then, which indicates that the management of the area in the past 20 years has been effective.
The RSEI reveals that areas classified as normal and bad are largely concentrated in the south-west. Between 2000 and 2005, there was a noticeable decrease in the proportion of bad areas, from 15.1%, as illustrated in Figure 5f, down to 13.92%, as shown in Figure 5g. This suggests a positive trend, reflecting an improvement in the regional ecological conditions over the specified timeframe. Nevertheless, the share of undesirable areas rose to 20.11% by 2010, as depicted in Figure 5h, and has since remained relatively constant. This stabilization suggests a period of ecological decline in the region commencing from 2010. The areas rated as best and good for land quality are predominantly situated in the south-eastern part of the region. In 2010, the best category constituted the smallest proportion at 0.02%, as indicated in Figure 5h. By 2015, the proportion of the best land category reached its zenith at 6.67%, as illustrated in Figure 5i. Conversely, the good category hit its lowest point in 2010 with 27.58%, according to Figure 5h, yet peaked in 2005 at 43%, as documented in Figure 5g. The higher RSEI values in the south-eastern part of the QTP indicate that the local ecological conditions are better; this can be attributed to the lower altitude in the region possibly. The area experiences limited human activity and boasts favorable hydrothermal conditions that are conducive to vegetation growth. The good growth of vegetation can in turn nourish water and better feed back to the ecosystem [50]. In general, the ecological quality of the QTP has been improving over the past 20 years, but deterioration exists in some areas.

3.2.3. RSEI Spatial Trends

Figure 5 shows that there are large spatial and temporal differences in the RSEI (Figure 5a–j), and changes in ecological conditions cannot be obtained by analyzing only its multi-period status. Therefore, GIS pro was introduced to treat long time series of remote sensing images as trend analysis. Consequently, the spatio-temporal variations of the RESI were delineated and are depicted in Figure 6.
From the spatial trend of the RSEI in Figure 6a, the QTP can be divided into three parts, i.e., western, central, and eastern. The western region encompasses the largest expanse, and, over the last two decades, the RSEI has exhibited a consistent increasing trend while most of the eastern region showed a decreasing trend. Based on the contents of Table 3, it can be seen that the areas where the RSEI has shown a very significant increase account for 1.55%, predominantly located in the western region. In contrast, areas with non-significant growth constitute 36.62% and are found in the same regions as those with significant increases. The sum of the area with a non-significant decrease and highly significant decrease in RSEI was about 61.73%, which was mainly located in the eastern region. Overall, the ecological environment quality of the QTP is exhibiting a trend of divergence. From the perspective of the significance of the results, most areas do not show significant changes, while a minority of regions in the west have demonstrated a clear upward trend in the RSEI.
During the study period, the area around Qinghai Lake in the north-eastern part emerged as the region with the most severe decline in the RSEI within the entire study area, which is part of the Qaidam Basin with high surface temperatures and scarce surface and groundwater resources. Groundwater resources are scarce; thus, vegetation is scarce, and most of the land is barren [51], resulting in a declining trend of RSEI values in the region over the past 20 years. Furthermore, given the region’s ecological fragility, it necessitates the implementation of appropriate conservation measures to prevent further ecological degradation.

3.3. SPEI and RSEI Correlation at Different Scales

The annual correlation and significance of the SPEI (SPEI-1, SPEI3, SPEI-6, SPEI-9, SPEI-12, and SPEI-24) and RSEI at different time scales, based on the likelihood meta-reflecting the degree of the ecological quality of the QTP region affected by meteorological aridity, were studied.
The correlation and significance test results of different time-scale SPEIs (SPEI-1, SPEI3, SPEI-6, SPEI-9, SPEI-12, and SPEI-24) with RSEI based on the image metrics reflecting the degree of ecological environmental quality conditions affected by meteorological aridity are reflected in Figure 7. The results show that the spatial distribution of the correlation between SPEI and RSEI at 1, 3, 6, 9, 12, and 24-month time scales is relatively consistent, short-term aridity has less impact on ecological environment quality due to the presence of certain moisture in the soil, and vegetation can take up water from the soil through the root system, while long-term aridity has a greater impact on ecological environment quality.
Furthermore, Figure 7 shows that the 1, 3, 6, 9, 12, and 24-month time-scale SPEI and RSEI, which showed E-S-N-C and S-N-C, accounted for 18.83%, 21.35%, 18.18%, 19.88%, 17.12%, and 18.27% of the total study area, respectively, and the negatively correlated areas were more distributed in the south-western part of the study area. Additionally, Figure 7 shows that the area is most of the ecological quality classified as bad and worst, indicating that the ecological quality of these areas is negatively correlated with aridity conditions, which may be due to the fact that there is more glacier snow in the area. Furthermore, it shows that the occurrence of meteorological aridity will increase the glacier snow melt, which will be conducive to the growth of vegetation to promote the improvement of the quality of the ecological environment. The positively correlated regions are smaller and more scattered, and the positively correlated regions of the SPEI and RSEI at 1, 3, 6, 9, 12, and 24 time scales are mainly concentrated in the south-eastern, south-western, and northern regions of the QTP, where the sum of the S-P-C and the E-S-P-C accounted for 20.89%, 17.41%, 13.44%, 11.53%, 12.08%, and 15.56%, respectively, suggesting that the ecological integrity of these areas is significantly compromised by meteorological aridity.

3.4. Ecological Quality Sustainability Changes

Figure 8 shows the Hurst Index of the QTP. The mean value is 0.47, and the Hurst Index of the RSEI of the QTP is greater than 0.5 in 52.03% of the regions and less than 0.5 in 47.97% of the regions. This indicates that the future trend of ecological quality on the QTP is stronger in the same direction than in the opposite direction, suggesting that the future change in ecological quality on the QTP is weakly persistent, i.e., the trend of the RSEI change in most regions is the same as that in the past.
The spatial variation of the Regional Sustainable Environment Index (RSEI) and the Hurst Index provides insights into the future ecological trends (Figure 6). Regions showing sustained improvement in ecological environment quality are primarily located in the western part of the Qinghai–Tibet Plateau (QTP), while areas of sustained degradation are found in the eastern and south-eastern parts. The south-west of the study area exhibits a clear persistence in ecological environment quality, which is in stark contrast to the non-continuous characteristic observed in the south-eastern part, where human activities have had a significant impact.

4. Discussion

The evolution of ecosystem patterns in the QTP region has been caused by a variety of factors, with natural and anthropogenic factors being the main contributors to the changes, and the SPEI, which is a meteorological Aridity Index that takes into account the effects of the dispersal of signs, has been shown by many scholars to be suitable for characterizing the aridity in the study area [52,53,54]. It has been found that temperature is the most influential factor on the rate of change in ecological patterns in the QTP; therefore, the SPEI was chosen as a driver of ecological quality in this study to investigate the response of aridity events to changes in ecological quality in the QTP region. Six SPEI with different time scales were selected to reflect the aridity conditions in the QTP region on monthly, quarterly, and annual scales, etc. The study findings are largely in accordance with the conclusions from studies monitoring the ecological environment of the QTP using other aridity indices. This indicates that using the SPEI as an influencing factor to reflect changes in the ecological environment quality of the QTP region holds a high degree of credibility.
The QTP is a unique ecosystem where altitudinal gradients, temperature, and aridity lead to the slow growth and uneven distribution of natural vegetation. The spatial patterns observed in the RSEI indices are influenced by various geographical and climatic factors. The lower elevation in the south-east contributes to the favorable water and heat conditions that promote better vegetation growth, hence the higher NDVI values. As elevation increases towards the north-west, hydrothermal conditions become less conducive to vegetation development, resulting in lower NDVI values. WET values are higher in the south-east due to the influence of warm and humid monsoons, lower terrain, and higher precipitation, which facilitates water vapor accumulation and higher humidity. In contrast, the higher terrain, lower rainfall, and lower temperatures in the north-west result in lower WET values. NDBSI values are lower in the south-east due to the influence of the monsoon, which brings high rainfall, a humid climate, dense vegetation, and good ground cover. In the north-west, the dry westerly winds, low rainfall, and arid climate result in sparser vegetation and poorer surface cover, leading to higher NDBSI values.
As shown in the Results section, the distribution pattern of the RSEI is higher in the south-east, decreasing in the north-west, and is lowest in the Qaidam Basin. The reason for the low RSEI values in the Qaidam Basin region may be that the location is in the north-eastern part of the QTP, with low relief, arid climate, scanty precipitation, and scarce water resources, and the arid climatic conditions make the vegetation growth limited and the soil susceptible to erosion and salinization phenomena. Secondly, the long-term unsustainable development and utilization patterns influenced by human activities have led to serious problems such as land degradation and water resource pollution. In conclusion, the Qaidam Basin’s poor ecological environment is largely the result of a confluence of factors, including aridity, human-induced influences, and natural calamities. The higher RSEI values in the south-eastern part of the QTP indicate that the local ecological conditions are better; this is probably due to the lower altitude in the region. The area experiences less human interference, and the conditions of water and heat are favorable for the growth of vegetation. The good growth of vegetation in turn can contain water and better feedback to the ecosystem. Overall, the ecological quality of the QTP has been improving over the past 20 years, but there is deterioration in some areas, which is consistent with the general knowledge in this study field [53,55,56].
From the spatial scale, the negative correlation between the SPEI and RSEI in the QTP over the past two decades is more distributed in the south-western region, which is due to the fact that the warming and humidification of the climate of the QTP has led to a decrease in the Aridity Index of the region, and the climate is humid and the temperature is higher, which is favorable for the growth of vegetation, so the ecological quality of the environment in the south-west is improved. The positive correlation in the south-eastern part of the region is due to the fact that the eastern and south-eastern parts of the low-altitude region are topographically enclosed by high mountains. The climate in the region, shaped by the south-east monsoon, is characterized by abundant rainfall and extended periods of sunshine. This results in higher precipitation levels that surpass evaporation rates, leading to more effective soil moisture retention. The area benefits from a warmer environment that is propitious for the growth of vegetation. Prior research has established the correlation between climatic factors and the ecological quality, particularly in relation to the flourishing of vegetation in the south-east region. Overall, topographic and climatic factors led to better vegetation recovery, a higher RSEI, better environmental quality, and stronger response to aridity changes in the south-eastern QTP, while vegetation growth in the central and western regions was slow, and the response to aridity changes was weaker.
The persistence of ecological improvement in the south-west of the QTP can be attributed to the region’s susceptibility to climate change. With global warming, the temperature in this high-altitude area is expected to rise, potentially altering precipitation patterns. This climate warming could lead to the melting of alpine ice and snow, thereby increasing water resource availability, which is beneficial for ecological restoration and improvement in the region. Conversely, the non-continuous RSEI trend in the south-eastern part is likely a result of intensified human activities. These activities have disrupted the local ecological balance and biodiversity, contributing to ecosystem fragility and an increase in the Hurst Index value, indicating a less predictable and more erratic ecological state in the area.
Given the results presented previously, it can be concluded that the spatial distribution of the correlation between the SPEI and RSEI in the QTP region at the 1, 3, 6, 9, 12, and 24-month time scales is basically the same, and short-term aridity has a smaller impact on the quality of the ecological environment because there is a certain amount of water in the soil. Moreover, the vegetation can absorb the water from the soil through the root system, whereas long-term droughts have a larger impact on the quality of the ecological environment. This is because, under short-term aridity conditions, a certain amount of water may still remain in the soil, which provides the necessary water source for vegetation to absorb and sustain growth through the root system. This buffering effect of soil moisture mitigates the direct effects of short-term drought on vegetation. At the same time, vegetation may have certain adaptations to drought conditions. For example, some plants may develop deeper root systems or more efficient water uptake and utilization mechanisms in response to arid conditions. Such adaptations allow vegetation to withstand the effects of aridity in the short term.
In contrast, long-term aridity is cumulative and has a greater impact on ecological quality due to its longer duration. Climate change may lead to changes in precipitation patterns, affecting the frequency and intensity of aridity events. Prolonged aridity may be associated with long-term changes in regional climate patterns, which may have far-reaching impacts on ecosystems.
The QTP is a relatively complex ecosystem, with alpine rocks as the main terrain type in the central and western part of the region, which is susceptible to factors such as wind speed and sunshine, and thus the correlation between the RSEI and Aridity Index is relatively weak. Therefore, meteorological factors such as sunshine hours and average wind speed can be added as driving factors in the future study direction so that the changes of ecological environment quality in the QTP region can be analyzed more comprehensively. The study calculated the Aridity Index only from the perspective of meteorological data and did not consider hydrological, agricultural, or socio-economic factors. In order to further strengthen the rigor of the study, a comprehensive Aridity Index that integrates information from multiple sources can be used, which can more accurately monitor the occurrence of meteorological aridity and provide technical support for a more objective and comprehensive assessment of the extent and scope of the impact on the quality of the ecosystem and the reduction in the losses incurred.
In this study, five years of data were used for the M–K test and Hurst Index analyses, given the limitations of time and resolution access. This temporal resolution of the data may have had some impact on the precision and sensitivity of the analyses. Future research will aim to obtain annual data with higher temporal resolution for analysis with a view to obtaining more robust analytical results.
In addition, the composition of ecosystems in different study areas has different ecosystems with different natural characteristics, and different natural factors can be selected to construct ecological environment evaluation models in future studies, so it is suggested that the natural characteristics of land-use types can be combined to select appropriate ecological factors in ecological environment evaluation so as to construct a more diversified ecological environment evaluation model to better analyze the ecological environment. Overall, the use of the SPEI facilitates the evaluation and quantification of ecological quality on the QTP and improves our understanding of aridity change as it relates to the ecosystem.

5. Conclusions

In this study, the response of meteorological aridity to ecological environment quality was studied, and the temperature and precipitation data of the QTP from 2000 to 2020 were selected. By calculating the SPEI indices on the time scales of 1, 3, 6, 9, 12, and 24 months, the spatial and temporal distributions of aridity characteristics on the QTP, as well as the ecological quality of the QTP region and its trends from 2000 to 2020, were analyzed. The results were also analyzed by overlaying the results with the correlation and significance of the SPEI at different time scales, respectively, for evaluating the degree of the response of meteorological aridity to the quality of ecological environment. The main conclusions of the article are summarized as follows.
  • The Aridity Index of the QTP shows a differentiated spatial and temporal distribution pattern. At shorter time scales, such as 1, 3, 9, and 12 months, the distribution of the SPEI is higher in the north-east and part of the south-west and lower in the center. Monthly, quarterly, and annual scales of SPEI values show an increasing trend, indicating that the QTP has shown a non-significant trend of becoming drier in the last two decades.
  • The ecological environmental quality of the QTP has a strong spatial distribution pattern. The ecological environment quality in the south-east is significantly better, while the Qaidam Basin and the west have lower ecological environment quality indexes, but the overall trend of environmental quality is getting better. During the study period, the changes in ecological environmental quality in the QTP region were the result of a combination of natural and human factors. Therefore, in the future, the implementation of ecological projects in each region should be continued, with a focus on areas where the ecological quality has deteriorated and where it may deteriorate in the future in order to further optimize and improve the ecological protection scheme.
  • There is a correlation between aridity conditions and ecological environment quality in the QTP. The areas with significant positive correlation between the RSEI and SPEI are mainly concentrated in the south-east, south-west and north of the QTP, indicating that the ecological environment quality in these areas is more seriously affected by meteorological aridity. At the same time, the correlation between the SPEI and RSEI increased with the increase in aridity time scale.

Author Contributions

Conceptualization, Y.Y. and J.C.; methodology, Y.Y. and J.C.; software, Y.Y., Y.G. and J.C.; validation, S.W., X.H. and Y.Y.; formal analysis, Y.Y.; investigation, X.L. and J.C.; resources, Y.Y. and J.C.; data curation, Y.Y.; writing—original draft preparation, Y.Y.; writing—review and editing, Y.Y. and J.C.; visualization, Y.Y.; supervision, Y.Y. and Y.H.; project administration, S.W.; funding acquisition, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Project of Industry University Research cooperation between Tsinghua University and China Forestry Group Corporation on Forestry carbon sink development (ZLJT-THU2022110101) and the National Natural Science Foundation (72104118).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fava, F.; Vrieling, A. Earth Observation for Drought Risk Financing in Pastoral Systems of Sub-Saharan Africa. Curr. Opin. Environ. Sustain. 2021, 48, 44–52. [Google Scholar] [CrossRef]
  2. Dikshit, A.; Pradhan, B.; Huete, A. An Improved SPEI Drought Forecasting Approach Using the Long Short-Term Memory Neural Network. J. Environ. Manag. 2021, 283, 237–248. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, F.; Wang, X.; Huang, J.; Ye, J.; Zhong, J. Analysis of Spatial-Temporal Characteristics of Drought and Flood in the Zhejiang Province Based on GIS. Bull. Sci. Technol. 2009, 14, 178–190. [Google Scholar]
  4. Wang, F.; Wang, Z.; Yang, H.; Zhao, Y.; Zhang, Z.; Li, Z.; Hussain, Z. Copula-Based Drought Analysis Using Standardized Precipitation Evapotranspiration Index: A Case Study in the Yellow River Basin, China. Water 2019, 11, 1298. [Google Scholar] [CrossRef]
  5. Chen, Q.; Timmermans, J.; Wen, W.; van Bodegom, P.M. Ecosystems Threatened by Intensified Drought with Divergent Vulnerability. Remote Sens. Environ. 2023, 289, 113512. [Google Scholar] [CrossRef]
  6. Liu, J.; Ji, Y.-H.; Zhou, G.-S.; Zhou, L.; Lyu, X.-M.; Zhou, M.-Z. Temporal and Spatial Variations of Net Primary Productivity (NPP) and Its Climate Driving Effect in the Qinghai-Tibet Plateau, China from 2000 to 2020. J. Appl. Ecol. 2022, 33, 1533–1538. [Google Scholar]
  7. Cui, X.-L.; He, Y.-L.; Li, Z.-S. Spatial-Temporal Variation of Vegetation Water Use Efficiency and Its Relationship with Climate Factors over the Qinghai-Tibet Plateau, China. J. Appl. Ecol. 2022, 33, 1525–1532. [Google Scholar]
  8. Xia, M.; Jia, K.; Zhao, W.; Liu, S.; Wei, X.; Wang, B. Spatio-Temporal Changes of Ecological Vulnerability across the Qinghai-Tibetan Plateau. Ecol. Indic. 2021, 123, 320–333. [Google Scholar] [CrossRef]
  9. Wei, Y.; Lu, H.; Wang, J.; Wang, X.; Sun, J. Dual Influence of Climate Change and Anthropogenic Activities on the Spatiotemporal Vegetation Dynamics Over the Qinghai-Tibetan Plateau From 1981 to 2015. Earth Future 2022, 10, 467–481. [Google Scholar] [CrossRef]
  10. Cao, J.; Wu, E.; Wu, S.; Fan, R.; Xu, L.; Ning, K.; Li, Y.; Lu, R.; Xu, X.; Zhang, J. Spatiotemporal Dynamics of Ecological Condition in Qinghai-Tibet Plateau Based on Remotely Sensed Ecological Index. Remote Sens. 2022, 14, 4234. [Google Scholar] [CrossRef]
  11. Peng, X.; Zhang, S.; Peng, P.; Chen, A.; Li, Y.; Wang, J.; Bai, M. Unraveling the Ecological Tapestry: A Comprehensive Assessment of Changtang Nature Reserve’s Ecological and Environmental Using RSEI and GEE. Land 2023, 22, 2001. [Google Scholar] [CrossRef]
  12. Zhao, Z.; Li, T.; Zhang, Y.; Lu, D.; Wang, C.; Lu, Y.; Wu, X. Spatiotemporal Patterns and Driving Factors of Ecological Vulnerability on the Qinghai-Tibet Plateau Based on the Google Earth Engine. Remote Sens. 2022, 14, 5279. [Google Scholar] [CrossRef]
  13. McKee, T.B.; Doesken, N.J.; Kleist, J. The Relationship of Drought Frequency and Duration to Time Scales. In Proceedings of the 8th Conference on Applied Climatology, Boston, MA, USA, 17–22 January 1993; Volume 34, pp. 90–104. [Google Scholar]
  14. Changnon, S.A.; Easterling, W.E. Measuring Drought Impacts: The Illnois Cace. J. Am. Water Resour. Assoc. 1989, 364, 2350–2366. [Google Scholar]
  15. Skøien, J.O.; Blöschl, G.; Western, A.W. Characteristic Space Scales and Timescales in Hydrology. Water Resour. Res. 2003, 39, 756–777. [Google Scholar] [CrossRef]
  16. Hao, Z.; Singh, V.P. Drought Characterization from a Multivariate Perspective: A Review. J. Hydrol. 2015, 527, 668–678. [Google Scholar] [CrossRef]
  17. Stagge, J.H.; Tallaksen, L.M.; Gudmundsson, L.; Van Loon, A.F.; Stahl, K. Candidate Distributions for Climatological Drought Indices (SPI and SPEI). Int. J. Climatol. 2015, 553, 1289–1293. [Google Scholar] [CrossRef]
  18. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef]
  19. Zarei, A.R.; Mahmoudi, M.R. Assessment of the Effect of PET Calculation Method on the Standardized Precipitation Evapotranspiration Index (SPEI). Arab. J. Geosci. 2020, 13, 182. [Google Scholar] [CrossRef]
  20. Han, X.; Qiu, X.H. A Remote Sensing Index for Assessment of Regional Ecological Changes. China Environ. Sci. 2013, 33, 889–897. [Google Scholar]
  21. Peng, L.; Zhang, J.; Liang, E.; Hu, M. Evaluation of Natural Ecological Environment Change in Manasi River Basin Based on RSEI. J. Shihezi Univ. 2017, 35, 506–512. [Google Scholar]
  22. Sun, C.; Li, J.; Liu, Y.; Cao, L.; Zheng, J.; Yang, Z.; Ye, J.; Li, Y. Ecological Quality Assessment and Monitoring Using a Time-Series Remote Sensing-Based Ecological Index (Ts-RSEI). GISci. Remote Sens. 2022, 59, 1793–1816. [Google Scholar] [CrossRef]
  23. Cui, R.; Han, J.; Hu, Z. Assessment of Spatial Temporal Changes of Ecological Environment Quality: A Case Study in Huaibei City, China. Land 2022, 11, 257–267. [Google Scholar] [CrossRef]
  24. Kuang, Q.; Yuan, Q.; Han, J.; Leng, R.; Wang, Y.; Zhu, K.; Lin, S.; Ren, P. A Remote Sensing Monitoring Method for Alpine Grasslands Desertification in the Eastern Qinghai-Tibetan Plateau. J. Mt. Sci. 2020, 17, 15. [Google Scholar] [CrossRef]
  25. Wang, G.; Han, L. The Vegetation NDVI Variation Trend in Qinghai-Tibet Plateau and Its Response to Climate Change. In Proceedings of the 2012 2nd International Conference on Remote Sensing, Environment and Transportation Engineering, Nanjing, China, 1–3 June 2012; IEEE: Piscataway, NY, USA, 2012; Volume 33, pp. 423–451. [Google Scholar]
  26. Gao, Q.; Guo, Y.; Xu, H.; Ganjurjav, H.; Li, Y.; Wan, Y.; Qin, X.; Ma, X.; Liu, S. Climate Change and Its Impacts on Vegetation Distribution and Net Primary Productivity of the Alpine Ecosystem in the Qinghai-Tibetan Plateau. Sci. Total Environ. 2016, 554, 34–41. [Google Scholar] [CrossRef] [PubMed]
  27. Peng, S.; Ding, X.; Liu, W. 1 km Monthly Temperature and Precipitation Dataset for China from 1901 to 2017. Earth Syst. Sci. Data 2019, 11, 1931–1946. [Google Scholar] [CrossRef]
  28. Peng, S. 1-km Monthly Precipitation Dataset for China (1901–2020); National Tibetan Plateau Data Center: Beijing, China, 2020. [Google Scholar] [CrossRef]
  29. Peng, S. 1-km Monthly Mean Temperature Dataset for China (1901–2022); National Tibetan Plateau Data Center: Beijing, China, 2019. [Google Scholar] [CrossRef]
  30. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS; NASA Technical Reports Server (NTRS); NASA: Washington, DC, USA, 1974; Volume 521, pp. 321–333. [Google Scholar]
  31. Crist, E.P.; Laurin, R.; Colwell, J.E.; Kauth, R.J. Investigations of Vegetation and Soils Information Contained in Landsat Thematic Mapper and Multispectral Scanner Data; NASA: Washington, DC, USA, 1984; Volume 127, pp. 34–41. [Google Scholar]
  32. Baig, M.H.A.; Zhang, L.; Shuai, T.; Tong, Q. Derivation of a Tasselled Cap Transformation Based on Landsat 8 At-Satellite Reflectance. Remote Sens. Lett. 2014, 5, 423–431. [Google Scholar] [CrossRef]
  33. Rogers, R.H. Important Prospect: Third Thematic Conference on Remote Sensing for Marine and Coastal Environments, to Be Held in Seattle, Washington, USA, from 18–20 September 1995. Environ. Conserv. 1994, 21, 376. [Google Scholar] [CrossRef]
  34. Chen, S.; Wen, Z.; Zhang, S.; Huang, P.; Ma, M.; Zhou, X.; Liao, T.; Wu, S. Effects of Long-Term and Large-Scale Ecology Projects on Forest Dynamics in Yangtze River Basin, China. For. Ecol. Manag. 2021, 496, 119463. [Google Scholar] [CrossRef]
  35. Moreno, R.; Ojeda, N.; Azócar, J.; Venegas, C.; Inostroza, L. Application of NDVI for Identify Potentiality of the Urban Forest for the Design of a Green Corridors System in Intermediary Cities of Latin America: Case Study, Temuco, Chile. Urban For. Urban Green. 2020, 55, 126821. [Google Scholar] [CrossRef]
  36. Kumar, V.; Bharti, B.; Singh, H.P.; Topno, A.R. Assessing the Interrelation between NDVI and Climate Dependent Variables by Using Granger Causality Test and Vector Auto-Regressive Neural Network Model. Phys. Chem. Earth Parts A B C 2023, 131, 103428. [Google Scholar] [CrossRef]
  37. Kimosop, P. Characterization of Drought in the Kerio Valley Basin, Kenya Using the Standardized Precipitation Evapotranspiration Index: 1960−2016. Singap. J. Trop. Geogr. 2019, 40, 239–256. [Google Scholar] [CrossRef]
  38. Chashina, B.; Ramazanova, N.; Atasoy, E. Analysis of the State of Plant Cover of the Berezovka River Basin of the West-Kazakhstan Region Based on NDVI. J. Geogr. Environ. Manag. 2021, 651, 245–258. [Google Scholar] [CrossRef]
  39. Ding, J.; Tiyip, T. Study on Oasis Vegetation Ecological Landscape Change Based on NDVI. Geogr. Territ. Res. 2002, 18, 23–26. [Google Scholar]
  40. Wagenseil, H.; Samimi, C. Assessing Spatio-temporal Variations in Plant Phenology Using Fourier Analysis on NDVI Time Series: Results from a Dry Savannah Environment in Namibia. Int. J. Remote Sens. 2006, 27, 3455–3471. [Google Scholar] [CrossRef]
  41. Yang, X.; Meng, F.; Fu, P.; Zhang, Y.; Liu, Y. Spatiotemporal Change and Driving Factors of the Eco-Environment Quality in the Yangtze River Basin from 2001 to 2019. Ecol. Indic. 2021, 131, 108214. [Google Scholar] [CrossRef]
  42. Mann, H.B. Nonparametric Tests against Trend. Econom. J. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  43. Hamed, K.H.; Rao, A.R. A Modified Mann-Kendall Trend Test for Autocorrelated Data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  44. Hamed, K.H. Trend Detection in Hydrologic Data: The Mann–Kendall Trend Test under the Scaling Hypothesis. J. Hydrol. 2008, 349, 350–363. [Google Scholar] [CrossRef]
  45. Zhang, Y.; Yan, D.; Ping-an, H. Quantitative Analysis of the Relationship of Biology Species Using Pearson Correlation Coefficient. Comput. Eng. Appl. 2005, 33, 79–82. [Google Scholar]
  46. Zhou, Y.; Feng, X.; Enhedelihai; Luo, T.; Yang, X.; He, M. Statistical Analysis of Fracture Properties Based on Particle Swarm Optimization and Pearson Correlation Coefficient Method. Glob. Geol. 2015, 18, 41–48. [Google Scholar]
  47. Arif, C.; Setiawan, B.I.; Hasanah, N.A.I.; Mizoguchi, M. Estimating Greenhouse Gas Emissions from Irrigated Paddy Fields in Indonesia under Various Water Managements. In Proceedings of the 1st International Conference on Mechanical Electronic and Biosystem Engineering (MEBSE 2018), Bogor, Indonesia, 15–16 December 2018; IOP Publishing: Bristol, UK, 2019; Volume 557, p. 012034. [Google Scholar]
  48. Wang, X.; Yi, S.; Wu, Q.; Yang, K.; Ding, Y. The Role of Permafrost and Soil Water in Distribution of Alpine Grassland and Its NDVI Dynamics on the Qinghai-Tibetan Plateau. Glob. Planet. Chang. 2016, 147, 40–53. [Google Scholar] [CrossRef]
  49. Ge, G.; Shi, Z.; Yang, X.; Hao, Y.; Guo, H.; Kossi, F.; Xin, Z.; Wei, W.; Zhang, Z.; Zhang, X. Analysis of Precipitation Extremes in the Qinghai-Tibetan Plateau, China: Spatio-Temporal Characteristics and Topography Effects. Atmosphere 2017, 8, 127. [Google Scholar] [CrossRef]
  50. Li, X.R.; Jia, X.H.; Dong, G.R. Influence of Desertification on Vegetation Pattern Variations in the Cold Semi-Arid Grasslands of Qinghai-Tibet Plateau, North-West China. J. Arid Environ. 2006, 64, 505–552. [Google Scholar] [CrossRef]
  51. Chang, B.; He, K.-N.; Li, R.-J.; Sheng, Z.-P.; Wang, H. Linkage of Climatic Factors and Human Activities with Water Level Fluctuations in Qinghai Lake in the Northeastern Tibetan Plateau, China. Water 2017, 9, 552. [Google Scholar] [CrossRef]
  52. Wang, W.; Guo, B.; Zhang, Y.; Zhang, L.; Ji, M.; Xu, Y.; Zhang, X.; Zhang, Y. The Sensitivity of the SPEI to Potential Evapotranspiration and Precipitation at Multiple Timescales on the Huang-Huai-Hai Plain, China. Theor. Appl. Clim. 2021, 143, 87–99. [Google Scholar] [CrossRef]
  53. Llanes-Cárdenas, O.; Gaxiola-Hernández, A.; Estrella-Gastelum, R.D.; Norzagaray-Campos, M.; Troyo-Diéguez, E.; Pérez-González, E.; Ruiz-Guerrero, R.; de, J. Pellegrini Cervantes, M.; Pellegrini Cervantes, M. Variability and Factors of Influence of Extreme Wet and Dry Events in Northern Mexico. Atmosphere 2018, 9, 122. [Google Scholar] [CrossRef]
  54. Chen, S.; Zhang, L.; Liu, X.; Guo, M.; She, D. The Use of SPEI and TVDI to Assess Temporal-Spatial Variations in Drought Conditions in the Middle and Lower Reaches of the Yangtze River Basin, China. Adv. Meteorol. 2018, 2018, 9362041. [Google Scholar] [CrossRef]
  55. Sun, T.; Yang, Y.; Wang, Z.; Yong, Z.; Xiong, J.; Ma, G.; Li, J.; Liu, A. Spatiotemporal Variation of Ecological Environment Quality and Extreme Climate Drivers on the Qinghai-Tibetan Plateau. J. Mt. Sci. 2023, 20, 2282–2297. [Google Scholar] [CrossRef]
  56. Guo, B.; Zhou, Y.; Zhu, J.; Liu, W.; Wang, F.; Wang, L.; Yan, F.; Wang, F.; Yang, G.; Luo, W.; et al. Spatial Patterns of Ecosystem Vulnerability Changes during 2001–2011 in the Three-River Source Region of the Qinghai-Tibetan Plateau, China. J. Arid Land 2016, 8, 23–35. [Google Scholar] [CrossRef]
Figure 1. The Qinghai–Tibet Plateau covers an area of about 2.5 × 106 km2. With an average altitude of about 4500 m, it is the source of many of China’s great rivers.
Figure 1. The Qinghai–Tibet Plateau covers an area of about 2.5 × 106 km2. With an average altitude of about 4500 m, it is the source of many of China’s great rivers.
Remotesensing 16 03461 g001
Figure 2. Logical flow. First, the Qinghai–Tibet Plateau is selected as the study area, and the SPEI was calculated by obtaining the precipitation, temperature, and evaporation data. Second, the data required for constructing the RSEI were obtained from MODIS satellite data. Then, PCA analysis was carried out at GEE using NDVI, LST, WET, and NDBSI to obtain RSEI, and the correlation between RSEI and SPEI was also analyzed. Finally, the future trend of RSEI was analyzed. Sub-figure (af) the annual mean Standardized Precipitation Evapotranspiration Index (SPEI) values.
Figure 2. Logical flow. First, the Qinghai–Tibet Plateau is selected as the study area, and the SPEI was calculated by obtaining the precipitation, temperature, and evaporation data. Second, the data required for constructing the RSEI were obtained from MODIS satellite data. Then, PCA analysis was carried out at GEE using NDVI, LST, WET, and NDBSI to obtain RSEI, and the correlation between RSEI and SPEI was also analyzed. Finally, the future trend of RSEI was analyzed. Sub-figure (af) the annual mean Standardized Precipitation Evapotranspiration Index (SPEI) values.
Remotesensing 16 03461 g002
Figure 3. (af) denote the spatio-temporal characteristics of SPEI on the QTP from 2000 to 2020 at different time scales. Note: The gray area in the figure is the Hengduan Mountains; SPEI data are missing.
Figure 3. (af) denote the spatio-temporal characteristics of SPEI on the QTP from 2000 to 2020 at different time scales. Note: The gray area in the figure is the Hengduan Mountains; SPEI data are missing.
Remotesensing 16 03461 g003
Figure 4. (at) Spatial and temporal characteristics of NDVI, WET, LST, and NDBSI on the QTP, 2000–2020. Where (ae) denote the values of NDVI on the QTP for 2000–2020, (fj) denote the values of LST on the QTP for 2000–2020, (ko) denote the values of WET on the QTP for 2000–2020, and (pt) denote the values of NDBSI on the QTP for 2000–2020.
Figure 4. (at) Spatial and temporal characteristics of NDVI, WET, LST, and NDBSI on the QTP, 2000–2020. Where (ae) denote the values of NDVI on the QTP for 2000–2020, (fj) denote the values of LST on the QTP for 2000–2020, (ko) denote the values of WET on the QTP for 2000–2020, and (pt) denote the values of NDBSI on the QTP for 2000–2020.
Remotesensing 16 03461 g004
Figure 5. Spatial and temporal characteristics of RSEI on the Qinghai–Tibet Plateau in 2000, 2005, 2010, 2015, and 2020. (ae) show the temporal and spatial distribution characteristics of RSEI in QTP. (fj) show the total area (km2) of best, good, normal, bad, and worst levels in each studying period.
Figure 5. Spatial and temporal characteristics of RSEI on the Qinghai–Tibet Plateau in 2000, 2005, 2010, 2015, and 2020. (ae) show the temporal and spatial distribution characteristics of RSEI in QTP. (fj) show the total area (km2) of best, good, normal, bad, and worst levels in each studying period.
Remotesensing 16 03461 g005
Figure 6. (a,b) Trends of RSEI and their significance in the Tibetan Plateau, 2000–2020. (a) represents the spatial distribution of RSEI trend values. Figure 6b represents the Z-value of the trend test. (b) can be classified as significant decline, nonsignificant decline, nonsignificant increase, and significant increase depending on the value of Z.
Figure 6. (a,b) Trends of RSEI and their significance in the Tibetan Plateau, 2000–2020. (a) represents the spatial distribution of RSEI trend values. Figure 6b represents the Z-value of the trend test. (b) can be classified as significant decline, nonsignificant decline, nonsignificant increase, and significant increase depending on the value of Z.
Remotesensing 16 03461 g006
Figure 7. Spatial correlation and significance test of RSEI and SPEI on the Tibetan Plateau. N-S-P-C: Not significantly positively correlated. S-P-C: Significantly positive correlation. E-S-P-C: Extremely significant positive correlation. E-S-N-C: Extremely significant negative correlation. S-N-C: Significantly negative correlation. N-S-N-C: Not significantly negatively correlated. (af) show the correlation and significance distribution characteristics of RSEI and SPEI in QTP. (gl) show the total area (km2) of N-S-N-C, S-N-C, E-S-N-C, E-S-P-C, S-P-C, and N-S-P-C in each studying period.
Figure 7. Spatial correlation and significance test of RSEI and SPEI on the Tibetan Plateau. N-S-P-C: Not significantly positively correlated. S-P-C: Significantly positive correlation. E-S-P-C: Extremely significant positive correlation. E-S-N-C: Extremely significant negative correlation. S-N-C: Significantly negative correlation. N-S-N-C: Not significantly negatively correlated. (af) show the correlation and significance distribution characteristics of RSEI and SPEI in QTP. (gl) show the total area (km2) of N-S-N-C, S-N-C, E-S-N-C, E-S-P-C, S-P-C, and N-S-P-C in each studying period.
Remotesensing 16 03461 g007
Figure 8. Hurst Index distribution of RSEI on the Tibetan Plateau and its future persistence characteristics.
Figure 8. Hurst Index distribution of RSEI on the Tibetan Plateau and its future persistence characteristics.
Remotesensing 16 03461 g008
Table 1. Indices for computing RSEI.
Table 1. Indices for computing RSEI.
IndexSatelliteResolutionLevel
NDVIMOD13A11000 m16 dL3
LSTMOD11A21000 m8 dL3
WETMOD09A1500 m8 dL3
NDBSIMOD09A1500 m8 dL3
Table 2. Aridity classification based on Standardized Precipitation Evapotranspiration Index (SPEI).
Table 2. Aridity classification based on Standardized Precipitation Evapotranspiration Index (SPEI).
LevelAridity DegreeSPEI Value
1No aridity−0.5 < SPEI
2Mild aridity−1.0 < SPEI ≤ −0.5
3Moderate aridity−1.5 < SPEI ≤ −1.0
4Severe aridity−2.0 < SPEI ≤ −1.5
5Extreme ariditySPEI ≤ −2.0
Table 3. Statistics on the trend of remote sensing ecological indices within the QTP region.
Table 3. Statistics on the trend of remote sensing ecological indices within the QTP region.
Z ValueTrend of RSEIPercentage
−2.20 < Z ≤ −1.96Significant decline0.3%
−1.96 < Z ≤ 0Nonsignificant decline61.45%
0 < Z < 1.96Nonsignificant increase36.62%
1.96 ≤ Z < 2.20Significant increase1.55%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, Y.; Cao, J.; Gu, Y.; Huang, X.; Liu, X.; Hu, Y.; Wu, S. An Evaluation of Ecosystem Quality and Its Response to Aridity on the Qinghai–Tibet Plateau. Remote Sens. 2024, 16, 3461. https://doi.org/10.3390/rs16183461

AMA Style

Yan Y, Cao J, Gu Y, Huang X, Liu X, Hu Y, Wu S. An Evaluation of Ecosystem Quality and Its Response to Aridity on the Qinghai–Tibet Plateau. Remote Sensing. 2024; 16(18):3461. https://doi.org/10.3390/rs16183461

Chicago/Turabian Style

Yan, Yimeng, Jiaxi Cao, Yufan Gu, Xuening Huang, Xiaoxian Liu, Yue Hu, and Shuhong Wu. 2024. "An Evaluation of Ecosystem Quality and Its Response to Aridity on the Qinghai–Tibet Plateau" Remote Sensing 16, no. 18: 3461. https://doi.org/10.3390/rs16183461

APA Style

Yan, Y., Cao, J., Gu, Y., Huang, X., Liu, X., Hu, Y., & Wu, S. (2024). An Evaluation of Ecosystem Quality and Its Response to Aridity on the Qinghai–Tibet Plateau. Remote Sensing, 16(18), 3461. https://doi.org/10.3390/rs16183461

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