Next Article in Journal
Solution of Multi-Crew Depots Railway Crew Scheduling Problems: The Chinese High-Speed Railway Case
Previous Article in Journal
Positive Organizational Practices, Life Satisfaction, and Psychological Capital in the Public and Private Sectors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Economic Evaluation of Posidonia oceanica Ecosystem Services along the Italian Coast

1
Centro Euro-Mediterraneo sui Cambiamenti Climatici, Ocean Predictions and Applications Division, 73100 Lecce, Italy
2
Laboratory of Experimental Oceanology end Marine Ecology, Department of Ecological and Biological Sciences (DEB), Università degli Studi della Tuscia, 00053 Civitavecchia, Italy
3
Italian National Institute for Environmental Protection and Research (ISPRA), 00144 Rome, Italy
4
Department of Forest Environment and Resources, University of Tuscia, 01100 Viterbo, Italy
5
Port Authority System of the Central Northern Tyrrhenian Sea, 00053 Civitavecchia, Italy
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(1), 489; https://doi.org/10.3390/su14010489
Submission received: 10 November 2021 / Revised: 9 December 2021 / Accepted: 29 December 2021 / Published: 3 January 2022

Abstract

:
This study reports the quantification and analysis of the ecosystem services (ESs) value of Posidonia oceanica (L.) Delile (1813) meadows in Italy (Liguria, Tuscany, Latium, Campania, Calabria, Apulia, Sardinia, and Sicily regions). The ES evaluation method of P. oceanica meadows applied in this study was obtained from a previous study in which the site-specific approach for the definition of the benefits and services was applied. The distribution of P. oceanica and the ESs economic value have been managed through an open-source geographic information system (QGIS), focusing on five essential ESs: carbon sequestration, bioremediation, oxygen production, erosion protection, and food production. The average value of the ESs obtained on the Italian national scale is €21,660.5 ha−1 yr−1, which is comparable with the values reported in the international literature concerning P. oceanica ESs. The results of this study confirm that the economic evaluation of ESs is an essential tool for the management of the coastal marine environment, especially considering the modularity of the applied approach. The value of the total benefits, considering the entire extension of the national P. oceanica meadows, represents significant value with respect to the Italian gross domestic product, and the individual budgets of the considered Italian regions.

1. Introduction

About 30% of the Italian population lives along the coast, which has an extension of about 7500 km [1]. The strong anthropogenic pressure favors the highly vulnerable conditions of the coastal marine environment, which is characterized by multiple physical and ecological conditions that can support multiple natural resources uses. In this context, conflicts between various interests that affect both human uses and natural resources arise [2]. To address this issue, and minimize the mutual impact of coastal uses, management plans must include the conservation of ecological structures and processes involved in these uses, in addition to evaluating their usefulness in human activities [2].
The conservation status of a given ecosystem affects both the quality and quantity of its ecosystem services (ESs), which are defined as “The benefits people obtain from ecosystems. These include provisioning services such as food, water, timber, and fibre; regulating services that affect climate, floods, disease, wastes, and water quality; cultural services that provide recreational, aesthetic, and spiritual benefits; and supporting services such as soil formation, photosynthesis, and nutrient cycling” [3]. Originally, this concept considered only the services provided by nature, but over time it has also come to include socio-economic concepts and conservation objectives [4]. The assignment of economic value to ESs allows access to a broader information framework, which allows better management choices of resource allocation.
Italy is a hot spot of marine biodiversity at the habitat level. Understanding the habitat and population distributions along the Italian coasts, despite having increased considerably in recent years because of large-scale surveys, literature reviews, and expert consultations, is still inadequate for natural capital assessment and environmental accounting [5].
Among the Italian marine and coastal ecosystems, and in general, in the Mediterranean Sea basin, P. oceanica meadows are key coastal habitats, hosting about 25% of all species present in the Mediterranean Sea [6]. Nevertheless, extensive regression of P. oceanica meadows has been recorded across the Mediterranean Sea due to climate change and human activities. In Italy, a total regressed area of 34,472 ha has been documented along the continental coasts of Liguria, Tuscany, Latium, and Apulia regions over the last 20–30 years. In Sardinia, large areas of dead matte, amounting to 23,215 ha, were reported in the Gulfs of Cagliari, Olbia, and Asinara [7]. The assessment of the distribution and density of P. oceanica meadows can therefore be used as a tool both for coastal planning programs, and conservation and monitoring programs. ESs estimates of P. oceanica meadows provides fundamental information for sustainable management of the related resources.
Several examples of the application of the concept of natural capital and the estimation of ESs for the P. oceanica meadows are available. A recent and accurate compendium is the technical report of the European Union prepared by the United Nations Environment Program World Conservation Monitoring Centre (UNEP-WCMC) in collaboration with the Institute for European Environmental Policy (IEEP) and the Centre for Social and Economic Research on the Global Environment (CSERGE) at the University of East Anglia [8]. Furthermore, the works of Campagne et al. [9] and Vassallo et al. [10], together with those of Costanza and Marcelli [2,11], which are used as a starting point for the evaluation of P. oceanica ESs in this work, describe different methodological approaches, and are particularly useful for the comparative analysis in terms of the economic estimate of the services considered. Campagne et al. [9] provided an extensive review of the existing economic assessments on the benefits considered for P. oceanica on the French national territory, using the “transfer of value” method and, in some cases, with site-specific economic assessments. Vassallo et al. [10] proposed an approach called “Emergy to Money Ratio” in which the flow of matter and energy from various inputs into the ecosystem (in this specific case, P. oceanica) is estimated and regulated by the necessary solar energy required to maintain this process.
Hence, this work aimed to analyze and quantify the ESs economic value of the Italian P. oceanica meadows, using synoptic representations through an open-source geographic information system (QGIS) to evaluate the intervention and alternatives of investment more effectively in cases of environmental impact assessment, and facilitate the process of defining and choosing action priorities.

2. Materials and Methods

2.1. Study Area and Data Organization

The study area included the Italian coastal zones on which P. oceanica forms large and extensive meadows. For this reason, the small patches located in the Friuli region (northeast of Italy) are not considered in this work. The regions under study included Liguria, Tuscany, Latium, Campania, Calabria, Apulia, Sardinia, and Sicily. For each region, both the distribution of P. oceanica meadows present inside the sites of community importance and marine special protection areas (SCIs and SPAs, respectively), and the distributions of P. oceanica meadows outside the limits of the SCIs and marine SPAs were considered.
To evaluate the ecosystem services of P. oceanica meadows along the Italian coasts, a dataset that included 356 measurement sites was used (Figure 1). Each measurement site reported several parameters: (1) point ID; (2) site code; (3) geographical location of the sampling point; (4) region; (5) geographical coordinates; (6) absolute density of the meadow in m2; (7) standard deviation of the absolute density; (8) relative density (m2); (9) standard deviation of relative density; (10) meadows coverage (%); (11) depth (m); (12) description of the lower limit; (13) the average number of shoots per m−2; (14) average leaf area; (15) leaf area index (LAI); (16) average number of leaves; (17) average length of the rhizome; (18) foliar production (gDW yr−1); (19) average rhizome production (gDW yr−1); (20) sampling year; and (21) data source. These data were collected by the Italian Institute for Environmental Protection and Research (ISPRA) [1], as it is necessary for the European (EU) reports relating to the initial assessment and definition of the Good Environmental Status (GES), and identification of the Environmental Targets (articles 8–10 of the Marine Strategy Framework Directive [MSFD] 2008/56/EC; http://www.db-strategiamarina.isprambiente.it/app/#/, accessed on 28 December 2021) [12].
The P. oceanica meadows taken into consideration had a total extension of 358,896.5 ha. In detail, among the eight regions considered, the largest regional extension was reached in Sardinia (174,489.7 ha), whereas the smallest was in Liguria (5050.1 ha), as shown in Table 1. The extension of Posidonia oceanica meadows (ha) in the Italian regions under consideration reported in Table 1 was obtained by QGIS elaboration of the data collected by ISPRA.
To calculate the P. oceanica ESs on the Italian territory, the dataset was preliminarily analyzed to derive the distribution of the P. oceanica meadows for each Italian region, the spatial distribution of the different types of substrates on which the meadows insist, and the extent of the SCIs.
The regional context coincides with the administrative subdivision of the maritime districts, and is the most suitable for analyzing the distributions of marine SCIs and SPAs, as these are of regional competence.
The data from the 356 measuring sites located in the Italian regions were organized using an open-source geographic information system (QGIS), which was previously prepared considering the shapefiles of the SCI areas.

2.2. Estimation of P. oceanica Meadows ESs

The method applied for the evaluation of P. oceanica ESs was obtained from a report by [11] in which a site-specific approach for the definition of the benefits and services was applied [2]. This approach was used both to identify the P. oceanica meadows ESs, for the definition of the methods, and for the estimation of each single service. Five benefits were considered in this work for P. oceanica meadows: (1) carbon sequestration; (2) oxygen production; (3) bioremediation; (4) protection from erosion; (5) food production [2]. In detail, carbon sequestration was computed considering mean regional (Maritime Districts) values of biomass and relative density data; oxygen production was calculated considering the regional mean values of biomass and the number of shoots per m2; bioremediation was computed considering the mean regional values of the number of shoots per m2; protection from erosion was computed considering wave height, period, and direction on different points at regional level; and food production was computed considering the annual revenues of professional fishing for each Maritime district. The transfer of value was performed for each ES after considering what reported by Marcelli et al. [2], except for bioremediation, for which the operating cost of a wastewater plant was used. The methods for the economic evaluation of benefits are summarized in Table 2. The results of the estimate of the economic value of P. oceanica benefits are expressed in € ha−1 yr−1.
For the calculation of the values of some benefits, such as carbon sequestration, erosion protection, and food production, the spatial distribution of the different types of substrates on which the meadow stands were considered. Regarding the substrates, dead matte, degraded P. oceanica, P. oceanica on matte, P. oceanica on rock, P. oceanica on sand, and unspecified substrates were considered.

2.2.1. Carbon Sequestration

With respect to the method proposed by Marcelli et al. [2], litter and soil compartments of the P. oceanica meadows were also considered in the estimation of carbon sequestration value following the approach reported in the InVEST Coastal Blue Carbon model (Natural Capital Project, Stanford University; https://naturalcapitalproject.stanford.edu/, accessed on 9 November 2021) [13].
The soil includes rhizomes/dead roots and organic sediment (matte), and represents the main pool of organic carbon stocked by P. oceanica meadows. The carbon sequestered in this compartment (CS) depends on the matte thickness, as its distribution decreases exponentially with the depth [14]. Taking into account a reference value of about 79 kg Corg m−2 in about 5 m of core [14], and assuming a different average thickness based on the type of substrate (we have arbitrarily assigned 0.5 and 2 m thick for the matte on rock and sand, respectively), it is possible to roughly estimate a quantity of carbon equal to 39.5 and 53.44 kg Corg m−2 for Posidonia on rock and for the plants located on the sand, respectively.
Litter means that the leaves that detach from the plants become stranded on the coast, and form massive supra-littoral deposits called banquettes. In this case, the dry weight values (BL) of 102.126 g m−2 and 22.749 g m−2 for high- and low-density meadows, respectively, were used [15]. The shoot density limit value is 300 shoots m−2, which corresponds to a “sparse meadow” according to the classification proposed by Giraud [16].
The carbon stock in the above- and below-ground biomass of P. oceanica was calculated using the foliar (Be) and rhizomal (Bi) biomass data, respectively [2].
To convert the dry weight determined in the biomass and litter compartments of P. oceanica meadow into organic carbon content, a common value (Fc) of 0.336 [17] was used.
The total organic carbon (CT) value stocked in the P. oceanica meadows is given by Equation (1) shown below:
C T = ( B e × F c ) + ( B i × F c ) + ( B L × F c ) + C S
For the transfer of value, Equation (2) was then used:
C = 0.01 × C T × 24.7
in which C is the value of the carbon sequestration benefit for P. oceanica meadows (€ ha−1 yr−1), 0.01 is the conversion factor from gm−2 to tonne ha−1, and 24.7 is the average annual exchange value of 1 tonne of CO2 based on the Emission Trading Scheme, updated to 2019 [18].

2.2.2. O2 Production

For the calculation of the O2 supply by P. oceanica, we referred to what was reported by Marcelli et al. [2]. The net oxygen production is estimated using literature data attesting the quantity of 5 mmol of O2 per gram of dry weight of the biomass produced daily during the metabolic cycles of P. oceanica [19]. Here, O2 supply estimation is further refined considering the relative rate attributable to the meadow-associated communities, which is about 0.018 mol m−2 day−1, given the condition of autotrophy found in the epiphytic communities during the year [20]. The economic value of O2 production of P. oceanica was therefore obtained using the cost of industrial methods to produce oxygen, about 0.05 € kg−1 [21], and expressed as a value per annual ha.
The total O2 supply is affected by the O2 supply attributable to the biomass (Ob) described by Equation (3), and by the oxygen supply attributable to the epiphytic community—Oe; Equation (4):
O b = [ B × ( M a × 32 ) ] × 10 C
O e = ( M c × 32 ) × 10 C
in which B is the meadow biomass (g m−2 yr−1), Ma is the oxygen moles supplied in a year as obtained from the literature data of 0.03 mol m−2 day−1 multiplied by 365 days after considering the net primary production [19], Mc is the moles supplied in a year by the epiphytic communities (equal to 0.018 mol m−2 day−1 multiplied by 365 days) [20], 32 is the O2 molecular mass, 10 the conversion factor of grams into kilograms and m2 into ha, and C is the O2 cost of industrial production per kg (equal to €0.05). Finally, the total O2 supply is shown in Equation (5):
O t o t = O b + O e

2.2.3. Bioremediation

Bioremediation refers to the capability to perform pollutant degradation and remediation through metabolic processes. Specifically, this work considers the capability of P. oceanica meadows to remove considerable quantities of nitrogen (N) and phosphorus (P) from coastal waters and sediments. P. oceanica plants and organism complexes associated with them carry out processes of excess N and P removal, thus allowing portions of each to become available again for coastal biological communities. The removal of excess loads of nutrients available for photosynthetic organisms causes limitations in the proliferation of algae and bacteria, among other factors, exercising control over the onset of eutrophication processes.
To calculate the economic value of this benefit, data extracted from recent studies relating to the costs for wastewater treatment systems from municipal plants in the United States of America were used (https://www.epa.gov/sites/production/files/2015-04/documents/nutrient-economics-report-2015.pdf, accessed on 9 December 2021) [22]. In regards to the estimation of the quantity of P and N used by the metabolic processes from P. oceanica, we refer to Marcelli et al. [2], using the available data of the number of shoots/m2, and the literature data of N and P concentration as the minimum daily quantity required by the plant [23]. In detail, the values selected in this situation correspond to 0.09 mg shoots day−1 of N, and 0.01 mg shoots day−1 for P after considering non-flowering meadows. The literature values were then compared with the values of the number of shoots per m2 of each available measurement site, and expressed on a regional scale as the mass of P and N ha yr−1. The relationship is described in Equation (6):
B = { [ E × ( n × 10 4 ) × A ] [ G × m ] } × C
in which B is the average annual N and P bioremediation value per ha of N expressed in €, E is the minimum daily amount of N or P required from the meadows (respectively 0.09 mg shoots−1 day−1 of N and 0.01 mg shoots−1 day−1 of P), n is the number of shoots per m2, A represents the days in one year (365), G is 3.78 (conversion factor from gallons to L), m is 15 for N and 1.9 for P (average milligrams contained in 1 L of water entering the urban treatment plant, obtained according to the efficiency of the treatment), and C is the cost in euros for the operation and management of an urban wastewater treatment plant (0.018 and €0.01 for N and P, respectively, per L). It should be specified that the cost in this case is related to the operation and annual maintenance of existing plants, thus not considering the capital investment cost for the construction of new urban waste treatment plants. The lowest operating and management cost was chosen considering information from the report by the United States Environmental Protection Agency [22] after comparison with the average operating management cost of Italian plants as reported by Sicuro [24], in which an average cost per m3 of sludge disposal in small- and medium-sized urban treatment plants of approximately €1 was chosen.

2.2.4. Erosion Protection

The economic value of erosion protection by P. oceanica meadows is obtained by calculating the percentage of the wave attenuation (WAP) given by:
WAP = H c H 0 × 100
in which H0 is the incident wave height in the offshore area, and Hc is the resulting wave height after the passage on the P. oceanica meadow. H0 values were extracted in proximity (offshore area) of each sampling site from the Copernicus product “Mediterranean Sea Wave Reanalysis” (https://resources.marine.copernicus.eu/product-detail/MEDSEA_MULTIYEAR_WAV_006_012, accessed on 9 December 2021) [25], which includes the main wave parameters (significant wave height, peak period, and mean direction) as calculated by the Med-Waves modelling system [25] in the entire Mediterranean basin, with the spatial and temporal resolutions of 0.042° (about 4.5 km) and 1 h, respectively. The data covered the period between 1993 and 2020.
For the constant depth case, Hc could be approximated using the following formula [26,27]:
H c s = H 0 exp ( Y X )
in which X is the cross-shore extension of the P. oceanica meadow, and Y is the attenuation coefficient proposed by Infantes et al. [28]:
Y = C D a v N 3 π s i n h 3 ( k p l v ) + 3 s i n h ( k p l v ) ( s i n h 2 k p h + 2 k p h ) s i n h ( k p h ) k p H 0
in which kp is peak wave number (m−1), lv is the vegetation length (m), N is the number of shoots per unit area (m2), a1v is the plant surface area per height unit (m), h is the water depth in which the wave attenuation is calculated, and CD is the drag coefficient that describes the drag associated with the individual blades. The friction due to seagrasses can be determined in terms of the Keulegan–Carpenter (KC) number [27,29] as:
C D = 22.9 K C 1.09
for KC values between 15 and 425 [27].
The wave attenuation percentage was considered for each sampling site, considering the characteristics of the P. oceanica meadows in terms of shoots density, blade lengths and leaf areas, and the offshore wave climate given by significant wave height and peak period as computed over the last 27 years.
To obtain an economic estimate, the cost of obtaining the same result through the construction of a submerged barrier, which is commonly used in coastal defense works, was calculated for the estimate, and installations of submerged barrier modules of the same size were evaluated in a non-continuous way to allow a constant exchange of water [2]. The value can be estimated using the costs of module construction, and normalizing them to obtain an energy dissipation comparable to that determined by the meadow in the measurement site. The cost per m of a structure that allows an average energy dissipation of 20% is approximately €1432.08. This cost is extended to one ha, assuming a square area of 100 m2, and the cost is spread over a period for which the structure is guaranteed (approximately 20 years). The benefit is calculated with the following relationship:
E = [ 1432.08 × ( H r m s H r m s 0 × 100 ) 20 ] × 5
in which E is the value of the benefit in € ha−1 yr−1, €1432.08 is the cost per linear m of a collapsed barrier capable of guaranteeing wave attenuation of about 20%, (Hrms/Hrms0) × 100 is the percentage of wave attenuation obtained from the model used from the Copernicus website (http://marine.copernicus.eu/services-portfolio/access-to-products/, accessed on 9 December 2021) [30], 20 is the proportion for the calculation of the value relating to the submerged barrier (modular structures in reinforced concrete, that were laid and juxtaposed on the seabed along a continuous line, parallel to the coast, and at a distance of at least one hundred m from it. This setup dissipated the energy of the wave motion, and favored the sliding of the sand towards the shore and its return to limit the erosion of the coasts). A conversion factor of 5 is given by the ratio between an area of 1 ha and the average lifetime of a submerged barrier (20 years).

2.2.5. Food Production

Regarding the calculation of the benefit connected to food production, the method presented in Marcelli et al. [2] was extended nationwide. The historical data of P. oceanica were supplemented with data regarding the catches of commercial fish species updated to 2012. These were found from the publications of the Economic Research Institute for Fisheries and Aquaculture (IREPA) (www.irepa.org, accessed on 9 November 2021) [31]. Starting in 2004, IREPA has reported the production of the 49 most commercially relevant fish species, with an indication of the monthly data relating to fishing days, catches, and revenues by regions and fishing systems. The fishing fleets considered in the IREPA data are commercial trawling, purse seine, hydraulic dredging, small-scale fishing, and passive multipurpose fishing.
The regional scale is the most suitable for the analysis of data relating to the benefit of food production, since the estimate of this parameter associated with P. oceanica meadows refers to fishing. The benefit of food production is calculated with the following relationship:
P = R r E t E d E m + R v E t
in which P is the value of the benefit in € ha−1 yr−1, R r is the revenue in euros of the species identified as a resident, E t is the total extension in hectares of P. oceanica meadows, E d is the extension of degraded meadows, E m is the extension in ha of the dead matte substrate, and R v is the revenue in euros of the visiting species. The identification of the two types of species considered in the calculation of the food production benefit was performed by extending what was described by Marcelli et al. [2] for the eight Italian regions under study.

3. Results

The dataset was used to create distribution maps of the P. oceanica meadows using the QGIS environment for each Italian region under study using the Maritime Districts as limits. Figure 2 reports the P. oceanica meadow distribution, and details of the substrates present in the Sardinia region, as an example. The creation of synoptic maps allowed us to analyze the characteristics of the P. oceanica, and the distribution of the data available for each region in detail. Subsequently, the ESs values were calculated by applying the methods described in the previous chapter. The results of the assessment of the regional and national total masses of CO2 and O2, and average volumes of N and P are reported in Table 3 and Table 4. The economic values of P. oceanica along the Italian territory are reported in Table 5, and Figure 3 and Figure 4.
This work allowed us to calculate the masses per ha of CO2 and O2, and the volumes of the average N and P (Table 3) for each Italian region under study. The maximum CO2 value was found in the Campania region (2.5 tonne ha−1 yr−1), whereas the minimum value was in Latium (0.9 tonne ha−1 yr−1). The maximum O2 value was found in Latium (76.2 tonne ha−1 yr−1), whereas the minimum value was in Apulia (21.5 tonne ha−1 yr−1). The maximum average value of N and P (L ha−1 yr−1) was found in Campania (531,590.7), whereas the minimum was in Sardinia (345,417.9). Using these calculated values, it was possible to compute the total masses for CO2 and O2 (tonne), and the total volumes of N and P (m3) at national level for a total amount of 523,354.8 tonne yr−1 of CO2; 11,179,022 tonne yr−1 of O2; and 147,553,144.2 m3 yr−1 of N and P (Table 4).
Regarding the carbon sequestration, a national average value of €38.9 ha−1 yr−1 was calculated. The maximum value was found in the Campania region (€61.1 ha−1 yr−1), whereas the minimum value was found in the Latium region (€22.60 ha−1 yr−1). In general, slight variations could be observed between the considered regions (Table 5). The average value of O2 production for the Italian P. oceanica meadows was €1929.4 ha−1 yr−1; the maximum value of O2 production was calculated for the Latium region (€3809.9 ha−1 yr−1), whereas the Apulia region showed the minimum value (€1075 ha−1 yr−1), as shown in Table 5. The Italian average bioremediation value amounted to €6396.2 ha−1 yr−1, with the maximum value for the Campania region (€7442.3 ha−1 yr−1), and the minimum value calculated for Sardinia (€4835.9 ha−1 yr−1). Some variability between regions was found, especially between the coastal areas of the north of Campania and Sicily compared to Apulia, Calabria, and Sardinia (Table 5). The average erosion protection value was €8964.7 ha−1 yr−1; the Sardinia region showed the maximum value (€15,914.92 ha−1 yr−1), whereas the minimum value of erosion protection was found in Liguria (€3629.61 ha−1 yr−1), as shown in Table 5. In this case, a marked variability could be observed. Considering that the estimate was based on the exposure of the coast to wave directions, it was possible to identify a geographical effect attributable to both the islands (Sicily and Sardinia) and coastal development (Tuscany). As an example, Figure 5 shows the comparison between the wave attenuation values of the Sardinia and Tuscany regions.
For P. oceanica food production, the average Italian value amounts to €4331.3 ha−1 yr−1, with the maximum value found in Tuscany (€8130.65 ha−1 yr−1), and the minimum value found in Sardinia (€338.97 ha−1 yr−1). The observed variability could be attributed to the different relationship between the total regional extensions of the P. oceanica, and the actual number of fishing fleets. In the case of Sardinia, the very low value reflected the low number of existing fleets when compared with the enormous extension of the existing meadows (174,500 ha), which alone represented 48.6% of the national P. oceanica extension considered in this study, whereas the fishing fleet was 13.6% of the total (equal to 9495 boats).
The percentage breakdown of the economic value of the benefits provided by the P. oceanica meadows in the considered Italian regions is shown in Table 6. In Liguria, Latium, and Apulia, the benefit that produces the greatest economic value was bioremediation (40.44%, 35.83%, and 34.18%, respectively), whereas the minimum value was related to the carbon sequestration (0.22%, 0.11%, and 0.18%, respectively). Tuscany, Campania, Calabria, Sicily, and Sardinia showed that the greatest economic value was connected to erosion protection (42.77%, 31.54%, 40.21%, 49.80%, and 70.68%, respectively), whereas the minimum value was represented by carbon sequestration (0.15%, 0.25%, 0.21%, 0.17%, and 0.15%, respectively).

4. Discussion and Conclusions

In this work, we assessed the economic value of the P. oceanica meadows’ ESs on eight Italian regions by applying the transfer of value methodology, focusing on five essential ecosystem services (Table 2). Compared to other marine phanerogams present in the Mediterranean Sea (e.g., Cymodocea nodosa, Zostera marina, and Zostera noltii), P. oceanica offers the greatest contribution to the considered ESs, as has been reported for the carbon sequestration in the result of the Deliverable C1 of the LIFE BLUE NATURA (LIFE14CCM/ES/00957; http://life-bluenatura.eu/, accessed on 6 December 2021) [32].
Even though the values of the different P. oceanica ESs obtained in this study (Table 5) could vary by orders of magnitude, the Italian national average P. oceanica ESs value obtained in this study (€21,660.5 ha−1 yr−1) is in line with data from other literature studies regarding the ESs of P. oceanica [2,11]. The results reported by Costanza et al. [11], Marcelli et al. [2], and Campagne et al. [9] are comparable in that they all considered the methodological approach (market values basic benefit transfer) and the value of the estimate (Table 7). However, in the case of the study by Campagne et al. [9], the value per annual ha was significantly lower. This discrepancy is due to the different economic evaluations carried out by the authors, especially during the transfer of the value of the individual benefits. The value reported by Vassallo et al. [10] is much higher than the others (Table 6). In this case, the reason can be traced back to the selected methodology (Emergy to Money Ratio), which is a completely different parameter from those in the other studies.
Our results highlight a different distribution of economic values of benefits at the regional level (Table 5, Figure 4). The different distribution of the erosion protection value of Tuscany, Sicily, and Sardinia compared to other regions is evident (examples of Sardinia and Tuscany are reported in Figure 5). In these regions, coastal development implies a more intense wave attenuation effect by the P. oceanica meadows than in regions with less coastal development and less exposure (as in the case of Liguria and Apulia, shown in Table 8) [33].
A similar consideration is possible concerning the values of food production in which the distribution of the estimates can be attributed to the relationship between the number of fishing fleets available at the regional level, and the actual extension of the P. oceanica meadows. This consideration explains the very low value of food production in Sardinia, even though the Sardinian meadows make up almost 50% of the total area considered in this work (Table 9).
The values of O2 production have a higher percentage weight for the regions of Liguria and Latium, whereas they constitute a lower percentage of the total value for the other regions. Bioremediation, on the other hand, appears much more constant at the regional level, with peaks in Liguria, Latium, and Campania.
Regarding carbon sequestration, the results of this study were obtained using the 2019 average annual exchange value of 1 tonne of CO2 [18]. Consequently, the assessment for subsequent years must consider the high fluctuation of this value.
Carbon sequestration is the benefit that has by far the least weight in the total calculation in all regions. The very low values found that agree with what is present in the literature [2] can be, at least partially, justified by an underestimation of the Emission Trading Scheme in the sense that the trading market is probably still incapable of registering these levels.
The applied methodology demonstrates that the ESs economic evaluation can provide an essential tool for the management of the coastal marine environment after considering the modularity of the used approach. The considered benefits can be further expanded by considering, for example, the use of tourism, which is an important economic sector in Italy, representing about 15% of the national general government spending (GDP). Tourism, together with recreational activities, are very important in Italy. However, these services are very heterogeneous on the Italian national territory. Consequently, to consider the value of these services in the total economic evaluation, in-depth sector studies are required.
The total value of the benefits obtained considering the entire extension of P. oceanica meadows in Italy is approximately eight billion euros, a significant value (about 0.44%) concerning both the national GDP (1787.7 billion euros, http://data.istat.it/, accessed on 6 December 2021) and the individual budgets of the regions under study. This aspect underlines not only the importance in terms of the value of P. oceanica meadows in Mediterranean Sea, but also the possibility of the use of this methodology of ESs economic estimation as a comparison tool for the management of the coastal areas, for the resolution of existing conflicts of use, and for the analysis of anthropogenic impacts in terms of natural capital loss of value.
The results obtained in this work can be a starting point for future works concerning the value of ecosystem services offered by P. oceanica meadows on a national scale in Italy, focused on the analysis of potential threats to these habitats, which could cause a loss of ecosystem value. Based on the analysis of potential threats, it will be possible to define mitigation and compensation works for lost ecosystem services.

Author Contributions

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

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors thank the Environmental Office of the Port Authority System of the Central Northern Tyrrhenian Sea for the support given to the development of the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ISPRA. Tematiche in Primo Piano-Annuario dei Dati Ambientali; ISPRA: Roma, Italy, 2012; ISBN 978-88-448-0602-6. [Google Scholar]
  2. Marcelli, M.; Scanu, S.; Frattarelli, F.M.; Mancini, E.; Carli, F.M. A benthic zonation system as a fundamental tool for natural capital assessment in a marine environment: A case study in the Northern Tyrrhenian Sea, Italy. Sustainability 2018, 10, 3786. [Google Scholar] [CrossRef] [Green Version]
  3. Millennium Ecosystem Assessment (MEA). Ecosystems and Human Well-Being: Synthesis; Island Press: Washington, DC, USA, 2005; ISBN 1-59726-040-1. Available online: http://www.millenniumassessment.org/en/Reports.html (accessed on 15 October 2018).
  4. Daily, G.C. Nature’s Services: Societal Dependence on Natural Ecosystems; Yale University Press: New Haven, CT, USA, 1997; pp. 454–464. [Google Scholar]
  5. Fraschetti, S.; Guarnieri, G.; Bevilacqua, S.; Terlizzi, A.; Claudet, J.; Russo, G.F.; Boero, F. Conservation of Mediterranean habitats and biodiversity countdowns: What information do we really need? Aquat. Conserv. Mar. Freshw. Ecosyst. 2011, 21, 299–306. [Google Scholar] [CrossRef]
  6. Ministero dell’Ambiente e della Tutela del Territorio e del Mare (MATTM). Praterie a Fanerogame Marine. Piante con Fiori nel Mediterraneo–Quaderni Habitat; MATTM: Rome, Italy, 2008. [Google Scholar]
  7. Telesca, L.; Belluscio, A.; Criscoli, A.; Ardizzone, G.; Apostolaki, E.T.; Fraschetti, S.; Gristina, M.; Knittweis, L.; Martin, C.S.; Pergent, G.; et al. Seagrass meadows (Posidonia oceanica) distribution and trajectories of change. Sci. Rep. 2015, 5, 12505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Weatherdon, L.V.; Friedrich, L.A.; Martin, C.S.; King, S. Experimental Seagrass Ecosystem Accounts: A Pilot Study for One Component of Marine Ecosystem Accounts. In Report Prepared for European Commission; European Commision: Brussels, Belgium, 2017. [Google Scholar]
  9. Campagne, C.S.; Salles, J.M.; Boissery, P.; Deter, J. The seagrass Posidonia oceanica: Ecosystem services identification and economic evaluation of goods and benefits. Mar. Pollut. Bull. 2015, 97, 391–400. [Google Scholar] [CrossRef] [PubMed]
  10. Vassallo, P.; Paoli, C.; Rovere, A.; Montefalcone, M.; Morri, C.; Bianchi, C.N. The value of the seagrass Posidonia oceanica: A natural capital assessment. Mar. Pollut. Bull. 2013, 75, 157–167. [Google Scholar] [CrossRef] [PubMed]
  11. Costanza, R.; d’Arge, R.; De Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neil, 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]
  12. European Commission. Directive 2008/56/EC of the European Parliament and of the Council of 17 June 2008 establishing a framework for community action in the field of marine environmental policy (Marine Strategy Framework Directive). Off. J. Eur. Communities 2008, 164, 19–40. [Google Scholar]
  13. Natural Capital Project. Available online: https://naturalcapitalproject.stanford.edu/ (accessed on 9 November 2021).
  14. Serrano, O.; Mateo, M.A.; Renom, P.; Julia, R. Characterization of soils beneath a Posidonia oceanica meadow. Geoderma 2012, 185, 26–36. [Google Scholar] [CrossRef]
  15. Romero, J.; Pergent, G.; Pergent-Ma16rtini, C.; Mateo, M.A.; Regnier, C. The detritic compartment in a Posidonia oceanica meadow: Litter features, decomposition rates, and mineral stocks. Mar. Ecol. 1992, 13, 69–83. [Google Scholar] [CrossRef]
  16. Giraud, G. Essai de Classement des Herbiers de Posidonia Oceanica (Linné) Delile. Bot. Mar. 1997, 20, 487–492. [Google Scholar] [CrossRef]
  17. Duarte, C.M.; Sand-Jensen, K. Seagrass colonization: Patch formation and patch growth in Cymodocea nodosa. Mar. Ecol. Prog. Ser. 1990, 65, 193–200. [Google Scholar] [CrossRef]
  18. Emission Trading Scheme. 2019. Available online: https://www.eea.europa.eu/data-and-maps/dashboards/emissions-trading-viewer-1 (accessed on 10 November 2021).
  19. Duarte, M.C.; Marbà, N.; Gacia, E.; Fourqurean, W.J.; Beggins, J.; Barrón, C.; Apostolaki, T.E. Seagrass community metabolism: Assessing the carbon sink capacity of seagrass meadows. Glob. Biogeochem. Cycles 2010, 24. [Google Scholar] [CrossRef] [Green Version]
  20. Gazeau, F.; Duarte, M.C.; Gattuso, P.J.; Barrón, C.; Navarro, N.; Ruiz, S.; Prairie, T.Y.; Calleja, M.; Delille, B.; Frankignoulle, M.; et al. Whole-system metabolism and CO2 fluxes in a Mediterranean Bay dominated by seagrass beds (Palma Bay, NW Mediterranean). Biogeosciences 2005, 2, 43–60. [Google Scholar] [CrossRef] [Green Version]
  21. Li, J.; Ren, Z.; Zhou, Z. Ecosystem services and their values: A case study in the Qinba mountains of China. Ecol. Res. 2006, 21, 597–604. [Google Scholar] [CrossRef]
  22. USEPA (US Environmental Protection Agency). A Compilation of Cost Data Associated with the Impacts and Control of Nutrient Pollution; USEPA: Washington, DC, USA, 2015.
  23. Gobert, S.; Lejeune, P.; Lepoint, G.; Bouquegneau, J.M. C, N, P concentrations and requirements of flowering Posidonia oceanica. Hydrobiologia 2005, 533, 253. [Google Scholar] [CrossRef]
  24. Sicuro, G. Analisi Economico-Finanziaria di Piccoli Impianti di Depurazione Anche Soggetti ad Upgrading con M.B.R. Master’s Thesis, Ingegneria per L’Ambiente ed il Territorio, Dipartimento di Ingegneria Idraulica, Geotecnica ed Ambientale, Università Degli Studi di Napoli Federico II, Naples, Italy, 2012. [Google Scholar]
  25. Korres, G.; Ravdas, M.; Zacharioudaki, A. Mediterranean Sea Waves Analysis and Forecast (CMEMS MED-Waves); AGRIS: Rome, Italy, 2019. [Google Scholar]
  26. Kobayashi, N.; Raichle, A.W.; Asano, T. Wave attenuation by vegeta-tion. J. Waterw. Port Coast Ocean Eng. 1993, 119, 30–48. [Google Scholar] [CrossRef]
  27. Sánchez-González, J.F.; Sánchez-Rojas, V.; Memos, C.D. Wave attenuation due to Posidonia oceanica meadows. J. Hydraul. Res. 2011, 49, 503–514. [Google Scholar] [CrossRef]
  28. Infantes, E.; Orfila, A.; Simarro, G.; Terrados, J.; Luhar, M.; Nepf, H. Effect of a seagrass (Posidonia oceanica) meadow on wave propagation. Mar. Ecol. Prog. Ser. 2012, 456, 63–72. [Google Scholar] [CrossRef] [Green Version]
  29. Mendez, F.J.; Losada, I.J. An empirical model to estimate the propagation of random breaking and nonbreaking waves over vegetation fields. Coast Eng. 2004, 51, 103–118. [Google Scholar] [CrossRef]
  30. Copernicus Marine Service. Available online: http://marine.copernicus.eu/services-portfolio/access-to-products/ (accessed on 9 November 2021).
  31. Istituto di Ricerche Economiche per la Pesca e l’Acquacoltura (IREPA). Available online: www.irepa.org (accessed on 9 November 2021).
  32. Life Blue Natura. Available online: https://life-bluenatura.eu/en/the-project/ (accessed on 6 December 2021).
  33. Ferretti, O.; Barsanti, M.; Delbono, I.; Furia, S. Elementi di Gestione Costiera-Parte I Tipi Morfo-Sedimentologici dei Litorali Italiani; ENEA (Ente per le Nuove tecnologie l” Energia el” Ambiente): Rome, Italy, 2003. [Google Scholar]
Figure 1. Location of the measurement sites of Posidonia oceanica characteristics along the Italian coast (decimal degrees, World Geodetic System 84).
Figure 1. Location of the measurement sites of Posidonia oceanica characteristics along the Italian coast (decimal degrees, World Geodetic System 84).
Sustainability 14 00489 g001
Figure 2. Distribution of P. oceanica and sites of community interest (SCIs) in the Sardinia region, and details of the substrates present in some areas (AC) (decimal degrees, WGS 84).
Figure 2. Distribution of P. oceanica and sites of community interest (SCIs) in the Sardinia region, and details of the substrates present in some areas (AC) (decimal degrees, WGS 84).
Sustainability 14 00489 g002
Figure 3. Percentage breakdown of economic services (ESs) values at a regional level (decimal degrees, WGS 84).
Figure 3. Percentage breakdown of economic services (ESs) values at a regional level (decimal degrees, WGS 84).
Sustainability 14 00489 g003
Figure 4. ESs economic values for each considered Italian region.
Figure 4. ESs economic values for each considered Italian region.
Sustainability 14 00489 g004
Figure 5. Regional comparison of wave motion attenuation values due to P. oceanica meadows for the Sardinia region (a) and Tuscany region (b) level (decimal degrees, WGS 84).
Figure 5. Regional comparison of wave motion attenuation values due to P. oceanica meadows for the Sardinia region (a) and Tuscany region (b) level (decimal degrees, WGS 84).
Sustainability 14 00489 g005aSustainability 14 00489 g005b
Table 1. Extension of Posidonia oceanica meadows (ha) in the Italian regions under consideration.
Table 1. Extension of Posidonia oceanica meadows (ha) in the Italian regions under consideration.
Italian RegionExtension (ha)
Liguria5050.1
Tuscany31,587.3
Sardinia174,489.7
Sicily74,928.4
Latium21,563.6
Campania9684.7
Calabria7749.9
Apulia33,842.8
Table 2. Summary of the considered benefits, and the respective data used for the calculations.
Table 2. Summary of the considered benefits, and the respective data used for the calculations.
BenefitDataTransfer of Value
Carbon sequestrationBiomass, relative densityEmission trading scheme
Erosion protectionWave height, period, directionLinear meter cost of submerged barrier
BioremediationNumber of shoots m−2Operating cost of wastewater plant
Food productionAnnual revenues, substrates distributionAnnual revenues professional fishing
O2 productionBiomass, number of shoots m−2Cost of the industrial method of Oxygen production
Table 3. Regional masses of carbon dioxide and oxygen (CO2 and O2, respectively), and L of the average nitrogen and phosphorus (N and P, respectively) per ha each year derived from the benefits of P. oceanica meadows.
Table 3. Regional masses of carbon dioxide and oxygen (CO2 and O2, respectively), and L of the average nitrogen and phosphorus (N and P, respectively) per ha each year derived from the benefits of P. oceanica meadows.
RegionCO2 Sequestration
(Tonne ha−1 yr−1)
O2 Production
(Tonne ha−1 yr−1)
Bioremediation
(L ha−1 yr−1)
Liguria1.657.4529,814.3
Tuscany1.726.7492,614.3
Sardinia1.427.8345,417.9
Sicily1.629.3493,231.4
Latium0.976.2527,321.4
Campania2.545.4531,590.7
Calabria1.724.4356,615
Apulia1.221.5378,375.7
Table 4. Total regional and national masses of C and O, and m3 of the average N and P per year derived from the benefit of P. oceanica meadows.
Table 4. Total regional and national masses of C and O, and m3 of the average N and P per year derived from the benefit of P. oceanica meadows.
RegionP. oceanica
Extension (ha)
CO2
Sequestration
(Tonne yr−1)
O2 Production
(Tonne yr−1)
Bioremediation
(m3 yr−1)
Liguria5050.18080.2289,875.72,675,615.2
Tuscany31,587.353,698.4843,380.915,560,355.7
Sardinia174,489.7244,285.64,850,813.760,271,865.7
Sicily74,928.4119,885.42,195,402.136,957,039.6
Latium21,563.619,407.21,643,146.311,370,947.7
Campania9684.724,211.8439,685.45,148,296.5
Calabria7749.913,174.8189,097.62,763,730.6
Apulia33,842.840,611.4727,620.212,805,293.1
Italy358,896.5523,354.81,117,902.2147,553,144.2
Table 5. Regional ecosystem services (ESs) values.
Table 5. Regional ecosystem services (ESs) values.
RegionCO2
Sequestration
(€ ha−1 yr−1)
O2 Production
(€ ha−1 yr−1)
Bioremediation
(€ ha−1 yr−1)
Erosion
Protection
(€ ha−1 yr−1)
Food
Production
(€ ha−1 yr−1)
Liguria40.128697417.43629.614386.4
Tuscany41.881334.96896.612,259.88130.65
Sardinia34.51392.24835.915,914.92338.97
Sicily40.61462.76905.211,620.333307.08
Latium22.63809.97382.57349.552042.1
Campania61.122717442.37778.987107.43
Calabria41.771220.54992.67904.65499.21
Apulia28.610755297.352603838.94
Italy38.91929.46396.28964.74331.3
Table 6. Breakdown of the economic value of the benefits provided by the P. oceanica meadows in the Italian regions under study.
Table 6. Breakdown of the economic value of the benefits provided by the P. oceanica meadows in the Italian regions under study.
BenefitLiguriaTuscanyLatiumCampaniaCalabriaApuliaSicilySardinia
Carbon sequestration0.22%0.15%0.11%0.25%0.21%0.18%0.17%0.15%
O2 production15.64%4.66%18.49%9.21%6.21%6.94%6.27%6.18%
Bioremediation40.44%24.06%35.83%30.18%25.40%34.18%29.59%21.48%
Erosion protection19.79%42.77%35.67%31.54%40.21%33.94%49.80%70.68%
Food production23.91%28.37%9.91%28.82%27.97%24.77%14.17%1.51%
Table 7. Ecosystem services economic assessment comparison between different studies and methodologies.
Table 7. Ecosystem services economic assessment comparison between different studies and methodologies.
Literature ReferencesMethodologyValue (€ ha−1 yr−1)
Costanza et al. [11]Market values basic benefit transfer20,300
Campagne et al. [9]Market values basic benefit transfer283–513
Vassallo et al. [10]Environmental Emergy to Money Ratio1,720,000
Marcelli et al. [2]Market values basic benefit transfer12,614
Table 8. Regional coastal length.
Table 8. Regional coastal length.
RegionCoastal Length (km)
Liguria345.176
Tuscany581.8
Sardinia1708.9
Sicily1429.5
Latium339.9
Campania447.9
Calabria712.1
Apulia859.7
Table 9. Regional fishing boats number and ratio with P. oceanica meadow extension.
Table 9. Regional fishing boats number and ratio with P. oceanica meadow extension.
RegionNumber of Fishing BoatsTrawling Percentage (%)P. oceanica Extension (ha)Ratio
Liguria525575050.10.104
Tuscany60061.331,587.30.019
Sardinia129263.4174,489.70.007
Sicily294666.874,928.40.039
Latium58275.321,563.60.027
Campania112436.29684.70.116
Calabria854617749.90.110
Apulia15726733,842.80.046
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Scanu, S.; Piazzolla, D.; Bonamano, S.; Penna, M.; Piermattei, V.; Madonia, A.; Manfredi Frattarelli, F.; Mellini, S.; Dolce, T.; Valentini, R.; et al. Economic Evaluation of Posidonia oceanica Ecosystem Services along the Italian Coast. Sustainability 2022, 14, 489. https://doi.org/10.3390/su14010489

AMA Style

Scanu S, Piazzolla D, Bonamano S, Penna M, Piermattei V, Madonia A, Manfredi Frattarelli F, Mellini S, Dolce T, Valentini R, et al. Economic Evaluation of Posidonia oceanica Ecosystem Services along the Italian Coast. Sustainability. 2022; 14(1):489. https://doi.org/10.3390/su14010489

Chicago/Turabian Style

Scanu, Sergio, Daniele Piazzolla, Simone Bonamano, Marina Penna, Viviana Piermattei, Alice Madonia, Francesco Manfredi Frattarelli, Simone Mellini, Tania Dolce, Riccardo Valentini, and et al. 2022. "Economic Evaluation of Posidonia oceanica Ecosystem Services along the Italian Coast" Sustainability 14, no. 1: 489. https://doi.org/10.3390/su14010489

APA Style

Scanu, S., Piazzolla, D., Bonamano, S., Penna, M., Piermattei, V., Madonia, A., Manfredi Frattarelli, F., Mellini, S., Dolce, T., Valentini, R., Coppini, G., Fersini, G., & Marcelli, M. (2022). Economic Evaluation of Posidonia oceanica Ecosystem Services along the Italian Coast. Sustainability, 14(1), 489. https://doi.org/10.3390/su14010489

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