A multiscale analysis of ecosystem services supply in the NW Iberian Peninsula from a functional perspective

In recent years, the assessment of ecosystem services (ES) supply has been based on the use of Land Use/Land Cover (LULC) data as proxies for spatial representation of ecosystems. Nevertheless, some shortcomings of this method, such as uncertainties derived from generalization of the ecosystem types and assumptions of invariance across spatial scales, indicate the need for new approaches. Such approaches could be aimed at improving knowledge of the relationships between ecosystem services and landscape structure and the spatial characteristics of ES patterns. In this study, we propose an integrative approach that involves the generation and analysis of continuous maps representing the supply of five ES potentially related to the amount of biomass. Five remote sensing images of the Northwestern Iberian Peninsula, obtained with Landsat-5 TM, were used to generate a proxy for net primary production by combining the normalized difference vegetation index (NDVI) of each image to calculate a ΣNDVI index that could act as a potential indicator of some ecosystem services. This information was combined with three variables - terrain slope, population density and occurrence of protected areas - to produce spatial models for the five ES and eventually a series of five supply maps. Food, materials and energy provision services showed a clustered pattern, with high supply values in flat zones and areas with high population densities. In contrast, mass flow and climate regulation services were more widely distributed throughout the study area. The five ecosystem service patterns were analyzed at different scales by two methods: lacunarity and four term local quadrat variance (4TLQV) analysis. These methods revealed differences in the spatial pattern: lacunarity analysis was useful for detection of scale thresholds at the local level, whereas 4TLQV was more sensitive to scale thresholds at larger spatial levels. Thus, the variance analysis yielded higher values for larger windows sizes, particularly for provisioning services. The results demonstrated the suitability of the proposed approach for the spatially explicit modeling of ecosystem services, avoiding the uncertainty of other assessments such as those based on LULC data, and for the exploratory analysis of ES supply from a spatial point of view.