Modeling multiaquifer wells with MODFLOW.

Christopher J Neville; Matthew J Tonkin

doi:10.1111/j.1745-6584.2004.t01-9-.x

Home/Archive/PubMed

PubMed · 15584304

Modeling multiaquifer wells with MODFLOW.

Christopher J Neville·Matthew J Tonkin

View original source ↗Export metadata ↓

Abstract

Multiaquifer wells, i.e., wells that are open across more than one aquifer, can have a profound effect on the hydraulics of a ground water system. These wells change the physical system by establishing direct hydraulic links between otherwise isolated strata. Several methods are available to simulate multiaquifer wells in the context of comprehensive ground water flow simulators. In this paper, we review four methods to represent multiaquifer wells with the widely used code MODFLOW. These methods include a specialized code developed, but never formally released, by the U.S. Geological Survey (USGS), the Multi-Aquifer Well (MAW1) Package. An expanded implementation of the techniques in the MAW1 Package has been incorporated in the Multi-Node Well Package released recently by the USGS (Halford and Hanson 2002). We examine the performance of the methods in the context of a benchmarking study against the analytical solutions of Papadopulos (1966) and Sokol (1963). Our results demonstrate that results obtained with the MAW1 Package closely match exact solutions for pumping and nonpumping conditions, using both coarse and refined grids.

Explore related subjects

Geological Phenomena Geology Models, Theoretical Pressure Water Movements Water Supply Journal Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Christopher J Neville, Matthew J Tonkin. Modeling multiaquifer wells with MODFLOW.. https://doi.org/10.1111/j.1745-6584.2004.t01-9-.x

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2006

Applied tracer tests in fractured rock: Can we predict natural gradient solute transport more accurately than fracture and matrix parameters?

Applied tracer tests provide a means to estimate aquifer parameters in fractured rock. The traditional approach to analysing these tests has been using a single fracture model to find the parameter values that generate the best fit to the measured breakthrough curve. In many cases, the ultimate aim is to predict solute transport under the natural gradient. Usually, no confidence limits are placed on parameter values and the impact of parameter errors on predictions of solute transport is not discussed. The assumption inherent in this approach is that the parameters determined under forced conditions will enable prediction of solute transport under the natural gradient. This paper considers the parameter and prediction uncertainty that might arise from analysis of breakthrough curves obtained from forced gradient applied tracer tests. By adding noise to an exact solution for transport in a single fracture in a porous matrix we create multiple realisations of an initial breakthrough curve. A least squares fitting routine is used to obtain a fit to each realisation, yielding a range of parameter values rather than a single set of absolute values. The suite of parameters is then used to make predictions of solute transport under lower hydraulic gradients and the uncertainty of estimated parameters and subsequent predictions of solute transport is compared. The results of this study show that predictions of breakthrough curve characteristics (first inflection point time, peak arrival time and peak concentration) for groundwater flow speeds with orders of magnitude smaller than that at which a test is conducted can sometimes be determined even more accurately than the fracture and matrix parameters.

Geological Phenomena↗

PubMed2006

Integration of field measurements and reactive transport modelling to evaluate contaminant transport at a sulfide mine tailings impoundment.

Over a decade of field observations including geochemical, mineralogical and hydrological information are available on the generation of acid mine drainage from the Pistol Dam region of the P-area of Inco's tailings impoundment in Copper Cliff, Ontario. This work focuses on the integration and quantitative assessment of this data set using reactive transport modeling. The results of the reactive transport simulations are in general agreement with the field observations; however, exact agreement between the field and simulated results was not the objective of this study, and was not attained. Many factors contribute to the discrepancies between the field observations and simulation results including geochemical and hydrogeological complexities and necessary model simplifications. For example, fluctuating water levels observed at the site were averaged and described using a steady state flow system. In addition, the lack of representative thermodynamic and rate expression data contributed to the discrepancies between observations and simulation results, thus further research into the applicability of laboratory-derived thermodynamic and rate expression data to field conditions could minimize these discrepancies. Despite the discrepancies between the field observations and simulated results, integrating field observations with numerical modelling of the P-area tailings impoundment allowed for a more complete understanding of what affects the complex geochemical reactions.

Geological Phenomena↗

PubMed2006

Geochemical and biological controls on trace metal transport in an acid mine impacted watershed.

Water samples collected in an acid mine impacted watershed indicated that the concentrations of dissolved trace metals were diurnally influenced by mineral saturation, which is controlled primarily by pH and water temperature. Measurements taken suggested that these variations only occur at sample locations immediately downstream from the confluence of acidic and alkaline waters. It is at these locations where initial mineral precipitation occurred and where subtle changes in solubility were most affected, increasing trace metal removal when both the rate of photosynthesis (influencing pH in headwaters) and water temperature were at a maximum. The role of iron photoreduction (increased midday production of ferrous iron) on overall Cu, Mn, and Zn transport was also evaluated, but found to be inconclusive. Iron photoreduction may however influence adsorption and/or coprecipitation of trace metals through associated changes in oxidation state, solubility, and mineralogy of various iron colloids, which are produced upon the neutralization of acidic, metal enriched water. Furthermore, measured values of copper and zinc were compared to relative USEPA chronic criterion for exposure to continuous concentration (CCC) of metals by the calculation of a "toxicity unit" (TU). It was found that average values of both copper and zinc only exceeded the CCC (TU>1) in the acid mine-impacted Leona Creek. In general, zinc toxicity decreased while copper toxicity increased downstream of the confluence of the mine impacted Leona Creek and background Lion Creek (sampled at Lake Aliso), indicating a significant source of zinc in upstream, non mine-impacted samples.

Geological Phenomena↗