Analytical and numerical analysis of tides and salinities in estuaries; part I: tidal wave propagation in convergent estuaries
Leo C. van Rijn
Ocean Dynamics Theoretical, Computational and Observational Oceanography
ISSN 1616-7341 Volume 61 Number 11
Ocean Dynamics (2011) 61:1719-1741 DOI 10.1007/s10236-011-0453-0
1 23 Your article is protected by copyright and all rights are held exclusively by Springer- Verlag. This e-offprint is for personal use only and shall not be self-archived in electronic repositories. If you wish to self-archive your work, please use the accepted author’s version for posting to your own website or your institution’s repository. You may further deposit the accepted author’s version on a funder’s repository at a funder’s request, provided it is not made publicly available until 12 months after publication.
1 23 Author's personal copy
Ocean Dynamics (2011) 61:1719–1741 DOI 10.1007/s10236-011-0453-0
Analytical and numerical analysis of tides and salinities in estuaries; part I: tidal wave propagation in convergent estuaries
Leo C. van Rijn
Received: 3 February 2011 /Accepted: 9 June 2011 /Published online: 10 July 2011 # Springer-Verlag 2011
Abstract Analytical solutions of the momentum and closed-end channels is important in the most landward 1/3 energy equations for tidal flow are studied. Analytical length of the total channel length. In strongly convergent solutions are well known for prismatic channels but are less channels with a single forward propagating tidal wave, well known for converging channels. As most estuaries there is a phase lead of the horizontal and vertical tide close have a planform with converging channels, the attention in to 90o, mimicking a standing wave system (apparent this paper is fully focused on converging tidal channels. It standing wave). will be shown that the tidal range along converging channels can be described by relatively simple expressions Keywords Tidal propagation . Salinity. Convergent solving the energy and momentum equations (new estuary. Prismatic channel . Analytical model . Numerical approaches). The semi-analytical solution of the energy model equation includes quadratic (nonlinear) bottom friction. The analytical solution of the continuity and momentum equations is only possible for linearized bottom friction. 1 Introduction The linearized analytical solution is presented for sinusoidal tidal waves with and without reflection in strongly The shape of most alluvial estuaries is similar all over the convergent (funnel type) channels. Using these approaches, world, see Dyer (1997), McDowell and O’Connor (1977), simple and powerful tools (spreadsheet models) for tidal Savenije (2005), and Prandle (2009). The width and the analysis of amplified and damped tidal wave propagation in area of the cross-section reduce in upstream (landward) converging estuaries have been developed. The analytical direction with a river outlet at the end of the estuary solutions are compared with the results of numerical resulting in a converging (funnel shape) channel system solutions and with measured data of the Western Scheldt (see Fig. 1). The bottom of the tide-dominated section Estuary in the Netherlands, the Hooghly Estuary in India generally is fairly horizontal. Often, there is a mouth bar at and the Delaware Estuary in the USA. The analytical the entrance of the estuary. Tidal flats or islands may be solutions show surprisingly good agreement with measured present along the estuary (deltas). tidal ranges in these large-scale tidal systems. Convergence The tidal range in estuaries is affected by four dominant is found to be dominant in long and deep-converging processes: (a) inertia related to acceleration an deceleration channels resulting in an amplified tidal range, whereas effects, (b) amplification (or shoaling) due to the decrease bottom friction is generally dominant in shallow converging of the width and depth (convergence) in landward direction channels resulting in a damped tidal range. Reflection in and (c) damping due to bottom friction, and (d) partial reflection at landward end of the estuary, whereas bathy- Responsible Editor: Roger Proctor metric convergence produces continuous reflection along the entire length of the channel (see point b). * L. C. van Rijn ( ) The classical solution of the linearized mass and Deltares and University of Utrecht, Utrecht, The Netherlands momentum balance equations for a prismatic channel of e-mail: [email protected] constant depth and width is well-known (Hunt 1964; Author's personal copy
1720 Ocean Dynamics (2011) 61:1719–1741
Fig. 1 Tidal estuary (planform and longitudinal section) Q-river
SEA
Bo +x -x B tidal excursion = 10 to 20 km
HW High water H MSL mean sea level
LW Low water h
Mouth bar
Dronkers 1964, 2005; Ippen 1966; Le Blond 1978; Verspuy more, they have not given the full solution including the 1985; Prandle 2009). This solution for a prismatic channel precise damping coefficient and wave speed expressions for represents an exponentially damped sinusoidal wave which both amplified and damped converging channels. Hunt dies out gradually in a channel with an open end or is (1964) briefly presents his solution for a converging reflected in a channel with a closed landward end. In a channel and focusses on an application for the Thames frictionless system with depth ho, both the incoming and Estuary in England. The analytical model is found to give 0.5 reflected wave have a phase speed of co=(gho) and have very reasonable results fitting the friction coefficient. Hunt equal amplitudes resulting in a standing wave with a virtual shows that strongly convergent channels can produce a wave speed equal to infinity due to superposition of the single forward-propagating tidal wave with a phase lead of incoming and reflected wave. Including (linear) friction, the the horizontal and vertical tide close to 90°, mimicking a wave speed of each wave is smaller than the classical value standing wave system (apparent standing wave). A basic co (damped co-oscillation). According to Le Blond (1978), feature of this system is that the wave speed is much larger 0.5 frictional forces dominate accelerations over most of the than the classical value co=(gho) , in line with observa- tidal cycle and hence a diffusive type of solution (neglect- tions. For example, the observed speed of the tidal wave in ing accelerations) provides accurate results for most of the the amplified Western Scheldt Estuary in The Netherlands tydal cycle in a prismatic channel. Using the classical is between 11 and 15 m/s, whereas the classical value is of approach, the tidal wave propagation in a funnel-type about 10 m/s. estuary can only be considered by schematizing the channel Parker (1984) has given a particular solution for a into a series of prismatic subsections, each with its own converging tidal channel with a closed-end focusing on the constant width and depth, following Dronkers (1964) and tidal characteristics (only sinusoidal (M2) tide) of the many others. Unfortunately, this approach eliminates to a Delaware Estuary (USA). He shows that the solution based large extent the effects of convergence in width and depth on linear friction and exponential decreasing width yields on the complex wave number and thus on the wave speed very reasonable results for the Delaware Estuary fitting the (Jay 1991). A better approach is to represent the planform friction coefficient. of the estuary by a geometric function. When an exponen- Harleman (1966) also included the effect of width tial function with a single-length scale parameter (Lb)is convergence by combining Greens’ law and the expressions used, the linearized equations can still be solved analyti- for a prismatic channel. Predictive expressions for the cally and are of an elegant simplicity. friction coefficient and wave speed were not given. Instead, The analytical solution for a funnel-type channel with he used measured tidal data to determine the friction exponential width and constant depth is, however, not very coefficient and wave length. well known. Hunt was one of the first to explore analytical Godin (1988) and Prandle and Rahman (1980) have solutions for converging channels using exponential and addressed a channel with both converging width and depth. power functies to represent the width variations. Both Le They show that that the analytical solution can be Floch (1961) and Hunt (1964) have given solutions for formulated in terms of Bessel functions for tidal elevations exponentially converging channels with constant depth. and tidal velocities in open and closed channels. However, However, their equations are not very transparent. Further- the complex Bessel functions involved obscure any Author's personal copy
Ocean Dynamics (2011) 61:1719–1741 1721 immediate physical interpretation. Therefore, their results horizontal and vertical tide. Only bulk parameters are were illustrated in diagrammatic form (contours of ampli- considered; hence, the time effect is not resolved. Since tude and phase) for a high- and low-friction coefficient. reflection is not considered, their equations cannot deal Like Hunt, Jay (1991),basedonananalyticalpertur- with closed-end channels. Various approaches have been bation model of the momentum equation for convergent used to arrive at their four equations. According to the channels (including river flow and tidal flats), has shown authors, the combination of different approaches may that a single, incident tidal wave may mimic a standing introduce inconsistencies, which may limit the applicabil- wave by having an approximately 90° phase difference ity of the equations. This may not be a real problem as between the tidal velocities and tidal surface elevations long as measured data sets are available for calibration of and a very large wave speed without the presence of a the tidal parameters. reflected wave. The tidal wave behavior to lowest order is In this paper, analytical solutions of the momentum and dominated by friction and the rate channel convergence. energy equations for converging tidal channels are pre- This result is confirmed by the energy-based approach sented. The basic tidal equations are presented in Section 2. presented herein. Friedrichs and Aubrey (1994)have After that in Section 3, an energy-based approach is presented a first-order solution of tidal wave propagation introduced for converging tidal channels including quadrat- which retains and clarifies the most important properties of ic (nonlinear) bottom friction. It will be shown that the tidal tides in strongly convergent channels with both weak and range (H) can be simply related to the friction parameter strong friction. Their scaling analysis of the continuity and and the rate of channel convergence. Then, an analytical momentum equation clearly shows that the dominant solution of the linearized equations of continuity and effects are: friction, surface slope, and along-channel momentum for converging tidal channels will be presented gradients of the cross-sectional area (rate of convergence). (Section 4). The linearized analytical solution of the Local advective acceleration is much smaller than the equations of continuity and momentum is presented for other contributions. M2 tidal waves with and without reflection in strongly Lanzoni and Seminara (1998) have presented linear and convergent (funnel type) channels. This latter analytical nonlinear solutions for tidal propagation in weakly and solution involves three cases: case A=synchronous tidal strongly convergent channels by considering four limiting channel with constant tidal range, case B=amplified tidal cases defined by the relative strength of dissipation versus range, and case C=damped tidal range. Focusing on local inertia and convergence. In weakly dissipative semidiurnal M2 tides, the present approach applies to channels, the tidal propagation is essentially a weakly midlatitude estuaries. nonlinear problem. As channel convergence increases, the Both approaches (energy and momentum equations) distortion of the tidal wave is enhanced and both the tidal offer simple and powerful tools (spreadsheet models) for wave speed and height increase leading to ebb dominance. tidal analysis of amplified and damped tidal wave propa- In strongly dissipative channels, the tidal wave propagation gation in converging estuaries. The analytical solutions are is a strongly nonlinear process with strong distortion of the compared with the results (Section 5) of measured data of wave profile leading to flood dominance. They use a the Western Scheldt Estuary in the Netherlands, the nonlinear parabolic approximation of the full momentum Hooghly Estuary in India, and the Delaware Estuary in equation. the USA. The analytical solutions of the energy flux Prandle (2003a, b) has presented localized analytical equation and the linearized continuity–momentum equa- solutions for the propagation of a single tidal wave in tions show surprisingly good agreement with measured channels with strongly convergent triangular cross- tidal ranges in these large-scale tidal systems. Section 6 sections, neglecting the advective terms and linearizing presents comparison results of the analytical and numerical the friction term. The solutions apply at any location where models. the cross-sectional shape remains reasonably congruent and It is noted that the linearized solution cannot deal with the spatial gradient of tidal elevation amplitude is relatively the various sources of nonlinearity such as quadratic small (ideal or synchronous estuary). Analyzing the tidal friction, finite amplitude, variation of the water depth under characteristics of some 50 estuaries, he proposed an the crest and trough, effects of river flow, and effects of expression for the bed friction coefficient as function of tidal flats causing differences in wave speed and hence the mud content yielding a decreasing friction coefficient wave deformation (see Jay 1991). In Section 7, these effects for increasing mud content. are discussed for a series of schematized cases based on Finally, Savenije et al. (2008) have presented analytical detailed numerical solutions including all terms. One- solutions of the 1D hydrodynamic equations in a set of dimensional numerical models can be setup easily and four equations for the tidal amplitude, the peak tidal quickly and produce fairly accurate results if the geometry velocity, the wave speed, and the phase difference between and topography is resolved in sufficient detail. Analytical Author's personal copy
1722 Ocean Dynamics (2011) 61:1719–1741 models can only deal with schematized cases, but offer the continuity and motion for depth-averaged flow in a channel advantage of simplicity and transparency. Simple spread- with a rectangular cross-section are: sheet solutions can be made for a quick scan of the b@h @Q parameters involved. The usefulness of the present analyt- þ ¼ 0 ð1Þ @t @x ical approaches is that the effects of human interventions such as channel deepening and widening can be assessed quickly (quick scans). The simple analytical tidal models 1 @Q g @h QjQj can be easily combined with salt intrusion models (see part þ þ ¼ 0 ð2Þ A @t @x C2A2R II), sediment transport models, ecological models, etc. for a quick first analysis of the problems involved. Based on this, in which, A=bho=area of cross-section, b=surface width, the parameter range can be narrowed down substantially so ho=depth to mean sea level (MSL), R=hydraulic radius that the minimum number of detailed numerical model runs ðÞffi ho if b ho , η=water level to mean sea level, Q= need to be made. discharge=Au, u=cross-section-averaged velocity, and C= Chézy coefficient (constant). Similar equations can be derived for compound channels 2 Basic equations of continuity and momentum (Van Rijn 2011).
A clear understanding of the relative importance of the terms of the equations of continuity and momentum can be 3 Energy-based approach obtained by scaling analysis. Based on measured data from three very different tidal systems with depths in the range of The principle of tidal wave amplification defined as the 5–10 m and decreasing widths and cross-sectional areas, increase of the wave height due to the gradual change of the Friedrichs (1993) and Friedrichs and Aubrey (1994) have geometry of the system (depth and width) can be easily convincingly shown that the e-folding length scale Lb of the understood by considering the wave energy flux equation, width (width b=bo exp(−x/Lb)) and the e-folding length scale which is known as Green’s law (1837). This phenomenon is LA of the cross-sectional area (cross-sectional area A=Ao exp also known as wave shoaling or wave funneling. The total (−x/LA)) are of the same order of magnitude and are much energy of a sinusoidal tidal wave per unit length is equal to 2 smaller than the length scale LU of the horizontal velocity E=0.125 ρgbH with b=width of channel and H=wave variations. This implies that the nonlinear convective accel- height. The propagation velocity of a sinusoidal wave is 0.5 eration term (u@u=@x; u ¼ cross-section-averaged velocity) given by co=(gho) with ho=water depth. is negligibly small in the major part of the estuary with Assuming that there is no reflection and no loss of exception of the region close to the landward boundary energy (due to bottom friction), the energy flux (see where a river (open end) or a zero tidal flow boundary Eq. 3) is constant resulting in: Eoco ¼ Excx or Hx=Ho ¼ À0:5 À0:25 (closed end) may be present. In most estuaries with a loose ðÞbx=bo ðÞhx=ho . Thus, the tidal wave height Hx sediment bed, the peak horizontal velocity is almost constant increases for decreasing width and depth. The wave length (of the order of 1 m/s) along the channel due to L=coT will decrease as co will decrease for decreasing depth morphological adjustment, which causes the tidal velocity resulting in a shorter and higher wave. to be rather uniform in space. If horizontal velocity gradients The variation of the tidal range H in a real estuary with will exist, then the channels will not be stable resulting in net exponentially decreasing width and depth can be derived from deposition in regions of relatively low bed–shear stresses and an overall energy-based approach including bottom friction. net erosion in regions with relatively high bed–shear The width and depth are assumed to vary exponentially. stresses. For the channels considered by Friedrichs and Waves transfer energy in horizontal direction during Aubrey (1994), the along-channel gradients in discharge wave propagation. To be able to propagate itself, a wave are dominated by along-channel gradients in cross- must transfer energy to the fluid in rest in front of the wave. sectional area as regards the continuity equation. They The rate at which energy is transferred from one section to also show that the local acceleration term @u=@t is another section is known as energy flux F. Basically, it is relatively small compared to the other terms. The two the time-averaged work done by the dynamic pressure force most important terms of the momentum equation are the per unit time during a wave period, like in short-wave pressure term (water surface slope term) and the theory. Wave energy is extracted from the system through frictional term. the work done by the bed shear stress (causing wave Neglecting the convective acceleration term (u@u=@x, damping). It is realized that the representation of a tidal very small in tidal channels; see Friedrichs and Aubrey wave by a single progressive sinusoidal wave is a crude 1994, Jay 1991, and Van Rijn 1993, 2011), the equations of schematization of a real distorted tidal wave (including Author's personal copy
Ocean Dynamics (2011) 61:1719–1741 1723 reflection) in an estuary. Therefore, this method only yields difference between horizontal and vertical tide, f=8g/C2= a first-order description of the bulk parameters involved friction coefficient, C=Chézy coefficient, and x=horizontal (tidal range H) for conditions without much reflection. coordinate (positive in landward direction). The energy flux balance for depth-integrated tidal flow Neglecting friction (f=0), Eq. 5 yields: dH=dx ¼ 0:5ðÞbþg x reads, as: 0:5ðÞb þ g HorH=Ho ¼ e . The tidal wave height increases exponentially for exponentially decreasing width ð Þ= þ ¼ ð Þ d bF dx bDw 0 3 and depth in line with Green’s law. Equations 5 and 6 include two terms: the effect of depth and width convergence (funneling) and the effect of = þ = þ ¼ ð Þ b dF dx Fdb dx bDw 0 4 quadratic (nonlinear) bottom friction. Amplification will occur if the convergence (first) term is largest. Damping with F =energy flux per unit width per unit time and Dw= will occur if the friction (second) term is largest. To produce energy dissipation per unit area and time due to bottom the tidal range along the channel, it can be simply solved in friction, b=width of estuary channel, x=horizontal coordi- b a spreadsheet when u, φ, ho, f (or C) and the boundary nate (positive in landward direction, x=0=mouth). The unit b condition H=Ho at x=0 are known. The u and φ 3 of each term is kilogram per m per cubic second (kg m/s ). parameters will be taken from the linearized solution of The energy flux can be expressed as (see Appendix I): the equations of continuity and momentum (see Table 1 of ¼ : 2 ¼ ¼ : 2 F 0 125 rgH co Eco with: E 0 125 rgH=energy Section 4) for converging channels. of a wave per unit area, co=wave propagation velocity. In the case of a compound channel (consisting of main Thus, the energy flux per unit width of a progressive wave channel and tidal flats), the depth ho should be replaced by in deep water is equal to the energy of a wave per unit the wave propagation depth heff ¼ A=bs ¼ ahho (Van Rijn length of the wave and the wave propagation velocity 2011). Furthermore, the friction coefficient will be larger (similar to the expression used in short wave theory). and the wave speed will be smaller in a compound channel. The instantaneous energy dissipation per unit area and Savenije (2005) has given an expression for the tidal – time due to work done by the bed shear stresses can be range which is based on an analytical solution of the ¼ =ðÞðÞ= 2 b3 expressed (see Appendix I): Dw 4 3p r g C u simplified momentum equation for a converging estuary The width and depth of the estuary channel is repre- using a Lagrangian coordinate system. The momentum sented as: equation is applied at times of high water (HW) and low water (LW) and then subtracted to determine dH/dx. His b ¼ boexpðÞ¼À bx boexpðÞÀx=Lb and thus db=dx ¼Àbb expðÞ¼ÀÀ bx b b final result reads as: o ¼ ðÞ¼À ðÞÀ = h hoexp gx hoexp x Lh and thus dH=dxðÞ¼ 1 þ ! bH À f =ðÞH=h bu=c cos 8 ð7Þ dh=dx ¼ÀghoexpðÞ¼ÀÀgx g h with c=wave speed, bu=peak tidal current velocity, φ= with, b =width at entrance x=0, β=1/L =width convergence o b phase lead angle (constant), h=mean depth, β=1/L =width coefficient, γ=1/L =depth convergence coefficient, L = b h b convergence coefficient, f /=g/C2=friction coefficient and convergence length scale for width, L =convergence length h a ¼ gH= 2cbu cos8 . The assumption of Savenije (2005): scale for depth. The length scale L is of the order of 10– b dbu/dx=(bu/H) dH/dx leads to the α term, which is not 50 km for most estuaries. The length scale L is much larger h present in Eq. 6. Savenije et al. (2008) have improved their than L for most estuaries as the depth generally is fairly b method by presenting a set of four equations for the constant or very weakly decreasing in landward direction. variables: H, bu, c, and φ. Using these expressions (see Appendix I), the solution of Equation 6 can be used to estimate the length of a Eq. 4 yields: (synchronous) estuary with a constant tidal range (dH/dx=0) 2 ¼ = ¼ ðÞ2 2 8 = 3 f bu yielding b 1 Lb fH c cos 6p gho H . dH=dx ¼ 0:5ðÞb þ g H À ð5Þ 2 8 Using c =gho, it follows that the converging length scale is: 3pgh cos 2 Lb ¼ 6pho =ðÞfHcos 8 for a synchronous estuary. As- or as γ=0: suming that the length L of the synchronous part of an alluvial 2 2 2 8 estuary is about L ffi 2Lb ffi 12 p ho =ðÞfHcos 8 .Using = ¼ : À fHc cos ð Þ dH dx 0 5 bH 3 6 f≅0.025 and cosφ≅0.7, yields that L≅2,000 (ho/H) ho. 12pghðÞ0 The ratio H/ho is about 0.1–0.4inthemouthofmostestuaries ≅ – – with, H=tidal range, β=1/Lb=width convergence coefficient and thus L (5,000 20,000)ho.Usingho=1 20 m, the length b (Lb=e-folding length scale), ho=depth (constant), u=peak of funnel-shaped, synchronous, alluvial estuaries will be in tidal velocity along channel, c=local wave speed, φ=phase the range of 5–400 km. Prandle (2004) has found 1724 Table 1 Analytical solutions for prismatic and converging channels (with and without reflection)
Type of wave Prismatic channels Converging channels
Excluding reflection at landward end (channel open at end) ¼ b ½Àmx ½ðÞÀ ¼ b ½À"x ½ðÞÀ hx;t ho e cos wt kx hx;t ho e cos wt kx b Àmx b À"x ux;t ¼ uo½e cosðÞwt À kx þ 8 ux;t ¼ uo½e cosðÞwt À kx þ 8
with with b ¼ÀðÞb = ðÞ= ½8 " ¼À : þ uo ho ho w k cos 0 5b m b ¼ÀðÞb = ðÞ= ½8 μ=friction parameter uo ho ho w k cos k=wave number c ¼ w=k; tan 8 ¼ sinhi8= cos 8 ¼ ðÞ0:5b þ m =k; 0:5 φ=phase lead sin 8 ¼ ðÞ0:5bhiþ m = ðÞ0:5b þ m 2 þ k2 0:5 2 tan φ=μ/k cos 8 ¼ðkÞ= ðÞ0:5b þ m þ k2 ; μ x=horizontal coordinate; positive in landward =friction parameter Author's direction k=wave number φ=phase lead (between horizontal and vertical tide)
β=1/Lb =convergence parameter x=horizontal coordinate; positive in landward direction h h personal ¼ : ðÞÀ1 ÀmðÞxÀL ð À ¼ : b ðÞÀ1 ðÞ"1 xþmL ð À Including reflection at landward end (channel closed at end) hx;t 0 5ho fA e cos wt hx;t 0 5ho fA e cos wt i i kxðÞÞþÀ L emðÞxÀL cosðÞwt þ kxðÞÀ L kxðÞÞþÀ L eðÞ"2 xÀmL cosðÞwt þ kxðÞÀ L À À0:5 ¼ : ðÞb = ðÞ1 2 þ 2 À0:5 ux;t 0 5w ho ho fA k m ¼ : ðb = Þð ÞÀ1ð 2 þ 2Þ h ux;t 0 5w ho ho fA k m copy ÀmðÞxÀL ðÞÀ ðÞþÀ 8 ð" Þ þ e cos wt kx L ½e 1 x mL cosðwt À kðx À LÞþ8Þ i ðÞÀ ð"2ÞxÀmL ~ Àem x L cosðÞwt þ kxðÞþÀ L 8 Àe cosðwt þ kðx À LÞþ8Þ with with
fA=amplification/damping factor= ε1 =0.5 β−μ =[cos2(kL)+sinh2(μL)]0.5 ε =0.5β+μ 2 61:1719 (2011) Dynamics Ocean L=channel length fA=amplification/damping factor= x=horizontal coordinate; positive in landward =[cos2(kL)+sinh2(μL)]0.5 direction L=channel length, x=horizontal coordinate, positive in landward direction – 1741 Author's personal copy
Ocean Dynamics (2011) 61:1719–1741 1725
1:25 0:5 L ffi 2; 500ðÞho =H . He presents a set of 50 estuaries The water surface slope can be expressed as: @h=@x ¼ : Àmx 2 2 0 5 (UK and USA) with lengths of 2–200 km and depths of 1– bn0ðÞe ðÞm þ k cosðÞwt À kx þ y with: tan y ¼ k=m 20 m. or cotany ¼ m=k. Since tan 8 ¼ m=k (see Table 1), it follows that cotan==tanφ or y ¼ 900 À 8.Thus,the velocity has a phase lag of y À 8ðÞor 90o À 28 with respect 4 Analytical and numerical solutions of equations to the surface slope ∂η/∂x.Thisphaselagvariesfromsmall of continuity and momentum values (10° or 0.3 hours) in shallow, friction-dominated conditions (quasi-steady flow; velocity reacts directly to a 4.1 Prismatic channels decreasing water level in the direction of the flow) to about 90° in deep water (velocity responds delayed to water Neglecting the convective acceleration term u@u=@x , the surface slope). equations of continuity and momentum for a wide, Figure 2 shows the velocity amplitude and the phase prismatic channel with rectangular cross-section and con- lead angle φ1 (in hours) as a function of the water level stant mean depth ðÞh ¼ ho þ h can be linearized by amplitude at x=0 m for three water depth values h=5, 10, replacing the quadratic friction term by a linear term (based and 20 m in an open-end system (without reflection). The on the Lorentz (1922, 1926) method), as follows: bed roughness (ks) value is 0.05 m. The peak velocities (flood) at the mouth are in the range of 0.5–1.5 m/s for the @h @ þ h0 u ¼ ð Þ smallest water depth of 5 m and in the range of 0.35–2.5 m/s @ @ 0 8 t x for the largest depth of 20 m. The effect of the water depth is rather small. Roughly, the peak velocities (which are very @u g @h þ þ m bu ¼ 0 ð9Þ realistic values for semi-diurnal estuaries) are in the range of @t @x 0.4–1.7 m/s for for bn0 =0.5–3 m. The phase lead angles are in which, u=cross-section-averaged velocity, u=depth- between 1 and 1.4 hours for h=5 m and between 0.25 and b averaged velocity, u=amplitude of close section-averaged tidal 1hourforh=20 m. The phase lead angle decreases with current velocity, m=friction coefficient= 8g bu = 3pC2R , increasing depth (less bottom friction).
R=hydraulic radius (R≅ho and u≅u for a very wide channel For a closed channel at one end, the solution represents a b>>ho). Using the Lorentz method, the total linear friction standing wave system characterized by incident and reflected over the tidal cycle is the same as that for quadratic friction. It waves of equal amplitude which individually propagate at 0.5 is noted that the friction coefficient depends on the peak speed c (smaller than the classical value co=(gho) ), see velocity (usually taken at x=0). The classical analytical Table 1. Resonance phenomena do occur when the channel solution of these two equations with and without reflection length is of the order of 0.25 L with L=tidal wave length. at the head of the channel is given in Table 1 (Dronkers 1964, 2005;Verspuy1985; Van Rijn 1993, 2011). The solution is an 4.2 Convergent (funnel type) channels exponentially damped wave due to bottom friction. The phase lead φ of the velocity with respect to the water level can be The topography of most estuary channels is sufficiently expressed as tanφ=μ/k with μ=damping coefficient, and k= complex that any realistic representation requires the division wave number. The wave speed c is always smaller than the of the channel into a series of prismatic channels each with its 0.5 classical value co=(gho) . own constant width b. The classical analytical approach is to
Fig. 2 Peak tidal velocity at 3 mouth and phase lead angle as Velocity amplitude (h = 5 m) Velocity amplitude (h = 10 m) function of water level ampli- 2.5 Velocity amplitude (h = 20 m) Phase lead (h = 5 m) tude at mouth (0.5H) and water Phase lead (h = 10 m) depth (negative flood velocity is Phase lead (h = 20 m) 2 plotted as positive value); pris- ks = 0.05 m matic channel without reflection 1.5
1 phase lead (hours)
0.5 Tidal velocity amplitude (m/s) and and amplitude (m/s) Tidal velocity
0 0123456 Water level amplitude (m) Author's personal copy
1726 Ocean Dynamics (2011) 61:1719–1741 apply the equations for a prismatic channel to each channel @ @ section from seaward to landward direction using the output h þ ho u À " ¼ ð Þ @ @ u ho 0 12 of the previous section as input for the next section (Dronkers t x 1964). Unfortunately, this method eliminates the large effect of The momentum equation can be simplified to: channel convergence on the complex wave number and thus @u g @ h onthewavespeedc (Jay 1991). Furthermore, this segment þ þ m u ¼ 0 ð13Þ @t @x type of approach is laborious and offers no real advantage b 2 b compared with a direct 1D numerical solution method. with m ¼ nA ¼ nbho ¼ 8g u =ðÞffi3pC R 8g u = Therefore, another approach is required which takes the 2 ð3pC hoÞ¼ Lorentz friction parameter (dimension 1/s), u= channel convergence properly into account. This can be cross-section-averaged velocity, bu=characteristic peak velocity done by using an exponential function for the channel (average value over traject), C=Chézy=coefficient, R= width allowing an analytical solution of the equations of hydraulic radius. The river velocity ur ¼ Qr=A can easily be continuity and momentum. The width of the channel is included. This leads to additional friction and more tidal ¼ Àbx ¼ = ¼ described by: b boe with b 1 Lb convergence damping. coefficient, Lb=converging length scale (constant in time). In the case of a compound channel consisting of a The analytical solutions with and without reflection are main channel and tidal flats, it may be assumed that the given in Table 1 (Van Rijn 2011). flow over the tidal flats is of minor importance and only The basic assumptions are: bottom is assumed to be contributes to the tidal storage. The discharge is horizontal, channel depth ho to MSL is assumed to be conveyed through the main channel. Inclusion of tidal ¼ À"x " ¼ constant in space and time, width is b boe with flats in a 1D approach requires that the equations of = ¼ 1 Lb convergence coefficient, Lb=converging length continuity and momentum are modified to account for @ =@ ¼ scale, convective acceleration u u x 0 is neglected; the different widths of the main channel and the flats linearized friction is used and x-coordinate is positive in (Van Rijn 2011). In Eq. 11, the multiplier of ∂η/∂tshould landward direction and negative in seaward direction (flood be bs (=surface width) rather than b resulting in co=(αb velocity is positive and ebb velocity is negative). 0.5 gho) with ab ¼ bc=bs ¼ width coefficient (in the range The friction term is linearized as follows (Van Rijn 2011): 0.5–1), bc=width of main channel, and bs=surface width. @ @ The transfer of momentum from the main flow to the flow 1 Q þ g h þ ¼ ð Þ @ @ nQ 0 10 over the tidal flats can be seen as additional drag exerted A t x on the main flow (by shear stresses in the side planes in which, n is a constant friction factor, n ¼ 8g Qb between the main channel and the tidal flats). This effect =ðÞ3pC2A2R . can be included crudely by increasing the friction in the The mass balance equation can be expressed as ðA ¼ main channel through a smaller effective depth heff=αhho b h and h ¼ ho þ hÞ: with αh in the range of 0.8–1 (to be used in the friction b @h bh @u ub @h uh @b term). The width of the mainchannelandthesurface þ o þ þ o ¼ 0 ð11Þ @t @x @x @x widtharebothassumedtovaryexponentiallyoverthe length of the channel (α =constant). @ =@ ¼À À bx ¼À b The gradient of the width is b x b bo e b b. The (laborious) analytical solution based on complex ∂η ∂ Neglecting the term ub / x, the mass balance equation functions is given by Van Rijn (2011) and is summarized in becomes: Table 1. The general solutions for k and μ are: "# : 0:5 2 0 5 À0:5 2 2 2 k ¼ð2Þ ko þ1 À ðÞ0:5 b=ko þÀ1 þ ðÞ0:5 b=ko þ ðÞm=w ð14Þ
"# : 0:5 2 0 5 À0:5 2 2 2 2 ¼ð2Þ ko À1 þ ðÞ0:5"=ko þÀ1 þ ðÞ0:5"=ko þ ðÞm=w ð15Þ
0.5 with ko ¼ w=co and co=(gho) Author's personal copy
Ocean Dynamics (2011) 61:1719–1741 1727 2 This involves three cases (for details, see Van Rijn 2011): 2 ¼ ðÞmw = 2 co ðÞ1=k 2 can also be expressed as : ð Þ Case A: Amplitude of water level and velocity are constant hi20 0:5 " ¼ 0:5 0:5 2 2 2 0:5 (ideal estuary or synchronous) H=Ho 2 ¼ 0:5 ðÞw=co ! À1 þ 1 þ m = w ! 0.5 2w=co ¼ 2ko or 1=Lb ¼ 2ko and c=co=(gho) ; convergence and friction are about equal. The actual propagation velocity c is: Case B: Amplification is dominant (hypersynchronous) L w c H>H with " 32w=c or 1=L > 2k and c>c ; c ¼ ¼ ¼ o ð21Þ o o b o o T k ðÞ: ! 0:5ðÞþ 0:5 inertia and convergence dominate bottom friction 0 5 1 a2 : hi : 2 2 2 0 5 0:5 0 5 with a2 ¼ ðÞ1 þ m =ðÞw ! k ¼ 0:50:5ðÞw=c !0:5 À1 þ 1 þ m2= w2!2 o This latter expression yields a wave propagation velocity with ! ¼ 0:25"2ðÞc =w 2 À 1 which is smaller than that of a frictionless wave (c