The “cubage”technique (“kubatuur” in Dutch) is a relatively simple technique to calculate the hydrodynamics in an estuary (Smets, 1996; Plancke et aL, 2011). It requires only topobathymetric data and water levels at different stations to calculate discharges and crosssectional averaged flow velocities (see Figure 3) starting from the conservation of mass formula. In contrast to numerical process models, this technique does not require any flow resistance (“roughness”) coefficients to calibrate the model. Where the cubage technique uses mass conservation (which is an exact relationship) the error of the method is solely related to the error in water level and topobathymetric data. Hereafter some mathematical background on this technique is given.
Output parameters
The “cubage”technique calculates discharges in time between every 2 consecutive crosssections in the estuary. From this parameter following additional parameters can be derived:
(a) Maximum ebb and flood discharges
(b) Ebb and flood discharges averaged over one tidal cycle
(c) Ebb and flood tidal volumes
(d) Evolution in time of sectionaveraged velocity
(e) Maximum ebb and flood velocity
(f) Ebb and flood velocities averaged over one tidal cycle
ul>
Selected input data and implementation
Within the TIDEproject it was chosen to perform the cubage calculation for mean tidal conditions (only input data available for mean tidal conditions). In a first step the mean tidal parameters (HW, LW, time of rising and time of falling) were derived from measurements for the present situation for a station at the mouth, a station in the middle of the estuary, and a station at the upstream boundary of the estuary (see Table 3). Where the cubage technique requires a continuous time series of water levels, one period (or 2 periods for the Elbe and Humber) with several tidal cycles was chosen from the continuous measurements (see ‘selected tides’ Table 3) for which the low and high waters have the best agreement with the MHWL and MLWL values for the present situation for all 3 selected stations. For the Elbe and Humber 2 time periods were selected (one with low and one with high river discharge) because of the important differences in river discharge at the upestuary boundary of these estuaries. From the selected time period, water level measurements for all tidal stations were selected as input for the cubage calculation (location of all tidal stations see Figure 9 (Scheldt), Figure 10 (Elbe), Figure 11 (Weser) and Figure 12 (Humber)) .
In a next step cross sections were defined perpendicular to the thalweg. The thalweg is defined as the longitudinal profile linking the deepest points of each cross section. The distance between two sections was chosen at about 500 m: this resolution was on the one hand sufficient to represent the topobathymetry, while on the other hand the number of cross section remained acceptable. These cross sections were exported in a MIKE11format. Subsequently the different cross sections were imported in MIKE11 and the necessary parameters for each cross section were derived automatically within this software:
 Wet cross section area at different heights
 Width at different heights
At the upestuarine boundaries the “cubage” technique requires the implementation of a fresh water discharge. Based on the daily values of the present situation, two characteristic fresh water discharges were calculated: P(5%) representing a “summer” condition and P(95%) representing a “winter” condition.
Cubage calculations with a winter/summer condition were within this study performed for 3 of the 4 estuaries (Elbe, Weser and Humber) (see Plancke et al., 2012a,b,c). For the calculation of the high/low discharge period of the Elbe and Humber, appropriate tides were selected in the high/low discharge period with implementation of respectively the P(95%) and P(5%) freshwater discharges. For the Weser, appropriate tides were selected only during one time period of the year. The high/low discharge variation was again simulated by implementation of respectively the P(95%) and P(5%) freshwater discharges. For the Scheldt, an earlier cubage calculation with implementation of a mean freshwater discharge (P(50%)) was used (Plancke et al., 2011). Here, no distinction was thus made between winter and summer condition, due to the limited effect of the discharge on the water level.
More detailed information on the cubage technique is given by Plancke et al. (2012a,b,c).
Table 3 – Input data used for the cubage calculations of each estuary

Scheldt 
Elbe 
Weser 
Humber 
Topobathymetry 
2001 
2006 
2009 
2005 
HW and LW 
19912000 
20012010 
20012010 
2005 
Stations used for selection of mean tide 
Vlissingen
Antwerpen
SintAmands 
Groβer Vogelsand
Kollmar
Wehr Geesthacht 
LT Alte Weser
Nordenham
Grosse Weserbrücke 
Spurn Point
South Ferriby
Naburn Lock (Ouse)
Carlton on Trent (Trent) 
Selected tides 
2122/06/2009 
5/03/2006 (high Q)
25/05/2006 (low Q) 
78/06/2010 
56/01/2005 (high Q)
34/05/2005 (low Q) 
Freshwater discharge 
20002010
(P50) 
20012010
(P95 and P5) 
20012010
(P95 and P5) 
2010
(P95 and P5) 
Tidal damping scale
Introduction
The tidal penetration in an estuary is influenced by several factors, with the most important factors being the funnel shape of the estuary, leading to an amplification of the tidal range, and the friction within the estuary, leading to a decrease in tidal range. The relative importance of these 2 factors is presented in the tidal damping scale (), which is a simple analytical expression to describe tidal amplification or damping (Savenije, 2001). Tidal damping (1/ is defined as:
With:
ß = tidal damping scale [m]
b = width convergence length [m]; derived from following regression:
f’ = adjusted friction factor =
v = amplitude of flow velocity (average of maximum ebb and maximum flood velocity at mean tidal
conditions) [m/s]
e = phase lag between HW and slack HW
c = celerity of tidal wave (= sqrt(g.h)) [m/s]
h = averaged water depth over tidal cycle (= (hHW + hLW)/2) [m]
h = tidal amplitude (= tidal range/2) [m]
C = Chézy roughness coefficient [m1/2/s]
g = gravitation constant (9,81 m/s²)
In the case that the funneling of the estuary is more important than the friction in the estuary, 1/b will be positive, and tidal amplification can be expected. If friction becomes more important, 1/b will be negative and tidal damping can be expected.
Assumptions
To calculate 1/β the following assumptions were made:
 Data were averaged over 5 km blocks along the thalweg
 Only the estuary was regarded and not the mouth area (based on mouthgeo, for definition see §3.4)
 Data with a ratio TR/h > 0.80 were excluded from the analysis (TR =tidal range, h= averaged water depth over tidal cycle). Above this threshold, tidal damping calculations become less accurate (Savenije, 1998). For the Scheldt, 5 data points were excluded (in total 25 km), for the Humber 11 datapoints (in total 55 km, a consequence of the shallowness of the estuary). For Elbe and Weser no data points were excluded.
 Chézy roughness coefficient C = 50; phase lag difference ε = 50’ (see sensitivity analysis below)
Estuary convergence
In the tidal damping equation (see Equation (1)), the estuary convergence is described by the width convergence length (b). This parameter is determined based on the following equation:
With:
B = the average crosssectional and tidalaverage width (m)
B^{0} = the width at the estuary mouth (m)
x = the distance from the mouth
The slope of the relationship between ln(B) and x then returns 1/b (see Appendix A). The larger 1/b, the more convergent the estuary is. Classically, the b value is determined for the whole estuary (Savenije, 1998), however we observed some important breaks in ln(B) for 3 out of 4 estuaries (Scheldt, Elbe and Humber). The best example is the HumberOuse estuary where the Humber is clearly much more convergent than the Ouse tributary (see Figure A 4). When defining breaks, we assumed that estuary convergence has mainly an effect on tidal amplification/damping over a large distance along the estuary (> 20 km). Smaller, more local variations in ln(B) were not considered. The found 1/b values are listed in Table 4.
Sensitivity analysis
For the Scheldt estuary a sensitivity analysis was performed on the tidal damping scale (Figure 5 ) in relation to the following parameters:
 Chézy roughness coefficient: 40 – 50 – 60
 Phase lag HW – slack HW: 30’ – 50’ – 70’
These values were selected since they are common values in analysis for the Scheldt estuary (Savenije, 2001, 2005; van Rijn, 2011; Winterwerp, 2012)
It was decided to use C = 50 and ε = 50’ to calculate 1/b based on the following criteria:
 There should be a good agreement between 1/b and the change in tidal range: if 1/b > 0 we can expect an increase in tidal range, if 1/b < 0 we can expect a decrease in tidal range
 The selected values for C and ε should be realistic values. A Chézy coefficient of 50 is a commonly used value in estuarine modeling (Savenije, 2001, 2005; van Rijn, 2011; Winterwerp, 2012). However, it should be pointed out that Chézy coefficients may vary in time for one estuary, and in between different estuaries (from 3080, see Winterwerp (2012)). Under mean tidal conditions, a phase lag difference of 50 minutes is representative for estuaries with semidiurnal tides (Savenije, 2001)
For the calculation of 1/b , some of the parameters are variable along the length axis of the estuary, others or constant. An overview is given in Table 4.
Table 4 – Overview of the parameters used to calculate 1/b
Parameter 
Scheldt 
Elbe 
Weser 
Humber 





1/b [m1] 
2.42*E5 (TIDE km 052)
3.74*E5 (TIDE km 52157)
3.39*E5 (mean) 
1.16*E5 (TIDE km 024)
3.17*E5 (TIDE km 24149)
2.75*E5 (mean) 
2.79*E5 
2.17*E5 (TIDE km 044)
5.4*E5 (TIDE km 44114)
4.18*E5 (mean) 
v [m/s] 
Calculated for each “cubage” crosssection (based on cubage) 
e [minutes] 
50 
50 
50 
50 
c [m/s] 
Calculated for each “cubage” crosssection 
h [m] 
Calculated for each “cubage” crosssection 
h [m] 
Calculated for each “cubage” crosssection (based on measured data) 
C [m1/2/s] 
50 
50 
50 
50 
Sediment fluxes
Introduction
Based on the derived tidal discharges calculated by applying the cubage technique (§3.2.1, see output parameters (b)) and the observed SPM values (see §3.1.4), sediment fluxes can be calculated according to the following equation:
With:
SF = sediment flux per unit cross section (kg/m²/s)
Cs = concentration suspended particle matter (kg/m³)
Q = tidal discharge over one tidal cycle at mean tidal conditions, positive from upestuary to downestuary (ebb = +  flood = ) (m³/s)
A = wet section at mean tidal conditions during half tide (= mean of wet section at low water and high water) (m²)
Implementation
Tidal discharges over a tidal cycle, representative for mean tidal conditions, were calculated at every water level station (§3.1.2). For the Scheldt these tidal discharges were calculated under conditions of mean riverine discharges (see Plancke et al., 2011). For the Elbe, Weser and Humber tidal discharges were calculated under conditions of high and low riverine discharges (see Plancke et al., 2012a,b,c). If the netto discharge over a tidal cycle is positive than the water is transported in the ebb direction, if the netto discharge over a tidal cycle is negative than the water is transported in the flood direction. SPM values at the water level stations were derived by use of linear interpolation (SPM stations not at the same locations as the water level stations, see Figure 9 till Figure 12). Based on the available SPM data and calculated tidal discharges (see Table 5), sediment fluxes were calculated according to equation (3). The sediment fluxes in the flood direction are hereby positive and the sediment fluxes in the ebb direction are negative (see equation (3)). Sediment fluxes under mean riverine discharges for the Elbe, Weser and Humber were determined by averaging the calculated sediment fluxes under high and low riverine discharge.
For each water level stations a netto sediment flux was calculated, as the sum of the different sediment fluxes per time step over a full tidal cycle.
Table 5 – Overview of the input parameters used to calculate sediment fluxes. (R = riverine discharge)

SPM 
Tidal discharges 
Scheldt 
surface and watercolumn 
mean R 
Elbe 
surface 
low and high R 
Weser 
surface 
low and high R 
Humber 
surface 
low and high R 
It should be mentioned that this methodology is only a simple approximation of the real sediment fluxes. The methodology uses a constant SPM over the tidal cycle (due to limited data availability), while in reality the SPM will vary strongly over a tidal cycle (Figure 6). Therefore the calculated fluxes should be seen as a first approximation that will be used in the evaluation of topic 3.
With:
ρ = density of water = 1000 kg/m³
V_{flood} = flood volume at mean tidal conditions (m³)
v_{flood} = mean flow velocity during flood at mean tidal conditions (m/s)
g = gravitation constant (9.81 m²/s)
MHW = mean high water level (m)
MLW = mean low water level (m)
Q = freshwater discharge (m³/s)
T = duration of a tidal cycle (s)
A_{ebb} = wet section at low water (m²)
h_{topo} = thalweg elevation relative to MLW at the most downstream cubage crosssection (m)
Implementation
Energy fluxes were calculated for every defined “cubage” crosssection. At these crosssections, the necessary input parameters for the energy calculations, are the output parameters of the cubage calculation (see (§3.2.1). As for the cubage calculations of the Elbe, Weser and Humber, energy fluxes were calculated for a winter and summer condition (using respectively P(95%) and P(5%) values for the freshwater discharge). Within this study, an earlier cubage calculation was used for the Scheldt using mean freshwater discharge (Plancke et al., 2011). Consequently, calculated energy fluxes are only representative for mean freshwater discharges. To make energy fluxes comparable within the interestuarine comparison (§5.3), mean fluxes were calculated for the Elbe, Weser and Humber by averaging the winter and summer energy fluxes.
Residence time
Introduction
Different simplified methods can be used to calculate residence times. A comparison between the different methodologies was made by Ides (2011b). It was found that for estuaries with a strong salinity gradient, the most appropriate method to calculate residence times is the fractional fresh water method (Dyer, 1973). As all estuaries within the TIDE project have a strong salinity gradient, the fractional fresh water method was applied, in which the residence time is defined as:
With:
T = residence time (days)
f = the fractional fresh water concentration = (S_{s} – S_{n})/S_{s}
S_{s} = salinity at the most downstream cubage crosssection (PSU)
S_{n} = salinity in a given segment of the estuary (PSU)
V = volume of the estuary segment (m³)
R = river discharge (m³/s)
Implementation
The defined estuary segments were based on the cubage crosssections (see Figure 3). For each estuary segment a fractional fresh water concentration was attributed based on the corresponding salinity profile. The segment volumes were calculated based on the mean crosssectional area multiplied by the segment length. The mean crosssectional area is hereby calculated as the average of the crosssectional areas at MHWL and MLWL at both ends of the segment. River discharges for the Elbe and Weser are equal to the freshwater discharge of the main channel. Since the discharges of the tributaries are approximately 100 times smaller than the main channel discharge (see also §3.1.5), the river discharges along the estuary segments are constant. For the Scheldt and Humber this is not the case. There, each tributary has an important contribution to the river discharge (§3.1.5). As a consequence, the river discharge for the different estuary segments increases to the mouth by adding the discharges of the different tributaries. Residence times along the estuaries are eventually found by summing up the calculated residence times of the individual estuary segments.