Next Article in Journal
Application of Risk Analysis in the Screening of Flood Disaster Hot Spots and Adaptation Strategies
Next Article in Special Issue
Ecosystem Services Valuation for the Sustainable Land Use Management by Nature-Based Solution (NbS) in the Common Agricultural Policy Actions: A Case Study on the Foglia River Basin (Marche Region, Italy)
Previous Article in Journal
Value and Quality in the Dialectics between Human and Urban Capital of the City Networks on the Land District Scale
Previous Article in Special Issue
Green Infrastructure, Urbanization, and Ecosystem Services: The Main Challenges for Russia’s Largest Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Resolution Land Cover Integrating Copernicus Products: A 2012–2020 Map of Italy

1
Italian Institute for Environmental Protection and Research (ISPRA), Via Vitaliano Brancati 48, 00144 Rome, Italy
2
Department of Innovation in Biology, Agri-Food and Forest Systems (DIBAF), University of Tuscia, Via San Camillo de Lellis SNC, 01100 Viterbo, Italy
*
Author to whom correspondence should be addressed.
Land 2022, 11(1), 35; https://doi.org/10.3390/land11010035
Submission received: 28 November 2021 / Revised: 16 December 2021 / Accepted: 20 December 2021 / Published: 26 December 2021
(This article belongs to the Special Issue Dynamics of Urbanization and Ecosystem Services Provision)

Abstract

:
The study involved an in-depth analysis of the main land cover and land use data available nationwide for the Italian territory, in order to produce a reliable cartography for the evaluation of ecosystem services. In detail, data from the land monitoring service of the Copernicus Programme were taken into consideration, while at national level the National Land Consumption Map and some regional land cover and land use maps were analysed. The classification systems were standardized with respect to the European specifications of the EAGLE Group and the data were integrated to produce a land cover map in raster format with a spatial resolution of 10 m. The map was validated and compared with the CORINE Land Cover, showing a significant geometric and thematic improvement, useful for a more detailed and reliable evaluation of ecosystem services. In detail, the map was used to estimate the variation in carbon storage capacity in Italy for the period 2012–2020, linked to the increase in land consumption

1. Introduction

Territory is a source of resources such as food, biomass and raw materials and provides essential ecosystem services that support production functions, regulate natural cycles, provide cultural and spiritual benefits [1].
The European Environment Agency introduced the concept of “land system”, which defines the territory as a set of terrestrial components that includes all the processes and activities related to its anthropic use [2,3]. This concept considers the territory as an integrated system [4] that combines everything related to land use and land cover. The study of land cover and land use is essential to understand the causes and effects of human-determined changes [5], seeing as now changes and transformations that occur within the land system lead to consequences for the well-being of humans and the environment at the local, regional and global levels.

1.1. Estimation and Monitoring of Ecosystem Services

The management of territory is fundamental, since its transformation alters the environmental processes and related ecosystem services [6], which are the benefits that humans obtain directly or indirectly from terrestrial ecosystems [7].
Since 2005, the concept of ecosystem services has been placed on the political agenda [8,9,10] thanks also to the Millennium Ecosystem Assessment (MEA). The MEA is an international research project launched in 2001 with the aim of evaluating the consequences of ecosystem changes on human well-being and establishing actions to improve the conservation and sustainable use of ecosystems and their contribution to health.
The Millennium Ecosystem Assessment [11] provides a classification of ecosystem services which consist of four groups:
Provisioning services, which provide products obtained from ecosystems, such as food, raw materials and water.
Regulating services, i.e., the benefits provided through ecosystem processes, such as carbon storage, erosion control.
Cultural services, which represent the nonmaterial benefits that people obtain through spiritual enrichment, cognitive development and aesthetic experience.
Supporting services, which constitute a “transversal” category that supports the production of other services, providing living space for plants and animals or maintaining genetic diversity. They differ from other categories since their impact on people is indirect or is visible after a very long period.
The Economics of Ecosystems and Biodiversity [12] provides a classification aligned with MEA, while the Common International Classification of Ecosystem Services (CICES) excludes “support services” and renames the category of “regulation services” as “regulation and maintenance services”, including habitat maintenance [13,14].
In 1997, Costanza [7] published one of the first assessments of ecosystem services, which was then updated in 2014 [15]. Other application examples for more limited areas are readily available in the scientific literature: assessment on pollination can be found both from a biophysical [16,17] and economic [18,19,20] point of view. The United Kingdom was among the first European country to draw up a complete official report in line with the Millennium Ecosystem Assessment [21]. Rabe et al. [22] developed a network of national indicators on ecosystem services in Germany using CORINE Land Cover data for the analysis of nine ecosystem services divided into three classification categories (according to CICES). Spain published an official report in 2016 [23] related to twelve ecosystem services. Outside the European context, there are applications on a national scale in China [24], South Africa [25] and South Korea [26].
Soil has fundamental functions for nature and humankind and is the source of many ecosystem services [27,28,29]:
Fertility: the nutrient cycle ensures fertility in the soil and, at the same time, the release of nutrients necessary for plant growth;
Filter and reserve: the soil can act as a filter against pollutants and can store large quantities of water, useful for plants and for the mitigation of floods;
Structural: soils represent the support for plants, animals and infrastructures;
Climate regulation: the soil, in addition to being the largest carbon sink, regulates the emission of important greenhouse gases (N2O and CH4);
Conservation of biodiversity: soils are an immense reservoir of biodiversity. They represent the habitat for thousands of species capable of preventing the action of parasites or facilitating waste disposal;
Resource: soils can be an important source of supply of raw materials.
All soils perform their functions at the same time (food production, water purification, carbon sequestration, etc.) in a different way according to land use and pedogenetic characteristics. For example, the rate of carbon sequestration and water purification is higher in a natural area than in an agricultural one, which however has greater production capacity.
Carbon sequestration and storage is an important regulatory service linked to the attitude of ecosystems to fix greenhouse gases. This service contributes to climate regulation and is fundamental in defining adaptation strategies to climate change [30,31]. The capacity to store carbon depends, among other things, on land use and cover and on the climate [32]. Land use, land use change and forestry (LULUCF) activities can act as sources of emissions or store carbon by acting as sinks. In particular, natural and seminatural forest ecosystems have the highest carbon sequestration potential. Once natural land is urbanized or degraded, it loses its ability to retain carbon which, as a result, is emitted into the atmosphere [33]. Urban expansion, land consumption, deforestation and forest degradation limit the ability of natural areas to store carbon and have contributed to these emissions by releasing carbon stored in forests, vegetation and soil [6,34,35].

1.2. Land Monitoring

Land cover and land use are strongly related and for many applications both information categories are required [36]. To meet different monitoring needs, data with different characteristics from a spatial, temporal and thematic point of view were introduced.
In this respect, different initiatives have been developed. The purpose of the Copernicus program is to collect information on the earth’s surface and organize it according to criteria that allow to compare different data, to exchange data between EU countries and to increase the number of users. The Copernicus Land Monitoring Service (CLMS) allows researchers to obtain geographic information on soils and on numerous variables related to them (such as the state of the vegetation or the water cycle), supporting applications in a wide variety of sectors, such as territorial planning, management of water resources and forests, agriculture and food security. CORINE Land Cover is one of the main products belonging to CLMS. It has guaranteed information for the whole European territory since 1990, with 44 land cover and land use classes and geometric detail of 25 hectares.
Recently, data with higher spatial and thematic details have been introduced in the context of the CLMS Local Component. It aims at providing detailed information on critical areas from an environmental point of view, which require specific and detailed monitoring. Currently, this Copernicus component offers land cover and land use maps in vector format, with high spatial resolution and a 6-year update frequency for four categories of areas. Urban Atlas refers to the CLC classification system, describing with higher detail the land cover and land use characteristics of urban areas, while Riparian Zone and Natura 2000 use the ecosystem types defined in the Mapping Assessment of Ecosystems and their Services (MAES) [37], which are based on the CLC classes too.
These aforementioned data adopt classification systems based on different combinations of land cover and land use classes that are difficult to compare and integrate with those of other data. In order to coordinate data flows from a thematic point of view, the EAGLE group (EIONET Action Group on Land monitoring in Europe) was created. It aims at defining a conceptual methodology to describe land cover and land use information in a consistent data model. EAGLE is not a classification system but a tool to describe classes of a given classification system by tracing them to the segments related to the three categories. This allows to better understand the characteristics, the overlaps and the possible conversions between different classification systems and provides a basis to define new ones. The EAGLE model aims at separating the land cover and land use components through data modelling systems applicable at different scales and in different contexts, while maintaining compatibility with existing databases.
The EAGLE data model is based on the definition of three blocks, called “categories”:
  • Land cover components (LCC), which refer to the definition of "land cover" provided by the INSPIRE directive 2007/2/CE. The LCCs are mutually exclusive and exhaustive and can be used as a modelling element to semantically describe a class definition or to map landscape;
  • Land use attributes (LUA), that follow in principle the Hierarchical INSPIRE Land Use Classification System (HILUCS), with some changes to fit the purpose of the EAGLE concept. The LUA are attached to the land cover unit;
  • Landscape characteristics (CH), which describe further details of the land cover components. The first level distinguishes “land management”, “spatial pattern”, “crop type", “mining product type”, “ecosystem types”, “height zone”, “(bio-)physical characteristics”, “general parameters”, “status” and “temporal” parameters. This enhances the integration between national activities and European land monitoring initiatives encouraging a bottom-up approach in data production.
The problem of interoperability and non-homogeneity between data is also evident at a national level. The National Land Consumption Map offers national coverage, with annual update and EAGLE compliant classification system, while most of the data available at the regional level are inconsistent, not updated and difficult to relate to each other.
Despite the large amount and variety of land cover and land use data available at national and European level, currently CLC is the only product capable of supporting an assessment of ecosystem services on a national scale [38], since it guarantees the mapping of the entire national territory and has a thematic detail suitable for the purpose. However, the low spatial resolution and the presence of mixed classes reduce the reliability of the assessments based on them.
In this sense, the first objective of this research concerns the development of a methodology that makes the main Copernicus and national land cover and land use data comparable and integrable, in order to obtain a product with national coverage that allows to overcome the limits of the CLC in terms of classification system and geometric detail.
Furthermore, the activity refers to an EAGLE compliant classification system with a thematic detail useful for conducting an assessment of ecosystem services, with particular reference to the variation of carbon stocks. This change was assessed with respect to theincrease in land consumption between 2012 and 2020.

2. Materials and Methods

2.1. Overview

The methodology presented in this study integrates Copernicus and national data for the production of a land cover map capable of supporting the ecosystem services assessment. Data were reclassified according to an EAGLE compliant classification system and merged into a 10 m resolution land cover map of Italy. The map was used to assess the loss of carbon storage capacity for the period 2012–2020, associated with land consumption (Figure 1).

2.2. Study Area

The analysis was carried out for the entire Italian territory (Figure 2), which covers 301,338 km2. The country is composed mainly of hills (41.0%) and mountains (35.0%), while the remaining 23.0% of the territory is covered by plains. To the north is the mountain range of the Alps, which exceeds 4,000 m in altitude. In this area, the alpine climate prevails, with high rainfall with a maximum of 2,500–3,500 mm. In the peninsular area there is the Apennine mountain range, reaching its highest peak in Abruzzo with Gran Sasso (2,912 m) and characterized by a continental climate. The coastal area has a Mediterranean climate, with average annual rainfall that reaches a minimum of 500 mm in Apulia and Molise.
Land cover is characterized by forest in mountain areas, with conifer concentration in alpine areas. Crops and most of the urbanized areas are concentrated in the plains and along the coast.

2.3. Land Cover Classification System

The activities described in this paper refers to a sixteen-class classification system. The classes are defined in accordance with the EAGLE group specifications [39,40,41] and are organized into five levels (Table 1).
The classification system is based on previous activities of the working group [41] and improved to maintain the thematic detail offered by the Copernicus and national input data. The first three classes coincide up to the third level of detail with Eagle concept land cover components. Wetland class and the fourth and fifth classification levels are based on EAGLE characteristics (LCH) definitions.
  • Abiotic non-vegetated: The class includes any unvegetated surfaces. At the second classification level, the class is subdivided between man-made artificial structures (artificial abiotic surfaces) and natural material surfaces (natural abiotic surfaces).
Artificial abiotic surfaces include impervious surfaces and reversible land consumption, according to the definition of ISPRA National Land Consumption Map [42]. Reversible land consumption indicates areas where the natural cover has been removed or replaced due to anthropogenic interventions, such as soil compaction or excavation, forming a non-impermeable and undeveloped surface. The main difference from the EAGLE model is the inclusion of quarries and extraction sites in this class, as detailed by De Fioravante et al. [41]. Natural abiotic surfaces are any kind of material that remains in its natural consistence or form, either with or without anthropogenic influence. The latter class presents a third-level distinction between consolidated unvegetated natural surfaces (bare rocks, cliffs) and unconsolidated unvegetated natural surfaces (beaches, dunes, sands).
2.
Biotic vegetated: The class includes any vegetated surfaces, with or without anthropogenic influence. At the second classification level, woody and herbaceous vegetation are distinguished.
The third classification level of woody vegetation divides trees and shrubs. Trees are then classified at the fourth level among broad-leaved, needle-leaved and permanent crops. Broad-leaved and needle-leaved have a fifth level based on CORINE Land Cover, while permanent crops distinguish olive groves, orchards and woody plantation. Shrubs are distinguished on the fourth level in vineyards and natural areas.
Herbaceous vegetation is divided into periodic and permanent. The periodic herbaceous class corresponds to the managed areas, subdivided at the fourth level into pastures and arable land, while the natural herbaceous is related to natural unmanaged grassland.
3.
Water surfaces: The class includes natural or artificial solid and liquid water. The second classification level distinguishes water bodies from permanent snow and ice. Water bodies includes liquid water regardless of shape, position, salinity and origin. Permanent snow and ice includes accumulations that persist throughout the year regardless of seasonal variations.
4.
Wetlands: The class does not have a direct correspondence with the EAGLE LCC, as it is considered an LCH. The class was however included to maintain the information content offered by the input data. In detail, a definition was adopted aligned with the CORINE Land Cover, including in the class the inland wetlands (inland marshes and peat bogs) and coastal wetlands (salt marshes, salines, intertidal flats) while lagoons and estuaries are associated to water bodies.

2.4. Selection of Input Data

The map is based on the integration of national and European data. At the European level, the main CLMS data for 2012 were considered, which is the reference year of most of the available data.
At the national level, the National Land Consumption Map (LCM) was used for the identification of artificial abiotic surfaces and some regional land cover and land use maps supported the characterization of agricultural areas.

2.4.1. CLMS Data

Data part of the local and pan-European components of the CLMS were considered. Regarding the local component, Urban Atlas (UA), Riparian Zones (RZ) and Natura 2000 (N2000) data were used (Table 2), while for the pan-European component, reference was made to CORINE Land Cover (CLC) and the high resolution layers (HRL) (Table 3).

2.4.2. National Data

LCM was used as national data [41,43,44,45]. It is a 10 m raster available for the entire Italian territory. In addition, regional LC/LU maps were selected for Apulia, Latium, Abruzzo, Veneto, Liguria, Lombardy and Basilicata, in order to increase the spatial and thematic detail in areas not covered by the Copernicus local component data.

2.5. Production of a National Land Cover Map Based on Copernius and National Data

2.5.1. Data Pre-Processing

Data were projected in the European Terrestrial Reference System 1989 (ETRS-89) and in Lambert Azimuthal Equal Area (LAEA) projection. Then they were converted into 10 m resolution raster, aligned with LCM and HRLs data (the only two data already in raster format).

2.5.2. Map Production

Map production involves the steps described below. From a geometric point of view, priority was given to Copernicus local products (UA, RZ, N2000), which offer better spatial resolution than other available land cover/land use data. In the areas not covered by the Local data, the CLC and the regional maps available for 2012 (derived from CLC) have been inserted.
From a thematic point of view, the input data were homogenized with respect to the classification system of Table 1. For the mixed classes, Copernicus HRL data were used to distinguish the woody component from herbaceous vegetation, then reference was made to the definitions of such mixed classes to distinguish natural areas from agricultural ones. The CLC data made it possible to attribute a detailed prevalent forest categories to areas classified as broad-leaved and needle-leaved, while the LCM locates the consumed land.

Reclassification of Copernicus UA, RZ, N2000 Data and Creation of the Basic Mosaic

Local component data show the highest thematic and geometric detail, therefore they were used as a basis for the map. They were first reclassified according with the classification system defined in Table 1, while temporary codes have been introduced for the classes without a direct correspondence (Table 4) and for the artificial surfaces and mixed forests. These areas have been assigned a land cover class according to the procedure described below, starting from the information provided by the HRL, LCM and CLC data.
Reclassification was carried out considering the correspondences of Table A1 for N2000, of Table A2 for UA and of Table A3 for RZ, which are reported in Appendix A.
The three reclassified data were merged. Where more data were present at the same time, priority was given to UA, as the higher spatial resolution allows for a more detailed description of the complex pattern that characterizes urban areas. Regarding RZ and N2000, the comparison between the two data in the overlapping areas showed greater detail in RZ’s description of the territory, which was considered a priority over N2000.

Reclassification and Introduction of CLC Data and Regional Maps

The CLC in 10 m raster format has been reclassified considering the correspondences of Table A4 (see Appendix A). CLC was then integrated with regional data. In detail, the polygons related to the following land cover classes (Table 1) were exported from the regional maps: 121, 122, 21131, 21132, 2121, 2122, 2212 and 222. They have been converted to raster and superimposed on the CLC in order to allow the detection of patches smaller than the CLC MMU. These data were then added to the mosaic in the areas not covered by UA, RZ and N2000.

Use of HRLs for the Classification of Uncertain Classes (Temporary Code 91–99)

The HRLs DLT, GRA and WAW related to 2012 have been mosaicized and used to assign a land cover class to the uncertain classes of Table 4.
In detail, the mixed land cover classes were compared with the HRLs DLT and GRA from the geometric point of view and with reference to the class definitions, in order to identify the data matches (Table 5). This comparison made it possible to distinguish the arboreal component (identified by HRL Forest) from the herbaceous and the unvegetated areas (identified by HRL Grassland) in the mixed classes. To distinguish natural vegetation from the agricultural areas, reference was made to the initial definition of mixed classes.
HRL WAW was used to refine water bodies. In areas covered by high resolution data (UA, RZ, N2000) the water bodies of the HRL WAW are classified as 122 “unconsolidated natural abiotic surfaces”. In areas covered only with CLC, water bodies mapped by the HRL were added, as they are smaller than the CLC minimum mapping unit (MMU).

Use of HRLs to Classify Broad-Leaved and Needle-Leaved in the Mixed Forest (Temporary Code 999)

The mixed forest pixels coming from UA, RZ and N2000 were distinguished in broad-leaved and needle-leaved through HRL DLT. In the mixed forest pixels where there was no correspondence with the HRL, the class was attributed according to a proximity criterion, starting from the Euclidean allocation made on HRL DLT.

Use of the Fourth CLC Level for the Attribution of a More Detailed Prevalent Class to Broad-Leaved and Needle-Leaved Pixel

Broad-leaved and needle-leaved woodlands were reclassified to the fifth classification level (not shown in Table 1) on the basis of the fourth CLC classification level available for the Italian territory. In the forest pixels without direct correspondence with the fourth CLC level, the more detailed class was assigned according to a proximity criterion, starting from the Euclidean allocation conducted on broad-leaved and needle-leaved CLC polygons.

Inclusion of LCM and Assignment of the Class to Temporary Codes 998

The LCM for 2012 and 2020 were superimposed on the map obtained in the previous steps, since it is the most detailed available data. Pixels classified as 998 that do not fall within the LCM are related to urban areas of the Copernicus data without artificial land cover. They have therefore been attributed to the “permanent herbaceous” class.

2.5.3. Accuracy Assessment

The 16-class map for 2012 was validated. Accuracy assessment consists of a first phase of quality control conducted through a systematic visual search for macroscopic errors. A quantitative accuracy assessment was then performed through the photointerpretation of a sample of points, which were then compared with the values of the land cover map at the same locations. The sample size was assessed using the methodology proposed by Olofsson [46], which is widely adopted in literature [47,48,49].
The sample size (n) is calculated starting from the areas of each class and from the definition of a first attempt user accuracy, using the following equation [47,48,49,50]:
( W i S i ) 2 [ S ( O ^ ) ] 2
where:
  • Wi—is area proportion of each classes in the considered map
  • Ui—user accuracy of class i. A conservative scenario was assumed, considering Ui = 0.6 for all classes.
  • Si—standard deviation of stratum i, Si = √(Ui(1 − Ui)) [50]. Considering Ui = 0.6, it turns out Si = 0.49 for all classes.
  • S(Ô)—is the target standard error for overall accuracy. It was assumed to be 0.01 as suggested by Olofsson [46], which corresponds to a confidence interval of 1%.
  • A sample of size 2400 was obtained (Table 6).
The 2400 points were distributed among the classes considering the average value between equal and area-proportional distribution (see https://fromgistors.blogspot.com/2019/09/Accuracy-Assessment-of-Land-Cover-Classification.html, accessed on 25 December 2021) (Table 7).
A stratified random sampling was conducted to identify on each class the corresponding number of points calculated in Table 7. The points were photointerpreted with very high resolution images, considering 2012 as the reference year.

2.5.4. Ecosystem Services—Carbon Storage Capacity Assessment

The land cover map was used as input data for the assessment of ecosystem services, and in particular for the carbon storage capacity. The analysis is based on a simplified scheme that considers the amount of carbon stock constant over time. The variation of this ecosystem service refers to two versions of the land cover map, which exploit the National Land Consumption Map for 2012 and 2020. In this sense, the variation in carbon storage capacity must be understood as a reduction linked to the increase in land consumption. Actually, while urbanization improves human social and economic well-being, on the other hand it has a negative impact on human ecological well-being, which is closely related to the level of ecosystem services and which puts urban development and human well-being at risk in the future [6].
As input data for the carbon stock, many different bibliographic sources were used and integrated, such as the National Inventory of Forests and Forest Carbon Tanks (INFC) and the recent map of organic carbon created as part of the activities of the Global Soil Partnership [51]. Specific coefficients were used to identify the contribution deriving from the different pools [52,53,54]. There are four main pools of carbon in nature [55], recognized and classified by the Intergovernmental Panel on Climate Change [56], which are analyzed for each portion of the territory and each type of land cover:

Above-Ground Biomass (AGB)

It includes all the tissues of plant organisms outside the soil (such as stems, branches, leaves, seeds, etc.). The volume is calculated as:
AGB = a * GSV + b * GSV * ec * GSV
where:
GSV = growing stock volume
a, b, c = specific coefficients for each forest type [53]
To switch from biomass to the fraction of stored carbon, the values are multiplied by 0.5

Below-Ground Biomass (BGB)

It includes the root system of plants. The volume is calculated as [54]:
BGB = GSV * BEF * WBD * R
where:
  • GSV = growing stock volume
  • BEF = biomass expansion factor
  • WBD = wood basic density
  • R = crown/roots ratio, tabulated for the different species [52,54].
To switch from biomass to the fraction of stored carbon, the values are multiplied by 0.5

The Carbon Contained in the Dead Organic Substance (DOS)

The pool includes the necromass, the woody plant residues, the litter, the finer residues not yet decomposed. As regards the epigeal biomass, specific multiplicative coefficients are considered to be applied to the values obtained from the calculation shown above, for example 0.20 for evergreen plants and 0.14 for deciduous trees [56].
Specific formulas for each species present in the bibliography were used for the litter [52,54].

Soil Carbon

The pool includes organic and mineral layers including up to a depth of 30 cm. It is evaluated starting from the data produced by CREA-ABP, CNR-Ibimet, regions and some universities as part of the Global Soil Partnership/FAO initiative [51] as an Italian contribution to the Global Soil Organic Carbon map. The map offers the values of the carbon contained in the soil in raster format with a resolution of 1 km.
In detail, for the forest cover areas, data from the National Inventory of Forests and Forest Carbon Tanks (INFC) were used [57]. The inventory provides differentiated values both by region and according to the different plant species, with reference to the classes of the CLC fourth classification level.
For the other land cover classes, estimates from the literature were used: the pool values for artificial areas were considered zero while for the other natural and agricultural surfaces the literature values reported in Table 8 were used [34].
With reference to the carbon values of permanent crops, the values of Table 9 were considered [58].

3. Results

3.1. Land Cover Map and Accuracy Assessment

The map of Figure 3 was obtained by applying the procedure described in the previous chapter.
Accuracy assessment was conducted on the map of Figure 3 and provides the results shown in Figure 4 and summarized in Table 10. The map has an overall accuracy of 90%. The omission error is less than 20% in all classes and slightly higher than this value for permanent snow and ice and wetlands.
The commission error is just over 20% in two of the 16 classes, while it is less than 3% for the artificial abiotic surfaces, beaches, dunes and sands, olive groves and permanent snow and ice.
Analyzing the first classification level of the 2012 land cover map, over 88% of the national surface is vegetated, followed by abiotic surfaces (9.83%) and water bodies and wetlands (1.46 and 0.16%). At the second classification level, in the abiotic class the artificial component prevails (6.98%) followed by bare rocks. Regarding vegetation, woody and herbaceous vegetation show comparable surfaces (45.7% and 42.7% respectively). In the woody class the arboreal component prevails, which occupies 38.91% of the national surface (Table 11).
Considering the maximum thematic detail, class 2212 (arable land) prevail, which occupies 30.47% of the national territory, followed by broad-leaved trees (25.90%). All other classes occupy less than 10% of the territory and 11 out of 16 classes less than 5% (Table 12).

3.2. Estimation of the Carbon Storage Capacity

The application to the 2012 map of the methodology for calculating the carbon stocks obtained the results shown in Figure 5, which show a strong concentration of the carbon stocks in the alpine and mountain areas, while the value is significantly reduced in the agricultural areas of the Padana plain and in particular in Sicily and in the Tavoliere delle Puglie.
In Italy, 2,898,672 tons of stored carbon (stock) were lost due to land consumption between 2012 and 2020 (Table 13). In detail, this value relates to transformations from natural to artificial land cover, excluding restorations and changes between other different land cover classes. Analysing the results on a regional scale, almost a quarter of the total losses are concentrated in the Veneto (384,537 tons, equal to 13.27% of the total) and Lombardy (319,666 tons, 11.03% of the total), while each of the other regions is affected by less than 10% of the changes. The minor losses are in Aosta Valley (13,206 tons), Molise (21,242 tons) and Liguria (25,237 tons), which together host less than 2% of the changes.

4. Discussion

The effectiveness of ecosystem services assessment methodologies is linked to the availability of spatial data that allow an accurate and reliable description of the territory [59,60].
The CLC is still one of the few data able to guarantee information on land cover and land use for the entire national territory and is a reference for studies on a national scale. However, the new monitoring needs in terms of updating frequency and spatial resolution have revealed the limits of the project. The increase in geometric detail guaranteed by the introduction of high resolution data from the Copernicus local component and national data made it possible to describe the territory more accurately than the CLC in critical areas, such as urban areas, riparian areas and protected areas, although important portions of the national territory are still covered only by CLC.
Figure 6 shows a more detailed representation of green urban areas, which is relevant if we consider that the ecosystem functions and services that urban soils are able to offer are often neglected [61].
The increase in geometric detail also makes it possible to improve the representation of small urban agglomerations, roads and riparian areas (Figure 6, Figure 7, Figure 8 and Figure 9).
Anyway, all the Copernicus data used for this research present critical issues from a thematic point of view, since their classification system is based on land cover and land use attributes with numerous mixed classes, such as the CLC class 2.4 (heterogeneous agricultural areas) and 3.2 (shrub and/or herbaceous vegetation associations).
The proposed methodology allows a more detailed description of the mixed classes, (Figure 8 and Figure 9) and makes it possible to distinguish the arboreal component from the herbaceous one, the agricultural areas from the natural ones, arable land from permanent crops. That allows the carbon storage capacity of the different classes to be analysed separately and more accurately. Although aspects such as the diameter of the trunk, the density of the trees or the characteristics of the undergrowth are not considered, the map distinguishes the tree species in a detailed way, maintaining the fourth classification level of CLC, which is available only for the Italian territory [62].
An aspect that will require further development is the possibility of defining the correspondence between the classes of Copernicus data and a classification system oriented to the description of habitats, which would be more functional for conducting studies on ecosystem services. In this sense it would be necessary to integrate ancillary data not easily available on a national scale [63].
The data available at the time of the research limited the study on carbon stocks only to variations caused by increased land consumption. The availability of new Copernicus data updated to 2018 will allow the production of land cover maps capable of evaluating the variations of ecosystem services associated with other land cover changes occurred between 2012 and 2018 [60]. However the update frequency remains too low for numerous monitoring activities. Actually, the LCM is updated annually, while Copernicus data every 6 years and the maps available at regional level in Italy are often based on CLC data and are updated in a few cases (Lombardy, updated to 2018, Tuscany 2016, Liguria 2018, Latium 2016), while other maps are less up to date (2012 for Sicily and Apulia, 2013 for Abruzzo and Basilicata, 2006 for Calabria and 2009 for Campania).
Initiatives such as the next “CLC Plus” are expected to be decisive in the near future, guaranteeing the introduction of updated and interoperable products, more suitable for carrying out the monitoring activities necessary to meet institutional needs. ISPRA is conducting other research activities in this direction [41,43,44,45,64], through the definition of a land cover classification methodology for the production of maps with Sentinel resolution, annual update frequency and EAGLE compliant classification system, capable of providing updated and reliable products for monitoring on a national scale, which can be integrated with the activities of the “CLC Plus” and the National Strategic Plan for the Space Economy.

5. Conclusions

Since the 19th century, anthropogenic activities have led to a significant increase in the level of carbon dioxide in the atmosphere [35] and negatively affected the regeneration capacity and balance of ecosystems. Urban expansion, deforestation and forest degradation have contributed to these emissions by releasing the carbon stored naturally in forests, vegetation and soil [34,35]. This type of carbon is added to greenhouse gas emissions related to industries and energy production and to the products of impure combustion. Terrestrial ecosystems are able to sequester as much carbon as is currently in the atmosphere but over the course of the century terrestrial biosphere is likely to become a net source of carbon due to factors connected with climate change, pollution and the over-exploitation of resources that will alter the structure, reduce biodiversity and perturb functioning of most ecosystems, and compromise the services they currently provide.
The land cover and land use changes involve landscape fragmentation, reduction of biodiversity and loss of green areas important for carbon accumulation and more generally for the provision of ecosystem services. Current conservation practices are generally poorly prepared to adapt to this level of change, and effective adaptation responses are likely to be costly to implement.
The monitoring of carbon stock accounting is an institutional duty enshrined in the Kyoto protocol and the Paris agreements and an important driver in defining adaptation strategies to climate change [33]. Effective monitoring strategies of land cover and land use changes are essential for studying the phenomenon. In this sense, the methodology presented in this paper represents a step forward for large-scale assessments of ecosystem services more relevant to reality, since compared to the CLC it provides products for the entire national territory with greater geometric detail and a better description of mixed areas.
The methodology is easily applicable in other territorial areas, since it is based on Copernicus data available for many European countries, furthermore the use of an EAGLE compliant classification system makes the methodology easily adaptable to the specific availability of national data.
The main limitation of the methodology concerns the low update frequency of the input data, which limited the monitoring of ecosystem services to variations related to land consumption, which is the only data updated annually for the Italian territory.
A first future development concerns the updating of the map using the new Copernicus Local and Pan-European products for 2018. This implementation will allow the evaluation of the variations in the carbon stocks associated not only with land consumption but also with other land cover changes.
The products of the application of the methodology are also in continuity with other activities carried out by the working group and can constitute a useful support tool for the development of land cover classification methodologies with high update frequency for the satisfaction of the institutional needs envisaged by the new “Space Economy” Strategic Plan and for the creation of “Istances” in the new CLC Plus Project. In this sense, an important added value of this research is linked to the suitability of products with respect to present and future national and European initiatives and standards in the field of remote sensing.

Author Contributions

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

Funding

This research was funded by the Italian Institute for Environmental Protection and Research (ISPRA) structural funds.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data presented in this study are available on request from the corresponding author. The data are not publicly available because they are part of ongoing research.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Below are the tables containing the correspondences referred to for the first reclassification of the data N2000 (Table A1), UA (Table A2), RZ (Table A3) and CORINE Land Cover (Table A4). In this first reclassification there are the temporary codes of Table 3.
Table A1. Reclassification table for N2000.
Table A1. Reclassification table for N2000.
N2000Class NameLCN2000Class NameLC
1110Urban fabric (predominantly public and private units)9984212Seminatural grassland without woody plants (C.C.D. ≤ 30%)93
1120Industrial, commercial and military units9984220Alpine and subalpine natural grassland22200
1210Road networks and associated land9985110Heathland and Moorland22200
1220Railways and associated land9985120Other scrub land21220
1230Port areas and associated land9985200Sclerophyllous vegetation21220
1240Airports and associated land9986100Sparsely vegetated areas96
1310Mineral extraction, dump and construction sites9986210Beaches and dunes12200
1320Land without current use9986220River banks12200
1400Green urban, sports and leisure facilities9986310Bare rocks and rock debris12100
2110Arable land221206320Burnt areas (except burnt forest)97
2120Greenhouses221206330Glaciers and perpetual snow32000
2210Vineyards, fruit trees and berry plantations997100Inland marshes40000
2220Olive groves211327210Exploited peat bog40000
2310Annual crops associated with permanent crops927220Unexploited peat bog40000
2320Complex cultivation patterns928110Coastal salt marshes40000
2330Land principally occupied by agriculture with significant areas of natural vegetation938120Salines40000
2340Agroforestry948130Intertidal flats40000
3110Natural and seminatural broad-leaved forest21118210Coastal lagoons31000
3120Highly artificial broad-leaved plantations21118220Estuaries31000
3210Natural and seminatural coniferous forest21129110Interconnected water courses31000
3220Highly artificial coniferous plantations21129120Highly modified water courses and canals31000
3310Natural and seminatural mixed forest9999130Separated water bodies belonging to the river system31000
3320Highly artificial mixed plantations999210Natural water bodies31000
3410Transitional woodland and scrub959220Artificial standing water bodies31000
3420Lines of trees and scrub959230Intensively managed fish ponds31000
3500Damaged forest959240Standing water bodies of extractive industrial sites31000
4100Managed grassland2211010000Sea and ocean31000
4211Seminatural grassland with woody plants (C.C.D. ≥ 30%)93
Table A2. Reclassification table for UA.
Table A2. Reclassification table for UA.
UAClass NameLCUAClass NameLC
11100Continuous Urban Fabric (>80%)99813400Land without current use998
11210Discontin. Dense Urban Fabric (50–80%)99814100Green urban areas91
11220Discontin. Medium Density Urban Fabric (30–50%)99814200Sports and leisure facilities998
11230Discontinuous Low Density Urban Fabric (10–30%)99821000Arable land (annual crops)22120
11240Discontinuous Very Low Density Urban Fabric (<10%)99822000Permanent crops (vineyards, fruit trees, olive groves)99
11300Isolated Structures99823000Pastures22110
12100Industrial, commercial, public, military and private units99824000Complex and mixed cultivation patterns92
12210Fast transit roads and associated land99825000Orchards at the fringe of urban classes21131
12220Other roads and associated land99831000Forests999
12230Railways and associated land99832000Herbaceous vegetation associations (natural grassland, moors...)22200
12300Port areas99833000Open spaces with little or no vegetations (beaches, dunes, bare rocks, glaciers)96
12400Airports99840000Wetland40000
13100Mineral extraction and dump site99850000Water bodies31000
13300Construction sites998
Table A3. Reclassification table for RZ.
Table A3. Reclassification table for RZ.
RZClass NameLCRZClass NameLC
11110Continuous Urban Fabric (IMD ≥ 80%)99841000Managed grassland22110
11120Dense Urban Fabric (IMD ≥ 30–80%)99842100Seminatural grassland93
11130Low Density Fabric (IMD < 30%)99842200Alpine and subalpine natural grassland22200
11200Industrial, commercial and military units99851100Heathland and Moorland22200
12100Road networks and associated land99851200Other scrub land21220
12200Railways and associated land99852000Sclerophyllous vegetation21220
12300Port areas and associated land99861000Sparsely vegetated areas96
12400Airports and associated land99862100Beaches and dunes12200
13100Mineral extraction, dump and construction sites99862200River banks12200
13200Land without current use99863100Bare rocks and rock debris12100
14000Green urban, sports and leisure facilities99863200Burnt areas (except burnt forest)97
21100Arable land2212063300Glaciers and perpetual snow32000
21200Greenhouses2212071000Inland marshes40000
22100Vineyards, fruit trees and berry plantations9972100Exploited peat bog40000
22200Olive groves2113272200Unexploited peat bog40000
23100Annual crops associated with permanent crops9281100Coastal salt marshes40000
23200Complex cultivation patterns9281200Salines40000
23300Land principally occupied by agriculture with significant areas of natural vegetation9381300Intertidal flats40000
23400Agroforestry9482100Coastal lagoons31000
31100Natural and seminatural broad-leaved forest211182200Estuaries31000
31200Highly artificial broad-leaved plantations211191100Interconnected water courses31000
32100Natural and seminatural coniferous forest211291200Highly modified water courses and canals31000
32200Highly artificial coniferous plantations211291300Separated water bodies belonging to the river system31000
33100Natural and seminatural mixed forest99992100Natural water bodies31000
33200Highly artificial mixed plantations99992200Artificial standing water bodies31000
34100Transitional woodland and scrub9592300Intensively managed fish ponds31000
34200Lines of trees and scrub9592400Standing water bodies of extractive industrial sites31000
35000Damaged forest9510000Sea and ocean31000
Table A4. Reclassification table for CLC.
Table A4. Reclassification table for CLC.
CLCClass NameLCCLCClass NameLC
111Continuous urban fabric99831113117Broad-leaved forest2111
112Discontinuous urban fabric99831213125Needle-leaved forest2112
121Industrial or commercial units99831313132Mixed forest999
1211Photovoltaic fiends9983211Continuous natural grasslands22200
122Road and rail networks and associated land9983212Discontinuous natural grasslands22200
123Port areas998322Moors and heathland21220
124Airports9983231High Mediterranean scrub21220
131Mineral extraction sites9983232Low Mediterranean scrub and the garrigue21220
132Dump sites998324Transitional woodland shrub95
133Construction sites9983241Forest harvesting95
141Green urban areas91331Beaches, dunes, sands12200
142Sport and leisure facilities998332Bare rocks12100
2111Intensive nonirrigated arable land22120333Sparsely vegetated areas96
2112Extensive nonirrigated arable land22120334Burnt areas97
212Permanently irrigated land22120335Glaciers and perpetual snow32000
213Rice fields22120411Inland marshes40000
221Vineyards21210412Peat bogs40000
222Fruit trees and berry plantations21131421Salt marshes40000
223Olive groves21132422Salines40000
224Woody plantation21133423Intertidal flats40000
2241New woody plantation21133511Water courses31000
231Pastures22110512Water bodies31000
241Annual crops associated with permanent crops92521Coastal lagoons31000
242Complex cultivation patterns92522Estuaries31000
243Land principally occupied by agriculture, with significant areas of natural vegetation93523Sea and ocean31000
244Agroforestry94

References

  1. Munafò, M. Consumo di Suolo, Dinamiche Territoriali e Servizi Ecosistemici—Edizione 2018; ISPRA: Rome, Italy, 2018. [Google Scholar]
  2. Verburg, P.H.; Erb, K.-H.; Mertz, O.; Espindola, G. Land System Science: Between global challenges and local realities. Curr. Opin. Environ. Sustain. 2013, 5, 433–437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pedroli, G.B.M.; Meiner, A. Landscapes in Transition: An Account of 25 Years of Land Cover Change in Europe; European Environment Agency: København, Denmark, 2017. [Google Scholar]
  4. EEA. Land Systems at European Level—Analytical Assessment Framework. 2018. Briefing no. 10/2018. pp. 1–9. Available online: https://www.eea.europa.eu/themes/landuse/land-systems (accessed on 28 November 2021).
  5. Anderson, J.R. Land use and land cover changes. A framework for monitoring. J. Res. Geol. Surv. 1977, 5, 143–153. [Google Scholar]
  6. Liu, R.; Dong, X.; Wang, X.; Zhang, P.; Liu, M.; Zhang, Y. Study on the relationship among the urbanization process, ecosystem services and human well-being in an arid region in the context of carbon flow: Taking the Manas river basin as an example. Ecol. Indic. 2021, 132, 108248. [Google Scholar] [CrossRef]
  7. Costanza, R.; D’Arge, R.; de Groot, R.; Farberll, S.; Grassot, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’Neilltt, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 387115, 253–260. [Google Scholar] [CrossRef]
  8. Gómez-Baggethun, E.; De Groot, R.; Lomas, P.L.; Montes, C. The history of ecosystem services in economic theory and practice: From early notions to markets and payment schemes. Ecol. Econ. 2010, 69, 1209–1218. [Google Scholar] [CrossRef]
  9. Danley, B.; Widmark, C. Evaluating conceptual definitions of ecosystem services and their implications. Ecol. Econ. 2016, 126, 132–138. [Google Scholar] [CrossRef]
  10. Teague, A.; Russell, M.; Harvey, J.; Dantin, D.; Nestlerode, J.; Alvarez, F. A spatially-explicit technique for evaluation of alternative scenarios in the context of ecosystem goods and services. Ecosyst. Serv. 2016, 20, 15–29. [Google Scholar] [CrossRef]
  11. Corvalan, C.; Hales, S.; McMichael, A.; Butler, C.; Campbell-Lendrum, D.; Confalonieri, U.; Leitner, K.; Lewis, N.; Patz, J.; Polson, K.; et al. A Report of the Millennium Ecosystem Assessment. In Ecosystems and Human Well-Being; Island Press: Washington, DC, USA, 2005; Volume 5. [Google Scholar]
  12. Xepapadeas, A. The Economics of Ecosystems and Biodiversity: Ecological and Economic Foundations; Kumar, P., Ed.; Earthscan: London, UK; Washington, DC, USA, 2011; Volume 16, pp. 239–242. ISBN 978-1-84971-212-5 (HB). [Google Scholar]
  13. La Notte, A.; D’Amato, D.; Mäkinen, H.; Paracchini, M.L.; Liquete, C.; Egoh, B.; Geneletti, D.; Crossman, N.D. Ecosystem services classification: A systems ecology perspective of the cascade framework. Ecol. Indic. 2017, 74, 392–402. [Google Scholar] [CrossRef]
  14. Potschin, M.; Haines-Young, R. Defining and measuring ecosystem services. In Routledge Handbook of Ecosystem Services; Potschin, M., Haines-Young, R., Fish, R., Turner, R.K., Eds.; Routledge: London, UK; New York, NY, USA, 2016; pp. 25–44. [Google Scholar]
  15. 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]
  16. Zulian, G.; Maes, J.; Paracchini, M.L. Linking land cover data and crop yields for mapping and assessment of pollination services in Europe. Land 2013, 2, 472–492. [Google Scholar] [CrossRef] [Green Version]
  17. Groff, S.C.; Loftin, C.S.; Drummond, F.; Bushmann, S.; McGill, B. Parameterization of the InVEST crop pollination model to spatially predict abundance of wild blueberry (Vaccinium angustifolium Aiton) native bee pollinators in Maine, USA. Environ. Model. Softw. 2016, 79, 1–9. [Google Scholar] [CrossRef] [Green Version]
  18. Leonhardt, S.D.; Gallai, N.; Garibaldi, L.A.; Kuhlmann, M.; Klein, A.-M. Economic gain, stability of pollination and bee diversity decrease from southern to northern Europe. Basic Appl. Ecol. 2013, 14, 461–471. [Google Scholar] [CrossRef]
  19. Lautenbach, S.; Seppelt, R.; Liebscher, J.; Dormann, C.F. Spatial and temporal trends of global pollination benefit. PLoS ONE 2012, 7, e35954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Hanley, N.; Breeze, T.D.; Ellis, C.; Goulson, D. Measuring the economic value of pollination services: Principles, evidence and knowledge gaps. Ecosyst. Serv. 2015, 14, 124–132. [Google Scholar] [CrossRef] [Green Version]
  21. UK National Ecosystem Assessment. The UK National Ecosystem Assessment: Synthesis of the Key Findings; UNEP-WCMC: Cambridge, UK, 2011. [Google Scholar]
  22. Rabe, S.-E.; Koellner, T.; Marzelli, S.; Schumacher, P.; Grêt-Regamey, A. National ecosystem services mapping at multiple scales. The German exemplar. Ecol. Indic. 2016, 70, 357–372. [Google Scholar] [CrossRef]
  23. Santos-Martín, F.; García Llorente, M.; Quintas-Soriano, C.; Zorrilla-Miras, P.; Martín-Lopez, B.; Loureiro, M.; Benayas, J.; Montes, M. Spanish National Ecosystem Assessment: Socio-economic valuation of ecosystem services in Spain. Synthesis of the key findings; Biodiversity Foundation of the Spanish Ministry of Agriculture, Food and Environment: Madrid, Spain, 2016.
  24. Xie, G.; Zhang, C.; Zhen, L.; Zhang, L. Dynamic changes in the value of China’s ecosystem services. Ecosyst. Serv. 2017, 26, 146–154. [Google Scholar] [CrossRef]
  25. Turpie, J.K.; Forsythe, K.J.; Knowles, A.; Blignaut, J.; Letley, G. Mapping and valuation of South Africa’s ecosystem services: A local perspective. Ecosyst. Serv. 2017, 27, 179–192. [Google Scholar] [CrossRef]
  26. Choi, H.-A.; Song, C.; Lee, W.-K.; Jeon, S.; Gu, J.H. Integrated approaches for national ecosystem assessment in South Korea. KSCE J. Civ. Eng. 2018, 22, 1634–1641. [Google Scholar] [CrossRef]
  27. Daily, G.C.; Matson, P.A.; Vitousek, P.M. Ecosystem services supplied by soil. In Nature’s Services: Societal Dependence on Natural Ecosystems; Daily, G.C., Ed.; Island Press: Washington, DC, USA, 1997; pp. 113–132. [Google Scholar]
  28. Wall, D.H.; Bardgett, R.D.; Covich, A.P.; Snelgrove, P.V.R. The Need for Understanding How Biodiversity and Ecosystem Functioning Affect Ecosystem Services in Soils and Sediments; Island Press: Washington, DC, USA, 2004. [Google Scholar]
  29. 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]
  30. Fryer, J.; Williams, I.D. Regional carbon stock assessment and the potential effects of land cover change. Sci. Total Environ. 2021, 775, 145815. [Google Scholar] [CrossRef] [PubMed]
  31. Mathew, I.; Shimelis, H.; Mutema, M.; Minasny, B.; Chaplot, V. Crops for increasing soil organic carbon stocks–A global meta analysis. Geoderma 2020, 367, 114230. [Google Scholar] [CrossRef]
  32. Clerici, N.; Cote-Navarro, F.; Escobedo, F.J.; Rubiano, K.; Villegas, J.C. Spatio-temporal and cumulative effects of land use-land cover and climate change on two ecosystem services in the Colombian Andes. Sci. Total Environ. 2019, 685, 1181–1192. [Google Scholar] [CrossRef]
  33. Savaresi, A.; Perugini, L.; Chiriacò, M.V. Making sense of the LULUCF Regulation: Much ado about nothing? Rev. Eur. Comp. Int. Environ. Law 2020, 29, 212–220. [Google Scholar] [CrossRef]
  34. Sallustio, L.; Quatrini, V.; Geneletti, D.; Corona, P.; Marchetti, M. Assessing land take by urban development and its impact on carbon storage: Findings from two case studies in Italy. Environ. Impact Assess. Rev. 2015, 54, 80–90. [Google Scholar] [CrossRef] [Green Version]
  35. Vashum, K.T.; Jayakumar, S. Methods to estimate above-ground biomass and carbon stock in natural forests-a review. J. Ecosyst. Ecography 2012, 2, 1–7. [Google Scholar] [CrossRef]
  36. Sallustio, L.; Munafò, M.; Riitano, N.; Lasserre, B.; Fattorini, L.; Marchetti, M. Integration of land use and land cover inventories for landscape management and planning in Italy. Environ. Monit. Assess. 2016, 188, 48. [Google Scholar] [CrossRef]
  37. Maes, J.; Teller, A.; Erhard, M.; Liquete, C.; Braat, L.; Berry, P.; Egoh, B.; Puydarrieus, P.; Fiorina, C.; Santos, F.; et al. An Analytical Framework for Ecosystem Assessments under Action 5 of the EU Biodiversity Strategy to 2020; Publications Office of the European Union: Luxembourg, 2013; ISBN 9789279293696. [Google Scholar]
  38. Cabral, P.; Feger, C.; Levrel, H.; Chambolle, M.; Basque, D. Assessing the impact of land-cover changes on ecosystem services: A first step toward integrative planning in Bordeaux, France. Ecosyst. Serv. 2016, 22, 318–327. [Google Scholar] [CrossRef] [Green Version]
  39. Arnold, S.; Kosztra, B.; Banko, G.; Smith, G.; Hazeu, G.; Bock, M. The EAGLE concept—A vision of a future European Land Monitoring Framework. In Proceedings of the 33th EARSeL Symposium towards Horizon 2020, Matera, Italy, 3–6 June 2013; pp. 551–568. [Google Scholar]
  40. Kleeschulte, S. Technical Specifications for Implementation of a New Land-Monitoring Concept Based on EAGLE. D3: Draft Design Concept and CLC-Backbone, CLC-Core Technical Specifications, Including Requirements Review, 3rd ed.; European Environment Agency: Copenhagen, Denmark, 2019; p. 78. [Google Scholar]
  41. De Fioravante, P.; Luti, T.; Cavalli, A.; Giuliani, C.; Dichicco, P.; Marchetti, M.; Chirici, G.; Congedo, L.; Munafò, M. Multispectral Sentinel-2 and SAR Sentinel-1 Integration for Automatic Land Cover Classification. Land 2021, 10, 611. [Google Scholar] [CrossRef]
  42. Munafò, M. Consumo di Suolo, Dinamiche Territoriali e Servizi Ecosistemici. Edizione 2020; ISPRA: Rome, Italy, 2020. [Google Scholar]
  43. Munafò, M. Consumo di Suolo, Dinamiche Territoriali e Servizi Ecosistemici Edizione 2021 Rapporto ISPRA SNPA; ISPRA: Rome, Italy, 2021; ISBN 9788844810597. [Google Scholar]
  44. Luti, T.; De Fioravante, P.; Marinosci, I.; Strollo, A.; Riitano, N.; Falanga, V.; Mariani, L.; Congedo, L.; Munafò, M. Land Consumption Monitoring with SAR Data and Multispectral Indices. Remote Sens. 2021, 13, 1586. [Google Scholar] [CrossRef]
  45. Strollo, A.; Smiraglia, D.; Bruno, R.; Assennato, F.; Congedo, L.; De Fioravante, P.; Giuliani, C.; Marinosci, I.; Riitano, N.; Munafò, M. Land consumption in Italy. J. Maps 2020, 16, 113–123. [Google Scholar] [CrossRef]
  46. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  47. FAO Map Accuracy Assessment and Area Estimation: A Practical Guide. National Forest Monitoring and Assessment Working Paper. 2016. 46/E. p. 69. Available online: https://www.fao.org/publications/card/fr/c/e5ea45b8-3fd7-4692-ba29-fae7b140d07e/ (accessed on 28 November 2021).
  48. Stehman, S.V.; Czaplewski, R.L. Design and Analysis for Thematic Map Accuracy Assessment—An application of satellite imagery. Remote Sens. Environ. 1998, 64, 331–344. [Google Scholar] [CrossRef]
  49. Stehman, S.V.; Foody, G.M. Key issues in rigorous accuracy assessment of land cover products. Remote Sens. Environ. 2019, 231, 111199. [Google Scholar] [CrossRef]
  50. Cochran, W.G.; William, G. Sampling Techniques; John Wiley& Sons: New York, NY, USA, 1977. [Google Scholar]
  51. Lupia, F. MARSALa-A Model-Based Irrigation Water Consumption Estimation at Farm Level; INEA: Brussels, Belgium, 2013. [Google Scholar]
  52. Vitullo, M.; De Lauretis, R.; Federici, S. La Contabilità del Carbonio Contenuto Nelle Foreste Italiane; Agenzia per la Protezione dell’ambiente e per i Servizi Tecnici: Rome, Italy, 2007. [Google Scholar]
  53. Di Cosmo, L.; Gasparini, P.; Tabacchi, G. A national-scale, stand-level model to predict total above-ground tree biomass from growing stock volume. For. Ecol. Manag. 2016, 361, 269–276. [Google Scholar] [CrossRef]
  54. Romano, D.; Arcarese, C.; Bernetti, A.; Caputo, A.; Contaldi, M.; De Lauretis, R.; Di Cristofaro, E.; Gagna, A.; Gonella, B.; Taurino, E.; et al. Italian Greenhouse Gas Inventory 1990–2015. National Inventory Report 2017. 2017. Available online: https://www.isprambiente.gov.it/en/publications/reports/italian-greenhouse-gas-inventory-1990-2015.-national-inventory-report-2017 (accessed on 28 November 2021).
  55. Li, S.; Liu, Y.; Yang, H.; Yu, X.; Zhang, Y.; Wang, C. Integrating ecosystem services modeling into effectiveness assessment of national protected areas in a typical arid region in China. J. Environ. Manag. 2021, 297, 113408. [Google Scholar] [CrossRef]
  56. Penman, J.; Gytarsky, M.; Hiraishi, T.; Krug, T.; Kruger, D.; Pipatti, R.; Buendia, L.; Miwa, K.; Ngara, T.; Tanabe, K. Good practice guidance for land use, land-use change and forestry. Good Pract. Guid. L. Use Land-Use Chang. For. 2003. Available online: https://www.ipcc.ch/publication/good-practice-guidance-for-land-use-land-use-change-and-forestry/ (accessed on 28 November 2021).
  57. Marchetti, M.; Sallustio, L.; Ottaviano, M.; Barbati, A.; Corona, P.; Tognetti, R.; Zavattero, L.; Capotorti, G. Carbon sequestration by forests in the National Parks of Italy. Plant Biosyst. Int. J. Deal. Asp. Plant Biol. 2012, 146, 1001–1011. [Google Scholar] [CrossRef]
  58. Canaveira, P.; Manso, S.; Pellis, G.; Perugini, L.; De Angelis, P.; Neves, R.; Papale, D.; Paulino, J.; Pereira, T.; Pina, A. Biomass Data on Cropland and Grassland in the Mediterranean Region. Final Report for Action A4 of Project MediNet. 2018. Available online: http://www.lifemedinet.com/ (accessed on 28 November 2021).
  59. Burke, T.; Whyatt, J.D.; Rowland, C.; Blackburn, G.A.; Abbatt, J. The influence of land cover data on farm-scale valuations of natural capital. Ecosyst. Serv. 2020, 42, 101065. [Google Scholar] [CrossRef]
  60. Tomlinson, S.J.; Dragosits, U.; Levy, P.E.; Thomson, A.M.; Moxley, J. Quantifying gross vs. net agricultural land use change in Great Britain using the Integrated Administration and Control System. Sci. Total Environ. 2018, 628, 1234–1248. [Google Scholar] [PubMed]
  61. Morel, J.L.; Chenu, C.; Lorenz, K. Ecosystem services provided by soils of urban, industrial, traffic, mining, and military areas (SUITMAs). J. Soils Sediments 2015, 15, 1659–1666. [Google Scholar] [CrossRef]
  62. Sambucini, V.; Marinosci, I.; Bonora, N.; Chirci, G. La Realizzazione in Italia del Progetto CORINE Land Cover 2006; ISPRA: Rome, Italy, 2010; ISBN 9788844804770. [Google Scholar]
  63. Tomaselli, V.; Dimopoulos, P.; Marangi, C.; Kallimanis, A.S.; Adamo, M.; Tarantino, C.; Panitsa, M.; Terzi, M.; Veronico, G.; Lovergine, F. Translating land cover/land use classifications to habitat taxonomies for landscape monitoring: A Mediterranean assessment. Landsc. Ecol. 2013, 28, 905–930. [Google Scholar] [CrossRef] [Green Version]
  64. Spadoni, G.L.; Cavalli, A.; Congedo, L.; Munafò, M. Analysis of Normalized Difference Vegetation Index (NDVI) multi-temporal series for the production of forest cartography. Remote Sens. Appl. Soc. Environ. 2020, 20, 100419. [Google Scholar] [CrossRef]
Figure 1. Workflow of the methodology for the production of the land cover map and evaluation of ecosystem services. The mosaic between Urban Atlas, Riparian Zones, and Natura 2000 was created (projected and rasterized at 10 m), then supplemented with the regional maps and CORINE Land Cover in the areas not covered by these three data. The distinction of forest categories and mixed classes was based on HRLs, while the National Land Consumption Map for 2012 and 2020 made it possible to identify artificial abiotic surfaces and their variation. The map allowed the calculation of three of the four pools considered for the estimation of the variation in carbon storage capacity.
Figure 1. Workflow of the methodology for the production of the land cover map and evaluation of ecosystem services. The mosaic between Urban Atlas, Riparian Zones, and Natura 2000 was created (projected and rasterized at 10 m), then supplemented with the regional maps and CORINE Land Cover in the areas not covered by these three data. The distinction of forest categories and mixed classes was based on HRLs, while the National Land Consumption Map for 2012 and 2020 made it possible to identify artificial abiotic surfaces and their variation. The map allowed the calculation of three of the four pools considered for the estimation of the variation in carbon storage capacity.
Land 11 00035 g001
Figure 2. Study area—Italy.
Figure 2. Study area—Italy.
Land 11 00035 g002
Figure 3. Land cover map of Italy for 2012.
Figure 3. Land cover map of Italy for 2012.
Land 11 00035 g003
Figure 4. Result of the accuracy assessment with reference to the sample of photo-interpreted points.
Figure 4. Result of the accuracy assessment with reference to the sample of photo-interpreted points.
Land 11 00035 g004
Figure 5. Carbon stock (2012).
Figure 5. Carbon stock (2012).
Land 11 00035 g005
Figure 6. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on urban area.
Figure 6. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on urban area.
Land 11 00035 g006
Figure 7. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on water courses.
Figure 7. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on water courses.
Land 11 00035 g007
Figure 8. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on small patches of natural vegetation and agricultural areas.
Figure 8. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on small patches of natural vegetation and agricultural areas.
Land 11 00035 g008
Figure 9. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on heterogeneous agricultural areas.
Figure 9. Comparison between land cover map for 2012 (a), CORINE Land Cover (b) and satellite image (c), focusing on heterogeneous agricultural areas.
Land 11 00035 g009
Table 1. Land cover classification system.
Table 1. Land cover classification system.
Land Cover
I LevelII LevelIII LevelIV LevelV Level
1Abiotic
Non-vegetated
surfaces
11Artificial abiotic
12Natural abiotic121Consolidated (bare rocks, cliffs)
122Unconsolidated (beaches, dunes, sands)
2Biotic vegetated surfaces21Woody
vegetation
211Trees2111Broad-leaved
2112Needle-leaved
2113Permanent crops21131Orchards
21132Olive groves
21133Wood plantations
212Shrubs2121Vineyards
2122Shrubland
22Herbaceous
vegetation
221Periodically2211Pastures
2212Arable land
222Permanent
3Water surfaces31Water bodies
32Permanent snow and ice
4Wetlands
Table 2. Copernicus Land Monitoring Service—local component data.
Table 2. Copernicus Land Monitoring Service—local component data.
Urban AtlasRiparian ZonesNatura 2000
Data typeVectorVectorVector
Classes275655
MMU0.25 ha (class 1)
1 ha (class 2–5)
0.5 ha0.5 ha
Table 3. Copernicus Land Monitoring Service—Pan-European component data.
Table 3. Copernicus Land Monitoring Service—Pan-European component data.
High Resolution LayersCORINE Land Cover
Data typeRasterVector
Classes444
MMUPixel 10 × 10 m25 ha (status)
5 ha (changes)
Table 4. Temporary codes for mixed land cover classes, artificial surfaces and mixed forests.
Table 4. Temporary codes for mixed land cover classes, artificial surfaces and mixed forests.
Temporary CodesUARZN2000CLC
91Green urban areas 14100 141
92Annual crops associated with permanent crops 231002310241
Complex cultivation patterns24000232002320242
93Land principally occupied by agriculture with significant areas of natural vegetation 233002330243
94Agroforestry 234002340244
95Transitional woodland and scrub 341003410324
Damaged forest 3500
96Sparsely vegetated areas33000610006100333
97Burnt areas 632006320334
99Permanent crops (vineyards, fruit trees, olive groves)22000221002210
998Urban areas and artificial surfaces11100–13400, 1420011110–140001110–1400111–133, 142
999Mixed forest3100033100 3131–3132
Table 5. Attribution of classes to temporary codes using HRL. DLT = dominant leaf type, GRA = grassland.
Table 5. Attribution of classes to temporary codes using HRL. DLT = dominant leaf type, GRA = grassland.
Temporary CodeArea Covered by DLTArea Covered by GRANo GRA and No DLT Area
9121131Orchards2211Pastures2121Vineyards
922212Arable land
932111
or
2112
Broad-leaved
or
Needle-leaved
2211Pastures2212Arable land
942212Arable land
952122Shrubland
96222Permanent
97222Permanent121Natural abiotic surfaces
9921131Orchards2121Vineyards
Table 6. Calculation of sample size.
Table 6. Calculation of sample size.
Land Cover ClassArea (ha)WiWi*Si
Artificial abiotic surfaces2,102,2880.0700.034
Consolidated (bare rocks, cliffs)802,5390.0270.013
Unconsolidated (beaches, dunes, sands)59,3250.0020.001
Broad-leaved7,805,2170.2590.127
Needle-leaved2,167,2610.0720.035
Orchards619,7590.0210.010
Olive plantations1,088,0990.0360.018
Wood plantations46,3560.0020.001
Vineyards689,7460.0230.011
Shrubland1,375,1250.0460.022
Pastures1,372,7900.0460.022
Arable land9,183,5190.3050.149
Permanent herbaceous2,337,8260.0780.038
Water bodies402,8300.0130.007
Permanent snow and ice37,0000.0010.001
Wetlands50,2950.0020.001
Total30,139,975 1
Total number of samples2400
Table 7. Sample allocation.
Table 7. Sample allocation.
ClassesAllocation
EqualProportionalFinal
Artificial abiotic surfaces150167159
Consolidated(bare rocks, cliffs)15064107
Unconsolidated (beaches, dunes, sands)1505100
Broad-leaved150622386
Needle-leaved150173162
Orchards15049100
Olive plantations15087119
Wood plantations150477
Vineyards15055103
Shrubland150109130
Pastures150109130
Arable land150731441
Permanent herbaceous150186168
Water bodies1503291
Permanent snow and ice150377
Wetlands150477
Total240024002427
Table 8. Carbon content for land cover classes.
Table 8. Carbon content for land cover classes.
LC ClassAGBBGBDOSTotal
M g   C h a
Agricultural areas5--58.1
Orchards10--62.1
Woody plantations28.555.251.7599.45
Pastures---78.9
Other natural areas3.05--69.95
Urban areas----
Sparsely vegetated areas----
Table 9. Carbon stock values for some agricultural classes.
Table 9. Carbon stock values for some agricultural classes.
LC ClassAGBBGB
M g   C h a
Olive trees9.12.6
Vineyards5.54.4
Fruit trees8.35.6
Table 10. Overall, user and producer accuracy of the 2012 land cover map.
Table 10. Overall, user and producer accuracy of the 2012 land cover map.
Overall Accuracy
0.83
Class nameUser’s
accuracy
Producer’s
accuracy
Artificial abiotic surfaces0.980.83
Consolidated (bare rocks scree, cliffs)0.890.82
Unconsolidated (beaches, dunes, sands)0.970.85
Broad-leaved0.920.91
Needle-leaved0.880.96
Orchards0.820.86
Olive plantation0.980.94
Wood plantations0.900.84
Vineyards0.790.89
Shrubland0.840.83
Pastures0.880.93
Arable land0.920.93
Permanent herbaceous0.760.83
Water bodies0.900.95
Permanent snow and ice0.990.79
Wetlands0.880.79
Table 11. Land cover classes 2012 (first and second classification level).
Table 11. Land cover classes 2012 (first and second classification level).
ha% Total% Class
Abiotic Surfaces2,964,1519.83
Artificial abiotic surfaces802,5392.6627.07
Natural abiotic surfaces59,3250.202.00
Bioticvegetated26,685,69688.54
Woody vegetation13,791,56145.7651.68
Herbaceous vegetation12,894,13542.7848.32
Water surfaces439,8301.46
Water bodies402,8301.3491.59
Permanent snow and ice37,0000.128.41
Wetlands50,2950.17
Table 12. Land cover classes 2012 (up to fifth classification level).
Table 12. Land cover classes 2012 (up to fifth classification level).
LC CodeClass Nameha%
11000Artificial abiotic surfaces2,102,2886.98
12100Consolidated (bare rocks scree, cliffs)802,5392.66
12200Unconsolidated (beaches, dunes, sands)59,3250.20
21110Broad-leaved7,805,21725.90
21120Needle-leaved2,167,2617.19
21131Orchards619,7592.06
21132Olive plantations1,088,0993.61
21133Wood plantations46,3560.15
21210Vineyards689,7462.29
21220Shrubland1,375,1254.56
22110Pastures1,372,7904.55
22120Arable land9,183,51930.47
22200Permanent herbaceous2,337,8267.76
31000Water bodies402,8301.34
32000Permanent snow and ice37,0000.12
40000Wetlands50,2950.17
Total 30,139,972100.00
Table 13. Variation in carbon storage capacity (2012–2020).
Table 13. Variation in carbon storage capacity (2012–2020).
RegionC Storage Variation (2012–2020)
(Tons)(%)
Abruzzo86.0342.97
Basilicata52.0861.80
Calabria102.1383.52
Campania200.6616.92
Emilia-Romagna274.2059.46
Friuli-Venezia Giulia91.2223.15
Latium203.8187.03
Liguria25.2370.87
Lombardy319.66711.03
Marche81.7462.82
Molise21.2430.73
Piedmont203.7257.03
Apulia210.7237.27
Sardinia88.8913.07
Sicily225.9727.80
Tuscany120.1434.14
Trentino-South Tyrol134.4704.64
Umbria58.9472.03
Aosta Valley13.2070.46
Veneto384.53713.27
Italy2.898.672100.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Fioravante, P.; Strollo, A.; Assennato, F.; Marinosci, I.; Congedo, L.; Munafò, M. High Resolution Land Cover Integrating Copernicus Products: A 2012–2020 Map of Italy. Land 2022, 11, 35. https://doi.org/10.3390/land11010035

AMA Style

De Fioravante P, Strollo A, Assennato F, Marinosci I, Congedo L, Munafò M. High Resolution Land Cover Integrating Copernicus Products: A 2012–2020 Map of Italy. Land. 2022; 11(1):35. https://doi.org/10.3390/land11010035

Chicago/Turabian Style

De Fioravante, Paolo, Andrea Strollo, Francesca Assennato, Ines Marinosci, Luca Congedo, and Michele Munafò. 2022. "High Resolution Land Cover Integrating Copernicus Products: A 2012–2020 Map of Italy" Land 11, no. 1: 35. https://doi.org/10.3390/land11010035

APA Style

De Fioravante, P., Strollo, A., Assennato, F., Marinosci, I., Congedo, L., & Munafò, M. (2022). High Resolution Land Cover Integrating Copernicus Products: A 2012–2020 Map of Italy. Land, 11(1), 35. https://doi.org/10.3390/land11010035

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