Multifunctionality at what scale? A landscape multifunctionality assessment for the European Union under conditions of land use change

The provision of multiple ecosystem services (ES) within a landscape is commonly referred to as landscape multifunctionality. Modifying landscapes to increase multifunctionality and reduce trade-offs with concurrent services bears the potential to enhance sustainability in human-dominated landscapes. Assessing landscape multifunctionality is thus crucial for land management and planning, but lack of a clear definition and operationalization of multifunctionality impedes comparisons of different study results. We want to address how elements of the study design affect results of multifunctionality assessments. Furthermore, we want to quantify future multifunctionality in the European Union (EU) and indicate the role of land use change and land use diversity on multifunctionality. We analyzed diverging scenarios depicting land use change in the EU between 2000 and 2040 for their effects on landscape multifunctionality. We tested different multifunctionality indicators at various spatial scales based on the modelling of 12 ES and biodiversity indicators. Particularly the analysis scale determines the interpretation of landscape multifunctionality. Coldspots identified by different indicators are in higher agreement than hotspots. We could not confirm links between land use diversity and landscape multifunctionality. While, at EU scale, multifunctionality slightly increases in future scenarios, agricultural intensification and (peri-)urban growth pose large threats to multifunctional landscapes. The choice of indicator and analysis scale strongly determine possible interpretations of the results. Rather than focusing on the impacts of land use change on multifunctionality, it is recommended to base land use policy on the impacts of location-specific change on ES supply and demands.


Introduction
Landscapes provide diverse functions and services. Human-dominated landscapes, such as agricultural land and urban agglomerations, are often managed to serve a particular purpose. However, besides the intended purpose, multiple other services may be provided. The diversity and abundance of different Electronic supplementary material The online version of this article (doi:10.1007/s10980-016-0459-6) contains supplementary material, which is available to authorized users. landscape functions and services within the same landscape are commonly referred to as landscape multifunctionality. Multifunctional landscapes are often thought to be related to high levels of biodiversity (Otte et al. 2007;Pasari et al. 2013). Land use change may affect multifunctionality. Modifying a landscape to maximize the supply of a particular function often involves the depletion of other functions (Goldstein et al. 2012). These trade-offs are particularly strong between regulating services, such as water purification, and provisioning services, such as crop production (Bennett et al. 2009;Raudsepp-Hearne et al. 2010;Howe et al. 2014). In case studies, it has been shown that appropriate land management decisions are able to counteract these trade-offs and facilitate multifunctional landscapes (Paletto et al. 2012;Bradford et al. 2014;Schindler et al. 2014). Modifying landscapes to increase their multifunctionality, therefore, bears the potential to enhance sustainability in human-dominated landscapes (Selman 2009;Waldhardt et al. 2010). Increasing multifunctionality is a proclaimed goal of landscape planning in many EU member states (Helming et al. 2008;van Zanten et al. 2014;Galler 2015). Assessment tools to quantify ecosystem services (ES) and to assess multifunctionality have become popular to determine current and future multifunctionality (Maes et al. 2012;Rodríguez-Loinaz et al. 2015).
The concept of landscape multifunctionality itself is not new (e.g., Beese 1996;Altieri 2000). Conceptual frameworks have been developed in time, and assessments at various scales have been published (e.g., Maes et al. 2012;Paletto et al. 2012;Plieninger et al. 2013b). Most recently, the analysis of time series and the effect of land use change on multifunctionality gained increasing attention (Elmhagen et al. 2015;Rodríguez-Loinaz et al. 2015).
Despite the concept of multifunctionality being elaborated over the last decade, there is still no clear definition as to what the key functions and services are that a particular landscape has to provide to be considered 'multifunctional'. Across the literature, there is no one single way to quantify multifunctionality. Often, the characterization and quantification of multifunctionality is based on metrics, such as diversity indicators that originate from biodiversity research and landscape ecology (Crouzat et al. 2015;Queiroz et al. 2015). The lack of a clear definition and operationalization of multifunctionality in ES assessments leads to problems in interpreting and comparing different studies.
The quantification of multifunctionality is dependent on the spatial scale of analysis. Multifunctionality can be found at the farm-, landscape or regional level (McGranahan 2014). Based on the research question, the chosen grain size of analysis influences the detected type of multifunctionality. For example, analyses on the grid-level (e.g. 1 km) can elucidate the local multifunctionality, while regional analyses may illustrate benefits retrieved from interacting ecosystems (Mastrangelo et al. 2014). Yet, at the regional level, multifunctionality may as well be achieved by diverse monofunctional landscape units together. The scale-effects of multifunctionality are poorly covered in the literature (Mastrangelo et al. 2014). Many assessments are performed solely on arbitrarily chosen scales, such as municipalities or Nomenclature of Territorial Units for Statistics (NUTS) regions (Kienast et al. 2009;Queiroz et al. 2015), or at the level of ecosystem or spatial units (e.g., grid level, Burkhard et al. 2009). Assessments range from local assessments (often relying on field data)  to continental assessments (primarily based on modeled data, e.g. Maes et al. 2014).
Most ES, and thus also multifunctionality, are susceptible to changes in land use and land use intensity (Dobson et al. 2006;Kremen et al. 2007;Metzger et al. 2008;Allan et al. 2015). In Europe, changes such as (peri-)urban growth, afforestation, abandonment of agricultural land, and changes in the land use intensity of agriculture profoundly changed European landscapes over the past decades (Fuchs et al. 2015;Gingrich et al. 2015). Particularly polarization of (rural) land use that consists of opposing trends of land abandonment at marginal locations in combination with intensification of land use at more suitable sites, potentially exerts strong impacts on multifunctional landscapes (Plieninger et al. 2013a).
We hypothesize that land use changes across Europe will affect multifunctionality in different ways and that the ongoing polarization of land use is leading to a decrease in multifunctionality, particularly at local levels.
The objective of this study is to analyze spatial and temporal patterns of multifunctionality in the European Union (EU) and to assess the impacts of land use change on multifunctionality. In analyzing the patterns of multifunctionality we test and compare various ways to quantify multifunctionality and analyze the effect of spatial scale on multifunctionality indicators. To analyze the effect of future land use change on multifunctionality we use a set of ES at the extent of the EU, available for the year 2000, and two diverging future scenarios for the year 2040. We discuss the findings in terms of the methodology and the ambitions of maintaining multi-functional land use.

Quantification of multifunctionality
In the literature, many approaches to quantify multifunctionality are presented. Each of the presented indicators focusses on particular characteristics of multifunctionality. To quantify multifunctionality in our study, we selected four multifunctionality indicators that were recently presented and applied in the literature and that focus on different aspects of multifunctionality (Table 1): Number of ES that exceed particular thresholds (T 30 -T 70 ), the summed ES supply MESLI, Shannon's H 0 (SH), and Simpson's reciprocal index (SRI) (Raudsepp-Hearne et al. 2010;Pasari et al. 2013;Plieninger et al. 2013b;Rodríguez-Loinaz et al. 2015). In all multifunctionality calculations, we applied the same weight to all ES and biodiversity indicators. In the following section, the various approaches to quantify multifunctionality with a single indicator are presented in more detail.
A basic approach, used for example in field studies that rely on survey data, characterizes multifunctionality by the presence or absence of multiple services (Plieninger et al. 2013b), giving an indication of the service richness at a site. However, there is wide agreement that individual services are only beneficial if their supply meets a certain critical level (Raudsepp-Hearne et al. 2010). Therefore, other approaches incorporate this condition and only count those services that exceed this critical level as a contribution to multifunctionality. For many services, such as water quality or carbon sequestration, critical levels can be derived from regulation standards or biogeophysical thresholds.
For the majority of services, however, no such universal thresholds are known or easy to incorporate. To circumvent this issue and to assure that a certain level of service supply is reached, it is possible to derive thresholds from cumulative distribution functions or quantiles of ES supply or biodiversity in multifunctionality assessments (Pasari et al. 2013;Allan et al. 2015). Here, only ES that exceed a threshold set at a particular qquantile (e.g., 0.3, 0.4, …, 0.7) at a given site are included to quantify the multifunctionality of a site. In determining multifunctionality, the number of ES that exceed the chosen threshold are weighted equally, and variations in the level of ES supply are not accounted for.
An approach that aims to integrate the quantity of service supply is employed by Rodríguez-Loinaz et al. (2015). Their multiple ecosystem services landscape index (MESLI) sums ES scores normalized to lower and upper supply benchmarks that can conceptually be represented by thresholds and target levels of ES. In the absence of clear-defined thresholds, however, these levels are set to minimum and maximum observed values following where ES represents ES i supply at location j, and ES i min and ES i max the minimum and maximum supply of ES i . The MESLI is very similar to the total ecosystem service value (TESV) applied by Maes et al. (2012). Both of these approaches, while valuable in discerning coldspots and hotspots of multiple ES supply, may bear the weakness that the resulting multifunctionality indicator may be skewed towards locations that feature few services in abundance, while disregarding the service richness that is commonly associated with multifunctionality.
Diversity indicators that aim to include both the service richness and their supply were originally introduced as a metric for biodiversity (Pielou 1966;Tramer 1969;Hill 1973). One of the most widely employed indicators is the Shannon's H 0 . Shannon's H' quantifies multifunctionality by taking both the abundance of ES and their provision at a given location into account (Eq. 2).
where H 0 is Shannon's H 0 , and p ij represents the relative supply of ES i at location j relative to the total ES supply at location j. The natural logarithm in the equation effects a relative reduction of the weight of very highly supplied ES compared to low supplied ES.
Another diversity indicator is the SRI. The SRI (Eq. 3) is a diversity indicator commonly used in biodiversity research but with some applications in ES research (Raudsepp-Hearne et al. 2010). Simpson's reciprocal index reflects the evenness of ES supply and gives more weight to service indicators that are highly supplied. It relates the relative supply of a particular service at a site to the total number of services.
where SRI is Simpson's reciprocal index, n is the quantity of ES i at location j, and N is the number of ES considered.

Ecosystem service models
We used a set of indicators to quantify twelve ES and biodiversity indicators that are contributing to the multifunctionality of European landscapes ( Table 2). The indicators represent provisioning (n = 2), regulating (n = 4), cultural (n = 3) and supporting ES and biodiversity (n = 3) at EU extent (excluding Croatia). All indicators are mapped at 1 km 2 for the year 2000 and the year 2040 for two diverging future scenarios of simulated land use change. The ES indicators we used for this study are calculated from a range of biogeophysical, geographic and socioeconomic data, including information on land cover, and land management intensity. The models employed are based on simplified representations of the ecosystem processes underlying ES and biodiversity. In that sense, these models go beyond simply mapping ES based on the occurrence of land use types as commonly applied (Burkhard et al. 2009). When ES would be simply related to land use types, multifunctionality would just be dependent on the number of services provided by a particular land use type. The capacity of a land use type to provide services very much depends on the geophysical and socio-economic context as well as on the landscape structure, that is accounted for in many of the ES models employed. In the following, ES and biodiversity will be referred to as ecosystem services. Detailed descriptions of the modelling of ES indicators can be found in the Supplementary Material. Multifunctionality assessments require normalization of ES maps to a common scale. Not all indicators are normally distributed, and some feature infrequent extreme values. To exclude extreme values from the analysis, we transformed the data by means of a 90% winsorization that set all values below the 5th percentile and above the 95th percentile to these thresholds. To standardize between the reference year and the two future scenarios, for each ES, the thresholds used for the winsorization were computed accounting for all values from the three maps. To prevent underestimating these thresholds, only grid cells that can provide a particular ES were included in the calculations. Then, each ES map was normalized from 0 to 1. To be able to assess multifunctionality at different spatial scales, we not only considered the original grid level of the ES maps, but also calculated for each ES map focal means of the normalized maps for six differently-sized square moving windows with variable radiuses (1, 2, 5, 10, 25 and 50 km); for disambiguation, the analysis-scale on the original importance in planning processes and regional policy application within the EU. River catchments were also included to reflect an ecological rather than administrative unit. Computing the different ES indicators for various spatial scales and at several time steps resulted in a total of 360 ES maps used in this study.

Land cover and land use intensity
Land cover in the excluding Croatia (EU) in 2000 was based on CORINE land cover (CLC), spatially and thematically aggregated to 1 km 2 and 16 categories (Verburg and Overmars 2009). Land use and land management indicators are dynamically simulated for the 2000-2040 period based on a series of models at a high spatial resolution of 1 km 2 for two diverging scenarios (Stürck et al. 2015). These models account for global macro-economic developments, European policy, land use history, and local geophysical and socio-economic conditions. The modeling chain is described in more detail in Verburg et al. (2013) and Stürck et al. (2015). The land use and land management scenarios were input to the ES models and the multifunctionality indicators.
We simulated two diverging scenario storylines that resemble those employed by the IPCC Special Report on Energy Scenarios (SRES) A1 and B2 (Nakicenovic and Swart 2000), but were elaborated to represent the EU context. The A1 storyline involves strong economic growth in a globalizing world with strong population growth, no incentives for climate change adaptation and weak regulation of land use change. The Common Agriculture Policy (CAP) is abolished in this scenario, leading to a discontinuation in EU subsidies and programs that strongly effect agricultural land use change in this scenario (e.g. through land abandonment in less suitable or marginal lands).
The B2 storyline represents a fragmented world with regional markets and modest economic growth, modest population growth, and regional implementation of environmental objectives. Under B2 conditions, the CAP is maintained.

Analyses of multifunctionality
We calculated multifunctionality indicators based on ES and biodiversity indicators to answer the following questions: (1) How much does the choice of multifunctionality indicator and scale of analysis influence multifunctionality results? (2) Does land use diversity significantly affect multifunctionality at a European  (2016) All indicators are available at EU extent at 1 km 2 resolution for the year 2000 and the future scenarios V-A1 and V-B2 Landscape Ecol scale? (3) How does multifunctionality in the EU change as result of the land use change scenarios? Each of the four multifunctionality indicators presented in ''Quantification of multifunctionality'' section is used to quantify multifunctionality in the reference year 2000 and for the two scenarios in 2040. Considering that four different multifunctionality indicators are included and the sensitivity of five functional thresholds is tested, this resulted in 8 multifunctionality maps per year and scale of analysis ( Fig. 1).

Comparing multifunctionality indicators
To gain insights into ways in which the choice of multifunctionality indicators influences results, we first retrieved Pearson's r between multifunctionality indicators performed for a random sample of 1000 points across the EU. For analysis at the scale of administrative units or watersheds, we used the entire dataset to calculate correlations.
Secondly, we compared EU-wide aggregated averages of each indicator at each scale of analysis with the map comparison statistic (MCS, Schulp et al. 2014a) to evaluate the spatial correspondence of the different indicators. The MCS measures the correspondence between any two multifunctionality maps, where an MCS of 0 indicates perfect agreement and an MCS of 1 indicates complete disagreement between two maps. To calculate the MCS for each pairing of multifunctionality maps, we normalized the multifunctionality indicators to a range between 0 and 1.
Finally, we compared the maps in terms of the identified hotspots and coldspots of multifunctionality across the EU. To do so, we classified the top and bottom 20% of grid cells (or zones) as indicated by each multifunctionality indicator as multifunctionality hotspot or coldspot. At each analysis scale, we overlaid these hot-and coldspots to find locations of high and low agreement between the indicators used in this study.

The relationship between land use diversity and multifunctionality
To assess the relationship between land use diversity and multifunctionality, we used Shannon's equitability (E H ) to characterize the diversity of land use types in the landscape (Eq. 4).
Where E H is Shannon's equitability, p is the % cover of land use class i within a fixed radius around location j or within a zone j, and N the number of land use classes. This indicator ranges from 0 to 1, where 1 represents complete evenness in the distribution of land use classes. We calculated E H for different analysis scales, including NUTS2 and NUTS3 regions, river catchments, and differently sized square moving windows analogous to the multifunctionality analysis scales. Based on a random sample of 1000 points, we calculated Pearson's r for correlations between multifunctionality indicators and land use diversity (E H ). Again, we used the entire dataset for analysis scales covering administrative units and river catchments. In a second step, we nuanced the analysis by looking at the same relationship separately for different landscape types. To do so, we selected, at each scale of analysis, landscapes that contain, respectively, [50% agriculture, [50% forests, [50% built-up, and[50% (semi-)natural land. We subset the landscape types to two groups. The first group displayed an E H greater than the 80% quantile, that we use here as a cut-of value to address land use diversity hotspots. The second group had an E H lower than 20% quantile that we used to define land use diversity coldspots. Whether multifunctionality in these two groups was significantly different (p value \ 0.05) from each other was assessed with a t test for each landscape type.
Land use change scenarios (2000 -2040) Ecosystem service models

Multifunctionality and land use change
To assess the effect of land use change on multifunctionality, we calculated EU-wide aggregated change according to each multifunctionality indicator at each analysis scale for both scenarios of land use change. Then, we visualized gains and losses in multifunctionality in the EU based on the relative change of multifunctionality between the years 2000 and 2040. Changes in multifunctionality were assessed at locations where different land use change trajectories take place. An overview of land use change trajectories (such as specific land conversions or land use intensity changes) as described in Stürck et al. (2015) is provided in the Supplementary Material.

The role of multifunctionality indicators
We quantified the effect of the scale of analysis on EUscale average values of multifunctionality (Table 3;  Table S2). Figure 2 contains multifunctionality maps for a selection of multifunctionality indicators at various spatial scales. The various indicators respond differently to the scale of analysis: where T 30 and MESLI are least sensitive to the scale of analysis, SRI and Shannon's H 0 increase with increasing scale of analysis (Fig. 2). Multifunctionality decreases for T 50 , T 60 and T 70 with increasing analysis scale.
At the finest resolution, the majority of multifunctionality indicators suggest high multifunctionality in Scandinavia, Spain, Eastern France, South-Western Germany, and Romania. Consistent low levels are found in England, Northern Italy, and Northern France. At larger scales, dissimilarities between the maps increase. For example, at NUTS2 level, with SRI and T 30 , maxima are found in central Europe, while at T 60 , Scandinavia and southern European regions stand out. In contrast, at NUTS2 level, Shannon's H 0 and MESLI both level out intraregional differences that can be found at smaller analysis scales, resulting in a more homogeneous pattern of multifunctionality across Europe.
To assess the overall agreement between the different multifunctionality indicators, we calculated Pearson's r for all pairings of multifunctionality indicators for each time step at each scale (Table 4). All results are provided in Table S3.
Shannon's H 0 and SRI showed correlations higher than 0.93 consistently throughout scales and time steps. The agreement of Shannon's H 0 with other indicators is significantly lower and agreement decreases with analysis scale. Negative correlations, as seen for example at NUTS2 level, indicate contrasting multifunctionality estimations between two indicators. Across all scales, the MESLI showed most frequently the highest agreement with the other indicators.
For all pairings of multifunctionality maps in each year at each analysis scale, we calculated the MCS. The MCS increased with increasing analysis scale, indicating that agreement in the spatial pattern of multifunctionality indicators decreased with increasing analysis scale (Table 5). All results are provided in Table S4.
Finally, we compared hotspots and coldspots (top and bottom 20% of grid cells) as identified by the various multifunctionality indicators (Fig. 3). Hotspots of multifunctionality (Fig. 3a-c) are generally in lower agreement between indicators across all scales than coldspots (Fig. 3d-f). The agreement of hotspot locations was particularly high at smaller analysis scales (radius B 10 km) and between similar indicators. Correspondence in multifunctionality hotspots declines with increasing scale of analysis.
In general, hotspots are associated with mountainous, forested lands on the grid level, most markedly in Germany, Austria, Romania and Scandinavia. Coldspots of multifunctionality are frequently found in landscapes dominated by agriculture, most markedly in Italy, Northern France and England. Also other homogenous landscapes, for example those dominated by pastures and (semi-)natural grasslands in Ireland and Scotland, show a high agreement of coldspots across all scales (Fig. 3d-f).

The relationship between land use diversity and multifunctionality
We quantified the relationship between land use diversity and multifunctionality. Correlations between a land use diversity indicator E H and multifunctionality indicators at various scales suggested only a weak relationship between land use diversity and multifunctionality (Table 6) with correlation coefficients Landscape Ecol smaller than 0.1. With increasing analysis scale, even weak negative correlations are common, indicating negative relationships between land use diversity and multifunctionality. Figure 4 displays a transect in the United Kingdom, in which scale-dependent land use diversity and multifunctionality according to several indicators are presented. Clearly, the different indicators feature partly opposing relationships with land use diversity. Only the MESLI seems to indicate lower multifunctionality where land use diversity drops along the transect.
We also tested the effect of land use diversity on multifunctionality for different landscape types based on the landscape composition, i.e. agriculture-, forest-, nature-, and urban-dominated landscapes (C50% of coverage within the analysis window). Results for agriculture-and forest-dominated landscapes are indicated in Fig. 5a and b, respectively.
The t test results indicate that in most cases, landscapes that differ in their land use diversity also differ significantly in their multifunctionality, though the relationship is less clear on larger analysis scales. We found that the relationship between land use diversity and multifunctionality differs for different landscape types. For example, homogenous agricultural landscapes (i.e., displaying a low E H ) yield significantly lower multifunctionality than more mosaic agricultural landscapes (Fig. 5a). While this effect is consistent across scales, it does not apply to the functional thresholds (T 30 to T 70 ). Multifunctionality in forest-dominated landscapes shows a different relationship to land use diversity, where, instead, land use diversity coldspots are associated with higher multifunctionality scores for the MESLI and T 40 (Fig. 5b).
For the tested analysis scales, the average multifunctionality indicator in land use diversity hotspots and coldspots differed significantly from the average multifunctionality for forest-dominated landscapes. However, no consistent relationship between land use diversity and multifunctionality can be observed. All results are detailed in Table S5.

The impact of land use change for future landscape multifunctionality
We quantified the effect of the land change scenarios on multifunctionality. Figure 6 exemplarily shows changes in multifunctionality between the year 2000 and 2040 for the V-B2 scenario. Results for the V-A1 scenario can be found in Figure S1. Land use change in both future scenarios suggests slight increases in the overall multifunctionality according to indicators defined by a threshold above the 40%-quantile (T 40 -T 70 ) and the MESLI. The increase in multifunctionality is consistent throughout scales of analysis. Shannon's H 0 and the SRI, however, slightly decreased in time. The overall change in multifunctionality is smaller when the analysis is made at a broader spatial scale (Fig. 6). This indicates that at the regional scale, local losses in multifunctionality are compensated for at intraregional level. For the V-B2 scenario, two consistent, regional trends can be distinguished. Agricultural dominated landscapes in Northern France and Southern UK, as well as regions with strong urban growth in the Netherlands, Belgium and Germany, consistently lose multifunctionality according to all tested indicators and across all scales of analysis. Gains in multifunctionality emerge in the northern UK, Scandinavia, Northern Spain, Italy, and Greece. This development is linked to reforestation and deintensification trends in agriculture.
The impact of individual land use changes on multifunctionality was assessed by calculating the   (Table 7). Results for the V-A1 scenario can be found in Table S6. We visualized the impact on individual ES using radar diagrams for each of the typical land change trajectories (Fig. 7) and assessed the significance of change between the start and end year using t tests (p value \ 0.05). Land changes that negatively affect multifunctionality in both scenarios are peri-urban and urban growth, and the intensification of pastures. Particularly locations in the UK, Belgium and the Netherlands are affected by these trajectories. In Poland, for example, losses of multifunctionality can be associated to cropland intensification. Positive impacts of land change on multifunctionality can be Landscape Ecol associated to pasture de-intensification, cropland deintensification, and conversion of cropland to pasture (Fig. 7). Other land change impacts vary more depending on the ES considered. For example, recultivation of green space increases the supply of agriculture-related services, but entails losses in regulating services and biodiversity. Expansion of large-scale nature leads to a broader supply of regulating and cultural services, but results in the loss of productive functions. These changes are reflected ambiguously by the different multifunctionality indicators (Table 7). Land use change, particularly broad conversions to another land system, result in large trade-offs between ES. Land abandonment entails losses in agricultural dependent services, but increases nearly all other ES. However, also the landscape in which land abandonment occurs is important. In the proximity of large-scale nature, MESLI increases slightly, while otherwise, land abandonment results in slight multifunctionality losses. De-intensification in agricultural systems is associated with subtle gains in multifunctionality without drastic losses in any ES provided. Radar diagrams depicting changes in ES supply by land Landscape Ecol change trajectory for the V-A1 scenario can be found in Figure S2.

Limitations and possibilities of quantifying multifunctionality
In the literature, there is no consistency in the methods used to quantify multifunctionality of land use. Therefore, we tested various approaches to quantify multifunctionality developed or employed in recent studies. Each of the indicators highlights specific aspects of multifunctionality. Quantile-based thresholds help to find locations that provide a vast number of ES that are supplied above a certain value. In this form, it is a simple, straight-forward tool that translates into a comparable multifunctionality indicator. However, the higher the chosen threshold, fewer locations provide multiple ES at that quantity (Table 2). High values of such an indicator tend to point towards locations that provide few services abundantly, while the status of the other services is entirely disregarded. In contrast, when the threshold is set at a low level, the severe depletion of one service could be easily offset by comparatively small gains of another service (Table 7), which hampers the informative value of this indicator.
The MESLI indicator, in turn, is strongly affected by the level of supply of individual services. Because it takes into account the entire distribution of ES in one site, it is an apt tool to gain an overview of the overall level of service provision. It can accurately depict the level of change in multifunctionality, for example as a result of land use change. However, due to the lack of minimum thresholds for many ES, it is necessary to investigate the individual services to ascertain whether high multifunctionality values are a result of few highly supplied services, or if the services are more evenly supplied. When interpreting change, this indicator does not differentiate between losses in an ES that has a high supply or a low supply (leading to a decrease in service richness).
Shannon's H 0 limits the weight of highly supplied ES within the multifunctionality analysis, and thus, accentuates the presence of rare, or lowly supplied, ES that support the overall ES diversity. Our results indicate that Shannon's H 0 has particularly low values in homogenous landscapes. As such, it is a highly scale-dependent indicator that requires careful interpretation.
Simpson's Reciprocal Index accentuates the occurrence of very highly supplied ES within a landscape. As it takes into account both the ES supply and the evenness of the supply across ES, it has shown to be particularly pronounced at agricultural fringes (Fig. 8).
While each of the indicators has its virtue, they also entail a range of limitations, and the results of our assessment show that it is of utmost importance to bear the multifunctionality indicator's characteristics in mind while interpreting the results. Relying on a single indicator may lead to results that only address a single element of multifunctionality. Particularly when analyzing change in multifunctionality, it is shown that  It is questionable if there is one singular indicator used in recent literature that is apt to solely base land management decisions upon, or that could replace the thorough investigation of the individual ES that a location comprises. As pointed out by Hansen and Pauleit (2014), solely aiming to improve a landscape's multifunctionality score bears the risk of disregarding the depletion of single ES that have little effect on multifunctionality, but that are fundamental to ecosystem well-being. In particular, in the case of tele-connected effects of ES (Liu et al. 2013), the contribution of individual ES to multifunctionality is convoluted and almost impossible to discern with approaches that are currently employed. We found that the scale of analysis significantly impacts the spatial configuration and overall magnitude of multifunctionality, yet each multifunctionality indicator showed an individual response (Table 3). T 30 as a quantile-based indicator with a low threshold ''profits'' from larger analysis scales because multiple different ecosystems in the landscape may synchronously provide different services, while at higher spatial resolutions, ES are separated in different units of analysis. The higher the chosen threshold, however, the smaller the likelihood for distinct ecosystems to feature a wide range of ES that also surpass the chosen threshold. For Shannon's H 0 and SRI, however, thresholds do not apply, and thus both, on average, increase with increasing scale of analysis. The MESLI, focusing on the sum of normalized ES supply, is fairly consistent over various spatial scales.
Land use change in the EU and its role for multifunctional landscapes In this study we have chosen for a wide range of ES and biodiversity indicators to contribute to multifunctionality because in many landscapes, both ES supply and biodiversity conservation are important objectives that do not necessarily coincide (Bullock et al. 2011;Mace et al. 2012). While ecosystem services are central to human well-being, biodiversity is not only inherently valuable, but also affects the underlying ecosystem service functioning (Cardinale et al. 2012). ES and biodiversity both represent the diversity of functions in the landscape that we tried to integrate in our approach and therefore, we have chosen to base multifunctionality on both ES and biodiversity indicators. We simulated land use change based on two scenarios that represent different political, socioeconomic and technological pathways. These pathways influence the frequency of particular land change trajectories and their spatial occurrence across the EU. In the V-A1 scenario, the abolition of the Common Agricultural Policy causes high ratios of land abandonment, particularly in Southern Europe. This, of Landscape Ecol course, entails drastic effects on agricultural production, but also farm-related services and agro-biodiversity (Rey Benayas 2007), and is, in our results, predominantly associated to multifunctionality loss across the EU (Table 7). However, as shown in Fig. 7a, land abandonment can also be considered as a chance by increasing regulating and cultural services as well as habitat functions and biodiversity (megafauna habitat) in context of the wider landscape configuration (Cerqueira et al. 2015). These locationdependent results are indicative for the locationspecific impacts of particular land use changes on multifunctionality across the EU territory. We also found that locations of urban growth and peri-urban growth are often characterized by comparatively low multifunctionality scores. These trajectories are bound to the vicinity of larger urban agglomerations, and therefore, also reflect the currently low level of multifunctionality in the neighborhood of agglomerations and the rural-urban fringe (e.g., Fig. 7b).
Overall, the results of the scenarios show a spatially variable palette of multifunctionality responses. The small aggregated changes at EU scale do not represent the much larger changes within specific regions and localities that are a result of the many land use changes simulated under the scenarios. Assessments of location-specific impacts on the various dimensions of multifunctionality therefore remain important.

Methodological limitations of the study
This study has relied on the application of a wide range of existing modelling tools to assess ES and land use change at the EU scale. Although these methods have been applied before, they all contain an inherent range of uncertainty, often described in the source publications of the individual methods. Particularly the EUwide methods for ES quantification show large variations when compared to other models (Schulp et al. 2014a).  Multifunctionality results should, therefore, be interpreted as being indicative while absolute values may be subject to uncertainty. Another source of uncertainty derives from the selection of indicators we chose. Not all ES can be represented in the analysis due to a lack of appropriate quantification methods and/or a limited relevance of some ES. We chose a set of indicators that covers a wide range of ES types, such as provisioning, regulating, cultural and supporting services, and biodiversity, to address multifunctionality in a broad sense and to avoid the resulting multifunctionality indicators being biased towards a particular service type. Yet, the choice of ES included in the analysis largely impacts the resulting multifunctionality, and changes or additions to the set of indicators will definitely alter the results.
Particularly in assessing land use change impacts on multifunctionality, the choice of ES will affect the result as land use change often results in trade-offs between different ES types, such as provisional and regulating ES.
Furthermore, it is important to understand what the respective ES models are able to quantify, i.e. potential supply, or use (ES flows). In our case, only few models (e.g., pollination) capture ES flows, while others do not take the use or demand side into account.
Without having knowledge about the demand side, however, and without explicitly taking spatio-temporal trade-offs between ES into account when modelling ES supply (or use), multifunctionality maps as prepared here cannot be used to indicate the sustainability of a given land use.

Implications for land management and further research
In many policy documents, the achievement of multifunctional land use is mentioned without specification of which aspects of multifunctionality, such as service richness or diversity across ES types, are targeted. Implementation of land management and management of land change require a more specific understanding of the scales at which multifunctionality should be achieved and what aspects of multifunctionality are given priority (Raudsepp-Hearne and Peterson 2016). Such will allow an evaluation of alternatives, ex-ante assessment of land use options as well as monitoring.
Our results show that high levels of land use diversity do not necessarily correspond to high levels of multifunctionality. Different land use types do supply different ES, but, particularly when the level of  (1) 3.87 (-2) 3.39 (-5) 1.52 (-21) 3.93 (-11) 1.70 (6) 5.07 (6) Landscape Ecol supply is accounted for, diversity does not straightforwardly translate into high levels of multifunctionality. The local context and the composition and configuration of the landscape are important determinants of its functioning that cannot be captured in simplified indicators. When re-designing landscape management, rather than simply targeting high land use diversity or multifunctionality, the functioning of the ecosystems in relation to specific service demands should be accounted for. Similar considerations hold for the design of land use policies to guide land use change. In scenarios, multiple land change trajectories occur over the EU territory at the same time, each having its own impact on multifunctionality in often opposing ways. Rather than prohibiting certain land use changes at large, it is more beneficial to guide and steer land use changes towards those locations where most beneficial impacts are expected. The analyses made in this paper provide a tool to identify such locations and measure the multiple aspects of multifunctionality related to land use change.

Conclusions
Our study has indicated that in quantifying landscape multifunctionality at EU level, the choice of indicator and scale of analysis can strongly affect the results. There is no single best multifunctionality indicator. Choosing a suitable indicator depends on the research question. It might require multiple indicators to address the different aspects of multifunctionality. We have shown that particularly the scale of analysis affects the results, and thus, interpretation of multifunctionality indicators is dependent of the scale of analysis. Particularly analyses at coarse scales, such as the administrative NUTS2 regions, may lump ES provided by different ecosystems and landscapes occurring within the administrative unit rather than referring to the multifunctionality of single land use units or landscapes. In the analysis at EU level, but also for different landscape types, we could not find a generally higher multifunctionality at locations providing higher levels  of land use diversity. The relationship between multifunctionality and land use diversity is clearly more complex and location dependent. Similar considerations hold for the diverse impacts of modeled land change in the EU on multifunctionality across scales and indicators. Large impacts in locations and regions are cancelled out at higher levels of aggregation due to the different impacts of different land change trajectories and the trade-offs between ES upon land use change. Therefore, rather than focusing on the impact of land use change on aggregated multifunctionality indicators, it is recommended to base land use policy on the location-specific impacts of land use changes on ecosystem services and functions, and societal demand for them.