part: Baseflow Separation

In USGS-R/DVstats: Functions to manipulate daily-values data

Baseflow Separation

Description

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

Usage

part(Flow, Dates, Start = NULL, End = NULL, da, STAID = "Unknown")

Arguments

Flow	the daily streamflow to be separated missing values are not permitted within the time specified by Start and end.
Dates	the date for each x, should be of class "Date." Missing values are not permitted.
Start	the start date for the analysis, can be either a character string or class "Date."
End	the end date for the analysis, can be either a character string or class "Date."
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(Flow, datetime, da=113,
STAID="01491000"))
ChopPart

## End(Not run)

USGS-R/DVstats documentation built on Oct. 11, 2022, 6:03 a.m.