Glover model of stream depletion fraction

This function estimates stream depletion at time `t` as a fraction of pumping from an individual pumping well. See Glover (1954).

The function requires variables x1, K, D, V, t. These variables can be specified as columns of df, or as named variables in the function call. If df is specified, the named variables are ignored.

get_stream_depletion_fraction(
  df,
  x1 = NULL,
  K = NULL,
  D = NULL,
  V = NULL,
  t = NULL
)

Arguments

df

data.frame with columns specifying all parameters

x1

Distance between well and river

K

Saturated hydraulic conductivity

D

Depth of aquifer

V

Drainable porosity of aquifer

t

Time from pumping onset at which to calculate stream depletion fraction

Details

Glover model of stream depletion, including image well

Examples

# Reproduce example from Glover
library(units)
x1 <- set_units(c(1, 5, 10) * 1e3, "ft")
D <- set_units(100, "ft")
K <- set_units(0.001, "ft/sec")
t <- set_units(5, "year")
V <- 0.2 # unitless

# Specifying parameters as numeric or vector inputs
stream_depletion_fraction <- get_stream_depletion_fraction(x1 = x1, K = K, D = D, V = V, t = t)
stream_depletion_fraction
#> [1] 0.9365474 0.6905933 0.4259739

# Specifying parameters as named data.frame columns
library(tibble) # simplifies specifying data.frames with units objects
df <- tibble(x1 = x1, K = K, D = D, V = V, t = t)
stream_depletion_fraction <- get_stream_depletion_fraction(df)
stream_depletion_fraction
#> [1] 0.9365474 0.6905933 0.4259739