streamDepletr (version 0.1.1)

intermittent_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, ...)

Value

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

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.

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)