Last data update: 2014.03.03

R: Baseflow Separation
BaseflowSeparationR Documentation

Baseflow Separation

Description

This function reads a streamflow dataset and produces a baseflow dataset. It can be run using 1, 2 or 3 passes

Usage

BaseflowSeparation(streamflow, filter_parameter = 0.925, passes = 3)

Arguments

streamflow

A vector containing streamflow values

filter_parameter

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 
>