Next Article in Journal
Impact of Typical Land Use Expansion Induced by Ecological Restoration and Protection Projects on Landscape Patterns
Previous Article in Journal
Correction: Yan, S.; Growe, A. Regional Planning, Land-Use Management, and Governance in German Metropolitan Regions—The Case of Rhine–Neckar Metropolitan Region. Land 2022, 11, 2088
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Cost–Benefit Analysis for the Economic Evaluation of Ecosystem Services Lost Due to Erosion in a Mediterranean River Basin

1
Department of Soil, Plant and Food Sciences (DISSPA), University of Bari Aldo Moro, Via Amendola 165/a, 70126 Bari, Italy
2
Reagri s.r.l., Via Chiatona 62, 74016 Massafra, Italy
3
Green Agency C.C., Via A. De Curtis 14/g, 70023 Gioia del Colle, Italy
4
Faculty of Forest Sciences, Agricultural University, Rruga Paisi Vodica 1025, Koder-Kamez, 1025 Tirane, Albania
5
Department of Environment and Natural Resources, Faculty of Agriculture and Environment, University of Tirana, Rruga Paisi Vodica 1025, Koder-Kamez, 1025 Tirane, Albania
*
Author to whom correspondence should be addressed.
Land 2024, 13(9), 1512; https://doi.org/10.3390/land13091512
Submission received: 24 July 2024 / Revised: 10 September 2024 / Accepted: 14 September 2024 / Published: 18 September 2024

Abstract

:
Soil degradation in Europe is mostly represented by soil erosion that, as a sediment production mechanism, is the main environmental threat to many watersheds, including the Bovilla watershed (Tirana), useful for the supply of drinking water to the city, and therefore, the care of water quality is of particular interest. The soil erosion of the Bovilla watershed was monitored in a work set up in June 2017. Following this work, this research updates the previous data on soil loss at the Bovilla watershed in t/ha/year to September 2019 and focuses on the identification and monetary evaluation of the ecosystem services (ESs) linked to soil erosion (loss of carbon, loss of mineral elements, habitat quality, crop productivity, and sustainable tourism suitability). Then, we applied the replacement cost analysis to test the economic convenience and suggest the adoption of sustainable land management practices (e.g., reforestation) able to improve the quality water in this watershed. The study carried out demonstrates that the values of soil lost due to erosion vary depending on the type of land use (ranging from average values of 120.32 t/ha for bare land to values of 8.16 t/ha for wooded areas). Furthermore, from the application of monetary methods for the evaluation of some ecosystem services linked to erosion (loss of carbonaceous and mineral elements, habitat quality, productivity, suitability for sustainable tourism), it clearly emerges that the value of the productivity of agricultural crops varies from EUR 0 to 35,320.50/ha and that the service represents a more incisive service than the previous ones, so much so as to make the conversion of some agricultural land with high productivity values into wooded areas economically disadvantageous. The data from this study may help to develop Bovilla watershed management strategies for erosion and pollution control and sediment remediation mainly in agricultural lands. A program of measures can be effective for controlling soil erosion, but it must be implemented over long time frames, and it requires relevant investments from the public and private sectors, also with a view to increase the allocation of economic values of monetary compensation aimed at those who decide to start forestation projects on highly productive soils.

1. Introduction

Soil is a complex system within which a series of physical and biochemical processes are located that are vital for the characteristics of soil including the numerous species that live inside guaranteeing biodiversity and health for our globe [1], allowing the soil/plant system to produce biomass, monitor nutrient balance, and permit chemical recycling and water storage [2]. As for the environment in a global sense, soil can also be considered a natural capital considered as a resource capable of expressing, with goods and services, its profitability also in the future [3]; therefore, in addition to the biomass produced, ecosystems are able to perform a variety of services, defined by the literature of interest as ecosystem services (ESs), capable of providing public goods [4].
Anthropogenic pressure and above all human activities have exerted pressure leading to the degradation of soil functions, highlighting processes such as erosion, a decline in organic matter; local and diffuse contamination; sealing; compaction; a decline in biodiversity; salinization; floods; and landslides [5].
There is rich scientific literature that highlights how the erosive process is related to the loss of ecosystem services [6,7,8] and the strategies to be adopted to combat it, in response to the CAP and the 2030 Agenda SDGs [9]. The European Commission also defines measures to protect and restore soils and ensure that they are used sustainably in the approval document of the “Soil Strategy for 2030”. The main objective is to ensure that, by 2050, all member states of the European Union avoid consuming land (zero net land take) and ensure that their soils are “healthy” through concrete actions, to be implemented by 2030. The EU Soil Strategy also aims to ensure the following by 2050:
  • All European soils are healthy and more resilient and can continue to provide their fundamental services (ecosystem services);
  • Net land consumption is reduced to zero, and soil pollution is brought back to levels that are not harmful to people’s health or ecosystems;
  • The soils are protected and managed in a sustainable way, also restoring those currently degraded.
The document also indicates a series of fixed points and actions to be implemented in the coming years, such as the following:
  • The presentation of a European law for the protection of soil health, which incorporates all the contents of the strategy;
  • Sustainable soil management, i.e., a European soil management practice, promoted through specific actions of the Community Agricultural Policy, aimed at sharing and developing the best agronomic management practices, and through free soil analysis agricultural campaigns;
  • Encourage the accumulation of organic carbon in soils to mitigate the effects of climate change, including through legislative actions that protect and constrain wetlands and organic soils;
  • The prevention of desertification through the development of a common methodology to evaluate its level and prevent land degradation [10].
In recent years, the evolution of lifestyles, food needs, and, consequently, agricultural techniques and the exacerbation of climate change have caused an increasingly negative impact on the soil [11], highlighting the inability of soils to continue to guarantee food security and conditioning their ability to maintain the quality and quantity of ESs supplied, requiring more sustainable soil management [12,13]. In fact, the sustainability of the management of these ecosystems (forests, agricultural, and urban areas) passes from the essential relationship that exists between the natural soil capital and the ESs that it is able to provide, such as, the contrast to the degradation and impoverishment of soil fertility in forest systems, crop rotations, and minimum tillage techniques in agricultural areas, together with the ability of non-degraded soils to preserve and accumulate organic matter, water, and micronutrients, regulate greenhouse gas emissions, and increase resistance to pests and diseases [14]. On the other hand, all these excellent capacities are threatened by the adoption of highly intensive cultivation systems and practices to meet the ever-increasing food needs of humans and animals: this behavior has the consequence of reducing the ability of natural soil capital to provide positive externalities and to be essential for the sustainability of ecosystems and communities [15]. Non-sustainable land uses are always linked to disservices (loss of biodiversity, acidification erosion, nutrient losses, etc.) that affect soil productivity and, consequently, cultivation quality and yields, affecting the profitability of agricultural activity [16,17].
According to the EU, these essential soil functions are under several degradation pressures due to human activities. These degradation pressures include erosion; a decline in organic matter; local and diffuse contamination; sealing; compaction; a decline in biodiversity; salinization; floods; and landslides [5].
One of the most dangerous causes of soil degradation is erosion. Soil erosion causes not only the loss of land but also water pollution. Many factors are the basis of soil erosion: the geological structure, climate, morphology, lack of vegetation cover, and structure of land use. Solid and organic matter are transported through surface and ground water and are the main cause of a reduced volume of water in the basin and at the same time are the main pollution factors for the water and the environment.
Scholars have carried out research on modelling soil erosion rates under different climate and land use conditions [18,19].
The impoverishment of the soil, and therefore of the ESs related to it, inevitably causes negative impacts on productivity; the above could lead to risks of extremely high restoration costs [16,17] and therefore further degradation and economic devaluation.
The economic value of soil biodiversity stems from the idea that it can be valued as a natural capital asset, from which a flow of soil ecosystem services is produced [20,21]. Furthermore, each of these negative externalities also affects social well-being; in fact, a reduction in soil biodiversity can generate effects on the cost of environmental mitigation and on the increase in production costs, which are essential to counteract the loss of productivity. The adoption of all these sustainable techniques can generate long-term positive impacts by strengthening the role and effects of ESs [22,23].
ESs have been classified using different schematizations [4,24,25], which have in common the division into four categories: production/provisioning services, regulating services, habitat/supporting services, and information/cultural services, as reported in Table 1. Furthermore, in the literature, there are other more specific classifications both in soil-specific [11,26,27,28,29,30,31,32,33,34,35] and agroecosystem contexts [36,37,38,39,40] (Table 1).
The concept of payments for ecosystem services (PESs) requires that landowners and land users be compensated for providing a service to society, such as water flow regulation, the seizure of carbon, biodiversity conservation, etc. [41], adopting practices to ensure ecosystem conservation and rehabilitation [42].
Here, the focus is on shedding light on the role of the economic valuation of ecosystem services linked to soil loss in a study area represented by the Bovilla watershed to allow such values to be made explicit to society in general and policy-making in particular.
The Bovilla watershed constitutes the main reserve source of drinking water for the city of Tirana, and part of the water is also used for the irrigation of crops cultivated in that area. A report from the Albania National Forest Inventory from 2004 showed that the land cover/use is changing especially in the sloping and hilly areas in favor of cultivated areas. This may mean that this change occurred in less favorable environmental conditions (e.g., shallow soils, slopes, difficult to access, etc.). As a result of improper land use practices in the watershed, soil erosion has become an inhibiting factor for the sustainable development of agriculture and the conservation of natural resources including water resources.
Soil degradation is an issue of growing concern in Europe because 12% of the total European land area has been affected by water erosion and 4% by wind erosion [43]. The studies conducted in this area have investigated soil erosion through the application of the RUSLE model at a basin scale and a timely soil loss monitoring activity [42]. In fact, twenty-five plots (≈65–100 m2) were installed under different slope (10–20%, 21–40%, 41–60%, and >60%) and land cover (forest, pasture, agricultural, orchard) conditions and were used to measure the amount of sediment and runoff from June 2017 to May 2018 in previous work [44], and this was updated until September 2019 in the present study. This research focuses on identifying and evaluating the possibility of establishing schemes for payments for ecosystem services (PESs) in the Bovilla watershed (Tirane, Albania) to promote the adoption of sustainable land management practices that will result in improved quality water in this watershed.
This research aims to provide an update on the previous works in the Bovilla [44] watershed in order to (1) update the data of soil loss in/t/ha/year of the Bovilla watershed from June 2018 till September 2019, (2) monetarily estimate the ecosystem services linked to the erosion of the Bovilla basin, and (3) suggest an implementation of PESs in particular for the sub-watersheds Ranxe, Vilez, and Zall-Bastar that were selected as critical erosion areas and the main hotspots [45], which contribute to about 65% of the sediment yield of the watershed. The data from this study may help to develop Bovilla watershed management strategies for erosion and pollution control and sediment remediation mainly in agricultural lands. To achieve these research objectives, using the first results on Bovilla soil loss and the erosion map published by the authors of [44], we identified the main ecosystem services directly and indirectly linked to soil erosion, we assessed and monetarily quantified these ecosystem services, and finally, we compared the soil loss to the basin restoration strategies (e.g., reforestation) through a cost–benefit analysis.

2. Materials and Methods

2.1. Study Area

The study area is in the north-eastern part of Tirana city (Figure 1), with a width of 95 km2, belonging to the Tirana (52.2 km2) and Kruja (39.8 km2) municipalities.
The Bovilla reservoir is 15 km away from Tirana city with a maximum water filling capacity of 80 million m3. The total area of the Bovilla reservoir is 4.6 km2, and the maximal depth is 53 m. The potable water volume of the Bovilla reservoir is estimated at about 55 mln m3/year, and the rest might be used for irrigation. The Bovilla watershed comprises six sub-watersheds (Table 2).
The Bovilla watershed is part of the Dajti National Park according to the Decision of the Council of Ministers No. 402 of 21 June 2006. The dam of the Bovilla reservoir has a height of 91 m as well as a length of 130 m and an elevation of about 321 m above sea level (a.s.l.) and is in the mouth of the Terkuza river. The climate is typical of the Mediterranean with relatively cold and wet winters and warm and dry summers. The area has a wide altitudinal range, directly affecting the climate variability. The mean annual temperature is 13.9 °C, the rainfall is 1718.6 mm/year, and the potential evaporation is 916 mm. Rainfall data are characterized by a non-uniform monthly distribution, where most of the precipitation falls during the autumn and winter months, being the main source of water supply for the hydrologic network in the Bovilla watershed.
Each watershed has its own runoff pattern and erosion rate which depends on the (i) geology, (ii) vegetation, (iii) land use, and (iv) topography. According to erodibility, the rocks of the Bovilla watershed can be classified as (i) high-eroded rocks (Quaternary deposits, 6.8 km2 or 7% of the watershed), (ii) medium-eroded rocks (Flysch and Molasse, 70.2 km2 or 71.6% of the watershed), and (iii) low-eroded rocks (Limestone and Dolomite, 21 km2 or 22% of the watershed). The area is characterized by a strongly variable terrain with alternating mountainous, hilly, and flat areas. Based on the watershed DEM, around 80% of the Bovilla watershed has a slope from 15 to 30 degrees. The elevation of the Bovilla watershed varies from 300 to 1800 m a.s.l., and 80% of the watershed area has an altitude from 300 to 1000 m a.s.l. Some pedological properties of the main soils are provided in Table 3.
According to their texture, the soils studied can be divided into two broad soil groups, namely light silty soils (5 out of the 8 studied soils were silt loam) and medium soils (3 out of the 8 studied soils were silty clay loam). In general, the soils of the first group have comparably low organic matter contents and low aggregate stability. Therefore, these soils have a high risk of erosion where runoff occurs. The soils of the second group have a high content of silt and consequently are not as stable. They are susceptible to erosion, as the content of organic matter is low.
The main land use classes identified in the study area were (i) broadleaved forests, (ii) riparian vegetation along rivers and water streams, (iii) agriculture land, and (iv) pastures. Besides the identification of vegetation cover, we have estimated the respective area for each land use class within the Bovilla watershed (Table 4) and carried out their mapping by means of ArcGis software, version 9.5 (Figure A1). Based on the site conditions identified within the Bovilla watershed, we argue that the study area is particularly sensitive to soil erosion due to deforestation in the high-elevation areas and the presence of sparse vegetation areas.
In the previous study [44], areas with a high erosion risk (hotspot areas) were determined using the Revised Universal Soil Loss Equation (RUSLE) and GIS.

2.2. Determination of Soil Loss

Sedimentation in the Bovilla reservoir has caused many problems with water quality, a decrease in storage capacity, and water supply for the population as a result of turbidity rising. In order to update the data on soil loss at the Bovilla watershed, we started from a previous work [44] by using the same methodology. In that previous work, the sub-watersheds Ranxe, Vilez, and Zall-Bastar were selected as critical erosion areas, which contribute to about 65% of the sediment yield of the watershed, and in these critical areas, it is necessary to apply PESs programs [44].
Then, 25 experimental sample plots with regular and irregular shapes, located in the villages of Zall Bastar and Mner Sperm, were installed near the closing section of the selected sub-basins to monitor the soil loss caused by water erosion (weight of the soil derived from an area ranging from 65 to 100 m2). The land use (forest, pasture, agricultural, bare lands, degraded forest land, and orchard) and the slopes (0–20%, 21–40%, 41–60%, and >60%) were the criteria to locate the hotspot areas for soil erosion. A mixture of all sediments collected during monitoring during the two-year test was examined in the laboratory, determining the exact weight of the sediments.
Each experimental erosion plot was coded in the field, and then, respective information regarding the geographic position, elevation, land use, slope, and exposure was recorded. For the coding of sample plots, the following rule was used: e.g., ZB1F/M1F, where the first letter(s) denotes the village name—Zall Bastar (ZB) or Mner (M); the number denotes the sample plot number; and the last letter (F) denotes the land use, such as “forest (F)”. The area of each experimental sample plot was measured by means of a highly accurate Trimble GPS device. This information is summarized in Table A1. The authors of this research [44] elaborated on these data for the period from June 2017 to May 2018, while we update the study to September 2019. The map of the spatial distribution of soil erosion for the period June 2017–May 2018 [44] is reported in Figure A1 and clearly demonstrates that soil loss of less than 2.46 t/ha/year (tolerable soil loss in EU countries) covers only 96% of the study area, and the rest is dominated by higher erosion values (over 2.46 t/ha).

2.3. Identification of Ecosystem Services Correlated to Soil Erosion

We selected five ecosystem services linked to soil loss: the loss of carbon (regulating services), loss of mineral elements (K and P—supporting functions), habitat quality (supporting functions), crop productivity (provisioning services), and sustainable tourism suitability (cultural ecosystem services) (Figure 2).

2.4. Application of Monetary Evaluation for Ecosystem Services Linked to Erosion

The economic evaluation of the ecosystem services was standardized in European currency (EUR); the value of each ecosystem service was estimated for each single plot and then summed. The soil quality was monitored with respect to K, P, and organic matter. For each plot, soil samples were collected and sent to the laboratory to analyze each parameter using the standard methods reported in Table A2. These data were used to quantify the monetary value of the loss of carbon and loss of mineral elements. The five ecosystem services and the related estimation methods are described as follows:
  • Loss of carbon (regulating services): The starting data are represented by the values expressed in average t/year (three-year period 2017–2019) of sediment loss due to the erosive process. Starting from these data, it was possible to quantify the percentage of organic substances [46] and therefore the content in kg of organic matter/ha. The transition from kg of organic matter per kg of carbon was carried out using the value of 58% (application of the inverse formula of the Van Bemmelen factor = 1.74). The transformation of C (kg) into CO2 was carried out with the application of the stoichiometric ratio resulting from the oxidative process of 12 kg of carbon, therefore equal to 3.66 (12 kg of C + 32 kg of O2 = 44 kg of CO2). For the conversion into EUR, the prices (EUR/t CO2) [47] for the year 2023 were used and compared with the average data for 2023 for the last 52 weeks [48] = 60.86 + 98.43 (average 79.45 EUR/t CO2).
  • Loss of mineral elements (K and P—supporting functions): The content of K and P was evaluated from an economic point of view using the prices of mineral fertilizers (monopotassium phosphate and potassium nitrate), which contain a higher content of K and P, respectively. Starting from the values expressed in average t/year (three-year period, 2017–2019) of sediment loss due to the erosive process, the potassium and phosphorus contents were calculated (obtained from laboratory analyses for the 25 plots using the Mehlich 3-Extrable Elements Soil Analysis method), expressed in mg/100 g of analyzed soil, subsequently converted into kg/ha. The result expressed in kg was multiplied by the average trading values (EUR), obtained from market analyses relating to monopotassium phosphate and potassium nitrate. The indicative amount is around 3–4 EUR/kg. It should be noted that this study estimates the loss of P and K by analyzing the soil and not the sediment. It is known that there is usually an enrichment of nutrients in the sediments compared to the soils (sediment enrichment). Therefore, in this study, there is an underestimation of the concentrations of P and K and, therefore, of the performance in terms of the loss of the specific ecosystem service.
  • Economic assessment of habitat quality (supporting functions): In this study, the habitats present in the basin were evaluated in terms of “habitat quality”, and a monetary value was subsequently attributed to them. “Habitat quality” was assessed using the Invest software (version 3.14.2), developed by the University of Stanford—Natural Capital Project [49]. The Invest software considers “habitat quality” as the ability of the ecosystem to provide conditions appropriate for individual and population persistence, and it is considered a continuous variable in the model, ranging from low to medium to high, based on the resources available for survival, reproduction, and population persistence, respectively [50]. Habitats with a high quality are relatively intact and have a structure and function within the range of historic variability. “Habitat quality” depends on a habitat’s proximity to human land uses and the intensity of these land uses. Generally, “habitat quality” is degraded as the intensity of nearby land use increases [51,52,53]. The model runs using raster data, where each cell in the raster is assigned a Land Use Land Cover (LULC) class, which can be a natural (unmanaged) class or a managed class. LULC types can be given at any level of classification detail. Besides a map of LULC and data that relate LULC to habitat suitability, the model also requires data on the habitat threat density and its effects on habitat quality. In general, we consider threats to be human-modified LULC types that cause habitat fragmentation, edge, and degradation in neighboring habitats. The following inputs were used for the case study: (i) the land use map (current landcover map) and (ii) the threat maps. All maps used in the Invest software have a cell size of 20 m. The characteristics of “habitat” or “no-habitat” were attributed to each LULC type. As regards the threats, this study considered the presence of roads, the presence of crops and agricultural activities, the presence of urbanized areas, and the proximity to contexts affected by erosion.
  • The result of the processing is a habitat quality map in which each cell is given a score ranging from 0 (low HQ) to 1 (high HQ). The input maps and the resulting “habitat quality map” are shown in the following figures (Figure 3, Figure A2, Figure A3, Figure A4, Figure A5 and Figure A6). The HQ values obtained (range 0–1) were then multiplied by the economic value given by the authors of [54] which assigns an economic amount to each area with different land use. The assigned value, expressed in 2007 USD/ha/year, was converted into EUR 2018 (intermediate year of this work) through the purchasing power parities (PPPs) [55], subsequently adjusted to 2018 by consulting the data available from the Italian Statistics Institute [56].
  • Crop productivity (provisioning services): The economic value of the agricultural production obtainable from the plots under investigation was quantified through the application of the Italian Ministerial Decree no. 236,537 of 5 May 2023 that established the identification of the standard values for plant production to determine the value of the average annual production and the maximum values that can be insurable on the subsidized market and to join mutual funds (EUR/ha). For forestry crops, however, the formula 0.5*d2*H (Bouvard formula) was applied. The diameters, heights, and number of plants per hectare divided by age were reported in the Bovilla database simple plot—forest info (not published). The total volume (q/ha) was multiplied by 6 EUR/q (for coniferous trees) and by 14 EUR/q (for broadleaved trees, average purchase and sale value of coniferous and broadleaved forests expressed in EUR).
  • Sustainable tourism suitability (cultural ecosystem services): For the evaluation of this ecosystem service, we considered the “sustainable tourism suitability” (STS) defined as the potential of the territory to offer services for sustainable tourism. To estimate the performance of this service within the park territory, we started from the following assumption: a territory, or parts of it, can have a performance in terms of tourism suitability depending on the presence or absence of some territorial elements: roads, hiking trails, habitat/naturalistic areas, restaurants, bar, hotels, sports facilities, and historical–cultural elements. A total of 4 macro-groups of territorial elements were defined and are shown in Table 5.
  • To evaluate the spatial distribution of the ecosystem service, we then proceeded to evaluate the spatial distribution of each characteristic by carrying out a density analysis in a GIS environment. The density map of each feature was normalized (Figure A7, Figure A8, Figure A9 and Figure A10). An overlay analysis was then performed to obtain a performance map of STS for which each cell (cell size of 20 m) was attributed with a value within the range of 0–4 (low–high STS) as a result of the sum of the spatial density values of the individual territorial elements (Figure 4).
  • The values of direct tourist flows within the “Parku Kombetar Mali i Dajtit” National Park provided by the National Agency for Protected Areas were used. The average annual flow (three-year period, 2021–2023) was 126,057 units that used the cableway to reach the protected area (within which the Bovilla Basin lies), departing from the city of Tirana. Official investigations show that the price of the cableway ticket is around 2 EUR/each. Therefore, the economic value in terms of tourist attractiveness of the entire protected area, considering only the value of the public transport ticket, is equal to 252,115 EUR/year. The extent of the entire protected area is 29,193.27 ha; thus, each single area can be standardized to the value of EUR 8.63.

2.5. Application of the Replacement Cost Method

We used the replacement cost method to transform bare ground, orchards, and pastures into forest soil to compare the result with the total value of the previous five ecosystem services. The transformation value (i.e., the replacement cost) of the agricultural land, pastures, and bare land under investigation was economically compared with the creation of an initial afforestation plant. The determination of the transformation value (EUR/ha) was obtained with two different methods. First, the economic values of the transformation of agricultural land into forestation interventions were used, in compliance with the Rural Development Programme (RDP) of the Puglia Region, which is the main funding instrument within the European Agricultural Fund for Rural Development (FEASR) through which the Puglia regional government promotes key interventions for the development of rural areas. Subsequently, the second method was adopted by applying the unit values (EUR/ha) of the cost items taken from the price list of the Puglia Region, attached to the Decision of the Regional Council (DGR) LLP/DEL/2023/00012. In both methods, the economic values of income loss and the maintenance values of the forestation systems were added to the values obtained, multiplied by 12 years (see the notice of the first afforestation of agricultural land—Measure 8.1 of the RDP of Puglia Region 2014–2020) [57].

3. Results

3.1. Determination of Soil Loss

The amount of soil loss for the whole monitoring period was 2705.6 t/ha, while the mean soil loss at the sample plot level was 108.22 t. The soil loss varied between land use and years (Table 6).
The mean soil loss varied between land use categories. Over the monitoring period, the lowest values of soil loss measured in forest lands ranged from 5.24 t/ha (2019) to 12.61 t/ha (2018), while the maximum values were reached in bare lands, from 56.72 t/ha (2019) to 154.76 t/ha (2017) (Table 7).
On average, the mean soil loss over the monitoring period (2017–2019) was 36.07 t/ha, ranging from 8.16 t/ha (forest) to 120.32 t/ha (bare land). We noticed that the soil loss per 1 hectare varied between sample plots, and the largest values were reached in 2019 (Figure 5). As expected, different amounts of runoff and sediment were generated from soils with different land use statuses. The trend of soil loss was as follows: bare land > orchard > degraded forest > agricultural land > pasture > forest. The soil loss varied between land uses from 8.16 t/ha in forest lands to 120.32 t/ha in bare lands. The soil loss proportionally increased when the slope rose from 20% to more than 61%.

3.2. Use of Monetary Methods for the Evaluation of Ecosystem Services Linked to Erosion

The value of the five analyzed ecosystem services (loss of carbon, loss of mineral elements K and P, habitat quality, crop production, and cultural ecosystem services) is reported in Table A3. From the obtained results, we observed the following:
  • The economic value of carbon loss (CO2/ha) varies from 0.01 to 348,94.00 EUR/t;
  • The economic value of the loss of chemical elements varies from 0 to 1641.19 EUR/ha for K and from 0 to 224.22 EUR/ha for P;
  • The economic value of the habitats varies from 0 to approximately 5000.00 EUR/ha;
  • The economic value of crop production varies from 0 to 35,320.50 EUR/ha; thus, the economic weight/ha of this service is more incisive than that of the previous ones.

4. Discussion

4.1. Replacement Cost

The value of the calculated replacement cost is reported in Table A4. The value of the restoration cost granted by the Puglia Region (estimated only for agricultural land, meadows, and pastures) is around 10,500.00 EUR/ha considering the contribution values provided by the Rural Development Program [57], while it is around EUR 18,192.00 if the estimate is made only with data from the Puglia Region price list for the year 2023 [58]. The already existing forests were not investigated.

4.2. Comparison

Comparing the economic values of restoration (afforestation) with the data derived from the application of the cost–benefit analysis, the following key points emerge:
  • It is always economically convenient to change the intended use of the land from pasture to forest (where possible and not in conflict with current environmental legislation), from agricultural land (corn or fodder) to forest, and from bare land to forest (where possible and not in conflict with current legislation). There is economic convenience if the values derived from the cost–benefit analysis are compared with the economic values of replacement cost from the Puglia Region price list (Table A5, value “Replacement cost”), while there is not economic convenience if the replacement cost values derive from the application of Measure 8.1 of the Rural Development Program of the Puglia Region [57].
  • It is possible to compare the values described above for a period of 12 years (minimum period to allow a forestry seedbed to become young and stable, to grant the economic prize for forest maintenance and compensate the loss of income from the previous crop within the RDP of Puglia Region). Also, in this case, the considerations of convenience reported in point 1 are confirmed.
  • It is more convenient from an economic point of view to transform bare land or pasture to forest (where in compliance with the legislation), since the bare land and pastures are almost non-productive. It is less convenient if the conversion takes place from productive agricultural land, even simply arable land (Table A5, value “Economic Convenience EUR/ha”), because the value of the ecosystem service of production (supply of products) is calculated with the MASAF 2023 Standard Value [59].
  • It is also underlined that for the ZB4O orchard, for the M1O orchard, and for the ZB2A arable land, the convenience of carrying out the intervention is uneconomical even considering the replacement cost values and the values of ordinary maintenance that are higher than the price list of the Puglia Region (Table A5).
  • It is not economically convenient to transform highly productive land (vineyards, plum orchards, horticultural crops, etc.), which is characterized by ecosystem services of very high economic value and not appropriately compensated by the replacement costs, maintenance prize, and income loss. The maintenance prize should be raised from EUR 2846.00 to at least EUR 7566.20 [57], as described in Table A5 (column “Total cost”). This is due to more frequent extreme weather events that drastically reduce the rooting of new plants, especially in the summer period, with the need to carry out at least six irrigations per year.

5. Conclusions

This research updated the previous data on soil loss at the Bovilla watershed to September 2019 and focused on the identification and monetary evaluation of the ecosystem services (ESs) linked to soil erosion. Agriculture is a key driver of soil erosion and soil contamination by fertilizers. The presence of soil erosion in the Bovilla watershed caused by agriculture is closely related to farming practices (land ploughing, irrigation, crop rotations, crop management, etc.), which might be considered as indicators of pressures on watershed erosion. Therefore, agriculture management should be monitored in the Bovilla watershed as an important factor influencing soil erosion; measures should be implemented to find suitable remedial measures and encourage farmers to adopt soil conservation techniques.
The evaluation of the on-site and off-site impacts of soil erosion is closely related to the pedoclimatic characteristics and land uses. This evaluation is not only difficult but also expensive, and specific data are often not available. Indeed, to overcome the lack of data, the evaluation of the ecosystem services was performed by gathering data from official documents of the Rural Development Programme of Puglia Region (Southern Italy) and from the Italian Institute of Statistics. This was conducted because Puglia and some areas of Albania have similar pedoclimatic conditions and land uses. Therefore, future research should be focused on collecting data directly from the Bovilla watershed.
Land use and management result from human activities, and they are the most important factors that influence and control accelerated soil erosion. Land cover, especially for forest land, has been altered to degraded forest land, affecting the erosion rates in the Bovilla watershed. The main factors in land cover change in the study area are illegal logging and forest fires. Therefore, the promotion of funding schemes in the framework of payments for ecosystem services is essential to support activities including reforestation, regeneration, dam construction, and so on.
The economic consequences of soil erosion in the Bovilla watershed are multifold. The loss of topsoil is associated with a reduction in soil productivity associated with a reduction in arable crop yield (e.g., wheat and corn), forage, orchards, and vineyards, which are crucial to ensure the food security of the local population.
For the specific conditions of Albanian forests, their management should be focused on the principles of sustainable mountain development. In this way, two different goals can be pursued: (1) to improve the natural capital and quality of life of upland populations and (2) to ensure ecosystem stability also in the urbanized territories downstream of the watershed basin. This means a greater involvement of local populations in planning and implementing management strategies and programs, often with more direct benefits to mountain populations.
Broadleaved forests, coniferous forests, and herbaceous vegetation provide good-quality fresh water. To increase the water quality parameters in the Bovilla watershed, it is crucial to enhance the green surface, promote reforestation, and protect soil from erosion. We are aware that the correct use of land cannot be linked exclusively to forestation or reforestation interventions and should be guaranteed by the adoption of sustainable agronomic techniques for food production [60].
Nevertheless, national and international policies aimed at the maintenance and expansion of habitats and natural territories to the detriment of agriculture, especially in areas with intensive agriculture, should make a better allocation of resources in future programming, rewarding as a priority, and with much greater contributions to current, the execution of afforestation works, especially in terms of compensation for income loss from previous crops.
Overall, this research highlighted that the combination of the monetary assessment of ESs linked to soil erosion and mapping on a spatial scale can help decision-makers to develop sustainable territorial planning and policies; quantifying the value of ESs could become fundamentally important in the decision-making processes related to land use.

Author Contributions

Conceptualization, G.O.P., A.D.B., C.A. and R.R.; methodology, G.R.R., V.L., E.T. and R.R.; software, V.L.; validation, G.R.R., G.O.P. and R.R.; formal analysis, M.C. and V.L.; investigation, G.R.R.; resources, M.C.; data curation, E.T. and R.K.; writing—original draft preparation, G.R.R. and M.C.; writing—review and editing, G.O.P., R.K., A.D.B., C.A. and R.R.; visualization, G.O.P.; supervision, A.D.B., C.A. and R.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no specific grant from funding agencies in the public, commercial, or non-profit sectors.

Data Availability Statement

The data presented in this study are available on request from the corresponding author (the data are not publicly available due to privacy).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. General information about the sample plots in Bovilla watershed.
Table A1. General information about the sample plots in Bovilla watershed.
No.Sample Plot CodeLocationElevation
a.s.l (m)
Slope Area
(m2)
Soil Type
NEDegree%
1ZB1AN 41°26′31.71″E 19°54′55.79″3252.8569.3Fluvic Kambisols
2ZB2AN 41°26′34.71″E 19°54′57.47″327.62.54.689.5Fluvic Kambisols
3ZB3AN 41°26′38.18″E 19°54′53″321.5713.2145.3Regosols Eutric
4ZB1FN 41°26′42.45″E 19°54′51.23″32918.434.672.9Haplic Kambisol
5ZB1BN 41°26′41.57″E 19°54′59.25″35731.86242.5Haplic Leptosols
6ZB2BN 41°26′38.50″E 19°55′01.79″35835.77246Haplic Leptosols
7ZB3BN 41°26′38.76″E 19°55′01.91″36036.57441.9Haplic Leptosols
8ZB3FN 41°26′37.62″E 19°55′03.04″3632374466.9Haplic Leptosols
9ZB4FN 41°26′38.40″E 19°55′01.79″35332.66474.5Haplic Kambisols
10ZB4BN 41°26′36.63″E 19°55′02.05″35023.34366.7Haplic Leptosols
11ZB2PN 41°26′35.97″E 19°54′58.94″34019.83670Kambisols Eutric
12ZB3PN 41°26′36.97″E 19°54′57.89″33613.624.285Kambisols Eutric
13ZB4PN 41°26′36.23″E 19°55′00.00″34119.83665.6Kambisols Eutric
14ZB5PN 41°26′36.64″E 19°55′00.98″3467.41385.5Kambisols Eutric
15ZB4ON 41°26′37.52″E 19°54′58.55″34215.82872.3Kambisols Eutric
16ZB2FN 41°26′35.70″E 19°54′59.23″33731.86259.3Kambisols Eutric
17M1ON 41°27′45.42″E 19°55′51.39″49120.24174.9Haplic Luvisols
18M1PN 41°27′44.86″E 19°55′51.97″474295519.4Haplic Luvisols
19M1FN 41°27′47.20″E 19°55′53.49″47321.64063.6Haplic Luvisols
20M1DFN 41°27′47.16″E 19°55′54.35″46717.83273.8Kambisols Eutric
21M2PN 41°27′47.34″E 19°55′52.22″48421.74072.9Haplic Luvisols
22M1BN 41°27′54.53″E 19°55′47.63″507.525.648159.2Kambisols Eutric
23M2FN 41°27′53.72″E 19°55′46.83″49715.42872.9Haplic Leptosols
24M3PN 41°27′45.14″E 19°55′53.18″46832.26378.5Kambisols Eutric
25M4PN 41°27′46.25″E 19°55′54.15″46533.06539.4Haplic Luvisols
Figure A1. Land cover map of Bovilla watershed.
Figure A1. Land cover map of Bovilla watershed.
Land 13 01512 g0a1
Table A2. Methods and standards used for chemical compound analysis.
Table A2. Methods and standards used for chemical compound analysis.
Sample Plot CodeK
(Mehlich 3-Extrable Elements Soil Analysis)
(kg/ha)
P
(Mehlich 3-Extrable Elements Soil Analysis)
(kg/ha)
Organic Matter Converted into CO2 1 (kg/ha)
ZB5P124.8822.287384.509
ZB4P94.3981.729289.138
ZB3P123.1792.256385.254
ZB4O820.593140.1374183.873
ZB2P222.3074.072688.056
ZB4F100.9632.088768.651
ZB3F122.9602.54411.441
ZB2F117.0222.1434.075
ZB1P101.7322.10410.868
ZB3A−306.43213.57410.401
ZB2A495.35221.94310.571
ZB1A543.46024.07410.741
ZB4B462.35821.1450.679
ZB3B502.35119.1350.424
ZB2B92.2841.94730.382
ZB1B94.8743.9970.148
M1B267.0733.3732.717
M2F7.0097.0094.818
M1O295.574.9657.854
M1P317.654.8403.545
M1F125.810.7594.818
M1DF203.4674.9064.670
M2P110.3551.685.455
M3P111.8561.7045.307
M4P155.4312.3685.349
1 Application of the inverse formula of the Van Bemmelen factor: organic C content equal to 58% of the weight of the organic substance; transition from organic C to CO2 by applying the stoichiometric formula 12 kg of C + 32 kg of O2 = 44 kg of CO2, from which it can be deduced that the ratio between C and CO2 is 3.66 (values expressed in kg, subsequently indicated as t for the purpose of a monetary evaluation).
Figure A2. Land use.
Figure A2. Land use.
Land 13 01512 g0a2
Figure A3. Erosion map.
Figure A3. Erosion map.
Land 13 01512 g0a3
Figure A4. Roads.
Figure A4. Roads.
Land 13 01512 g0a4
Figure A5. Urban areas.
Figure A5. Urban areas.
Land 13 01512 g0a5
Figure A6. Agricultural areas.
Figure A6. Agricultural areas.
Land 13 01512 g0a6
Figure A7. Road density map.
Figure A7. Road density map.
Land 13 01512 g0a7
Figure A8. Network of paths density map.
Figure A8. Network of paths density map.
Land 13 01512 g0a8
Figure A9. Habitat quality.
Figure A9. Habitat quality.
Land 13 01512 g0a9
Figure A10. Point of interest density map.
Figure A10. Point of interest density map.
Land 13 01512 g0a10
Table A3. The value of the five ecosystem services investigated in the plots.
Table A3. The value of the five ecosystem services investigated in the plots.
Sample Plot CodeLand Use Classes
(Corine Land Cover 2012)
I
(Regulating)
II
(Regulating)
III
(Provisioning)
IV
(Supporting)
V
(Cultural)
CO2
(EUR/t)
Habitat Quality (EUR/ha)Agricultural and Forestry Production
(EUR/ha)
K
(EUR/ha)
P
(EUR/ha)
CES
(EUR/ha)
ZB5PPasture32.07 2892.09 326.00 249.76 3.66 5.57
ZB4PPasture24.11 2930.39 326.00 188.80 2.77 5.51
ZB3PPasture32.13 2923.44 326.00 246.36 3.61 5.42
ZB4OOrchard348.94 0.0035,320.501641.19224.222.20
ZB2PPasture57.38 2931.38326.00444.616.525.45
ZB4FConifer forests64.11 0.002600.88201.933.342.45
ZB3FConifer forests0.95 2377.922600.88245.924.076.00
ZB2FOak forest with pasture0.34 2203.63326.00234.04 3.43 5.46
ZB1PConifer forests0.91 2524.452701.08203.463.376.39
ZB3AAgricultural land (corn)0.87 0.003948.00612.8621.723.11
ZB2AAgricultural land (corn and beans)0.88 0.00 5757.50 990.7135.11 2.20
ZB1AForage0.90 0.002086.001086.9238.52 2.20
ZB4BBare lands0.06 0.00 0.000.00 0.005.89
ZB3BBare lands0.04 0.00 0.000.000.002.49
ZB2BBare lands0.03 0.00 0.000.000.002.45
ZB1BBare lands0.01 0.00 0.000.000.006.39
M1BBare earth with renaturalization0.23 0.00 0.00534.155.401.47
M2FBroadleaved forests0.40 0.00 1147.9314.0211.221.48
M1OOrchard0.66 5156.14 32,843.330.000.005.60
M1PCultivated meadow (alfalfa)0.30 3867.062114.00635.317.745.64
M1FBroadleaved forests0.40 2900.631147.93251.63 1.215.67
M1DFDegraded forest0.39 2905.11653.520.000.005.67
M2PCultivated meadow (alfalfa)0.45 3840.832114.00220.712.695.67
M3PPasture—meadow (alfalfa)0.44 3863.92770.00223.712.735.60
M4PPasture—meadow (alfalfa)0.45 3862.22770.00310.863.795.60
Table A4. Comparison between replacement cost values and authorized costs reported in the Rural Development Programme (RDP) of Puglia Region, with detail of the values of maintenance prize and income loss for each individual plot.
Table A4. Comparison between replacement cost values and authorized costs reported in the Rural Development Programme (RDP) of Puglia Region, with detail of the values of maintenance prize and income loss for each individual plot.
Sample Plot CodeSub-Basin NameLand Use Classes
(Corine Land Cover 2012)
Replacement Cost
(EUR/ha)
Authorized Costs (RDP Puglia)
(EUR/ha)
Loss of Income
(EUR/ha)
Loss of Income for 12 Years
(EUR/ha)
Maintenance Prize
for 12 Years (EUR/ha)
Total Cost 1
(EUR/ha)
ZB5PZall BastarPasture18,192.0010,500.000.00 0.00 34,154.4052,346.40
ZB4PZall BastarPasture18,192.0010,500.000.00 0.00 34,154.4052,346.40
ZB3PZall BastarPasture18,192.00 10,500.000.00 0.00 34,154.4052,346.40
ZB4OZall BastarOrchard18,192.00 10,500.00950.0011,400.00 34,154.4063,746.40
ZB2PZall BastarPasture18,192.00 10,500.000.000.0034,154.4052,346.40
ZB4FZall BastarConifer forests0.00 0.00 0.00 0.00 0.00 0.00
ZB3FZall BastarConifer forests0.00 0.00 0.00 0.00 0.00 0.00
ZB2FZall BastarOak forest with pasture0.00 0.00 0.00 0.00 0.00 0.00
ZB1PZall BastarConifer forests0.00 0.00 0.00 0.00 0.00 0.00
ZB3AZall BastarAgricultural land (corn)18,192.00 10,500.00700.008400.0034,154.4060,746.40
ZB2AZall BastarAgricultural land (corn and beans)18,192.00 10,500.00700.00 8400.0034,154.4060,746.40
ZB1AZall BastarForage18,192.0010,500.0064.00768.0034,154.4053,114.40
ZB4BZall BastarBare lands18,192.0010,500.000.000.0034,154.4052,346.40
ZB3BZall BastarBare lands18,192.0010,500.000.000.0034,154.4052,346.40
ZB2BZall BastarBare lands18,192.0010,500.000.000.0034,154.4052,346.40
ZB1BZall BastarBare lands18,192.0010,500.000.000.0034,154.4052,346.40
M1BMner i SipermBare earth with renaturalization18,192.0010,500.000.000.000.0018,192.00
M2FMner i SipermBroadleaved forests0.00 0.00 0.000.000.000.00
M1OMner i SipermOrchard18,192.0010,500.00950.0011,400.0034,154.4063,746.40
M1PMner i SipermCultivated meadow (alfalfa)18,192.0010,500.00700.008400.0034,154.4060,746.40
M1FMner i SipermBroadleaved forests0.000.000.000.000.000.00
M1DFMner i SipermDegraded forest0.000.000.000.000.000.00
M2PMner i SipermCultivated meadow (alfalfa)18,192.0010,500.00700.008400.0034,154.4060,746.40
M3PMner i SipermPasture—meadow (alfalfa)18,192.0010,500.00700.008400.0034,154.4060,746.40
M4PMner i SipermPasture—meadow (alfalfa)18,192.0010,500.00700.008400.0034,154.4060,746.40
1 Obtained by summing the replacement cost, the income loss, and the forest maintenance prize for 12 years.
Table A5. Financial convenience in the transformation of agricultural land and pastures into forests (cost–benefit analysis and replacement cost).
Table A5. Financial convenience in the transformation of agricultural land and pastures into forests (cost–benefit analysis and replacement cost).
Sample Plot CodeLand Use Classes (Corine Land Cover 2012)Total Tons of CO2 (EUR)Monetary Revaluation of Habitat Quality 2007–2018
(EUR)
Agric. and Forestry Production
(EUR/ha)
K
(EUR/ha)
P
(EUR/ha)
CES
(EUR/ha)
Total ES 1 (EUR/ha)Total
for 12
Years (EUR/ha)
Replacement Cost
(EUR/ha)
Lost Income
(EUR/ha)
Lost Income for 12 Years
(EUR/ha)
Maintenance
(EUR/ha)
Total Cost 2
(EUR/ha)
Economic Convenience (EUR/ha)
ZB5PPasture32.072982.09326.00 249.76 3.66 5.57 3509.1610,230.01 18,192.000.00 0.00 34,154.4052,346.4042,116.39
ZB4PPasture24.112930.39326.00 188.80 2.77 5.51 3477.58 9430.5318,192.000.00 0.00 34,154.4052,346.4042,915.87
ZB3PPasture32.132923.44326.00 246.36 3.61 5.423537.0610,220.64 18,192.00 0.00 0.00 34,154.4052,346.4042,125.76
ZB4OOrchard348.940.0035,320.501641.19224.222.2037,540.37450,418.10 18,192.00 950.0011,400.00 34,154.4063,746.40−386,671.70
ZB2PPasture57.382931.38326.00444.616.525.453771.3412,945.55 18,192.00 0.000.0034,154.4052,346.4039,400.85
ZB4FConifer forests64.110.002600.88201.933.342.452872.7034,443.06 0.00 0.00 0.00 0.00 0.00 no
ZB3FConifer forests0.952377.922600.88245.924.076.00 5235.7436,599.82 0.00 0.00 0.00 0.00 0.00 no
ZB2FOak forest with pasture0.342203.63326.00234.04 3.43 5.462772.908969.40 0.00 0.00 0.00 0.00 0.00 no
ZB1PConifer forests0.912524.452701.08203.463.376.395439.6637,430.27 0.00 0.00 0.00 0.00 0.00 no
ZB3AAgric. land (corn)0.870.003948.00612.8621.723.11 4586.5655,001.41 18,192.00 700.008400.0034,154.4060,746.405744.99
ZB2AAgric. land (corn and beans)0.880.005757.50 990.7135.11 2.206786.4081,410.36 18,192.00 700.00 8400.0034,154.4060,746.40−20,663.96
ZB1AForage0.900.002086.001086.9238.52 2.203214.5438,548.04 18,192.0064.00768.0034,154.4053,114.4014,566.36
ZB4BBare lands0.060.000.000.00 0.005.895.95 0.68 18,192.000.000.0034,154.4052,346.4052,345.72
ZB3BBare lands0.040.000.000.000.002.492.53 0.42 18,192.000.000.0034,154.4052,346.4052,345.98
ZB2BBare lands0.030.000.000.000.002.452.48 0.38 18,192.000.000.0034,154.4052,346.4052,346.02
ZB1BBare lands0.010.000.000.000.006.396.40 0.15 18,192.000.000.0034,154.4052,346.4052,346.25
M1BBare earth with renaturalization0.230.000.00534.155.401.47541.246477.26 18,192.000.000.000.0018,192.0011,714.74
M2FBroadleaved forests0.400.001147.9314.0211.221.481175.05 14,082.79 0.00 0.000.000.000.00no
M1OOrchard0.665156.4432,843.330.000.005.6038,005.73 399,283.96 18,192.00950.0011,400.0034,154.4063,746.40−335,537.56
M1PCultivated meadow (alfalfa)0.303867.062114.00635.317.745.646630.0536,955.30 18,192.00700.008400.0034,154.4060,746.4023,791.10
M1FBroadleaved forests0.402900.631147.93251.63 1.215.674307.47 19,714.71 0.000.000.000.000.00no
M1DFDegraded forest0.392905.11653.520.000.005.673564.6910,752.03 0.000.000.000.000.00no
M2PCultivated meadow (alfalfa)0.453480.832114.00220.712.695.676184.36 31,895.10 18,192.00700.008400.0034,154.4060,746.4028,851.30
M3PPasture—meadow (alfalfa)0.443863.92770.00223.712.735.604866.40 15,826.51 18,192.00700.008400.0034,154.4060,746.4044,919.89
M4PPasture—meadow (alfalfa)0.453862.22770.00310.863.795.604952.92 16,883.41 18,192.00700.008400.0034,154.4060,746.4043,862.99
1 Obtained by summing the values of CO2 (EUR/t), the monetary revaluation of habitat quality (EUR/ha), agricultural and forestry production (EUR/ha), K, P (EUR/ha), and CES (sustainable tourism suitability) (EUR/ha). 2 Obtained by summing the replacement cost, the income loss, and the maintenance prize for 12 years.

References

  1. Blum, A. Drought Resistance, Water-Use Efficiency, and Yield Potential—Are They Compatible, Dissonant, or Mutually Exclusive? Aust. J. Agric. Res. 2005, 56, 1159. [Google Scholar] [CrossRef]
  2. Blum, W.E.H. Functions of Soil for Society and the Environment. Rev. Environ. Sci. Biotechnol. 2005, 4, 75–79. [Google Scholar] [CrossRef]
  3. Costanza, R.; d’Arge, R.; De Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neill, R.V.; Paruelo, J.; et al. The Value of the World’s Ecosystem Services and Natural Capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  4. Leemans, R.; De Groot, R.S. Millennium Ecosystem Assessment, Ecosystem and Human Well-Being: A Framework for Assessment; Island Press: Washington, DC, USA, 2005. [Google Scholar]
  5. European Commission. Thematic Strategy for Soil Protection; European Commission: Brussels, Belgium, 2006. [Google Scholar]
  6. Steinhoff-Knopp, B.; Kuhn, T.K.; Burkhard, B. The Impact of Soil Erosion on Soil-Related Ecosystem Services: Development and Testing a Scenario-Based Assessment Approach. Environ. Monit. Assess. 2021, 193, 274. [Google Scholar] [CrossRef]
  7. Guerra, C.A.; Maes, J.; Geijzendorffer, I.; Metzger, M.J. An Assessment of Soil Erosion Prevention by Vegetation in Mediterranean Europe: Current Trends of Ecosystem Service Provision. Ecol. Indic. 2016, 60, 213–222. [Google Scholar] [CrossRef]
  8. Sartori, M.; Philippidis, G.; Ferrari, E.; Borrelli, P.; Lugato, E.; Montanarella, L.; Panagos, P. A Linkage between the Biophysical and the Economic: Assessing the Global Market Impacts of Soil Erosion. Land Use Policy 2019, 86, 299–312. [Google Scholar] [CrossRef]
  9. Panagos, P.; Katsoyiannis, A. Soil Erosion Modelling: The New Challenges as the Result of Policy Developments in Europe. Environ. Res. 2019, 172, 470–474. [Google Scholar] [CrossRef]
  10. EU Soil Strategy for 2030—European Commission. Available online: https://environment.ec.europa.eu/publications/eu-soil-strategy-2030_en (accessed on 10 September 2024).
  11. Barrios, E. Soil Biota, Ecosystem Services and Land Productivity. Ecol. Econ. 2007, 64, 269–285. [Google Scholar] [CrossRef]
  12. Brevik, E.C.; Calzolari, C.; Miller, B.A.; Pereira, P.; Kabala, C.; Baumgarten, A.; Jordán, A. Soil Mapping, Classification, and Pedologic Modeling: History and Future Directions. Geoderma 2016, 264, 256–274. [Google Scholar] [CrossRef]
  13. Adhikari, K.; Hartemink, A.E. Linking Soils to Ecosystem Services—A Global Review. Geoderma 2016, 262, 101–111. [Google Scholar] [CrossRef]
  14. Balbi, S.; Prado, A.D.; Gallejones, P.; Geevan, C.P.; Pardo, G.; Pérez-Miñana, E.; Manrique, R.; Hernandez-Santiago, C.; Villa, F. Modeling Trade-Offs among Ecosystem Services in Agricultural Production Systems. Environ. Model. Softw. 2015, 72, 314–326. [Google Scholar] [CrossRef]
  15. FAO. The Future of Food and Agriculture—Trends and Challenges; Food and Agriculture Organization of the United Nations: Rome, Italy, 2017. [Google Scholar]
  16. Pereira, P.; Brevik, E.C.; Muñoz-Rojas, M.; Miller, B.A.; Smetanova, A.; Depellegrin, D.; Misiune, I.; Novara, A.; Cerdà, A. Soil Mapping and Processes Modeling for Sustainable Land Management. In Soil Mapping and Process Modeling for Sustainable Land Use Management; Elsevier: Amsterdam, The Netherlands, 2017; pp. 29–60. ISBN 978-0-12-805200-6. [Google Scholar]
  17. Robinson, D.A.; Fraser, I.; Dominati, E.J.; Davíðsdóttir, B.; Jónsson, J.O.G.; Jones, L.; Jones, S.B.; Tuller, M.; Lebron, I.; Bristow, K.L.; et al. On the Value of Soil Resources in the Context of Natural Capital and Ecosystem Service Delivery. Soil Sci. Soc. Am. J. 2014, 78, 685–700. [Google Scholar] [CrossRef]
  18. Ricci, G.F.; Jeong, J.; De Girolamo, A.M.; Gentile, F. Effectiveness and Feasibility of Different Management Practices to Reduce Soil Erosion in an Agricultural Watershed. Land Use Policy 2020, 90, 104306. [Google Scholar] [CrossRef]
  19. Ricci, G.F.; D’Ambrosio, E.; De Girolamo, A.M.; Gentile, F. Efficiency and Feasibility of Best Management Practices to Reduce Nutrient Loads in an Agricultural River Basin. Agric. Water Manag. 2022, 259, 107241. [Google Scholar] [CrossRef]
  20. Turner, R.K.; Daily, G.C. The Ecosystem Services Framework and Natural Capital Conservation. Environ. Resour. Econ. 2008, 39, 25–35. [Google Scholar] [CrossRef]
  21. Francis, R.A. Peter Kareiva, Heather Tallis, Taylor H. Ricketts, Gretchen C. Daily, and Stephen Polasky (eds), Natural Capital: Theory and Practice of Mapping Ecosystem Services. Oxford: Oxford University Press, 2011; 365 pp.: 9780199588992, £75 (hbk). Prog. Phys. Geogr. Earth Environ. 2011, 35, 701–704. [Google Scholar] [CrossRef]
  22. Kragt, M.E.; Robertson, M.J. Quantifying Ecosystem Services Trade-Offs from Agricultural Practices. Ecol. Econ. 2014, 102, 147–157. [Google Scholar] [CrossRef]
  23. Pascual, U.; Termansen, M.; Hedlund, K.; Brussaard, L.; Faber, J.H.; Foudi, S.; Lemanceau, P.; Jørgensen, S.L. On the Value of Soil Biodiversity and Ecosystem Services. Ecosyst. Serv. 2015, 15, 11–18. [Google Scholar] [CrossRef]
  24. De Groot, R.S.; Wilson, M.A.; Boumans, R.M.J. A Typology for the Classification, Description and Valuation of Ecosystem Functions, Goods and Services. Ecol. Econ. 2002, 41, 393–408. [Google Scholar] [CrossRef]
  25. Kumar, P. (Ed). The Economics of Ecosystems and Biodiversity: Ecological and Economic Foundations; Routledge: London, UK, 2012; ISBN 978-1-136-53880-3. [Google Scholar]
  26. Susan, S.; Karlen, D.L.; Cambardella, C.A. The Soil Management Assessment Framework: A Quantitative Soil Quality Evaluation Method. Soil Sci. Soc. Am. J. 2004, 68, 1945–1962. [Google Scholar] [CrossRef]
  27. Dominati, E.; Patterson, M.; Mackay, A. A Framework for Classifying and Quantifying the Natural Capital and Ecosystem Services of Soils. Ecol. Econ. 2010, 69, 1858–1868. [Google Scholar] [CrossRef]
  28. Faber, J.H.; Van Wensem, J. Elaborations on the Use of the Ecosystem Services Concept for Application in Ecological Risk Assessment for Soils. Sci. Total Environ. 2012, 415, 3–8. [Google Scholar] [CrossRef] [PubMed]
  29. Haygarth, P.M.; Ritz, K. The Future of Soils and Land Use in the UK: Soil Systems for the Provision of Land-Based Ecosystem Services. Land Use Policy 2009, 26, S187–S197. [Google Scholar] [CrossRef]
  30. Lavelle, P.; Decaëns, T.; Aubert, M.; Barot, S.; Blouin, M.; Bureau, F.; Margerie, P.; Mora, P.; Rossi, J.-P. Soil Invertebrates and Ecosystem Services. Eur. J. Soil Biol. 2006, 42, S3–S15. [Google Scholar] [CrossRef]
  31. Palm, C.; Sanchez, P.; Ahamed, S.; Awiti, A. Soils: A Contemporary Perspective. Annu. Rev. Environ. Resour. 2007, 32, 99–129. [Google Scholar] [CrossRef]
  32. Smit, E.; Bakker, P.A.H.M.; Bergmans, H.; Bloem, J.; Griffiths, B.S.; Rutgers, M.; Sanvido, O.; Singh, B.K.; Van Veen, H.; Wilhelm, R.; et al. General Surveillance of the Soil Ecosystem: An Approach to Monitoring Unexpected Adverse Effects of GMO’s. Ecol. Indic. 2012, 14, 107–113. [Google Scholar] [CrossRef]
  33. Wall, D.H.; Bardgett, R.; Covich, A.; Snelgrove, P. The Need for Understanding How Biodiversity and Ecosystem Functioning Affect Ecosystem Services in Soils and Sediments. Sustaining Biodiversity and Ecosystem Services in Soils and Sediments; Wall, D.H., Ed.; Island Press: Washington, DC, USA, 2004; Volume 64. [Google Scholar]
  34. Warr, B.; Ayres, R. Accounting for Soils: Towards an Integrated Sustainability and Productivity Assessment for Soils; Centre for the Management of Environmental and Social Responsibility (CMER): Fontainbleau, France, 2004. [Google Scholar]
  35. Weber, J.L. Accounting for Soil in the SEEA; European Environment Agency: Rome, Italy, 2007. [Google Scholar]
  36. Porter, J.; Costanza, R.; Sandhu, H.; Sigsgaard, L.; Wratten, S. The Value of Producing Food, Energy, and Ecosystem Services within an Agro-Ecosystem. AMBIO J. Hum. Environ. 2009, 38, 186–193. [Google Scholar] [CrossRef]
  37. Sandhu, H.S.; Wratten, S.D.; Cullen, R. Organic Agriculture and Ecosystem Services. Environ. Sci. Policy 2010, 13, 1–7. [Google Scholar] [CrossRef]
  38. Sandhu, H.S.; Wratten, S.D.; Cullen, R. The Role of Supporting Ecosystem Services in Conventional and Organic Arable Farmland. Ecol. Complex. 2010, 7, 302–310. [Google Scholar] [CrossRef]
  39. Swinton, S.M.; Lupi, F.; Robertson, G.P.; Hamilton, S.K. Ecosystem Services and Agriculture: Cultivating Agricultural Ecosystems for Diverse Benefits. Ecol. Econ. 2007, 64, 245–252. [Google Scholar] [CrossRef]
  40. Zhang, W.; Ricketts, T.H.; Kremen, C.; Carney, K.; Swinton, S.M. Ecosystem Services and Dis-Services to Agriculture. Ecol. Econ. 2007, 64, 253–260. [Google Scholar] [CrossRef]
  41. Baloshi, V.; Gjoka, V.; Collaku, N.; Toromani, E. The Possibility of Establishing a Scheme for Payments for Ecosystem Services at the Bovilla Watershed (Tirana). Albanian J. Agric. Sci. 2019, 18, 50–57. [Google Scholar]
  42. Wunder, S. Payments for Environmental Services: Some Nuts and Bolts; Center for International Forestry Research (CIFOR): Bogor, Indonesia, 2005. [Google Scholar]
  43. Van der Born, G.J.; De Haan, B.J.; Pearce, D.W.; Howarth, A. Technical Report on Soil Degradation in Europe: An Integrated Economic and Environmental Assessment; RIVM: Bilthoven, The Netherlands, 2000. [Google Scholar]
  44. Baloshi, V.; Gjoka, F.; Çollaku, N.; Toromani, E. Determination of Soil Loss by Erosion in Different Land Covers Categories and Slope Classes in Bovilla Watershed, Tirana, Albania. Int. J. Environ. Ecol. Eng. 2019, 13, 57–61. [Google Scholar] [CrossRef]
  45. Baloshi, V.; Gjoka, V.; Collaku, N.; Toromani, E. Relationship between Water Quality and Land Use in Bovilla Watershed, Tirana Region (Albania). Fresenius Environ. Bull. 2019, 28, 4435–4441. [Google Scholar]
  46. ISO 14235:1998; Soil Quality—Determination of Organic Carbon by Sulfochromic Oxidation. ISO: Geneva, Switzerland, 1998. Available online: https://www.iso.org/standard/23140.html (accessed on 20 June 2024).
  47. Prezzi CO2—Sendeco2. Available online: https://www.sendeco2.com/it/prezzi-co2 (accessed on 20 June 2024).
  48. Investing.com—Stock Market Quotes & Financial News. Available online: https://www.investing.com/ (accessed on 20 June 2024).
  49. Habitat Quality|The Natural Capital Project. Available online: https://naturalcapitalproject.stanford.edu/invest/habitat-quality (accessed on 20 June 2024).
  50. Hall, L.S.; Krausan, P.R.; Morrison, M.L. The Habitat Concept and a Plea for Standard Terminology. Wildl. Soc. Bull. 1997, 25, 173–182. [Google Scholar]
  51. Nellemann, C.; Thomsen, M.G. Long-Term Changes in Forest Growth: Potential Effects of Nitrogen Deposition and Acidification. Water Air Soil Pollut. 2001, 128, 197–205. [Google Scholar] [CrossRef]
  52. McKinney, M.L. Urbanization, Biodiversity, and Conservation. BioScience 2002, 52, 883. [Google Scholar] [CrossRef]
  53. Forman, R.T.T. (Ed.) Road Ecology: Science and Solutions; Island Press: Washington, DC, USA, 2003; ISBN 978-1-55963-932-3. [Google Scholar]
  54. Costanza, R.; De Groot, R.; Sutton, P.; Van Der Ploeg, S.; Anderson, S.J.; Kubiszewski, I.; Farber, S.; Turner, R.K. Changes in the Global Value of Ecosystem Services. Glob. Environ. Chang. 2014, 26, 152–158. [Google Scholar] [CrossRef]
  55. Conversion Rates—Purchasing Power Parities (PPP)—OECD Data. Available online: https://data.oecd.org/conversion/purchasing-power-parities-ppp.htm (accessed on 20 June 2024).
  56. Istat—Rivaluta. Available online: https://rivaluta.istat.it/ (accessed on 20 June 2024).
  57. Bando Sottomisura 8.1—2022—PSR Puglia. Available online: https://psr.regione.puglia.it/bando-sottomisura-8.1-2022?redirect=%2Fsottomisura-8.1 (accessed on 22 July 2024).
  58. Elenco Prezzi Opere Pubbliche. 2024. Available online: https://www.regione.puglia.it/web/territorio-mobilita-e-infrastrutture/elenco-prezzi-opere-pubbliche-2024 (accessed on 22 July 2024).
  59. Masaf—Risultato Della Ricerca. Available online: https://www.politicheagricole.it/flex/cm/FixedPages/Common/Search.v3.php/L/IT/s/1 (accessed on 22 July 2024).
  60. Madenoğlu, S.; Pınar, M.Ö.; Şahin, S.; Erpul, G. Sustainable Land Management for Mitigating Soil Erosion at the Catchment Scale. Türkiye Tarımsal Araştırmalar Derg. 2024, 11, 176–190. [Google Scholar] [CrossRef]
Figure 1. Location and physical map of the Bovilla watershed.
Figure 1. Location and physical map of the Bovilla watershed.
Land 13 01512 g001
Figure 2. Economic evaluation methods of ecosystem services (authors’ elaboration of [42,43,44]).
Figure 2. Economic evaluation methods of ecosystem services (authors’ elaboration of [42,43,44]).
Land 13 01512 g002
Figure 3. Habitat quality.
Figure 3. Habitat quality.
Land 13 01512 g003
Figure 4. Sustainable tourism suitability (STS, score range 0–4 = road network density map score range 0–1 + network of paths density map score range 0–1 + habitat quality map score range 0–1 + point of interest density map score range 0–1).
Figure 4. Sustainable tourism suitability (STS, score range 0–4 = road network density map score range 0–1 + network of paths density map score range 0–1 + habitat quality map score range 0–1 + point of interest density map score range 0–1).
Land 13 01512 g004
Figure 5. Soil loss (t/ha) by sample plots over the monitoring period.
Figure 5. Soil loss (t/ha) by sample plots over the monitoring period.
Land 13 01512 g005
Table 1. Soil ES classification quoted from the Millenium Ecosystem Assessment. (Source: authors’ elaboration of [4]).
Table 1. Soil ES classification quoted from the Millenium Ecosystem Assessment. (Source: authors’ elaboration of [4]).
Description and MeaningCategory
Products people obtain from ecosystems, such as food, fuel, fibers, and fresh water.Provisioning services
The benefits people obtain from the regulation of ecosystem processes, including air quality maintenance, climate regulation, erosion control, regulation of human diseases, and water purification.Regulating services
The non-material benefits people obtain from ecosystems through spiritual enrichment, cognitive development, reflection, recreation, and aesthetic experiences.Cultural services
Services that are necessary for the production of all other ESs, such as primary production, production of oxygen, and soil formation.Supporting functions
Table 2. Description of each sub-watershed in the study area.
Table 2. Description of each sub-watershed in the study area.
Sub-WatershedArea (ha)Slope (degree)Elevation (m) a.s.l.
Min Max
Bovilla reservoir345.781.70314.60356.29
Ranxe Bruz-Zall1626.2620.45314.601249.95
Bruz Mal1422.9522.99319.821623.64
Mner i Siperm967.4116.86314.601180.65
Vilez2343.9320.00321.571520.81
Zall Bastar1979.4218.11335.581606.49
Total9515.24
Table 3. Chemical composition of soils in the study area.
Table 3. Chemical composition of soils in the study area.
Sample
Code
Soil
Type
Laboratory Analysis
pH H2OpH CaCl2K AvailableP AvailableOrganic Matter (%)SandSiltClay
mg/100 g %
M1BLeptosol8.167.389.500.121.2815.7649.9032.10
ZB1AKambisol7.606.9522.801.015.0612.4665.6621.89
ZB3FLeptosol8.167.3914.500.305.3916.4252.7930.79
ZB3PKambisol6.665.8613.100.251.9318.8052.9128.90
ZB4OKambisol7.426.9021.903.745.2616.7051.6527.43
M1FPhaeozem7.767.0711.600.072.2721.1751.8626.97
M2PKambisol7.837.0610.500.162.5716.0348.6135.35
M2FKambisol8.137.3980.951.262.2719.2149.2233.30
Table 4. Land use classes in Bovilla watershed.
Table 4. Land use classes in Bovilla watershed.
Land Use Classes (Corine Land Cover 2012)Area (ha)
Pasture302.531
Complex cultivated areas 453.654
Lands covered with crops and natural vegetation1247.31
Broadleaved forests2423.62
Conifer forests 149.153
Mixed forests33.012
Herbaceous vegetation517.519
Sclerophyll vegetation 1143.24
Transition areas from forests to shrub lands1964.31
Others130.525
Bare lands0.00
Sparse vegetation areas656.366
Aquatic areas (rivers, streams, lakes)436.144
Total area9457.384
Table 5. Methods and standards used for chemical compound analysis.
Table 5. Methods and standards used for chemical compound analysis.
Territorial ElementDescription
Road networkAll roads in the park
Network of pathsThe network of pedestrian roads and hiking trails
Habitat qualityThe spatial distribution of elements with high naturalistic value within the park
Points of interestRestaurants, bars, hotels, sports facilities, and historical–cultural elements
Table 6. Sum of soil loss (t/ha) by land use and year.
Table 6. Sum of soil loss (t/ha) by land use and year.
Land UseTotal Soil Loss by Year (t/ha)
201720182019
Agricultural land57.5990.3729.05
Bare land773.78747.47283.61
Degraded forest25.2735.2411.53
Forest33.1363.0326.20
Orchard81.6379.3019.36
Pasture154.81136.0958.15
Total1126.211151.50427.90
Table 7. Mean soil loss (t/ha) by land use and year.
Table 7. Mean soil loss (t/ha) by land use and year.
Land UseMean Soil Loss by Year (t/ha)Mean Soil Loss (2017–2019)
t/ha
201720182019
Agricultural land19.2030.129.6819.67
Bare land154.76149.4956.72120.32
Degraded forest25.2735.2411.5324.01
Forest6.6312.615.248.16
Orchard40.8239.659.6830.05
Pasture17.2015.126.4612.93
Total45.0446.0617.1236.07
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

Romanazzi, G.R.; Ottomano Palmisano, G.; Cioffi, M.; Leronni, V.; Toromani, E.; Koto, R.; De Boni, A.; Acciani, C.; Roma, R. A Cost–Benefit Analysis for the Economic Evaluation of Ecosystem Services Lost Due to Erosion in a Mediterranean River Basin. Land 2024, 13, 1512. https://doi.org/10.3390/land13091512

AMA Style

Romanazzi GR, Ottomano Palmisano G, Cioffi M, Leronni V, Toromani E, Koto R, De Boni A, Acciani C, Roma R. A Cost–Benefit Analysis for the Economic Evaluation of Ecosystem Services Lost Due to Erosion in a Mediterranean River Basin. Land. 2024; 13(9):1512. https://doi.org/10.3390/land13091512

Chicago/Turabian Style

Romanazzi, Giuliano Rocco, Giovanni Ottomano Palmisano, Marilisa Cioffi, Vincenzo Leronni, Ervin Toromani, Romina Koto, Annalisa De Boni, Claudio Acciani, and Rocco Roma. 2024. "A Cost–Benefit Analysis for the Economic Evaluation of Ecosystem Services Lost Due to Erosion in a Mediterranean River Basin" Land 13, no. 9: 1512. https://doi.org/10.3390/land13091512

APA Style

Romanazzi, G. R., Ottomano Palmisano, G., Cioffi, M., Leronni, V., Toromani, E., Koto, R., De Boni, A., Acciani, C., & Roma, R. (2024). A Cost–Benefit Analysis for the Economic Evaluation of Ecosystem Services Lost Due to Erosion in a Mediterranean River Basin. Land, 13(9), 1512. https://doi.org/10.3390/land13091512

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