Coastal lowlands and estuaries are subjected to increasing flood risks during storm surges due to global and regional changes. Tidal wetlands are increasingly valued as effective natural buffers for storm surges by dissipating wave energy and providing flood water storage. While previous studies focused on flood wave attenuation within and behind wetlands, this study focuses on the effects of estuarine wetland properties on the attenuation of a storm tide that propagates along the length of an estuary. Wetland properties including elevation, surface area, and location within the estuary were investigated using a numerical model of the Scheldt estuary (Belgium, SW Netherlands). For a spring tide lower wetland elevations result in more attenuation of high water levels along the estuary, while for a higher storm tide higher elevations provide more attenuation compared to lower wetland elevations. For spring and storm tide a larger wetland surface area results in a better attenuation along the estuary up to a threshold wetland size for which larger wetlands do not further contribute to more attenuation. Finally a wetland of the same size and elevation, but located more upstream in the estuary, can store a larger proportion of the local flood volume and therefore has a larger attenuating effect on upstream high water levels. With this paper we aim to contribute towards a better understanding and wider implementation of ecosystem-based adaptation to increasing estuarine flood risks associated with storms.
Deltas and estuaries worldwide are subjected to an increasing risk of flooding disasters. Recent examples of storm surge flood disasters include Hurricane Katrina in New Orleans (2005), Hurricane Sandy in New York (2012) and Typhoon Haiyan in the central Philippines (2013). Global climate warming accelerates sea level rise (Church et al., 2013) and there is growing consensus that it increases the intensity of storms and associated storm surges (Emanuel, 2005; Webster et al., 2005; Knutson et al., 2010; Lin et al., 2012; Grinsted et al., 2013). These global change effects combined with regional changes, including growing coastal populations and land subsidence, expose low-lying deltas and estuaries to increasing storm surge flood risks (Hanson et al., 2011; Hallegatte et al., 2013). It is estimated that in 2005 around 40 million people and USD 3000 billion of assets were located in flood-prone coastal cities, and that these numbers could increase to 150 million people and USD 35 000 billion of assets by 2070 (Hanson et al., 2011).
For low-lying coasts, deltas and estuaries, the presence of wetland ecosystems, such as tidal marshes and mangroves, is considered to play an important role in the landward attenuation of storm surge flood waves (Temmerman et al., 2013). An increasing number of studies shows that these wetland ecosystems can reduce the height of storm surges due to their dense vegetation and the friction they exert on landward propagating flood waves (Costanza et al., 2008; Das and Vincent, 2009; Krauss et al., 2009; Wamsley et al., 2010; Gedan et al., 2011; Temmerman et al., 2012; Zhang et al., 2012; Sheng et al., 2012; Liu et al., 2013; Barbier et al., 2013; Hu et al., 2015). Unfortunately these wetland ecosystems are and have been under a lot of pressure due to urban, agricultural and industrial expansion. An estimated 25 % of the world's intertidal estuarine habitat has been lost due to land reclamation (French, 1997). Over the last 3 decades globally, 50 % of salt marsh and 35 % of mangrove ecosystems have been lost or have degraded due to human activity (Alongi, 2002; Millenium Ecosystem Assessment, 2005). With the loss of these wetland ecosystems, their protective function is also lost, making low-lying coasts, deltas and estuaries more vulnerable to flood risks. Although traditional hard engineering solutions, like dikes and storm surge barriers, are widely perceived as the ultimate solution to combat flood risks, the combination with ecosystem-based flood defense – i.e. the conservation and restoration of tidal wetlands for flood attenuation – is likely to be more sustainable and cost-effective in the long term for many critical deltas and estuaries worldwide (Temmerman et al., 2013). Wetland ecosystems can effectively keep up with sea level rise by natural sediment accretion (Kirwan and Temmerman, 2009; Kirwan et al., 2010), which makes them in the long term more sustainable than static engineering structures. Apart from storm surge attenuation, wetlands also effectively attenuate wind waves and associated shoreline erosion, even during extreme wave conditions (Möller et al., 2014). They further deliver additional ecosystem services like carbon sequestration, water quality regulation, nutrient cycling, biological production, and many others (Barbier et al., 2011). Only a limited number of places around the world are implementing that ecosystem-based flood defense on a large scale, such as the conversion of reclaimed land back into wetlands in the UK (Pethick, 2002; French, 2006) and Belgium (Cox et al., 2006; Maris et al., 2007; Vandenbruwaene et al., 2011), and more in-depth research is needed in order to stimulate and support more widespread implementation (Temmerman et al., 2013). Harmonizing hard engineering solutions, where necessary, with ecosystem-based solutions, where possible, is a challenge for future development and sustainability of estuaries and coasts.
Existing studies on the degree of storm surge height reduction by tidal
wetlands rely on fragmentary observations of peak water levels within and
behind tidal wetlands (e.g. Lovelace, 1994; McGee et al., 2006; Krauss et
al., 2009; Wamsley et al., 2010) and on hydrodynamic model simulations (e.g. Resio
and Westerink, 2008; Loder et al., 2009; Wamsley et al., 2009, 2010; Zhang et
al., 2012; Temmerman et al., 2012; Sheng et al., 2012; Liu et al., 2013; Hu et
al., 2015). These studies reported that peak water levels can be reduced by 5
to 50 cm km
Whereas the existing studies discussed above focus on
In this paper we present a 2-D hydrodynamic model that has been calibrated
and validated for tides and storm tides in the Scheldt estuary (Belgium and
SW Netherlands), and use the model to investigate the effect of scenarios
with different intertidal wetland properties, including wetland elevation,
size and location along the estuary, on
The Scheldt estuary is situated in the Netherlands and Belgium (Fig. 1a). The connection with the North Sea is formed by the Vlakte van de Raan, a shallow water area with several deeper channels. The part of the estuary from the mouth until the Dutch/Belgian border (located at 67 km from the mouth, measured along the thalweg) is called Western Scheldt and is characterized by different ebb and flood channels surrounding large intertidal sand and mud flats. The part further upstream from the border until Ghent (located at 170 km from the mouth) (Fig. 1b) is called Sea Scheldt and is characterized by a single channel bordered by much smaller intertidal flats and marshes.
The estuary mouth near Vlissingen (km 2) (Fig. 1b) is approximately 5 km
wide and flood enters twice a day with an average flood volume of 1.04 Gm
The Scheldt estuary occupies an area of 35 000 ha of which 9000 ha (26 %)
is intertidal area (yellow area in Fig. 1b) and 2800 ha (8 %) are
vegetated marshes located above mean high water level (MHWL) (dark green
areas in Fig. 1b) (Vandenbruwaene et al., 2013). Saeftinghe (51.33
In the Scheldt estuary large-scale historic wetland reclamation – about 100 000 ha of intertidal area has been reclaimed over the last 1000 years of which about 15 000 ha since 1800 – has largely contributed to increasing water levels in the remaining estuarine channels. Nowadays tidal marsh restoration on formerly reclaimed land is brought into large-scale practice as an essential part of the flood defense system (Meire et al., 2014). Saeftinghe is the largest remaining intertidal area today with 3000 ha.
A 2-D hydrodynamic model of the Scheldt estuary and tidal tributaries was
made with the finite element model software TELEMAC 2-D (v6p0). The
hydrodynamics are solved by the depth-averaged Saint-Venant equations. The
triangular irregular mesh consists of 42 632 nodes and 75 889 elements. Its
resolution ranges from 300 m in the mouth area to 150 m in the Western
Scheldt and 50 m near Antwerp (km 91) to 5 m at the upstream end of the
model. The bathymetry of 2009 (20 by 20 m resolution from mouth to the
border and 1 by 1 m further upstream, obtained from single- and multibeam
sonar measurements) was interpolated onto the mesh with the bathymetric
depth being negative below and positive above the reference plane, i.e. NAP,
the Dutch reference level which is about local mean sea level. A time series
of measured water levels with 10 min time interval was used as boundary
condition on the seaward boundary (Fig. 1b). Upstream there are eight
boundaries with an imposed daily averaged discharge (Fig. 1b). The timestep
was set to 6 s. A constant viscosity model was used with a value of
10
Since bottom friction is regarded as the main sink for momentum this
parameter was used for calibration of the model. Telemac uses the
Chézy roughness formulation, with the Chézy roughness coefficient,
Storm surge generation by wind and atmospheric pressure effects were not explicitly taken into account in the model domain, assuming that wind and atmospheric pressure effects on water levels are mainly generated on the North Sea (e.g. Gaslikova et al., 2013), which was not included in the model domain. Measured storm surges are implicitly included in the imposed measured water levels at the seaward boundary of the model domain (Fig. 1b) (similarly as in Temmerman et al., 2012), and assuming that additional wind and atmospheric pressure effects on water levels within the estuarine model domain are limited because of the limited fetch length and limited water depths in such a narrow and shallow estuarine setting. Two real storm surges were included in the boundary conditions and were used for model validation. The first storm dated from 10 February 2009 and the second one from 6 December 2013. The storm return periods were calculated for the measured water levels at the mouth of the estuary according to Willems et al. (2007) and resulted in 2.4 and 20 years, respectively. Due to the extra wind upset high water levels at the seaward boundary of the model increased to 0.74 and 1.54 m, for both storms respectively, above the astronomical calculated high water levels. For both storms time series of water level measurements were available for the entire estuary and so modelled water levels could be compared with measured ones, revealing an acceptable agreement between both (see model validation in Sect. 3.1). The 2013 storm was used in the scenario simulations.
Schematic summary of the scenarios (see explanations in the text).
The colours used in the small maps indicate the subtidal area (
In this paper we focus on the effect of the largest currently remaining
intertidal area of Saeftinghe (3000 ha; Fig. 1b) on the storage of storm
water and attenuation of the tidal wave travelling through the estuary. We
differentiate between realistic (
For the realistic scenarios the measured bathymetries of 1931, 1963 and 2009 were used (see Wang and Temmerman, 2013). These scenarios are called P.R.1, P.R.2 and P.R.3 respectively, where P refers to platform elevation and R to realistic topography. The topography of 2009 was used in the calibrated and validated model and this scenario, P.R.3, will be used as a reference against which all other realistic scenarios will be compared. In 1931 this area was mostly only non-vegetated tidal flat with an average elevation of 0.85 m NAP (see small figure in Fig. 2). Through the years this area trapped a lot of sediments and a vegetated marsh started to develop (see Wang and Temmerman, 2013). In 1963 the average elevation was 1.40 m NAP and in 2009 it was 1.85 m NAP. In 2009 Saeftinghe developed into a vegetated marsh area with large non-vegetated drainage channels.
For the schematic scenarios a uniformly flat platform without drainage channels was implemented. The first three levels coincide with the average elevations of 1931 (0.85 m NAP), 1963 (1.40 m NAP) and 2009 (1.85 m NAP) respectively. A fourth schematic scenario was added with an even higher platform elevation of 2.5 m NAP. These scenarios will be called P.S.1, P.S.2, P.S.3 and P.S.4, respectively, where P stands again for platform elevation and S for schematic topography (Fig. 2). Scenario P.S.3 with a platform elevation that is equal to the average wetland elevation in the realistic scenario P.R.3 will be used as a reference against which all schematic scenarios will be compared.
Here we test the impact of the size of a wetland on the tidal propagation within the estuary. For the realistic scenarios Saeftinghe was excluded from the estuary (S.R.1; where S stands for surface area), 2000 ha in the back of the wetland were excluded by creating a dike in the bathymetry leaving about 1000 ha intertidal area (S.R.2), 1000 ha was excluded in the back leaving about 2000 ha of intertidal area (S.R.3) and the scenario with full size area of about 3000 ha (S.R.4) is actually the same as the reference (P.R.3) (Fig. 2).
For the schematic scenarios a flat platform elevation of 2 m NAP was used. The Saeftinghe area was replaced by a larger rectangle at the same location. This area was made bigger for successive scenarios: 1000 ha (S.S.1), 1500 ha (S.S.2), 2000 ha (S.S.3), 2500 ha (S.S.4), 3000 ha (S.S.5) and 3500 ha (S.S.6) (Fig. 2).
Peak water levels (
We hypothesize that intertidal storage of a certain amount of water in the mouth area of an estuary will have different effects on tidal and storm surge attenuation than the same intertidal storage volume but located further upstream. To study this effect, different scenarios were simulated with a wetland area at different locations along the length of the estuary. These scenarios are all with schematic flat wetland. The first one is with a 3000 ha large square intertidal area near the mouth of the estuary, around km 10 in the estuary (L.S.1; where L stands for location). Platform elevation was kept the same as in the schematic scenarios of surface area, i.e. 2 m NAP surface elevation. The second scenario is with the same area on the left bank around km 28 in the estuary (L.S.2) (Fig. 2). Like the Saeftinghe area these new locations are also situated next to the main channel. For the third scenario a schematized flat version of the Saeftinghe area was used (L.S.3). Further upstream the total discharge through the main channel decreases and so the size of the flooding area was also decreased. For the fifth scenario a 1500 ha large area was placed near km 85 on the right bank (L.S.5). To compare the results of this smaller sized area an extra scenario with a downscaled size of the Saeftinghe area (1500 ha) was made (L.S.4).
The hydrodynamic model was validated by simulating the whole year 2009. The model represents the measured tides reasonably well: for most tidal gauge stations along the estuary (Table 1), the average difference between all measured and modelled water levels is limited from 0.04 to 0.07 m (with a RMSE of 0.06 to 0.09 m), except for the three most inland tidal gauge stations, where the average difference is 0.12 to 0.15 m and RMSE is 0.15 to 0.18 m. Phase differences between the moments of high and low water stayed below 20 min for all stations. For storm tide attenuation the most important variable is the high water level along the estuary. Without special calibration for storm tides, the model performed fairly well in reproducing the two storm tides in the estuary, with average differences between measured and modelled high water levels between 0.01 and 0.22 m for the 2009 storm, and between 0.06 and 0.20 m for the 2013 storm (Table 1).
Average and RMSE (for the whole year 2009) of difference in
measured and modeled water level (
High water levels rise starting from the estuary mouth until km 120 along the estuary due to convergence, which results in flood wave amplification (Fig. 3). Further upstream, the effects of friction will dominate more than the convergence and the high water levels start decreasing again, resulting in flood wave damping. This trend is the same for neap, spring and storm tides.
Scenarios with realistic wetland elevation relative to reference scenario P.R.3 (given by horizontal zero line). Scenarios given for spring and storm tide (2013). The green area shows the location of the Saeftinghe wetland along the estuary.
Scenarios with schematic wetland elevations relative to reference scenario P.S.3. Scenarios given for spring tide and storm tide (2013). The green area in the figure shows the location of the Saeftinghe wetland in the model.
All scenarios are compared with the reference scenario (P.R.3 or P.S.3 for realistic or schematic scenarios respectively). This means that for each scenario the high water levels along the estuary are calculated relative to the high water levels of the reference scenario (Scenario – Reference). Negative values mean that a scenario results in lower high water levels compared to the reference scenario, and vice versa. The coloured area in the figures indicates the location of the wetland depending on its elevation.
For the scenarios with different realistic wetland elevations, Fig. 4 shows
the differences in high water levels relative to the reference scenario for
a spring tide and the 2013 storm. For a spring tide, the lower the wetland
elevation, the more the water levels are attenuated along the estuary with a
maximum difference of almost 0.2 m upstream of the wetland area (Fig. 4). For
the storm tide, the effects are more complex and differences between the
different scenarios are much smaller than for the spring tide. Upstream
differences are small (
Water surface slopes for scenarios with different schematic
wetland elevations (P.S.1, P.S.2, P.S.3 and P.S.4) for a spring tide between
middle of the main estuarine channel and the wetland edge (distance was
about 1000 m). Positive slopes occur when the water floods from the main
channel onto the wetland and negative slopes when the water drains again
from the wetland towards the main channel. As a reference, the water level
fluctuation in the main estuarine channel (thick black line) is given with
values on the right
The same trend as with the realistic scenarios for a spring tide is seen in
Fig. 5 for the schematic scenarios: the lower the wetland elevation, the
lower the water levels become along the estuary. For the lowest elevation
scenario (P.S.1) the maximum attenuation of high water levels occurs just
downstream of the wetland and amounts to 0.18 m (as compared to the
reference scenario P.S.3 with wetland elevation
For the storm tide the effects are again more complex and the difference between all scenarios is again much smaller compared to the spring tide. Scenarios P.S.1 and P.S.2 with lower wetland areas have higher water levels (up to 0.08 m) especially just downstream of the wetland and upstream along most of the estuary. For the highest elevation (P.S.4) the high water levels first increase (up to 0.07 m) downstream of the wetland. Near the wetland area the high water levels drop 0.14 m and stay low upstream of the wetland compared to the other scenarios. Overall, upstream of the wetland the storm tide is more attenuated when the marsh has a higher elevation, which is the opposite effect to the spring tide, which is more attenuated when the wetland has a lower elevation.
Water surface slopes for scenarios P.S.1, P.S.2, P.S.3 and P.S.4 for
the storm tide of 1994 between middle of the main estuarine channel and the
wetland edge (distance was about 1000 m). Positive slopes occur when the
water floods from the main channel onto the wetland and negative slopes when
the water drains again from the wetland towards the main channel. As a
reference, the water level fluctuation in the main estuarine channel (thick
black line) is given with values on the right
Scenarios for different realistic surface areas relative to reference scenario S.R.4. Scenarios given for a spring tide and storm tide (2013). The green area in the figure shows the location of the Saeftinghe wetland in the model.
Scenarios for different schematic surface areas relative to reference scenario P.S.3. Scenarios given for a spring tide and the storm tide of 2013. For spring tide all scenarios gave nearly the same result. The green area in the figure shows the location of the Saeftinghe wetland in the model.
Water surface slope between main estuarine channel and wetland edge (distance between them was about 1000 m) simulated for the storm tide of 2013 for different schematic scenarios with different wetland surface areas.
In order to further examine the differences in attenuation between the
spring and storm tide, we analysed the time evolution of the water surface
slope between the main estuarine channel and the adjacent edge of the
wetland, which is the driving factor for water fluxes from the main
estuarine channel onto the wetland (a higher positive slope indicates a
faster water flux from the channel towards the wetland and vice versa). For
the schematic scenarios the slopes for the spring and storm tide are plotted
as a function of time in Figs. 6 and 7, respectively. For a spring tide,
Fig. 5 shows that the lowest platform elevation (P.S.1) resulted in the strongest
tidal attenuation along the estuary. Looking at the slopes in Fig. 6, the
scenario with the lowest platform elevation (P.S.1) is the first for which
water starts to flow onto the wetland (
For the storm tide of 2013 the highest platform elevation (P.S.4) results in the strongest attenuation of the storm tidal wave at the location of the wetland and upstream along the estuary (Fig. 5). In Fig. 7 the slopes for the different schematic platform elevations are plotted for the storm of 2013. The lowest platform (P.S.1) starts to store water first, but the slope reaches a lower maximum value and it decreases faster to zero than for the higher platform scenarios (P.S.2, P.S.3, P.S.4; see Fig. 7). The higher the platform the later it starts storing water, but the higher the maximum slope and the longer it can keep storing water (Fig. 7).
The wetland surface area has a clear effect on tidal and storm surge
attenuation along the estuary: the larger the wetland surface area, the
larger the attenuation of the tidal wave along the estuary, both for the
spring and storm tide (Fig. 8). The effect of the wetland surface area on
attenuation is also bigger for the storm tide than for the spring tide. At
spring tide the loss of 1000 ha (
For a spring tide all scenarios gave similar results (Fig. 9). An increase
in high water levels starts downstream of the large wetland area growing
towards an increase of 0.17 m upstream and it starts decreasing again from
km 85 on. For the storm tide of 2013 the smallest surface area scenario (S.S.1)
results in the highest increase in high water levels (up to 0.18 m)
compared to the reference scenario. The second scenario (S.S.2) also has
higher water levels (but only up to 0.12 m) compared to the reference. The
scenario with surface area of 2000 ha (
Scenarios with different wetland (with surface area
Figure 10 presents the water surface slopes between the main estuarine channel
and the wetland for the schematic surface area scenarios for the storm tide (2013).
Because the wetland elevation is the same for all scenarios, the
water starts flooding onto the wetland (
Scenarios with different wetland (with surface area
For the scenarios with a 3000 ha wetland area (scenarios L.S.1, L.S.2 and L.S.3)
the attenuation of the tidal wave upstream of the wetland area is
stronger if the wetland is located more upstream along the estuary, both for
a spring and storm tide (Fig. 11). The difference in high water levels
between L.S.1 (with wetland at 15–20 km from estuary mouth) and L.S.2 (at
26–32 km) is small (
For the two scenarios with a 1500 ha wetland area, the scenario L.S.5 is
compared with L.S.4, used as reference (
For locations more downstream in the estuary, the total flood volume that passes through a cross-section of the estuary is larger (illustrated for the storm tide in Table 2, fourth column). Consequently, wetlands of the same size and elevation, but located more downstream in the estuary, can store a smaller percentage of the total flood volume that passes through the estuary at that location (Table 2, sixth column). Additionally, the difference in water volume stored on the wetland (Table 2, fifth column) between the first three scenarios (with equal size and elevation but different locations along the estuary) is explained by the amplification of the high water levels in upstream direction along the estuary (Fig. 3). This results in longer periods of flooding and thus in larger volumes stored on the wetlands upstream of the estuary. Table 2 shows that the wetland of 1500 ha located most upstream (L.S.5) is most effective in storing the largest proportion of water that passes locally trough the estuarine channel.
For the schematic scenarios with different wetland locations, the storm tidal flood volumes (storm 2013) that pass through the main estuarine channel next to the wetland and the volumes stored on top of the wetland are given. The last column shows the flood volumes stored onto the wetland as the percentage of the volume that passes through the main channel.
The role of wetlands for attenuation of storm surges is over the last few
years demonstrated by an increasing number of field-based studies
(e.g. Costanza et al., 2008; Das and Vincent, 2009; Kraus et al., 2009; Gedan et
al., 2011) and hydrodynamic modelling studies (e.g. Wamsley et al., 2010;
Temmerman et al., 2012; Zhang et al., 2012; Sheng et al., 2012; Liu et al.,
2013; Hu et al., 2015). All these studies focussed on the attenuation of
storm surges that propagate through coastal wetlands, which we call here
In this paper we focus on the effect of different wetland characteristics, including wetland elevation, wetland surface area, and wetland location along the estuary, on the attenuation of spring and storm tides along the length of the estuary. In summary, for a spring tide lower wetland elevations result in more attenuation of high water levels along the estuary, while for a higher storm tide higher wetland elevations provide more attenuation compared to lower wetland elevations. For spring and storm tide a larger wetland surface area results in a better attenuation along the estuary. The schematic scenarios suggest that the increase in attenuation is non-linearly related to the increase in size and that for wetlands larger than a certain threshold there is no significant increase in attenuation anymore with increasing size. Finally a wetland of the same size and elevation, but located more upstream in the estuary, can store a larger proportion of the local flood volume and therefore has a larger attenuating effect on upstream high water levels.
For both the realistic and the schematic scenarios at spring tide large
differences in high water levels along the estuary are shown in Figs. 4 and 5.
For spring tide the lower the elevation and the lower the high water levels,
the better the
In contrast with a spring tide, a storm tide brings much more water to our
wetland in our model setup. Now, depending on the elevation the wetland can
be completely filled with water. A lower elevation causes an early start of
the water storage (Fig. 7), but as the wetland gets completely filled with
water the slope starts to decrease and falls back to zero. This happens much
earlier than with higher elevated wetlands. For higher elevations, the
wetland flooding starts later. The wetland's storage capacity is not filled
completely and thus higher water surface slopes are reached during flooding
of a higher elevated wetland. The available storage capacity keeps on
filling to a time much closer to the high water level in the estuarine
channel (Fig. 7). Hence there is an optimum wetland elevation for different
tides depending on the maximum high water level. For the storm tide we
modeled, if the wetland elevation is too low, the storage capacity is
completely used, resulting in lower slopes between channel and wetland and
slower filling rates. If the elevation is too high, filling the storage area
starts too late and the full capacity is not used. In that case the slopes
are high and the rate at which water flows onto the wetland is high, but the
duration is too short, leading to a sub-optimal use of the storage capacity.
Loder et al. (2009) stated that lower wetland elevation facilitated the
landward propagation of coastal storm surges. Their focus is on
The realistic scenarios (S.R.1, S.R.2, S.R.3 and S.R.4) showed that for
increasing wetland surface areas the
The realistic scenarios (S.R.2, S.R.3 and S.R.4) have large drainage
channels, which are wider and deeper near their connection with the main
estuarine channel. This means that cutting of a large part at the landward
side of the wetland also results in a decrease of the average elevation of
the wetland site. So for these scenarios a change in surface area
also results in a change in average elevation (effects discussed in previous
section). The schematic scenarios (S.S.1 to S.S.6) with a flat wetland
topography show only the effect of the surface area. It is clear that for a
certain storm tide with specific maximum high water level there is a
limitation on the time that water can flow onto the tidal wetland, and that
this time of flooding continues longer as the wetland surface area is
larger. Figure 10 shows that for smaller wetland areas (like S.S.1) the water
surface slope decreases to zero far before the high water level in the main
estuarine channel is reached. This is because small wetland areas are more
rapidly filled with flood water, hence decreasing the slope and discharge
towards the wetland before the peak water level in the main estuarine
channels is reached. Accordingly, smaller wetland areas (like S.S.1) are less
effective in attenuating peak water levels along the estuary during a storm
tide (Fig. 9). If the wetland area is large enough (like S.S.2, S.S.3 and
S.S.4) the water surface slope and hence the discharge towards the wetland
remain high, until the water level in the main estuarine channel will fall
again, decreasing the slope and thus the discharge. Our results imply that
the higher the incoming storm tide, the longer the water surface slope
towards the wetland can stay positive and the larger the wetland surface
area needs to be to maintain this slope and to attenuate the storm tide most
effectively. In this respect if we keep the width (i.e. the length of
connection between wetland and estuarine channel) of the wetland fixed, like
in our scenarios, increasing the surface area means increasing the depth
(i.e. the length, perpendicular to the estuarine channel and towards the
back of the wetland). Increasing the depth means that water needs a longer
time to reach the back of the wetland and so to use all the available
storage capacity. If we keep the depth fixed and increase the surface area,
it means that the wetland area is widened and this will give different
results for the
An increase in intertidal surface area can have large effects on hydrodynamics up- and downstream of the affected intertidal area by increasing the tidal prism. Townend and Pethick (2002) showed this with a conceptual and idealized model. They increased the tidal prism by managed retreat of dikes and levees and show that this will increase the size of the estuarine channels to accommodate for the increased water demand by the added intertidal area (see also Fries et al., 2014), possibly changing the characteristics of the system (from flood to ebb dominancy). Townend and Pethick conclude that increasing the intertidal surface area will decrease tidal water levels and decrease flood risk. We come to the same conclusion in this paper and share their concern that interfering in the tidal prism will possibly change the morphodynamics of the estuary. By choosing elevations of 2 m NAP for our schematic scenarios interference with a normal spring-neap tidal cycle will be minimized, i.e. it will not affect the tidal prism for most tides (like the spring tide scenarios in Fig. 9). Therefore the effects on the morphodynamics of the estuary (which is governed by regular neap to spring tides instead of exceptional storm tides) will be minimal.
Surface area was also found to be an important parameter in previous studies on storm surge attenuation within foreshore wetlands. The larger these areas, the more high water levels reduce during storm surges within and behind coastal wetlands (e.g. Loder et al., 2009; Wamsley et al., 2010). In these studies the surface area is more a measure for the length that friction can act on the flood wave, whereas in our study the surface area is a parameter determining the total storage capacity of a wetland.
Table 2 shows a decrease in total storm flood volumes in the main estuarine channel in upstream direction. As a consequence, the more upstream the less flood water needs to be stored to have a similar effect on the tidal wave attenuation along the estuary. Furthermore, the high water levels are amplified in upstream direction by the funnel shape of the estuary (until the location where the effect of friction becomes dominant over convergence; Fig. 3) so that more upstream located wetlands with the same elevation will be flooded by higher water levels and thus during a longer time. As a result, the degree of attenuation of a spring tide is different for the wetlands with different locations (Figs. 11 and 12). Table 2 also shows that more upstream smaller surface areas can store higher percentages of the total storm flood flow at the specific locations. So the higher the percentage of flood water stored onto the wetland the larger the attenuation. The results show further a spatial extent of high water level reduction along the estuary. Figure 12 shows clearly for the most upstream wetland at storm tide that the effect of attenuation is situated mainly upstream. This means that when we want to have a specific flood wave attenuation at a specific reach along the estuary, the location where we want to add extra wetland surface area is very important. If the wetland is located too far downstream from this location it will need to have a larger surface area to create the same amount of attenuation than a wetland that is located near our area of interest.
In real estuaries all these parameters interfere with each other in a complex way and when discussing the impact of one of them, we cannot ignore the others. For example if the wetland elevation is increased the water will enter the wetland at a later time and so the duration of flow onto the wetland will be shorter and thus less surface area (storage capacity) is needed. Further, more characteristics are important, like the connection between wetland and the estuarine channel. The longer this connection length, the more water can access the wetland. This was not taken into account in this paper, but will be a topic for further research.
This paper presented the results of a hydrodynamic model study on the
along-estuary propagation of a spring and storm tide and the attenuation effect of
different geometric characteristics of a tidal wetland, i.e. the wetland
elevation, surface area, and location along the estuary. Based on our model
results the following conclusions can be drawn:
For a spring tide lower wetland elevations result in more attenuation of
high water levels along the estuary, while for a higher storm tide higher
elevations provide more attenuation compared to lower wetland elevations. The
wetland elevation relative to the height of the tidal wave determines the start
and duration of flow from the estuarine channel to the wetland. We found that
there is an optimal wetland elevation relative to the height of the tidal wave,
for which the water fluxes from the estuarine channel reach a maximal rate and
duration so that the storage capacity of the wetland is most optimally used. Larger wetland surface areas result in more flood wave attenuation along
the estuary. For a specific wetland elevation and location along the estuary,
the amount of flood wave attenuation does not further increase with increasing
wetland area above a certain wetland depth, i.e. the distance from the estuarine
channel, because time is lacking to fill and use the storage capacity of that
additional space. The more downstream in an estuary, the larger the flood volumes in the
main estuarine channel will be and the larger the wetland surface area will
have to be to have a significant effect on upstream tidal and storm tide
attenuation along the estuary. A wetland of the same size will store a smaller
percentage of the total flood volume at a more downstream location, than at a
more upstream location, and therefore will have a smaller effect on attenuation
of high water levels upstream along the estuary.
In conclusion, this paper provides a first insight in the separate effects
of geometric properties of estuarine tidal wetland on tidal and storm tide
attenuation along an estuary. Further research is necessary for example to
account for the combined effects of wetland geometric properties, the
effects of multiple wetlands located at different positions along an
estuary, and effects of storm properties, such as duration and intensity, on
flood attenuation along an estuary, in order to support a better
understanding and wider implementation of nature-based adaptation to
increasing flood risks associated with storms (Temmerman et al., 2013).
The authors like to thank the Antwerp Port Authority and the Dehousse scholarship of the Antwerp University for the financial support. We like to acknowledge the technical support of the Flanders Hydraulics Research institute. We like to thank Rijkswaterstaat for all measurements and data sharing. Finally we would like to thank the developers of the Telemac system and Blue Kenue software. Edited by: J. Brown Reviewed by: C. Ibanez and another anonymous referee