  Groundwater flow in Washington DC greatly influences the surface water
quality in urban areas. The current methods of flow estimation, based on
Darcy's Law and the groundwater flow equation, can be described by the
diffusion equation (the transient flow) and the Laplace equation (the
steady-state flow). The Laplace equation is a simplification of the diffusion
equation under the condition that the aquifer has a recharging boundary. The
practical way of calculation is to use numerical methods to solve these
equations. The most popular system is called MODFLOW, which was developed by
USGS. MODFLOW is based on the finite-difference method in rectangular Cartesian
coordinates. MODFLOW can be viewed as a "quasi 3D" simulation since it only
deals with the vertical average (no z-direction derivative). Flow calculations
between the 2D horizontal layers use the concept of leakage. In this project,
we have established a mathematical model based on gradually varied functions
for groundwater data volume reconstruction. These functions do not rely on the
rectangular Cartesian coordinate system. A gradually varied function can be
defined in a general graph or network. Gradually varied functions are suitable
for arbitrarily shaped aquifers. Two types of models are designed and
implemented for real data processing: (1) the gradually varied model for
individual (time) groundwater flow data, (2) the gradually varied model for
sequential (time) groundwater flow data. In application, we also established a
MySQL database to support the related research. The advantage of the gradually
varied fitting and its related method does not need the strictly defined
boundary condition as it is required in MODFLOW.
