Variably saturated groundwater modelling for optimizing managed aquifer recharge using trench infiltration

Spreading‐basin methods have resulted in more than 130 million cubic metres of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and (or) clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water‐table depth, alternate wet/dry periods, and number of parallel trenches. Modelling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water‐table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings. Copyright © 2014 John Wiley & Sons, Ltd.     

A new adiabatic kernel for the MC2 model

Abstract

Traditional semi‐implicit formulations of nonhydrostatic compressible models may not be stable in the presence of steep terrain when pressure gradient terms are split and lagged in time. If all pressure gradient terms and the divergence are treated implicitly, the resulting wave equation for the pressure contains off‐diagonal cross‐derivative terms leading to a highly nonsymmetric linear system of equations. In this paper we present a more implicit formulation of the Mesoscale Compressible Community (MC2) model employing a Generalized Minimal Residual (GMRES) Krylov iterative solver and a more efficient semi‐Lagrangian advection scheme. Open boundaries now permit exact upwind interpolation and the ability to reproduce simulations to machine precision is illustrated for one‐way nesting at equivalent resolution. Numerical simulations of hydrostatic and nonhydrostatic mountain waves demonstrate the stability and accuracy of the new adiabatic kernel. The computational efficiency of the model is reported for 1D Jacobi and 3D Alternating Direction Implicit (ADI) line relaxation preconditioned implemented with a parallel data transposition strategy.