Quantitative Evaluation of the Eco-Environment in a Coalfield Based on Multi-Temporal Remote Sensing Imagery: A Case Study of Yuxian, China

With the exploitation of coalfields, the eco-environment around the coalfields can become badly damaged. To address this issue, “mine greening” has been proposed by the Ministry of Land and Resources of China. The sustainable development of mine environments has now become one of the most prominent issues in China. In this study, we aimed to make use of Landsat 7 ETM+ and Landsat 8 OLI images obtained between 2005 and 2016 to analyze the eco-environment in a coalfield. Land cover was implemented as the basic evaluation factor to establish the evaluation model for the eco-environment. Analysis and investigation of the eco-environment in the Yuxian coalfield was conducted using a novel evaluation model, based on the biological abundance index, vegetation coverage index, water density index, and natural geographical factors. The weight of each indicator was determined by an analytic hierarchy process. Meanwhile, we also used the classic ecological footprint to calculate the ecological carrying capacity in order to verify the effectiveness of the evaluation model. Results showed that the eco-environment index illustrated a slowly increasing tendency over the study period, and the ecological quality could be considered as “good”. The results of the evaluation model showed a strong correlation with the ecological carrying capacity with a correlation coefficient of 0.9734. In conclusion, the evaluation method is a supplement to the time-series quantitative evaluation of the eco-environment, and also helps us to explore the eco-environment in the mining area.


Introduction
China has been experiencing rapid urban development with the economic reform taking place over the past 20 years. Coal has become one of the most important power sources and accounts for about 70% of China's primary energy consumption [1]. The original balance of the rock stratum is inevitably changed after the exploitation of coal resources. The collapse of the Earth's surface caused by mining activities can trigger a series of environmental disasters [2]. To address this issue, "mine greening" has been proposed by the Ministry of Land and Resources of China [3], which is a new idea for mineral resources management. Huang [4] constructed the evaluation index system for "mine greening" based on legitimate mining, efficient utilization, cleaner production and standardized

Study Area
The Shanhou and Cuijiazhai coalfield is located between 114 • 24 40"E-114 • 32 30"E and 39 • 50 54"N-39 • 58 25"N in Yuxian County, Zhangjiakou City, Hebei Province, China, shown in Figure 1. This area belongs to the warm temperate continental monsoon climate zone, and the annual precipitation is between 380 and 683 mm. The monsoon climate generates four distinct seasons in Yuxian County. A cool summer, changeable autumn, and lower temperature in winter are the characteristics of the climate in Yuxian County. The temperature in winter can occasionally reach −20 • C. Coal industry is the most important industry in Yuxian County, resulted from many superiorities of this coalfield, such as the shallow burial depth, low slope, and low water inflow in the interior of the coalfield. However, some problems have arisen with the exploitation of the coalfield, including ground fissures and collapse (see Figure 2). The mining exploitation has also increased the likelihood of geological disasters. Moreover, sustained precipitation may further aggravate the geological disasters and cause casualties or property loss.  Figure 2). The mining exploitation has also increased the likelihood of geological disasters. Moreover, sustained precipitation may further aggravate the geological disasters and cause casualties or property loss.

Data Sources and Processing
A complementary set of land-cover observations were provided by the Operational Land Imager (OLI) and Enhanced Thematic Mapper Plus (ETM+) sensors. The Landsat Earth resources satellite system was the first system designed to provide near global coverage of the Earth's surface on a regular and predictable basis [33]. The 30-m spatial resolution of these sensors makes it possible to extract surface information such as land cover, vegetation distribution, and urban impervious  Figure 2). The mining exploitation has also increased the likelihood of geological disasters. Moreover, sustained precipitation may further aggravate the geological disasters and cause casualties or property loss.

Data Sources and Processing
A complementary set of land-cover observations were provided by the Operational Land Imager (OLI) and Enhanced Thematic Mapper Plus (ETM+) sensors. The Landsat Earth resources satellite system was the first system designed to provide near global coverage of the Earth's surface on a regular and predictable basis [33]. The 30-m spatial resolution of these sensors makes it possible to extract surface information such as land cover, vegetation distribution, and urban impervious

Data Sources and Processing
A complementary set of land-cover observations were provided by the Operational Land Imager (OLI) and Enhanced Thematic Mapper Plus (ETM+) sensors. The Landsat Earth resources satellite system was the first system designed to provide near global coverage of the Earth's surface on a regular and predictable basis [33]. The 30-m spatial resolution of these sensors makes it possible to extract surface information such as land cover, vegetation distribution, and urban impervious surfaces. With the data continuity of Landsat over the past decades, Landsat data based time-series analysis allows for an effective characterization of the temporal and spatial variation of the eco-environment [34][35][36]. The images used in this study were downloaded from the U. S. Geological Survey (USGS) (http://earthexplorer.usgs.gov/) and the Geospatial Data Cloud (http://www.gscloud.cn/). The acquisition period of these data was from 2005 to 2016. Multi-temporal analysis of Landsat imagery presents challenges in terms of data availability. The L1TP products were utilized in our work and the atmospheric correction was done by FLAASH model which has been demonstrated to be sufficiently accurate for the image classification [37][38][39][40].
The images in 2010 were deficient because of the heavy cloud cover, so we skipped this year. Finally, we collected 44 scenes, covering the March, July, September, and December of each year. Given the data availability, we replaced the data which were not available with data from adjacent months. The data list is listed in Table 1. Because the original remote sensing data had only been processed with rough radiometric correction and geometric correction, it was necessary to preprocess the data, including clipping, re-sampling, and projection conversion. The Landsat calibration module and the Fast Line-of sight Atmosphere Analysis of Spectral Hypercube (FLAASH) algorithm provided in ENVI5.1 software (Harris Geospatial Solutions, Broomfield, CO, USA) was used to achieve atmospheric correction. The geometric registration root-mean-square (RMS) deviation was found to be less than 0.5 pixels, when compared with field investigation. Finally, the images were converted to the World Geodetic System-1984 (WGS-84) coordinate system and Universal Transverse Mercator (UTM) projection. The ETM+ sensor Scan Line Corrector (SLC) malfunctioned in May 2003, so the captured images contained missing data strips after 2003. We therefore applied the integrated algorithm proposed by Zeng et al. [41] to repair the missing data strips.
In addition, auxiliary data were also integrated into the environmental quality evaluation. Meteorological data were provided by the China Meteorological Data Service Center of the China Meteorological Administration (http://data.cma.cn/). Statistical data were acquired from Statistical Yearbooks of China's counties and cities, which were downloaded from the National Bureau of Statistics of China (http://www.stats.gov.cn/tjsj/).

Construction of the Environmental Quality Evaluation Index System
The evaluation of environmental quality is a multi-disciplinary and multi-technology task, which combines ecology, statistics, earth observation, etc. After comprehensive consideration of the data availability and the natural conditions of the Yuxian County coalfield, the index system was established based on the "Technical specifications of eco-environment evaluation (for trial implementation)" (HJ/T192-2006). Specifically, the individual indices were made up of the BAI, VCI, WDI, and NGF [42,43]. The main influencing factors of these four major aspects were cultivated land, water body, grassland, building land, bare land, digital elevation model (DEM), precipitation, vegetation condition, and surface runoff area, the first five of which were obtained by classification of Landsat imagery.
The impact of mining activities on the surrounding environment is mainly manifested in the change of land cover. On the other hand, land cover is an important environmental variable and a key input parameter for ecological models [34]. The ROI is divided into training set and test set when the supervised classification is performed. The overall accuracy of the classification in this study was over 85% verified by the test set, which can satisfy the demand of eco-environment evaluation. The accuracy assessment on 2016 as the example is listed in Table 2. In addition, vegetation coverage is also a key parameter for environmental evaluation, as vegetation plays an important role in soil and water conservation, environmental purification, and oxygen exchange. Due to the fact that abnormal precipitation and the topographic gradient are both significant factors which can cause geological disasters in mining areas, both of these factors were selected as representative indices. The weight is a quantitative coefficient of the importance distribution of the impact factor of the research object [44]. AHP is a multi-objective decision analysis method which combines quantitative analysis with qualitative analysis [31]. It can be combined with fuzzy mathematics to establish a membership function for a correlation analysis. When the indicators are few in number and a high evaluation accuracy is required, AHP can be regarded as a good choice. The assigned weight is normalized on the basis of Saaty's scale, considering two themes and classes at a time, on the basis of their relative importance, to determine the eco-environmental index (EI). Thereafter, pair-wise comparison matrices of the weights assigned to the different thematic layers and their individual classes are constructed using Saaty's AHP, and the weights are normalized by an eigenvector approach. The consistency ratio (CR) is then calculated to examine the normalized weights of the various thematic layers and their individual classes, as per the recommendation of Saaty [28]. The following steps are carried out to compute the CR of the various thematic layers and their individual classes.
(1) Hierarchy Model Construction The HJ/T192-2006 was promulgated by the national environment protection bureau of China which specified the eco-environmental assessment indexes and the corresponding computing methods. In HJ/T192-2006, EI is utilized to evaluate the regional eco-environment quality with a range of 0 to 100 comprehensively. Combining the indexes given by HJ/T192-2006 and the field investigation, we used the BAI, VCI, WDI and NGF as the indicators of criterion layer to established the hierarchy model. This lays an important foundation for the environmental quality evaluation in the coalfield. The hierarchy model is shown in Figure 3. The weight is a quantitative coefficient of the importance distribution of the impact factor of the research object [44]. AHP is a multi-objective decision analysis method which combines quantitative analysis with qualitative analysis [31]. It can be combined with fuzzy mathematics to establish a membership function for a correlation analysis. When the indicators are few in number and a high evaluation accuracy is required, AHP can be regarded as a good choice. The assigned weight is normalized on the basis of Saaty's scale, considering two themes and classes at a time, on the basis of their relative importance, to determine the eco-environmental index (EI). Thereafter, pair-wise comparison matrices of the weights assigned to the different thematic layers and their individual classes are constructed using Saaty's AHP, and the weights are normalized by an eigenvector approach. The consistency ratio (CR) is then calculated to examine the normalized weights of the various thematic layers and their individual classes, as per the recommendation of Saaty [28]. The following steps are carried out to compute the CR of the various thematic layers and their individual classes.
(1) Hierarchy Model Construction The HJ/T192-2006 was promulgated by the national environment protection bureau of China which specified the eco-environmental assessment indexes and the corresponding computing methods. In HJ/T192-2006, EI is utilized to evaluate the regional eco-environment quality with a range of 0 to 100 comprehensively. Combining the indexes given by HJ/T192-2006 and the field investigation, we used the BAI, VCI, WDI and NGF as the indicators of criterion layer to established the hierarchy model. This lays an important foundation for the environmental quality evaluation in the coalfield. The hierarchy model is shown in Figure 3.

Ibai
Water density index

WDI = φwdi × W · Iwdi
Natural geographical factors Biological abundance index Eco-environmental Index

EI = W · [ BAI,VCI,NGF,WDI]
Target Layer Criterion Layer Index Layer The computing methods of all indicators in the index layer are shown in Table 3, the area of different land use is obtained from the post-classification statistics. The weight and the normalization coefficient of each indicator are introduced in next section. The computing methods of all indicators in the index layer are shown in Table 3, the area of different land use is obtained from the post-classification statistics. The weight W and the normalization coefficient of each indicator ϕ are introduced in next section. Pair-wise comparison matrices are used to compare the importance of two indices. The importance of the pair-wise comparison judgment is shown in Table 4. The pair-wise comparison matrix A is built as follows: where n represents the number of evaluation criteria considered, each a ij of the matrix A represents the importance of the ith criterion relative to the jth criterion. If A > 1, then the ithcriterion is more important than the jth criterion, while if a ij < 1, then the ith criterion is less important than the jth criterion. Moreover, the a ij and a ji satisfy the reciprocal constraint. a ij denotes the normalized result by column. The upper half and lower half of the matrix diagonal are reciprocal to each other, therefore matrix A is a reciprocal matrix. Table 4. The importance of the pair-wise comparison judgment [45].

1
Equally preferred 3 Moderately preferred 5 Strongly preferred 7 Very strongly preferred 9 Extremely strongly preferred 2, 4, 6, 8 The adjacent middle value judgment inversion Comparison of factor i to j is b, while the factor j to i comparison scale is 1/b (3) Calculation of the Weight Vector The eigenvalue and the eigenvector are calculated as follows: where W is the criteria weight vector, W is utilized to calculate the average eigenvalue of the pair-wise comparison matrix λ max .

(4) Consistency Inspection of the Judgment Matrix
For the judgment matrix, a consistency method is applied to test the reliability. If the pair-wise comparison matrix deviates from the consistency, the reliability of the results will be lower. To judge the uncertainty, Saaty's consistency index (CI) is used, which is calculated using Equation (5): where n is the number of indices. CR is a measurement of consistency of the pair-wise comparison matrix, and is calculated using Equation (6): where RI is the ratio index. The standard values of RI are shown in Table 5. The CR is acceptable if CR ≤ 0.1; otherwise, we re-evaluate the corresponding weights to avoid inconsistency. Practically, 0.00001 was utilized as substitute for zero in the calculation of CR.

Ecological Carrying Capacity Calculation Based on the Ecological Footprint
The ecological footprint refers to the various types of resources meeting the daily consumption of a certain human population and the area of ecologically productive land necessary for assimilating various types of associated domestic waste [46,47]. The ecological carrying capacity is calculated from the perspective of land supply based on quantification of a group of indices regarding the area of ecologically productive land [48]. The computational equation for the ecological carrying capacity is as follows [49]: where ECC represents the total ecological carrying capacity; N represents the total population of the research area; ecc is the per capita ecological carrying capacity; 12% represents the 12% land area deducted from the ecological supply for protection of the biological productivity; γ is the equivalence factor; y is the yield factor; i is the area of biologically productive land required; and a i is the per capita area of the productive land of i.
The equivalence factor and the yield factor have been given by WWF on a big-scale and are not suitable on the regional scale. We utilized the equivalence factor and the yield factor for Hebei province which are calculated based on NPP [50], because the provincial ecological footprint model based on NPP can reflect the actual land productivity in the region of a province or city. In this way, the ecological footprint can be of practical value in medium-and small-scale regions. The relevant parameters are shown in Table 6.

Comprehensive Evaluation Model Construction
Considering the reality of a coalfield, we adopted the comprehensive indices approach to calculate the final score of environmental quality. This approach not only is based on a clear principle, but also emphasizes the integrality and objectivity of the evaluation. The "integral" refers that the hierarchical structure is integral, each layer has its own calculational criterion and between layers the connection weights are under control by CR. The "objective" refers that the indicators are objective [51]. Based on the four main indicators and the weight sorting supported by HJ/T192-2006, the evaluation model of environmental quality is as follows: where EI is the eco-environmental index, w cli denotes the weight of the ith indicator in criterion layer. V BAI , V VCI , V NGF and V WDI denote the values of the indicators in criterion layer respectively. The average value of EI in the four quarters of each year is used as the final EI of the year to exclude the influence of contingency. During the calculation of the V BAI , V VCI , V NGF and V WDI , the normalization coefficient of each indicator ϕ is utilized to guarantee the score of each indicator is less than or equal to 100: where w il−bai , w il−vci w il−ng f and w il−wdi represent the weights of index layer for BAI, VCI, NGF and WDI respectively. Taking the maximum value among the middle values of each indicator, the normalization coefficient of each evaluation indicator is calculated. The score of each indicator is less than or equal to 100, and the final EI value is also between 0 and 100. The normalization coefficient ϕ BAI , ϕ VCI , ϕ NGF and ϕ WDI are calculated using Equation (13): where ϕ i is the normalization coefficient of each indicator; and I max represents the maximum value among the middle values of each indicator. The I max of each indicator is as follows: where time represents the different period. Finally, the ecological carrying capacity based on the ecological footprint is applied in order to verify the effectiveness of the evaluation model. The flow chart of quantitative evaluation of the ecological environment is shown in Figure 4. = max (17) where represents the different period. Finally, the ecological carrying capacity based on the ecological footprint is applied in order to verify the effectiveness of the evaluation model. The flow chart of quantitative evaluation of the ecological environment is shown in Figure 4.

Index Analysis
The eco-environment evaluation model for the mining area is based on HJ/T192-2006. The weights of criterion layer and index layer are calculated to make the evaluation model achieve the most suitable trade-off using Equations (1)-(6) and the results are shown in Table 7.

Index Analysis
The eco-environment evaluation model for the mining area is based on HJ/T192-2006. The weights of criterion layer and index layer are calculated to make the evaluation model achieve the most suitable trade-off using Equations (1)-(6) and the results are shown in Table 7. Combined with the BAI, VCI, WDI, and NGF, the EI is calculated based on AHP. The scores of the different indices are shown in Figure 5 and Table 8.
Since the integration of the state's coal mining enterprises in Yuxian County in 2008, small-scale coal enterprises have been closed down or merged with large state-owned coal enterprises. The BAI and the VCI of the mining area showed slow growth on the whole ( Figure 5). Vegetation growth was restrained in the year of 2014 in the mining area. Meanwhile, the WDI was also relatively low. The WDI was mainly used to characterize the abundance of water resources in study area, and was thus affected by precipitation and groundwater. The WDI fluctuated between 73-77, and the variation range was relatively small. a negative correlation between EI and energy consumption index. When the energy consumption index becomes higher, the EI will present a downward trend in the next two years, vice versa. The Cuijiazhai mining area and the Shanhou mining area are the part of the Kailuan Mining Bureau. Technological innovation and improvement not only promote economic development, but also reduce energy consumption and has important significance for the restoration and management of the eco-environment.  The correlation coefficients were calculated between the EI and the sub-indices in order to find the driving factors of the EI. The correlation coefficients are shown in Figure 7.  The NGF was influenced by vegetation coverage, DEM, precipitation, etc. Overall, the NGF was mainly influenced by the slope of the study area, because of the mining causing ground collapse and landslides. An undulating terrain formed as a result of these events. The NGF was mainly influenced by DEM and precipitation from 2008 to 2012. In 2014, the NGF was affected by precipitation and vegetation coverage. The BAI and the VCI show similar trends to the EI in Figure 5.
Since the "7.14" coal mine disaster that occurred in Yuxian in 2010, which claimed the lives of 35 coal mine workers, China has begun to integrate coal mining enterprises. Shutting down small-scale coal mining enterprises, restructuring, and technological innovation have been put into effect. Figure 6 shows the output energy consumption of Yuxian before rectification and after rectification. The output energy consumption of each mineral company, especially the Kailuan Group, decreased with time. Energy consumption index represents that how much energy would be consumed for unit production value. The bigger value the index represents, the more coal would be consumed, which bring on more environment pollution. Based on Figures 5 and 6, it could be found that there is a negative correlation between EI and energy consumption index. When the energy consumption index becomes higher, the EI will present a downward trend in the next two years, vice versa. The Cuijiazhai mining area and the Shanhou mining area are the part of the Kailuan Mining Bureau. Technological innovation and improvement not only promote economic development, but also reduce energy consumption and has important significance for the restoration and management of the eco-environment.  The correlation coefficients were calculated between the EI and the sub-indices in order to find the driving factors of the EI. The correlation coefficients are shown in Figure 7. The R 2 values between the EI and BAI/VCI are 0.9951 and 0.9911, respectively, which indicates that the BAI and the VCI are both driving factors of the EI. However, the BAI and the VCI were calculated by land cover in the study area, and there is sufficient evidence to say that environmental quality is mainly determined by land-cover types in the Yuxian coalfield. The next decisive factor is NGF, with R2 equal to 0.7060. The correlation between the EI and the WDI is not significant, with R 2 equal 0.1729, and the reason for this is the precipitation. The ground runoff fluctuates with the change of precipitation. The precipitation does not have a fixed regular pattern in the warm temperate continental monsoon climate, so the ground runoff is also irregular.

Environmental Quality Analysis and Validation
The EI has shown an increasing trend, as a whole, and this phenomenon has occurred since 2008. The R 2 values between the EI and BAI/VCI are 0.9951 and 0.9911, respectively, which indicates that the BAI and the VCI are both driving factors of the EI. However, the BAI and the VCI were calculated by land cover in the study area, and there is sufficient evidence to say that environmental quality is mainly determined by land-cover types in the Yuxian coalfield. The next decisive factor is NGF, with R2 equal to 0.7060. The correlation between the EI and the WDI is not significant, with R 2 equal 0.1729, and the reason for this is the precipitation. The ground runoff fluctuates with the change of precipitation. The precipitation does not have a fixed regular pattern in the warm temperate continental monsoon climate, so the ground runoff is also irregular.

Environmental Quality Analysis and Validation
The EI has shown an increasing trend, as a whole, and this phenomenon has occurred since 2008. However, there was an exception, the EI declined in 2014. The reason for this may be the growth of vegetation was blocked in 2014. Wang et al. [52] analyzed the temporal and spatial characteristics of the NPP of the Yuxian coalfield between 2013 and 2015, based on high-resolution remote sensing data, and their results showed that the NPP value was the lowest in 2014.
The final results are divided into different levels to evaluate the quality of the eco-environment. The five levels of eco-environmental quality and four levels of change degree are divided based on HJ/T192-2006. The details are shown in Tables 9 and 10. According to the rank of the eco-environmental quality, the environmental condition of the study area is 'good'. Table 9. Grade of eco-environmental quality [53].  There is no clear change in the ecological environment If 2 < ∆EI ≤ 5, the eco-environment is slightly better, if −2 > ∆EI ≥ −5, the eco-environment is slightly worse.
On the other hand, the ecological carrying capacity calculated based on the ecological footprint was used to verify the result of the evaluation model. Based on the land utilization in the mining area, we obtained data such as the area of cultivated land, area of grassland, area of building land, area of productive land, area of water, and area of bare land, from 2005 to 2016, which is shown in Table 11. These data were divided by the total population of the study area and then multiplied by the corresponding balance factors and yield factors. The computational results were accumulated to obtain the per capita ecological carrying capacity of the coalfield in the study area. Based on the suggestions for computation of the ecological carrying capacity made by the World Commission on Environment and Development [54], 12% should be deducted from the ecological carrying capacity as the productive area for biological protection.
The ecological carrying capacity, which is shown in Table 12, Figures 8 and 9, shows a similar trend to the EI. Because of the ground collapse caused by the mining, a large amount of land was fully utilized. The reason for the low ecological carrying capacity is the large population and the mining activities in the area. The per capita ecological carrying capacity increased from 0.0413 (hm 2 /cap) to 0.0604 (hm 2 /cap) between 2005 and 2016, after deducting 12% of the productive area for biological protection. This demonstrates that the ecological carrying capacity is low. The EI increased from 66.62 to 70.71 between 2015 and 2016, which was relatively large. It is shown in Table 12 that the eco-environment quality was slightly improved between 2015 and 2016. According to the reports of government work in Yuxian County, the total grain output of Yuxian County was 55,624 ton in 2010 and it almost tripled in 2016, reaching 145,800 ton. Besides, since Zhangjiakou City, which is the superior administrative agency of Yuxian County, was elected as the host city of XXIV Olympic Winter Games in July 2015, a series of greening actions began to be implemented and seventy-five square kilometers of fine greening project was completed in 2016. This can also serve as evidence that the improvement of EI in 2016 was relatively large. obtain the per capita ecological carrying capacity of the coalfield in the study area. Based on the suggestions for computation of the ecological carrying capacity made by the World Commission on Environment and Development [54], 12% should be deducted from the ecological carrying capacity as the productive area for biological protection. The ecological carrying capacity, which is shown in Table 12, Figures 8 and 9, shows a similar trend to the EI. Because of the ground collapse caused by the mining, a large amount of land was fully utilized. The reason for the low ecological carrying capacity is the large population and the mining activities in the area. The per capita ecological carrying capacity increased from 0.0413 (hm 2 /cap) to 0.0604 (hm 2 /cap) between 2005 and 2016, after deducting 12% of the productive area for biological protection. This demonstrates that the ecological carrying capacity is low. The EI increased from 66.62 to 70.71 between 2015 and 2016, which was relatively large. It is shown in Table 12 that the eco-environment quality was slightly improved between 2015 and 2016. According to the reports of government work in Yuxian County, the total grain output of Yuxian County was 55,624 ton in 2010 and it almost tripled in 2016, reaching 145,800 ton. Besides, since Zhangjiakou City, which is the superior administrative agency of Yuxian County, was elected as the host city of XXIV Olympic Winter Games in July 2015, a series of greening actions began to be implemented and seventy-five square kilometers of fine greening project was completed in 2016. This can also serve as evidence that the improvement of EI in 2016 was relatively large.
The correlation coefficients between ecological carrying capacity and the EI were calculated in order to verify the effectiveness of the evaluation model. The fitting curve of the EI and ECC, with high R 2 values, indicates that the EI and ECC are closely correlated. This phenomenon also illustrates that the results of the novel evaluation model based on multi-temporal remote sensing images and auxiliary data are acceptable.
Although environmental quality has shown an improving trend from 2005 to 2016, the protection and restoration of the eco-environment in the Yuxian coalfield is still an urgent problem that needs to be solved. The eco-environment of the coalfield had been seriously affected by the largescale and heavy exploitation of the coal resources during the past decades. This damage has a longterm influence and is difficult to recover. The EI has become better in recently years because of the introduction of artificial greening project and the adequate precipitation. Tiny fractures in rock that can store water were collapsed during periods of coal mining. And then the surface vegetation growth would rely on the external water sources and human disturbance. When WDI dropped in 2014, EI also declined. The aquifer restoration requires long-term efforts. In recent years, the government and relevant departments have become aware of the importance of ecological restoration of coal mining areas, and some measures have been taken to control and improve the situation. However, if these problems are not paid enough attention, environmental pressure will accumulate because of the exploitation of the coal resources. Moreover, species diversity will decrease and the human living environment may be threatened in the mining area. Ultimately, the eco-environment will not be able to adapt to the development of the economy and society.

Conclusions
This paper provides a quantitative evaluation of the eco-environmental quality of the Yuxian coalfield, based on land-cover and statistical data from 2005 to 2016, using remote sensing and a multi-criteria decision-making technique. The land cover was utilized as the basic indicator for ecoenvironmental analysis of mining area with a small number of auxiliary indicators. What's more, the weights are calculated to make the evaluation model achieve the most suitable trade-off by AHP. The mine-areas-oriented model is more suitable than the recommended indicators and weights given by HJ/T 192-2006. To verify the validity of our model, the ecological carrying capacity based on the ecological footprint is introduced. The eco-environmental quality of the Yuxian coalfield is not in an optimal state of harmonious development, based on the computed results. Although the EI has shown a slowly increasing trend, the value of the EI is low, as a whole, and the level of ecoenvironmental quality is only 'good'. The BAI and the VCI are driving factors which have a good correlation with the EI. It can be seen from Table 9 that the per capita ecological carrying capacity increased to 0.0604 (hm 2 /cap) in 2016 from 0.0413 (hm 2 /cap) in 2005, which is consistent with the EI, and the R 2 is 0.9734. By combining features extracted from remoting sensing and auxiliary data, the comprehensive evaluation model of the mining area in Yuxian County was established, achieving a quantitative evaluation of the eco-environment of the mining area of Yuxian County.
In conclusion, the research results could provide a basis for the management of the ecoenvironment of the Yuxian mining area, and a reference and scientific basis for achieving sustainable development of the eco-environment in the coalfield. The evaluation method is a supplement to the time-series quantitative evaluation of the eco-environment, and also helps us to explore the ecoenvironment in the mining area.
Author Contributions: The research was designed and performed by X.W. and K.T. The data was collected and analyzed by Y.C., K.X. and J.D. The paper was written by X.W. and K.T., K.X., and finally checked and revised by X.W., K.T. K.X., Y.C. and J.D. All authors read and approved the final manuscript. The correlation coefficients between ecological carrying capacity and the EI were calculated in order to verify the effectiveness of the evaluation model. The fitting curve of the EI and ECC, with high R 2 values, indicates that the EI and ECC are closely correlated. This phenomenon also illustrates that the results of the novel evaluation model based on multi-temporal remote sensing images and auxiliary data are acceptable.
Although environmental quality has shown an improving trend from 2005 to 2016, the protection and restoration of the eco-environment in the Yuxian coalfield is still an urgent problem that needs to be solved. The eco-environment of the coalfield had been seriously affected by the large-scale and heavy exploitation of the coal resources during the past decades. This damage has a long-term influence and is difficult to recover. The EI has become better in recently years because of the introduction of artificial greening project and the adequate precipitation. Tiny fractures in rock that can store water were collapsed during periods of coal mining. And then the surface vegetation growth would rely on the external water sources and human disturbance. When WDI dropped in 2014, EI also declined. The aquifer restoration requires long-term efforts. In recent years, the government and relevant departments have become aware of the importance of ecological restoration of coal mining areas, and some measures have been taken to control and improve the situation. However, if these problems are not paid enough attention, environmental pressure will accumulate because of the exploitation of the coal resources. Moreover, species diversity will decrease and the human living environment may be threatened in the mining area. Ultimately, the eco-environment will not be able to adapt to the development of the economy and society.

Conclusions
This paper provides a quantitative evaluation of the eco-environmental quality of the Yuxian coalfield, based on land-cover and statistical data from 2005 to 2016, using remote sensing and a multi-criteria decision-making technique. The land cover was utilized as the basic indicator for eco-environmental analysis of mining area with a small number of auxiliary indicators. What's more, the weights are calculated to make the evaluation model achieve the most suitable trade-off by AHP. The mine-areas-oriented model is more suitable than the recommended indicators and weights given by HJ/T 192-2006. To verify the validity of our model, the ecological carrying capacity based on the ecological footprint is introduced. The eco-environmental quality of the Yuxian coalfield is not in an optimal state of harmonious development, based on the computed results. Although the EI has shown a slowly increasing trend, the value of the EI is low, as a whole, and the level of eco-environmental quality is only 'good'. The BAI and the VCI are driving factors which have a good correlation with the EI. It can be seen from Table 9 that the per capita ecological carrying capacity increased to 0.0604 (hm 2 /cap) in 2016 from 0.0413 (hm 2 /cap) in 2005, which is consistent with the EI, and the R 2 is 0.9734. By combining features extracted from remoting sensing and auxiliary data, the comprehensive evaluation model of the mining area in Yuxian County was established, achieving a quantitative evaluation of the eco-environment of the mining area of Yuxian County.
In conclusion, the research results could provide a basis for the management of the eco-environment of the Yuxian mining area, and a reference and scientific basis for achieving sustainable development of the eco-environment in the coalfield. The evaluation method is a supplement to the time-series quantitative evaluation of the eco-environment, and also helps us to explore the eco-environment in the mining area.
Author Contributions: The research was designed and performed by X.W. and K.T. The data was collected and analyzed by Y.C., K.X. and J.D. The paper was written by X.W. and K.T., K.X., and finally checked and revised by X.W., K.T. K.X., Y.C. and J.D. All authors read and approved the final manuscript.