Applicability of Earth Observation for Identifying Small-Scale Mining Footprints in a Wet Tropical Region

The unpredictable climate in wet tropical regions along with the spatial resolution limitations of some satellite imageries make detecting and mapping artisanal and small-scale mining (ASM) challenging. The objective of this study was to test the utility of Pleiades and SPOT imagery with an object-based support vector machine (OB-SVM) classifier for the multi-temporal remote sensing of ASM and other land cover including a large-scale mine in the Didipio catchment in the Philippines. Historical spatial data on location and type of ASM mines were collected from the field and were utilized as training data for the OB-SVM classifier. The classification had an overall accuracy between 87% and 89% for the three different images—Pleiades-1A for the 2013 and 2014 images and SPOT-6 for the 2016 image. The main land use features, particularly the Didipio large-scale mine, were well identified by the OB-SVM classifier, however there were greater commission errors for the mapping of small-scale mines. The lack of consistency in their shape and their small area relative to pixel sizes meant they were often not distinguished from other land clearance types (i.e., open land). To accurately estimate the total area of each land cover class, we calculated bias-adjusted surface areas based on misclassification values. The analysis showed an increase in small-scale mining areas from 91,000 m2—or 0.2% of the total catchment area—in March 2013 to 121,000 m2—or 0.3%—in May 2014, and then a decrease to 39,000 m2—or 0.1%—in January 2016.


Introduction
Small-scale mines are defined as a mine that operates using rudimentary mining and milling methods (minimal mechanization and minimal, or absence of, structures that control mine and mill spillages) at the level of individuals, families or cooperatives, and temporarily makes a specific mining claim or in other cases may be illegal [1,2].Small-scale mines can be major contributors to local economies and livelihoods [3,4] but in many cases, this comes with significant environmental impacts due to unsuitable locations, rudimentary mining methods, and lack of effective regulation [5,6].The diffuse and remote locations of the mine sites-and their relatively short and unpredictable life span-impede the systematic documentation of small-scale mine development and footprints [7].This is challenging for accurately surveying land use, the accurate attribution of impacts and the targeting of regulatory resources.The difficulty of effectively monitoring small-scale mining using traditional surveys calls for supplementary approaches based on remote sensing.The application of remote sensing from satellites, in conjunction with supervised techniques for classifying land use and land cover (LULC), is attractive for several well-established reasons.Remote sensing provides the opportunity to observe changes in LULC over inaccessible regions.By combining images from different passes of one or more satellites, the evolution of LULC can be modelled.For example, Landsat now provides 45 years of historical records from 1972 [8].Continual improvements in spatial, spectral, radiometric and temporal resolutions are permitting greater levels of information about the type, timing and extent of changes to be obtained.
Small-scale mining, however, introduces some challenges to applying remote sensing and associated land use classification techniques.The surface footprints for smaller mines may only be tens of meters, and less for underground mines or mines in early stages of development, reducing the applicability of Landsat and other freely accessible products with moderate spatial resolutions.Another challenge is that the visible footprint of the mines-mainly vegetation and top-soil clearance -is similar to many other land cover and land use changes such as clearance for agriculture or forestry and is thus likely to have similar spectral responses.Moreover, mining areas can be confused with recently burned forest [9].Finally, small-scale mining is a dynamic land use and mines can be established and abandoned within a period of months therefore they may be difficult to capture with, bi-annual or annual images for example.These challenges are especially problematic in wet tropical regions, where cloud cover is likely to prevent the use of a large proportion of available images, and regular rainfall-driven land-slides may confound the identification of mines.
A key challenge for the remote sensing of mining land cover is that several forms of land use have similar spectral signatures [10].Using pixel-based land cover classification approaches based only on differences in spectral signatures can therefore result in misclassifications.This can be exacerbated if the distributions of reflectance for each class are not normally distributed [11].The limitations of pixel-based classification approaches have given rise to development of object-based segmentation techniques [12], which are collectively known as geographic-object-based image analysis (GEOBIA) [12,13].Using GEOBIA, pixels with homogenous properties are merged to produce multiple segments (objects) where pixels are grouped based on not only their spectral properties but also the spatial and textural characteristics of the segments to which they potentially belong [13].
A recent advancement in GEOBIA classification methods has been the application of Support Vector Machines (SVM) to image objects instead of pixels using a technique known as object-based support vector machine (OB-SVM) [14].Recent examples in the literature show that SVMs used with both pixels [13] and image objects [14,15] can outperform other classification algorithms (e.g., Decision Tree, Maximum Likelihood, Nearest Neighbor and Neural Networks).The SVM classifier is not based directly on differences between the statistical distribution of attributes of separate classes [15,16] but instead it uses non-parametric machine learning algorithms that determine the optimal boundaries among classes [11].The SVM classifier is noteworthy for its ability to separate complex classes with limited training data [17].In addition, the overall accuracy of an OB-SVM can be improved by including other spatial data such as Digital Elevation Model (DEM) and Vegetation Indices (VI) [14].
This research explored the utility of remote sensing for capturing land use change dynamics in wet tropical regions subject to transient small-scale mining using a case-study in the Didipio catchment region, Philippines.We applied the OB-SVM classifier to map land cover change for three different scenes between 2013 and 2016.Imagery was acquired from Pleiades-1A for the 2013 and 2014 scenes and from SPOT-6 for the 2016 scene.Historical spatial data on location and type of ASM mines were also collected from the field and were utilized as training data for the OB-SVM classifier.This was supplemented by manual image interpretation and ancillary spatial data for the other land covers in the region which included forests, rivers and the nearby Didipio large-scale mining operation.We conclude by discussing the utility of the methods and imagery in the context of the remote sensing of ASM in wet tropical regions.

Study Area
The Didipio catchment is situated in the Cagayan region of northeast Philippines.The anthropogenic land use features in this area include open land (deforested area), cultivated areas, small-scale gold mining and large-scale gold mining.Prior to this study, there was no publicly available documentation of land use in this region.Mining has become a principal land use since 2013, when operations started at the large-scale Didipio mine.This is composed of the open-pit mine and other mining structures such as haulage roads, stockpile areas and motor pool open wash bays with an estimated (as of January 2016) surface area of 2.23 km 2 (Figure 1).The mine expanded from phase 1 to phase 2 during the period between 2012 and 2013.The area has seven principal rivers of interest (Figure 1).The anthropogenic use of the Alimit, lower Didipio, Dupit and Surong rivers is primarily to support agriculture, whereas the Camgat, Camgat-Surong, Dinauyan and upper Didipio rivers are available for general community purposes but are also the focus of small-scale mining activity.

Study Area
The Didipio catchment is situated in the Cagayan region of northeast Philippines.The anthropogenic land use features in this area include open land (deforested area), cultivated areas, small-scale gold mining and large-scale gold mining.Prior to this study, there was no publicly available documentation of land use in this region.Mining has become a principal land use since 2013, when operations started at the large-scale Didipio mine.This is composed of the open-pit mine and other mining structures such as haulage roads, stockpile areas and motor pool open wash bays with an estimated (as of January 2016) surface area of 2.23 km 2 (Figure 1).The mine expanded from phase 1 to phase 2 during the period between 2012 and 2013.The area has seven principal rivers of interest (Figure 1).The anthropogenic use of the Alimit, lower Didipio, Dupit and Surong rivers is primarily to support agriculture, whereas the Camgat, Camgat-Surong, Dinauyan and upper Didipio rivers are available for general community purposes but are also the focus of small-scale mining activity.On the other hand, the availability of cloud-free satellite data is generally only possible in the dry season between January and April, although a high percentage of cloud cover persists for most of that period.Landsat-8 Operational Land Imager (OLI), a satellite with a 16-day acquisition frequency, was tested for cloud-free image availability of the catchment between 2013 and 2016.Over the four years, only five to eight of its images were suitable for image analysis, while only two full images of delineated Didipio catchment were unobstructed by atmospheric effects.The remainder of the images are contaminated either by clouds, haze, or shadows.
Since 2012, the large-scale mine operator has kept a database of small-scale mine locations, based on ground surveys every three months.Each site is labelled as either an active or inactive mine.An inactive mine includes operational mines where there happened to be no activity on the day of On the other hand, the availability of cloud-free satellite data is generally only possible in the dry season between January and April, although a high percentage of cloud cover persists for most of that period.Landsat-8 Operational Land Imager (OLI), a satellite with a 16-day acquisition frequency, was tested for cloud-free image availability of the catchment between 2013 and 2016.Over the four years, only five to eight of its images were suitable for image analysis, while only two full images of delineated Didipio catchment were unobstructed by atmospheric effects.The remainder of the images are contaminated either by clouds, haze, or shadows.
Since 2012, the large-scale mine operator has kept a database of small-scale mine locations, based on ground surveys every three months.Each site is labelled as either an active or inactive mine.An inactive mine includes operational mines where there happened to be no activity on the day of survey, as well as seemingly abandoned mines.Supplementary information, such as the mining technique and/or presence of processing equipment and earth-moving vehicles has also been recorded.This database was supplemented by ground surveys by the authors in November 2015 and February 2016.Figure 2 shows the evolution of the location and type of small-scale mines; and photographs for sample sites are presented in Figure 3.
Remote Sens. 2017, 9, 945 4 of 22 survey, as well as seemingly abandoned mines.Supplementary information, such as the mining technique and/or presence of processing equipment and earth-moving vehicles has also been recorded.This database was supplemented by ground surveys by the authors in November 2015 and February 2016.Figure 2 shows the evolution of the location and type of small-scale mines; and photographs for sample sites are presented in Figure 3.In March 2012, the majority of small-scale mining sites were near the rivers.During this period, excavating and panning the river sediments was the principal mining method (Figure 4).Out of the  survey, as well as seemingly abandoned mines.Supplementary information, such as the mining technique and/or presence of processing equipment and earth-moving vehicles has also been recorded.This database was supplemented by ground surveys by the authors in November 2015 and February 2016.Figure 2 shows the evolution of the location and type of small-scale mines; and photographs for sample sites are presented in Figure 3.In March 2012, the majority of small-scale mining sites were near the rivers.During this period, excavating and panning the river sediments was the principal mining method (Figure 4).Out of the In March 2012, the majority of small-scale mining sites were near the rivers.During this period, excavating and panning the river sediments was the principal mining method (Figure 4).Out of the 43 panning sites, 26 were located in the upper and 8 in the lower Didipio River, while 9 were situated in smaller tributary creeks of the Dinauyan River (Figure 1).In addition, small-scale underground mining began to develop with tunnels operating perpendicular to these creeks, which became part of the Didipio mine phase 1 during the same year.
Remote Sens. 2017, 9, 945 5 of 22 43 panning sites, 26 were located in the upper and 8 in the lower Didipio River, while 9 were situated in smaller tributary creeks of the Dinauyan River (Figure 1).In addition, small-scale underground mining began to develop with tunnels operating perpendicular to these creeks, which became part of the Didipio mine phase 1 during the same year.In 2013, hydraulic mining with assistance of hydraulic excavators was introduced in order to increase production yield.These mining practices became dominant in sites operating near the Dinauyan River and near the junction of the Camgat and Surong rivers.The development of smallscale underground mines remained common throughout 2013, with 4 newly constructed tunnels beside the lower Didipio River.All of the surveyed small-scale mines were between 50 and 100 m width in terms of land clearance and between 400 m 2 and 14,000 m 2 surface area.Those mines located on exposed rocks (left of lower Didipio River) appear similar to the open-pit mine in the images, whereas the mines near the rivers had a diverse appearance of sediments and topsoil (Figure 2).
By 2014, the majority of the small-scale mines had moved from their initial locations to higher elevations.Along the Camgat River, newly established larger scale mines (but still within the definition of small-scale mines) were observed.Meanwhile, the development of new underground mines slowed.In 2015, most of the mine tunnels had been abandoned but small-scale mines on the surface, particularly near the Camgat River, were still operating (Figure 5).In 2013, hydraulic mining with assistance of hydraulic excavators was introduced in order to increase production yield.These mining practices became dominant in sites operating near the Dinauyan River and near the junction of the Camgat and Surong rivers.The development of small-scale underground mines remained common throughout 2013, with 4 newly constructed tunnels beside the lower Didipio River.All of the surveyed small-scale mines were between 50 and 100 m width in terms of land clearance and between 400 m 2 and 14,000 m 2 surface area.Those mines located on exposed rocks (left of lower Didipio River) appear similar to the open-pit mine in the images, whereas the mines near the rivers had a diverse appearance of sediments and topsoil (Figure 2).
By 2014, the majority of the small-scale mines had moved from their initial locations to higher elevations.Along the Camgat River, newly established larger scale mines (but still within the definition of small-scale mines) were observed.Meanwhile, the development of new underground mines slowed.In 2015, most of the mine tunnels had been abandoned but small-scale mines on the surface, particularly near the Camgat River, were still operating (Figure 5).

Satellite Imagery
The size of small-scale mines was the primary consideration in selecting the source satellites with visible footprint widths as low as 50 m.Landsat 8 was discounted due to its low spatial resolution (~30 × 30 m pixel size).The availability of satellite images was also a limitation due to cloud cover.Two satellite data sources were selected, which addressed both these limitations-SPOT-6 and Pleiades-1A, with pixel sizes of 2 m and 6 m, respectively (Table 1).The Pleiades and SPOT images originate from pairs of identical satellites (SPOT-6 and 7; and Pleiades-1A and 1B), which revolve in the same orbital plane [18], making it possible to secure coherent images with identical spatial extents every one to two days, increasing the probability of obtaining cloud-free images.In addition, both sources provide three spectral channels in the visible spectrum, a single channel of near-infrared and a panchromatic band.For the Didipio catchment, two images with near-zero cloud cover were captured by Pleiades-1A on 20 March 2013 and 5 May 2014, while one image with near-zero cloud cover was captured by SPOT-6 on 6 January 2016 (Table 1).Out of all the SPOT and Pleiades images from January 2012 to January 2016 with acceptable deviations from nadir, these images were the only ones identified with less than 2% cloud cover.Although images with more than 2% cloud coverage show some small-scale mining areas, the direct cloud areas combined with the cloud shadows were deemed to be too contaminating.Furthermore, the clouds tend to be at higher elevations and therefore the sampling errors introduced are not random.Likewise, Pleiades-1A and SPOT-6 images with low deviation of sensor angle from nadir (e.g., less than 25 degrees) were prioritized as they have less geometric errors to be corrected during orthorectification process [19,20].

Radiometric Calibration of Sensors
The radiometric calibration of the images was performed to: (1) normalize brightness values and (2) convert the data into Top-of-the-Atmosphere (ToA) reflectance.The calibration was conducted in two steps.The original geometrically corrected images provided by the supplier were calibrated to spectral radiance (L λ ) using (Equation ( 1)) [21][22][23] and gain and bias values obtained from the supplementary metadata for each acquisition.The output of Equation (1) for each image was used in Equation (2) to derive the ToA reflectance [21,22].The band-specific exoatmospheric irradiance at a given sun elevation angle (θ SE ) and earth-sun distance (ESUN) during the acquisition of each image were also provided from the metadata.The earth-sun distances were manually calculated using Equation (3) [24].The Julian date in Equation (3) corresponds to the number of days between Universal time on 1 January 4713 BC and the capture time of the satellite image [25].

Atmospheric Correction
An atmospheric correction was then applied to the calibrated satellite images to reduce the effects of atmospheric absorption and scattering [24].Atmospheric correction was applied using the Quick Atmospheric Correction (QUAC) algorithm provided in the ENVI software [26][27][28].QUAC is an approximate, in-scene atmospheric approach that normalizes images based on the statistical properties of object spectra found within an image [29].The QUAC algorithm can be applied when there are at least 10 diverse objects in the image [30].The average reflectance of each object was used as a baseline reference.Prior to correction, objects with high reflectance values such as clouds were manually masked and highly structured regions in the image (e.g., forest areas) were selectively masked by the algorithm.This was followed by the linear transformation of spectral signatures by the QUAC algorithm through plotting their reflectance against a set of universal reference spectra of image objects (e.g., road networks, bare soil, river networks, roofs) from the software's spectral library [30].The resulting outputs from the application of the atmospheric correction were surface reflectance values as opposed to the ToA reflectance values produced in the previous step.

Enhancement of Spatial Resolution through Image Fusion and Supplementary Pre-Processing Techniques
To increase the spatial resolution, image fusion with a pan-sharpening algorithm was applied to reduce the pixel size of the original multispectral images to that of the panchromatic band (0.50 m for Pleiades-1A and 1.5 m for SPOT-6).The nearest-neighbor diffuse (NND) algorithm was used from the ENVI software for this task.Unlike other pan-sharpening algorithms, this method processes spectral bands separately.NND operates on a per pixel basis, thus providing more accurate spectral information, which can improve land cover classification [31].Furthermore, a mosaicing technique was applied to sections of the 2014 Pleiades-1A image with clouds over fixed land covers (e.g., forest).This section of the reference image had been clipped, and the reference image was merged with the cloud-free equivalent of this section from the 2013 Pleiades-1A image.For the concluding steps, a bilinear resampling technique was used for all the images.The spatial resolution of the Pleiades-1A images was reduced to 1 m to remove pixel noise and to increase the classification processing speed.In contrast, the pansharpened SPOT-6 image (1.5 m) was also resampled to 1 m.To smoothen the intensity contrast of image object boundaries, a 3 × 3 pixels median kernel was applied for all the satellite images.The final products were satellite images with 1 m spatial resolutions.

Classification Method
For the image classification, a two-step process was used (Figure 6).The example-based feature extraction workflow by ENVI was used for the application of OB-SVM, with segmentation parameters developed for each image.In the first step, each image was segmented.The scale of segmentation was configured using the edge algorithm.This algorithm locates the outermost edge between image objects, which corresponds to adjacent pixels with the highest intensity contrast [32].The segmentation work flow first creates small objects/segments which are then combined through the application of the full lambda schedule algorithm, which merges objects based on similarity of their spatial and spectral properties [33].Before segmentation, a kernel size of 19 by 19 pixels was configured to calculate the values of the textural attributes of the segments.Since the pixel size of the final image is 1 m, the kernel resembles a 19 × 19 m box or 361 m 2 .This kernel size is sufficient to capture the textural difference (e.g., footprints) of small-scale mines from other land clearances.Although the segmentation parameters were configured in similar ways for all images, the SPOT-6 image had larger and fewer segments than Pleiades-1A due to the difference in spatial resolution between these two satellite sensors.
Remote Sens. 2017, 9, 945 9 of 22 configured to calculate the values of the textural attributes of the segments.Since the pixel size of the final image is 1 m, the kernel resembles a 19 × 19 m box or 361 m 2 .This kernel size is sufficient to capture the textural difference (e.g., footprints) of small-scale mines from other land clearances.
Although the segmentation parameters were configured in similar ways for all images, the SPOT-6 image had larger and fewer segments than Pleiades-1A due to the difference in spatial resolution between these two satellite sensors.After image segmentation, a proportion of segments were manually selected and classified to create training data for OB-SVM (Table 2).We classified 10 land cover types ranging from forest to small-scale mining.The classification was based on the baseline land cover map developed from the ecological assessment of the Didipio Project, the small-scale mining ground-truthing and manual interpretation of the images.Table 2 describes the number of training segments for each land cover type.Approximately, 8% and 10% of the catchment areas were used as training data for the Pleiades-1A and SPOT-6 images, respectively.Due to differences in segment/object sizes for the SPOT-6 image, a lower number of segments was used as training however the total training area was larger at 10% of the catchment.There were also differences in the number of segments used for training data between land cover types.For example, as there are several forms of tree canopies present in the catchment, small segments were generated for forest classes and thus more segments were assigned.In addition, general expansion of various forms of land clearance (e.g., open land, large-scale mine) as well as the tailings storage facility were observed in the 2014 image, therefore more segments were selected to represent these land cover types in the training data.For the small-scale mine (SSM) class, the ground surveys that best coincided with the image acquisition dates were used to develop training data.After image segmentation, a proportion of segments were manually selected and classified to create training data for OB-SVM (Table 2).We classified 10 land cover types ranging from forest to small-scale mining.The classification was based on the baseline land cover map developed from the ecological assessment of the Didipio Project, the small-scale mining ground-truthing and manual interpretation of the images.Table 2 describes the number of training segments for each land cover type.Approximately, 8% and 10% of the catchment areas were used as training data for the Pleiades-1A and SPOT-6 images, respectively.Due to differences in segment/object sizes for the SPOT-6 image, a lower number of segments was used as training however the total training area was larger at 10% of the catchment.There were also differences in the number of segments used for training data between land cover types.For example, as there are several forms of tree canopies present in the catchment, small segments were generated for forest classes and thus more segments were assigned.In addition, general expansion of various forms of land clearance (e.g., open land, large-scale mine) as well as the tailings storage facility were observed in the 2014 image, therefore more segments were selected to represent these land cover types in the training data.For the small-scale mine (SSM) class, the ground surveys that best coincided with the image acquisition dates were used to develop training data.The OB-SVM classifier was then applied to the three scenes and ancillary data using the training data described above.Ancillary data in the form of Interferometric Synthetic Aperture Radar (IFSAR) Digital Elevation Model (DEM) with a vertical resolution of 5 m (resampled to 1 m) and NDVI calculated from the satellite imagery were used.The IFSAR DEM of the catchment was acquired between March and July 2013.Mine pits may be detected through DEMs of sufficient vertical resolution.This may aid the OB-SVM classifier to delineate mining areas from adjacent land covers (e.g., grassland), which have lower and more regular slopes.The DEM was used to form a composite image with the four spectral channels of each image and their respective NDVI maps through image stacking.These multiple spatial data were used as inputs for the calculation of the spatial, spectral and textural attributes of training data.The training data were used by the OB-SVM classifier to generate the thematic map of land use classes.
In the application of an OB-SVM, a range of kernels can be used for separating classes.Known kernels commonly used in remote sensing are polynomial and radial basis functions kernels, which perform well for classes that are nonlinearly separable [11].For this research, a nonlinear radial basis function kernel of gamma (γ) = 0.01 was used for the Pleiades-1A image in 2013 based on previous literature and testing.For example, Huang et al. [11] tested values of γ from 0.1 to 1 and found γ = 0.1 gave the optimum boundaries among classes with less misclassification rate, whereas γ = 0.01 was used by Lin et al. [15] in classifying 10 LULC features.However, a γ of 0.05 was selected for the 2014 Pleiades-1A and 2016 SPOT-6, as our testing found it to have the best result for these images.

Accuracy Assessment
Seventy-five test pixels for each LULC class were sampled from the map of classification results and compared with the same pixels from the baseline assessment map.Stratified sampling was used, as implemented in the ENVI software.This sample size meets the guidelines provided from previous studies [34,35].The training data pixels were excluded during sampling of test pixels.Correct and incorrectly classified pixels were used to generate a confusion matrix for each image.The Overall Accuracy of classification of each image (Equation ( 5), Producer's Accuracy (error of omission; Equation ( 6)) and User's Accuracy (error of commission; Equation (7)) were then calculated [35].
where We also assessed biases in classification.For each class j, the proportion of the total catchment area wrongly classified as i was estimated as P ij (Equation ( 8)).These proportions were used to adjust the original classification results (A Total ) to give an improved estimate of the actual surface area for this class (A j ) (Equation ( 9)) [36].The standard error for the proportion of each class was determined (Equation ( 10)) to generate the specific standard error of every bias-adjusted surface area (Equation ( 11)).Furthermore, to obtain the margin of error of a class, the standard error of a bias-adjusted surface area is multiplied by the z-score of the standard distribution at 95% confidence that is 1.96 or simply equivalent to two.
S A j = S P j •A Total (11) where P ij = weight of misclassified test pixels in each row; t ij = number of misclassified test pixels in each column; t i = number of assigned test pixels in each row; W i = ratio of the calculated class area and total catchment area; A j = Bias-adjusted surface area; S(P j ) = Standard error.

Comparison of Spectral Signatures across LULC Classes
Comparing the spectral signatures for the training data (Figure 7) gives some insight into the challenge of classification.Unlike the vegetated areas (e.g., forest, grassland, cultivated area), which have distinct reflectance values in the near infrared channel (825 nm), several classes (Figure 7 Plot 1B, 2B, 3B) in which soil is exposed have similar reflectance across all spectral channels.These classes include: open land made up of stripping and burning of vegetation; small-scale mines in which deeper soil layers were brought to the surface; roads mainly made up of compacted sand and gravel; and the mine pit benches.Adding to the challenges, the stream network comprises of clear water, highly turbid water and riverbanks, and results in a spectral signature similar to the forms of bare soil.The application of NDVI maps helped separate vegetated areas from bare soil and the variation of reflectance in the near infrared region provides some information on the potential presence of patches of vegetation.For instance, the open land in the 2013 image, despite having the lowest overall reflectance in the visible spectrum, has the highest surface reflectance in the near infrared channel.However, this was unhelpful for distinguishing small-scale mines from other forms of bare soil.
The results in Figure 7, therefore, illustrate that the differences in spectral signatures are small between several classes of land and spectral attributes alone are unlikely to be sufficient to separate the small-scale mines from other LULC classes that are or appear to be stripped land.This illustrates the importance of GEOBIA, which also considers the spatial and textural properties of image objects.

Accuracy Assessment of Classification for Pleiades-1A and SPOT-6 Images
The application of the Pleiades-1A and SPOT-6 images and the OB-SVM classifier to locate small-scale mines had acceptable performance in terms of the Producer's Accuracy (PA) values.The PA values for the SSM class were 90%, 93% and 96% for the 2013, 2014 and 2016 images, respectively (Table 3).The high PA values illustrate that almost all of the small-scale mines were correctly classified.However, lower values of User's Accuracy (UA) were obtained for the SSM class: 76%, 69% and 69% for the 2013, 2014 and 2016 images, respectively (Table 3).The lower UA values illustrate that a significant number of pixels were misclassified as SSM.This was due to the pixels classified as SSM which included irregular segments of vacant house lots, large-scale mining, cultivated areas and open land.Still, the overall accuracy of the OB-SVM classifier applied to these images can be considered satisfactory, with OA values of 89%, 87% and 89% for the 2013, 2014 and 2016 images, respectively.
The misclassification of other LULCs as SSM may be attributed to several factors.Firstly, Table 3 shows that 4% to 13% of test pixels in the 2013, 2014 and 2016 images that were actually large-scale mining were misclassified as SSM.This is likely to be due to spectral and textural similarities of the SSM class with the construction of the in-pit ramps and benches between 2013 and 2016.The application of the IFSAR DEM as ancillary data during image classification only contributed in discriminating the sections of large-scale mine that were already present at the time of acquisition of the IFSAR DEM.The IFSAR DEM provided the information on the ground depression difference between the large-scale mine and the small-scale mines as well as their distribution throughout the catchment.The large-scale open pit mine was also among the well-classified LULCs, with UA values of 95%, 97% and 100%.Its size, unique shape and textural consistency generated distinguishable segments that were easily recognized by the OB-SVM classifier.
Moreover, in the 2013 and 2016 images, 8% of the cultivated area test pixels were misclassified as SSM.In SSM, stripping of ground cover is selective, and difficult-to-access ore bodies are left partially stripped and less undisturbed.On several occasions, the extraction of only a single section of the ore was observed over an extended period.Therefore, although deeper soil layers are exposed, these areas are often appearing as superficial disturbance that is easily misclassified as cultivated area or open land.In addition, both SSM and open land classes include irregularly shape segments.This kind of misclassification was also encountered in 2014 image, but with lower frequency.This is likely to be because the higher spatial resolution of this image allowed the textural features of small-scale mines to be better distinguished.
The resulting UA for the SSM class may be improved if sufficient cloud-free images could be obtained to capture the growing season, when cultivated areas and open lands have minimal stripping activities and therefore are less likely to be misclassified as SSM.The efficacy of the OB-SVM classifier may be examined in other catchments less affected by cloud cover.In addition, the Didipio catchment was challenging due to presence of several types of SSM with limited coverage.The spatial distribution of the ten LULC classes is presented in Figure 8 with areas listed in Table 4.Although there are misclassifications among the classes (Table 4)-if it is assumed that biases identified for the test pixels represent the whole catchment-these classified areas can be bias-adjusted with given margins of error (Table 4) [36].The areas referred to in the discussion below are the bias-adjusted values.Forest was the dominant natural landscape and, according to the OB-SVM classifier, it occupied 60% of the total catchment on average over the three years, mostly in the outer, higher altitude regions of the catchment.Meanwhile, grasslands were classified as being the next most dominant, representing an average of 23% of the area, concentrated at lower altitudes surrounding the large-scale mine.For stripped areas, the dominant LULCs according to the classifier were open land, large-scale mining and the road system/vacant lots, with average areas of 1.94, 1.95 and 0.68 km 2 .The cultivated areas were constant at approximately 1.3 km 2 .The area classified as SSM was on average 0.098 km 2 representing only 0.25% of the total catchment area, the smallest LULC footprint among the ten LULC classes.However, the small-scale mining activity may involve, as well as the mines themselves, LULCs included under the open land, road system/vacant lot classes.
Over the three years, there were significant changes in the mining activity.The areas classified as large-scale mine and the tailings storage facility (TSF) increased by 0.77 km 2 and 0.36 km 2 , respectively.The area classified as SSM increased slightly from 0.091 to 0.121 km 2 between the 2013 and 2014 images.However, there was a decrease in 2016 by 0.039 km 2 .Based on that result, the estimated area coverages of SSM were 0.091 km 2 , 0.121 km 2 and 0.082 km 2 for the 2013, 2014 and 2016 images.The decline in small-scale mining in the 2016 image is counter-intuitive if it is considered that all abandoned mines should be classed as SSM.However, abandoned mines that have revegetated were not included in the ground truth surveys, and if they were included in the training data, they will be recognized under other classes.Uncertainty in the change in SSM coverage, based on the estimated margin of errors (Table 4), is relatively high and it is possible that the change from 2013 to 2014 is an artefact of misclassification errors.If it is real, this temporal change in SSM class coverage can be said to be significant at local scale only.
The generated thematic maps, apart from evaluating the SSM coverage, show how sediment pits and panning areas are more likely to develop once active whereas underground mines have more static footprints.In addition, these maps also permit judgment about the type of small-scale mining and how it varies throughout the catchment.For instance, the proximity to a river of areas classed as SSM as well as the presence of small ponds within these areas indicate sediment pits (Figure 5).From this, it is inferred that the small-scale mines neighboring the Camgat River, and those situated north of Camgat-Surong River besides the mine facilities/camp site are mostly sediment pits.This is supported by the ground surveys, where pumps and pipes were observed in these areas.Areas in the upper and lower Didipio River that were classified as SSM were identified in the ground surveys as well as in the images as scoured river banks associated with panning.In contrast, areas classed as SSM to the south of the upper Didipio River (Figure 2c), which were ground-truthed as tunnels, are almost 800 m away from the nearest river, and a decreasing surface area was observed within three years.Despite the identified dynamics, overall the areas classified as SSM remained in the same parts of the catchment.It may be concluded that high resolution satellite images may be useful, not only for identifying areas that are likely to include small-scale mines, but also for distinguishing between the main types of small-scale mine.

Figure 1 .
Figure 1.Major features of the Didipio catchment area in the Philippines.The land cover is shown only within the Didipio catchment area, while the extent of the figure corresponds to the boundary of the Pleiades-1A and SPOT-6 images.

Figure 1 .
Figure 1.Major features of the Didipio catchment area in the Philippines.The land cover is shown only within the Didipio catchment area, while the extent of the figure corresponds to the boundary of the Pleiades-1A and SPOT-6 images.

Figure 4 .
Figure 4. Active small-scale mining sites derived from field data from 2012 to 2015 categorized by mining method.

Figure 4 .
Figure 4. Active small-scale mining sites derived from field data from 2012 to 2015 categorized by mining method.

Figure 5 .
Figure 5. Locations and types of active small-scale mine sites based on ground surveys, superimposed on remote-sensed images.2010 Satellite Image courtesy of the DigitalGlobe Foundation; Pleiades: © CNES (2013, 2014), Distribution Airbus DS/Spot Image; SPOT: © AIRBUS DS (2016).Note: 2010 data used due to lack of alternative cloud-free image in 2012.

Figure 5 .
Figure 5. Locations and types of active small-scale mine sites based on ground surveys, superimposed on remote-sensed images.2010 Satellite Image courtesy of the DigitalGlobe Foundation; Pleiades: © CNES (2013, 2014), Distribution Airbus DS/Spot Image; SPOT: © AIRBUS DS (2016).Note: 2010 data used due to lack of alternative cloud-free image in 2012.

Figure 6 .
Figure 6.Flowchart of image pre-processing and image analysis (classification and post-classification).

Figure 6 .
Figure 6.Flowchart of image pre-processing and image analysis (classification and post-classification).
OA = Overall accuracy; PA = producer's accuracy; UA = User's accuracy; n = Number of classes; R i = Number of correctly classified test pixels in class i; N i = Number of test pixels that are classified to be in class i (=75); M i = Number of test pixels that are in class I in the baseline map.

Figure 7 . 6 Figure 7 .
Figure 7. Plots with suffixes A, B and C describe variation in spectral signatures of vegetated areas, stripped lands and water bodies, respectively, from the three different images where prefixes 1, 2 and 3 correspond to 2013 Pleiades-1A, 2014 Pleiades-1A and 2016 SPOT-6 images.

Table 1 .
Acquisition dates and characteristics of the three scenes(2013, 2014, 2016)and corresponding sensor characteristics.

Table 2 .
Descriptions of the land use and land cover classes and size of training data in each image.

Table 3 .
Cont.Coverage and Spatial Distribution of Areas Classed as SSM Relative to other LULC Classes

Table 4 .
Classified areas and Bias-adjusted areas of each LULC class.