uh: Unit hydrograph estimation

Description Usage Arguments Details Value Examples

View source: R/uh.R

Description

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

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
uh(hl, ...)

## Default S3 method:
uh(hl, uc, deltat, ...)

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

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

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

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).

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.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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$time,"h"), cumsum(uh1$prob), type = "b",
xlab = "Travel~time", ylab = "Cumulative~distribution")

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

transfR documentation built on Sept. 8, 2021, 9:08 a.m.