baseflowB function

Baseflow removal (after Ladson et al)

Baseflow removal (after Ladson et al)

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

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

Returns

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

Details

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

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

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

hydroEvents packageRead PDF manual