Date of Award

Summer 1996

Document Type


Degree Name

Doctor of Philosophy (PhD)


Civil/Environmental Engineering

Committee Director

A. Osman Akan

Committee Member

David R. Basco

Committee Member

Resit Unal

Committee Member

Dayanand N. Naik


A comprehensive conjunctive-use management model is developed. The dynamics of flow and solute transport processes in connected surface water groundwater systems are integrated by a dual programming management model. The governing aquifer flow parameters and streamflows are treated as stochastic random processes. Multiple realizations of the random field are generated and are explicitly incorporated in a non-linear optimization model along with other system, environmental, and management constraints. To facilitate management of large aquifer systems, a linked simulation-optimization approach is used. The simulation program generates the response matrices for flow and transport processes. The management model then determines optimal well discharges and optimal surface water diversion rates. Further, the model determines optimal concentration injection rates for recharge wells and optimal concentration disposal into surface water bodies.

Implicit finite difference method is used in modeling the two-dimensional, unsteady groundwater flow and transport process. Iterative alternating direction implicit method is used in its solution. Leaky aquifer, evapotranspiration, aerial recharge, and induced infiltration can be accounted for. Advection, diffusion, and dispersion of conservative and non-conservative substances are considered in the transport model. In modeling advective component four schemes were investigated. It was found that the quadratic upstream interpolation method is the best.

Implicit finite difference method is used in modeling the one-dimensional, unsteady surface water flow and solute transport processes. The surface water groundwater interaction, stream mass balance, initial and boundary conditions are input to the management model as system constraints.

The individual components of the model developed are tested by comparing the numerical results obtained with analytical solutions and other generic numerical models where available. The model developed is demonstrated for a large aquifer. The considered hypothetical aquifer in model application is based on Yorktown-Eastover aquifer characteristics of Southeastern Virginia Aquifer system.

The comprehensive model developed in this study can be used for field applications of large aquifer systems over a longer management horizon.