bf_part: Baseflow Separation

In smwesten-usgs/recharge: Estimating Recharge with HYSEP, PART, and Recursive Digital Filter methods.

Description

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

Usage

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

## 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 March 8, 2021, 11:51 a.m.