intermittent_pumping: Streamflow depletion for an intermittent pumping schedule...

In szipper/streamDepletr: Estimate Streamflow Depletion Due to Groundwater Pumping

Streamflow depletion for an intermittent pumping schedule using superposition.

Description

Streamflow depletion for an intermittent pumping schedule using superposition.

Usage

intermittent_pumping(t, starts, stops, rates, method = "glover", d, S, Tr, ...)

Arguments

t	times you want output for [T]
starts	vector of times to start pumping [T] (must be same length as stops and rates)
stops	vector of times pumping stops [T] (must be same length as starts and rates)
rates	vector of pumping rates [L3/T] (must be same length as starts and stops)
method	analytical solution to use (options= 'glover', 'hunt'). Defaults to 'glover'.
d	distance from well to stream [L]
S	aquifer storage coefficient (specific yield if unconfined; storativity if confined)
Tr	aquifer transmissivity [L2/T]
...	any other inputs required for your method of choice; for example, hunt needs lmda (streambed conductance)

Details

This function superimposes wells and image wells to calculate a timeseries of streamflow depletion. Unlike the streamflow depletion models (e.g. glover, hunt) this is not fractional depletion (Qf) because there can be different pumping rates at different times.

Value

A numeric of Qs, streamflow depletion rate [L3/T].

References

Jenkins, C.T. (1968). Techniques for Computing Rate and Volume of Stream Depletion by Wells. Ground Water 6(2): 37-46. doi:10.1111/j.1745-6584.1968.tb01641.x

Examples

Qs <- intermittent_pumping(t = seq(0, 1000, 5),
 starts = seq(0, 900, 10), stops = seq(9, 909, 10), rates = seq(1, 1000, length.out=91),
 method = "hunt", d = 100, S = 0.1, Tr = 100, lmda = 10)

Qs <- intermittent_pumping(t = seq(0, 1000, 5),
 starts = seq(0, 900, 10), stops = seq(9, 909, 10), rates = seq(1, 1000, length.out=91),
 method = "hunt", d = 100, S = 0.1, Tr = 100, lmda = 100000, lmda_max = 10)

szipper/streamDepletr documentation built on Feb. 2, 2024, 4:10 a.m.