Next Article in Journal
Assessing the Nonlinear Relationship between Land Cover Change and PM10 Concentration Change in China
Previous Article in Journal
Conservation Management Practices for Biodiversity Preservation in Urban Informal Green Spaces: Lessons from Central European City
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ecosystem Services Supply from Peri-Urban Watersheds in Greece: Soil Conservation and Water Retention

by
Stefanos Stefanidis
1,*,
Nikolaos Proutsos
2,
Vasileios Alexandridis
3 and
Giorgos Mallinis
3
1
Forest Research Institute, Hellenic Agricultural Organization “DIMITRA”, Vassilika, 57006 Thessaloniki, Greece
2
Institute of Mediterranean Forest Ecosystems, Hellenic Agricultural Organization “DIMITRA”, Terma Alkmanos, 11528 Athens, Greece
3
Laboratory of Photogrammetry and Remote Sensing (PERS Lab), School of Rural and Surveying Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Land 2024, 13(6), 765; https://doi.org/10.3390/land13060765
Submission received: 30 April 2024 / Revised: 27 May 2024 / Accepted: 28 May 2024 / Published: 29 May 2024

Abstract

:
The main objective of this study is to investigate the supply of ecosystem services (ESs) within peri-urban watersheds (PUWs) across Greece, utilizing earth observation (EO) data and empirical models. To achieve these goals, several open-access geospatial datasets were utilized and processed within a GIS environment. Specifically, the supply of soil conservation and water retention services was assessed at the watershed scale. The results indicates that the supply soil conservation service withing the examined PUWs range from 44.41 t ha−1 y−1 to 441.33 t ha−1 y−1 with an average value of 161.99 t ha−1 y−1. Water retention services vary from 35.09 to 154.63 mm within the PUWs, with an average of 91.45 mm. Additionally, the variation in ES values across distinct altitudinal zones and ecosystem types provides useful insights for implementing protection and management measures. It is evident from the analysis that the soil conservation service gradually increases with elevation, with alpine and subalpine areas exhibiting the highest values. Moreover, mountainous and semi-mountainous regions demonstrate higher values compared to the estimated average soil conservation services of the entire study area. Furthermore, the provision of water retention services varies across different altitudinal zones. Specifically, the mountainous and semi-mountainous areas exhibit the highest values, followed by the subalpine and hilly regions, while lower values are observed in the lowland and alpine regions. The analysis also reveals distinct patterns in ecosystem services across various ecosystem types within the PUWs. Woodland and forest, along with heathland and shrubs, demonstrate higher values in terms of both soil conservation and water retention services. Our findings contribute to understanding the dynamics of ESs in PUWs, highlighting their significance for sustainable land management, and informing policy decisions aimed at preserving ecosystem health and resilience.

1. Introduction

Since 2008, more than half of the global population has been living in urban regions, and this figure is expected to rise to 68% by 2050 according to average projections [1]. This unprecedented urban development witnessed over the last century and the forthcoming urbanization prospects in both developed and developing countries represent a phenomenon with profound economic, social, and environmental implications at regional to global scales [2]. While urbanization can positively influence economic development and productivity, the associated functional, morphological, and structural changes often impose increased pressure on adjacent natural resources, leading to adverse effects on environmental quality, ecosystem functions, and sustainability [2,3].
At the landscape level, the primary challenge lies in identifying the most effective approach to managing various land cover options [4]. The growth of urban regions leads to substantial changes in landscapes, transforming natural ecosystems into impermeable urban zones at an accelerated rate [5,6]. Despite the growing recognition of the environmental impacts of such transformations, unsustainable land-use practices persist. Land management decisions, influenced by policy and environmental planning, play a pivotal role in shaping land cover, and directly impact ecosystem services (ES) [7,8]. The Millennium Ecosystem Assessment (MEA), a major UN-sponsored effort to analyze the impact of human actions on ecosystems and human well-being, identified four major categories of ecosystem services: provisioning, regulating, cultural services, and supporting services [9].
Peri-urban landscapes, characterized by their adjacency to urban centers and the blending of urban and rural areas, face heightened vulnerability to the pressures of urban expansion and land-use alterations. The conversion of peri-urban regions from natural to anthropogenic landscapes not only compromises their ecological integrity but also undermines the provision of essential ESs crucial for human well-being and environmental sustainability. For a considerable period, the lack of awareness regarding the significance of ESs has led to numerous challenges, resulting in the excessive exploitation of natural resources, and affecting ecological security on the regional and global scales [10]. As urban sprawl encroaches upon these areas, the demand for green infrastructure, including peri-urban forests and open green spaces, intensifies.
Green infrastructures play a vital role in mitigating the adverse effects of urbanization by providing habitats for biodiversity, regulating microclimates, and enhancing the overall quality of life for residents [11,12,13]. Their significance spans from local to global scales, contributing substantially to the achievement of numerous Sustainable Development Goals outlined by the UN for 2030. They also offer a spectrum of material and non-material benefits to urban and peri-urban dwellers, encompassing carbon sequestration and climate regulation [14], atmospheric pollutant removal and air quality improvement [15,16], water flow regulation [17], purification [18], provision of wood and non-wood products [14], recreational and cultural opportunities [19], educational and science opportunities [20], renewable energy production [21], disaster risk reduction [22], and soil conservation [23]. Particularly noteworthy is their role in promoting physical and mental well-being, as evidenced during crises such as the COVID-19 pandemic [24,25].
Among these services, soil conservation and water retention capacity are of paramount importance, especially in the Mediterranean basin. This area is acknowledged as a focal point for both climate change and biodiversity [26,27], and it is anticipated to encounter heightened challenges in water-related services due to the climate crisis. The Food and Agricultural Organization (FAO) of the United Nations (UN) developed a comprehensive and practical international Forests and Water Agenda to guide future action [28]. While ESs are universally important, across all landscapes and at different scales, they are especially crucial in mountainous watersheds [29]. Moreover, water-related services should be studied at the watershed scale, following the trend of the EU directives (e.g., 2007/60/EC, 2000/60/EC) for integrated water resources management. Watershed areas are crucial units for studying ES dynamics due to their unique ecological, economic, and geographic characteristics. They provide a nuanced understanding of the supply and demand dynamics of ecosystem services, characterized by intricate interactions between natural processes and human activities. Scholars address these challenges by evaluating ecosystem services using watershed-scale approaches [30,31,32]. However, most studies are conducted as case studies rather than on a broader scale. In recent years, the expanding accessibility of earth observation (EO) data coupled with the widespread utilization of geospatial technologies enables the extensive evaluation of ESs on a large scale [33].
This study aims to evaluate the supply of soil conservation and water retention capacity ES in peri-urban watersheds (PUWs) across Greece using EO data. More specifically, this study provides insight into the following: (i) hotspot identification of the upstream PUWs for ES provision and (ii) exploration of the pattern between altitudinal zones and ecosystem types with ES values. By delving into these complex dynamics within watershed contexts, the research contributes to a deeper understanding of sustainable land management practices and conservation strategies in peri-urban environments.

2. Materials and Methods

2.1. Study Area

The study was conduct in selected PUWs, located upstream of major cities across Greece. The location of these Greek cities can be found in Figure 1, while the characteristics of the PUWs are detailed in Table A1 (Appendix A). Specifically, the study examined 54 watersheds upstream of 14 cities.
To define the study area, a set of open-access geospatial data were utilized. Spatial boundaries of cities were derived from data provided by the Organisation for Economic Co-operation and Development (OECD), accessible at https://www.oecd.org/regional/regional-statistics/functional-urban-areas.htm (accessed on 10 January 2024). The OECD, in cooperation with the EU, has developed a harmonized definition of cities. According to this definition, a city is a grouping of local administrative units (LAUs) in European countries, such as municipalities or local authorities, where at least 50% of the population resides in an urban center. An urban area is defined by a grouping of adjacent grid cells covering a single square kilometer, with a population density of at least 1500 residents per square kilometer and a total population of no less than 50,000 individuals. To this end, fourteen cities were identified in Greece. Regarding the population data from the aforementioned source, two out of fourteen cities are characterized as metropolitan (with population more than 1 million inhabitants) named Athens and Thessaloniki, while the rest of the cities are rather smaller, and their population range from 46,800 to 156,000 inhabitants. The population of each city can be seen in detail in Table A1.
In terms of watershed, the thematic layers were retrieved from the Hellenic Ministry of the Environment and Energy (MEEN), accessible at https://geodata.gov.gr/en/dataset/oreines-lekanes-aporroes-2es-taxes (accessed on 10 January 2024). This dataset includes the mountainous watersheds of the Greek territory, as described in the torrent registry system compiled for the entire territory and used from local authorities as a basis for watershed management. Subsequently, the watersheds that intersect the borders of the cities were identified, thus determining PUWs that formed the study area of the present research (Table A1). The area of the PUWs included in this study ranged from 6.38 km2 to 406.51 km2.

2.2. Soil Conservation Service

The assessment of ecosystem services related to soil conservation relied upon the widely recognized empirical model of the Universal Soil Loss Equation (USLE) combined with some features of the updated model’s revised form [34]. To accomplish this, potential and actual erosion rates were calculated using Equations (1) and (2), respectively. Furthermore, the soil conservation service, representing the disparity between potential and actual soil erosion rates was computed as detailed in Equation (3).
Potential soil erosion rate:
E p = R × K × L S
Actual soil erosion rate:
E a = R × K × L S × C × P
Soil conservation service:
E c = E P E a
where Ep and Ea represent the potential and actual rates of soil erosion, respectively (t ha−1 y−1), while Ec denotes the soil conservation service (t ha−1 y−1). The variable R stands for the rainfall erosivity factor (MJ mm ha h−1 y−1), K denotes the soil’s erodibility factor (t ha h ha−1 MJ−1 mm−1), LS signifies the topographic factor (dimensionless), which is determined by slope length (L) and slope steepness (S), C indicates the cover management factor (dimensionless), and P represents the support practice factor (dimensionless).
The R factor, which serves as the model’s climate component, was obtained from a newly created dataset by leveraging the Rainfall Erosivity Database at European Scale (REDES) in conjunction with UERRA regional reanalysis rainfall data spanning the period 1990–2018 [35]. This dataset, furnished by the European Soil Data Center (ESDAC) with spatial resolution 1 km, is accessible at https://esdac.jrc.ec.europa.eu/content/rainfall-erosivity-european-union-and-switzerland (accessed on 10 January 2024).
The soil erodibility (K-factor) represents the soil’s resistance to the erosive effects of rainfall, influenced by various soil properties such as physical, chemical, and mineralogical characteristics, as well as their interactions. In this research, the soil erodibility factor is determined by employing following equation provided by Wischmeier and Smith [36]:
Κ = [ 2.1 × 10 4 Μ 1.14 12 O Μ + 3.25 s 2 + 2.5 p 3   100 × 0.1317 ]
where M represents the grain size parameter, OM denotes organic matter content (%), s represents soil structure, and p signifies permeability. The calculation of M relies on the on the percentages of the “silt + very fine sand” fraction (Ps%) and clay fraction (Pc%) as depicted in the equation below:
M = P S ( 100 P C )
In this study, soil properties used in the calculations were obtained from SoilGrids, a state-of-the-art digital soil mapping system developed by the International Soil Reference and Information Centre (ISRIC—World Soil Information), accessible at https://soilgrids.org/ (accessed on 10 January 2024). SoilGrids employs advanced machine learning techniques along with a vast dataset comprising soil profile data (WoSIS) and environmental layers to provide precise, high-resolution soil property information at a resolution of 250 m [37]. This comprehensive resource enables the accurate characterization of soil attributes, including fractions of sand, silt, and clay, as well as organic matter content and descriptors related to soil structure and permeability.
Τhe topographic factor (LS) is composed of individual factors, slope length (L) and steepness (S), which collectively describe their combined effect on soil loss. A crucial parameter for calculating the LS-factor is the digital elevation model (DEM). In this study, we utilized the newly developed Forest And Buildings removed DEM (FABDEM). FABDEM is a global elevation map that eliminates biases from building and tree heights in the Copernicus GLO-30 Digital Elevation Model (DEM) using machine learning techniques [38]. FABDEM emerges as the preferred choice due to its superior accuracy, particularly in areas with dense vegetation, presenting a significant enhancement over existing open-access DEM products with similar resolutions (30 m) [39].
The calculation of the LS-factor was conducted using advanced algorithms from the System for Automated Geoscientific Analyses (SAGA) within the latest version of QGIS (v3.22) software. This software package integrates the multiple flow algorithm, enabling comprehensive terrain analysis [40].
The formula introduced by McCool et al. [41] was utilized to derive the S-factor, relying on the slope gradient (ϑ) expressed in degrees:
S = 10.8 × s i n ϑ + 0.03 , ϑ < 0.09 16.8 × s i n ϑ 0.5 , ϑ > 0.09
Furthermore, the L-factor was determined utilizing the recommended equation by Desmet and Govers [42], which has demonstrated reliability in regions characterized by intricate terrain [40]. This approach, incorporating the concept of unit-contributing area to account for variations in slope steepness across the area, is detailed by Equation (7):
L = A i , j + D 2 m + 1 A i , j i n m + 1 D m + 2 × x i , j m × 22.13 m
In this method, Ai,j,–in represents the contributing area (m2) at the inlet of the grid pixel denoted by (i,j). D stands for the size of the grid pixel in meters. The variable xi,j denotes the combined sine and cosine values of the aspect direction (αi,j) of the grid pixel, calculated as xi,j = sin αi,j + cos αi,j. The coefficient m is associated with the ratio β of rill to inter-rill erosion. Its values range from 0 to 1, while ϑ represents the slope angle in degrees. The m coefficient is calculated as follows:
m = β β + 1
β = s i n ϑ 0.0896 0.56 + 3 × s i n ϑ 0.8
The cover management factor (C) delineates the impact of land cover on erosion. Herein, the estimation of the C-factor relied on a national-scale ecosystem type product developed as a part of the LIFE-IP 4 Natura project [43], based on land use/cover (LULC) classification. This dataset, specifically tailored for ecosystem services mapping and assessment, served as the basis for our analysis and was supplied by the Natural Environment and Climate Change Agency (N.E.C.C.A.), the primary governmental body overseeing the National System of Governance of Protected Areas (PAs) in Greece. To determine the C-factor, appropriate values were assigned to each LULC based on the literature [44,45]).
Finally, the impact of support practices (P-factor), such as contour farming, stone walls, terracing, and grass margins on erosion mitigation, was evaluated using data retrieved from the European Soil Data Centre (ESDAC), accessible at https://esdac.jrc.ec.europa.eu/content/support-practices-factor-p-factor-eu (accessed on 10 January 2024), providing comprehensive coverage for the entire Europe [46].
The raster calculator tool was employed to combine various factors, resulting in the derivation of soil conservation service values, with the subsequent exportation of statistics at the watershed level per city. Additionally, zonal statistics were conducted to illustrate the outcomes across distinct altitudinal zones and ecosystem types.

2.3. Water Retention Service

The water storage capacity plays a pivotal role in hydrology, influencing both peak discharge and water supply. The assessment of water retention ecosystem service was based on the Soil Conservation Service-Curve Number (SCS-CN) approach [47], initially formulated by the SCS under the U.S. Department of Agriculture (now referred to as the Natural Resources Conservation Service—NRCS). According to this approach, the maximum water retention capacity (S) is directly influenced by the runoff coefficient (CN), as expressed by the following equation:
S = 25400 C N 254
In the SCS-CN method, there is an inverse relationship between the maximum water retention capacity and the runoff coefficient. Soils with higher retention capacities exhibit lower runoff coefficients, while those with lower capacities generate runoff more rapidly.
To calculate this service, we estimated the maximum water retention capacity based on the current LULC (SC) and compared it to a scenario without any vegetation (SnoVeg). This scenario typically excludes vegetation parameters from the Curve Number (CN) estimation, following the methodology of the erosion conservation service. Therefore, the water retention ES was defined as the difference between these two states. The empirical equation developed in the framework of the Deucalion Project [48] was used to calculate the CN for average humidity conditions. The updated evaluation of the CN parameter has improved the accuracy and spatial representation of CN across the entire watershed [49].
The CN value is determined by considering the predominant factors such as water permeability, land cover, and drainage capacity of the watershed, utilizing the following empirical equations [48,49]:
Current CN:
C N C = 10 + 9 × i P E R M + 6 × i V E G + 3 × i S L O P E
CN without vegetation:
C N n o V e g = 10 + 9 × i P E R M + 30 + 3 × i S L O P E
where iPERM represents water permeability, iVEG denotes vegetation density, and iSLOPE indicates drainage capability, each variable ranging from 1 to 5 based on literature [48,49].
To determine the iPERM values, the parent material of the studied Permeable Urban Watersheds (PUWs) was identified using the national scale soil map of Greece. The vector data of the parent material are accessible through the Hellenic Ministry of the Environment and Energy (MEEN) website at http://mapsportal.ypen.gr/maps/289 (accessed on 10 January 2024). Subsequently, each parent material was assigned an appropriate value based on the relevant literature [49]. A ranking system ranging from 1 to 5 was then employed, where an index of one corresponds to substrata with very high permeability (e.g., karst formations), while an index of five denotes substrata with very low permeability (e.g., dense rocks).
Vegetation classes are delineated based on land characteristics such as retention mechanisms, soil roughness, and filtration capacity, influenced notably by root zone growth. Leveraging the LULC product and the associated ecosystem types established within the LIFE-IP 4 Natura project, as described above [45], the vegetation type within the designated area is first classified as dense, moderate, undergrowth, sparse, or absent. Then, a grading scale from 1 to 5 is implemented, with a score of 1 indicating regions with dense vegetation (e.g., evergreen forests), while an index of 5 indicates bare soil [49].
The initial classification of the drainage capability encompasses categories such as negligible, low, moderate, high, and very high, followed by the assignment of ranks ranging from 1 to 5. These rankings are determined solely by considering five terrain slope categories, as this information can be readily obtained through standard DEM processing. A rank of one is designated for nearly flat areas, while a rank of five corresponds to slopes exceeding 30% [49].
Similarly with the soil conservation, water retention service values were derived using raster calculator tool and then exported as statistics at the watershed level per city. Additionally, zonal statistics were conducted to visualize the outcomes across different altitudinal zones and ecosystem types.

2.4. Overall Workflow

Various geospatial datasets were acquired and processed to fulfill the objectives of this study, including topographic, climatic, land cover, pedological, and geomorphological data. With the aid of the open-source Quantum GIS (QGIS) software, the datasets were structured into GIS thematic layers. The study area was delineated, and factors influencing soil conservation and water retention were identified. The output pixel size for ecosystem service (ES) analysis was set to 100 m, with the coordinate system set as the standard European Coordinate Reference System, utilizing the European Terrestrial Reference System 1989 (ETRS89) datum and the Lambert Azimuthal Equal Area (LAEA) projection (EPSG: 3035). The selection of a 100 m pixel scale was considered optimal, allowing for the adjustment of the C-factor layer, which operates at a similar resolution, in response to policy interventions affecting land use, as highlighted in previous large-scale assessments [50]. In summary, the methodology’s entire workflow is depicted in the accompanying figure (Figure 2).

3. Results

The results demonstrate the critical role of PUWs in providing ESs that contribute to sustainable development and human welfare. The average soil conservation supply in the PUWs located upstream of the examined cities in Greece was found to be 165.54 t ha−1 y−1 with a range from 44.41 t ha−1 y−1 in Dendropotamos watershed (Thessaloniki) to 441.33 t ha−1 y−1 in Nedonas watershed (Kalamata). Figure 3 depicts the average soil conservation service along with the standard deviation (SD) yielded by the PUWs for the studied cities.
The analysis identified the upstream PUWs of Kalamata (309.03 t ha−1 y−1, SD: 146,37 t ha−1 y−1), Ioannina (260.20 t ha−1 y−1, SD: 122.42 t ha−1 y−1), and Xanthi (245.84 t ha−1 y−1, SD: 58.53 t ha−1 y−1) as the top hotspots for soil conservation supply, while the PUWs upstream of Chania, Athens, and Thessaloniki were found to yield the lowest soil conservation services reaching 99.35 (62.00), 84.59 (SD: 34.75), and 64.19 (SD: 18.05) t ha−1 y−1, respectively. The spatial distributions of the soil conservation service in selected PUWs are presented in the following figure (Figure 4).
Moreover, the water retention service of the PUWs was assessed as an additional regulating ecosystem service within the PUWs, revealing a range from 35.09 mm in Vrilision watershed (Athens) to 155.48 mm in Xirila watershed (Kalamata), with an average of 91.45 mm. In detail, the average water retention supply along with the SD values of the PUWs upstream of the cities under study is illustrated in Figure 5, while their spatial distribution can be seen in Figure 6.
The PUWs demonstrated varying levels of water retention service across the examined cities. The higher amount of water retention service was provided by the PUWs over Kalamata (120.27 mm, SD: 40.35 mm), Trikala (112.3 mm), and Ioannina (103.23 mm, SD: 20.70 mm). On the contrary, the lowest values were found in the PWUs over Chania (82.1 mm, SD: 7.28 mm), Thessaloniki (76.41 mm, SD: 21.24 mm), and Irakleio (62.38 mm, SD: 11.93 mm). Figure 5 provides the spatial distribution of the water retention service within the selected PUWs.
Subsequently, an investigation was carried out to discern the diverse values of ecosystem services across distinct altitudinal zones and ecosystem types. The analysis of the altitudinal zones revealed a significant dominance of lowland areas (0–200 m) and hilly areas (201–500 m) within the PUWs, accounting for 34.3% and 32.4% of the total area, respectively. As the elevation increased, there was a gradual decrease in coverage. Semi-mountainous areas (501–1000 m) represented 21.8% of the total area, while mountainous areas (1001–1500 m) comprised 9.2%. In contrast, subalpine and alpine areas presented relatively small percentages (<2%) in the overall landscape. The average soil conservation and water retention capacity services in each altitudinal are presented in Figure 7.
Based on the above graph, it is evident that the soil conservation service gradually increases as elevation increases. The alpine and subalpine areas exhibit the highest values, while the mountainous and semi-mountainous regions demonstrate higher values in contrast to the calculated mean soil conservation services across the entirety of the study region. Additionally, the provision of water retention services varies across different altitudinal zones. Specifically, the mountainous and semi-mountainous areas exhibit the highest values, followed by the subalpine and hilly regions. Conversely, lower values are observed in the lowland and alpine regions.
Afterwards, the variation in ecosystem service supply between different terrestrial ecosystem types were examined based on the classification scheme for the ecosystem type mapping and assessment (MAES) level 2 [51]. Croplands, heathland and shrubs, and woodland and forest occupy roughly equal percentages (about 30%) in the research area, while grasslands and sparsely vegetated areas occupy rather smaller percentages, 5.2% and 1.9%, respectively. Also, the percentage coverage of wetlands is negligible (0.11%). The graphical representation of the average soil conservation and water retention services within each MAES category is depicted in the figure below (Figure 8).
The analysis of the provided data reveals distinct patterns in ecosystem services across various ecosystem types within the PUWs. Woodland and forest, along with heathland and shrubs, demonstrate higher values in terms of both soil conservation and water retention services. In contrast, croplands and wetlands exhibit lower values in these services. It is worth mentioning that grasslands and areas with sparse vegetation surpass the anticipated average soil conservation services of the study region.

4. Discussion

The ES concept was developed with the intention of integrating it into decisionmakers’ agendas, employing various methods focused on environmental and natural resource management. However, its limited incorporation thus far can be attributed to the absence of institutional capacity and the credibility of the generated information [52]. Due to the ongoing increase in the availability of EO data, the established utilization of geospatial technologies, and the use of empirical models, conducting the large-scale assessments of ES has become feasible [33,53]. In pursuit of this objective, national-scale approaches have been employed to evaluate the provision of numerous ecosystem services for conservation planning [51,54].
In recent decades, demographic trends and urban development have intensified the pressures on ecosystems near city boundaries, thereby presenting challenges in estimating ESs for both urban and peri-urban landscapes. This direct linkage between ecosystem services and the advancement in implementing the United Nations Sustainable Development Goals (SDGs) underscores the significance of comprehending and managing these dynamic relationships [55,56]. Nevertheless, limited studies have been carried out on the scale of peri-urban watersheds. As a part of such an approach, An et al. [33] evaluated temporal changes in soil conservation services across four large basins worldwide, providing a valuable conceptual framework.
This research utilized the SCS-CN and RUSLE modeling approaches to estimate the ESs provided by the Greek PUWs. Both models have been widely applied, especially in ungauged watersheds. However, their uncertainties need to be acknowledged to contextualize the findings appropriately. For instance, both methods rely on parameters like land cover, topography, and soil properties, which are derived from Earth Observation (EO) data. However, EO data come with inherent limitations in spatial resolution and temporal coverage, which can impact the accuracy of soil conservation and water retention estimates. Moreover, these parameters might not capture site-specific conditions accurately. In particular, the largest uncertainties in the SCS-SN model, regarding maximum water capacity, arise from the lack of knowledge about the previous soil moisture conditions [57]. On the other hand, when it comes to erosion estimation, the uncertainties in the RUSLE model stem from the fact that critical processes, such as sediment transport, deposition, and routing within hydrographic networks, are not considered [58]. Nonetheless, comparing erosion rates and maximum water retention capacity between the current LULC and a scenario without vegetation effectively quantifies the supply of the related ESs.
Our analysis was performed in the upstream watersheds of the largest Greek cities. It delves into the pivotal role of PUWs in Greece as key contributors to essential ecosystem services, focusing on soil conservation and water retention capacity. These services play a critical role in fostering sustainable urban development and enhancing human well-being. The spatial distribution of soil conservation and water retention services highlights distinct hotspots and areas with lower contributions. The findings reveal the varying degrees of soil conservation and water retention services provided by different PUWs across Greek cities. Soil conservation services, crucial for preventing soil erosion and sedimentation, exhibit substantial spatial disparities. Noteworthy is the identification of Kalamata, Ioannina, and Xanthi as hotspots for robust soil conservation supply, offering potential models for effective soil preservation strategies. Conversely, regions like Athens, Thessaloniki, and Chania face challenges in soil conservation, underlining the complexities associated with urbanized landscapes. Water retention capacity, crucial for urban resilience against flooding and water scarcity, also varies across PUWs. Kalamata, Trikala, and Ioannina emerge as leaders in this regard, indicating their potential to manage stormwater runoff and mitigate flood risks. On the other hand, areas like Chania, Thessaloniki, and Irakleio demonstrate lower water retention services, underscoring their vulnerability to water-related urban challenges.
Altitudinal analysis elucidates the relationship between elevation and ecosystem services, revealing how higher elevations tend to enhance soil conservation and water retention capacities. This emphasizes the importance of alpine and subalpine areas as valuable contributors to erosion control and flood management. The relationship between average soil conservation services and elevation has been evaluated in similar studies. It was found that these services increase with elevation up to 1000 m, but decrease for higher altitudes, despite the non-linear behavior observed in the latter case [33]. The investigation into different terrestrial ecosystem types reveals significant variations in service supply. Woodland and forest, along with heathland and shrubs, emerge as prominent contributors to both soil conservation and water retention. These findings align with established ecological principles, as vegetation cover and root systems in forests and shrublands play a pivotal role in stabilizing soils and regulating water flow [59]. In contrast, croplands and wetlands display relatively lower levels of these ecosystem services, necessitating context-specific management interventions to enhance their contributions. The study accentuates the need for integrating green infrastructure and sustainable land management practices into urban planning to optimize ecosystem services.
The insights gained from this study have far-reaching implications for urban planning and policy development. By recognizing and leveraging the varying capacities of different PUWs, cities can formulate site-based strategies in an effort to improve ES supply. Despite the importance of the soil conservation and water retention ESs, they are closely linked to other regulating services, such as carbon sequestration and climate regulation. By preventing soil erosion, these services help maintain soil fertility and structure, which in turn supports plant growth and carbon storage. This interconnection highlights the role of PUWs in mitigating climate change and enhancing ecosystem resilience. Moreover, the preservation of soil and water resources contributes to provisioning services by sustaining agricultural productivity and ensuring a stable water supply. Healthy soils and efficient water retention are essential for crop growth, which directly supports food security and local economies. These services also reduce the need for artificial fertilizers and irrigation, leading to more sustainable agricultural practices.
Incorporating these findings into urban development plans can lead to more sustainable land use, improved disaster risk reduction, and enhanced urban livability. While this study provides a comprehensive analysis of soil conservation and water retention capacity in Greek PUWs, there are avenues for further research. Studying the economic value of these ESs can offer a deeper insight into their tangible benefits for society. Moreover, exploring the potential synergies and trade-offs among different ESs within PUWs could provide valuable guidance for making more holistic land-use decisions.

5. Conclusions

In conclusion, this study underscores the critical importance of PUWs as essential suppliers of ESs that contribute to sustainable urban development and human well-being. The analysis of soil conservation and water retention capacity across diverse landscapes and urban contexts offers valuable insights into the potential of these natural systems to enhance urban resilience.
By recognizing the spatial disparities in ES supply, cities can strategically allocate resources and prioritize conservation efforts. The identification of hotspots and areas with lower contributions provides a roadmap for targeted interventions that can maximize the benefits of ecosystem services for both current and future generations.
The integration of ecosystem service considerations into urban planning and policy can pave the way for more resilient and livable cities. Embracing the principles of sustainable land management, green infrastructure, and nature-based solutions can facilitate a harmonious coexistence between urbanization and ecological functioning. As cities continue to grapple with the challenges of rapid urbanization, the findings of this study offer valuable guidance for decisionmakers, planners, and researchers seeking to create urban environments that are both sustainable and supportive of human well-being. By harnessing the potential of peri-urban watersheds and their ecosystem services, we can build cities that thrive in the face of environmental uncertainties while fostering a healthier and more vibrant urban future. Hence, there is a necessity for a unified directive that addresses the management of watersheds and the services they offer, particularly in areas adjacent to urban and peri-urban zones.

Author Contributions

Conceptualization, G.M.; methodology, S.S. and G.M.; software, S.S. and V.A.; formal analysis, S.S. and N.P.; investigation, S.S.; data curation, S.S.; writing—original draft preparation, S.S. and N.P.; writing—review and editing, S.S., N.P. and G.M.; visualization, S.S. and N.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the European Commission LIFE Integrated Project, LIFE-IP 4 NATURA “Integrated Actions for the Conservation and Management of Natura (2000) sites, species, habitats and ecosystems in Greece”, Grant Number: LIFE 16 IPE/GR/000002, and the Green Fund, Ministry of Environment and Energy (Greece).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Upstream watersheds for the fourteen cities under study.
Table A1. Upstream watersheds for the fourteen cities under study.
IdCityPopulationWatershed NameWatershed Area (km2)
1Athens3,608,000Asopou150.56
2Filis113.08
3Group of Gullies East Hymettus197.35
4Group of Gullies West Hymettus82.31
5Kifissou406.51
6Lauriou251.69
7Martathonos293.63
8Rafinas161.74
9Prelesion71.56
10Chania62,040Chanion91.90
11Soudas156.44
12Ioannina64,500Kranoulas67.90
13Marmaron58.08
14Irakleio138,600Gazanou294.79
15Giofiriou175.12
16Karterou230.66
17Kalamata52,500Eleochoriou54.53
18Nedona132.18
19Xirila114.54
20Katerini57,200Gerakari55.48
21Mavroneriou613.29
22Tsiamanti45.75
23Kavala46,800Chalkerou56.47
24Group of torrents between Iraklitsa-Kavala54.19
25Kokkinoxomatos22.32
26Larisa144,000Sikouriou99.90
27Patras156,000Charadrou20.54
28Elekistrtas12.90
29Finikos127.84
30Glafkou76.75
31Milichou25.83
32Romanou12.80
33Selemnou17.75
34Serres57,000Agion Anargiron80.63
35Anatolikou Christou6.38
36Anatolikou Lefkonos21.21
37Ditikou Christou27.45
38Ditikou Lefkonos10.90
39Eleonos26.76
40Eptamilon11.79
41Kamenikion38.23
42Thessaloniki1,062,000Dendropotamou99.66
43Pilaias34.64
44Thermis51.52
45Trikala51,000Lithaiou211.28
46Volos77,000Aligarorematos20.07
47Anavrou10.40
48Arkoudorematos19.77
49Krausidona30.52
50Seskouliti32.97
51Xiria73.95
52Xanthi60,000Kimerion48.62
53Kosinthou240.63
54Laspopotamou53.09

References

  1. United Nations Department of Economic and Social Affairs & Population Division. World Urbanization Prospects: The 2018 Revision; United Nations Department of Economic and Social Affairs & Population Division: New York, NY, USA, 2019; p. 799. Available online: https://population.un.org/wup/Publications (accessed on 30 January 2024).
  2. Chen, M.; Zhang, H.; Liu, W.; Zhang, W. The global pattern of urbanization and economic growth: Evidence from the last three decades. PLoS ONE 2014, 9, e103799. [Google Scholar] [CrossRef]
  3. García-Nieto, A.P.; Geijzendorffer, I.R.; Baró, F.; Roche, P.K.; Bonneau, A.; Cramer, W. Impacts of urbanization around Mediterranean cities: Changes in ecosystem service supply. Ecol. Indic. 2018, 91, 589–606. [Google Scholar] [CrossRef]
  4. Bai, Y.; Wong, C.P.; Jiang, B.; Hughes, A.C.; Wang, M.; Wang, Q. Developing China’s Ecological Redline Policy using ecosystem services assessments for land use planning. Nat. Commun. 2018, 9, 3034. [Google Scholar] [CrossRef]
  5. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global Consequences of Land Use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef]
  6. Mallinis, G.; Boutsis, N.; Arianoutsou, M. Monitoring land use/land cover transformations from 1945 to 2007 in two peri-urban mountainous areas of Athens metropolitan area, Greece. Sci. Total Environ. 2014, 490, 262–278. [Google Scholar] [CrossRef] [PubMed]
  7. Euliss, N.H.; Smith, L.M.; Liu, S.; Feng, M.; Mushet, D.M.; Auch, R.F.; Loveland, T.R. The need for simultaneous evaluation of ecosystem services and land use change. Environ. Sci. Technol 2010, 44, 7761–7763. [Google Scholar] [CrossRef]
  8. Maring, L.; Blauw, M. Asset management to support urban land and subsurface management. Sci. Total Environ. 2018, 615, 390–397. [Google Scholar] [CrossRef]
  9. MEA. Ecosystems and Human Well-Being: Biodiversity Synthesis; World Resources Institute: Washington, DC, USA, 2005. [Google Scholar]
  10. Ouyang, Z.; Song, C.; Zheng, H.; Polasky, S.; Xiao, Y.; Bateman, I.J.; Liu, J.; Ruckelshaus, M.; Shi, F.; Xiao, Y.; et al. Using gross ecosystem product (GEP) to value nature in decision making. Proc. Natl. Acad. Sci. USA 2020, 117, 14593–14601. [Google Scholar] [CrossRef] [PubMed]
  11. Dobbs, C.; Eleuterio, A.A.; Amaya, J.D.; Montoya, J.; Kendal, D. The benefits of urban and peri-urban forestry. Unasylva 2018, 69, 22–29. [Google Scholar]
  12. Proutsos, N.D.; Solomou, A.D.; Petropoulou, M.; Chatzipavlis, N.E. Micrometeorological and Hydraulic Properties of an Urban Green Space on a Warm Summer Day in a Mediterranean City (Attica–Greece). Land 2022, 11, 2042. [Google Scholar] [CrossRef]
  13. Livesley, S.J.; Escobedo, F.J.; Morgenroth, J. The biodiversity of urban and peri-urban forests and the diverse ecosystem services they provide as socio-ecological systems. Forests 2016, 7, 291. [Google Scholar] [CrossRef]
  14. Berglihn, E.C.; Gómez-Baggethun, E. Ecosystem services from urban forests: The case of Oslomarka, Norway. Ecosyst. Serv. 2021, 51, 101358. [Google Scholar] [CrossRef]
  15. Baró, F.; Gómez-Baggethun, E. Assessing the potential of regulating ecosystem services as nature-based solutions in urban areas. In Nature-Based Solutions to Climate Change Adaptation in Urban Areas: Linkages between Science, Policy and Practice; Springer: Berlin/Heidelberg, Germany, 2017; pp. 139–158. [Google Scholar]
  16. García-Gómez, H.; Aguillaume, L.; Izquieta-Rojano, S.; Valino, F.; Avila, A.; Elustondo, D.; Santamaría, J.M.; Alastuey, A.; Calvete-Sogo, H.; González-Fernández, I.; et al. Atmospheric pollutants in peri-urban forests of Quercus ilex: Evidence of pollution abatement and threats for vegetation. Environ. Sci. Pollut. Res. 2016, 23, 6400–6413. [Google Scholar] [CrossRef] [PubMed]
  17. Seitz, J.; Escobedo, F. Urban Forests in Florida: Trees Control Stormwater Runoff and Improve Water Quality. Edis 2008, 5, 1–4. [Google Scholar] [CrossRef]
  18. Ramírez-Agudelo, N.A.; Porcar Anento, R.; Villares, M.; Roca, E. Nature-based solutions for water management in peri-urban areas: Barriers and lessons learned from implementation experiences. Sustainability 2020, 12, 9799. [Google Scholar] [CrossRef]
  19. Paletto, A.; Guerrini, S.; De Meo, I. Exploring visitors’ perceptions of silvicultural treatments to increase the destination attractiveness of peri-urban forests: A case study in Tuscany Region (Italy). Urban For. Urban Green. 2017, 27, 314–323. [Google Scholar] [CrossRef]
  20. Dou, Y.; Zhen, L.; De Groot, R.; Du, B.; Yu, X. Assessing the importance of cultural ecosystem services in urban areas of Beijing municipality. Ecosyst. Serv. 2017, 24, 79–90. [Google Scholar] [CrossRef]
  21. Kraxner, F.; Aoki, K.; Kindermann, G.; Leduc, S.; Albrecht, F.; Liu, J.; Yamagata, Y. Bioenergy and the city—What can urban forests contribute? Appl. Energy 2016, 165, 990–1003. [Google Scholar] [CrossRef]
  22. Cariñanos, P.; Calaza, P.; Hiemstra, J.; Pearlmutter, D.; Vilhar, U. The role of urban and peri-urban forests in reducing risks and managing disasters. Unasylva 2018, 69, 53–58. [Google Scholar]
  23. Li, Z.; Ning, K.; Chen, J.; Liu, C.; Wang, D.; Nie, X.; Hu, X.; Wang, L.; Wang, T. Soil and water conservation effects driven by the implementation of ecological restoration projects: Evidence from the red soil hilly region of China in the last three decades. J. Clean. Prod. 2020, 260, 121109. [Google Scholar] [CrossRef]
  24. Grima, N.; Corcoran, W.; Hill-James, C.; Langton, B.; Sommer, H.; Fisher, B. The importance of urban natural areas and urban ecosystem services during the COVID-19 pandemic. PLoS ONE 2020, 15, e0243344. [Google Scholar] [CrossRef] [PubMed]
  25. Kolimenakis, A.; Solomou, A.D.; Proutsos, N.; Avramidou, E.V.; Korakaki, E.; Karetsos, G.; Kontogianni, A.B.; Kontos, K.; Georgiadis, C.; Maroulis, G.; et al. Public Perceptions of the Socioeconomic Importance of Urban Green Areas in the Era of COVID-19: A Case Study of a Nationwide Survey in Greece. Land 2022, 11, 2290. [Google Scholar] [CrossRef]
  26. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; Da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  27. Diffenbaugh, N.S.; Giorgi, F. Climate change hotspots in the CMIP5 global climate model ensemble. Clim. Change 2012, 114, 13–822. [Google Scholar] [CrossRef] [PubMed]
  28. Food and Agriculture Organization of the United Nations. Forests and Water: International Momentum and Action; FAO: Rome, Italy, 2013. [Google Scholar]
  29. Springgay, E. Headwater Catchments: Foundation Pillars for Ecosystem Services. In Ecosystem Services of Headwater Catchments, 1st ed.; Křeček, J., Haigh, M., Hofer, T., Kubin, E., Promper, C., Eds.; Capital Publishing Company: New Delhi, India, 2017; pp. 3–6. [Google Scholar]
  30. Meng, Q.; Zhang, L.; Wei, H.; Cai, E.; Xue, D.; Liu, M. Linking Ecosystem Service Supply–Demand Risks and Regional Spatial Management in the Yize River Basin, Central China. Land 2021, 10, 843. [Google Scholar] [CrossRef]
  31. Bekri, E.S.; Kokkoris, I.P.; Skuras, D.; Hein, L.; Dimopoulos, P. Ecosystem accounting for water resources at the catchment scale, a case study for the Peloponnisos, Greece. Ecosyst. Serv. 2024, 65, 101586. [Google Scholar] [CrossRef]
  32. Vatitsi, K.; Ioannidou, N.; Mirli, A.; Siachalou, S.; Kagalou, I.; Latinopoulos, D.; Mallinis, G. LULC Change Effects on Environmental Quality and Ecosystem Services Using EO Data in Two Rural River Basins in Thrace, Greece. Land 2023, 12, 1140. [Google Scholar] [CrossRef]
  33. An, Y.; Zhao, W.; Li, C.; Ferreira, C.S.S. Temporal changes on soil conservation services in large basins across the world. Catena 2022, 209, 105793. [Google Scholar] [CrossRef]
  34. Renard, K.G.; Foster, G.R.; Weesies, G.A.; Porter, J.P. RUSLE: Revised universal soil loss equation. J. Soil Water Conserv. 1991, 46, 30–33. [Google Scholar]
  35. Bezak, N.; Ballabio, C.; Mikoš, M.; Petan, S.; Borrelli, P.; Panagos, P. Reconstruction of past rainfall erosivity and trend detection based on the REDES database and reanalysis rainfall. J. Hydrol. 2020, 590, 125372. [Google Scholar] [CrossRef]
  36. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses, a Guide to Conservation Planning; U.S. Department of Agriculture: Washington, DC, USA, 1978; Volume 537, p. 62.
  37. Hengl, T.; Mendes de Jesus, J.; Heuvelink, G.B.; Ruiperez Gonzalez, M.; Kilibarda, M.; Blagotić, A.; Shangguan, W.; Wright, M.N.; Geng, X.; Bauer-Marschallinger, B.; et al. SoilGrids250m: Global gridded soil information based on machine learning. PLoS ONE 2017, 12, e0169748. [Google Scholar] [CrossRef] [PubMed]
  38. Hawker, L.; Uhe, P.; Paulo, L.; Sosa, J.; Savage, J.; Sampson, C.; Neal, J. A 30 m global map of elevation with forests and buildings removed. Environ. Res. Lett. 2022, 17, 024016. [Google Scholar] [CrossRef]
  39. Meadows, M.; Jones, S.; Reinke, K. Vertical accuracy assessment of freely available global DEMs (FABDEM, Copernicus DEM, NASADEM, AW3D30 and SRTM) in flood-prone environments. Int. J. Digit. Earth 2024, 17, 2308734. [Google Scholar] [CrossRef]
  40. Pilesjö, P.; Hasan, A. A Triangular Form-based Multiple Flow Algorithm to Estimate Overland Flow Distribution and Accumulation on a Digital Elevation Model. Trans. GIS 2014, 18, 108–124. [Google Scholar] [CrossRef]
  41. McCool, D.K.; Foster, G.R.; Mutchler, C.K.; Meyer, L.D. Revised Slope Length Factor for the Universal Soil Loss Equation. Trans. ASAE 1989, 30, 1387–1396. [Google Scholar] [CrossRef]
  42. Desmet, P.J.J.; Govers, G. A GIS procedure for automatically calculating the USLE LS factor on topographically complex landscape units. J. Soil Water Conserv. 1996, 51, 427–433. [Google Scholar]
  43. Verde, N.; Kokkoris, I.P.; Georgiadis, C.; Kaimaris, D.; Dimopoulos, P.; Mitsopoulos, I.; Mallinis, G. National scale land cover classification for ecosystem services mapping and assessment, using multitemporal Copernicus EO data and google earth engine. Remote Sens. 2020, 12, 3303. [Google Scholar] [CrossRef]
  44. Panagos, P.; Borrelli, P.; Meusburger, K.; Alewell, C.; Lugato, E.; Montanarella, L. Estimating the soil erosion cover-management factor at the European scale. Land Use Policy 2015, 48, 38–50. [Google Scholar] [CrossRef]
  45. Efthimiou, N.; Psomiadis, E.; Papanikolaou, I.; Soulis, K.X.; Borrelli, P.; Panagos, P. A new high resolution object-oriented approach to define the spatiotemporal dynamics of the cover-management factor in soil erosion modelling. Catena 2022, 213, 106149. [Google Scholar] [CrossRef]
  46. Panagos, P.; Borrelli, P.; Meusburger, K.; van der Zanden, E.H.; Poesen, J.; Alewell, C. Modelling the effect of support practices (P-factor) on the reduction of soil erosion by water at European scale. Environ. Sci. Policy 2015, 51, 23–34. [Google Scholar] [CrossRef]
  47. SCS (Soil Conservation Service). SCS National Engineering Handbook, Section 4: Hydrology; Soil Conservation Service: Washington, DC, USA, 1972; pp. 1–44.
  48. Efstratiadis, A.; Koukouvinos, A.; Michaelidi, E.; Galiouna, E.; Tzouka, K.; Koussis, A.D.; Mamassis, N.; Koutsoyiannis, D. Description of Regional Approaches for the Estimation of Characteristic Hydrological Quantities, DEUCALION—Assessment of Flood Flows in Greece under Conditions of Hydroclimatic Variability: Development of Physically-Established Conceptual-Probabilistic; Department of Water Resources and Environmental Engineering—National Technical University of Athens, National Observatory of Athens: Athens, Greece, 2014. [Google Scholar]
  49. Savvidou, E.; Efstratiadis, A.; Koussis, A.D.; Koukouvinos, A.; Skarlatos, D. The curve number concept as a driver for delineating hydrological response units. Water 2018, 10, 194. [Google Scholar] [CrossRef]
  50. Panagos, P.; Borrelli, P.; Poesen, J.; Ballabio, C.; Lugato, E.; Meusburger, K.; Montanarella, L.; Alewell, C. The new assessment of soil loss by water erosion in Europe. Environ. Sci. Policy 2015, 54, 438–447. [Google Scholar] [CrossRef]
  51. Kokkoris, I.P.; Drakou, E.G.; Maes, J.; Dimopoulos, P. Ecosystem services supply in protected mountains of Greece: Setting the baseline for conservation management. Int. J. Biodivers. Sci. Ecosyst. Serv. Manag. 2018, 14, 45–59. [Google Scholar] [CrossRef]
  52. Ruckelshaus, M.; McKenzie, E.; Tallis, H.; Guerry, A.; Daily, G.; Kareiva, P.; Polasky, S.; Ricketts, T.; Bhagabati, N.; Wood, S.A.; et al. Notes from the field: Lessons learned from using ecosystem service approaches to inform real-world decisions. Ecol. Econ. 2015, 115, 11–21. [Google Scholar] [CrossRef]
  53. Cong, W.; Sun, X.; Guo, H.; Shan, R. Comparison of the SWAT and Invest models to determine hydrological ecosystem service spatial patterns, priorities and trade-offs in a complex basin. Ecol. Indic. 2020, 112, 106089. [Google Scholar] [CrossRef]
  54. Mitchell, M.G.; Schuster, R.; Jacob, A.L.; Hanna, D.E.; Dallaire, C.O.; Raudsepp-Hearne, C.; Bennett, E.; Lehner, B.; Chan, K.M. Identifying key ecosystem service providing areas to inform national-scale conservation planning. Environ. Res. Lett. 2021, 16, 014038. [Google Scholar] [CrossRef]
  55. Maes, M.J.; Jones, K.E.; Toledano, M.B.; Milligan, B. Mapping synergies and trade-offs between urban ecosystems and the sustainable development goals. Environ. Sci. Policy 2019, 93, 181–188. [Google Scholar] [CrossRef]
  56. Rozas-Vásquez, D.; Spyra, M.; Jorquera, F.; Molina, S.; Caló, N.C. Ecosystem Services Supply from Peri-Urban Landscapes and Their Contribution to the Sustainable Development Goals: A Global Perspective. Land 2022, 11, 2006. [Google Scholar] [CrossRef]
  57. Hosseini, S.M.; Mahjouri, N. Sensitivity and fuzzy uncertainty analyses in the determination of SCS-CN parameters from rainfall–runoff data. Hydrol. Sci. J. 2018, 63, 457–473. [Google Scholar] [CrossRef]
  58. Kumar, M.; Sahu, A.P.; Sahoo, N.; Dash, S.S.; Raul, S.K.; Panigrahi, B. Global-scale application of the RUSLE model: A comprehensive review. Hydrol. Sci. J. 2022, 67, 806–830. [Google Scholar] [CrossRef]
  59. Qiu, D.; Xu, R.; Wu, C.; Mu, X.; Zhao, G.; Gao, P. Vegetation restoration improves soil hydrological properties by regulating soil physicochemical properties in the Loess Plateau, China. J. Hydrol. 2022, 609, 127730. [Google Scholar] [CrossRef]
Figure 1. Location map of the major Greek cities.
Figure 1. Location map of the major Greek cities.
Land 13 00765 g001
Figure 2. The overall workflow of the methodology.
Figure 2. The overall workflow of the methodology.
Land 13 00765 g002
Figure 3. Average and SD values of soil conservation service provided by PUWs upstream of Greek cities.
Figure 3. Average and SD values of soil conservation service provided by PUWs upstream of Greek cities.
Land 13 00765 g003
Figure 4. Spatial distribution of soil conservation service in the examined PUWs of Greece.
Figure 4. Spatial distribution of soil conservation service in the examined PUWs of Greece.
Land 13 00765 g004
Figure 5. Average and SD values of the water retention service provided by PUWs upstream of Greek cities.
Figure 5. Average and SD values of the water retention service provided by PUWs upstream of Greek cities.
Land 13 00765 g005
Figure 6. Spatial distribution of water retention service in the examined PUWs of Greece.
Figure 6. Spatial distribution of water retention service in the examined PUWs of Greece.
Land 13 00765 g006
Figure 7. Average ES supply per altitudinal zone.
Figure 7. Average ES supply per altitudinal zone.
Land 13 00765 g007
Figure 8. Average ES supply per MAES category.
Figure 8. Average ES supply per MAES category.
Land 13 00765 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Stefanidis, S.; Proutsos, N.; Alexandridis, V.; Mallinis, G. Ecosystem Services Supply from Peri-Urban Watersheds in Greece: Soil Conservation and Water Retention. Land 2024, 13, 765. https://doi.org/10.3390/land13060765

AMA Style

Stefanidis S, Proutsos N, Alexandridis V, Mallinis G. Ecosystem Services Supply from Peri-Urban Watersheds in Greece: Soil Conservation and Water Retention. Land. 2024; 13(6):765. https://doi.org/10.3390/land13060765

Chicago/Turabian Style

Stefanidis, Stefanos, Nikolaos Proutsos, Vasileios Alexandridis, and Giorgos Mallinis. 2024. "Ecosystem Services Supply from Peri-Urban Watersheds in Greece: Soil Conservation and Water Retention" Land 13, no. 6: 765. https://doi.org/10.3390/land13060765

APA Style

Stefanidis, S., Proutsos, N., Alexandridis, V., & Mallinis, G. (2024). Ecosystem Services Supply from Peri-Urban Watersheds in Greece: Soil Conservation and Water Retention. Land, 13(6), 765. https://doi.org/10.3390/land13060765

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop