Evaluation of the Ecological Environment Quality of the Kuye River Source Basin Using the Remote Sensing Ecological Index

Landsat remote sensing images obtained from 2000, 2005, 2010, 2015, and 2020 were analyzed. The normalized vegetation index (NDVI), moisture index (WET), land surface temperature (LST), and normalized building-soil index (NDBSI) were extracted based on the four aspects of greenness, humidity, heat, and dryness. The Remote Sensing Ecological Index (RSEI) was calculated using principal component analysis to quantitatively analyze and dynamically monitor and evaluate the ecological environment changes in the Kuye River Basin over the past 20 years. From the perspective of spatial and temporal distribution, the ecological and environmental quality of Kuye River Basin had a downward trend from 2000 to 2020. The overall RSEI grade was medium or poor, and the average RSEI decreased. The proportion of excellent and good grade watershed areas decreased, whereas that of medium, low, and poor grade watershed areas increased over the study period. Spatially, RSEI decreased gradually from southeast to northwest. The degraded areas were mainly distributed in urban areas with frequent human activities. Conversely, the superior eco-environmental quality areas were mainly distributed in eastern sections of the watershed. Compared with 2000, the eco-environmental quality of the Yulin urban area and Shenmu County in the southern section of the watershed are worsening.


Introduction
The level of ecological environment quality is related to human production [1], economic and social activities, and interactions of the three factors which can restrict and affect the organic whole [2]. The ecological environment quality not only affects the quality of life of the human inhabitants, but also directly affects the stable sustainable development of the economy [3]. Ecological environmental quality assessment monitors the ecological environmental quality by investigating and analyzing ecological environmental quality of actual situations [4]. The process utilizes scientific experimental methods, data, evaluation indicators, and an evaluation system to analyze and assess changes in ecological environmental quality [5]. With the development of science and technology, human material civilizations have reached unprecedented heights but remain increasingly associated with serious ecological environment issues [6]. Humans excessively extract natural resources in the blind pursuit of their economic interests [7]. The contradiction between economic development and the allocation of resources is increasing and poses a serious threat to the survival and development of human beings [8]. Globally, people are paying increasing attention to the ecological environment. The concept of "sustainable development" was introduced at the United Nations Conference on Development and Environment in Brazil in 1992 [9]. It provided clear recognition of the need to protect the environment while pursuing economic development. Protecting the environment and achieving sustainable and 2020 to extract the NDVI, WET, LST, and NDBSI based on the four aspects of greenness, humidity, heat, and dryness, respectively. RSEI provides the ecological environmental changes in the Kuye River Basin over the past 20 years using quantitative analysis with dynamic monitoring. Our research objective is to use remote sensing and GIS methods to extract ecological environment assessment indicators, to track and monitor the ecological environment, which can effectively grasp the ecological environment quality of the watershed, improve the quality of decision-making, and quickly prevent the excessive deterioration of the ecological environment.

Study Area
The basin is located 38°28'-39°52' north and 109°00'-110°52' east. Kuye River is the second largest tributary of the middle reaches of the Yellow River. It originates from Badingou, Dongsheng City, Inner Mongolia Autonomous Region. The Kyye River has a total length of 242 km and a total basin area of 8 706 km 2 . It flows through Ejin Horo Qi and Fugu, Shaanxi Province and into the Yellow River. The topography of the Kuye River Basin is complex and diverse. The terrain is high in the northwest and low in the southeast. The surface vegetation is scarce, the forest and grass coverage is low, the underlying surface conditions are complex, and there are desert steppe, meadow, artificial irrigation, and agricultural vegetation. The climate type of the basin is semi-arid continental monsoon in the north temperate zone. The annual average precipitation is 419 mm, which varies greatly over a year. The annual maximum precipitation mainly occurs in July and August, which accounts for 55-65% of the total annual precipitation. The rainfall is spatially distributed in the southeast and the northwest. The annual average temperature is 8.2 °C, with the highest temperatures in July ( Figure 1).

Date Source and Date Processing
In 2021, the United States geological survey (USGS) used the GEE platform (https:// developers.google.com/earth-engine/datasets/catalog/landsat, accessed on 27 September 2022) to create a new Landsat Collection Level 2 dataset containing positively-corrected Landsat surface reflectance and surface temperature bands [26]. In our study area, the bands better reflect the ecological environment conditions in the vegetation growing season. We selected images from the dataset with cloud cover less than 20% in the years of 2000, 2005, 2010, 2015, and 2020. We used the cloud removal algorithm provided in GEE platform to remove the cloud and its shadow using QA_PIXEL of the pixel quality band and synthesized the multi-layer image into a multi-band image with good quality using the median value. To ensure sufficient image quantity and quality, we selected only images with the cloud cover below 20%. Due to the large time scale, the data sets used different satellites/sensors in different years. Landsat 7 ETM images provided the data for 2000 and 2005, Landsat 5 TM images were used for the 2010 data, and Landsat 8 Oli/TIRS images were used for the 2015 and 2020 data.
The runoff and sediment transport data from the Kuye River Basin from 1956-2009 were obtained from the Yellow River Hydrological Yearbook, while the same data for 2010-2020 were obtained from the Yellow River sediment Bulletin. The precipitation data for 1956-2020 were obtained from the surface daily meteorological observation data of four meteorological stations (Dongsheng, Hequ, Yulin, and Xingxian) from the China Meteorological Data Network (http://data.cma.cn, accessed on 27 September 2022). The spatial interpolation used the mean values for the basin to analyze the climate change trends.

RSEI Calculation
The functions of principal component transform of remote sensing image are image compression, image denoising, image enhancement, image fusion, and feature extraction. Principal component analysis (PCA) is used to reduce the dimension of the index bands of greenness, humidity, heat, and dryness to obtain the irrelevant linear component data, and then extract the components with large contribution rate and discard the redundant component data with little information and analyze the component data with large weight. In our study, the dimension of the remote sensing data is degraded by compressing the component data of the four indicators, removing redundant data, enhancing the remote sensing image, and obtaining the eigenvalues of the component indicators after principal component analysis.
The RSEI contains four indices-greenness, wetness, dryness, and heat. These are represented by the NDVI, WET, NDBSI, and LST, respectively [27]. They are obtained by remote sensing and then coupled in a PCA. The greenness index (NDVI) was calculated using the NDVI. NDVI is an important biophysical parameter reflecting the state of surface vegetation, which is closely related to plant biomass, vegetation coverage, and leaf area index. In the field of ecological environment, due to the expansion of ecosystem research results at the landscape patch level, NDVI has become the connection point of spatial scale expansion. Vegetation Index can be used to reflect the change in vegetation cover and land use, to realize the study of ecological environment quality. Vegetation connects soil, atmosphere, and water in the terrestrial ecosystem, regulates the material and energy exchange between the interior and exterior of the ecosystem, and has an important influence on the global energy cycle and the biochemical ring of substances. The vegetation has obvious interannual and seasonal changes. When the vegetation cover changes in the region, it will have a direct impact on the ecological environment quality of the whole region. The wetness index (WET) was calculated using a tasseled cap transformation of the remote sensing images. The dryness (NDBSI) was obtained using the bare land index and the built-up area index as average values. The heat index (LST) was obtained by calculating the actual temperature of the ground surface using thermal infrared. The RSEI is defined as follows: where NIR and R are the reflectance of the near -infrared and red bands, respectively. 2 Wetness The tasseled cap transformation algorithm can compress and remove redundant data. Its greenness, humidity, and brightness components are directly related to the physical parameters of the surface [25], and widely used in ecological environment monitoring. The moisture component of the Tasseled Cap Transformation Algorithm is closely related to the soil moisture and vegetation. In this study, the moisture component of the silk hat transformation provides the humidity index to reflect the surface humidity in the basin. When calculating the humidity component of the silk hat transformation, the coefficients of each band of the different Landsat sensor images were different, which are expressed as follows: where Wet (TM), wet (ETM+) and wet (Oli) are the moisture components of the Tasseled Cap Transformation Algorithm of the Landsat TM, ETM+, and Oli images, respectively. ρ blue , ρ green , ρ red , ρ NIR , ρ SWIR1 , and ρ SWIR2 are the surface reflectance of Landsat TM, ETM+, and Oli images in blue, green, red, near infrared, short wave infrared 1, and short-wave infrared 2 bands, respectively. 3 Dryness Soil drying causes serious pollution in the regional ecological environment. The dryer the soil, the more harm to the environment. In addition to bare soil (such as rock, sandy land, bare ground), urban building areas also cause soil "drying" [28]. Therefore, the dryness index (NBDSI) was developed using the bare soil index (IS) and building index (IBI), which can be calculated using the following formula [29]: IBI = 2ρ swirl /ρ swirl + ρ nir − ρ nir /(ρ nir + ρ red ) + ρ green /ρ green + ρ swirl 2ρ swirl /ρ swirl + ρ nir + ρ nir /(ρ nir + ρ red ) + ρ green /ρ green + ρ swirl (7) where Landsat TM, ETM+, and Oli correspond to the reflectance in each band. SI is the bare earth index, and IBI is the building index. 4 Heat (LST) The atmospheric correction method was used for the inversion of LST [30]. After selecting the thermal infrared band of the image for the radiometric calibration, the blackbody radiation brightness was obtained [31], and the LST was calculated using the Planck function [11], with the following formula [32]: where L λ is the brightness value of the thermal infrared radiation; B (LST) is blackbody radiation brightness; LST is the real surface temperature; Epsilon for Surface emissivity is τ; L↓ and L↑ are the atmospheric profile parameters.
Since the four indicators have their own dimensions and each indicator has an unbalanced weight, there was no way to directly collate them into the principal component analysis. Therefore, the data was normalized into [0, 1] prior to performing the PCA. The following formula for each indicator was as follows [33]: where NI i is the result of the image standardization; I i is the numerical pixel value of the image. I max and I min are the maximum and minimum values of the image pixel, respectively; PCA is the principal component analysis, and f is the set of functions used to perform the PCA; The maximum and minimum values of RSEI 0 are expressed as RSEI 0max and RSEI 0min , respectively.

Spatial Autocorrelation Analysis
Spatial autocorrelation is an important index to investigate whether the attribute value of an element significantly correlates with the attribute value of its adjacent space. This reveals the correlation of the attribute eigenvalues between spatial reference units and their adjacent spatial units. We used both Moran's index (global spatial autocorrelation) and the local spatial correlation index (Lisa) to analyze the spatial correlation of the ecoenvironmental quality. The global spatial autocorrelation analysis investigated the overall spatial distribution characteristics of the regional RSEI which was represented by Global Moran's I using the following formulas [34]: (15) where x i is the attribute value of position i; N is the total number of grids in the study area; W ij is the weight of the matrix, which represents the relationship of the spatial objects, i = 1, 2, 3..., N, j = 1, 2, 3..., M, W ij = 0 when i and j are adjacent. The Moran's I ranges from approximately +1 (positive spatial autocorrelation) to −1 (negative autocorrelation), with zero indicating no spatial autocorrelation. LISA can reflect the local correlation of the spatial distribution of an object's attributes. Even when I in the global Moran's I is 0, there may be local clustering. Therefore, it is necessary to combine the LISA and Global Moran's I analyses. The formula was expressed as follows: where LISA represents the spatial clustering of similar values (high value or low value) around the spatial unit, and a negative LISA represents the spatial clustering between dissimilar values. The results of the local spatial autocorrelation analysis generated a LISA cluster, which presented five spatial distribution types: high-high (H-H) cluster, low-low (L-L) cluster, high-low (H-L) outliers, low-high (L-H) outliers and insignificant values. H-H was the presence of a high value surrounded by a high value, and L-L was the presence of a low value surrounded by a low value. H-L indicates a high-value anomaly and L-H indicates a low-value anomaly. Non-significant values indicate the attribute values are close to randomly distributed.

RSEI Flow Chart
RSEI uses a principal component transformation by aggregating the main information provided by the NDVI, WET, NDBSI, and LST indicators into the first principal component so that they can be described using a single indicator [35]. This method integrates the information closely related to the four indicators and has advantages in constructing a RSEI. The AM, ETM, and image data stored in the GEE cloud platform were used in the spatiotemporal distribution maps of the RSEI in Kuye River Basin in 2001, 2005, 2010, 2015, and 2020 and generated using the PCA. The Moran's index (global spatial auto-correlation) and LISA (local indicator of spatial association) analyzed the spatial correlations for the eco-environmental quality. The driving factors of the eco-environmental quality changes were analyzed [17]. The lower the RSEI value, the worse the ecological conditions. The RSEI was further normalized between 0 and 1, in which 0 indicated an extremely poor status and 1 indicated an excellent status. Then, following previous studies, the normalized RSEI was reclassified into five levels (excellent: 0.8-1.0; good: 0.6-0.8; moderate: 0.4-0.6; fair: 0.2-0.4; poor: 0-0.2) using an equal interval of 0.2. A detailed workflow of this study is provided in Figure 2.

RSEI Flow Chart
RSEI uses a principal component transformation by aggregating the main information provided by the NDVI, WET, NDBSI, and LST indicators into the first principal component so that they can be described using a single indicator [35]. This method integrates the information closely related to the four indicators and has advantages in constructing a RSEI. The AM, ETM, and image data stored in the GEE cloud platform were used in the spatiotemporal distribution maps of the RSEI in Kuye River Basin in 2001, 2005, 2010, 2015, and 2020 and generated using the PCA. The Moran's index (global spatial auto-correlation) and LISA (local indicator of spatial association) analyzed the spatial correlations for the eco-environmental quality. The driving factors of the eco-environmental quality changes were analyzed [17]. The lower the RSEI value, the worse the ecological conditions. The RSEI was further normalized between 0 and 1, in which 0 indicated an extremely poor status and 1 indicated an excellent status. Then, following previous studies, the normalized RSEI was reclassified into five levels (excellent: 0.8-1.0; good: 0.6-0.8; moderate: 0.4-0.6; fair: 0.2-0.4; poor: 0-0.2) using an equal interval of 0.2. A detailed workflow of this study is provided in Figure 2.   Table 1 states the contribution rate of the eigenvalue of the first principal component of the RSEI in the fifth period is above 85%, and the contribution rate of each index has the same positive and negative distribution in the first principal component. Specifically, the contribution rate of NDVI (representing vegetation coverage) and WET (representing environmental humidity) in the first principal component is positive. LST (representing land surface temperature) and NDBSI (representing hardening degree of buildings and bare soil) were negative values. This result indicates NDVI and WET played a promoting role in improving the ecological environment. Conversely, LST and NDBSI played a restraining role in improving the ecological environment. This result was consistent with the reality of the effects of the four indicators on the ecological environment.  In 2000, 2005, 2010, 2015, and 2020, the contribution of the dryness index in PC1 was −0.851, −0.853, −0.861, −0.879, and −0.891, respectively. This demonstrates a gradually increasing trend. The cause is the urbanization of Kuye River Basin which accelerated after 2000, increasing the area of build-up land. The dryness index reflects the change in the construction land in the study area. This phenomenon demonstrates that the urban planning and construction of the Kuye River Basin has obvious influences on the urban eco-environmental quality (especially the increased construction land area caused by urban expansion) which has become the key factor affecting the eco-environmental quality of the Kuye River Basin.

NDVI
The higher the greenness value, the higher the vegetation coverage rate and the better the ecological environment quality [1]. Conversely, the lower the greenness value, the lower the vegetation coverage and the worse the ecological environment quality [17,27]. The greenness index in 2000, 2005, 2010, 2015, and 2020 after the standardization treatment are between 0 and 1. As demonstrated in Figure 3, the overall distribution characteristics of the NDVI in the five different years are basically the same, with the characteristics of low in the northwest and high in the southeast, and the NDVI changes significantly in the southeast. In 2000, the whole Kuye River Basin was dominated by a high value area (blue), the NDVI high value area was mainly distributed in the south, and the NDVI low value area was mainly distributed in the north. In 2005, when compared with 2000, the low value range increased significantly (light brown), indicating that the vegetation coverage of the region was decreasing. In 2010, when compared with 2005, the NDVI value in northern China decreased slightly. In 2015, when compared with 2010, the NDVI value in southern China decreased slightly. In 2020, the NDVI value decreased significantly in the south when compared to 2015.

WET
The higher the humidity value, the higher the surface water content, the better the ecological environment quality, and vice versa [31]. The humidity index in 2000, 2005, 2010, 2015, and 2020 after the standardization treatment is between 0 and 1. As demonstrated in Table 2, from 2000-2020, the humidity in Kuye River Basin initially decreased before increasing, but the overall trend was decreasing. This indicates the increase in vegetation coverage increased the demand for water, resulting in a decrease in the ground humidity.

NDBSI
The higher the NDBSI value, the more serious the soil degradation in the region and the worse the ecological environment quality. Conversely, a low NDBSI value relates to a high ecological quality [36]. The values of the dryness index in 2000, 2005, 2010, 2015, and 2020 after the standardization treatment is between 0 and 1. As demonstrated in Table 2, the mean values of the NDBSI index of the five stages were 0.282, 0.321, 0.332, 0.423, and 0.528, respectively. Therefore, as land construction in the Kuye River Basin increased the soil deteriorated sharply from 2000 to 2020 as economic development and urbanization accelerated. In addition, high intensity groundwater extraction led to decreases in surface water in the Kuye River Basin, an increase in the dryness index in some areas, and the deterioration of the ecological environment quality.

LST
In Figure 3, green and red represent areas with high and low heat values, respectively. The higher the heat value, the worse the ecological quality within the region. Table 2 Figure 3 graphically demonstrates the five-stage LST distribution map, the high temperature area (in red) encompasses the main urban area and areas of exposed soil, while the non-urban areas and the river possess an appropriate temperature (in green). Figure 4 provides the spatial and temporal distribution map of the ecological environmental quality in Kuye River Basin from 2000-2020. Red, white, yellow, light blue, and dark blue represent excellent, good, medium, poor, and bad ecological grades, respectively. The area ( Figure 5) and proportion (Table 3) of each ecological level were calculated using the five RSEI maps from 2000, 2005, 2010, 2015, and 2020. As demonstrated in Figure 5, the area of moderate and poor RSEI increased over time. Conversely, the area of good and excellent RSEI decreased, indicating that the overall ecological environment quality of the Kuye River Basin was declining. Table 3  In Figure 3, green and red represent areas with high and low heat values, respectively. The higher the heat value, the worse the ecological quality within the region. Table  2 shows the regions with high LST values in Kuye River Basin exhibited an upward trend from 2000 to 2020. The LST mean values in the five periods were 19.521, 20.139, 20.333, 20.536, and 20.602, respectively. Figure 3 graphically demonstrates the five-stage LST distribution map, the high temperature area (in red) encompasses the main urban area and areas of exposed soil, while the non-urban areas and the river possess an appropriate temperature (in green). Figure 4 provides the spatial and temporal distribution map of the ecological environmental quality in Kuye River Basin from 2000-2020. Red, white, yellow, light blue, and dark blue represent excellent, good, medium, poor, and bad ecological grades, respectively. The area ( Figure 5) and proportion (Table 3) of each ecological level were calculated using the five RSEI maps from 2000, 2005, 2010, 2015, and 2020. As demonstrated in Figure  5, the area of moderate and poor RSEI increased over time. Conversely, the area of good and excellent RSEI decreased, indicating that the overall ecological environment quality of the Kuye River Basin was declining. Table 3      From 2000-2005, the area with high ecological and environmental quality in Kuye River Basin decreased significantly, from 66.98% (GE) to 27.49%. Firstly, the local climate conditions caused large areas of vegetation to die, resulting in decreased vegetation coverage. Second, the rapid economic development, population density, and per capita GDP growth increased the human activities, damaging the surface ecological environment. In 2005-2010, the area with high river basin ecological environment quality increased, this was due to the development of soil and water conservation engineering. The active promotion of the ecological construction engineering by river basin ecological environment construction, organization of local people in controlling desertification, development of agriculture in arid areas, solving the problems of food and clothing for the farmers, and adjusting the rural industrial structure. Ecological awareness drives the adjustment and optimization of land use structure, develops the river basin economy, and attempts to promote improvements in the ecological environment quality. From 2015-2020, the areas with high ecological and environmental quality in the Kuye River Basin decreased significantly. During this period, the rapid urbanization of the Kuye River Basin led to the deterioration of the ecological environment quality. The environment quality was especially impacted by the rapid urban population growth which involved the reclamation of new land, destroying the original soil structure. Urbanization affects the soil quality such as  From 2000-2005, the area with high ecological and environmental quality in Kuye River Basin decreased significantly, from 66.98% (GE) to 27.49%. Firstly, the local climate conditions caused large areas of vegetation to die, resulting in decreased vegetation coverage. Second, the rapid economic development, population density, and per capita GDP growth increased the human activities, damaging the surface ecological environment. In 2005-2010, the area with high river basin ecological environment quality increased, this was due to the development of soil and water conservation engineering. The active promotion of the ecological construction engineering by river basin ecological environment construction, organization of local people in controlling desertification, development of agriculture in arid areas, solving the problems of food and clothing for the farmers, and adjusting the rural industrial structure. Ecological awareness drives the adjustment and optimization of land use structure, develops the river basin economy, and attempts to promote improvements in the ecological environment quality. From 2015-2020, the areas with high ecological and environmental quality in the Kuye River Basin decreased significantly. During this period, the rapid urbanization of the Kuye River Basin led to the deterioration of the ecological environment quality. The environment quality was especially impacted by the rapid urban population growth which involved the reclamation of new land, destroying the original soil structure. Urbanization affects the soil quality such as surface runoff, soil water infiltration, and soil water content, resulting in the reduction of surface vegetation coverage and potentially causing excessive soil and water loss. There are many coal mining enterprises in Kuye River Basin. Coal mining requires a lot of resources and energy, and discharges pollution. Coal mining lowers the underground water level, which reduces or ceases river water flow. The natural ecosystem of the Kuye River Basin is fragile, the cultivated land in the basin has sharply reduced. Land desertification is becoming increasingly serious, the forest resources have sharply reduced and are now limited. The freshwater resources are seriously limited, and the biodiversity is constantly decreasing. Therefore, the natural environment in sections of the basin have been seriously damaged, and the ecological environment quality has been significantly decreased. Figure 4 demonstrates that in 2000 the RSEI level increased gradually from upstream to downstream and reduced gradually from east to west. The ecological environment quality is better downstream and is superior in the west. The precipitation diminishes from southeast to northwest. The temperature increased from east to west. The vegetation coverage was better in the lower reaches than in the upper reaches. From 2005 to 2015, excellent eco-environmental quality areas were mainly distributed in the eastern section of the watershed. Conversely, the eco-environmental quality in the Yulin urban area and Shenmu County in the southern section of the watershed showed a worsening trend when compared to the results in 2000. This change is mainly because coal mining in these regions caused serious environmental pollution and damage to the watershed. In 2020, the area of low ecological and environmental quality reached 91.8% (BPM). The low-quality areas were mainly distributed in the northwestern region of Kuye River Basin, which is the desert area of the southern edge of the Mu Us desert and covered by yellow sand. In the west, sandstorms cause serious erosion, and the area belongs to an ecologically fragile district. The RSEI level of this area is poor. The regional natural ecological environment conditions are relatively unsatisfactory with an unreasonable ecological structure, ecological environment destruction, a weak ecological function and resilience system, and a very difficult area for undertaking ecological restoration.  19.891, and in 2020, they were 0.255, 7.75, respectively, with higher Z values than the critical level of 1.96. At the significance level of 5%, the analysis revealed the RSEI was not a completely random state, with a significant correlation in space. RSEI positively correlated with spatial information. The remote sensing ecological index increases in areas with large spatial distributions of ecological environmental quality. Conversely, the remote sensing ecological index decreases if the spatial distribution of ecological environmental quality is small. From 2000 to 2020, Moran's index had a downward trend, the degree of spatial autocorrelation became weaker, and the spatial distribution dispersed. This result was consistent with the spatial and temporal variation characteristics of the ecological environment quality level results.

Local Indicators of Spatial Autocorrelation of RSEI
To understand the spatial and temporal distribution of ecological environmental quality, we analyzed the local spatial correlation patterns of RSEI using LISA clustering. The LISA clustering diagram (Figure 7) revealed H-H clustering areas were mainly distributed in the east of the watershed. L-L was mainly distributed in the northwest of the watershed, as well as along the river and coal mining development zones. These areas are densely populated, especially the L-L cluster area which increased significantly from 2015 to 2020, due to the rapid urban growth and the increase in industrial and mining land. This development destroys the original soil structure, reduces the surface vegetation coverage, and potentially increases soil and water loss. Our results support previous studies. The area of built-up land in Kuye River Basin increased from 78.29 hm 2 in 1980 to 877.25 hm 2 in 2020. The built-up area has the highest growth rate among all the land use types, indicating that human activities rapidly expanded causing built-up land and generated a higher demand on the environment.

Local Indicators of Spatial Autocorrelation of RSEI
To understand the spatial and temporal distribution of ecological environmental quality, we analyzed the local spatial correlation patterns of RSEI using LISA clustering. to 2020, due to the rapid urban growth and the increase in industrial and mining land. This development destroys the original soil structure, reduces the surface vegetation coverage, and potentially increases soil and water loss. Our results support previous studies. The area of built-up land in Kuye River Basin increased from 78.29 hm 2 in 1980 to 877.25 hm 2 in 2020. The built-up area has the highest growth rate among all the land use types, indicating that human activities rapidly expanded causing built-up land and generated a higher demand on the environment.

The Influence of Natural Factors on the Ecological Environment
Good research results can indeed be obtained by using a variety of remote sensing data for ecological environmental quality assessment, but the time scale of these data is not long, so it is impossible to study the ecological environmental quality change of Kuye

The Influence of Natural Factors on the Ecological Environment
Good research results can indeed be obtained by using a variety of remote sensing data for ecological environmental quality assessment, but the time scale of these data is not long, so it is impossible to study the ecological environmental quality change of Kuye River Basin in 2000. Compared with other satellites, Landsat series satellite image data have large relative data storage, long time series and relatively moderate 30 m resolution. More importantly, GEE (Google Earth Engine) platform has level-2 remote sensing data synthesized by TOA of Landsat series released by USGS. Cloud interference and noise can be effectively removed, and seamless image stitching can be achieved to reduce the lack of image information.
Greenness, humidity, heat, and dryness are important components of an ecological environment. The ecological index is based on the main information from the four indexes and better reflects the ecological environment quality of Kuye River Basin [37]. Each index of the ecological environment interacts and influences each other. When the temperature rises, the activity of some microorganisms increases, and the activity of some microorganisms decreases which changes the types and quantities of microorganisms, resulting in changes in soil oxygen content and fertility, RSEI index, and finally ecological environment quality. When rainfall increases, there is sufficient water in the basin, river runoff increases, plants grow, soil moisture content is rich, and soil structure is stable. Meanwhile, increasing rainfall causes soil erosion to intensify. There are serious impacts of soil erosion such as, soil internal structure changes, and finally RSEI is affected. Figure 8 demonstrates the annual precipitation change rate in the Kuye River Basin was 0.375 mm/A (p > 0.1). The average annual precipitation had an increasing trend that was not statistically significant, and the average annual temperature increased at a rate of 0.026°C/a in Kuye River Basin (p > 0.1). The precipitation changes did not differ statistically. However, the average annual temperature showed a statistically insignificant but increasing trend. Water and heat are two of the most important natural conditions for promoting vegetation growth and vegetation cover, which affect the ecological processes of plant communities and the plant growth cycle. Precipitation changes are directly related to NDVI and WET in this region, they can improve the quality of the local ecological environment [14]. Therefore, water and heat conditions play an important role in the ecological environment of the Kuye River. However, in the study area, both temperature and rainfall showed no significant increasing trend, while NDVI and WET showed a decreasing trend, indicating other factors are affecting NDVI and WET. NDBSI (which is negatively correlated with the ecological environment) contributes the most to RSEI, which corresponds to the climate warming characteristics in the Kuye River Basin as global climate change progresses. The negative correlation between LST and RSEI corresponds with the obvious heat island effect in the basin. rainfall showed no significant increasing trend, while NDVI and WET showed a decreasing trend, indicating other factors are affecting NDVI and WET. NDBSI (which is negatively correlated with the ecological environment) contributes the most to RSEI, which corresponds to the climate warming characteristics in the Kuye River Basin as global climate change progresses. The negative correlation between LST and RSEI corresponds with the obvious heat island effect in the basin.

The Influence of Human Activities on Ecological Environment Change
Human activity refers to the various livelihood activities [38]. From 2000 to 2020, the population density in Kuye River Basin showed a straight upward trend, with the population increasing from 500,000 in 2000 to 1 million in 2020, and the per capita GDP had a

The Influence of Human Activities on Ecological Environment Change
Human activity refers to the various livelihood activities [38]. From 2000 to 2020, the population density in Kuye River Basin showed a straight upward trend, with the population increasing from 500,000 in 2000 to 1 million in 2020, and the per capita GDP had a "J-shaped" growth trend. However, the basin area is fixed, and the ecological carrying capacity is limited, so the population and economic growth have gradually approached the environmental limit. Population growth not only causes a higher consumption of resources, but also enlarges the ecological disturbance. The increased population leads to frequent human activities. To meet human demands, it is necessary to accelerate the exploitation of environmental resources, forcing the city to continuously expand outward, changing other land types into construction land, and ultimately changing the ratio of the land type structure, affecting the ecological environment quality.
Population growth and economic development requires human beings to consume resources for faster development, which leads to vegetation destruction, changes in land use type, and the decline of the ecological environment quality. The Kuye rivers region is mainly in the arid areas which were inundated for a long time. As the population increased in the basin, the river was used for irrigation and the groundwater was continuously exploited. As a result, the river surface runoff is significantly reduced (Figure 9), and a reduced surface runoff can lead to a reduced water vapor content, reducing the soil moisture. In the 1980s, a large coal reserve was discovered in the Kuye River Basin, accounting for 20% of China's coal reserves. Local people began to exploit the coal, and the volume of mined coal is continuously increasing. The average annual volume of mined coal increased from approximately 1 × 107 t in 1997 to 2.5 × 108 t in 2017. The underlying surface conditions have changed, affecting the discharge of surface runoff. In addition, most coal mining enterprises are located in mountainous areas, and coal mining seriously damages the surface vegetation, resulting in a significant decrease in vegetation coverage, a significant reduction in river runoff, and often interrupts water flow. The large-scale afforestation in the Kuye River Basin also negatively impacts the ecological environment. One of the main manifestations of this is the likelihood of increased drought events in dry years. Kuye river basin is attempting to return the farmland to forest which should dramatically change the watershed vegetation quantity, increasing the frequency and intensity of drought events. An increased vegetation coverage may result in a decline in the soil moisture over a long time scale. This will bring adverse effects on the returning farmland to forest References project, causing dwarf vegetation with stunted growth and a reduced resistance to drought, causing massive vegetation death. The internal surface of Kuye River Basin is unsustainable, with serious land desertification in a very fragile ecosystem. This causes the speed of improvement to have difficulty in restricting the speed of deterioration of the ecological environment quality. In addition, in rural areas with a high dependence on natural resources, the ecological environmental protection and economic and social welfare are complementary. Only by eliminating poverty can we effectively improve the environmental quality of the caves and wild rivers, to avoid achieving development at the expense of environment.

Conclusions
In this study, we used a PCA coupling NDVI, WET, NBDSI, and LST on the GEE platform to construct a macroscopic and efficient ecological and environmental quality monitoring method based on RSEI. We quantitatively analyzed the spatial and temporal distribution of the ecological and environmental quality in the Kuye River Basin. We obtained regional ecological environment quality information quickly providing an effective and simple means to reveal the dynamic changes of the ecological environment in the Kuye River Basin. Our analysis objectively and quantitatively investigated time and space and overcame the shortcomings of complex ecological evaluation indexes and cumbersome processes used in the past. The results show that the RSEI evaluation index can effectively reflect the ecological environment quality of the Kuye River Basin with high applicability. NDVI and WET had a positive effect on the ecological environment quality of the Kuye River Basin, while NDBSI and LST had a negative effect on the ecological environment quality. The inhibitory effect of NDBSI and LST on the ecological environment were significantly greater than the promoting effect of NDVI and WET on the ecological

Conclusions
In this study, we used a PCA coupling NDVI, WET, NBDSI, and LST on the GEE platform to construct a macroscopic and efficient ecological and environmental quality monitoring method based on RSEI. We quantitatively analyzed the spatial and temporal distribution of the ecological and environmental quality in the Kuye River Basin. We obtained regional ecological environment quality information quickly providing an effective and simple means to reveal the dynamic changes of the ecological environment in the Kuye River Basin. Our analysis objectively and quantitatively investigated time and space and overcame the shortcomings of complex ecological evaluation indexes and cumbersome processes used in the past. The results show that the RSEI evaluation index can effectively reflect the ecological environment quality of the Kuye River Basin with high applicability. NDVI and WET had a positive effect on the ecological environment quality of the Kuye River Basin, while NDBSI and LST had a negative effect on the ecological environment quality. The inhibitory effect of NDBSI and LST on the ecological environment were significantly greater than the promoting effect of NDVI and WET on the ecological environment.
The WET had a greater promoting effect on the ecological quality of the study area and had a better effect on the improvements in the ecological quality of the study area. The NDBSI index had a great inhibitory effect on the changes in ecological quality, indicating that urban expansion activities will lead to a significant reduction in the regional ecological quality.
Since 2000, the ecological environment quality of the Kuye River Basin has fluctuated and declined. Although a series of measures such as industrial structure adjustment, sand control, and grass planting have achieved initial results, reflecting the ecological environment has recovered to a certain degree. In some areas of the basin, however, the land desertification is serious, the vegetation is sparse, short and low coverage rate, and the ecosystem is fragile. The speed of improvement of ecological environment quality cannot catch up with the speed of deterioration of ecological environment quality. With the implementation of the basin development strategy, several energy and chemical bases have been built in the basin, the rapid expansion of cities and the rapid growth of population have rendered the basin's ecological environment diverse, volatile, and fragile.
From the perspective of spatial and temporal distribution, the ecological and environmental quality of Kuye River Basin had a downward trend from 2000 to 2020. The overall ecological quality level was mainly moderate or poor. The mean value of the RSEI decreased from 0.595 in 2000 to 0.387 in 2020. The sum proportion of excellent RSEI grade areas decreased from 66.98% in 2000 to 3.2% in 2020. The sum proportion of medium, low, and poor RSEI grade areas increased from 33.02% in 2000 to 91.8% in 2020. The spatial distribution of RSEI decreased gradually from southeast to northwest. The degraded areas were mainly distributed in urban areas with frequent human activities. The excellent RSEI grade areas were mainly distributed in the eastern section of the watershed. The ecological environmental quality of the Yulin urban area and Shenmu County in the southern section of the watershed had a worsening trend compared with 2000.
To improve the ecological environment of the Kuye River Basin in the future, we should continue to strengthen the supervision of coal enterprises in the Kuye River Basin. All coal enterprises should deal with pollution at the source, and this should be combined with a series of industrial upgrading measures to improve the coal industry impacts on the environment. We will continue to strengthen the implementation of ecological remediation projects, strengthen the publicity and protection of the ecological environment, actively promote the transformation and development of mining areas, and strive to build a resourcesaving and environmentally friendly development strategy. In the development planning of the Kuye River Basin, it is necessary to focus on green development by implementing scientifically supported projects returning farmland to forest or grassland. It is important to coordinate development with the ecological environment and land use if we are to achieve sustainable development in the basin.
Author Contributions: The research design, date collection, methodology and formal analysis were completed by Q.L. and X.M. The manuscript was written by Q.L. and F.Y. All authors have read and agreed to the published version of the manuscript.