Mapping Annual Land Disturbance and Reclamation in a Surface Coal Mining Region Using Google Earth Engine and the LandTrendr Algorithm: A Case Study of the Shengli Coalfield in Inner Mongolia, China

: The development and utilization of mining resources are basic requirements for social and economic development. Both open-pit mining and underground mining have impacts on land, ecology, and the environment. Of these, open-pit mining is considered to have the greatest impact due to the drastic changes wrought on the original landform and the disturbance to vegetation. As awareness of environmental protection has grown, land reclamation has been included in the mining process. In this study, we used the Shengli Coalfield in the eastern steppe region of Inner Mongolia to demonstrate a mining and reclamation monitoring process. We combined the Google Earth Engine platform with time series Landsat images and the LandTrendr algorithm to identify and monitor mining disturbances to grassland and land reclamation in open-pit mining areas of the coalfield between 2003 and 2019. Pixel-based trajectories were used to reconstruct the temporal evolution of vegetation, and sequential Landsat archive data were used to achieve accurate measures of disturbances to vegetation. The results show that: (1) the proposed method can be used to determine the years in which vegetation disturbance and recovery occurred with accuracies of 86.53% and 78.57%, respectively; (2) mining in the Shengli mining area resulted in the conversion of 89.98 km 2 of land from grassland, water, etc., to barren earth, and only 23.54 km 2 was reclaimed, for a reclamation rate of 26.16%; and (3) the method proposed in this paper can achieve fast, efficient identification of surface mining land disturbances and reclamation, and has the potential to be applied to other similar areas.


Introduction
Rapid economic and social development have driven increasing demand for mineral resources. However, as mining activities occupy and also damage land resources, the exploitation of these resources often has a significant impact on the ecological environment [1], leading to various secondary environmental problems such as land loss, land occupation, vegetation degradation, and soil erosion [2]. Therefore, mining and urbanization are currently regarded as the economic activities that have the most destructive impacts on natural ecosystems. These activities have various negative environmental consequences, and the governments of many countries require that areas degraded by mining be recovered. Open-pit mining, ore mining, and placement of soil cover (overburden dumps) following strip mining can cause changes in land cover and land use. Land-cover changes around large surface mines have already caused serious environmental landscape degradation in parts of the USA [3]. Sustainable mining and green mining both require continuous monitoring of these changes to identify the long-term impacts on the environment and land cover. Monitoring such changes can then provide basic security measures as well as basic data for use in formulating land reclamation and ecological restoration strategies. In addition, quantifying the spatiotemporal changes in an area as a result of surface mining is particularly important for understanding the effects of mining activities and for assessing their ecological, environmental, and socioeconomic impacts [4].
Disturbances to land and vegetation caused by open-pit mining are key concerns for researchers and land management agencies. There have been numerous studies of the changes caused by mining at the micro-level, including changes to the physical and chemical properties of land [5], vegetation disturbance [6], and ground subsidence [7]. Before the advent of widespread remote-sensing techniques, field investigation had been an effective means of obtaining information about groundlevel damage and surface monitoring prior to reclamation. However, due to the historical legacy of past mining, field surveys and studies cannot determine quantitative changes, and are rather slow and expensive.
Advances in satellite-to-ground observation technology have made remote sensing an effective method of measuring the nature and extent of vegetation [8]. The earliest application of remotesensing data was in comparisons between data for different years. Satellite remote-sensing detection can be performed visually or by using computer-aided interpretation via automatic pattern recognition technology [9,10]. Remote-sensing technology can monitor and evaluate vulnerable landscapes, natural resources, ecosystems, and biodiversity conservation areas cost-effectively. In recent decades, extensive studies have been conducted on the use of low-and medium-resolution optical satellite data, including Landsat and MODerate-resolution Imaging Spectroradiometer (MODIS), to monitor the impacts of mining activities on ecosystems [11][12][13]. In addition, radar satellite data have been used to observe the deformations resulting from surface subsidence in mining areas [14].
The methods for monitoring mining disturbance and reclamation activities include various remote sensing classification methods. Several studies have measured coal mining surface disturbance using vegetation indices, such as the Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), etc., which can typically be applied to monitoring via remote sensing. For example, a time series NDVI has been used to monitor mining-related grassland vegetation disturbances [15,16]. Jia et al. [17] and Li et al. [18] both studied the land damage and restoration processes associated with grassland open-pit and rare-earth mining using the Single Spectrum Analysis (SSA)-Mann Kendall and CART classification on multi-source time series NDVI, respectively. Multivariate change detection has been used to measure the expansion and spread of gold mining on the ground during two time periods [19]. Machine learning algorithms such as support vector machines and random forests have also been used to study land use classification in mining areas [20].
Monitoring studies have evolved from the use of several images to yearly time series of images. More recent studies have focused on using data with various levels of resolution to analyze and identify the spatiotemporal process of vegetation disturbance [21] and evaluate the effects of land reclamation [22,23]. This type of research is currently focused on monitoring the changes in time series remote-sensing data. For example, Townsend et al. [24] used the Landsat time series between 1976 and 2006 to detect changes in Central Appalachian surface mining and reclamation areas. Time series remote-sensing data can also be used to analyze the impact of surface mining on cultivated land and food security [25]. However, this method of manual interpretation and classification using multi-period images produces a heavy workload and readily accumulates classification errors.
Another problem is that the yearly image from a particular time node in a specific year can often only represent the vegetation and land use on the surface in a transient state. It cannot recognize characteristics that might change over a period of time. This makes monitoring events with nonlinear or highly spatiotemporal heterogeneity difficult. Often, annual data are available only for different time points in each year due to factors such as satellite transit time and weather conditions, including cloud cover [26]. In many studies in which satellite images for the vegetation growth period have been selected, these data are usually for a specific time period between May and September [27]. This, however, challenges the credibility and comparability of the monitoring results.
In recent years, the development of cloud computing technology has made possible the use of available, massive remote-sensing satellite transit data for detecting changes in fine-grained, hightemporal-resolution data. The Google Earth Engine (GEE) launched by the American multinational technology company Google LLC has collected common remote-sensing data sets such as MODIS, Landsat, and Sentinel [28], and can use online or offline programming to obtain and process shared data. Cloud computing analyzes and processes remote-sensing data, thereby avoiding the tedious processes of data downloading and preprocessing required by traditional remote-sensing analysis models. The GEE platform offers strong data analysis and processing capabilities, and is gradually being applied to areas such as urban expansion [29], dynamic monitoring of ecological quality [30,31], and the monitoring of vegetation coverage [32].
Related algorithm-based technologies for time series remote-sensing change monitoring have also led to rapid progress, and the use of time series remote-sensing data analysis has increased in recent years. Current vegetation disturbance detection methods are now quite mature, with many excellent algorithms [33]. They are used to monitor changes in land cover by analyzing images as temporal sequences. These time series algorithms have enabled time series fitting based on the spectral index of vegetation, using methods such as the best index slope extraction (BISE) algorithm, time series decomposition methods such as wavelet changes, breakpoint detection algorithms based on time series segmentation such as the vegetation change tracker (VCT) model [34] and the LandTrendr algorithms [35], and multiclass times series curve detection algorithms such as the continuous change detection and classification (CCDC) [36].
Based on the premise that vegetation index patterns change over time along particular trajectories, mines can be differentiated from undisturbed natural areas because the trajectories of reclaimed mines differentiate them from urban disturbances [37]. Long-term sequence change detection methods are classified into four main types: spectral variables, image classification, spectral trajectory-based analysis, and data fusion methods [38]. The currently popular Landsat-based time series algorithms come mainly from North America. These include the Landsat interference and trend monitoring algorithm LandTrendr (for Landsat-based detection of trends in disturbance and recovery), which uses spectral trajectories developed by Oregon State University and the United States Forest Service, first proposed by Kennedy et al. [35]. It is used primarily in forest change detection, but has gradually been extended to areas of water and to cities. There are also other fields of application, such as the use of water probability to generate surface water change maps [39]. The advantages of LandTrendr are that it relies on Landsat archive data that can be accessed for free, and it is simple to operate. It can provide useful data for environmental assessment and remediation, including raster maps of bands that represent change trends, identification of points in time when abrupt changes occur, and analysis of various changing pixel attributes to promote better environmental assessment. The LandTrendr algorithm has been used to detect and analyze the timing of damage caused by surface mining [40].
In this study, in order to identify the spatiotemporal processes of surface mining and land reclamation more efficiently and accurately, an attempt was made to synthesize year-round images for the study period and the study area using the GEE cloud computing platform and the LandTrendr algorithm. The process can be summarized as data preprocessing, NDVI extraction, LandTrendr processing of the yearly trajectories, and mapping of the years of disturbance and reclamation. The Landsat series of remote-sensing data were used to generate an annual NDVI index that characterizes the annual changes in pixels by constructing long-term annual pixel trajectories.
Timely reclamation following disturbances to land and vegetation is the key to ensuring green mining and sustainable mining [41]. At the national level, it is also hoped that mining and reclamation can be conducted concurrently [42]. However, most mining companies are principally concerned about the efficiency of their mining operations, and so neglect land reclamation and ecological restoration. Therefore, accurate and timely monitoring of areas of damaged and reclaimed land is particularly important. The method described in this paper achieves rapid and efficient extraction of the spatiotemporal processes of open-pit coal mining and reclamation, thus providing basic data for safe and sustainable mining production and land reclamation. The results of this study may also come to play a crucial role in monitoring and evaluating the implementation of land reclamation.

Study Area
Open-pit mining is a driving force for changes in land cover, leading to the loss of natural vegetation and topsoil as well as disruption of ecosystem service flows. In the 21st century, mining has gradually decreased in North America and Europe, but has expanded in South America (including Chile and Brazil) and Asia (including India and China). In the latter regions, the proportion of open-pit mining sites has increased significantly [43]. Surface mining is mining on the ground in which the soil and rock over the mineral deposits are removed. This contrasts with underground mining, in which the overlying rock is left in place and the mineral is removed via shafts or tunnels. The area chosen for this study lies on the eastern plateau of Inner Mongolia, where surface mining has severely disturbed the grassland vegetation. In open-pit mining, the process of disturbing the surface includes stripping of overlying strata, coal-seam mining, solid backfilling, return of top soil, and reclamation.
The study area consisted of the Shengli Coalfield, which is located in Zhenglanqi, to the northwest of Xilingol City, in the Inner Mongolia Autonomous Region (43°54'15"-44°13'52"N, 115°24'26"-116°26'30"E) of China ( Figure 1). The terrain of the study area is a gentle undulating piedmont plain, while the climate is that of a semi-arid steppe, with an average temperature of 1.7 °C, a maximum temperature of 38.3 °C, and a minimum temperature of −42.4 °C. The rainfall there averages 294.7 mm a year; the maximum rainfall is 481.0 mm, the minimum rainfall is 146.7 mm, and the average annual evaporation is 1794.6 mm. The area is usually windy in the spring, with wind speeds of 2.1-8.4 m·s −1 and an average annual wind speed of 3.5 m·s −1 . The soil types there are chestnut calcareous soil, meadow calcareous soil, and meadow soil. Desertified, gravelly calcareous soil has been formed due to grassland degradation. The soil organic matter content is low and the soil fertility is poor. There are two coal-bearing sections in the Shengli Coalfield: the lower coal section of the Cretaceous Bayanhua Group Xilin Formation and the Shengli Formation coal-bearing section. There are seven coal mining groups, with an average burial depth of 200-500 m and an average thickness of 8-60 m. In 1989, the proven and retained reserves were 15.932 billion tons and 15.931 billion tons, respectively.
Under the temperate continental climate of the region, this vulnerable natural ecosystem is typical of temperate semi-humid and semi-arid steppes. The area is an important source of green agricultural and livestock products in China. There are 13,786 km 2 of high-quality natural grasslands available in the region. The average vegetation growth season is 141 days, which includes the period from May through September [44]. Vegetation coverage is usually maintained at 40% throughout the year [45].

Overview of the Methodology
Our methodology is based on the premise that open-pit mining and reclamation activity will cause abrupt changes to vegetation, which can be shown by the trajectory of the vegetation indices. As shown in Figure 2a, there are three main vegetation time series trajectories of pixels in the study area, which represent areas that are undisturbed, disturbed and reclaimed, and disturbed without reclamation. Time series analysis was used to find breakpoints in the trajectories, thereby indicating sudden changes: surface mining including strip mining, open-pit mining, and mountaintop removal mining. In the surface mining process, soil and rock overlying the mineral deposits (the overburden) are removed, and the vegetation is disturbed in this process. During overburden removal, substantial vegetation is disturbed by both stripping and dumping (second curve in Figure 2a). The dumps are subsequently reclaimed and revegetated, and the open-pit is covered by topsoil and revegetated following mining activities (third curve in Figure 2a). The aforementioned process includes vegetation disturbance and reclamation. According to our field study, natural recovery usually takes more than three years; thus, these sudden changes within two years can be considered to be caused by human activities.
The LandTrendr algorithm was used to find breakpoints; a simple concept map based on the work of Kennedy et al. [35] is shown in Figure 2b. This process of line segmentation was guided by the curve itself, and then the actual breakpoints among the potential breakpoints were decided according to the parameter settings. The entire process is shown in Figure 3 and has two main components: (i) disturbance to the vegetation of the surface mining area and detection of reclamation and (ii) accuracy verification. First, Landsat surface reflection data were preprocessed by masking snow and cloud cover. Then the 2003-2019 sequence image collection was processed to form the annual 95th percentile NDVI. The selection of Landsat images for the entire period avoided manual selection and shadow accumulation errors. The annual 95th percentile NDVI was then used to reduce vegetation index extraction errors caused by the selection of specific time nodes and to facilitate the analysis of this vegetation disturbance and reclamation process in a more efficient and time-saving manner. Finally, spatiotemporal information about sudden changes in the NDVI time series was gained from using the LandTrendr algorithm. In the accuracy verification, we picked sample points from areas where changes were detected and validated them with the images obtained from Google Earth.

Data Acquisition and Preprocessing
This study selected all of the Landsat surface reflectance (SR) images on the GEE platform from 1 January 2001 through 30 December 2019, including images from Landsat 5, 7, and 8. The Landsat SR data from the GEE platform had already had the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) atmosphere correction applied. The number of SR images available for the study area for the period from 2001 through 2019 is shown in Figure 4a. Using the quality assessment (QA) band to indicate the possibilities of cloud, cloud shadow, snow, etc., in each pixel, the cloud, cloud shadow, and snow pixels were masked. The number of remaining valid images in each pixel of the study area is shown in Figure 4b. In this way, the synthesis of good observations produced a collection of cloud-free SR images for each year [46] for subsequent processing. The collection of good observations was then ready to generate a collection of 20 images (2003-2019) with a specified index indicating changes that required observation. A previous study by Karan showed that NDVI was the best indicator for detecting changes in coal mine vegetation [47]. Thus, we used the NDVI to represent vegetation growth. The NDVI was calculated as follows: where and are the surface reflection values in the red and infrared bands, respectively. The annual phenology curves of vegetated pixels differ from year to year, but some characteristics of the curve, such as peaks, can indicate the growth of vegetation. We then used the 95th percentile NDVI to represent the growth of vegetation for that year. Finally, the 2001-2019 image collection comprised of 20 images was ready for time series fitting analysis using the LandTrendr algorithm.

Identification of Disturbance and Reclamation Years
Mining activity includes land cover removal, overburden dumps, and construction that causes the natural environment to abruptly change to barren land or a structure [48]. Various reclamation strategies can be used to seek recovery of vegetation and restoration of ecological function [49]. These reclamation processes usually require months or even years, while mining itself results in relatively fast, severe losses of vegetation. Reclamation leads to faster and greater vegetation recovery than natural processes driven by climate, and is thus suitable for measurement via interannual vegetation changes. As shown in Figure 2, the mining disturbance date of a pixel occurred when an abrupt minimum NDVI was observed, and the reclamation date occurred when vegetation recovery suddenly increased and stabilized.
The LandTrendr algorithm, with its GEE version available online according to Kennedy et al. [50], uses an iterative calculation that includes continuous determination of breakpoints and fitting based on the time series trajectory of a single 30 m × 30 m Landsat image pixel. Breaks are set according to particular parameter values, such as identified mutation types, the maximum number of breakpoints, etc. Different parameters did not affect the result; thus, default parameters were used to simplify the workflow. Each pixel in the resulting image had a matrix indicating whether there had been an abrupt change and, if so, when it changed. By choosing vegetation loss, information about the spatiotemporal distribution of vegetation disturbance occurring in the mining area was gained. Pixels for a sudden increase in vegetation were also obtained. These points were determined and are shown on the map, along with the time of the abrupt change. The mining and reclamation years were determined based on the assumptions that vegetation loss represented mining and that vegetation gain represented reclamation. Thus, maps of vegetation damage and vegetation restoration caused by mining and reclamation, respectively, were generated.

Data Post-processing, Accuracy Assessment, and Validation
The results of our study indicate that disturbances occurred in the period from 2004-2019 and reclamation occurred in the period from 2007-2019. Thus, we randomly selected 50 points in the disturbance area and 30 points in the reclamation area for each year. This produced a total of 800 disturbance and 390 reclamation sample points. Then, we used the high-resolution image data on Google Earth to perform interactive visual calibrations and determine the disturbance and recovery years of each sample point. By comparing the sample label with the algorithmic identification results, the user, producer, overall accuracy, and kappa coefficients for mining disturbance and reclamation monitoring were calculated, as described in the next section.

Disturbance and Reclamation Accuracy Assessment
There were 26,153 identified reclamation pixels, of which 25,905 (99.05%) were disturbed. The overall accuracies of vegetation disturbance and reclamation date (year) identification were 0.84 and 0.81, and the kappa coefficients were 0.83 and 0.79, respectively ( Figure 5). Thus, the overall disturbance and reclamation recognition accuracies were found to be quite good.

Spatiotemporal Characteristics of Vegetation Changes
The spatiotemporal mapping of the mining disturbance and reclamation years is shown in Figures 6 and 7. There are three main mining areas in the Shengli Coalfield, which are circled in red in the figures. Figures 6 and 7 each include a year map of human activity (Figures 6a and 7a), an enlarged map of the three main mining areas (Figures 6b and 7b), and a time series of a region representing the typical tracks of mining and reclamation, respectively, with high-resolution Google Earth images below the trajectory (Figures 6c and 7c).
According to the results of our study, large-scale mining began after 2007, while most recovery processes occurred after 2011. Reclamation and mining activities both show spatial patterns of starting on one side of the mine and gradually advancing.   Figure 8a is a histogram representing the area disturbed and reclaimed each year in the study area. It is worth noting that mining activity has decreased since 2014 and reclamation activity has become more concentrated since 2011. In Figure 8b, the cumulative mining and reclamation areas are shown in the cumulative area chart, and the yearly reclamation rate is drawn as a line chart, with a specific number noted above the curve. The area of cumulative mining minus that of reclamation reaches a dynamic balance after 2010 due to reclamation commencing. A significant rise in the reclamation rate during 2010-2011 is evident, and the reclamation rate stabilizes at 25% after 2014. We further investigated the mining and reclamation time of each pixel and counted the number of pixels in different time periods. According to the statistics presented in Table 1, in which the reclamation year is given horizontally and the vegetation disturbance year (caused by mining activity) is given vertically, 99,982 pixels were disturbed between 2003 and 2019, accounting for approximately 90 km 2 . During the same period, 26,153 pixels were monitored for reclamation, accounting for approximately 24 km 2 . Of these, 25,948 pixels were previously detected as disturbed, with a drastic decline in vegetation coverage, while 205 pixels were not. The latter group may be reclamation pixels that did not experience any severe disturbance. It is evident that the earlier the mining started, the higher the proportion of land reclamation. For example, 90% of the area subjected to vegetation loss caused by mining activities in 2004 has been recovered, while less than 10% of the area disturbed after 2012 has been repaired.  Note: UP indicates the number of unchanged pixels, i.e., those that underwent mining without reclamation, while RM indicates the percentage of pixels that were reclaimed in a later period after being mined in a particular year.
The information in Table 1 can be visualized using Figure 9, which shows a pixel reclamation matrix for various years. The following characteristics are most noticeable: those pixels disturbed during the period from 2007-2011 have been reclaimed in large numbers in the period from 2011-2014. The reclamation chart of the region forms a triangle in the picture. In addition, reclamation performed in 2011 covers mining from 2004 through 2008. Finally, it can be seen that the pixels mined in 2004 have been reclaimed for many years until 2019.

Discussion
The grasslands of Inner Mongolia are rich in coal. Their delicate ecosystems have been threatened by increased mining activity as a result of the depletion of the coal resources of other provinces following years of mining [51]. The developing Shengli Coalfield in Xilingol was chosen as the study area for examination of yearly mining and reclamation activities, an area that has been the subject of numerous studies, varying from the impacts of mining and reclamation on vegetation dynamics to their impacts on the landscape [52]. Utilizing the GEE platform and the LandTrendr algorithm, our methodology was used to determine the spatiotemporal distributions of mining and reclamation rapidly, accurately, and with high spatial resolution (30 m). One benefit of applying this methodology to open-pit mining is that it is an unsupervised process, unlike time series filtering and classification, which require sample training [53] and integrated clustering of pixel classification [54]. It is also worth noting that the use of SR data and operations on image collections enable one to examine mining and reclamation situations where grassland is the original vegetation. Other disturbance and reclamation detection projects have used similar methods in forested areas [40,55].
If the year map on which pixels detect disturbances and reclamation appears to be clustered rather than discrete, this indicates that there are massive mining or reclamation activities. This method can be applied to the open-pit coal mines in the steppe area and can detect the mining and reclamation years over long periods of time.

Applicability of the Method
The mining of open pits is often accompanied by the destruction of vegetation, which causes a sudden reduction in the vegetation index. Therefore, time series change monitoring and the LandTrendr algorithm can be used to determine the year when the quantity of vegetation changed abruptly. The NDVI is more robust than other vegetation indices, such as the EVI in grasslands, and the 95th percentile of the annual NDVI can reflect vegetation growth in that year quite well [56]. It can be used broadly in areas that exhibited vegetation coverage before open-pit mining began. This algorithm is based on the GEE web platform, which has almost no computing requirements and is easy to transfer.

Comparison with Existing Products
Using the information on vegetation disturbance years presented in this paper, it is possible to generate a map of vegetation changes for any time period in the study area. In addition, comparing multiphase images is a traditional method for detecting vegetation changes in open-pit mines. We used the National Land Cover Database (NLCD) data set to generate vegetation change maps from 2010 through 2018 [57]. It should be noted that NLCD is a product that is interpreted manually and cannot be used as a factual basis. For direct comparison, in Figure 6c, we chose a rectangular area of 4 km × 5 km. The results show that the area of reduced vegetation from NLCD products is greater than the results presented in this paper. Moreover, the image spots of the NLCD products are more concentrated, mainly because NLCD is an object-oriented manual interpretation product and the algorithm used in this study is pixel-based. Thus, it is difficult to avoid the salt-and-pepper noise phenomenon [58]. However, our results detected more details: the vegetation in the black ellipse in Figure 10 has not been significantly reduced, but NLCD mistakenly identified it as vegetation. In addition, the blue ellipse in Figure 10 shows that the vegetation disturbance boundary identified by results in this study is more similar to the original image of the boundary of topsoil stripping due to mining. For this reason, the method of detecting vegetation disturbance put forward in this paper is superior to the traditional comparative method.

Method Defects and Generalization
The defects in and prospects for generalization of this methodology are discussed primarily from three angles: applicable objectives, algorithmic operability, and accuracy. This method can be applied to detect mining and reclamation activity in open-pit mines. As for underground mining, subsidence caused by mining does not usually disturb the vegetation, and thus is hard to monitor from the surface. With regard to algorithmic operability, adjusting the algorithm in the GEE platform is relatively simple. However, pixel-based algorithms usually have accuracy problems. As shown in Table 1, there are 41 pixels that were incorrectly identified, with their reclamation year earlier than their mining year. A few pixels outside the main mining area were detected as disturbed or reclaimed. We can erase such pixels via eight-field smoothing. In addition, there are some undetected pixels. These types of errors are common in pixel-based algorithms due to salt-and-pepper noise and other causes. Moreover, we used this method in a specific mining region, which means that the mining boundary was identified prior to analysis. Thus, how to implement the proposed method for an area with an unknown mining boundary and avoid the impacts resulting from urbanization and related vegetation disturbance requires more work to be done.

Implications for Ecological Monitoring and Reclamation
The methodology that we have pioneered can identify effectively areas newly subject to reclamation and mining, as well as the net area that is destroyed (the cumulative reclamation area minus the cumulative mining area) each year. Figure 8b shows that the areas of cumulative damage by mining changed little after 2010, as the area already disturbed but not reclaimed tended to vary slowly, with few fluctuations. This reflects the effectiveness of reclamation and the dynamic balance that can be maintained by conducting reclamation while continuing mining activities. This serves to protect the environment and maintain the original ecosystem functions. The results extracted by the algorithm indicate clear monitoring of mining disturbance and reclamation. This algorithm can play a role in monitoring vast areas for mine reclamation in the future. At the enterprise level, the use of this algorithm for rapid and extensive nonmanual monitoring is low-cost, and can be combined with direct observation of land reclaimed from mining. It may enable a quantitative understanding of the effects of mining and reclamation, and allow the exploration of more efficient reclamation methods. Nevertheless, further specific analyses need to be performed in order to explore the mechanisms and impacts of mining activities, such as comparisons with the data for annual coal production [59] and various analyses of vegetation dynamics that address the impacts of mining on the health of ecosystems [60].

Conclusions
China's coal consumption is substantial, and the output of coal mines is indispensable to economic development. Pressure from environmental and energy needs, however, means that coal mines must be transformed. Open-pit mines produce extensive surface damage, and application of soil reconstruction and vegetation reclamation is especially important. To this end, we proposed a method of rapid, largescale mining and reclamation monitoring to improve oversight of the changes resulting from open-pit mining. This can help us to understand the natural vegetation recovery mechanisms in disturbed areas without reclamation, as well as the reclamation status of reclaimed areas.
The current GEE-based algorithms are fast and effective, and are increasingly the preferred means of performing remote-sensing research. The method described in this paper monitored mining and reclamation in the Shengli open-pit mining area effectively. The disadvantage of this pixel-based method is that it lacks neighborhood information, and thus is not helpful for identifying sources of sudden interference. Reasonable explanations can be obtained only in limited scenarios, such as those of the mining and reclamation activities described in this paper. In the future, we may integrate object-oriented mining and reclamation monitoring methods and compare the results with pixelbased algorithms.