Next Article in Journal
Sea Surface Temperature Retrieval from the First Korean Geostationary Satellite COMS Data: Validation and Error Assessment
Next Article in Special Issue
sUAS-Based Remote Sensing of River Discharge Using Thermal Particle Image Velocimetry and Bathymetric Lidar
Previous Article in Journal
Assessment of Satellite-Derived Precipitation Products for the Beijing Region
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Empirical Assessment Tool for Bathymetry, Flow Velocity and Salinity in Estuaries Based on Tidal Amplitude and Remotely-Sensed Imagery

by
Jasper R. F. W. Leuven
*,
Steye L. Verhoeve
,
Wout M. Van Dijk
,
Sanja Selaković
and
Maarten G. Kleinhans
Faculty of Geosciences, Utrecht University, Princetonlaan 8A, 3584 CB Utrecht, The Netherlands
*
Author to whom correspondence should be addressed.
Remote Sens. 2018, 10(12), 1915; https://doi.org/10.3390/rs10121915
Submission received: 8 November 2018 / Revised: 22 November 2018 / Accepted: 23 November 2018 / Published: 30 November 2018

Abstract

:
Hydromorphological data for many estuaries worldwide is scarce and usually limited to offshore tidal amplitude and remotely-sensed imagery. In many projects, information about morphology and intertidal area is needed to assess the effects of human interventions and rising sea-level on the natural depth distribution and on changing habitats. Habitat area depends on the spatial pattern of intertidal area, inundation time, peak flow velocities and salinity. While numerical models can reproduce these spatial patterns fairly well, their data need and computational costs are high and for each case a new model must be developed. Here, we present a Python tool that includes a comprehensive set of relations that predicts the hydrodynamics, bed elevation and the patterns of channels and bars in mere seconds. Predictions are based on a combination of empirical relations derived from natural estuaries, including a novel predictor for cross-sectional depth distributions, which is dependent on the along-channel width profile. Flow velocity, an important habitat characteristic, is calculated with a new correlation between depth below high water level and peak tidal flow velocity, which was based on spatial numerical modelling. Salinity is calculated from estuarine geometry and flow conditions. The tool only requires an along-channel width profile and tidal amplitude, making it useful for quick assessments, for example of potential habitat in ecology, when only remotely-sensed imagery is available.

Graphical Abstract

1. Introduction

Estuaries are characterised by fresh water inflow at the landward boundary and an open connection to the sea. Within these boundaries, tidal flows form dynamic patterns of channels and bars (Figure 1a). Many of these systems are managed to balance the needs of flood safety, access to harbours and ecological quality. The main channels in large systems such as the Western Scheldt, Elbe and Yangtze are dredged for access to harbours [1], while intertidal area on bars and estuarine shorelines forms valuable ecological habitat [2]. Depth, inundation time, flow velocity and salinity are the prime abiotic factors that affect living organisms in estuarine environments [2,3,4]. However, the general lack of information on these abiotic factors limits the prediction and sustainable management of habitat area in most estuaries [5]. Additionally, the effects of human interventions and rising sea-level on the hydrodynamic conditions, equilibrium morphology and, therefore, also on ecology are largely unknown. Unfortunately, data to study these effects are generally only available for a few economically exploited estuaries that are already under great pressure by human influence. Data are lacking for most other systems around the world. Here, we explore to what degree it is possible to globally estimate estuarine characteristics, such as bed elevation, inundation duration, flow velocity, and salinity based on limited but widely publicly available data: remotely-sensed imagery and tidal range at the estuary mouth. To aid application and further investigation, we make the toolbox available.
Currently, several alternative methods are available to study the equilibrium morphology and hydrodynamic conditions of estuaries: in situ measurements, numerical modelling [6,7,8], idealised form models [9], stability analyses [10,11,12] and stability relations [13,14,15,16]. Only numerical models, for example Delft3D, can predict the spatial patterns of bed elevation, flow velocity and salinity with considerable accuracy at the spatial detail that is required for the purpose of habitat estimation. However, morphological models are computationally expensive and require calibration and detailed measurement of boundary and initial conditions. In situ measurements are an alternative, but they are expensive and equally time-consuming. For example, Rijkswaterstaat typically covers 2 km 2 /day in the Western Scheldt. While linear stability analysis may predict typical estuarine bar patterns, it is not straightforwardly applicable to fully developed bars and estuaries with irregular planforms.
For converging ideal estuaries, predictive relations are already available between width-averaged bed elevation and hydrodynamic conditions along the system [9,17,18] because of the imposed constraint that the tidal damping by friction is balanced by the gain in energy by convergence. While these concepts are useful when applied to end-members, they leave the more common estuaries with irregular planforms unexplored. We recently found that the concept of an ideal estuary can be used to predict bed elevation and bar patterns in estuaries with irregular planform shapes as a function of the deviation from the ideal shape [19,20]. However, we still lack empirical or theoretical relations that describe flow velocity in these non-ideal systems, while this is an essential input variable for salinity predictions [21,22,23] and habitat suitability [24,25]. Based on observations in the Western Scheldt, we hypothesise that a relation exists between the depth and tidal flow velocity amplitude (Figure 1a,b). Here, we derive a predictor for flow velocity as a function of bed elevation to complement the available relations for bed elevations, bar patterns and salinity in estuaries with irregular planforms.
Here, existing empirical relations and theoretical concepts were selected and combined in a fast-to-apply tool that predicts equilibrium morphology and hydrodynamic conditions. To assess the applicability, we test it against a range of estuaries and other end-member systems, characterised by varying degrees of tidal and river influence [9,26] as well as varying degrees of adaptation to boundary conditions [9,27] and varying degrees of bar presence (Figure 2). These systems are a river-dominated delta branch [26], a tide-dominated delta branch (an ideal estuary cf. Savenije [18]), an alluvial estuary filled with tidal bars [19,26], a relatively wide estuarine valley partially filled with bars [26,28], and a tidal basin lacking fluvial input [29,30,31] (Figure 2). In the discussion, we show applications for management of estuaries, palaeogeographic reconstructions and quantifying habitat area based on species-specific preferences.

2. Methods

Below, first, the general approach and minimum data required to run the tool are described. Then, the general structure and functions used in the tool are presented. Finally, it is described how tool output was compared with measured data and numerical modelling for two estuaries filled with bars and three end-member systems (Figure 2), including the Western Scheldt and Columbia River. A general flowchart of the tool is given in Figure 3 and a list of symbols used is given in Appendix A Table A1. The full Python code of the tool and instructions on how to use it are available in the supplementary material and accessible on Github. The tool output comprises data exported in spreadsheet format and maps of bed elevation, inundation duration, peak flow velocity and salinity throughout the estuary.
General assumptions underlying the entire method are as follows. There are no significant geological constraints on the depth along the estuary. At the landward and seaward boundary, we assume hydraulic geometry. Moreover, there is no significant effect of the offshore wave climate on the morphology of the mouth and bed elevations in the estuary. The along-channel tidal range is simplified as an along-channel profile that is constant or linearly increasing or decreasing. For calculations of inundation duration, a sinusoidal M2 tide is assumed.

2.1. Tool Approach

In general, the level of predictable detail depends on the available data. Data that are generally available for estuaries from all over the world are remotely-sensed imagery and a typical tidal range at the mouth, which means that our predictive relations should be based on and derivable from this information only. From remotely-sensed imagery, the along-channel width profile can be obtained (Section 2.2). It was previously found that the along-channel width profile correlates with the cross-sectional hypsometry [20]. The deviation from a perfectly converging planform translates in the occurrence of intertidal bars and therefore also in the shape of the hypsometric curve [33,34,35]. However, to scale dimensionless hypsometry to vertical bed elevation distributions, an along-channel depth profile is required. Estuary depth generally decreases linearly to near-linearly from the estuary mouth to the upstream tidal limit citepSavenije2015, enabling channel-depth estimation at any location in the estuary when the depth at the mouth of the estuary and upstream river are known. For the landward boundary, hydraulic geometry relations for rivers predict channel depth from width or bankfull discharge [36,37], while the geometry of the estuary mouth is related to local tidal prism [13,15,17,38,39]. The typical tidal range at the mouth in combination with the along-channel width profile gives an estimate of the local tidal prism. This means that bed elevations and inundation duration are predictable from only the along-channel width profile and tidal amplitude profile, which is the minimum input data required to run the tool.
In order to be able to use additional data when available, we allow some variables to override estimations by the above relations (Figure 3). Specifically, measured maximum mouth depth and maximum depth at the upstream boundary make their estimation based on tidal prism and hydraulic geometry redundant. A specified depth and width at the upstream boundary implicitly impose a user defined hydraulic geometry. Similarly, measured fresh water discharge may override the estimation of discharge based on channel width, which replaces the high uncertainty of the empirical relation with the lower measurement uncertainty. Finally, salinity at the mouth and landward boundary may be specified instead of the default values of 35 ppt at the mouth and 0 ppt at the landward boundary.

2.2. Bed Elevation Prediction

From remotely-sensed imagery, the first step is to obtain the along-channel width profile. Multiple options are available ranging from manually measuring the width at equally spaced transects to software tools that calculate river widths from remotely-sensed imagery (e.g., RivWidth [40,41]). Here, the following approach was used (following Leuven et al. [19]): first, we digitised the estuary planform in GIS software. Second, the polygon was imported in GIS software and the channel centreline was automatically extracted (e.g., Polygon To Centerline, ArcGIS function). Last, transects were drawn equally spaced on the centreline (e.g., Generate Transects Along Lines, ArcGIS function) and cropped with the planform polygon. The length of successive transects resulted in the along-channel width profile.

2.2.1. Estuary Shape Predicts Cross-Sectional Hypsometry

Estuaries typically narrow from the tidal inlet to the upstream river because the tidal wave generally dampens out in landward direction such that the combined discharge, or flow energy, of the river and the tides reduces. The estuary convergence is often characterised by a convergence length of width [18,42,43]. In ideal estuaries, the energy loss due to friction is balanced by channel convergence such that the tidal energy per unit of width remains constant up to the estuary head [17,44,45]. However, natural estuaries adapted in varying degrees to their initial and boundary conditions [9,27], resulting in planform shapes that deviate from the converging shape expected according to the theory for the ideal estuary state. The deviation in width from a converging shape was found to be a fairly accurate predictor of bar pattern and cross-sectional bed elevation distributions [19,20].
The convergence length ( L W ) is the length over which the channel width reduces by a factor 2.72 (e) and is calculated as:
L W = s ln W r W m
where W m is the width of the estuary mouth, W r is width of the river at the most landward boundary of the estuary and s is the along-channel distance measured between the seaward and landward boundary. We followed the same guidelines on the selection of the landward and seaward boundary as in Leuven et al. [20], given on p. 766.
The ideal width ( W i d e a l ( x ) ) and excess width ( W e x c e s s ( x ) ) are assumed to be related to this convergence length and the local width ( W ( x ) ) of the estuary:
W i d e a l ( x ) = W m e x L W
W e x c e s s ( x ) = W ( x ) W i d e a l ( x )
Here, x is the streamwise coordinate measured from the mouth along a centreline and y is the coordinate perpendicular to the centreline.
Along-channel variations in excess width translate into bed elevation with cross-sectional hypsometric curves. Hypsometric curves describe bed elevations with a cumulative profile. Multiple empirical relations have been proposed for the hypsometric shape of (partially) submerged bodies [34,35,46,47] and terrestrial landscapes [33]. The original formulation of Strahler [33] (Equation (5)) appeared capable to describe hypsometries that occur in estuaries [20]. Two parameters in the Strahler [33] equation can be tuned to modify the hypsometric shape: r sets the slope of curvature at the inflection point and z sets the concavity of the function, with lower values representing a more convex cross-sectional bed profile and higher values representing a more concave profile.
Reasonably accurate predictions are obtained in estuaries with r set to a constant value of 0.5, while the along-channel variation of best fitting z strongly depends on the along-channel variation in width [20]. Moreover, at locations where the estuary is much wider than expected from an ideal shape, bars are typically more abundant, the hypsometry is more convex and z-values are higher. The along-channel variation in concavity (z) is given by:
z ( x ) = 1 . 4 W i d e a l ( x ) W ( x ) 1.2
The resulting predictions are accurate within a factor 2 of the measured value for intertidal and subtidal area in multiple estuaries of various size and character. This means that it is possible to predict the hypsometric shape per cross-section based on the estuary shape and, in this study, its predictive capacity for end-member cases is assessed as well.
The along-channel variation in concavity of the hypsometric shape ( z ( x ) ) can subsequently be used in the general hypsometric curve, formulated by Strahler [33] as:
h z ( x , y ) = r r 1 z ( x ) 1 ( 1 r ) y + r 1 z ( x )
in which h z ( x , y ) is the bed elevation above which proportion y of the width profile occurs. Alternatively said, h z ( x , y ) is the proportion of total section height and y the proportion of section width. r is set to a constant value of 0.5 and z ( x ) is given by Equation (4).
The bed elevation ( h z ( x , y ) ) and the width fraction (y) in Equation (5) are dimensionless, where h z ( x , y ) is normalised between the local high water level (HWL) and the maximum estuary depth for that cross-section. Values for y are dimensionalised with the local estuary width. Thus, to transform the predicted hypsometry to dimensional depth distributions, an along-channel profile of the maximum depth is required, which is estimated from the channel depth at the mouth and the depth at the landward boundary as explained in the next subsection.

2.2.2. Tidal Prism and River Discharge Predict Channel Geometry at the Mouth and Landward River

To calculate dimensional hypsometric profiles, an along-channel depth profile is required. Since the along-channel depth profile often shows a linear or almost linear profile [19,42], we fitted a linear along-channel profile on the maximum depth at the mouth of the estuary and the maximum depth at the landward river. When these values are known, the user can specify them, which overrides the estimation. When unspecified, the maximum and average depth at the mouth and at the upstream river are estimated with hydraulic geometry, using respectively the tidal prism and river discharge.
Hydraulic geometry is an empirical construct and should be used with care. However, it is often the best estimate in the absence of observed data. The resulting predictions are dependent on the specific region or river type for which they were developed. Here, we implemented, and used, the equation Hey and Thorne [37] as an example. Nevertheless, users of the tool can specify a measured channel depth and river discharge as input when available because it overrides the calculation with Equations (6) and (7) and therefore implicitly impose a hydraulic geometry. Alternatively, the user can calculate these values with a hydraulic geometry equation suitable for the type of river considered or implement their own equation in the source code.
When unknown, river discharge is estimated from the channel width at the landward boundary with hydraulic geometry. Here, we used Hey and Thorne [37], rewritten as:
Q b = W r 3.67 1 / 0.45
in which Q b is the bankfull discharge and W r is the width at the upstream river, measured at a location with almost zero tidal influence.
Subsequently, the width-averaged depth at the landward boundary is estimated from bankfull discharge as:
h r ¯ = 0.33 Q b 0.35
The maximum depth at the landward boundary is a function of the average depth and the geometric shape of the channel. Therefore, the following equation was adopted:
h m a x , r = s r h r ¯
in which s r is the shape factor of the cross-section at the landward river, where 1 is a perfect rectangular and 2 a V-shaped cross-section. If unknown, a value of 1.85 is assumed for the landward boundary, based on geometry of five estuary mouths (Appendix A Figure A1).
The average depth at the mouth ( h m ¯ ) is calculated as the cross-sectional area at the mouth ( A c s a , m , Equation (10)) divided by the local width ( W m ), in which the cross-sectional area is estimated with the tidal prism. The local tidal prism is given as:
P ( x ) = x x + E 2 W ( x ) a ( x ) + Q b 4 t
in which P ( x ) is the local tidal prism resulting from the tidal amplitude ( a ( x ) ) and surface area upstream of location x. The contribution of river discharge ( Q r ) to the local tidal prism is approximated as the bankfull discharge ( Q b ) divided by 4. E is the tidal excursion length, which is the distance a water particle travels over half a tidal cycle. Since E is dependent on flow velocities, which are unknown at this point in the assessment tool, E has to be estimated at this point. The tidal excursion length is estimated by multiplying a typical tidal flow velocity with t, which is the duration of half a tidal period. According to Savenije [18], the order of magnitude of the tidal excursion length can be estimated even in ungauged estuaries because the peak velocity amplitude in alluvial estuaries from all over the world is remarkably similar, being approximately 1 ms 1 . Here, we use this typical flow velocity to obtain an estimate of the tidal excursion length, which will be evaluated in the results. In Leuven et al. [19], the upper limit for summation of the tidal prism was the upstream estuary boundary, but given that the input of the tool may comprise estuaries much longer than the tidal excursion length, the summation is here limited to the tidal excursion length.
Subsequently, the cross-sectional area at the mouth is related to the tidal prism. Earlier approaches derived these relations for tidal inlets or specific tidal systems [13,15,38,39,48,49,50]. Recently, we derived an empirical relation for multiple cross-sections within 35 estuaries (supplementary material of [19]). This equation is used to estimate the cross-sectional area at the mouth:
A c s a ( x ) = 0.13 × 10 3 P ( x )
in which A c s a ( x ) is the along-channel cross-sectional area.
Now, the average depth at the mouth is calculated as:
h m ¯ = A c s a , m W m
and converted to a maximum depth at the mouth with
h m a x , m = s m h m ¯
in which s m is the shape factor of the cross-section at mouth. If unknown, a value of 1.65 is suggested (Appendix A Figure A1).

2.3. Inundation Duration Prediction

Dimensional depth distributions directly translate in inundation durations assuming a simple harmonic tide. Dimensional depth profiles are obtained from scaling the cross-sectional hypsometric profiles with the along-channel profile of maximum depth. Bed elevations below the low water level (below −1 times the amplitude) are always submerged and are therefore assigned a relative inundation duration of 1. For the intertidal zone, the inundation duration was calculated as:
I ( x , y ) = 0.5 0.5 sin 0.5 π a ( x ) h z ( x , y )

2.4. Bar Pattern Prediction: Bar Width and Braiding Index

The along-channel variation in bar pattern correlates with along-channel variations in width. Moreover, the summed width of bars is equal to the excess width and thus Braiding Index ( B I ( x ) ) can subsequently be calculated by dividing the excess width ( W e x c e s s ( x ) ) by the predicted bar width ( w b a r ( x ) ) [19]:
w b a r ( x ) = 0.39 W ( x ) 0.92
w b a r s s u m m e d ( x ) = W e x c e s s ( x )
B I ( x ) = w b a r s s u m m e d ( x ) w b a r ( x )

2.5. Flow Velocity Prediction

A relation for typical flow conditions as a function of estuarine geometry was still lacking. Therefore, we explored this relation in a calibrated hydrodynamic model of Rijkswaterstaat (SCALWEST, [32]) for the bathymetry of the Western Scheldt in 2009. We hypothesised that a relation would exist between bed elevation and typical flow velocity (Figure 1b). Peak tidal flow velocity ( u p e a k ) indeed showed a good correlation ( R 2 = 0.73 ) with depth below high water level (Figure 4):
u p e a k ( x , y ) = 0.10 + 0.91 10 log ( h ( x , y ) )
Typical average flow velocity over half a tidal cycle was estimated depending on the deviation from the maximum along-channel tidal amplitude as follows (Appendix A Figure A5):
u a v g ( x , y ) = u p e a k ( x , y ) 1 a ( x ) max ( a ( x ) ) 2 π

2.6. Salinity Prediction

Salinity along the estuary depends on the freshwater discharge at the upstream boundary and the degree of mixing by the tides. Three along-channel salinity predictors were implemented in the assessment tool [21,22,23] and the quality of these predictors for the cases studied here will be described in the results.
Savenije [21] proposed a one-dimensional model that uses predictive equations for mixing and dispersion, such that it can be applied to systems for which limited information is available. The resulting model has been tested on a set of 15 estuaries. Gisen et al. [23] revised the predictive equation of Savenije [21], calibrated the model on 20 estuaries and validated it with another 10 estuaries. Brockway et al. [22] proposes a solution for advection-diffusion equations in converging estuaries and implements an empirical relation for longitudinal mixing based on data from the Incomati Estuary (Mozambique).
The output maps in the tool are based on the averaged prediction of Savenije [21] and Gisen et al. [23]. Nevertheless, the output for Brockway et al. [22] is also provided. The along-channel salinity profile of Brockway et al. [22] is typically exponentially decreasing to a salinity of 0, while the other two are generally more s-shaped with a strong gradient in the middle of the curve. Taking the average of all three predictors would thus not make sense because their characteristic shape would smoothen out. First, the latter is described because it has a simpler solution and less constitutive variables.
The salinity profile of Brockway et al. [22] ( S b ( x ) ) is given as:
S b ( x ) = S m exp Q r L A c s a 1 K A c s a , m e ( L A c s a 1 x ) 1
in which S m is the salinity at the mouth of the estuary (35 ppt), A c s a , m is the cross-sectional area at the mouth, L A c s a is the convergence length of the cross-sectional area and K is an empirical equation for the longitudinal mixing coefficient based on average river discharge ( Q r ) and tidal range ( a ( x ) ), being K = 0.28 Q r + 13 a ( x ) ¯ .
The salinity profile of Savenije [21] ( S s ( x ) ) is given as:
S s ( x ) = ( S m S r ) D ( x ) D m 1 / K + S r
in which S r is the salinity at the fresh water river (0 ppt), K is the longitudinal mixing (Van der Burgh) coefficient (typically 0.07), D ( x ) is the dispersion coefficient and D m is the dispersion at the mouth of the estuary (typically 7000 m 2 s 1 ). The ratio between the along-channel dispersion and dispersion at the mouth is given in Savenije [21] as:
D ( x ) D m = 1 K L A c s a Q r D m A c s a , m e x L A c s a 1
with
D m = 220 40 h m ¯ L A c s a N r max ( u p e a k ( m o u t h , y ) ) E
in which
N r = Q r t / P ( m ) ρ w max ( u p e a k ( m o u t h , y ) ) 2 / Δ ρ g h m ¯
where ρ w is the density of water, Δ ρ is the density difference between salt and fresh water and g is the gravitational acceleration (9.81 ms 2 ).
The salinity profile of Gisen et al. [23] ( S g ( x ) ) is equal to the profile of Savenije [21] for the part landward of the inflection point, which is the part used in this study, except for the equation for K. The longitudinal mixing coefficient (K) in Savenije [21] ( K s ) is given as:
K s = 0.16 × 10 6 h m ¯ 0.69 g 1.12 t 2.24 ( 2 a ( m ) ) 0.59 L W 1.1 W m 0.13
while, in Gisen et al. [23], it ( K g ) is given as:
K g = 151.35 × 10 6 W r 0.3 ( 2 a ( m ) ) 0.13 t 0.97 W m 0.30 C m 0.18 [ max ( u p e a k ( m o u t h , y ) ) ] 0.71 L W 0.11 h m ¯ 0.15 r s 0.84
in which a ( m ) is the tidal amplitude at the mouth, C m is the Chezy roughness, estimated in the assessment tool as 42 m 0.5 s 1 , which is adjustable in the tool code, max ( u p e a k ( m o u t h , y ) ) is the peak tidal flow velocity in the cross-section at the mouth, and r s is the storage ratio, calculated as the width of intertidal area at the mouth divided by the width of the subtidal area at the mouth. The minimum value of K is 0 and the maximum value is 1.

2.7. Tool Validation

For purposes of illustration and validation, the tool was applied to the Columbia River estuary (USA), which is an example of a relatively wide valley that is partially filled with bars [26], and the Western Scheldt (NL), which is an alluvial estuary with large deviations from a converging along-channel width profile and filled with tidal bars [19,26] (Figure 2c,d). For validation, channel depth at the landward boundary was calculated in the tool with measured river discharge. Predicted velocities are validated on the Columbia River only because the regression was based on data for the Western Scheldt. Bed elevation and salinity predictions are validated on both systems. The predicted flow velocities from the tool were compared with results from numerical modelling for the Columbia River estuary [7,28] and from a calibrated hydrodynamic model of Rijkswaterstaat for the Western Scheldt (SCALWEST, [32]). Results from numerical modelling were used because field measurements of flow velocities do not provide full coverage in time and space. Salinity data were obtained from Jay and Smith [51] for the Columbia River and from de Brauwere et al. [52] and Vroom et al. [53] for the Western Scheldt. Habitat preferences for potential habitat maps were obtained from the Marine Life Information Network (MarLIN), which is a database on the biology of species and the ecology of habitats.

2.8. Tool Application to End-Member Cases

We tested the tool’s performance against end-member systems. The dominant factors controlling the tool output are (i) the along-channel width profile, (ii) river discharge at the landward boundary and (iii) tidal range at the seaward boundary. Given these three factors, we designed the three most extreme end-member cases (Figure 2a,b,e), which were tested in addition to two estuarine cases with variable width profile (the Western Scheldt and Columbia River) (Figure 2c,d). The first end-member is an estuary with a perfectly converging shape from the mouth to the upstream river, i.e., negligible width variation. This case was based on the dimensions of the Tran De branch of the Mekong Delta (Vietnam). The second end-member is a river with a constant channel width, i.e., no tidal range at the seaward boundary and negligible width variation. This case was based on the dimensions of a distributary of the Saskatchewan River, near Cumberland Lake (Saskatchewan, Canada). The third end-member is a tidal basin, i.e., no river discharge and landward increasing width instead of decreasing. This case was based on dimensions of the Ameland inlet and represented by a narrow mouth, an increasing channel width in landward direction and an infinitely small river at the landward boundary. It should be noted that, in the last case, along-channel width is measured perpendicular to the tidal channels, which means that the transects for an idealised tidal basin are semicircles. The predictions for the Tran De branch were compared with measured data [54,55,56] and the predictions for the tidal basin were compared with numerical model results [31,57].

3. Results

3.1. Tool Output

The tool output comprises maps of predicted bed elevation with respect to the high water level (HWL), relative inundation duration, flow velocity and salinity (Figure 5 and Figure 6). In general, the bed deepens in seaward direction. At locations where the estuary is relatively wide, average bed elevations are higher and at locations where the estuary width is close to the ideal width, bed elevations are lower (Equations (4) and (5)). Zones where the excess width is large are the locations where stretches of intertidal area are wider and where inundation durations are less than 100%. Peak flow velocities increase in seaward direction with maxima around 1.8 ms 1 at the mouth. Peak flow velocities in the intertidal area are generally below 0.8 ms 1 . The salinity gradually decreases in landward directions and approaches river salinity at the landward boundary.
Before, we had to estimate a typical peak flow velocity because it was required for the estimation of the tidal excursion length (Equation (9)). At that point, a value of 1 ms 1 was assumed [18]. This assumption only partly influences the output flow velocity, which depends on bed elevation, which in turn is controlled by the estuary planform shape, the tidal range and the tidal excursion length. Nevertheless, output maps justify the assumption of 1 ms 1 for the purpose of calculating the tidal excursion length. The average peak tidal flow velocity is 1 ms 1 for the Western Scheldt and 1.2 ms 1 for the Columbia River when averaged over the surface area of a typical tidal excursion length (Figure 5 and Figure 6). Median values of peak tidal flow velocity per transect are also remarkably close to the typical value proposed by Savenije [18] along the entire estuary (Figure 7a,c).
Additional output (Appendix A Figure A3 and Figure A4) consists of predictions of bar width and braiding index per cross-section, maximum peak flow velocity and maximum average flow velocity per transect. Furthermore, four zones are calculated based on elevation with respect to the tidal range: subtidal deep (<−2a), subtidal low, intertidal low (more than 50% of time submerged) and intertidal high (less than 50% of time submerged). For each of these zones, the width of the zone, average depth and average and peak tidal flow velocities are calculated as output.

3.2. Tool Validation

3.2.1. Bed Elevations

To validate the predicted bed elevations, median and minimum values per transect were compared with measurements from bathymetry. Median and minimum bed elevations are generally within a factor 2 of the measured values from bathymetry (Figure 7a,b). Median bed elevations deviate the most in the landward part for the Western Scheldt (Figure 7a). In this section, a larger portion of the width is dredged, which causes the median depth to fall within the deeper dredged channel. The deviation between predicted and measured values lacks a trend with a position along the estuary (as indicated by the colours in Figure 7a,b). The largest outliers for the Columbia River are found at the mouth where the predicted median bed elevation is much lower than measured from bathymetry (Figure 7b). This is caused by the small embayments enclosed by groynes, located on the sides of the mouth area. In this zone, a larger proportion of high bed elevations causes a higher median bed elevation.

3.2.2. Flow Velocity

The range of predicted flow velocities is generally in the same range as the modelled values for both the peak and average flow velocities (Figure 7c,d). However, the along-channel variation in minimum and maximum flow velocities is larger for the models than for the tool output (Figure 8). We explain this by a simplification of some morphological elements in the tool. For example, the presence and location of scours and sills, which cause variation in the along-channel maximum depth, is not predicted by the tool. However, they do occur in natural systems and therefore also cause variations in along-channel flow velocity pattern. In the tool, the profile of maximum depth is linear along-channel, which results in lower variations in the maximum flow velocity. Nevertheless, the velocity predictions do capture the along-channel variability caused by the large-scale morphological variation in bed elevation, as illustrated with the median values per transect (Figure 8).
The tool does not reproduce the highest observed flow velocities in the Columbia River, which are well above 2 ms 1 in the first 15 km from the mouth. To obtain velocities above 2 ms 1 in the tool, bed elevation of at least 120 m below high water level are required (Equation (17)), which generally do not occur in natural systems. This means that the tool underpredicts flow velocities in cases where the cross-sectional area is relatively small compared to the local tidal prism, which might also suggest disequilibrium. Moreover, the regression between flow velocity and depth, simplifies the increasing scatter for larger depths and flow velocities, which is observed for the Western Scheldt (Figure 4).
Predicted velocity against modelled velocity shows a trend for the Columbia River with modelled velocities below 1 ms 1 being overpredicted by the tool and velocities above 1 ms 1 being underpredicted (Figure 7d). The degree of misprediction does not correlate with position along the estuary. The only explanation for the deviation is a different depth dependency of peak flow velocity caused by a different ratio between tidal prism and cross-sectional area.
Therefore, an alternative approach for calculating flow velocities was based on the local tidal prism and the cross-sectional channel geometry, i.e., continuity. For each cross-section, a typical velocity is calculated at the average depth of the cross-section (i.e., P ( x ) / ( T h ( x ) ¯ ) ). Nevertheless, the quality of these predictions appeared to be less (Figure 9) for two reasons: (1) at the upstream boundary, velocities are mainly sensitive to river discharge. Here, local tidal prism approaches river discharge because it is obtained from integrating over the distance of a tidal excursion length upstream (Equation (9)) plus river discharge. In the case of the Columbia River, this leads to overprediction of velocities because river discharge is large, while, for the Western Scheldt, it would lead to underpredictions. (2) Flow velocities at the intertidal area and deepest parts still have to be extrapolated from the average prediction creating uncertainty in the most important parts for respectively ecology and shipping.
Following this alternative approach, we anticipated that there would be a relation between the width variation and the misprediction based on depth only (Equation (17)). It was expected that a local increase in width would lead to flow expansion and therefore typically lower flow velocities than for a cross-section with similar tidal prism but lower width. However, the lack of trend between velocity deviation and width measures (Appendix A Figure A2) indicates that velocity predictions will not improve from including width in the regression.
The alternative approach shows that a prediction of flow velocities depending on tidal prism and depth does not reproduce modelled flow velocities better (Figure 9a,b,e,f). A predictive relation between flow velocity and depth is thus of value. This urges the need to compile a dataset with flow velocities and water depths of multiple systems. From that data, a relation that is better applicable to more systems could be derived, but this data is not yet available. Nevertheless, the present relation predicted flow velocity within a factor 2 of the modelled values even in the most extreme cases, which means it is of value for an estimate when limited data is available.
Predicted flow velocity and local depth result in a local tidal prism that is in the same range as the modelled values (Figure 10). Local tidal prism obtained from integrating a tidal range over the full length of the estuary is typically 1.5–2 times larger at the mouth than the tidal prism obtained from integrating over a typical tidal excursion length. This explains why using the latter is required to obtain accurate predictions of the depth at the estuary mouth and additionally is less sensitive to the length of the estuary.

3.2.3. Salinity

Predicted salinities are within the range of measured values [51] and values obtained with a numerical model with shallow-water and tracer-transport equations [52]. For the Western Scheldt, local salinity can vary between 4–8 ppt over a tidal cycle, with averages of 28 ppt at 10 km, 22 ppt at 30 km and 16 ppt at 45 km from the mouth [52]. In the Columbia River estuary, temporal salinity variation is even larger [51], with local variations of 20–25 ppt in the most seaward 15 km. Predicted salinities fall within the measured range, but salinities are overpredicted in the most landward reach, where measured values approach 0 ppt from approximately 45 km from the mouth [51].
The salinity predictor of Gisen and Savenije [50] is most accurate for the salinities observed in the Western Scheldt (Figure 7e), which was a system used for validation in Gisen and Savenije [50] and not for calibration. This is the only salinity predictor for which predicted salinity is within a factor 1.5 of the measured value for the entire along-channel section. Predictors of Savenije [21] and Brockway et al. [22] show large deviations at the landward side, where the predicted salinity gradient is respectively too weak and too strong compared to measured values (Figure 7e).

3.3. Applicability to End-Member Systems: A River, an Ideal Estuary, and a Tidal Basin

In the case of a river and ideal estuary, the local channel width will always be equal to the ideal channel width as obtained from a maximum fitting converging shape. This means that the shape of the cross-sectional hypsometry is concave and constant along-channel (Equations (4) and (5), Figure 11). The concave hypsometry results in a relatively constant along-channel bed elevation profile (Figure 11). Only the along-channel variation in maximum depth and increase in channel width in seaward direction cause a slight along-channel variation. The implication of an along-channel constant concave hypsometry is that intertidal area is very small and inundation duration is 100% for the largest part of the estuary. Typical flow velocities are also along-channel constant at approximately 1–1.2 ms 1 . A strong salinity gradient is predicted at the mouth of the estuary, with salinities dropping to ≈5 ppt at 30 km from the mouth.
Ideal estuaries are said to be characterised by an along-channel constant velocity amplitude and depth [18,54] and these features are reproduced by the tool predictions. Moreover, the predictions correspond to observations in the Tran De branch of the Mekong: observed peak flow velocity is approximately 1 ms 1 along the entire branch [55,56] and salinity intrusion reaches up to 45–55 km from the mouth [54]. In addition, the along-channel predicted salinities fall within the range of observed salinities [54]. The only feature not completely reproduced is the along-channel constant depth. When the maximum depth at the mouth is calculated from tidal prism, predicted depth at the mouth is a factor 1.5 larger, 5 m deeper, than observed.
The prediction for depth at the mouth leads to a misprediction at the most downstream boundary for the river (Appendix A Figure A6). In that case, maximum depth at the landward boundary is controlled by river discharge (to the power of ≈0.35), the width of the river without tidal influence (to the power of ≈0.77) or if both of these are unspecified, the channel width at the landward boundary. Depth at the mouth is calculated with the tidal prism, which in the case of a river is equal to the river discharge times the duration of half a tidal cycle. Therefore, depth at the mouth is also a function of river discharge or landward river width, but to the power of ≈1. This means that the predicted channel depth at the mouth will always be deeper than at the upstream boundary, especially in the end-member case of a straight river.
An ideal tidal basin is characterised by a strong landward increase in width, while the ideal width approaches 0 at the landward end. This means that the ratio of W i d e a l ( x ) to W ( x ) strongly decreases in landward direction, resulting in convex hypsometry (Equation (4)) and therefore large amounts of intertidal area within the basin (Figure A7). In this end-member case, the depth at the mouth is set by the surface area of the tidal basin multiplied with the tidal range. In the case of a tidal basin based on the Ameland inlet, the tool predicts a maximum depth at the mouth of 20 m below the high water level. This seems appropriate when compared with modelled [31,57] and measured values [29] (Figure A8). Flow velocities at the mouth are in the order of 1 ms 1 , which is comparable with measured and modelled values for tidal basins in equilibrium [57,58].

4. Discussion

4.1. Sensitivity of Tool Output to Tidal Range, River Discharge, and along-Channel Width Profile

The sensitivity of the input tidal range, river discharge, and along-channel width profile was tested on the following parameters: (1) the ratio between intertidal and subtidal area; (2) the ratio between low and high dynamic environment; and (3) peak flow velocity at the estuary mouth (Figure 12). Low dynamic environment, which forms potential habitat, is defined as intertidal area where flow velocities are below 0.8 ms 1 . All other areas are classified as a high dynamic area.
While most sensitive dependent variables are increasing or decreasing monotonously with the independent variable, the ratio of low and high dynamic intertidal area as a function of tidal amplitude shows an optimum. An increase in tidal amplitude increases the ratio between intertidal and subtidal area (Figure 12a) because the increased tidal range increases intertidal area. However, increases in tidal amplitude also generate a proportional increase in tidal prism. A larger tidal prism at the mouth will increase the predicted cross-sectional area and because the width is fixed; this will increase predicted depth at the mouth. Because maximum depth along the estuary is dependent on the maximum depth at the mouth and the maximum depth at the landward river, this increases the maximum depth along the entire estuary. The increased depths subsequently increase the predicted flow velocities (Figure 12g). This means that two opposing mechanisms affect the ratio between low dynamic and high dynamic environment: increased tidal range increases intertidal area which forms a more low-dynamic environment. In contrast, the deepening by increase in tidal prism increases flow velocities and therefore reduces the low dynamic area. Therefore, the combined effect results in an optimum curve for the ratio between low dynamic (velocity < 0.8 ms 1 ) intertidal area and high dynamic area as a function of tidal range (Figure 12d).
An increase in the river discharge results in less intertidal area, less low dynamic area and larger peak flow velocities (Figure 12b,e,h). After all, increased river discharge increases the tidal prism, which increases channel depths.
The ratio between ideal width and excess width affects the distribution of bed elevations per cross-section and the tidal prism. In estuaries that adapt to the ideal shape, the excess width becomes lower. This results in more concave cross-sectional hypsometry and therefore narrower stretches of intertidal flats (Figure 12c). The tidal prism decreases as a result of reduced estuary area, resulting in lower maximum depths and therefore also lower flow peak velocities (Figure 12i).
A comparison between predicted and measured or modelled values showed that bed elevations, flow velocity and salinity can be predicted within a factor two of the measured or modelled values (Figure 7). The largest deviations between predicted and modelled flow velocities occur when they peak above 2 ms 1 in natural systems. Largest mispredictions for salinity occur at the transition from brackish to fresh water, where salinity is also the most dynamic over the seasons. Bed elevation predictions are most uncertain when the tidal range is low for the depth at the mouth, which is also sensitive to offshore wave climate that we have not considered here. This implies that the tool is capable to provide a first-order estimation for estuaries with limited data, but that it cannot replace measurements or numerical modelling.

4.2. Applications

4.2.1. Illustration of Application for Ecological Assessment

The major abiotic factors that determine habitat area of living organisms in estuaries are depth, flow velocity, salinity and inundation time [2,3,4,59]. Habitat suitability models (e.g., HABITAT, [24]) can translate these abiotic factors into a characteristic habitat suitability index, which for example has been applied to study the effect of climate change and dam construction on potential habitat in rivers [25]. Multiple case studies of individual systems or individual species are also available for estuaries [60,61,62,63,64], but application to more systems worldwide and application to scenarios of future climate change and human engineering are hampered by a lack of data.
The tool presented in this study can provide typical abiotic factors that are input for habitat suitability models within minutes. Therefore, the output can lead to a first-order estimation of the potential habitat (Figure 13) based on species specific preferences (Table 1). The potential habitat area based on tool output generally overpredicts the measured habitat. First, the tool always predicts a zone with intertidal area because the predicted hypsometry is stretched between the high water line and maximum channel depth. Second, other constraining biotic and abiotic factors were excluded, such as the competition between species and the water quality. Lastly, the settling and growth conditions may deviate from the conditions under which species occur in a later life stage. For example, seedlings of salt marsh species generally only settle when peak flow is below 0.25 ms 1 and, in their initial life stage, they may only sustain flows up to 0.5 ms 1 [65]. Nevertheless, the application of the tool for ecological assessment is of value because it is capable of providing abiotic factors of individual systems and scenarios within minutes, rather than days to months for numerical modelling and in situ measurements.

4.2.2. Application for Management of Estuaries

The main channels in estuaries are often used as shipping fairways, while the intertidal area forms important ecological habitat [2]. Given a certain channel width and depth required for shipping, the empirical predictions of bed elevations that result from this study can estimate volumes that need to be dredged and typical locations where this is necessary. While, for many of the present-day shipping fairways, bathymetry is available, this application may be relevant for quick assessments by consultancy in developing countries with growing economies for the construction of fairways, harbours or studies on land loss in these areas [66,67].
Another implication from this study is that the wider an estuary is relative to the fitted converging shape, the wider the zone of intertidal area will be. Bed elevations and basin hypsometry form an indication of possible habitat composition [35] and low-dynamical intertidal areas are generally labelled as valuable areas [68] because they are suitable for settling and feeding. The output can indicate which areas satisfy these conditions (e.g., suitable flow conditions) and which areas can be transformed into habitat by disposal of dredged sediment. For example, in the Western Scheldt, dredged sediment has been disposed on the edges of tidal flats to maintain intertidal area [68], which aids in achieving European Union Natura 2000 regulations [69]. This approach requires bathymetry, which may not be available at locations where new fairways or harbours are planned. The predictive tool can provide a first-order estimation in these cases.

4.2.3. Application in Palaeogeographic Reconstruction

Our empirical method is capable of making quantitative estimates of the bar dimensions, the number of bars and channels (Braiding Index), channel dimensions, flow conditions and salinities. The tool is applicable for palaeogeographical reconstructions [70,71] where the shape of tidal channel belts or estuaries has been identified, but wherein bar configuration is typically not recognisable and typical flow conditions are unknown. If the along-channel width profile is known and the tidal range can be estimated, the tool presented in this study can be applied. This application may be particularly relevant for high-stand estuaries that are filled with sediment. Besides the forward approach of palaeogeographical reconstructions based on estuary outline, the tool and equations also allow a reverse approach. A single measurement or preferably a few measurements of bar or channel dimensions [12,72,73,74,75,76] would be enough to reconstruct typical channel width and local tidal prism and therefore also typical flow conditions.
Additionally, the combination with ecological habitat analyses may improve the understanding of the evolution of estuaries. It is known that vegetation succession stabilises tidal sand bars [65,77,78,79], but it remains an open question to what extent this succession confines the width of estuaries on the longer term. Additionally, it is unknown what the equilibrium shape of an estuary will be under a given set of stabilising and destabilising eco-engineering species.

5. Conclusions

A comprehensive set of empirical relations allows estimating bathymetry, flow velocity, salinity and inundation durations for estuaries for which limited data are available. Morphological characteristics are based on the concept that the degree to which an estuary deviates from a maximum fitting converging shape is reflected in the locations where intertidal bars are found, which also determines the bed elevation distributions. This study provides a new correlation between depth below high water level and peak tidal flow velocity, based on numerical model outcomes, which complements the bed elevation predictions and allows for salinity predictions. The resulting predictions are within a factor two of the measured or modelled values, but generally better. Largest deviations between predicted and modelled or measured occur for flow velocities when they peak above 2 ms 1 , for salinity at the transition from brackish to fresh water and for bed elevations when the tidal range is low for the depth at the mouth. The end-member cases of river-dominated delta branches and tide-dominated delta branches, i.e., rivers and ideal estuaries, are most sensitive to depth at both boundaries. In these cases, when possible, it should be measured in a single cross-section and specified as input to prevent misprediction.
We conclude that the tool is suitable for alluvial estuaries filled with tidal bars, for relatively wide valleys that partially filled with bars and for tidal basins. The results are usable to quantify valuable ecological habitat area, to estimate dredging volumes in management of estuaries or to make palaeographic reconstructions when limited data are available.

Supplementary Materials

The following is available at https://www.mdpi.com/2072-4292/10/12/1915/s1, A zip file with tool code and instructions to use the tool. The ZIP contains: (1) ’Instructions.pdf’ with instructions on how to use the tool; (2) ’Input_variables.xls’, which is file where input parameters need to be specified; (3) ’.csv’ files with examples of along-channel width profiles; (4) ’Model_v1_0.py’, which is the tool code that needs to be run in Python; (5) ’Examples’ folder with examples of the output for the Western Scheldt.

Author Contributions

Contributions were in following proportions to conception and design, tool programming, data collection and analysis, conclusions and manuscript preparation: J.R.F.W.L. (65%, 50%, 60%, 85%), S.L.V. (5%, 50%, 20%, 5%), W.M.v.D. (5%, 0%, 10%, 5%), S.S. (5%, 0%, 5%, 0%), and M.G.K. (20%, 0%, 5%, 5%).

Funding

J.R.F.W.L., W.M.v.D. and M.G.K. were supported by the Dutch Technology Foundation TTW (grant Vici 016.140.316/13710 to M.G.K., which is part of the Netherlands Organisation for Scientific Research (NWO), and is partly funded by the Ministry of Economic Affairs). S.L.V. was supported by Future Deltas, Utrecht University (seed grant to J.R.F.W.L.). S.S. was supported by an ERC Consolidator grant (647570) to M.G.K.

Acknowledgments

This work is part of the PhD research of J.R.F.W.L. and MSc research of S.L.V. We thank three anonymous reviewers for their comments on the manuscript. We thank Rijkswaterstaat for providing bathymetry of the Western Scheldt and Edwin Elias for providing the model of the Columbia River. Discussion with Sepehr Eslami Arab, Barend van Maanen and Tjeerd Bouma helped to improve the manuscript. The code of the tool is made available open access online with the manuscript and through GitHub (https://github.com/JasperLeuven/EstuarineMorphologyEstimator). Other data sources have been referenced in the text.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Supplementary Figures and Table

Figure A1. Cross-sectional profiles at the mouth (left) and landward river boundary (right) for five estuaries. The maximum depth (solid red line) and an approximation of the average water depth (dashed line) are indicated. The ratio between maximum depth and average depth, which is used in Equations (12) and (8), is given in the figures. Average ratio for the geometry of the estuary mouth is 1.85 and for the landward river boundary 1.65.
Figure A1. Cross-sectional profiles at the mouth (left) and landward river boundary (right) for five estuaries. The maximum depth (solid red line) and an approximation of the average water depth (dashed line) are indicated. The ratio between maximum depth and average depth, which is used in Equations (12) and (8), is given in the figures. Average ratio for the geometry of the estuary mouth is 1.85 and for the landward river boundary 1.65.
Remotesensing 10 01915 g0a1
Figure A2. Measured velocity from the model divided by predicted velocity from the tool against three expressions of excess width. Lack of trend between velocity deviation and width measures indicates that velocity predictions will not improve from including width in the regression. (a) ideal width divided by local width, (b) excess width divided by ideal width and (c) along-channel derivative of excess width divided by ideal width. Colours in the plot become more yellow with increasing point density.
Figure A2. Measured velocity from the model divided by predicted velocity from the tool against three expressions of excess width. Lack of trend between velocity deviation and width measures indicates that velocity predictions will not improve from including width in the regression. (a) ideal width divided by local width, (b) excess width divided by ideal width and (c) along-channel derivative of excess width divided by ideal width. Colours in the plot become more yellow with increasing point density.
Remotesensing 10 01915 g0a2
Figure A3. Tool output for the Western Scheldt. (top) predicted along-channel bar width, (middle) predicted along-channel braiding index and (bottom) the cross-sectional maximum and average of the predicted along-channel peak tidal flow velocity and the cross-sectional average of the mean flow velocity over half a tidal cycle.
Figure A3. Tool output for the Western Scheldt. (top) predicted along-channel bar width, (middle) predicted along-channel braiding index and (bottom) the cross-sectional maximum and average of the predicted along-channel peak tidal flow velocity and the cross-sectional average of the mean flow velocity over half a tidal cycle.
Remotesensing 10 01915 g0a3
Figure A4. Tool output for the Western Scheldt. (Top) along-channel width of subtidal, intertidal low and intertial high zones; (center) average depth of subtidal, intertidal low and intertial high zones; (bottom) average and maximum flow velocity of subtidal, intertidal low and intertial high zones.
Figure A4. Tool output for the Western Scheldt. (Top) along-channel width of subtidal, intertidal low and intertial high zones; (center) average depth of subtidal, intertidal low and intertial high zones; (bottom) average and maximum flow velocity of subtidal, intertidal low and intertial high zones.
Remotesensing 10 01915 g0a4aRemotesensing 10 01915 g0a4b
Figure A5. (Left) maximum ebb velocity as a function of mean ebb velocity; (Right) maximum ebb velocity as a function of mean ebb velocity. Histograms indicate occurrence within bins and colours in the plot become darker with increasing density. Drawn lines indicate the relation x = 2 / π y , which is the expected average velocity when calculated as the average of half a sine function with an amplitude equal to the peak velocity.
Figure A5. (Left) maximum ebb velocity as a function of mean ebb velocity; (Right) maximum ebb velocity as a function of mean ebb velocity. Histograms indicate occurrence within bins and colours in the plot become darker with increasing density. Drawn lines indicate the relation x = 2 / π y , which is the expected average velocity when calculated as the average of half a sine function with an amplitude equal to the peak velocity.
Remotesensing 10 01915 g0a5
Figure A6. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for a river based on dimensions of a distributary of the Saskatchewan River, Cumberland Lake (Saskatchewan (Canada), coordinates: 54 04 N 102 22 W) (Figure 2a). In the left panels (ad) the boundary condition for maximum depth at the mouth and the landward end are both based on hydraulic geometry for rivers (Equations (6)–(8)). In the right panels (eh), the boundary condition for maximum depth at the mouth is based on the tidal prism relation (Equations (10) and (12)) and boundary condition at the landward end is based on hydraulic geometry for rivers (Equations (8)–(6)). See Bolla Pittaluga et al. [80] for measurements in this river.
Figure A6. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for a river based on dimensions of a distributary of the Saskatchewan River, Cumberland Lake (Saskatchewan (Canada), coordinates: 54 04 N 102 22 W) (Figure 2a). In the left panels (ad) the boundary condition for maximum depth at the mouth and the landward end are both based on hydraulic geometry for rivers (Equations (6)–(8)). In the right panels (eh), the boundary condition for maximum depth at the mouth is based on the tidal prism relation (Equations (10) and (12)) and boundary condition at the landward end is based on hydraulic geometry for rivers (Equations (8)–(6)). See Bolla Pittaluga et al. [80] for measurements in this river.
Remotesensing 10 01915 g0a6
Figure A7. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for a tidal basin with dimensions based on the Ameland inlet.
Figure A7. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for a tidal basin with dimensions based on the Ameland inlet.
Remotesensing 10 01915 g0a7aRemotesensing 10 01915 g0a7b
Figure A8. The blue solid line indicates the resulting basin hypsometry from the prediction for a tidal basin with dimensions based on the Ameland inlet (Figure A7) plotted in Figure 7 of van Maanen et al. [31].
Figure A8. The blue solid line indicates the resulting basin hypsometry from the prediction for a tidal basin with dimensions based on the Ameland inlet (Figure A7) plotted in Figure 7 of van Maanen et al. [31].
Remotesensing 10 01915 g0a8
Table A1. List of symbols used.
Table A1. List of symbols used.
SymbolUnitsVariable
A c s a , m [m 2 ]Cross-sectional area at the mouth
A c s a ( x ) [m 2 ]Cross-sectional area at coordinate x
a ( m ) [m]Tidal range at the mouth
a ( x ) [m]Tidal range at the mouth
B I ( x ) [-]Braiding Index at coordinate x
C m [m 0.5 s 1 ]Chezy roughness at the estuary mouth
D m [m 2 s 1 ]Dispersion coefficient at the mouth of the estuary
D ( x ) [m 2 s 1 ]Dispersion coefficient at coordinate x
Δ ρ [-]Density difference between salt and fresh water
E[m]Tidal excursion length
g[ms 2 ]Gravitational acceleration
h m a x , r [m]Maximum channel depth at the landward boundary
h r ¯ [m]Average channel depth at the landward boundary
h m a x , m [m]Maximum channel depth at the estuary mouth
h m ¯ [m]Average channel depth at the estuary mouth
h z ( x , y ) [-]Dimensionless bed elevation at coordinate (x,y)
h ( x , y ) [m]Dimensional bed elevation at coordinate (x,y)
I ( x , y ) [× t]Inundation duration
K[m 2 s 1 ]Empirical equation for the longitudinal mixing coefficient
K s [-]Longitudinal Van der Burgh mixing coefficient in Savenije [21]
K g [-]Longitudinal Van der Burgh mixing coefficient in Gisen et al. [23]
L W [m]Width convergence length
L A c s a [m]Cross-sectional area convergence length
N r [-]Estuarine Richardson number
P ( x ) [m 3 ]Local tidal prism at coordinate x
P ( m ) [m 3 ]Local tidal prism at the estuary mouth
Q b [m 3 s 1 ]Bankfull river discharge
Q r [m 3 s 1 ]Average river discharge
r[-]Coefficient in the Strahler [33] equation
r s [-]Storage ratio (intertidal area/subtidal area)
ρ w [kgm 3 ]Water density
S m [ppt]Salinity at the estuary mouth
S r [ppt]Salinity at the landward boundary
S b ( x ) [ppt]Salinity at coordinate x [22]
S g ( x ) [ppt]Salinity at coordinate x [23]
S s ( x ) [ppt]Salinity at coordinate x [21]
s[s]Distance between the mouth and landward boundary
s m [-]Shape factor of the channel at the estuary mouth
s r [-]Shape factor of the channel at the landward boundary
t[s]Duration of half a tidal cycle
u a v g ( x , y ) [ms 1 ]Tidal average flow velocity at coordinate (x,y)
u p e a k ( x , y ) [ms 1 ]Peak tidal flow velocity at coordinate (x,y)
u p e a k ( m o u t h , y ) [ms 1 ]Peak tidal flow velocity at the estuary mouth
W e x c e s s ( x ) [m]Excess width at coordinate x
W i d e a l ( x ) [m]Ideal width at coordinate x
W m [m]Width at the estuary mouth
W r [m]Width at the landward boundary
W ( x ) [m]Local width at coordinate x
w b a r ( x ) [m]Predicted bar width at coordinate x
w b a r s s u m m e d ( x ) [m]Summed width of bars at coordinate x
x[m]Streamwise coordinate measured from the mouth along the centreline
y[m]Coordinate perpendicular to the centreline
z[-]Coefficient in the Strahler [33] equation
z ( x ) [-]Value of z-coefficient at coordinate x

References

  1. De Vriend, H.J.; Wang, Z.B.; Ysebaert, T.; Herman, P.M.; Ding, P. Eco-morphological problems in the Yangtze Estuary and the Western Scheldt. Wetlands 2011, 31, 1033–1042. [Google Scholar] [CrossRef]
  2. Bouma, H.; de Jong, D.; Twisk, F.; Wolfstein, K. A Dutch Ecotope System for Coastal Waters (ZES. 1). To Map the Potential Occurence of Ecological Communities in Dutch Coastal And Transitional Waters; Technical Report, Rijkswaterstaat, Report RIKZ/2005.024; Rijksinstituut voor Kust en Zee/RIKZ: Middelburg, The Netherlands, 2005. [Google Scholar]
  3. De Jong, D. Ecotopes in the Dutch Marine Tidal Waters: A Proposal for a Classification of Ecotopes and a Method to Map Them; Technical Report, Rijkswaterstaat, Report RIKZ/99.017; Rijksinstituut voor Kust en Zee/RIKZ: Middelburg, The Netherlands, 1999. [Google Scholar]
  4. Gurnell, A.M.; Bertoldi, W.; Corenblit, D. Changing river channels: The roles of hydrological processes, plants and pioneer fluvial landforms in humid temperate, mixed load, gravel bed rivers. Earth-Sci. Rev. 2012, 111, 129–141. [Google Scholar] [CrossRef]
  5. Cozzoli, F.; Smolders, S.; Eelkema, M.; Ysebaert, T.; Escaravage, V.; Temmerman, S.; Meire, P.; Herman, P.M.; Bouma, T.J. A modeling approach to assess coastal management effects on benthic habitat quality: A case study on coastal defense and navigability. Estuar. Coast. Shelf Sci. 2017, 184, 67–82. [Google Scholar] [CrossRef] [Green Version]
  6. Van der Wegen, M.; Roelvink, J. Reproduction of estuarine bathymetry by means of a process-based model: Western Scheldt case study, the Netherlands. Geomorphology 2012, 179, 152–167. [Google Scholar] [CrossRef]
  7. Elias, E.P.; Gelfenbaum, G.; Van der Westhuysen, A.J. Validation of a coupled wave-flow model in a high-energy setting: The mouth of the Columbia River. J. Geophys. Res. Oceans 2012, 117, C09011. [Google Scholar] [CrossRef]
  8. Braat, L.; van Kessel, T.; Leuven, J.R.F.W.; Kleinhans, M.G. Effects of mud supply on large-scale estuary morphology and development over centuries to millennia. Earth Surf. Dyn. 2017, 5, 617–652. [Google Scholar] [CrossRef] [Green Version]
  9. Townend, I. The estimation of estuary dimensions using a simplified form model and the exogenous controls. Earth Surf. Process. Landf. 2012, 37, 1573–1583. [Google Scholar] [CrossRef]
  10. Seminara, G.; Tubino, M. Sand bars in tidal channels. Part 1. Free bars. J. Fluid Mech. 2001, 440, 49–74. [Google Scholar] [CrossRef]
  11. Schramkowski, G.; Schuttelaars, H.; De Swart, H. The effect of geometry and bottom friction on local bed forms in a tidal embayment. Cont. Shelf Res. 2002, 22, 1821–1833. [Google Scholar] [CrossRef] [Green Version]
  12. Leuven, J.R.F.W.; Kleinhans, M.G.; Weisscher, S.A.H.; van der Vegt, M. Tidal sand bar dimensions and shapes in estuaries. Earth-Sci. Rev. 2016, 161, 204–233. [Google Scholar] [CrossRef]
  13. O’Brien, M.P. Equilibrium flow areas of inlets on sandy coasts. In Journal of the Waterways and Harbors Division: Proceedings of the American Society of Civil Engineers (ASCE); American Society of Civil Engineers: New York, NY, USA, 1969; pp. 43–52. [Google Scholar] [CrossRef]
  14. Wang, Z.; Jeuken, C.; De Vriend, H. Tidal Asymmetry and Residual Sediment Transport in Estuaries; Delft Hydraulics Report Z2749; Deltares (WL): Delft, The Netherlands, 1999. [Google Scholar]
  15. Lanzoni, S.; D’Alpaos, A. On funneling of tidal channels. J. Geophys. Res. Earth Surf. 2015, 120, 433–452. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, Z.; Coco, G.; Townend, I.; Gong, Z.; Wang, Z.; Zhang, C. On the stability relationships between tidal asymmetry and morphologies of tidal basins and estuaries. Earth Surf. Process. Landf. 2018, 43, 1943–1959. [Google Scholar] [CrossRef]
  17. Langbein, W. The hydraulic geometry of a shallow estuary. Hydrol. Sci. J. 1963, 8, 84–94. [Google Scholar] [CrossRef]
  18. Savenije, H.H. Salinity and Tides in Alluvial Estuaries; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  19. Leuven, J.R.F.W.; Haas, T.; Braat, L.; Kleinhans, M.G. Topographic forcing of tidal sand bar patterns for irregular estuary planforms. Earth Surf. Process. Landf. 2018, 43, 172–186. [Google Scholar] [CrossRef]
  20. Leuven, J.R.F.W.; Selaković, S.; Kleinhans, M.G. Morphology of bar-built estuaries: empirical relation between planform shape and depth distribution. Earth Surf. Dyn. 2018, 6, 763–778. [Google Scholar] [CrossRef]
  21. Savenije, H.H. Predictive model for salt intrusion in estuaries. J. Hydrol. 1993, 148, 203–218. [Google Scholar] [CrossRef]
  22. Brockway, R.; Bowers, D.; Hoguane, A.; Dove, V.; Vassele, V. A note on salt intrusion in funnel-shaped estuaries: Application to the Incomati estuary, Mozambique. Estuar. Coast. Shelf Sci. 2006, 66, 1–5. [Google Scholar] [CrossRef]
  23. Gisen, J.; Savenije, H.; Nijzink, R. Revised predictive equations for salt intrusion modelling in estuaries. Hydrol. Earth Syst. Sci. 2015, 19, 2791–2803. [Google Scholar] [CrossRef]
  24. Haasnoot, M.; van de Wolfshaar, K. Combining a conceptual framework and a spatial analysis tool, HABITAT, to support the implementation of river basin management plans. Int. J. River Basin Manag. 2009, 7, 295–311. [Google Scholar] [CrossRef]
  25. Van Oorschot, M.; Kleinhans, M.; Buijse, T.; Geerling, G.; Middelkoop, H. Combined effects of climate change and dam construction on riverine ecosystems. Ecol. Eng. 2018, 120, 329–344. [Google Scholar] [CrossRef]
  26. Dalrymple, R.W.; Zaitlin, B.A.; Boyd, R. Estuarine facies models: Conceptual basis and stratigraphic implications: Perspective. J. Sediment. Res. 1992, 62. [Google Scholar] [CrossRef]
  27. De Haas, T.; Pierik, H.; van der Spek, A.; Cohen, K.; van Maanen, B.; Kleinhans, M. Holocene evolution of tidal systems in The Netherlands: Effects of rivers, coastal boundary conditions, eco-engineering species, inherited relief and human interference. Earth-Sci. Rev. 2017, 177, 139–163. [Google Scholar] [CrossRef]
  28. Stevens, A.; Gelfenbaum, G.; MacMahan, J.; Reniers, A.; Elias, E.; Sherwood, C.; Carlson, E. Oceanographic Measurements and Hydrodynamic Modeling of the Mouth of the Columbia River, Oregon and Washington, 2013; Technical Report; U.S. Geological Survey: Reston, VA, USA, 2017. [CrossRef]
  29. Marciano, R.; Wang, Z.B.; Hibma, A.; de Vriend, H.J.; Defina, A. Modeling of channel patterns in short tidal basins. J. Geophys. Res. Earth Surf. 2005, 110, F01001. [Google Scholar] [CrossRef]
  30. Wang, Z.B.; Hoekstra, P.; Burchard, H.; Ridderinkhof, H.; De Swart, H.E.; Stive, M.J.F. Morphodynamics of the Wadden Sea and its barrier island system. Ocean Coast. Manag. 2012, 68, 39–57. [Google Scholar] [CrossRef]
  31. van Maanen, B.; Coco, G.; Bryan, K.R. Modelling the effects of tidal range and initial bathymetry on the morphological evolution of tidal embayments. Geomorphology 2013, 191, 23–34. [Google Scholar] [CrossRef]
  32. Verbeek, H.; Van der Male, K.; Jansen, M. Het SCALWEST-model (in Dutch); Technical Report, Technical Report RIKZ/OS/2000.814; Rijksinstituut voor Kust en Zee/RIKZ: Middelburg, The Netherlands, 2000. [Google Scholar]
  33. Strahler, A.N. Hypsometric (area-altitude) analysis of erosional topography. Geol. Soc. Am. Bull. 1952, 63, 1117–1142. [Google Scholar] [CrossRef]
  34. Boon, J.D.; Byrne, R.J. On basin hyposmetry and the morphodynamic response of coastal inlet systems. Mar. Geol. 1981, 40, 27–48. [Google Scholar] [CrossRef]
  35. Townend, I.H. Hypsometry of estuaries, creeks and breached sea wall sites. Proc. Inst. Civ. Eng.-Marit. Eng. 2008, 161, 23–32. [Google Scholar] [CrossRef]
  36. Leopold, L.B.; Maddock, T., Jr. The Hydraulic Geometry of Stream Channels and Some Physiographic Implications; Technical Report, U.S. Geological Survey, Professional Paper 252; U.S. Geological Survey: Reston, VA, USA, 1953.
  37. Hey, R.D.; Thorne, C.R. Stable channels with mobile gravel beds. J. Hydraul. Eng. 1986, 112, 671–689. [Google Scholar] [CrossRef]
  38. Eysink, W. Morphologic response of tidal basins to changes. Coast. Eng. Proc. 1990, 1, 1948–1961. [Google Scholar] [CrossRef]
  39. Friedrichs, C.T. Stability shear stress and equilibrium cross-sectional geometry of sheltered tidal channels. J. Coast. Res. 1995, 11, 1062–1074. [Google Scholar]
  40. Pavelsky, T.M.; Smith, L.C. RivWidth: A software tool for the calculation of river widths from remotely sensed imagery. IEEE Geosci. Remote Sens. Lett. 2008, 5, 70–73. [Google Scholar] [CrossRef]
  41. Donchyts, G.; Schellekens, J.; Winsemius, H.; Eisemann, E.; van de Giesen, N. A 30 m resolution surface water mask including estimation of positional and thematic differences using landsat 8, srtm and openstreetmap: A case study in the Murray-Darling Basin, Australia. Remote Sens. 2016, 8, 386. [Google Scholar] [CrossRef]
  42. Savenije, H.H. Prediction in ungauged estuaries: An integrated theory. Water Resour. Res. 2015, 51, 2464–2476. [Google Scholar] [CrossRef] [Green Version]
  43. Davies, G.; Woodroffe, C.D. Tidal estuary width convergence: Theory and form in North Australian estuaries. Earth Surf. Process. Landf. 2010, 35, 737–749. [Google Scholar] [CrossRef]
  44. Pillsbury, G. Tidal Hydraulics, rev. ed.; US Army Corps of Engineers: Washington, DC, USA, 1956.
  45. Dronkers, J. Convergence of estuarine channels. Cont. Shelf Res. 2017, 144, 120–133. [Google Scholar] [CrossRef]
  46. Wang, Z.; Jeuken, M.; Gerritsen, H.; De Vriend, H.; Kornman, B. Morphology and asymmetry of the vertical tide in the Westerschelde estuary. Cont. Shelf Res. 2002, 22, 2599–2609. [Google Scholar] [CrossRef]
  47. Toffolon, M.; Crosato, A. Developing macroscale indicators for estuarine morphology: The case of the Scheldt Estuary. J. Coast. Res. 2007, 23, 195–212. [Google Scholar] [CrossRef]
  48. Jarrett, J.T. Tidal Prism-Inlet Area Relationships; Technical Report, DTIC Document WES-GITI-3; DTIC: Fort Belvoir, VA, USA, 1976.
  49. Shigemura, T. Tidal prism–throat area relationships of the bays of Japan. Shore Beach 1980, 48, 30–35. [Google Scholar]
  50. Gisen, J.I.A.; Savenije, H.H. Estimating bankfull discharge and depth in ungauged estuaries. Water Resour. Res. 2015, 51, 2298–2316. [Google Scholar] [CrossRef] [Green Version]
  51. Jay, D.A.; Smith, J.D. Circulation, density distribution and neap-spring transitions in the Columbia River Estuary. Prog. Oceanogr. 1990, 25, 81–112. [Google Scholar] [CrossRef]
  52. De Brauwere, A.; De Brye, B.; Blaise, S.; Deleersnijder, E. Residence time, exposure time and connectivity in the Scheldt Estuary. J. Mar. Syst. 2011, 84, 85–95. [Google Scholar] [CrossRef]
  53. Vroom, J.; De Vet, L.; Van der Werf, J. Validatie Waterbeweging Delft3D-NeVla model Westerscheldemonding; Technical Report, Rapport 1210301-001-ZKS-0001; Deltares: The Nederland, 2015. [Google Scholar]
  54. Nguyen, A.; Savenije, H. Salt intrusion in multi-channel estuaries: A case study in the Mekong Delta, Vietnam. Hydrol. Earth Syst. Sci. Discuss. 2006, 10, 743–754. [Google Scholar] [CrossRef]
  55. Nowacki, D.J.; Ogston, A.S.; Nittrouer, C.A.; Fricke, A.T.; Van, P.D.T. Sediment dynamics in the lower Mekong River: Transition from tidal river to estuary. J. Geophys. Res. Oceans 2015, 120, 6363–6383. [Google Scholar] [CrossRef] [Green Version]
  56. Xing, F.; Meselhe, E.; Allison, M.; Weathers, H., III. Analysis and numerical modeling of the flow and sand dynamics in the lower Song Hau channel, Mekong Delta. Cont. Shelf Res. 2017, 147, 62–77. [Google Scholar] [CrossRef]
  57. Van Maanen, B.; Coco, G.; Bryan, K.R. A numerical model to simulate the formation and subsequent evolution of tidal channel networks. Aust. J. Civ. Eng. 2011, 9, 61–72. [Google Scholar] [CrossRef]
  58. Dissanayake, D.; Roelvink, J.; Van der Wegen, M. Modelled channel patterns in a schematized tidal inlet. Coast. Eng. 2009, 56, 1069–1083. [Google Scholar] [CrossRef]
  59. Dame, R. Estuaries. In Encyclopedia of Ecology; Jörgensen, S.E., Fath, B.D., Eds.; Academic Press: Oxford, UK, 2008; pp. 1407–1413. [Google Scholar]
  60. Kemp, W.M.; Batleson, R.; Bergstrom, P.; Carter, V.; Gallegos, C.L.; Hunley, W.; Karrh, L.; Koch, E.W.; Landwehr, J.M.; Moore, K.A.; et al. Habitat requirements for submerged aquatic vegetation in Chesapeake Bay: Water quality, light regime, and physical-chemical factors. Estuaries 2004, 27, 363–377. [Google Scholar] [CrossRef] [Green Version]
  61. Vinagre, C.; Fonseca, V.; Cabral, H.; Costa, M.J. Habitat suitability index models for the juvenile soles, Solea solea and Solea senegalensis, in the Tagus estuary: Defining variables for species management. Fish. Res. 2006, 82, 140–149. [Google Scholar] [CrossRef]
  62. Barnes, T.; Volety, A.; Chartier, K.; Mazzotti, F.; Pearlstine, L. A habitat suitability index model for the eastern oyster (Crassostrea virginica), a tool for restoration of the Caloosahatchee Estuary, Florida. J. Shellfish Res. 2007, 26, 949–959. [Google Scholar] [CrossRef]
  63. Degraer, S.; Verfaillie, E.; Willems, W.; Adriaens, E.; Vincx, M.; Van Lancker, V. Habitat suitability modelling as a mapping tool for macrobenthic communities: An example from the Belgian part of the North Sea. Cont. Shelf Res. 2008, 28, 369–379. [Google Scholar] [CrossRef] [Green Version]
  64. Feyrer, F.; Newman, K.; Nobriga, M.; Sommer, T. Modeling the effects of future outflow on the abiotic habitat of an imperiled estuarine fish. Estuar. Coasts 2011, 34, 120–128. [Google Scholar] [CrossRef]
  65. Lokhorst, I.; Braat, L.; Leuven, J.R.F.W.; Baar, A.W.; van Oorschot, M.; Selaković, S.; Kleinhans, M.G. Morphological effects of vegetation on the fluvial-tidal transition in Holocene estuaries. Earth Surf. Dyn. Discuss. 2018, 2018, 1–28. [Google Scholar] [CrossRef] [Green Version]
  66. Lotze, H.K.; Lenihan, H.S.; Bourque, B.J.; Bradbury, R.H.; Cooke, R.G.; Kay, M.C.; Kidwell, S.M.; Kirby, M.X.; Peterson, C.H.; Jackson, J.B. Depletion, degradation, and recovery potential of estuaries and coastal seas. Science 2006, 312, 1806–1809. [Google Scholar] [CrossRef] [PubMed]
  67. Yap, W.Y.; Lam, J.S.L. 80 million-twenty-foot-equivalent-unit container port? Sustainability issues in port and coastal development. Ocean Coast. Manag. 2013, 71, 13–25. [Google Scholar] [CrossRef]
  68. Depreiter, D.; Sas, M.; Beirinckx, K.; Liek, G.J. Flexible Disposal Strategy: Monitoring as a Key to Understanding And Steering Environmental Responses to Dredging and Disposal in the Scheldt Estuary; Technical Report; International Marine & Dredging Consultants: Antwerp, Belgium, 2012. [Google Scholar]
  69. Vikolainen, V.; Bressers, H.; Lulofs, K. A shift toward building with nature in the dredging and port development industries: Managerial implications for projects in or near natura 2000 areas. Environ. Manag. 2014, 54, 3–13. [Google Scholar] [CrossRef] [PubMed]
  70. Vos, P.; de Koning, J.; Van Eerden, R. Landscape history of the Oer-IJ tidal system, Noord-Holland (The Netherlands). Neth. J. Geosci. 2015, 94, 295–332. [Google Scholar] [CrossRef]
  71. Pierik, H.; Cohen, K.; Stouthamer, E. A new GIS approach for reconstructing and mapping dynamic late Holocene coastal plain palaeogeography. Geomorphology 2016. [Google Scholar] [CrossRef]
  72. Wood, L.J. Predicting tidal sand reservoir architecture using data from modern and ancient depositional systems, integration of outcrop and modern analogs in reservoir modeling. AAPG Mem. 2004, 80, 45–66. [Google Scholar]
  73. Yoshida, S.; Johnson, H.D.; Pye, K.; Dixon, R.J. Transgressive changes from tidal estuarine to marine embayment depositional systems: The Lower Cretaceous Woburn Sands of southern England and comparison with Holocene analogs. AAPG Bull. 2004, 88, 1433–1460. [Google Scholar] [CrossRef]
  74. Jackson, M.D.; Yoshida, S.; Muggeridge, A.H.; Johnson, H.D. Three-dimensional reservoir characterization and flow simulation of heterolithic tidal sandstones. AAPG Bull. 2005, 89, 507–528. [Google Scholar] [CrossRef]
  75. Steel, R.J.; Plink-Bjorklund, P.; Aschoff, J. Tidal deposits of the Campanian Western Interior Seaway, Wyoming, Utah and Colorado, USA. In Principles of Tidal Sedimentology; Springer: Berlin, Germany, 2012; pp. 437–471. [Google Scholar] [CrossRef]
  76. Saha, S.; Burley, S.D.; Banerjee, S.; Ghosh, A.; Saraswati, P.K. The morphology and evolution of tidal sand bodies in the macrotidal Gulf of Khambhat, western India. Mar. Pet. Geol. 2016, 77, 714–730. [Google Scholar] [CrossRef]
  77. Gray, A.J. Saltmarshes: Morphodynamics, Conservation and Engineering Significance; Allen, J., Pye, K., Eds.; Cambridge University Press: Cambridge, UK, 1992; pp. 63–79. [Google Scholar]
  78. Thomson, A.; Huiskes, A.; Cox, R.; Wadsworth, R.; Boorman, L. Short-term vegetation succession and erosion identified by airborne remote sensing of Westerschelde salt marshes, The Netherlands. Int. J. Remote Sens. 2004, 25, 4151–4176. [Google Scholar] [CrossRef]
  79. Van der Wal, D.; Wielemaker-Van den Dool, A.; Herman, P.M. Spatial patterns, rates and mechanisms of saltmarsh cycles (Westerschelde, The Netherlands). Estuar. Coast. Shelf Sci. 2008, 76, 357–368. [Google Scholar] [CrossRef]
  80. Bolla Pittaluga, M.; Coco, G.; Kleinhans, M.G. A unified framework for stability of channel bifurcations in gravel and sand fluvial systems. Geophys. Res. Lett. 2015, 42, 7521–7536. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) bathymetry of the Western Scheldt in 2012 from the mouth (left) to the Dutch-Belgian border (right); (b) peak tidal flow velocities from a calibrated hydrodynamic model of Rijkswaterstaat (SCALWEST, [32]) for the bathymetry of the Western Scheldt in 2009 and (c) ecotope map of the Western Scheldt (2012). Bathymetry, hydrodynamic output and ecotope map were obtained from Rijkswaterstaat.
Figure 1. (a) bathymetry of the Western Scheldt in 2012 from the mouth (left) to the Dutch-Belgian border (right); (b) peak tidal flow velocities from a calibrated hydrodynamic model of Rijkswaterstaat (SCALWEST, [32]) for the bathymetry of the Western Scheldt in 2009 and (c) ecotope map of the Western Scheldt (2012). Bathymetry, hydrodynamic output and ecotope map were obtained from Rijkswaterstaat.
Remotesensing 10 01915 g001
Figure 2. Examples of the systems used in this study. (a) a delta branch of a distributary of the Saskatchewan River (Canada); (b) the Tran De branch of the Mekong Delta (Vietnam); (c) the Western Scheldt estuary (The Netherlands); (d) the Columbia River estuary (USA); (e) the Ameland inlet and tidal basin (The Netherlands).
Figure 2. Examples of the systems used in this study. (a) a delta branch of a distributary of the Saskatchewan River (Canada); (b) the Tran De branch of the Mekong Delta (Vietnam); (c) the Western Scheldt estuary (The Netherlands); (d) the Columbia River estuary (USA); (e) the Ameland inlet and tidal basin (The Netherlands).
Remotesensing 10 01915 g002
Figure 3. Flowchart of the calculation steps in the tool.
Figure 3. Flowchart of the calculation steps in the tool.
Remotesensing 10 01915 g003
Figure 4. Relation between peak tidal flow velocity and depth below high water level from the Western Scheldt model (Figure 1a,b). Histograms indicate occurrence within bins and colours in the plot become more saturated with increasing point density.
Figure 4. Relation between peak tidal flow velocity and depth below high water level from the Western Scheldt model (Figure 1a,b). Histograms indicate occurrence within bins and colours in the plot become more saturated with increasing point density.
Remotesensing 10 01915 g004
Figure 5. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for the Western Scheldt.
Figure 5. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for the Western Scheldt.
Remotesensing 10 01915 g005
Figure 6. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for the Columbia River.
Figure 6. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for the Columbia River.
Remotesensing 10 01915 g006
Figure 7. (a,b) predicted bed elevations against measured bed elevations taken from bathymetry for the Western Scheldt (a) and Columbia River (b); (c,d) predicted flow velocity from the tool against modelled flow velocity for numerical models of the Western Scheldt (c) and Columbia River (d). Filled circles indicate median values per cross-section, open circles indicate maximum values per cross-section. Dashed lines indicate a spread of a factor 2 around the x = y line. Colours indicate the position along the estuary, with blue being the mouth and yellow being the most landward location; (e) predicted salinity for the three predictors against measured salinity in the Western Scheldt. Solid lines indicate mean measured salinity per transect against predicted salinity and dashed lines indicate maximum and minimum measured values.
Figure 7. (a,b) predicted bed elevations against measured bed elevations taken from bathymetry for the Western Scheldt (a) and Columbia River (b); (c,d) predicted flow velocity from the tool against modelled flow velocity for numerical models of the Western Scheldt (c) and Columbia River (d). Filled circles indicate median values per cross-section, open circles indicate maximum values per cross-section. Dashed lines indicate a spread of a factor 2 around the x = y line. Colours indicate the position along the estuary, with blue being the mouth and yellow being the most landward location; (e) predicted salinity for the three predictors against measured salinity in the Western Scheldt. Solid lines indicate mean measured salinity per transect against predicted salinity and dashed lines indicate maximum and minimum measured values.
Remotesensing 10 01915 g007
Figure 8. Comparison between modelled [red] and predicted [blue] flow velocities along the Western Scheldt (a,b) and the Columbia River (c,d). Shaded area shows the range of peak velocities (a,c) and range of average velocities over half a tidal cycle (b,d) bounded by the 5 and 95 percentile values. The solid line indicates the median and dashed line maximum value per transect.
Figure 8. Comparison between modelled [red] and predicted [blue] flow velocities along the Western Scheldt (a,b) and the Columbia River (c,d). Shaded area shows the range of peak velocities (a,c) and range of average velocities over half a tidal cycle (b,d) bounded by the 5 and 95 percentile values. The solid line indicates the median and dashed line maximum value per transect.
Remotesensing 10 01915 g008
Figure 9. Mean velocity over a tidal cycle (left) and peak tidal flow velocity (right) as a function of local tidal prism (a,b), bed elevation (c,d) and numerically modelled values (e,f) for the Columbia River.
Figure 9. Mean velocity over a tidal cycle (left) and peak tidal flow velocity (right) as a function of local tidal prism (a,b), bed elevation (c,d) and numerically modelled values (e,f) for the Columbia River.
Remotesensing 10 01915 g009
Figure 10. Resulting predictions of local tidal prism for the Western Scheldt (a) and Columbia River (b). Red lines are model results, blue lines are tool output.
Figure 10. Resulting predictions of local tidal prism for the Western Scheldt (a) and Columbia River (b). Red lines are model results, blue lines are tool output.
Remotesensing 10 01915 g010
Figure 11. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for an ideal estuary based on dimensions of the Tran De branch of the Mekong river.
Figure 11. Resulting predictions of bed elevation with respect to the high water level, inundation duration, peak flow velocity and salinity for an ideal estuary based on dimensions of the Tran De branch of the Mekong river.
Remotesensing 10 01915 g011
Figure 12. Sensitivity of tool output for the Western Scheldt to input tidal amplitude, river discharge and width profile. The ratio between intertidal and subtidal area increases with tidal amplitude (a) and decreases with river discharge and adaptation to ideal shape (b,c). The ratio between low dynamic (velocity < 0.8 ms 1 ) intertidal area and high dynamic area has an optimum curve for tidal amplitude (d) and decreases with river discharge and adaptation to ideal shape (e,f). Peak flow velocity at the mouth of the estuary increases with tidal amplitude and river discharge (g,h) and decreases with adaptation to ideal shape (i).
Figure 12. Sensitivity of tool output for the Western Scheldt to input tidal amplitude, river discharge and width profile. The ratio between intertidal and subtidal area increases with tidal amplitude (a) and decreases with river discharge and adaptation to ideal shape (b,c). The ratio between low dynamic (velocity < 0.8 ms 1 ) intertidal area and high dynamic area has an optimum curve for tidal amplitude (d) and decreases with river discharge and adaptation to ideal shape (e,f). Peak flow velocity at the mouth of the estuary increases with tidal amplitude and river discharge (g,h) and decreases with adaptation to ideal shape (i).
Remotesensing 10 01915 g012
Figure 13. Potential habitat area (in green) as part of total width for the species in Table 1.
Figure 13. Potential habitat area (in green) as part of total width for the species in Table 1.
Remotesensing 10 01915 g013
Table 1. Habitat preferences of species used in this study based on MarLIN database.
Table 1. Habitat preferences of species used in this study based on MarLIN database.
NameDepth Range (-)Flow Velocity (ms 1 )Salinity (ppt)
Mud shrimp (Corophium volutator)Intertidal<0.52–40
Blow lug (Arenicola marina)Intertidal<4.018–40
Salt marsh species 1 (Salicornia)Intertidal above Mean Sea Level<0.518–40
Salt marsh species 2 (Spartina anglica)Intertidal above 0.5 a <1.018–40

Share and Cite

MDPI and ACS Style

Leuven, J.R.F.W.; Verhoeve, S.L.; Van Dijk, W.M.; Selaković, S.; Kleinhans, M.G. Empirical Assessment Tool for Bathymetry, Flow Velocity and Salinity in Estuaries Based on Tidal Amplitude and Remotely-Sensed Imagery. Remote Sens. 2018, 10, 1915. https://doi.org/10.3390/rs10121915

AMA Style

Leuven JRFW, Verhoeve SL, Van Dijk WM, Selaković S, Kleinhans MG. Empirical Assessment Tool for Bathymetry, Flow Velocity and Salinity in Estuaries Based on Tidal Amplitude and Remotely-Sensed Imagery. Remote Sensing. 2018; 10(12):1915. https://doi.org/10.3390/rs10121915

Chicago/Turabian Style

Leuven, Jasper R. F. W., Steye L. Verhoeve, Wout M. Van Dijk, Sanja Selaković, and Maarten G. Kleinhans. 2018. "Empirical Assessment Tool for Bathymetry, Flow Velocity and Salinity in Estuaries Based on Tidal Amplitude and Remotely-Sensed Imagery" Remote Sensing 10, no. 12: 1915. https://doi.org/10.3390/rs10121915

APA Style

Leuven, J. R. F. W., Verhoeve, S. L., Van Dijk, W. M., Selaković, S., & Kleinhans, M. G. (2018). Empirical Assessment Tool for Bathymetry, Flow Velocity and Salinity in Estuaries Based on Tidal Amplitude and Remotely-Sensed Imagery. Remote Sensing, 10(12), 1915. https://doi.org/10.3390/rs10121915

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop