pastas.recharge.Peterson

class pastas.recharge.Peterson

Recharge to the groundwater calculated based on Peterson and Western [2014].

The water balance for the unsaturated zone reservoir is written as:

\[\frac{dS}{dt} = P_e - E_a - R\]

where the fluxes $P_e$, $E_a$ and $R$ are calculated as:

\[P_e = P \left(1 - \frac{S}{\hat{S_{cap}}}\right)^\alpha\]

\[E_a = E_p \left(\frac{S}{\hat{S_{cap}}}\right)^\gamma\]

\[R = \hat{k_{sat}}\left(\frac{S}{\hat{S_{cap}}}\right)^{\hat{\beta}}\]

with the parameters:

\[\hat{S_{cap}} = 10^{S_{cap}}; \hat{k_{sat}} = 10^{k_{sat}}; \hat{\beta} = 10^{\beta}\]

Note that the method currently uses forward Euler method to solve the ODE so significant water balance errors can occur.

Methods

get_init_parameters(→ pandas.DataFrame)	Method to obtain the initial parameters.
simulate(...)	Simulate the recharge flux.
get_recharge(→ tuple[pastas.typing.ArrayLike, ...)	Internal method used for the recharge calculation.