Sentinel-2a Msi and Landsat 8 Oli Provide Data Continuity for Geological Remote Sensing

Sentinel-2A MSI is the Landsat-like spatial resolution (10–60 m) super-spectral instrument of the European Space Agency (ESA), aimed at additional data continuity for global land surface monitoring with Landsat and Satellite Pour l'Observation de la Terre (SPOT) missions. Several simulation studies have been conducted in the last several years to show the potential of Sentinel-2A MSI (MultiSpectral Instrument). Now that real data are available, the first confirmations of this potential and comparisons with other operational systems are being made. This paper aims at evaluating Sentinel-2A MSI band ratio products that are relevant for geological remote sensing. A Sentinel-2A MSI and a Landsat 8 OLI (Operational Land Imager) scene were processed from their respective levels L1C and L1T to level L2A (bottom of atmosphere reflectance). Then, three band ratios originally defined for Landsat TM (Thematic Mapper) were used to map mineralogy associated with a hydrothermal alteration system in southeast Spain. The results obtained with Sentinel-2A MSI were compared with those obtained with Landsat 8 OLI and a simulated Sentinel-2A MSI dataset that was used before actual data were released. Results show that the images appear similar to the human eye having a correlation of approximately 0.8 and higher, but that the associated data ranges differ significantly. The resulting products are also compared to a published geologic map of the study area, and it is shown that the resulting maps correspond with the conceptual geologic model of the epithermal deposit.


Introduction
The European Commission (EC) and European Space Agency (ESA) established the European Earth Observation programme Copernicus for monitoring of the environment.The space component of this programme consists of five missions: Sentinel 1-5 [1].These missions are to provide routine observations for operational Copernicus services and data continuity for already operational satellite systems [2].Sentinel-2 carries a super-spectral imager with 13 bands covering the Visible and Near InfraRed (VNIR) and ShortWave InfraRed (SWIR) wavelength region.The spatial resolution of these bands is 10-60 m and has a coverage between −56 and +84 degrees latitude with a 290 km swath width.The minimum revisit time at the equator is 10 days for Sentinel-2A, which will decrease to five days when both planned imagers are operational [2].When combined, the Sentinel-2A MSI (MultiSpectral Instrument) and Landsat 8 OLI (Operational Land Imager) sensors will provide a 10-30 m multi-spectral global coverage approximately every three days [3].The Sentinel-2A satellite was launched on 23 June 2015 and the first scenes were delivered a few days later [4].
The potential of Sentinel-2A MSI encompasses a wide range of applications [5].Based on simulated datasets, several data exploitation methods have been published in the preparation for this mission, with the majority comprising vegetation related studies: Leaf Area Index (LAI) [6][7][8]; chlorophyll and nitrogen [9]; biophysical parameters [10,11]; red edge position [9,12]; and vegetation state [13].Outside the green realm are evaluation of water quality [14] and mapping of coral reefs [15].The use of Sentinel-2A MSI for geological remote sensing was first evaluated by Van der Meer et al. [16], who compared the capabilities of Sentinel-2A MSI to the ASTER (Advanced Spaceborne Thermal Emission and Reflectance Radiometer) imager onboard the Terra satellite for mapping surface mineralogy associated with hydrothermal alteration systems.
In the 1980s, the Landsat MultiSpectral Scanner (MSS) was deployed for the first time to derive iron oxide maps [17].With the advent of the Landsat Thematic Mapper (TM), band ratios were used to separate argillic from non-argillic materials (using TM bands 5/7) and for mapping ferric/ferrous oxides (using TM bands 3/1).Since the year 2000, the ASTER sensor [18][19][20] has been the work horse of geologic remote sensing for mapping surface mineralogy.Cudahy and Hewson [21] published several ASTER band ratios that serve as proxies for mineralogy; some of these bands overlap with the bands of Sentinel-2A MSI (Figure 1).Van der Meer et al. [16] then also concluded that Sentinel-2A MSI bands can be used in band ratios that were originally defined with Landsat 5 TM or ASTER and are expected to provide data continuity for geological remote sensing.Mielke et al. [22] evaluated several current and next-generation satellite sensors, including Sentinel-2A MSI, for mapping the iron absorption feature depth in a band-ratio approach.They concluded that sensors like Sentinel-2A MSI are suitable for detection of Gossan with their method, and that they potentially can save costs when used for target detection prior to hyperspectral data acquisition.Van der Werff and van der Meer [23] demonstrated the capabilities of Sentinel-2A MSI for geological remote sensing by mapping iron oxide mineralogy.This was done by introducing a modelling approach on the 0.9 µm iron absorption feature using Sentinel-2A MSI's near-infrared (NIR) bands 6-9.This was, however, all done prior to the availability of real Sentinel-2A MSI data.Now that actual Sentinel-2A MSI data are available, several studies have been undertaken to confirm prior findings.The majority comprises vegetation related studies, such as crop and tree species classification [24]; forest crowns [25,26]; chlorophyll and nitrogen [27] and biomass [28], as well as applications involving water, such as mapping water bodies [29]; water quality [30] and lake shore habitats [31].
This study aims to confirm the usability of Sentinel-2A MSI for geological remote sensing, as was laid out by Van der Meer et al. [16].In their paper, a Sentinel-2A MSI dataset was simulated from a 2004 airborne hyperspectral image that covered a hydrothermal alteration area [32] near the former mining town of Rodalquilar, in the Sierra del Cabo de Gata in southeast Spain (Figure 2).They compared their simulated Sentinel-2A MSI dataset to an ASTER scene for mapping alteration mineralogy using band ratios published for ASTER [33] and Landsat 5 TM [34].In May 2016, a Sentinel-2A MSI and a Landsat 8 OLI image were acquired over the same area, in comparable weather conditions and season as when the hyperspectral dataset was acquired in 2004.This study aims to compare results of the Sentinel-2A MSI and Landsat 8 OLI imagery when applied for mapping alteration mineralogy, using band ratios developed for Landsat 5 TM [34].As the geology and environmental conditions of the area have not really changed in the last 12 years, the results of 2016 are also compared with the simulated dataset dating from the year 2004, as used by Van der Meer et al. [16].(A)) and the location of the image dataset used in this study (inset (B)).The Cabo de Gata volcanic field consists of calc-alkaline volcanic rocks (andesites and rhyolites) of late Tertiary age which have been extensively altered to form an assemblage of metamorphic minerals from high to low temperature as: silica, alunite, kaolinite, montmorillonite and chlorite.The geology, geochemistry and mineralization of the area is documented [35], while several remote sensing studies have been conducted at this site [36][37][38][39].Image modified after Arribas et al. [35] and Rytuba et al. [32] and used with permission from Van der Meer et al. [16].

Method
The datasets analysed in this study are a Sentinel-2A MSI L1C image, a Landsat 8 OLI L1T image and a simulated Sentinel-2A MSI dataset.The latter is based on a geometrically corrected hyperspectral dataset that was previously processed to surface reflectance.
The Landsat 8 OLI L1T dataset (12-bit radiometric quantization, 190 km swath width, ID: LC81990352016145LGN00, 24 May 2016 at 10:44 UTC, solar azimuth 125 • , solar elevation 66 • ) was radiometrically processed in IDL ENVI version 5.3 (Harris Corporation, Melbourne, FL, USA) [44] by applying the supplied gain and offset values.The L2 dataset was processed on demand by the U.S. Geological Survey Science Research and Development to surface reflectance, also being the hemispherical-directional reflectance factor (HDRF).The dataset was then subset to the study area delimited with Sentinel-2A MSI (leaving 269 × 135 pixels).The simulated Sentinel-2A MSI dataset used by Van der Meer et al. [16] is derived from a mosaic of HyMAP airborne hyperspectral images.The spectral coverage of HyMAP is nearly continuous in the VNIR and SWIR wavelength regions with small gaps in the middle of the 1.4 and 1.9 µm atmospheric water bands.The wavelength coverage, band width and spectral sampling interval for the 4 detectors are as follows: VIS 0.45-0.89µm, 15-16 nm, 15 nm; NIR 0.89-1.35µm, 15-16 nm, 15 nm; SWIR1 1.40-1.80µm, 15-16 nm, 13 nm; SWIR2 1.95-2.48µm, 18-20 nm, 17 nm.The spatial configuration of the instrument gives an Instantaneous Field Of View (IFOV) of 2.5 mrad along track and 2.0 mrad across track.The resulting scene that was used as a source for the simulated Sentinel-2A MSI dataset had originally a nominal pixel size of 4 × 4 m.
The HyMAP data (16 bit radiometric quantization) were acquired on 18 May 2004 by the German Aerospace Center (DLR) and processed to reflectance-at-surface values at a 4 m spatial sampling interval.This included geometric correction with PARGE (ReSe Applications Schläpfer, Zürich, Switzerland) [45] and a 5 m interval Digital Elevation Model (DEM), followed by an atmospheric correction with the ATCOR 4 (Atmospheric CORrection) software (ReSe Applications Schläpfer, Zürich, Switzerland) [46].The HyMAP data were then spectrally convolved to the lower resolutions of the Sentinel-2A MSI bands using spectral response functions provided by ESA (see also [47]).After the spectral resampling, the data were spatially resampled from the original 4 m pixel resolution to the 20 m resolution of Sentinel-2A MSI bands.This was done following Kruse and Perry [48] using pixel aggregation (neighbourhood averaging).The subset of the study area was eventually covered by 402 × 201 pixels of 20 × 20 m.
For a statistical comparison, all three datasets were once more resampled to a 10 m grid covering 800 × 400 pixels.The geometric correction that comes with the images is not sufficiently accurate to allow a direct comparison.Visual inspection showed an offset of multiple pixels.The geometric correction was, however, found to be sufficiently precise, as hardly any internal deformation was observed.For this reason, it was decided not to follow an image-to-image matching approach as, for example, published by Yan et al. [3].Instead, the low accuracy was manually corrected by shifting the Landsat 8 OLI image three pixels north (y-direction).The simulated dataset, based on a parametric geocoding, was shifted four pixels east (x-direction) and four pixels south (y-direction) in order to overlay the Sentinel-2A MSI image.The resulting datasets were subset to only include the pixels encompassed in all three datasets, leaving 790 × 390 pixels for comparative analysis.The weather conditions in the days before each data acquisition are in Table 1.Band ratios serve as proxies for mineral assemblages or individual mineral groups.A summary of Landsat 5 TM and ASTER band ratios suitable for use with Sentinel-2A MSI and Landsat 8 OLI is provided in Table 2.We focused on the band ratios originally defined for Landsat 5 TM by Sabins [34], which require bands 2 (Blue), 4 (Red), 6 (SWIR1) and 7 (SWIR2) of Landsat 8 OLI and bands 2, 4, 11 and 12 of Sentinel-2A MSI.In addition, Landsat 8 OLI band 5 (NIR) and Sentinel-2A MSI band 8a were needed for masking vegetation.First, the correlation between comparable bands of the Sentinel-2A MSI, Landsat 8 OLI and the simulated datasets were determined and interpreted with the help of scatterplots.Then, the Normalized Difference Vegetation Index (NDVI) was calculated for each dataset to prepare the masking of pixels with a vegetation cover.Although it was stated above that the environmental conditions of the area have not changed, it is likely that differences in vegetation cover do occur between 2004 and 2016.A vegetation index such as NDVI will, as a result, not be uniform for all scenes.It appeared that not only surface cover had changed over time, but also the range of values appeared to differ between the Sentinel-2A MSI and Landsat 8 OLI images acquired only three days apart.As no uniform thresholds could be defined, three vegetation masks were eventually created based on visual analysis, resulting in the following value ranges: 0-0.25 for Sentinel-2A MSI and Landsat 8 OLI and 0-0.4 for the simulated Sentinel-2A MSI dataset.The three resulting masks were subsequently fused and uniformly applied to all image datasets.Three of the band ratios defined in Table 2 were selected for mapping the alteration mineralogy around the town of Rodalquilar: A combination that highlights hydroxyl bearing alteration (TM 5/7, Landsat 8 OLI 6/7 and Sentinel-2A MSI 11/12), a combination that highlights iron oxides in general (TM 3/1, Landsat 8 OLI 4/2 and Sentinel-2A MSI 4/2), and a combination that highlights ferrous iron oxides specifically (TM 3/5, Landsat 8 OLI 4/6 and Sentinel-2A MSI 4/11).The correlation between these products were also analysed in scatterplots, which are shown in Section 3. Lastly, the mineral alteration map obtained with Sentinel-2A MSI is compared with a published geological map [50].

Results
The comparison of the three datasets is shown as images and as scatterplots.In the scatterplots, the sensors are indicated with S2A (Sentinel-2A MSI), LS8 (Landsat 8 OLI) and Sim (simulated Sentinel-2A MSI).Sentinel-2A MSI and Landsat 8 OLI results show both atmospherically corrected (L2A) datasets (indicated as bottom-of-atmosphere reflectance (BoA)) as well as the original data values of the L1T and L1C datasets (indicated as top-of-atmosphere reflectance (ToA)).The results of the simulated Sentinel-2A MSI dataset are only available in surface reflectance values.
Figure 4 shows the NDVI image results.Visually, a comparable spatial pattern can be seen in the Sentinel-2A MSI and Landsat 8 OLI results dating from May 2016.The simulated dataset, dating from May 2004, shows overall a similar pattern, but also shows a couple of differences that appear to be predominantly associated with agricultural fields (particularly in the northeast) and a mountain ridge (in the northwest).Although the patterns look similar in a 2%-98% image stretch, the range of NDVI values does differ considerably.The simulated dataset has the narrowest range (0.26-0.49), while Landsat 8 OLI has the widest range (0.06-0.61).The same is reflected in the scatterplots resulting from the NDVI products, which are shown in Figure 5.The correlation between Sentinel-2A MSI and Landsat 8 OLI is 0.85 for the bottom-of-atmosphere datasets (centre figure, S2a-BoA and LS8-BoA) and 0.82 for the top-of-atmosphere datasets (right figure, S2a-ToA and LS8-ToA).There is significantly more scatter between the bottom-of-atmosphere Sentinel-2A MSI dataset and the simulated Sentinel-2A MSI dataset (left figure, S2A-BoA and Sim).The high NDVI values that are associated with agricultural fields in the northeast of this scene even create a second population in the feature space, and makes the correlation drop to 0.40.   Figure 6 shows the image processing results for mapping mineral alteration with three band ratios for mapping mineralogy.When shown as colour composites, the three datasets are easily comparable to the human eye.Especially the ratio TM 5/7 (SWIR1/SWIR2, shown in red) shows similar patterns in all images.While the range of the Landsat 8 OLI ratio values is significantly higher (1.33-1.81)than both of the others, the data range for the two products based on the Sentinel-2A MSI band definitions (1.12-1.65 and 1.12-1.65)overlaps to a large extent.As a contrast, the spatial pattern for the ratio TM 3/1 (Red/Blue, shown in green) seems similar to the human eye for Landsat 8 OLI and Sentinel-2A MSI, but its presence in the simulated dataset is weaker.Between all three datasets, the data ranges also seem further apart.The ratio TM 3/5 (Red/SWIR1, shown in blue) shows similar data ranges for the two products based on Sentinel-2A MSI bands (0.40-0.72 and 0.40-0.72)and a significantly lower range for the Landsat 8 OLI product (0.27-0.43).
This result can be explained when looking at the scatterplots of these datasets in Figure 7.The correlation between the values of the Landsat 8 OLI and Sentinel-2A MSI bottom-of-atmosphere products is 0.95 for the TM 5/7 ratio (top row), 0.56 for the TM 3/1 ratio (centre) and 0.87 for the TM 3/5 ratio (bottom row).The correlation of the top-of-atmosphere reflectance products is higher than the bottom-of-atmosphere reflectance products (0.87 and 0.84 versus 0.56 and 0.78).The correlation between the Sentinel-2A MSI results of 2016 and the simulated dataset of 2004 (in the left column) is overall significantly lower, with correlations between 0.48 and 0.62.While the correlation is lower, it can however be observed that data ranges for band ratios 5/7 and 3/5 are similar, while between the Sentinel-2A MSI and Landsat 8 OLI bottom-of-atmosphere products, major differences appear in gain and offset.The Sentinel-2A MSI and Landsat 8 OLI top-of-atmosphere products suffer considerably less from this.Lastly, it can be observed that the two scatterplots that involve band 1 (blue) of the atmospherically corrected Sentinel-2A MSI products (center-left and centre-centre) both show a nonlinear trend.A comparison between the top-of-atmosphere and bottom-of-atmosphere products (Supplemental Materials) shows a correlations between 0.85 and 0.95, which does not reveal any abnormal behaviour for this particular band.The top row shows ratio TM 5/7, the centre row shows ratio TM 3/1 and the bottom row shows ratio TM 3/5.It can be observed that the correlation between the top-of-atmosphere reflectance products is higher than the correlation between the bottom-of-atmosphere reflectance products.
Comparing the Sentinel-2A MSI band ratio product with the geological map, in Figure 8, shows that pixels with a relatively high value for TM 5/7 (hydroxyl bearing alteration, in red) coincide with Dacitic and Rhyolithic pyroclasts, domes and outflows (units 8, 10 and 11).It also appears that large parts of these volcanic rock units fall within image pixels that are masked, and are apparently covered by vegetation.The iron oxides, depicted in green and blue, are generally within the areas that have been mapped by Arribas [50] as quaternary deposits (unit 21).

Discussion
Three band ratios, previously defined for Landsat 5 TM, were used for mapping mineralogy associated with a hydrothermal alteration system in southeast Spain.The datasets were a subset of a Landsat 8 OLI and a Sentinel-2A MSI scene acquired in 2016 and a simulated Sentinel-2A MSI scene that was acquired in 2004.These data cover a geologically varied terrain, with marine limestones in the western part of the image, hills and slopes formed by highly alterated volcanic tuffs containing sulfates and iron oxides in the centre part, and alluvial planes that contain eroded products from the surrounding hills.The surface cover of this area is complex.However, the similarity between the image products is, to the human eye, striking.Even when considering that the 12 years between the acquisition dates are not that relevant on a geological time scale, it must be concluded that, besides the vegetation cover in the agricultural fields, not much has changed in this area when it is observed with band ratios for depicting alteration mineralogy.
The similarity between the two multispectral scenes acquired in May 2016 is higher than the similarity with the dataset created from a hyperspectral scene acquired in May 2004.An explanation for this difference points towards land cover, weather and sensor differences.A small change in land cover might not easily be observed within single 10-20 m pixels.While the weather conditions of data acquisition both consisted of dry and sunny weather for several days in a row before data acquisition (Table 1), there may still remain a difference in how far the season (spring) was in either year.
The pre-processing done on the hyperspectral dataset involved a MODTRAN (MODerate resolution atmospheric TRANsmission) model to come to surface reflectance values.The correction applied to the Landsat 8 OLI dataset is based on the S6 atmospheric model, while the correction applied to the Sentinel-2A MSI dataset with the sen2cor software is based on "look-up tables compiled using an atmospheric radiative transfer model based on libRadtran1" [51].The results in Figure 7 reveal that the Sentinel-2A MSI bottom-of-atmosphere band 2 does not compare as well as the top-of-atmosphere product does.Nevertheless, a direct comparison between the bottom-of-atmosphere and top-of-atmosphere product (Supplemental Materials) does not reveal a problem in the atmospheric correction.
Toming et al. [30], in their work on mapping lake water quality parameters with Sentinel-2A MSI, concluded that the currently provided atmospheric correction is in a need of improvement.However, most other studies [24][25][26][27][28][29]31] on initial results with Sentinel-2A MSI did not specifically mention this as an observation or a problem.
Overall, the obtained data products look similar and are likely to lead, when visually interpreted by one and the same end-user, to similar outcomes.This might well change when the products are interpreted or compared based on the values of the data products.Although the correlation between Landsat 8 OLI and Sentinel-2A MSI for this area and datasets was found to be in the range of 0.56-0.97,there are significant differences in the data ranges.In a (semi-)automated processing chain, an additional effort will be needed to make these data products repeatable in time and transportable between sensors.In our study, this became already apparent in the pre-processing of the data, when a masking of vegetation with NDVI had to be carried out manually due to the significant different data ranges that were obtained.This is, however, a common feature of remote sensing, and is not specific for the Sentinel-2A MSI or Landsat 8 OLI imagers in particular.
We have compared Sentinel-2A MSI with Landsat 8 OLI and used band ratios designed for Landsat family broadband multispectral sensors.ASTER band ratios have partially been ported to Sentinel-2A MSI by Van der Meer et al. [16], but these were not explored in detail in this study.However, this does not undermine the importance of Sentinel-2A MSI to also provide data continuity for ASTER.Despite the ASTER sensor going defunct in 2009 and the availability of a commercial high-spatial resolution sensor [52] that has four of ASTER's five SWIR bands, the use of ASTER for geological studies is still widespread.In any way, having a more frequent data acquisition is a welcome benefit.
The five to ten-day minimum revisit time of Sentinel-2A MSI is not only important for environmental studies that require frequent data: geological remote sensing would also benefit from frequent data.An example is the monitoring of geological activity associated with volcanoes, tectonics and geothermal systems, but also secondary effects related to earth sciences, such as land cover changes associated with natural degradation or natural hazards.However, the major leap that data continuity with operational satellite sensors such as Sentinel-2A MSI and Landsat 8 OLI can bring is the generation of physically consistent, reproducible and seamless data products, such as the seamless maps of surface mineralogy published for the continent Australia [53].The future of geological remote sensing will incorporate the generation of such continental scale maps, and data continuity in time and space is therefore essential.

Conclusions
Three band ratios based on ESA's Sentinel-2A MSI were used to map mineralogy associated with a hydrothermal alteration system in southeastSpain.The results were compared with similar products derived from Landsat 8 OLI and a simulated Sentinel-2A MSI dataset used in Van der Meer et al. [16].Although the surface cover of this area is complex in geological terms, a clear similarity between the image products can be observed with the human eye.When the products are, however, compared based on data values, it appears that particular band combinations have a considerable difference in gain and/or offset, despite also having an average correlation of approximately 0.8.Any specific problem with the atmospheric correction of Sentinel-2A MSI or Landsat 8 OLI could, however, not be found.
The resulting image products demonstrate a good correspondence between Sentinel-2A MSI and Landsat 8 OLI VNIR and SWIR bands.A number of band ratios originally proposed for Landsat 5 TM correlate favourably when ported to Sentinel-2A MSI.The resulting products are also compared to a published geologic map of the study area, and it is shown that the resulting maps support the existing conceptual geologic model of the epithermal deposit.The future of geological remote sensing will also incorporate the generation of continental scale maps that somehow have to be consistent in time and space.Sensor cross and inter-calibration as well as pre-processing methods, therefore, remain an important issue.Although the repeatability and portability of the Sentinel-2A MSI and Landsat 8 OLI data products is still a work in progress, this study shows that the Sentinel-2A MSI mission can provide data continuity for Landsat 8 OLI when studying mineralogy at the Earth surface.

Figure 2 .
Figure 2.Location of the former mining town of Rodalquilar in the volcanic belt of southeast Spain (inset (A)) and the location of the image dataset used in this study (inset (B)).The Cabo de Gata volcanic field consists of calc-alkaline volcanic rocks (andesites and rhyolites) of late Tertiary age which have been extensively altered to form an assemblage of metamorphic minerals from high to low temperature as: silica, alunite, kaolinite, montmorillonite and chlorite.The geology, geochemistry and mineralization of the area is documented[35], while several remote sensing studies have been conducted at this site[36][37][38][39].Image modified after Arribas et al.[35] and Rytuba et al.[32] and used with permission from Van der Meer et al.[16].

Figure 3 .
Figure 3. True colour composite of Sentinel-2A MSI bands 4 (in red), 3 (in green) and 2 (in blue).The image covers the 404 × 203 pixels subset of the study area around the mining town of Rodalquilar, Cabo de Gata, Spain.

Figure 4 .
Figure 4. Normalized Difference Vegetation Index (NDVI) image result of the simulated Sentinel-2A MSI data cube dating from 2004 (top), the Sentinel-2A MSI dating from 2016 (centre), and the Landsat 8 OLI (bottom).The images shown are histogram stretched from 2%-98%.The town of Rodalquilar in the centre has relatively low NDVI values, while the hills north and south of it have the highest values.Other high values, such as those associated with the fields in the northeastern part of the simulated image, seem to be related to agricultural activity and not to geology.

Figure 5 .
Figure 5. Scatterplots of the Normalized Difference Vegetation Index (NDVI) made with Landsat 8 OLI, Sentinel-2A MSI and a simulated Sentinel-2A MSI data cube.The figure on the left shows the scatter between the 2016 Sentinel-2A MSI bottom-of atmosphere (S2A-BoA) NDVI product and the 2004 simulated dataset (Sim).Note the relatively high values in the simulated datasets that are associated with the agricultural fields in the northeast of the scene.The centre figure shows the scatter between the 2016 Sentinel-2A MSI and Landsat 8 OLI bottom-of-atmosphere (S2A-BoA and LS8-BoA) products.The figure on the right shows the scatter between the Sentinel-2A MSI and Landsat 8 OLI top-of-atmosphere (S2A-ToA and LS8-ToA) products.

Figure 6 .
Figure 6.Sabins' band ratio colour composites of a simulated Sentinel-2A MSI data cube (top), Sentinel-2A MSI (centre) and Landsat 8 OLI (bottom).Shown are ratios TM 5/7 (ShortWave InfraRed 1/ShortWave InfraRed 2) in red, TM 3/1 (Red/Blue) in green and TM 3/5 (Red/ShortWave InfraRed 1) in blue.The patterns appear similar to the human eye, but differences can be observed in the ratio TM 3/1 between the simulated dataset of 2004 and the Landsat 8 OLI and Sentinel-2A MSI datasets of 2016.The colour ramps show that data ranges differ.

Figure 7 .
Figure7.Scatterplots of Sabins' band ratios made with Landsat 8 OLI, Sentinel-2A MSI and a simulated Sentinel-2A MSI data cube.The top row shows ratio TM 5/7, the centre row shows ratio TM 3/1 and the bottom row shows ratio TM 3/5.It can be observed that the correlation between the top-of-atmosphere reflectance products is higher than the correlation between the bottom-of-atmosphere reflectance products.

AFigure 8 .
Figure 8.A colour composite made with Sentinel-2A MSI (A) consisting of band ratios TM 5/7 in red, 3/1 in green and 3/5 in blue.The patterns show a similarity with the volcanic rock units in the Rodalquilar area in a subset of a published geological map [50] (B).

Table 1 .
[49]her conditions in Almeria, Spain, in the week before the data acquisitions in 2004 and 2016 (indicated in grey).Data source is available online[49].