Next Article in Journal
Evaluation of Gridded Precipitation Datasets in Malaysia
Next Article in Special Issue
Progress in the Reconstruction of Terrain Relief Before Extraction of Rock Materials—The Case of Liban Quarry, Poland
Previous Article in Journal
Time Series Remote Sensing Data-Based Identification of the Dominant Factor for Inland Lake Surface Area Change: Anthropogenic Activities or Natural Events?
Previous Article in Special Issue
Biomass Estimation for Semiarid Vegetation and Mine Rehabilitation Using Worldview-3 and Sentinel-1 SAR Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Filipe Silveira Nascimento
1,2,
Markus Gastauer
2,
Pedro Walfir M. Souza-Filho
2,3,*,
Wilson R. Nascimento, Jr.
2,
Diogo C. Santos
2 and
Marlene F. Costa
4
1
Vale S.A. Mina de Águas Claras (MAC), Nova Lima, Minas Gerais 34006-270, Brazil
2
Instituto Tecnológico Vale. Belém, Pará 66055-090, Brazil
3
Geosciences Institute, Universidade Federal do Pará. Belém, Pará 66075-110, Brazil
4
Vale S.A. Gerência de Meio Ambiente Corredor Norte, São Luís, Maranhão 65085-582, Brazil
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(4), 611; https://doi.org/10.3390/rs12040611
Submission received: 30 October 2019 / Revised: 27 January 2020 / Accepted: 30 January 2020 / Published: 12 February 2020

Abstract

:
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.

Graphical Abstract

1. 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 an enterprise starts. Then, international good practices indicate opportunities for remediation, i.e., reverse the residual impacts, and these practices are considered before offsets are outlined to address unavoidable impacts [9]. The four pillars of the mitigation hierarchy process increase environmental sustainability by avoiding net losses or generating positive impacts [10,11] as well as by increasing social acceptance of mining [12].
Remediation includes revegetation, i.e., greening, and the progressive rehabilitation of biodiversity, ecosystem structure and ecosystem function in degraded, damaged or destroyed mine lands [2]. Although considerable progress can be achieved in short time periods [13], such mine land rehabilitation remains challenging, especially in tropical ecosystems, as issues related to species selection, biological invasions, and monitoring the effectiveness of rehabilitation activities are not fully resolved [14]. Uncertain rehabilitation trajectories or unclear outcomes of rehabilitation activities [15] require constant monitoring to effectively regulate mining [16]. In addition, information about the status of environmental rehabilitation is derived from ground data [13]. This monitoring requires accurate quantification of the temporal and spatial extents of revegetation and rehabilitation sites within the mines [17].
An increasing demand for accurate and timely information on the nature and extent of land use and land cover changes in and around open-cast mines highlights the importance of remote sensing methods, which can provide accurate tools to track these changes over time [18,19,20,21,22]. Remote sensing covers large areas and has a higher temporal frequency and lower costs than ground-based investigations. Thus, remote sensing may play a fundamental role in the environmental monitoring of rehabilitation activities in mining areas [1,23] given that geographic information systems, satellite images and digital classification systems are available for the automated detection of land use and land cover changes [24,25,26]. Specifically, the automated analysis of time series, as proposed by the Landsat-based detection of trends in disturbance and recovery (LandTrendr) approach [27], enables the successful tracking of revegetation dynamics within mines based on vegetation indices [28]. However, because these vegetation indices by themselves are unable to discriminate between natural vegetation and revegetated areas, the separation between classes is based on only historical evidence. Furthermore, the LandTrendr approach is restricted to Landsat imagery with a low spatial resolution (30 × 30 m) and thus is unable to adequately map revegetation activities on steep benches. More recently, high-resolution satellite images have been used to capture land use change dynamics [19,29,30] and mine reclamation patterns [31] via fine-scale land mining interpretation. Hence, we adopt a geographic object-based image analysis (GEOBIA) approach to detect small-scale vegetation suppression and revegetation activities in open-cast mines from high spatial resolution image datasets captured by spaceborne remote sensors.
The objective of this study was to provide a methodology to track land cover and land use changes in mines, including the spatial and temporal dynamics of revegetated mine lands, with a special focus on the recognition of small revegetated areas on steep benches. To do so, we quantified the amount of natural vegetation converted to mining areas as well as the amount of mine land revegetated between 2011 and 2015 in the largest iron ore open-cast mining complex in the world, which is situated in the Carajás National Forest (CNF), in eastern Amazon, Brazil.

2. Materials and Methods

2.1. Study Site

The study area includes the iron ore open-cast mines from the N4-N5 complex situated in the CNF in the watershed of the Itacaiúnas River, eastern Amazon (Figure 1). The study area has a mean temperature of approximately 25 °C and annual precipitation between 1900 and 2000 mm [32]. The climate of the region is tropical warm, with rainy summers and dry hot winters. The CNF is a protected area dominated by dense and open evergreen forests as well as by semideciduous submontane and montane forest formations that cover the hillsides and lower portions of the landscape. In two mountain ridges within the CNF, Serra Sul and Serra Norte, patches of banded iron formations outcrop on the mountain tops. These patches are covered by hyperdiverse, endemic savanna vegetation [33,34] that is locally called canga. These cangas cover the largest iron ore reserves in the world [35] and represent the major business in the region. Since the 1980s, the world’s largest iron ore mining complex has been in operation, and it extracts 120 Mt of iron ore annually [36].
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.

2.2. Remote Sensing Data Sources

High-resolution images from three different sensors (2011, 2013 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.

2.3. 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.

2.4. 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].

2.5. 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.

2.6. 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 Table 2 and Table 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.

2.7. 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 indicated the land use changes from 2011 to 2012, 2012 to 2013, 2013 to 2015, and 2011 to 2015 (Figure 2).

2.8. 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.

2.9. 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, nij. Hence, pij was substituted by p ^ ij and applied for simple random and stratified sampling [53] using the following equation:
p ^ i j =   W i n i j n i
where Wi 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, and p ^ i j is the sum over class k. The estimated area of class k is:
A ^ k =   p ^ k   ×   A t o t a l
where A is the total map area. The standard error of the estimated area S( A ^ k ) is given by:
S ( A ^ k )   =   S (   p ^ k   )   ×   A t o t a l
An approximate 95% confidence interval is obtained as A ^ k ± 1.96 × S( A ^ k ). For details about the matrix nomenclature, consult Olofsson et al. [53].

3. Results

3.1. 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 A ^ m i n e = p ^ m i n e × A t o t a l = 0.184 × 10,953.7 = 1928.74 ha. The mapped area of mine (Amine,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).

3.2. 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.
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).
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).

4. Discussion

4.1. 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].

4.2. 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.

5. 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.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/12/4/611/s1, Table S1. Confusion matrix for image classification: (a) 2011 Ikonos; (b) 2012 GeoEye; (c) 2013 Ikonos; and (d) 2015 WorldView-3. The matrix shows the number of verification points, omission (OE) and commission (CE) errors, user (UA) and producer’s (PA) accuracy, Kappa index per class, overall accuracy (OA). Table S2. Normalized confusion matrices: (a) 2011 Ikonos; (b) 2012 GeoEye; (c) 2013 Ikonos; and (d) 2015 WorldView-3. The matrix shows the number of verification points, omission (OE) and commission (CE) errors, user (UA) and producer’s (PA) accuracy, Kappa index per class, overall accuracy (OA), overall disagreement (OD); allocation disagreement (AD); and quantity disagreement (QD).

Author Contributions

Conceptualization, P.W.M.S.-F.; field work collection, F.S.N., W.R.N.J., D.C.S., P.W.M.S.-F., M.F.C.; digital image processing, GIS and formal analysis, accuracy assessments, F.S.N., W.R.N.J., D.C.S., P.W.M.S.-F., writing–original draft preparation, M.G., F.S.N., P.W.M.S.-F., review and editing, P.W.M.S.-F., M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank the Vale Institute of Technology (ITV) for its management support. The authors are also grateful to the Vale Mining Company for logistical support. PWMSF was supported by CNPq through research scholarships. This project was carried out in the Carajás National Forest with permission from IBAMA (SISBIO 35594-2).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Koruyan, K.; Deliormanli, A.H.; Karaca, Z.; Momayez, M.; Lu, H.; Yalçin, E. Remote sensing in management of mining land and proximate habitat. J. S. Afr. Inst. Min. Metall. 2012, 112, 667–672. [Google Scholar]
  2. Gastauer, M.; Souza Filho, P.W.M.; Ramos, S.J.; Caldeira, C.F.; Silva, J.R.; Siqueira, J.O.; Furtini Neto, A.E. Mine land rehabilitation in Brazil: Goals and techniques in the context of legal requirements. Ambio 2018, 48, 74–88. [Google Scholar] [CrossRef]
  3. Matlaba, V.J.; Maneschy, M.C.; Filipe dos Santos, J.; Mota, J.A. Socioeconomic dynamics of a mining town in Amazon: A case study from Canaã dos Carajás, Brazil. Miner. Econ. 2019, 32, 75–90. [Google Scholar] [CrossRef]
  4. Santos, J.F.; Maia, M.B.; Maneschy, M.C.; Matlaba, V.; Mota, J.A. Redes sociais ao longo da estrada de ferro Carajás na Amazônia Oriental. Finisterra 2018, 109, 149–166. [Google Scholar]
  5. Skirycz, A.; Castilho, A.; Chaparro, C.; Carvalho, N.; Tzotzos, G.; Siqueira, J.O. Canga biodiversity, a matter of mining. Front. Plant Sci. 2014, 5, 1–9. [Google Scholar] [CrossRef] [Green Version]
  6. Bisone, S.; Chatain, V.; Blanc, D.; Gautier, M.; Bayard, R.; Sanchez, F.; Gourdon, R. Geochemical characterization and modeling of arsenic behavior in a highly contaminated mining soil. Environ. Earth Sci. 2016, 75, 69621. [Google Scholar] [CrossRef]
  7. International Finance Corporation. Performance Standard 6: Biodiversity Conservation and Sustainable Management of Natural Resources; IFC: Washington, DC, USA, 2012. [Google Scholar]
  8. Maron, M.; Brownlie, S.; Bull, J.W.; Evans, M.C.; von Hase, A.; Quétier, F.; Watson, J.E.M.; Gordon, A. The many meanings of no net loss in environmental policy. Nat. Sustain. 2018, 1, 19–27. [Google Scholar] [CrossRef] [Green Version]
  9. Ekstrom, J.; Bennun, L.; Mitchell, R. A Cross-Sector Guide for Implementing the Mitigation Hierarchy; Cambridge: London, UK, 2015. [Google Scholar]
  10. Bull, J.W.; Suttle, K.B.; Gordon, A.; Singh, N.J.; Milner-Gulland, E.J. Biodiversity offsets in theory and practice. Oryx 2013, 47, 369–380. [Google Scholar] [CrossRef] [Green Version]
  11. Bull, J.W.; Gordon, A.; Watson, J.E.M.; Maron, M. Seeking convergence on the key concepts in ‘no net loss’ policy. J. Appl. Ecol. 2016, 53, 1686–1693. [Google Scholar] [CrossRef]
  12. Moffat, K.; Zhang, A. The paths to social licence to operate: An integrative model explaining community acceptance of mining. Resour. Policy 2014, 39, 61–70. [Google Scholar] [CrossRef] [Green Version]
  13. Gastauer, M.; Caldeira, C.F.; Ramos, S.J.; Trevelin, L.C.; Santiago, A.; Aurélio, M.; Carneiro, C.; Coelho, F.; Silva, R.; Souza-Filho, P.W.M.; et al. Integrating Environmental Variables by Multivariate Ordination enables the Reliable Estimation of Mineland Rehabilitation Status. J. Environ. Manag. 2019, 256, 109894. [Google Scholar] [CrossRef]
  14. Gastauer, M.; Silva, J.R.; Caldeira Junior, C.F.; Ramos, S.J.; Souza Filho, P.W.M.; Furtini Neto, A.E.; Siqueira, J.O. Mine land rehabilitation: Modern ecological approaches for more sustainable mining. J. Clean. Prod. 2018, 172, 1409–1422. [Google Scholar] [CrossRef]
  15. Perring, M.P.; Standish, R.J.; Hobbs, R.J. Incorporating novelty and novel ecosystems into restoration planning and practice in the 21st century. Ecol. Process. 2013, 2, 1–8. [Google Scholar] [CrossRef] [Green Version]
  16. Padmanaban, R.; Bhowmik, A.; Cabral, P. A Remote Sensing Approach to Environmental Monitoring in a Reclaimed Mine Area. ISPRS Int. J. Geo Inf. 2017, 6, 401. [Google Scholar] [CrossRef] [Green Version]
  17. Townsend, P.A.; Helmers, D.P.; Kingdon, C.C.; McNeil, B.E.; de Beurs, K.M.; Eshleman, K.N. Changes in the extent of surface mining and reclamation in the Central Appalachians detected using a 1976-2006 Landsat time series. Remote Sens. Environ. 2009, 113, 62–72. [Google Scholar] [CrossRef]
  18. Asner, G.P.; Llactayo, W.; Tupayachi, R.; Luna, E.R. Elevated rates of gold mining in the Amazon revealed through high-resolution monitoring. Proc. Natl. Acad. Sci. USA 2013, 110, 18454–18459. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, W.; Li, X.; He, H.; Wang, L. A Review of Fine-Scale Land Use and Land Cover Classification in Open-Pit Mining Areas by Remote Sensing Techniques. Remote Sens. 2018, 10, 15. [Google Scholar] [CrossRef] [Green Version]
  20. Sonter, L.J.; Barrett, D.J.; Soares-Filho, B.S.; Moran, C.J. Global demand for steel drives extensive land-use change in Brazil’s Iron Quadrangle. Glob. Environ. Chang. 2014, 26, 63–72. [Google Scholar] [CrossRef]
  21. Souza-Filho, P.W.M.; de Souza, E.B.; Silva Júnior, R.O.; Nascimento, W.R.; Versiani de Mendonça, B.R.; Guimarães, J.T.F.; Dall’Agnol, R.; Siqueira, J.O. Four decades of land-cover, land-use and hydroclimatology changes in the Itacaiúnas River watershed, southeastern Amazon. J. Environ. Manag. 2016, 167, 175–184. [Google Scholar] [CrossRef]
  22. Souza-Filho, W.P.; Nascimento, R.W.; Santos, C.D.; Weber, J.E.; Silva, O.R.; Siqueira, O.J. A GEOBIA Approach for Multitemporal Land-Cover and Land-Use Change Analysis in a Tropical Watershed in the Southeastern Amazon. Remote Sens. 2018, 10, 1683. [Google Scholar] [CrossRef] [Green Version]
  23. Karan, S.K.; Samadder, S.R.; Maiti, S.K. Assessment of the capability of remote sensing and GIS techniques for monitoring reclamation success in coal mine degraded lands. J. Environ. Manag. 2016, 182, 272–283. [Google Scholar] [CrossRef] [PubMed]
  24. Demirel, N.; Emil, M.K.; Duzgun, H.S. Landuse change detection in a surface coal mine using multi-temporal high-resolution satellite images. Int. J. Min. Reclam. Environ. 2011, 25, 342–349. [Google Scholar] [CrossRef]
  25. Vogelmann, J.E.; Gallant, A.L.; Shi, H.; Zhu, Z. Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data. Remote Sens. Environ. 2016, 185, 258–270. [Google Scholar] [CrossRef] [Green Version]
  26. Ma, B.; Chen, Y.; Zhang, S.; Li, X. Remote Sensing Extraction Method of Tailings Ponds in Ultra-Low-Grade Iron Mining Area Based on Spectral Characteristics and Texture Entropy. Entropy 2018, 20, 345. [Google Scholar] [CrossRef] [Green Version]
  27. Kennedy, R.E.; Yang, Z.; Cohen, W.B. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 1. LandTrendr—Temporal segmentation algorithms. Remote Sens. Environ. 2010, 114, 2897–2910. [Google Scholar] [CrossRef]
  28. Yang, Y.; Erskine, P.D.; Lechner, A.M.; Mulligan, D.; Zhang, S.; Wang, Z. Detecting the dynamics of vegetation disturbance and recovery in surface mining area via Landsat imagery and LandTrendr algorithm. J. Clean. Prod. 2018, 178, 353–362. [Google Scholar] [CrossRef]
  29. Isidro, M.C.; McIntyre, N.; Lechner, M.A.; Callow, I. Applicability of Earth Observation for Identifying Small-Scale Mining Footprints in a Wet Tropical Region. Remote Sens. 2017, 9, 945. [Google Scholar] [CrossRef] [Green Version]
  30. Demirel, N.; Emil, M.K.; Duzgun, H.S. Surface coal mine area monitoring using multi-temporal high-resolution satellite imagery. Int. J. Coal Geol. 2011, 86, 3–11. [Google Scholar] [CrossRef]
  31. Maxwell, A.E.; Warner, T.A.; Strager, M.P.; Pal, M. Combining RapidEye Satellite Imagery and Lidar for Mapping of Mining and Mine Reclamation. Photogramm. Eng. Remote Sens. 2014, 2, 179–189. [Google Scholar] [CrossRef]
  32. Piló, L.B.; Auler, A.S.; Martins, F. Carajás National Forest: Iron Ore Plateaus and Caves in Southeastern Amazon. In Landscapes and Landforms of Brazil; Vieira, B.C., Salgado, A.A.R., Santos, L.J.C., Eds.; Springer: Dordrecht, The Netherlands, 2015; pp. 273–283. ISBN 978-94-017-8023-0. [Google Scholar]
  33. Viana, P.L.; Mota, N.F.D.O.; Gil, A.D.S.B.; Salino, A.; Zappi, D.C.; Harley, R.M.; Ilkiu-Borges, A.L.; Secco, R.D.S.; Almeida, T.E.; Watanabe, M.T.C.; et al. Flora das cangas da Serra dos Carajás, Pará, Brasil: História, área de estudos e metodologia. Rodriguésia 2016, 67, 1107–1124. [Google Scholar] [CrossRef]
  34. de Mota, N.F.O.; Silva, L.V.C.; Martins, F.D.; Viana, P.L. Vegetação sobre sistemas ferruginosos da Serra dos Carajás. In Geossistemas ferruginosos do Brasil: Áreas prioritárias para conservação da diversidade geológica e biológica, patrimônio cultural e serviços ambientais; Carmo, F.F., Kamino, L.H., Eds.; 3iEditora: Belo Horizonte, Brazil, 2015; pp. 289–315. [Google Scholar]
  35. Resende, N.P. Carajás: Memória da descoberta; Editora Gráfica Stampa: Rio de Janeiro, Brazil, 2009. [Google Scholar]
  36. Paradella, W.R.; Ferretti, A.; Mura, J.C.; Colombo, D.; Gama, F.F.; Tamburini, A.; Santos, A.R.; Novali, F.; Galo, M.; Camargo, P.O.; et al. Mapping surface deformation in open pit iron mines of Carajás Province (Amazon Region) using an integrated SAR analysis. Eng. Geol. 2015, 193, 61–78. [Google Scholar] [CrossRef] [Green Version]
  37. PCI. Geomatica Image processing with Focus. In Geomatica I Course Guide, Version 0.2; PCI Geomatica Enterprises: Markham, ON, Canada, 2015; pp. 59–79. [Google Scholar]
  38. PCI. Geomatica Geomatica Orthoengine: Course Exercises; PCI Geomatica Enterprises: Markham, ON, Canada, 2015; p. 172. [Google Scholar]
  39. Tarpley, J.D.; Schneider, S.R.; Money, R.L. Global Vegetation Indices from the NOAA-7 Meteorological Satellite. J. Clim. Appl. Meteorol. 1984, 23, 491–494. [Google Scholar] [CrossRef]
  40. McFEETERS, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  41. Blaschke, T. Object based image analysis for remote sensing. ISPRS J. Photogramm. Remote Sens. 2010, 65, 2–16. [Google Scholar] [CrossRef] [Green Version]
  42. Blaschke, T.; Hay, G.J.; Kelly, M.; Lang, S.; Hofmann, P.; Addink, E.; Queiroz Feitosa, R.; van der Meer, F.; van der Werff, H.; van Coillie, F.; et al. Geographic Object-Based Image Analysis—Towards a new paradigm. ISPRS J. Photogramm. Remote Sens. 2014, 87, 180–191. [Google Scholar] [CrossRef] [Green Version]
  43. Desclée, B.; Bogaert, P.; Defourny, P. Forest change detection by statistical object-based method. Remote Sens. Environ. 2006, 102, 1–11. [Google Scholar] [CrossRef]
  44. Duveiller, G.; Defourny, P.; Desclée, B.; Mayaux, P. Deforestation in Central Africa: Estimates at regional, national and landscape levels by advanced processing of systematically-distributed Landsat extracts. Remote Sens. Environ. 2008, 112, 1969–1981. [Google Scholar] [CrossRef]
  45. Souza-Filho, P.W.M.; Giannini, T.C.; Jaffé, R.; Giulietti, A.M.; Santos, D.C.; Nascimento, W.R., Jr.; Guimarães, J.T.F.; Costa, M.F.; Imperatriz-Fonseca, V.L.; Siqueira, J.O. Mapping and quantification of ferruginous outcrop savannas in the Brazilian Amazon: A challenge for biodiversity conservation. PLoS ONE 2019, 14, e0211095. [Google Scholar] [CrossRef] [Green Version]
  46. Raši, R.; Bodart, C.; Stibig, H.-J.; Eva, H.; Beuchle, R.; Carboni, S.; Simonetti, D.; Achard, F. An automated approach for segmenting and classifying a large sample of multi-date Landsat imagery for pan-tropical forest monitoring. Remote Sens. Environ. 2011, 115, 3659–3669. [Google Scholar] [CrossRef]
  47. Baatz, M.; Schäpe, A. Multiresolution Segmentation: An Optimization Approach for High Quality Multi-Scale Image Segmentation. 2000. Available online: http://www.agit.at/papers/2000/baatz_FP_12.pdf (accessed on 23 December 2012).
  48. Rejaur Rahman, M.; Saha, S.K. Multi-resolution segmentation for object-based classification and accuracy assessment of land use/land cover classification using remotely sensed data. J. Indian Soc. Remote Sens. 2008, 36, 189–201. [Google Scholar] [CrossRef]
  49. Lu, D.; Batistella, M.; Li, G.; Moran, E.; Hetrick, S.; Freitas, C.D.C.; Dutra, L.V.; Sant’Anna, S.J.S. Land use/cover classification in the Brazilian Amazon using satellite images. Pesqui. Agropecu. Bras. 2012, 47, 10–1590. [Google Scholar] [CrossRef] [PubMed]
  50. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices, 2nd ed.; Taylor & Francis: London, UK, 2008. [Google Scholar]
  51. Pontius, R.G., Jr.; Millones, M. Death to Kappa: Birth of quantity disagreement and allocation disagreement for accuracy assessment. Int. J. Remote Sens. 2011, 32, 4407–4429. [Google Scholar] [CrossRef]
  52. Gallego, F.J. Remote sensing and land cover area estimation. Int. J. Remote Sens. 2004, 25, 3019–3047. [Google Scholar] [CrossRef]
  53. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  54. Puertas, O.L.; Brenning, A.; Meza, F.J. Balancing misclassification errors of land cover classification maps using support vector machines and Landsat imagery in the Maipo river basin (Central Chile, 1975–2010). Remote Sens. Environ. 2013, 137, 112–123. [Google Scholar] [CrossRef]
  55. Johansen, K.; Erskine, P.D.; McCabe, M.F. Using Unmanned Aerial Vehicles to assess the rehabilitation performance of open cut coal mines. J. Clean. Prod. 2019, 209, 819–833. [Google Scholar] [CrossRef]
  56. Kimball, S.; Lulow, M.; Sorenson, Q.; Balazs, K.; Fang, Y.-C.; Davis, S.J.; O’Connell, M.; Huxman, T.E. Cost-effective ecological restoration. Restor. Ecol. 2015, 23, 800–810. [Google Scholar] [CrossRef]
  57. Chazdon, R.L. Tropical forest regeneration. Bol. Mus. Para. Emílio Goeldi. Cienc. Nat 2012, 7, 195–218. [Google Scholar]
  58. Chen, J.; Li, K.; Chang, K.-J.; Sofia, G.; Tarolli, P. Open-pit mining geomorphic feature characterisation. Int. J. Appl. Earth Obs. Geoinf. 2015, 42, 76–86. [Google Scholar] [CrossRef]
Figure 1. Location map of the N4-N5 mining complex indicating mine pits (N) and waste piles (WP) in relation to the Itacaiunas River watershed, southeastern Amazon, Brazil, South America. 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).
Figure 1. Location map of the N4-N5 mining complex indicating mine pits (N) and waste piles (WP) in relation to the Itacaiunas River watershed, southeastern Amazon, Brazil, South America. 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).
Remotesensing 12 00611 g001
Figure 2. Flowchart including the main steps of digital image processing and geographic object-based image analysis (GEOBIA).
Figure 2. Flowchart including the main steps of digital image processing and geographic object-based image analysis (GEOBIA).
Remotesensing 12 00611 g002
Figure 3. Land-cover and land-use classes detected in the N4-N5 mining complex and possible transitions among them. Bold lines indicate most frequent transitions observed in this study.
Figure 3. Land-cover and land-use classes detected in the N4-N5 mining complex and possible transitions among them. Bold lines indicate most frequent transitions observed in this study.
Remotesensing 12 00611 g003
Figure 4. Kappa index (a), allocation and quantity disagreements (b), user accuracy (c) and producer accuracy (d) per class for 2011, 2012, 2013, and 2015 images in the N4-N5 Carajás mining complex. F is forest, C is canga, W is water, M is mine and R is revegetated area.
Figure 4. Kappa index (a), allocation and quantity disagreements (b), user accuracy (c) and producer accuracy (d) per class for 2011, 2012, 2013, and 2015 images in the N4-N5 Carajás mining complex. F is forest, C is canga, W is water, M is mine and R is revegetated area.
Remotesensing 12 00611 g004
Figure 5. Estimated area ± margin of error normalized by mapped area.
Figure 5. Estimated area ± margin of error normalized by mapped area.
Remotesensing 12 00611 g005
Figure 6. Classified high-resolution 2011 Ikonos, 2012 GeoEye, 2013 Ikonos and 2015 WorldView-3 images illustrating the temporal and spatial changes in forest, canga, open-cast mine and revegetated (REV) land cover 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.
Figure 6. Classified high-resolution 2011 Ikonos, 2012 GeoEye, 2013 Ikonos and 2015 WorldView-3 images illustrating the temporal and spatial changes in forest, canga, open-cast mine and revegetated (REV) land cover 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.
Remotesensing 12 00611 g006
Figure 7. Land use classes in the N4-N5 Carajás mining complex between 2011 and 2015. REV indicates revegetated areas.
Figure 7. Land use classes in the N4-N5 Carajás mining complex between 2011 and 2015. REV indicates revegetated areas.
Remotesensing 12 00611 g007
Figure 8. Conversion (in ha) of major land cover and land use classes between 2011 and 2015. (a) Unchanged canga, forest, mine and revegetated areas (REV). (b) Conversion of forest to mine, canga to mine, mine to REV, and REV to mine.
Figure 8. Conversion (in ha) of major land cover and land use classes between 2011 and 2015. (a) Unchanged canga, forest, mine and revegetated areas (REV). (b) Conversion of forest to mine, canga to mine, mine to REV, and REV to mine.
Remotesensing 12 00611 g008
Figure 9. Distribution of LULC changes from the GEOBIA during the periods 2011–2012, 2012–2013, 2013–2015, and 2011–2015 using the “from-to” change detection approach.
Figure 9. Distribution of LULC changes from the GEOBIA during the periods 2011–2012, 2012–2013, 2013–2015, and 2011–2015 using the “from-to” change detection approach.
Remotesensing 12 00611 g009aRemotesensing 12 00611 g009b
Figure 10. Increase in revegetated area (REV) on three selected waste piles (South 4 [S 4], West [W] and Northwest 2 [NW 2]) from the N4–N5 mining complex from Carajás, Brazil, between 2011 and 2015. WorldView 3 satellite image in color-composite 1R2G3B acquired in 2015.
Figure 10. Increase in revegetated area (REV) on three selected waste piles (South 4 [S 4], West [W] and Northwest 2 [NW 2]) from the N4–N5 mining complex from Carajás, Brazil, between 2011 and 2015. WorldView 3 satellite image in color-composite 1R2G3B acquired in 2015.
Remotesensing 12 00611 g010
Table 1. General characteristics of WorldView-3, GeoEye and IKONOS images acquired during the dry season (August) used in this study.
Table 1. General characteristics of WorldView-3, GeoEye and IKONOS images acquired during the dry season (August) used in this study.
Satellite
Acquisition date
WorldView-3
1 August 2015
GeoEye
1 July 2012
Ikonos
23 May 2011; 22 July 2013
Spectral Resolution
Coastal400–450 nm------
Blue450–510 nm450–520 nm450–520 nm
Green510–580 nm520–600 nm520–600 nm
Yellow585–625 nm------
Red630–690 nm625–695 nm630–690 nm
Red Edge705–745 nm------
Near Infrared 1770–895 nm760–900 nm760–900 nm
Near Infrared 2860–1040 nm
Panchromatic450–800 nm450–900 nm450–900 nm
Spatial Resolution
Panchromatic0.3 m0.5 m1 m
Multispectral1.24 m2 m4 m
Radiometric Quantification11 bits per pixel11 bits per pixel11 bits per pixel
Scene Size13.1 km15.2 km11.3 km
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.
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.
ClassLayerIkonos 2011GeoEye 2012Ikonos 2013WorldView 2015
ForestsB1: Red-0–1.9--
B2: Green2.1–5.5 *2.1–5.5 *2.1–5.5 *1.7–38 *
B3: Blue2.8–6.3 *1.8–6.3 *2.5–7.3 *0.9–6.3 *
B4: Infrared----
B5: NDVI0.85–10.76–10.65–10.81–1
B6: NDWI----
B7: DTM----
B8: SM----
CangasB1: Red1.3–2.71.3–2.71.3–3.41.3–2.5
B2: Green3.8–5.75 *3.8–5.75 *3.8–5.75 *3.1–5.75 *
B3: Blue4.5–8.4 *4.5–7.4 *4.5–9.3 *4.5–7.4 *
B4: Infrared18–26 *10–18.5 *10–18 *10–18.5 *
B5: NDVI0.78–0.88--0.64–0.8
B6: NDWI−0.7–−0.4--−0.47–−0.33
B7: DTM---562–700
B8: SM0–17.5 *0–17.1 *0–17.5 *0–23 *
Complementary cangas (threshold condition: objects adjoining canga edges)B1: Red0.85–3.3 *0.85–3.3 *0.85–3.6 *0.85–3.3 *
B2: Green3–6.5 *3–6.5 *3–6.5 *3–6.5 *
B3: Blue3.1–9.6 *3.1–9.6 *3.1–9.9 *3.1–9.6 *
B4: Infrared----
B5: NDVI0.74–0.940.72–0.940.57–0.940.67–0.94
B6: NDWI----
B7: DTM-0–7000–7000–722
B8: SM----
Mining areasB1: Red0.27–7.7 *0.27–9.6 *0.27–10 *0.27–9.6 *
B2: Green1.6–13.5 *1.6–14.6 *1.6–17 *1.6–14.6 *
B3: Blue6–254.3–474.3–473.5–47
B4: Infrared----
B5: NDVI----
B6: NDWI−9.5–0.35 *−9.5–0.35 *−9.5–0.35 *−9.5–0.35 *
B7: DTM----
B8: SM----
Water bodiesB1: Red----
B2: Green----
B3: Blue----
B4: Infrared----
B5: NDVI----
B6: NDWI0.1–1 *0.1–1 *0.1–1 *0.1–1 *
B7: DTM----
B8: SM----
* indicates open intervals; otherwise, intervals are closed.
Table 3. Layers, functions and thresholds used to classify microsegments from mining areas in the classes of forests, revegetated sites and rehabilitated sites. 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.
Table 3. Layers, functions and thresholds used to classify microsegments from mining areas in the classes of forests, revegetated sites and rehabilitated sites. 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.
ClassLayersIkonos 2011GeoEye 2012Ikonos 2013WorldView 2015
Complementary forests 1B1: Red----
B2: Green2.1–5.5 *2.1–5.5 *2.1–5.5 *2.1–5.5 *
B3: Blue2.8–6 *2.8–6 *2.8–7 *2.8–6 *
B4: Infrared----
B5: NDVI0.87–1 *0.87–1 *0.73–1 *0.87–1 *
B6: NDWI----
B7: DTM----
B8: SM----
Revegetated and rehabilitated sitesB1: Red1.5–5 *1.4–5 *1–5 *1–5 *
B2: Green4.5–11 *3.3–11 *2–11 *2–11 *
B3: Blue5.5–14.2 *3.6–14.2 *2.3–14.2 *2.3–14.2 *
B4: Infrared----
B5: NDVI0.7–0.90.65–0.920.56–0.920.6–0.92
B6: NDWI−1–−0.3 *--−1–−0.3 *
B7: DTM----
B8: SM----
Complementary forests 2B1: Red----
B2: Green----
B3: Blue----
B4: Infrared----
B5: NDVI0.78–10.78–10.78–10.78–1
B6: NDWI----
B7: DTM----
B8: SM----
* indicates open intervals; otherwise, intervals are closed.
Table 4. Error matrix expressed in terms of classification of the land change map (A) and estimated proportions of areas (B). W = mapped area proportion, UM = unchanged mine, M-REV = from mine to revegetation, C-M = from canga to mine, UC = unchanged canga, F-M = from forest to mine, UF = unchanged forest, UREV = unchanged revegetation, REV-M = from revegetation to mine, ME = Margin of error, and NME = Normalized margin of Error.
Table 4. Error matrix expressed in terms of classification of the land change map (A) and estimated proportions of areas (B). W = mapped area proportion, UM = unchanged mine, M-REV = from mine to revegetation, C-M = from canga to mine, UC = unchanged canga, F-M = from forest to mine, UF = unchanged forest, UREV = unchanged revegetation, REV-M = from revegetation to mine, ME = Margin of error, and NME = Normalized margin of Error.
(A) Error Matrix of Classification of the Land Change Map (Line) against the Reference Data (Column) for the Sample Sites
AreaClassesUMM-REVC-MUCF-MUFUREVREV-MTotals
1928.74UM1631000010165
377.54M-REV42700001032
353.66C-M10300000031
671.71UC00351010055
264.06F-M00101401117
6790.13UF20030617110633
378.08UREV10010556164
189.81REV-M40000241727
10953.7Totals1752834551462574191024
Producer’s accuracy93.196.488.292.7100.098.775.789.5
User’s accuracy98.884.496.892.782.497.587.563.0
Kappa per class0.990.840.970.920.820.940.870.62
Agreement163.0027.0030.0051.0014.00617.0056.0017.00975.0
By chance28.200.881.032.950.23386.354.630.50424.8
Overall accuracy =0.952Kappa index =0.918
(B) Error Matrix by Estimated Proportions of Areas
WClassesUMM-REVC-MUCF-MUFUREVREV-MTotals
0.176UM0.1740.0010.0000.0000.0000.0000.0010.0000.176
0.034M-REV0.0040.0290.0000.0000.0000.0000.0010.0000.034
0.032C-M0.0010.0000.0310.0000.0000.0000.0000.0000.032
0.061UC0.0000.0000.0030.0570.0000.0010.0000.0000.061
0.024F-M0.0000.0000.0010.0000.0200.0000.0010.0010.024
0.620UF0.0020.0000.0000.0030.0000.6040.0110.0000.620
0.035UREV0.0010.0000.0000.0010.0000.0030.0300.0010.035
0.017REV-M0.0030.0000.0000.0000.0000.0010.0030.0110.017
1.000Totals0.1840.0300.0360.0600.0200.6090.0470.0131.0
Producer’s accuracy94.496.586.894.2100.099.264.11884.789
User’s accuracy98.884.496.892.782.497.587.50062.963
Area (ha)2019.4330.2394.4660.9217.56674.3515.96140.95
ME (95%)71.953.455.560.249.392.091.648.0
Area (ha)2019 ± 72330 ± 53394 ± 56661 ± 60217 ± 496674 ± 92516 ± 92141 ± 48
Overall accuracy =0.96
Normalized area1.0470.8751.1150.9840.8240.9831.3650.743
NME0.0370.1410.1570.0900.1870.0140.2420.253
Standard error36.66327.25528.32530.70025.16646.95746.75324.481
Standard deviation0.1090.3630.1770.2600.3810.1570.3310.483

Share and Cite

MDPI and ACS Style

Nascimento, F.S.; Gastauer, M.; Souza-Filho, P.W.M.; Nascimento, W.R., Jr.; Santos, D.C.; Costa, M.F. Land Cover Changes in Open-Cast Mining Complexes Based on High-Resolution Remote Sensing Data. Remote Sens. 2020, 12, 611. https://doi.org/10.3390/rs12040611

AMA Style

Nascimento FS, Gastauer M, Souza-Filho PWM, Nascimento WR Jr., Santos DC, Costa MF. Land Cover Changes in Open-Cast Mining Complexes Based on High-Resolution Remote Sensing Data. Remote Sensing. 2020; 12(4):611. https://doi.org/10.3390/rs12040611

Chicago/Turabian Style

Nascimento, Filipe Silveira, Markus Gastauer, Pedro Walfir M. Souza-Filho, Wilson R. Nascimento, Jr., Diogo C. Santos, and Marlene F. Costa. 2020. "Land Cover Changes in Open-Cast Mining Complexes Based on High-Resolution Remote Sensing Data" Remote Sensing 12, no. 4: 611. https://doi.org/10.3390/rs12040611

APA Style

Nascimento, F. S., Gastauer, M., Souza-Filho, P. W. M., Nascimento, W. R., Jr., Santos, D. C., & Costa, M. F. (2020). Land Cover Changes in Open-Cast Mining Complexes Based on High-Resolution Remote Sensing Data. Remote Sensing, 12(4), 611. https://doi.org/10.3390/rs12040611

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop