Ecosystem services are inherently defined mutual interaction between ecological and social systems because only those ecosystem processes that contribute to the fulfillment of human needs, are defined as ecosystem services    . The ecosystem services approach offers the ability to explore the influence of land use and practices on natural capital stocks, on the processes that build and degrade these stocks, and on the flow of ecosystems services from the use of these stocks . Service-providing ecosystems are denoted as natural capital  . Ecosystem services linked to natural capital can be divided into three main categories (provisioning, regulating and cultural)  . Agroecosystem services are the multiple goods and services provided to humanity by nature in combination with additional anthropogenic inputs in agricultural systems . Concept of agroecosystem service combines environmental and socio-economic approach to the analysis and evaluation of natural capital. For practical use as well as legislative use of the concept of agroecosystem services in planning and prospective studies, qualitative and quantitative analysis and evaluation of agroecosystem linked to spatial visualization at the required level are essential . The concept of agroecosystem services receives increasing attention in recent years with an emphasis on national assessment. Following this trend, an integrated evaluation and mapping of agroecosystem services was carried out in the Slovak Republic .
Agricultural systems are intensely managed by humans and are more controlled and regulated than majority of the other ecosystems . Traditionally, agroecosystems have been considered primarily as sources of provisioning services, but more recently their contributions to the other types of ecosystem services have been recognized   , such as regulating and cultural services   . Recent conceptual works used the ecosystem services approach to highlight the importance of pedosphere to the human well being and prosperity  . We consider the agroecosystem not only as means of production but also as a part of the natural environment, where the pedosphere has functions other than food production  . The agroecosystem based on soil is multifunctional in all conditions, both in terms of its processes and functions and services . The cascade model developed by Haines-Young and Potschin  demonstrated how soil functions can contribute to ecosystem services. According to Greiner et al.  soil and their functions are critical to ensure the provision of ecosystem services. Nevertheless, agroecosystem services at the national level, were not assessed on the context of the Slovak Republic. In Slovakia, Juráni  , Bujnovský et al.  , Makovníková, Barančíková and Pálka  and Barančíková et al.  have been also dealing with assessment of soil functions that are the basis for the evaluation of agroecosystem services. These authors define a minimum collection of soil indicators for sufficient assessment of soil functions. These indicators form our basis for agroecosystem services evaluation.
While the concept of natural capital and economic services are widely accepted and their potential contribution to better environmental management is also acknowledged    , their practical applications such as distribution and mapping are still insufficient and limited. The models and their map presentations should reflect biophysical factors on one hand, but on the other hand they should be applicable within administrative units for better application of models in decision-making . Ecosystem services and natural capital are inherently spatial by nature although some services are easier to map then others. According to Costanza et al.  , landscape metrics quantifies physical landscape structures which themselves determine processes and functions. Landscape metrics has therefore to be considered as meaningful parameters together with the others in ES mapping and evaluation . Ecosystem services evaluation and mapping are useful for assessment of landscape capacities and potentials to supply ecosystem services  and to adapt the management to local conditions  . These maps are also important to assess spatial trade-offs among ES and synergies among multiple ES . Because the provision of ecosystem services depends on biophysical conditions and changes over space and time due to human induced land cover, land use and climatic changes   , the supply and demand of services may differ geographically . A number of recent studies have evaluated and mapped the supply of services at global  , continental  , national     or regional scales. For this reason, it is necessary to understand better where and what services are provided by a local area of land or region . The most common indicators for mapping ecosystem services are land cover, soils, vegetation and nutrient related indicators. According to a review of different approaches used to model ecosystem services, regulating services were the most commonly modeled, followed by provisioning and cultural services and supporting processes . Provisioning services strongly depend on the extent of managed land and the land use intensity. However, for quantification of productivity or food provision, further information, for example on soil quality or soil management, location or climatic zone, is essential to derive a valid estimation on food provision . Cultural ecosystem services are less in the foreground to be put on maps, because researchers must rely on proxies for their quantification . Ecosystem services models can vary from simple expert based scoring systems to complex ecological models cycles of carbon, nitrogen and water . Linking ecosystem services models to landscape coverage is one of the conditions for using these models to monitor changes in land use management, spatial planning, and implementation of the assessment of the potential of natural capital services in socio-economic planning within the region and landscape . Burghard and Maes  highlights the importance of the process of mapping ecosystem services that can be used to visualise impacts which are often considered invisible externalities of agriculture, both positive and negative.
The aim of this study was to assess and map agroecosystem services on national level and to evaluate agroecosystem services along the climatic gradient in Slovak Republic. The paper aims to describe the use of GIS techniques in creating a uniform spatial unit for agroecosystem services inventory and compare the spatial configuration, synergies and trade-offs of agroecosystem services along the climatic regions.
2. Material and Methods
2.1. Study Area
Slovakia is a land-locked country in Central Europe between latitudes 47˚ - 49˚N and longitudes 15˚ - 21˚E. The average rainfall in lowlands is about 600 millimeters per year, in midlands about 700 millimeters per year and the biggest average rainfall rate belongs to mountain areas―approximately 1500 millimeters per year. The daily average temperatures in winter are around freezing (0˚C or 32˚F), while in summer they are around 13/15˚C (55/59˚F) at night, and 25/27˚C (77/81˚F) during the day. Its terrain is mostly hilly, upland and mountainous in the central, north and north-eastern parts of the country where the permanent grasslands mainly occur . For this study, we used a classification of agro-climatic regions provided by the Information Service of the National Agricultural and Food Centre―Soil Science and Conservation Research Institute . In this classification, 11 agro-climatic regions were identified according to long-term average temperatures in January, average growing-season temperatures, daily average temperatures sums (T > 10˚C), the length of period with daily temperatures td > 5˚C and the climatic moisture indicator according Budyko calculated by Tomlain . For our purpose, the original vector layer with 11 categories were merged into 4 categories (climatic regions: very warm (VW), warm (W), moderate warm (MW) and moderate cool climatic region (MC)) and transfered to a raster with a resolution of 100 m (Figure 1).
Agroecosystems, which occupy 49.3% of the Slovak republic, were assessed. To estimate the surplus area of agroecosystems (arable land and permanent grasslands) in each of the climatic regions (Table 1) we used data provided by Land Parcel Identification System (LPIS).
The structure of the assessment is given by one ecosystem type―agricultural ecosystem and 6 ecosystem services delivered from this ecosystem (provisioning services, regulating services―regulation of water regime, regulation of water erosion, cleaning potencial of ecosystem, regulation of climate and cultural services-potencial of outdoor recreation). Ecosystem types are further classified into two ecosystem categories based on the management, arable land and permanent grassland.
To the primary geo-referenced data belongs the Digital database of soil profiles of Geochemical atlas of Slovakia (GchA-2965 localities on agricultural soil). Geochemical atlas database contains data of agrochemical soil properties and risk elements concentration determined from the samples collected during the national project “Geochemical atlas of soils of Slovakia” . Further, there are geo-referenced data, as a source of additional information (data of relevant soil properties) on the primary geo-referenced data. As secondary geo-referenced data was used Digital database of Soil monitoring of Slovakia (CMSP) represented by digital data archive of the “Partial monitoring system-Soil” as a part of the Complex Environment Monitoring of Slovakia . Soil indicators are included in the soil monitoring system in Slovakia according to the recommendation of the European Commission for united soil monitoring system in Europe . All monitored indicators are quantifiable.
Table 1. Surface areas: arable land and grassland in climatic regions (ha).
Figure 1. Categories of climatic regions in the Slovak Republic.
2.3. Methods of Assesssment and Mapping the Potential of Agroecosystem Service
Ecosystem services potential (capacity) has been characterized by Burghard et al.  as the hypothetical maximum yield of selected ecosystem services. Integrated approach towards the assessment of ecosystem services with the application of expert methods enables to link data on soil parameters, morphological and biophysical parameters of soil with the data on soil exploitation into one concept, i.e. the concept of agroecosystem services. Biophysical spatial quantified data is used by many authors to evaluate ecosystem services in their work   . For analysis and evaluation of potential of agroecosystem services we have created a mapping unit combining these input layers: slope topography, soil texture and usage of land in four climatic regions. The resulting layer is elaborated at the national level for the whole territory of the Slovak Republic using the methods and tools offered by geographic information systems (GIS). We used ESRI’s software package ArcGIS for Desktop Advanced version 10.3. Mapping units are also compatible with the spatial units in international database (Corine Land Cover) which carries information on the use of land. Within each unit, we calculated a weighted average of the potential of each agroecosystem service for the territory of the Slovak Republic, which is characteristic for the spatial aggregate of functional unit. Forouzangohar et al.  concluded that the best sampling and mapping strategy would be a regular grid scheme.
Moreover, the assessment and mapping methodology developed in this paper is replicable and could be applied by planners in case they are proficient in geographical information systems at different levels, from regional (district) to national level.
2.4. The Provisioning Services
The basis for analysing the potential for the provisioning agroecosystem services was a point value within a range of 0 - 100 that indicate a productive potential based on typological and production classification of agricultural soils of Slovakia  . Point value is the sum of points assigned to the average yield of crop production according to the soil type (0 - 60 points), slope and its aspect (0 - 15 points), stone content and soil depth (0 - 15 points) and soil texture (0 - 10 points). Values were categorised into five groups (range of 20 points) and the categories of provisioning service potencial are as follows: 1) very low potential (lower than 20 points), 2) low potential (20 - 40 points), 3) medium potential (41 - 60 points), 4) high potential (61 - 80 points), 5) very high potential (more than 80 points).
2.5. The Regulating Services
Potential of regulation of water regime (soil water storage) was obtained from maps and databases  (Bujnovský et al., 2009). Its values are given in mm and are determined on the basis of the value of retention water capacity recalculated to soil water storage in context with the soil depth. Values were categorised into five groups and the categories are as follows: 1) very low potential (<135 mm), 2) low potential (135 - 175 mm), 3) medium potential (175 - 215 mm), 4) high potential (216 - 275 mm), 5) very high potential (>275 mm).
2.6. Potential of Regulation of Soil Erosion, Regulation of Water Erosion
Regulation of water erosion was derived from maps and databases based on empirical model of the universal soil loss equation―USLE  (Wischmeier and Smith, 1978), Styk and Pálka, 2005). The relative ratio of the calculated values of soil loss and acceptable erosion expresses the degree of soil erosion endangerment (SEOP value,  ). Values were categorised into five categories: 1) very low potential (more than 2.60), 2) low potential (2.21 - 2.60), 3) medium potential (1.81 - 2.20), 4) high potential (1.40 - 1.80), 5) very high potential (less then 1.40).
2.7. Cleaning Potential of Agricultural Land Ecosystem
Cleaning potential of ecosystem of agricultural land depends on the actual soil contamination and potential of soil sorbents and was calculated as accumulative function:
Cleaning potential = Sorption potential of soil + Potential of total content of inorganic contaminants evaluated according to The Slovak Law 220/2004 Coll. (the method is mentioned detaily in our previous article . Point evaluation of Sorption potential of soil was calculated as a sum of quality factors (pH (0 - 4 points), (0 - 1 points) and quantity factors (Cox (0 - 1 points), H-depth of humus horizon (0 - 2 points)) according to the function: PS = F(pH) + F(Q46) + F(Cox) × F(H). The overall rating is determined as a sum of the soil contamination and Sorption potential of soil. The high soil contamination was evaluated by the high point value and present high risk (0 - 5 points). On the other hand, high soil Sorption potential results in low point value and decreases possible transport risk of harmful elements in soil. Sum of values were categorised into five groups and respective categories are as follows: 1) very low potential (more than 6.50 points), 2) low potential (5.51 - 6.50), 3) medium potential (4.51 - 5.50 points), 4) high potential (3.50 - 4.50 points), 5) very high potential (lower than 3.50 points).
2.8. Climate Regulation
In agroecosystems of agricultural land, soil organic matter represents the largest share of total organic carbon found in the soil. Agroecosystems contribute to climate regulation by sequestration of organic carbon in the soil. Soil organic carbon stock (SOCS) was calculated as a function:
Soil organic carbon stock (depth 0 - 30 cm) in t・ha−1 = 10 × (BD (0 - 10 cm) × SOC (0 - 10 cm) + BD (10 - 20 cm) × SOC (10 - 20 cm) + BD (20 - 30 cm) × SOC (20 - 30 cm)), BD―soil bulk density in g・cm−3 , SOC―soil organic matter content in % . The categories are as follows: 1) very low potential (lower than 58.00 t C ha−1), 2) low potential (58.00) 62.00 t C ha−1, 3) medium potential (62.01 - 67.00 t C ha−1), 4) high potential (67.01 - 72.00 t SOC ha−1) 5) very high potential (more than 72.00 t SOC ha−1).
2.9. Cultural Ecosystem Services
All agro-ecosystems are considered to be potential providers of such services. Agroecosystems also have the potential for providing cultural ecosystem services, particularly recreational activities linked to natural resources, such as hiking, biking, cross-country skiing. Recreation potential was evaluated through agroecosystem landscape components that have a specific link with summer, winter and year-round recreation. The capacity of ecosystems to provide recreational services depends on particular uniqueness of the site, its accessibility and the surrounding infrastructure. Point value is the sum of points assigned to the altitude (1 - 3 points), inclination (1 - 3 points), drainage and precipitation and temperature (1 - 4 points) and distance to the roads (1 - 4 points). The recreational potential for all these activities was calculated as the sum of potentials for individual recreational activities without added points (Natura 2000). These were added only to the final sum in order to prevent multiple evaluations of additional factors. This method is described in detail in our previous article  , where we compared the use of this model and SolVES model (according to Sherrouse et al.,  ) in the case study. The categories of agroecosystem to provide outdoor recreational activity are as follows: 1) very low potential (lower than 6.00 points), 2) low potential (6.01 - 9.00 points), 3) medium potential (9.01 - 12.00 points), 4) high potential (12.01 - 15.00 points), 5) very high potential (more than 15.00 points) and:low relevant capacity (from 2.01 to 3.09 points), medium relevant capacity (from 3.10 to 5.09 points), high relevant capacity (from 5.10 to 7.09 points), very high relevant capacity (higher than 7.10 points).
3. Results and Discussion
Evaluation of agroecosystem services linked to spatial visualization allows to optimize the management of agroecosystems, thereby to promote synergies between ecosystem functioning and the social dynamics of respective region. Explicit modelling of agroecosystem services is considered to be one of the main requirements for implementation of the concept of these services in institutional decision-making. Krkoška, Lorenzová et al.  and Kanianska et al.  identified land use and land cover as a major driving force leading to the differencies in agroecosystem services potential. Another important driving force affecting the distribution of agroecosystem services is climatic region  , that impacts distribution of these services as well as interaction between them   . The correlations between climatic region categories and agroecosystem services for arable land and grassland in the Slovak Republic shows in Table 2.
Table 2. Correlation between clima categories and agroecosystem services for arable land and grassland.
Significance labels: ***p < 0.001, **p < 0.01, *p < 0.05, ns: non significant.
Positive correlation coefficient indicates the positive effect of cold climate zone, while negative coefficient indicates the positive influence of warm climate zone. The warm and dry lowland region has higher production potential, water regime regulation, pollutant filtration, and soil drainage in comparison to slightly warm to cool regions. These results are consistent with the place of occurence of soil, its properties, processes and functions in the concept of agroecosystem services   .
The potential of the provisioning service of agricultural land is determined by its location in the landscape with the climatic conditions (temperature and precipitation) and it is a combination of abiotic, biotic, morphological and socio-economic factors (Figure 2). In Slovakia 29.14% of ecosystems of agricultural land has very high potential for provisioning services. They are mainly ecosystems of arable land, located in warm climate in the Danube basin. The high level of provisioning services (crop poroduction) indicates Spake et al.  in his work for widest lowland valleys in North Alps, where agricultural areas are overrepresented. The main limiting features of low potential of provisioning services are relatively cool climate and considerable sloppiness, limited depth or relatively high stoniness of soils in such ecosystem. In this category the predominant ecosystems are permanent grasslands. Our results show that the proportion of the category with higher potential of the provisioning service increases from the cool climate to very warm (Figure 2) (for arable land from cool climate with prevailing low potential to very warm with very high potential, for grassland from cool climate with lowpotential to very warm with high potential of provisioning services).
In Slovakia 27.47% of the area of agricultural ecosystems has very high potential for regulation of water regime (accumulation of water in the soil) (Figure 3).
They are mostly ecosystems of arable land located in Eastern Slovak Lowland, Danubian Upland, South-Slovak Basin and Košice Basin with heavy clay loam and clayey deep soils without skeleton. In Slovakia out of the total area of agricultural land the highest proportion (35.96%) have ecosystems with high potential for regulation of water regime (Figure 3). Ecosystems with low potential for water storage occupy 21.21% and they are located on deep to moderately deep, light soils without skeleton and on moderately heavy, slightly
Figure 2. The potential of the provisioning service of agricultural land (arable land―AL―Figure to the left, grassland―G―Figure to the right), percentage of the area of each category of potencial in four climatic regions. Explanations: the categoriesof potential: 1) very low potential, 2) low potential, 3) medium potential, 4) high potential, 5) very high potential.
Figure 3. The potential of regulation of water regime of agricultural land (arable land―AL―Figure to the left, grassland―G―Figure to the right), percentage of the area of each category of potencial in four climatic regions. Explanations: the categoriesof potential: 1) very low potential, 2) low potential, 3) medium potential, 4) high potential, 5) very high potential.
to moderately skeletal soils. The greatest influence on water storage potential in both ecosystems has climate, but the impact of soil texture is also significant. Potential of water regime regulation within both types of ecosystems increaces from the cooler to the warmer climate and the highest potential is most represented in areals with warm climate, where deep clay soils developed on clay sediments of former seas and lakes, as well as rivers (alluvial cones, aggradational levee) occur and are situated in the foothill parts of lowlands and in the basins .
Agroecosystems of arable soils have a high to very high potential for regulation of soil erosion, regulation of water erosion (92% out of the total area of arable land). Arable land is located mainly in flat areas where low risk of water erosion occurs (Figure 4). These values of soil erosion significantly correlate with the attribute of slope (r = −0.72, p = 0.0008). Another prerequisite for higher potential for regulation of water erosion on arable land is presence of deep soils, and consequently higher limit for acceptable soil loss. When considering the overall coverage of land by permanent grassland (land registered in LPIS as permanent grasslands), the potential for soil erosion control achieves very high levels (total area of permanent grasslands―100%, is classified to the category with the very high potential) . The potential to regulate soiltransport is increasing with warmer climate in case of arable land as well as permanent grassland. This is related to the occurrence of deep soils in lowland areas, where warm and very warm climatic regions predominate, and the limit for acceptable loss of soil is higher. In case of arable soils, there is a potential for regulation of soil erosion at approximately the same level in moderate cool (MC) and moderate warm (MW) climatic zone (Figure 4). In general, grassland potential to regulate soiltransport rises from the MC climatic zone to the VW climatic zone. The grassland potential is over 80% in the VW climatic zone (Figure 4).
Cleaning potential of ecosystems in agricultural land depends on the potential for contamination and potential of soil sorbents with high affinity to inorganic pollutants. Out of the total agriculrural land in Slovakia, 41.67% of ecosystems have very high potential for soil cleaning (immobilization of inorganic pollutants). They are mainly ecosystems of arable land with high carbonate content developed on loess, located in the Danube and the Eastern Slovak Lowlands without any anthropogenic and geochemical depositions (Figure 5). This is based on agrosystems with optimal soil parameters in relation to the ecosystem filtration service    . Ecosystems of arable land of low potential (41.12% of the area) are developed on fluvisols (along Váh River, Hron River and Bodrog River). Low cleaning potential of these ecosystems is due to a higher number of risk elements in alluvial sediments, anthropogenic deposition and low potential of soil sorbents (low pH, low content of carbonates, and low content of organic matter of lower quality). However, grassland ecosystems are also strongly involved in low category of cleaning potential. Permanent grasslands mostly use farmland ecosystems located at higher altitudes as well as steeper slopes, on soils with lower sorption potential, and on soils developed on substrates with higher content of risk elements. Areas occuring in moderately cold to moderately warm climatic regions, with steeper slopes, and a higher percentage of clay particles have very low and low potentials. Very warm climatic region, mild slope and medium content of clay fraction prevail in areas with high and very high filtration potential.
Figure 4. The potential of soil erosion in agricultural land regulation (arable land―AL―Figure to the left, grassland―G―Figure to the right), percentage of of each category area for potencial in four climatic regions. Explanations: the categoriesof potential: 1) very low potential, 2) low potential, 3) medium potential, 4) high potential, 5) very high potential.
Figure 5. The cleaning potential of agricultural land (arable land―AL―Figure to the left, grassland―G―Figure to the right), percentage of the area of each category of potencial in four climatic regions. Explanations: the categoriesof potential: 1) very low potential, 2) low potential, 3) medium potential, 4) high potential, 5) very high potential.
Carbon stored in ecosystems is an important indicator of regulation services potential  whose amount depends on land use and land management practices . In agroecosystems of agricultural land, soil organic matter represents the largest share of the total soil organic carbon stock. Agroecosystems contribute to climate regulation by organic carbon sequestration in soil. Soil organic carbon content can be used to represent the carbon sequestration and can be used as an indicator in climate regulation potential service . The results of percentage distribution of various categories of potential of climate regulations are significantly influenced by the ecosystems of arable land due to high share of area of these ecosystems in the total area of agricultural land. Out of the total area of agroecosystems for agricultural land up to 94.83% belong to the category of low potential of climate regulations (Figure 6). Agroecosystems of arable land located in lowlands are characterized by low potential for climate regulation and low average stocks of soil organic carbon. In higher altitudes, the average organic carbon stocks, and thus the potential of climate regulation is slightly rising.
The low potential of climate control for agroecosystems of arable soils is mentioned in Burghard et al.  work. The average stock of soil organic matter in arable land of Slovakia (at a depth 0 - 30 cm) ranges from 59 t C ha−1 (soils in altitudes from 0 to 300 m a.s.l.) to 67 t C ha−1 . Carbon sequestration in arable soils is lower compared to the grassland  within the same soil type, therefore, in case of agroecosystems of arable soils there are no categories of high and very high potential of climate control. Cambel and Souster  and Schnitzer  report that intensive soil management reduces the amount of organic matter in the soil . The average stock of soil organic matter in ecosystems of slovakian grasslands (0 - 30 cm deep) is ranging from 73 t C ha−1 (soil at altitudes from 0 to 300 m nm) to 86 t C ha−1 (land at an altitude above 600 m a. s. l.) . The high and very high category of climate regulation is up to 58.04% of the agricultural grasslands area. In case of grassland ecosystems, cambisols, which are the most widespread soil’s type in Slovakia  , are dominant soil type at higher altitudes (>300 m nm). Higher representation of the high category of climate control was also established in warm climatic zone (Figure 6). At elevated average temperatures, soil organic matter stock in arable soils decreased, as organic carbon inputs of plant residues and manure are not able to eliminate carbon losses due to faster mineralization  . Agroecosystems of arable soils in colder regions have a higher potential for climate control.
Agroecosystems also have the potential for providing cultural ecosystem services, particularly recreational activities linked to natural resources, such as hiking, biking, cross-country skiing. The capacity of ecosystems to provide recreational services depends on particular uniqueness of the site, its accessibility and the surrounding infrastructure. Agroecosystems of arable land are predominantly of very low to low potential of natural conditions for recreation. In Slovakia 53.82% of permanent grassland area has high and very high natural conditions potential for outdoor recreation. On the contrary to ecosystems of arable land, grassland agroecosystems are located close to the protected Natura 2000 sites in areas with steeper sloppiness and at higher altitude. Agroecosystems have predominantly very low to low potential of natural preconditions for recreation (Figure 7). The perception of the landscape in intensively and extensively farmed areas is also different, which, according to Martens  , may negatively affect the attractiveness of intensive agricultural regions for recreation. The overall assessment of the potential of natural preconditions for recreation in % of agricultural land is strongly influenced by the ecosystem of arable soils (Figure 7), which has higher share in total area of evaluated agroecosystems. According to Burkhard et al.  , agricultural use of grassland has predominantly medium potential to provide recreational services (outdoor activities).
Trade offs and Synergies between Agroecosystem Services
Ecosystem services are non-lineary linked and changes in one service can impact the others in positive or negative way   . Spake et al.  define synergies when multiple services are enhanced simultaneously and trade-offs when the provision of one service is reduced due to increased use of another service. For identifying and mapping agroecosystem services associations with regards to known land use (arable land or grassland), climatic zone PCA analysis and Spearmans pairwise correlation analysis were used to evaluate the relationships between ecosystem services  .
The first axis (Figure 8) of principal component analysis biplot for agroecosystem services of arable land represents a spatial trade-offs between climatic zones and provisioning services, regulation of water regime, pollutants filtration and regulation of soil erosion. The second axis PC2 describese a synergy between climatic zones and clima regulation and their trade-offs with potential of natural conditions for recreation. The first two components accounted for 70% of the total variation in agroecosystem services. The trade-offs between climatic zones and provisioning services and regulation of water regime are more pronounced in the case of grassland.
Positively correlated agroecosystem services are assumed to be synergistic, while negative correlations infer trade-offs   . The relationship of individual services expressed through correlation coefficients is given in Table 3.
In the case of arable land trade-offs between the potential provisioning services and potential of climate regulation as well as potential of natural conditions for recreation was determined. The synergistic effect is between the provisioning agroecosystem service, regulation of water regime and soil cleaning (immobilization of inorganic pollutants). The similar relationships have been established for permanent grassland. However, in this land use management, there was no trade-offs between the potential provisioning services and potential of climate regulation as well as potential of natural conditions for recreation. The synergistic effect of regulation of water regime and soil state in its work Lescourret et al. . Statistically significant synergistic effects are determined within each climatic area. Synergistic relationship between the potencial of erosion regulation and the potencial of water regime regulation as well as provisioning services has been established in all climatic regions. Only in moderate warm and very warm climatic region we have determined synergistic relationship between the potencial of provisioning services and the potencial of soil cleaning (immobilization of inorganic pollutants). In the moderate warm climate zone and warm climate zone, the potencial of provisioning services and the potencial of water regulation are also mutually beneficial.
Figure 6. The climate regulation potential of agricultural land (arable land―AL―Figure to the left, grassland―G―Figure to the right), percentage of area for each category of potencial in four climatic regions. Explanations: the categoriesof potential: 1) very low potential, 2) low potential, 3) medium potential, 4) high potential, 5) very high potential.
Figure 7. The potential for providing cultural ecosystem services of agricultural land, percentage of the area of each category of potencial in four climatic regions. Explanations: the categoriesof potential: 1) very low potential, 2) low potential, 3) medium potential, 4) high potential, 5) very high potential.
Figure 8. PCA analyse Biplot for arable land PCA analyse Biplot for grassland. Labels: PAESS―provisioning agroecosystem, RWR―Regulating water regime, RSC―regulating soil cleaning, RSE―regulating soil erosion, RC―regulating clima, CAESS―cultural agroecosystem services.
Table 3. Spearman correlations coefficients between agroecosystem services.
Significance labels: ***p < 0.001, **p < 0.01, *p < 0.05, ns: non significant.
This study applies a first assessment of agroecosystem services in the Slovak Republic and allows us to link the analysis of land use and differences of particular agroecosystem services in climatic regions. Evaluation of agroecosystem services linked to spatial visualization allows to optimize the management of agroecosystems, and thereby to promote synergies between ecosystem functioning and the social dynamics of the region. Provisioning, regulating and cultural ecosystem services of agricultural land are analyzed, modelled and evaluated in spatial grid scheme, which is replicable and could be applied by planners in case that they are proficient in geographical information systems. This proposed mapping system can also be used to assess agro ecosystem services in the regions and districts in the Slovak Republic that provide guidelines and limits for policy development on land management and land use changes at local and regional levels. Applying the agroecosystem service concept can help to show the effects of land use, climatic conditions as well as human interventions by qualitatively and quantitatively analyzing trade-offs between different services and by supporting the development of site-specific, more sustainable land use strategies. According to Burghard et al. 2013 and Krkoška, Lorencova et al. (2016), the mainstreaming of agroecosystem services into national policy and decision making needs to be futher supported by assessments based on local or national data and more accurate modelling approaches.
This study suggests that climate has the most significant impact on agroecosystem services. Warm, dry lowland region has a higher potential of provisioning services, regulation of water regime, filtration of pollutants and control of soil erosion in comparison to moderately warm and cold regions. In moderate cold region, more than 90% of the total arable land area has low potential for regulation water regime and cleaning potential (immobilization of risk elements). In moderate warm climatic region, there is a high share of categories of low and moderate potential of provisioning services and low and moderate potential of regulation of water regime. Majority of the total area of warm climatic region belongs to the categories of moderate to high potential of provisioning services and high potential of regulation of water regime. In this climatic zone categories of low potential of risk elements immobilization is present in more than 65% of the total arable land area. On the other hand, in very warm climatic zone, more than 89% of the total arable land area belongs to the category with very high cleaning potential (immobilization of risk elements). Potential of natural conditions for recreation is higher only in moderate cold and moderate warm climatic zones with higher proportion of grassland area agroecosystems and protected areas NATURA 2000.
Ecosystem services are non-lineary linked and changes in one service can impact the other in positive or negative way. High potential of provisioning service is linked to the high potential of regulation of water regime, pollutants filtration and soil erosion. The opposite trend has the potential provisioning services to potential of natural conditions for recreation. However, increasing of primary and secondary production of agroecosystems must be managed with regard to the sustainability of soil multifunctionality and also the sustainability of potential of agroecosystem to provide ecosystem services in their integrity. Agroecosystems management should always be oriented to the optimization of providing of the current needs within the sustainable use of agroecosystems.
The authors acknowledge the Slovak Research and Development Agency for the financial support via Project No. APVV-0098-12 “Analysis, modelling and evaluation of agroecosystem services” amd Project No. APVV-15-0160 “Elimination of degradation processes in soil by biodiversity restoring”.
 Nieto-Romero, M., Oteros-Rozas, E., Gonzáles, J.A. and Martin-López, B. (2013) Exploring the Knowledge Landscape of Ecosystem Services Assessments in Mediterranean Agroecosystems: Insights for Future Research. Environmental Science & Policy, 37, 21-133.
 Burkhard, B., Kandziorai, M.S. and Müller, F. (2014) Ecosystem Service Potentials, Flows and Demands—Concepts for Spatial Localisation, Indication and Quantification. Official Journal of the International Association for Landscape Ecology (IALE-D), 34, 1-32.
 Frélichová, J., Vackár, D., Pártl, A., Loucková, B., Harmácková, Z. and Lorencová, E. (2014) Integrated Assessment of Ecosystem Services in the Czech Republic. Ecosystem Services, 8, 110-117.
 Birghofer, K., Diehl, E., Andersson, J., Ekroos, J., Früh-Müller, A., Machnikowski, F., Mader, V.L., Nilsson, L., Sasaki, K., Rundlof, M., Wolters, V. and Smith, H.G. (2015) Ecosystem Services—Current Challenges and Opportunities for Ecological Research. Frontiers in Ecology and Evolution, 2.
 Dominati, E., Patterson, M. and Mackay, A. (2010) A Framework for Classifying and Quantifying the Natural Capital and Ecosystem Services of Soils. Ecological Economics, 69, 1858-1868.
 Makovníková, J., Pálka, B., Siráň, M., Kanianska, R., Kizeková, M. and Jadudová, J. (2017) Modeling and Evaluation Agroecosystem Services. (Modelovanie a Hodnotenie Agroekosystémovych Sluzieb). Belianum Vydavatelstvo Univerzity Mateja Bela, Banskej Bystrici.
 (2005) Ecosystems and Human Well-Being: Our Human Planet: Summary for Decision Makers. The Millennium Ecosystem Assessment Series, MEA (Millennium Ecosystem Assessment), Vol. 5, Island Press, Washington DC.
 Burkhard, B., Kroll, F., Nedkov, S. and Müller, F. (2012) Mapping Supply, Demand and Budgets of Ecosystem Services. Ecological Indicators, 21, 17-29.
 Daniel, T.C., Muhar, A., Arnberger, A., Aznar, O. and von der Dunk, A. (2012) Contributions of Cultural Services to the Ecosystem Services Agenda. PNAS, 109, 8812-8819.
 Forouzangohar, M., Crossman, N.D., Richard, J., MacEwan, R.J., DugalWallace, O. and Bennett, L.T. (2014) Ecosystem Services in Agricultural Landscapes: A Spatially Explicit Approach to Support Sustainable Soil Management. Scientific World Journal, 2014, Article ID: 483298.
 Dominati, E.J., Mackay, A., Lynch, B., Heath, N. and Millner, I. (2014) An Ecosystem Services Approach to the Quantification of Shallow Mass Movement Erosion and the Value of Soil Conservation Practices. Ecosystem Services, 9, 204-215.
 Coyle, C., Creamer, R.E., Schulte, R.P.O., O’Sullivan, L. and Jordan, P. (2016) A Functional Land Management Conceptual Framework under Soil Drainage and Land Use Scenarios. Environmental Science & Policy, 56, 39-48.
 Haines-Young, R., Potschin, M. and Kienast, F. (2012) Indicators of Ecosystem Service Potential at European Scales: Mapping Marginal Changes and Trade-Offs. Ecological Indicartos, 21, 39-53.
 Greiner, L., Keller, A., Gret-Regamey, A. and Papritz, A. (2017) Soil Function Assessment: Review of Methods for Quantifying the Contributions of Soils to Ecosystem Services. Land Use Policy, 69, 224-237.
 Juráni, B. (2005) Protection, Land Use and Human Impact (Ochrana, vyuzívanie pody a vplyv cloveka). In: Sobocká, J., Ed., Proseedings of Soil Science Conference, Societas Pedologica Slovaca and Soil Science and Conservation Research Institute, Bratislava, 178-180. (In Slovak)
 Bujnovsky, R., Vilcek, J., Blaas, G., Skalsky, R., Barancíková, G., Makovníková, J., Balkovic, J. and Pálka, B. (2011) Assessment of Capacities of Soil and Effects from Its Use (Hodnotenie kapacít pody a efektov z jej vyuzívania). Soil Science and Conservation Research Institute, Bratislava, 70 p. (In Slovak)
 Makovníková, J., Barancíková, G. and Pálka, B. (2007) Approach to the Assessment of Transport Risk of Inorganic Pollutants Based on the Immobilisation Capability of Soil. Plant, Soil and Environment, 53, 365-373.
 Burkhard, B., Crossman, N., Nedkov, S., Petz, K. and Alkemade, R. (2013) Mapping and Modelling Ecosystem Services for Science, Policy and Practice. Ecosystem Services, 4, 1-3.
 Scolozzi, R., Morri, E. and Santolini, R. (2012) Delphi-Based Change Assessment in Ecosystem Service Values to Support Strategic Spatial Planning in Italian Landscapes. Ecological Indicators, 21, 134-144.
 Makovníková, J., Kobza, J., Pálka, B., Malis, J., Kanianska, R. and Kizeková, M. (2016) An Approach to Mapping the Potential of Cultural Agroecosystem Services. Soil and Water Research, 11, 44-52.
 Crossman, N.D., Burkhard, B., Nedkov, S., Willemen, L., Petz, K., Palomo, I., Drakou, E.G., Martín-Lopez, B., McPhearson, T., Boyanova, K., Alkemade, R., Egoh, B., Dunbar, M. and Maes, J. (2013) A Blueprint for Mapping and Modelling Ecosystem Services. Ecosystem Services, 4, 4-14.
 Naidoo, R., Balmford, A., Costanza, R., Fisher, B., Green, R.E., Lehner, B., Malcolm, T.H. and Ricketts, T.H. (2008) Global Mapping of Ecosystem Services and Conservation Priorities. Proceedings of the National Academy of Sciences of the USA, 105, 9495-9500.
 Schulp, C.J.E., Lautenbach, S. and Verburg, P.H. (2014) Quantifying and Mapping Ecosystem Services: Demand and Supply of Pollination in the European Union. Ecological Indicators, 36, 131-141.
 Bateman, I.J., Mace, G.M., Fezzi, C., Atkinson, G. and Turner, K. (2011) Economic Analysis for Ecosystem Service Assessments. Environmental and Resource Economics, 48, 177-218.
 Honigová, I., Vackár, D., Lorencová, E., Melichar, J., Gotzl, M., Sonderegger, G., Ousková, V., Hosek, M. and Chobot, K. (2012) Survey on Grassland Ecosystem Services, Report to the European Topic Centre on Biological Diversity. Nature Conservation Agency of the Czech Republic, Prague, 78.
 Frélichová, J. and Fanta, J. (2015) Ecosystem Service Availability in View of Long-Term Land-Use Changes: A Regional Case Study in the Czech Republic. Ecosystem Health and Sustainabilty, 1, 1-15.
 Maes, J., Braat, L., Jax, K., Hutchins, M., et al. (2011) A Spatial Assessment of Ecosystem Services in Europe: Methods, Case Studies and Policy Analysis—Phase 1. PEER Report No. 3, Partnership for European Environmental Research, Ispra.
 Hamnett, R.G. (2006) Country Pasture/Forage Resources Profile-Slovakia.
 NAFC-SSCRI (2015) Distribution of Agricultural Land in Agro-Climatic Regions of Slovakia.
 Kobza, J., Barancíková, G., Dodok, R., Hrivňáková, K., Makovníková, J., Pálka, B., Pavlenda, P., Schlosserová, J., Styk, J. and Siráň, M. (2014) Soil Monitoring of the Slovak Republic (Monitoring pod SR). National Agriculture and Food Centre-Soil Science and Conservation Research Institute, Bratislava, 252. (In Slovak)
 Van Camp, B., Bujarrabal, A.R., Gentile, R.J.A., Jones, L., Montanarella, L., Olazabal, O. and Selvaradjou, S.K. (2004) Reports of the Technical Working Groups Established under the Thematic Strategy for Soil Protection. EUR 21319 EN/5, 872. Office for Official Publications of the European Communities, Luxembourg.
 Alam, M., Dupra, J. and Messier, C.H. (2016) A Framework towards a Composite Indicator for Urban Ecosystem Services. Ecological Indicators, 60, 38-44.
 Bujnovsky, R., Balkovic, J., Barancíková, G., Makovníková, J. and Vilcek, J. (2009) Assessment and Valuation of Ecological Functions of Agricultural Land in Slovakia. Soil Science and Conservation Research Institute, Bratislava, 72. (In Slovak)
 Majtaníková, J. (2011) Vyhodnotenie eróznej ohrozenosti vybraného podneho celku na území povodia Hostovského potoka. In: Krajinné inzinierstvo-súcasny stav a vyhlad do budúcnosti, Slovenská polnohospodárska univerzita, Nitra.
 Barancíková, G., Gutteková, M., Halas, J., Koco, S., Makovníková, J., Nováková, M., Skalsky, R., Tarasovicová, Z. and Vilcek, J. (2011) Soil Organic Carbon in the Agricultural Landscape-Modeling of Changes in Space and Time. Soil Science and Conservation Research Institute, Bratislava, 85 p. (In Slovak)
 Sherrouse, B.C., Clement, J.M. and Semmens, D.J. (2011) A GIS Application for Assessing, Mapping, and Quantifying the Social Values of Ecosystem Services. Applied Geography, 31, 748-760.
 Krkoska Lorencová, E., Harmácková, Z.V., Landová, L., Pártl, A. and Vackár, D. (2016) Assessing Impact of Land Use and Climate Change on Regulating Ecosystem Services in the Czech Republic. Ecosystem Health and Sustainability, 2, e01210.
 Kanianska, R., Jadudová, J., Makovníková, J. and Kizeková, M. (2016) Assessment of Relationships between Earthworms and Soil Abiotic and Biotic Factors as a Tool in Sustainable Agricultural. Sustainability, 8, 906.
 Montoya, J.M. and Raffaelli, D. (2010) Climate Change, Biotic Interactions and Ecosystem Services. Philosophical Transactions of the Royal Society B, 365, 2013-2018.
 Birkhofer, K. and Wolters, V. (2012) The Global Relationship between Climate, Net Primary Production and the Diet of Spiders. Global Ecology and Biogeography, 21, 100-108.
 Diehl, E., Sereda, E., Wolters, V. and Birkhofer, K. (2013) Effects of Predator Specialization, Host Plant and Climate on Biological Control of Aphids by Natural Enemies: A Meta-Analysis. Journal of Applied Ecology, 50, 262-270.
 Spake, R., et al. (2017) Unpacking Ecosystem Service Bundles: Towards Predictive Mapping of Synergies and Trade-Offs between Ecosystem Services. Global Environmental Change, 47, 37-50.
 Barancíková, G. and Makovníková, J. (2003) The Influence of Soil Humic Acid Quality on Sorption and Mobility of Heavy Metals. Plant, Soil and Environment, 49, 565-571.
 Siráň, M., Makovníková, J. and Barancíková, G. (2013) Monitoring of Soil Bulk Demsity-Basis for Determination of Organic Matter Stock in Soil. Proseedings of Soil Science and Conservation Research Institute, Bratislava, 166-173. (In Slovak)
 Campbell, C.A. and Souster, W. (1982) Loss of Organic Matter and Potentially Mineralizable Nitrogen from Saskatchewan Soils Due to Cropping. Canadian Journal of Soil Science, 62, 651-656.
 Schnitzer, M., Mcarthur, D.F.E., Schulten, H.R., Kozak, L.M. and Huang, P.M. (2006) Long-Term Cultivation Effects on the Quality and Quantity of Organic Mater in Selected Canadian Prairie Soils. Geoderma, 130, 141-156.
 Conant, R.T., Ryan, M.G., Agren, G.I., Birge, H.E., Davidson, E.A., Eliasson, P.E., Evans, S.E., Frey, S.D., Giardina, C.P., Hopkins, F.M., Hyvonen, R., Kirschbaum, M.U.F., Lavallee, J.M., Leifeld, J., Parton, W.J., Steinweg, J.M., Wallenstein, M.D., Wetterstedt, J.A.M. and Bradford, M.A. (2011) Temperature and Soil Organic Matter Decomposition Rates-Synthesis of Current Knowledge and a Way Forward. Global Change Biology, 17, 3392-3404.
 Martens, S.N., Breshears, D.D. and Meyer, C.W. (2000) Spatial Distributions of Understory Light along the Grassland/Forest Continuum: Effects of Cover, Height, and Spatial Pattern of Tree Canopies. Ecological Modeling, 126, 79-93.
 Felipe-Lucia, M.R., Comin, F.A. and Bennett, E.M. (2014) Interactions among Ecosystem Services across Land Uses in a Floodplain Agroecosystem. Ecology and Society, 19, Art. 20.
 Mahbubul, A., Duprasb, J. and Messier, C.H. (2016) A Framework towards a Composite Indicator for Urban Ecosystem Services. Ecological Indicators, 60, 38-44.
 Mouchet, M.A., Paracchini, M.L., Schulp, C.J.E., Stürck, J., Verkerk, P.J., Verburg, P.H. and Lavorel, S. (2017) Bundles of Ecosystem (Dis-)Services and Multifunctionality across European Landscapes. Ecological Indicators, 7, 23-28.
 Lescourret, F., Magda, D., Richard, G., Blondon, A.F.A., Bardy, M., Baudry, J., Doussan, I., Dumont, B., Lefèvre, F., Litrico, I., Clouaire, R.M., Montuelle, B., Pellerin, S., Plantegenest, M., Tancoigne, E., Thomas, A., Guyomard, H. and Soussana, L.F. (2015) A Social-Ecological Approach to Managing Multiple Agro-Ecosystem Services. Environmental Sustainability, 14, 68-75.