baseflowB {hydroEvents}

Baseflow removal (after Ladson et al)

Description

This function calculates baseflow using a recursive digital filter and is based on the implementation described in Ladson et al (2013).

Usage

baseflowB(q, alpha = 0.925, passes = 3, r = 30)

Arguments

q

The vector series of streamflow

alpha

Filter parameter

passes

Number of passes

r

number of points reflected at start and end of data set

Details

The reflected points act to resolve spin up issues and are removed before the baseflow is removed.

Value

A list of the baseflow and baseflow index at each timestep.

References

Ladson, A., Brown, R., Neal, B., & Nathan, R. (2013). A standard approach to baseflow separation using the Lyne and Hollick filter. Australian Journal of Water Resources, 17(1).

Examples

library(hydroEvents)
data(dataBassRiver)
alpha.list = c(0, 0.9, 0.925, 0.95, 0.98, 0.987)
BFI = numeric(length(alpha.list))
for (i in 1:length(alpha.list)) {
  bf = baseflowB(dataBassRiver, alpha = alpha.list[i])
  BFI[i] = sum(bf$bf)/sum(dataBassRiver)
}
print(cbind(alpha.list, BFI))