Calculate measure of central tendency and baseflow indices using the Lynne-Hollick filter
baseflows(flow.ts, a, n.reflected = 30, ts = "mean")
Dataframe with date and discharge data in columns named "Date" and "Q" respectively. Date must be in POSIX format (see ts.format). Missing values are ignored.
The alpha value used in the Lynne-Hollick filter for digital baseflow separation. Default value is 0.975
The number of days that are reflected at the start and end of the series to provide a burn in for the digital filter. Default value is 30. (See Ladson et al. 2013).
ts="mean" returns means for the entire time series ts="annual" returns annual averages. Note this function does not currently use hydrologic years even when defined ts="daily" returns complete series
A dataframe. See below for details. The original dataframe with appended columns "bf" and "bfi". See ts="annual" for details.
ts="mean"
The number of years of record in the series
proportion of non-missing observations
mean daily flow
median daily flow
mean baseflow volume
mean baseflow index
the record year
no of observations in each year
mean daily flow in each year
mean baseflow volume in each year
baseflow index for each year
baseflow index for each observation
baseflow index associated with each observation
Technically the LH filter cannot be calculated where there are missing data. Here the function removes missing values and is applied to a concatenated version of the time series. Missing dates are reinserted after the filter has been applied for the purpose of returning annual or daily series. The function further reports the number of missing values leaving the user to decide on the reliability of the baseflow estimates.
Ladson, A. R., R. Brown, B. Neal and R. Nathan (2013) A standard approach to baseflow separation using the Lyne and Hollick filter. Australian Journal of Water Resources 17(1): 173-18
Lynne, V., Hollick, M. (1979) Stochastic time-variable rainfall-runoff modelling. In: pp. 89-93 Institute of Engineers Australia National Conference. Perth.
# NOT RUN {
data(Acheron)
Acheron<-ts.format(Acheron)
baseflows(Acheron,a=0.975, ts="mean")
baseflows(Acheron,a=0.975, ts="annual")
head(baseflows(Acheron,a=0.975, ts="daily"))
# }
Run the code above in your browser using DataLab