Active–Passive Remote Sensing Evaluation of Ecological Environment Quality in Juye Mining Area, China

: The coal industry is a crucial component of China’s energy sector. However, the persistent exploitation of coal resources has gravely impacted the ecological environment. While the Remote Sensing Ecology Index (RSEI) is predominantly used for assessing ecological quality, its primary focus has been urban or aquatic environments. There is limited research focused on the evaluation of the ecological environment quality in mining areas. Moreover, the information regarding surface deformation caused by coal mining extraction is an essential factor in the ecological monitoring of mining areas. Therefore, this study proposed the Modified Remote Sensing Ecology Index (MRSEI). This enhanced model merges active and passive remote sensing techniques and incorporates a deformation factor (Surface Deformation Index, SDI) to provide a holistic evaluation of mining area ecologies. Furthermore, for comparative verification, we developed the Eco-environmental Quality Index (EQI) model by selecting 12 ecological parameters and employing a hierarchical analysis. The Juye mining area in Shandong Province was selected as the region of study. MRSEI results from 2015 to 2021 indicate a decline in the ecological quality of the Juye mining area, with MRSEI values registering at 0.691, 0.644, and 0.617. The EQI model mirrors this decreasing trend over the same period. Despite MRSEI using fewer indicators, its assessments align closely with the multi-indicator EQI method. This validates the accuracy of the MRSEI method, providing reliable technical support for the monitoring and evaluation of ecological environment quality in mining areas.


Introduction
The coal industry plays a pivotal role in China's energy sector, significantly contributing to the nation's economic and social advancement.However, coal mining processes can inflict substantial harm on the ecological environment.For instance, it often leads to the depletion of surface and groundwater levels, diminishing surface vegetation cover, and crop failure.Additionally, it predisposes the rock layers to deformations, precipitating various geological disasters such as surface subsidence, sinkholes, waterlogging, and landslides, among others [1][2][3][4][5][6][7].Hence, it is imperative to bolster the efforts in monitoring and assessing the mining environment to thwart environmental degradation and ecological damage.Enhancing these measures will foster the recovery and rejuvenation of the ecological environment, facilitating a harmonious coexistence between economic progression and ecological conservation.
Traditionally, ecological environment monitoring in mining areas has heavily relied on ground surveys and sampling.These methods, while essential, come with a set of limitations, such as the inability to retrieve historical data on ecological and environmental elements, as well as being time-consuming, labor-intensive, and costly.Moreover, they are not conducive to regular, extensive area measurements [8][9][10].The advent and progression of satellite remote sensing technology, coupled with advancements in computer technology, have considerably enriched the efficacy and precision of environmental monitoring methods.Enhanced quality, spatial and temporal resolution, and spectral resolution of remote sensing data have been achieved.These technologies offer multiple platforms, higher accuracy, macroscopic insights, and extensive coverage, effectively overcoming the limitations associated with traditional survey and sampling methods, thereby furnishing robust support for environmental monitoring and evaluation.Among various methodologies, the Remote Sensing Ecology Index (RSEI) stands out as a predominant tool for evaluating ecological environment quality [11][12][13][14][15]. Initiated by Xu Hanqiu in 2013 [16], RSEI incorporates four key ecological indicators: greenness (NDVI), humidity (WET), dryness (NDBSI), and heat (LST).Utilizing Principal Component Analysis (PCA) facilitates a comprehensive evaluation of the ecological environment quality in urban areas.In recent years, scholars have adeptly applied RSEI or Modified RSEI in ecological environment monitoring across various cities and waters, yielding commendable and insightful evaluation outcomes [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32].
The existing studies predominantly utilize the Remote Sensing Ecology Index (RSEI) for evaluating ecological quality in urban and water areas, with limited application in mining regions.Moreover, coal mining will cause surface subsidence and geologic disasters, which are very harmful to the ecological environment, so the surface deformation information is an indispensable and important factor for ecological monitoring in mining areas [33,34].Currently, there is a notable absence of comprehensive methods that integrate optical and radar remote sensing data for these evaluations.Addressing this, our study introduces a Modified Remote Sensing Ecology Index (MRSEI), innovatively incorporating both active and passive remote sensing techniques.An additional Surface Deformation Index (SDI) is also integrated into the traditional RSEI model, enhancing its applicability in assessing the ecological impacts within mining areas.Additionally, seven ecological factors were extracted from meteorological and statistical data, including PM2.5, net primary productivity, land use intensity, annual precipitation, average temperature, population density, and per capita GDP.These factors were integrated with five ecological indices from the MRSEI: NDVI, WET, NDBSI, LST, and SDI.Utilizing the Analytical Hierarchy Process (AHP), we calculated the weights of each index, enabling the construction of the EQI model.This model was instrumental in comparing and validating the MRSEI method.A comprehensive analysis was performed to explore the spatial and temporal variations of ecological elements in mining areas.
In this research, the Juye Coalfield, encompassing an area of 1210 square kilometers, was meticulously chosen as the focal study area.The applicability of the method was validated by comprehensively utilizing six scenes of Landsat 8 remote sensing imagery and fifty-eight scenes of Sentinel 1A imagery.The subsequent sections are organized as follows: Section 2 provides a detailed description of the study area and dataset; Section 3 meticulously introduces the MRSEI and the EQI models; Section 4 elucidates the findings of the research; Section 5 engages in a discussion analyzing the strengths and weaknesses of the MRSEI and EQI; and Section 6 concludes by summarizing the key outcomes of the study.

Study Area
The Juye Mining Area is strategically situated in the southwestern part of Shandong Province, nestled within the boundaries of Heze and Jining cities.It lies south of the Huaihe River Basin, an area characterized by a rich network of rivers and a complex geological structure predominated by grid-pattern faults running predominantly east-west and north-south.This mining locale comprises seven substantial minefields: Yuncheng, Guotun, Zhaolou, Longgu, Wanfu, Pengzhuang, and Liangbaosi Coal Mines, collectively spreading over an area of 1210 km², as shown in Figure 1 [35].Remarkably, the Juye Mining Area boasts a significant production capacity, yielding approximately 18 million tons of coal annually.Since its inception in 2007, 16 years of operational history have unfortunately led to notable ecological degradation in the area.The repercussions of extensive mining activities manifest in various environmental detriments such as groundwater contamination, substantial land and ground subsidence, notable vegetation damage, and pervasive heavy metal pollution.These adverse ecological impacts have notably compromised the quality of life and health of the local communities and pose significant challenges to national ecological and environmental sustainability objectives.Therefore, it is essential to monitor the ecological environment quality of the Juye mining area, providing necessary support and assistance for subsequent ecological restoration efforts.

Datasets
The extraction of index data for NDVI, WET, NDBSI, and LST was accomplished using six scenes of Landsat 8 remote sensing imagery with less than 5% cloud cover, spanning from May 2015 to May 2021.Detailed parameters of these images are presented in Table 1.
For the SDI index, data were extracted from fifty-eight scenes of Sentinel-1A images in IW mode, collected between March 2015 and December 2021.A comprehensive list of this data is displayed in Table 2.The DEM data utilized SRTM data with a horizontal resolution of 90 m, which is essential for removing the terrain phase in the time-series differential interferometry.ESA POD precision orbit determination ephemeris data were employed for orbit information correction.Additionally, vector administrative boundaries, vector boundaries of the Juye Mining area, county statistical yearbooks of Shandong Province, and meteorological data from 2015 to 2021 were utilized to augment the research on ecological environment monitoring in the mining area.

Methods
This paper addresses the deficiency of comprehensive monitoring and evaluation methods that combine optical remote sensing data and radar remote sensing data in the assessment of the ecological environment in mining areas.A MRSEI model is proposed, incorporating unique deformation factors endemic to the ecological environment of mining areas, based on the foundation of the Remote Sensing Ecology Index (RSEI).The AHP is utilized to construct a model grounded in the EQI, facilitating a comprehensive analysis of the spatiotemporal variation characteristics of ecological elements within mining areas.
Initially, data from Landsat 8 and Sentinel 1A are processed to extract essential ecological indicators such as greenness, humidity, dryness, heat, and deformation factors.A principal component analysis is applied to establish the MRSEI model.Subsequently, leveraging statistical and meteorological data, seven ecological elements are extracted: PM2.5, net primary productivity, land use intensity, annual precipitation, average temperature, population density, and per capita GDP.These elements are then combined with five ecological indicators from the MRSEI: NDVI, WET, NDBSI, LST, and SDI.Lastly, a comparative analysis is conducted to discern the relative merits and demerits of the two methodologies, enabling a comprehensive evaluation of the ecological environment quality in the Juye mining area.The specific technical approach adopted in this research is depicted in Figure 2.

Modified Remote Sensing Ecological Index Model (MRSEI)
The prevailing RSEI is primarily adapted to monitor and evaluate urban environments.However, due to the intrinsic differences between mining and urban areas, the mining processes inevitably lead to surface subsidence or even collapse in mining regions.Consequently, the applicability of RSEI in monitoring mining areas is somewhat limited.Considering the specific conditions of the Juye mining area, this study modifies four ecological indicators initially included in the RSEI model, incorporating surface deformation indicators, and formulates the MRSEI to assess the ecological environmental quality of the Juye mining area.

Masked Waters
The water body has high reflectivity and low radiation temperature and is greatly different from vegetation and soil.In order to avoid the impact of the water body on the calculation results, the Modification of Normalized Difference Water Index (MNDWI) is used [36].Create a mask file to remove the water body.The formula is as follows: For regions with average water resources, the range of water bodies is typically between [0.25, 1].Therefore, this article sets 0.25 as the water body threshold, encompassing the majority of water bodies.

Extraction of Traditional Remote Sensing Ecological Indicators
Traditional remote sensing ecological indicators include the greenness index (NDVI) [37], humidity index (WET) [38], dryness index (NDDBSI) [39], and heat index (LST) [40].The principle is as follows: Extracting the greenness index (NDVI).The greenness index reflects the growth status of surface vegetation and can effectively distinguish between vegetation and bare land.The greenness index uses the Normalized Difference Vegetation Index (NDVI), and its formula is: where NIR  and red  represent the reflectivity of the red band and the near-infrared band.The value of NDVI is usually between [−1,1].The higher the NDVI value, the better the vegetation coverage.

Extracting the humidity index (WET).
The humidity index is used to monitor surface humidity and evaluate the water and heat balance in a certain area.This paper uses the tassel cap transformation method to obtain the humidity index.The calculation formula is: Extracting the dryness index (NDBSI).Normalized Difference Build and Soil Index (NDBSI) is an index that reflects the distribution of surface bare soil and buildings.It is composed of Soil Index (SI) and Index-based Built-up Index (IBI) synthesis; the calculation formula is: ) / ( ) Extracting the heat index (LST).The heat index is an index that reflects the surface temperature.By studying the distribution and changes of the surface temperature, it provides a scientific basis for thermal environment monitoring and ecological environment assessment.This article uses the radiative transfer equation method to invert the surface temperature.The formula is as follows: where 10 L is the radiance value of the 10th band of the Landsat TIRS image;  is the landmark specific radiance; () B LST is the blackbody radiance value; L  , L  are the upward and downward radiances of the atmosphere respectively;  is the atmospheric transmittance;

Extraction of Deformation Indicators
The procedure for extracting deformation indicators is illustrated in Figure 3. Data preprocessing.Calibrate and position the original Sentinel-1A data to obtain SLC data.Geocode the DEM to obtain the DEM data in the radar coordinate system, which is used to remove the terrain phase.Add precision orbit parameter files to accurately position and Correction.Due to the limited number of Sentinel-1A images in 2015, the Stacking-InSAR technique was employed to extract deformation information.For the years 2018 and 2021, the SB-InSAR method was utilized for this purpose.
Image registration.The super master image is selected through overall correlation coefficient analysis.All secondary images are registered and resampled to the super primary image.
Stacking-InSAR Processing.The Stacking-InSAR experiment used a total of 10 scenes of Sentinel-1A images from two orbits covering the Juye mining area, with the time span from 30 July 2015 to 3 November 2015.The Goldstein adaptive filtering method was used to filter the registered images, and the minimum cost flow method was used for phase unwrapping.Finally, Stacking-InSAR was used to obtain the deformation rate in the LOS direction of the Juye mining area.
SB-InSAR Processing.This SB-InSAR experiment uses a total of 48 Sentinel-1A images from an orbit covering the Juye mining area.The time span is from 3 January 2018 to 29 December 2018 and from 11 January 2021 to 13 December 2021.(1) Small baseline interference atlas generation.On the premise of ensuring that there are no isolated points in the spatiotemporal baseline network, the vertical baseline and time baseline thresholds are set to construct a small baseline interference atlas.Use DEM for simulation and remove terrain contribution.(2) Surface deformation rate inversion.In the StaMPS software environment, high coherence points are selected based on coherence, 3D phase unwrapping is performed on the extracted high coherence points, and the orbit residual error, atmospheric delay error, and DEM error are estimated and corrected, and finally, the least squares estimation is used Invert surface deformation rate information in the LOS direction of high coherence points.
Extract the SDI indicator.The surface settlement rate of LOS in the Juye mining area is the deformation index, and the formula is as follows: ( ) SDI is the deformation index, and () DX refers to the deformation rate of the pixel in the LOS direction.

Construction of Modified Remote Sensing Ecological Index Model
Principal Component Analysis.The original data is converted into a new set of independent variables through linear transformation.These new variables are called principal components.When the contribution rate of the principal component is greater than 85%, it is considered valid data [16].If the contribution rate of the principal component does not exceed 85%, it is recommended to retain more principal components at this time to better retain the information in the original variables.In this experiment, the contribution rate of the first principal component is less than 85%, so the first two principal components are retained.The formula for principal component analysis is as follows.
  , , , X is a matrix with n rows and m columns.
Perform a linear transformation of the formula and compute the contribution rate of each principal component.
where i F represents the linear transformation formula; R is the contribution rate; m is the first few principal components; and p is the number of principal components.
Indicator normalization.Since the scale of the five factors is not uniform, performing PCA directly could result in a disproportionate weighting of the indicators.Therefore, it is necessary to normalize the five indicators before conducting the principal component analysis.The linear normalization formula is as follows: Indicator synthesis.The traditional RSEI method only uses the first principal component to construct the RSEI model and conduct research on regional ecological environment assessment.Due to the poor overall ecological environment of the Juye mining area, using the first principal component cannot ensure that most of the image information is reasonably utilized.Therefore, in order to include more effective information, this paper retains the first two principal components to construct MRSEI.The formula is:  ) ( ) where MRSEI is the modified remote sensing ecological index, the range is [0,1] , the larger the value of MRSEI , the better the quality of the ecological environment, 0,min

RSEI
is the minimum value of

RSEI
is the maximum value of 0 RSEI 的最小值 .

Selection of Evaluation Indicator Factors
The impact factors of changes in ecological and environmental quality are very complex.Referring to the "Technical Specifications for Ecological and Environmental Condition Assessment", 12 indicators are selected [41].These 12 indicators are divided into three layers: the target layer, the element layer, and the indicator layer (Table 3).

Analytical Hierarchy Process to Establish Evaluation Index Weights
Establishment of a hierarchical structure.The selected 12 factor indicators are stratified, as shown in Table 3.
Constructing the judgment matrix.The order of the relative importance of each evaluation indicator is ranked according to the 9-point ratio, specifically divided into nine levels, respectively: 1, 3, 5, 7, 9, 7, 5, 3, 1.The judgment matrix D of the evaluation indicators is constructed according to the order, and the specific differences in the degree of importance between the elements are shown in Table 4.
Create a judgment matrix based on the above Table 4: a , a , a , , a , ,1 where A is judgment matrix; anj represents the indicator of the n-th and j-th.Calculate the maximum eigenvalue: where max  is the maximum eigenvalue; A is the judgment matrix; n represents the num- ber of indicators; i  represents the i-th eigenvalue.In order to ensure the reliability of the calculation results, CR consistency testing needs to be done.The calculation formula is as follows: where CR is the consistency result; CI is a consistency indicator; RI is a randomized consistency indicator; max  is the maximum eigenvalue; n represents the number of indica- tors.CR < 0.1, yielding the result that this judgment matrix has a more reliable consistency, where RI was obtained from Table 5. Hierarchical general ordering.The result of the total ranking of levels reflects the relative importance between the evaluation indicators and is a key step in determining the evaluation results of the AHP method (Table 6).The formula for checking the consistency of the total hierarchical ordering is: where CRtotal is the consistency of the total hierarchical ordering; CIj is the j-th consistency indicator; RIj is the j-th randomized consistency indicator; aj represents the total hierarchical result of the indicator layer on the target layer.Determination of the final weight for each indicator.The judgment matrix is tested for consistency in the overall ranking of the hierarchy and then averaged and normalized to obtain the final weight value.

Subsubsection
Evaluation of ecological environment quality should comprehensively consider the contribution of the above factors, assigning values to each evaluation index and multiplying them by the weights to be accumulated, and the accumulated result is the EQI.The formula is as follows: In the formula, i V is the normalized value of the i -th evaluation index, w i is the weight of the i -th index, and n is the number of evaluation indicators.Since LST, NDBSI, SDI, population density, per capita GDP, and PM2.5 are inversely related to the merits of ecological quality, their normalized values need to be transformed to the opposite number, i.e., 1− i V , in the arithmetic process.However, since the dimensions of each ecological environment parameter are different, it is necessary to normalize them by removing the dimensions.

Results and Analysis
In this study, reliable results for the greenness, humidity, dryness, and heat indices were achieved through meticulous radiometric and atmospheric corrections.The average deformation rate and standard deviation in the non-deformation areas were 2.47 mm/a and 1.17 mm/a, respectively.These indicators provide a solid foundation for the credible construction of the MRSEI model.

Principal Component Analysis Results
Table 7 shows the results of the principal component eigenvalues.The table shows the first principal ingredient contributes more than 55% and has the highest contribution in 2015, 2018, and 2021.This is followed by the second principal ingredient, both of which exceed 15%.The traditional Remote Sensing Ecological Index (RSEI) proposes that the first principal ingredient can be used to establish the RSEI when it can reach 85% or higher.However, because of the specificity of the Juye mining area, it is not sufficient to use only the first principal ingredient to represent the quality of the ecological environment.So, we use the first principal ingredient combined with the second principal ingredient to construct MRSEI.  8. From the table, the correlation coefficients of NDVI and WET of the first principal ingredient in each year of the Juye mining area are positive and show a positive correlation with MRSEI, which has a certain effect on the ecological environment.The NDBSI, LST, and SDI correlation coefficients are all negative, showing a negative correlation with the MRSEI, which is inhibitory to the ecosystem.NDVI has the greatest impact on the process of ecological quality change in the Juye mining area.Improving the quality of the ecosystem in mining areas requires an increase in vegetation cover and an increase in the vegetation index cover.From the negative correlation index, the influence of the NDBSI index is more obvious in the process of ecological environment quality coercion in the mining area, and the SDI index gradually increases with time and occupies a larger and larger proportion.

Analysis of MRSEI Model Results
Table 9 and Figure 4 present the statistical values and the spatio-temporal distribution of the MRSEI index in the Juye mining area, respectively.The MRSEI values are distributed within an interval ranging from 0 to 1, where higher values correspond to a better ecological environment.As depicted in Table 9, the mean MRSEI values for the Juye mining area exhibit a declining trend from 2015 to 2021, with average values recorded as 0.691, 0.644, and 0.617 over the period.Specifically, there was a decrease of 6.80% in the mean MRSEI from 2015 to 2018, and a further reduction of 4.19% was observed between 2018 and 2021.Figure 4 illustrates that, from 2015 to 2021, areas within the Juye mining region exhibiting poorer ecological quality were predominantly situated in industrial and mining lands, followed by town construction areas and settlements.Conversely, areas characterized by higher MRSEI values-indicative of better ecological quality-were mainly those with a distribution of woodlands, grasslands, and croplands.

Evaluation of the Ecological Environment Quality Based on the MRSEI Model
Based on the regional division, the Juye mining area can be divided into seven regions (Table 10).As shown in Table 10, the mean MRSEI values were calculated for each mining area to analyze the changes in the ecological quality of the mining area.
From 2015 to 2021, the average MRSEI values of the Liangbaosi and Pengzhuang mining areas were higher than the overall average values in the Juye mining area.In these two mining areas, the MRSEI index exhibited a trend of a slight increase followed by a decline.According to visual interpretation results, it is observable that these areas have a smaller expanse of construction and industrial land, resulting in more abundant vegetation cover and superior ecological quality within the Juye mining district.Conversely, the mean MRSEI values for the Yuncheng and Guotun mining areas were below the Juye mining area's average.MRSEI values in these regions display a decreasing trend, with all ecological indicators showing signs of deterioration.This decline is attributed to an increase in bare land due to extensive mining operations.For the Zhaolou, Longgu, and Wanfu mining areas, the MRSEI index generally exhibits a fluctuating decline, indicating that the ecological quality in these regions is of a medium level compared to other parts of the Juye mining area.To more directly represent the fundamental state of the ecological environment, the MRSEI values were categorized into five grades based on equal intervals, following the Technical Specification for the Evaluation of Ecological Environment Condition issued by the State in 2015 (see Table 11).This categorization yielded the MRSEI grading map of the Juye Mining Area (refer to Figure 5).Analysis of Figure 5 reveals that regions with very poor and poor grades of ecological environment quality are primarily concentrated in industrial and construction lands.Medium-grade areas are predominantly located on the periphery of industrial and urban construction zones, while areas with less premium and premium grades are mainly distributed in cultivated, forested, and grassland areas.Over the period from 2015 to 2021, it is evident that regions with poor grades of the ecological environment are expanding, signifying a worsening ecological condition, primarily centered around the Yuncheng coal mine, Guotun coal mine construction zone, and Longgu coal mine subsidence area.Statistical analysis of the area and percentage for each MRSEI grade was conducted, as detailed in Table 12.From 2015 to 2021, areas rated as very poor initially expanded from 2.93% to 8.21%, subsequently decreasing to 3.99%.The areas categorized as poor, medium, and slightly above average observed a consistent increase, while areas classified as premium experienced a continuous decline.This suggests that in the Juye mining area, the quality of the ecological environment has been perpetually diminishing due to ongoing mining activities.

Evaluation of the Ecological Environment Quality Based on the EQI Model
The results of the EQI model evaluation indicators are displayed in Table 13.These findings led to the creation of an ecological environment quality assessment map (Figure 6) and the corresponding statistical index values for the Juye mining area, which are detailed in Table 14.As depicted in Figure 6, areas with poor ecological conditions are primarily concentrated in the construction zones of the Yuncheng and Guotun mines, as well as the subsided and waterlogged areas of the Longgu coal mine.There are also sporadic dots in other areas, mainly within towns and rural settlements.The spatial distribution results are similar to those obtained through the MRSEI model.In Table 14, the mean values for the years 2015, 2018, and 2021 are 0.508, 0.436, and 0.464, respectively.Between 2015 and 2018, there was an initial decrease of 14.17%, followed by an increase of 6.42% from 2018 to 2021.The overall trend indicates a decline, with the magnitude of the downward trend being more significant than the upward trend.This is primarily due to the smaller baseline of the EQI.Following the guidelines of the National Technical Specification for the Evaluation of Ecological Environmental Conditions, the EQI is categorized into five levels within the range of 0 to 1, as outlined in Table 11.The areas designated with each ecological level were quantified by area, and the results are presented in Table 15 and Figure 7.The findings reveal that in 2015, the largest percentage of ecological quality was rated as less premium, while in 2018 and 2021, the medium grade constituted the largest percentage, accounting for 50% or more.The second-highest percentages of the area were allocated to the medium grade in 2015 at 23.32%, the less poor grade in 2018 at 18.75%, and the less premium grade in 2021 at 33%.The combined percentages of both the less poor and poor grades exhibit an increasing and then decreasing trend.These areas are primarily characterized by land cover types such as urban and industrial areas, and the ecological quality is not conducive to human habitation.In contrast, the area designated as premium grade decreased from about 10% in 2015 to less than 1% in 2018 and 2021.
Comparing the grading statistics of the MRSEI and EQI methods reveals differences in ecological assessment.In 2015, MRSEI identified a higher percentage of areas as premium and less premium grades, totaling 67%, while the EQI method classified the majority as less premium and medium grades at 76%.By 2018, MRSEI still categorized a significant portion as less premium and premium, totaling 59%, showing an increase in the less premium grade compared to 2015.Meanwhile, EQI focused on medium and less poor grades, totaling 76%.In 2021, MRSEI results shifted, highlighting the less premium and medium grades, comprising 75% of the total, showing a shift away from the dominance of premium and less premium grades seen in previous years.EQI, on the other hand, reverted to prioritizing medium and less premium grades.The MRSEI results suggest a better ecological ranking predominantly in the central part of the Juye mining area, characterized by smaller urban and industrial expanses compared to the more prominent cultivated lands and forested and grassland areas.Between 2015 and 2021, the encroachment of urban and industrial areas upon cultivated lands and natural habitats has led to land degradation and a continuous decline in ecological quality.EQI results corroborate the MRSEI findings, confirming a general trend of ecological deterioration in the area.Thus, both methods concur on the declining trajectory of ecological quality in the region.

Discussion
This section comprehensively compares and analyzes the MRSEI and EQI methods, detailing the advantages and disadvantages of each.
In terms of the temporal distribution, it can be inferred from the MRSEI method that the ecological environment has experienced a declining trend.Similarly, the EQI method also indicates a general trend of ecological quality deterioration.The downward trend in EQI is more pronounced than that in MRSEI, with a greater magnitude.This discrepancy can be attributed to the smaller baseline of EQI.Nonetheless, the absolute value of changes in both methods follows a similar trend, falling within the range of 4% to 15% or less.
In terms of spatial distribution, both methods exhibit similar patterns.Regions with poor ecological conditions are concentrated in the construction areas of the Yuncheng and Guotun mines, as well as in the subsided and waterlogged areas of the Longgu coal mine.In other areas, sporadic points are observed, primarily within towns and rural settlements.However, the MRSEI method demonstrates more pronounced discrepancies in the data points compared to the EQI method, which is more coherent.This difference may be attributed to the greater number of indicators used in the EQI method.
The spatial and temporal distribution maps of the ecological environment in the Juye mining area reveal similarities between the MRSEI and multi-indicator AHP methods, even with MRSEI employing fewer indicators.When comparing these methods, it appears that the MRSEI method can effectively reflect the regional ecological environmental status to a certain extent.It rationally links key ecological elements with remote sensing characteristics, utilizing remote sensing factors such as vegetation cover, surface humidity, aridity, surface temperature, and surface deformation as evaluation indices.This approach facilitates a convenient and effective estimation of the regional ecological environment's quality, enabling the execution of ecological environment quality assessments purely based on remotely sensed data.The MRSEI model indicators are stable and readily accessible, possessing a degree of applicability.Parameters can be adjusted as necessary to accommodate varying research requirements, making the inclusion of surface deformation factors more adaptable to the mining area's modeling prerequisites.However, these indicators can be sensitive to changes within the study area, influencing the corresponding estimation models' accuracy.In terms of vegetation parameters, the model heavily relies on NDVI due to a limited selection of indicators.Furthermore, the model predominantly depends on remote sensing indicators in its evaluation criteria.A broader selection of indicators, incorporating an array of ecological parameters, could better represent the mining area's ecological environment.
The EQI method features a richer set of evaluation indicators and a clear hierarchical structure, making it accessible even to non-specialists.The model is stable enough, and, with the support of NDVI, the selection of each evaluation indicator is more flexible and holds certain ecological significance.However, the process of collecting indicators is somewhat cumbersome, and there is a high level of subjectivity involved in the selection and weighting of the indicators in the model.

Conclusions
This study focuses on the Juye mining area in southwest Shandong Province, utilizing Landsat 8 OLI and Sentinel-1A data from 2015, 2018, and 2021 for basic data processing and quantitative analysis.Five indicators-NDVI, WET, NDBSI, LST, and SDIwere extracted and utilized in establishing the MRSEI through principal component analysis.Additionally, the EQI was constructed using ecological factors such as PM2.5, NDVI, NPP, population density, and GDP per capita.The aim is to explore and analyze the changing patterns and influencing factors of ecological and environmental quality in the Juye mining area from 2015 to 2021.The following conclusions were derived: (1) The image is obtained by normalizing five indicators: the greenness indicator, humidity indicator, dryness indicator, heat indicator, and deformation indicator.These indicators undergo principal component transformation.The cumulative contributions of the first and second principal components of the MRSEI exceed 85%.This indicates that PC1 and PC2 capture the majority of information from the five indicators.Specifically, the greenness, humidity, and dryness indicators have positive impacts on ecological, environmental quality and the environment.In contrast, the heat and deformation indicators negatively affect ecological environmental quality.
(2) Analysis of Ecological Environment Quality in the Juye Mining Area Based on the MRSEI Model.From 2015 to 2021, the MRSEI values in the Juye mining area displayed a decreasing trend, with values of 0.691 in 2015, 0.644 in 2018, and 0.617 in 2021.In 2015 and 2018, the areas were primarily graded as less premium and premium in terms of ecological environmental quality.However, by 2021, the dominant grades shifted to less premium and medium, indicating a deterioration in the ecological environment quality over this period.Liangbaosi and Pengzhuang coal mines exhibit the highest MRSEI values, reflecting the best ecological environment quality within the area.Conversely, the Guotun coal mine area has the lowest ecological environment quality within the Juye mining area.
(3) Analysis of Ecological Environment Quality in the Juye Mining Area Based on the EQI Model.From 2015 to 2021, the EQI values in the Juye mining area displayed a decreasing trend.The dominant ecological environment quality grades in each year are medium and less premium, collectively accounting for over 70% of the study area.Both the poor and very poor grades show an initial increase followed by a decrease in area proportion, with their main land cover types being building land and industrial and mining land.The ecological quality of these areas is not conducive to human life.Additionally, the area of excellent grades sharply declined from about 10% in 2015 to less than 1% in 2018 and 2021.
(4) Comparison reveals that the MRSEI method is simple, objective, and easy to implement but heavily reliant on remote sensing indicators.On the other hand, the EQI method, with its rich evaluation indicators and clear hierarchical structure, demands higher standards for indicators and involves a certain degree of subjective judgment.

Figure 1 .
Figure 1.Schematic diagram of the geographical location of the study area.(a) The study area location.(b) The mining area location.

Figure 2 .
Figure 2. Workflow of the proposed method.
of the green band and the mid-infrared band 1.
of the blue band, green band, red band, mid-infrared band 1, and mid-infrared band 2.


are the reflectivity of the blue band, green band, red band, and mid-infrared band 1.The value of NDBSI is between [−1, 1].The higher the value of NDBSI, the larger the area of buildings and bare soil areas.

2 K 1 K  , 2 K
 is the calibration parameter, for the 10th band of the Landsat image TIRS sensor has  values of 480.89 and 1201.14, and LST which is the surface temperature.

Figure 3 .
Figure 3. Flowchart for extraction of deformation indicators.
eigenvalue of the principal component; where i NI is the result of the indicator value; i I is the original value of the corresponding indicator pixel i ; and max I and min I are the maximum and minimum values of the cor- responding indicator, respectively.

where 0 RSEI 2 PC
is the initial remote sensing ecological index; 1 PC represents the first principal component; represents the second principal component; p represents the eigenvalue of 1 PC ; and q represents the eigenvalue of 2 PC .In order to facilitate the measurement and comparison of the size of each indicator, linear normalization 0 RSEI is performed to obtain MRSEI.The formula is as follows:

Figure 4 .
Figure 4. Spatial and temporal distribution of MRSEI in the Juye mining area.

Figure 5 .
Figure 5. Grading of MRSEI in the Juye mining area.

Figure 6 .
Figure 6.Spatial and temporal distribution of ecological environment quality index in the Juye mine area.

Figure 7 .
Figure 7. Grading of EQI in the Juye mining area.

Author Contributions:
Conceptualization, Y.C. and Z.S.; methodology, H.L.; validation, Y.C., Z.S., and H.L.; formal analysis, Y.C; investigation, H.L; writing-original draft preparation, Z.S.; writing-review and editing, Y.C. and Z.S.; visualization, H.C., and Q.L.; supervision, Y.C.; project administration, Y.C.; funding acquisition, Y.C.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by the National Natural Science Foundation of China, grant number 42171312, the Open Fund of Key Laboratory of Geographic Information Science (Ministry of Education), East China Normal University, Grant No. KLGIS2023A02.Data Availability Statement: Publicly available datasets were analyzed in this study.

Table 3 .
Ecological environment quality evaluation index system.

Table 4 .
Weight factor judgment table.

Table 5 .
Randomized consistency indicator reference table.

Table 6 .
Total weight factor judgment table.

Table 7 .
Statistical table of the results of the eigenvalues of the principal components.

Table 8 .
Principal component analysis table for various indicators.

Table 9 .
Statistical results of the MRSEI.

Table 10 .
Mean MRSEI value for the mining area.

Table 11 .
Ecological environmental quality grading standards.

Table 12 .
Statistical results of MRSEI grading in the Juye mine area.

Table 13 .
Results of the EQI model evaluation indicators.

Table 14 .
Statistical results of the EQI.

Table 15 .
Statistical results of EQI grading in the Juye mine area.