Rivers and estuaries are flanked by floodplains built by mud and vegetation. Floodplains affect channel dynamics and the overall system's pattern through apparent cohesion in the channel banks and through filling of accommodation space and hydraulic resistance. For rivers, effects of mud, vegetation and the combination are thought to stabilise the banks and narrow the channel. However, the thinness of estuarine floodplain, comprised of salt marsh and mudflats, compared to channel depth raises questions about the possible effects of floodplain as constraints on estuary dimensions. To test these effects, we created three estuaries in a tidal flume: one with recruitment events of two live vegetation species, one with mud and a control with neither. Both vegetation and mud reduced channel migration and bank erosion and stabilised channels and bars. Effects of vegetation include local flow velocity reduction and concentration of flow into the channels, while flow velocities remained higher over mudflats. On the other hand, the lower reach of the muddy estuary showed more reduced channel migration than the vegetated estuary. The main system-wide effect of mudflats and salt marsh is to reduce the tidal prism over time from upstream to downstream. The landward reach of the estuary narrows and fills progressively, particularly for the muddy estuary, which effectively shortens the tidally influenced reach and also reduces the tidal energy in the seaward reach and mouth area. As such, estuaries with sufficient sediment supply are limited in size by tidal prism reduction through floodplain formation.
Introduction
The size and shape of natural estuaries are potentially modified by life, as they are often flanked by salt marsh and mudflats
. In rivers, floodplains formed by vegetation development and mud deposition affect the channel dimensions and channel–bar
pattern for review, see. Often, floodplain destruction by migrating and displacing channels is balanced by floodplain
formation. It is plausible that floodplains in estuaries, comprised of salt marsh and mudflats, have similar effects. However, there is in principle
no limit to the available flow discharge from the sea, unlike in rivers where long-term discharge frequency and magnitude are determined by hinterland
characteristics and climate. This raises the question of what limits the widening of estuaries, particularly how vegetation and mud sedimentation affect
the channel pattern in estuaries as far as such estuaries are not laterally constrained by valley walls. The objective is to understand the different
mechanisms by which vegetation or mud affect local channel dynamics and system-scale effects on the channel–bar pattern. Two alternative hypotheses
for effects of dynamic floodplain on the dimensions of estuaries are proposed here from observations and mechanisms in rivers.
The first hypothesis relates to the conservation of tidal energy in the landward direction. In some systems, floodplain formation along estuaries may
have led to relatively deep channels allowing further landward tidal penetration than in shallow estuaries. Estuaries that formed as part of
tidally dominated deltas such as the Mekong, Mahakam, Yangtze and Rhine developed narrow, deep channels surrounded by high floodplains
, as did some late Holocene ingressive estuaries . The channel depth and
convergence possibly cause the tides to propagate far landwards with backwater effects up to the delta apex.
The alternative, opposing, hypothesis relates to landward tidal energy reduction. The reduction of intertidal area by accreting floodplains with
vegetation reduces the tidal energy, especially on the fluvial–tidal transition, and pushes the tidal limit in the seaward direction. Indeed, on
longer timescales, floodplain formation along estuaries was concurrent with a reduction in tidal discharge, landward tidal penetration and overall
filling that reduced tidal system extent and depth over some portion of the Holocene . Modelling showed
that mudflat sedimentation alone is a sufficient condition for net infilling , while absence of fine sediment is known
to lead to drowning or keeps estuaries unfilled .
The expansion or filling of estuaries depends not only on mudflat and salt marsh formation but also on the morphodynamic interactions between river
discharge and tides and the channels and bars. While significant river discharge causes net seaward sediment transport
, sea-level rise, subsidence and shipping fairway deepening cause flood dominance that promotes sediment import
. Reconstructions of Holocene coastal plains show that salt marsh and mangrove development likely
amplified vertical accretion in agreement with numerical modelling
. The direct link of tidal prism and tidal asymmetry with the area occupied by bars and
floodplain is potentially another mechanism through which salt marsh and mudflat expansion affect the dimensions of estuarine systems. The underlying
cause for the reducing planform estuary dimensions may be the reduction of tidal prism by
sediment filling of intertidal storage area and vegetation causing hydraulic resistance and amplifying vertical accretion. In other words,
hypothetically, there is a direct effect of mud and vegetation by filling and an additional, indirect effect on the tidal dynamics that may determine
net sediment import or export.
System-scale effects of estuarine floodplain formation are relevant in a great number off settings. Intertidal mudflats and supratidal marshes,
riparian forest, mires and swamps in the estuarine landscape are important ecosystems along the river continuum and in the coastal zone
e.g.. Vegetation protects dikes against wave attack and enhanced sedimentation may
eventually provide fertile agriculture soils and may protect urban areas against flooding . The channelised parts function as
shipping fairways to provide access to the world's major ports . Understanding floodplain-forming processes is important in view
of the projected sea-level rise acceleration, ongoing attempts to build a natural elevated buffer around estuaries and the creation of intertidal
basins along estuaries to reduce flood levels. Whether this is a viable pathway depends on the effects of such floodplain dynamics on the large-scale
estuarine flow and morphodynamics e.g..
Effects of mud or vegetation in isolation on the entire estuary cannot be inferred well from observations in nature, because mud and vegetation often
settle in the same locations and because estuaries with either one or the other do not exist. On the other hand, numerical models are useful tools to
study their combined effects in rivers and estuaries e.g.. However, such models are also sensitive to choices
of a host of parameterisations of flow resistance, sediment transport , mud characteristics and vegetation
characteristics . Landscape experiments, as presented here, may complement numerical modelling by real physical and biological
processes at the cost of some known scale effects . Recent advances in tidal experiments allow for the formation of estuarine
systems with similarities in planform pattern of channel and bar morphology and mudflats and the long-term development as well as response to dredging
and sediment disposal . These experiments focussed on partially filled, multi-channel estuaries with
mid-channel bars that have also been studied in data and models
. Here, we also focus on this type of estuary (examples in
Fig. ).
Methods
Three estuaries were formed in experiments in the Metronome, a flume of 20 m by 3 m and 0.4 m deep that tilts periodically to
drive reversing tidal currents . The tilting causes similar sediment mobility of coarse sand to nature
. The tidal forcing was kept simple with one primary tidal component only. In all experiments, we supplied a minor river
discharge and applied monochromatic waves by a paddle in front of the ebb delta. The first experiment had only sand, while the second and third had
live vegetation or mud to compare the two styles of floodplain formation against a control without floodplain. The experiment with live vegetation
reported here is novel and is based on vegetation experiments in and . Low-density sediment was supplied as a
mud simulant; this experiment is also reported in . Furthermore, to avoid practical difficulties of flow measurements, we apply a
numerical flow model to all experiments , modifying the roughness where vegetation is present. In this section, the
experimental methods, vegetation treatment and numerical modelling methods are described.
Experimental methods
Experiments were conducted in the Metronome, which tilts with a sinusoidal motion over the short middle axis with a slope amplitude of
0.0068 mm-1 at a period of 40 s. This drives tidal currents with a depth-averaged velocity amplitude in the channels of about
0.3 ms-1.
The initial bed was screed flat at millimetre accuracy. The initial convergent channel of 0.2–1.0 m width and 0.03 m depth was carved
into a sand bed of 0.07 m thick, while the sea was kept free of sand between 17.8–20 m. The initial water depth above the channel
bed, set by the seaward weir, was 0.025 m. We used poorly sorted sand with a D10 of 0.33 mm, a D50 of 0.57 mm and a
D90 of 1.2 mm. Flume settings were chosen such that sediment mobility, expressed as the Shields number, was up to 0.5
.
The landward boundary condition is a river discharge of 0.1 Ls-1 supplied for half the tidal cycle when the flume was tilted seaward. The
river discharge in isolation did not mobilise the sediment but was essential to maintain a long estuary see online supplement
in. The seaward boundary condition is a broad-crested weir with constant water supplied from a basin to have a constant head
condition. The weir moves up and down in the opposite phase of the tilting motion such that the sea remains approximately horizontal. This prevents
surges into the estuary and drawdown that would cause incision in the ebb delta and ensures that the tidal flow is entirely forced by the
periodically varying gradient. The weir amplitude was reduced linearly with the progradation of the delta towards the boundary. Furthermore,
monochromatic waves of 0.5 s period and 0.01 m height were generated by a horizontal paddle during the landward tilting half of the
tidal cycle. Waves in isolation hardly mobilised the sediment but in combination with tidal currents caused coastal transport diffusion to round the
delta during progradation and form beach ridges. Testing and analyses of scaling of waves is described in the supplement to .
Bathymetry was measured every 1000 cycles for slight deviations in timing, see on the dry bed by stereophotography and structure
from motion by Agisoft PhotoScan, constrained by control points at millimetre accuracy also see. Vegetation and mud
were recognised in the overhead imagery through colour also see.
Floodplain formation
One experiment had regular vegetation recruitment events through periodic supply of vegetation seeds of two species followed by 4 d of rest to
allow sprouting above the still water surface. The seeding was started after 4500 cycles to allow time for formation of bars where vegetation could
settle and was done every 2000 cycles. Vegetation was added by seeding at the river boundary, similarly to the hydrochorous seed distribution in the
river experiment of for scaling of vegetation, also see. After growth tests, we selected two species
that grow above the water surface to simulate effects of salt marsh species by enhanced flow resistance and capturing of suspended sediment
. Veronica beccabunga forms 10 mm tall plants that grow at the waterline in dense elongated patches and above the
waterline in sparse cover, and Lotus pedunculatus grows up to 20 mm tall at somewhat higher elevations above still water, usually with
sparse cover and sometimes in tussocks. The combination of species that settle in somewhat different zones is expected to lead to a larger vegetation
cover than a single species, but since all species have similar hydraulic roughness parameters , the differences in their
effects on the hydrodynamics may be negligible.
Our first scaling consideration for the selection of these species was geometric: roots in natural estuaries have lengths of a fraction of the main
channel depth, and our smallest plants still have relatively large roots. Our second scaling consideration dynamic: bank erosion and channel incision
should reduce. Earlier experiments with vegetation (pilots conducted for ; also see ) showed that more extensive
and interlocking root systems can completely fixate systems. In view of this experimental difficulty, we chose the smallest plants, which have
measurable bank erosion reduction effects as shown in in bespoke bank erosion tests at the scale of the experiments. The
third scaling consideration was also dynamic: vegetation should cause hydraulic resistance. As long as the stems penetrate the water surface and there
is sufficient stem density, the vegetation has a strong measurable hydraulic resistance effect . Unlike roots that are small
relative to channel depth in large natural estuaries, the vegetation settles at such high elevations that its effect on shallow flow can be large, as
shown in numerical modelling of meandering rivers with riparian vegetation and of estuaries with salt marsh vegetation
. However, the small scale of the 20 m flume does not allow for wave-driven salt marsh cliff erosion.
All seeds were soaked for 24 h to prevent floating and speed up germination. We released batches of seeds at the upstream boundary after every
dry bed photograph. After 10 tidal cycles for initial wetting, 12.5 g or about 10 000 seeds of Lotus pedunculatus were
supplied. Another 25 cycles later, 3.75 g or 15 000 seeds of Veronica beccabunga were allowed to disperse for 35 more tidal cycles
and then left to germinate for four days without tilting but with river discharge. These 70 tidal cycles were subtracted from the 1000 cycle run
before the next dry bed photography session so that spacing between the DEMs was exactly 1000 cycles. Laboratory conditions were around
300 lux light intensity by daylight-toned luminescent tube at all times as in the controlled vegetation experiments
of, a water temperature of 20 ∘C and a room temperature of about 17–20 ∘C. As found in earlier
experiments with vegetation, these conditions do not limit sprouting , whereas water depth and unrooting are limiting or
terminating growth . The plant species stopped growing after 5–7 d as the sand and tap water were free of
nutrients. Chlorine was added to the water as pest control to prevent algae, fungi and bacterial growth, and the bed of the vegetated experiment was
further treated with anti-algae spray at 8000, 10 500 and 12 500 tidal cycles. The duration of the vegetated experiment was about 2 months.
The other experiment with floodplain had a supply of suspended sediment from the start sufficient to simulate the formation of mudflats. Mud was
simulated by a supply of crushed nutshell at the river boundary. This experiment has already been published earlier and is
reanalysed here for different aspects. The nutshell had a diameter of 0.2 mm and a dry density of 1350 kgm-3. The large grain
size prevents mud infiltration in the sand bed. The nutshell was kept in suspension in a mixing tank and was supplied at a dry volume of about
1 mL per tidal cycle, such that a total volume of 0.013 m3 (0.009 m3 without 30 % pore space) was added to the experiment. Given a 20 s river discharge duration per tidal cycle, the mud concentration in the river influx was about 405 mgL-1 or 810 mg per tidal cycle. The velocity required to suspend the mud is much lower than that to suspend the sand, so that the mud can deposit in much shallower flow .
Numerical flow model
The numerical flow model Nays2D solves the shallow water equations to predict two-dimensional, depth-averaged flow velocity, water surface elevation
and water depth in small-scale systems such as flumes and was adapted to drive the flow by tilting as in the Metronome
. Previous results showed good similarity between modelled flow dynamics, flow velocity and depth measured by surface
particle imaging velocimetry and water colour-derived depth in a narrow estuary in the Metronome and in a meandering experiment
. This code was applied to the measured bathymetries and applied boundary conditions of the three experiments. The output of
water depth and flow velocity was calculated at a 2.5 cm by 2.5 cm resolution. For the sand and the mud experiments, the measured bed
surface elevation is taken with a spatially constant surface Manning roughness of 0.02 sm-1/6.
For the vegetated experiment, spatially variable roughness is applied depending on the estimated vegetation density as
in. Vegetation was filtered off the bed surface elevation map and assigned the same roughness as the other two experiments in
unvegetated areas. Vegetation roughness was assumed to be related to plant stem density and diameter. Earlier experiments demonstrated that, under
similar conditions, the flow velocity halved and the roughness doubled for a vegetation density of 2 plantscm-2. Here, the roughness is calculated from estimated stem density and known drag coefficient and stem thickness using the
relation for emergent vegetation, which was in good agreement with data for the species and conditions of the experiments
and has been used successfully in numerical models in rivers and estuaries
. A Manning coefficient of 0.1675 sm-1/6 was calculated for the densest vegetation with 20 stemscm-2
and a stem diameter of 0.5 mm. Based on this, a roughness map was created for each time step based on observed vegetation density.
Bathymetry after 13 000 tidal cycles (see movie in the Supplement). (a) Control experiment with only sand. (b) Estuary with two vegetation species, with cover indicated in green. For visualisation, the binary map of vegetation was dilated by a kernel of 4 mm in MATLAB. (c) Estuary with mud simulant, with cover indicated in black low mud experiment from.
ResultsMorphodynamics pattern and development
All three experiments developed a convergent, multi-channel estuary with mid-channel bars and a prograding ebb delta (Fig. , movie in
the Supplement). All three estuaries widened by bank erosion from the initial, monotonously converging estuary to form channels and
bars. As the original 0.03 m high sand bed flanking the estuaries was never flooded, vegetation or mud was confined in the reworked area. The
vegetation and the mud first settled on the most upstream mid-channel bars and the shore-connected bars (Fig. b and c). Seeds and mud
were also observed in suspension in the channels and deposited seaward of the ebb delta (without sprouting).
Time series of estuary depth below initial bed level, subdivided into an (a) upstream part (1–6 m), a (b) middle part (6–12 m) and a (c) downstream part (12–18 m). The 5th, 50th and 90th percentiles of depth represent the shoal depth, the median depth and the main channel depth.
The distribution of bed elevations generally broadened in the experiments (Fig. ). Initially, the carved channel
rapidly shallowed by deposition of sediment eroded from the banks. After about 1000 cycles, the channel–bar morphology had formed and the slower
process of estuary widening dominated the trend in bed elevations. Especially in the middle reach (Fig. b), the
muddy shoals increased several millimetres in elevation, while the channels deepened. All three estuaries gained a broader bed elevation distribution
with deeper channels and higher shoals in the middle reach, as well as in the upstream reach for the vegetated estuary. In the other reaches, the bed
elevation distributions did not broaden much. In the upstream reach (Fig. a), the channels became shallower for the
sandy estuary and stayed constant for the muddy estuary, while the vegetated estuary developed slightly deeper channels. In the downstream reach
(Fig. c), bed elevations increased particularly in the vegetated estuary, while the sandy estuary remained the
deepest.
Morphodynamics characterised by cumulative absolute bed level change over the duration of the experiments with (a) sand only, (b) vegetation and (c) mud. To avoid overrepresentation by the early stages, only data of every 1000 cycles were used. The absolute bed level change is calculated from the moment that an area became part of the estuary for the first time.
The channel–bar patterns and bed elevation distributions were caused by the morphodynamics of channel erosion and migration, estuarine bank erosion,
and bar formation and accretion. In the upstream and middle reaches of the vegetated and muddy estuaries, these processes reduced considerably
compared to the control experiment (Fig. ). As a result, the upstream 7 m of the estuaries with floodplain
remained narrower, especially of the muddy estuary. The downstream half of the estuaries (10–16 m) widened similarly but upstream the muddy
estuary showed fewer channels (Fig. ) and less lateral channel migration (Fig. ) than the
vegetated estuary and the control. The vegetated estuary was intermediate both in channel carving and migration
(Figs. and ).
The differences in dynamics are also visible in the virtual stratigraphy of the cross sections (Fig. ). These were constructed
from subsequent bed elevation maps, where eroded sediment is removed and accretion is assigned to the time step of the overlying bed surface elevation
for method, see. The cross sections of the sandy control show rapidly varying ages of deposits, indicating perpetual
avulsion of channels. For example, on the 15 m cross section, young (yellow) deposits overlie old (blue) deposits in several places with
3000–6000 tidal cycle hiatuses between them (Fig. ). The vegetated estuary (vegetation map shown in
Fig. ) tends to channel migration (more gradual lateral shifts in colour in Fig. ), while the muddy
estuary had predominantly vertically accreting bars (more gradual vertical shifts in colour in Fig. ), showing that the mud was
more effective at reducing bank erosion and lateral channel migration than the vegetation. This is surprising as the rooting depth of the vegetation
is approximately the same as the vegetation height and the typical channel depth , while the mud layer thickness is an order
of magnitude smaller .
Morphodynamics characterised by the age of the largest depth reached over the duration of the experiments with (a) sand only, (b) vegetation and (c) mud. To avoid overrepresentation by the early stages, only data of every 1000 cycles were used.
Cross sections with age of deposit, showing the preserved part of the morphodynamics in the experiments with (a) sand only, (b) vegetation and (c) mud for four positions along the flume. Bed level data of every 1000 cycles were used. The vertical exaggeration is a factor of 50.
Hydrodynamics: flow and tidal prism
The hydrodynamics were characterised by numerical flow modelling for the entire tidal cycle on all measured bathymetries. At the end of the
experiments, the path of maximum ebb velocity during the tidal cycle mainly follows a single channel in all three experiments, while the maximum flood
velocity is more distributed over the width in the downstream half of the estuaries (Fig. ). The velocities in the upper reach
are generally lower in the muddy estuary than in the sandy estuary, which is consistent with the lower width and the raised bars due to mud
deposition. The vegetated estuary shows a striking difference with the other two: the higher friction of vegetated bars (compare
Fig. and Fig. ) not only strongly reduces flow velocity over the bars but also focuses the flow into
the channels, even during the flood phase when the flow enters the system through the relatively wide mouth unconfined by vegetation.
Example of Manning roughness values for the vegetated estuary specified as initial condition in the Nays2D flow model at 13 000 tidal cycles. The Manning roughness of the unvegetated bed is set to 0.02 in the model.
Maximum flow velocity calculated with the flow model for the end of the experiment during ebb (left) and during flood (right) in the experiments with (a) sand only, (b) vegetation and (c) mud.
Distributions of maximum flow velocity calculated with the flow model for the middle (top; a, b) and end (bottom; c, d) of the experiment during flood (left; a, c) and during ebb (right; b, d). The selected depths ranges are shoal areas (0 > depth >-0.015 m) and channels (-0.025 > depth >-0.050 m) for consistency with Fig. b.
To quantify the velocity patterns in the different experiments, the probability distributions of the flow velocities were calculated for the
shallowest and the deepest parts in middle section, halfway and at the end of the experiments (Fig. , consistent with the bed
elevation percentiles shown in Fig. b). Comparison between 6000 cycles and 13 000 cycles shows a reduction of flow
velocity over the shoals (peaks of dashed lines shift to the left from Fig. a–c and b–d), especially for the vegetated
experiment. Comparison between the experiments at 6000 cycles shows that the ebb and flood flow in the channels is lower and broader distributed in
the muddy estuary (Fig. a and b), although the differences disappear towards the end of the experiments. The vegetated
experiment, on the other hand, has a narrower distribution of high channel velocities around 6000 cycles and lower velocities over the shoals around
13 000 cycles when vegetation has settled more generally. Clearly, the vegetation baffles the flow by resistance on the shoals, whereas the mud
reduces the flow in the channels.
Tidal prism along the experimental estuaries calculated with the flow model after (a) 6000 and (b) 13 000 tidal cycles. The estuary mouth is located at 17.8 m.
The tidal prism here is calculated not only at the mouth but at every cross section along the estuary to resolve more local effects of vegetation and
mud as in. The tidal prism increased in the seaward direction and is the smallest for the experiments with mud
(Fig. ). After 13 000 tidal cycles, the tidal prism in the mouth area had been reduced most for the muddy experiment compared to
6000 cycles. In the upstream half, the change was small though the tidal prism of the sandy experiment continued to increase. This is consistent with
the rapid initial channel shallowing of the upstream reach (Fig. ) and concurrent infilling
(Fig. ). The spatial variation in tidal prism shows the expected seaward increase, but has a superimposed secondary pattern of
around 6 and 14 m. A similar pattern in the tidal prism for the initial bathymetry (dashed line in Fig. ) shows
that this is likely due to a deviation of tidal flow generation by tilting from what is generally observed in natural estuaries. Specifically, in
somewhat more open and less filled reaches of the estuaries, the tilting was visually observed to drive periodically reversing flow regardless of the
connection with the mouth. Nevertheless, the landward penetration of tides reduced over time. The most upstream flood flow velocity was slightly
increased in the muddy estuary as the flow was concentrated in a single channel and flanked by high mudflats, but the tidal prism at that point was
already smaller than in the other estuaries.
Time series of tidal prism calculated along the estuaries at distances of (a) 6 m, (b) 12 m and close to the mouth at (c) 17.5 m.
Development of vegetation cover as a fraction of estuary width along the estuary.
The tidal prism changed through time in different ways along the estuaries (Fig. ). The vegetated estuary shows a pattern of
gradual increase of tidal prism until about 8000 cycles followed by a decrease, unlike the muddy estuary. This can also be seen in the higher channel
velocities in the vegetated experiment (Fig a and b). The decrease is due to the gradual increase in vegetation cover as the
recruitment events added up and the particularly rapid expansion of vegetation after about 8000 cycles (Fig. ). The tidal prism
reduction in the vegetated estuary is approaching that of the muddy estuary after 13 000 cycles (Fig. ).
In the upstream reach, the tidal prism of the muddy estuary decreased rapidly in the beginning and then remained about constant, while that of the
sandy and vegetated estuaries decreased almost linearly over time. In the middle reach, the tidal prism in the muddy and vegetated estuaries stay
about constant after initially increasing, while that of the sandy estuary continues to increase nearly until the end. Close to the mouth, the tidal
prism is the collective result of tides along the entire estuary (excluding a secondary effect of locally generated tides due to tilting) and the
differences in magnitude are largest. Here, the tidal prism in the muddy estuary is about 70 % of that in the sandy estuary, and the tidal prism
in the vegetated estuary is about 80 % of the sandy control. The muddy and sandy estuaries initially show a rapid increase of tidal prism but then
a gradual, though limited, decline. This is consistent with the declining depth in the downstream reach
(Fig. c). The ongoing decline at the end of the experiments shows that complete convergence to some equilibrium has
not yet taken place, but the slowing of the change and the opposite trends in the middle and downstream reaches suggest that convergence is nearly
reached.
The reduced tidal prism in the muddy estuary is not merely lower than in the other two experiments because a certain volume of mud was supplied that
filled space: the total added volume was 0.013 m3, which is only a small fraction of the difference in tidal prism between the sandy and the
muddy estuary (about 0.5 m3 during nearly the entire experiment; see Fig. ). As also shown in
Figs. and , mud had a system-wide morphological effect on the estuary development. Regardless of
the limited filling, the mud reduced the overall channel dynamics.
It is unclear whether the muddy estuary would have filled further, given the mud supply, or whether the vegetated estuary would have been covered more
extensively by vegetation with the ongoing seed distribution events. While the tidal prism at the end of the experiments was nearly constant, the
erosion of the outer estuary banks, composed of pure sand, and the expansion of the ebb delta could slowly continue.
Discussion
The experimental results suggest that vegetation and mud have similar, but not the same, effects on the landward reduction of tidal energy. While mud
fills the accommodation space to reduce lateral channel mobility and overall depth, the vegetation reduces the tidal flow through vegetation-enhanced
hydraulic resistance on the bars. Surprisingly, the enhanced bar accretion in the muddy estuary reduced the flow over the bars less than the
vegetation despite the blockage of flow by mudflats growing up to the water surface. This is all the more surprising as the shallow experimental flow
was probably laminar for a greater part of the tidal cycle which would increase friction
more.
Equally surprising is the much lower lateral channel mobility in the muddy estuary than in the vegetated estuary, despite the fact that the mud
deposits were much thinner than the channels , which allows for unhindered undercutting into the non-cohesive sand underlying the
mud. Moreover, the mud simulant, nutshell, was not found to be so cohesive as to reduce bank erosion, even after tens of days Fig. 13
in. In contrast, the vegetation rooting could be as deep as the channels but this did not reduce the mobility
as much as the mud. The growth of vegetation in nature, even at supratidal level, causes high hydraulic resistance during high tides. This leads to
reduced and often negligible flow on the bars and strong focusing of the flow in the channels as also found in experiments and models of river systems
. This happens particularly in upstream reaches of the estuaries where the
tidal dynamics are reduced sufficiently for vegetation to settle and where estuarine bars and mudflats can accrete to the high intertidal and
supratidal levels required for vegetation to settle . Possibly, the difference between the vegetated and
muddy estuary occurred because the cover of vegetation is smaller than that of the mud, and given more time for recruitment, or more
inundation-resistant species, the vegetated estuary might have developed more similar to the muddy estuary.
A morphological effect of filling of intertidal space by vegetation and by mud is the reduction of flow shear stress over tidal bars. This, combined
with the increased resistance against erosion, diminishes the likelihood of cross-cutting bars by channels that would otherwise lead to channel
braiding as in rivers . The data suggest a transition from rapid channel avulsion in the control experiment to more gradual
migration in the vegetated estuary and even reduced migration in the muddy estuary. This filling effect is paralleled in river morphodynamics, where
initiation of braiding through chute cutoffs and unhindered bank erosion is inhibited by the vegetation and the suspended sediment on the
floodplain. Both vegetation and suspended sediment reduce excess shear stress, regardless of whether the sediment is cohesive
or non-cohesive . Note that this inner-bend effect of floodplain on the channel pattern, unlike the usually
assumed outer-bend effect of floodplain on eroding bank stability, works both in rivers and estuaries without rooting depth or cohesive sediment
deposit thickness that approach the channel depth .
The reduced channel-cutting tendency is important for overall floodplain formation because channel migration is far more effective in eroding bars and
removing floodplain by undercutting than direct overflow during flooding. This was also evident from the comparison between the vegetated and muddy
estuaries with the control. In field data, the tendency to stabilise channels and reduce bar cross cutting was also observed
, but the effects of increasing sand bar height, mud sedimentation and salt marsh expansion occurred
simultaneously and could not be separated, unlike in the experiments.
Moreover, channels in this kind of estuary are rather fixated in position by topographic forcing around the mid-channel bars ,
by embankments and by dredging , which nearly completely removes the tendency of floodplain destruction by lateral bank
undercutting. This was also observed in earlier experiments . As floodplain accretion continues unhindered,
cohesive sediment and plants also reduce potential erosion by overflow . Where the more protected flanks of the estuary fill
up with fine sediment and vegetation, the tidal system is effectively constrained and reducing in dimensions as observed in reconstructions
and models . This contrasts with idealised modelling and with
experimental conditions with cohesionless sediment only so that tidal systems continue to expand . The only factor
reducing the tidal penetration in such estuaries is tidal damping by the bottom friction, which increases with reducing depth and shallow bar area,
and tidal damping by large flood storage areas, which were absent in the experiments. Friction increased more rapidly in the muddy estuary due to the
upstream depth reduction and in the vegetated estuary due to the large relative surface area with vegetation (Fig. ).
The two local effects outlined above, namely space filling and flow resistance, could both lead to lateral constraints and tidal channel deepening,
which would, according to the first hypothesis, enhance tidal energy in the upstream estuary and cause further landward tidal penetration. However,
the long-term development of tidal prism in the experiments demonstrates a positive feedback between floodplain formation, channel morphodynamics and
planform dimensions. The filling of intertidal space, on the fluvial–tidal transition and further seaward, causes a reduction in the tidal prism at
the mouth as hypothesised in . This is especially the case for the muddy estuary where accommodation space is filled with mud and
less so for the vegetated area where flow over vegetated bars is reduced. Regardless of the widening by cohesionless bank erosion, the tidal prism at
the mouth reduced due to vegetation and mud and, in the sandy estuary, also due to shallowing. In this sense, the effect of floodplain formation on an
estuary is different from that on a river, where the flow discharge is locally unchanged, as it is externally imposed. As a result, narrowing of the
channel in rivers leads to deepening at the same time, which is not the case in estuaries where the tidal discharge reduces.
The long-term effects of a reducing tidal prism on the system are considerable. Where a large intertidal area with delayed outflow would have enhanced
ebb dominance and sediment export, its reduction due to floodplain formation may invert this tendency and cause import
. In other words, mud deposition and vegetation development could cause large-scale filling of
estuaries. This direct link between sedimentation and vegetation growth and estuary dimensions dominates over the effect of floodplain observed in
rivers because the total flow discharge is modified in estuaries as opposed to rivers. As filling progresses, this zone is expected to shift seawards
as well, but testing this requires experiments with initially unfilled estuaries and control on the landward sediment transport from the coastal
zone. Experiments conducted in the Metronome after the experiments reported in this paper indeed showed that mud and vegetation tend to fill up
unfilled lagoonal estuaries, given sufficient sediment supply from the coastal and fluvial boundaries for the vegetation to settle
. While those experiments show how coastal plains drowned in the early Holocene could fill, the present findings
explain how initial estuary widening, following catastrophic ingression , can be followed by constraints on widening, and even
filling, despite the easily erodible banks that would have led to unbridled braiding in the case of a river.
In nature, the calm hydrodynamic conditions in salt marshes could be conducive to organic matter production and storage leading to organogenic rise of
the bed surface and mangrove , especially once removed further from the active estuary where
precipitation leads to freshwater conditions. The latter is more likely to occur upstream in the estuary, which is consistent with the observation
that Holocene tidal systems filled in the seaward direction . On the other hand, past numerical modelling
showed that the reduced flow through the canopy also reduces the sedimentation within the vegetation . This could not be
tested well in the experiments, because the flow through the experimental vegetation becomes laminar and would not suspend sediment
. Regardless, the collective effect of the different mechanisms of mudflat and salt marsh formation is to constrain estuary
dimensions and, given sediment input, the forming of land.
Conclusions
Salt marsh and mudflat development in natural multi-channel estuaries reduces system dimensions and dynamics similarly to floodplains in rivers, but
vegetation and mud have subtly different effects. While both reduce the overall tidal prism along the estuary and the landward penetration of tides,
vegetation effectively concentrates tidal flow in channels, while mud effectively reduces lateral channel migration. Both likely contributed to the
lateral constraining of estuaries initially formed by catastrophic ingression. Salt marsh expansion and mudflat sedimentation also have indirect
effects on the tidal dynamics, which lead to a positive feedback of enhanced filling through the reduction in tidal prism, an effect absent in rivers
where flow discharge is mostly independent of the upstream morphology.
Data availability
Bathymetries (in MATLAB and in NetCDF4 formats) and numerical flow model input and output data are made available at Yoda 10.24416/UU01-R3ZUC9.
Video supplement
This includes the video of overhead imagery for all three experiments (10.24416/UU01-R3ZUC9; ).
Author contributions
MGK conceived the study, built the laboratory facility and wrote the manuscript. LR made most of the figures. LR, IRL and LB conducted the experiments supervised by MGK. SAHW and LR set up and conducted the numerical models with help of SAHW, and all authors contributed to idea development and manuscript preparation.
Competing interests
The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
William Kearney, associate editor Claire Masteller and an anonymous reviewer are gratefully acknowledged for helpful comments and questions for clarification of methods. Technical support of the experiments by the lab technicians and help during experiments by Jasper Leuven are cordially acknowledged. Discussion with Harm Jan Pierik and Tjalling de Haas helped to improve the discussion.
Financial support
This research has been supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (grant no. Vici 016.140.316/13710) and the European Research Council, H2020 European Research Council (ESTUARIES (grant no. 647570)).
Review statement
This paper was edited by Claire Masteller and reviewed by William Kearney and one anonymous referee.
References
Ashmore, P.:
How do gravel-bed rivers braid?, Can. J. Earth Sci., 28, 326–341, 1991.Baar, A. W., Boechat Albernaz, M., van Dijk, W. M., and Kleinhans, M. G.:
Critical dependence of morphodynamic models of fluvial and tidal systems on empirical downslope sediment transport, Nat. Commun., 10, 4903, 10.1038/s41467-019-12753-x, 2019.
Baptist, M., Babobic, V., Rodriguez Uthurburu, J., Keijzer, M., Uittenbogaard, R., Mynett, A., and Verwey, A.:
On inducing equations for vegetation resistance, J. Hydraul. Res., 45, 1–16, 2006.Boechat Albernaz, M., Roelofs, L., Pierik, H. J., and Kleinhans, M. G.:
Complementing scale experiments of rivers and estuaries with numerically modelled hydrodyanmics, Earth Surf. Proc. Land., 45, 3824–3841, 10.1002/esp.5003, 2020.Bouma, T., van Belzen, J., Balke, T., Zhu, Z., Airoldi, L., Blight, A., Davies, A., Galvan, C., Hawkins, S., Hoggart, S., Lara, J., Losada, I., Maza, M., Ondiviela, B., Skov, M., Strain, E., Thompson, R., Yang, S., Zanuttigh, B., Zhang, L., and Herman, P.:
Identifying knowledge gaps hampering application of intertidal habitats in coastal protection: Opportunities and steps to take, Coast. Eng., 87, 147–157, 10.1016/j.coastaleng.2013.11.014, 2014.Braat, L., van Kessel, T., Leuven, J. R. F. W., and Kleinhans, M. G.:
Effects of mud supply on large-scale estuary morphology and development over centuries to millennia, Earth Surf. Dynam., 5, 617–652, 10.5194/esurf-5-617-2017, 2017.Braat, L., Leuven, J., Lokhorst, I., and Kleinhans, M.:
Effects of estuarine mudflat formation on tidal prism and large-scale morphology in experiments, Earth Surf. Proc. Land., 44, 417–432, 10.1002/esp.4504, 2019.
Braudrick, C., Dietrich, W., Leverich, G., and Sklar, L.:
Experimental evidence for the conditions necessary to sustain meandering in coarse-bedded rivers, P. Natl. Acad. Sci. USA, 106, 936–941, 2009.Brown, J. and Davies, A.:
Flood/ebb tidal asymmetry in a shallow sandy estuary and the impact on net sand transport, Geomorphology, 114, 431–439, 10.1016/j.geomorph.2009.08.006, 2010.Brückner, M., Braat, L., Schwarz, C., and Kleinhans, M.:
What came first, mud or biostabilizers? Elucidating interacting effects in a coupled model of mud, saltmarsh, microphytobenthos and estuarine morphology, Water Resour. Res., 56, e2019WR026945, 10.1029/2019WR026945, 2020.Brückner, M. Z., McMahon, W. J., and Kleinhans, M. G.:
Muddying the waters: modeling the effects of early land plants in paleozoic estuaries, Palaios, 36, 173–181, 10.2110/palo.2020.073, 2021.
Dalrymple, R. and Choi, K.:
Morphologic and facies trends through the fluvial-marine transition in tide-dominated depositional systems: A schematic framework for environmental and sequence-stratigraphic interpretation, Earth-Sci. Rev., 81, 135–174, 2007.de Haas, T., Pierik, H., van der Spek, A., Cohen, K., van Maanen, B., and Kleinhans, M.: Long-term evolution of tidal systems: effects of rivers, coastal boundary conditions, eco-engineering species, inherited relief and human interference, Earth-Sci. Rev., 177, 139–163, 10.1016/j.earscirev.2017.10.006, 2018.de Haas, T., van der Valk, L., Cohen, K., Pierik, H., Weisscher, S., Hijma, M., van der Spek, A., and Kleinhans, M.:
Long-term evolution of the Old Rhine estuary: Unravelling effects of changing boundary conditions and inherited landscape, Depositional Record, 5, 84–108, 10.1002/dep2.56, 2019.Dronkers, J.: Convergence of estuarine channels, Cont. Shelf Res., 144, 120–133, 10.1016/j.csr.2017.06.012, 2017.FitzGerald, D. M. and Hughes, Z.:
Marsh Processes and Their Response to Climate Change and Sea-Level Rise, Annu. Rev. Earth Pl. Sc., 47, 481–517, 10.1146/annurev-earth-082517-010255, 2019.
Friedrichs, C. and Aubrey, D.:
Non-linear tidal distortion in shallow well-mixed estuaries: a synthesis, Estuar. Coast. Shelf S., 27, 521–545, 1988.
Gran, K. and Paola, C.:
Riparian vegetation controls on braided stream dynamics, Water Resour. Res., 37, 3275–3283, 2001.Kirwan, M., Temmerman, S., Skeehan, E., Guntenspergen, G., and Fagherazzi, S.:
Overestimation of marsh vulnerability to sea level rise, Nat. Clim. Change, 6, 253–260, 10.1038/NCLIMATE2909, 2016.Kleinhans, M., de Vries, B., Braat, L., and van Oorschot, M.:
Living landscapes: muddy and vegetated floodplain effects on fluvial pattern in an incised river, Earth Surf. Proc. Land., 43, 1618–1632, 10.1002/esp.4437, 2018.Kleinhans, M. G.: Sorting out river channel patterns, Prog. Phys. Geog., 34, 287–326, 10.1177/0309133310365300, 2010.Kleinhans, M. G., van Dijk, W., van de Lageweg, W., Hoyal, D., Markies, H., van Maarseveen, M., Roosendaal, C., van Weesep, W., van Breemen, D., Hoendervoogt, R., and Cheshier, N.:
Quantifiable effectiveness of experimental scaling of river- and delta morphodynamics and stratigraphy, Earth-Sci. Rev., 133, 43–61, 10.1016/j.earscirev.2014.03.001, 2014.Kleinhans, M. G., Braudrick, C., van Dijk, W., van de Lageweg, W., Teske, R., and van Oorschot, M.:
Swiftness of biomorphodynamics in Lilliput- to Giant-sized rivers and deltas, Geomorphology, 244, 56–73, 10.1016/j.geomorph.2015.04.022, 2015a.Kleinhans, M. G., Terwisscha van Scheltinga, R., van der Vegt, M., and Markies, H.:
Turning the tide: growth and dynamics of a tidal basin and inlet in experiments, J. Geophys. Res.-Earth, 120, 95–119, 10.1002/2014JF003127, 2015b.Kleinhans, M. G., van der Vegt, M., Leuven, J., Braat, L., Markies, H., Simmelink, A., Roosendaal, C., van Eijk, A., Vrijbergen, P., and van Maarseveen, M.:
Turning the tide: comparison of tidal flow by periodic sea level fluctuation and by periodic bed tilting in scaled landscape experiments of estuaries, Earth Surf. Dynam., 5, 731–756, 10.5194/esurf-5-731-2017, 2017.Kleinhans, M. G., Roeelofs, L., Weisscher, S. A. H., Lokhorst, I. R., and Braat, L.: Data supplement to “Estuarine morphodynamics and development modified by floodplain formation”, YODA [data set],
10.24416/UU01-R3ZUC9, 2021.Leuven, J., de Haas, T., Braat, L., and Kleinhans, M. G.:
Topographic forcing of tidal sandbar patterns for irregular estuary planforms, Earth Surf. Proc. Land., 43, 172–186, 10.1002/esp.4166, 2018a.Leuven, J., Verhoeve, S., Van Dijk, W., Selaković, S., and Kleinhans, M.:
Empirical Assessment Tool for Bathymetry, Flow Velocity and Salinity in Estuaries Based on Tidal Amplitude and Remotely-Sensed Imagery, Remote Sens.-Basel, 10, 10.3390/rs10121915, 2018b.
Leuven, J. R., Pierik, H. J., van der Vegt, M., Bouma, T. J., and Kleinhans, M.:
Sea-level-rise-induced threats depend on the size of tide-influenced estuaries worldwide, Nat. Clim. Change, 9, 986–992, 2019.Leuven, J. R. F. W., Braat, L., van Dijk, W. M., de Haas, T., van Onselen, E. P., Ruessink, B. G., and Kleinhans, M. G.:
Growing Forced Bars Determine Nonideal Estuary Planform, J. Geophys. Res.-Earth, 123, 2971–2992, 10.1029/2018JF004718, 2018.Lokhorst, I. R., Braat, L., Leuven, J. R. F. W., Baar, A. W., van Oorschot, M., Selaković, S., and Kleinhans, M. G.:
Morphological effects of vegetation on the tidal–fluvial transition in Holocene estuaries, Earth Surf. Dynam., 6, 883–901, 10.5194/esurf-6-883-2018, 2018.Lokhorst, I. R., de Lange, S. I., van Buiten, G., Selakovic, S., and Kleinhans, M. G.:
Species selection and assessment of eco-engineering effects of seedlings for biogeomorphological landscape experiments, Earth Surf. Proc. Land., 44, 2922–2935, 10.1002/esp.4702, 2019.Robins, P. and Davies, A.:
Morphological controls in sandy estuaries: the influence of tidal flats and bathymetry on sediment transport, Ocean Dynam., 60, 503–517, 10.1007/s10236-010-0268-4, 2010.Savenije, H.:
Prediction in ungauged estuaries: an integrated theory, Water Resour. Res., 51, 2464–2476, 10.1002/2015WR016936, 2015.Stark, J., Smolders, S., Meire, P., and Temmerman, S.:
Impact of intertidal area characteristics on estuarine tidal hydrodynamics: A modelling study for the Scheldt Estuary, Estuar. Coast. Shelf S., 198, 138–155, 10.1016/j.ecss.2017.09.004, 2017.Swinkels, C., Jeuken, C., Wang, Z., and Nicholls, R.:
Presence of connecting channels in the Western Scheldt Estuary, J. Coastal Res., 25, 627–640, 10.2112/06-0719.1, 2009.Tal, M. and Paola, C.: Effects of vegetation on channel morphodynamics: results and insights from laboratory experiments, Earth Surf. Proc. Land., 35, 1014–1028, 10.1002/esp.1908, 2010.Tamura, T., Saito, Y., Bateman, M., Nguyen, V., Ta, T., and Matsumoto, D.:
Luminescence dating of beach ridges for characterizing multi-decadal to centennial deltaic shoreline changes during Late Holocene, Mekong River delta, Mar. Geol., 326–328, 140–153, 10.1016/j.margeo.2012.08.004, 2012.van de Lageweg, W., van Dijk, W., and Kleinhans, M.:
Channel belt architecture formed by an experimental meandering river, Sedimentology, 60, 840–859, 10.1111/j.1365-3091.2012.01365.x, 2013.
van der Spek, A.:
Tidal asymmetry and long-term evolution of Holocene tidal basins in The Netherlands: simulation of palaeo-tides in the Schelde estuary, Mar. Geol., 141, 71–90, 1997.van der Wegen, M.: Numerical modeling of the impact of sea level rise on tidal basin morphodynamics, J. Geophys. Res. Earth Surf., 118, 447–460, 10.1002/jgrf.20034, 2013.van Dijk, W., Teske, R., van de Lageweg, W., and Kleinhans, M. G.:
Effects of vegetation distribution on experimental river channel dynamics, Water Resour. Res., 49, 7558–7574, 10.1002/2013WR013574, 2013a.van Dijk, W. M., Van de Lageweg, W. I., and Kleinhans, M. G.:
Formation of a cohesive floodplain in a dynamic experimental meandering river, Earth Surf. Proc. Land., 38, 1550–1565, 10.1002/esp.3400, 2013b.van Dijk, W. M., Cox, J. R., Leuven, J. R., Cleveringa, J., Taal, M., Hiatt, M. R., Sonke, W., Verbeek, K., Speckmann, B., and Kleinhans, M. G.:
The vulnerability of tidal flats and multi-channel estuaries to dredging and disposal, Anthropocene Coasts, 4, 36–60, 10.1139/anc-2020-0006, 2021.van Oorschot, M., Kleinhans, M., Geerling, G., and Middelkoop, H.:
Distinct patterns of interaction between vegetation and morphodynamics, Earth Surf. Proc. Land., 41, 791–808, 10.1002/esp.3864, 2016.
Vos, P. and van Kesteren, W.:
The long-term evolution of intertidal mudflats in the northern Netherlands during the Holocene; natural and anthropogenic processes, Cont. Shelf Res., 20, 1687–1710, 2000.Wang, Z., van Maren, D., Ding, P., Yang, S., van Prooijen, B., de Vet, P., Winterwerp, J., de Vriend, H., Stive, M., and He, Q.:
Human impacts on morphodynamic thresholds in estuarine systems, Cont. Shelf Res., 111, 174–183, 10.1016/j.csr.2015.08.009, 2015.Weisscher, S., Van den Hoven, K., Pierik, H., and Kleinhans, M.:
Building and raising land: mud and vegetation effects in infilling estuaries, J. Geophys. Res.-Earth, 127, e2021JF006298, 10.1029/2021JF006298, 2022.Weisscher, S. A. H., Shimizu, Y., and Kleinhans, M. G.: Upstream perturbation and floodplain formation effects on chute-cutoff-dominated meandering river pattern and dynamics, Earth Surf. Proc. Land., 44, 2156–2169, 10.1002/esp.4638, 2019.Weisscher, S. A. H., Boechat-Albernaz, M., Leuven, J. R. F. W., Van Dijk, W. M., Shimizu, Y., and Kleinhans, M. G.:
Complementing scale experiments of rivers and estuaries with numerically modelled hydrodynamics, Earth Surf. Dynam., 8, 955–972, 10.5194/esurf-8-955-2020, 2020.Whitfield, A., Elliott, M., Basset, A., Blaber, S., and West, R.:
Paradigms in estuarine ecology – A review of the Remane diagram with a suggested revised model for estuaries, Estuar. Coast. Shelf S., 97, 78–90, 10.1016/j.ecss.2011.11.026, 2012.Woodroffe, C., Rogers, K., McKee, K., Lovelock, C., Mendelssohn, I., and Saintilan, N.:
Mangrove Sedimentation and Response to Relative Sea-Level Rise, Annu. Rev. Mar. Sci., 8, 243–266, 10.1146/annurev-marine-122414-034025, 2016.Ysebaert, T., van der Hoek, D.-J., Wortelboer, R., Wijsman, J. W., Tangelder, M., and Nolte, A.:
Management options for restoring estuarine dynamics and implications for ecosystems: A quantitative approach for the Southwest Delta in the Netherlands, Ocean Coast. Manage., 121, 33–48, 10.1016/j.ocecoaman.2015.11.005, 2016.Zhou, Z., Coco, G., Jimenez, M., Olabarrieta, M., van der Wegen, M., and Townend, I.:
Morphodynamics of river-influenced back-barrier tidal basins: The role of landscape and hydrodynamic settings, Water Resour. Res., 50, 9514–9535, 10.1002/2014WR015891, 2014.