intermittent_pumping {streamDepletr}

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)