Land Cover Changes in Open-Cast Mining Complexes Based on High-Resolution Remote Sensing Data

Remote sensing technologies can play a fundamental role in the environmental assessment of open-cast mining and the accurate quantification of mine land rehabilitation efforts. Here, we developed a systematic geographic object-based image analysis (GEOBIA) approach to map the amount of revegetated area and quantify the land use changes in open-cast mines in the Carajás region in the eastern Amazon, Brazil. Based on high-resolution satellite images from 2011 to 2015 from different sensors (GeoEye, WorldView-3 and IKONOS), we quantified forests, cangas (natural metalliferous savanna ecosystems), mine land, revegetated areas and water bodies. Based on the GEOBIA approach, threshold values were established to discriminate land cover classes using spectral bands, the normalized difference vegetation index (NDVI), normalized difference water index (NDWI) and a light detection and range sensor (LiDAR) digital terrain model and slope map. The overall accuracy was higher than 90%, and the kappa indices varied between 0.82 and 0.88. During the observation period, the mining complex expanded, which led to the conversion of canga and forest vegetation to mine land. At the same time, the amount of revegetated area increased. Thus, we conclude that our approach is capable of providing consistent information regarding land cover changes in mines, with a special focus on the amount of revegetation necessary to fulfill environmental liabilities.


Introduction
The societal and environmental impacts of mining activities have become a great focus of public interest [1]. On the one hand, societies expect the production of inexpensive raw materials to drive the economy. The exploitation of mineral reserves should be embedded in efficient and responsible management plans that have a strong focus on the development of neighboring communities and sustainable land use. This balancing act has led to the emergence of environmental regulations in many countries, including Brazil [2], and a series of self-commitments by the mining industry [3,4].
The transformation of natural ecosystems during open-cast mining reduces the amount of wildlife habitat and endangers populations of rare, endemic or threatened species [5,6]. To reduce these environmental changes, the mitigation hierarchy offers a useful form of guidance [7,8]. According to this framework, measures should be taken to avoid and minimize any potential impacts before WorldView-3 satellite image in color-composite 1R2G3B acquired in 2015. The red line represents the boundaries of mining features; WP represents a waste pile; and N represents the mine pits (N is derived from Serra Norte, and the number refers to the original canga body, numbered consecutively for each mountain range in the CNF).

Remote Sensing Data Sources
High-resolution images from three different sensors ( IKONOS, 2012 GeoEye and 2015 WorldView-3) were acquired and used in this research. Table 1 provides detailed information about this remote sensing dataset. All images were acquired in August of each year (i.e., the dry season) to minimize the cloud cover and the spectral differences in vegetation cover as a function of seasonality along the analyzed time series. The spatial resolution differed among the images, i.e., 4 m for IKONOS, 2 m for GeoEye and 1.24 m for WorldView-3. Additionally, a digital terrain model was used to map the land cover and land use classes in the N4-N5 mining complex. Figure 2 illustrates the main steps of the digital image processing and GEOBIA proposed in this study. During mining operation, the original vegetation, i.e., forests and cangas, is logged, and topsoil and mining waste are removed and deposited in waste piles to guarantee access to high-quality iron ore. Benches in active mine pits are revegetated to reduce soil erosion and increase their stabilization. Furthermore, benches from permanent structures, such as waste piles, are subjected to permanent environmental rehabilitation to reduce the overall impact of mining enterprises on natural resources [14]. In both cases, similar mixtures of native and fast-growing and nonnative species are applied together with organic and inorganic fertilizers and fixing material via hydroseeding after bench preparation using biotextiles and/or manual trench digging. The species mixtures used for temporary revegetation in mining pits do not include tree species, as trees may destabilize the landscape by falling; however, the establishment of trees or shrubs from forest and canga ecosystems is implemented during the environmental rehabilitation of waste pile benches.

Remote Sensing Data Sources
High-resolution images from three different sensors ( IKONOS, 2012 GeoEye and 2015 WorldView-3) were acquired and used in this research. Table 1 provides detailed information about this remote sensing dataset. All images were acquired in August of each year (i.e., the dry season) to Remote Sens. 2020, 12, 611 4 of 21 minimize the cloud cover and the spectral differences in vegetation cover as a function of seasonality along the analyzed time series. The spatial resolution differed among the images, i.e., 4 m for IKONOS, 2 m for GeoEye and 1.24 m for WorldView-3. Additionally, a digital terrain model was used to map the land cover and land use classes in the N4-N5 mining complex. Figure 2 illustrates the main steps of the digital image processing and GEOBIA proposed in this study.

Atmospheric Corrections and Orthorectification
Conversions from digital number to ground reflectance data were carried out in the Atmospheric Correction (ATCOR) module of PCI Geomatica 2016 software [37]. Specifically, the ATCOR wizard was applied to calibrate the radiometric models considering the entire scene and to

Atmospheric Corrections and Orthorectification
Conversions from digital number to ground reflectance data were carried out in the Atmospheric Correction (ATCOR) module of PCI Geomatica 2016 software [37]. Specifically, the ATCOR wizard was applied to calibrate the radiometric models considering the entire scene and to generate the reflectance images. As images may suffer from deformation due to sensor inclination and altitude variations, the images were orthorectified after atmospheric correction. This orthorectification, including georeferencing of the imagery, was carried out in the OrthoEngine Tool of PCI Geomatica [38] using the digital terrain model as the basis for the coordinates (latitude, longitude) and altitude. The root mean square error was approximately one pixel in size. To map the land cover as accurately as possible, the pixel size of GeoEye was 2 m, whereas the IKONOS and WorldView-3 images were resampled to 2 × 2 m during the orthorectification process. Hence, the IKONOS image had the original reflectance preserved, while only WorldView-3 reflectance values needed to be resampled from nearest neighbor interpolation. In the next step, images were mosaicked to cover the complete study site.

Elaboration of Remote Sensing Indices
In addition to the original multispectral bands of the orbital sensors (B1-red, B2-green, B3-blue, and B4-infrared) and the light detection and ranging sensor (LiDAR, see next section for details), a digital terrain model derived from LiDAR was used for image interpretation. We also computed the normalized difference vegetation index-NDVI [39] and normalized difference water index-NDWI [39,40] (Figure 2). The WorldView-3 image presents two spectral bands in the infrared region; we used the spectral band with wavelengths between 770 and 895 nm, which is similar to the near infrared band from further sensors used in this study. The NDVI highlights the presence of vegetation, while the NDWI separates water surfaces from all further land classes. For both indices, higher values indicate the presence of more photosynthetically active vegetation [40].

LiDAR Data Processing
An overflight with a full waveform LiDAR was carried out in June 2012, with six elevation points per square meter. From that, the digital terrain model (DTM) was produced, and ground returns were separated from vegetation returns using the GroundFilter DTM by Triangular Irregular Network interpolation to generate a bare earth surface using the raster calculator in ArcMap. A digital slope map (SM) was created using geoprocessing tools. Both maps were generated in ArcGIS 10.3 software (ESRI) at 2 m spatial resolution.
Using the high-resolution satellite images alone, the cangas were often difficult to distinguish from old mining areas (mining areas that remained mining areas throughout the observation period), such as weathered slopes or rocky outcrops within the mine, which presented similar reflectance parameters as the canga areas. Hence, the inclusion of LiDAR (DTM and SM) data resulting from a single overflight increased the accuracy of the classification in these cases.

GEOBIA: Image Segmentation, Multilayer Calibration and Hierarchical Classification
Segmentation is the division of a scene into homogeneous parts to extract image objects [41]. The segmentation of the N4-N5 iron ore open-cast mine scenes was carried out by a multiresolution segmentation algorithm based on the homogeneity definition [42]. A four-date segmentation was conducted from all of the ground reflectance bands from 2011 and 2013 IKONOS, 2012 GeoEye and 2015 WorldView-3 with the same resampled pixel size (2 m). To detect land cover and land use changes, we carried out a segmentation process from four separate single-date images, as suggested by Descleé et al. [43] and Duveiller et al. [44]. Multidate segmentation allowed for the comparison of four single images based on objects with the same geometry, delineating spatially and spectrally consistent segments and avoiding misclassification [45]. Hence, it was possible to estimate the most accurate and rapid process, reducing additional processing efforts required to outline polygons [46]. The segmentation approach was developed on two levels using eCognition software (Trimble, München, Germany). First, the macrosegmentation aimed to segregate large features of the image using a less-refined segmentation (scale: 100, compactness: 0.5, and form: 0.1, using weight 1 for all 8 bands). In the second step, a more refined microsegmentation was performed (scale: 30, compactness: 0.5, and form: 0.1, using weight 25 for the NDVI band, 15 for the infrared band, and 1 for the other bands) (Figure 2). Different weights were chosen from a trial and error (heuristics) approach to enhance specific objects that were effectively differentiated in predetermined spectral bands and indices [22,47].
To calibrate the threshold values for the classification of segments, 35 ground control points were used for each class. This procedure was separately carried out for each segmentation step and resulted in the definition of the threshold values shown in Tables 2 and 3. A hierarchical classification was carried out from fuzzy logic [48]. Based on these threshold layer values, macrosegments were classified into forests, cangas, mining areas and water bodies using the threshold values shown in Table 2. Objects located in the boundary between canga and forest classes (i.e., canga-forest transitions) presented different spectral and topographic characteristics. After their identification, these objects were pooled into the canga class. The classification of microsegments in the mining area class enables the differentiation between revegetated and non-revegetated mining areas (Table 3). This refinement was necessary to identify revegetated and rehabilitated areas within mines and to separate them from misclassified forests, which were defined as complementary forests 1 and forests 2 and were grouped with the forest class during the first segmentation. Moreover, the DTM and SM derived from LiDAR data were important for discriminating the canga class from the mining areas. Table 2. Layers, functions and thresholds used to classify macrosegments in the classes of forests, cangas, mining areas and water bodies. NDVI represents the normalized difference vegetation index, NDWI represents the normalized difference water index, DTM represents the digital terrain model and SM represents the digital slope map.

Class
indicates open intervals; otherwise, intervals are closed.

Detection of Land Cover and Open-Cast Mine Changes
In the next step, the polygons were vectorized to calculate the areas of the different classes for each year [22]. To do so, the classifications obtained for each year were exported as shapefiles (.shp) and loaded into ArcGIS 10.3 software (ESRI, Readlands, CA, USA) before the area (in hectares) was calculated using the geometry calculation function present in the layer attributes table. For the analysis Remote Sens. 2020, 12, 611 8 of 21 of temporal changes, a conceptual tree of the possible typologies and their transitions was elaborated ( Figure 3). Changes between classes during the observation period were detected by comparison among polygons from maps of different dates. Each map was compared to its previous map by subtraction to generate a new thematic map resulting from this differentiation and a table of changes between classes [49]. The "dissolve" geoprocessing function was used for each date of the scene and used to unify polygons from each layer by class. Then, the polygons were compared between subsequent dates using the "Intersect" tool in ArcGIS to produce maps that

Detection of Land Cover and Open-Cast Mine Changes
In the next step, the polygons were vectorized to calculate the areas of the different classes for each year [22]. To do so, the classifications obtained for each year were exported as shapefiles (.shp) and loaded into ArcGIS 10.3 software (ESRI, Readlands, CA, USA) before the area (in hectares) was calculated using the geometry calculation function present in the layer attributes table. For the analysis of temporal changes, a conceptual tree of the possible typologies and their transitions was elaborated ( Figure 3). Changes between classes during the observation period were detected by comparison among polygons from maps of different dates. Each map was compared to its previous map by subtraction to generate a new thematic map resulting from this differentiation and a table of changes between classes [49]. The "dissolve" geoprocessing function was used for each date of the scene and used to unify polygons from each layer by class. Then, the polygons were compared between subsequent dates using the "Intersect" tool in ArcGIS to produce maps that

Classification Accuracy Assessment
To evaluate the accuracy of the GEOBIA classification, polygons were used as the sampling units [50]. It is important to mention that ground control points were not collected in the field. Hence, in each scene, 1024 calibration points were randomly stratified per class in the high-resolution images. This method guaranteed that the number of points was proportional to the size of the land use and land cover class in each scene. For these points, the GEOBIA classification was compared to actual land use and land cover to generate a confusion matrix (Table S1). The actual vegetation was evaluated by manual image analysis by two independent reviewers who were familiar with the scene. Then, the relation of omission and commission errors, the producer's and the user's accuracy, the kappa index per class, the overall kappa index, overall accuracy [50], quantity disagreement (QD) and allocation disagreement (AD) indices [51] were computed for each multitemporal mosaic.

Accuracy Assessment of Land Change
To assess the accuracy and estimate the area of the land change map, 1024 check points were selected to validate the land change map from 2011 to 2015. The objective of this analysis followed best practices to assess the accuracy of the change classification and estimate the area of change in terms of a classification error matrix [52,53]. The error matrix was used to cross-tabulate the land

Classification Accuracy Assessment
To evaluate the accuracy of the GEOBIA classification, polygons were used as the sampling units [50]. It is important to mention that ground control points were not collected in the field. Hence, in each scene, 1024 calibration points were randomly stratified per class in the high-resolution images. This method guaranteed that the number of points was proportional to the size of the land use and land cover class in each scene. For these points, the GEOBIA classification was compared to actual land use and land cover to generate a confusion matrix (Table S1). The actual vegetation was evaluated by manual image analysis by two independent reviewers who were familiar with the scene. Then, the relation of omission and commission errors, the producer's and the user's accuracy, the kappa index per class, the overall kappa index, overall accuracy [50], quantity disagreement (QD) and allocation disagreement (AD) indices [51] were computed for each multitemporal mosaic.

Accuracy Assessment of Land Change
To assess the accuracy and estimate the area of the land change map, 1024 check points were selected to validate the land change map from 2011 to 2015. The objective of this analysis followed best practices to assess the accuracy of the change classification and estimate the area of change in terms of a classification error matrix [52,53]. The error matrix was used to cross-tabulate the land change class labels allocated by the classification of the high-resolution images against the reference check points collected at the sample sites. Accuracy parameters derived from a sample error matrix of unchanged classes (forest, canga, mine, revegetation and water) and changed classes (forest to mine, canga to mine, mine to revegetation, and revegetation to mine) included overall accuracy, user accuracy of class i (the proportion of the area mapped as class i that was referenced as class i), and producer accuracy of class j (the proportion of the area referenced as class j that was mapped as class j) as proposed by [53].
To assess the accuracy of the land change map, the accuracy was reported in terms of estimated area proportions (p ij ) and not in terms of sample counts, n ij . Hence, p ij was substituted byp ij and applied for simple random and stratified sampling [53] using the following equation: where W i is the proportion of the area mapped as class i, n i is the sum of the mapped class at row (i, k) in the error matrix, andp ij is the sum over class k. The estimated area of class k is: where A is the total map area. The standard error of the estimated area S(Â k ) is given by: An approximate 95% confidence interval is obtained asÂ k ± 1.96 × S(Â k ). For details about the matrix nomenclature, consult Olofsson et al. [53].

High-Resolution Satellite Image Accuracy Assessment and Estimated Area of Land Change
Based on the evaluation of 1024 points per scene, the overall accuracy of the final land cover and land use maps varied between 90% and 94%; the highest value was obtained in 2011, and the lowest value was obtained in 2012. Therefore, the overall disagreement is very low, and the AD and QD are lower than 8% and 5%, respectively (Tables S1 and S2). The highest kappa index per class was found for forests and cangas in all scenes, while the lowest values were detected for water (except in 2011) and rehabilitated areas (Figure 4a). The AD and QD per class were lower than 5.1% and 2.75%, respectively (Figure 4b). The largest user's and producer's accuracy values were found in forest and mining areas (Figure 4c,d). The greatest omission and commission errors of classification were associated with the misclassification of revegetation and water areas (Table 4).
AD allowed for the evaluation of classification as spatial mismatches during map comparisons, thus contributing to the detection of false transitions. QD was important for computing areal differences in classes among maps [22]. Our results showed that the AD of each class was frequently higher than the QD except for the forest and canga classes in 2011 and canga class in 2013, indicating that the area of the classes tended to be accurate for this purpose. However, the land change detection between different maps tended to be less reliable than that of the computed areas; hence, we expected that none of the results would be significantly affected.
The accuracy assessment and estimated area of land change were calculated from Equation (2) and Equation (3). For example, the estimated area of the mine class wasÂ mine =p mine × A total = 0.184 × 10,953.7 = 1928.74 ha. The mapped area of mine (A mine,1 ) was 2019.44 ha; thus, the mine area was overestimated by 90.7 ha. Hence, a degradation area was identified with a 95% confidence interval equal to 2019.4 ± 71.9 ha.
The overall accuracy of the estimated area of land change was 96% ( Table 4). The user accuracy varied between 63% (from revegetation to mine) and 99% (unchanged mine), resulting in a slight underestimation of the area changed from revegetation to mine and an overestimation of the unchanged rehabilitation area ( Figure 5).    underestimation of the area changed from revegetation to mine and an overestimation of the unchanged rehabilitation area ( Figure 5).

Analysis of the Spatial-Temporal Distribution of Land Cover and Land Use Classes
The multiresolution classification with the GEOBIA approach effectively classified the high-resolution satellite images into land cover and open-cast mine classes. Figure 6 illustrates the multitemporal maps throughout the study site for the years 2011, 2012, 2013 and 2015, showing the reduction in forest and canga areas and the expansion of mining and revegetated areas.

Analysis of the Spatial-Temporal Distribution of Land Cover and Land Use Classes
The multiresolution classification with the GEOBIA approach effectively classified the high-resolution satellite images into land cover and open-cast mine classes. Figure 6  Between 2011 and 2015, the natural vegetation, i.e., cangas and forests, was reduced by 687 ha, while the mining and revegetated areas increased by 420 ha and 279 ha, respectively (Figure 7). Throughout the observation period, forests remained the largest unchanged class in the study site (Figure 8a).   (Figure 8a). During this period, the largest conversion of forest to mining area that was observed (207 ha) occurred due to the logging of the S2 waste pile area (Figure 8b). At the same time, the area of canga converted to mine area was lower than that in the following periods (42 ha). Furthermore, 171 ha of mine land was revegetated; however, 165 ha of this revegetated area was reconverted to mine land (Figure 8b).
In the subsequent observation period (2012 to 2013), the areas of unchanged mine and revegetated land increased by 10% and 15%, respectively (Figure 8b). The conversion of land from forest to mine was reduced by more than 50%, while 116 ha of canga was converted to mining area, i.e., an increase of more than 100% compared to the previous period (2011 to 2012). Although 199 ha of mine land was revegetated between 2012 and 2013, the amount of revegetated area increased only marginally because 149 ha of revegetated area was reconverted to mine land (Figure 8b).
Between 2013 and 2015, the unchanged mine area declined in response to the expanding revegetated area, reaching the highest values observed in this study (Figure 8a). This pattern is in response to the establishment of spontaneous vegetation on the S2 waste pile (Figure 1). The conversion of forest to mine increments remained stable (90 ha), while the conversion of canga to mine reached the highest value (212 ha) in response to the installation of the N4 WS mine (Figure 8b). The reuse of revegetated areas by mining activities reached the lowest value of the entire study period (129 ha) (Figure 8b).  This finding indicates a clear tendency of a decrease in canga area and its surrounding forest. Throughout the period of change (2011-2015), we observed that the mining area expanded mainly over canga (354 ha), followed by forest (264 ha). The revegetation of the mining area increased over time (378 ha), while its reuse by mining decreased (Figure 6). The reformulation of the northern (N) waste pile and the suppression of smaller portions of temporary revegetation in the mining area were responsible for the transformation of revegetated sites in mining areas. Figure 9 illustrates the From 2011 to 2012, unchanged canga occupied an area of 983 ha, while unchanged mine and revegetation areas occupied 2135 ha and 403 ha, respectively (Figure 8a). During this period, the largest conversion of forest to mining area that was observed (207 ha) occurred due to the logging of the S2 waste pile area (Figure 8b). At the same time, the area of canga converted to mine area was lower than that in the following periods (42 ha). Furthermore, 171 ha of mine land was revegetated; however, 165 ha of this revegetated area was reconverted to mine land (Figure 8b).
In the subsequent observation period (2012 to 2013), the areas of unchanged mine and revegetated land increased by 10% and 15%, respectively (Figure 8b). The conversion of land from forest to mine was reduced by more than 50%, while 116 ha of canga was converted to mining area, i.e., an increase of more than 100% compared to the previous period (2011 to 2012). Although 199 ha of mine land was revegetated between 2012 and 2013, the amount of revegetated area increased only marginally because 149 ha of revegetated area was reconverted to mine land (Figure 8b).
Between 2013 and 2015, the unchanged mine area declined in response to the expanding revegetated area, reaching the highest values observed in this study (Figure 8a). This pattern is in response to the establishment of spontaneous vegetation on the S2 waste pile (Figure 1). The conversion of forest to mine increments remained stable (90 ha), while the conversion of canga to mine reached the highest value (212 ha) in response to the installation of the N4 WS mine (Figure 8b). The reuse of revegetated areas by mining activities reached the lowest value of the entire study period (129 ha) (Figure 8b).
This finding indicates a clear tendency of a decrease in canga area and its surrounding forest. Throughout the period of change (2011-2015), we observed that the mining area expanded mainly over canga (354 ha), followed by forest (264 ha). The revegetation of the mining area increased over time (378 ha), while its reuse by mining decreased (Figure 6). The reformulation of the northern (N) waste pile and the suppression of smaller portions of temporary revegetation in the mining area were responsible for the transformation of revegetated sites in mining areas. Figure 9 illustrates the unchanged land cover and land use (LCLU) classes and the "from-to" change detection classes based on a bitemporal 2011-2015 mosaic image analysis.
Nevertheless, the proportion of revegetated areas in relation to mining areas increased from 24.63% in 2011 to 31.55% in 2015. This increase was partially due to the establishment of spontaneous vegetation at waste pile S2, which, despite being logged, was still not used during the observation period. When we zoomed in on the selected permanent structures, the annual/biannual increases became clearly visible in the following waste piles: Northeast 2-NW2, West-W and South 4-S4 ( Figure 10).  This finding indicates a clear tendency of a decrease in canga area and its surrounding forest. Throughout the period of change (2011-2015), we observed that the mining area expanded mainly over canga (354 ha), followed by forest (264 ha). The revegetation of the mining area increased over time (378 ha), while its reuse by mining decreased (Figure 6). The reformulation of the northern (N) waste pile and the suppression of smaller portions of temporary revegetation in the mining area were responsible for the transformation of revegetated sites in mining areas. Figure 9 illustrates the unchanged land cover and land use (LCLU) classes and the "from-to" change detection classes based on a bitemporal 2011-2015 mosaic image analysis. Nevertheless, the proportion of revegetated areas in relation to mining areas increased from 24.63% in 2011 to 31.55% in 2015. This increase was partially due to the establishment of spontaneous vegetation at waste pile S2, which, despite being logged, was still not used during the observation period. When we zoomed in on the selected permanent structures, the annual/biannual increases became clearly visible in the following waste piles: Northeast 2-NW2, West-W and South 4-S4 ( Figure 10).

Assessment of the High-Resolution Satellite Image Accuracy and the GEOBIA Approach
The two-step segmentation of the scenes and the utilization of eight different layers, including spectral bands, elevation data and SMs, allowed for the mapping of the spatial land use dynamics of the N4-N5 mining complex from the CNF between 2011 and 2015 with acceptable accuracy. The less detailed scale of the first segmentation step permitted the separation of the entire scene in terms of the areas of natural vegetation, i.e., canga and forest vegetation, and mining areas (including mining and revegetated areas). From the level 1 of segmentation, the mining area class was segregated into small-scale (>0.2 ha) revegetated and/or forest areas within the mine. The second segmentation presented a level of detail that was approximately three times greater than that of the first segmentation, allowing the identification of small and medium structures within the mine; this enhanced level of detail enabled us to map the revegetated areas within the mine operations. Thus, despite the use of different sensors that required individual atmospheric corrections and individual calibration, the resolution of our mapping was higher than that obtained from similar approaches [28]. The multispectral classification further allowed the separation of revegetated areas from areas covered by natural vegetation from the time series without preliminary knowledge. The GEOBIA approach presented here showed a straightforward method of mapping the revegetation inside mining areas, and it avoided misclassifications with other vegetation cover that occurred in the neighborhood.
Overall, the kappa indices and accuracies of the classification indicate marginal confusion; these mismatches occurred principally between small revegetated areas and areas of exposed soil in mining sites, and these results confirmed the difficulties in identifying and classifying small features [54]. These limits are set by the spatial resolution of the scenes. However, large revegetated areas from structures such as waste piles were well classified, as the spectral characteristics differentiated them from the surrounding mining areas or natural vegetation. In relation to allocation and quantity disagreement, AD was greater than QD most of the time, indicating that the area of the classes tended to be more accurate for this purpose. However, the change detection between different years tended to be less reliable than the computed areas [22]. The accuracy assessment of land change maps has contributed to estimates of the uncertainty of the area estimates [53], providing more efficient and rigorous methods to estimate the accuracy and area of land change in iron ore mining sites. The major limitation of the GEOBIA approach was determining an appropriate scale in the two levels of segmentation to improve the classification results. New developments must be made to try to calculate the segmentation scale based on the pixel size and the area of the main target to be mapped.
The individual calibration of scenes originating from different sensors resulted in threshold values that differed slightly between images. Nevertheless, the spectral differences in all scenes between land use classes, including the differences between natural vegetation and revegetated areas, enabled the detection of emerging structures such as areas of vegetation logging, i.e., land conversion from natural ecosystems to mine lands, and emerging rehabilitation sites, i.e., conversion of mine lands to revegetated areas. Thus, these spectral differences were sufficient to track the land use changes and accurately quantify the effectively revegetated areas. However, the assessment of revegetation quality still requires ground monitoring [13]. Specifically, the two-step segmentation of the high-resolution images enabled the identification of small revegetated structures, such as strait slopes, indicating the suitability of the proposed methodology for this study [22].
The evaluation of satellite images by segmentation and classification via eCognition as proposed here represents an effective tool to monitor the spatial extent of revegetation activities within open-cast mines. This technology may be helpful for monitoring the revegetation of mine lands conducted by mining companies and thus could become an important tool for environmental monitoring by official licensing agencies. Furthermore, the approach proposed here is suitable for monitoring revegetation activities; the intersection and quantification of revegetated areas with areas planted by operations contributes to the identification of revegetation constraints, e.g., inadequate planting periods or substrates disabling plant survival, thus contributing to more effective mine land revegetation.
The reduction in natural vegetation and the increase in mining area show that the N4-N5 Carajás mining complex is a fully expanding mine and indicates that the iron ore reserves are not yet exhausted. The amount of expansion (absolute as well as relative) of revegetated areas within the mining complex shows that the mine land has been revegetated, which is the first step to address environmental liabilities by mining companies. Ground points are necessary to assess the quality of the revegetation, i.e., to assess whether revegetated areas are resilient ecosystems as required by legislation [14]. The combination of ground investigations with additional remote investigations of the scene may result in approaches that can be used to map the vegetation structure [55] or species composition of revegetated areas, which may enable the tracking of the advance of revegetated sites towards rehabilitated ecosystems [56].

Revegetation Analysis from GEOBIA using High-Resolution Satellite Data
Revegetation activities in the N4-N5 mining complex are carried out for the temporary revegetation of benches to protect them from soil erosion and reduce dust loading during periods in which they are not in use. Furthermore, permanent revegetation is carried out on permanent structures, such as waste piles, to rehabilitate the biodiversity and ecosystem functioning of mine lands, reducing the overall impact of mining operations. Our classification was unable to differentiate between temporary and permanent rehabilitation areas, which is not surprising because similar seed mixtures and rehabilitation techniques, such as hydroseeding, biotextiles or manual trench digging, were applied for both purposes. Nevertheless, the high amount of revegetated area that was transformed to mining area (i.e., one-third of all revegetated areas mapped in 2011) indicates that a high amount of mine land greening is composed of temporary revegetation. In contrast, approximately two-thirds of the revegetated area mapped in 2011 may correspond to permanent revegetation, and this land is intended to rehabilitate mine land to improve long-term biodiversity and ecosystem services and reduce the overall impact of mining operations.
Our classification revealed that only a very small amount (0.42 ha) of revegetated area was confounded with natural forest vegetation during the observation period (less than 0.05% of the entire rehabilitated area). However, this figure shows that the revegetated mining structures maintained their spectral profiles during the observation period, which is expected because environmental rehabilitation in tropical ecosystems may span centuries [57], while the studied revegetated areas are, at most, decades old.
Assessing the surface mining area through high-resolution satellite images has a high potential to become an operational tool to monitor and evaluate the dynamics of land cover and land use changes in open-cast mining complexes. It is important to emphasize that high spatial resolution satellite imageries have been commercially available since 1999 with the IKONOS system. Therefore, high-resolution satellites at one meter and below have been operating for more than twenty years, thus providing a substantial amount of spatial and spectral information that is extremely useful for the recognition of the geomorphic features of mines [58] and analysis of different trajectory types of land cover and land use in mining areas [30]. Thus, the GEOBIA approach presented in this paper can be used to study other sites of open-cast mines. Enabling the accuracy assessment and the estimation of the area of land change can become a good practice for monitoring these complex man-made environments. Further methodological advances should focus on the recognition and discrimination of different stages of rehabilitation in mining areas (e.g., herbaceous, shrub and forest cover) from high-resolution satellite and unmanned aerial vehicle systems to remotely track environmental advances of revegetated areas.

Conclusions
Our classification permitted the mapping of land cover and open-cast mine changes in the N4-N5 mining complex with sufficient accuracy. After the atmospheric correction of images, the generation of the digital terrain and slope models, and calibration of the threshold values for individual scenes, the proposed methodology was shown to be sufficiently robust to monitor land use changes in mining sites, offering a powerful tool for all stakeholders. Specifically, the remote sensing of mine land revegetation dynamics was able to monitor the occurrence of revegetation activities and the attendance of the mitigation hierarchy as required by national environmental laws. It is important to highlight the importance of the proposed two-step segmentation and multispectral classification to accurately separate the land cover and land use classes in the mining scene. Combined with ground data, additional remote investigations of the scene should be used to develop classification approaches for revegetated sites according to their structure or species composition, which may be useful for tracking advances in rehabilitating sites towards natural ecosystems. Moreover, the use of multitemporal high spatial resolution satellite images can become an operational tool for the monitoring of cumulative changes in land use related to mining operations.