The value recommended by Nathan and McMahon (1990) is 0.925, however, the user might want to play with this value (0.9-0.95)
passes
The number of times you want the filter to pass over your data. 1-3
Value
This will return a 2 column data frame with nrow = length of input streamflow data. The first column contains
baseflow, while the second contains quickflow, both in the same units as the input.
Author(s)
Josephine Archibald
References
Lyne, V. D. and M. Hollick (1979). Stochastic time-variable rainfall-runoff modelling. Hydrology and Water Resources Symposium, Perth, Institution of Engineers, Australia.
Nathan, R. J. and T. A. McMahon (1990). "Evaluation of automated techniques for base flow and recession analysis." Water Resources Research 26(7): 1465-1473.
Examples
######### Look at a dataset for Owasco Lake in NY:
data(OwascoInlet)
summary(OwascoInlet)
## Get an approximation for baseflow using a 3 pass filter:
bfs<-BaseflowSeparation(OwascoInlet$Streamflow_m3s, passes=3)
## You can check out how this looks with the hydrograph function:
hydrograph(input=OwascoInlet,streamflow2=bfs[,1])
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(EcoHydRology)
Loading required package: operators
Attaching package: 'operators'
The following objects are masked from 'package:base':
options, strrep
Loading required package: topmodel
Loading required package: DEoptim
DEoptim package
Differential Evolution algorithm in R
Authors: D. Ardia, K. Mullen, B. Peterson and J. Ulrich
Loading required package: XML
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/EcoHydRology/BaseflowSeparation.Rd_%03d_medium.png", width=480, height=480)
> ### Name: BaseflowSeparation
> ### Title: Baseflow Separation
> ### Aliases: BaseflowSeparation
>
> ### ** Examples
>
> ######### Look at a dataset for Owasco Lake in NY:
> data(OwascoInlet)
> summary(OwascoInlet)
date P_mm Streamflow_m3s baseflow_m3s
Min. :2009-04-22 Min. : 0.000 Min. : 0.310 Min. : 0.000
1st Qu.:2009-11-29 1st Qu.: 0.000 1st Qu.: 1.590 1st Qu.: 1.329
Median :2010-07-09 Median : 0.000 Median : 3.060 Median : 2.516
Mean :2010-07-17 Mean : 3.121 Mean : 4.647 Mean : 3.212
3rd Qu.:2011-03-16 3rd Qu.: 2.540 3rd Qu.: 5.077 3rd Qu.: 3.960
Max. :2011-10-26 Max. :112.522 Max. :56.350 Max. :15.737
Tmax_C Tmin_C
Min. :-12.220 Min. :-26.110
1st Qu.: 8.193 1st Qu.: -1.670
Median : 18.330 Median : 6.670
Mean : 15.977 Mean : 5.059
3rd Qu.: 25.000 3rd Qu.: 12.220
Max. : 36.670 Max. : 23.330
>
> ## Get an approximation for baseflow using a 3 pass filter:
> bfs<-BaseflowSeparation(OwascoInlet$Streamflow_m3s, passes=3)
>
> ## You can check out how this looks with the hydrograph function:
> hydrograph(input=OwascoInlet,streamflow2=bfs[,1])
>
>
>
>
>
> dev.off()
null device
1
>