Impact of Environmental Factors on the Spectral Characteristics of Lava Surfaces: Field Spectrometry of Basaltic Lava Flows on Tenerife, Canary Islands, Spain

We report on spectral reflectance measurements of basaltic lava flows on Tenerife Island, Spain. Lava flow surfaces of different ages, surface roughness and elevations were systematically measured using a field spectroradiometer operating in the range of 350–2500 nm. Surface roughness, oxidation and lichen coverage were documented at each measured site. Spectral properties vary with age and morphology of lava. Pre-historical lavas with no biological coverage show a prominent increase in spectral reflectance in the 400–760 nm range and a decrease in the 2140–2210 nm range. P¯ ahoehoe surfaces have higher reflectance values than 'a'¯ a ones and attain a maximum reflectance at wavelengths < 760 nm. Lichen-covered lavas are characterized by multiple lichen-related absorption and reflection features. We demonstrate that oxidation and lichen growth are two major factors controlling spectra of Tenerife lava surfaces and, therefore, propose an oxidation index and a lichen index to quantify surface alterations of lava flows: (1) the oxidation index is based on the increase of the slope of the spectral profile from blue to red as the field-observed oxidation level strengthens; and (2) the lichen index is based on the spectral reflectance in the 1660–1725 nm range, which proves to be highly correlated with lichen coverage documented in the field. The two spectral indices are applied to Landsat ETM+ and Hyperion imagery of the study area for mapping oxidation and lichen coverage on lava surfaces, respectively. Hyperion is shown to be capable of discriminating different volcanic surfaces, i.e., tephra vs. lava and oxidized lava vs. lichen-covered lava. Our study highlights the value of field spectroscopic measurements to aid interpretation of lava flow characterization using satellite images and of the effects of environmental factors on lava surface evolution over time, and, therefore, has the potential to contribute to the mapping as well as dating of lava surfaces.


Introduction
Volcanic eruptions pose a persistent threat to humans and their activities.Amongst the most frequent volcanic hazards are lava flows, e.g., [1,2].Lava emplacement is an important constructive geological process that creates and repaves the earth's crust and contributes to reshaping natural landscapes, but also has immediate and long lasting impacts on human activities.In order to assess the hazards and long term impacts posed by lava flows, it is vital to understand aspects such as the return period of effusive eruptions, to map the areas covered by eruptions in the past and to characterize the evolution of lava flow surfaces after emplacement.
Remote sensing techniques have a significant role to play in studying volcanic surfaces, their emplacement processes and environmental impacts.Current remote sensing studies of lava primarily focus on monitoring and mapping, e.g., [3,4].Thermal observations are often used to monitor potential and ongoing volcanic activities by detecting and characterizing hot spots, e.g., [5][6][7][8][9][10].Satellite data acquired in the visible to near infrared, thermal infrared ranges as well as radar signals have been proven useful in mapping lava flows [11][12][13].Other studies highlighted the possibility to discriminate weathered lavas and even to relatively date them [14][15][16].All these applications rely on material-specific spectral responses of volcanic surfaces, which may vary both spatially and temporally.Improving the understanding of the spectra of lava surfaces has the potential to help to better interpret the spectra of volcanic areas recorded on satellite images and, therefore, to facilitate the realization of our long-term goals, i.e., mapping and dating lava flows.Whereas several studies have investigated the factors controlling the contrasted characteristics of lava in the thermal infrared and radar wavelength ranges [14,[17][18][19], only limited attention has been dedicated to unravelling the factors controlling the large spectral contrasts observed within lava flow fields in the visible to near-infrared range.There is therefore a need to build a spectral library of lava surfaces, to find out how they behave spectrally, to characterize how this spectral behavior varies in space and time and to investigate which flow characteristics or environmental factors influence these variations.
Characterization of surface spectral reflectance by satellite remote sensing is constrained by the spectral range and resolution (i.e., number of spectral bands) as well as by the spatial resolution of the imagery.Whereas multispectral imagery can be acquired at very high spatial resolution (e.g., Pleiades, 0.5-2 m; [20]), spatial resolution of hyperspectral satellite data remains low (e.g., Hyperion, 30 m; [21,22]) and spectral mixing is thus a major issue.The spectral reflectance of lava of different compositions has also been documented using laboratory spectrometry with decimeter-size samples [15].For accessible volcanic terrains, field spectrometry offers a useful alternative approach for characterizing the spectral reflectance of contrasted lava surfaces and for documenting its spatial variation at different spatial scales [21].
Hyperspectral analysis of lava flows, especially in the visible and shortwave infrared regions, is rare and limited to a few volcanoes.Spectral reflectance of Hawaiian 'a'ā lavas have been collected in the laboratory and used to discriminate lava flows of different age groups [15].Volcanic surfaces at Mt. Etna have been documented based on field spectrometry, including the reflectance spectra of dated, fresh, altered and oxidized lava surfaces [21,23,24].All previous studies, however, were based on a limited number of spectroscopic measurements, not representing the spatial variability of spectral reflectance within and between lava flow surfaces.Earlier work also did not systematically investigate the environmental factors and surface conditions influencing the measured reflectance spectra.
As lava surfaces are long exposed to the environment, the question is whether and how the spectra of lava flows change with time.Studies of Nyamuragira volcano, for example, have shown that in humid, tropical climates, vegetation colonization has a significant impact on the spectra of aging lava flows by increasing the spectral response in the near infrared band [20,22].The spectral reflectance of lava flows is also affected by lichen cover [25,26] and chemical weathering [26][27][28][29][30]. Compared to lava spectra, lichens are characterized by higher reflectance with several minima and maxima at specific wavelengths (e.g., [31][32][33]).Lichen cover on lava would therefore alter the pattern of lava spectra [24,33].Exposure to the environment also promotes chemical weathering of lava.This is a complex process controlled by multiple factors such as climate, relief, lithology, vegetation cover, lava age as well as the surface roughness resulting from lava flow emplacement dynamics (e.g., [28,30]).Although oxidation has been reported to increase the slope of the reflectance curve in the visible part of the spectrum over time [15], how lava spectra are progressively influenced by other surface characteristics and environmental conditions (e.g., surface roughness [8,14,16]), and to what extent, and how this understanding could contribute in improving automatic mapping and relative dating of lava flow fields are important questions that remain to be investigated.
In order to contribute to the general objective of using spatio-temporal spectral variations of lava surfaces for improved mapping and relative dating, this study assesses the control of environmental factors on the spectral characteristics of lava surfaces.Specific objectives are: (1) to collect a large dataset of lava surface spectra using field spectroscopy, representing the spatial variability of spectral reflectance within and between lava flow surfaces; (2) to investigate the in situ environmental factors and surface conditions influencing the measured reflectance spectra, particularly age specific aspects such as degree of oxidation and lichen coverage, and time-independent surface roughness; and (3) to compare and analyze the correlation between field and satellite based reflectance measurements of lava surfaces.

Study Area
Tenerife is the largest among the seven volcanic islands of the Canaries located in the Atlantic Ocean off the west coast of Morocco, Africa (inset of Figure 1a).Tenerife has a sub-tropical climate with high insolation throughout the year with average temperatures varying from 19 ˝C to 25 ˝C.As a result of trade winds and topography of the island, the northern slopes at 600-1800 m a.s.l.receive most of the rainfall (annual precipitation > 500 mm), and the south and west flanks are much warmer and drier (<300 mm) [34], leading to less vegetation cover.Moreover, the central part of the island hosts the Las Cañadas caldera and Teide stratovolcano--which are located above 2000 m a.s.l.--the lower temperatures at these elevations prevent vegetation growth [35].
As a complex volcanic island, Tenerife is characterized by both effusive and explosive eruptions [36].Historical volcanic activity has been characterized by basaltic lava flows on the volcano flanks, including the Boca Cangrejo (1492), Fasnia (1705), Garachico (1706) and Chinyero (1909) lava flows, and also the Chahorra (1798) flow in the south of the caldera ( [36]; Figure 1).The latter covers several pre-historical, undated, basaltic flow units displaying a higher degree of weathering (Figure 1b).The historical basaltic lavas have either a basanite or phonolotic tephrite composition with limited crystal content [36,37].
Two field campaigns were conducted in November 2013 and September 2014, during which 43 lava flow sites of different ages and three tephra surface sites (mostly lapilli with <10% bombs) were documented with field spectrometry (Figure 1; Table A1).Measurement sites were selected in order to optimally represent the spatial variability of lava surfaces of different ages, with specific attention given to different elevations under different climatic conditions.Understandably, sites are selected only if they are accessible.Other criteria include surface homogeneity and a minimum distance of 30 m away from roads or lava flow edges.).For more information on these eruptions and how they were dated, the reader is referred to [36] and references therein.

Field Spectral Measurements
The spectroradiometer used for collecting field hyperspectral data is a portable FieldSpec 3, manufactured by US-based Analytical Spectral Devices (ASD).This instrument measures targets within the 350-2500 nm spectral range using three spectrometers, namely a VNIR (visible and near infrared) detector for 350-1000 nm and two SWIR (shortwave infrared) ones, SWIR1 for ~1000-1850 nm and SWIR2 for 1700-2500 nm [39].
The procedure for obtaining the surface reflectance of lava flows using the ASD spectroradiometer is shown in Figure 2. Twenty spots, evenly distributed within a 30-meter-side square, were measured around each site center (Figure 2b).The size of the site area corresponds to the pixel size of commonly used satellite data such as Landsat and Hyperion.The pistol grip was held at a height of 1.2 m with bare fore optic.With the inherent 25-degree field of view (FOV), this results in the scanning of a surface area of roughly 50 ˆ50 cm 2 .The reflectance at each spot was measured four times (each time averaging 40 scans), resulting in a total of 80 spectral curves for each site.At the center of the selected site, a spectralon panel was stabilized on a tripod for calibration at the start and after every five spots' measurement.All measurements were conducted within 2.5 h before and after solar noon to avoid shadow and low sunlight angle effects.This measurement strategy ensured that we obtained a representative average reflectance spectrum for each site, while also documenting the local variability within the site.This example shows that the surface of the spot scanned by the spectroradiometer pistol is 40% covered by lichens.

Quantification of Weathering
Lavas are highly vesicular and porous volcanic rocks.Surface weathering is influenced by climatic conditions and by the texture of the lava flow, which will vary between pāhoehoe and 'a'ā (e.g., [28,40]).Lava on Tenerife has experienced different forms of surface alteration, including chemical weathering (mostly oxidation), biological weathering (mostly lichen colonization) and vegetation growth (Figure 3).Recovery of vegetation (ferns, pine trees and bush species, e.g., spartocytisus supranubius and erysimum scoparium) is limited to undated flows within the Teide National Park and the lower part of the Garachico lava flow (due to a local reforestation program), and thus was not taken into account during our field measurements.Fruticose lichens (mostly stereocaulon vesuvianum) grow predominantly on the shady sides of 'a'ā lavas, especially on Garachico and occasionally on Boca Cangrejo and Chinyero.Lichens are capable of weathering lava rocks biologically and, therefore, have been previously studied from various perspectives [25,41,42].The degree of lava surface alteration was characterized in a semi-quantitative way by placing a 50 ˆ50 cm 2 plastic net consisting of 25 equal-sized grids on the lava surface and calculating the percentage of grids filled by oxidized surface and lichens (Figure 2c).The size of the plastic net is equal to the FOV of the fore optic, which allows associating the lava spectrum to the degree of oxidation and the amount of lichens.Whereas the identification of the lichen cover fraction is unambiguous, the proportion of oxidized surface should be rather interpreted as a qualitative indicator, as it is based on visual interpretation of the transition from grey/black for fresh lava surfaces, to light brown or reddish for oxidized ones.The visual assessment of the "redness" of the rock should obviously correlate with the reflectance measured in the red wavelengths, even though the plastic net method assesses the proportion of oxidized surface rather than the intensity.This independent characterization helps to assess the consistency of the field inspection of surface weathering and the spectral measurement.It also helps to understand whether surface oxidation is associated with other spectral features.
The fractions of lichen-covered and oxidized lava surfaces were visually estimated for each measured spot and averaged at the site scale for the 29 sites measured during the September 2014 campaign.For the 18 sites measured in November 2013, the level of lichen coverage and oxidation were estimated at the site scale based on field observations and photos.

Surface Roughness
Surface roughness defined by the micro-topographic variations is an important characteristic of lava surfaces, controlled by its emplacement dynamics (e.g., [18,19,43]).Roughness controls the lava surface porosity and will therefore influence the type and rate of surface weathering [44,45].Lava surface roughness has been previously correlated with emission in the thermal infrared range [14] and backscatter of radar signals [8,46].We here wish to investigate whether lava flow surface roughness also influences the spectral reflectance in the VNIR and SWIR ranges.To characterize the surface roughness at the scale of each site, the chain method was applied [47].A 30-meter-long chain was placed in E-W and S-N directions along the lava surface, in such a way that the chain follows the surface's irregularities.The measured horizontal distance (Table A1) covered by the chain decreases as surface roughness increases.The ratio of the chain length to the mean horizontal distance along the two transects was used to characterize the average surface roughness of each site: this ratio is 1 for perfectly smooth surfaces and increases with surface roughness.

Processing of Field Spectra
Field spectra acquired from a spectroradiometer are often noisy and need to be processed prior to further analysis.Typical processing practice includes jump correction, noise removal and smoothing [48].
A small jump at the VNIR-to-SWIR1 spectrometer transition at 1000 nm is observed for spectral data acquired in the field.This spectral discontinuity is attributed largely to a slight offset of FOVs of the fiber optics for each spectrometer and the fact that the signal-to-noise ratio (SNR) is lowest at connecting wavelengths [38].To correct the jump, we shifted the reflectances at the 350-1000 nm range and beyond 1000 nm half the difference toward each other [49].
The ~1400 nm and ~1900 nm and >2350 nm ranges are affected by high noise for all acquired field spectra, which is attributed to strong absorption of incident downwelling radiation by atmospheric water vapor [50].In our case, the affected ranges 1351-1459 nm, 1771-2029 nm and 2351-2500 nm were excluded from spectral analysis.Afterwards, the data was smoothed by averaging seven adjacent points to build a lava spectral library.

Spectral Indices
In order to link the ground-based reflectance spectra with lava surface oxidation and lichen coverage observed in the field, we tested spectral indices of three different types: subtraction (ref2 ´ref1), ratio (ref2/ref1) and normalized difference ((ref2 ´ref1)/(ref2 + ref1)).Both ref1 and ref2 in the formulations are the reflectances at a specific wavelength or the average reflectances in wavelength ranges, resampled over an interval of 5 nm.The wavelengths or wavelength ranges of interest are either observed from the field spectral data or taken from the literature (see Results).
By plotting the proposed spectral indices against field-observed lava surface alteration, the wavelengths and the type of spectral indices that best characterize the observed alterations can be determined.

Satellite Remote Sensing Data
Satellite remote sensing data of Tenerife that were used in this research include a Landsat ETM+ image acquired on 23 January 2001 and a Hyperion image acquired on 20 September 2003, both freely obtained through USGS Earth Explorer site.The Landsat ETM+ image covers almost the entire island except for its southernmost part while the Hyperion image only covers the western flank of Tenerife, including the Garachico, Boca Cangrejo and Chinyero lava flows.ETM+ is a multispectral radiometer mounted on the seventh Landsat mission, capable of providing imaging information of the Earth's surface in eight bands at a spatial resolution of 30 m (VNIR and SWIR bands b1-b5 and b7) [51].Hyperion is a hyperspectral imaging sensor mounted on Earth Observing-1 Satellite, scanning the Earth's surface at a spatial resolution of 30 m in 220 unique spectral bands within the range of 400-2500 nm [52].The spectral range of this spaceborne instrument is almost the same as that of the field spectroradiometer.
The spatial resolution of the two image datasets is equal to the size of the sites measured in the field, enabling a comparison of the field spectral data and the satellite image data.The interval of over 10 years between the images and field data is not seen as an obstacle to comparison, as the lava surfaces are more than 100 years old and are expected to have evolved slowly.
Both images were atmospherically corrected and georeferenced to a WebGIS-based 1:5000 topographic map of Tenerife (2009) [53] with a planimetric error <10 m.A topographic correction (C method) [54,55] was also tested with the Landsat ETM+ image but the dark surfaces in shadow zones show too little spectral contrast to return useful information and have to be masked out.The field spectra were resampled at the wavelength interval of ETM+ and Hyperion (i.e., 10 nm), and the Hyperion spectra were smoothed to diminish noise.

Lava Flows
The common features of their reflectance spectra are illustrated in Figure 4. Based on the shape of their spectral curves, lava sites on the island can be categorized into three groups.(f) absolute reflectances for selected sites.Grey background behind spectral curves are standard deviations of reflectance at each wavelength for sites C11, F03T and G04, measured from the standard deviation of the 20 spots measurement per site.Wavelengths corresponding to significant spectral features are indicated by dashed lines.The sites whose names end with a "P" or "T" were measured on pāhoehoe lava surfaces and tephra fields, respectively.All other sites were measured on 'a'ā lava surfaces.Spectral ranges centered around 1400, 1900 nm and beyond 2350 nm are removed due to high noise.
The first group includes all the sites measured on the Chahorra, Fasnia and the undated lava flows, as well as the B01P of Boca Cangrejo (Figure 4a,b).All these sites are characterized by elevations >1470 m a.s.l.(referred to as the high elevation sites).A noticeable reflectance maximum is observed around 750-770 nm for 'a'ā surfaces while this maximum shifts to shorter wavelengths for pāhoehoe surfaces (C03P and B01P).The steepest part of the reflectance curves in the visible range of the spectrum is located at 500-550 nm.In the NIR, most lava surfaces' spectra show a decrease except for a few prehistorical lava sites (U01P, U04 and U10).For wavelengths beyond 1500 nm, the reflectance curves are very flat except for a noticeable decrease in reflectance from 2130-2150 nm to 2200-2220 nm for all lava sites.The latter feature is most clear on prehistorical lava flows.
The second group consists of all the sites measured on the Garachico lava flow (Figure 4c), which are characterized by low elevations (<1200 m a.s.l., referred to as the low elevation sites) and significant lichen coverage.After a sharp increase in reflectance from around 390 nm and to about 680 nm, the reflectance spectra of lichen-covered lavas show a continuous increase in the VNIR range, reaching a maximum at 1300-1340 nm and a slope discontinuity at 1180-1200 nm.A well-defined spectral feature marked by a maximum at 1650-1670 nm associated with a minimum at 1720-1730 nm is observed for all the Garachico sites.A second spectral feature is observed at wavelengths >2000 nm, with a maximum at 2220-2235 nm and two minima at 2085-2105 nm and 2300-2315 nm.
The third group includes the sites measured on the Boca Cangrejo (except B01P) and Chinyero lava flows (Figure 4d), which are characterized by intermediate elevations (1200-1470 m a.s.l., referred to as the intermediate elevation sites) and limited lichen coverage and chemical weathering.The increase in reflectance up to around 800 nm is similar for all sites on 'a'ā surfaces of the two lava flows.The profiles vary in slope toward longer wavelengths.Sites B02, Y01 and Y02 exhibit a prominent peak around 1125-1150 nm and pronounced spikes in the SWIR2 range, particularly beyond 2000 nm, which are very likely due to noise as a result of unfavorable weather conditions.

Tephra
As discussed in the methodology, tephra refers mainly to lapilli sized pyroclasts (2-64 mm in diameter).Reflectance spectra of tephra are different from those of lavas (Figure 4e).Tephra show a local maximum around 700-750 nm but their reflectances are generally low, less than 5% (Figure 4f), with limited variation.Their flat spectra increase in the visible range and decrease in the NIR range (750-1000 nm), followed by a gradual increase in the SWIR range.

Comparison with Hyperspectral Image Data
Figure 5 shows the field-derived and Hyperion spectral curves for eight of the 18 sites covered by the Hyperion image.The Hyperion reflectance in the violet to blue range is often negative after atmospheric correction, due to the low SNR of the sensor at these wavelengths, and is therefore deleted before comparison.The comparison shows that: (1) field spectral reflectances are generally higher than Hyperion spectral reflectances except for the tephra surface (G01T); (2) large differences in reflectance are in the SWIR region for the Garachico lava sites; (3) the two spectral curves have similar shapes, notably for the Garachico lava sites; (4) the Hyperion and field derived spectral curves always show correlations >90%, except for the very flat spectral curve of the tephra surface; and (5) noise in the NIR and SWIR of Hyperion spectra tends to weaken spectral features or to create false ones, which may imply that derivation of spectral indices is less reliable.More importantly, the Hyperion data show clear differences between the various surfaces: the reflectance of tephra surfaces is up to 5%, which is lower than lava surfaces; oxidized lava surfaces tend to have flatter spectra than lichen-covered ones, which show reflectance higher than 10% in the NIR and clear lichen-related features, e.g., the decrease in reflectance between 1660 nm and 1725 nm.

Oxidation
Based on the spectral signatures of oxidized lava surfaces and existing literatures, the following wavelengths and wavelength ranges were used for defining the oxidation indices:

‚
The slopes of the reflectance curves between 400 and ~800 nm have been reported to show a systematic increase with the relative ages of lava flows as well as the ferric iron content [15].

‚
The ratio of the Landsat red band to blue band has been reported as useful for measuring the red coloration of rocks [56].We therefore considered the wavelength ranges 450-520 nm and 630-690 nm, corresponding to Landsat band 1 (blue) and band 3 (red), respectively.The central wavelengths (485 nm and 660 nm) of the two ranges were also tested.

‚
We also considered the 450-495 nm and 620-750 nm ranges, corresponding to the entire blue and red ranges of the electromagnetic spectrum, respectively, as well as their central wavelengths (475 nm and 680 nm).
Linear relationships between different types of indices for each of the above wavelength definitions and the observed lava oxidation were analyzed.The sites of the Garachico lava flow were excluded from this analysis, as there is significant lichen coverage.The regression models obtained for the different indices at site and spot scales prove to be similar for the different wavelength definitions, except for 400 nm and 800 nm wavelengths (Figure 6 and Table A2).At the site scale, subtraction oxidation indices have clearly higher adjusted R 2 values than ratio and normalized difference indices.For the 450-520 nm and 630-690 nm ranges, equivalent to Landsat bands 1 and 3, there are no differences in whether the entire range or the central wavelength is used for defining the oxidation indices.However, use of the central wavelengths for the 450-495 nm and 620-750 nm ranges, the entire red and blue ranges of the electromagnetic spectrum, slightly improves the adjusted R 2 values for all oxidation indices at the site scale.The results show that the subtraction oxidation indices using the 450-520 nm and 630-690 nm ranges (Landsat ETM+), and the central wavelengths of the entire red and blue spectrum best characterize the relationship between spectral signatures of lava surfaces and their surface oxidation.Contrasted trends between the recent (Chahorra and Chinyero) and old flows (Boca Cangrejo and the undated) suggest that that field-based estimates of oxidized surface might be biased by the overall color of the lava surface, leading to overestimation of oxidization for fresh flows and underestimation for older ones.

Lichen Cover
Based on data for the lava sites of Garachico, Boca Cangrejo and Chinyero where lichens were observed, several spectral indices were tested to estimate lichen coverage at both site and spot scales.The focus is on the spectral features of lichen-covered surfaces discussed above, in particular the increase from 680 nm to 1320 nm and the decreases from 1660 nm to 1725 nm and from 2230 nm to 2300 nm.The relationships between the subtraction, ratio and normalized difference indices using the specific wavelengths and the measured lichen coverage were found to be mostly linear except for the ratio and normalized difference indices using 680 nm and 1320 nm (Figure 7 and Table A3).Results show that indices using 1660 nm and 1725 nm are in general better than those using the other wavelengths.Although the indices using 2230 nm and 2300 nm have a very strong correlation with the lichen coverage, their capacity is overshadowed by the use of wavelengths in the SWIR region which is easily affected by atmospheric conditions.More interestingly, there are no differences between the ratio and normalized difference indices using 1660 nm and 1725 nm.As normalization helps to minimize the effects of variable irradiance (illumination) levels and to limit the index values from ´1 to 1 [57,58], we considered that the normalized difference index using 1660 nm and 1725 nm is better in characterizing the relationship between spectral signatures of lava surfaces and their surface lichen coverage.

Effect of Surface Roughness
To understand how surface roughness can be related to reflectance, we focused on the correlation between the roughness index with both the reflectance and the first derivate of reflectance at each wavelength (Figure 8).Only sites of the Chahorra lava flow which are characterized by limited surface oxidation and absence of lichen were considered for this analysis.Correlation between roughness and reflectance is generally low for all wavelengths, particularly in the NIR and SWIR regions, with higher negative correlations (<´0.5)observed for wavelengths < 500 nm, i.e., the violet and blue of the visible region.When considering roughness and the first derivative of reflectance, we observe that the correlation coefficient varies strongly with wavelengths (Figure 8b), yet for some parts of the spectrum high correlations are obtained.The highest negative correlation of ´0.806 occurs at 1244 nm and it implies that a flatter reflectance curve at 1244 nm is indicative of a smoother surface, while a steeper decrease in reflectance in this spectral range is associated with rougher lava surfaces.The correlation values in the SWIR are however highly variable, suggesting a low reliability of this estimate.On the other hand, a lower, but consistent, positive correlation is obtained between the roughness index and the first derivate of the spectra at the red to NIR transition (~720 nm).

Interpretation and Discussion
Understanding the spectral behaviors of different volcanic terrains could help to contrast and map different volcanic products using satellite imagery [21,24].With the spectra of lava and vegetation endmembers, spectral unmixing has previously been used to map lava surfaces of different ages and vegetation colonization [20].A thematic map of tephra distribution of Mt.Etna was created by combining the Hyperion spectra and field spectra of tephra [21].Another potential application of spectral data would be relatively dating lava flows, which requires taking more factors (e.g., oxidation, vegetation, lichen cover) into account in the modelling.Although such applications fall beyond the scope of our current contribution, highlighting the contrasted spectral features of these different volcanic surfaces is a first step towards this improved mapping of volcanic terrains.

Interpretation of Lava Spectral Curves
In general, Tenerife's lavas show the typical spectra of basalt (see C11 & U03 vs. basalt in Figure 9).Lava flows with no lichen but subjected to different levels of surface oxidation (the high elevation group) show a spectral reflection feature centered at 760 nm, which was previously identified at 800 nm for Hawaiian lava flows [15].The general increase in the visible range can be attributed to the charge transfer in ferric ion (Fe 3+ ) ( [15,59]; Figure 9).The three absorption features in the VNIR range of hematite, centered at 530 nm, 630 nm and 880 nm [60][61][62], are weakly expressed in Tenerife lavas (Figure 9).The declining reflectance in the range beyond 760 nm is explained by ferrous iron (Fe 2+ ) electronic effects.The weak absorption feature around 2210 nm, which was also tested for its association with the observed surface oxidation, is likely owing to formation of clay [15].Comparison between field-derived spectra for Tenerife lavas with lab-measurement of samples for basalt, hematite and lichen.The basalt (igneous-mafic-fine-basal2) and hematite (oxide-none-coarse-o01a) spectra are obtained from the ASTER spectral library [63].The lichen spectrum shows typical spectral features of Stereocaulon [33].Spectral features for lichen and lavas, and for hematite, are indicated by dashed lines and arrows, respectively.
When comparing lava flows of different ages, we observed that for older undated flows there is a higher reflectance at 760 nm and also a steeper ascending gradient from 400 nm to 760 nm (Figures 1a and 4a,b,f).This agrees well with the observation made of the seven Hawaiian 'a'ā lava flows of various age groups [15].In addition, the decrease from 2130-2150 nm to 2200-2220 nm is much stronger for the undated flows than for the Chahorra flow (erupted in 1798), which might be explained by accumulation of clay on the lava surfaces over time.
The lava surfaces which are characterized by significant lichen coverage (the low elevation sites, Figure 4c) have a spectral reflectance that is strongly influenced by the lichen (Figure 9).Published lichen reflectance datasets highlight minimum values at 660-675 nm (agreeing with chlorophyll absorption maxima at 642 nm and 662 nm), 1160-1200 nm (probably because of absorption by water) and 2065-2110 nm (probably due to absorption by cellulose), and maximum values at 1265-1385 nm, 1655-1675 nm and 2230-2250 nm (e.g., [33]; Figure 9), as observed for the Garachico lavas.In addition, the 1720-1730 nm absorption in our lichen-covered lava spectra (likely because of the presence of cellulose) agrees with the minimum at 1720 nm for lichen-affected rock spectra [32] and at 1730-1790 nm measured for lichen spectra [33], in spite of a slight shift to the left of the latter (Figure 9).The absorption at 2300-2315 nm is characteristic of some lichen species [33] or might be explained by the presence of clay minerals containing magnesium hydroxide (Mg(OH) 2 ) [64].
Lavas of the intermediate elevation sites (Figure 4d) seem to be characterized by spectral features found at both high and low elevation lava sites.They show a maximum reflectance at about 760 nm and a decrease from 2140 nm to 2210 nm, like reflectance recorded at high elevation sites, and an increase in the NIR range, which was also observed in the low elevation sites.Accordingly, we may conclude that: (1) lava spectra at high elevations are mainly affected by chemical weathering; (2) lava spectra at low elevations mostly controlled by lichen and its associated biological weathering; and (3) lava spectra at intermediate elevations are jointly influenced by both effects.
Tephra (Figure 4e) measured on Tenerife range from 3 mm to 10 mm in diameter and from dark grey to black in color.Reflectances of the tephra surfaces are low and the recorded spectra are flat, which are similar to those described for the tephra on Mt.Etna [21].These featureless spectra are significantly different from the lava surfaces measured on Tenerife.

Comparison with Spaceborn Sensor Data
Comparison between Tenerife's field and Hyperion spectra shows that the differences of tephra spectra are more prominent than that of lava spectra, agreeing with the findings in the study of Mt.Etna [21].Our results, showing larger reflectance gaps in the SWIR region than in the visible part of the spectrum, confirm the report of much higher errors in the Hyperion SWIR range as compared to its visible range [24].This is attributed to the fact that photons at longer wavelengths carry less energy than at shorter wavelengths, and thus more likely suffer from atmospheric disturbance.In addition, the Hyperion sensor has low SNRs (VNIR SNR < 200 and SWIR SNR < 100): it is thus more difficult to acquire spectral information of low reflectance materials [24,65], such as lava in volcanic terrains.Furthermore, the occurrence of negative reflectance values for wavelengths shorter than 500 nm after atmospheric correction also reduces the applicability of some of the proposed spectral indices.In spite of these shortcomings, the Hyperion spectra still present well-marked spectral features that are similar to what is observed in the field spectral data.This allows differentiating volcanic surfaces, e.g., lapillus from lava, and oxidized lava from lichen-covered lava.

Effect of Weathering
The spectral evolution of lava flow surfaces associated with iron oxidation or lichen cover is reflected by contrasted features in the spectra of lavas at different wavelengths (Figure 3 and Table A1).These spectral features are dependent not only on time but also on climate factors controlling the type and rate of alteration.Quantitative analysis was performed on remote sensing and ground-based data to establish the link between field-observed surface alterations and spectral features, which could be used to characterize lava surfaces.

Oxidation Index and Its Application to Landsat ETM+ Image
Although our field observation of oxidization using the net method is only semi-quantitative, the observations based on surface color changes, from fresh dark black to weathered reddish-brown, seem to be a reasonable proxy for ferric iron content of lava surface [15,66].The contrasted relationships obtained for recent and older flows in Figure 5 suggest, however, that the field-based oxidation estimate might have been influenced by the local context and created a bias when comparing different flows.The subtraction oxidation index using the wavelength ranges corresponding to the Landsat ETM+ red and blue bands was shown to be useful for estimating the oxidation level of lava surfaces.
Mapping the level of oxidation across an entire flow field is of interest to discriminate flow units of different ages or to identify the contrast in weathering which is associated with the primary oxidation during the emplacement of varying weathering rate.The subtraction oxidation index using the red (630-690 nm) and blue (450-520 nm) bands of the Landsat ETM+ is highly correlated with the same oxidation indices extracted from the field spectra (Figure 10).This suggests that the spectral characteristics generated from a Landsat image pixel are properly validated by the field measurements.Figure 10 highlights the large contrast in oxidation level between the Chahorra and the undated flows located in the caldera.The age variation is, however, not systematic: the Chinyero and Boca Cangrejo lava flows, situated at lower elevations, display a higher level of oxidation than the Chahorra, irrespective of their age.The correspondence of field-and satellite-derived oxidation indices suggests that it is possible to map the oxidation of lava surfaces from multispectral imagery (Figure 11), although the oxidation index has not been validated with independent chemical analyses but only with visual description of the surface.The map of the oxidation index derived from the Landsat ETM+ data shows a clear boundary between the Chahorra lava flow and the surrounding undated lava surfaces, as well as zones with various levels of oxidation within the pre-historical lava flows.Although the Chahorra and its surrounding lava surfaces are part of the same environment, they display remarkably different oxidation levels that in this case correlate positively with the ages of these lavas.The oxidation index is also sensitive to slope orientations, especially in shadow zones (Figure 11b).In an attempt to minimize the topographic effect, the ratio and normalization oxidation indices using the Landsat ETM+ bands were tested but they did not succeed in reducing the effect of shadows.Further research is required to investigate whether part of this effect might be related to real oxidation variation with slope gradient and orientation, or how a slope-independent oxidation index could be established within these low reflectance surfaces.

Lichen Index and Its Application to the Hyperion Image
The lichen index defined with field measurements has been calculated for the 17 lava sites covered by the Hyperion image.Figure 12 illustrates the linear relationship between lichen coverage estimates derived from field spectra and from smoothed Hyperion spectra.Only four sites with very low lichen coverage (B04,Y03, Y04, and Y07) do not follow this trend.These four outliers cannot be attributed to site heterogeneity-as their sites are all located on homogeneous lava surfaces-but is very likely due to the low SNRs of the Hyperion sensor, particularly in the SWIR wavelengths.The spectral feature in the rather narrow 1660-1725 nm range for defining the lichen index might therefore be less clear in the spectral profile, especially for low lichen coverage.As lichen presence was very limited on the Boca Cangrejo and Chinyero lava flows due to their younger age, or higher elevations and/or drier climate than the Chahorra, the mapping of lichen coverage using Hyperion was applied to the Garachico flow only.
In order to minimize non-lichen effects, vegetated and urbanized areas on the Garachico lava flow were excluded.For comparison, the lichen coverage was modelled and mapped from the Hyperion image using a calibrated lichen index (Figure 13d) based on the relationship between field As lichen presence was very limited on the Boca Cangrejo and Chinyero lava flows due to their younger age, or higher elevations and/or drier climate than the Chahorra, the mapping of lichen coverage using Hyperion was applied to the Garachico flow only.
In order to minimize non-lichen effects, vegetated and urbanized areas on the Garachico lava flow were excluded.For comparison, the lichen coverage was modelled and mapped from the Hyperion image using a calibrated lichen index (Figure 13d) based on the relationship between field spectra and Hyperion derived indices shown in Figure 12.The lichen coverage may also be estimated (Figure 13e) using the high positive correlation between lichen coverage estimated in the field and elevation (Figure A strong linear relationship is obtained when considering only the Garachico sites, because in this case, interacting factors such as differences in age and varying slope orientation are excluded.The lichen estimation with the calibrated Hyperion lichen index (Figure 13d) shows that: (1) lichen coverage decreases with increasing elevation from 550 m to 1200 m a.s.l; (2) lichen coverage is lower on the steep slopes below 550 m a.s.l.; (3) lava surfaces close to the eruptive cone, where lava flow is partially covered by terphra, have lower lichen coverage; and (4) an oblique linear pattern within the index map an artefact in the Hyperion image.Figure 13e also shows that an extrapolation of the lichen coverage based purely on elevation is over-simplistic, as the elevation-lichen relationship in Figure 13a is probably only valid for the 1000-1200 m elevation range, where field measurements were acquired.
Spatial contrast in lichen coverage can be explained by surface conditions and climate differences (Figure 13b).Despite of high rainfall, the lava areas close to the volcanic cone (>1200 m a.s.l.) have low lichen coverage, which is due largely to the prominent presence of tephra covering the lava surface.The significant amount of lichen at the intermediate elevations (550-1200 m. a.s.l.) is likely attributable to the favorable temperature and precipitation conditions, which have not yet allowed a full vegetation recovery.As for the low lichen coverage at elevations <550 m, this can be mostly attributed to the steep and dissected topography as well as potential human impact on the northwestern coast of Tenerife.This spatially-explicit estimation of lichens, as demonstrated in this study, creates new opportunities for characterizing the spatio-temporal pattern of biological colonization of newly emplaced lava flows.

Effect of Surface Roughness
In Figure 4a, it is shown that the wavelength corresponding to the reflectance peak in the visible part of the spectrum for the lava surfaces affected only by oxidation is shifted to shorter wavelengths for pāhoehoe lavas (B01P and C03P), as compared to 'a'ā flows.This is confirmed by the moderate positive correlation between the slope of the reflectance curve and the roughness index measured in the field for the 640-730 nm wavelengths (Figure 8b).
In addition, the absolute reflectance in the visible region is higher for the pāhoehoe sites than for the 'a'ā sites, especially in the shortest wavelengths (Figure 4f).This is confirmed by the negative correlation between the absolute reflectance in the visible region (particularly in the violet to blue (350-500 nm) range) and lava surface roughness (Figure 8a).Fresh pāhoehoe lava surfaces (sites B01P and C03P) are very smooth, have a blue tint to some extent, and are shiny, thus reflecting more energy over the entire spectral range but more specifically in the violet to blue range (Figures 2a and 3h).Crystallization has been reported to alter the VNIR and mid-infrared spectral range of basalt (e.g., [67]) but the crystal content of Tenerife lavas is generally low [36,37] so the effect is limited.Loss of the quenched shiny surface through alteration might reduce this contrast between pāhoehoe and 'a'ā surfaces, but sufficient pāhoehoe surfaces could not be found on different flows to test this systematically.Finally, a high negative correlation of more than 0.8 was found between the first-order derivative of reflectance and the roughness (Figure 8b) but the interpretation of this result is not clear.
One might question the reliability of the roughness measured by the chain method.Unlike soil, lava surfaces (particularly 'a'ā) are rather rough, with scattered blocks of different sizes up to one meter, which makes surface roughness vary greatly from site to site, with also high directional variation and huge differences in scales (from millimeter to meter scale).This measurement strategy is however the only practicable method that could be adopted in such a physically challenging environment: its consistent application makes our results comparable with each other.

Conclusions
This study presents the results of spectral analyses of in situ spectral data of basaltic lava surfaces on Tenerife Island.Our main conclusions are as follows:

‚
Chemical weathering in the form of oxidation and biological weathering as a result of lichen growth can affect the spectra of lavas considerably.

‚
Bare lava surfaces with no lichen show low and flat reflectance spectra, increasing in the visible and decreasing in the infrared part of the spectrum.This trend tends to be more expressed for older lava flows which are affected by significant surface oxidation.Particularly, the slope of the reflectance curve from the blue to the red range of the spectrum increases as oxidation becomes stronger, a spectral feature which is useful for defining an oxidation index that is able to relatively quantify the oxidation degree of lava surfaces.

‚
Lavas by lichens have spectral features that are similar to those of lichens.Wavelengths (1660 nm and 1725 nm) sensitive to lichen colonization have been identified and used to develop an NDVI-like lichen index.A strong linear relationship between this index and lichen coverage recorded in the field has been observed, which might be used for spatial-temporal analysis of lichen cover development.
‚ Surface roughness has a limited but significant influence on the spectra of lavas.Smoother lava surfaces have higher reflectance than rough ones, especially at short wavelengths, and their maximum reflectance value is found at shorter wavelengths.Roughness is therefore positively correlated with the first-order derivative of reflectance at 720 nm.It requires further exploration to find out whether the visible and NIR bands of multispectral and hyperspectral images allow discriminating pāhoehoe from 'a'ā lava surfaces.
Our study presents the first systematic spectroscopic survey of lava flows coupled with field-based characterization of surface oxidation, lichen coverage and roughness.Spectral properties of lava surfaces are quantitatively analyzed and important spectra-controlling factors are identified, contributing to the knowledge of spectral evolution of lava.
We also illustrate the use of the oxidation and lichen indices generated from Landsat ETM+ and Hyperion imagery, respectively, for mapping oxidation and lichen coverage on lava surfaces.These approaches could help in the future in mapping complex lava flow field.Further work on quantifying the degree of chemical alteration by oxidation or other processes using laboratory analysis might be beneficial to better understand how these processes affect spectral signatures of lava surfaces.Spatial quantification of chemical and biological weathering of lava flows using satellite imagery is useful for contrasting environmental and time controls on lava spectra evolution.It therefore might contribute to using lava flow spectra as a reliable indicator of relative lava flow age.

Figure 1 .
Figure 1.Location of studied lava flows on Tenerife Island and distribution of field spectroscopic measurements sites.A shaded relief image derived from a 25 m resolution DEM is overlaid with contour lines with an interval of 50 m and GIS data of roads and municipal boundaries (GIS data are freely available at [38]).A detailed description of these sites is given in Table A1.(a) Overview of the island; (b) Chahorra flow and surrounding sites on undated flows; (c) Garachico flow; (d) Chinyero (top) and Boca Cangrejo (bottom) flows; (e) Fasnia (top) and Siete Fuentes (bottom) lava flows and tephra deposit (sites F03T and F04T).For more information on these eruptions and how they were dated, the reader is referred to[36] and references therein.

Figure 2 .
Figure 2. Illustration of the field measurement strategy: (a) spectral measurement using an ASD field spectroradiometer, an IBM laptop, a spectralon panel fixed on a tripod and a handheld GPS placed on the ground; (b) distribution of spots (green squares) around the center of the site on the lava flow surface; (c) quantifying the surface oxidation and lichen coverage with the plastic net method.This example shows that the surface of the spot scanned by the spectroradiometer pistol is 40% covered by lichens.

Figure 3 .
Figure 3. Illustration of different lava surfaces on Tenerife: (a) fresh 'a'ā lavas on the Chahorra flow at site C05; (b) a large view of the Chahorra lava flow at site C07; (c) undated oxidized 'a'ā lavas in caldera at site U10; (d) a large view of the undated surface at site U10; (e) comparison between the Chahorra lava (dark color in the background) and the undated lava (foreground); (f) partly oxidized 'a'ā lavas with few lichens on Chinyero at site Y07; (g) lichen-covered 'a'ā lavas on Garachico at site G06; (h) patch of fresh pāhoehoe lavas on Chahorra at site C03P; and (i) undated oxidized pāhoehoe lavas in caldera at site U06P.A 15 cm long pencil in (a), (c), (f) and (h) and humans in (d) and (i) are placed for visual comparison of scale.

Figure 4 .
Figure 4. Field reflectance spectra measured at different volcanic sites on Tenerife (vertically offset for clarity).Each curve is the average of ~80 measurements in 20 spots within a 30 ˆ30 m 2 square.(a) Lava sites on the Chahorra (C), Fasnia (F) and Boca Cangrejo (B) lava flows; (b) lava sites on prehistorical lava flows in the caldera (U); (c) lava sites on the Garachico (G) lava flow; (d) lava sites on Boca Cangrejo (B) and Chinyero (Y) lava flows; (e) tephra deposit sites on the Garachico and Fasnia cones;(f) absolute reflectances for selected sites.Grey background behind spectral curves are standard deviations of reflectance at each wavelength for sites C11, F03T and G04, measured from the standard deviation of the 20 spots measurement per site.Wavelengths corresponding to significant spectral features are indicated by dashed lines.The sites whose names end with a "P" or "T" were measured on pāhoehoe lava surfaces and tephra fields, respectively.All other sites were measured on 'a'ā lava surfaces.Spectral ranges centered around 1400, 1900 nm and beyond 2350 nm are removed due to high noise.

Figure 5 .
Figure 5. (a) Hyperion image of Tenerife acquired on 10 September 2003, showing the extent of lava flows and measured sites indicated by their names; (b) Graphs comparing field spectra resampled to the wavelength of the Hyperion sensor (red curves) and smoothed Hyperion spectra (green curves) for selected sites including one tephra site on the Garachico cone (G01T), three Garachico lava sites characterized by lichen cover (G05-G07), one Boca Cangrejo lava site (B03) and three Chinyero lava sites characterized mainly by intermediate level of surface oxidation (Y04, Y05 and Y07).The correlations between the two spectral curves are provided below the site labels.

Figure 6 .
Figure 6.Relationship between lava surface oxidation estimated in the field and different oxidation indices (see Table A2) at site and spot scales.(a) (ref 800 ´ref 400 ) vs. oxidation at site scale; (b) (ref 800 ´ref 400 ) vs. oxidation at spot scale; (c) ((ref 630´690 ´ref 450´520 )/(ref 630´690 + ref 450´520 )) vs. oxidation at site scale; (d) ((ref 630´690 ´ref 450´520 )/(ref 630´690 + ref 450´520 )) vs. oxidation at spot scale; (e) (ref 630´690 ´ref 450´520 ) vs. oxidation at site scale; (f) (ref 630´690 ´ref 450´520 ) vs. oxidation at spot scale; (g) (ref 685 ´ref 475 ) vs. oxidation at site scale; (h) (ref 685 ´ref 475 ) vs. oxidation at spot scale.At site scale (a,c,e and g), data for 35 sites is plotted including lava sites measured on Boca Cangrejo, Chahorra, Chinyero, Fasnia and undated lava flows in 2013 and 2014; at spot scale (b,d,f and h), there are 415 spots belonging to the sites measured on Boca Cangrejo, Chahorra, Chinyero and undated lava flows in 2014.Only sites and spots with no or limited lichen coverage were used to constrain the relationships.Standard deviations were calculated for the field oxidation observations and spectral indices for each site.The solid line and dashed line at different slopes represent the trend lines for the young and old lava sites, respectively.

Figure 7 .
Figure 7. Relationship between lichen coverage observed in the field and the proposed lichen indices at site and spot scales.(a) (ref 1660 /ref 1725 ) vs. lichen coverage at site scale; (b) (ref 1660 /ref 1725 ) vs. lichen coverage at spot scale; (c) ((ref 1660 ´ref 1725 )/(ref 1660 + ref 1725 )) vs. lichen coverage at site scale; (d) ((ref 1660 ´ref 1725 )/(ref 1660 + ref 1725 )) vs. lichen coverage at spot scale.The site scale graphs include the 19 lava sites measured on the Garachico, Boca Cangrejo and Chinyero lava flows in 2013 and 2014, and the spot graph includes 260 spots belonging to the sites measured in 2014 only.Standard deviations for lichen coverage and spectral indices were calculated based on the measurements at the spot scale for each site.In order to enhance the visibility, one tenth of the standard deviations were shown for the spectral indices.The large variations for the spectral indices are because 1660 nm and 1725 nm are in the SWIR range that is easily affected by atmospheric conditions and errors propagate through multiple mathematical operations.

Figure 8 .
Figure 8.(a) Correlation between surface roughness and reflectance at each wavelength; (b) Correlation between surface roughness and the first derivate of reflectance at each wavelength.

Figure 9 .
Figure 9.Comparison between field-derived spectra for Tenerife lavas with lab-measurement of samples for basalt, hematite and lichen.The basalt (igneous-mafic-fine-basal2) and hematite (oxide-none-coarse-o01a) spectra are obtained from the ASTER spectral library[63].The lichen spectrum shows typical spectral features of Stereocaulon[33].Spectral features for lichen and lavas, and for hematite, are indicated by dashed lines and arrows, respectively.

Figure 10 .
Figure 10. between the subtraction oxidation index (ref 630´690 ´ref 450´520 ) derived from the Landsat ETM+ image bands 1 and 3 and from the field spectral data, for 33 lava sites measured on the Boca Cangrejo, Chahorra, Chinyero and prehistorical lava flows.The grey dashed line is the 1:1 line.

Figure 11 .
Figure 11.Chahorra lava flow (1798) shown on the Landsat ETM+ image displayed in true color (R:G:B = b3:b2:b1); (b) Oxidation mapped from the Landsat ETM+ image using the relationship between oxidized surface and oxidation index shown in Figure 7c, calibrated based on the relationship between ETM+ and field spectra derived indices shown in Figure 10.The Chahorra lava flow is discriminated from its surrounding undated lava surfaces due to the contrast in oxidation level.Areas affected by steep topographic terrains and variable slope orientations are masked out (shown in grey) in the figure.

18 Figure 11 .
Figure 11.(a) Chahorra lava flow (1798) shown on the Landsat ETM+ image displayed in true color (R:G:B = b3:b2:b1); (b) Oxidation mapped from the Landsat ETM+ image using the relationship between oxidized surface and oxidation index shown in Figure 7c, calibrated based on the relationship between ETM+ and field spectra derived indices shown in Figure 10.The Chahorra lava flow is discriminated from its surrounding undated lava surfaces due to the contrast in oxidation level.Areas affected by steep topographic terrains and variable slope orientations are masked out (shown in grey) in the figure.

Figure 12 .
Figure 12.Comparison between Hyperion-derived lichen index and field-derived lichen index ((ref − ref )/(ref + ref )).The linear regression model was defined without the four outliers indicated by black circles.The grey dashed line is the 1:1 line.

Figure 12 .
Figure 12.Comparison between Hyperion-derived lichen index and field-derived lichen index (pref 1660 ´ref 1725 q { pref 1660 `ref 1725 q).The linear regression model was defined without the four outliers indicated by black circles.The grey dashed line is the 1:1 line.

Figure 13 .
Figure 13.Lichen colonization on lava surfaces.(a) Relationship between field-derived spectra of lichen-covered lava sites on the Boca Cangrejo, Chinyero and Garachico lava flows and their elevations.Lichens tend to grow more at elevations <1200 m a.s.l.The black solid fitting line is plotted only for the Garachico sites and the grey dashed line is fitted to all the displayed sites; (b) Scatterplot of elevation vs. temperature and rainfall for the Garachico lava flow area.The temperature and rainfall maps were obtained from the Spanish Meteorological Agency (AEMET) [34]; (c) Lichen coverage observed in the field displayed with their point size and color representing the value of lichen coverage on the lava surface; (d) Lichen coverage derived from the calibrated Hyperion lichen index (Figure 12); (e) Lichen coverage based on the relationship between lichen coverage and elevation found in Figure 13a (the linear regression model only with the Garachico sites).

Table A1 .
of sites measured on Tenerife.