Assessment of State Transition Dynamics of Coastal Wetlands in Northern Venice Lagoon, Italy

: Coastal wetlands represent particularly valuable natural resources, characterized by the interaction between their geomorphological and biological components. Their adaptation to the changing conditions depends on the rate and extent of spatial and temporal processes and their response is still not fully understood. This work aims at detecting and improving the understanding of the transition dynamics on eco-geomorphological structures in a coastal wetland ecosystem. The approach could support sustainable habitat management improving the detection and optimizing the offer of Earth Observation (EO) products for coastal system monitoring. Such course of action will strengthen evidence-based policy making, surface biophysical data sovereignty and the Space Data downstream sector through remote sensing techniques thanks to the capability of investigating larger scale and short-to-long-term dynamics. The selected case study is the Lido basin (Venice Lagoon, Italy). Our methodology offers a support in the framework of nature-based solutions, allowing the identiﬁcation of ecosystem-level indicators of the surface biophysical properties inﬂuencing stability and evolution of intertidal ﬂats on which a conceptual model is implemented. Landsat satellite imagery is used to delineate the spatial and temporal variability of the main vegetation and sediment typologies in 1990–2011. Within this period, speciﬁc anthropic activities were carried out for morphological restoration and ﬂood protection interventions. Speciﬁcally, the lower saltmarsh shows its more fragmented part in the Baccan islet, a residual sandy spit in front of the Lido inlet. The area covered by Sarcocornia-Limonium , that triggers sediment deposition, has ﬂuctuated yearly, from a minimum coverage of 13% to a maximum of 50%. The second decade (2001–2009) is identiﬁed as the period with major changes of halophytic and Algae-Bioﬁlm cover typologies distribution. The power law and related thresholds, representing the patch size frequency distribution, is an indicator of the ecosystem state transition dynamics. The approach, based on multi-temporal and spatial EO analysis, is scalable elsewhere, from regional to local-to-global scale, considering the variability of climate data and anthropogenic activities. The present research also supports sustainable habitat management, improving the detection, and optimizing the offer of EO products for coastal system monitoring.


Introduction
Coastal wetlands represent complex ecosystems prone to continuous changes. These changes are especially standing out nowadays due to the existing feedback loop between between traditional and innovative remote sensing platforms (e.g., satellites and drones) allow the near real-time monitoring of natural and human-induced changes in tidal flows. The blended monitoring approach (traditional with innovative) allows expanding the model to a variety of complex systems, from marine to marine-coastal (subtidal belt) and transitional (intertidal belt) systems, each subject to different prevailing external forcing.
Complex systems, where small-scale interactions scale up to large-scale changes, are characterized by a range of small and large patches that determines the spatial structure of the system. The Power Law Distribution (PLD), a statistical indication of the scaledependent distributions, offers a unique tool to explore physical mechanisms and selforganization of sediments and vegetation spatial patterns [38][39][40]. Particularly, PLD of patch sizes can be interpreted as a sign of self-organization and it allows characterizing the spatial processes occurring at the boundary of the existing patches, determined by their spatial extension or contraction [27,30,41,42].
According to recent studies [40,43,44], it is fundamental to understand the non-linear dynamics observed in the power law tail, whereas the deviation in the patches' size distribution is crucial to forecast the system response to environmental changes. It evidences the tuning of morphological restoration (mudflat-saltmarsh anthropogenic reconstruction) to the subsequent complex feedbacks between the vegetation-sediment distribution and external drivers. Recent studies in lake systems, oceans, and forests show how the gradual change is in some cases interrupted by sudden changes that determine the transition to a new state, when a threshold value in environmental conditions is exceeded [45].
Temporal trends' characterization, as extent and size of patches' patterns, is further fundamental to discriminate the changes in both space and time directions. The Empirical Orthogonal Function (EOF) [26] is a valuable tool, largely used in climatology, oceanography, and Earth Observation (EO) disciplines to rank spatial patterns of variability, their time variation and the importance of each pattern on the basis of variance [46]. Particularly, it allows to isolate the major phases of changes through the EOFs' expansion coefficients along the time series and to assign a measure of the importance of each pattern.
The integration of field observations and quantitative tools, such as Earth Observation (EO) enable an accurate insight into the wide range of spatial scales and mechanisms. Satellite remote sensing allows accurate and repeatable data acquisitions over large areas [8,14,[47][48][49][50], therefore, it is a key data source, especially due to nowadays with easy access to open satellite data [51,52]. Given the fact that small perturbations cause the instability of the system and rapid changes in the landscape, EO is an interesting tool to study the scale variance and invariance of bio-geomorphological process [27,32,[52][53][54][55][56][57][58].
EO-based studies can highlight the conditions under which variations in patch size distributions are important determinants, e.g., for animals living in patchy intertidal landscapes like many migratory and non-migratory ornithic species [50,59]. Adaptation in the species' population size may occur responding to habitat loss (e.g., reclamation or flood) and to human recreational or productive activities. Seasonal hydrologic variability, as well as artificially induced water-level fluctuations, has been proven to threaten the reproductive success of some bird species in different climatic zones of the world [60]. Tidal flats are significantly impacted by rising sea levels and intensifying storms under climate change [61]. For example, given that most tidal flats in the Yellow Sea are bounded by rock walls and coastal development, ecosystem migration in response to rising sea levels cannot occur, which lead to inundation of the ecosystem, reducing the key process of regular tidal inundation, changes in distribution and declines in extent [62]. Therefore, the spatial and temporal trends of changes in intertidal areas within a lagoon can be fundamental for conservation and resilience studies [63].
This paper presents a comprehensive method developed through the analysis of multi-source open data applied to a shallow tidal basin of saltmarshes and mudflats in the Northern part of the Venice Lagoon (North-East of Italy), the Lido basin. The research questions addressed are built on the conceptual model of mudflats and their transition dynamics, and have the goal to improve the understanding of eco-geomorphological struc-Sustainability 2021, 13, 4102 4 of 24 ture changes using vegetation and biofilm cover classes. The spatial structure dynamics and the way spatial patterns manifest themselves over time (i.e., the PLD and the EOF analyses) are used to detect the influence of external drivers. In the study area, particular emphasis has been given to bio-stabilization processes linked to morphological restoration as a major thematic area. The spatial and temporal ecosystem-level indicators allow identifying which kind of processes influence the behavior of the spatial structures' response to environmental changes, leading to possible evidence of criticalities. This is the base for sustainable management policies driven by the adoption of naturebased solutions. Using tradeoffs between exploitation and conservation, especially for utilization of sediment, water and biological resources, is key to guaranteeing the maintenance of ecosystem services provision.

Study Area: The Lido Basin in the Venice Lagoon
The study area is in one of the most valuable coastal transitional ecosystems in the Mediterranean Sea: The Venice Lagoon. The Lagoon is suffering of ecosystem degradation. It is located along the North-Western coastline of the Adriatic Sea, bounded by the Sile River to the North and the Brenta River to the South. It is connected to the sea by three inlets of 500 m to 900 m wide and to 15 m deep. Venice Lagoon is divided into three basins, named after the inlets: Lido in the North (276 km 2 ), Malamocco in the middle (112 km 2 ) and Chioggia in the South (111 km 2 ).
The Lagoon is the largest wetland in Italy, with a surface of 550 km 2 (Table 1). Venice Lagoon is a microtidal environment in which the neap tide is 30 cm, the mean tide is 55 cm and the spring tide is 110 cm [19]. Tidal waves are semidiurnal and asymmetric in form. Two high tides (HWL) and low tides (LWL) occur every day. HWL are increased by Scirocco wind, low atmospheric pressure and river discharge. Acqua Alta events can increase the HWL (~1 m) during spring tides and southerly winds. The maximum HWL reached up to 192 cm on 4 November 1966. Tidal wave is slowed down when entering the Lagoon and propagates within the Lagoon different delays. The tide height is also lowered within the Lagoon due to the friction exerted by the bottom [64]. The complexity of its landscape results from the interaction between natural processes and human interventions. Human activities have modified the environment of the Lagoon in several ways during the last five centuries. A drastic change has occurred in the hydrological regime of the rivers flowing into the Lagoon. This modern morphology has been created by a variety of human activities, also during the Anthropocene, and is classified into six categories (Table 1).
During the time interval considered in the present study (two decades between 1990-2011), several interventions were carried out within the Venice Lagoon as part of a Morphological Master Plan of the Consorzio Venezia Nuova (CVN). Since the 1990s, several dredging and saltmarsh restoration were carried out. In the Northern part of the Venice Lagoon, saltmarsh restoration was carried out in the period between 2000 and 2002. Then, since 2007, the realization of hard infrastructures at the three inlets for the Lagoon took place. The study area has thus been always subject to natural siltation mitigated by different types of anthropogenic interventions such as dredging, protection against high water, river diversion, and morphological restoration [66]. Recently, the Experimental Electromechanical Module (MOSE) project was installed across the three inlets that link the Lagoon to the Adriatic Sea, to protect the city of Venice from flooding [67]. The mobile barriers consist of 79 elements placed perpendicularly to the main three inlets and represent a significant anthropogenic structure that influence the hydrodynamics, sediment transport, and the habitat distribution of the Lagoon [67,68].
The MOSE project has been designed to protect the city of Venice against the Acqua Alta phenomenon, but during its construction, the winter seasons in 2008-2010 were characterized by numerous episodes of high tide and heavy rainfall. 2008-2010 is characterized as anomalous in comparison to the last 30 years and this change in the frequency of inundation led to eco-morphological changes.
The selected case study (16 km 2 ) is within the Lido basin, represented by temporarily exposed mudbanks and barene and velme (local name for the saltmarshes and mudflats) ( Figure 1). The basin (276 km 2 ; 50% of Venice Lagoon surface) comprises the islet of Baccan (a residual sandy spit) and part of the San Felice marshy wetlands, incised by meandering tidal networks. Venice Lagoon shows unique ecological characteristics at a national level and in the context of the Mediterranean basin. In fact, it is included in the list of wetlands of international importance [69] and the European Community, with the financial instrument for the environment (LIFE), is granting a wide variety of funds for nature conservation and restoration programs.
Venice Lagoon is regulated by special laws, thereby defining a benchmark for pioneering legislation of coastal sustainable development to balance anthropogenic impact and natural capital conservation. The Lido basin area was selected because of its ecological and biological value in terms of habitat conservation. The Northern part of the Venice Lagoon holds several listed habitats of the Habitats and the Birds Directives (1992/43/EEC and 2009/147/EEC). In addition, the Directive 2000/60/EC of the European Parliament and of the Council, establishes a Framework for Community Action in the Field of Water Policy to protect aquatic environments by making it compulsory to achieve a good environmental status for the coastal waters. In Italy, a series of bioindicator species have been selected to establish the ecological status of marine systems.
Mainly halophytic species colonized the basin, e.g., Spartina maritima, Salicornia veneta, Limonium narbonense, Sarcocornia fruticosa, and Juncus spp. [8,47]. The continuous changes in the sediment stability of the Northern Adriatic wetland system have led to fluctuations of the population of migratory ornithic species, such as the Little Tern (Sterna albifrons Pallas). Italy hosts the largest population in Europe of the Little Tern and about 85% of the population lives along the coast of the Northern Adriatic. Specifically, a post-breeding colony of thousands of individuals have historically sited in the Venice Lagoon between mid-June and early September, in correspondence to the existence of an intertidal sandy bar that is not submerged during the average tide [19,50].
The glossary used to describe the coastal tidal system metrics in the following sections is graphically displayed in Figure 2.

Climate Settings
Climate contributes to the unusualness of the landscape and biodiversity of the Venice Lagoon. The Northern Adriatic is a slightly more humid sector than usual in the Mediterranean Sea as a transition between the Mediterranean and the Atlantic climate conditions. Data obtained for the Cavallino (Treporti) station (ISPRA service for the Venice Lagoon) show that precipitation peaks characterize the study area in spring and autumn. Higher values are obtained in September and October (98-104 mm) and lower from January to March (42-49 mm). The annual mean temperature is 14 • C, with a relatively hot and humid summer (Table 2) [70]. Table 2. The monthly cumulative precipitation (P) and the monthly temperature (T) in Venice Lagoon. Data obtained for the Cavallino (Treporti) station for the period of 1990-2011. Source: Italian Institute for Environmental Protection and Research (ISPRA), freely available online [70].

Monthly Mean
Annual

Characterization of Dominant Processes
Tides and winds (Bora and Scirocco, from the North-and South-East, respectively) are the main driving forces that affect the hydrological and sedimentary dynamics of the Lagoon [71,72]. Land subsidence also contributes to the elevation loss with respect to the mean sea level. In fact, most of the Northern Adriatic coastland lays below the mean sea level. Ground displacements of this coastland are caused by the individual and/or combined action of natural and anthropogenic factors (i.e., eustacy and groundwater pumping, respectively). In the Northern sector of the Venice Lagoon, this process continues with sinking rate on saltmarshes and intertidal environments ranging from 5 to 10 mm/year [73,74]. The short-term sedimentation and surface elevation change of wetlands in the Lagoon are due to changing water levels which have been measured at several stations with varying sediment input/availability and wave energy [75].
Mainly medium-fine sediments compose the Lagoon: silty sand, sandy silt, and silt. The grain size distribution reflects the environmental hydrodynamic energy. The coarsest sediments (sands) are distributed along the inlets and main tidal channels [19]. Both the tidal regime and the wave action reach their maximum in these areas, where sediment transportation is based on traction, rolling, and saltation. Instead, the finer sediments (silts and clays) concentrate along the inland areas, where a low flow energy deposits the fine matter that remained suspended in the water [76].

Characterization of Tidal Regime
The environment of the Lagoon is characterized by the tidal forcing with a first harmonic at 12 h with an average tidal range of about 60 cm, and about 1 m excursion in spring. The tide propagates along the deep narrow channels onto the tidal flats and tidal marshes [77]. The frequency distribution of high-extreme tide levels over 110 cm increased in the last decade, especially in 2008-2010. In the same period, Venice Lagoon was affected by the highest mean sea level, which reached 40.5 cm in 2010 ( Figure 3). This value lays into the highest observed in the last 120 years [70].

Materials
The developed methodology required different kind of inputs, namely EO data, field spectral data, and the annual tide level and mean sea level for the period of 1990-2011 ( Figure 4). From the United States Geological Survey Earth Explorer, 18 cloud-free Landsat images (level 1T, Top-Of-Atmosphere radiance) covering the period 1991-2011 were acquired. The images are acquired by the Landsat-4, Landsat-5, and Landsat-7 satellite platforms (Table A1)   The selected images for the case study are those preferably captured in July due to the higher vegetation vigorousness. With a pixel spatial resolution of 30 × 30 m, the entire surface selected as study area is about 15.804 km 2 . The processed data cube has an extent of 17,560 pixels in terms of surface and it is made by a stack of 108 layers of information (18 dates × 6 spectral bands each) for a total of 1,896,480 cells. Values of the tide level, recorded at the EO images acquisition time (Table A1), are used to represent the transition between the permanently emerged and submerged environment to mask the extension of the analysis. This allows avoiding detecting spatiotemporal variations that rely on tidal fluctuations, thereby reducing uncertainties.
The field spectral library has been collected in June 2013 to validate the spectral endmembers selected on the EO imagery and to characterize the main sediment, biofilm and vegetation cover and the mixtures of them (Figure 1).

Methodology
The implementation of tools for monitoring for sustainable development need to detect variables in the spatial and temporal dimensions. The methodology proposed to identify ecosystem-level indicators for state transition dynamics in the case study consists of two main processes: (i) image field-based analysis (Linear Spectral Mixture Analysis, LSMA), based on cover typologies classification and validation using field spectral sampling [53], and (ii) a spatiotemporal analysis applying the EOF and the PLD [27] (Figure 5).

Image Field-Based Analysis
The image field-based analysis includes a sub-pixel processing technique, i.e., the Linear Spectral Mixture Analysis (LSMA), to obtain the spatial distribution of vegetation typologies and different wet and dry sediment transitions. The selection of the endmembers, one of the main phases of LSMA process to achieve cover fractional abundance maps at a pixel level (30 × 30 m in the case study with the Landsat images), is validated using the spectral library obtained by field radiometry.
Field campaigns are performed for better understanding the sediments' variability, considering the influence of the biofilms covering the tidal bare substrates, and to collect the spectral library used to classify the vegetation typologies and validate the endmember collection [50,53].
The field spectral library acquired in June 2013 identifies the main vegetation and sediment typologies, e.g., Sarcocornia, Limonium, Juncus, and Algae-Biofilm. First, the field spectral library is resampled to the Landsat satellite spectral resolution. Then, it is compared with the spectral signatures of the vegetation and sediment endmembers previously selected. Subsequently, the field spectral signatures are used to simulate a reflectance spectrum (ρ(λ) linear ) as similar as possible to the reference spectral profile of the selected endmembers, according to Equation (1).
This equation was based on finding a linear relation between the reflectance values of the field spectral signature of each endmember (ρ em n ) and the equivalent endmember percentage ( f em n ). The difference (∆ρ(λ)) between the simulated spectral signature (ρ(λ) 0 ) and the reference one (ρ(λ)linear) must be as minimum as possible (Equation (2)). To achieve this goal, the Microsoft Excel Solver was used. This tool allows finding the values of certain cells in a spreadsheet that optimize (maximize or minimize) a certain objective. Thus, the optimization model comprises three components: 1.

2.
The changing cells: the cells of the spreadsheet that can be adjusted to optimize the target cell. In this case, the endmember fraction ( f em n ). 3.
The constraints: the restrictions placed on the changing cells.

Spatiotemporal Analysis
The main aim of analyzing the spatiotemporal trends of this coastal wetland is identifying the tipping points that are likely inducible to a critical state, as well as the main climatic, hydrodynamic, and morphological variables that influence and intensify this behavior.
Then, the permanently emerged areas, characterized by the presence of the four cover typologies, were assayed by means of its temporal, i.e., the EOF, and spatial, i.e., the PLD, distribution in order to characterize the variation trends over time.

Temporal Trends: The Empirical Orthogonal Function
To characterize the temporal trends in vegetation and sediment patterns, time series of the fractions cover (obtained from LSMA) are analyzed to discriminate the magnitude of the interannual changes. To this end, the EOF [13,26,55,[78][79][80] is applied. This method decomposes the time series dataset in terms of orthogonal functions, allowing to determine the spatial patterns of change and their variation over time. In addition, it provides a quantitative measure of the contribution of each of these patterns to the overall multitemporal change.
The outputs are a dimensionless map of change in the study area and the expansion coefficients, which represent the development of the phenomena during the study period. The time trend of an EOF shows how the vegetation and sediment fractions cover oscillate in time. The analysis of the expansion coefficients leads to the identification of peaks, both positive and negative, in the pattern change over the time series.

Spatial Trends: The Power Law Distribution
To characterize the spatial patterns of vegetation and sediment, the fractions' cover (obtained from LSMA) is analyzed in terms of patch size frequency distribution. The aim is to test whether the frequency distribution could be fit by a power law (Equation (3)) to detect the presence of scale-invariant patterns [81].
where P is the probability, x is the patch size and α the scaling exponent of the distribution [26][27][28]82]. The patch size is here forth considered as a non-linearity threshold that might describe changes in the main climatic, hydrodynamic, and morphological variables when a deviation in the power law tail occurs. Thus, the latter indicates the shift from a scale-invariant, self-organized, vegetation pattern to a scale-variant one. The mathematical simplicity of this statistical distribution allows describing complex natural phenomena by defining just one parameter, the scaling exponent, which might represent the presence and the level of transformation/disturbance experienced by the eco-morphological system.

Endmembers' Validation: The Field Spectral Library
Six endmembers have been distinguished through their spectral signature looking at the axes of the Landsat multi-temporal scatterplot, namely sediment, shallow water, two types of vegetation cover (with low and high reflectance), Algae-Biofilm, and water. The spectral signatures of the endmembers are validated using the field spectral library collected in 2013 ( Figure 6). Vegetation with low reflectance corresponds to the upper saltmarsh and, more specifically, to the transition between Sarcocornia-Limonium and Juncus. On the other hand, vegetation with high reflectance values corresponds to the lower saltmarsh, specifically to Sarcocornia-Limonium and Sarcocornia (Figure 7). The identification of the different vegetation typologies (i.e., Sarcocornia, Limonium, Juncus, and Algae) is possible through the correlation between the pixel spectral signature and the field data. Above 2500 spectral bands were resampled according to the satellite sensor spectral features (i.e., the spectral bands of the sensors on board Landsat-4, Landsat-5, and Landsat-7). This procedure assures the representativeness of the endmembers in the study area. The endmembers' validation provided well-fitted correlation for both Juncus (R 2 > 0.90) and Sarcocornia-Limonium (R 2 > 0.80) (Figure 7).

Classification of Spatial Patterns
Once the endmembers are defined and validated through the field spectral library, the LSMA allows for detecting their fractional cover abundance at a pixel level (30 × 30 m). Specifically, the fraction extraction and its yearly distribution in the period covered by the Landsat images (1990 to 2011) has been analyzed for Algae-Biofilm, dry sand, Sarcocornia-Limonium, and Juncus mixture due to their representativeness spatiotemporal variation in the study area (

Rate of Variation in the Spatial Patterns' Distribution
The annual rate of variation in the spatial patterns' distribution has been quantified for all four cover typologies and compared with the mean sea level. The analysis assesses the amount of change using 1990 as reference year rather than the absolute quantity. A significant reduction of the vegetation distribution is detected after 1994 (Figure 9). This fact could be related to the increased mean sea level (changing from 22.7 cm in 1994 to 28.2 cm in 1995, and with the highest peak in the last decade of 40.5 cm in 2010) (Figure 9), and especially the increased frequency of tide levels over 110 cm (going from 125 in 1994 to 799 in 1995, with the highest peak of 3657 also in 2010) (Figure 3).

Spatiotemporal Analysis
Based on the vegetation cover typologies previously identified through the satellite image classification, the spatiotemporal trends are analyzed within the Lido basin in the period of 1990-2011. The aim is clearly identifying the areas that are more likely inducible to a critical state. This way, the natural (climatic, hydrodynamic, or morphological) as well as the human-induced variables that influence and intensify this changing behavior can also be identified over the studied period.  These above peaks are frequently natural-(climatic) or human-induced change events. Specially, the maps of EOF expansion coefficient 2 (EOF 2) allow to indicate that changes in the upper saltmarsh and sediment stability are more significant in the Central and Southern part than in the Northern part ( Figure 11). EOF 1 shows the most obvious variations ( Figure 10). This is the time course of the most typical changes. Instead, EOF 2 allowed detecting the time course of less typical changes either natural or human-induced ( Figure 11).
The heterogeneous nature of the investigated site suggests that EOF changes are related to both natural and anthropogenic factors. So, considering the artificial anthropogenic restoration of the Lagoon morphology, the EO data-based pattern analysis shows the spatial and temporal variability, highlighting changes of state over time that have to be interpreted using the PLD.  The heterogeneous nature of the investigated site suggests that EOF changes are related to both natural and anthropogenic factors. So, considering the artificial anthropogenic restoration of the Lagoon morphology, the EO data-based pattern analysis shows

Spatial Trends: The Power Law Distribution (PLD)
Changes from power laws identify stochastic conditions. The PLD allows for analyzing the interaction among the vegetation and sediment patch size. A frequent deviation in the tail of the distributions is related to the larger upper saltmarsh's patches ( Figure 12). Disturbance-recovery processes result in predictable patterns over a range of spatial scale. Despite the patterns' irregularity, the graph shows that patches of all sizes are characterized by a predictable PLD. Figure 12 shows the PLD results for the years of higher interest in terms of linearity deviations, namely 1992namely , 2002namely , 2010namely , and 2011namely . Firstly, 1992 shows an important increase in the frequency of tide levels over 110 cm (Figure 3). The PLD of the lower saltmarsh shows a linear trend, related to more resilient and adaptive species to tide variations (pulsating, seasonal and interannual events). Instead, the upper saltmarsh PLD follows a non-linear trend of the largest patches. This fact can be related to climatic conditions, i.e., frequent tide levels exceeding 110 cm. As for Algae-Biofilm, significant variations occur. This element of the ecosystem increases its spatial expansion as a significant part of the surface subjected to the variation of the tide increases, thereby undergoing a colonization of the diatoms, which begins to build the mudflat and contributes to vertical accretion. Additionally, in 1992, there is a higher frequency of the smaller patches than in previous years. The PLD of the sediment shows a linear trend with little significant variations but the inexistence of large patches.
The period until 2002 is characterized by a constant increasing trend of the mean sea level and a variable frequency of tide levels over 110 cm (Figure 3). The PLD of the lower saltmarsh perfectly follows a linear trend. Instead, the upper saltmarsh shows a significant deviation in its tail for larger patches. The year 2010 is characterized by an increase in the trend of sea level measurements, as well as a significant increase in the frequency of tide levels over 110 cm (over 2500). The winter periods 2009-2010 and 2008-2009 were characterized by numerous episodes of high tide and heavy rainfall and are anomalous in comparison to the last 30 years. The lower saltmarsh shows its more fragmented part in the Baccan islet. There is no detour. The trend of the PLD in 2010 is like that of the upper saltmarsh in which the frequency size of the patches is linear but some intermediate size classes are missing. Previous years are characterized by a more pronounced variation of the size frequencies. Instead, the period of 2009-2010 undergoes a more gradual variation. As for the Algae-Biofilm, the size frequency distribution follows a linear trend with a deviation towards larger patches. The exposed part is greater, probably due to the decreasing tide level in 2009, and the larger patches are evident. Finally, the sediment PLD shows less frequency of larger patches. This fact is related to the coverage of biofilm. The graph shows the lack of intermediate classes and never reaches large patches.
2011 is the last studied year of the time series. It is significant because of the sea and tide levels. The lower saltmarsh shows linearity on its PLD trend and a higher frequency of larger sizes. The distribution of patch sizes is gradual and well distributed with a greater density of the intermediate classes. As for the upper saltmarsh, the PLD shows a high deviation for larger patches. The Algae-Biofilm has a similar trend to the lower saltmarsh.
As for the sediment, the PLD of 2004 is clearly related to the changes displayed by the EOF (Figure 10, year 2004). The sediment PLD greatly increases the slope and the absence of larger patches highlights an unstable phase (Figure 12, bottom right).

Discussion
The analysis and implemented conceptual models allow evaluating spatial patterns (i.e., small to large patches) and the time range within these patterns fluctuate before a change of state. The analysis of threshold values makes it also possible to understand the stability of these systems with respect to disturbances and risks associated with climate change or anthropogenic alterations.
Open satellite data analysis, namely the LSMA, the EOF, and the PLD have been used to spatiotemporally characterize the saltmarsh in the Lido basin in the Venice Lagoon. The LSMA allows obtaining a preliminary classification of the entire study area with the extrapolation of vegetation, sediment, and biofilm cover typologies, while the EOF characterizes vegetation and sediment interannual trends.
The observed peaks in the EOF maps and the deviation in the tails of the PLD, specifically in the period of 2001-2004 and 2006-2009, allowed detecting the main changes in the system.
Specifically, in the EOF results ( Figures 10 and 11), it is possible to see the effect of artificial mudflats reconstruction, protected by wooden piles along their perimeter (Figure 11c). Jetties, confined disposal facilities, dredging operations, and mobile barriers on the sea-bottom were carried out together with morphological restoration of intertidal areas in the mentioned years.
The PLD deviation points out changes in environmental conditions at saltmarsh edge ( Figure 12). The deviation in the tail of the distributions related to the largest patches (the more resistant in terms of stability) is observed in the periods of major variance in tidal levels' frequency and this represents a critical state change influenced by hydrodynamics of the area.
The upper saltmarsh shows linearity along the temporal series, even with the inexistence of the major patch classes, i.e., there is no change in the slope but a minimal deviation of the tail. As for the Algae-Biofilm, the PLD of 2002 shows the deviation of the tail towards the larger patches. Conversely, in the same and in the following years, the sediment shows a linear distribution, never reaching the largest patches. This fact could explain an extension of the stable mudflat surface. In 2004, the sediment shows an increase of the intermediate classes whereas almost an absence of the larger ones and a double slope in the trend line. Small patches are gradually progressing until the intermediate sizes. The decrease in the tide level can be associated to the larger percentage of sediment presence (18%).
The only strongly aggregated point is found in the Baccan area that is located in front of the Lido inlet and thereby more sensitive to the changes in the hydrodynamic regime. The results suggest that a feedback mechanism, mediated by microphytobenthos' influence, causes sediment bio-stabilization. This may explain the exceptional preservation of the intertidal mudflats in the Northern part of the study area [19,72,85]. Such a mechanism [18,19] can be detected by the EO methods for the three steps as follows ( Figure 13):

•
Step A: Subtidal areas, which are colonized by seagrass, are better protected against erosion compared to unvegetated areas. The ability of vegetation to reduce the applied bed shear stress enhances deposition of suspended sediment.

•
Step B: Vertical accretion can occur over time; bottom sediments are elevated from a subtidal to an intertidal environment, where time of exposure is higher.

•
Step C: Micro-organisms (diatoms and cyanobacteria), which are more abundant in intertidal environments, start colonizing the substrata producing Extracellular Polymeric Substances (EPS), particularly during summer and intertidal exposure, and increasing the erosion threshold and sediment stability (i.e., resistance against sediment waves and tidal currents). This feedback mechanism can be correlated to human interventions at the main three inlets of Venice Lagoon, where the construction of jetties was necessary for the placement of MOSE's barriers on the seafloor. In the last 30 years, morphological structures were realized with the reuse of about 19.5 × 106 m 3 of sediments from dredging channels. At the beginning-early 1990s, re-implantation of sea grass (Zoostera marina on sandy subtidal areas and Zoostera noolti) was tested because vegetation cover favors siltation and morphological restoration of shoals and intertidal flats. Then, dredging of sediments from navigation channels and their use to increase the elevation of velme (mudflat) and barene (saltmarshes) was the most common approach since 2000.
It does not surprise that in the period of 1989-2015, the presence of ornithic species has been targeted to understand how the artificial barren islands were populated. For example, for the little tern, even considering the natural fluctuating behavior of the species in the Lagoon during the nesting periods, a preference for the artificial barrens has been recorded with a complete abandonment of the natural barrens [86].
Considering the feedback mechanism proposed by [19], which explains how microorganisms contribute to the seasonal stabilization of sediments in the Northern intertidal mudflats of the Venice Lagoon, the change of state highlighted in the PLD results could imply a passage from one ecosystem regime to another, and therefore between basins of attraction-or stable states. The environment can take on different states driven by external changes that lead to the transition from one to another. The system opposes these changes until it reaches a threshold value, the deviation from linearity, specifically in 1992, 2002, 2010, and 2011, where internal relations begin to break. It is in fact maintained by internal relations between species, but at the same time, these dynamics can interact with forcing that act on larger scales. Therefore, the saltmarsh vegetation patterns, as well as their irregular spatial structure and its self-organization, could be used as an indicator of the ecosystem stability, either inducing steady states or pushing the ecosystem towards transition state dynamics [27,51,72,85]. The spatiotemporal trends can show the modulation of the system response to external drivers over time and support the estimation of the system resistance or resilience stability to disturbances. There is evidence of different inclinations of the PLD as well as different degrees of deviations over vegetation typologies and sediment over time. Detecting and understanding these trends is important to contribute to the description and forecasting of wetland ecosystems.
It could also be considered how the method can be applied to the MOSE management for the bio-stabilization processes during prolonged periods of still-water condition within the Lagoon. Will the closure of the mobile barriers influence the sedimentation and sediment stabilization processes and the morphological changes?
Effects of global change on development and sustainability of Venice Lagoon management is taken in great consideration. In 2012, CVN investigated the possibility to use the MOSE not only to protect the city of Venice from Acqua Alta, but also to partially open the floodgates in order to exploit nature-based energy (tidal cycles) to force or slow down flows within the Lagoon by using and generating different water level between the Northern and Southern part.
Changes in tidal height and frequency or the reclamation of lands over the intertidal zones can create unpredictable environments where there may be changes in biogeophysical variables. Studies regarding the adaptation of water level changes and the response of ecosystem structure (i.e., the biological communities) and functions (i.e., the state transition from bare tidal flat to saltmarshes) are widely developed, i.e., the Eastern Asian coastal ecoregion [61]. Key threatening processes often include damming and sea level rise or modification of the large systems inducing degradation, fragmentation and decline in the areal extent of tidal flats, which can lead to increased risk for the biodiversity [87], making available interesting analogies with the uniqueness of the MOSE effects in the Venice Lagoon unless the spatial scales and the climate region.
With this study, we are aiming at observing and discussing the bridge between the short-term (decadal) and long-term (3 decades) ecosystem-level transitions driven by the summed effects of sea level and human-driven changes of the Lagoon hydraulics. It is a contribution to studies on the adaptive capacities of biological communities and it contributes to better understanding of ecosystem state transitions.
Findings of the present paper can support coastal wetlands restoration programs that aim to increase the resilience of systems considering different scenarios of sea level rise and the natural capital conservation (i.e., the value of being a Ramsar and a Natura 2000 site). The system of adaptive management enables adjustments to be made at a site level but also touching the regional level, particularly in response to the impacts of naturalor human-driven events affecting water-sediment-vegetation patterns. It is important to highlight that adaptation initiatives need to be analyzed on both short-and long-term and both wide to fine scale.
Sustainable management and conservation agendas commonly include the requirements for mapping and monitoring environmental patterns and processes because, without monitoring and data collection, stakeholders are unable to implement conservation and short-and long-term protection actions, which are the fundamentals for planning anthropogenic activities, risk recovery, and prevention actions [88].

Conclusions
The present manuscript presents an approach, mainly based on satellite imagery acquired from Landsat, that allows identifying ecosystem-level indicators, and assessing the state transition dynamics in the Lido basin in the North of Venice Lagoon (North-East of Italy) in the period of 1990-2011. The methodology is based on multi-temporal and spatial EO-based analysis (the EOF and the PLD), thereby allowing its application elsewhere. The results of the spatial pattern analysis (PLD) manifest the influence of external drivers, mainly derived by anthropogenic activities. PLD's thresholds prove to be an effective method to detect changes in the size frequency distribution of different cover classes (upper and lower saltmarsh, Algae-Biofilm and sediment). Finally, the work has derived a conceptual model on bio-stabilization processes linked to morphological restoration.
The studied wetland can be considered an example of what is happening in many other wetlands due to sea level rise effects or hydrodynamic regime due to management. The EO-based analysis, coupled with climate data, measures of changes in erosion patterns, and water quality, can make different areas comparable. Moreover, biologically, short-and long-term repetition of environmental parameters' change have great impact on the aquatic species' adaptability. Therefore, the current approach can be further developed monitoring the response and resilience of environment to climate changes and/or human impacts on fishes and other species, the sediment transport and water quality (i.e., Anthropocene evidence). This would allow continuing the improvement of the understanding of the ecological response, scaling up from ecosystem-level to landscape-level indicators.
Considering the development of the current use of Copernicus data and products, and the European Commission's six high-priority candidate missions to address EU's policies and gaps in Copernicus users' needs, the presented operational method, tested on Landsat data, could be further expanded. Specifically, considering within the six upcoming missions the analysis of data from CHIME (Copernicus Hyperspectral Imaging Mission), the method applied in this study will optimize, as a monitoring tool, the set of results already analyzed in this context. Author Contributions: Conceptualization, A.T., E.V. and S.C.; data curation, E.V. and L.P.; methodology, A.T. and E.V.; Formal analysis and validation, E.V., L.P. and M.R.; All authors have contributed to the investigation, visualization, writing, review and editing; supervision, E.V and S.C.; project management and funding acquisition, A.T. All authors have read and agreed to the published version of the manuscript. Acknowledgments: Special thanks to the ISPRA groups of Chioggia, Venezia and Ozzano for sharing their knowledge, expertise, and for the fieldwork support. Thanks to Giuseppina Persichillo for the dataset collection and preliminary processing of EO data.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Satellite Data Used
The following table comprises the satellite data used, acquired by the sensors on-board Landsat satellite platforms, for the studied period of 1990-2011, to analyze the studied area of Lido basin in the Venice Lagoon.