bf_part: Baseflow Separation

Description Usage Arguments Value Note References Examples

View source: R/part.R

Description

Extract baseflow from a daily streamflow record using the method described by Rutledge (1998).

Usage

1
bf_part(date, discharge, da, STAID = "Unknown")

Arguments

date

the date for each x, should be of class "Date." Missing values are not permitted.

discharge

the daily streamflow to be separated missing values are not permitted within the time specified by Start and end.

da

the drainage area of the basin in square miles.

STAID

the station identifier for the data.

Value

an object of class "baseflow" and inherits class "data.frame" of the selected data, a data frame of the baseflow information, and other information about the analysis.

Note

The estimates from this routine will occasionally differ slightly from the original FORTRAN code from Rutledge (1998). It is expected that those small differences are due to rounding differences affecting comparison between numbers that differ by very small amounts. There are also differences in the initial baseflow estimates between this version and the original version. Those differences are due to slightly different initialization routines.

The estimate of baseflow from this routine is computed from a linear interpolation of the Part1 and Part2 estimates rather than the curvilinear interpolation that was proposed without detials in Rutledge (1998).

References

Rutledge, A.T., 1998, Computer programs for describing the recession of ground-water discharge and for estimating mean ground-water recharge and discharge from streamflow records—Update: U.S. geological Survey Water-Resources Investigations Report 98-4148. 43 p.

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
library(smwrData)
data(ChoptankFlow)
# Process by calendar year as that is the retrieval range
ChopPart <- with(ChoptankFlow, part(discharge, datetime, da=113,
STAID="01491000"))
ChopPart

## End(Not run)

smwesten-usgs/recharge documentation built on Jan. 27, 2020, 3:08 p.m.