Front. Mar. Sci. Frontiers in Marine Science Front. Mar. Sci. 2296-7745 Frontiers Media S.A. 10.3389/fmars.2023.1198536 Marine Science Original Research Erosion and accretion of salt marsh in extremely shallow water stages ChenDezhi 1 2 TangJieping 3 XingFei 4 * ChengJun 5 LiMingliang 6 ZhangYiyi 7 ShiBenwei 4 ShiLianqiang 2 WangYa Ping 4 8 * 1Key Laboratory of Tropical Marine Ecosystem and Bioresource, Fourth Institute of Oceanography, Ministry of Natural Resources, Beihai, China 2Guangxi Key Laboratory of Beibu Gulf Marine Resources, Environment and Sustainable Development, Fourth Institute of Oceanography, Ministry of Natural Resources, Beihai, China 3School of Electronic and Information Engineering, Guangdong Ocean University, Zhanjiang, China 4State Key Laboratory of Estuarine and Coastal Research, East China Normal University, Shanghai, China 5Department of Environmental & Sustainability Sciences, Kean University, Union, NJ, United States 6State Department of Environmental Geology, Geological Survey of Jiangsu, Nanjing, China 7Research Department of Tidal Flat, Tidal Flat Research Center of SOA (Jiangsu), Nanjing, China 8School of Geography and Ocean Science, Nanjing University, Nanjing, China

Edited by: Yifei Zhao, Nanjing Normal University, China

Reviewed by: Aijun Wang, State Oceanic Administration, China; Guoxiang Wu, Ocean University of China, China

*Correspondence: Ya Ping Wang, ypwang@nju.edu.cn; Fei Xing, fxing@sklec.ecnu.edu.cn

10 05 2023 2023 10 1198536 01 04 2023 24 04 2023 Copyright © 2023 Chen, Tang, Xing, Cheng, Li, Zhang, Shi, Shi and Wang 2023 Chen, Tang, Xing, Cheng, Li, Zhang, Shi, Shi and Wang

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Salt marshes, which commonly exist on the upper tidal flat, provide a natural barrier against sea level rise and coastal storm. The extremely shallow water stages (water depth< 0.2 m), including the initial stage of flood tides and the last stage of ebb tides, can induce a significant impact on sediment dynamics of saltmarshes and associated tidal flats, despite lasting for only a short time (around 10 min), which has been less studied. In this study, two parallel field sites were established to quantify erosion-accretion processes and morphological changes during extremely shallow water stages in salt marshes within Doulonggang tidal flat along the Jiangsu coast. Our results revealed that obvious accretion occurred during extremely shallow water stages, with a total deposition amount of +33.8 mm in vegetated areas and +20.8 mm in unvegetated areas. In contrast, erosion dominated during deep water stages, with a total erosion amount of -22.3 mm at the vegetated site and -32.7 mm at the unvegetated site. The magnitude of bed-level change during extremely shallow water stages was 7~8 times greater than that during deep water stages, even though the duration of extremely shallow water stages was only about 14~15% of the entire tidal cycle. Furthermore, strong winds significantly impacted deposition during extremely shallow water stages compared to calm weather. During the strong wind period, the average bed level change rate reached +0.15 mm/min and +0.12 mm/min in the vegetated and unvegetated areas, respectively. This is significantly higher than the +0.05 mm/min and +0.01 mm/min during the calm weather period. These results reveal that extremely shallow water stages have substantial impacts on sedimentary processes, which are vital for the maintenance of tidal flat systems.

 erosion accretion tidal flat salt marsh extremely shallow water stages shear stress section-in-acceptance Coastal Ocean Processes Introduction

As an important ecological type of coastal wetland, salt marshes play a key role in providing tidal flat ecosystem services and coastal protection. Its formation and development are affected by a series of physical processes and biological factors, including climate, shoreline shape, waves, tides, sediment sources, sediment characteristics, sea level change, and plant cover (Jacobson and Jacobson, 1989; Allen, 1996; Cahoon et al., 1996; Allen, 2000; Davidson-Arnott et al., 2002). The sedimentary dynamic process is the main controlling factor of salt marsh evolution, which determines the evolution of salt marsh morphology and the response of vegetation to sediment transport (Fagherazzi and Mariotti, 2012; Fagherazzi et al., 2013; Ganju et al., 2015; Ganju et al., 2017).

The extremely shallow water stages, which are defined as stages of water depth (h) less than 0.2 m occurring at the beginning of flood tide and at the end of the ebb tide, exhibit different hydrodynamic characteristics compared to the deep water stages (h > 0.2 m) (Shi et al., 2019). When tides propagate landward to the shallow tidal flat, a velocity surge appears along with the tidal front due to the gentle slope and bottom friction, and this surge can last for several minutes (Xu et al., 1994; Gao, 2010; Nowacki and Ogston, 2013). The extremely shallow water stages play an important role in controlling bottom hydrodynamics and further affect sediment transport and geomorphic evolution. The short-term strong hydrodynamics lead to a large volume of sediment resuspension (Fagherazzi and Mariotti, 2012; Zhang et al., 2021), accompanied by strong erosion and accumulation fluxes(Gao, 2010). The maximum current velocity and suspended sediment concentration (SSC) occurred during extremely shallow water stages (Zhang et al., 2016a), causing large and rapid topography changes (Nowacki and Ogston, 2013; Shi et al., 2017a; Shi et al., 2019). For example, Zhang et al. (2016b) found that bed shear stress is large enough to resuspend and transport a large amount of sediment during extremely shallow water stages, which resulted in severe scouring during the flood stages. Shi et al. (2019) found that bed shear stress during the extremely shallow water stage was twice of that during the deep water stages during the flood tides, which resulted in extensive erosion. During the ebb stages, the shear stress during extremely shallow water stages was only half of that during the deep water stages, resulting in large accretion.

Due to the difficulties in measurements within extremely shallow water environments, few continuous field investigations of hydrodynamic and sediment transport on salt marshes had been conducted, and the impact of extremely shallow water stages in salt marshes on morphological evolution has often been less studied. Numerical models also neglect the process of water and sediment movement during extremely shallow water stages. For example, the Delft3d model uses a water depth of about 0.1 m as the critical value for dry and wet cell determination as default [WL| Delft Hydraulics, 2010 (Eds.)], resulting in high uncertainties of morphological changes in extremely shallow water environments. However, sediment movement and morphological changes are very active during extreme shallow water stages, particularly in the presence of extensive salt marshes. Therefore, it is of importance to investigate sediment dynamics under extremely shallow water conditions to improve the accuracy of prediction regarding the morphological change of tidal flat systems.

This paper aims to investigate the differences in hydrodynamic conditions, SSCs, and bed erosion–accretion processes between extremely shallow water stages and deep water stages in salt marshes. This is crucial for understanding the long-term morphological evolution mechanism of salt marshes in high-turbidity intertidal zones.

 Study area

The study area is located in the open intertidal zone on the Doulong coast of Jiangsu, China. The intertidal flat is situated between the abandoned Yellow River mouth and Changjiang (Yangtze) River estuary (Figure 1). It has an average width of 8-10 km and a slope of 0.05%, with the offshore Radial Sand Ridges acting as the primary source of sediments to the system (Fan et al., 2016). This area has experienced continuous growth of intertidal flats in the past several decades (Wang et al., 2012). The Doulong coast is characterized as meso-to-macrotidal with irregular semidiurnal tides and an average tidal range of 3.68 m. Tidal currents are southward-dominated during the flood phase and northward-dominated during the ebb phase. The winter monsoon mainly originates from the north and northeast, with an annual average speed of 4.21 m/s, while the summer monsoon mainly originates from the south and southeast, with an annual average speed of 2.76 m/s. The northward wave is prevalent, with 85% of the frequency of waves less than 1 m.

Location of the Doulonggang intertidal flat (A, B) enlarged view of the area (red box), (C) zoomed-in view showing the location of observation sites, (D) instruments deployed at SM site (ADV, OBS, EMCM, SBE26plus), (E) instruments deployed at TF site (ADV, OBS, EMCM, RBR-solo D).

The surficial sediments on the tidal flat become finer landward, dominated by silt and fine sand. The exotic plant Spartina alterniflora occupies the upper tidal flat (elevation above 0.4 m). The SM site was situated within the vegetated areas while the TF site was within the unvegetated areas. The distance between the two sites was around 22 m. No distinct tidal creeks were observed in the upper tidal flat.

 Methods Field data collection

The field observation was conducted from November 6th to November 15th, 2018, and all in situ measurements were synchronized at both sites. Acoustic Doppler Velocimeters (6.0 MHz vector current meter, Nortek AS, Norway) were deployed in a down-looking configuration with the transmit transducer to measure 3D velocity at 16 Hz for 256 sec every burst (4096 points per 5-min time series). The ADV probe was installed at a height of 20 cm above the seabed, and the intratidal bed-level changes were indicated by the distance between the probe and the bed surface recorded by the ADV. The accuracy of bed-level changes was ±1 mm based on laboratory tests in previous studies (Hosseini et al., 2006; Shi et al., 2015; Shi et al., 2017b). To capture the hydrodynamics during extremely shallow water stages, an electromagnetic current meter (EMCM) was used because this instrument can precisely measure velocity within a small volume using a small probe (diameter of ~3 cm). To measure the 2D near-bed current velocity for extremely shallow water conditions, the EMCM was deployed at a height of 0.05 m above the bed and operated at a burst period of 30 s and sampling frequency of 2 Hz (Figures 1D, E).

Wave parameters such as wave height and wave period were measured using self-logging sensors, including the SBE-26plus Seagauge (Sea-Bird Electronics, Washington, USA) and RBRsolo |wave (RBR Ltd., Ottawa, Canada). The instruments were horizontally placed on the sediment surface, with their pressure sensors located 5 cm above the bed surface (Figures 1D, E). To obtain high-frequency water level measurements, the pressure data were recorded at 4 Hz for a duration of 256 seconds, resulting in a total of 1024 measurements per burst.

The optical backscatter (OBS) sensors (OBS-3A, D&A Instrument Company, Washington, USA) were used to measure turbidity in the water column every five minutes. The probe of OBS was positioned 5 cm above the bed surface. The suspended sediment concentration (SSC) values were obtained by converting the turbidity values measured by the OBSs through calibration using in situ water samples (Figure 2).

Calibration curves used to convert optical turbidity (T; NTU) recorded by OBSs to suspended sediment concentration (SSC) (C; g/L). (A) calibration curves for the TF site; (B) calibration curves for the SM site.

 Data processing and calculation Calculation of bed shear stress

The bed shear stress is a critical parameter controlling the erosion and deposition of sediments. The total bed shear stress due to waves and currents ( τ c w , N / m 2 ) is derived from the components of bed shear stress due to waves alone ( τ w , N / m 2 ) and currents alone ( τ c , N / m 2 ).

The shear stress generated by waves is closely related to the wave orbital velocity ( U w , m/s). When waves propagate onto the tidal flat, they are transformed due to local topography. The wave orbital velocity is calculated using the high-frequency velocity recorded by the ADV (Wiberg and Sherwood, 2008).

U w = 2 ∑ i S u v , i Δ f

where S u v = S u u + S v v is the combined horizontal spectrum of eastern and northern velocity. This method removes the influence of turbulence on wave parameters, thereby eliminating the overestimation of orbital velocity that would be induced by using the linear wave theory (Xiong et al., 2018).

However, the ADV probes were exposed to air and stopped working during extremely shallow water stages when the water depth was<0.2 m. Therefore, we used the SBE-26plus Seagage and RBRsolo | wave to observe waves during shallow water stages applying the linear wave theory, expressed as follows:

U w = π H T s i n h ( 2 π h / L )

where H ( m ) is the wave height, T ( s ) is the wave period, L ( m ) ( = ( g T 2 2 π ) tanh ( 2 π h L ) ) is the wavelength, g ( = 9.8   m / s 2 ) is the gravity acceleration, and h ( m ) is water depth.

The bed shear stress induced by waves ( τ w ) can be calculated from wave parameters as (Soulsby, 1995):

τ w = 1 2 ρ f w U w 2

where ρ is the fluid density (= 1030 kg/m3), f w is the wave friction factor which depends on the hydraulic regime. f w can be calculated as follows (Soulsby, 1997):

f w = { 2 R w − 0.5 ,           R w ≤ 10 5   ( laminar flow ) 0.0521 R w − 0.187 ,   R w > 10 5   ( smooth turbulent ) 0.237 r − 0.52 ,             ( rough turbulent )

where   R w   ( = U w A ν ) is the wave Reynolds number, r   ( = A k s ) is the relative roughness, v ( m / s ) is the kinematic viscosity of seawater, A (= U W T 2 π , m ) is the semi-orbital excursion, and k s (= 2.5 d 50 ,   d 50   is the median grain size) is the Nikuradse roughness (Fredsøe, 1984).

The bed shear stress induced by current ( τ c )   is calculated from friction velocity ( u * ) (Soulsby, 1995):

τ c = ρ w   u * 2

The logarithmic velocity profile (LP method) is used to describe the velocity structure in the bottom boundary layer within a weak hydrodynamic environment (Soulsby and Dyer, 1981; Grant and Madsen, 1986; Shi et al., 2019), expressed as (Dyer, 1986; Whitehouse et al., 1999)

u = u * k ln ( z z 0 )

where u ( m / s )   is the velocity at height z ( m ) above the bed, measured by EMCM during extremely shallow water stages and ADV during deep water stages. K ( = 0.4 ) is the von Karman’s constant, and Z 0 ( = k s 30 , m ) is the bed roughness length related to the Nikuradse grain roughness (Whitehouse et al., 1999).

The bed shear stress due to the combined action of waves and currents ( τ c w ) is determined using the proposed models by Soulsby and Clarke (2005), expressed as:

τ c w = ( τ m + τ w | c o s φ c w | ) 2 + ( τ w | s i n φ c w | ) 2

where φ c w   is the angle between waves and currents, and the average total shear stress   τ m is calculated as:

τ m = τ c [ 1 + 1.2 ( τ w τ c + τ w ) 3.2 ] Cumulative bed level changes during extremely shallow water stages

The cumulative bed level changes of extremely shallow water stages for two consecutive tidal cycles were obtained using the following method:

Δ D i _ shallow  =   D i _ e n d   –   D i + 1   _ i n i t i a l

where D i e n d ( m m ) is the relative distance between the ADV probe and the bed surface at the last effective burst during the ebb tide stage, D i + 1   _ i n i t i a l is the distance between the ADV probe and the bed surface for the first effective burst of the next tide. A positive value indicates deposition and a negative value indicates erosion.

 Critical shear stress for erosion ( <inline-formula> <mml:math display="inline" id="im37"> <mml:mrow> <mml:msub> <mml:mi>τ</mml:mi> <mml:mrow> <mml:mi>c</mml:mi> <mml:mi>e</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>)

The critical shear stress for erosion ( τ c e ) was estimated following the approach recommended by Andersen et al. (2007) and Shi et al. (2015), which has been validated in tidal flat environments. According to their theory, the values of τ c w recorded at the point of erosion initiation (i.e., when there was an noticeable decrease in bed level elevation and a significant increase in SSC) are considered as the critical shear stress for erosion ( τ c e ).

The bed level started to drop during the flooding tide of T3 when τ c w reached a high value of 0.24 N / m 2 (Figures 3D, F) at the TF site, whereas at the SM site, erosion started when τ c w was about 0.17 N / m 2 during the flooding tide of T9 (Figures 3D, F), as determined from the bed-level changes obtained using the high-resolution ADV. These values (0.24 and 0.17 N / m 2 ) are considered as the τ c e values for the TF and SM sites, respectively.

Time series of (A) water depth, (B) bed shear stress due to currents ( τ c ); (C) bed shear stress due to waves ( τ w ); (D) bed shear stress due to currents and waves ( τ c w ); (E) suspended sediment concentration (SSC); (F) bed level changes (positive value represents deposition, the negative value represents erosion, and the single point represents the value of each burst.

 Index of agreement

The Index of agreement ( I n d e x ) is often used to evaluate the difference between the same physical quantity from two different methods by showing the degree of similarity quantitatively:

I n d e x =   1 − ∑ ​ ( x − y ¯ ) 2 ∑ ​ ( | x − y ¯ | + | y − y ¯ | ) 2

Where x and y are the two datasets involved in the comparison, and the I n d e x ranges from 0 to 1. The larger the value of the I n d e x , the higher the degree of similarity between the two datasets. I n d e x =1 indicates that the two datasets are completely consistent. In this study, the index of agreement is used to compare the wave orbital velocity calculated using different methods.

 Results Hydrodynamics

During our field measurements, wind patterns showed significant temporal variations at both sites, which were classified into two periods. The first period, composed of T1~T3 and T9~T11, was characterized by strong northerly onshore winds with speeds ranging from 3.4 to 10.5 m/s and averaging at 7.6 m/s (Figure 4A). These relatively strong onshore winds generated large waves, and the average significant wave height reached a maximum of 0.43 m and 0.33 m at the two sites. The average significant wave heights during extremely shallow water stages were 0.03 m and 0.05 m at the SM and TF sites, respectively (Figure 4C). The tide-average maximum significant wave height was 0.27 m at the SM site, 25% lower than that at the TF site (0.36 m) due to the damping effect of vegetation.

Time series of (A) wind speed and direction. Wind data were recorded every hour; (B) water depth; (C) significant wave height; (D) current velocity at the TF, and the SM sites.

The second period, composed of T4~T8 and T12~T16, was characterized by relatively weak winds ranging from 0.85~8.6 m/s, with an average of 4.6 m/s. During this period, the tide-averaged maximum wave heights were much lower, with an average of 0.12 m and 0.19 m at the SM and TF sites, respectively. Since wave height was limited by water depth in tidal flat environments (Figure 3A), the average significant wave heights during extremely shallow water stages were 0.02 m and 0.03 m at the SM and TF sites, respectively (Figure 4C).

The tidal current near the bottom (at the height of 0.05 m above the bed) at both sites was characterized by a rotational flow pattern (Figure 4D). High current velocities normally occurred during the flood stages of extremely shallow water stages, followed by low velocities during deep water stages, and the minimum current speed occurred during slack water. At the SM site, the tide-averaged current velocity near the bed during extremely shallow water stages ranged from 0.01 to 0.22 m/s (average value of 0.05 m/s), which was much higher than that during deep water stages (0.01–0.09 m/s; the average value of 0.02 m/s). Similarly, the velocities at the TF site during extremely shallow water stages (ranging from 0.01 to 0.29 m/s with an average of 0.06 m/s) were higher than those during deep water stages (ranging from 0.01 to 0.23 m/s with an average of 0.05 m/s).

 Erosion and accretion processes Shear stress

The wave orbital velocity calculated by the two different methods were shown in Figure 5. Linear wave theory was used to calculate U w using the data recorded by the SBE26plus and RBRsolo |wave, while the spectral method was used to estimate U w using ADV recorded data. For both sites, the time series of U w obtained using different methods show similar variation patterns during tidal cycles. At the SM site, the two methods provide close U w values with an I n d e x value reaching 0.91(Table 1). However, at the TF site, agreement in the magnitude of U w is weaker. The Uw values calculated from the two methods were more deviant during deep water stages than during shallow water stages ( I n d e x = 0.75), demonstrating that the two methods produce more consistent results in shallow water environments. For convenience, the spectral method is used in the calculation of U w , except during extremely shallow water stages when ADV stopped working and the linear wave theory is used.

Time-series of (A) wave orbital velocity calculated by different methods at the SM site (B) wave orbital velocity calculated by different methods at the TF site.

Index of agreement for wave orbital velocity using different methods.

Site Method Linear Spectrum
SM Linear 1.00 0.91
Spectrum 0.91 1.00
TF Linear 1.00 0.75
Spectrum 0.75 1.00

The values of τ c , τ w and τ c w are shown in Figures 3B–D and statistically summarized in Tables 2, 3. Generally speaking, for all tidal cycles, the average values of τ c during flood tides (0.06 N / m 2 ) in extremely shallow water stages were higher than those in ebb tides (0.02 N / m 2 ) and those in deep water stages because of the flood surge (Zhang et al., 2021) at the SM site. The differences in τ c were more pronounced on the unvegetated mudflat. The average value during shallow water stages in flood tides (0.18 N / m 2 ) was six times higher than that in ebb tides (0.03 N / m 2 ) and slightly higher than that in deep water stages (0.15 N / m 2 ). The τ w values during extremely shallow water stages were very small due to the shallow water environment. The average values during deep water stages (0.10 N / m 2 at TF site, 0.07 N / m 2 at SM site) were two times higher than those during extremely shallow water stages (0.05 N / m 2 at TF site, 0.03 N / m 2 at SM site).

Comparison of bed shear stress due to waves ( τ w ), currents ( τ c ), and combined current–wave action ( τ c w ), and SSC at water depths (h) of<0.2m (Extremely Shallow Water Stages) and >0.2m (deep water stages) at the SM site.

Tides τ c on average (N/m2) τ w on average (N/m2) τ c w on average (N/m2) S S C on average (kg/m3)
Extremely shallow water stages Deep water stages Extremely shallow water stages Deep water stages Extremely shallow water stages Deep water stages Extremely shallow water stages Deep water stages
Flood Ebb AVG Flood Ebb AVG Flood Ebb AVG Flood Ebb AVG
T1 0.09 0.02 0.04 0.11 0.02 0.04 0.03 0.11 0.10 0.05 0.07 0.21 0.30 0.11 0.16 1.68
T2 0.08 0.01 0.05 0.09 0.01 0.01 0.01 0.10 0.09 0.02 0.06 0.19 0.64 0.05 0.30 1.40
T3 0.10 0.03 0.06 0.08 0.03 0.04 0.03 0.12 0.09 0.07 0.08 0.20 0.88 0.10 0.40 1.92
T4 0.11 0.02 0.06 0.04 0.03 0.03 0.03 0.07 0.14 0.05 0.09 0.12 0.49 0.07 0.19 0.86
T5 0.11 0.01 0.06 0.04 0.05 0.03 0.04 0.05 0.14 0.04 0.09 0.09 0.69 0.24 0.29 1.24
T6 0.03 0.01 0.01 0.02 0.03 0.03 0.03 0.02 0.05 0.05 0.05 0.03 0.28 0.14 0.17 2.31
T7 0.02 0.02 0.02 0.01 0.03 0.03 0.03 0.03 0.05 0.05 0.05 0.04 0.83 0.08 0.26 0.91
T8 0.02 0.03 0.02 0.03 0.03 0.03 0.03 0.02 0.05 0.06 0.05 0.04 0.18 0.06 0.11 0.52
T9 0.07 0.05 0.05 0.04 0.03 0.03 0.03 0.08 0.10 0.08 0.08 0.12 0.80 0.14 0.20 2.81
T10 0.06 0.05 0.05 0.05 0.07 0.02 0.05 0.08 0.13 0.07 0.10 0.13 0.46 0.06 0.10 0.87
T11 0.01 0.02 0.02 0.10 0.05 0.04 0.05 0.14 0.07 0.07 0.07 0.23 1.65 0.06 0.28 2.03
T12 0.04 0.03 0.03 0.04 0.03 0.04 0.03 0.04 0.08 0.07 0.07 0.08 0.36 0.12 0.13 0.51
T13 0.06 0.04 0.04 0.05 0.01 0.04 0.03 0.09 0.07 0.08 0.06 0.15 0.73 0.16 0.25 1.07
T14 0.02 0.01 0.01 0.05 0.03 0.02 0.03 0.01 0.05 0.04 0.04 0.05 0.05 0.05 0.05 0.15
T15 0.03 0.02 0.03 0.04 0.04 0.03 0.04 0.07 0.08 0.05 0.06 0.12 0.96 0.59 0.64 3.06
T16 0.04 0.01 0.02 0.03 0.03 0.02 0.02 0.03 0.06 0.03 0.04 0.06 0.24 0.07 0.14 0.84
AVG 0.06 0.02 0.04 0.05 0.03 0.03 0.03 0.07 0.08 0.06 0.07 0.12 0.60 0.13 0.23 1.39

Comparison of bed shear stress due to waves ( τ w ), currents ( τ c ), and combined current–wave action ( τ c w ), and SSC at water depths (h) of<0.2m (Extremely Shallow Water Stages) and >0.2m (deep water stages) at the TF site.

Tides τ c on average (N/m2) τ w on average (N/m2) τ c w on average (N/m2) S S C on average (kg/m3)
Extremely shallow water stages Deep water stages Extremely shallow water stages Deep water stages Extremely shallow water stages Deep water stages Extremely shallow water stages Deep water stages
Flood Ebb AVG Flood Ebb AVG Flood Ebb AVG Flood Ebb AVG
T1 0.18 0.04 0.09 0.14 0.13 0.08 0.10 0.15 0.27 0.11 0.15 0.29 0.53 0.25 0.30 3.07
T2 0.21 0.02 0.09 0.11 0.03 0.08 0.06 0.15 0.23 0.11 0.16 0.26 2.11 0.21 0.75 3.65
T3 0.20 0.01 0.11 0.14 0.03 0.05 0.04 0.16 0.22 0.07 0.15 0.31 2.13 0.24 0.66 3.59
T4 0.14 0.02 0.06 0.13 0.09 0.07 0.08 0.11 0.22 0.10 0.14 0.25 1.23 0.64 0.78 3.23
T5 0.21 0.02 0.12 0.14 0.05 0.07 0.06 0.09 0.26 0.09 0.18 0.22 1.27 0.96 1.03 2.85
T6 0.12 0.02 0.06 0.05 0.02 0.04 0.03 0.03 0.14 0.06 0.09 0.08 1.89 1.08 1.14 3.58
T7 0.15 0.04 0.09 0.10 0.06 0.06 0.06 0.04 0.21 0.09 0.15 0.14 2.01 0.67 1.02 1.26
T8 0.15 0.04 0.09 0.27 0.02 0.03 0.02 0.02 0.16 0.07 0.11 0.28 0.74 0.74 0.74 1.55
T9 0.22 0.06 0.12 0.11 0.02 0.07 0.05 0.11 0.23 0.12 0.17 0.21 2.52 0.55 1.14 5.82
T10 0.19 0.03 0.11 0.16 0.08 0.02 0.05 0.11 0.27 0.05 0.16 0.27 0.56 0.56 0.56 4.87
T11 0.25 0.03 0.10 0.15 0.03 0.07 0.05 0.18 0.27 0.10 0.15 0.32 2.05 0.24 0.64 5.09
T12 0.23 0.03 0.11 0.17 0.02 0.04 0.03 0.07 0.25 0.07 0.14 0.23 1.26 0.31 0.67 1.60
T13 0.23 0.01 0.10 0.11 0.05 0.03 0.03 0.14 0.27 0.04 0.13 0.24 1.39 0.56 0.93 2.66
T14 0.10 0.02 0.04 0.25 0.05 0.08 0.07 0.02 0.12 0.09 0.11 0.27 0.84 0.13 0.47 0.77
T15 0.23 0.02 0.10 0.12 0.02 0.04 0.03 0.11 0.25 0.07 0.13 0.23 2.07 1.78 1.82 5.34
T16 0.11 0.03 0.06 0.20 0.06 0.04 0.05 0.06 0.17 0.06 0.11 0.25 0.82 0.73 0.77 1.95
AVG 0.18 0.03 0.09 0.15 0.05 0.05 0.05 0.10 0.22 0.08 0.14 0.24 1.46 0.60 0.84 3.18

The maximum τ c w within a tidal cycle usually appeared during deep water stages (Figure 3D), which was much larger than that during extremely shallow water stages. Correspondingly, the average τ c w during the flood stages (0.22 N / m 2 at TF site, 0.08 N / m 2 at SM site) and ebb stages (0.08 N / m 2 at TF site, 0.06 N / m 2 at SM site) during extremely shallow water stages was much lower than the τ c e values (0.24 N / m 2 at TF site, 0.17 N / m 2 at SM site). Therefore, effective resuspension cannot be initiated.

 Suspended sediment concentration

The SSC near the bed had a similar temporal variation pattern to τ c w (Figures 3D, E), with peak values usually occurring at the highest water level. The average SSC during flood stages was greater than that during ebb stages. The average SSC during deep water stages (3.18 kg / m 3 at the TF site, 1.39 kg / m 3 at the SM site) was much higher than that during extremely shallow water stages at both sites (Figure 3E; Tables 2, 3).

 Duration of extremely shallow water stages and deep water stages

Due to continuous emergence and submergence, salt marshes experience temporal-spatial variations in water depth, leading to the frequent occurrence of extremely shallow water conditions (twice every tide). Compared with the relatively long duration of the deep water stages (224 min at the TF site, 233 min at the SM site), the duration of extremely shallow water tides during the flood tides (10 min at the TF site, 9 min at the SM site) was only half of that during the corresponding ebb stages (19 min at the TF site, 23 min at the SM site), and the total duration of extremely shallow water stages accounted for 14~15% of the entire tidal cycle (Tables 4, 5).

Comparison of duration of deep (> 0.2 m) and extremely shallow water stages (< 0.2 m) within different tidal cycles and the rates of bed-level change (mm/min) for different stages at the TF site.

Tide Entire tidal cycle(min) deep water stages extremely shallow water stages
Duration(min) % Rate (mm/min) Duration(min) % Rate (mm/min)
Flood Ebb Flood Ebb Flood Ebb Total Flood Ebb Total* Flood Ebb
T1 327 124 167 38% 51% -0.05 -0.01 -0.02 7 29 33 2% 9% 0.11
T2 289 120 153 42% 53% 0.04 -0.03 0.00 5 11 14 1% 4% 0.36
T3 291 121 152 42% 53% -0.15 -0.03 -0.08 5 13 22 1% 4% 0.14
T4 193 62 102 32% 53% -0.03 -0.01 -0.01 9 20 31 5% 10% 0.02
T5 210 74 108 35% 51% 0.01 0.00 0.00 11 17 26 5% 8% -0.02
T6 145 37 82 26% 57% 0.04 0.02 0.02 9 17 31 6% 12% 0.00
T7 238 107 102 45% 43% -0.02 0.02 0.00 14 15 29 6% 6% -0.01
T8 159 57 67 36% 42% 0.05 -0.03 0.01 14 21 30 9% 13% 0.03
T9 293 112 152 38% 52% -0.04 0.01 -0.01 9 20 29 3% 7% 0.00
T10 244.5 77 145 31% 59% 0.04 -0.02 0.00 9 14 21 4% 6% 0.15
T11 288.5 122 147 42% 51% -0.02 0.01 0.00 7 13 22 2% 4% 0.13
T12 153 52 77 34% 50% -0.01 0.01 0.00 9 15 27 6% 10% -0.03
T13 235 95 112 40% 48% -0.03 0.02 0.00 12 16 27 5% 7% 0.09
T14 88.5 17 32 19% 36% 0.00 -0.01 0.00 11 29 38 12% 32% -0.02
T15 220.5 102 87 46% 39% -0.01 0.00 0.00 9 23 40 4% 10% 0.05
T16 209 95 72 45% 34% 0.01 -0.02 0.00 17 25 / 8% 12% /
AVG 224 86 110 37% 48% -0.01 0.00 -0.01 10 19 28 5% 10% 0.07

*The total duration of extremely shallow water stages was the ebb duration of the previous tide plus the flood duration of the next tide during extremely shallow water stages.

Comparison of duration of deep (> 0.2 m) and extremely shallow water stages (< 0.2 m) within different tidal cycles and the rates of bed-level change (mm/min) for different stages at the SM site.

Tide Entire tidal cycle(min) deep water stages extremely shallow water stages
Duration(min) % Rate (mm/min) Duration(min) % Rate (mm/min)
Flood Ebb Flood Ebb Flood Ebb Total Flood Ebb Total* Flood Ebb
T1 328 125 167 38% 51% -0.06 0.02 -0.01 5 31 37 2% 9% 0.13
T2 282 122 138 43% 49% -0.04 0.00 -0.02 6 16 22 2% 6% 0.30
T3 282 122 138 43% 49% -0.05 0.00 -0.02 6 16 23 2% 6% 0.13
T4 211 70 108 33% 51% -0.01 0.03 0.01 7 26 37 3% 12% 0.06
T5 227 75 108 33% 48% 0.00 0.01 0.00 11 33 43 5% 15% 0.00
T6 169 45 90 27% 53% 0.03 -0.02 0.00 10 24 37 6% 14% 0.24
T7 231 107 90 46% 39% -0.03 0.02 -0.01 13 21 33 6% 9% 0.00
T8 183 65 82 35% 45% 0.01 0.00 0.00 12 24 32 7% 13% 0.00
T9 276 117 136 42% 49% -0.04 0.01 -0.01 8 15 23 3% 5% -0.07
T10 260 82 150 32% 58% -0.02 0.00 -0.01 8 20 27 3% 8% 0.05
T11 286 125 136 44% 48% -0.01 -0.01 -0.01 7 18 26 2% 6% 0.17
T12 191 57 98 30% 51% -0.04 0.01 -0.01 8 28 35 4% 15% -0.02
T13 237 100 109 42% 46% -0.03 0.02 0.00 7 21 32 3% 9% 0.08
T14 120 23 42 19% 35% -0.07 0.00 -0.02 11 24 32 9% 20% 0.03
T15 230 107 94 47% 41% -0.02 0.02 0.00 8 21 31 3% 9% 0.06
T16 213 97 82 45% 39% 0.01 0.01 0.01 10 24 / 5% 11% /
AVG 233 90 111 38% 47% -0.02 0.01 -0.01 9 23 31 4% 10% 0.08

* The total duration of extremely shallow water stages was the ebb duration of the previous tide plus the flood duration of the next tide during extremely shallow water stages

 Bed level changes

The data of bed level changes measured by ADV during the period with effective data (h > 0.2 m) are shown in Figure 3F. Previous studies have shown that the high sediment concentration water layer near the seabed and fluid mud can affect the quality of data obtained from ADV measurements (Sahin et al., 2012; Mehta et al., 2014). During our observational period, fluid mud was not detected, and the high SSCs only occurred during deep water stages. The low-quality data were manually excluded to obtain high-quality realistic bed-level changes. The result showed that the measured bed level changes decreased during T1~T3, which was defined as the ‘erosion phase’ and then increased slowly, defined as the ‘accretion phase’ for the back-siltation during T4~T16 at the TF site. At the SM site, deposition dominated during T1–T8, and erosion occurred continuously from T9 to T16. The cumulative vertical bed-level change was +11.5 mm at the SM site and -13.5 mm at the TF site at the end of T16.

Accumulation occurred during extremely shallow water stages with a total deposition amount of +33.8 mm at the SM site and +20.8 mm at the TF site with an average bed-level change rate of 0.08 mm/min at the SM site and 0.07 mm/min at the TF site. On the other hand, erosion dominated the two sites during deep water stages, with a total erosion thickness of -22.3 mm at the SM site and -32.7 mm at the TF site with an average bed-level change rate of -0.01 mm/min at both sites.

Wind also had a significant impact on determining the deposition rate during extremely shallow water stages. Strong wind induced high waves and faster current speeds, moving more sediment to the upper flat zone. Our observation showed that the tide cycle-averaged deposition during the strong wind period was + 3.1 mm at the SM site and + 2.9 mm at the TF site, respectively. This was several times higher than that during the calm weather period (+ 1.7 mm at the SM site and + 0.4 mm at the TF site, respectively) during extremely shallow water stages. The average bed-level deposition rate reached +0.15 mm/min at the SM site and +0.12 mm/min at the TF site during the strong wind period, much higher than that during the calm weather period (+0.05 mm/min at the SM site and +0.01 mm/min at the TF site).

 Discussion Geomorphological changes during extremely shallow water stages in salt marshes

Tidal current and corresponding sediment transport patterns have shaped the unique geomorphic features of the mudflat. In previous studies, surges at the beginning of the flood period were as large as tens of centimeters per second (Gao, 2010; Shi et al., 2019; Zhang et al., 2021), leading to large values of bed shear stress (reaching 1.5 N/m2 at the beginning of the flood tide) (Zhang et al., 2016a), which is an order of magnitude higher than the critical shear stress found on intertidal flats. Such a phenomenon also occurred in our study area. A monitoring site (CT02) was set seawards in the intertidal zone, which was 600 m away from the SM and TF site. When the tidal flat was exposed to the air, a ruler was used to measure the distance between the ADV and seabed (the red points in Figure 6E) to reflect the erosion and accretion results during the extremely shallow water stage at the CT02 site. As shown in Figure 6, surges were commonly observed during extremely shallow water stages. Affected by the high flow velocity (0.2 m/s), the shear stress greatly increased, with an average value of 0.42 N/m2 during the flood tide. The short-term strong hydrodynamics led to a large volume of bottom sediment resuspension, and the SSC reached the first peak, accompanied by strong erosion fluxes with an average magnitude of -21.9 mm per tide. These results were similar to those of previous studies (Xu et al., 1994; Gao, 2010; Hughes, 2012; Zhang et al., 2016a; Shi et al., 2019).

Time series of hydrodynamics and sediment transport-related parameters at site CT02 (A) water depth and significant wave height, (B) current velocity, (C) suspended sediment concentration (SSC), (D) bed shear stress due to currents ( τ c ), waves ( τ w ), and current and wave interactions ( τ c w ); (E) bed level changes (the first point in T1 was set to zero), blue bars indicate the flood tide during extremely shallow water stages, and green bars indicate the ebb tide during extremely shallow water stages.

However, the mechanism by which the surge was formed and how it influences sediment transport in the salt marsh was still uncovered, where hydrodynamic characters were quite different from those on the mudflat. The existence of vegetation increases surface roughness, reduces flow velocity, and influences the generation of turbulence and energy dissipation, weakening the hydrodynamic force within the salt marsh (Leonard and Luther, 1995; Leonard et al., 2002; Neumeier and Ciavola, 2004; Neumeier and Amos, 2006). In this study, the flood tide was relatively slow when it propagated to the salt marsh, and the flow velocity was only 0.05 m/s in the vegetated areas and 0.06 m/s in the unvegetated areas during extremely shallow water stages. Breaking water was not observed at those sites. Additionally, limited by the shallow water depth, waves were very small, resulting in very weak hydrodynamic force during extremely shallow water stages. As a result, the mean τ c w value ranged from 0.04~0.1 N / m 2 in the vegetated areas and 0.09~0.18 N / m 2 in neighboring unvegetated areas, much lower than the critical bed shear stress (0.24 N / m 2 in the unvegetated areas and 0.17 N / m 2 in neighboring vegetated areas). Therefore, the tidal front water was too weak to generate effective resuspension, and the local sediments in the salt marsh were not stirred at this stage. Meanwhile, the tidal front water with high SSC, caused by sediment erosion of the lower tidal flat, was carried to the upper flat along the tidal flat profile (Gao, 2010). In this area, the low-sloped morphology and vegetation lead to high bed resistance, causing the settling of suspended sediment and a decrease in suspended sediment concentration during this period. These sediments deposited in the salt marsh with a mean value of 2.3 mm per tide and 1.4 mm per tide in the unvegetated areas during extremely shallow water stages (Table 6), and the water column was relatively clear for a short period. Compared with the unvegetated areas, the vegetated areas had a higher flow reduction rate and sediment capture effectiveness (Leonard and Croft, 2006; Yang and Nepf, 2019; Chen et al., 2020), resulting in more deposition. The cumulative net deposition thickness during extremely shallow water stages of the whole tidal cycle in the vegetated areas was + 33.8 mm and + 20.8 mm in the unvegetated areas.

Statistics of the bed level changes (unit: mm) at the TF and the SM sites.

Tide SM TF
Deep water stages Extremely shallow water stages Deep water stages Extremely shallow water stages
1 -3 4.9 -8.1 3.5
2 -5.1 6.5 -0.4 5
3 -6.5 2.9 -22.6 3
4 2.1 2.2 -2.7 0.6
5 0.8 -0.1 0.6 -0.4
6 -0.3 8.7 3.3 0
7 -1.6 0.1 -0.3 -0.2
8 0.8 0.1 1.3 0.8
9 -2.8 -1.6 -2.7 0
10 -1.4 1.4 0.2 3
11 -2.7 4.3 -1.1 2.7
12 -1.4 -0.8 0.2 -0.8
13 -0.1 2.4 -0.2 2.3
14 -1.9 0.9 -0.2 -0.8
15 -0.8 1.9 -0.8 2.1
16 1.6 / 0.8 /
Total -22.3 33.8 -32.7 20.8

After the extremely shallow water stage, flooding tides and waves became stronger, and τ w greatly increased to a mean value of 0.07 N / m 2 (vegetated areas) and 0.10 N / m 2 (unvegetated areas) during the deep water stages of different tidal cycles (Tables 2, 3). The mean τ c w also increased to higher values ranging between 0.05 to 0.23 N / m 2 (vegetated areas) and 0.09~0.32 N / m 2 (unvegetated areas), exceeding the local critical bed shear stress and causing strong sediment resuspension. Accordingly, SSC increased (Figure 3E), and sediment was transported away by tidal currents. Erosion was more obvious in the unvegetated areas because of the relatively strong hydrodynamic force and lack of vegetation protection. Our results showed that the net erosion amount during deep water stages was -32.7 mm in the unvegetated areas and -22.3 mm in the vegetated areas. The distinct sediment dynamics between the two sites resulted in a net deposition in the vegetated areas and net erosion in the unvegetated areas, thus influencing the morphological evolution of the tidal flat system.

Although the duration of extremely shallow water stages in the salt marsh was only a few minutes, accounting for about 14-15% of the whole tidal cycle, the rate of bed level changes during these stages was 7~8 times greater than that during deep water stages at both sites. This resulted in a profound impact on sediment transport, highlighting the importance of extremely shallow water stages in the replenishment of sediments and the maintenance of salt marshes.

Strong winds can significantly contribute to the development of salt marshes (French and Spencer, 1993; Schuerch et al., 2013), according to previous studies. Storms can generate strong waves in a short period of time and directly affect the front of the tidal flat. This process is crucial for understanding the long-term morphological evolution mechanism of salt marshes in high-turbidity intertidal zone. Abundant sediments are brought into the salt marsh from the subtidal zone under the influence of strong onshore winds during extremely shallow water stages. Vegetation causes sediment adhesion and deposition (Shi et al., 2000; Leonard and Croft, 2006), and the average SSC increased from 1.35 g/L (during calm weather) to 1.65 g/L (during rough weather) in unvegetated areas, while it increased from 0.48 g/L to 0.79 g/L in the vegetated areas during extremely shallow water stages in the initial flood stages (Figure 3E). Consequently, the average bed level changes rate reached +0.15 mm/min in the vegetated areas and +0.12 mm/min in the unvegetated areas, much higher than that during calm weather. Our results indicate that weather conditions are also crucial factors in determining sediment transport patterns within salt marshes during extremely shallow water stages.

 Linking extremely shallow water stages with microtopography in salt marshes

Many studies have found that the formation of ripples is the result of local turbulence occurring at the interface between the seabed and water (Coleman and Melville, 1994; Bartholdy et al., 2005; Bose and Dey, 2012). Previous field observations showed a high erosion rate and strong turbulence lead to an increase in ripple height, which can cause the formation of microtopography, such as sand ripples and grooves (Zhang et al., 2016a), and further promote the formation of large geomorphic units, such as tidal creek (Zhou et al., 2014). During our observations, sand ripples approximately 1-2 cm in height and 5-10 cm in length were found on the tidal flat. The wave crest can be flattened by the sheet flow (Gao, 2010) and evolved into a ‘flat bed’, which may occur during extremely shallow water stages.

In this study, the salt marsh was observed to be dominated by strong accumulation during extremely shallow water stages, resulting in a lack of ripples at the two sites. Abundant accumulation filled the wave trough, and the relatively strong deposition rate and weak turbulence led to a reduction in ripple height (Shi et al., 2017a). Therefore, the topography was relatively flat in the observed salt marsh area (Figure 7). Additionally, salt marshes reduced the resuspension rate of sediment by decreasing flow velocity and stabilizing sediment. The surface sediment within 10 cm was re-mixed by the burrowing organisms within 4 to 6 hours, thereby also flattening current-formed ripple patterns. Overall, the hydrodynamic processes within salt marshes during extremely shallow water stages limit the generation of micro topography.

Photographs of ripples at the (A) SM site and (B) TF site.

As discussed above, strong deposition occurred during extremely shallow water stages, while erosion occurred during deep water stages in salt marshes, resulting in significant bed-level changes. We conclude that although extremely shallow water stages are transient, this period substantially influences sediment transport within tidal cycles, and plays a significant role in the formation and evolution of microtopography on tidal flats.

 Conclusions

Extremely shallow water stages had an important impact on the maintenance of salt marshes by inducing significant sediment accretion in vegetated areas. This study explored the morphological dynamics and quantified the sediment erosion–accretion processes during extremely shallow water stages through integrated near-bed field measurements within salt marshes. The main conclusions are as follows:

(1) Current and wave induced shear stress was less than the critical shear stress for erosion ( τ c e ) during extremely shallow water stages but larger than τ c e during deep water stages. These changes led to deposition during extremely shallow water stages and erosion during deep water stages in almost all observed tidal cycles. In addition, gradual accumulation processes played a very important role in replenishing sediments in salt marshes.

(2) The deposition rate was more significant under windy weather during extremely shallow water stages. The average bed level changes rate reached +0.15 mm/min in the vegetated areas and +0.12 mm/min in the unvegetated areas, much higher than that during calm weather (+0.05 mm/min in the vegetated areas and +0.01 mm/min in unvegetated areas).

(3) Although extremely shallow water stages made up only 14~15% of the entire tidal cycle, the tide-average bed level change rate during extremely shallow water stages was 7 times faster than that during deep water stages.

(4) The deposition during extremely shallow water stages filled the wave trough, resulting in the relatively flat topography in salt marsh and restricting the generation of microtopography within salt marshes.

 Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

 Author contributions

The contributions made by each of the authors are listed as follows: 1) YW & FX put forward the idea, designed the experiments and funded the study. 2) DC processed the main measurements/experiments data and completed the major sections of the manuscript. 3) JT helped processing data. 4) JC, ML,YZ & BS reviewed this article and made suggestions to improve it. All authors contributed to the article and approved the submitted version.

 Acknowledgments

This research was supported by the following grants: the Jiangsu Special Program for Science and Technology Innovation (JSZRHYKJ202106), the Youth Foundation of Guangxi Zhuang Autonomous Region (2022GXNSFBA035566), innovation Driven Development Foundation of Guangxi (AD22080035), the Young Scientists Fund of the National Natural Science Foundation of China (42006149), the Key Laboratory of Coastal Salt Marsh Ecosystems and Resources, Ministry of Natural Resources (KLCSMERMNR2021108). We thank Gao C, Pan YP, Lan TF, Lu T for their assistance in the field work and laboratory measurements.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 References Allen J. R. L. (1996). Geomorphology and sedimentology of estuaries. Sedimentary Geol. 105, 110–111. doi: 10.1016/0037-0738(95)00149-2 Allen J. R. L. (2000). Morphodynamics of Holocene salt marshes: a review sketch from the Atlantic and southern north Sea coasts of Europe. Quaternary Sci. Rev. 18, 1155–1231. doi: 10.1016/S0277-3791(99)00034-7 Andersen T. J. Fredsoe J. Pejrup M. (2007). In situ estimation of erosion and deposition thresholds by acoustic Doppler velocimeter (ADV). Estuarine Coast. Shelf Sci. 75, 327–336. doi: 10.1016/j.ecss.2007.04.039 Bartholdy J. Flemming B. W. Bartholomä A. Ernstsen V. B. (2005). Flow and grain size control of depth-independent simple subaqueous dunes. J. Geophysical Research: Earth Surface 110 (F4). doi: 10.1029/2004JF000183 Bose S. K. Dey S. (2012). Instability theory of sand ripples formed by turbulent shear flows. J. Hydraulic Eng. 138, 752–756. doi: 10.1061/(ASCE)HY.1943-7900.0000524 Cahoon D. Lynch J. Powell A. (1996). Marsh vertical accretion in a southern California estuary, U.S.A. Estuar. Coast. Shelf Sci. - Estuar. Coast. SHELF Sci. 43, 19–32. doi: 10.1006/ecss.1996.0055 Chen D. Li M. Zhang Y. Zhang L. Tang J. Wu H. . (2020). Effects of diatoms on erosion and accretion processes in saltmarsh inferred from field observations of hydrodynamic and sedimentary processes. Ecohydrology 13, 2246. doi: 10.1002/eco.2246 Coleman S. E. Melville B. W. (1994). Bed-form development. J. Hydraulic Eng. 120, 544–560. doi: 10.1061/(ASCE)0733-9429(1994)120:5(544) Davidson-Arnott R. G. D. van Proosdij D. Ollerhead J. Schostak L. (2002). Hydrodynamics and sedimentation in salt marshes: examples from a macrotidal marsh, bay of fundy. Geomorphology 48, 209–231. doi: 10.1016/S0169-555X(02)00182-4 Dyer K. R. (1986). Coastal and estuarine sediment dynamics. Chichester, Sussex (UK): John Wiley Sons, 1986, 358. Fagherazzi S. Mariotti G. (2012). Mudflat runnels: evidence and importance of very shallow flows in intertidal morphodynamics. Geophysical Res. Lett. 39 (14). doi: 10.1029/2012GL052542 Fagherazzi S. Mariotti G. Wiberg P. L. McGlathery K. J. (2013). Marsh collapse does not require Sea level rise. Oceanography 26, 70–77. doi: 10.5670/oceanog.2013.47 Fan X. Tao J. Zeng Z. Coco G. Zhang C. (2016). Mechanisms underlying the regional morphological differences between the northern and southern radial sand ridges along the jiangsu coast, China. Mar. Geol. 371, 1–17. doi: 10.1016/j.margeo.2015.10.019 Fredsøe J. (1984). Turbulent boundary layer in wave-current motion. J. Hydraulic Eng. 110, 1103–1120. doi: 10.1061/(ASCE)0733-9429(1984)110:8(1103) French J. R. Spencer T. (1993). Dynamics of sedimentation in a tide-dominated backbarrier salt marsh, Norfolk, UK. Mar. Geol. 110, 315–331. doi: 10.1016/0025-3227(93)90091-9 Ganju N. K. Defne Z. Kirwan M. L. Fagherazzi S. D'Alpaos A. Carniello L. (2017). Spatially integrative metrics reveal hidden vulnerability of microtidal salt marshes. Nat. Commun. 8, 14156. doi: 10.1038/ncomms14156 Ganju N. K. Kroeger K. D. Kirwan M. L. Dickhudt P. J. Guntenspergen G. R. Cahoon D. R. (2015). Sediment transport based metrics of wetland stability. Geophysical Res. Lett. 42, 7992–8000. doi: 10.1002/2015GL065980 Gao S. (2010). Extremely shallow water benthic boundary layer processes and the resultant sedimentological and morphological characteristic(in Chinese with English abstract). Acta Sedimentologica Sin. 28, 926–932. doi: 10.14027/j.cnki.cjxb.2010.05.005 Grant W. D. Madsen O. S. (1986). The continental-shelf bottom boundary layer. Annu. Rev. Fluid Mechanics 18, 265–305. doi: 10.1146/annurev.fl.18.010186.001405 Hosseini S. A. Shamsai A. Ataie-Ashtiani B. (2006). Synchronous measurements of the velocity and concentration in low density turbidity currents using an acoustic Doppler velocimeter. Flow Measurement Instrumentation 17, 59–68. doi: 10.1016/j.flowmeasinst.2005.05.002 Hughes Z. J. (2012). Tidal channels on tidal flats and marshes. In: Davis R. Jr Dalrymple R. (eds) Principles of Tidal Sedimentology. 269–300. doi: 10.1007/978-94-007-0123-6_11 Jacobson H. A. Jacobson G. L. Jr. (1989). Variability of vegetation in tidal marshes of Maine, U.S.A. Can. Joumal Bot. 67 (1), 230–238. doi: 10.1139/b89-03 Leonard L. A. Croft A. L. (2006). The effect of standing biomass on flow velocity and turbulence in spartina alterniflora canopies. Estuar. Coast. Shelf Sci. 69, 325–336. doi: 10.1016/j.ecss.2006.05.004 Leonard L. A. Luther M. E. (1995). Flow hydrodynamics in tidal marsh canopies. Limnology Oceanogr. 40, 1474–1484. doi: 10.4319/lo.1995.40.8.1474 Leonard L. A. Wren P. A. Beavers R. L. (2002). Flow dynamics and sedimentation in spartina alterniflora and phragmites australis marshes of the Chesapeake bay. Wetlands 22, 415–424. doi: 10.1672/0277-5212(2002)022[0415:FDASIS]2.0.CO;2 Mehta A. J. Samsami F. Khare Y. P. Sahin C. (2014). Fluid mud properties in nautical depth estimation. J. Waterway Port Coastal Ocean Eng. 140, 210–222. doi: 10.1061/(ASCE)WW.1943-5460.0000228 Neumeier U. R. S. Amos C. L. (2006). The influence of vegetation on turbulence and flow velocities in European salt-marshes. Sedimentology 53, 259–277. doi: 10.1111/j.1365-3091.2006.00772.x Neumeier U. Ciavola P. (2004). Flow resistance and associated sedimentary processes in a spartina maritima salt-marsh. J. Coast. Res. 20, 435–447. doi: 10.2112/1551-5036(2004)020[0435:FRAASP]2.0.CO;2 Nowacki D. J. Ogston A. S. (2013). Water and sediment transport of channel-flat systems in a mesotidal mudflat: willapa bay, Washington. Continental Shelf Res. 60, S111–S124. doi: 10.1016/j.csr.2012.07.019 Sahin C. Safak I. Sheremet A. Mehta A. J. (2012). Observations on cohesive bed reworking by waves: atchafalaya shelf, Louisiana. J. Geophysical Research: Oceans 117 (C9). doi: 10.1029/2011JC007821 Schuerch M. Vafeidis A. Slawig T. Temmerman S. (2013). Modeling the influence of changing storm patterns on the ability of a salt marsh to keep pace with sea level. J. Geophysical Res. Earth Surface 118, 84–96. doi: 10.1029/2012JF002471 Shi B. Cooper J. R. Li J. Yang Y. Yang S. L. Luo F. . (2019). Hydrodynamics, erosion and accretion of intertidal mudflats in extremely shallow waters. J. Hydrology 573, 31–39. doi: 10.1016/j.jhydrol.2019.03.065 Shi B. Cooper J. R. Pratolongo P. D. Gao S. Bouma T. J. Li G. . (2017a). Erosion and accretion on a mudflat: the importance of very shallow-water effects. J. Geophysical Res. Oceans 122 (12), 9476–9499. doi: 10.1002/2016JC012316 Shi Z. Hamilton L. J. Wolanski E. (2000). Near-bed currents and suspended sediment transport in salt marsh canopies. J. Coast. Res. 16, 909–914. Available at: http://www.jstor.org/stable/4300101. Shi B. Wang Y. P. Yang Y. Li M. Li P. Ni W. . (2015). Determination of critical shear stresses for erosion and deposition based on In situ measurements of currents and waves over an intertidal mudflat. J. Coast. Res. 31, 1344–1356. doi: 10.2112/JCOASTRES-D-14-00239.1 Shi B. W. Yang S. L. Wang Y. P. Li G. C. Li M. L. Li P. . (2017b). Role of wind in erosion-accretion cycles on an estuarine mudflat. J. Geophysical Res. Oceans 122 (1), 193–206. doi: 10.1002/2016JC011902 Soulsby R. L. (1995). Bed shear-stresses due to combined waves and currents. Adv. Coast. Morphodynamics, 4–20. Soulsby R. (1997). Dynamics of marine sands: a manual for practical applications. Oceanographic Lit. Rev. 9 (44), 947. Soulsby R. Clarke S. (2005). Bed shear-stresses under combined waves and currents on smooth and rough beds. Technical Report. HR Wallingford, UK. Soulsby R. L. Dyer K. R. (1981). The form of the near-bed velocity profile in a tidally accelerating flow. J. Geophys. Res. 86, 8067–8074. doi: 10.1029/JC086iC09p08067 Wang Y. P. Gao S. Jia J. Thompson C. E. L. Gao J. Yang Y. (2012). Sediment transport over an accretional intertidal flat with influences of reclamation, jiangsu coast, China. Mar. Geol. 291-294, 147–161. doi: 10.1016/j.margeo.2011.01.004 Whitehouse R. J. S. Soulsby R. L. Roberts W. (1999). Dynamics of estuarine muds: a manual for practical applications. Thomas Telford, London UK. Wiberg P. L. Sherwood C. R. (2008). Calculating wave-generated bottom orbital velocities from surface-wave parameters. Comput. Geosciences 34, 1243–1262. doi: 10.1016/j.cageo.2008.02.010 WL| Delft Hydraulics . (2010). User manual Delft3D-FLOW (Delft, The Netherlands: Deltares). Xiong J. Wang Y. P. Gao S. Du J. Yang Y. Tang J. . (2018). On estimation of coastal wave parameters and wave-induced shear stresses. Limnology Oceanogr.: Methods 16, 594–606. doi: 10.1002/lom3.10271 Xu Y. Wang B. Zhang k. (1994). On flood surge over muddy tidal flat along the shanghai coast and its mechanism(in chinese). Geographical Res. 13, 60–68. doi: 10.11821/yj1994030007 Yang J. Q. Nepf H. M. (2019). Impact of vegetation on bed load transport rate and bedform characteristics. Water Resour. Res. 55, 6109–6124. doi: 10.1029/2018WR024404 Zhang Q. Gong Z. Zhang C. Lacy J. Jaffe B. Xu B. . (2021). The role of surges during periods of very shallow water on sediment transport over tidal flats. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.599799 Zhang Q. Gong Z. Zhang C. Townend I. Jin C. Li H. (2016a). Velocity and sediment surge: what do we see at times of very shallow water on intertidal mudflats? Continental Shelf Res. 113, 10–20. doi: 10.1016/j.csr.2015.12.003 Zhang Q. Gong Z. Zhang C. Zhou Z. Townend I. (2016b). Hydraulic and sediment dynamics at times of very shallow water on intertidal mudflats: the contribution of waves. J. Coast. Res. 75, 507–511. doi: 10.2112/SI75-102.1 Zhou Z. Olabarrieta M. Stefanon L. D'Alpaos A. Carniello L. Coco G. (2014). A comparative study of physical and numerical modeling of tidal network ontogeny. J. Geophysical Research: Earth Surface 119, 892–912. doi: 10.1002/2014JF003092