Aquifer Modeling by Numerical Methods Applied to an Arizona Groundwater Basin
KeywordsAquifers -- Arizona -- Mathematical models.
Groundwater -- Arizona -- Mathematical models.
Aquifers -- Mathematical models.
Groundwater -- Mathematical models.
MetadataShow full item record
RightsCopyright © Arizona Board of Regents
Collection InformationThis title from the Hydrology & Water Resources Technical Reports collection is made available by the Department of Hydrology & Atmospheric Sciences and the University Libraries, University of Arizona. If you have questions about titles in this collection, please contact firstname.lastname@example.org.
AbstractFLUMP, a recently developed mixed explicit -implicit finite -element program, was calibrated against a data base obtained from a portion of the Tucson Basin aquifer, Arizona, and represents its first application to a real -world problem. Two previous models for the same region were constructed (an electric analog and a finite -difference model) in which calibration was based on prescribed flux boundary conditions along stream courses and mountain fronts. These fluxes are not directly measured and estimates are subject to large uncertainties. In contrast, boundary conditions used in the calibration of FLUMP were prescribed hydraulic heads obtained from direct measurement. At prescribed head boundaries FLUMP computed time - varying fluxes representing subsurface lateral flow and recharge along streams. FLUMP correctly calculated fluctuations in recharge along the Santa Cruz River due to fluctuations in storm runoff and sewage effluent release rates. FLUMP also provided valuable insight into distributions of recharge, discharge, and subsurface flow in the study area.Properties of FLUMP were compared with those of two other programs in current use: ISOQUAD, a finite -element program developed by Pinder and Frind (1972), and a finite- difference program developed by the U.S. Geological Survey (Trescott, et al., 1976). It appears that FLUMP can handle a larger class of problems than the other two programs, including those in which the boundary conditions and aquifer parameters vary arbitrarily with time and /or head. FLUMP also has the ability to solve explicitly when accuracy requires small time steps, the ability to solve explicitely in certain parts of the flow region while solving implicitly in other parts, flexibility in mesh design and numbering of nodes, computation of internal as well as external fluxes, and global as well as local mass balance checks at each time step.
Series/Report no.Technical Reports on Natural Resource Systems, No. 32
Showing items related by title, author, creator and subject.
Modeling System Reliability For Digital Preservation: Model Modification and Four-Copy Model StudyHan, Yan; Chan, Chi Pak (British Library, 2008)Research has been studied to evaluate the reliability of storage media and the reliability of a computer backup system. In this paper, we use the Continuous Time Markov Chain to model and analyze the reliability of a computer backup system. We propose a modified model from that of the Constantopoulos, Doerr and Petraki . We analyze the difference, show computational results, and propose new input parameters (e.g. time to repair) for the model from our experience. Further we developed a four-copy data model to test if it fulfills the sample reliability rate set by the RLG-NARA. The modeling process can be applied to construct models for computer preservation systems using different storage media. The reliability of constructed models can be calculated so that preservation institutions can have quantitative data to decide their preservation strategies.
Dendroclimatic Analysis Using Thornwaite-Mather-Type Evapotranspiration Models: A Bridge Between Dendroevology and Forest Simulation ModelsLeBlanc, David; Terrell, Mark; Department of Biology, Ball State University, Muncie, IN (Tree-Ring Society, 2001)The objective of this study was to document correlations between radial growth of white oak (Quercus alba L.) at 128 sites in the eastern US and variables related to early growing season site water balance, including the ratio of actual to potential evapotranspiration (AE/PE) computed based on the procedure described by Thornthwaite and Mather (1957). White oak radial growth was strongly correlated with all measures of early growing season water balance, but was most consistently and strongly correlated with Palmer Drought Severity Index (PDI) and AE /PE ratio computed using a modified Thornthwaite-Mather model. We propose that dendroecological analyses of tree growth responses to climate based on AE /PE variables could provide empirical data useful for improving climate response algorithms in forest simulation models. This change of standard practice could also improve biological interpretations derived from such dendroecological analyses.
The effects of surface roughness and stress on lattice gas models using kinetic Monte Carlo modelingJackson, Kenneth A.; Bentz, Daniel N. (The University of Arizona., 2003)This thesis reports kinetic Monte Carlo computer simulations using a lattice gas model conducted on a variety of systems. These studies may be divided into two main categories: rod eutectics and related surface roughening, and surface morphology changes due to local stresses. The first grouping is a study of irregular rod eutectic systems. Simulations of directional solidification of rod eutectic systems were conducted using a model similar to the spin one Ising model. Growth of the rods was initiated from columns of pure A atoms embedded in a matrix of B atoms. The growth characteristics of the eutectic depend on the location of the surface roughening transition for the two phases. The surface roughening transition was determined using fluctuation dissipation theory, where the kinetic behavior of the interface is related to a characteristic time of fluctuations about an equilibrium position. These times were determined by time correlations. Results show a sharp transition in the kinetic behavior of the interface as a function of Jackson's alpha factor. This is the first time this method has been used to locate the surface roughening transition. An applied temperature gradient supplied the restoring force for the interface. The roughening behavior of binary alloys was also examined and compared to pure component systems. The second set of simulations reported here examine the effects of local stress on surface morphology. The weakening of bonds due to the dislocation stress field was studied as the origin of the formation of etch-pits at dislocations. Atoms from a diamond cubic lattice were irreversibly removed with a probability which depends on an local surface configuration as well as on the local stress developed from its physical location with respect to a dislocation in the lattice. In accordance with experimental observations, both faceted and non-faceted dislocation etch-pits have been observed. Simulations of crystal surfaces near equilibrium have reproduced direct experimental results using atomic force microscopy (AFM). The probe tip interacts with the shape and the motion of step edges, and the motion of a step is retarded in the vicinity of the tip.