Identification of Multi-Style Hydrothermal Alteration Using Integrated Compositional and Topographic Remote Sensing Datasets

The western part of the island of Milos, Greece has undergone widespread, intense alteration associated with a range of mineralization, including seafloor Mn-Fe-Ba, sub seafloor Pb-Zn-Ag, and epithermal Au-Ag. The surrounding country rocks are a mixture of submarine and subaerial calc-alkaline volcanic rocks ranging from basaltic andesite to rhyolite in composition, but are predominantly andesites and dacites. The current surface spatial distribution of the alteration mineralogy is a function not only of the original hydrothermal, but also subsequent tectonic and erosional processes. The high relief and the excellent rock exposure provide ideal conditions to evaluate the potential of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite remote sensing data to identify and differentiate the different styles of alteration mineralisation. Laboratory spectral reflectance and calculated emittance measurements of field samples, supported by XRD analysis and field mapping, were used to support the analysis. Band ratio and spectral matching techniques were applied to the shortwave-infrared (SWIR) reflectance and thermal-infrared (TIR) emissivity imagery separately and were then integrated with topographic data. The band ratio and spectral matching approaches produced similar results in both the SWIR and TIR imagery. In the SWIR imagery, the advanced argillic, argillic and hydrous silica alteration zones were clearly identifiable, while in the TIR imagery, the silicic and advanced argillic alteration zones, along with the country rock, were differentiable. The integrated mineralogical–topographic datasets provided an enhanced understanding of the spatial and altitude distribution of the alteration zones when combined with conceptual models of their genesis, which provides a methodology for the differentiation of the multiple styles of alteration.


Introduction
Milos island is a young (<3 Ma) and recently emergent submarine volcano that records a range of contemporaneous mineralisation styles, from seafloor Mn-Fe-Ba [1], to sub seafloor Pb-Zn-Ag mineralisation [2], to epithermal Au-Ag mineralization [3,4].It also shares similarities with active submarine mineralization associated with the volcanic centres of Palinuro [5] and Kolumbo [6] in the Aeolian and Aegean arcs (Figure 1A).Moreover, this shallow submarine volcanic-hydrothermal environment is where large gold and silver volcanic massive sulfide deposits such as Eskay Creek are thought to form [7]. environment is where large gold and silver volcanic massive sulfide deposits such as Eskay Creek are thought to form [7]. A sequence of submarine and subaerial volcanism, followed by hydrothermal alteration, significant tectonic activity, and subsequent erosion, has resulted in an area of high relief (Figure 2).The combination of the high spatial compositional heterogeneity of the mineralisation, the scale of the study area (100 km 2 ), and the high relief pose a significant challenge to field-based studies.The excellent rock exposure and the very dry climate, with resultant lack of vegetation (Figure 3), makes Milos an ideal location to study the complex interactions between these different mineralization styles and the efficacy of remote sensing approaches to detect a range of mineralization styles, with applications to mineral exploration for gold-rich volcanic massive sulfide systems.
This study evaluates the utility of band ratio and spectral matching techniques to Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) SWIR reflectance and thermal emissivity integrated with topographic data to identify and differentiate the multiple styles of alteration in the western part of the island of Milos.A sequence of submarine and subaerial volcanism, followed by hydrothermal alteration, significant tectonic activity, and subsequent erosion, has resulted in an area of high relief (Figure 2).The combination of the high spatial compositional heterogeneity of the mineralisation, the scale of the study area (100 km 2 ), and the high relief pose a significant challenge to field-based studies.The excellent rock exposure and the very dry climate, with resultant lack of vegetation (Figure 3), makes Milos an ideal location to study the complex interactions between these different mineralization styles and the efficacy of remote sensing approaches to detect a range of mineralization styles, with applications to mineral exploration for gold-rich volcanic massive sulfide systems.
This study evaluates the utility of band ratio and spectral matching techniques to Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) SWIR reflectance and thermal emissivity integrated with topographic data to identify and differentiate the multiple styles of alteration in the western part of the island of Milos.

Geological Setting of the Study Area
The island of Milos, located in the Aegean arc (Figure 1A), is associated with a zone of Pliocene to modern volcanism related to active backarc extension caused by the subduction of the African plate beneath the Aegean microplate [3].It is built on continental crust, comprises seven major volcanic centers, and is located 120-250 km north of the Hellenic Trench (Figure 1A).The volcanic rocks are calc-alkaline with localized high-K variants and range from basalt to rhyolite in composition, but are predominantly andesites and dacites [22].Present-day hydrothermal activity

Geological Setting of the Study Area
The island of Milos, located in the Aegean arc (Figure 1A), is associated with a zone of Pliocene to modern volcanism related to active backarc extension caused by the subduction of the African plate beneath the Aegean microplate [3].It is built on continental crust, comprises seven major volcanic centers, and is located 120-250 km north of the Hellenic Trench (Figure 1A).The volcanic rocks are calc-alkaline with localized high-K variants and range from basalt to rhyolite in composition, but are predominantly andesites and dacites [22].Present-day hydrothermal activity

Geological Setting of the Study Area
The island of Milos, located in the Aegean arc (Figure 1A), is associated with a zone of Pliocene to modern volcanism related to active backarc extension caused by the subduction of the African plate beneath the Aegean microplate [3].It is built on continental crust, comprises seven major volcanic centers, and is located 120-250 km north of the Hellenic Trench (Figure 1A).The volcanic rocks are calc-alkaline with localized high-K variants and range from basalt to rhyolite in composition, but are predominantly andesites and dacites [22].Present-day hydrothermal activity comprises both low-enthalpy (Aegina, Sousaki, Methana, Greece) and high-enthalpy (Milos, Nisyros, Santorini, Greece) systems.Late Pliocene submarine and late Pleistocene to Holocene subaerial volcanic rocks overlie Mesozoic metamorphic basement and upper Miocene-lower Pliocene marine sedimentary rocks, recording a transition from a shallow submarine (200 m) to subaerial volcanic setting (Figure 1B).Emergence probably occurred in the middle to late Pleistocene [22].Plutonic rocks are not known on Milos and have only been reported as ignimbrite hosted granitic xenoliths from the nearby islet of Kimolos [22].
The oldest submarine volcanic rocks occur on western Milos and host low sulfidation gold-silver mineralization-the Profitis Ilias (5.5 Mt at 4.4 g/t gold and 43 g/t silver) and Chondro Vouno (1.2 Mt at 1 g/t gold and 124 g/t silver) deposits (Figure 1B), which constitute a 20 km 2 epithermal system with a combined resource of 0.8 Moz gold and 12.4 Moz silver.Mineralization is hosted by silicified, sericitized, and adularized rhyolitic tuffs and ignimbrites and consists of crustiform or colloform-banded quartz (or chalcedony), adularia, barite, and sericite veins related to quartz, adularia, sericite (or illite) wall-rock alteration.Mineralized veins extend to depths of at least 300 m below the current surface.Hypogene metallic minerals include pyrite, galena, chalcopyrite, sphalerite, marcasite, tetrahedrite, native gold, electrum and rare tellurides.Sulfides mainly occur deep in the system as disseminated stockwork mineralization of unknown vertical extent.Elevated gold values are concentrated above the base-metal zone and are spatially related to boiling [3,4].
Milos island hosts several hydrothermal mineral deposits, and the host volcanic rocks have undergone intense hydrothermal alteration.Alfieris, et al. [2] state that alteration zones are similar to those described in high-and intermediate/low-sulfidation epithermal mineralizing systems.On the basis of alunite grain size and habit, they also distinguish hypogene and steam heated advanced argillic alteration.However, this distinction cannot be determined from multi-spectral satellite datasets such as ASTER.Therefore, for the purposes of this study, we are only examining the distribution of advanced argillic alteration and do not make any reference to its origin.In the north of western Milos, they identify NE-trending silicic alteration outcrops.These grade laterally and vertically into advanced argillic, argillic, and propylitic alteration and fresh rock.In addition, there is localized phyllic (sericitic) alteration.Silicic alteration comprises fine-grained quartz, TiO 2 polymorphs, alunite, and pyrite, and occurs as vuggy and massive-silica types.Advanced argillic alteration is typified by intergrowths of alunite, quartz, chalcedony, opal-CT, kaolinite, halite, dickite, pyrophyllite, diaspore, pyrite, chalcedony, and jarosite.Towards the south of western Milos are blanket-like zones of advanced argillic alteration.Locally, in the northern part, this style of alteration overprints the hypogene advanced argillic alteration [2].Argillic alteration comprises illite ˘kaolinite ˘smectite ˘sericite.Its main occurrence is marginal to the main centers of advanced argillic alteration.However, it can also occur in fault zones [2].In terms of abundance of individual alteration minerals, kaolinite is generally more widespread in the advanced argillic alteration associated with sub-seafloor Pb-Zn-Ag mineralization in the Triades-Galana-Kondaros area (Location 2, Figure 1B), in contrast, illite is more abundant in the middle and southern part of the study area, and is associated with epithermal Au-Ag mineralization in the Chondro Vouno-Profitis Ilias area (Location 3, Figure 1B).Phyllic alteration is restricted to the area around Galana, and is typified by the presence of sericite, quartz, and pyrite.Adularia-sericite alteration-typical of epithermal Au mineralization-envelopes quartz veins at Profitis Ilias, Chondro Vouno, Galena, and Triades [2].

Field Mapping and XRD Analyses
Two sets of field and XRD datasets were utilized in this study.An extensive field campaign was carried out by the authors in May 2010, and a representative set of field samples was acquired (Figure 2).XRD analyses were carried out on these of samples, with a subset of the results presented in Table 1.In addition, Alfieris et al. [2] present a map of the alteration zones of western Milos produced using field mapping and XRD analysis [2], Figure 4.

VNIR-SWIR Reflectance Spectroscopy
VNIR-SWIR (λ = 0.4´2.5 µm) reflectance spectra were collected using an Analytical Spectral Devices (ASD) FieldSpec Pro FR spectroradiometer.The ASD has a spectral range of 0.35-2.5 µm.Spectra were collected at sampling intervals of 1.4 nm at 0.35-1.0µm and 2 nm at 1.0-2.5 µm, with the following full width half maximum (FWHM) spectral resolutions: 6 nm on average at VNIR wavelengths (though up to 3 nm at ~0.7 µm) and 11 nm at SWIR wavelengths.Absolute reflectance values were calculated by calibrating each batch of ~20 repeat measurements against a Spectralon white reference panel that was also measured at the end of each batch to monitor for instrumental drift.Individual measurements sample an ellipse approximately 10 ˆ5 mm in size.In order to maximize the quality of VIS-SWIR spectra, ~50 ˆ50 mm flat surfaces were cut into the samples and polished coarsely in order to ensure good optical coupling with the ASD contact probe.The effect of compositional heterogeneity on sample spectra was mitigated by averaging over up to 94 evenly spaced repeat measurements taken from across sample surfaces.The typical precision of individual spectra estimated from repeat measurements of an identical field of view is <0.5.

ASTER Satellite Remote Sensing Dataset
ASTER datasets covering the SWIR and thermal-infrared (TIR) waveranges (Table 2) were used in the study [22].An ASTER SWIR surface reflectance data product (AST-07XT) was downloaded from the Earth Remote Sensing Data Analysis Center in Japan.The AST-07XT imagery was used, as it already has the CROSSTALK correction applied and has been found to be more reliable than previous ASTER products for discriminating hydrothermal alteration minerals and mineral groups without the use of additional spectral data from the site for calibration [18,23].The ASTER spectral emittance data used in this study was the standard product (AST_05), which was produced at the EROS Data Center using a MODTRAN-based atmospheric correction and the Temperature-Emissivity Separation TES) algorithm developed by [24].The spectral emittance data were co-registered with the SWIR data by applying a 6ˆpixel duplication of the TIR data.

ASTER Data Processing
Band Ratios (Table 3), [18][19][20]25] and spectral matching [26] analyses were applied to the SWIR and TIR datasets.The ASTER satellite has 14 spectral bands offering many permutations of Band Ratios [18][19][20]25], and spectral matching [26] analyses were applied to the SWIR and TIR datasets.The spectral matching algorithm implemented in this study was the Constrained Energy Minimisation (CEM) technique [26].It performs a matched filtering (MF) of multispectral/ hyperspectral images, which minimises the response of background materials by projecting each pixel vector onto a subspace-which is orthogonal to the background spectra-and then maximises the response of the endmembers of interest by comparing the residual pixels to each of the reference spectra [26].It was applied to both the SWIR and TIR ASTER images separately.The spectral matching approach can utilise a range of input spectral profile sources.Acquisition of a representative set of field spectra of the appropriate alteration zones should provide the optimum mapping approach; however, in practice, this is often very difficult to achieve due to sampling-related issues.Alternative approaches are using library spectra (e.g., the ASTER spectral library [27]) or spectral endmembers retrieved from the imagery.There are two pre-processing stages required to retrieve the spectral endmembers from the image dataset.Firstly, a Minimum Noise Function (MNF) transformation must be applied to determine the inherent dimensionality of the image data and to reduce the computational requirements for subsequent processing [28,29].Secondly, a Pixel Purity Index (PPI) is then applied using the MNF result as the input data, which extracts the image endmembers [28,29].

Spectral Reflectance
Laboratory spectra covering the visible, near-infrared, and shortwave infrared waveranges of minerals characteristic of hydrothermally altered rocks are shown in Figure 5.The spectral matching algorithm implemented in this study was the Constrained Energy Minimisation (CEM) technique [26].It performs a matched filtering (MF) of multispectral/hyperspectral images, which minimises the response of background materials by projecting each pixel vector onto a subspace-which is orthogonal to the background spectra-and then maximises the response of the endmembers of interest by comparing the residual pixels to each of the reference spectra [26].It was applied to both the SWIR and TIR ASTER images separately.The spectral matching approach can utilise a range of input spectral profile sources.Acquisition of a representative set of field spectra of the appropriate alteration zones should provide the optimum mapping approach; however, in practice, this is often very difficult to achieve due to samplingrelated issues.Alternative approaches are using library spectra (e.g., the ASTER spectral library [27]) or spectral endmembers retrieved from the imagery.There are two pre-processing stages required to retrieve the spectral endmembers from the image dataset.Firstly, a Minimum Noise Function (MNF) transformation must be applied to determine the inherent dimensionality of the image data and to reduce the computational requirements for subsequent processing [28,29].Secondly, a Pixel Purity Index (PPI) is then applied using the MNF result as the input data, which extracts the image endmembers [28,29].

Spectral Reflectance
Laboratory spectra covering the visible, near-infrared, and shortwave infrared waveranges of minerals characteristic of hydrothermally altered rocks are shown in Figure 5. Laboratory spectra of samples collected on Milos representative of the silicic, advanced argillic, argillic, propylitic, and adularia-sericitic alteration zones are shown in Figure 6A.The reflectance spectra of samples MI-05 and MI-27 from the study area indicates that they contain minerals that are typical of advanced argillic alteration, attributed to AL-OH absorption with the kaolinite features at 2.17 and 2.2 μm being very prominent.The reflectance spectra of samples MI-22 and MI-28 from the study area exhibit an intense absorption feature at 2.17 μm, along with secondary features located near 2.32 and 1.760 μm and the doublet O-H feature in the 1.5 μm region, which are diagnostic of alunite.Although kaolinite is present, the intense alunite features obscure the kaolinite features at 2.17 and 2.2 μm.Examination of the reflectance spectra of sample MI-22 and MI-28 from the study Laboratory spectra of samples collected on Milos representative of the silicic, advanced argillic, argillic, propylitic, and adularia-sericitic alteration zones are shown in Figure 6A.The reflectance spectra of samples MI-05 and MI-27 from the study area indicates that they contain minerals that are typical of advanced argillic alteration, attributed to AL-OH absorption with the kaolinite features at 2.17 and 2.2 µm being very prominent.The reflectance spectra of samples MI-22 and MI-28 from the study area exhibit an intense absorption feature at 2.17 µm, along with secondary features located near 2.32 and 1.760 µm and the doublet O-H feature in the 1.5 µm region, which are diagnostic of alunite.Although kaolinite is present, the intense alunite features obscure the kaolinite features at 2.17
The XRD analyses show that silicate-based minerals make up a significant proportion of the rock samples collected on Milos.Using the results from the XRD analyses, the emission spectra from the Jet Propulsion Laboratory, and the Arizona State University spectral libraries, emission spectra for the sampled rocks were generated (Figure 7B).The emission spectra from the sample representative of the silicic alteration zone has a very strong double absorption feature between 8.4 and 8.95 μm produced by the almost pure quartz mineralogy.The emissivity spectra from the advanced argillic zone has two distinctive absorption features at 9.2 and 9.9 μm.The argillic and the propylitic emission spectra are similar, with a broad absorption feature from 8.95 to 10 μm.There is a distinguishable shift of the absorption feature to longer wavelengths moving from the silic to the advanced argillicargillic-propylitic alteration zones.
The XRD analyses show that silicate-based minerals make up a significant proportion of the rock samples collected on Milos.Using the results from the XRD analyses, the emission spectra from the Jet Propulsion Laboratory, and the Arizona State University spectral libraries, emission spectra for the sampled rocks were generated (Figure 7B).The emission spectra from the sample representative of the silicic alteration zone has a very strong double absorption feature between 8.4 and 8.95 µm produced by the almost pure quartz mineralogy.The emissivity spectra from the advanced argillic zone has two distinctive absorption features at 9.2 and 9.9 µm.The argillic and the propylitic emission spectra are similar, with a broad absorption feature from 8.95 to 10 µm.There is a distinguishable shift of the absorption feature to longer wavelengths moving from the silic to the advanced argillic-argillic-propylitic alteration zones.1).

Results of ASTER Image Analysis
The accuracy of the mineralogical maps output from the analysis of the ASTER SWIR and TIR datasets were assessed using the alteration zones map [2] and XRD analyses of rock samples.

SWIR Data
The resampled ASTER laboratory spectra of the field samples retain definition of the most intense absorption features, particularly the intense absorption feature centered near 2.205 μm (Figure 6B).The different alteration zones are differentiable despite the significant decrease in spectral resolution.The argillic spectrum has a single differentiable feature centered on 2.205 μm, with the shoulder at 2.165 μm significantly lower than the shoulder at 2.260 μm.The adularia-sericite spectrum is also centered on 2.205 μm, but the higher wavelength shoulder is located at 2.330 μm.The advanced argillic spectrum has a pronounced absorption feature which extends from 2.165 to 2.330 μm with a distinct gradient between 2.205 and 2.260 μm.The silicified sinter (silicic) spectrum has a distinct absorption feature extending to 2.330 μm.Both the adularia-sericitic and the slightly altered dacite samples only have one very weak absorption feature centered at 2.216 μm The spectral endmembers representing rocks containing spectrally-dominant alunite-kaolinite, illite-smectite-sericite, and hydrous silica identified using the PPI processing are shown in (Figure 8).Only three spectral distinct mineralogical endmembers could be resolved from the imagery.The spectra representative of the advanced argillic alteration zone shows the distinctive, intensive, broad absorption feature between bands five and six are characteristic of significant alunite content with secondary kaolinite.The spectra representative of the argillic alteration zone shows the distinctive (but relatively weak) absorption feature in band six, while the spectra representative of hydrous silica shows a distinctive spectra edge between bands five and nine.The spectral responses from sericite, smectite, and illite in the propylitic and adularia-sericitic alteration zones were not sufficient to differentiate from the argillic zone.The overall shapes of the ASTER-derived spectra are different from the field spectra resampled to the ASTER bandwidths because of the effects of several factors that influence the spectra, but are obviated in the laboratory.They include the presence of variable mineral mixtures, grain-size variations, residual atmospheric absorption features, and ASTER calibration errors [30].1).

Results of ASTER Image Analysis
The accuracy of the mineralogical maps output from the analysis of the ASTER SWIR and TIR datasets were assessed using the alteration zones map [2] and XRD analyses of rock samples.

SWIR Data
The resampled ASTER laboratory spectra of the field samples retain definition of the most intense absorption features, particularly the intense absorption feature centered near 2.205 µm (Figure 6B).The different alteration zones are differentiable despite the significant decrease in spectral resolution.The argillic spectrum has a single differentiable feature centered on 2.205 µm, with the shoulder at 2.165 µm significantly lower than the shoulder at 2.260 µm.The adularia-sericite spectrum is also centered on 2.205 µm, but the higher wavelength shoulder is located at 2.330 µm.The advanced argillic spectrum has a pronounced absorption feature which extends from 2.165 to 2.330 µm with a distinct gradient between 2.205 and 2.260 µm.The silicified sinter (silicic) spectrum has a distinct absorption feature extending to 2.330 µm.Both the adularia-sericitic and the slightly altered dacite samples only have one very weak absorption feature centered at 2.216 µm The spectral endmembers representing rocks containing spectrally-dominant alunite-kaolinite, illite-smectite-sericite, and hydrous silica identified using the PPI processing are shown in (Figure 8).Only three spectral distinct mineralogical endmembers could be resolved from the imagery.The spectra representative of the advanced argillic alteration zone shows the distinctive, intensive, broad absorption feature between bands five and six are characteristic of significant alunite content with secondary kaolinite.The spectra representative of the argillic alteration zone shows the distinctive (but relatively weak) absorption feature in band six, while the spectra representative of hydrous silica shows a distinctive spectra edge between bands five and nine.The spectral responses from sericite, smectite, and illite in the propylitic and adularia-sericitic alteration zones were not sufficient to differentiate from the argillic zone.The overall shapes of the ASTER-derived spectra are different from the field spectra resampled to the ASTER bandwidths because of the effects of several factors that influence the spectra, but are obviated in the laboratory.They include the presence of variable mineral mixtures, grain-size variations, residual atmospheric absorption features, and ASTER calibration errors [30].

Results
The results of the band ratio and matched analysis of the SWIR image are shown in Figure 9.In both the band ratio and MF images, the degree of match is indicated by higher numbers representing good matches, and vice versa [28,29].Pixels corresponding to the highest degree of match for each class were transferred digitally to the aerial photography mosaic and assigned a colour.Determination of the threshold range to represent each class was based on the coherence of the spatial distribution of pixels and lack of widely scattered pixels.Pixels that are not colour-coded are unclassified.The overall distribution of altered rocks is very similar in both the band ratio and MF images.
The advanced argillic altered rocks are distributed unevenly across the study area with a discontinuous belt trending northwest to southeast across the central section of western Milos.There are also a number of significant advanced argillic discrete occurrences in the northern part of the study area.The distribution of the argillic alteration rocks are mainly located in a belt trending northwest to southeast across the centre of the study area.
There are additional small occurrences found peripherally to some of the advanced argillic occurrences.The hydrous silica index resolves the pozzolan quarries on the southern coast of the study area very clearly.In addition, there are a few very small, discrete occurrences in the centre and north of the study area.The samples M-04, MI-05, MI-25, and MI-27 are located in the mapped argillic altered zone, supporting the accuracy of the ASTER-derived maps of advanced argillic alteration.Sample MI-26 is indicative of slightly altered/country rock which again supports the accuracy of the ASTER-derived mineralogy maps.

Results
The results of the band ratio and matched analysis of the SWIR image are shown in Figure 9.In both the band ratio and MF images, the degree of match is indicated by higher numbers representing good matches, and vice versa [28,29].Pixels corresponding to the highest degree of match for each class were transferred digitally to the aerial photography mosaic and assigned a colour.Determination of the threshold range to represent each class was based on the coherence of the spatial distribution of pixels and lack of widely scattered pixels.Pixels that are not colour-coded are unclassified.The overall distribution of altered rocks is very similar in both the band ratio and MF images.
The advanced argillic altered rocks are distributed unevenly across the study area with a discontinuous belt trending northwest to southeast across the central section of western Milos.There are also a number of significant advanced argillic discrete occurrences in the northern part of the study area.The distribution of the argillic alteration rocks are mainly located in a belt trending northwest to southeast across the centre of the study area.
There are additional small occurrences found peripherally to some of the advanced argillic occurrences.The hydrous silica index resolves the pozzolan quarries on the southern coast of the study area very clearly.In addition, there are a few very small, discrete occurrences in the centre and north of the study area.The samples M-04, MI-05, MI-25, and MI-27 are located in the mapped argillic altered zone, supporting the accuracy of the ASTER-derived maps of advanced argillic alteration.Sample MI-26 is indicative of slightly altered/country rock which again supports the accuracy of the ASTER-derived mineralogy maps.

TIR Data
Resampling the high-resolution laboratory spectra to the five ASTER TIR bandpasses results in a general loss of definition of specific absorption features, except for the main quartz reststrahlen feature located in ASTER band 12 (Figure 10A).Only three spectral endmembers were identified from the PPI analysis of the TIR wavelength imagery representing the silicic (quartz), advanced argillic (alunite-kaolinite) and country rock (Figure 10B).The spectra representative of the silicic alteration zone shows the distinctive, intense, broad emissivity minimum between bands 10 and 12 characteristic of a very high quartz content.The spectra representative of the advanced argillic alteration zone shows a distinctive emissivity minimum in band 12 but with higher values in bands 10 and 11.The spectra representative of the country rock also shows an emissivity minimum between bands 11 and 12 with a weak emissivity feature in band 14.Spectral endmembers representing the argillic and propylitic alteration zones were not identifiable.The coarse spatial resolution of the ASTER TIR data can have a significant impact on the capability to map surface mineralogy.Despite the effects due to the size of the field-of-view and surface and atmospheric effects, the generated emissivity spectra of MI-21 and MI-05 are similar to the ASTER-derived emissivity spectra.
The results of the Quartz Index (QI) and MF of the TIR imagery are shown in Figure 11, overlain by the alteration map [2].The overall distribution of silicic rocks is very similar in both the QI and MF images.The silicic altered rocks are distributed unevenly across the study area, with a discontinuous belt trending across the central section of western Milos and a significant proportion aligned across the Profitis Ilias ridge.There are also a number of discrete occurrences in the central and northern parts of the study area.The distribution of the argillic alteration rocks in the MF image are associated peripherally to the silicic altered rocks.The distribution of the country rocks correlates closely with the distribution mapped by [2].

TIR Data
Resampling the high-resolution laboratory spectra to the five ASTER TIR bandpasses results in a general loss of definition of specific absorption features, except for the main quartz reststrahlen feature located in ASTER band 12 (Figure 10A).Only three spectral endmembers were identified from the PPI analysis of the TIR wavelength imagery representing the silicic (quartz), advanced argillic (alunite-kaolinite) and country rock (Figure 10B).The spectra representative of the silicic alteration zone shows the distinctive, intense, broad emissivity minimum between bands 10 and 12 characteristic of a very high quartz content.The spectra representative of the advanced argillic alteration zone shows a distinctive emissivity minimum in band 12 but with higher values in bands 10 and 11.The spectra representative of the country rock also shows an emissivity minimum between bands 11 and 12 with a weak emissivity feature in band 14.Spectral endmembers representing the argillic and propylitic alteration zones were not identifiable.The coarse spatial resolution of the ASTER TIR data can have a significant impact on the capability to map surface mineralogy.Despite the effects due to the size of the field-of-view and surface and atmospheric effects, the generated emissivity spectra of MI-21 and MI-05 are similar to the ASTER-derived emissivity spectra.
The results of the Quartz Index (QI) and MF of the TIR imagery are shown in Figure 11, overlain by the alteration map [2].The overall distribution of silicic rocks is very similar in both the QI and MF images.The silicic altered rocks are distributed unevenly across the study area, with a discontinuous belt trending across the central section of western Milos and a significant proportion aligned across the Profitis Ilias ridge.There are also a number of discrete occurrences in the central and northern parts of the study area.The distribution of the argillic alteration rocks in the MF image are associated peripherally to the silicic altered rocks.The distribution of the country rocks correlates closely with the distribution mapped by [2].

Integration of ASTER Mineralogical Maps with Topography
Integration of the mineralogical composition maps produced from the analysis of the SWIR and TIR datasets with a Digital Elevation Model (DEM) provides additional information into the spatial and altitudinal distribution of the altered rocks, which enhances understanding of the threedimensional extent of the alteration zones and their association with the different styles (depths) of alteration, providing a methodology to differentiate the different styles of the multiple alteration events that have occurred on western Milos.
The discontinuous belt of alteration minerals across the centre of the study area is aligned along a ridge with the majority of the silicic alteration zone being location on Profitis Ilias.The spatially discrete alteration occurrences located in the centre and north of the study area are located at significantly lower altitudes.By utilising the conceptual models of the different styles of mineralisation, we can propose differentiating these alteration occurrences on the basis of their current altitude, with the occurrences located at the higher altitude (>250 m) being associated with a low sulfidation epithermal Au-Ag mineralization (Figure 12A) and those at a lower altitude (<250m) being associated with with advanced argillic alteration (Figure 12B) [3].

Integration of ASTER Mineralogical Maps with Topography
Integration of the mineralogical composition maps produced from the analysis of the SWIR and TIR datasets with a Digital Elevation Model (DEM) provides additional information into the spatial and altitudinal distribution of the altered rocks, which enhances understanding of the three-dimensional extent of the alteration zones and their association with the different styles (depths) of alteration, providing a methodology to differentiate the different styles of the multiple alteration events that have occurred on western Milos.
The discontinuous belt of alteration minerals across the centre of the study area is aligned along a ridge with the majority of the silicic alteration zone being location on Profitis Ilias.The spatially discrete alteration occurrences located in the centre and north of the study area are located at significantly lower altitudes.By utilising the conceptual models of the different styles of mineralisation, we can propose differentiating these alteration occurrences on the basis of their current altitude, with the occurrences located at the higher altitude (>250 m) being associated with a low sulfidation epithermal Au-Ag mineralization (Figure 12A) and those at a lower altitude (<250m) being associated with with advanced argillic alteration (Figure 12B) [3].

Discussion
The analysis of the spectral reflectance data in the SWIR imagery clearly outlined the alteration zones from the country rocks.Despite the significant reduction in spectral resolution in the ASTER SWIR data, the hydrous silica, argillic, and advanced argillic altered rocks were clearly identifiable and differentiable.The laboratory spectra of the propylitic, phyllic, and adularia-sericite alteration zones showed that the absorption features were weak, which could explain why these alteration zones were not clearly identifiable in the ASTER data.The spatial distribution of the identified altered rocks were similar in the band ratio and MF analysis, with a discontinuous belt of altered rock running across the centre of the study area clearly visible and with small, discrete occurrence in the central and northern parts of western Milos.Advanced argillic outcrops were very clearly identifiable because of the high alunite and natro-alunite concentrations and the excellent rock exposure.The argillic alteration zones are usually smaller in extent and are often found peripheral to the advanced argillic outcrops.
Although ASTER TIR data has the potential to overcome a significant limitation of the SWIR data-the identification of silicate minerals-mixtures of quartz with other minerals can also affect spectral shape.The presence of these minerals in sufficient concentrations will decrease the prominence of the diagnostic quartz emissivity peak at band 11 or create an emissivity absorption feature at those wavelengths.The presence of montmorillonite, kaolinite, or muscovite (among other minerals) will result in a substantial decrease in emissivity between bands 10 and 12. ASTER SWIR [26,27] can be used to identify many rock-forming minerals that will affect the TIR spectral response.The absence of an ASTER band near 9.8 μm-resulting from intense atmospheric absorption-is an important limitation in discriminating alunite and kaolinite from quartz.

Discussion
The analysis of the spectral reflectance data in the SWIR imagery clearly outlined the alteration zones from the country rocks.Despite the significant reduction in spectral resolution in the ASTER SWIR data, the hydrous silica, argillic, and advanced argillic altered rocks were clearly identifiable and differentiable.The laboratory spectra of the propylitic, phyllic, and adularia-sericite alteration zones showed that the absorption features were weak, which could explain why these alteration zones were not clearly identifiable in the ASTER data.The spatial distribution of the identified altered rocks were similar in the band ratio and MF analysis, with a discontinuous belt of altered rock running across the centre of the study area clearly visible and with small, discrete occurrence in the central and northern parts of western Milos.Advanced argillic outcrops were very clearly identifiable because of the high alunite and natro-alunite concentrations and the excellent rock exposure.The argillic alteration zones are usually smaller in extent and are often found peripheral to the advanced argillic outcrops.
Although ASTER TIR data has the potential to overcome a significant limitation of the SWIR data-the identification of silicate minerals-mixtures of quartz with other minerals can also affect spectral shape.The presence of these minerals in sufficient concentrations will decrease the prominence of the diagnostic quartz emissivity peak at band 11 or create an emissivity absorption feature at those wavelengths.The presence of montmorillonite, kaolinite, or muscovite (among other minerals) will result in a substantial decrease in emissivity between bands 10 and 12. ASTER SWIR [26,27] can be used to identify many rock-forming minerals that will affect the TIR spectral response.The absence of an ASTER band near 9.8 µm-resulting from intense atmospheric absorption-is an important limitation in discriminating alunite and kaolinite from quartz.
The results of the band ratios and MF analysis gave similar results, which matched closely to the alteration zones mapped by [2] and the XRD sample compositions.Spectral matching techniques should in theory produce more accurate results, particularly if field spectra from the study area are available for use.However, in practice, issues relating to sampling in terms of how representative the ground targets are relative to the ASTER image can cause significant problems.While band ratios are a relatively imprecise method of quantifying absorption feature shape and can show significant confusion between surface materials with similar spectral properties, they can often produce better results than the matched-filter technique.
Integration of topographic data with surface composition data can provide additional insight into the nature of alteration processes.The western part of Milos hosts a wide range of hydrothermal alteration styles which were generated at a variety of depths.The present-day topographic distribution of these alteration events has been further affected by tectonic and erosional processes, such that, epithermal Ag-Au mineralisation now outcrops in a discrete higher altitudinal zone, whereas as the advanced argillic alteration occurs at significantly lower altitudes.

Conclusions
This study has demonstrated that ASTER SWIR and TIR imagery have the capability to resolve multiple styles of alteration mineralisation.The TIR datasets are essential in order to resolve non-hydrous silicates associated with hydrothermal alteration and thus fully characterise all the alteration zones associated with hydrothermal alteration processes.The integration of topographic data, combined with a clear understanding of conceptual models of mineralisation, provides a powerful additional tool to help identify and differentiate multiple styles of alteration.

Figure 2 .
Figure 2. Topography of western Milos (in meters) with sample locations.

Figure 3 .
Figure 3. Topography and vegetation of western Milos-looking northeast to Profitis Ilias.

Figure 3 .
Figure 3. Topography and vegetation of western Milos-looking northeast to Profitis Ilias.

Figure 3 .
Figure 3. Topography and vegetation of western Milos-looking northeast to Profitis Ilias.

Figure 6 .
Figure 6.(A): Laboratory spectra of field samples representing the alteration zones.Numbers in labels indicate field site (Table 1); (B): spectra representative of the alteration zones (shown in Figure 3A) resampled to the nine ASTER VNIR + SWIR bandwidths.ASTER band wavelength position shown at top.

Figure 6 .
Figure 6.(A): Laboratory spectra of field samples representing the alteration zones.Numbers in labels indicate field site (Table 1); (B): spectra representative of the alteration zones (shown in Figure 3A) resampled to the nine ASTER VNIR + SWIR bandwidths.ASTER band wavelength position shown at top.

Figure 8 .
Figure 8. ASTER image spectra used as reference in match-filter processing of SWIR data to map advanced argillic, argillic, and hydrous silica classes which are shown in Figure 9. ASTER band wavelength position shown at top.ASTER band wavelength position shown at top.

Figure 8 .
Figure 8. ASTER image spectra used as reference in match-filter processing of SWIR data to map advanced argillic, argillic, and hydrous silica classes which are shown in Figure 9. ASTER band wavelength position shown at top.ASTER band wavelength position shown at top.

Figure 10 .
Figure 10.(A): Calculated emissivity spectra of samples shown in Figure 7B resampled to the five ASTER TIR bandwidths (Table 1); (B): ASTER image spectra used as reference in matched-filter processing of TIR data.ASTER band wavelength position shown at top.

Figure 10 .Figure 11 .
Figure 10.(A): Calculated emissivity spectra of samples shown in Figure 7B resampled to the five ASTER TIR bandwidths (Table 1); (B): ASTER image spectra used as reference in matched-filter processing of TIR data.ASTER band wavelength position shown at top.

Figure 11 .
Figure 11.Results of (A) Quartz Index and (B) Matched Filter analysis of TIR data with map of alteration zones [2] overlain showing the distribution of three classes of hydrothermally altered rocks superposed on aerial photography grey-tone mosaic of the study area: green-advanced argillic; red-silicic; and blue-country rock.

Figure 12 .
Figure 12. Results of Matched Filter analysis of SWIR and TIR data superposed on aerial photography grey-tone mosaic of the study area: green-advanced argillic; red-silicic and blue-argillic; (A) alteration mineralisation distribution above 250 m altitude; (B) alteration mineralisation distribution below 250 m altitude.

Figure 12 .
Figure 12. Results of Matched Filter analysis of SWIR and TIR data superposed on aerial photography grey-tone mosaic of the study area: green-advanced argillic; red-silicic and blue-argillic; (A) alteration mineralisation distribution above 250 m altitude; (B) alteration mineralisation distribution below 250 m altitude.

Table 1 .
Description of the field sites, samples collected, and the results of the XRD analyses.

Table 1 .
Description of the field sites, samples collected, and the results of the XRD analyses.