High-Resolution Monitoring of Tidal Systems Using UAV: A Case Study on Poplar Island, MD (USA)

: Tidal processes regulating sediment accretion rates and vegetated platform erosion in tidal systems strongly affect salt marsh evolution. A balance between erosion and deposition in a restored salt marsh is crucial for analyzing restoration strategies to be adopted within a natural context. Marsh morphology is also coupled with tidal mudﬂats and channel networks and this makes micro-tidal systems crucial for a detailed assessment of restoration interventions. Here, we present a methodological approach for monitoring channel morphodynamics and vegetation variations over a time frame of six years in a low tidal energy salt marsh of the Paul S. Sarbanes Ecosystem Restoration Project at Poplar Island (Maryland, USA). The project is a restoration site where sediment dredged from the shipping channels in the upper Chesapeake Bay is used to restore a tidal marsh habitat in mid-Chesapeake Bay. Aerial surveys with an Unmanned Aerial Vehicle (UAV) have been performed for the high-resolution mapping of a small tidal system. Flight missions were planned to obtain a Ground Sample Distance (GSD) of 2 cm. Structure-from-Motion (SfM) and Multi-View-Stereo (MVS) algorithms have been used to reconstruct the 3D geometry of the site. The mapping of channel morphology and an elevation assessment on the mudﬂat were performed using orthomosaics, Digital Terrain Models (DTMs) and GNSS survey. The results highlight that the workﬂow adopted in this pilot work is suitable to assess the geomorphological evolution over time in a micro-tidal system. However, issues were encountered for salt marsh due to the presence of dense vegetation. The UAV-based photogrammetry approach with GNSS RTK ground surveys can hence be replicated in similar sites all over the world to evaluate restoration interventions and to develop new strategies for a better management of existing shorelines.


Introduction
The impact of human activities combined with meteorological events [1,2] such as hurricanes and storm surges [3] or sea level rise (SLR) [4][5][6], affects the sediment fluxes of coastal regions in a significant manner. The presence of vegetation along the shoreline [7,8] contributes in mitigating erosion of the most impacted areas. In fact, such erosion is mainly caused by a high shear stress on the bottom of un-vegetated tidal creeks, often associated with extreme events. The friction of dense vegetation reduces erosion, particularly when the tides are elevated, and the drag created by 2 m-tall marsh plants can be critical in slowing currents [9]. However, the factors mentioned above also impact on the capability of vegetated salt marshes to survive climate change [10][11][12]. Due to the important role Several approaches have been proposed in the literature [28] to evaluate and analyze shoreline changes. The use of satellite images [29] is suitable for large areas and advanced techniques such as a cluster-based segmentation may be applied [30]. However, the accurate detection of small systems requires high resolutions with a Ground Sample Distance (GSD) of a few centimeters. Digital Elevation Models (DEMs) represent primary tools in order to both detecting and characterizing morphology of such systems [31,32]. Different techniques exist for collecting data for DEM generation: (i) direct ground point surveys, (ii) Terrestrial Laser Scanning (TLS), (iii) aerial Light Detection And Ranging (LiDAR) systems and (iv) aerial photogrammetry from Unmanned Aerial Vehicles (UAVs) [33,34]. Direct methods are very time-consuming: (i) requires the survey of each point of interest, whereas (ii) implies post-processing operations to aligning and registering point clouds with the additional task of a data filtering. Aerial LiDAR is similar to (ii), essentially with the difference of a non-terrestrial data acquisition. UAV photogrammetry combines the flexibility of drones in carrying almost every kind of device, in particular cameras, with the Structure-from-Motion (SfM) and Multi-View-Stereo (MVS) approach that allows one to reconstruct the 3D geometry of objects. The density of photogrammetric dense point clouds is comparable to the one of TLS and LiDAR [35][36][37] and the high spatial resolution is also matched by high model accuracy [38]. For such reasons, UAV photogrammetry has been widely and successfully used in a multitude of applications: coastal environments [39][40][41], such as the monitoring of dune systems [31,42,43] and beach topography [44][45][46], river morphology characterization [47], assessment of glacier dynamics [48,49], landslide monitoring [50,51] and also architectural heritage [52,53]. The high level of accuracy achievable with DEMs derived from UAV imageries has been investigated and confirmed [37,54], as well as the reproducibility of the results [55,56].
In this work, we present the use of UAV photogrammetry for assessing the behavior of a micro-tidal system in terms of evolution, with regard to the restored salt marsh, the mudflat and the channel network. In fact, the potential of UAVs for monitoring the morphodynamics of mudflats has been investigated by Jaud et al. [57] in the Seine Estuary, finding out that UAV datasets enable, for instance, multi-scale approaches, but a complete micro-tidal system has not been investigated yet. In particular, the methodology proposed here has been applied focusing on a small cell (Cell 1B, abbreviated as C-1B in the following) [58], which is located in the middle of the overall tidal system described above and hence C-1B is connected with both Cells 1A and 1C. In fact, C-1B experiences all of the dynamics of CB shorelines. For instance, between 2 March and 10 March 2018 Winter Storm Riley caused a water level range of up to 2 m. This was recorded by tide gauges in Cambridge (Figure 2) along the Choptank River and in Annapolis. Both the tide gauges are about 30 km away from Poplar Island. Considering that the total (maximum) predicted (Tide predictions have been retrieved from NOAA (National Oceanic and Atmospheric Administration) at https://tidesandcurrents.noaa.gov accessed on 3 June 2019. Predictions reported here does not account for local meteorological effects.) tidal range during the winter storm (without storm effects) is approximately 0.7 m, this was a typical Nor'easter. In order to evaluate the evolution of the cell under SLR and meteorological events, a monitoring over time has therefore been started.
The aim of the proposed study is to apply a simple and flexible survey procedure to help and improve tidal system monitoring, also for providing basis for high-resolution numerical modeling of small tidal systems, exploring the potential of UAVs in this field. Moreover, the study site selected for the methodology development is one of the first worldwide restoration project where an interdisciplinary marsh design has been tested. Due to the small size of both the restoration project and the features of interest, traditional satellite remote sensing (e.g., Landsat 8, Sentinel 2 imagery) does not have a sufficient ground resolution and an innovative monitoring method is needed for tracking geomorphological and ecological changes on a short time scale and with a quick time frame. Originality consists of applying this approach to a small cell (micro-tidal system) in order to assess the evolution over time investigating the feasibility of the overall procedure and evaluating the actual level of accuracy attainable.

Materials and Methods
In this section, we present the survey methodology ( Figure 3) based on UAV photogrammetry together with ground GNSS RTK survey that was adopted for the reconstruction and characterization of the 3D geometry of C-1B. The morphology of the different parts of the cell (illustrated in Figure 4) was directly assessed through a comparison over time to identify possible erosion processes that are in progress. In order to achieve this result, we tested a methodological approach to collect geospatial information about the cell morphology: acquiring aerial imagery from UAV to evaluate the distribution of the vegetation and to generate Digital Terrain Models (DTMs) of the emerged areas, performing a GNSS RTK survey to integrate missing data in submerged channel areas and to collect the coordinates of Ground Control Points (GCPs). The procedures presented below can hence be replicated in practice in any similar environment all over the world.

Image Acquisition and Related In-Situ Operations
Cell mapping with a centimeter-level ground resolution can be performed using a camera mounted on a small UAV and flying at an appropriate altitude. In this methodological research, we used a DJI Phantom 3 Professional (DJI-P3P) multicopter equipped with the DJI FC300X camera for image acquisitions. In fact, this is a versatile, lightweight and easy-to-use drone. Aircraft and camera specifications [59] are reported in Table 1. The theo-retical Ground Sample Distance (GSD), that represents the pixel size on the ground, can be computed as: GSD = pixel size · flight altitude focal length (1) and hence depends on the actual flight altitude of the UAV. In order to obtain a GSD of ≈2 cm, we planned image acquisition at a 40 m altitude: in this way, the theoretical GSD is 1.8 cm. The first aerial survey was performed in April 2018 and a total amount of three flight missions were planned and executed for the survey of the entire extent of the cell. A second survey was also carried out one year later (in April 2019) using the same DJI-P3P drone. In this case, four different flight missions were necessary for the mapping of the cell, due to the need for limiting the flight time to about 10-12 min. Finally, a third aerial, the last for this methodological investigation, survey was also performed in November 2019 with the vegetation close to the end of its seasonal cycle. Photo density was found to be sufficient (with more than 9 images per ground point) to ensure a high number of reprojections of each ground point on the acquired images for each image dataset we processed. Flight plan parameters are shown in Table 2. With regard to water levels (tide conditions), the survey in April 2018 was performed in high tide (approximate water level of about +0.10 to +0.20 m NAVD 88), while the April 2019 and November 2019 in mid-low tide (approximate water level of about −0.10 to 0.00 m NAVD 88). In order to perform a model georeferencing with a final centimeter-level accuracy, a set of GCPs and Check Points (CPs) were deployed every time. Each GCP/CP consisted of a target of about 60 cm × 60 cm, clearly visible on the acquired images. Coordinates of the GCPs were collected using two Topcon HiPer V geodetic receivers in Real-Time Kinematic (RTK) mode. The base station was set up on a benchmark at the edge of the cell (see Figure 5), while the rover was used to record the coordinates of the center of each target with a 30-s Stop-and-Go. Figure 5 illustrates the location of GCPs and CPs for the survey repetitions of April 2018 and April 2019. In addition to GCPs and CPs, non-signalized Validation Points (VPs) were also surveyed in RTK mode for an accuracy assessment of the DTMs. All the coordinates were recorded in NAD 83(2011) reference system, while ellipsoidal heights were converted into orthometric elevations referred to the NAVD 88 vertical datum by applying the GEOID12B separation model. It is worth noting that the RTK mode used in this work can be replaced with a Network RTK (NRTK) in most situations. However, RTK mode can be easily used in any kind of environment, simply monumenting a local base benchmark.

Photogrammetric Reconstruction
UAV-acquired images were used to reconstruct the three-dimensional geometry of the cell by using an SfM-MVS approach. The image datasets acquired during aerial flights (see Section 2.1) were processed in Agisoft PhotoScan/Metashape Professional software. The overall process, from the stage of the acquisition planning (desired GSD, flight altitude and number of missions) to the generation of final products of the photogrammetric workflow, is depicted in Figure 6. The reconstruction of the geometry using SfM-MVS essentially consists in (i) detecting features on the images and describing them using, for instance, scale invariant feature transform (SIFT), then (ii) matching the features across the different images looking for correspondences and filtering out outliers and finally (iii) performing a Bundle Block Adjustment (BBA), also taking into account the camera model estimation (interior orientation parameters through a self-calibration). In order to increase the accuracy of the sparse point cloud made up of all the tie points detected on the images and successfully matched, those points with a low reprojection accuracy (threshold used in Agisoft PhotoScan/Metashape: 0.3) or a high reconstruction uncertainty (threshold used in Agisoft PhotoScan/Metashape: 10) were removed prior to continuing the data processing. In these initial stages, the GCPs were introduced (i.e., specified on every image and provided with their own surveyed coordinates). In this way, the photogrammetric model was georeferenced. Finally, non linearities in camera modeling were considered through a global optimization of the interior orientation parameters: focal length f , principal point offset c x , c y , affinity and non-orthogonality (skew) coefficients b 1 , b 2 , radial distortion coefficients k 1 − k 3 , tangential distortion coefficients p 1 − p 4 [60]. Processing parameters used in the software are reported in Table 3. Table 3. Processing parameters.

Alignment
Camera model Normal Accuracy High Key point limit 200,000 Tie point limit 20,000

Depth Maps and Dense Point Cloud
Quality High Depth filtering Aggressive The assessment of the actual modeling accuracy was carried out by using the CP dataset: the comparison between model-derived coordinates and surveyed coordinates provided an indicator of the overall level of accuracy of the 3D model. The residuals computed for the aerial survey performed in April 2018, for both GCPs and CPs, are shown in Table 4. Values obtained for the survey in April 2019 and November 2019 were found to be similar.
Dense point clouds were reconstructed within the same software. In this work, dense clouds were made up of 55,000,000 to 170,000,000 points for the overall extent of the cell, leading to an average point density spacing from a little less than 400 points/m 2 to almost 1200 points/m 2 . A further procedure to filtering out the vegetation was also applied using a built-in function of Agisoft PhotoScan/Metashape. This method splits the entire model into a grid of cells with a user-defined size, then assumes that the lowest point within each cell (5 × 5 m) belongs to the ground and finally classifies all the remaining points on the basis of angle and distance thresholds [60]. From the classified ground points, a georeferenced Digital Terrain Model (DTM) was generated and exported. Similarly, a centimeter-level orthomosaic was also generated and exported from each photogrammetric model. In addition, a VP-based validation of the DTM was also performed prior to analyzing morphology changes. Despite the application of a slope detection algorithm to classify ground points in Agisoft PhotoScan/Metashape, the DTM did not reach a centimeter-level accuracy in vegetated regions (i.e., salt marsh). Figure 7 shows the location of VPs and the difference between the GNSS elevation and the April 2019 DTM elevation (H GNSS − H DTM ). The results, divided per region (i.e., salt marsh, mudflat and channel edges) are reported in Table 5. The vertical error ranged from a few centimeters to a maximum of about 40 cm into the salt marsh. This represents a common and expected issue in marshes with dense vegetation [61][62][63]. In un-vegetated zones the differences were found to be lower, with RMSE values comparable to the expected accuracy of the DTM. Similar values were found for the other aerial surveys in this work. Hence, a centimeter-level vertical accuracy was obtained in un-vegetated regions (i.e., channel edges and mudflat). Figure 8 illustrates the location of a wider dataset of VPs on the mudflat and the distribution of elevation differences computed again as H GNSS − H DTM for each VP. The mean value of −0.8 cm and the standard deviation of 2 cm confirmed the level of accuracy mentioned above. Also, the standard deviation is still comparable with the 3D Root Mean Square Error (RMSE) of GCPs and CPs (see Table 4). Due to the high uncertainty of the DTM in vegetated regions, in this work we did not process elevation data of the salt marsh platform to track geomorphological changes.

Initial as-Built Dataset
Data collected by aerial flight missions and GNSS RTK was compared to the initial status of C-1B (in 2012) provided by the MES. The initial profile of a total of six cross-sections (see Figure 5) was hence reconstructed using the elevation data. Three cross-sections were selected along a main channel and the remaining were located on the secondary channels. Also, an as-built point dataset was used to evaluate elevation variations on the mudflat using a discrete (point-by-point) approach. Vegetation planting design was also provided to us.

Vegetation Boundary for Evaluating Channel Width
The high-resolution orthomosaic generated from UAV-acquired images [41,64] and the GNSS survey data were then used for evaluating the morphology of the channels. In this work, we evaluated the channel width for each monitored cross-section through two different approaches. Since the possibility to assess the channel width from a DTM depends strongly on the water level (tide conditions) at the time of each survey, here we propose an alternative method based on vegetation boundary to assess lateral variations as first approach.
The first aerial survey was performed in high tide conditions (April 2018), with a water level W L N AVD 88 ≥ 5-10 cm and hence with the channel banks almost completely submerged. This situation should be avoided in general. However, in practice, the possibility to acquire aerial images is often too much related to wind speed in the middle of the bay, thus waiting for low tide could be almost not feasible in most situations. This led us to consider the use of vegetation boundary on the salt marsh to evaluate lateral erosion along the scarps as well as the channel width. In fact, due to the vegetation that is present on the salt marsh platform and that is suddenly interrupted whenever these scarps are present, the erosion can be recognizable with almost any water level. It is also worth noting that the mapping using this approach was repeated twice in similar conditions of the seasonal growth of the vegetation (in April 2018 and April 2019, respectively). This ensured the best level of repeatability of the mapping operations based on UAV orthomosaics for channel width evaluation. A third repetition was finally conducted on November 2019 data.
The detection of the erosion scarp mentioned above consisted in a manual mapping carried out on the orthomosaics with a GSD of about 2 cm. Operator error magnitude (σ op ) can be estimated as a few GSDs (up to σ op = 10 cm) on each channel bank. Channel width error is hence σ w = 2 · σ op ≈ 0.15 m. Such a value appears to be quite high, but it is worth noting that vegetation loss is often at a level of clods, since the roots entrap a wide portion of soil. For this reason, we did not consider as significant width variations less than 2%. This threshold corresponds to a variation of σ w = 0.15 m over a total width of w = 7.5 m. Then, we computed the variation rate per year twice for the period 2018-2019 as: where the factor 12 19 takes into account the actual number of months. Also, we computed the average ∆w = ∆w 1 +∆w 2 2 .

GNSS RTK Survey: Cross-Section Profiles and Channel Width
In order to reconstruct the profile of cross-sections and integrate the aerial surveys, a direct GNSS RTK survey can be carried out. This method was used for surveying crosssection points located within the channels. This made it possible to acquire additional information in submerged regions to assess (i) the evolution of the bedlevel elevation and (ii) the shape of lateral slopes for further analyzing the effect of real shear stresses, making available, at the same time, (iii) a second approach to evaluate a more accurate channel width. Due to the high time consumption of this task, however, the GNSS RTK survey can be performed with a longer time frame (1-2 years). In general, the temporal resolution depends on the hydro-dynamic conditions and marsh type of the monitored site. The availability of as-built data allowed us to reconstruct the initial shape of the profiles for a comparison over time. The georeferencing of all of the surveys within the same NAD 83(2011) horizontal datum is a crucial aspect to ensure a consistent comparison of the cross-sections with respect to each channel axis, without any significant lateral offset. For the same reason, elevations were referred to the NAVD 88 vertical datum in this work.
In this research, we assumed a reference elevation at a water level W L N AVD 88 = 0 m. In fact, this elevation represents a good indicator for channel width in C-1B (see Figure 9) with respect to the predicted tide levels reported in Figure 2.  Table 7, surveyed in 2018 (GNSS) and 2012 (as-built).

Elevation Variations on the Mudflat
In addition to mapping channel edges and estimating channel widths, an assessment of the elevation variations on the mudflat can be carried out. This region experiences sedimentation and erosion whose effects are not reduced by the vegetation. In order to compare the morphology of the cell with its initial one (in 2012), the elevations extracted from the DTM were considered. In particular, we performed an analysis using the point dataset described in Section 2.3. An average 2018-2019 elevation (H) and its variation from the time of as-built data (i.e., 2012) was computed. This assumption is feasible for our purpose of evaluating the elevation change over ca. 6 years due to the centimeterlevel accuracy of DTMs and sedimentation/erosion rates on the mudflat not exceeding this threshold over a single year and also allowed us to exclude possible outliers from the comparison.

Vegetation Boundary for Evaluating Channel Width
As mentioned in Section 2.4, channel width was estimated for six cross-sections. Channel width estimates based on the vegetation boundary are reported in Table 6. Also, the variation rate ∆w per year is reported together with the width estimates. The location of each cross-section (with abbreviations) is shown in Figure 5. With regard to the main channel, MC-B seems to show a moderate enlargement. On the contrary, MC-C and MC-D can be assumed as stable considering the uncertainty of channel width using the vegetation boundary. The secondary channel cross-sections, instead, present different situations. SC-C shows a moderate enlargement in 1-year period, significantly extending itself from April 2019 to November 2019. SC-J and SC-F did not show significant variations. Both MC-B and SC-C enlargements were due to a loss of some clods of muddy soil with vegetation on the right side of the respective cross-section. This was confirmed by a visual inspection of the orthomosaics.

GNSS RTK Survey: Cross-Section Profiles and Channel Width
Channel width was also estimated in correspondence of a water level assumed as a reference, as previously mentioned in Section 2.5. Table 7 reports the results we obtained with this approach. It is worth noting that the values reported in Tables 6 and 7 cannot be compared since the approach we adopted was different (i.e., vegetation boundary vs. elevation). We can note that all of the cross-sections are wider than the time they were built. On average, the variation rate is up to ≈+6 % per year. Looking at the profiles (see Figure 9) directly surveyed in GNSS RTK mode and compared to the initial (and almost trapezoidal) shape derived from the as built data (2012), we can note that the shape is deeply changed, with a rounding also due to the properties of the fine and cohesive sediment which the cell is made up of. Both MC and SC cross-sections highlighted a lateral erosion on left and right sides as well as a deposition on the channel bed. In particular, the profiles of the MC cross-sections (MC-B, MC-C and MC-D) showed a deposition up to ≈40 cm in proximity of the center of the channel. Similarly, a value of about 20 cm has been generally detected for the secondary channels. Lateral erosion trends have been found to be more spatially variable across the monitored cross-sections. Most of the variations we found took place at elevations that are regularly exposed to the air and submerged during tide cycles (intertidal). Additionally, due to the low velocities in the channelized area, the cross-sections showed deposition at the bottom of the tidal creeks. An important role is might played by the sediment dynamics which is locally removed from lateral banks and deposited right at the center of the channel.

Elevation Variations on the Mudflat
The elevation variation on the mudflat, computed using H as reported in Section 2.6, is shown in Figure 10. Points presenting elevation losses are located close to the channel edges. This further confirms the presence of lateral erosion along the channels since the construction of C-1B (in 2012). The average elevation loss is 4.5 cm for the moat channel and 3.6 cm for the secondary channel. Points located near the vegetation on the mudflat showed a significant elevation increase (≈+7 cm on average). This is due to sediment trapping and organic matter production [65].
Conversely, the inner region of the mudflat is characterized by lower elevation changes and seems to be essentially stable in the last period. However, the mudflat area shows local signs of erosion which might highlight elevated variability in all the area. Rates of elevation change detected in this region (+0.5 cm on average) have the same magnitude of the model uncertainties that were evaluated through the quality assessment of photogrammetric model and DTM (i.e., using CPs and VPs). Such variations cannot be considered as significant.

Discussion
The results of the methodology adopted in this first case study on restored salt marshes located in a micro-tidal system allowed us to reconstruct the three-dimensional geometry of the cell (i.e., C-1B) with a centimeter-level ground resolution with the possibility of monitoring geomorphological variations with respect to previous conditions. The quality assessment of photogrammetric models is performed with CPs for the validation and confirmed the centimeter-level accuracy of the 3D reconstruction obtained in this work. The subsequent introduction of VPs (i.e., non-signalized points) is useful to evaluate the quality of the DTM and highlighted that only un-vegetated regions reached a centimeterlevel accurate elevation in this research, since the vegetation in the salt marsh was not adequately filtered by ground point classification algorithms. This is an expected issue whenever the vegetation is very dense. However, the overall survey accuracy of about 2 cm, due to both the GNSS survey in RTK mode and the GSD of the aerial datasets, allows one to evaluate the progress of channel morphodynamics and vegetation variations. Difficulties related to the presence of the vegetation could be overcome using, for instance, a different approach based on UAV-LiDAR systems to obtain vegetation characteristics [66]. In fact, a new algorithm has been applied to high-resolution UAV-based LiDAR point clouds for estimating ground elevation, vegetation height and vegetation density in coastal salt marshes, without the necessity of additional datasets (such as direct ground GNSS surveys). This method proved to produce maps of ground elevation and vegetation characteristics (of Spartina alterniflora) with a resolution of 40 × 40 cm 2 [67].
In general, two different approaches can be adopted for evaluating channel width variations in correspondence of selected cross-sections. In this work, we evaluated the channel width for a total of six cross-sections using (i) the high-resolution orthomosaics through an image-based mapping of vegetation boundary and (ii) the direct survey of the entire profile with a GNSS RTK approach. Width estimates based on GNSS RTK have been also compared to the as-built survey of the cell (2012) since both the elevation data have been referred to the NAVD 88 vertical datum for a consistent comparison. Using vegetation boundary for estimating channel width and monitoring lateral erosion trends proved to be a quick approach characterized by a low level of accuracy (σ w ≈ 0.15 m). However, the land loss due to the entrapping of Spartina alterniflora roots is often at a level of clods, and hence recognizable. Moreover, this approach can be also used in high tide conditions if the vegetation is not submerged. In order to integrate aerial surveys with data collected in the submerged regions of the channels, cross-sections can be also surveyed with GNSS RTK to reconstruct the corresponding profiles. This method is similar to the technique commonly used to integrate topography and bathymetry data, in which the vertical reference [68] is still crucial to ensure an accurate and consistent reconstruction of the overall morphology. However, this is a time-consuming task and its repetition can be performed less frequently than the mapping based on orthomosaics. Here, the comparison of cross-section profiles with as-built ones highlighted a deposition of sediments on the bed of both main and secondary channels. A general increase of about 20 cm has been found in the middle of the channels, with a maximum up to ≈40 cm. At the same time, the additional information provided by GNSS data has shown how the deposition on the bed of the channels mentioned above is also matched by a lateral erosion of the submerged slopes. This change of shape is correlated with the real shear stresses caused by ordinary tidal cycles [69]. In addition, the pronounced curvature of the channels can enhance the impact of the incoming and outgoing fluxes on certain channel banks. The channel width trends found are essentially consistent with orthomosaic-based mapping carried out on a 1-year reference period.
An elevation assessment has been also performed for the mudflat, with a focus on the recent stabilization. Its innermost region has been found to be stable and detected elevation variations were within the range of modeling errors. Conversely, points assessed on the edges of the mudflat, close to the channel edges, lost elevation. This fact further confirms the general lateral erosion highlighted by the monitoring over time performed in correspondence of the cross-sections. A moderate elevation increase was found near the vegetation on the mudflat. This marsh platform accretion represents the response to SLR [18]. For tidal flats it is important to highlight that the presence of water ponds and/or water layers can interfere with the DTM retrieval. Flying when the water level is increasing can help to mitigate and prevent this issue, leaving time to small water ponds to disappear if possible. Conversely, flying when the water level is decreasing gives the possibility to survey more emerged areas within the channels and on the channel banks. Hence, the optimal solution for performing UAV surveys should be carefully evaluated on the basis of which kind of data has to be collected.
Even if we still have a partial understanding of the forces driving tidal marsh resilience, the methodological imprint presented in this study, along with its preliminary results obtained in a micro-tidal system on Poplar Island (USA), establishes that channel design and vegetation growth are key processes for salt marsh restoration and their assessment over time is crucial. Belliard et al. [70] state that the impacts of vegetation on sediment deposition might be enough for marshes to withstand low amounts of SLR, however, numerical models forecast wide amount of marsh loss in tidal system with moderate to high levels of SLR (>10 mm/yr) and hence restoration interventions should be deeply investigated with a short time frame to establish the most efficient way to restore salt marshes.
Therefore, it is important that numerical modeling, tidal system monitoring and restoration actions improve rapidly enough to suggest actions to counteract the most severe SLR projections. Fleri et al. [27] showed that vegetation characteristics and channel stability are integrally linked.
Recently, restoration projects by ecological engineers in Black Water National Refuge in Maryland (USA) are focusing on channel digging in natural marshes as a solution for dredging swamped area under erosion. Those recent interventions show the crucial role of channels network, and additionally, demonstrate how excavating new channels as an extension of a natural tidal creek will connect the flooded parts of the marsh to the bay draining the entrained water.

Conclusions
This work highlights how the use of UAV+SfM-MVS and GNSS RTK surveying allows the monitoring over time of a micro-tidal system. The very high spatial resolution of orthomosaics and DTMs, especially if further combined with a high temporal resolution in terms of small-time lapse between survey repetitions, enables the assessment of cell morphology with a centimeter-level ground resolution and accuracy on tidal flats, while issues were encountered for salt marshes due to the presence of vegetation. Channel width, that is an indicator of lateral erosion, can be assessed with an uncertainty of ≈0.15 m from orthomosaics, enabling a quick mapping at almost any water level (tide conditions). More accurate data can be collected through a GNSS RTK survey of channels, but this task is very time-consuming and hence can be used less frequently.
Because of the low operative costs in using UAVs, survey repetitions can be made whenever needed and they will provide additional data to assess the behavior of tidal systems in the future. This information could be also considered by managers as a reliable tool in marsh creation and restoration projects. In fact, from a geomorphological point of view, restored marsh and channel systems rapidly develop up to the order of months, requiring high resolution and frequency of surveys. In particular, our study site well shows interesting results which might help management and design of future restoration projects. The methodology presented in this work can hence be replicated in any micro-tidal system in order to assess the response of restoration interventions to ordinary tide cycles and extreme events, both for geomorphology and vegetation tracking.
Further improvements of the research presented in this work is likely represented by the use of an aircraft with on-board RTK [71] that limits the need for a wide set of GCPs. The same DJI Phantom 3 Professional we used for the aerial surveys in this work is already available in a more recent version that integrates such a technology (DJI Phantom 4 RTK). Moreover, the advantages of GNSS kinematic modes are also available in post-processing [72] with a final centimeter-level accuracy and the acquisition of oblique images [73] leads to higher vertical accuracies. This can be a crucial aspect for improving and developing a precise monitoring of the mudflat. The combination of on-board RTK with oblique images is also able to provide consistent estimates not only for the exterior orientation, but also with regard to the interior orientation [74,75]. In fact, the accurate estimation of the focal length is crucial for achieving the best vertical accuracies [76], in particular when using UAV image datasets.
Finally, numerical simulations of hydro-dynamics and sediment transportation within such a small, but quite complex, tidal system, based on high-resolution elevation models will explain the driving factors of the morphological evolution enabling a better understanding of restoration interventions of salt marshes for a better management of coastal environments.

Abbreviations
The following abbreviations are used in this manuscript: