Reﬁning the Tiered Approach for Mapping and Assessing Ecosystem Services at the Local Scale: A Case Study in a Rural Landscape in Northern Germany

: Spatially explicit assessments of ecosystem services (ES) potentials are a key component in supporting a sustainable land use management. The ES matrix method is a commonly used approach as it allows for a comparably fast, comprehensible and accessible ES assessment. As it is often based on land use / land cover data (LULC) with no spatial variability, a main critique is that the results fail to assess spatial variability at landscape levels, which limits the reliability of the outputs for spatial planning applications. By using the case study area of Bornhöved in northern Germany, we analyzed three assessment methods that combine expert judgments, LULC data with di ﬀ erent resolutions and ecosystem condition indicators, in order to ﬁnd the required resolution and data for ES assessment and mapping at a local scale. To quantify map discrepancies, we used the structural similarity index (SSIM) and analyzed the di ﬀ erences in local mean, variance and covariance between the maps. We found that using di ﬀ erent spatial resolutions led to a relatively small di ﬀ erence in the outcomes, in which regulation and maintenance services are more a ﬀ ected than the other services categories. For most regulation, maintenance and cultural ES, our results indicate that assessments based only on LULC proxies are not suitable for a local quantitative assessment of ES, as they cannot su ﬃ ciently cover the spatial heterogeneity of ES capacities that arise from di ﬀ erent ecosystem conditions.


Introduction
Ecosystem services (ES) are defined as the benefits that ecosystems provide in support of human well-being [1]. ES assessment and mapping have become popular approaches for sustainable planning and management of natural resources, for instance in EU policies [2,3]. ES approaches can also help to preserve and enhance landscape multifunctionality [4,5]. Different methodologies have been developed for assessing and mapping ES (for a review see [5]), that can be classified in three mdifferent tiers [6]. The appropriate ES assessment methods or tier depend on the purpose of the assessment (e.g., research • How does a higher spatial resolution with more information on small landscape elements affect the results of an ES-matrix assessment on a local scale? • Does the integration of ecosystem condition information add value to the ES assessment and can patterns between different ES and ES categories be detected? • What conclusions can be drawn for practical applications in landscape management?

Case Study Area (CSA)
The CSA is located in the region of the Bornhöved Lake District in the province of Schleswig-Holstein in northern Germany ( Figure 1). It was the focus of a former integrative ecological study project and is a long-term ecological research (LTER) site. The present study area is bigger than the former LTER site and includes the administrative borders of the constitutive municipalities, resulting in a CSA of 147 km 2 . Agroecosystems dominate the landscape in a catchment area of five glacially formed and consecutively connected lakes. The CSA is composed of relatively small and diverse ecosystems, as well as a high proportion of hedgerows, shaping a highly diverse and fragmented landscape (see Figure 1 and [27]). The former LTER site was used in several ES assessments focusing on one single ES [28] or a group of ES [7][8][9]26].

Ecosystem Services Assessments
To guide through the following section, we introduce the workflow of the different levels of assessment developed in this study. Figure 2 illustrates this workflow and how the different datasets are integrated in the different levels of assessment. We prepared two geospatial maps using three different datasets: the CLC datasets for the first assessment level, the ATKIS and InVeKoS datasets were combined for the second and third assessments. We used the LULC classes from the three datasets to build a unique ES matrix based on a common typology and solicited local experts. For the first and second level assessment, we applied the resulting matrix to the different geospatial datasets to map ES potential in the CSA. For the third level assessment, we added ecosystem condition indicators to the analysis: ES potentials were calculated based on the environmental condition scores of each ecosystem and the ES scores from the matrix. Finally, we statistically compared the different levels of assessment and the differences in the resulting ES potential maps.

Ecosystem Services Assessments
To guide through the following section, we introduce the workflow of the different levels of assessment developed in this study. Figure 2 illustrates this workflow and how the different datasets are integrated in the different levels of assessment. We prepared two geospatial maps using three different datasets: the CLC datasets for the first assessment level, the ATKIS and InVeKoS datasets were combined for the second and third assessments. We used the LULC classes from the three datasets to build a unique ES matrix based on a common typology and solicited local experts. For the first and second level assessment, we applied the resulting matrix to the different geospatial datasets to map ES potential in the CSA. For the third level assessment, we added ecosystem condition indicators to the analysis: ES potentials were calculated based on the environmental condition scores of each ecosystem and the ES scores from the matrix. Finally, we statistically compared the different levels of assessment and the differences in the resulting ES potential maps.

Geospatial Datasets
The LULC typology used for the first level assessment is derived from the CLC 2018 dataset (https://land.copernicus.eu/pan-european/corine-land-cover/clc2018).
The geospatial dataset used in the second level assessment is composed of ATKIS (version 2012) and improved by using InVeKoS data (version 2010). ATKIS is a national standardized topographic information system provided by the German Federal Surveying Authorities. InVeKoS is a control system established by the European Commission to control payments of income support from the The LULC typology used for the first level assessment is derived from the CLC 2018 dataset (https://land.copernicus.eu/pan-european/corine-land-cover/clc2018).
The geospatial dataset used in the second level assessment is composed of ATKIS (version 2012) and improved by using InVeKoS data (version 2010). ATKIS is a national standardized topographic information system provided by the German Federal Surveying Authorities. InVeKoS is a control system established by the European Commission to control payments of income support from the European Common Agricultural Policy to farmers. Both datasets feature an object-based comprehensive description of the Earth's surface. Topographical and LULC features are modeled as points, polylines or polygons. The spatial resolution of these two datasets depends on the feature classes and has a MMU between 0.1 and 1 hectare [29] compared to the 25 hectares of CLC.
The two datasets differ in represented LULC classes: "transitional woodland-shrub" and "land principally occupied by agriculture, with significant areas of natural vegetation" from the CLC dataset are not part of the ATKIS/InVeKoS classification systems. The CLC class "transitional woodland/shrub" corresponds to "areas representing natural development of forest formations, consisting of young plants of broad-leaved and coniferous species, with herbaceous vegetation and dispersed solitary adult trees. Transitional process can be for instance natural succession on abandoned agricultural land, regeneration of forest after damages of various origin (e.g., storm, avalanche), stages of forest degeneration caused by natural or anthropogenic stress factors (e.g., drought, pollution), reforestation after clearcutting, afforestation on formerly non-forested natural or semi-natural areas, etc." (https://land.copernicus.eu/user-corner/technical-library/corine-land-cov er-nomenclature-guidelines/html/index-clc-324.html). The CLC class "land principally occupied by agriculture, with significant areas of natural vegetation" is defined as "areas principally occupied by agriculture, interspersed with significant natural or semi-natural areas (including forests, shrubs, wetlands, water bodies, mineral outcrops) in a mosaic pattern." (https://land.copernicus.eu/user-cor ner/technical-library/corine-land-cover-nomenclature-guidelines/html/index-clc-243.html). Table 1 and Figure 3 show all LULC classes and their respective proportions for each LULC dataset are shown in Table 1. The main differences between CLC and ATKIS/InVeKoS are found in the higher number of LULC classes and also in total surface covered by "arable land", by landscape elements (LULC subclass which regroups the classes "hedgerows", "tree rows", "copse, thickets", "single trees", "field borders") as well as by artificial surfaces.

ES Matrix Approach
The ES matrix approach by Burkhard et al. [7][8][9] is based on a look-up table where the capacity of each geospatial unit (e.g., LULC types) to supply ES is quantified. In this study, we assessed the ES potential of the different LULC classes. We did not consider the actual use or flow of ES, nor the demand for ES (see Syrbe et al. [30] for definitions of ES potential, stock, use and demand). ES potentials are scored on a scale from zero (no relevant potential) to five (very high relevant potential). The methodology is also described in Burkhard and Maes [31] and in Campagne and Roche [13]. We adhered to the following steps: (1) identification of the relevant ES and ecosystems to be assessed, (2) selection of an expert panel, (3) collection of expert scorings and (4) compilation of the final matrix values.
The list of ES in the matrix is based on the common international classification of ecosystem services (CICES V5.1) [32] and covers all three main CICES sections, i.e., regulation and maintenance ES, provisioning ES and cultural ES. Each ES is briefly explained in the matrix in the supplement materials ( Table 1). The list of ecosystems corresponds to the CLC, ATKIS and InVeKoS LULC classes in the CSA (see Table 1 "Corresponding ATKIS/InVeKoS classes"). The experts had a wide range of ES and ecosystems expertise, had researched the study area and were familiar with ES concepts.
Müller et al. [33] developed an ES assessment matrix for CLC land cover classes within the German Federal State of Schleswig-Holstein, to which more than 55 experts contributed. We used this matrix for the first ES assessment, as all the selected ES and LULC were scored in this matrix. For the second assessment, we completed the matrix with all the missing LULC, i.e., the LULC classes within the ATKIS/InVeKoS datasets that were not represented in CLC classes. These classes were assessed through a new evaluation process, in which several preliminary discussion rounds were organized with four experts in order to identify in the missing values of the matrix. Other experts

ES Matrix Approach
The ES matrix approach by Burkhard et al. [7][8][9] is based on a look-up table where the capacity of each geospatial unit (e.g., LULC types) to supply ES is quantified. In this study, we assessed the ES potential of the different LULC classes. We did not consider the actual use or flow of ES, nor the demand for ES (see Syrbe et al. [30] for definitions of ES potential, stock, use and demand). ES potentials are scored on a scale from zero (no relevant potential) to five (very high relevant potential). The methodology is also described in Burkhard and Maes [31] and in Campagne and Roche [13]. We adhered to the following steps: (1) identification of the relevant ES and ecosystems to be assessed, (2) selection of an expert panel, (3) collection of expert scorings and (4) compilation of the final matrix values.
The list of ES in the matrix is based on the common international classification of ecosystem services (CICES V5.1) [32] and covers all three main CICES sections, i.e., regulation and maintenance ES, provisioning ES and cultural ES. Each ES is briefly explained in the matrix in the supplement materials ( Table 1). The list of ecosystems corresponds to the CLC, ATKIS and InVeKoS LULC classes in the CSA (see Table 1 "Corresponding ATKIS/InVeKoS classes"). The experts had a wide range of (ATKIS/InVeKoS) of the total area, from which 67.2% (CLC) or 60.6% (ATKIS/InVeKoS) is dedicated to arable land. Forests covers only 9.2% (CLC) or 10.7% (ATKIS/InVeKoS) of the case study area. 6.4 (CLC) to 3.8% (ATKIS/InVeKoS) of the area is covered by artificial surfaces. In the CLC datasets, 0.1% of the area is covered by transitional woodland-shrub which has no equivalent in the other dataset. Landscape elements represent 2.9% of the total area in the ATKIS/InVeKoS datasets.

Selecting Suitable Indicators
The MAES working group of the European Commission (http://biodiversity.europa.eu/maes) aims at providing support to the EU member states in regard to ES mapping and assessment [34]. The fifth MAES report presents a list of key indicators for mapping and assessing ecosystem condition of the main ecosystem types at the European level [35]. From this report, we selected ecosystem condition key indicators for which we could find local data or easily develop an equivalent that was spatially explicit and sensitive to local land use change. Whereas MAES provides one list of indicators per ecosystem type, the indicators selected for this study are not ecosystem type-dependent Land 2020, 9, 348 8 of 23 and are relevant for artificial areas, agroecosystems, forests as well as other natural and semi-natural ecosystems. Three indicators were selected: the landscape fragmentation index, a Red List Index for ecosystems and the soil and organic carbon (SOC) ( Table 2).

Landscape Fragmentation Index
Landscape fragmentation results from the transformation of large natural and semi-natural ecosystems into smaller and isolated, fragmented ecosystems [36]. This process can for instance be caused by the construction of transport networks and built-up areas, which have greatly expanded over the last 50 years. Accessibility is a prerequisite for the exploitation and supply of a wide range of ES. Up to a given density, there is a positive relationship between road density and the potential use of provisional and cultural services, as it generally depends on access to the ES supply areas. However a high fragmentation has a negative impact on agricultural and timber production, as small land parcels tend to have lower yields [36]. Transportation networks also have negative impacts on cultural ES, above all for recreational use and landscape aesthetics. Fragmentation seriously affects regulation ES by, among other, species movement impediment and habitat reduction [37], water cycle disruption [38] and erosion phenomenon exacerbation [38]. A detailed list of the effects of fragmentation on ecosystem conditions and services can be found in the report of the European Environment Agency (EEA) [36].
The EEA applied the method of effective mesh density to measure the degree of landscape fragmentation in the EU member states [36]. This landscape metric unit is based on the probability that two points chosen randomly in an area are connected and not separated by any natural or man-made barriers. The effective mesh size represents the degree of fragmentation of a landscape: the smaller the value, the more fragmented the landscape. The datasets used to calculate the mesh are CLC (2006) for artificial areas such as urban and industrial units but also waterbodies, TeleAtlas Multinet for the roads and railroads, WorldClim and Nordregio for topographic and climatic barriers [36]. The resulting map presents landscape fragmentation at 1 km 2 resolution.

Local Red List Index
The Red List Index (RLI) is one prominent indicator in the Convention for Biological Diversity (CBD) indicator sets for biodiversity conservation [39], but has until now not been used for ES assessments. RLI is a key condition indicator in the fifth MAES report that indicates "conservation status and trends of species of community interest" and "conservation status and trends of habitats of community interest" [35] and that aims at supporting the evaluation of the efficiency of conservation policies. The RLI measures status and trends in the overall extinction risk of species by using weight scores based on the Red List status of each assessed species, which range from 0 (least concern) to 5 (extinct/extinct in the wild). To date, global RLIs have been calculated for birds, mammals, amphibians, corals and cycads [40]. They serve as indicators to evaluate progress towards meeting the CBD "2010 biodiversity target" and help to develop a better understanding of which taxa, regions or ecosystems are declining or improving globally. Global RLIs and Red List assessments are hardly sensitive to local land-use management and decisions, so we instead used local extinction risk assessments lists provided by the Land 2020, 9, 348 9 of 23 State Agency for Nature and Environment of Schleswig-Holstein [41][42][43]). Moreover, as we needed a spatially explicit indicator, we calculated local RLIs for each ecosystem and not for each species group, following the method by Juslén et al. [44]. Each species was associated to one or more ecosystems based on their habitats described in national or Schleswig-Holstein Red List assessments. A total of 214 species of birds, 44 mammals and 18 amphibians and reptiles were included in the assessment.
The RLI values were calculated following the method developed by Butchart et al. [40]. We multiplied the number of taxa in each ecosystem by the category weight (0 for least concern, 1 for near threatened, 2 for vulnerable, 3 for endangered, 4 for critically endangered and 5 for extinct/extinct in the wild). These products are summed up and then divided by the number of taxa multiplied by the maximum weight 5, the "maximum possible denominator", which corresponds to the "worst case scenario". To obtain the RLI value, this sum is subtracted from 1. The resulting index varies between zero and one: the lower the value, the closer the set of taxa living in the ecosystem is heading towards extinction. The Red List status of the pool of species living in an ecosystem is used as a proxy for this ecosystem condition. To make this indicator sensitive to land management change, we included the RLI value trends as well as the values from the last Red List assessment. The "RLI status" corresponds to the values from the last Red List assessment. The "RLI trend" corresponds to the difference between the last and the first Red List assessment values (i.e., RLI status-RLI 1990 (or 1995)) for each ecosystem.

Soil Organic Carbon (SOC)
Soil is a key component of terrestrial ecosystems and plays an essential role in ecosystem functioning [45]. Despite its importance, relatively little emphasis has been placed on soil in studies that map ES, mainly because our understanding on how soils support ES delivery is incomplete [46,47]. Soil quality and its capacity to provide a wide range of ES depend on key properties such as soil type and texture, soil organic carbon content, pH, soil porosity and structure [48][49][50]. Soil organic carbon (SOC) is mainly determined by the balance between net primary production and the rate of decomposition of organic material, which both depend on land use, climate and microorganism action. A literature review showed that SOC can be (directly or indirectly) linked to a broad range of ES, in all three ES CICES categories [48]. As no map was readily available at the desired local scale, we developed a SOC map based on the data and method from the Federal Institute for Geosciences and Natural Resources' (BGR) soil information system (FISBo), which provides soil information at land and national scale [51,52]. The method from the BGR allows to determine the SOC from soil characteristics, climatic data and land cover types.

First and Second Level
Ecosystem mapping was carried out for the CSA using CLC (first assessment level) and ATKIS/InVeKoS (second assessment level) datasets. All spatial analyses were executed using ArcGis 10.6.1. ES potential maps were obtained by joining each LULC map with the ES potential matrix. Exemplary ES potential maps (crop production, pollination and recreation and tourism) are presented in the results, the other maps can be found in the Supplementary Materials.

Third Level
For the third level assessment, we assume that good ecosystem condition is prerequisite for a high supply of ES [35] and therefore first combined the three selected indicators into one ecosystem condition map ( Figure S1). However, most provisioning ES, such as agricultural production (crop, livestock, biomass for energy, fodder and fiber production) but also timber production, often arise much more from human interventions than from ecological properties and functions of the ecological system [53] (but see Balvanera et al. [54] and Harrison et al. [55]). Some ecosystem properties have a direct positive impact on agricultural production, but it can be generally substituted by human actions. For example, a high SOC maximizes crop and pasture biomass. A low SOC can be however compensated by adding off-farm organic residues. As we cannot assume that the potential of agricultural and timber production services increases with ecosystem conditions, we did not include them into the third assessment level.
A second assumption, in the absence of further evidence [19], was that each indicator equally contributes to the ecosystem condition of each ecosystem type except aquatic ecosystems for which the SOC is not relevant, following the Equation (1): where EC is the aggregated indicator value for the ecosystem condition, FI the fragmentation index, RLI status the RLI for Schleswig-Holstein for 2014 and RLI trend the improvement or degradation of the RLI over the years, SOC the soil organic content. All indicators are scaled down to the interval [0,1], to give all the indicators an equivalent weight. The EC values are also scaled into the range [0,1]. EC values from zero (the ecosystem is in a very bad condition) and one (the ecosystem is in a very good condition). As the SOC is not relevant for aquatic ecosystems, we did not include this indicator in the EC of water bodies nor water courses. For these ecosystems, Equation (2) is defined as follows: The ecosystem capacity to deliver ecosystem services is determined by the ecosystem condition, following Equation (3): where ES potential is the calculated ecosystem services potential resulting from the ecosystem condition (EC) and the ES potential scores (ES) from the matrix.

Influence of Typology and Resolution
The CLC and ATKIS/InVeKoS datasets use different typologies to describe some of the LULC classes. The differences in the ES maps resulting from the different LULC datasets are therefore not only due to different resolutions, but also to a difference in typologies. To differentiate these two effects, we compared the two maps using the Union tool from ArcGIS. The resulting map is a map containing each input feature from the two LULC datasets, associated with its attributes. The resulting attribute table shows how the two datasets overlap and this allowed to determine which CLC LULC classes corresponds to each ATKIS/InVeKoS class. Using this attribute table, we then calculated the percentage of each CLC class, which was classified (or not) to the corresponding ATKIS/InVeKoS class (see Table 1). To assess the effect of spatial resolution, we assessed the percentage of LULC classes that were misclassified according to the size of each feature.

Statistical Maps Comparison
We first compared the maps resulting from the first and second assessments to evaluate the improvement in mapping precision. We then compared the maps resulting from the second and third assessments to estimate the mapping accuracy discrepancy. We used the structural similarity index (SSIM index) from Jones et al. [56] for visual interpretation and statistical analysis. This index provides three distinctive measures of spatial heterogeneity in map structures: similarity in means (SIM), similarity in variance (SIV) and similarity in pattern (SIP) of spatial covariance. SIM is the ratio of twice the product of the local means to their summed squares. It ranges from 0 to 1, in which 0 characterizes two maps for which means are dissimilar (in this case, different modelled ES local abundances), 1 if means are similar (similar modelled ES local abundances). The SIV is the ratio of twice the product of the local standard deviations to their summed variances. SIV ranges from 0 to 1, in which 0 is indicated if variances are dissimilar (ES are not similarly distributed), 1 if variances are similar (ES are similarly distributed). The SIP is the ratio of the local covariance to the product of the local standard deviation. It ranges between −1 and 1, in which −1 indicates a negative spatial correlation, 1 a positive spatial correlation and 0 indicates that the distributions are independent. The SSIM index compares output values of maps at pixel level and accounts for the values of neighboring pixels through a window of 3 × 3 pixels. The method allows for a spatially explicit difference analysis between two maps, including spatially explicit differences in means, variances and covariances, taking spatial dependencies between neighboring cells into account. All statistical analyses were conducted using the software R [57] and based on the script provided by Jones et al. [56].

First Level Assessment
For the first level assessment, the ES potential supply was evaluated and mapped using the ES assessment matrix and CLC data. Figure 4 shows the spatial distribution of three ES potentials in the CSA: crop production, pollination and recreation and tourism services. A relative scale ranging from zero (no relevant capacity) to five (very high relevant capacity) was used. As the study area is mainly composed of "arable land" (see Table 1 for a description of land use classes share within the CSA), the maps show an overall high crop production potential, a low pollination potential and a moderate recreation and tourism potential for the area. High and very high pollination potentials are patchily distributed and correspond to "broad-leaved forest", "coniferous forest", "Pastures" and "land principally occupied by agriculture, with significant areas of natural vegetation". High and very high recreation and tourism potentials correspond to "broad-leaved forest", "coniferous forest" and "water bodies".   Figure 5 shows the spatial distribution of three ES in the CSA: crop production, pollination and recreation and tourism services. These maps are based on the ATKIS/InVeKoS datasets and show a more fragmented ES potential than the maps from the first level assessment, due to the inclusion of small habitats and landscape elements. Yet the resulting maps show the same global patterns ( Figure  5): fewer areas sustain crop production, but have an overall high crop production potential, a low pollination potential and a moderate recreation and tourism potential.  Figure 5 shows the spatial distribution of three ES in the CSA: crop production, pollination and recreation and tourism services. These maps are based on the ATKIS/InVeKoS datasets and show a more fragmented ES potential than the maps from the first level assessment, due to the inclusion of small habitats and landscape elements. Yet the resulting maps show the same global patterns Land 2020, 9, 348 12 of 23 ( Figure 5): fewer areas sustain crop production, but have an overall high crop production potential, a low pollination potential and a moderate recreation and tourism potential. 5.0: very high potential. Figure 5 shows the spatial distribution of three ES in the CSA: crop production, pollination and recreation and tourism services. These maps are based on the ATKIS/InVeKoS datasets and show a more fragmented ES potential than the maps from the first level assessment, due to the inclusion of small habitats and landscape elements. Yet the resulting maps show the same global patterns ( Figure  5): fewer areas sustain crop production, but have an overall high crop production potential, a low pollination potential and a moderate recreation and tourism potential.

Third Level Assessment
The final ES potential maps ( Figure 6) were obtained by combining the ecosystem condition index map ( Figure S1 in the Supplementary Materials) and the ES potential maps from the second level of assessment. Only those ES which have a positive relationship with EC were included in this

Third Level Assessment
The final ES potential maps ( Figure 6) were obtained by combining the ecosystem condition index map ( Figure S1 in the Supplementary Materials) and the ES potential maps from the second level of assessment. Only those ES which have a positive relationship with EC were included in this assessment (see Section 2.2.3), so for instance crop production was excluded. As none of the ecosystems in the CSA is in a very good condition (EC values are between 0.052 and 0.514, whereas the maximum achievable value is 1), the overall potentials of the CSA to supply ES are greatly reduced.

CORINE Land Cover and ATKIS/InVeKoS Typology Differences
Overall, there is a relatively high correspondence between the two datasets: around 70% of the CSA area is classified in the same LULC class in the two datasets. The main difference between CLC and ATKIS/InVeKoS datasets comes from the different resolutions for the LULC classes in both datasets, as LULC classes, which are mainly composed of large elements in the ATKIS/InVeKoS

CORINE Land Cover and ATKIS/InVeKoS Typology Differences
Overall, there is a relatively high correspondence between the two datasets: around 70% of the CSA area is classified in the same LULC class in the two datasets. The main difference between CLC and ATKIS/InVeKoS datasets comes from the different resolutions for the LULC classes in both datasets, as LULC classes, which are mainly composed of large elements in the ATKIS/InVeKoS dataset, have a high matching level with CLC LULC classes: this is the case for "water bodies" (99%), "non-irrigated arable land" (81%), "industrial or commercial units" (64%) and "coniferous forest" (63%) classes. The LULC classes "discontinuous urban fabric", "mineral extraction sites", "pastures" and "broad-leaved forest" have however a matching level around 50% with a perfect correspondence for elements bigger than the CLC MMU.
There is a clear typology mismatch between the two datasets for the LULC classes that can be found only in the CLC dataset: The CLC LULC class "land principally occupied by agriculture, with significant areas of natural vegetation" mainly corresponds to arable land (27%) or pastures (42%) in the ATKIS/InVeKoS LULC classification, whereas "transitional woodland/shrub" is mainly classified as arable land (52%) or mixed forest (37%). Unlike the CLC classification, the ATKIS/InVeKoS dataset does not differentiate conventional arable land or pastures from agricultural land with natural areas. The ATKIS/InVeKoS classification considers "transitional woodland/shrub" areas as arable land or mixed forest, but does not consider afforestation/deforestation processes unlike the CLC classification (see Section 2.

for a definition of CLC classes)
ES scores were recalculated for each CLC LULC class according to the percentage of each composing ATKIS/InVeKoS LULC class and their ES matrix scores (for instance, if the CLC LULC Class 1 was composed of 10% ATKIS/InVeKoS Class 1 (ES score 2) and 90% of the ATKIS/InVeKoS Class 2 (ES score 0), then the CLC LULC Class 1 has a ES score of 0.2 in the recalculated matrix). The ES scores from the original and the recalculated matrix were then compared. The LULC class "land principally occupied by agriculture, with significant areas of natural vegetation" showed a low variance between the ES scores from the initial matrix and the recalculated matrix (0.6 point) but showed a high variance for the CLC class "transitional woodlands, sclerophyllous vegetation" (1.6 point). This explains some variations in ES mapping, however only for a relatively small fraction of the CSA (the LULC classes "land principally occupied by agriculture, with significant areas of natural vegetation" and "transitional woodlands, sclerophyllous vegetation" represent less than 1% of the area).
For those LULC classes that were not included in the CLC dataset, the landscape elements correspond in the CLC map to the "non-irrigated arable land" (77%) and "pastures" (17%). This indicates that landscape elements are not considered in the CLC map, as they do not match the "land principally occupied by agriculture, with significant areas of natural vegetation" CLC LULC class. We can therefore assume that this LULC class is not integrated in the CLC map because the MMU is too coarse.

Statistical Comparison of First and Second Level Assessment
The resulting maps have similar ES potentials local abundance, variance and are highly correlated with a high to very high similarity in mean values (73.85 to 99.80% of the CSA has a SIM > 0.5), a high similarity in their variances (83.24 to 99.19% of the CSA has a SIV > 0.5) and a high positive spatial correlation (99.61 to 100% of the CSA has a SIP > 0) ( Table S2 in the Supplementary Materials). Different patterns can be observed for the different ES and ES categories: the higher discrepancy between ES maps are for nutrient regulation, crop ( Figure 7) and fiber productions, water purification and global climate regulation services, whereas minerals and fish provisioning have the lowest discrepancy. The differences between crop and fiber supply maps are mainly due to higher dissimilarity in local abundance, whereas the nutrient regulation maps differ in local abundance and variance. Cultural services tend to have a higher SIM than other services maps. Most of the regulation and maintenance services such as pollination, nutrient regulation and global climate regulation, but also cultural services and wild food production maps have a relatively lower SIV, whereas provisioning between ES maps are for nutrient regulation, crop ( Figure 7) and fiber productions, water purification and global climate regulation services, whereas minerals and fish provisioning have the lowest discrepancy. The differences between crop and fiber supply maps are mainly due to higher dissimilarity in local abundance, whereas the nutrient regulation maps differ in local abundance and variance. Cultural services tend to have a higher SIM than other services maps. Most of the regulation and maintenance services such as pollination, nutrient regulation and global climate regulation, but also cultural services and wild food production maps have a relatively lower SIV, whereas provisioning services are characterized by a high SIV (see Figure 7 for exemplary spatial agreement maps showing crop production services).

Statistical Comparison Second and Third Level Assessment
For most of the ES maps, the two mapping methods show highly different local ES potentials, however the variances and covariances have a low to a very high similarity in means (21.71 to 99.11% of the CSA has a SIM > 0.5), but very high similarity in variance (89.09 to 99.82% of the CSA has a SIV > 0.5) and an almost perfect spatial correlation (99.96 to 100% of the CSA has a SIP > 0) for all the ES (Table S3 in the Supplementary Materials). The correlation between the maps is to be expected as the map of the third level assessment is calculated based on the other map. We observe different patterns for each ES category: cultural services have the lowest SSIM values, followed by regulation services (except nutrient regulation) and provisioning services, which have the highest SSIMs. Overall, regulation and maintenance services tend to have a lower SIV than other ES categories (except for groundwater recharge and air quality regulation).
The spatial agreement between the maps depends on the number of ecosystems that deliver each ES, their sizes, the state of the ecosystem and the proportion of the area that potentially deliver each ES. ES that rely on few and in good condition LULC classes tend to have a higher SIM: fish provisioning has the higher SIM values and is sustained by few LULC classes (water bodies and water courses), which represent a small fraction of the CSA and tend to have a better ecological condition than other LULC classes. Nutrient regulation is supplied by a broader range of LULC types and has a high SIM (Figure 8). This may be due to the fact that this ES is mainly sustained by ecosystems that do not represent a high share of the total CSA (landscape elements and forests).
ES. ES that rely on few and in good condition LULC classes tend to have a higher SIM: fish provisioning has the higher SIM values and is sustained by few LULC classes (water bodies and water courses), which represent a small fraction of the CSA and tend to have a better ecological condition than other LULC classes. Nutrient regulation is supplied by a broader range of LULC types and has a high SIM (Figure 8). This may be due to the fact that this ES is mainly sustained by ecosystems that do not represent a high share of the total CSA (landscape elements and forests).

Discrepancies in ES Maps Based on Different Spatial Resolutions
The results of the statistical map comparisons show that using coarse spatial resolution leads to a relatively small, but still nonnegligible differences. Most of the provisioning services and regulation and maintenance services maps are more affected by the spatial resolution than the cultural services maps. The differences in SSIM can be explained by the overestimation or underestimation of different ecosystems in the CLC dataset and by the capacity of these ecosystems to provide multiple or single

Discrepancies in ES Maps Based on Different Spatial Resolutions
The results of the statistical map comparisons show that using coarse spatial resolution leads to a relatively small, but still nonnegligible differences. Most of the provisioning services and regulation and maintenance services maps are more affected by the spatial resolution than the cultural services maps. The differences in SSIM can be explained by the overestimation or underestimation of different ecosystems in the CLC dataset and by the capacity of these ecosystems to provide multiple or single ES. Arable lands mainly deliver one provisioning ES: crop production. Arable lands are overestimated in the CLC dataset, and increasing the spatial resolution reduces the total arable land area. Hence, the capacity of the CSA to produce crops diminishes in the second assessment. Landscape elements (such as hedgerows and field borders) can provide a wide range of ES, particularly maintenance and regulation services. Landscape elements only appear in high resolution LULC data; therefore, maintenance and regulation services are not appropriately mapped when using CLC as input data. As landscape elements are relatively small, the two mapping approaches using different LULC datasets still have a high similarity in variance. Cultural services are generated by more diverse LULC classes and are therefore relatively less affected by the spatial resolution.
For an ES assessment at the scale of the CSA, using CLC or more detailed LULC data would provide outputs with a low degree of structural differences, particularly for cultural services. There is however no unique appropriate resolution and our conclusions may be valid only for a similar scale of assessment. Kandziora et al. [26] showed for instance that the ATKIS/InVeKoS dataset was required for a reliable ES assessment of provisioning services at the plot scale. Similarly, Rioux et al. [58] showed that a spatial resolution of 5 m or finer is required to adequately map the ES provided by small greenspaces in urban areas. This would be especially true for highly fragmented landscapes, with a high share of linear elements and small habitats, which cannot be captured by coarse spatial resolution data such as CLC [59]. On the other hand, a high level of accuracy can be misleading when working at a larger scale and the results may not be serviceable for or understood by decision makers.

Discrepancies in ES Maps Based on Different Data Complexity
The results of the statistical map comparisons show that including ecosystem condition indicators in ES assessments has a high influence on the resulting maps for most of the ES, mainly because of the dissimilarities in mean values, however with a low variance discrepancy and a high spatial correlation. This suggests that LULC proxies, even with a high spatial resolution, do not sufficiently capture spatial heterogeneity in ES potential: for most of the considered ES potentials, local abundance would change considerably when accounting for ecosystem condition. These results concur with the studies by Eigenbrod et al. [25] and van der Biest et al. [22] who question the appropriateness of using LULC data as ES delivery proxies for different assessment scales and for different mapping purposes (but see Roche and Campagne [13] for contrasting results). LULC proxies are not appropriate for a local quantitative assessment i.e., to assess how much ES can be delivered at a specific location, under the assumption that it substantially depends on ecosystem condition. However, the high similarity in variances implies that the maps have similar spatial ES clusters and therefore that mapping ES hot spots and cold spots with the two different LULC datasets would obtain similar results.

Validation
We compared statistical maps to measure the degree of spatial conformity between comparable models with different spatial resolutions and data complexities. This does not constitute a measure of the overall accuracy of the assessment. Assessing model accuracy generally implies a validation using independent data such as primary data, which is generally feasible only for some provisioning services such as agricultural or timber productions [60], but see Sinclair et al. [61] and Chabert and Sarthou [62]. However, models adapted to a local context, with a higher spatial resolution or taking site-specific characteristics into account, generally provide a more precise assessment of ES processes, and their outputs are more reliable [63]. Moreover, if we assume that the capacity of an ecosystem to deliver ES depends on its condition, the results including EC indicators should be closer to the true distribution of ES potential.

Sources of Uncertainties of the Method
ES assessment at the landscape scale is subject to intrinsic uncertainties, which originate from the complexities of landscape and ecosystem processes. First, land cover classifications can imply uncertainties, mainly because of scale mismatches between landscape structure and the spatial data [21]. Even though the accuracy of spatial data has significantly increase with the development of remote sensing technologies, LULC classifications still contains landscape simplifications (as every map), generalizations or even misclassifications [64]. The use of higher spatial resolution data can help to reduce data uncertainties, so in this sense, the use of ATKIS/InVeKoS datasets is more reliable than the CLC dataset.
Second, uncertainties are inherent to the ES assessment field and our still insufficient knowledge on the relationships between ecosystems and ES supply [17][18][19]. For many ES, it is difficult to allocate ES delivery to specific LULC classes, above all for regulation and cultural services, as they are often provided by a combination of LULC with specific characteristics. Even if the importance of a good ecosystem condition for the delivery of ES is widely recognized [24,55,65], the effects of ecosystem condition parameters on ecosystem functioning and the basic mechanisms are still poorly understood [4,66]. van Oudenhoven et al. [67], Kandziora et al. [68], Balvanera et al. [18] and Harrison et al. [55] provide overviews of relationships between different EC/ES indicators and (some) ES, but they do not reach a consensus. Moreover, according to Braat and Brink [69] and Cardinale et al. [70], there may not be a linear relationship between biodiversity, EC and ES. As there is no clear and quantifiable information, it was not possible to adjust the importance of each indicator for each ES.
Lastly, the method entails uncertainties inherent to expert-based methods. Resulting ES assessments are strongly dependent on the level of knowledge of experts which can be subjective in scoring procedures [21]. For instance, less well-known ecosystems and non-directly used services may be underestimated by experts [14,71]. The number and selection of experts are critical and can reduce the level of uncertainty of the model [13,21]. The small size of the CSA considerably reduced the number of potentially relevant experts, as the task necessitates local knowledge on ecosystems and related ES.

Data Availability
The feasibility of the different levels of assessment generally depends on the data availability for the study area. Among the data used in this study, only CLC and landscape fragmentation data are directly available and usable for all EU member states. Land cover datasets with a high resolution (covering for landscape elements) may or may not be available for other countries or regions. In Germany, this information can however not be accessed freely, due to data protection (InVeKoS) or commercial (ATKIS) reasons. For most of the ecosystem condition indicators recommended by the MES report, the needed data were not available for our CSA. Another issue is the scale for which the different indicators have been developed: most of them are not adapted to the local scale and their spatial resolutions do not allow for a local differentiation, and therefore spatially explicit assessment. Even if the availability of remote-sensing data has increased (e.g., Sentinel-2 based data products), data accessibility and availability is still a major issue for ES assessment. A typical example is data on ecosystem pressures such as the use of pesticides, fertilizers and biomass uptake that would greatly improve the index. The data should be made available and used in models to increase the accuracy of results.

Selecting Adapted Ecosystem Condition Indicators
The three ecosystem condition indicators used for this study were chosen because they were recommended in the fifth MAES report (see Maes et al. [35]). These indicators were developed to support national ecosystem assessments and were selected due to their spatial explicitness and sensitivity to change. Indicators that are suitable at the national scale may not necessarily be suitable for a sub-national or a local assessment: the spatial resolution may not be adaptable, and the indicators may be not sensitive to local land management and/or ecosystem condition changes. From the MAES list of recommended key indicators, we selected only relevant and complementary indicators, as a relatively large number of the indicators proposed in the report can be considered redundant and/or not adaptable to the local scale. The three ecosystem condition indicators selected for this study are still not completely independent and to some extent contain redundant information. For example, as biodiversity is usually positively affected by SOC values and negatively by landscape fragmentation, RLIs probably already include some information on the first two indicators. Despite their redundancy, the chosen indicators provide different information on ecosystem condition and their combination provides a relatively good estimation of the condition of different ecosystems. Another issue is that the different indicators are strongly linked to LULC: the different RLIs were calculated per ecosystem type and the SOC and fragmentation values are also largely a consequence of local land use. The outcomes of the second and third level assessments could have been more contradictory if the ecosystem condition indicators would have been measured differently, e.g., based on field data or if we had selected other ecosystem condition indicators.

Conclusions
Our study addresses some important issues in regard to ES assessments at a local level, where the quality of data is often critical. To our knowledge, this is one of the first studies to assess and map ES at the local scale based on ecosystem condition indicators recommended by the fifth MES report.
Land 2020, 9, 348 18 of 23 In contrast to most of the previously published work, a broad range of ES was considered, and we compared the effects of different assessment methods on individual ES and on ES categories.
One of the main challenges for local ES assessments and mapping is finding and applying indicators, related quantification methods and data for comprehensible, reliable and feasible assessments. A reliable analysis should be as detailed as necessary and possible but is however often hampered by the lack of appropriate data for the quantification of each service. The context in which the assessment is implemented should guide the selection of the ES assessment methods and data. Yet a critical aspect for mapping accuracy is the capacity of the used model to capture spatial heterogeneity of ES potentials at the desired level. This paper compared three assessment methods that combine expert judgments, LULC data and ecosystem condition indicators in order to find the required resolution and data for an assessment in a local-scale case study.
Our results suggest that a high spatial resolution including landscape elements (such as hedgerows and field borders) does not drastically change the output of the spatial assessment at the local scale. However, by introducing LULC classes spatial variability, the use of EC indicators in the mapping can lead to different ES potential spatial patterns. This mainly depends on the number and characteristics of the chosen LULC classes, but also the landscape patterns of the studied area: map discrepancies tend to increase with the variety and size of LULC that are linked to ES in the matrix and the proportion of small elements in the studied area.
The results of this study show that expert-based methods combined with local data on ecosystem condition can be a feasible and efficient method for a local assessment of the capacity of a landscape to provide a wide range of ES. Nevertheless, the ES assessment method needs to be improved by exploring the different relationships between ES and the different EC indicators, but also by adding other essential ecosystem condition indicators and socioecological factors. Knowledge of the causal relationships between ES and ecosystem conditions needs to be better documented in different ecosystem types. This should also increase the accuracy of ES models.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-445X/9/10/348/s1, Table S1: ES assessment matrix with expert scores. Figure S1: Distribution of the ecosystem condition index values within the CSA. Figure S2: Crop (human nutrition) production service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S3: Biomass for energy service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure  S4: Crop (fodder) production service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S5: Livestock production service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S6: Timber service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S7: Fibers service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure  S8: Wood fuel service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S9: Wild food service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S10: Fish and seafood service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S11: Beach wrack and flotsam organic material service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S12: Ornamental service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S13: Drinking water service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S14: Abiotic energy service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S15: Mineral production service potential maps using CLC (left) and ATKIS/InVeKoS (right) land cover data and the ES potential matrix. Figure S16: Groundwater recharge service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S17: Local climate regulation service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S18: Global climate regulation service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S19: Flood protection service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S20: Air quality service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S21: Erosion regulation (wind) service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S22: Erosion regulation (water) service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S23: Nutrient regulation service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S24: Water purification service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S25: Pest and disease control service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S26: Pollination service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S27: Recreation and tourism service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S28: Landscape esthetic and inspiration service potential maps using CLC land cover data and the ES potential matrix(left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S29: Knowledge systems service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S30: Cultural heritage service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S31: Regional identity service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Figure S32: Natural heritage service potential maps using CLC land cover data and the ES potential matrix (left), ATKIS/InVeKoS land cover data and the ES potential matrix (middle), ATKIS/InVeKoS land cover data, the ES potential matrix and ecosystem condition indicators (right). Table S2: Spatial agreement between ES potential maps based on CLC and ATKIS/InVeKoS datasets for three ES. Table S3: Spatial agreement between ES potential maps based on ATKIS/InVeKoS dataset only and ATKIS/InVeKoS combined with ecosystem condition indicators for three ES. Funding: This study was part of IMAGINE project (BiodivERsA) supported for Germany by the BMBF-Federal Ministry for Research and Education.