accumulate_runoff_muskingum_cunge: Apply Muskingum-Cunge river routing

In mkkallio/hydrostreamer: A Package for Interpolating Distributed Runoff Products on to Explicit River Segments

Apply Muskingum-Cunge river routing

Description

The function implements Muskingum routing scheme where the storage parameter k is computed using user input wave celerity, and the length of a river segment. Using Muskingum for runoff data with time interval day may cause instability in the output. If the interval is too high, it is advised to use another routing algorithm.

Usage

accumulate_runoff_muskingum_cunge(
  HS,
  manning,
  slope,
  channel_width = "Moody",
  q_ref = NULL,
  celerity = NULL,
  verbose = FALSE
)

Arguments

HS	A HS object.
manning	Manning's roughness coefficient. Either a constant, a vector of values corresponding to each river segment, or the name of the column of HS where the values are stored.
slope	River bed slopes. Either a constant, a vector of values corresponding to each river segment, or the name of the column of HS where the values are stored.
channel_width	Channel widths. See details.
q_ref	Reference discharge for which channel width computations are based on. The default is NULL, which means that q_ref is automatically computed from inflow timeseries.
celerity	If desired, (constant) wave celerities for each river segment can be given. Default NULL means that wave celerity is automatically computed for each river segment and each timestep.
verbose	Whether or not to print progress information. Defaults to FALSE.

Details

Muskingum hydrologic routing method consists of two parameters; K and x. K corresponds to travel time through a river segments, and can be estimated using Wave celerity and length of the flow line. x is an empirical parameter for which the value can be found through trial and error. hydrostreamer currently does not offer automated calibration of x.

Channel widths can be given as a vector, however by default channel width is computed from a power-law relationship from reference discahrge (q_ref). If the length(channel_width) == 2, the first parameter is alfa and the second is beta in alfa*q_ref^beta. "Allen" and "Moody" correspond to optimised parameters.

Value

Returns the input object HS) with an added list column discharge_ts containing routed discharge estimates for each river segment.

mkkallio/hydrostreamer documentation built on Oct. 14, 2023, 9:38 p.m.