,
Kompetenzzentrum Wasser Berlin gGmbH [cph]| Title: | Calculation of Groundwater Mounding Beneath an Infiltration Basin |
|---|---|
| Description: | Calculation groundwater mounding beneath an infiltration basin based on the Hantush (1967) equation (<doi:10.1029/WR003i001p00227>). The correct implementation is shown with a verification example based on a USGS report (page 25, <https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf#page=35>). |
| Authors: | Michael Rustler [aut, cre] ,
Kompetenzzentrum Wasser Berlin gGmbH [cph] |
| Maintainer: | Michael Rustler <[email protected]> |
| License: | MIT + file LICENSE |
| Version: | 0.3.0 |
| Built: | 2024-10-29 03:16:14 UTC |
| Source: | https://github.com/KWB-R/kwb.hantush |
Hantush equation base properties
baseProperties(time = 10, basinWidth = 10, basinLength = 10, infiltrationRate = 0.5, horizConductivity = 10, iniHead = 10, specificYield = 0.2, numberTimeSteps = 150)baseProperties(time = 10, basinWidth = 10, basinLength = 10, infiltrationRate = 0.5, horizConductivity = 10, iniHead = 10, specificYield = 0.2, numberTimeSteps = 150)
time |
time elapsed since recharge began (T), (Default: 1.5 days) |
basinWidth |
half width of the recharge basin (L), (Default: 10 m) |
basinLength |
half length of the recharge basin (L), (Default: 10 m) |
infiltrationRate |
recharge (infiltration) rate (L/T), (Default: 0.5 m/d) |
horizConductivity |
horizontal hydraulic conductivity (L/T), (Default: 10 m/d) |
iniHead |
initial head (height of the water table above the base of the aquifer);(L), (Default: 10) |
specificYield |
specific yield (Default: 0.2) |
numberTimeSteps |
number of time steps to be used for average aquifer thickness calculation (Default: 150) |
Base properties for Hantush equation
p.22, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
USGS verification example
baseProps_ex1()baseProps_ex1()
p.23, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
USGS verification example: compare R results for water level increase with other models and calculate statistical goodness of fit values (e.g. RMSE, PBIAS, NSE)
compareModelResults(conf = example1())compareModelResults(conf = example1())
conf |
list as retrieved by example1() |
data.frame with R results and other model including all goodness of fit criteria calculated with gof() of package hydrogof
Table 5, p.25, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
example1 for the USGS example parameterisation with distances
Error function
erf(x)erf(x)
x |
x |
Error function result
https://stat.ethz.ch/R-manual/R-devel/library/stats/html/Normal.html
USGS verification example: model parameterisation (multiple distances)
example1(x = c(0, 0.3, 3.3, 6.6, 10, 20, 25, 30, 40, 50, 75, 100, 150, 200), baseProps = baseProps_ex1(), dbg = FALSE)example1(x = c(0, 0.3, 3.3, 6.6, 10, 20, 25, 30, 40, 50, 75, 100, 150, 200), baseProps = baseProps_ex1(), dbg = FALSE)
x |
distance from the center of the recharge basin in the x direction (L) |
baseProps |
basic model properties as retrieved by baseProps_ex1() |
dbg |
If True additional messages on debug messages |
p.23, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
res <- example1() if (FALSE) { #### Head for each time step (defined with parameter "numberTimeSteps) xyplot(head ~ x | as.factor(sprintf("%f days", timeSteps)), data = res$timeSteps, type = "b", las = 1, as.table = TRUE ) #### Head at end of simulation plot(head ~ x, data = res$simTime, type = "b", las = 1 ) } #### Water level increase at end of simulation & compare to alternative models modelComparison <- compareModelResults(conf = res) plotModelComparison(modelComparison = modelComparison)res <- example1() if (FALSE) { #### Head for each time step (defined with parameter "numberTimeSteps) xyplot(head ~ x | as.factor(sprintf("%f days", timeSteps)), data = res$timeSteps, type = "b", las = 1, as.table = TRUE ) #### Head at end of simulation plot(head ~ x, data = res$simTime, type = "b", las = 1 ) } #### Water level increase at end of simulation & compare to alternative models modelComparison <- compareModelResults(conf = res) plotModelComparison(modelComparison = modelComparison)
USGS verification example: get model comparison table
getModelComparisonTable()getModelComparisonTable()
data.frame with water level increase of different model approaches
Table 5, p.25, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
modelComparison <- getModelComparisonTable()modelComparison <- getModelComparisonTable()
Hantush equation
hantush(x = 0, y = 0, baseProps = baseProperties(), dbg = TRUE)hantush(x = 0, y = 0, baseProps = baseProperties(), dbg = TRUE)
x |
distance from the center of the recharge basin in the x direction (L) |
y |
distance from the center of the recharge basin in the y direction (L) |
baseProps |
basic model properties as retrieved by baseProperties() |
dbg |
If True additional messages on integration quality of function hantushSstar are printed on screen |
Head at a given time after recharge begins
p.22, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
baseProperties for basic model properties
Hantush distance: for multiple coordinates
hantushDistances(x = 0:10, y = rep(0, length(x)), baseProps = baseProperties(), dbg = TRUE)hantushDistances(x = 0:10, y = rep(0, length(x)), baseProps = baseProperties(), dbg = TRUE)
x |
vector with distances from the center of the recharge basin in the x direction (L), (Default: each meter between 0-100m) |
y |
vector with distances from the center of the recharge basin in the y direction (L), (Default: 0 times length of x) |
baseProps |
basic model properties as retrieved by baseProperties() |
dbg |
If True additional debug messages are printed on screen |
Head at a given time after recharge begins
hantush for parameterizing the Hantush equation
Hantush distances & base properties: allows input of vector of x,y coordinates and also a vector for one of the base properties
hantushDistancesBaseProps(x = seq(0, 200, 5), y = rep(0, length(x)), baseProps = baseProperties(time = 2^(0:6), infiltrationRate = 1, basinWidth = 10, basinLength = 50, horizConductivity = 10, iniHead = 10, specificYield = 0.2), dbg = FALSE)hantushDistancesBaseProps(x = seq(0, 200, 5), y = rep(0, length(x)), baseProps = baseProperties(time = 2^(0:6), infiltrationRate = 1, basinWidth = 10, basinLength = 50, horizConductivity = 10, iniHead = 10, specificYield = 0.2), dbg = FALSE)
x |
vector with distances from the center of the recharge basin in the x direction (L), (Default: every 5 meter between 0-200m) |
y |
vector with distances from the center of the recharge basin in the y direction (L), (Default: 0 times length of x) |
baseProps |
as retrieved by baseProperties(), but one property is allowed to be a vector (e.g. infiltrationRate = c(1,2,4)) |
dbg |
If True additional debug messages are printed on screen |
List with sublists "dat" (x,y,head & WLincrease), "changedBaseProp.Name" (name of base property with multiple values) and "baseProps" (complete base parameterisation)
baseProperties for basic model properties
baseProps <- baseProperties( time = 2^(0:6), infiltrationRate = 1, basinWidth = 10, basinLength = 50, horizConductivity = 10, iniHead = 10, specificYield = 0.2, numberTimeSteps = 10 ) res <- hantushDistancesBaseProps(baseProps = baseProps) cols <- length(unique(res$dat[[res$changedBaseProp.Name]])) mainTxt <- sprintf("Changed baseProperty: %s", res$changedBaseProp.Name) xyplot(WLincrease ~ x, groups = res$dat[[res$changedBaseProp.Name]], data = res$dat, type = "b", auto.key = list(columns = cols), main = mainTxt )baseProps <- baseProperties( time = 2^(0:6), infiltrationRate = 1, basinWidth = 10, basinLength = 50, horizConductivity = 10, iniHead = 10, specificYield = 0.2, numberTimeSteps = 10 ) res <- hantushDistancesBaseProps(baseProps = baseProps) cols <- length(unique(res$dat[[res$changedBaseProp.Name]])) mainTxt <- sprintf("Changed baseProperty: %s", res$changedBaseProp.Name) xyplot(WLincrease ~ x, groups = res$dat[[res$changedBaseProp.Name]], data = res$dat, type = "b", auto.key = list(columns = cols), main = mainTxt )
Helper function hantushS
hantushS(x, alpha, beta)hantushS(x, alpha, beta)
x |
distance between 0 and half length of recharge basin |
alpha |
alpha |
beta |
beta |
Hantush star
p.22, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
Hantush function Sstar
hantushSstar(alpha, beta, dbg)hantushSstar(alpha, beta, dbg)
alpha |
alpha |
beta |
beta |
dbg |
If True additional messages on integration quality of function hantushSstar are printed on screen |
Hantush Sstar result
p.22, https://pubs.usgs.gov/sir/2010/5102/support/sir2010-5102.pdf
USGS verification example: plot model comparison results
plotModelComparison(modelComparison = compareModelResults(), title = "", ...)plotModelComparison(modelComparison = compareModelResults(), title = "", ...)
modelComparison |
data.frame as retrieved by compareModelResults(), Default: (compareModelResults()) |
title |
to be used as title above plot (Default: "") |
... |
further arguments passed to function xyplot() |
model comparison goodness of fit plots
compareModelResults for comparison with USGS benchmark models
### Plot model comparison with title "Model comparison" and one plot for ### each page plotModelComparison( title = "Model comparison", layout = c(1, 1) )### Plot model comparison with title "Model comparison" and one plot for ### each page plotModelComparison( title = "Model comparison", layout = c(1, 1) )