baseflow | R Documentation |
Given a stream flow hydrograph of flows (regular time series), the base flow is separated. The minima of a period (default block.len = 5)
is calculated and turning points are identified. At turning points the base flow equals the actual flow, in between, linear interpolation is carried out.
baseflow(x, tp.factor = 0.9, block.len = 5)
x |
numeric vector containing flows |
tp.factor |
numeric vector of length one. Towards high flows, allow the
central value of three consecutive minima only to be of a factor
|
block.len |
numeric vector of length one. |
A numeric vector of length(x)
. It contains NA
s as until the
first turning point, the base flow 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.
## reproducing Tallaksen and van Lanen (2004) ## Example 5.3 Base Flow Index" data(ray) ray <- as.xts(ray) # calculate base flow 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.