Multiple Endmember Spectral Mixture Analysis (MESMA) Applied to the Study of Habitat Diversity in the Fine-Grained Landscapes of the Cantabrian Mountains

Heterogeneous and patchy landscapes where vegetation and abiotic factors vary at small spatial scale (fine-grained landscapes) represent a challenge for habitat diversity mapping using remote sensing imagery. In this context, techniques of spectral mixture analysis may have an advantage over traditional methods of land cover classification because they allow to decompose the spectral signature of a mixed pixel into several endmembers and their respective abundances. In this work, we present the application of Multiple Endmember Spectral Mixture Analysis (MESMA) to quantify habitat diversity and assess the compositional turnover at different spatial scales in the fine-grained landscapes of the Cantabrian Mountains (northwestern Iberian Peninsula). A Landsat-8 OLI scene and high-resolution orthophotographs (25 cm) were used to build a region-specific spectral library of the main types of habitats in this region (arboreal vegetation; shrubby vegetation; herbaceous vegetation; rocks-soil and water bodies). We optimized the spectral library with the Iterative Endmember Selection (IES) method and we applied MESMA to unmix the Landsat scene into five fraction images representing the five defined habitats (root mean square error, RMSE <= 0.025 in 99.45% of the pixels). The fraction images were validated by linear regressions using 250 reference plots from the orthophotographs and then used to calculate habitat diversity at the pixel (alpha-diversity: 30 x 30 m), landscape (gamma-diversity: 1 x 1 km) and regional (epsilon-diversity: 110 x 33 km) scales and the compositional turnover (beta- and delta-diversity) according to Simpson's diversity index. Richness and evenness were also computed. Results showed that fraction images were highly related to reference data (R-2 >= 0.73 and RMSE <= 0.18). In general, our findings indicated that habitat diversity was highly dependent on the spatial scale, with values for the Simpson index ranging from 0.20 +/- 0.22 for alpha-diversity to 0.60 +/- 0.09 for gamma-diversity and 0.72 +/- 0.11 for epsilon-diversity. Accordingly, we found beta-diversity to be higher than delta-diversity. This work contributes to advance in the estimation of ecological diversity in complex landscapes, showing the potential of MESMA to quantify habitat diversity in a comprehensive way using Landsat imagery.

Automated summary

Using Multiple Endmember Spectral Mixture Analysis (MESMA), this study quantified habitat diversity in fine-grained landscapes of the Cantabrian Mountains using Landsat imagery. The results demonstrated a significant influence of spatial scale on habitat diversity, with fraction images highly correlated with reference data, showcasing the potential of MESMA in comprehensive habitat diversity quantification with Landsat imagery.