Non-Steady Mixing Cell Model for Assessing Ground Water Fluxes and Storage Capacity in Complex Aquifers

Monday, April 20, 2009: 10:50 a.m.
Coronado I (Hilton Tucson El Conquistador Golf & Tennis Resort )
Eilon M. Adar , Zuckerberg Inst. for Water Research, Ben Gurion University of the Negev, Sede Boqer Campus, Israel
Nimrod Halamish , Zuckerberg Inst. for Water Research, Ben Gurion University of the Negev, Sede Boqer Campus, Israel
Shaul Sorek , Zuckerberg Inst. for Water Research, Ben Gurion University of the Negev, Sede Boqer Campus, Israel
This paper presents a non-steady mixing cell approach that enables quantitative assessment of groundwater fluxes and storage capacity in a complex aquifer and transient hydro-geological systems. The algorithm proposed in this study is based on a multi-compartmental flow model incorporated with environmental and anthropogenic dissolved minerals in a transient hydrological system.
The flow model was built according to the principle of division into mixing cells known as “multi-compartmental model”. The sub-division of the aquifers into discrete sub-hydrological systems for every selected time period creates hydrological units with relatively homogeneous and fixed values of every hydrological parameter without gradients of hydraulic heads or solutes within the compartments. The model assumes mixing of solutes in each compartment due to water contribution from various sources with differing solute compositions resulting a unique hydrochemical composition in each cell for specific time.
All potential sources of recharge into each cell must be identified and characterized by its hydro-chemical composition. The flow pattern is constructed such that any hydrological connection between cells will have a unique flow direction. Fluxes into and out of cells are proposed along the downstream hydraulic head gradients.
The mathematical algorithm is based on set of balance equations for the flux of water and for the associated fluxes of solutes for each cell in the modelled aquifer over a given time period. The Least Squares mathematical algorithm E04NCF (by NAG®) is proposed for solving the water and solutes balance expressions in a multi-compartmental flow model with spatial and temporal variations of groundwater fluxes. It is capable of accommodating several orders of magnitude differences among the unknown parameters while limiting the solution to only positive values of water fluxes.
The paper demonstrates the use of the mathematical algorithm for a quantitative assessment of groundwater fluxes and recharge in the Central Arava Valley, Israel.