Learn R Programming

transfR (version 1.1.4)

uh: Unit hydrograph estimation

Description

Estimate the unit hydrograph from a sample of hydraulic lengths and a streamflow velocity.

Usage

uh(hl, ...)

# S3 method for default uh(hl, uc, deltat, ...)

# S3 method for units uh(hl, uc, deltat, ...)

# S3 method for stars uh(hl, ...)

# S3 method for transfR uh(hl, verbose = TRUE, ...)

Value

A data.frame with vectors of class units, or if hl is a transfR object, the same transfR object incremented by the "uh" attribute.

Arguments

hl

hydraulic length of class stars, matrix, vector or transfR. If no unit is provided, hl is assumed to be in [m].

...

further arguments passed to or from other methods

uc

streamflow velocity. If no unit is provided, uc is assumed to be in [m/s].

deltat

time step of the time series. If no unit is provided, deltat is assumed to be in [min].

verbose

boolean indicating if information messages should be written to the console

Details

The function estimates the unit hydrograph from geomorphometric information. A travel time to the outlet is estimated by assuming an average streamflow velocity (uc) within the river network and by applying uc over the sample of hydraulic lengths (hl). The unit hydrograph is the probability distribution of this travel time to the outlet given at each time step (deltat).

Examples

Run this code
data(Oudon)
uh1 <- uh(hl=Oudon$hl[[1]], uc=units::set_units(0.5,"m/s"),
deltat=units::set_units(1,"h"))
plot(units::set_units(uh1$max_time,"h"), cumsum(uh1$prob), type = "b",
xlab = "Travel~time", ylab = "Probability~of~non-exceedance")

object <- as_transfr(st = Oudon$obs, hl = Oudon$hl)
object <- velocity(object)
object <- uh(object)
plot(object, i = 1, attribute = c("uh"))

Run the code above in your browser using DataLab