Changing landscapes: habitat monitoring and land transformation in a long-time used Mediterranean coastal wetland

Zone Umide della Capitanata, located in the northeastern part of the Apulia Region, is one of the most extensive coastal wetlands of the Italian peninsula and one of the largest components of the Mediterranean wetland system. Despite its high ecological importance, this site has been undergoing a variety of pressures intensified in recent decades. This study analyzes and evaluates the changes occurred in this area between 2010 and 2020. Land cover and habitat maps were performed by photointerpretation and on-site surveys, and classified according to the FAO-LCCS and EUNIS taxonomies, respectively. To focus on local dynamics, four subset areas were analyzed separately. A set of landscape metrics was computed to analyze the landscape structure. The anthropogenic pressures affecting the study area were described through the Driving Forces-Pressures-State-Impact-Response (DPSIR) conceptual framework. Numerous changes were identified, deriving from both finest thematic redefinition and anthropogenic pressures. Both class conversions and class modifications were identified and quantified by means of transition matrices. Most of the observed conversions were borne by classes belonging to saltmarshes and to coastal dune systems. In particular, landscape configuration of coastal dune classes was well highlighted by a set of specific landscape metrics. Agriculture practices and changes in water flow pattern turned out to be the main driving forces exerting pressures on these natural systems. Significant differences were found between the four subsets under analysis, thus indicating that different management strategies lead to different levels of conservation.


Introduction
Coastal environments have a fundamental role in keeping and enhancing a wide range of ecosystem services, such as improvement in water quality, equilibration of the water cycle, carbon sequestration, natural habitat for migratory birds, recreation (Maltby and Acreman 2011;Mehvar et al. 2018). Coastal wetlands and sandy dune systems are characterized by high biological diversity, which includes several habitats of European Community interest, as well as plant and animal species listed in the Annexes to 92/43/EEC Habitat Directive (HD). These environments are also characterized by a great fragility due to different threat factors. Land claim, agricultural intensification, hydrological modifications, urbanization, coastal erosion, and introduction of alien species are the main drivers of change (Davy et al. 2009;Janssen et al. 2016). In addition, climate changes, such as increase in temperatures and/or decrease in precipitation during the spring and summer seasons (Giorgi and Bi 2005), and sea level rise (Janssen et al. 2016), may affect coastal plant communities and habitats, determining changes in their distribution and extension. In the Mediterranean region, coastal wetlands and dune systems have progressively undergone increasing pressures, with resulting habitat degradation, loss and fragmentation (Perennou et al. 2012;Bazzichetto et al. 2020;Margiotta et al. 2020;Rodríguez-Santalla and Navarro 2021). Moreover, complex, diverse and highly fragmented landscape patterns have been determined by a variety of historical land uses in whole the Mediterranean (Naveh 1990;Ruiz and Sanz-Sanchez 2020).
In this scenario, constant and consistent monitoring procedures, as well as awareness of the ongoing socio-economic dynamics, are fundamental to implement effective management policies and conservation strategies. Assessing changes in landscape ecological elements through land-cover and habitat monitoring over long time periods is essential to: (a) understand the drivers of the temporal changes; (b) make provisions for future trends; (c) design appropriate conservation policies (Fahrig 2003;Fisher and Lindenmayer 2007;Bunce et al. 2008;Lengyel et al. 2008;Nagendra et al. 2014).
The representation of the detected changes, that may occur in terms of both conversion and modification, is a crucial step and is closely depending on the Land Cover (LC) or habitat taxonomy used, that is to say, to its structure, semantics and especially to its level of thematic resolution (Jansen and Di Gregorio 2002;Bajocco et al. 2012;Tomaselli et al. 2016Tomaselli et al. , 2021. Quantifying the landscape's spatial structure provides an understanding of the underlying impact on ecological processes and of the effects that changing patterns have on ecosystem services provision (Turner et al. 2007). In this framework, Landscape Metrics (LM) are an essential tool to analyze and evaluate landscape mosaics and spatial arrangement of the landscape structure (Turner et al. 2001;Uuemaa et al. 2013). The application of LM in coastal areas can be relevant because these landscapes are prone to rapid transformations, due to numerous anthropogenic activities, with natural areas converted into other land uses or in other natural types, and vice-versa. Many analyses have been carried out to identify core sets of metrics for the assessment of landscape heterogeneity and the relationships with anthropogenic activities (Botequilha-Leitão and Ahren 2002; Nagendra et al. 2004;Schindler et al. 2008;Plexida et al. 2014). The analysis of a landscape pattern includes two main aspects: (a) composition, relating to the abundance and variety of patch types in the landscape; (b) configuration, relating to spatial arrangement and context of the patches (Riitters et al. 1995).
The definition of appropriate conservation policies requires a precise understanding of the driving forces and pressures influencing the environment. The Driver-Pressure-State-Impact-Response (DPSIR) model, developed by the European Environment Agency (EEA 1995) as a casual framework for describing the interactions between society and environment, has been widely used as analytical approach for determining pressures and impacts especially in coastal environments. DPSIR model is considered as an additional tool for policy makers concerning environmental resources management (Lozoya et al. 2011;Gari et al. 2015;Bruno et al. 2020).
In this paper, we focused on Capitanata (FG), one of the largest components of the Mediterranean wetland system that has historically been subject to various types of land use and that, despite its high ecological importance, has been undergoing a variety of pressures, which have intensified in recent decades and particularly in recent years.
The aim of this work is to detect, analyze, and evaluate the changes occurred in the study area between 2010 and 2020, in order to: (a) identify the main changes occurred over this time-frame, in a landscape that has been undergoing rapid transformations; (b) evaluate the effectiveness of selected landscape metrics in analyzing the landscape changes; (c) correlate driving forces, pressures and threats to the identified changes.

Study area
The site "Zone Umide della Capitanata e Paludi presso il Golfo di Manfredonia" [Site of Community Importance (SCI) IT9110005; Special Protection Area (SPA) IT9110038], is located in the north-eastern part of the Apulia Region (SE Italy). It includes three protected areas: the State Natural Reserves "Palude Frattarolo" and "Saline di Margherita di Savoia", and the Oasis "Lago Salso", enclosed in the Gargano National Park (Fig. 1). This site, extended more than 14.000 ha, is one of the most extensive wetlands of the Italian peninsula and one of the largest components of the Mediterranean wetland system, classified as Ramsar site and Important Bird Area (IBA 230 M).
The Capitanata's landscape consists of a system of lagoons, with brackish to salt water (depending on the specific water regime), fresh water wetlands, swamps and ponds, surrounded by cultivated field and farmlands. The agricultural landscape is characterized by intensive horticultural crop fields, locally called "arenili", widely extended along the coast and mainly used to cultivate onions, potatoes and carrots; towards the inland, by graminoid crops, in some case interspersed with tree crops (olive groves, vineyards, orchards).
The long-term exploitation of the site, mainly for agricultural purposes, has led to a progressive reduction and fragmentation of the original natural habitats. The natural vegetation is mostly represented by halophytic shrub, annual pioneer communities, reed thickets, and rush and sedges communities, typical of salt marshes. Along the sandy coast, fragments of dune vegetation are still present, although strongly reduced and altered due to anthropic activities and coastal erosion (Tomaselli and Sciandrello 2017). According to the erosion risk map of Apulian coasts (Bruno et al. 2020), the Capitanata falls within low (northern part) to high (southern part) erosion risk level.
Despite the whole of these negative processes, the study area still represents an important site for avifauna, especially water bird species, which find nesting sites or rest areas during the annual migrations.
In consideration of the important conservation value of these coastal wetlands and their delicate and fragile ecological equilibrium, two Life+ Nature and Biodiversity projects have been carried on in the last years: "Sipontine Wetlands"-Conservation actions of habitats in the coastal wetlands of SCI "Zone Umide della Capitanata" (www. lifez oneum ide. it), and "Conservation Activities for Priority Avifauna in the Lago Salso Oasis" (www. lifel agosa lso. it).
LC and habitat mapping LC and habitat mapping in 2020 were performed by means of photointerpretation and on-site surveys, and compared to an analogous set of maps referring to 2010. The thematic maps were produced in Arc-Gis 10.2 by digitizing color orthophotos, produced between 2019 and 2020 within the POR-PUGLIA project COHECO (www. coheco. it), in three different months (February, June and October). The same criteria of a previous mapping carried out in 2010 were followed, taking care of spatial co-registration between the two maps, with appropriate corrections and modifications, in order to avoid mismatching problems or errors in the estimation of changes. First, natural and semi-natural landscape elements were described as vegetation types defined on the base of phytosociological units, according to the Zurich-Montpellier method (Braun-Blanquet 1964). We adopted a representation scale 1:5000, which allowed representing the studied landscapes with a 2 m resolution. Next, vegetation units were reclassified in habitat types and subsequently in LC classes. Habitat mapping was performed using the EUropean Nature Information System (EUNIS) (Davies et al. 2004) classification scheme (levels III and IV), which is considered an effective standardizing tool for habitat classification in the European Union (EU) (Ichter et al. 2014). LC classes were defined on the base of the FAO-Land Cover Classification System (LCCS) taxonomy (Di Gregorio and Jansen 2005) whose potentialities in mapping Mediterranean coastal wetlands have been explored in previous papers. Indeed, by comparing the effectiveness of different LC and habitat taxonomies in monitoring coastal wetlands emerged that, for long-term habitat monitoring and change detection, the coupling EUNIS and LCCS is highly recommended (Adamo et al. 2014(Adamo et al. , 2016Tomaselli et al. 2013Tomaselli et al. , 2016Tomaselli et al. , 2021Gavish et al. 2018). The output maps were validated by in-field campaigns carried out in 2020 and 2021. Information on vegetation composition and structure, as well as 1 3 Vol:. (1234567890) Fig. 1 Geographical location of the SCI "Zone Umide della Capitanata e Paludi presso il Golfo di Manfredonia" agricultural practices or land use, was gathered, geocoded by GPS and integrated into a GIS geo-database for an accurate and detailed definition of some types.

Changes
Magnitude of Changes (MCs) in class area occurred between 2010 and 2020 for each habitat class was calculated by using the following formula (Abbas 2013;Abbas et al. 2018): where i is the habitat class considered and, for the case under study, T1 and T2 correspond to 2010 and 2020, respectively. CA represents the Class Area recorded for each class.
To describe the conversion size of habitat types in different periods, the transition matrix approach was used (Tomaselli et al. 2021). The transition matrix indicates the amount of different habitat types that remain unchanged and change in the study period. Based on the transition matrix, the following parameters were calculated for each habitat type: (1) the percentage of the 2010 habitat area that experienced a change (losses); (2) the percentage of the 2020 habitat area that resulted from a change (gains). For example, the gains or losses of habitat type i stand for the other habitat types having changed into i or habitat type i having been converted into other habitat types, respectively.
In order to compare in a more detailed way the changes occurred in some specific areas of the study site, we chose to analyze the transformations within "Palude Frattarolo" (PF), "Saline di Margherita di Savoia" (SMS), and in the "Lago Salso" wetland (LS), because subject to rapid transformations (PF and LS), and in order to compare the ongoing dynamics in protected areas with different management type (Fig. 1). To analyze the effects of the hydrogeological process of subsidence (Caldara et al. 2013), we added also the coastal area between Cervaro and Carapelle rivers, called "Ippocampo" (I) (Fig. 2).

Landscape metrics
Basing on an analysis of the literature on monitoring programs in Mediterranean coastal ecosystems, among the most widely used and effective in MPS, LPI and MESH result to be good metrics for assessing the patchiness of landscapes, SHAPE and ED for evaluating landscape complexity, defining habitat network and functional features (Almeida et al. 2016). These LM were implemented using the LecoS-Land cover statistics plugin (https:// plugi ns. qgis. org/ plugi ns/ LecoS/) of the open-source Geographic Information System (GIS) software QGIS (https:// www. qgis. org/ it/ site/).
In coastal landscapes, the assessment (measurement) of the "elongatedness" of patches, especially for dune systems, is a crucial issue to detect the integrity of these environments (Botequilha-Leitão and Ahren 2002; Tomaselli et al. 2012). Here we introduced a new metric, extracted by using eCognition Developer 8.9, the Mean Length/Width (MLW) (https:// docs. ecogn ition. com/). For the evaluation of the adjacency between dune systems habitats, useful to evaluate in what measure the standard zonation has been altered, we introduced the "Mean Rel. Border to" (MRBT) https:// docs. ecogn ition. com/).

Pressures and driving factors
In order to describe the anthropogenic pressures affecting the study area in an organic scheme, we resorted to the DPSIR conceptual framework (Gabrielsen and Bosch 2003). Within this frame, "driving forces" (causes) generate "pressures" on the environment, modifying its "state" (physical, chemical, and biological), leading to "impacts" on ecosystems (structure and function) and eventually to "responses" (policy). To relativize the main categories, we referred to the Unified Classifications of Direct Threats and Conservation Measures Actions (IUCN-CMP 2012a) for "drivers of changes", and to IUCN-CMP Classification of Stresses (IUCN-CMP 2012b) and to Nagendra et al. (2014) for "broad impact (stress) category". "Broad impact category", "specific type of impact" and "short description of the impact" refer to changes observed in the period of observation; "direct threat" (proximate pressure) and "underlying factors" (drivers of change) have been identified based on in field observations, and/or through interviews. In this latter respect, representatives of local authorities or management bodies, as well as members of the local communities (e.g., farmers, stakeholders, etc.) were interviewed using a semistructured interview approach, covering issues related to: the implementation of conservation/ restoration strategies (e.g., LIFE projects); water management and water supply systems; cultivation systems and agricultural practices (included crop rotation and irrigation); fire frequency.

Results
Habitat and LC maps of the whole site Figure 3 shows the output EUNIS maps obtained in 2010 and in 2020. LCCS maps are in Online resource 2. The complete list of habitat types, in relation to LC classes, are reported in Online resource 3.
A high degree of landscape heterogeneity characterizes the site, due to both an effective natural diversity of the biotope and a variety of land uses. In the following sections, further considerations on landscape composition and on the degree of fragmentation of the site are provided.

Landscape composition
The overall landscape composition of the area in 2020 results as following. Cultivated areas (Croplands) are the predominant type in the whole site, covering the 39% of the SCI, with arable lands (I1.1-intensive unmixed crops and I1.2-mixed crops of market gardens and horticulture) making up the most part (38%). The class I1.1 is mainly distributed landwards while I1.2 extends along the sandy coast, the socalled "arenili". The next dominant landscape type, which covers about 37% of the whole area, is represented by natural and artificial water bodies (Table 1), with water from fresh to salty, and including the intricate system of canals and pools used to drain the cultivated areas. The coastal lagoons are highly present with the dominant class X02 (31%), principally represented by the salines of Margherita di Savoia.
As regards the natural vegetation, the most extensive and representative types are the helophytic communities (9.2%), in which the class C3.2 prevails with 6.5%, and halophytic shrubs and annual herbaceous communities of saline to hypersaline environments (7.4%) in which the most relevant class is A2.526 with 6.1% (Table 2).
Although the coastline of the site is over 30 km long, the classes of the sandy coast system (classes of group B are shown in Online resource 4), including both vegetated and not vegetated areas, cover only the 0.70% of the whole area.
This general outcome does not differ much from the 2010, but this is true only if considering coarse categories. If going into detailed description of single classes, numerous changes can be described.

Changes
Changes observed over the period 2010-2020 in the study area are of two types: (a) conversion from one class to another (inter-class changes); (b) modifications within a specific class (intra-class changes). In the case of class conversion, they may be: (a) real changes; (b) conversions due to a change in thematic resolutions (from broad to more detailed classes).

Inter-class changes (class conversion)
Habitat maps from 2010 and 2020 were analyzed to obtain the habitat Transition Matrix (TM). The analysis of the TM (Fig. 4) revealed an overall percentage of areal changes equal to 3.99% (about 564 ha).
The Magnitude of Change (MC) in class area (ha) along with gains and losses (%) occurred during 2010-2020 are reported in Table 3.
The first evident change is that two classes, included in the map 2010, were completely converted: A2.5 (Coastal saltmarshes and saline reedbeds) and J2.7 (Rural construction and demolition sites). The change of J2.7 is due to the correction of a misleading labelling, whereas, in the case of A2.5 it is partly due to a thematic redefinition of the mapping product and partly to a real change. In particular, the class A2.5 represents a very coarse habitat class that, in 2010, was used to include, in some cases, both the perennial and the annual halophilous vegetation of saltmarshes. In the 2020 map, an increasing thematic resolution has been applied and class A2.5 has been replaced with the more detailed classes A2.526 (Mediterranean Most of the observed conversions is borne by classes belonging to the Eunis high level categories A (marine habitats, directly or indirectly connected to the marine waters, included saltmarshes and constructed marine saline habitats), and B (coastal habitats, including coastal dunes and beaches) (Davies et al. 2004).
Group B-the whole group, corresponding to the sand coast system, shows deep transformations due to processes of coastal erosion, in some traits, and accretion in other traits, resulting in a general reshaping of the coastline. In particular: -The class B1.1 (Sand beach driftlines) shows a drastic transformation with gains and losses Analyzing the TM (Fig. 4) and Table 3, other significant changes are: Class E1.61 reduces its surface (MC = − 199 ha) with the conversion into arable lands (I1.1 31.55%; I1.2 4.03%).
As regards all other classes belonging to groups G, I, J and X, they do not indicate major transitions, only small oscillations, mainly between cultivated typologies or other land use types. LCCS allows describing the invasion of helophytes in the Tamarix community (A12/A1A4A11B3-A12B14-Open ((70-60)-40%) Medium to High Shrubs), using stratification classifiers (A12/ A1A4A11B3F2F4F7G4-A12B9F8G11-Open ((70-60)-40%) Medium to High Shrub land with closed medium to tall herbaceous vegetation).

Data subset results
The analyses of the four data subsets corresponding to "Palude Frattarolo" (PF), "Lago Salso" (LS), "Ippocampo" (I) and "Saline di Margherita di Savoia" (SMS) lead to the following results. The Eunis maps of the four areas are shown in Online resource 5. The MC, gains and losses in class area occurred during 2010-2020 are reported in Online resource 6.
From the observation of the TMs (Online resource 7), the total change percentages of habitat surfaces are: PF 31.9%, LS 19.0%, I 9.4% and SMS 1.6%. The principal and more interesting results are discussed hereafter.
The increase of C1.3 water bodies (MC = 61 ha) is due to two factors: the natural phenomenon of subsidence that have involved some areas of Capitanata during the last two decades (Caldara et al. 2013), and the interventions carried out with the Life+ project "Conservation Activities for Priority Avifauna in the Lago Salso Oasis". Specifically, the Life actions consisted in the opening of ponds in the reedbeds and the creation of a 90-ha basin, in the southern part of the Oasis, with variable water levels (Fig. 8). The significant expansion of Juncus subulatus vegetation ((MC = 89.7 ha) is strictly linked to the subsidence.
In the site I a general decrease of saltmarshes (A2.526, MC = − 25.3 ha, with losses of 53.5%) has been observed, due, on one hand, to the conversion in arable lands and (Fig. 9), and on the other hand, to the transformation in coastal lagoons (19.4%). This latter is due to the effects of the subsidence, a trend widely generalized in the site that is causing also an increment of coastal lagoons (X02 with gains of 30.7%; X03 with gains of 67.5%) (Fig. 10),  the conversion of the classes of sandy coast system (16.3%) and arable or fallow lands in saltmarshes (Fig. 11), and the transformation of fallow lands in Suaeda vera J.F. Gmel. vegetation (class A2.516 with 37.5%).
The site SMS presents a percentage of change smaller than the general one. In detail, the changes are due to thematic redefinition or occur at the border of the site, where there is the contact with croplands.
Changes through landscape metrics CA data relative to all habitat classes are summarized in Fig. 12, data about Eunis classes of group A and B in Fig. 13; data have been partially discussed in Landscape general composition.
MLW is a measure of elongatedness of a patch; in the case of the sand dune systems (habitat types B1.1, B1.2, B1.31and B1.4) it is a proxy of the integrity of the vegetation strips forming this environment, and indicative of their level of fragmentation. The most striking outcome is for class B1.31 (Fig. 15), shifting from 2.8 to 10.4; this is in accordance with CA that has an important increment. No particular change is recorded for class B1.1, while for class B1.2 the MLW decreases from 13.5 to 6.3, indicating, along with a significative reduction, a fragmentation process, mainly due to coastal erosion.
MRBT is a measure to quantify the length of the edges in contact between two adjacent classes and indicative for the standard zonation. In the case of class B1.31 (Fig. 16A), an adjacency with B1.1 is expected, but the value is low in 2010 and becomes much lower in 2020, while it increases the contact with B1.2. As regards class B1.4 (Fig. 16B), a contact with B1.31 is expected, but no adjacency is recorded, neither in 2010 nor in 2020, while it borders on class B1.2 for a large part. This output indicates a deep alteration of the standard zonation, even worst in 2020.

Impacts
In Table 4 the impacts observed on natural and seminatural classes during the period of observation, limited to those with highest impact (e.g., in terms of % of surface area affected), have been reported. Most of the observed impacts (14 out of 20) fall within the broad impact category "Land cover/habitat conversion" (1.1 Ecosystem Conversion); a minor portion has been classified as "Land cover/habitat modification" (1.2 Ecosystem Degradation and/or 2.3 Indirect Species Effect) and as "Habitat fragmentation and change in landscape connectivity" (1.3 Indirect Ecosystem Effects). Actually, this last category concerns a higher number of habitat types, but here we reported those that have been affected by fragmentation and change in landscape connectivity to such an extent that processes and functions may result to be altered. This is the case of the habitat types of the coastal dune systems, which resulted more or less all affected by this process along the entire length of the coastline. Among the "underlying factors", "Agriculture" (Agricultural practices intensification) seems to be the main driver of change, along with "Water Management/Use" (Changing water flow patterns from their natural range of variation; Abstraction of ground water) and "Other Ecosystem Modifications"

Discussions and conclusions
During the period of observation, a wide range of changes has been detected, in part due to an improving in thematic resolution, but in a large extent, attributable to anthropogenic pressures. Most of the observed conversions is borne by classes belonging to the Eunis high-level categories A and B, corresponding to saltmarshes and to coastal dune systems respectively, which are the main types of natural ecosystems in the site.
In the case of salt marshes, the main pressures turned out to be intensification and expansion of agricultural areas, changes in land management, and modifications of the hydrological regime. In this last case, the underlying factors lead back to agricultural activities (e.g., uptake of surface and groundwater), land use conversions, but also to direct interventions to the natural environments. Subsidence, that plays a key role in some areas, may be largely linked back to water caption, driven by agricultural intensification. Similar processes, consisting of temporal changes in vegetation components, and related to water management practices, have been also observed in other In particular, agriculture has been identified as the economic sector with the highest impact on wetlands, both directly (land reclamation) and indirectly (water management practices) (MedWet Secretariat 2016).
In the case of the sand dune system, this ecosystem is long-time severely compromised, subject to severe pressures by both land (intensive agriculture, touristic facilities) and sea (coastal erosion). Some dune habitats have been lost long ago (e.g., Ammophila arenaria communities-Eunis B1.31 white dunes, Juniperus macrocarpa communities-Eunis B1.631 Dune prickly juniper thickets), while others, especially those corresponding to the first zones of the standard zonation, result severely altered in distribution pattern, structure and functions. During the period of observation, a rapid dynamic of the coastline and of the related habitat types has been observed, with processes of erosion and accretion and often with rearrangement of the natural vegetation. Among the drivers producing such processes, both natural (the reduction of sediment flow from the Ofanto river), and  anthropogenic causes (the realization of infrastructures, such as the Margherita di Savoia harbor and the groynes along the coastline in the southern part of the study area) have been identified, especially in terms of accelerating erosion rates (Caldara et al. 2008). This is in accordance with what observed in many other coastal areas in the Mediterranean, where decrease in sedimentary contribution, sea level rise, and subsidence, are the main driving factors of coastal retreat (Rodríguez-Santalla and Navarro 2021).
The Mediterranean basin represents one of the most responsive region (primary hot spot) to climate change caused by human activities (Giorgi 2006). The principal threats that affect the Adriatic portion of the basin are increase in temperatures, decrease of precipitation during the spring and summer seasons, sea level rise. These factors may affect in different ways coastal plant communities. Nevertheless, in the case of our analysis and considering the 10-year timespan, the anthropogenic direct factors act at a faster rate in the determination of changes in land cover of coastal habitats.
As regards landscape composition and configuration, landscape metrics indicate a reduction in landscape complexity and a more even distribution and higher patchiness for some habitat types of group A. The outcomes for habitat types of group B indicate a deep alteration in the standard zonation of the whole stretch of sandy coast, and this alteration has become even more severe in 2020. MLW and MRBT turned out to be very effective metrics to assess the integrity of habitat types of group B (Coastal dune systems) and of their standard zonation. In the case of class B1.31, the significant increment of CA along with MLW indicates a positive trend with the reestablishment of this habitat along large stretches of coastline; the creation of some artificial barriers built at the edges between cultivated fields ("arenili") and the beach, has determined the accumulation of sand with following colonization by the typical vegetation of the shifting dunes (e.g., Thinopyrum junceum). It is worth noting that the re-establishment of this habitat is supported by human intervention, although the purpose was aimed at the protection of cultivated fields. MRBT may be very useful in evaluating the integrity of the standard zonation, as it gives a concrete measure of the contact between the different vegetation belts.
The separate analysis of the sub-set areas revealed the way in which different management strategies and practices may affect landscape composition and nature conservation. In particular, in the protected areas PF and LS there are the most significant total change percentages, much higher than the percentage of the whole study area. These are the areas where human interventions and changes in management, mostly involving agriculture practices and water flow patterns, have significantly affected the natural systems in the last decade; specifically, it is the case of water extraction in surrounding agricultural areas, and also of some management practices provided in restoration programs (e.g., creation of drainage channels that have negatively impacted on the hydric balance of the wetland area).
A system that has been undergoing such rapid transformations requires well-targeted policies and effective interventions. Coastal zones are very sensitive and dynamic systems, whose management planning requires in-depth knowledge of the natural dynamics as well as a full understanding of the social, economic, and political context (Damiani et al. 2002). With a view to identifying appropriate management and conservation measures, the Integrated Coastal Zone Management (ICZM) is a strategic reference point. ICZM is a widely accepted approach for sustainable management of the coastal environment and consists of the legal and institutional framework necessary to ensure that development and management plans for coastal zones are integrated with environmental goals (McKenna et al. 2008;Dronkers 2022).
From the observed trends for both drivers/pressures and impacts obtained from the DPSIR analysis, and according to the key principles of the ICZM, some alternative recommendations could be provided: 1. Creation of buffer zones (limited to sensitive areas, such as those in spatial contact to agricultural areas subject to intensive exploitation), may contribute to mitigate the effects of agriculture on water quality and regime; in order to prevent conflicts with stakeholders, regulated and controlled uses can be allowed in these areas, as well as in adjacent areas-buffer zones could also be considered as effective means to mitigate the effects of sea level rising and coastal erosion; 2. Higher control and management of surrounding activities; 3. Promotion of pro-environmental activities that could have a lower environmental impact on the overall system quality, such as eco sustainable activities, organic farming and, above all, proper management of water resources. In relation to the latter, crucial point, among the ICZM specific objectives, "Limitation of soil subsidence" ("Enhancing sustainability and ecosystem services") is included, with the actions provided: (a) regulations for groundwater extraction and drainage; (b) groundwater management and alternative water supply (incl. recharge of aquifers). Moreover, within the specific object "Land-use" ("Economic development"), the actions provided are: (a) ban on urban development in sensitive zones; (b) coastal zone water management plan and implementation (Clark 1997;Fujita et al. 2013;Dronkers 2022).
In order to achieve a comprehensive and effective management plan, this framework should be implemented at different (both local and national) administrative levels, and an effective coordination is needed among various organizations.