2.0.0b10
catchment modelling framework
Loading...
Searching...
No Matches
Is hydrology kinematic?

... is the title of a paper by Singh (2002), where he lists numerous examples from all scales of hydrology where a kinematic wave function has been found to be a valid descriptor for hydrological processes. Anyway, numerous more or less conceptional models in hydrology use kinematic wave like formulations to calculate fluxes. While some authors are happy to call every 1st order flux calculation using some exponent a kinematic wave, other authors like to see the term reserved for more "wave like" phenomens in hydrology.

Hence a number of different names exist for the same concept:

  • Kinematic wave
  • Power law equation (eg. DHSVM, SUPERFLEX)
  • curve shape parameter (eg. MACRO, HBVlight)

A simplified equation for power law equations

A generalized form for some approaches is:

\[ q(V) = Q_0 \cdot \theta^\beta\ (1) \]

where:

  • \(\theta\) is a dimensonless variable representing a water content related state of a storage
  • \(\beta\) is a curve shape parameter
  • \(Q_0\) is the flux in m³/day at unity state \(\theta=1\)

For cmf, using the finite volume approach, we need a consistent function to relate the stored water volume \(V\) in a water storage to the variable \(\theta\). So we define:

\[ \theta = \frac{V - V_{res}}{V_0}\ (2) \]

where:

  • \(V\) is the water volume stored in a storage in \(m^3\)
  • \(V_{res}\) is an immobile water volume in the storage, that does not take place in this connection. It may be ommitted and set to \( 0 m^3\).
  • \(V_0\) is a reference volume, where, for any \(\beta\) the runoff from the water storage is \(Q_0\).

The famous linear storage equation can be seen as a special case of this approach where \(\beta=1\). In that case, we can rearrange the formula above to:

\[ q(V) = \frac{Q_0}{V_0} \cdot (V-V_{res})\ (3) \]

Here the quotient \(\frac{Q_0}{V_0}\) becomes the inverse of the storage residence time \(t_r^{-1}\), which simplifies (3) to

\[ q(V) = \frac 1 {t_r [days]} \cdot (V-V_{res})\ (4) \]

for the linear case

Implementation in cmf

Due to the different meanings of parameters between the linear and non-linear case, cmf provides now 2 distinct connections for the linear and the non-linear case.

LinearStorageConnection

ls = cmf.LinearStorageConnection(source, target, residencetime=3.2, residual=0)

The linear storage connection LinearStorageConnection implements (4) and takes the residence time in days and a residual water content in m³ as parameters. While the residence time must be given, the residual water content can be omitted and set to 0. Depending on system size and structure, residence times between minutes for surface water storages or fast drainages up to hundreds of years for large aquifers are suitable values.

PowerLawConnection

plc = cmf.PowerLawConnection(source, target, Q0=3.0, V0=5.3, beta=2.0, residual=0)

The power law connection (PowerLawConnection) needs the parameters Q0, V0 and beta as given in the generalized form of the power law equation. The reference flow \(Q_0 [\frac{m^3}{days}]\) occurs when the source water storage contains the reference volume \(V_0\) as mobile water. The curve shape parameter \(\beta\) bends the curve for mobile water volumes between 0 and \(V0\) down, when \(\beta>1\). The curve is bend up for \(\beta<1\). If \(V - V_{res}>V_0\), the effect is directly turned around, as shown in the figure.

The residual water content \(V_{res} [m^3]\) shifts the storage/runoff response to the right.

A residual water content makes sense eg. for percolation in simplified models where percolation starts only above field capacity, where the residual water content equals field capacity of the top soil, or it can be used as the volume behind a dam. To find a justification for \(\beta\) is more challenging: \(\beta>1\) describes systems, where flow is over linear slow for low storages and gets much faster for a large storage. Unsaturated flow can be modelled that way, if one takes an analogy with the Brooks-Corey retention curve, betas from 4 to 15 are reasonable. The standard beta for the "real" kinematic wave from the Manning-Strickler equation for sheet flow in \(\frac 5 3\). If you find a real world usage for \(\beta<1\), please inform the authors of cmf. However, if \(\beta<0.3\) it is reported, that numerical instabilities will arise.