Mapping Cultural Ecosystem Services Enables Better Informed Nature Protection and Landscape Management

: Cultural ecosystem services (CES) have speciﬁc richness and diversity provision patterns related to particular landscape features and land cover forms. Studies of their spatial distribution, however, are quite rare in the Slovak Republic and surrounding countries. This paper links land cover information based on an ecosystem services (ES) matrix, ﬁeld survey data and GIS method to assess CES supply in two selected Slovak regions. Our main focus is on the ecologically more valuable ‘hot-spots’ where socio-cultural values accumulate. We determined their spatial distribution, and our comparison with lower cultural value areas conﬁrmed that mountainous landscapes have the highest capacity to provide CES. This especially applies to the landscapes under National Park protection. While Slovak forests, rocks and water areas also form essential ecosystems for overall CES provision, the lowest overall capacity is in areas with residential buildings, construction, industrial and other artiﬁcial habitats. Finally, a comparison of our results with the National Ecosystem Assessment indicates that our detailed CES assessment will be more e ﬀ ective in supporting future participatory planning and management processes.


Introduction
Ecosystem services (ES) is a globally recognized concept which highlights the material and non-material benefits humans gain from ecosystems [1]. These services are combined in the following classes: (1) provisioning ES to provide services such as food and water, (2) regulating ES to monitor climate stability and disease control, (3) supporting ES to cover pollination and nutrient cycling and (4) cultural ES to provide the "non-material benefits people obtain from ecosystems through spiritual enrichment, cognitive development, reflection, recreation and aesthetic experience" [2,3]. In practice, cultural ecosystem services (CES) express appreciation of aesthetic beauty in natural scenery and inspiration for culture, art and design. They also provide a sense of place, belonging and spiritual and religious motivation. This leads to opportunities for environmental science, education, recreation and economic tourism [4].
Research into ES mapping and valuation has significantly increased in the last decade. The literature analysis [5] indicates that CES in Denmark and Spain are mostly devoted to national assessments [6,7]. However, most countries use ecosystem, land use, Corine Land Cover maps and ES assessment indicators derived from natural-environmental databases. For example, Finland, Luxembourg and the United Kingdom use the most elaborate indicator systems [8][9][10]. The ES Figure 1 depicts the two study areas assessed in the land-locked Slovak Republic in Central Europe. The Terchovská valley micro-region is in the Žilina and Dolný Kubín County in the country's north-west. This region has 17 municipalities with 33,987 ha total area and 31,699 inhabitants. The western part is mostly industrial because it borders the Žilina regional capital city with its KIA car factory and related companies. In contrast, the eastern area covers municipalities with preserved nature and historical landscape structures. These latter include dispersed settlements and terrace fields created in the 16th and 17th centuries. The protected Mala Fatra National Park is situated in the south-east of the study area, and an edge of the Kysuce Landscape protected area is in the northern region. There is also a wide range of habitats in this small area. Both these protected areas harbor valuable endemic species and glacial plant and animal relics. The Horný Liptov Region in north Slovakia has 15 municipalities and 41,414 inhabitants in its 40,560 ha. The central part is mostly industrial because it is close to Liptovsky Mikuláš city and surrounding municipalities. The industry there is mostly timber production, milk products and tourism. This region is one of the most attractive in Slovakia and it has unique natural assets and cultural monuments. The northern part contains the Tatra National Park (TANAP). This national park is a unique protected area with several endemic and protected species. The southern area has the Low Tatras National Park which is the largest in Slovakia. Its greatest attractions include numerous karstic phenomena and caves. The typical features of the Low Tatras National Park are long valleys, deep canyons and steep rocks coated by dwarf pines and grasslands.
Unfortunately, there has been extensive exploitation of Slovakia's forests. Large forest areas are rapidly becoming fragmented from timber harvesting and massive tourism. They have also suffered natural disasters. In November 2004, the Tatra National Park was decimated by a calamitous windstorm that destroyed 12,000 ha of the forest at 700-1350 m altitude. Our study area covers 3237 ha of the affected area and this accounts for almost 27% of the park's total area. This windstorm and consequent logging of damaged timber have radically changed the natural conditions for all ecosystem services [33].

Data and Methodology
We compared our regional results with the national CES assessment [32]. We chose the same three ecosystem services and the same valuation scale for the assessment that the national assessment used. The valuation methods are also similar. The main differences are our data accuracy and participatory approach. The national approach depended on the availability of national spatial data and statistics and our regional approach is based on local surveys. The primary data for our research was assembled from field research and structured stakeholder interviews. This was then supplemented by the CES matrix and land cover data. The collected field data comprised historical landscape structures, sacred, war and historical monuments and tourist infrastructure placement. This was followed by our important structured interviews, with the mayors of each municipality participating in the mapping. The Horný Liptov Region in north Slovakia has 15 municipalities and 41,414 inhabitants in its 40,560 ha. The central part is mostly industrial because it is close to Liptovsky Mikuláš city and surrounding municipalities. The industry there is mostly timber production, milk products and tourism. This region is one of the most attractive in Slovakia and it has unique natural assets and cultural monuments. The northern part contains the Tatra National Park (TANAP). This national park is a unique protected area with several endemic and protected species. The southern area has the Low Tatras National Park which is the largest in Slovakia. Its greatest attractions include numerous karstic phenomena and caves. The typical features of the Low Tatras National Park are long valleys, deep canyons and steep rocks coated by dwarf pines and grasslands.
Unfortunately, there has been extensive exploitation of Slovakia's forests. Large forest areas are rapidly becoming fragmented from timber harvesting and massive tourism. They have also suffered natural disasters. In November 2004, the Tatra National Park was decimated by a calamitous windstorm that destroyed 12,000 ha of the forest at 700-1350 m altitude. Our study area covers 3237 ha of the affected area and this accounts for almost 27% of the park's total area. This windstorm and consequent logging of damaged timber have radically changed the natural conditions for all ecosystem services [33].

Data and Methodology
We compared our regional results with the national CES assessment [32]. We chose the same three ecosystem services and the same valuation scale for the assessment that the national assessment used. The valuation methods are also similar. The main differences are our data accuracy and participatory approach. The national approach depended on the availability of national spatial data and statistics and our regional approach is based on local surveys. The primary data for our research was assembled from field research and structured stakeholder interviews. This was then supplemented by the CES matrix and land cover data. The collected field data comprised historical landscape structures, sacred, war and historical monuments and tourist infrastructure placement. This was followed by our important structured interviews, with the mayors of each municipality participating in the mapping.

Interviews and Participatory Mapping with Stakeholders
Participatory mapping enables stakeholders to express how and where they benefit from ecosystem services. This knowledge can then be used to understand their values in the area and possible conflicts between these services. It also has the potential to reveal perceptions of ecosystem services under different future management scenarios, although this has not been explored [34,35].
The 49 interviewed stakeholders included municipal mayors and their deputies, scientists working in national parks and the directors of regional tourism organizations. We conducted detailed interviews and requested respondents to identify positive and negative elements in the research areas (Table 1). We mapped the positive and negative landscape features also using GPS in the field and transformed the raw data into GIS. The locations were then transferred into QGIS 3.8 and used for later calculations [36].

Land Cover Maps and Cultural Ecosystem Services' Matrix
The provision of ecosystem services is significantly influenced by land cover and land use [12]. Therefore, we created land cover spatial layers in QGIS 3.8 for each study area using Petrovič et al.'s [24] landscape classes. Google and the 'Mapy.cz open map server' added the most current information and data [37,38]. We then verified the land cover classes in core areas of both regions during field research. Our study area had 45 land use classes in 6 groups. These groups were Forest and Woodland vegetation, Grasslands and Shrublands, Agricultural areas, Bare rock and Screes, Surface water and Settlements and Built-up Areas.
The CES matrix created by expert estimations formed the basis for the reclassification of land cover maps and in accordance with Burkhard et al.'s "ES matrix" [12]. We used this matrix as the base value and adapted it to the Slovak land use classes [39]. Each element provided three values: one for each class of CES. These expert's values were then verified with local stakeholders. In accordance with [12], each land cover element was then assigned to a 0-5 scaled land use class which measured its capacity to deliver CES. The zero value identified the lowest capacity and the value of "5" was allocated to the highest capacity. These values then provided the basis for further refinement in the following steps.

Mapping Important CES Features
To map the important CES features, we transferred our field research data to the QGIS map as spatial point, line and polygon layers. We then listed the positive and negative elements that affect ecosystem capacity in order to provide the three assessed CES (Table 1). Finally, we used these features to refine the values determined by land cover classes (see the next sub-section).

Data Processing and Evaluation
An expert estimation method assigned values to the positive and negative elements. The values depended on the extent of these elements' effect on ecosystem capacity to produce the selected CES. The values for each element ranged from 1 for low influence to 3 for high impact. We made the resulting layers of positive and negative elements as a sum of values of all partial elements. Finally, we reclassify these values to the scale (1.0-2.0) in case of positive elements and (0.5-1.0) in case of negative elements.
Spatial analysis by QGIS employed vector and raster data for the calculation. All data were converted to raster shape by the rasterize tool with 25 × 25 m resolution after the selected vector elements were redesigned by 'v.buffer'. The three layers of positive and negative elements and land cover layer with matrix data were then set for each CES category, and GRASS GIS [40] was used for all other calculations.
We then used 'r.patch' for layer-overlay. We combined the layers using 'r.mapcal' technology so that the final map was based on matrix values for the land cover layer.
Finally, we used multiplication coefficients to calculate all possible values. The matrix value of each element in a layer was multiplied by the resulting values of positive and negative element layers. The highest positive value adjusted the basic value by the coefficient of 2.0 and the highest negative value adjusted by 0.5. The coefficient of 1.0 means that no positive or negative features were identified.
Those calculations were made for each of the three CES, and the 'r.stat and r.univa' tools determined the values of the layers and their occurrence for our statistical analyses. We recorded the data using the 'r.univa' tool and interpreted it with the 'r.mapcalc' tool. We then converted the data to a percentage scale by the Ixi = (Xi − Xmin/Xmax − Xmin) * 100 formula, so that the final values of the landscape capacity for CES provision fit the 0-100 range. Finally, we reclassified the data into "hot-spots" with higher than 85% value and "cold-spots" with under 10% value. The defined hot-spot and cold-spot areas were calculated by the 'r.stat' zonal mechanism which uses category or object-oriented statistics.

Results
The specific land use land cover categories (LULC) are crucial in determining the attractiveness of the landscape for CES. The CES also informs us of the appropriate nature protection and landscape management required for a specific region. This particularly applies to recreational and nature conservation areas. The material-technical facilities and the attractiveness of cultural-historical structures are then secondary criteria in assessing the area's recreational value. Abiotic factors including morphometric relief parameters also confer importance. Our results indicate the landscape benefits to tourism from the natural conditions and preserved cultural-historical buildings, the cultural establishments and cultural and sports facilities. The resultant maps show the spatial distribution of cultural ecosystem services supply. The higher percentages indicate higher CES potential. The landscape's capacity for CES provision increases according to the area's importance for sustained nature and landscape conservation.

Recreation and Tourism
Only 1.6% of the study area contains hot-spots with very high capacity to provide eco-friendly recreation value. These hot-spots are also important for maintaining nature protection and landscape management. The areas are located in the Small Fatras, High Tatras and Low Tatras National Parks (Figure 2). The greatest number of recreation facilities are all in these areas. These include tourist trails, cycling routes, educational trails, climbing, skiing and cross-country skiing areas. Features that contribute to the attractiveness of these landscapes for visitors include the high natural value of preserved forests, alpine meadows with rich biodiversity, glacier lakes, geo-morphological localities and many wonderful observation points. However, there are limits imposed in high-value protected areas. These include seasonal closures and inaccessibility to some national park areas due to fauna and flora protection.

Landscape Aesthetics
Aesthetics are expressed as subtle, individual emotions, attitudes and perceptions of beauty. This applies to benefits accruing from landscape natural areas. It also includes their capacity to provide different kinds of scenic beauty but still maintain appropriate landscape management. These Additional hot-spots for recreation and tourism are related to cultural heritage. Historical landscape structures in dispersed settlements with preserved architecture and terrace fields are mainly used as recreational sites by urban residents [41][42][43]. The identified dispersed settlements are mostly located in the eastern parts of the Terchovská valley micro-region where Wallachian colonization occurred. The Horný Liptov region also has many war memorials which often form part of educational trails. These emphasize the objects themselves and the entire beauty of this historical area. Local castles are very important sites and significant historical objects greatly admired and appreciated by visitors. These are complemented by the interesting chateaus, bell towers, chapels and crucifixes in the natural environment.
We then selected areas up to 10% in value to highlight the cold-spots with very low capacity to provide CES (Figure 2a,b). These cover a larger area than hot-spots, with 14.5% of the entire study area. The proximity to negative landscape features decreases the attractiveness of a region. Negative features include industrial and built-up areas, opencast mining areas and highways and railways. These have little or no potential to develop recreation and tourism in their immediate vicinity. Further negative areas are found in forests affected by thunder-storms. These have unfavorable conditions for CES provision. From a long-term perspective, these areas can return to their natural state and then have very high potential to provide CES. However, their actual capacity to provide opportunities for recreation and tourism remains very low.
The most cold-spots in the Terchovská valley micro-region are located in the western, more industrial part of the study area (Figure 2a). Others are scattered in the north and east with intensive timber harvesting. However, those in the Horný Liptov region mostly surround the Liptovský Mikuláš city industrial sites. In addition, the High Tatras and Low Tatras National Park cold-spots are caused by timber harvesting initiated to limit bark beetle infestation (Figure 2b).

Landscape Aesthetics
Aesthetics are expressed as subtle, individual emotions, attitudes and perceptions of beauty. This applies to benefits accruing from landscape natural areas. It also includes their capacity to provide different kinds of scenic beauty but still maintain appropriate landscape management. These include the panoramic views and photogenic spots provided by Slovak water bodies, historical landscape structures, rock formations and wonderful mountains. In addition, landscape capacity relies on expert-based land-use ranking by ES matrix. The eco-service supply is then based on the visible properties of observation points, especially those along recreation routes.
Our results indicate that the landscape capacity for this ES is substantially higher in the core and buffer zones of protected regions than in surrounding areas ( Figure 3). Hot-spots provide very high capacity to provide aesthetic value, but these are present in only 4.1% of the study area (Figure 3a,b). They include the water bodies, rock formations, forests and grassland in national parks and historical landscape structures and cultural heritage sites. In contrast, landscape cold-spots cover 12.9% of the total area, and they have low to very low capacity to provide aesthetic value. (Figure 3a,b). Cold-spots in the study area includes some agricultural land, built-up areas, industrial and commercial units, dumpsites and mining and forest calamity areas. This contrasts sharply with the National ES assessment which does not acknowledge any of these negative CES values. The national assessment uses different methodologies and emphasis and it has a lower level of detail.

Natural and Cultural Heritage
We used a limited number of indicators to assess natural and cultural heritage. Our results identify landscapes with selected areas' cultural and natural heritage hot-spots and cold-spots (Figure 4). Only 0.8% of the total study area has a very high potential to provide this CES.

Discussion and Conclusions
Spatially explicit studies are rare in the cultural ecosystem service literature. There is little documentation on how specific landscapes and their features and ecosystem properties contribute to public perception. This absence limits possible benefits for landscape planning and nature protection management [44,45]. Our approach combines analysis and expert estimation of the landscape's capacity for the provision of recreation facilities and its aesthetic beauty. It also includes the values The most valuable parts in Figure 4 comparison are the protected areas' positive elements of water bodies, rock formations, historical monument and historical landscape structures. Figure 4a shows that the hot-spot density of cultural and natural heritage is higher in the Terchovská valley micro-region. It highlights that 0.8% of this area has a higher number of dispersed settlements and nature protected areas and less negative sites. Figure 4 also shows the natural and cultural heritage cold-spots which cover 11.0% of the two areas. These cold-spots particularly affect the capacity of large-block agricultural fields to supply natural and cultural heritage ES. They include the dump-sites and damaged forest from calamitous wind-storm areas not included in the national assessment. Finally, our closer examination of the spatial relationships between heritage values, nature protection, land use and landscape management will help to improve the understanding of cultural ecosystem services in both the scientific community and the general public. This includes the historic influence on landscapes and the enhanced communal sense of place and nature protection.

Discussion and Conclusions
Spatially explicit studies are rare in the cultural ecosystem service literature. There is little documentation on how specific landscapes and their features and ecosystem properties contribute to public perception. This absence limits possible benefits for landscape planning and nature protection management [44,45]. Our approach combines analysis and expert estimation of the landscape's capacity for the provision of recreation facilities and its aesthetic beauty. It also includes the values of the positive and negative landscape elements and utilizes important land use and land cover maps (LULC). This approach most closely resembles national ecosystem services assessments in Slovakia [32] and Russia [13]. However, we provide more detailed data from land cover, and we include the positive and negative elements that influence the selected cultural ecosystem services (CES).
Our mapping of the ecosystem potential for tourism and aesthetic value defines the extent to which Slovakia can potentially supply outdoor recreation services in selected areas. This potential is highly dependent on the spatial distribution of nature protection areas and its limits. It also depends on the quality of natural and semi-natural vegetation. Our results indicate high overall potential, with almost a third of the territory falling in the highest recreational capacity classes with several hot-spot areas. In contrast, the negative elements give uneven spatial distribution. In addition, 38% of the study area has cold-spots in regions with residential buildings, industrial and other artificial habitats, and these have very low CES capacity.
The results highlight the high tourist and aesthetic value of the traditionally used landscapes in the Terchovská valley micro-region dispersed settlements and the nature protected Tatra National Park. Coniferous forests, mixed and deciduous forests and rocks, screes and water areas are essential ecosystem types for national CES provision. However, these landscapes and ecosystems are currently at risk because of the labor-intensive care required and their low economic return. They are also vulnerable to future degradation, especially in remote municipal areas. Therefore, improved management is essential to ensure their future existence and negate the poor spatial extent of traditionally used landscapes. This will also help restore large damaged areas in nature parks and offset massive tourism expansion. This paper presents the first comparison of Slovak national ecosystem services assessments [32] with a more detailed approach to recreation and landscape aesthetics. These assessments highlight that limitations in the present methods must be counteracted with improved tourist and recreational facilities. These include tourist rest-areas and seating-benches and more accessible pathways and scenic views. In addition to the current land use and land cover detail, future LULC maps should address changes in land use which reflect the low CES values of increased industrial zones, transportation, built-up areas, deforestation and in-growths.
We identified inconsistencies in the comparison of our results and the national ecosystem services assessment [32]. Our approach and results are more critical and stricter than the national assessment. We used specific land use maps and considered the accurate spatial distribution of recreational and cultural-historical attractions in both study areas. Our participatory mapping provides a great amount of information for use in CES mapping at the local and regional scale. A further significant difference between these approaches relates to the mountain areas. From the national perspective, it appears that larger mountain forest areas have a higher capacity to provide CES than our research showed. This is because national assessment uses less detailed and up to date data. Our results confirm that it is imperative to use actual and more detailed data for landscape planning and management at the municipal and regional levels than the data which is available at the national level.
Our approach provides a framework for CES assessment and identification of the areas where the supply of services should be restored and special management implemented. The framework is also adapted to respond to over-reaching national policy and nature-protection goals. The approach aims to promote sustainable regional tourism planning and the cartographic outputs provide additional information on the spatial patterns of CES supply. These patterns particularly include the hot-spots and cold-spots and their altitudinal distribution throughout different regions, as well as the mountain landscapes with their high CES provision.
In addition to other benefits, greater CES knowledge can help both landscape planning and regional development strategies. Our approach to landscape planning reveals values that stakeholders attach to landscapes. These especially include social values often ignored by biophysical and economic ecosystem service assessments. These values were not previously considered in most Slovak local and regional economic and social development planning documents. However, they can provide a significant incentive to owning, using, managing and protecting land. For example, the facilities and activities stimulated by CES in our study areas include the growth of private and public nature reserves, tourism facilities, hobby farms, residential houses and country cottages.
People already perceive and use cultural services but they should understand that the potential of these services is related to the natural and socio-economic environment. Our approach identifies the regions that can promote new activities, such as building new tourism facilities. Alternative landscape protection and maintenance of endangered and protected areas can also be prioritized. Finally, the implementation of these results in the previously mentioned development documents could promote planning and future development in each municipality.
Enterprising improvement can also be initiated by incorporating data from the online Public Participatory GIS (PPGIS) survey. This examines inhabitants' CES perceptions and it can assess the actual regional CES supply and demand and the participant consumption of these ecosystem services. Finland [8] has already implemented this approach, and has arguably the most elaborate indicator system. The proposed methodology enabled a greater analytic scale, where information on recreation facilities can be included as elements that attract and promote the flow of visitors [41].
Recent scientific developments have also acknowledged social media potential for collecting data on recreation activities and landscape aesthetic beauty. Here, photographs uploaded to image-sharing websites provide information on the visited sites' spatial distribution, visitor appreciation and locations of the most accessible sites [46][47][48][49][50][51][52][53]. This information can be used in modeling the density of infrastructure for recreation and valuable aesthetic sites. It can also monitor areas that need extra protection to maintain the study regions' natural and cultural beauty.
Principally, our results provide the most important information base for design, regional implementation and evaluation of land use and nature protection policies.
In addition, traditional landscape management can benefit both the state and society with sustainable and environmentally-friendly agro-tourism. Finally, the landscape can be redesigned so that it is more accessible for experiencing all the cultural ecosystem services available in these study areas.