Given a stream flow hydrograph of flows (regular time series), the baseflow is separated. The minima of a period (default
block.len = 5) is calculated and turning points are identified. At turning points the baseflow equals the actual flow, in between, linear interpolation is carried out.
baseflow(x, tp.factor = 0.9, block.len = 5)
numeric vector containing flows
numeric vector of length one. Towards high flows, allow the
central value of three consecutive minima only to be of a factor
numeric vector of length one.
A numeric vector of
length(x). It contains
NAs as until the
first turning point, the baseflow cannot be determined.
Tallaksen, L. M. and Van Lanen, H. A. J. 2004 Hydrological Drought: Processes and Estimation Methods for Streamflow and Groundwater. Developments in Water Science 48, Amsterdam: Elsevier.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
## reproducing Tallaksen and van Lanen (2004) ## Example 5.3 Base Flow Index" data(ray) ray <- as.xts(ray) # calculate baseflow and plot it ray$baseflow <- baseflow(ray$discharge) ray96 <- ray[format(time(ray), "%Y") == "1996", ] plot(ray96$discharge, type = "l") lines(ray96$baseflow, col = 2) # aggregated base flows for river Ray # these are mean flow totals per day, not per year as written # in Tallaksen and van Lanen (2004) round(colSums(ray96[, c("discharge", "baseflow")]), 2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.