# Morphological evolution of creek networks in 10 restored coastal wetlands in the UK

Clementine Chirol, I D Haigh, N. Pontee, C E L Thompson, S L Gallop

► To cite this version:

Clementine Chirol, I D Haigh, N. Pontee, C E L Thompson, S L Gallop. Morphological evolution of creek networks in 10 restored coastal wetlands in the UK. Scientific Data, 2022, 9 (1), pp.144. 10.1038/s41597-022-01199-4. hal-04276319

HAL Id: hal-04276319

https://hal.inrae.fr/hal-04276319

Submitted on 8 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
scientific **data**

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# Morphological evolution of creek networks in 10 restored coastal wetlands in the UK

C. Chirol<sup>1,2</sup>✉, I. D. Haigh<sup>2</sup>, N. Pontee<sup>2,3</sup>, C. E. L Thompson<sup>2,4</sup> & S. L. Gallop<sup>5,6</sup>

Coastal wetlands provide crucial ecosystem services including flood protection and carbon storage, but are being lost rapidly worldwide to the combined effects of sea-level rise, erosion and coastal urbanisation. Managed Realignment (MR) aims to mitigate for these losses by restoring reclaimed land to tidal influence. Data of creek evolution is critical to assess the performance of design strategies and improve design and implementation practices. This data descriptor provides a dataset of the horizontal morphological evolution of creek systems from various initial conditions in 10 MR schemes across the UK. Using a semi-automated workflow, morphological creek parameters were extracted from 52 lidar datasets at 1 m horizontal resolution spanning 2 to 20 years post-breach. This constitutes the most comprehensive systematic monitoring of MR creek morphology to date. The dataset will assist future MR design and provide baseline morphological information for ecological and biogeochemical surveying.

## Background & Summary

Increasing focus on nature-based solutions has led to a rise in coastal wetland restoration projects worldwide<sup>1–3</sup>, to the point where the UN has declared 2021–2030 the Decade of Restoration. Coastal wetlands are being degraded at a rapid pace<sup>4–6</sup>, resulting in the loss of critical ecosystem benefits including biodiversity<sup>7</sup>, flood protection<sup>8</sup>, nursery habitats for juvenile fish<sup>9</sup>, pollutant filtering<sup>10</sup>, and carbon storage<sup>9</sup>. Managed Realignment (MR) schemes, in which embanked lands are opened to tidal influence to create new saltmarsh and mudflat habitats (Fig. 1a), are important sources of experimentation and innovation for coastal restoration strategies and have brought to light design aspects that require further scientific guidance<sup>11–13</sup>. Notably, creek networks play a crucial role in the distribution of water, sediment, nutrients and seeds through the marsh<sup>14–16</sup>, and it is thought that excavating initial channels may speed up creek network development, thus helping to support saltmarsh functioning and meet target conditions<sup>11,17–19</sup>. Of particular interest is the capacity of the initial morphology of MR creeks to develop horizontally and improve their distribution over the marsh.

However, even though a large variety of initial creek designs have been used in the UK and around the world<sup>11,20</sup>, integrated, standardised data sets over time for multiple sites are lacking if we are to compare the impact on creek evolution of different design strategies. This is mainly due to the computational challenge of mapping and parametrising complex features like creek networks. This data gap hinders our understanding of how both the physical and ecological aspects of MR schemes evolve, and hence the formulation of general design and implementation guidelines.

In this data descriptor, we apply a creek mapping algorithm<sup>21</sup> to 10 MR sites to monitor creek network evolution after implementation (Fig. 1c). The algorithm uses freely-available elevation maps collected from an airborne sensor (lidar) to extract relevant volumetric and topological information (Fig. 1b), and allows for the systematic comparison of non-engineered and highly engineered MR sites. The algorithm was tested on 13 natural saltmarshes that have been morphologically stable for over 20 years, where stability was defined as when geomorphic changes within the creek network within a time period are lower than the data resolution, and so

<footnote n="1"><sup>1</sup>LSE, University of Lorraine, Nancy, 54000, France. <sup>2</sup>Ocean and Earth Sciences, National Oceanography Centre Southampton, University of Southampton, Southampton, SO14 3ZH, UK. <sup>3</sup>Jacobs, One Glass Wharf, The West Wing, Bristol, BS2 0EL, UK. <sup>4</sup>Channel Coastal Observatory, National Oceanography Centre, Southampton, SO14 3ZH, UK. <sup>5</sup>School of Science, University of Waikato, Tauranga, 3110, New Zealand. <sup>6</sup>Environmental Research Institute, University of Waikato, Hamilton, 3240, New Zealand. ✉e-mail: clementine.chirol@univ-lorraine.fr</footnote>
a: creek design

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

b: creek mapping

<figure><img src="image_6.png" /></figure>

c: creek parameterisation

<figure><img src="image_7.png" /><figcaption>**Fig. 1** Post-implementation monitoring schemes of creek networks in managed realignment schemes as provided by this data descriptor. (a) Initial creek pattern following construction; (b) Semi-automated creek mapping from lidar using elevation and slope thresholds; (c) Creek parametrisation and morphological evolution of creek systems inferred from consecutive lidar datasets.</figcaption></figure>

are considered to have reached a mature state²¹. The results were consistent with those found by manual mapping and ground-truthing in the field.

The new dataset provides crucial information on the baseline (pre-breach) morphology of coastal wetland restoration schemes, and how they might evolve in the future. The outputs of this study are twofold: first, topographic maps showing the extent and topology of the creek network for reuse in further studies; second, result tables of creek network morphometric parameters and a series of interpretative figures. The dataset will assist future managed realignment design by showing how current creek systems have evolved since implementation for different environmental conditions and levels of engineering, as well as providing baseline morphological data for further research. Knowledge of the morphology and spatial distribution of creeks is particularly relevant for ecological or biogeochemical surveys. This dataset will therefore support the study of various saltmarsh processes such as headward erosion²², seed dispersal²³, development of vegetation patches²⁴, tidal flow and sediment transport²⁵, marsh flow attenuation²⁶, and carbon storage²⁷, and facilitate the interpretation of local observations at the scale of the entire saltmarsh. There is a strong need to conduct such surveys because current MR schemes often fail to provide the same ecosystem services as their natural counterparts²⁸,²⁹, and do not replicate the full plant diversity of natural wetlands³⁰. If not addressed properly, this could lower their value as wetland loss mitigation schemes.

Methods

The extraction of MR creek evolution datasets involved three main stages, explained in the details below:

**Stage 1: Site selection and raw data acquisition.** The first stage consisted of collating raw datasets that provide geomorphological and tidal information for MR sites, as well as background information on the context and design choices of these MR sites to inform site selection.

**Lidar data.** Lidar is an airborne sensor that measures distances by transmitting a laser pulse and calculating the return time of the reflected beam. Lidar is particularly efficient for the geomorphological monitoring of coastal environments³¹,³². The use of lidar is increasing due to its accessibility in many countries, high coverage (e.g., all of the UK territory), high horizontal resolution (up to 0.25 m horizontally), and the possibility to compare datasets over several years to infer evolution rates.
The lidar data used in this data descriptor were collected by the Environment Agency (EA) and accessed via their data services platform³³. The datasets were interpolated from a point cloud into a raster using a nearest neighbour interpolation to Euclidean distance method. Contrary to other interpolation methods, nearest neighbor preserves variations in the data such as small (~1 m) channels which would otherwise be smoothed over³¹, making this method particularly suited to representing creek networks³⁴. The freely available data comes in two forms: a raw Digital Surface Model (DSM) interpolated from the points of first laser beam return, and a Digital Terrain Model (DTM) where the above-ground features such as buildings and vegetation have been filtered out to obtain a bare-Earth model using EA proprietary algorithms³⁵.

Horizontal resolutions vary between 0.25 m and 2 m. Due to the fractal behaviour of creek networks³⁶, their morphological complexity increases with the horizontal resolution. Thus, in order to compare the morphological evolution of the site at different years, all datasets need to be at the same resolution. The most common available horizontal resolution, 1 m, is therefore used as a standard for this study. The vertical resolution of lidar data is estimated at 0.15 m or lower based on GPS ground-truthing surveys, as described in the EA's lidar quality control reports³³.

**Tide data.** In order to establish tidal forcing parameters and to estimate the tidal prism, mean predicted tidal levels were obtained for 582 standard (tidal data tabulated) and secondary (tidal data calculated from the standard ports) ports, as provided by the Admiralty Tide Tables 2014³⁷. The mean tidal levels at the 13 natural saltmarshes and 10 MR sites were then interpolated from the weighted mean of the surrounding ports' values (up to 30 km away). Tidal levels were converted from Chart Datum to Ordnance Datum Newlyn, using conversion values in the official Tide Tables and from the National Tidal and Sea Level Facility (https://ntslf.org/). The key tidal levels that were extracted are virtually stationary within the 18.6 year nodal-cycle. Therefore, the effect of sea-level rise was considered negligible in this study which spans a maximum of 20 years of creek evolution. Observed rates of sea-level rise around the English Channel in the 20th century range between 0.8 and 2.3 mm/yr³⁸, and have been estimated at 1.4 ± 0.2 mm/yr around the UK³⁹, so the change in tidal elevation over 20 years is unlikely to exceed 0.05 m, which lies below the vertical resolution of lidar. The following mean tidal levels are interpolated for each site: Highest Astronomical Tide (HAT), Mean High Water Spring (MHWS), Mean High Water Neap (MHWN), Mean Water Spring (MWS), Mean Low Water Neap (MLWN), and Mean Low Water Spring (MLWS).

**Site selection.** This study considers 10 MR schemes, implemented between 1995 and 2014 around the coast of England (Fig. 2). Sites were selected based on four criteria the first of which was data availability, i.e. the number of lidar datasets and case studies available. The second criterion was the sites' settings, which should capture a range of initial external conditions, i.e. the tidal range, size of scheme, location, implementation date, land use history, etc. The third criterion was the general scheme design, i.e. number and size of breaches, initial site elevation within the tidal frame, targeted habitats. The fourth criterion was the creek network design, i.e. the absence of initial creeks (strategy 1), excavation of a creek system from a natural template (strategy 2), or excavation of artificial creeks in the absence of a natural template (strategy 3). The context for each MR scheme was taken from the academic and grey literature (see all references in Table 1), and from an online database referencing coastal habitat creation schemes in North-Western Europe (ABPmer OMReg⁴⁰). The schemes span tidal ranges from 4.1–11 m, granulometries from clay to coarse silt / sand, suspended sediment concentrations from 60–1000 mg/L, mean elevations from 0.9–5.5 m above Mean Water Spring, and scheme areas between 0.1 and 3.61 km².

For each site, we collected 3 to 8 DSMs, surveyed and flown between 2002 and 2016 and spanning 2 to 20 years post-breach depending on data availability. No 1 m resolution lidar dataset could be found before 2002. Furthermore, changes in lidar technology mean that older datasets tend to be noisier. Therefore, while lidar datasets of 2 m horizontal resolution were generally omitted, an exception was made for Allfleet in 2007 in order to quantify the creek network close to the implementation date, and because the data was less noisy than the 2002 1 m resolution dataset taken at Freiston. A summary of selection criteria, lidar data availability and external references is given in Table 1 and the geographical repartition of the selected MR sites is given in Fig. 2.

**Stage 2: Preprocessing.** The second stage consisted of a series of preprocessing protocols applied using ArcGIS 10.2.2 to turn the lidar datasets obtained from the EA into the correct format of input files for the creek parametrisation algorithm.

Preprocessing included merging mosaics into a single dataset, interpolating to 1 m horizontal resolution, interpolating data gaps to the values of the nearest neighbours according to Euclidean distance (lidar data are generally collected at low tide when most creeks are drained, but remnant water within ponds for instance may lead to gaps in the dataset), cropping to the saltmarsh area, and extracting elevation and slope maps. The landward limit of the saltmarsh area is defined by the local HAT level which delimits the intertidal zone from the land, and the seaward limit as the mouth of the entry channel, consistent with Steel (1996); both are generally constrained by flood defences in managed realignment schemes. Creek edges tend to be more visible on the slope map than on the curvature map for the selected datasets, so the slope was chosen as a threshold parameter unlike previous studies⁴¹.

The preprocessing steps are minimalistic to provide monitoring tools that are easily reusable by coastal habitat restoration project designers and researchers for future saltmarsh monitoring efforts. The creek parameters are detected from freely available lidar DSM that have undergone minimal preprocessing, as is likely to be the case for most MR monitoring work performed by environmental agencies, contracted consulting companies, or by research projects in non-geomorphological disciplines that use creek morphology as a baseline. The two outputs of this preprocessing stage, an elevation map and a slope map in degrees, both at a horizontal resolution of 1 m, are converted into text files and exported to Matlab for the processing stage to extract relevant creek morphological parameters.
<figure><img src="image_8.png" /><figcaption>**Fig. 2** Location and aerial photography of the 10 MR sites considered. HAT: Highest Astronomical Tide interpolated along the British coastline from Admiralty Tide Tables 2014 mean predicted tidal levels. Red dashed lines: catchment area delimited by the remaining seawalls and HAT levels. Red arrows: breach areas.</figcaption></figure>

**Stage 3: Creek mapping, parameters extraction and visualisation.** In the third stage, the mean tidal levels and elevation and slope maps obtained in Stage 2 are used as input parameters for a coastal wetland creek parametrisation tool developed in 2014–2018 in collaboration with Jacobs. This tool outputs an Excel table of creek network morphometric characteristics and several figures. The algorithm, written in Matlab R2015a, is made available within the data repository along with the datasets. Readers can also refer a previously published methodology paper<sup>21</sup> which details the algorithm's functioning and validation process. A summary of the algorithm's running principles is given below.
<table><caption>**Table 1.** Selection criteria and data availability for each site. Creek design strategies considered: 1 = absence of initial creeks; 2 = excavation of a creek system from a natural template; 3 = excavation of artificial creeks in the absence of a natural template.</caption><thead><tr><th></th><th>Latitude/ Longitude</th><th>Site opening date</th><th>EA lidar datasets suitable for use between 1995-2016</th><th>Creek design strategy</th><th>References</th></tr></thead><tbody><tr><td>Abbots Hall</td><td>52.7846<br>0.8455</td><td>2002</td><td>2002; 2008; 2013; 2014; 2015</td><td>3</td><td>Blott and Pye<sup>57</sup> Shepherd<sup>58</sup> ABPmer White Papers<sup>59</sup> Scott et al.<sup>60</sup> Pendle<sup>61</sup> Brooks et al.<sup>45</sup></td></tr><tr><td>Alkborough</td><td>53.6914<br>-0.6773</td><td>2006</td><td>2007;2010; 2012; 2015</td><td>3</td><td>Manson et al.<sup>62</sup> Costa<sup>63</sup> Medlock et al.<sup>64</sup> Morris<sup>12</sup> Pendle et al.<sup>61</sup> Luisetti et al.<sup>9</sup> Pontee (in Esteves)<sup>65</sup> Pontee<sup>11</sup></td></tr><tr><td>Allfleet</td><td>52.6163<br>0.8364</td><td>2006</td><td>2007; 2011; 2013; 2015</td><td>3</td><td>ABPmer White Papers<sup>66</sup> Dixon<sup>67</sup> Medlock et al.<sup>64</sup> Pendle<sup>61</sup></td></tr><tr><td>Chowder Ness</td><td>53.6916<br>-0.4815</td><td>2006</td><td>2007; 2009; 2010; 2011; 2012; 2013; 2015; 2016</td><td>1</td><td>ABPmer White Papers<sup>68</sup> Luisetti<sup>9</sup> Costa<sup>63</sup> Morris<sup>12</sup> Pendle<sup>61</sup></td></tr><tr><td>Freiston</td><td>52.9646<br>0.0923</td><td>2002</td><td>2002; 2009; 2011; 2013; 2014</td><td>3</td><td>Friess<sup>69</sup> Friess et al.<sup>70</sup> Frost et al.<sup>71</sup> Hampshire<sup>72</sup> Nottage et al.<sup>73</sup> Rotman<sup>74</sup> Symonds et al.<sup>75</sup> Symonds et al.<sup>76</sup></td></tr><tr><td>Hesketh Out Marsh West (HOMW)</td><td>53.7203<br>-2.8876</td><td>2008</td><td>2009; 2010; 2011; 2014</td><td>2</td><td>Hampshire<sup>72</sup> Pontee (in Esteves)<sup>65</sup> Tovey et al.<sup>77</sup></td></tr><tr><td>Paull Holme Strays</td><td>53.7082<br>-0.2193</td><td>2003</td><td>2007; 2010; 2012; 2013; 2014</td><td>3</td><td>Clapp<sup>78</sup> Costa<sup>63</sup> Edwards et al.<sup>79</sup> Luisetti<sup>9</sup> Mazik et al.<sup>80</sup> Morris<sup>12</sup> Pendle<sup>61</sup> Pontee et al.<sup>81</sup></td></tr><tr><td>Steart</td><td>52.1983<br>-3.0506</td><td>2014</td><td>2014; 2015; 2016</td><td>3</td><td>Burgess et al.<sup>82</sup> Townsend<sup>83</sup> Wright et al.<sup>84</sup></td></tr><tr><td>Tollesbury</td><td>52.7673<br>0.8402</td><td>1995</td><td>2002; 2009; 2012; 2013; 2014; 2015; 2016</td><td>3</td><td>Atkinson et al.<sup>85</sup> Frost et al.<sup>71</sup> Garbutt et al.<sup>86</sup> Luisetti<sup>9</sup> Paramor et al.<sup>87</sup> Pendle<sup>61</sup> Reading et al.<sup>88</sup> Shepherd et al.<sup>58</sup> Steel et al.<sup>89</sup> Watts et al.<sup>90</sup></td></tr><tr><td>Welwick</td><td>53.6471<br>0.0095</td><td>2006</td><td>2007; 2009; 2010; 2011; 2012; 2013; 2014</td><td>1</td><td>ABPmer White Papers<sup>91</sup> Costa<sup>63</sup> Frost et al.<sup>71</sup> Hampshire<sup>72</sup> Luisetti<sup>9</sup> Morris<sup>12</sup> Pontee (in Esteves)<sup>65</sup> Pontee et al.<sup>81</sup></td></tr></tbody></table>

The creek algorithm is based on the threshold method: a creek network is defined as a connected feature which lies lower than the rest of the saltmarsh, and whose edges are delimited by a steeper slope (Fig. 3, Step 1). Once a raw creek logical mask has been detected, noise is filtered out by removing all connected elements smaller than a number of pixels defined by the user, while fragmented terminal channels are reconnected to the creek system using the shortest Euclidean distance as the repair path (Fig. 3, Step 2). Morphological thinning is then applied to shrink the creek network to a skeleton corresponding to the centerline of the channels (Fig. 3, Step 3). The topology of the creek network, expressed quantitatively in terms of reverse Strahler order<sup>21</sup>, is automatically computed by using a pruning process and assigning an order iteratively to all creek segments (Fig. 3, Steps 4 and 5). Compared to traditional Strahler ordering, reverse Strahler order ensures that the entry channel is always classified as the first order. The following parameters are calculated for each creek segment: sinuous length, straight length, sinuosity ratio, junction angle, and cross-sectional width, depth and area (Fig. 3, Step 6). The overall size and distribution of the creek system is given by the drainage density (the total channel length divided by the studied marsh area), overmarsh path length (OPL, the mean distance to the creek system at all points within the marsh), main channel length (the longest channel connected to the largest outlet), total channel length, number of creeks, total mouth cross-sectional area (sum of all outlets' cross-sectional areas), main channel mouth depth, planform area, creek volume, sinuosity ratio and the main channel gradient.

The algorithm is faster and less subjective than manual mapping, and interactive interfaces make it easy to use by researchers and stakeholders. Furthermore, the short running time at the pixel size used for the study allows for a number of sensitivity tests to refine the threshold parameters. The algorithm therefore allows the extraction of a systematic dataset that is comparable and can be produced quickly across multiple sites and years. Finally, it is well suited for the comparison of creek evolution in non-engineered and highly engineered sites, where creeks shaped by human intervention are poorly detected by the flow accumulation method<sup>42</sup>. A summary of the data types and processing protocols used in this study is provided in Table 2. The algorithm gives as outputs a number of morphometric parameters for the creek network (Table 3), following previous recommendation guidelines for creek design<sup>43</sup>.

## Data Records

The database presented herein consists in a set of tables, figures and text files containing creek network morphological characteristics for each of the 10 MR sites. The data and the algorithm used to generate it are accessible via ePrints (https://eprints.soton.ac.uk/434946/)⁴⁴ and via the Channel Coastal Observatory website (https://coastalmonitoring.org/ccoresources/education/chirol/). The complete list of data made available in the repository is as follows (summarised in Table 3):

Supporting Material SM1 is a table listing the initial conditions, design choices and available lidar datasets for each MR scheme considered (Stage 1). The table contains for each scheme: Implementation date, Latitude/Longitude, Catchment area, Time embanked (with references), Policy context, Site aims, Number of breaches (if MR), Initial creek system, Design details, Other references, available lidar data at 1 m resolution. The background information has been collected from various sources as listed in Table 1.
<figure><img src="image_9.png" /><figcaption>**Fig. 3** Creek parametrisation algorithm workflow. The 6 processing steps are grouped into three phases: creek detection (step 1); creek repair (steps 2 and 3) and parameter extraction (steps 4 to 6). The steps where user inputs (UIs) are necessary are marked as UI 1 to 4. The algorithm's functioning and validation process is detailed in a separate publication²¹.</figcaption></figure>

Supporting Material SM2 is a table listing the tidal data obtained for each site, including: Highest Astronomical Tide, Mean High Water Spring, Mean Low Water Spring, Mean Spring Tidal Range, Mean High Water Neap, Mean Low Water Neap and Mean Neap Tidal Range (**Stage 1**).

Supporting Material SM3 provides 52 text files containing XYZ data that cover the extent and elevation of UK managed realignment schemes, for all study sites and all available years, after implementing preprocessing protocols 1–4 (**Stage 2**).
<table><caption>Table 2. Summary of procedures for each data type.</caption><thead><tr><th>Type</th><th>Lidar DSM</th><th>Tide data</th></tr><tr><th>Source</th><th>Environment Agency</th><th>Admiralty Tide Tables 2014</th></tr></thead><tbody><tr><td>Preprocessing protocol 1</td><td>Mosaic merging (Merge, Data Management, ArcGIS)</td><td>Mean predicted tidal levels obtained for 582 standard (tabulated) and secondary ports (calculated from the standard ports)</td></tr><tr><td>Preprocessing protocol 2</td><td>Interpolate to 1 m horizontal resolution (Aggregate, Spatial Analyst, ArcGIS)</td><td></td></tr><tr><td>Preprocessing protocol 3</td><td>Gap interpolation (Nibble, Spatial Analyst, ArcGIS)</td><td></td></tr><tr><td>Preprocessing protocol 4</td><td>Cropping to relevant area delimited by HAT (Polygon feature, ArcGIS)</td><td></td></tr><tr><td>Preprocessing protocol 5</td><td>Slope maps acquisition (Slope, Data Analyst, ArcGIS)</td><td></td></tr><tr><td>Preprocessing protocol 6</td><td>Convert to ASCII (Raster to ASCII, Conversion, ArcGIS)</td><td></td></tr><tr><td>Processing protocol 1</td><td>Creek detection (Fig. 3, Phase 1)</td><td>Interpolation from the weighted mean of the surrounding ports' values up to 30 km away</td></tr><tr><td>Processing protocol 2</td><td>Creek repair (Fig. 3, Phase 2)</td><td></td></tr><tr><td>Processing protocol 3</td><td>Parameters extraction (Fig. 3, Phase 3)</td><td></td></tr><tr><td>Outputs</td><td>Creek extent mask<br>Creek order skeleton maps<br>Creek morphological parameters</td><td>Interpolated mean tidal levels (HAT, MHWS, MHWN, MWS, MLWN, MLWS)</td></tr></tbody></table>

<table><caption>Table 3. Summary of data record (access: https://eprints.soton.ac.uk/434946/).</caption><thead><tr><th>Name</th><th>Output Type</th><th>Provenance</th><th>Experimental manipulations</th><th>Content</th></tr></thead><tbody><tr><td>SM1</td><td>MR information</td><td>Academic and grey literature, online database ABPmer OMReg</td><td>Qualitative data collation from various sources</td><td>Time embanked, Policy context, Site aims, Number of breaches, Initial creek system, Design details, References</td></tr><tr><td>SM2</td><td>Mean tidal levels</td><td>Admiralty Tide Table 2014</td><td>Tide data preprocessing protocol 1, processing protocol 1</td><td>HAT, MHWS, MHWN, MWS, MLWN, MLWS</td></tr><tr><td>SM3</td><td>XYZ area maps</td><td>Environment Agency</td><td>Lidar data preprocessing protocols 1-4, processing protocol 1</td><td>textfiles containing XYZ data that cover the area and topography of the whole managed realignment schemes</td></tr><tr><td>SM4</td><td>XYZ creek maps</td><td>Environment Agency</td><td>Lidar data preprocessing protocols 1-6, processing protocols 1-2</td><td>textfiles containing XYZ data that cover the area and topography of creek networks</td></tr><tr><td>SM5</td><td>XYC creek skeleton maps</td><td>Environment Agency</td><td>Lidar data preprocessing protocols 1-6, processing protocols 1-3</td><td>textfiles containing XYC data that give the skeletonized area and Reverse Strahler order of creek networks</td></tr><tr><td>SM6</td><td>Morphometric parameters and figures</td><td>Environment Agency, Admiralty Tide Table 2014</td><td>Lidar data preprocessing protocols 1-6, processing protocols 1-3</td><td>PDF file containing:<br>-table of creek network morphometric parameters (Reverse Strahler order, Number of creeks per order, Total Length, Mean Length, Bifurcation Ratio, Sinuosity Ratio, Mean junction angle, Mean channel width, Mean channel depth, Mean cross-sectional area, A/D, W/D, Drainage density, Total Channel Length, Catchment area, Mean elevation above MWS, Overmarsh Path Length)<br>-Figures of cross-sectional area evolution of the largest entry channel mouth for all MR sites and all available years.<br>-Figures of MR creek morphometric parameters evolution per RS order over the years, plotted against the 95% spread of natural creek parameters.<br>-Figures showing the evolution of the MR creek extent, unchanneled length and overmarsh path length over the years for all sites<br>-Figures showing the elevation gains and losses of all MR sites between the first and the last year considered, correlated to initial site elevation and distance to creeks.</td></tr><tr><td>CHIROL CREEK EXTRACTION DEMO</td><td>Custom code used to generate the data</td><td>CHIROL-CREEK-ALGORITHM is a semi-automated coastal wetland creek parametrisation tool developed in 2018 by Clementine Chirol as part of her PhD at the University of Southampton, from 2014 to 2018, in collaboration with Jacobs. This tool uses a semi-automated elevation and slope thresholds method to detect a creek network, and gives as outputs an Excel table of morphometric characteristics and several figures. The scripts are run on MatLAB2015a.</td><td></td><td></td></tr></tbody></table>
Supporting Material SM4 provides 52 text files containing XYZ data that cover the extent and elevation of creek networks as detected by lidar and our developed algorithm for all study sites and all available years (Stage 3).

Supporting Material SM5 provides 52 text files containing XYZ data that cover the skeletonized creek network area as detected by lidar and our developed algorithm for all study sites and all available years, and gives the Reverse Strahler order of each branch (Stage 3).

Supporting Material SM6 provides tables and supporting figures that describe the morphological evolution of creek networks in MR schemes following implementation (Stage 3). They cover a broad range of parameters to facilitate their reuse for future studies. These are concatenated in a PDF file, and numbered Appendix A to E:

* Appendix A is a table listing all results from the creek network morphometric analysis performed at 10 MR schemes in the UK, utilising the semi-automated creek parametrisation algorithm for each available lidar dataset.

* Appendix B consists of 10 figures showing the cross-sectional area evolution of the largest entry channel mouth for all MR sites and all available years.

* Appendix C consists of 10 figures showing the evolution of MR creek morphometric parameters per Reverse Strahler order over the years, plotted against the 95% spread of natural creek parameters.

* Appendix D consists of 10 figures showing the evolution of the MR creek extent, unchanneled length and overmarsh path length over the years for all sites

* Appendix E consists of 10 figures showing the elevation gains and losses of all MR sites between the first and the last year considered, and correlates them to the initial site elevation and distance to creeks to analyse creek forming processes. The effect of creek proximity on the creek network is evaluated up to 20 m away from the creeks to reduce the impact of site elevation or of multiple creek influence on the results. Positive elevation changes are correlated to the most recent creek extent and negative elevation changes to the initial creek extent.

Finally, **CHIROL_CREEK_EXTRACTION_DEMO** contains the code used to generate the dataset, in the form of a list of Matlab functions and scripts. The associated **READ_ME** file contains a utilisation guide detailing the order in which the three primary Matlab functions should be called.

## Technical Validation

The various validation tests described below were applied to the lidar datasets and to the creek algorithm results.

**Lidar data validation.** The following quality control procedures were undertaken on the lidar datasets to justify carrying out a multi-annual monitoring of coastal wetlands using lidar:

1. Tests of the repeatability of lidar data collection over different years;

2. Investigations of the potential systematic bias caused by the presence of low-lying saltmarsh vegetation in the DSM and DTM; and

3. Comparison of the respective efficiency of the DSM and DTM at detecting the marsh surface

**Repeatability of lidar.** The repeatability of lidar data was tested for 5 different datasets collected near HOMW in 2007, 2009, 2010, 2011 and 2014. Elevation values were found using the nearest neighbour interpolation method along a low-lying unmoving feature, in this case a road (Appendix A1). The largest vertical repeatability error estimated from the mean standard deviation for the five points is +/-0.06 m, which is less than the 0.15 m vertical resolution of lidar.

**Effect of vegetation.** The effects of saltmarsh vegetation can lead to an overestimation of the marsh surface elevation by lidar and a systematic bias in the dataset. This effect needs to be quantified. Also, the capacity of the EA proprietary vegetation removal algorithms at minimising this bias needs investigating. Real-Time Kinematic (RTK) GPS data were collected with an accuracy <2 cm and a precision <1.5 cm in August 2015 in a MR scheme at Tollesbury and in an adjacent natural mature saltmarsh (Fig. 4). The marsh surface elevation data relative to ODM were collected across a 50 x 50 m section following the systematic grid sampling method described by Brooks et al.⁴⁵ and expanded by Lawrence (2018)⁴⁶. 36 focal sampling points are selected every 10 m, with 4 additional sampling points at 1 m N, S, E and W of each focal point. For 8 randomly selected focal points, 20 additional points are taken at 2 and 3 m N, S, E and W of each focal point. Finally, 25 additional points are selected within the grid at 5 m N, S, E and W from the focal points.

The method is used to study ground heterogeneity over length scales ranging from 1 to 80 m, and is well applicable to lidar ground-truthing. The closest in time lidar data available are from the 22nd February 2015. Contrary to the GPS data lidar may detect clumps of vegetation. Some marsh accretion may have occurred between February and August 2015; however, with accretion rates in Tollesbury MR of 2.3 cm/yr to 2.89 cm/yr as estimated by the present and a previous study⁴⁷, the difference is unlikely to be detected by RTK GPS or lidar. Most of the differences between the two datasets should thus be due to the vegetation present in February 2015.

Most data points (97% using the DSM and 96% using the DTM) fall within the limit of agreement (Fig. 5). A positive systematic bias of 0.09 m was found between the lidar and GPS data, for both the DSM and the DTM: lidar tends to overestimate the elevation values by 0.09 m. This value fits with the expected elevation of dense saltmarsh canopy in the UK⁴⁸. The persistence of this positive bias in the DTM shows that the EA vegetation
<figure><img src="image_10.png" /><figcaption>**Fig. 4** Layout of RTK-GPS points used for validation of lidar data. The RTK GPS points were collected in August 2015 in a natural marsh (NAT) and in the Tollesbury MR schemes (MR), overlain over lidar DSM collected in February 2015. RTK GPS dataset courtesy of Peter Lawrence.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 5** Pairwise difference between RTK GPS and lidar data in saltmarshes. The differences between lidar (a-c: DSM and d-f: DTM) and RTK-GPS data are investigated at a natural (NAT) and at an artificial (MR) marsh in Tollesbury. The sum of both locations is marked under TOT. The three red lines show the limits of agreement (2x standard deviation) and the mean value of the differences. The dashed green line shows the ideal mean difference if there is no bias between the two methods.</figcaption></figure>

removal algorithms fail to systematically remove the vegetation cover. This is expected when the vegetation cover is lower than the vertical resolution of lidar; however previous studies have found lidar DTM to record dense clumps of reeds as high as 2 m⁴⁹. The Root Mean Square Error of the lidar versus GPS elevation data is higher for the natural saltmarsh (0.13 m) than for the MR scheme (0.07 m), probably due to the higher ground heterogeneity⁴⁵ and higher vegetation diversity found in natural saltmarshes³⁰.
Finally, some fluctuation can be expected in the mean vegetation cover between summer and winter⁵⁰. This should be visible as fluctuations in the mean marsh elevation data from year to year. If interannual fluctuations in the mean marsh elevation data are small, and if the elevation changes are linear, then it can be assumed that the changes are due to accretion rather than vegetation growth and die-back. Elevation changes data as detected by lidar should also be confronted with field monitoring data of accretion rates when available in the literature.

**DSM/DTM comparison.** Compared to the uncorrected DSM, the EA proprietary DTM generation algorithms aim to remove first return points corresponding to vegetation or standing water, and to replace them with points of latter arrival corresponding to the ground. In order to explore the efficiency of those algorithms, the DTM was subtracted from the raw DSM in HOMW, were trees were present in 2009, for 5 different years. Most of the differences between the models occur along and inside the creek network, where the elevation is higher in the DSM than in the DTM. Some of the larger positive values correspond to trees, isolated or in a line along one of the main channels. According to the standard deviation values, discrepancies of about 0.26 m can be expected between the DSM and the DTM.

However, overall the amount of correction differs greatly between years: it is suspected that the algorithm itself has been updated over time, but that the newest algorithm has not been systematically applied to the older datasets, leading to inconsistencies. Some tiles in the 2007 dataset also seem to have undergone no correction at all, showing that the algorithms are not always systematically applied, in accordance with previous findings⁵¹.

Uncertainties concerning the functioning and consistency of application of EA algorithms make the DTM datasets unreliable for monitoring the morphological evolution of coastal wetlands. Furthermore, since the DTM datasets do not remove the systematic bias caused by low saltmarsh vegetation as seen previously, this study uses the raw DSM. This approach assumes that the low vegetation cover characteristic of saltmarsh areas is unlikely to mask the creek network or significantly affect the detection of creek edges, as based on aerial photography and field observations, plants develop on the creek banks but rarely within the creeks themselves. This choice might lead to underestimation of the channel depth, if the laser is reflected by residual water within the creek during low tide. It could also lead to overestimation of the saltmarsh elevation if the vegetation cover is detected as the ground level, and could limit the monitoring of accretion rates when performing the MR evolution analysis over several years⁵².

Some open-source algorithms attempt to correct vegetation from lidar data, but 1) they rely on precise knowledge of the vegetation distribution on site, obtained through regular field surveying, which defeats the purpose of lidar as a quick and cheap monitoring method over large areas, and 2) the correction factors based on local dominant vegetation height creates unrealistic “steps” in the dataset, which may complicate the detection of creek networks⁵³. Overall, vegetation removal has been a major issue in previous lidar-based saltmarsh monitoring studies⁵²: even state of the art lidar sensors fail to penetrate the saltmarsh canopy⁵⁴. This study provides an opportunity to estimate the efficiency of the uncorrected DSM at detecting evolving creek systems.

**Creek algorithm validation.** In order to verify its accuracy, the creek parametrisation algorithm was tested on 13 mature natural British saltmarshes. Creek morphological parameters were obtained independently from manual mapping and field surveying⁵⁵. We extracted creek parameters for those sites from recent lidar datasets (2014–2016) using the newly developed algorithm. Supposing that the natural saltmarshes are at equilibrium and that the effect of sea-level rise on marsh extent is negligible over the considered timescale, the algorithm should yield similar creek parameters. An example of comparison between the two creek mapping methods is given in Fig. 6, while the complete validation process is detailed in Chirol *et al.*²¹.

In order to assess the agreement between the two pairs of readings, detect outliers and visualise systematic biases, the differences were plotted on a Bland and Altman diagram (Fig. 7). Most of the differences between parameters lie within the limits of agreement of 2\* standard deviation. The errors of omission in creek number compared with Steel's results generally increase for the smaller creeks of higher Reverse Strahler Orders (Fig. 7a). No visible bias was found in the detection of mean creek length (Fig. 7b) and junction angle (Fig. 7d). Data from both methods were in good agreement for the creek numbers (negative mean difference of 7 creeks with the algorithm, Fig. 7a), the cross-sectional areas (mean negative difference of 1.1 m², Fig. 7f) and the bifurcation ratio (mean positive difference of 0.3, Fig. 7c). The depth of channels measured using lidar were ~ 0.4 m shallower than Steel's (1996) field validated results (Fig. 7h), probably due to the presence of residual water at the bottom of creeks: this is a limitation of using near infrared lidar data which cannot penetrate water⁵⁶. Furthremore, creek width is overestimated when adjacent creeks are detected as one channel due to the resolution of the dataset, leading to a mean difference of 2 m (Fig. 7g). Width overestimation and depth underestimation leads to a width/depth ratio overestimation of 6.3 compared to Steel's (1996) results (Fig. 7i). Even though the values fell within those expected of intertidal creek networks, between 5 and 34⁴³, no correlation could be found between the width/depth ratio and the reverse Stralher order. However, in the case of the cross-sectional area and the mean width value (given by the area/depth), the depth underestimation had a much lower impact, and the results are close to Steel's (1996)⁵⁵, with higher values for first reverse Strahler order channels, and a positive mean difference of 1.9 for the mean width/depth ratio (Fig. 7j).

Outliers in creek length for the entry channels and largest tributaries (first and second reverse Strahler orders) can generally be explained by the differences in detected creek branching. Indeed, the bifurcation ratio varies between ~1 and 2 depending on the creek extraction method chosen (Fig. 7c): those variations in detected branching change the distribution of creek orders, with a potential knock-on effect on the rest of the creek system (Fig. 8). A good example of this is the detected creek length at Longton marsh (Fig. 7b), where the first reverse Strahler order length is significantly underestimated while creek length is overestimated for the second order: subtle differences in the detection of high order creeks led to vastly different characteristics of the
<figure><img src="image_12.png" /><figcaption>**Fig. 6** Creek area and skeleton extraction result comparison for three saltmarshes using manual and semi-automated mapping. The creek networks in the natural saltmarshes Grange, Hen Hafod and Longton were manually mapped by Steel in 1996⁵⁵, then mapped semi-automatically using lidar data collected in 2014 and 2015. (a) Creek network skeleton manually extracted by Steel. (b) creek mask. (c) creek skeleton (b and c overlain over Steel's manual extraction results).</figcaption></figure>

entry channel. However, this problem does not affect the characteristics of the whole system such as the total channel length, drainage density, total creek volume and OPL.

Overall, there is a good agreement between the field-validated creek morphological parameters and those extracted by the algorithm, showing that this method applied to 1 m horizontal and 0.15 m vertical resolution is adequate to capture the evolution of the majority of active creeks. The chosen resolution for the lidar data also means that the likelihood of detecting ephemeral creeks is low. Most of the differences are linked to limits in ground detection by lidar, due to the presence of remnant water within creeks or of low-lying dense vegetation. Small differences in creek detection can have a knock-on effect on creek ordering, but the general creek
<figure><img src="image_13.png" /><figcaption>**Fig. 7** Morphological parameter differences for each reverse Strahler order (lidar extraction results minus Steel results⁵⁵). (a) number of creeks; (b) mean length; (c) bifurcation ratio; (d) junction angle; (e) sinuosity ratio; (f) cross-sectional area; (g) creek width; (h) creek depth; (i) width/depth ratio from top of creek; (j) area/depth² ratio (mean width/depth ratio). Red lines show the limits of agreement (2* standard deviation) surrounding the mean value of the differences. The dashed green line shows the ideal mean difference of 0 if there is no bias between the two methods.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Fig. 8** Illustration of how subtle differences in creek detection can significantly impact branching and creek order throughout the system (valid for both Strahler and Reverse Strahler ordering). Detection method **a** led to one 4th Reverse Strahler (RS) order creek to be omitted (red dashed circle) compared with method **b**, with a knock-on effect on the rest of the creek system, notably the length of the large RS order 1 creek. The detected entry channel (RS order 1) is twice as long and sinuous in **b** than in **a**.</figcaption></figure>

<table><caption>**Table 4.** Creek morphological parameters used in the study and data processing uncertainty mean calculated as the standard deviation of morphological parameters detected by the algorithm at Hesketh Out Marsh West (HOMW) when the elevation thresholds are changed by +/- 0.15 m.</caption><thead><tr><td>Parameter (tested on HOMW 2009–2014)</td><td>Symbol</td><td>Mean standard deviation (% of mean value)</td></tr></thead><tbody><tr><td>Mean elevation above MWS (m)</td><td>MWS</td><td>N/A</td></tr><tr><td>Drainage density (km/km²)</td><td>DD</td><td>0.5 (5%)</td></tr><tr><td>Overmarsh path length (m)</td><td>OPL</td><td>5.08 (11%)</td></tr><tr><td>Main channel length (m)</td><td>MCL</td><td>36.8 (3%)</td></tr><tr><td>Total channel length (m)</td><td>TCL</td><td>800 (5%)</td></tr><tr><td>Number of creeks (no unit)</td><td>NB</td><td>21.2 (9%)</td></tr><tr><td>Total mouth cross-sectional area (m²)</td><td>CSA</td><td>4.98 (4%)</td></tr><tr><td>Main channel mouth depth (m)</td><td>D</td><td>0.33 (12%)</td></tr><tr><td>Planform area (m²)</td><td>PA</td><td>3.49*10⁴ (20%)</td></tr><tr><td>Undermarsh tidal prism (creek volume) (m³)</td><td>TP</td><td>2.02*10⁴ (17%)</td></tr><tr><td>Sinuosity ratio (no unit)</td><td>SR</td><td>0.07 (6%)</td></tr><tr><td>Main channel gradient (°)</td><td>MCG</td><td>0.05 (4%)</td></tr></tbody></table>

characteristics such as total channel length, volume and OPL remain unaffected. The uncertainty linked to creek detection can be quantified using the standard deviation of morphological parameter values when the elevation thresholds vary by +/- 0.15 m. This uncertainty was tested for Hesketh Out Marsh West (HOMW), one of the largest and more complex creek networks considered, using 4 lidar datasets taken between 2009 and 2014 (Table 4).

## Code availability

The custom code used to produce the dataset, written in Matlab R2015a, is available within the data repository via https://eprints.soton.ac.uk/434946/ (see CHIROL_CREEK_EXTRACTION_DEMO zip folder and utilisation guide README_Creek_Extraction_Algorithm).

Readers can also refer a previously published methodology paper²¹ which details the algorithm's functioning and validation process.

Received: 21 November 2019; Accepted: 11 February 2022;

Published online: 01 April 2022

## References

1. Liu, L., Wang, Hjun & Yue, Q. China's coastal wetlands: Ecological challenges, restoration, and management suggestions. *Reg. Stud. Mar. Sci.* **37**, 101337 (2020).

2. Stewart-Sinclair, P. J. et al. Blue Restoration – Building Confidence and Overcoming Barriers. *Front. Mar. Sci.* **7** (2020).

3. Adam, P. *Salt marsh restoration. Coastal Wetlands: An Integrated Ecosystem Approach*. https://doi.org/10.1016/B978-0-444-63893-9.00023-X (2018).

4. Davidson, N. C. How much wetland has the world lost? Long-term and recent trends in global wetland area. *Mar. Freshw. Res.* **65**, 934–941 (2014).

5. Gu, J. et al. Losses of salt marsh in China: Trends, threats and management. *Estuar. Coast. Shelf Sci.* **214**, 98–109 (2018).
6. Alber, M., Swenson, E. M., Adamowicz, S. C. & Mendelssohn, I. A. Salt Marsh Dieback: An overview of recent events in the US. *Estuar. Coast. Shelf Sci.* **80**, 1–11 (2008).

7. Bobbink, R., Beltman, B., Verhoeven, J. T. A. & Whigham, D. F. Wetlands functioning in relation to biodiversity conservation and restoration. In *Wetlands: Functioning, Biodiversity Conservation, and Restoration* 1–12 (2006).

8. Gedan, K. B., Kirwan, M. L., Wolanski, E., Barbier, E. B. & Silliman, B. R. The present and future role of coastal wetland vegetation in protecting shorelines: Answering recent challenges to the paradigm. *Climatic Change* **106**, 7–29 (2011).

9. Luisetti, T. et al. Coastal Zone Ecosystem Services: From science to values and decision making; a case study. *Sci. Total Environ.* **493C**, 682–693 (2014).

10. Mudd, S. M. The life and death of salt marshes in response to anthropogenic disturbance of sediment supply. *Geology* **39**, 511–512 (2011).

11. Pontee, N. *Hydromorphological considerations for the restorations of intertidal habitats. Coastal Sediments Conference* 2015 (2015).

12. Morris, R. K. A. Management Managed realignment as a tool for compensatory habitat creation - A re-appraisal. *Ocean Coast. Manag.* **73**, 11 (2013).

13. Townend, I. H. Breach design for managed realignment sites. *Proc. ICE - Marit. Eng.* **161**, 9–21 (2008).

14. Sanderson, E. W., Ustin, S. L. & Foin, T. C. The influence of tidal channels on the distribution of salt marsh plant species in Petaluma Marsh, CA, USA. *Plant Ecol.* **146**, 29–41 (2000).

15. Morzaria-Luna, H., Callaway, J. C., Sullivan, G. & Zedler, J. B. Relationship between topographic heterogeneity and vegetation patterns in a Californian salt marsh. *J. Veg. Sci.* **14**, 523–530 (2004).

16. Zheng, Z., Zhou, Y., Tian, B. & Ding, X. The spatial relationship between salt marsh vegetation patterns, soil elevation and tidal channels using remote sensing at Chongming Dongtan Nature Reserve, China. *Acta Oceanol. Sin.* **35**, 26–34 (2016).

17. Crooks, S., Schutten, J., Sheern, G. D., Pye, K. & Davy, A. J. Drainage and elevation as factors in the restoration of salt marsh in Britain. *Restor. Ecol.* **10**, 591–602 (2002).

18. Wallace, K. J., Callaway, J. C. & Zedler, J. B. Evolution of tidal creek networks in a high sedimentation environment: A 5-year experiment at Tijuana Estuary, California. *Estuaries* **28**, 795–811 (2005).

19. Wu, Y. et al. The size and distribution of tidal creeks affects salt marsh restoration. *J. Environ. Manage.* **259**, 110070 (2020).

20. Faber, P. W. & Salt, P. Marsh Restoration Experience in San Francisco Bay. *J. Coast. Res.* **27**, 203–211 (2001).

21. Chirol, C., Haigh, I. D., Pontee, N., Thompson, C. E. & Gallop, S. L. Parametrizing tidal creek morphology in mature saltmarshes using semi-automated extraction from lidar. *Remote Sens. Environ.* **209**, 291–311 (2018).

22. Wu, F. et al. Disturbance is complicated: Headward-eroding saltmarsh creeks produce multiple responses and recovery trajectories. *Limnol. Oceanogr.* 1–15, https://doi.org/10.1002/lno.11867 (2021).

23. Shi, W., Shao, D., Gualtieri, C., Purnama, A. & Cui, B. Modelling long-distance floating seed dispersal in salt marsh tidal channels. *Ecohydrology* **13**, 1–13 (2020).

24. Taramelli, A., Valentini, E., Cornacchia, L., Monbaliu, J. & Sabbe, K. Indications of Dynamic Effects on Scaling Relationships Between Channel Sinuosity and Vegetation Patch Size Across a Salt Marsh Platform. *J. Geophys. Res. Earth Surf.* **123**, 2714–2731 (2018).

25. Ortals, C., Cordero, O., Valle-Levinson, A. & Angelini, C. Flows, Transport, and Effective Drag in Intertidal Salt Marsh Creeks. *J. Geophys. Res. Ocean.* **126**, 1–18 (2021).

26. Kiesel, J. et al. Effective design of managed realignment schemes can reduce coastal flood risks. *Estuar. Coast. Shelf Sci.* **242** (2020).

27. Tan, L. S. et al. The roles of vegetation, tide and sediment in the variability of carbon in the salt marsh dominated tidal creeks. *Estuar. Coast. Shelf Sci.* **239**, 106752 (2020).

28. Esteves, L. S. Is managed realignment a sustainable long-term coastal management approach? *J. Coast. Res.* 933–938, https://doi.org/10.2112/SI65-158.1 (2013).

29. Brady, A. F. & Boda, C. S. How do we know if managed realignment for coastal habitat compensation is successful? Insights from the implementation of the EU Birds and Habitats Directive in England. *Ocean Coast. Manag.* **143**, 164–174 (2017).

30. Mossman, H. L., Davy, A. J. & Grant, A. Does managed coastal realignment create saltmarshes with 'equivalent biological characteristics' to natural reference sites? *J. Appl. Ecol.* **49**, 1446–1456 (2012).

31. Lohani, B. & Mason, D. C. Application of airborne scanning laser altimetry to the study of tidal channel geomorphology. *ISPRS J. Photogramm. Remote Sens.* **56**, 100–120 (2001).

32. Mason, D. C., Scott, T. R. & Wang, H. Extraction of tidal channel networks from airborne scanning laser altimetry. *ISPRS J. Photogramm. Remote Sens.* **61**, 67–83 (2006).

33. Environment Agency. Defra Data Services Platform. (2019).

34. Rapinel, S., Clément, B., Nabucet, J. & Cudennec, C. Ditch network extraction and hydrogeomorphological characterization using LiDAR-derived DTM in wetlands. *Hydrol. Res.* **46**, 276–290 (2015).

35. Liu, X. Airborne LiDAR for DEM generation: some critical issues. *Prog. Phys. Geogr.* **32**, 31–49 (2008).

36. Vlaswinkel, B. M. & Cantelli, A. Geometric characteristics and evolution of a tidal channel network in experimental setting. *Earth Surf. Process. Landforms* **36**, 739–752 (2011).

37. Admiralty. Admiralty Tide Tables. Volume 1, 2014. Admiralty Charts and Publications, Chart NP201. United Kingdom Hydrographic Office, Taunton. (2014).

38. Haigh, I., Nicholls, R. & Wells, N. Mean sea level trends around the English Channel over the 20th century and their wider context. *Cont. Shelf Res.* **29**, 2083–2098 (2009).

39. Woodworth, P. L., Teferle, F. N., Bingley, R. M., Shennan, I. & Williams, S. D. P. Trends in UK mean sea level revisited. *Geophys. J. Int.* **176**, 19–30 (2009).

40. ABPmer. OMReg Habitat Creation Scheme Database - ABPmer. (2018).

41. Fagherazzi, S. et al. Tidal networks 1. Automatic network extraction and preliminary scaling features from digital terrain maps. *Water Resour. Res.* **35**, 3891–3904 (1999).

42. Sullivan, J. C., Torres, R. & Garrett, A. Intertidal Creeks and Overmarsh Circulation in a Small Salt Marsh Basin. *J. Geophys. Res. Earth Surf.* **124**, 447–463 (2019).

43. Zeff, M. L. Salt Marsh Tidal Channel Morphometry: Applications for Wetland Creation and Restoration. *Restor. Ecol.* **7**, 205–211 (1999).

44. Chirol, C., Haigh, I. D., Pontee, N., Thompson, C. E. & Gallop, S. L. Dataset: Morphological evolution of creek networks in 10 restored coastal wetlands in the UK. *University of Southampton* https://doi.org/10.5258/SOTON/D1114 (2019).

45. Brooks, K. L., Mossman, H. L., Chitty, J. L. & Grant, A. Limited Vegetation Development on a Created Salt Marsh Associated with Over-Consolidated Sediments and Lack of Topographic Heterogeneity. *Estuaries and Coasts* **38**, 325–336 (2015).

46. Lawrence, P. J. How to create a saltmarsh: Understanding the roles of topography, redox and nutrient dynamics A thesis submitted in partially fulfilment of the requirement of the Manchester Metropolitan University for the degree of Doctor of Philosophy Faculty of Scien. (2018).

47. Garbutt, R. A., Reading, C. J., Wolters, M., Gray, A. J. & Rothery, P. Monitoring the development of intertidal habitats on former agricultural land after the managed realignment of coastal defences at Tollesbury, Essex, UK. *Mar. Pollut. Bull.* **53**, 155–164 (2006).

48. Neumeier, U. Quantification of vertical density variations of salt-marsh vegetation. *Estuar. Coast. Shelf Sci.* **63**, 489–496 (2005).

49. Pontee, N. *The use of LiDAR to assess siltation within the Alkborough Scheme Humber Estuary Flood Risk Management Strategy.* **44** (2014).
50. Boorman, L. A. & Ashton, C. The productivity of salt marsh vegetation at Tollesbury, Essex, and Stiffkey, Norfolk, England. *Mangroves and Salt Marshes* 1, 113–126 (1997).

51. Leung, S. F. *Analysis of siltation rates in managed realignment schemes using lidar data.* (2017).

52. Wang, C. et al. Separation of ground and low vegetation signatures in LiDAR measurements of salt-marsh environments. *IEEE Trans. Geosci. Remote Sens.* 47, 2014–2023 (2009).

53. Hladik, C., Schalles, J. & Alber, M. Remote Sensing of Environment Salt marsh elevation and habitat mapping using hyperspectral and LIDAR data. *Remote Sens. Environ.* 139, 318–330 (2013).

54. Hladik, C. & Alber, M. Accuracy assessment and correction of a LIDAR-derived salt marsh digital elevation model. *Remote Sens. Environ.* 121, 224–235 (2012).

55. Steel, T. J. The morphology and development of representative British saltmarsh creek networks. (1996).

56. Brzank, A., Heipke, C., Goepfert, J. & Soergel, U. Aspects of generating precise digital terrain models in the Wadden Sea from lidar-water classification and structure line extraction. *ISPRS J. Photogramm. Remote Sens.* 63, 510–528 (2008).

57. Blott, S. J. & Pye, K. Application of lidar digital terrain modelling to predict intertidal habitat development at a managed retreat site: Abbotts Hall, Essex, UK. *Earth Surf. Process. Landforms* 29, 893–905 (2004).

58. Shepherd, D. et al. Modelling the effects and economics of managed realignment on the cycling and storage of nutrients, carbon and sediments in the Blackwater estuary UK. *Estuar. Coast. Shelf Sci.* 73, 355–367 (2007).

59. ABP mer. *Case Study on the Abbotts Hall Managed Realignment Scheme (England).* (2011).

60. Scott, C. R., Armstrong, S., Townend, I. H., Dixon, M. & Everard, M. Lessons learned from 20 years of managed realignment and regulated tidal exchange in the UK. 1–10, https://doi.org/10.1680/iczm.57494.365 (2011).

61. Pendle, M. & Wallingford, H. R. *Estuarine and coastal managed realignment sites in England selected case studies.* (2013).

62. Mason, W. L., McKay, H. M., Weatherall, A., Connolly, T. & Harrison, A. J. The effects of whole-tree harvesting on three sites in upland Britain on the growth of Sitka spruce over ten years. *Forestry* 85, 111–123 (2012).

63. Costa, S. S. A study of morphological changes within managed realignment sites along the Humber estuary (UK). (2013).

64. Medlock, A. J. M. & Vaux, A. G. C. Colonization of UK Coastal Realignment Sites by Mosquitoes: Implications for Design, Management, and Public Health Colonization of UK coastal realignment sites by mosquitoes: implications for design, management, and public health. *J. Vector Ecol.* 38, 53–62 (2013).

65. Pontee, N. Long-term sustainability of managed realignment. in *Esteves, L.S. Managed realignment: A viable long-term coastal management strategy?* 1–15 (2014).

66. ABP mer. *Case Study on the Wallasea (North) Managed Realignment Scheme (England).* (2010).

67. Dixon, M., Morris, R. K. A., Colclough, S., Birchenough, A. & Scott, C. R. Managed realignment–lessons from Wallasea, UK. *Proc. ICE - Marit. Eng.* 161, 61–71 (2008).

68. ABP mer. *Case Study on the Chowder Ness Managed Realignment Scheme (England).* (2011).

69. Friess, D. A. et al. Remote sensing of geomorphological and ecological change in response to saltmarsh managed realignment, The Wash, UK. *Int. J. Appl. Earth Obs. Geoinf.* 18, 57–68 (2012).

70. Friess, D. a. et al. Coastal saltmarsh managed realignment drives rapid breach inlet and external creek evolution, Freiston Shore (UK). *Geomorphology* 208, 22–33 (2014).

71. Frost, N., Hull, S. & Pontee, N. *Habitat quality measures and monitoring protocols. Technical report* (2005).

72. Hampshire, S. F. Evolution and design of tidal creeks within Managed Realignment sites: a study at Freiston Shore, Hampshire. (2011).

73. Nottage, A. S. & Robertson, P. A. *Salt Marsh Creation Handbook*. (Royal Society for the Protection of Birds, 2005).

74. Rotman, R., Naylor, L., McDonnell, R. & MacNiocail, C. Sediment transport on the Freiston Shore managed realignment site: An investigation using environmental magnetism. *Geomorphology* 100, 241–255 (2008).

75. Symonds, A. M. & Collins, M. B. The establishment and degeneration of a temporary creek system in response to managed coastal realignment: The Wash. *UK. Earth Surf. Process. Landforms* 32, 1783–1796 (2007).

76. Symonds, A. M. Impacts of coastal realignment on intertidal sediment dynamics: Freiston Shore, The Wash. (2006).

77. Tovey, E. L., Pontee, N. I. & Harvey, R. Managed Realignment at Hesketh Out Marsh West. *Eng. Sustain.* 162, 223–228 (2009).

78. Clapp, J. Managed realignment in the Humber estuary: factors influencing sedimentation. (2009).

79. Edwards, A. M. C. & Winn, P. S. J. The Humber Estuary, Eastern England: Strategic planning of flood defences and habitats. *Mar. Pollut. Bull.* 53, 165–174 (2006).

80. Mazik, K. et al. Managed realignment as compensation for the loss of intertidal mudflat: A short term solution to a long term problem? *Estuar. Coast. Shelf Sci.* 90, 11–20 (2010).

81. Pontee, N. I., Hull, S. C. & Moore, J. R. Banked realignment: a case study from the Humber Estuary. *UK. Eng. Sustain.* 159, 99–108 (2006).

82. Burgess, K., Pontee, N., Wilson, T., Lee, S. C. & Cox, R. *Steart Coastal management Project: Engineering Challenges in a Hyper-tidal Environment*. (2013).

83. Townsend, D. Sediment response to an artificial creek network. *Plymouth University Thesis* (2014).

84. Wright, A., Shipton, J., Carroll, B. & Armstrong, S. *Lessons for designing Managed Realignment sites along hyper tidal estuaries-a case study on The Bristol Port Company's Steart Habitat Creation Scheme*. ftp://cpc1-heme10-2-0-cust487.9-1.cable.virginm.net/disk1/share/Files/BeeWorkDocs/Site visits,conferences,etc/ICECoastalManagement2011/Draftpapers/94_draft.pdf (2011).

85. Atkinson, P. W. et al. Managed realignment in the UK - the first 5 years of colonization by birds. *Ibis* (Lond. 1859). 146, 101–110 (2004).

86. Garbutt, R., Reading, C. & Wolters, M. Monitoring the development of intertidal habitats on former agricultural land after the managed realignment of coastal defences at Tollesbury, Essex, UK. *Mar. Pollut. Bull.* 53, 155–164 (2006).

87. Paramor, Oa. L. & Hughes, R. G. Effects of the invertebrate infauna on early saltmarsh plant colonisation of managed realignment areas in south-east England. *Mar. Ecol. Prog. Ser.* 303, 61–71 (2005).

88. Reading, C. J. et al. *Managed realignment at Tollesbury: R&D Technical Report FD 1922/TR.* (2008).

89. Steel, T. J. & Pye, K. The development of saltmarsh tidal creek networks: evidence from the UK. In *Canadian Coastal Conference 267–280 (1997).*

90. Watts, C. W., Tolhurst, T. J., Black, K. S. & Whitmore, A. P. *In situ* measurements of erosion shear stress and geotechnical shear strength of the intertidal sediments of the experimental managed realignment scheme at Tollesbury, Essex, UK. *Estuar. Coast. Shelf Sci.* 58, 611–620 (2003).

91. ABP mer. *Case Study on the Welwick Managed Realignment Scheme (England).* (2011).

## Acknowledgements

Funding for this project came from NERC, internal University of Southampton funding sources, and Jacobs. We thank the Environment Agency for providing the freely available lidar data that made this research possible. We thank Dr Peter Lawrence for sharing the RTK GPS dataset used in the Technical Validation section. We thank the University of Southampton Institutional Research Repository ePrints Soton for hosting and preserving the datasets and the Channel Coastal Observatory for hosting a link to the dataset.
## Author contributions

Dr Clementine Chirol collected the data, developed the creek mapping algorithm, interpreted the results and wrote the data descriptor. Dr Ivan D. Haigh co-designed the research project and secured funding, helped with data interpretation as first PhD supervisor, provided codes for the creek mapping algorithm and reviewed the data descriptor. Dr Nigel Pontee co-designed the research project, secured funding, helped with data interpretation and reviewed the data descriptor. Dr Charlie Thompson helped with data interpretation and liaising with the Channel Coastal Observatory, and reviewed the data descriptor. Dr Shari Gallop co-designed the research project, helped secure funding, assisted with data interpretation and reviewed the data descriptor.

## Competing interests

The authors declare no competing interests

## Additional information

Correspondence and requests for materials should be addressed to C.C.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_15.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022