aquiferRouting {WRSS} | R Documentation |
base function for aquifer simulation
Description
Given a sort of demand(s), aquiferRouting
function simulates a lumped and simple model of an unconfined aquifer under an optional givn recharge time series, rechargeTS
, and specific yield, Sy
.
Usage
aquiferRouting(demand, priority = NA, area, volume,
rechargeTS = NA, leakageFraction = NA,
initialStorage = NA, Sy, simulation)
Arguments
demand |
(optional) A matrix: is column-wise matrix of demands, at which the rows present demands for each monthly time step and columns are for different individual demand sites (MCM). |
priority |
(optional) A vector: is a vector of priorities associated to |
area |
The area of aquifer (Km^2) |
volume |
The aquifer volume (MCM) |
rechargeTS |
(optional) A vector : a vector of water flowing into the aquifer (MCM) |
leakageFraction |
(optional) The leakage coeffcient of aquifer storage. The leakage is computed as the product of |
initialStorage |
(optional) The initial volume of aquifer at the first step of the simulation (MCM). If missing, the function iterates to carry over the aquifer |
Sy |
Specific yield (default: 0.1) |
simulation |
A list: |
Value
the aquiferRouting
function returns a list of objects as bellow:
-
release
: a matrix of release(s) equivalant to eachdemand
(MCM) -
leakage
: a vector of leakage time series (MCM) -
storage
: a vector of storage time series (MCM)
Author(s)
Rezgar Arabzadeh
References
Mart nez-Santos, P., and J. M. Andreu. "Lumped and distributed approaches to model natural recharge in semiarid karst aquifers." Journal of hydrology 388.3 (2010): 389-398.
See Also
Examples
area <-200
leakageFraction<-0.01
Sy <-0.15
volume <-20000
priority <-c(3,1,1,2)
rechargeTS <-rnorm(120,60,8)
demand <-matrix(rnorm(480,10,3),120)
simulation <-list(start='2000-01-01',end='2009-12-29',interval='month')
res<-
aquiferRouting(demand =demand ,
priority =priority ,
area =area ,
volume =volume ,
rechargeTS =rechargeTS ,
leakageFraction=leakageFraction,
Sy =Sy ,
simulation =simulation)
plot(res$storage$storage,ylab='Storage (MCM)',xlab='time steps(month)',type='o')