Transport of pollutants in aquifer largely depends on groundwater flow which is governed by aquifer hydraulic parameters. Determination of these parameters and associated groundwater modeling become essential for adequate remediation of contaminate groundwater. The objective of this paper is to analyze groundwater flow and determine the optimum hydraulic parameters by performing groundwater modeling based on sensitivity analysis for unconfined sandy gavel aquifer constructed in a laboratory scale under various boundary condition. Results revealed that the simulated drawdown was lower than the observed drawdown irrespective of boundary conditions. and specific yield (S$_{y}$) had less effect on the grondwater flow than permeability (K) in the aquifer. Water balance analysis showed that the measured drawdown in neighboring observation wells during pumping was higher than either simulated or recovered water table. The indicated that a difference might exist in the water tables between aquifer and wells. The difference was investigated by time domain reflectometry (TDR) measurements on water contents in the region of water table and capillary fringe, and explained by a delayed response of water table during gravitational drainage as the water table was lowered as a result of pumping.g.