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 km
2, 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).
where
PET is the potential evapotranspiration;
is the mean monthly air temperature;
is the Annual Heat Index; and
A is a function of
,
= 0.49 + 0.179
− 0.0000771
+ 0.000000675
. Potential evapotranspiration in this study uses the data sets that have been obtained (
Table 1).
is the difference between precipitation and potential evapotranspiration. The normalization of the sequence
can be calculated by Equations (4)–(7):
where
is the probability distribution function (log-logistic); G is the factorial function; and
,
, and
are the probability-weighted moments of
, computed as Equations (8) and (9):
where
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:
when the cumulative probability
≤ 0.5:
where the
SPEI is the Standard Precipitation Evapotranspiration Index.
,
,
,
,
and
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):
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
and
are the spectral reflectance in the infrared and near-infrared bands, respectively.
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:
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:
where
(
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).
where
,
,
,
and
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):
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 . The test statistic follows a normal distribution, indicated by “”. The “” term refers to the variance calculation, which is an essential part of the test’s methodology. Lastly, “” 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:
At a given confidence level α, the original hypothesis is invalid if || ≥ 1 − α/2. If > 0, the extent of fluctuation within the time series expands; if < 0, the range of variation of the time series decreases, then it indicates a decrease in the range of variation. In general, if || ≥ 2.58, the trend of change is deemed to be of extreme significance. If || ≤ 2.58, the trend is deemed significantly pronounced. If 1.96 ≤ || < 2.58, the trend is regarded as statistically significant. If 0 ≤ || < 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:
In this context,
r signifies the correlation coefficient, while
and
denote the two variables under consideration, and
and
represent the average values of the respective variables.
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:
For the ratio , a Hurst phenomenon exists in the time series if the following relationship exists . H is the Hurst exponent, which can be obtained by using least-squares fitting according to . 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.
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.