Modeling of Contaminant Movement Near Pumping Wells: Saturated-Unsaturated Flow with Particle Tracking

A transient axisymmetric saturated-unsaturated numerical flow model was coupled with a particle tracking model to investigate the movement of contaminants when a shallow unconfined aquifer is pumped at a constant rate. The particle tracking model keeps track of locations and masses of solutes in the aquifer, and the time of capture by the well. At the end of each time-step the flow model solves the Richard's equation for the hydraulic head distribution from which elemental velocities are calculated. Solutes are then displaced for a period equivalent to the time-step using both the magnitude and direction of the elemental velocities. Numerical experiments were performed to investigate effluent concentrations in wells with screens of different length and in different positions relative to zones of stratified contamination. At early times of pumping the effluent concentrations were similar to the concentrations adjacent to the well screen, but at late times, the concentrations approached the vertically averaged concentration in the aquifer. Time to attain the vertically averaged concentration was determined by the well geometry, initial location of the contaminant plume in relation to the well screen, and hydraulic properties of the aquifer. The results are consistent with the hydraulics of flow to a pumping well and of particular importance, they demonstrate that short-term pump tests could give erroneous design concentrations for pump-and-treat systems. The model provides a means of quantifying arrival times and mixing ratios. It could therefore provide a useful means of designing production wells in aquifers with stratified contamination and more efficient recovery systems for aquifer remediation.     

A Model of Regional Ground-Water Flow in Secondary-Permeability Terrane

The ground-water flow system in the Lower Susquehanna River Basin in Pennsylvania and Maryland can be considered as one complex unconfined aquifer in which secondary porosity and permeability are the dominant influences on the occurrence and flow of ground water. The degree of development of secondary porosity and permeability in the various lithologies of the lower basin determines the aquifer characteristics of each lithology. Based on qualitative evidence, the use of a porous-media model was assumed to be appropriate on a regional scale and a finite-difference ground-water flow model was constructed for the lower basin. The conceptual model of ground-water flow in the lower basin incorporates the major features of the flow system. Through the use of two layers, 21 hydrogeologic units, and five topographic settings, the conceptual model was systematically reduced to arrive at a simplified conceptual model. Further reduction produced a numerical model representation of the conceptual model, in which the essential features of the lower-basin flow system were quantified for input into the numerical model. The model was calibrated under both steady-state and transient conditions, and was used to evaluate the water-supply potential of the 21 hydrogeologic units. The carbonate units have the greatest potential for ground-water development and the Triassic sedimentary and crystalline units have the least potential. A total ground-water yield potential of about 900 million gallons per day could be obtained from the lower basin with a consequent 50-percent reduction of base flow in streams.     

Modeling Transport in Transient Ground-Water Flow: An Unacknowledged Approximation

Abstract. During unsteady or transient ground-water flow, the fluid mass per unit volume of aquifer changes as the potentiometric head changes, and solute transport is affected by this change in fluid storage. Three widely applied numerical models of two-dimensional transport partially account for the effects of transient flow by removing terms corresponding to the fluid continuity equation from the transport equation, resulting in a simpler governing equation. However, fluid-storage terms remaining in the transport equation that change during transient flow are, in certain cases, held constant in time in these models. For the case of increasing heads, this approximation, which is unacknowledged in these models'documentation, leads to transport velocities that are too high, and increased concentration at fluid and solute sources. If heads are dropping in time, computed transport velocities are too low. Using parameters that somewhat exaggerate the effects of this approximation, an example numerical simulation indicates solute travel time error of about 14 percent but only minor errors due to incorrect dilution volume. For horizontal flow and transport models that assume fluid density is constant, the product of porosity and aquifer thickness changes in time: initial porosity times initial thickness plus the change in head times the storage coefficient. This formula reduces to the saturated thickness in unconfined aquifers if porosity is assumed to be constant and equal to specific yield. The computational cost of this more accurate representation is insignificant and is easily incorporated in numerical models of solute transport.