3.1 Tidal waves

The tide causes the water surface elevation $z_{s}$ and discharge $Q$ to periodically oscillate over time $t$ , which suggests separating them into the components $z_{j}=\text{Re}\{z_{j}\}$ and $Q_{j}=\text{Re}\{Q_{j}\}$ with frequencies $\unicode[STIX]{x1D714}_{j}$ (Godin Reference Godin1991b ),

(3.1a ) $$\begin{eqnarray}\displaystyle z_{s}(t,x) & = & \displaystyle \mathop{\sum }_{j=0}^{\infty }z_{j}(t,x)=z_{0}(x)+\mathop{\sum }_{j=1}^{\infty }\text{Re}\{z_{j}(x)\exp (\text{i}\unicode[STIX]{x1D714}_{j}t)\},\end{eqnarray}$$

(3.1b ) $$\begin{eqnarray}\displaystyle Q(t,x) & = & \displaystyle \mathop{\sum }_{j=0}^{\infty }Q(t,x)=Q_{0}+\mathop{\sum }_{j=1}^{\infty }\text{Re}\{Q_{j}(x)\exp (\text{i}\unicode[STIX]{x1D714}_{j}t)\}.\end{eqnarray}$$

$x$

$j$

$(y)_{j}=(1/T)\int _{T}y(\cos \unicode[STIX]{x1D714}_{j}t+\text{i}\sin \unicode[STIX]{x1D714}_{j}t)\,\text{d}t$

$T$

The surface elevation of each frequency component has a distinct amplitude $|z_{j}|$ and phase $\unicode[STIX]{x1D711}_{jz}=\arctan (\text{Im}(z_{j})/\text{Re}\{z_{j}\})$ . The astronomical tide consists of an infinite number of constituents (Pugh Reference Pugh1987). Their frequencies are integer combinations of basic frequencies derived from the orbits of the celestial bodies (Doodson Reference Doodson1921; Cartwright & Tayler Reference Cartwright and Tayler1971; Souchay, Mathis & Tokieda Reference Souchay, Mathis and Tokieda2012). Several dozen constituents are required to accurately predict ocean tides, of which many constituents are of similar frequency and magnitude. Long time series are required to separate these constituents from each other.

River discharge not only determines the means $z_{0}$ and $Q_{0}$ ( $\unicode[STIX]{x1D714}_{0}=0$ ) but also modulates the tide. River discharge varies in an irregular manner over much shorter periods than necessary for a meaningful harmonic analysis. Therefore, we consider the tide for successive periods of just one tidal cycle and decompose it into a Fourier series, where the frequencies of the components are integer multiples of a single fundamental frequency, $\unicode[STIX]{x1D714}_{j}=j\unicode[STIX]{x1D714}_{1}$ . The Fourier components decay rapidly in amplitude, which allows a meaningful truncation of the series to just a few components. These components are referred to as tidal species and effectively lump tidal constituents of similar frequencies together (Kukulka & Jay Reference Kukulka and Jay2003b ; Guo et al. Reference Guo, van der Wegen, Jay, Matte, Wang, Roelvink and He2015). Alternatively to species, the tidal wave can be interpreted as a periodic function of arbitrary shape that is described by low water and high water (Savenije Reference Savenije2001; Savenije et al. Reference Savenije, Toffolon, Haas and Veling2008). This approach is supported by the observation that tidal waves travel upstream individually after each other and that the discharge of large rivers changes little over the time it takes a single wave to travel upstream. The incoming tide is thus roughly represented by a single frequency component that has an amplitude equal to half the tidal range. The range of the incoming tide changes from one tidal cycle to the next, most notably over to the spring-neap cycle. The river tide has therefore be predicted for each cycle individually, depending on the incoming tide and the river flow. The amplitude and phase of a wave change as a wave propagates up river, depending on the cross-section geometry and river discharge. For convenience, this is expressed in the form of the admittance $|z_{j}(x)|/|z_{j}(0)|$ and phase difference $\unicode[STIX]{x1D711}_{jz}(x)-\unicode[STIX]{x1D711}_{jz}(0)$ , where $x$ is the distance from the river mouth. The remainder of this section develops the theory of tidal wave propagation. It builds on previous works by Godin (Reference Godin1985) and Jay (Reference Jay1991). This section advances the theory on how tides propagate along rivers with varying cross-section geometry. The theory is held general and covers both mild depth and width convergence. Section 4 then analyses the backwater effect caused by a varying river discharge and a sloping bed.

3.2 Shallow-water equations

The flow in open channels is described by the one-dimensional shallow-water equations (Cunge, Holly & Verwey Reference Cunge, Holly and Verwey1980; Savenije Reference Savenije2012). These are the equation of continuity,

(3.2a ) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}A}{\unicode[STIX]{x2202}t}+\frac{\unicode[STIX]{x2202}Q}{\unicode[STIX]{x2202}x}=0,\end{eqnarray}$$

as well as the equation of motion,

(3.2b ) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}Q}{\unicode[STIX]{x2202}t}+\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}x}\left(\frac{Q^{2}}{A}\right)+\frac{1}{2}\frac{g}{w}\frac{\unicode[STIX]{x2202}A^{2}}{\unicode[STIX]{x2202}x}=-gA\frac{\unicode[STIX]{x2202}z_{b}}{\unicode[STIX]{x2202}x}+g\frac{A^{2}}{w^{2}}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}-c_{d}w\frac{Q|Q|}{A^{2}},\end{eqnarray}$$

where $A$ is the cross-sectional area, $Q$ the discharge, $w$ the channel width, $z_{b}$ the bed level, $g$ the gravitational acceleration and $c_{d}$ the drag coefficient. We analyse here only the case of a straight non-meandering channel that has a rectangular cross-section, i.e. no intertidal areas (figure 1), so that the depth $h=A/w$ and the surface elevation $z_{s}=z_{b}+h$ . The terms on the right-hand side in (3.2b ) represent the forces acting on the flow per unit distance along the channel. The forces determine how the tidal wave changes while propagating up river. Tidal flats are not taken into account, as intertidal areas are small in rivers. The reader is referred to Speer & Aubrey (Reference Speer and Aubrey1985), Jay (Reference Jay1991), Friedrichs & Madsen (Reference Friedrichs and Madsen1992), Savenije et al. (Reference Savenije, Toffolon, Haas and Veling2008) for the treatment of intertidal storage, and for tidal propagation along channels of arbitrary cross-sections, to Li & Valle-Levinson (Reference Li and Valle-Levinson1999). We assume that the channel is wide enough so that the hydraulic radius is well approximated by the water depth and narrow enough for Rossby circulation to be relatively small. We also neglect spatio-temporal variation of the drag coefficient $c_{d}$ between high and low river flow as well as between flood and ebb flow.

3.3 Wave equation

As the tide is a periodic function, it is purposeful to decompose the shallow-water equations into their frequency components. The equations are coupled by the interaction of the species due to the nonlinear terms. To transform the shallow-water equations into the wave equation, we consider the case where the tidal amplitude is small compared to the tidally averaged water depth $h_{0}=(1/T_{1})\int _{0}^{T_{1}}h\,\text{d}t$ ; $T_{1}=2\unicode[STIX]{x03C0}/\unicode[STIX]{x1D714}_{1}$ is the tidal period. We neglect the small effect of nonlinearity in $1/h$ , which has been discussed in the literature (Godin Reference Godin1985). We also neglect the advective acceleration term $(\unicode[STIX]{x2202}/\unicode[STIX]{x2202}x)(Q^{2}/A)$ because its magnitude is small (Savenije Reference Savenije2012). This holds as long as $h_{0}>|z_{1}|$ , which is the case as long as the bed slope is moderate or the river flow is strong.

For the mean flow $\unicode[STIX]{x1D714}_{j}=0$ , continuity is trivial ( $\unicode[STIX]{x2202}Q_{0}/\unicode[STIX]{x2202}t=\unicode[STIX]{x2202}Q_{0}/\unicode[STIX]{x2202}x=0$ ), and the momentum equation simplifies to the backwater equation that determines the tidally averaged water level $z_{0}$ ,

(3.3) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}z_{0}}{\unicode[STIX]{x2202}x}+\frac{c_{d}w}{\unicode[STIX]{x03C0}gA_{0}^{3}}F_{0}=0,\quad \unicode[STIX]{x1D714}_{j}=0,\end{eqnarray}$$

where $A_{0}=wh_{0}$ is the tidally averaged cross-sectional area, and $F_{0}$ is the mean component of $(1/\unicode[STIX]{x03C0})F=|Q|Q$ , the signed square of the friction term.

The oscillatory components ( $\unicode[STIX]{x1D714}_{j}>0$ ) are determined by the wave equation. We obtain the wave equation by first differentiating the continuity equation in time and the momentum equation in space and then eliminating the surface elevation $z$ by combining the equations,

(3.4) $$\begin{eqnarray}\displaystyle -\frac{1}{gh_{0}w}\frac{\unicode[STIX]{x2202}^{2}Q}{\unicode[STIX]{x2202}t^{2}}+\frac{1}{w}\frac{\unicode[STIX]{x2202}^{2}Q}{\unicode[STIX]{x2202}x^{2}}+\frac{1}{w^{2}}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}\frac{\unicode[STIX]{x2202}Q}{\unicode[STIX]{x2202}x}-\frac{c_{d}}{g{h_{0}}^{3}w^{2}}\frac{\unicode[STIX]{x2202}Q|Q|}{\unicode[STIX]{x2202}t}=0. & & \displaystyle\end{eqnarray}$$

We approximate the signed square of the friction term with a quadratic Chebyshev polynomial (Dronkers Reference Dronkers1964),

(3.5) $$\begin{eqnarray}\displaystyle & \displaystyle \frac{1}{\unicode[STIX]{x03C0}}F\approx |Q|Q, & \displaystyle\end{eqnarray}$$

(3.6) $$\begin{eqnarray}\displaystyle & \displaystyle F=f_{0}Q_{hr}^{2}+f_{1}Q_{hr}Q+f_{2}Q^{2}, & \displaystyle\end{eqnarray}$$

where $Q_{hr}$ is half the tidal range. The complex conjugate is indicated by the asterisk: $f_{0,1,2}$ are coefficients that depend on the relative strength of the river and the tidal flow; $f_{0}$ is always small. When the river flow is low so that $Q_{0}<Q_{hr}$ , then $f_{1}=8/2$ and $f_{2}$ is small. When the river flow is strong so that $Q_{0}\geqslant Q_{hr}$ , then $f_{1}=0$ and $f_{2}$ equals $\unicode[STIX]{x03C0}$ . Appendix C gives the detailed expressions for $f_{0,1,2}$ .

The expansion of the discharge as a Fourier series (3.1b ), yields one equation for each frequency component. By continuity, the tidal discharge is proportional to its derivative with respect to time so that the wave equation reduces to a second-order ordinary differential equation (Ippen Reference Ippen1966). As the surface elevation has been eliminated, the system consists only of one equation per frequency component,

(3.7) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}^{2}Q_{j}}{\unicode[STIX]{x2202}t^{2}}+g\frac{A_{0}}{w^{2}}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}\frac{\unicode[STIX]{x2202}Q_{j}}{\unicode[STIX]{x2202}x}-g\frac{A_{0}}{w}\frac{\unicode[STIX]{x2202}^{2}Q_{j}}{\unicode[STIX]{x2202}x^{2}}+\frac{c_{d}w}{\unicode[STIX]{x03C0}A_{0}^{2}}F_{j}^{\prime }=0,\quad \unicode[STIX]{x1D714}_{j}>0,\end{eqnarray}$$

where $(1/\unicode[STIX]{x03C0})F_{j}^{\prime }$ are the frequency components of $(\unicode[STIX]{x2202}/\unicode[STIX]{x2202}t)(|Q|Q)$ .

With the Chebyshev approximation, the frequency components $F_{j}$ and $F_{j}^{\prime }$ for $\unicode[STIX]{x1D714}_{j}=0$ , $\unicode[STIX]{x1D714}_{j}=\unicode[STIX]{x1D714}_{1}$ and $\unicode[STIX]{x1D714}_{j}=2\unicode[STIX]{x1D714}_{1}$ are

(3.8a ) $$\begin{eqnarray}\displaystyle F_{0} & = & \displaystyle f_{0}Q_{hr}^{2}+f_{1}Q_{0}Q_{hr}+f_{2}(Q_{0}|Q_{0}|+{\textstyle \frac{1}{2}}(|Q_{1}|^{2}+|Q_{2}|^{2})),\end{eqnarray}$$

(3.8b ) $$\begin{eqnarray}\displaystyle F_{1}^{\prime } & = & \displaystyle \text{i}\unicode[STIX]{x1D714}_{1}((f_{1}Q_{hr}+2f_{2}Q_{0})Q_{1}+f_{2}Q_{2}Q_{1}^{\ast }),\end{eqnarray}$$

(3.8c ) $$\begin{eqnarray}\displaystyle F_{2}^{\prime } & = & \displaystyle \text{i}\unicode[STIX]{x1D714}_{2}((f_{1}Q_{hr}+2f_{2}Q_{0})Q_{2}+{\textstyle \frac{1}{2}}f_{2}Q_{1}^{2}),\end{eqnarray}$$

$\unicode[STIX]{x1D714}_{1}$

Further analysis is limited to two frequency components, representing the main tidal species. For the main tidal species, we use the shorthand notation

(3.9a ) $$\begin{eqnarray}c_{2}\frac{\unicode[STIX]{x2202}^{2}Q_{j}}{\unicode[STIX]{x2202}x^{2}}+c_{1}\frac{\unicode[STIX]{x2202}Q_{j}}{\unicode[STIX]{x2202}x}+c_{0}Q_{j}=0,\end{eqnarray}$$

with

(3.9b ) $$\begin{eqnarray}\displaystyle \frac{c_{1}}{c_{2}} & = & \displaystyle -\frac{1}{w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x},\end{eqnarray}$$

(3.9c ) $$\begin{eqnarray}\displaystyle \frac{c_{0}}{c_{2}} & = & \displaystyle \frac{\unicode[STIX]{x1D714}_{1}^{2}}{gh_{0}}-\frac{\text{i}\unicode[STIX]{x1D714}_{1}c_{d}}{\unicode[STIX]{x03C0}wgh_{0}^{3}}(f_{1}Q_{hr}+2f_{2}Q_{0}),\end{eqnarray}$$

where we consider the case in which the magnitude of the overtide is small so that its feedback on the main tidal species through $f_{2}Q_{2}Q_{1}^{\ast }$ can be neglected. As the frequency components are trigonometric functions in time (cf. (3.1b )), they are proportional to their derivative $\unicode[STIX]{x2202}z_{j}/\unicode[STIX]{x2202}t=\text{i}\unicode[STIX]{x1D714}_{j}z_{j}$ . The surface elevation amplitude of each component can thus be determined by differentiating the discharge along $x$ .

Substitution of the tidal average of the friction term (3.8a ) into the backwater equation (3.3) yields $h_{0}\approx z_{b}+(f_{1}Q_{0}Q_{hr}+f_{2}Q_{0}|Q_{0}|)x$ near the sea, which shows that the water surface slope increases linearly with the river discharge when the river discharge is low and quadratically when it is high. Conversely, the frequency component of the friction term that corresponds to the main tidal species (3.8b ) increases linearly with the tidal discharge when the tidal discharge is low and quadratically when it is high.

As the friction term is nonlinear, it couples the equations between the frequency components. The friction term damps and delays the tide. In addition, it generates components of higher frequency, the overtide (Parker Reference Parker1991). The overtide changes the shape of the tidal wave as it propagates up river (Parker Reference Parker1991). River flow forces an overtide with twice the frequency of the incoming tide so that high water is advanced and low water is delayed (Godin Reference Godin1999). The overtide is different in estuaries with wide tidal flats, where the falling limb of the tide lasts shorter than the rising limb (Friedrichs & Madsen Reference Friedrichs and Madsen1992).

Similarly, the friction term generates lower frequency components when the incoming tide contains components with close frequencies (LeBlond Reference LeBlond1979; Buschman et al. Reference Buschman, Hoitink, van der Vegt and Hoekstra2009). These modulate the daily mean water level over the spring-neap cycle. Subtidal variations of the surface elevation are captured by the $Q_{hr}$ -terms in (3.8a ). Modelling of subtidal harmonics is discussed in Kukulka & Jay (Reference Kukulka and Jay2003a ). The overtide and subtidal harmonics are small in magnitude so that we ignore their feedback on the main tidal component in further analysis.

3.4 Propagation of tidal waves

The discharge and tidal amplitude can be expressed as the product of the initial values $Q_{j}(0),z_{j}(0)$ at the river mouth and a complex admittance factor that we define as

(3.10a ) $$\begin{eqnarray}\displaystyle z_{j}(x)=z_{j}(0)\exp \left(-\text{i}\int _{0}^{x}k_{jz}\,\text{d}x^{\prime }\right), & & \displaystyle\end{eqnarray}$$

(3.10b ) $$\begin{eqnarray}\displaystyle Q_{j}(x)=Q_{j}(0)\exp \left(-\text{i}\int _{0}^{x}k_{jQ}\,\text{d}x^{\prime }\right). & & \displaystyle\end{eqnarray}$$

$\{z_{1},Q_{1}\}$

$k_{1Q}$

$k_{1z}$

(3.11a ) $$\begin{eqnarray}\displaystyle \frac{1}{Q_{1}}\frac{\unicode[STIX]{x2202}Q_{1}}{\unicode[STIX]{x2202}x} & = & \displaystyle -\text{i}k_{1Q},\end{eqnarray}$$

(3.11b ) $$\begin{eqnarray}\displaystyle \frac{1}{z_{1}}\frac{\unicode[STIX]{x2202}z_{1}}{\unicode[STIX]{x2202}x} & = & \displaystyle -\text{i}k_{1z}=\left(\frac{1}{k_{1Q}}\frac{\unicode[STIX]{x2202}k_{1Q}}{\unicode[STIX]{x2202}x}-\frac{1}{w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}-\text{i}k_{1Q}\right).\end{eqnarray}$$

$k_{1Q}=k_{1z}=k_{1}$

(3.12a ) $$\begin{eqnarray}\displaystyle z_{j}(t,x) & = & \displaystyle z_{j}(t,0)\exp (\text{i}(\unicode[STIX]{x1D714}_{1}t-k_{1}x)),\end{eqnarray}$$

(3.12b ) $$\begin{eqnarray}\displaystyle Q_{j}(t,x) & = & \displaystyle Q_{j}(t,0)\exp (\text{i}(\unicode[STIX]{x1D714}_{1}t-k_{1}x)).\end{eqnarray}$$

$\text{Re}\{\exp (\text{i}(\unicode[STIX]{x1D714}_{1}t-k_{1}x))\}=\exp (\text{Im}\{k_{1}\}x)\cos (\unicode[STIX]{x1D714}_{1}t-\text{Re}\{k_{1}\}x)$

$\text{Re}\{k_{1}\}$

$\text{Im}(k_{1})$

$k_{1}$

$1/100~\text{km}^{-1}$

$k_{1Q}$

$k_{1z}$

$k_{1Q}$

The wave equation (3.9a ) can be separated into two first-order ordinary differential equations when (3.11a ) is inserted into (3.9a ). This yields a Riccati equation, from which $k_{1Q}$ can be obtained,

(3.13a ) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}k_{1Q}}{\unicode[STIX]{x2202}x}=\text{i}k_{1Q}^{2}-\frac{c_{1}}{c_{2}}k_{1Q}-\text{i}\frac{c_{0}}{c_{2}}.\end{eqnarray}$$

Far upstream, the flow is uniform and $k_{1Q}$ does not change along the channel. The left-hand side of (3.13a ) is then zero so that $-\text{i}k_{1Q}$ is a root of the characteristic polynomial $c_{2}r^{2}+c_{1}r+c_{0}=0$ . The roots of the characteristic polynomial are

(3.13b ) $$\begin{eqnarray}r^{\pm }(x)=\frac{\text{i}}{2w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}\pm \sqrt{\frac{\unicode[STIX]{x1D714}^{2}}{gh}-\frac{1}{4w^{2}}\left(\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}\right)^{2}+\frac{\text{i}c_{d}\unicode[STIX]{x1D714}}{\unicode[STIX]{x03C0}wgh^{3}}(2f_{2}Q_{0}+f_{1}Q_{1})}.\end{eqnarray}$$

For a river of constant width and depth, the roots remain constant along the channel. In this case, the roots have well-known limits for the case of no river flow: $r^{2}=-\unicode[STIX]{x1D714}^{2}/gh+\text{i}(8/3\unicode[STIX]{x03C0})c_{d}\unicode[STIX]{x1D714}(Q_{hr}/gwh^{3})$ , (Lorentz Reference Lorentz1926) and when the river flow is strong, $r^{2}=2\text{i}\unicode[STIX]{x1D714}c_{d}(Q_{0}/gwh_{0}^{3})$ (Godin Reference Godin1985). When the width changes along the channel, the roots can still remain constant as long as the change is exponential.

Downstream, where both width and depth converge, $k_{1Q}$ changes along the river. Thus, $k_{1Q}$ can be determined by integrating the initial value problem (3.13a ) from upstream to downstream. In general, there is no closed-form solution to this initial value problem. Further simplifications are necessary to determine how the tide propagates up river.

3.5 Wave propagation along rivers with a gradually varying cross-section

The solution to the wave equation is the superposition of two waves. One wave travels upstream, and the other one travels downstream. These are analogous to the Riemann invariants of the shallow-water equations. In the case of constant coefficients, which holds in channels of constant cross-section,

(3.14a ) $$\begin{eqnarray}\displaystyle Q_{1}(x) & = & \displaystyle Q_{1}^{+}(0)\exp (r^{+}x)+Q_{1}^{-}(0)\exp (r^{-}x),\end{eqnarray}$$

(3.14b ) $$\begin{eqnarray}\displaystyle z_{1}(x) & = & \displaystyle z_{1}^{+}(0)\exp (r^{+}x)+z_{1}^{-}(0)\exp (r^{-}x),\end{eqnarray}$$

$r^{\pm }$

$k_{1}$

$k_{1}=+\text{i}r^{-}$

$Q_{1}^{+}(x)=0$

$z_{1}^{+}(x)=0$

$Q_{1}=Q_{1}^{-}$

$z_{1}=z_{1}^{-}$

The wave propagates as a pure exponentially damped sine. The rate at which it travels corresponds to the imaginary part, and at which it is damped, to the real part of the respective root. When the cross-section geometry varies along the channel, then the incoming wave is partially reflected. It follows from (3.13a ) that in this case, the wavenumber differs from the corresponding root of the characteristic polynomial. The coefficients $c_{\{0,1,2\}}$ vary as well, and the wave propagates as

(3.15a ) $$\begin{eqnarray}\displaystyle \frac{\unicode[STIX]{x2202}Q_{1}^{-}}{\unicode[STIX]{x2202}x} & = & \displaystyle \left(r^{-}+\underbrace{\frac{-1}{r^{-}-r^{+}}\frac{\unicode[STIX]{x2202}r^{-}}{\unicode[STIX]{x2202}x}}_{T^{-}}\right)Q_{1}^{-}+\underbrace{\frac{-1}{r^{-}-r^{+}}\frac{\unicode[STIX]{x2202}r^{+}}{\unicode[STIX]{x2202}x}}_{R^{+}}Q_{1}^{+},\end{eqnarray}$$

(3.15b ) $$\begin{eqnarray}\displaystyle \frac{\unicode[STIX]{x2202}Q_{1}^{+}}{\unicode[STIX]{x2202}x} & = & \displaystyle \underbrace{\frac{+1}{r^{-}-r^{+}}\frac{\unicode[STIX]{x2202}r^{-}}{\unicode[STIX]{x2202}x}}_{R^{-}}Q_{1}^{-}+\left(r^{+}+\underbrace{\frac{+1}{r^{-}-r^{+}}\frac{\unicode[STIX]{x2202}r^{+}}{\unicode[STIX]{x2202}x}}_{T^{+}}\right)Q_{1}^{+},\end{eqnarray}$$

$T^{-}$

$T^{+}$

$R^{-}$

$R^{+}$

When the cross-section geometry varies smoothly at a low rate, then the amplitude of the reflected wave is negligible so that $Q_{1}\approx Q_{1}^{-}$ and $k_{1}\approx -\text{i}(\!r^{-}+(1/(r^{-}-r^{+}))(\unicode[STIX]{x2202}r^{-}/\unicode[STIX]{x2202}x)\!)$ . Thus, even when the reflected wave is small, the incoming wave can change considerably by transmission. For infinitesimally small waves, $r^{\pm }$ does not depend on $Q_{1}$ , and (3.15a ) gives direct insight into the propagation of the tidal wave along a river with known geometry.

For the sake of illustration, consider the case where the width remains constant along the channel so that $r^{-}=-r^{+}$ . When the cross-section geometry changes smoothly at a low rate, (3.15a ) and (3.15b ) simplify to

(3.16a ) $$\begin{eqnarray}\displaystyle \frac{1}{Q_{1}}\frac{\unicode[STIX]{x2202}Q_{1}}{\unicode[STIX]{x2202}x} & = & \displaystyle r^{-}-\frac{1}{2r^{-}}\frac{\unicode[STIX]{x2202}r^{-}}{\unicode[STIX]{x2202}x},\end{eqnarray}$$

(3.16b ) $$\begin{eqnarray}\displaystyle \frac{1}{z_{1}}\frac{\unicode[STIX]{x2202}z_{1}}{\unicode[STIX]{x2202}x} & = & \displaystyle r^{-}+\frac{1}{2r^{-}}\frac{\unicode[STIX]{x2202}r^{-}}{\unicode[STIX]{x2202}x},\end{eqnarray}$$

$\unicode[STIX]{x2202}r^{\pm }/\unicode[STIX]{x2202}x$

Equations (3.15b ) and (3.15a ) show that a convergence of the cross-section has the opposite effect on the upstream travelling wave ( $Q^{-},z^{-}$ ) and reflected waves ( $Q^{+},z^{+}$ ), as the sign in front of $\unicode[STIX]{x2202}r^{\pm }/\unicode[STIX]{x2202}x$ is equal. In contrast, friction damps the incoming and outgoing waves at the same rate, as the sign in front $r^{\pm }$ is the opposite. For the same reason, equations (3.16a ) and (3.16b ) show that convergence likewise has the opposite effect on the discharge $Q_{1}$ and surface elevation $z_{1}$ , while they are also damped at the same rate.

The tide can be approximated by integrating the approximate wavenumber (3.10a ) along the river (3.16a ). The wavenumber thus corresponds to the sum of the negative root of the characteristic polynomial and the coefficient of transmission $k_{1Q}\approx -\text{i}(r^{-}+T^{-})$ . The admittance is consequently the product of two factors. The first accounts for the effect of gravity and friction, and the second for the effect of width and depth convergence. Gravity and friction always act on the tide, even if the cross-section geometry does not vary along the river (3.2b ). They form the zero-order terms that enter $r^{\pm }$ . There is only convergence when the cross-section geometry varies along the river, which is represented by the partial derivatives that enter $T^{\pm }$ and $R^{\pm }$ .

3.6 The effect of gravity and friction

When only gravity acts on the wave, i.e. when both friction and $\unicode[STIX]{x2202}w/\unicode[STIX]{x2202}x$ are zero, the wavenumber (3.16a ) is identical to $\text{i}$ -times the negative root of the characteristic polynomial (3.13b ) and simplifies to

(3.17) $$\begin{eqnarray}k_{1,0}=\frac{\unicode[STIX]{x1D714}_{1}}{\sqrt{gh}},\quad a=0.\end{eqnarray}$$

When both gravity and friction act, i.e. $\unicode[STIX]{x2202}w/\unicode[STIX]{x2202}x$ is zero, the wavenumber is (cf. Godin (Reference Godin1985))

(3.18) $$\begin{eqnarray}k_{1,a}=k_{1,0}\sqrt{1+2\text{i}a},\end{eqnarray}$$

where $a$ is measures the strength of friction (cf. (3.8b ))

(3.19) $$\begin{eqnarray}a=\frac{1}{\unicode[STIX]{x03C0}}\frac{c_{d}}{\unicode[STIX]{x1D714}wh_{0}^{2}}\left(f_{2}Q_{0}+\frac{1}{2}f_{1}Q_{hr}\right).\end{eqnarray}$$

The wave travels in the direction into which it is driven by gravity, and friction acts against it. The surface amplitude and discharge are thus affected in the same manner, $k_{1,0Q}=k_{1,0z}$ . The friction scale $a$ varies along the river and with the strength of the flow. A Puiseux series expansion with respect to the parameter $a$ reveals the effect of friction for low river flow

(3.20a ) $$\begin{eqnarray}k_{1,a}=k_{1,0}(1+\text{i}a),\quad a\rightarrow 0\end{eqnarray}$$

and high river flow, respectively,

(3.20b ) $$\begin{eqnarray}k_{1,a}=k_{1,0}(1+\text{i})\sqrt{a},\quad a\rightarrow \infty .\end{eqnarray}$$

The imaginary and real parts of (3.20a ) and (3.20b ) determine the rates of damping and phase change. Substitution of (3.17) and (3.19) into (3.20a ) reveals that when the flow is low, the friction damps the tide proportionally to the discharge and $h_{0}^{-5/2}$ , but it does not influence the phase. When the flow is strong (3.20b ), the friction determines both the rates of damping and phase change. Both rates approach the same value that is proportional to the square root of the discharge and $h_{0}^{-3/2}$ .

3.6.1 Low river flow

During periods when the river discharge is much smaller than the tidal discharge, the friction coefficients in (3.8a ) and (3.8b ) attain the values $f_{1}=3/8$ and $f_{2}=0$ . This is identical to the approximation by Lorentz (Terra, van de Berg & Maas Reference Terra, van de Berg and Maas2005). In this case,

(3.21) $$\begin{eqnarray}a=\frac{c_{d}}{\unicode[STIX]{x1D714}wh_{0}^{2}}\frac{8}{3\unicode[STIX]{x03C0}}Q_{1}\left(1-4\frac{Q_{0}^{2}}{|Q_{1}|^{2}}\right),\quad |Q_{1}|\gg |Q_{0}|.\end{eqnarray}$$

Damping is thus asymptotically insensitive to river discharge when the tide is strong. River discharge does not add noticeably to the damping as long as $Q_{0}<(4/3\unicode[STIX]{x03C0})|Q_{1}|$ (figure 5 b). The water depth increases linearly with the river discharge when the river discharge is low, as from (3.3) and (3.8a ), it follows that $h_{0}\approx h_{0}|_{Q_{0}=0}+x(8/3\unicode[STIX]{x03C0})(c_{d}w/gA_{0}^{3})Q_{hr}Q_{0}$ (figure 5 a). Damping can thus even decrease with the river discharge before a threshold is reached, as the linear increase in water depth can reduce the friction by a larger amount than it is increased by the square of the river discharge, as long as the river flow does not considerably increase the roughness.

Figure 5. Magnitudes of the frequency components of the friction term ((3.8a )–(3.8c ), bold), depending on the relative strength of river and tidal flow ( $|Q_{0}|/|Q_{1}|$ ), as well as their low flow asymptotes ( $|Q_{0}|/|Q_{1}|\rightarrow 0$ , dash-dotted) and high flow asymptotes ( $|Q_{0}|/|Q_{1}|\rightarrow \infty$ , dashed). Note that the discharge scale is identical to the velocity scale as $|u_{0}|/|u_{1}|=|Q_{0}|/|Q_{1}|$ .

3.6.2 Strong river flow

When the river discharge is so strong that the flow does not reverse over the tidal cycle, then the friction coefficients (3.8a ) and (3.8b ) obtain the values $f_{1}=0$ and $f_{2}=\unicode[STIX]{x03C0}$ .

(3.22) $$\begin{eqnarray}\displaystyle a=\frac{c_{d}}{\unicode[STIX]{x1D714}wh_{0}^{2}}2Q_{0},\quad |Q_{0}|>|Q_{1}|. & & \displaystyle\end{eqnarray}$$

The tide only contributes to the damping rate by modulating the water depth when the river discharge is large, which does not affect the first-order approximation of the damping.

3.7 The effect of width and depth convergence

Width and depth convergence modify the wavenumber by the term $\unicode[STIX]{x0394}k_{1}$ . When the cross-section geometry changes smoothly along the river, the term is,

(3.23a ) $$\begin{eqnarray}\displaystyle \unicode[STIX]{x0394}k_{1} & = & \displaystyle \frac{1}{4(\text{i}-a)}\left((1+3\text{i}a)\frac{1}{h_{0}}\frac{\unicode[STIX]{x2202}h_{0}}{\unicode[STIX]{x2202}x}+(2+3\text{i}a)\frac{1}{w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}\right),\end{eqnarray}$$

(3.23b ) $$\begin{eqnarray}\displaystyle k_{1z} & = & \displaystyle k_{1,a}+\unicode[STIX]{x0394}k_{1},\end{eqnarray}$$

(3.23c ) $$\begin{eqnarray}\displaystyle k_{1Q} & = & \displaystyle k_{1,a}-\unicode[STIX]{x0394}k_{1}.\end{eqnarray}$$

In contrast to damping, width and depth convergence have the opposite effect on the surface amplitude and the discharge. This corollary of Green’s law (Green Reference Green1838; Jay Reference Jay1991) thus also holds in the presence of friction. The sign of the convergence term also depends on the direction in which the wave travels and in which the cross-section changes. Convergence increases the tidal amplitude when the cross-section becomes narrower and shallower.

In the limit of low friction, the change in wavenumber with the rate of convergence is

(3.24) $$\begin{eqnarray}\unicode[STIX]{x0394}k_{1}=-\text{i}\frac{1}{4}\frac{1}{h_{0}}\frac{\unicode[STIX]{x2202}h_{0}}{\unicode[STIX]{x2202}x}-\text{i}\frac{1}{2}\frac{1}{w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x},\quad a\rightarrow 0.\end{eqnarray}$$

A relative change in width thus has a larger effect than a relative change in depth of the same magnitude. This is known as Green’s Law (Green Reference Green1838; Jay Reference Jay1991). When friction is strong, the rate of convergence is

(3.25) $$\begin{eqnarray}\unicode[STIX]{x0394}k_{1}=-\text{i}\frac{3}{4}\frac{1}{h_{0}}\frac{\unicode[STIX]{x2202}h_{0}}{\unicode[STIX]{x2202}x}-\text{i}\frac{3}{4}\frac{1}{w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x},\quad a\rightarrow \infty .\end{eqnarray}$$

Strong friction enhances the effect of convergence, and in contrast to low friction, the relative changes of width and depth have the same effect.

However, the river discharge influences the effect of convergence not only indirectly by increasing friction but also directly, as depth convergence decreases with increasing discharge as $\unicode[STIX]{x2202}h_{0}/\unicode[STIX]{x2202}x=\unicode[STIX]{x2202}z_{0}/\unicode[STIX]{x2202}x-\unicode[STIX]{x2202}z_{b}/\unicode[STIX]{x2202}x$ (figure 6 a). The effects of width and depth convergence on the admittance of the tide increase monotonically with the friction and thus with the river and tidal discharge (figure 6 b). However, the depth convergence itself decreases with the river discharge as well, as the upstream water level rises, so that the effect on the admittance has a maximum for intermediate river discharges. Both width and depth convergence primarily affect the amplitude. The rate of phase change is only affected when friction is intermediate.

Figure 6. (a) Relative depth convergence (black) and friction scale (red) at the river mouth depending on the river discharge for infinitesimal waves; (b) effect of depth convergence on the damping rate (black) and rate of phase change (red) depending on the river discharge (solid) or tidal discharge (dashed); damping rate depending on friction and depth convergence (c) as well as width convergence (d) as approximated by (3.23a ), red line shows critical convergence ( $\text{Im}(k)=0$ ), blue and black are asymptotes for high and low friction, respectively.

Figure 7. (a) Tidally averaged water level along the river, where the backwater drawdown lengths are indicated by dots; (b) backwater and drawdown length $L$ for various states of river flow, defined as the distance to point where water depth deviates no more than 5 % $(L_{95})$ and approximated by simplified relations.

4.1 Tidally averaged water level

The tidally averaged water level changes along a river depending on the river discharge. When the river discharge is low, it forms a backwater profile, and depth decreases in the upstream direction (figure 8 black and red). Far upstream, the water surface slope asymptotically approaches the bed slope, so that the tidally averaged water depth remains constant along the river.

There are no analytic solutions to the backwater profile (3.3), with the exception of those of the Bresse type (Vatankhah & Easa Reference Vatankhah and Easa2011), which cannot be integrated into a general solution of the tide because they swap the dependent and independent variables. For the analysis of the river tide, we thus linearize the backwater equation (3.3) and (3.8a ) at the river mouth,

(4.1) $$\begin{eqnarray}h_{0}(x)=h_{0}(0)+\left(\left.\frac{\unicode[STIX]{x2202}z_{0}}{\unicode[STIX]{x2202}x}\right|_{0}-\frac{\unicode[STIX]{x2202}z_{b}}{\unicode[STIX]{x2202}x}\right)x+O\left(\left.\frac{\unicode[STIX]{x2202}^{2}z_{0}}{\unicode[STIX]{x2202}x^{2}}\right|_{0}x^{2}\right)=\frac{c_{d}}{g}\frac{Q_{0}|Q_{0}|}{w^{2}h_{0}^{3}},\end{eqnarray}$$

for the reach between the river mouth and the point where the flow becomes approximately uniform $(x=|h_{0}(0)-h_{u}|((\unicode[STIX]{x2202}h_{0}/\unicode[STIX]{x2202}x)|_{0})^{-1})$ ; $h_{u}$ is the depth in the reach of uniform flow far upstream.

When the river is at normal flow $(Q=Q_{n})$ , then the water surface slope is identical to the bed slope, and the tidally averaged water depth does not change along the river (blue in figure 7 a). When the river discharge is below normal flow, then the river forms a backwater curve (black, red and green in figure 7 a). When the river discharge is above normal flow, then the water surface forms a drawdown curve so that the depth increases into the upstream direction (orange in figure 7 a). For extremely low river discharge, the point where the flow becomes uniform approaches the point where the bed reaches sea level $x=L_{0}=h_{0}(\unicode[STIX]{x2202}z_{b}/\unicode[STIX]{x2202}x)^{-1}$ . As long as the river is in a state of backwater, this point is located the closer to the river mouth, the higher the discharge is (figure 7 b). As the wavenumber depends on the water depth, the backwater profile strongly influences the admittance of the tide.

4.2 Admittance along the river

In the case of normal flow, the depth does not converge. The tidal amplitude decreases exponentially along the channel so that the tidal amplitude drops most rapidly near the sea and less rapidly farther upstream (figure 8 b). The wave propagates with constant celerity (figure 9 a). At normal flow, the frictional damping is strong so that the damping rate is proportional to the square root of the river discharge, as given by (3.20a ). This is a well-known relation for the propagation of a tidal wave along a river with constant depth (LeBlond Reference LeBlond1978; Godin Reference Godin1985; Jay Reference Jay1991; Jay & Flinchem Reference Jay and Flinchem1997).

Figure 8. (a) Admittance of the tide for the case where the depth varies along the river (bold) and for the case where the depth remains constant along the river (dashed); inflexion points are indicated by dots, maxima by chevrons; (b) the logarithmic $y$ -axis reveals the inversion of the stage–amplitude relation in the asymptotic reach; scenarios as in figure 7(a).

Figure 9. (a) Phase difference between the tidal surface elevation along the river with respect to that at the mouth; (b) phase lag between the tidal velocity and surface level variation along the river. Without backwater effects, the phase difference is zero during low flow, as expected for a progressive wave (black dashed). During high flow, the phase difference is apparently half-way between that of a progressive and a standing wave, although there is no reflection. With the backwater effect, the phase of the tide always appears to be between that of a progressive and a standing wave (solid lines); scenarios as in figure 8.

However, when the river is not at normal flow, then the depth changes along the river, and the tide propagates differently compared to the case where the depth remains constant along the river. At low flow (black and red lines in figure 8 a), the tidally averaged water surface forms a backwater profile. The tidal amplitude decreases only slowly near the river mouth, and more rapidly farther upstream (black and red in figure 8 b). The transition is marked by an inflexion point at which $\unicode[STIX]{x2202}^{2}|z_{1}|/\unicode[STIX]{x2202}x^{2}=0$ . As (3.11b ) implies $(1/|z_{1}|)(\unicode[STIX]{x2202}|z_{1}|/\unicode[STIX]{x2202}x)=\text{Im}(k_{1z})$ , this implies

(4.2) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}\text{Im}(k_{1z})}{\unicode[STIX]{x2202}x}-\text{Im}(k_{1z})^{2}=0.\end{eqnarray}$$

The inflexion point is located closer to the sea for higher river discharges, similar to the transition point where the asymptotic flow is reached. For discharges higher than a particular value, the inflexion point is not observed anymore. Only when the river discharge is above this threshold is the tide admitted similarly to the conventional case where the depth is assumed to be constant along the river.

For very low discharges (black), the tidal range even increases in the landward direction to a point where a maximum is reached. The tidal admittance is hypersynchronous. At this point, the imaginary part of the wavenumber is zero,

(4.3) $$\begin{eqnarray}\text{Im}(k_{1z})=0.\end{eqnarray}$$

The maximum is located closer to the sea than the inflexion point. This point is related to the concept of critical damping that is used to characterize tidally dominated estuaries (Jay Reference Jay1991). Along such estuaries, width convergence can cancel frictional damping. If the tidally averaged water depth remains constant along an estuary, it can be critically damped along its entire length so that the tidal amplitude neither increases nor decreases, in which case, an estuary is considered to be ‘ideal’. A fundamental difference between tide-dominated estuaries and tidal rivers with a sloping bed is thus that tidal rivers cannot represent ideal estuaries, as the damping can be critical only at one point, and not along the entire tidally influenced reach. The location of this point furthermore depends on the river discharge. Similar to the location of the inflexion point, the maximum is located closer to the sea for higher river discharges and vanishes when the discharge exceeds a certain threshold.

4.3 Damping and convergence rates

The rate of amplitude change along the river is the combined effect of frictional damping and convergence of the width and depth. For large depths, the damping rate is proportional to $h_{0}^{-5/2}$ (3.20a ), while depth convergence is proportional to $h_{0}^{-1}$ (3.23a ). At low flow, the river forms a backwater profile, and the depth increases towards the sea. The effects of both the frictional damping and convergence thus decrease towards the sea, but frictional damping decreases more rapidly so that close to the sea, the effect of convergence is larger than that of damping, and the amplitude of waves travelling upstream increases along the downstream part of the tidal river. In this reach, the imaginary part of $k_{1z}$ is thus positive (figure 10). Conversely, the water depth decreases into the upstream direction towards the point where the bed reaches sea level when the river discharge is low. For small depths, damping is proportional to $h_{0}^{-3/2}$ , while convergence is still proportional to $h_{0}^{-1}$ . The effects of both damping and convergence thus increase in the upstream direction, but damping increases more rapidly so that eventually frictional damping dominates, and the tidal amplitude is reduced in the upstream part of the tidal river. In this reach, the imaginary part of $k_{1z}$ is thus negative (figure 10). The depth only increases into the upstream direction when the river discharge is very large so that the water surface forms a drawdown curve. In that case, damping is very strong, and convergence has the opposite effect compared to the situation for low flow, so the tidal amplitude rapidly decreases.

Figure 10. (a) Imaginary (dashed) and real (solid) parts of the wavenumber $z_{1}$ , corresponding to the rate of damping and phase change, respectively; critical damping occurs at $\text{Im}(k_{1z})=0$ . (b) Relative strength of convergence compared to that of gravity and friction ( $\text{Im}(\unicode[STIX]{x0394}k_{1}/k_{1,0})$ , solid) and its approximation by (3.23a ) (dashed), critical damping occurs at $-1$ ; scenarios correspond to those in figure 8.

4.4 Asymptotic admittance

The influence of the river discharge on the propagation of the tide in the upstream reach where the flow becomes uniform has been studied by Godin (Reference Godin1985). However, this study does not pay attention to the influence of the river flow on the tidally averaged water depth. When the slope of the river bed is very small, the tidally averaged water depth does not strongly change with the river discharge. In such a case, the damping decreases with the river discharge so that the tide intrudes far upstream during low flow (figure 8 b). However, this is not the case when the slope of the river bed is not small. The admittance is always marginal beyond the point where the bed reaches sea level (figure 8 b). In the asymptotic reach (inset in figure 8), the tidal amplitude even slightly increases with the river discharge. Far upstream, the river flow becomes asymptotically uniform. The ordinary differential equation for the water surface elevation (3.3) simplifies to an algebraic equation referred to as the Chézy formula. The algebraic relation between discharge and surface elevation allows one to eliminate either the river discharge or the water depth from the expression for the wavenumber. Far upstream, damping is strong so that the wavenumber is approximated by (3.20b ),

(4.4a ) $$\begin{eqnarray}\displaystyle k_{1z} & = & \displaystyle -(1-\text{i})\sqrt{\frac{c_{d}\unicode[STIX]{x1D714}Q_{0}}{gwh_{0}^{3}}},\end{eqnarray}$$

(4.4b ) $$\begin{eqnarray}\displaystyle & = & \displaystyle -(1-\text{i})\sqrt{\frac{\unicode[STIX]{x1D714}w}{Q_{0}}\frac{\unicode[STIX]{x2202}z_{b}}{\unicode[STIX]{x2202}x}},\end{eqnarray}$$

(4.4c ) $$\begin{eqnarray}\displaystyle & = & \displaystyle -(1-\text{i})\sqrt[4]{\frac{c_{d}\unicode[STIX]{x1D714}^{2}}{gh_{0}^{3}}\frac{\unicode[STIX]{x2202}z_{b}}{\unicode[STIX]{x2202}x}},\quad x>\left(\frac{1}{h_{0}}\frac{\unicode[STIX]{x2202}z_{b}}{\unicode[STIX]{x2202}x}\right)^{-1}.\end{eqnarray}$$

Figure 11. Tidal admittance for three stations located at a tidal river subject to backwater effects (black) and as predicted by a conventional model that neglects the backwater (red); dashed lines show the admittance of a wave with infinitesimal amplitude, solid lines show the admittance of a wave with finite amplitude; without backwater effects, the admittance decays exponentially; with backwater effects, the admittance decreases hyperbolically along the river.

As intuition suggests, the damping becomes the stronger, the larger the slope, due to the higher flow velocity and shallower depth. The expression also shows that damping is indeed asymptotically inversely proportional to the square root of the river discharge so that in the asymptotic reach, the amplitude increases with the river discharge, which is the opposite of what is predicted by conventional models, where the depth remains constant along the river. The expression bears yet another surprise: neither the rate of frictional damping nor that of the phase change asymptotically depend on the drag coefficient for a particular river discharge. The tidal damping along the upstream reach has thus to be calibrated by adjusting other parameters, such as bed slope or channel geometry, rather than the drag coefficient, which is a common parameter for the calibration of hydrodynamic models. Similarly, the rate of damping and phase change are asymptotically independent of the water depth for a particular river discharge.

We point out that upstream propagation of tidal waves is not primarily inhibited by the pull of gravity down the slope, but rather by high friction caused by either the low water depth during low discharges, or by high flow velocities during high discharges. For small discharges, the region where the damping is low extends farther upstream; however, in the asymptotic reach beyond the point where the bed reaches sea level, the tidal amplitude always decreases along the channel. The transition between the reach where friction is low and the asymptotic reach where friction is high is thus more rapid during periods of low river flow than during periods of high river flow.

4.5 At-a-station admittance

The tidal amplitude that is observed at a gauging station depends on the river discharge, in analogy to the admittance of the tide along the river (§ 4.2). Figure 11 shows the admittance of the tide at three stations along a tidal river with a sloping bed in comparison to the admittance that is predicted by a conventional model, which neglects the backwater dynamics and assumes a constant depth along the river.

At a station near the sea, the tidal amplitude is always higher than that predicted by the conventional model (figure 11 a) except at normal flow, when the amplitudes are equal. At a station in the transition reach (figure 11 b), the admittance is higher during low flow and high flow but lower during intermediate river flow. Beyond the point where the river bed reaches sea level (figure 11 c), the admittance is much smaller during low flow and slightly larger during high flow than predicted without bed slope. In general, when the backwater dynamics is neglected, the admittance is predicted to decrease exponentially with the river discharge. The backwater introduces a hyperbolic factor to the admittance due to the gradual change in depth. The tidal amplitude drops less rapidly than predicted by the exponential approximation when the river discharge is low, and more rapidly when the river discharge is high. The hyperbolicity of the tidal admittance can be shown on hand of the propagation of the upstream travelling wave $\unicode[STIX]{x2202}z_{1}^{-}/\unicode[STIX]{x2202}x=-\text{i}(k_{1,a}-\unicode[STIX]{x0394}k_{1})z_{1}^{-}$ (3.23c ), when the reflected wave is neglected ( $z_{1}\approx z_{1}^{-}$ ). The admittance is found by integration $z_{1}(x)=z_{1}(0)\exp (\int _{0}^{x}-\text{i}(k_{1,a}-\unicode[STIX]{x0394}k)\,\text{d}x^{\prime })$ . The linearization of the backwater equation leads to $\unicode[STIX]{x0394}k_{1}=-\text{i}(p/(1+(1/h_{0}(0))(\unicode[STIX]{x2202}h_{0}/\unicode[STIX]{x2202}x)x))(\unicode[STIX]{x2202}h_{0}/\unicode[STIX]{x2202}x)$ , and an admittance of the form

(4.5a ) $$\begin{eqnarray}\displaystyle z_{1}(x) & {\approx} & \displaystyle \left(\frac{h_{0}(0)}{h_{0}(x)}\right)^{p}\exp \left(\int _{0}^{x}-\text{i}k_{1,a}\,\text{d}x^{\prime }\right)z_{1}(0)\end{eqnarray}$$

(4.5b ) $$\begin{eqnarray}\displaystyle & = & \displaystyle \frac{1}{\displaystyle \left(1+\frac{1}{h_{0}(0)}\frac{\unicode[STIX]{x2202}h}{\unicode[STIX]{x2202}x}x\right)^{p}}\exp \left(\int _{0}^{x}-\text{i}k_{1,a}\,\text{d}x^{\prime }\right)z_{1}(0),\end{eqnarray}$$

$p$

$1/4$

$3/4$

4.6 Reflected wave

The approximation (3.16a ) considers only the effect of transmission on the tide but not reflection. However, waves are partially reflected whenever they propagate along a channel where the cross-section changes over a distance that is shorter than its wavelength (Lighthill Reference Lighthill2001). The tide is thus gradually reflected, as the depth diminishes along the backwater affected reach of a tidal river.

When the tidal amplitude is small, the wave reflected back downstream can be approximated by the following double integral,

(4.6) $$\begin{eqnarray}Q_{1}^{+}(x)=\int _{x}^{\infty }R^{-}Q_{1}^{-}(x^{\prime })\exp \left(\int _{x^{\prime }}^{x}r^{+}+T^{+}\,\text{d}x^{\prime \prime }\right)\,\text{d}x.\end{eqnarray}$$

The wave that is reflected back downstream is again partially reflected so that the upstream and downstream travelling waves are linearly dependent on each other (Witting Reference Witting1981) ((3.15a ) and (3.15b )). The step (4.6) can be iterated to decompose the tide into a series that can be interpreted as recursively reflected waves (Wilmer III & Costa Reference Wilmer and Costa2008). The reflected waves decrease rapidly in magnitude, due to the damping and gradual change in depth ( $|R^{\pm }|\ll 1$ ).

A decomposition of the tidal wave into the incoming and the reflected parts shows that the magnitude of the reflected wave is indeed not entirely negligible (figure 12). The reflected wave increases the amplitude of the total wave. As the reflected wave leaves the river, the amplitude and phase of the tide at the river mouth slightly differ from that of the incoming wave (figure 12, inset). The tide at the river mouth thus depends on the river geometry and river discharge. This is important for the prescription of boundary conditions. It seems natural to specify the surface level $z_{s}(t,0)$ at the river mouth, as the astronomical tide can reliably be predicted (Egbert & Erofeeva Reference Egbert and Erofeeva2002). This approach is typically applied in one-dimensional models of tidal rivers that do not include the adjacent sea (Bolla Pittaluga et al. Reference Bolla Pittaluga, Tambroni, Canestrelli, Slingerland, Lanzoni and Seminara2015). In § 3.6, we apply this to approximate the propagation of the wave. However, this implicitly reduces the amplitude of the tide entering from the sea by the amplitude of the reflected wave and thus violates causality, as it implies that the astronomical tide depended on the estuarine geometry and river discharge. For the numerical computation of the tide in § 4, we therefore only prescribe a value for the incoming wave $\hat{z}_{1}^{-}$ at the seaward boundary, and let the seaward travelling wave $\hat{z}_{1}^{+}$ pass without reflection. As only one value is specified per frequency component at each boundary, we express it as a linear combination of the incoming and seaward travelling wave, according to (3.14b ). The mean water level at the mouth also varies over a Metonic cycle (Woodworth Reference Woodworth2017) and is due to the river discharge, but these variations are negligible.

Conversely, a trivial boundary condition of zero tidal amplitude or discharge causes reflection at the upstream boundary. Reflection at the upstream boundary can be reduced by placing the upstream boundary far from the sea (Godin & Martínez Reference Godin and Martínez1994; Cai et al. Reference Cai, Savenije and Toffolon2014). We avoid reflection by expressing the upstream boundary condition as a linear combination of the upstream and downstream travelling waves and set only the amplitude of the downstream travelling wave to zero as $Q_{1}^{+}=0$ .

Figure 12. (a) Admittance of the combined incoming and reflected waves ( $|z_{1}(x)|/|z_{1}^{-}(0)|$ , solid) and only of the upstream travelling wave ( $|z_{1}^{-}(x)|/|z_{1}^{-}(0)|$ , dash-dotted) along the river. (b) amplitude (solid) and phase (dashed) of the combined wave with respect to the incoming wave at the river mouth, depending on the river discharge. Scenarios as in figure 8.

4.7 Tidal propagation where the bed reaches sea level

When the river discharge is very low, the water is shallow near the point where the bed reaches sea level so that the relative change in depth along the river is large. The wave is therefore partially reflected according to (3.15a ). Locally, the reflected wave adds noticeably to the water surface oscillation, but as it travels back, it rapidly decreases in amplitude because the friction is high and the divergence of the cross-sectional area reduces the amplitude of the reflected wave. Near the head of the tidal river, the amplitude of the main tidal species can exceed that of the water depth. The overtide rises low water by changing the shape of the tide so that the river bed does not necessarily fall dry during periods of low river flow.

Without river discharge, the bed dries up beyond the point where the bed reaches sea level so that the tide cannot propagate much farther upstream. The model (3.9a ) forces the tidal amplitude to zero at the point where the bed reaches sea level. The phase of the tide changes at an ever-increasing rate towards the point where the depth reaches zero so that the phase is undefined at this point. This is the case for all waves of zero amplitude and is not an artefact of the model. In practice, a short reach upstream of the head of the tidal river can periodically flood and dry over the tidal cycle. The flow furthermore concentrates near the thalweg, and shallow parts of the cross-section are periodically flooded. A model that has a bed that is not flat across the river can better predict how the tide propagates near the head. However, the solution is sensitive to unevenness of the river bed, to seasonal variation of the sea level as well as to the residual river flow. Wetting and drying of the entire cross-section is also beyond the scope of the shallow-water equations (3.2a ) and (3.2b ), as those require the depth to be non-zero at all time. In any case, the tide does not propagate far beyond the point where the bed reaches sea level during periods without river flow, as it rapidly decays due to the shallow depth as well as storage on the higher parts of the bed that are periodically flooded. The case of the wave vanishing when running up the slope thus contrasts with the case where the wave runs against a wall, where the surface amplitude obtains a maximum, as well as the case with a wave propagation along a channel with a horizontal bed, where the wave propagates incessantly, without ever entirely decaying to zero. A minimum base flow ensures that the bed of the Kapuas never falls entirely dry. However, bars and dunes emerge during low flow beyond the point where the bed reaches sea level, which is slightly upstream of the city of Sanggau.

Figure 13. Magnitudes of tidal discharge (a) and tidal velocity (b), along a tidal river with a sloping bed (solid), and constant depth (dashed); river discharge indicated by dash-dotted line in (a); scenarios as in figure 8.

4.8 Tidal discharge

The tidal discharge decreases with increasing river discharge (figure 13 a). It diminishes gradually along the river so that no extrema or inflexion points occur.

As long as the river discharge is low, the tidally averaged water surface forms a backwater profile, and the tidal discharge at the river mouth decreases at a higher rate than that predicted without a sloping bed.

The most upstream point of flow reversal is located closer to the sea, compared to the conventional model. Similar to the tidal water surface amplitude, the tidal discharge amplitude is marginal beyond the point where the river bed reaches sea level. The higher rate at which the discharge is reduced is caused by depth convergence, which, recalling (3.23a ), has an effect on the discharge that is opposite to the effect on the tidal surface level amplitude. The lower initial discharge at the river mouth is also a direct consequence of depth convergence. Recalling the continuity (3.2a ), the tidal discharge of the main species is

(4.7) $$\begin{eqnarray}Q_{1}=\text{i}w\frac{\unicode[STIX]{x1D714}}{k_{1Q}}z_{1},\end{eqnarray}$$

which can be directly evaluated at the river mouth in combination with the approximate wavenumber (3.16a ). The tidal discharge at the river mouth is thus smaller because it is inversely proportional to the wavenumber $k_{1Q}$ , which itself is larger in magnitude, due to depth convergence.

4.9 Tidal velocity

The cross-section averaged flow velocity is directly proportional to the discharge $u=1/A_{0}Q$ and can thus also be explicitly evaluated at the river mouth. The tidal velocity amplitude decreases with the river discharge. However, the magnitude of the tidal velocity amplitude decreases along the river at a lower rate than the tidal discharge when the river is below normal flow, which follows from the chain rule,

(4.8) $$\begin{eqnarray}k_{1u}=-\text{i}\frac{1}{u_{1}}\frac{\unicode[STIX]{x2202}u_{1}}{\unicode[STIX]{x2202}x}=\text{i}\frac{1}{h_{0}}\frac{\unicode[STIX]{x2202}h_{0}}{\unicode[STIX]{x2202}x}+\text{i}\frac{1}{w}\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}x}+k_{1Q}.\end{eqnarray}$$

Width and depth convergence increase the magnitude of the velocity compared to the discharge, similar to the surface elevation amplitude, but the phase is similar to that of the discharge. An inflexion point marks the transition to the upstream reach of the tidal river, where the tidal velocity amplitude is rapidly reduced. The inflexion point of the velocity is located the closer to the river mouth, the higher the river discharge (figure 13), similarly to the inflection point of the tidal surface elevation amplitude. Overall, the magnitude of the tidal velocity amplitude decreases at a slightly higher rate than the magnitude of the tidal surface elevation, and the inflexion point of the velocity is located closer to the sea.

When the river flow is very low, the tidal velocity amplitudes are higher than that predicted by a conventional model along the downstream reach of a tidal river. The value is furthermore lower for intermediate discharges compared to that predicted without a sloping bed. The tidal excursion length near the sea is similar but decreases more rapidly towards upstream than that predicted without bed slope (not shown). As the salinity intrusion depends on the tidal excursion (Savenije Reference Savenije2012), this implies that a sloping bed shortens the distance that salinity intrudes up river, as long as the slope does not strongly change the mixing.

4.10 Phase lag

The damping of the tide is related to the phase lag $\unicode[STIX]{x0394}\unicode[STIX]{x1D711}_{1}=\arg (u_{1})-\arg (z_{1})$ between the oscillation of the water surface and velocity. When the tidally averaged water depth does not change along the river, then the damping rate and phase lag are directly proportional as long as the damping is low, which follows from (3.20a )

(4.9) $$\begin{eqnarray}\unicode[STIX]{x0394}\unicode[STIX]{x1D711}_{1}=\arctan (a)\approx a,\quad a<1,\end{eqnarray}$$

which coincides with the definition of Savenije (Reference Savenije2012) for the case of non-converging estuaries. For weak damping ( $a\rightarrow 0$ ), the velocity is in phase with the surface elevation, as expected for a progressive wave. For strong damping $(a\gg 1)$ , the phase lag is close to $\unicode[STIX]{x03C0}/4$ , according to (3.20b ). The phase lag is apparently in between that of a progressive and a standing wave, although the wave is not reflected.

Along rivers with a sloping bed, the tidal amplitudes of the surface level and discharge change at different rates so that damping is no longer proportional to the phase lag (figure 9 b). This is the case because friction is always strong in the upstream part of the tidal river and because a part of the wave is reflected gradually along the backwater affected reach. As a consequence, the phase lag is close to $\unicode[STIX]{x03C0}/4$ independent of the river flow and can even exceed $\unicode[STIX]{x03C0}/4$ in the upstream part of the tidal river. The phase lag thus does not linearly increase with the river discharge along the backwater affected reach. In the asymptotic flow reach, the phase lag even decreases with the river discharge, opposite to what occurs in a river with a horizontal bed.

Figure 14. (a) At a station, water level variation over a tidal cycle for $z_{1}(0)=0.1h_{0}(0)$ ; (b) along river amplitude of the even overtide; (c) along river subtidal water level set-up; solid lines show the tide in a backwater affected river, and dashed lines show the prediction of the tide when the backwater is neglected; scenarios as in figure 8.

4.11 Overtide generation

Friction not only damps the tide but also generates higher frequency components, referred to as the overtide. The principal overtide advances high water and delays low water so that the tidal wave changes its shape while travelling upstream (figure 14 a). The river discharge forces the principal overtide with twice the frequency and an amplitude proportional to the square of the main tidal species, which is evident from the frequency decomposition of the friction (3.8a ).

When the tidally averaged width and water depth do not change along the river, then the even overtide is determined by the equation

(4.10) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}^{2}Q_{2}}{\unicode[STIX]{x2202}x^{2}}-2\text{i}r_{2}Q_{2}=-\text{i}\frac{r_{1}^{2}}{w}\frac{Q_{1}^{2}}{Q_{0}},\end{eqnarray}$$

which is a simplification of (3.4) and (3.8c ) for the limit that the river flow is much stronger than the tidal flow. The wavenumber $k_{1}$ is given by (3.18); $r_{1}$ is its imaginary part, and $k_{2}=\sqrt{2}k_{1}$ as well as $r_{2}=\sqrt{2}r_{1}$ . This gives the amplitude of the overtide as

(4.11) $$\begin{eqnarray}z_{2}=\frac{\unicode[STIX]{x1D714}_{1}w}{8r_{1}Q_{0}}z_{1}(0)^{2}(\exp (-2\text{i}k_{1}x)-\exp (-1\text{i}k_{2}x)),\end{eqnarray}$$

for the case when the amplitude of the overtide is zero at the river mouth. Note that even if the surface elevation amplitude of the overtide is zero at the river mouth, the amplitude of the velocity is not, which is due to continuity (4.7). Equation (4.11) is the sum of two decaying exponentials, which add to zero at the river mouth. Near the mouth the principal overtide increases by 0.6 times the rate of that at which the amplitude of the main tidal species decreases, and far upstream, it decreases by 1.4 times the rate of that of the main tidal species.

In a river with a sloping bed, the principal overtide increases near the sea, similar to the case when the tidally averaged water depth remains constant. However, in the upstream part of the tidal river, the tide is damped more rapidly when the bed is sloping. Similar to the main tidal species, the amplitude becomes marginal beyond the point where the bed level reaches sea level (figure 14 b).

4.12 Water level set-up

The tidal flow raises the tidally averaged water level by increasing friction (3.3). The water level, therefore, oscillates with the tidal range over the spring-neap cycle. Such subtidal harmonics have been described by LeBlond (Reference LeBlond1979), Kukulka & Jay (Reference Kukulka and Jay2003a ), Buschman et al. (Reference Buschman, Hoitink, van der Vegt and Hoekstra2009) and Sassi & Hoitink (Reference Sassi and Hoitink2013). The water level set-up $z_{0}^{\prime }$ caused by the tide with respect to the water depth $h_{0}^{\prime }$ that is only due to river flow can be approximated by a perturbation expansion of the backwater equation (3.3). For dominant river flow, this is done by substituting $h_{0}^{\prime }+z_{0}^{\prime }$ for $h_{0}$ , expanding the equations as a series in $z_{0}^{\prime }$ and keeping only the terms of largest magnitude,

(4.12) $$\begin{eqnarray}\frac{\unicode[STIX]{x2202}z_{0}^{\prime }}{\unicode[STIX]{x2202}x}=\frac{c_{d}}{gw^{2}h_{0}^{4}}\left(\frac{1}{2}|Q_{1}|^{2}h_{0}-3Q_{0}^{2}z_{0}^{\prime }\right).\end{eqnarray}$$

The first term generates and the second term damps the water level set-up. For a river at normal flow, where the tidally averaged water depth remains constant along the channel, the water level set-up is,

(4.13) $$\begin{eqnarray}\displaystyle & \displaystyle z_{0}^{\prime }=\frac{h_{0}}{6}\frac{r_{0}}{r_{0}-r_{1}}\frac{|Q_{1}(0)|^{2}}{|Q_{0}|^{2}}(\exp (-2r_{0}x)-\exp (-2r_{1}x)), & \displaystyle\end{eqnarray}$$

(4.14) $$\begin{eqnarray}\displaystyle & \displaystyle r_{0}=\frac{1}{2}\frac{c_{d}Q_{0}^{2}}{gw^{2}h_{0}^{4}}. & \displaystyle\end{eqnarray}$$

A similar relation has recently been published by Tolkova (Reference Tolkova2018). Here, $r_{1}=\text{Im}(k_{1,a})$ is the damping rate of the main tidal species as given by (3.18), and $r_{0}$ is the damping rate of the mean water level set-up; $Q_{1}(0)$ is the discharge at the river mouth as given by (4.7). The water level set-up is the difference of two decaying exponentials, which are zero at the river mouth. The water level set-up rises along the river twice as rapidly as the main tidal species is damped but thereafter, decays at a lower rate than the main tidal species, because $r_{0}\ll r_{1}$ (figure 14 c).

In rivers with a sloping bed, the tidally averaged water level set-up is always marginal beyond the point where the bed reaches sea level (figure 14 c). The water level set-up has thus a pronounced local maximum. This maximum is lower and closer to the river mouth for higher river discharge. The rapid reduction of the water level set-up in the upstream part of the tidal river is, as the reduction of the tidal amplitude, caused by the variation of the water depth with the river discharge.

Far upstream, where the flow becomes asymptotically uniform, the damping rate can be expressed as a function of the water level only, similar to that of the main tidal species (4.4c )

(4.15) $$\begin{eqnarray}r_{0}=-\frac{1}{2}\frac{1}{h_{0}}\frac{\unicode[STIX]{x2202}z_{b}}{\unicode[STIX]{x2202}x},\quad x\rightarrow \infty .\end{eqnarray}$$

This rate is higher than that of the main tidal species (4.4c ), and similar to it, decreases with the river discharge.

The amplitude of the principal overtide and the subtidal water level set-up change along the channel in a similar manner because they are both forced by the part of the friction term that is proportional to $Q_{1}^{2}$ ; however, the overtide damps out more rapidly. When the main tidal species is diurnal, as in the Kapuas, then the overtide is superimposed on the semi-diurnal tide that enters the river from the sea so that the semi-diurnal tide along the river strongly depends on the river flow.