We discussed the possibilities and open questions concerning laboratory reflectance spectroscopy, spectroscopic measurements in the field, hyperspectral image data from spacecraft, and integration of multiscale data. Open questions included: (1) bulk-rock spectral complexity, which provides the geologic basis for every spectroscopic analysis; (2) criteria for laboratory and field spectra classification, as a tool for (3) end-member selection for image data classification; (4) peculiar spectral characteristics of Mount Etna basalts; and (5) effects of remote-sensing data quality. The last three items emerged during a multiscale survey on the Mount Etna volcano.
Our laboratory spectroscopic analyses, supported by specific petrographic analyses, showed the relationship between absorption-band frequency and spectrally active functional groups and the unexpected effects of bulk-rock composition on this relationship. We studied the muscovite Al-OH band in quartzite and micaschist and Fe2+ band in pyroxene-bearing cumulates.
Laboratory reflectance spectra of rocks were classified using the concept of spectro-facies. In the case of metamorphic rocks, the result was a tree-structure of rock spectral classes mainly based on the predominant vibrational processes. In the case of basalts, characterized by an overall similarity in their composition, the classes were determined on the basis of overall shape of the spectral curve and on electronic process intensities. Here, we report first results of multiscale data integration for the Mount Etna volcano. Etna rocks consist of basalts, with very low albedo and variable degrees of alteration, and recent lava flows are characterized by overall low reflectance in both ASTER (advanced spaceborne thermal emission and reflection radiometer) and Hyperion color-composite images. We carried out Spectral Angle Mapper (SAM) classification of Hyperion images, where individual field spectra represented suitable end members for classification of recent lava and pyroclastic deposits. We used field spectra linear combinations to classify mixed pixels and to approximate the classification of altered and oxidized effusive products. Only two laboratory spectral classes coincided with field spectra classes; laboratory spectra were mainly used for spectral features attribution. The overall spectral shape of some of these spectra is still under study. Noise level in Hyperion data precluded the identification of subtle diagnostic iron absorption bands.
Spectroscopy at various scales is used here to mean spectroscopic measurements in experimental conditions characterized by different observational distances. According to the instantaneous field of view (IFOV), the increase of sensor-target distance causes lower at-sensor irradiance with a consequent decrease of the signal-to-noise ratio. As a consequence, ground and/or spectral resolutions are degraded, and the probability of both spectral inhomogeneity within the pixel and ambiguity in spectral interpretation increases.
Variation of the observational distance and consequent data degradation also control the data resolution threshold for surface complexity analysis. Surface complexity manifests on different scales. The macro-complexity is due to topographic effects. It produces lightning and radiance variations among the pixels and is not considered in this paper. The meso-complexity occurs at subpixel scale, including surface texture and mixed pixels. Surface texture is the object of recent theoretical studies on the reflectance distribution function (e.g., Hapke, 2001). The reflectance of a mixed pixel is generally assumed to result in a linear combination of the spectral properties of the component materials, proportionally to their areal distribution within the pixel. In agreement with this assumption, meso-complexity is not further discussed here. The micro-complexity is intrinsic to the surface individual lithologies. It derives from the fact that the rocks are formed by many minerals with different chemistries intimately mixed together. In addition, several rock-forming minerals with similar chemistries have distinct crystal structures and a variety of grain sizes arranged into a number of rock textural patterns. All of these micro-complexity variables affect bulk-rock spectral properties at different observational scales.
The micro-complexity of rocks is commonly analyzed with laboratory spectroscopy. However, hyperspectral imaging spectrometers are supposed to be suitable for investigating rock micro-complexity from remote distances, in spite of some problems due to low signal-to-noise ratios. Imaging interferometers recently developed for remote-sensing purposes operate using the whole at sensor radiance (e.g., Meigs et al., 1997; Barducci et al., 2004), thus increasing both signal-to-noise ratios and spatial and spectral resolutions with respect to the optical dispersion systems. Such improved performances allow remotely sensed data to be reasonably compared to field and laboratory data, opening up new chances for resolving the micro-complexity of the surface.
The purpose of hyperspectral remote-data analysis is to identify and map different targets: (1) chemical species, such as Fe and Ti distribution on a planetary surface (e.g., Wilcox et al., 2005); (2) local abundances of minerals, for example, olivine and pyroxene amount in soils on Mars (e.g., Mustard et al., 2005); and (3) nature and areal distribution of rocks for local or regional studies (e.g., McSween et al., 2003). In the last case, all the factors contributing to bulk-rock spectral complexity should be taken into consideration in the choice of the end members for the remote-data classification. The most suitable end members should therefore be laboratory bulk-rock spectra, and appropriate criteria of analysis and classification of these spectra should be applied in order to synthesize and rank the diagnostic spectral characteristics.
On the other hand, accurate models for the interpretation of the bulk-rock spectral properties have not yet been fully established. Therefore, in a discussion about multiscale spectroscopy of bulk rocks, both the state of the art in understanding the spectral properties of a rock system as a whole, and the capability of remote sensing to resolve rock subtle spectral features should be considered (Sgavetti et al., 2006).
The purpose of this paper is to contribute to the general understanding of bulk-rock spectral properties for compositional analysis, classification, and mapping of data acquired in different experimental conditions. The discussion focuses on (1) the analysis of rock spectral variability, which is a component of the micro-complexity, in order to establish relationships between mineral chemistry and absorption bands in bulk-rock spectra, focusing on rock sequences that formed according to a given geologic process; this provides the geologic basis for every spectroscopic analysis, regardless of resolution and observational distance; (2) the criteria for bulk-rock spectral classification, which represents the necessary link between laboratory and remote-sensing spectroscopy; and (3) the problems arising when data from laboratory, in situ, and remote measurements are integrated; some of these problems emerged from the preliminary results of a multi-scale survey in a particularly critical area characterized by an extremely low reflectance.
The spectral range investigated here is the 0.3–2.5 micron comprising near ultraviolet, visible and near infrared (UV-VIS-NIR) regions. Within this interval, electronic transitions and vibrational overtones are particularly sensitive to the chemical composition of several rock-forming minerals (e.g., Hunt, 1980).
We measured the reflectance spectra of both slightly polished slab surfaces and particulate samples using a double-beam, double-monochromator spectrophotometer (Perkin-Elmer, Mod. λ19) equipped with a 60 mm integrating sphere internally covered by gold. The reflectance spectra were acquired between 0.35 and 2.50 µm, using a halogen lamp as radiation source. Calibration was performed with plates of Spectralon®. The measurements were performed in directional-hemispherical reflectance geometry, at room temperature and normal atmospheric pressure, with a spectral sampling of 1 nm. The spectrometer system incorporated a dual-background correction function, which automatically corrected the background absorptions. The resolution of the spectrophotometer was ±0.2 nm and ±0.8 nm for the UV/VIS and NIR ranges, respectively. The average spectra after 10 measurements per each sample were considered as representative of rock modal composition (Fig. 1A). The spectrometer was available at the spectroscopy laboratory of IFAC (Istituto di Fisica Applicata Nello Carrara), Consiglio Nazionale delle Ricerche (CNR), Florence, Italy.
In situ measurements were made using an ASD FieldSpec Pro FR portable spectrometer operating in the 0.35–2.50 µm spectral range (Sgavetti et al., 2003). Calibration was performed with plates of Spectralon®. The spectrometer was equipped with 1.4 m optic cable, 25° full angle cone of acceptance; the acquisition speed was 10 spectra/s. For each target, 1 s time of acquisition was adopted. Up to 5 acquisitions were averaged for each target within a site, and are referred to as “spectrum” in the following (Fig. 1B). Field spectra were acquired with nadir geometry on surfaces with inclination up to 10°, from an average distance of 70 cm. A set of spectra of pyroclastic deposits was measured from a distance of 30 cm using a tripod.
ASTER (advanced spaceborne thermal emission and reflection radiometer) and Hyperion daytime images used for remote spectral analyses were acquired on 19 July 2003 and were provided by the INGV (Istituto Nazionale di Geofisica e Vulcanologia), Rome division. ASTER is a multispectral imager on board the National Aeronautics and Space Administration's (NASA) Terra spacecraft (http://asterweb.jpl.nasa.gov/). It covers a wide spectral region with 14 bands from the visible to the thermal infrared with high spatial, spectral, and radiometric resolution. An additional backward-looking near-infrared band provides stereo coverage. Spatial resolution and spectral bandpass vary with wavelength 01(Table 1). Hyperion is a hyperspectral imaging spectrometer on board NASA's EO-1 spacecraft (http://eo1.gsfc.nasa.gov/Technology/Hyperion.html). Hyperion produces spectra using 226 channels covering the range 0.4–2.4 µm with a spatial resolution of 30 m per pixel 01(Table 1). The Mount Etna ASTER and Hyperion overflights were accomplished under the supervision of M.F. Buongiorno, INGV, Rome.
Images of the multispectral Thematic Mapper sensor (TM) onboard Landsat 5 spacecraft (http://landsat.usgs.gov/index.php) were also used for a portion of the metamorphic basement of central Madagascar. TM has six bands in the 0.4–2.5 µm range and one in the 10–12 µm spectral ranges, with variable spectral and spatial resolutions, as shown in 01Table 1.
MICRO-COMPLEXITY AND ROCK SPECTRAL VARIABILITY
Several parameters, such as mineral chemistry and structure, grain size, and rock texture, contribute to the surface micro-complexity and affect the spectral properties of rocks. Variation in these parameters can produce both systematic and nonsystematic rock spectral variability. In this section, the spectroscopic variability and the effects of bulk-rock parameters are analyzed for igneous and metamorphic rocks, which form large portions of Earth's crust.
Bulk-Rock Spectroscopic Variability in Genetically Related Rocks
The diagnostic absorption features occurring in mineral and rock spectra are due to energy-matter interactions within molecular groups. The frequency of the quantic transitions between energy levels is peculiar to each individual absorption process and can have a range of variation, depending on both the chemical composition and crystal structure of minerals.
Metamorphism results from physical transformations and chemical reactions involving the original protolith, strictly controlled by pressure, temperature, and time. Protolith chemical composition, grain size, texture, and the nature of the circulating fluids also control the chemistry, mineralogy, and texture of metamorphic rocks. A sequence of metamorphic mineral assemblages associated in space indicates the range of metamorphic grade through an area.
Igneous rock-forming minerals primarily derive from fractional crystallization of a melt produced by mantle fusion and/or crustal assimilation, according to known sequences of processes. Main input variables are the composition of melt, temperature, pressure, oxygen, water, and gas fugacities. The output consists of series of mineral phases, some of which are solid solutions between two end members. A magmatic rock is therefore a complex chemical system, constrained by specific physic-chemical conditions that affect its composition, crystal structure, crystallinity degree, and texture.
Rock petrology and geochemistry resulting from different geologic processes affect the rock spectral signature. Hence, the spectral variability is intrinsic to the involved rock-forming process, whatever it is. As a result, genetically related rocks can display systematic variations of spectral parameters as functions of systematic variations of petrographical and geochemical parameters.
Spectroscopic Variability of Metamorphic and Mafic IntrusiveRocks
It has been shown that vibrational overtones of the fundamental Al-OH bond dominate the wavelength interval between 2.0 and 2.5 µm in spectra of phyllosilicate from rock types belonging to different geologic settings (e.g., Swayze et al., 1992; Duke, 1994; Longhi et al., 2001, 2004). The variation of Al-OH overtone band frequency as a function of the Al/Fe ratio has been observed for micas with several different compositions (Post and Noble, 1993), including hydrothermal muscovite (Swayze et al., 1992). A negative correlation has been independently observed by different researchers between the Al-OH overtone band centers and the octahedral Al content of muscovite, occurring both in metamorphosed graywackes from Black Hills, South Dakota (Duke, 1994), and in suites of quartzites and micaschists from the pre-Paleozoic basement of central Madagascar (Longhi et al., 2000, 2004), as shown in Figure 2A. The systematic shift toward shorter wavelengths of Al-OH overtone band centers as the [VI]Al abundance increases in muscovite has been correlated with the increasing metamorphic grade. In one case, the metamorphic grade shifts from the biotite zone toward the sillimanite + K-feldspar zone (Black Hills graywacke; Duke, 1994); in the other, from the diopside zone toward the sillimanite + K-feldspar zone (central Madagascar quartzites and micaschists; Moine, 1974). The increasing grade of metamorphism is responsible for the progressive [VI]Al enrichment in muscovite, according to progressively reduced phengite substitution within Al-undersaturated rocks (Guidotti and Sassi, 1976).
In Figure 2A, the quartzite-micaschists and graywacke suites form distinct clusters in the [VI]Al abundance versus Al-OH band center space. This was attributed to differing bulk-rock compositions (Longhi et al., 2000). In more detail, two slightly different variation trends also characterize the quartzite and micaschist distributions. The bulk-rock effect was analyzed considering the product of [VI]Al (apfu: atoms per formula unit) times modal muscovite, which estimates the distribution of the spectroscopically active group within the rock. As shown in Figure 2B, the relationship between Al-OH overtone band centers and this distribution parameter is still negative (at a confidence level of 99%) in the case of quartzites, analogous to the relationship in Figure 2A. On the contrary, the two variables are not correlated for micaschists (Fig. 2C). This is consistent with the different bulk compositions, because the quartzites are composed of quartz with a small content of muscovite and tourmaline, while the micaschists consist of quartz, muscovite, biotite, and albite in various amounts.
Spectral variability was also observed for diagnostic absorptions due to electronic processes in iron-bearing minerals occurring in mafic rocks. Crystal field (CF) transitions in Fe2+, preferentially when located in distorted octahedral crystallographic sites, commonly occur in ferromagnesian minerals such as olivine and pyroxene (Burns, 1993). The wavelength position of the resulting absorption bands has been accurately measured by several researchers, using sets of minerals containing different amounts of Fe, and calibration curves have been developed. Calibration curves are useful tools for quantitative analysis of such minerals, especially for remote spectroscopy of rocks (Cloutis and Gaffey, 1991).
This well-demonstrated experimental basis was used here to study the spectral behavior of a suite of cumulates (mainly norites and anorthosites) belonging to the Bjerkreim-Sokndal Layered Intrusion (Norway). Their composition varied within each megacyclic unit and even within sublayers. The rock modal composition mainly included plagioclase, ilmenite, and Ca-poor pyroxene. The absorption bands diagnostic of pyroxene near 0.9 and 1.8 µm dominated the rock slab surface spectra. Accurate electron microprobe and Mössbauer analyses supported the attribution of the observed bands to Fe2+/M2 site crystal field transitions. Throughout the whole set of samples, absorption band minima shifted toward longer wavelengths with increasing Fe2+ content in octahedral M2 sites (Fig. 3A), which is consistent with the measurements derived from spectra of pyroxene separates from the same rock set (Pompilio, 2005a, 2005b) and in agreement with previously published data for pyroxenes (e.g., Cloutis and Gaffey, 1991). On the other hand, Fe2+ crystal field band minima showed a negative dependence on the surface distribution parameter of active Fe2+ and consequently on the pyroxene modal abundance in norites, in contrast to what was expected (Fig. 3B). This is due to the inverse relationship between the abundance of spectrally active group in pyroxene (Fe2+/M2 sites) and modal pyroxene in this suite of norites (Fig. 3C).
ROCK SPECTRAL CLASSIFICATION
The analysis and interpretation of rock spectral variability discussed in the previous section provide insights into characteristics and origin of the complex spectral properties of rocks. In order to manage this complexity properly and recognize it in remotely sensed data, appropriate methods of spectral classification within rock sets must be used.
Methodological Approaches for Rock Spectral Classification
Several methods for remote-sensing image spectra analysis and classification have been developed, based on various statistical and mathematical approaches. Two main approaches are generally followed: (1) classifying the unknown spectra by comparing them with either image-derived end members, or laboratory and field reference spectra (e.g., Mazer et al., 1988; Kruse et al., 1993); and (2) identifying and modeling absorption bands through: (a) absorption feature parameterization, with automatic detection of band positions, depths, and widths, after removal of the continuum (e.g., Clark and Roush, 1984; Kruse et al., 1988); (b) linear and nonlinear spectral mixing modeling, using either linear combinations of mineral spectra from reference libraries or theoretical calculations (e.g., Adams et al., 1986); and (c) curve-fitting by Gaussian and modified Gaussian distributions (Hawthorne and Waychunas, 1988; Sunshine et al., 1990, respectively), which provide the spectral decomposition in absorption bands actually related to absorption processes.
As the rocks are complex systems, classification criteria should combine the original Steno's facies concept, which refers to the external “appearance” of an object, with the appropriate genetic processes. In a rock spectrum, the basic “appearance” is represented by the absorption features, and the primary genetic processes are electronic transitions and vibrational motions within individual ions and cation-anion bonds (e.g., Hunt, 1980; Burns, 1993). As also discussed in the previous section, these processes are controlled by mineral chemistry and crystal structure and depend on physical and chemical conditions during rock formation, diagenesis, and weathering, according to the geologic context. Thus, a rock spectrum is the result of a complex chain of processes. To have an understanding of it, we need to both accurately disentangle the individual processes responsible for the spectral features, and infer the geologic significance involved in their association within the spectrum.
The empirical approach described here follows a number of steps. First, sets of absorption features occurring within the rock spectrum and attributable to absorption processes in individual minerals are isolated to form absorption patterns. Absorption patterns are thus related to mineral chemistry. Then the absorption patterns are combined into higher-order groups, which have been called spectrofacies (e.g., Fig. 4; Sgavetti et al., 1995; Ferrari et al., 1996). A spectrofacies represents an association of genetically related absorption patterns. It spectrally identifies a rock, i.e., a chemical system resulting from the physical-chemical conditions acting in a given geologic context at the time it formed. Within a spectrofacies, the spectral patterns are ranked according to the relative absorption intensities. The most prominent spectral patterns often represent the most abundant mineral in the rock. Occasionally, dominant patterns are due to modally subordinate phases, but again, they provide a description of the rock mineralogy, useful for the identification of rock types and geologic contexts.
Automatic or semiautomatic feature-finding algorithms included in software packages devoted to spectral analysis were used to identify and parameterize spectral bands (e.g., DIS-PEC, JPL, Grove et al., 1992; LABCALC®, Galactic Ind. Corp.; OriginLab® Corporation, www.originlab.com). These algorithms were able to identify band minima and calculate band parameters with the required accuracy. Then, the empirical, concept-driven procedure based on identification of spectrofacies was followed.
The spectrofacies criterion proved to be particularly apt for classifying laboratory rock spectra dominated by groups of absorption features mainly due to vibrational processes, as shown below for muscovite-bearing metamorphic rocks. Further criteria are required for classifying rock spectra dominated by electronic absorptions, as for mafic igneous rocks. In particular, the spectral classification of mafic effusive rocks, with some effects due to texture and syneffusive alteration, needs to be improved.
Spectral Classification of Metamorphic and Mafic Volcanic Rocks
Laboratory spectra measured on rock slabs were used for the spectral classification of a suite of rocks belonging to pre-Paleozoic metamorphic sequences of the central Madagascar basement (02Tables 2 and 033; Longhi et al., 2001) resulting in a tree-structure classification of spectrofacies. The first-order classes were based on dominant spectral patterns corresponding to the most abundant minerals or to the spectrally dominant species, and were represented by impure carbonate rocks, muscovite-bearing quartzites and micaschists, ferromagnesian rocks, and felsic rocks with abundant altered feldspar. Second-order classes, based on subordinate spectral patterns, described the spectral variability due to spectrally less active or less abundant phases. The spectral classes were used to interpret remotely acquired images of the area in which the rock samples were collected. Figure 5 is a spatial subset of TM images of central Madagascar. Band color composite 7, 4, 1 enhances a number of image units. Four of these units correspond to the first-order spectral classes listed in 03Table 3. They recur throughout the entire area investigated (not shown here), and the ground truth was provided by both the samples relative to the distinct units (Longhi et al., 2001) and existing lithologic maps (Moine, 1974). Higher-order classes, not discriminated in this image, are expected to be identified and mapped by hyperspectral sensors with proper spectral and spatial resolution and signal-to-noise ratio.
The spectral classification of basalts is relative to the Mount Etna (Italy) lava units ranging in age from 100 ka to the present (Romano, 1982). Spectral patterns, spectrofacies, and spectral classes were derived after the analysis of laboratory measurements on rock slab surfaces of lava samples (04Tables 4 and 055). Figure 6 shows the representative spectra for each class. The overall lava modal composition consists of plagioclase, clinopyroxene, olivine, and opaques, occurring as phenocrysts, embedded in a mostly microcrystalline and amorphous groundmass. In pyroclastic deposits, clinopyroxene is more abundant than plagioclase. Four first-order classes can be recognized on the basis of spectral patterns identified mainly according to the overall spectral shape. NIR absorption bands due to well-defined minerals (CFM spectral pattern in 04Table 4) are present in just one class of spectra. The very broad band centered in the NIR (WS and WM in 04Table 4), however, has a controversial interpretation. The slopes toward the infrared have positive or negative trends, probably depending on texture (Fig. 6). Second-order classes are based on several weak absorption bands due to electronic transitions in Fe-bearing minerals in the VIS region. Only the representative spectrofacies for each second-order class are reported in 05Table 5. The combined effect of texture, type, and degree of alteration, neutral dark components, and amorphous groundmass is still not fully understood. The basalt spectra are discussed in more detail in the following section.
INTEGRATING MULTISCALE DATA: PRELIMINARY RESULTS OF A SURVEY OF RECENT LAVA FLOWS ON MOUNT ETNA
Hyperion and ASTER satellite images of Mount Etna volcano were acquired in July 2003 (Buongiorno et al., 2006), after the intense activity at the end of 2002 and in the first months of 2003. During this eruptive phase, a number of lava flows were emplaced on both the northeast and south flanks of the volcano (Fig. 7). Contemporaneous remote-sensing and in situ measurements 06(Table 6) and sampling of different rock units were carried out.
Mount Etna is a middle Pleistocene to present volcanic complex. Lava composition has slightly changed through time, from intermediate to sodic, progressively high-K rock kindred, toward alkali-basalt and hawaiites (e.g., Romano, 1982; Cristofolini et al., 1987; Armienti et al., 1994).
Five main types of surface characterize the volcanic complex (Fig. 8). The most ancient lava flows are covered by thick and fresh vegetation, and hence they are not included in this discussion. The surface types considered here include: (1) historic lava flows (1600–1900 A.D.), covered by scattered bushes, lichens, and coatings of iron oxides and weathering-derived products (Fig. 8A); (2) pyroclastic deposits (tephra) consisting of ash and small-sized (centimetric) scoriae (Fig. 8B), which drape wide and flat regions within the huge ancient caldera near the top of the complex, the flanks of main craters, and the southern and eastern flanks of the volcano (Andronico et al., 2005); (3) the most recent lava flows (emplaced during 2001–2003 eruptions), characterized by piled, very coarse scoriae and ash (Fig. 8C); and (4) portions of recent lava flows consisting of reddish scoriae produced by syneffusive oxidation, which is extremely intense near the effusive centers (Fig. 8D). Factors such as ash concentration, size and shape of the scoriae, alteration, degree of weathering, and lichen coating all control the surface texture and roughness on various scales. Texture radiometric effects are expected to be superimposed on the spectral features of fresh and altered lavas.
Integrating Field and Satellite Hyperspectral Data: Image Data Classification
In situ measurements were made on lava flows emplaced during the 1792 and 2001–2003 eruptions, both on distal and proximal positions with respect to the effusive centers. Various pyroclastic deposits were also measured on the south and northeast flanks of the volcano 06(Table 6).
Four classes of field spectra were identified and are shown in Figure 8. Spectra of pyroclastic deposits (“pyroclastics”) and slightly oxidized lava (“lava”) (Figs. 8B and 8C) are comparable to each other in terms of overall reflectance and position of the reflectance peak. The main differences regard the slope toward long wavelengths and the position of minima of the very broad absorption in the near-infrared interval. Spectra of strongly oxidized lava, measured in the proximal portions of the most recent lava flows (Fig. 8D), are dominated by the signature of hematite, although features due to other components can also be recognized, as discussed in the following section. Spectra of historic lava (“altered lava” in Fig. 8A) show a shoulder at wavelengths longer than the reflectance maxima of previous spectral classes, and a flat or slightly convex profile in the 0.75 and 2.5 µm wavelength range. A weak absorption feature near 0.5 µm always occurs in all the spectral classes.
Hyperion data have 10 nm nominal spectral resolution and 30 m ground resolution, and are therefore more powerful than ASTER data for mineralogical mapping in the VIS-NIR range 01(Table 1). Hyperion images used in this paper were corrected for atmospheric effects and calibrated to reflectance (Buongiorno et al., 2006). A Hyperion color-composite image of the northeast sector of the volcano is shown in Figure 9A. Hyperion data still proved to be affected by nonremovable noise, as shown in Figures 10B–D and 11B–E. For this reason, the Spectral Angle Mapper (SAM) provided with the ENVI® software package appeared to be the most suitable classification method. SAM classification treats pixel and end-member spectra as vectors in a space with dimensionality equal to the number of bands, and determines the similarity between individual pixels and end-member spectra by calculating the angle between the vectors (e.g., Kruse et al., 1993). The end members for the classification of Mount Etna Hyperion image are represented by in situ measurements and linear combinations resulting from them (Figs. 9C). Class statistics are reported in 07Table 7. SAM accuracy is in general evaluated through confusion matrix analysis, by comparing the classification results with a ground-truth image or pixels from the image. However, the purposes of this paper, i.e., comparison of data acquired in different experimental settings, required a different evaluation approach, and SAM accuracy was evaluated by directly comparing classified data with the end-member field spectra used for the classification (Figs. 10 and 11; 08Table 8).
Materials with the composition of slightly oxidized lavas (lava) constitute most of the area (Fig. 9B). Materials represented by tephra spectrum (pyroclastics) concentrate around the major effusive centers and along the northeast fracture. Very recent lava flows (see detail in Fig. 10A) are classified as “lava” and “pyroclastics” along their axial zone, while their borders include a number of pixels classified as “pyroclastics/lava/oxidized” and as “altered” because of the presence of iron oxides and weathering products. The average spectra of three image classes are shown in Figures 10B–10D, where they are compared with the corresponding end members. The standard deviation for each class of spectra and the root mean square (RMS) error between class average spectra and relative end members are listed in 08Table 8. A quite good spectral correlation can be observed (RMS ≤ 0.005).
Volcanic products affected by different degrees of alteration are concentrically distributed around recent effusive centers. Figure 11A shows that this distribution is apparently also controlled by the downdip flow direction of the effusive products along the northeast slope of Mount Etna. Transitions from materials containing both iron oxides (oxidized) and weathering products (altered/lava, altered, and lava/oxidized), to mixtures of variously altered lava (lava, lava/oxidized) have been identified moving from the center outward. Except for the “lava/oxidized” class, for most of these materials, the correlation between average spectra of individual SAM classes and relative end members was poorer than for recent lava and tephra. In this case, average class spectra best compare with linear combinations of end-member field spectra not coinciding exactly with the end-member combinations used for the classification (Fig. 11B–E; 08Table 8).
Comparison Between Laboratory and Field Data
The comparison between laboratory and in situ measurements discussed in detail here regards recent lava (E25 in Fig. 6 and “lava” in Fig. 8C), tephra (PC5 in Fig. 6 and “pyroclastics” in Fig. 8B), and proximal oxidized lava (PC7 in Fig. 6 and “oxidized lava” in Fig. 8D). Laboratory spectra of historic lava (E19 in Fig. 6) are very similar to those of recent lava, whereas the corresponding in situ spectra (Fig. 8) are not comparable to each other. This is intrinsic to the laboratory operational method for accomplishing measurements on exposed, altered surfaces. According to alteration rim thickness and surface roughness, the alteration coating can be destroyed during rock slab preparation for laboratory measurements. The spectral features observed in laboratory reflectance spectra are interpreted through petrographic analyses.
In the NIR region, both laboratory and in situ measurements of recent lava are characterized by a slightly negative slope and a shallow concavity at wavelengths longer than 0.7 µm (Fig. 12A). We attribute this broad band to the superposition of a number of crystal field transitions in Fe2+ in mostly zoned Ca-rich pyroxene, olivine, and plagioclase (cf. “[SFO] WS VO NS” spectrofacies in 05Table 5 and spectrum E25 in Fig. 6). Alternatively, it can be interpreted as the effect of the porphyritic texture of basalt, which is responsible for the coexistence of phenocrysts embedded in a poorly crystalline, mostly glassy groundmass. Note that laboratory spectra broad-band minima are shifted toward shorter wavelengths with respect to the in situ spectrum (Fig. 12A).
In the same NIR region, the laboratory spectrum of tephra samples (Fig. 12B) has an overall flatter profile than in situ measurements, in which an absorption due to ferromagnesian minerals can possibly be recognized (SFM SFO [CFM] VO PS spectrofacies in 05Table 5 and PC5 in Fig. 6). In contrast, the in situ spectrum shows a broad band similar to that observed in recent lava spectra (WS spectral pattern in 04Table 4). Tephra spectra are characterized by a slightly positive slope, likely due to the particulate nature of this material.
In the VIS interval, laboratory spectra of recent lava slab surfaces show generally weak absorption bands at 0.473, 0.513, 0.572 µm (Fig. 12C; SFO pattern in 04Table 4; cf. E25 and E11 spectra in Fig. 6). These band minima can be attributed to spin-forbidden crystal field transitions in a number of cations occurring as trace elements within olivine and Ca-rich pyroxene, both fresh and oxidized (King and Ridley, 1987; Straub et al., 1991; Cloutis, 2002). In situ measurements of recent lava display a single absorption in this region, at 0.47 µm. However, the most significant difference between the two types of spectra is the opposite slope shape toward the UV, shown in Figure 12A.
In the same VIS range, a tephra laboratory spectrum shows two absorption bands at 0.40 and 0.52 µm (Fig. 12D), both assigned to spin-forbidden transitions in iron oxides (patterns SFO in 04Table 4). In the field spectrum, the two bands are shifted to 0.48 and 0.54 µm, respectively.
The several discrepancies observed in the laboratory and field VIS-NIR spectra of recent lava and tephra are further discussed in the following section, although they are the object of in progress studies.
The absorption features occurring within the vibrational overtone region are diagnostic of weathering products, also observed in the field.
Field and laboratory spectra of highly oxidized lava, measured in proximity to the effusive centers, are significantly similar to each other for both band positions and intensities (Fig. 13; CT CFH WM VO NS spectrofacies in 05Table 5 and PC7 in Figure 6; “oxidized lava” in Fig. 8D). They are characterized by a reflectance peak at 1.20–1.25 µm and three weak bands near 0.50 and at 0.90 µm superimposed on the strong charge transfer band in the UV (Figs. 13A, 13C, and 13D), and represent a spectral pattern diagnostic of hematite (cf. Fig. 13B). The near-infrared composite absorption feature in both spectra, although shifted to 1.6–2.0 µm, is still consistent with the presence of pyroxene (WM pattern in 04Table 4). A secondary absorption peak at 2.2 µm still relates to alteration products (VO pattern).
Two major problems arise from the analysis of the three data sets. The first is the noise level in the remote data, which required several sensor channels to remove, mostly in the short wavelength region. As a consequence, the observation of some of the iron absorption bands critical for main composition and alteration analyses was precluded. Nevertheless, a certain similarity was observed between remotely sensed and field data, and, actually, field spectra proved to be suitable end members for the classification of the remotely acquired spectra of Mount Etna basalts and pyroclastic deposits. A lower accuracy was obtained in the classification of variously altered and oxidized materials, likely consisting of intimately mixed minerals. They were probably not systematically measured in the field, emphasizing the importance of appropriate end members for any image data classification procedure.
The second problem is the peculiar shape of the laboratory basalt spectra, particularly in the VIS and toward the UV, which is in contrast to the field spectra of the same rocks. As a result, only two of the spectral classes determined through laboratory spectroscopy (represented by spectra PC5 and PC7 in Fig. 6) clearly correspond to field spectral classes (“pyroclastic” and “oxidized lava,” respectively in Figs. 8C and 8D). The CFM class of laboratory spectra in 05Table 5 also includes the ancient lava (E29 in Fig. 6), for which field spectroscopy is not available. The ancient lava spectrum shows charge transfer absorption bands in the UV (CT pattern in 04Table 4), which produce a peak of reflectance in the VIS. The pronounced negative slope at longer wavelengths is quite common in mafic rock spectra (Harloff and Arnold, 2001; Pompilio, 2005a). When considering the spectral classes represented by other spectra, such as E11, E6, E25, E19, the most impressive discrepancies with both the ancient lava and in situ spectra regard the absence of the CT pattern. These laboratory spectra are relative to basalts from both distal areas of recent flows, and historic lava flows.
In spite of the quite similar mineral and chemical composition, both differences and similarities in basalt laboratory spectra did not appear to correlate to variations in the petrographic parameters determined by standard analysis, such as the porphyritic and vesiculation indexes. One hypothesis could be that, for these extremely low-reflectance rocks, the spectroscopic effects of high amounts of opaque and amorphous phases are exceedingly relevant (Hapke, 2001), particularly if these phases are a significant constituent of the micro- and submicrocrystalline groundmass. Particularly complex interactions between compositional and geometric parameters, such as grain size, crystallinity, viewing geometry, and texture, can also be expected. Answering the questions raised by the set of rocks discussed here is of primary importance to bridge the gap between laboratory and both field and remotely acquired data, and is critical to obtain a geologically significant output from hyperspectral data analysis and classification.
Nevertheless, in situ spectra were interpreted using spectral patterns and spectrofacies identified in laboratory spectra, and some of the end members for remote data classification were selected starting from in situ spectra closely corresponding to spectral classes determined through laboratory spectral classification. The conceptual background for this procedure was that spectral patterns and spectrofacies identified in laboratory spectra provide synthetic descriptions of interaction processes in minerals and whole-rock spectral properties that can be identified also in data at different scales. The concept of spectral classifications described in this paper was primarily developed for outcropping rocks with negligible or limited surface weathering, i.e., under conditions that are unable to significantly mask the clean rock absorption features. In the case of intense surface weathering, spectra of exposed rock surfaces must be included in the classification, with consequent addition of new spectral patterns.
Within individual spectrofacies, some rock spectral features can systematically or nonsystematically vary in wavelength position or intensity. The systematic shift of specific absorption band minima is expected to be somewhat correlated with well-defined mineral chemistry variations, which in turn are related to the geologic processes responsible for the rock formation. The analysis of rock spectral variability therefore provides the geologic conceptual background for the interpretation of spectroscopic data at every scale. In the examples discussed in this paper, the wavelength shift of the diagnostic absorption bands was found to be either directly or inversely dependent on the abundance of spectroscopically active groups within the minerals. The geologic processes, either the fractional crystallization from a melt or the rock transformation by metamorphism, evolved according to those relationships. However, in rocks with a high content of minerals that can be spectrally active, the relationships between position of the absorption band minima and mineral abundance can be lacking or even reversed with respect to the band position versus spectrally active group relationship. This depends on the degree and type of the correlation existing between mineral and bulk chemistry, constrained by the evolutionary history of the geologic system. Therefore, the systematic analysis of the relationship between bulk-rock composition and spectrally active group should always be considered in bulk-rock spectral modeling.
When sets of multiscale spectroscopic data are not available for an area, the concept of spectrofacies can be used to classify mineral and bulk-rock spectra from available spectral libraries. This classification can be assisted by existing calibration curves and by improved understanding of the other factors affecting the bulk-rock spectra, such as texture and crystallinity. This can also provide the geologic rationale for the use of expert systems that run multifeature identification algorithms and decision-making processes (e.g., Clark et al., 2003), thus optimizing the lithologic analysis of remote-sensing data.
Spectroscopy at various scales has a direct application in remote data interpretation. A remote-sensing survey is an intrinsically multiscale survey, both when in situ measurements and sampling return for laboratory spectroscopy are performed concomitantly with the air/spaceborne sensor overpass, and in frontier exploration surveys, where, in the absence of ground data from the study area, remotely acquired spectra have to be interpreted through comparison with pre-existing laboratory spectral libraries. In both cases, understanding the remotely acquired spectra strongly depends on the accuracy of the relationships between spectroscopic features and rock chemical and mineralogical composition. Integrated spectroscopic and petrographic-mineralogical laboratory measurements are the primary tool to set up the scientific background for spectral features assignment, identify calibration curves, and also develop classification criteria, that, with proper scaling, can be applied to the different experimental conditions.
Spectroscopic calibration curves involving spectroscopically active groups can be established in sets of pure minerals and bulk-rock suites: two examples illustrated this issue relative to metamorphic rocks and cumulates. In both cases, the important role of the bulk-rock composition in controlling the parameters of the calibration curves and their applicability to quantitative determinations has been demonstrated. Understanding bulk-rock effects is critical, because the rock spectral variabilities record the rock-forming geologic processes, and therefore provide the geologic conceptual background for spectroscopic data interpretation.
Selection of end members for image data classification should include the concept of spectrofacies, i.e., a synthetic representation of rock spectral characteristics. We reported two examples of rock classification based on laboratory spectra, relative to metamorphic rocks and basalts. Different orders of classes were identified, potentially recognizable in remotely acquired images, according to the sensor resolution.
In the Mount Etna multiscale survey, individual field spectra and linear combinations of field spectra represented suitable end members for satellite hyperspectral data classification, except for the case of probable intimate mixtures of minerals. On the other hand, only a limited number of basalt spectral classes determined by laboratory spectroscopy were comparable to field spectral classes. The usefulness of laboratory spectra was, in this case, more related to absorption features assignment than as end members for image classification.
Basalts are rocks for which the relationship between spectral characteristics and compositional variation has not yet been systematically analyzed. Work in progress includes extended analyses of groundmass composition, crystallinity degree, vesiculation patterns and relationship with alteration, relative roles of phenocrysts, microphenocrysts, submicroscopic and amorphous phases, and the nature of alteration processes. As basalts are important constituents of Earth's surface and are inferred to be the primary components of other bodies of the solar system, understanding the spectral properties of these rocks on different scales is critical for building spectral libraries appropriate to the analysis of those surfaces.
Currently operative spaceborne hyperspectral sensors show serious problems in resolving absorption features in both the VIS and NIR ranges, at least for surfaces characterized by a very low reflectance. More generally, low signal-to-noise ratios preclude the detection of rock spectral variability, which is critical for rock identification from remote data. New technologies are therefore required for successful mineralogical and lithological surveys.
Field and Hyperion spectroscopic data acquisition was made in collaboration with INGV (Istituto Nazionale di Geofisica e Vulcanologia, Rome and Catania) and Valerio Tramutoli of the University of Catania. We are particularly thankful to Ted L. Roush for the critical reading of a first version of the manuscript. The useful comments and suggestions of two anonymous reviewers improved the manuscript. Financial support by ASI (Agenzia Spaziale Italiana) and MIUR (Ministero Italiano Università e Ricerca) is acknowledged.