Mapping Potentially Acid Generating Material on Abandoned Mine Lands Using Remotely Piloted Aerial Systems

: Weathering and transport of potentially acid generating material (PAGM) at abandoned mines can degrade downstream environments and contaminate water resources. Monitoring the thousands of abandoned mine lands (AMLs) for exposed PAGM using ﬁeld surveys is time intensive. Here, we explore the use of Remotely Piloted Aerial Systems (RPASs) as a complementary remote sensing platform to map the spatial and temporal changes of PAGM across a mine waste rock pile on an AML. We focus on testing the ability of established supervised and unsupervised classiﬁcation algorithms to map PAGM on imagery with very high spatial resolution, but low spectral sampling. At the Perry Canyon, NV, USA AML, we carried out six ﬂights over a 29-month period, using a RPAS equipped with a 5-band multispectral sensor measuring in the visible to near infrared (400–1000 nm). We built six different 3 cm resolution orthorectiﬁed reﬂectance maps, and our tests using supervised and unsupervised classiﬁcations revealed beneﬁts to each approach. Supervised classiﬁcation schemes allowed accurate mapping of classes that lacked published spectral libraries, such as acid mine drainage (AMD) and efﬂorescent mineral salts (EMS). The unsupervised method produced similar maps of PAGM, as compared to supervised schemes, but with little user input. Our classiﬁed multi-temporal maps, validated with multiple ﬁeld and lab-based methods, revealed persistent and slowly growing ‘hotspots’ of jarosite on the mine waste rock pile, whereas EMS exhibit more rapid ﬂuctuations in extent. The mapping methods we detail for a RPAS carrying a broadband multispectral sensor can be applied extensively to AMLs. Our methods show promise to increase the spatial and temporal coverage of accurate maps critical for environmental monitoring and reclamation efforts over AMLs.


Introduction
Abandoned Mine Lands (AMLs) can adversely impact the environmental quality of watersheds in various ways. When sulfide minerals, such as pyrite, are brought to the surface during mining, they can react with air and water. Sulfuric acid and ferrous iron can form to create acid mine drainage (AMD), which can contaminate nearby surface and groundwater. If AMLs and their associated waste rock piles are continually exposed to weathering, large concentrations of sulfates and heavy metals can be found in surface and groundwater for decades after mining has ceased [1,2].
Containment strategies for mine waste rock or other solid mining wastes, include capping the waste rock with a soil or engineered cover. These can avert exposure and erosion of ( Figure 1). Perry Canyon hosts a system of historical mines, some of which are adjacent to Perry Creek (Figure 1b), an ephemeral creek with multiple springs and mine water seeps. Mining in Perry Canyon began in the late 1800's. The Jones-Kincaid mine, our primary focus, was developed in the 1870's, with the main commodities being silver and copper [18]. The mine consisted of a 300 m adit, 150 m deep shaft, and a mine waste rock pile (Figure 1b) [18].
We chose to monitor AMD and PAGM at the Perry Canyon AML because it is a site with problems characteristic of those present at many of the thousands of waste rock piles and adits in the Western United States [19] (Figure 1a). It has been established by the Bureau of Land Management (BLM) and the University of Nevada, Reno as an important site for testing and evaluating various monitoring, evaluation, and treatment methods for AMLs [20][21][22]. In 1995, BLM and the Nevada Bureau of Mines and Geology (NBMG) identified the Jones-Kincaid Mine and the associated mine waste rock pile (Figure 1b), as being potentially hazardous due to the presence of AMD [23]. Specifically, waters draining from the abandoned mine were among the most acidic in the state, with elevated levels of toxic metal concentrations 2300 times the Federal safe drinking water standard [23,24]. The 2004 Final Perry Canyon Site Investigation Report [17] confirmed the presence of AMD flowing from the adit seep at the Jones-Kincaid mine waste rock pile. Efflorescent mineral salts (EMS) sampled at the adit seep and downstream of the waste rock, contained various concentrations of heavy metals [17]. Samples from the mine waste rock pile revealed sulfate concentrations greater than 8000 mg L −1 and various metal leachate concentrations, with the most elevated being copper, zinc, and iron. The acid generating potential was greater than 100 t CaC0 3 Kt −1 , with no acid neutralization potential detected [17]. In September 2006, the Universal Environmental Nevada (UEN) consultants completed the remediation activities on the Jones-Kincaid mine waste rock pile [25]. The completed remediation included capturing the mine waters to redirect them away from the waste pile via a drainage pipe and adit backfilling, re-contouring and capping of the waste pile, and seeding for re-vegetation of disturbed areas [25]. An approximately 7 cm thick soil cover was applied over the mine waste rock to limit erosion of waste rock, infiltration, and AMD formation via leaching [17,25]. Despite the completed remediation, seasonal AMD discharge continues, and EMS still form at various locations across the site. Erosion of the Jones-Kincaid waste rock cover has occurred since its installation [20], exposing 'hotspots' of yellow and red soil, and red coated minerals. Previous site characterization data, summarized below in the section on field validation activities, indicate that these 'hotspots' are secondary sulfate minerals, including jarosite, and iron oxide rich soils. All together these site characteristics and previous work make for a unique opportunity to test the efficacy of mapping and classifying AMD and PAGM with a broadband sensor atop a RPAS, using common remote sensing techniques.

Materials and Methods
We mapped the temporal change of PAGM on the Jones-Kincaid mine waste rock pile with reflectance maps, and three remote sensing techniques. To collect needed imagery, we deployed a lightweight 5-band multispectral sensor (450-850 nm) via RPAS, in six flights over 29 months. We tested multiple classification methods to identify surface materials of interest, including Maximum Likelihood (ML), Spectral Angle Mapper (SAM), and Band Math Ratios (BMRs). To validate the remote mapping along select transects, we (1) collected soil samples and analyzed them for general soil chemistry and mineral identification using XRD, (2) imaged the transects at mm-resolution with a hyperspectral sensor (400-1000 nm) on a custom-built motion-control system, (3) collected high resolution in silico field spectra with a hand held spectrometer (350-2500 nm).

Mapping Targets and Key Indicators of Potentially Acid Generating Material
Our study focused on classifying, mapping, and validating surface materials and AMD on the Jones-Kincaid mine waste rock pile. The surface materials included yellow soil (jarosite), white and blue soil (EMS), and red soils (iron oxides and other secondary iron rich coatings indicative of AMD). Yellow mine-rich soils are typically associated with jarosite [26], a secondary iron-bearing sulfate mineral which forms under acidic conditions. Jarosite has been used in many studies of active and inactive mines to locate areas of PAGM [16,26,27]. Red soils are typically associated with iron oxides (ferric oxides), and have the ability to sequester heavy metals and generate more acidic water [28,29]. Iron rich minerals have unique and identifiable spectral reflectance signatures in the visible to near-infrared (VNIR) [16], making them good candidates to identify areas of PAGM for this study. Jarosite has an absorption feature at approximately 450 nm and a downturn at 700 nm [30], and iron oxides have an upturn with a steep slope from 550 to 700 nm [31].

RPAS Surveys: Micasense 5-Band Multispectral Sensor
To map the temporal dispersion of AMD discharge into Perry Creek, and the time variable distribution of jarosite and iron oxides at the Jones-Kincaid mine waste rock pile, we used a Tarot 650 RPAS equipped with a 5-Band MicaSense RedEdge sensor. The five bands ranged in wavelengths from 475 to 840 nm (Table 1). We surveyed over a 29-month period, and survey dates included June 2018, September 2018, November 2018, March 2019, July 2019, and October 2020. Universal Ground Control Station (UGCS), version 3.0.1302 software (SPH Engineering, Riga, Latvia), was used to plan and execute programmed missions, flown at 80 m above ground level (AGL) with 80% front and side overlap. The missions were flown with a flight speed of 5 m s −1 with an orientation parallel to topographic contours. GPS measurements were recorded with a Trimble RTX R10 receiver (Trimble Inc., Sunnyvale, CA, USA) with better than 5 cm vertical accuracy. Up to 15 Ground Control Points (GCPs) were established that covered both the canyon bottom and slope up to the ridge line. For surveys which covered only the canyon bottom, approximately eight GCPs were necessary for accurate registration. The analysis presented herein focused on subsets of each time series that match the June 18 survey, which encompassed the canyon bottom and waste rock pile. Sensor specifications are shown in Table 2.
We aligned the visual images obtained during the RPAS flights using Structure from Motion (SfM) in Pix4D v.3.2 (Pix4D, Lausanne, Switzerland) [32]. Geometry of the scene, camera positions, and orientation were solved simultaneously using an iterative bundle adjustment procedure to align overlapping images based on the automatic extraction of features [33]. We manually identified GCPs on multiple images within each survey, which enabled an assessment of any non-linear structural errors in the SfM reconstruction. We used the aligned images to produce digital surface models (DSMs). From here, built-in software functionality converted raw data to reflectance by taking advantage of white reference calibration panel images taken in the field during each survey and measurements from a Downwelling Light Sensor (DLS) mounted on top of the drone. DLS data were recorded in the header file of each image. Utliizing the DSMs and reflectance calibration, we created reflectance maps for each band and each survey, where the value of each pixel indicates the reflectance of the object. Reflectance maps for each band were combined into a single file in ascending band order, from Blue to NIR (Figure 2), for each survey period. Processing details for each flight are presented in Table 3.   Abbreviations: GSD = ground sampling distance, GCPs = ground control points, MRE = mean reprojection error.

Nano-Hyperspectral 271-Band Sensor Transects
Hyperspectral data were acquired on 9 October 2018, during cloud-free conditions. Based on field inspection of the waste pile cover we picked four transects to sample a wide range of surface types. Transect 1 and 2 were just downstream from the adit (Figure 3). Transect 3 was at a higher elevation on the center of the waste pile and Transect 4 was at the far western end ( Figure 3). We used a Nano VNIR Hyperspec camera from Headwall Photonics, MA, USA, which measured reflected light in the 400-1000 nm region. The camera was fitted with an 8 mm lens with a 30.4 degree field of view (FOV). It captured data in 271 bands with 12-bit radiometric resolution with a full width half maximum (FWHM) of 6 nm. The camera was a push broom sensor, thus in order to acquire spatial data, the camera needed to be moved in a smooth and controlled manner. Therefore, data were captured using a custom-built motioncontrol system [34] (Figure 3a) that points the instrument at nadir with an above ground height of 2.5 m. The exposure time was set such that the camera did not get oversaturated while imaging the 99%-reflectance standard (Labsphere Inc., North Sutton, NH, USA). The velocity of the camera moving on the system is a function of the distance between the camera and the ground (height) and the frame rate. The velocity and frame rate were chosen such that the resulting pixels of the image were square. The resulting images were 1.3 m wide, 5.5 m long, with a pixel scale of 2 mm (Figure 3, Transects 1-4). Raw data were converted to radiance in SpectralView (version 5.5.1., Headwall Photonics, MA, USA) using both a dark reference image to account for background noise and the sensor-specific calibration. Next, we converted the radiance to reflectance using the reflectance standard and the empirical line method [35]. We used five to six GCPs per transect, and GPS measurements were recorded with the same Trimble RTX R10 receiver used in the RPAS surveys. The reflectance cubes were then orthorectified in Arcmap 10.5.1 (ArcGIS, ESRI Inc, Redlands, CA, USA) with a first order polynomial fit and nearest-neighbor resampling. The final image cube was a 2 mm per pixel reflectance map ( Figure 3, Transects 1-4).
In order to improve the signal-to-noise ratio of the hyperspectral data, we re-sampled the 271 band data to 45 bands using a spectral binning algorithm that averages 5 spectral channels into one bin ( Figure 4). This reduced the FWHM but did not impact our ability to map the AMD materials of interest. We truncated the spectra at 850 nm, such that the final spectra also matched the wavelengths of the 5-band MicaSense multispectral sensor. For the validation hyperspectral transect data, given the larger number of channels and better definition of spectral shape, we used spectral angle mapper.

Mapping Methods
We built reflectance maps with high spatial resolution and relatively low spectral resolution, to explore how accurately and by which classification algorithm PAGM or AMD can be mapped. We compared two remote sensing supervised classification methods for mapping surface materials, Maximum Likelihood (ML) [36] and Spectral Angle Mapper (SAM) [37]. We then tested a rapid, unsupervised classification algorithm (Band Math Ratios (BMRs)), that has been used to map AMD with success [15]. Both ML and SAM required determination of initial classes or training areas. Our classes of interest included AMD, jarosite, iron oxides, and EMS. We included a soil training area in order to have a control class. We based the selection of training areas on initial site characterization from 2004, and field observation data from the RPAS flights, which confirmed the presence of AMD, iron oxides (red soils), jarosite (yellow soils) and EMSs. We explored two different Band Math Ratios (BMRs) to distinguish jarosite from iron oxides. For each image we masked vegetation using a Normalized Difference Vegetation Index (NDVI) and in some instances we restricted scene analysis to just the area of the mine waste rock pile.

Maximum Likelihood Classification
We first mapped temporal changes using the Maximum Likelihood classification (ML). We classified AMD, soil, jarosite, iron oxides, and EMS on each multispectral survey. ML is a supervised classification method based on the Bayes Theorem [36], where a discriminant function is used to assign a pixel to the class with the highest likelihood.
where i is the class, x = n-dimensional data (where n is the number of bands) and p(ω i ) = probability that class ω i occurs in the image and is assumed the same for all classes, | ∑ i | is the determinant of the covariance matrix of the data in class ω i , ∑ −1 i is its inverse matrix, and m i is the mean vector [36].
We defined separate training data for each individual survey using defined polygons which visually represented each class based on distinctive colors. In our case, yellow soil represented jarosite, red soil and red rocks represented iron oxides, orange colored drainage from the adit represented AMD, and white and blue soils in and around the adit and downstream of the Jones-Kincaid represented EMS. ML utilized the training data to estimate the means and variances of each class to estimate the probability and a frequency distribution for pixels with different training classes. A zone of overlap depicted the pixel values which were common to a specific class, and the pixel with the largest zone of overlap was chosen for the classification. We assigned different thresholds for each class and each survey, which are outlined in Table 4.
We attained an accuracy assessment post classification by assigning 50 random pixels to each class (predicted values), for every classified survey. We then manually inspected each random pixel and visually confirmed whether the point was correctly classified. We amended inaccurate classifications by manually assigning the correct classification to the random pixel (truth values), resulting in two datasets which consisted of predicted values and truth values. We then used this data in ENVI (version 5.5) to return confusion matrices computed from the input truth values and predicted values. The confusion matrices were used to calculate the overall accuracy (%) and the kappa coefficient. The kappa coefficient (κ) measured the agreement between predicted values and truth values, where 1 represents perfect agreement, and 0 represents no agreement. The kappa coefficient is computed as follows: where i is the class number, N is the classified values compared to truth values, m i,i is the number of values which belong to the truth class i that have also been classified as i, C i is the total number of predicted values in class i, and G i is the total number of truth values in class i. Table 4. Thresholds for Maximum Likelihood, where a -indicates that no threshold was used.

Spectral Angle Mapper
We applied Spectral Angle Mapper (SAM) [37] to classify soil, jarosite, and iron oxides, and utilized existing endmember reference spectra derived from spectral libraries sourced from the United States Geological Survey (USGS) database [31] in ENVI (version 5.5) (Figure 5a). The USGS endmember spectrum for iron oxides was an acid mine drainage assemblage, which is a mixture of secondary iron mineral precipitates coating a rock surface [31]. Taking this into consideration, we combined the AMD and iron oxide (AMD\iron oxide) into one class for SAM and labeled the class as "iron oxide". The USGS endmember reference spectra were both higher spectral resolution and covered a longer wavelength range than either the multispectral or hyperspectral sensors. We used ENVI to resample the USGS data in to the same 45 bands ( Figure 4) as shown in (Figure 5b) or 5 multispectral channels (Figure 5c). We provided the wavelengths and FWHM values for each band, used a Guassian model with FWHM spacings to resample the spectra of the USGS endmembers. We combined the four Jarosite species from the USGS into one jarosite class (Figure 5b,c). We removed the EMS from the SAM classification because there is no appropriate USGS reference spectrum for this particular class. Here we restricted scene analysis to just the area of the waste rock pile and part of Perry Creek. These scene-derived endmembers show less pronounced shapes than the library minerals, but exhibit similar trends to the re-sampled USGS data in (c). Iron oxide has an upturn and steep slope from 550 to 650 nm, while jarosite has a downturn from 700 to 850 nm. The shape of the AMD profile exhibits both trends. SAM classification is based on a comparison of the spectral image with an endmember reference spectrum. SAM treats the reference spectrum and the image spectrum as vectors, where dimensionality is equal to the number of bands, and calculates the inner (dot) product between these vectors which results in an angle.
where α is the angle between the reference and the image spectrum, X is the image spectrum, and Y is the reference spectrum. The angle α presents a variation anywhere between 0-90 degrees, and smaller angles represent closer matches to the reference spectrum.
We applied the same threshold angle of 0.1 for all surveys, and as a result, any unclassified pixel had a threshold angle greater than 0.1. The result was a classified image which shows the best match to a given reflectance material for each pixel, in each survey. We followed similar post classification steps as outlined in the previous section using a unified threshold angle of 0.1 for the SAM classification. Each random point was manually inspected using the spectral profiles to confirm whether the point was correctly classified.

Band Math Ratios
For our final mapping method, we implemented band math ratios (BMRs) to classify jarosite and iron oxides on similar scene areas used for SAM. BMRs enhanced the spectral differences between bands by dividing one spectral band by another, producing an image that provided relative band intensities. We used both the reference spectrum derived from the ML polygons and the USGS profiles to choose our band ratios (Figure 5c,d). We applied red green ( 668 560 nm) to capture the upturn representing iron oxides, and red N IR ( 668 840 nm) to capture the downturn representing jarosite. The sharp downturn for jarosite in the USGS profiles began at 725 nm (Figure 5c), whereas the weaker downturn for jarosite in the ML polygons began at 700 nm ( Figure 5d). The ML spectra were derived using defined polygons which visually represented each class based on distinctive colors, whereas the USGS spectra were derived from pure minerals. As a result, the ML polygons did not contain pure jarosite minerals, resulting in weaker features. We applied both rededge N IR ( 725 840 ) and red N IR ( 668 840 ) to map jarosite, and we ultimately found no difference between the two ratios. The downturn of the AMD profile matched that of the jarosite downturn, and the upturn matched that of the iron oxide profiles (Figure 5c,d). AMD was only present in a small number of pixels in three surveys along the adit drain, so we used the classes which have the most pixels, jarosite and iron oxides. We masked vegetation in the images using the Normalized Difference Vegetation Index (NDVI) and threshold regions of interest (ROIs) in ENVI. We then created ROIs using two different unified thresholds for jarosite and iron oxides to represent the mapped classes. We applied a threshold of 1.05-1.1 for jarosite, and 1.7-2.3 for iron oxides on each survey.

Validation: Field Data and Soil Sample Analysis
We collected seven samples from the four hyperspectral transects previously described in Section 2.4. The samples were outsourced to AuTec Innovative Extractive Solutions Ltd., Vancouver, BC, Canada, where they were analyzed with XRD, and peaks were determined using HighScore Plus software by Panalytical withthe the Crystallography Open Database (COD) [38]. Sample preparation and XRD details are listed in Appendix A (Tables A1 and A2).
Based on our study of the six RPAS flights we selected six locations for measurement using a portable field spectrometer and additional sample collection and analysis. Sample locations were selected based on the presence of AMD materials identified in the RPAS surveys. During the O20 RPAS flight we used a Spectral Evolution RS-5400 Portable Spectroradiometer (350-2500 nm), using both the 5 degree FOV lens and the sun as a source and the contact probe with internal halogen source. For interpretation of the field spectra we used standard spectral libraries available in ENVI and from the USGS [31].
Due to the persistent mapping of jarosite in one location, we dug a 0.3 × 0.3 m trench in order to understand whether this was forming at the surface or continuously exposing a jarosite hot spot. We collected yellow soil samples on the surface, at 0.2 and 0.3 m depth. We dug a similar trench 3.6 m west of this location and collected brown soil samples at 0.2 and 0.27 m depth. We also collected blue EMS samples and red soil samples at other locations. The XRD instrument and COD database [38] was available as part of the Chemistry Shared Instrument Lab at UNR. XRD details on these samples are listed in Appendix A (Table A2).

Results
We first discuss classification results for the RPAS surveys for each time increment and each classification method in the 29-month period of observation. We then outline our validation results from XRD on collected soil samples, the October 2018 hyperspectral rail surveys, and the 2020 field studies. We identified several locations where jarosite was mapped persistently, including Transects 1-3 (1,2,3), two zones of yellow soil (Y, site of the trench, and T2), and southeast of the adit (E). The trench site adjacent to Y with no yellow soil was labeled S. Iron oxides were consistently mapped at site A, and the EMS and red soil at B and E, respectively, were only mapped with the ML (as described in Section 2.4.1). The sample locations can be found on the June 18 maps in Figures 6-8, and in detail on Figures 9 and 10.

Maximum Likelihood
The overall accuracy for the total 29-month period using Maximum Likelihood was 72.5%, with a Kappa Coefficient of 0.6502 ( Figure 6). As noted in the individual surveys ( Figure 6), accuracies ranged from a low of 61% to a high of 87%. Jarosite hotspots were consistently mapped in Y, T2, E and  (Tables A3-A8).

Spectral Angle Mapper
The overall accuracy for the total 29-month period using spectral angle mapper was 80%, with a Kappa Coefficient of 0.7104 and ranges from 70% to 94% accuracy (Figure 7). Jarosite hotspots were again consistently mapped in Y, T2, E and

Band Math Ratios
Because band math ratios is a thresholded value, rather than a pixel classification, we could not perform an accuracy assessment similar to ML or SAM. However, we visually saw the same jarosite and iron oxide patterns from the previous classification methods at Y, T2, E, A, and Transects 1-3 (1,2,3) (Figure 8). BMRs did not map iron oxides on the west end of E in any survey, and AMD was mapped as jarosite in September 2018, March 2019, and July 2019.

Validation: X-ray Diffraction Soil Analysis
Soil samples were collected on the Jones-Kincaid mine waste pile in October 2018 and again in October 2020. XRD analysis on the October 18 soil samples detected varying modal abundances of jarosite at 58, 19.1, and 14.6 wt% at Transects 1-3, respectively, and a minor abundance of pyrite at 0.6 wt% at Transect 2 (Table A1). Trace phase grouping definitions, mineral identification results, and the ideal chemical formulas are listed in Appendix A (Table A1). The XRD ground surveys in October 20 were used to confirm the presence of jarosite and iron oxides in areas where they were mapped. Although various mineral assemblages exist within these soil samples, we focused solely on the presence of jarosite and iron oxides. Jarosite was detected in the red soil at E, and at Y on the surface, and at 0.2 m and 0.3 m. We did not detect any jarosite peaks in S, and hematite, an iron oxide, was detected in the red soil at E. XRD patterns for all October 20 samples had three strong peaks which were indicative of quartz. XRD profiles and 2θ peak values for all soil samples are shown in Appendix D ( Figure A1). Unfortunately, pentahydrate was not in the COD database, and therefore we relied on the spectral profiles taken in the field for identification.

Validation: Hyperspectral Transects
We performed validation of the September 2018 RPAS survey using the SAM classification on the hyperspectral data from Transects 1-3. We chose Transects 1-3 because they had confirmed modal abundances of jarosite (Table A1), and we validated the September 2018 RPAS survey because it was the closest survey in time to the October 2018 hyperspectral data collection. We represented the validation using the September 2018 BMR image, although each classification analysis mapped jarosite in these three transects, not only in September 2018, but in every other multispectral survey (Figure 9). By representing this one RPAS survey we then confirmed the locations in all surveys.

Validation: Field Portable Spectroradiometer
Additional validation included a field portable remote sensing spectroradiometer at locations Y, T2, B, A, and E ( Figure 10). These data exhibited the jarosite absorption feature at 450 nm in all areas with yellow soil (Y, T2, and E) ( Figure 11). The spectral profiles for the three samples taken at A showed an upturn and steep slope from 500 to 700 nm, indicative of an iron oxide ( Figure 11). The red soil at E displayed an absorption feature at 450 nm and a bumpy slope from 550 to 700 nm (Figure 11), which is similar to the USGS spectra of an iron oxide paint pigment [31]. The blue EMS at B had a wide absorption feature at approximately 800 nm (Figure 11), indicative of the copper sulfate salt pentahydrate [39]. Pentahydrate is a common bright blue copper sulfate salt [40] which forms as an EMS from waters containing iron and sulfate, the products of pyrite oxidation [2].

Discussion
Our goal was to monitor PAGM and AMD at the Perry Canyon, NV AML using an RPAS equipped with a 5-band multispectral sensor. We investigated two supervised classification methods, maximum likelihood (ML) and spectral angle mapper (SAM), and one unsupervised method, band math ratios (BMRs). We tested how accurately and efficiently each method could be used to classify and map PAGM and AMD. To our knowledge, these techniques have not been used in conjunction with high spatial resolution and low spectral sampling reflectance maps to monitor AMLs. We then used various validation methods to test the accuracy of the established supervised and unsupervised classification algorithms. We created multi-temporal maps of PAGM and AMD, and make recommendations for future use of these techniques on AMLs. These methods are not a replacement for ground measurements; however, we showed that a RPAS combined with broadband multi-channel sensors have the potential to monitor high priority AMLs and guide in silico work.

Maximum Likelihood
We tested three different classification methods, with the first being Maximum Likelihood (ML). Although ML is a commonly used classification scheme, for this study we found implementation across multiple surveys to be time consuming with a high degree of inaccuracy. We were unable to use the same training data for each flight such that training areas had to be independently selected for each survey. We also had to use different thresholds for each class because ML classifies every pixel in the image, so thresholds were set to minimize classification inaccuracies. For example, ML over-classified iron oxides in June 18, even with a threshold of 0.3, and only classified them to greater than 70% accuracy in September and November 2018 because the threshold was 0.3. The advantage of ML is that we were able to map AMD and efflorescent mineral salts (EMS) as their own classes. However, ML only classified AMD in March 19, and did not classify it in September 2018 and July 2019, although it was present. This is an indication that it may be difficult to map AMD dispersion from small adit seeps as its own class, even when using cm-scale reflectance maps. EMS was over-classified in July 2019, March 2020 and October 2020 without using a threshold ( Figure 6). In general, EMS are white in the reflectance images due to the salt content ( Figure 2). However, flying in very sunny conditions, as was the case with most of our surveys, also produces high albedo areas in the data. These high albedo areas are evident in the true color images of each survey ( Figure 2) and could be a reason EMS was consistently misclassified. This is an indication that it may be difficult to distinguish between efflorescent mineral salts and high albedo areas when using a broadband sensor in conjunction with the ML pixel classification. It is unlikely that applying a threshold to this class would have mitigated this issue. For the ML classification, we did not mask all surfaces outside of the Jones-Kincaid waste rock pile and Perry Creek. Although this may have helped increase our accuracy, we do not recommend maximum likelihood due to the time commitment and complicated implementation of this method across multiple surveys.

Spectral Angle Mapper
The second supervised classification method we tested was spectral angle mapper (SAM), which performed better and was easier to implement relative to ML. Utliizing the USGS database for SAM was quite simple for common acid generating materials such as iron oxides and jarosite, and a unified threshold (0.1) made the workflow seamless. The use of different thresholds might have provided more accurate maps, but at the loss of consistency across time. The major disadvantage to this method was that we could not map AMD as its own class, nor did we attempt to map EMS. With regard to EMS, there was not an appropriate mineral reference in the USGS database. With regard to AMD, there were not large regions of this material, although combining the AMD and iron oxide classes still captured the regions with potential acid generating materials (PAGM). This is evident in March 2019 and September 2018, where Iron oxides\AMD was accurately mapped discharging from the adit. Considering that AMD was not classified in September 2018 using ML, the best practice may be to combine AMD with iron oxides when using a broadband sensor. SAM under-predicted jarosite hotspots in November 18, and overpredicted jarosite in September 2018 (Figures 7 and 8). The over-prediction could have been caused by higher reflectance values in the September 18 image, and in this case, a lower threshold may have fixed the misclassification. The shadows in the November 2018 image (Figure 2) and soil moisture caused by typical November rainstorms in the region, likely caused the reflectance to be lower. As clouds and soil moisture will cause a reduction in reflectance, this lead to attenuation in the SAM accuracy for November 18. We recommend that all flights take place prior to any storm events, during cloud-free conditions. SAM under-predicted iron oxides on the west side of the E location in all surveys, and instead consistently mapped patches of jarosite. We confirmed hematite in the soil at the west end of E using both field spectroscopy and XRD (Figures 11b and A1). However, the spectral profile for the red soil at E (Figure 11b) had an absorption feature at 450 nm, and a downturn at 700 nm, similar to the jarosite profiles (Figure 11a). Because we used a 5-band multispectral sensor, SAM mapped these areas as only jarosite. Our results suggest that the soil at E was a mixture of both jarosite and hematite, and that a broadband sensor combined with SAM will only map one class, even when two are present.

Band Math Ratios
The only unsupervised classification we tested was band math ratios (BMRs), which we found performed as well or better than SAM and ML without the added complication of assessing individual pixels or additional classes. We determined uniform thresholds by creating polygons that best fit our visual inspection of the imagery, and we masked vegetation in the images using NDVI. Although BMR did not allow for an accuracy assessment, as it depends on relative band intensities, it is the easiest to implement, and no soil control class was needed. We simplified band ratio determination by utilizing the reference spectrum from both the USGS profiles and the ML polygon profiles. Similar to SAM, we were unable to map AMD as its own class, and we did not attempt to map EMSs. AMD was incorporated into the jarosite and iron oxide endmembers, and was successfully mapped as jarosite in September 2018, March 2019, and July 2019 ( Figure 8). As this was the only method which mapped AMD where it was present, we suggest that, when using band math ratios, incorporating AMD into another endmember may be the best approach to map AMD in small adit seeps. This method under-predicted iron oxides at E, again suggesting that broadband sensors, in conjunction with spectral endmember classification schemes, may only map one class even when two are present. The shadows and possible soil moisture in the N18 image also proved to be an issue when classifying jarosite using BMR, yet the higher reflectance in the S18 image did not appear to effect the classification. This is another reason we believe that the BMR method is the simplest and best approach for mapping PAGM with a broadband sensor, as it is not as sensitive to images with higher reflectance values.

Validation Using Hyperspectral Sensors, Ground Surveys, and Soil Analysis
The hyperspectral surveys evaluated using spectral angle mapper, combined with XRD on selected soils, provided a robust ground-truthing method. Hyperspectral sensors have more bands with narrower bandwidths, enabling them to better detect AMD and PAGM than the broadband 5-channel MicaSense sensor. In the hyperspectral data, the absorption feature of jarosite at 450 nm was clearly identified in the spectral profiles (Figure 4), allowing confident mapping of the distribution of this material (Figure 9). Future use of a RPAS equipped with a hyperspectral sensor could be used to survey the entire field area. The algorithms would most likely not require absorbing AMD into other endmembers, as was necessary with the broadband sensor. In addition, classifying EMS with the appropriate field collected spectra, would most likely be possible. This method of mapping AML is a promising area of future work once this type of sensor is more readily available to fly large areas on a RPAS.
We performed an extra validation step using a handheld spectrometer in areas along the Jones-Kincaid mine waste rock pile. Ground-truthing using the handheld spectrometer clearly identified the spectral signatures of jarosite and iron oxides at locations where they were mapped using the broadband multispectral sensor (Figure 10). Although we were unable to map EMS using spectral angle mapper and band math ratios, and it proved difficult to map them accurately with maximum likelihood, their presence indicates that EMS enriched with heavy metals were continuing to actively discharge into Perry Creek. The continued presence of EMS may indicate that further remediation may need to be implemented, and that an analysis of heavy metals in the soil and in the EMS may need to be investigated to gain further insight into the contamination of Perry Creek. Future studies may be able to map EMS using field collected spectra as an endmember for SAM. However, we did not test the field spectra we collected in our study.

Interpretation of Multi-Temporal Maps of PAGM and AMD
The persistently mapped 'hotspots' of jarosite on the Jones-Kincaid mine waste cover are an indication that jarosite is either forming on the surface, or is being exposed due to cover erosion over time. We mapped the same jarosite 'hotsposts' throughout each RPAS survey and with each classification technique. Jarosite was consistently mapped at Y, T2, E, and Transects 1-3. We confirmed jarosite at Y and Transects 1-3 using soil samples and XRD from the October 18 and October 20 ground surveys (Table A1 and Figure A1), and using field spectroscopy at Y, T2, and E (Figures 10 and 11a). The waste rock cover has been eroding over time during large storm events, suggesting that jarosite was being exposed as a result of erosion at Y, T2, and E, and not necessarily forming on the surface. The XRD results from the trench at Y confirmed this, as jarosite was present from the surface to at least a depth of 0.3 m. Jarosite was not detected at the adjacent trench site S, which implied that jarosite does not exist everywhere on the waste pile, but only in certain areas. This was consistent with our classification maps, as jarosite was not mapped everywhere on the waste rock. Based on the high levels of sulfates in the 2004 site characterization report and our validation data, we believe that jarosite was most likely already present on the mine waste pile prior to cover installation, and is being exposed via erosion in Y, T2, and E. However, the XRD results from the O18 ground surveys show that jarosite was present on the surface of Transects 1-3. Jarosite may have been forming on the surface due to adit seepage at Transects 1-2, and groundwater seepage at Transect 3. The pyrite at T2, although a small amount, was indicative of the potential for active AMD formation, which can lead to the formation of secondary sulfates such as jarosite. The pipe installation during remediation was laid in order to redirect mine waters away from waste dump material. The waste rock did not appear to be eroding at Transect 3 based on previous change detection studies, and the presence of jarosite in this transect could be a result of various scenarios. Jarosite at Transect 3 could mean that groundwater with mobilized sulfates is seeping into the waste rock and oxidizing on the surface, or that the pipe installed during remediation is not efficiently redirecting mine waters away from waste dump material, which can also cause the formation of jarosite on the surface. In either case, jarosite was consistently mapped in all surveys in Transect 3.
Iron oxides were more difficult to accurately map using the three methods. We did, however, map iron oxides at A in each survey, with each method. The iron oxide rocks at A were stable throughout each survey, and had a distinct upturned slope (Figure 11b). This is an indication that iron oxide rocks can be mapped accurately using a 5-band multispectral sensor in conjunction with both supervised and unsupervised classification schemes. The iron rich soils at E are most likely a mixture of jarosite and hematite, and it is likely that these are also consistent 'hotspots'. The iron rich soils were only mapped as distinct classes us the ML pixel based classification, and were mapped as jarosite using SAM and BMRs. We believe that classifying these mixed soils as jarosite sill reveals areas of PAGM on the waste rock. Therefore, we maintain that the BMR method is still the most efficent at mapping PAGM across AML mine waste piles, even in the case of mixed soils.

Recommendations for Use on Other Abandoned Mine Lands
Planning RPAS missions based on weather conditions and necessary time increments, is an important step for mapping temporal changes across multiple surveys efficiently. The cloud cover and possible soil moisture in the N18 survey lead to an under-prediction of jarosite across all methods. Multispectral sensors capture data using the light energy reflected off of objects in the environment. We recommend that flights take place within two and a half hours of local solar noon to reduce deep shadows, which can significantly affect results. The only exception would be flying in very sunny conditions with the sun directly overhead. This can produce high albedo areas which, in our case, likely caused the misclassification of EMS. The 3 to 6 months cadence of RPAS flights may not be necessary, as 'hotspots' of PAGM did not rapidly change throughout the 29-month period of our surveys. However, we did map an increase in PAGM at Y, T2, E, and Transects 1-3 (1,2,3) from June 2018 to July 2019, and again from July 2019 to October 2020, indicating that yearly time-steps would be sufficient to map any changes. As erosion in arid climates generally occurs just after a large rainstorm event, it may be the best practice to plan surveys post heavy rainfall events, when the soil has dried.
Accuracy assessments for classified maps can provide an estimation as to how well a chosen method worked. However, the accuracy assessment for ML and SAM was only performed using 50 random pixels, representing less than 0.1% of the total pixels classified. While the randomized approach should be statistically representative of the scene overall, a more comprehensive assessment would need to include many more pixels for each survey that would have to be manually evaluated. This would have been extremely time consuming and made for a difficult workflow. Here we chose not to evaluate every pixel, as the goal was to make a smooth and simple workflow which can be applied to broadband sensors. Therefore, band math ratios is the most straightforward method to eliminate the labourintensive accuracy assessment, without compromising the validity of the classification.
The Jones-Kinaid mine waste pile in Perry Canyon, NV was a unique AML to fly a broadband sensor and test different supervised and unsupervised classification methods. Based on the site characterization data, we knew we were looking for sulfates and iron rich minerals associated with AMD. The yellow and red soils provided us two minerals (jarosite and iron oxides) associated with PAGM, that we could map and classify using a broadband sensor and various techniques. However, there are many other sulfates associated with PAGM and AMD, and not all of them are yellow or red in color. The advantage of using a broadband multispectral sensor is that many iron bearing sulfate minerals commonly associated with PAGM, exhibit the same spectral features within the 400-850 nm range. In this paper, we showed that jarosite and iron oxides can be mapped accurately with RPAS surveys using a multispectral sensor combined with the band math ratio classification. The specific band ratios we used also have the potential to map and monitor AMLs in which no previous site characterization data is available. In this regard, we found band math ratios to be the easiest and most efficient temporal mapping method for PAGM on AMLs when using a remotely piloted aerial system equipped with a 5-band multispectral sensor.

Conclusions
Temporal monitoring of PAGM and AMD in Perry Canyon using a 5-band multispectral sensor provided an alternative approach to field measurements. In particular, for small area AMLs, a RPAS with a 5-band multispectral sensor provided an efficient way to map the temporal changes of AMD and PAGM, while increasing the spatial coverage and temporal resolution over traditional field methods. When using this type of broadband sensor, we recommend that anything in the image other than the area of interest be masked using NDVI, in order to improve classification accuracy. We mapped PAGM on the Jones-Kincaid mine waste rock pile, including jarosite and iron oxides\AMD. These were readily characterized with the multispectral sensor using maximum likelihood (ML), spectral angle mapper (SAM), and band math ratios (BMRs), with consistent results between the three methods. The EMS material was not accurately mapped due to confusion with other high albedo materials, or the lack of appropriate spectral reference data. However, the persistence of EMS on the waste rock does indicate that further remediation and an analysis of heavy metals may need to be investigated to gain further insight into the contamination of Perry Creek. BMRs combined iron oxides, AMD, and jarosite, but collectively these classes mapped the distribution of PAGM at this site. Although each classification was efficient and accurate in different ways, for future monitoring with frequent surveys, BMR was the simplest and easiest method to map areas of PAGM using a 5-band multispectral sensor. Our multi-temporal survey classifications indicated change in PAGM exposure is occuring slowly, such that 3-6 months time-steps may not be necessary when mapping temporal changes for this particular site. Although these methods are not a replacement for ground measurements, they can be used to guide in silico work to highest priority regions for detailed investigation and remediation. RPASs have demonstrated potential as a means to map high priority AMLs, using even broadband, multi-channel sensors like the MicaSense. Future development, such as incorporating hyperspectral sensors, is expected to improve accurate identification and mapping of many different types of PAGM. Although these methods are not a replacement for ground measurements, they can be used to guide in silico work to highest priority regions for detailed investigation and remediation. RPAS have demonstrated potential as a means to map high priority AMLs using even simple multi-channel cameras like the MicaSense. Future development, such as incorporating hyperspectral sensors, is expected to improve accurate identification and mapping of many different types of PAGM.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript:  Table A1. XRD Results for the October 2018 ground survey, including mineral identification, chemical formulas, and modal percent weight. Trace phase grouping definitions were based on weight percentage ranges, where greater than 30 were considered major phases, less than 30 and greater than 10 were considered moderate, less than 10 and greater than 2 were considered minor, and less than 2 were considered trace phases. When peak to background ratios were low, the detection limit was 2 wt%, and when peak to background ratio was high, the detection limit could be as low as 0.5 wt%. Each sample was ground for five minutes in a McCrone Micronizing Mill using reagent alcohol to reduce particles to between 5 and 10 µm.

Appendix B. Maximum Likelihood Confusion Matrices
We created confusion matrices for the Maximum Likelihood and Spectral Angle Mapper classifications, where columns represent true classes, and rows represent the predictions. All correct classifications are represented along the upper-left to lower-right diagonal.