Search USGS⌕ Search

SEARCH · Search USGS

Results for “Ground Water Modeling Investigation Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources
USGS2005

Hydrogeology and trichloroethene contamination in the sea-level aquifer beneath the Logistics Center, Fort Lewis, Washington

The U.S. Army disposed of waste trichloroethene (TCE) and other materials in the East Gate Disposal Yard near the Logistics Center on Fort Lewis, Washington, from the 1940s to the early 1970s. As a result, ground water contaminated with primarily TCE extends more than 3 miles downgradient from the East Gate Disposal Yard. The site is underlain by a complex and heterogeneous sequence of glacial and non-glacial deposits that have been broadly categorized into an upper and a lower aquifer (the latter referred to as the sea-level aquifer). TCE contamination was detected in both aquifers. This report describes an investigation by the U.S. Geological Survey (USGS) of the source, migration, and attenuation of TCE in the sea-level aquifer. A refined conceptual model for ground-water flow and contaminant migration into and through the sea-level aquifer was developed in large part from interpretation of environmental tracer data. The tracers used included stable isotopes of oxygen (18O), hydrogen (2H), and carbon (13C); the radioactive hydrogen isotope tritium (3H); common ions and redox-related analytes; chlorofluorocarbons; and sulfur hexafluoride. Tracer and TCE concentrations were determined for samples collected by the USGS from 37 wells and two surface-water sites in American Lake during 1999-2000. Ground-water levels were measured by the USGS in more than 40 wells during 2000-01, and were combined with measurements by the U.S. Army and others to create potentiometric-surface maps. Localized ground-water flow features were identified that are of particular relevance to the migration of TCE in the study area. A ridge of ground water beneath American Lake diverts the flow of TCE-contaminated ground water in the sea-level aquifer to the west around the southern end of the lake. Tracer data provided clear evidence that American Lake is a significant source of recharge to the sea-level aquifer that has created that ridge of ground water. High ground-water altitudes at locations north and northeast of the Logistics Center combined with the ridge beneath American Lake prevent TCE contaminated water beneath the Logistics Center from migrating toward municipal water-supply wells northeast of the site. The 1999-2000 TCE concentrations measured by the USGS at older wells screened in the sea-level aquifer were similar to those measured since 1995, but the known downgradient extent of the TCE contamination expanded nearly 2 miles after the Army installed and sampled new wells during 2003-04. Concentrations of TCE in the sea-level aquifer were consistently highest in the upper part of the aquifer throughout the plume, although TCE has spread throughout much of the thickness of the aquifer in the downgradient portions of the plume. Environmental tracer data indicated that the primary pathway for contaminant migration into the sea-level aquifer is through the previously identified confining unit window, an area where the predominately fine-grained confining unit is relatively coarse grained and more permeable. Other less substantial pathways for contaminant migration also were identified near the East Gate Disposal Yard and the I-5 pump-and-treat facilities. Those areas are near active pumping wells and ground-water reintroduction facilities, but there is no evidence that the contaminant migration was caused or enhanced by those activities. Within the sea-level aquifer, TCE concentrations continue to migrate westward in the flow field strongly influenced by ground-water recharge from American Lake. Historical data are not available to definitively determine if the 5-?g/L leading edge of the current TCE plume is stable or if it is still moving downgradient. However, an evaluation of the available data combined with TCE traveltime estimates indicates that the peak TCE concentrations in the sea-level aquifer may have not yet reached the wells near the currently defined leading edge of the plume. Hypothetically, the 5-?g/L leading edge

Scientific Investigations Report↗
USGS2004

Documentation of the Santa Clara Valley regional ground-water/surface-water flow model, Santa Clara Valley, California

The Santa Clara Valley is a long, narrow trough extending about 35 miles southeast from the southern end of San Francisco Bay where the regional alluvial-aquifer system has been a major source of water. Intensive agricultural and urban development throughout the 20th century and related ground-water development resulted in ground-water-level declines of more than 200 feet and land subsidence of as much as 12.7 feet between the early 1900s and the mid-1960s. Since the 1960s, Santa Clara Valley Water District has imported surface water to meet growing demands and reduce dependence on ground-water supplies. This importation of water has resulted in a sustained recovery of the ground-water flow system. To help support effective management of the ground-water resources, a regional ground-water/surface-water flow model was developed. This model simulates the flow of ground water and surface water, changes in ground-water storage, and related effects such as land subsidence. A numerical ground-water/surface-water flow model of the Santa Clara Valley subbasin of the Santa Clara Valley was developed as part of a cooperative investigation with the Santa Clara Valley Water District. The model better defines the geohydrologic framework of the regional flow system and better delineates the supply and demand components that affect the inflows to and outflows from the regional ground-water flow system. Development of the model includes revisions to the previous ground-water flow model that upgraded the temporal and spatial discretization, added source-specific inflows and outflows, simulated additional flow features such as land subsidence and multi-aquifer wellbore flow, and extended the period of simulation through September 1999. The transient-state model was calibrated to historical surface-water and ground-water data for the period 1970–99 and to historical subsidence for the period 1983–99. The regional ground-water flow system consists of multiple aquifers that are grouped into upper- and lower-aquifer systems. Ground-water inflow occurs as natural recharge in the form of streamflow infiltration and areal infiltration of precipitation along stream channels, artificial recharge from infiltration of imported water at recharge ponds and along selected stream channels, and leakage along selected transmission pipelines. Ground-water outflow occurs as evapotranspiration, stream base flow, discharge through pumpage from wells, and subsurface flow to the San Francisco Bay. The geohydrologic framework of the regional ground-water flow system was represented as six model layers. The hydraulic properties were redefined on the basis of cell-based lithologic properties that were delineated in terms of aggregate thicknesses of coarse-grained, fine-grained, and mixed textural categories. The regional aquifer systems also are dissected by several laterally extensive faults that may form at least partial barriers to the lateral flow of ground water. The spatial extent of the ground-water flow model was extended and refined to cover the entire Santa Clara Valley, including the Evergreen subregion. The temporal discretization was refined and the period of simulation was extended to 1970–99. The model was upgraded to MODFLOW-2000 (MF2K) and was calibrated to fit historical ground-water levels, streamflow, and land subsidence for the period 1970–99. The revised model slightly overestimates measured water levels with an root-mean-square error of -7.34 feet. The streamflow generally shows a good match on gaged creeks and rivers for flows greater than 1.2 cubic feet per second. The revised model also fits the measured deformation at the borehole extensometer site located near San Jose within 16 to 27 percent and the extensometer site near Sunnyvale within 3 percent of the maximum measured seasonal deformation for the deepest extensometers. The total ground-water inflow and outflow of about 225,500 acre-feet per year (acre-ft/yr) for the period 1970–89 and of about 205,300 acre-feet per year for the period for the period 1970–99 is comparable with that of the previous model, 207,200 acre-ft/yr for the period 1970–89. Overall the simulated net change in storage increased by about 189,500 acre-ft/yr for the entire period of simulation, which represents about one and a half years of the 1970–99 average pumping. The changes in ground-water flow and storage generally reflect the major climate cycles and the additional importation of water by Santa Clara Valley Water District, with the basin in recovery since the drought of the late 1980s and early 1990s. The average total recharge rate, from natural and artificial recharge and from streamflow infiltration for the revised model for the entire simulation period 1970–99, was about 157,100 acre-ft/yr, which represents about 59 percent of the inflow to the ground-water flow system. The average rate of artificial recharge of about 77,600 acre-ft/yr represents about 30 percent of the inflow to the ground-water flow system. The average pumpage for the entire 29.75-year simulation period is about 133,400 acre-ft/yr and represents about 69 percent of the outflow from the ground-water flow system. Most of the simulated recharge infiltrates and flows through the uppermost layers (i.e. model layers 1 and 3) of the aquifer system. Most of the water that flows to the deeper model layers is occurring through wellbores, with wellbore flow representing 19 percent of the total ground-water inflow between model layers.

California↗
USGS2010

Interpretation of Flow Logs from Nevada Test Site Boreholes to Estimate Hydraulic Conductivity Using Numerical Simulations Constrained by Single-Well Aquifer Tests

Hydraulic conductivities of volcanic and carbonate lithologic units at the Nevada Test Site were estimated from flow logs and aquifer-test data. Borehole flow and drawdown were integrated and interpreted using a radial, axisymmetric flow model, AnalyzeHOLE. This integrated approach is used because complex well completions and heterogeneous aquifers and confining units produce vertical flow in the annular space and aquifers adjacent to the wellbore. AnalyzeHOLE simulates vertical flow, in addition to horizontal flow, which accounts for converging flow toward screen ends and diverging flow toward transmissive intervals. Simulated aquifers and confining units uniformly are subdivided by depth into intervals in which the hydraulic conductivity is estimated with the Parameter ESTimation (PEST) software. Between 50 and 150 hydraulic-conductivity parameters were estimated by minimizing weighted differences between simulated and measured flow and drawdown. Transmissivity estimates from single-well or multiple-well aquifer tests were used to constrain estimates of hydraulic conductivity. The distribution of hydraulic conductivity within each lithology had a minimum variance because estimates were constrained with Tikhonov regularization. AnalyzeHOLE simulated hydraulic-conductivity estimates for lithologic units across screened and cased intervals are as much as 100 times less than those estimated using proportional flow-log analyses applied across screened intervals only. Smaller estimates of hydraulic conductivity for individual lithologic units are simulated because sections of the unit behind cased intervals of the wellbore are not assumed to be impermeable, and therefore, can contribute flow to the wellbore. Simulated hydraulic-conductivity estimates vary by more than three orders of magnitude across a lithologic unit, indicating a high degree of heterogeneity in volcanic and carbonate-rock units. The higher water transmitting potential of carbonate-rock units relative to volcanic-rock units is exemplified by the large difference in their estimated maximum hydraulic conductivity; 4,000 and 400 feet per day, respectively. Simulated minimum estimates of hydraulic conductivity are inexact and represent the lower detection limit of the method. Minimum thicknesses of lithologic intervals also were defined for comparing AnalyzeHOLE results to hydraulic properties in regional ground-water flow models.

Scientific Investigations Report↗
USGS1999

Age, double porosity, and simple reaction modifications for the MOC3D ground-water transport model

This report documents modifications for the MOC3D ground-water transport model to simulate (a) ground-water age transport; (b) double-porosity exchange; and (c) simple but flexible retardation, decay, and zero-order growth reactions. These modifications are incorporated in MOC3D version 3.0. MOC3D simulates the transport of a single solute using the method-of-characteristics numerical procedure. The age of ground water, that is the time since recharge to the saturated zone, can be simulated using the transport model with an additional source term of unit strength, corresponding to the rate of aging. The output concentrations of the model are in this case the ages at all locations in the model. Double porosity generally refers to a separate immobile-water phase within the aquifer that does not contribute to ground-water flow but can affect solute transport through diffusive exchange. The solute mass exchange rate between the flowing water in the aquifer and the immobile-water phase is the product of the concentration difference between the two phases and a linear exchange coefficient. Conceptually, double porosity can approximate the effects of dead-end pores in a granular porous media, or matrix diffusion in a fractured-rock aquifer. Options are provided for decay and zero-order growth reactions within the immobile-water phase. The simple reaction terms here extend the original model, which included decay and retardation. With these extensions, (a) the retardation factor can vary spatially within each model layer, (b) the decay rate coefficient can vary spatially within each model layer and can be different for the dissolved and sorbed phases, and (c) a zero-order growth reaction is added that can vary spatially and can be different in the dissolved and sorbed phases. The decay and growth reaction terms also can change in time to account for changing geochemical conditions during transport. The report includes a description of the theoretical basis of the model, a detailed description of input requirements and output options, and the results of model testing and evaluation. The model tests illustrate use of these modifications and demonstrate that accurate solutions can be obtained for these simple cases. Two test cases have no dispersion, illustrating the suitability of this method-of-characteristics model for simulation of advection-dominated transport in ground water.

Water-Resources Investigations Report↗
USGS2000

Hydrogeology and simulation of ground-water flow at Dover Air Force Base, Delaware

Dover Air Force Base in Kent County, Delaware, has many contaminated sites that are in active remediation. To assist in this remediation, a steady-state model of ground-water flow was developed to aid in understanding the hydrology of the system, and for use as a ground-watermanagement tool. This report describes the hydrology on which the model is based, a description of the model itself, and some applications of the model.Dover Air Force Base is underlain by unconsolidated sediments of the Atlantic Coastal Plain. The primary units that were investigated include the upper Calvert Formation and the overlying Columbia Formation. The uppermost sand unit in the Calvert Formation at Dover Air Force Base is the Frederica aquifer, which is the deepest unit investigated in this report. A confining unit of clayey silt in the upper Calvert Formation separates the Frederica aquifer from the lower surficial aquifer, which is the basal Columbia Formation. North and northwest of Dover Air Force Base, the Frederica aquifer subcrops beneath the Columbia Formation and the upper Calvert Formation confining unit is absent. The Calvert Formation dips to the southeast. The Columbia Formation consists predominately of sands, silts, and gravels, although in places there are clay layers that separate the surficial aquifer into an upper and lower surficial aquifer. The areal extent of these clay layers has been mapped by use of gamma logs. Long-term hydrographs reveal substantial changes in both seasonal and annual ground-water recharge. These variations in recharge are related to temporal changes in evaporation, transpiration, and precipitation. The hydrographs show areas where extensive silts and clays are present in the surficial aquifer. In these areas, the vertical gradient between water levels in wells screened above and below the clays can be as large as several feet, and local ground-water highs typically form during normal recharge conditions. When drought conditions persist, water drains off these highs and the vertical gradients decrease. At the south end of Dover Air Force Base, hydrographs of water levels in the Frederica aquifer show that off-Base pumping can cause the water levels to decline below sea level during part of the year.A 4-layer, steady-state numerical model of ground-water flow was developed for Dover Air Force Base and the surrounding area. The upper two layers represent the upper and lower surficial aquifers, which are in the Columbia Formation. In some areas of the model, a semi-confining unit is used to represent an intermittent clay layer between the upper and lower surficial aquifer. This semi-confining unit causes the local groundwater highs in the surficial aquifer. The third model layer represents the upper part of the Calvert Formation, a confining unit. The fourth model layer represents the Frederica aquifer. The model was calibrated to hydraulic heads and to ground-water discharge in Pipe Elm Branch, both of which were measured in September 1997. For the calibrated model, the root-mean-squared errors for the hydraulic heads and the ground-water discharge in the Pipe Elm Branch were 9 percent of the range of head and 3 percent of discharge, respectively. Heads simulated by use of the model were consistent with a map showing average water levels in the region. The U.S. Geological Survey?s MODPATH program was used to simulate ground-water-flow directions for several areas on the Base. This analysis showed the effects of the local groundwater highs. In these areas, ground water can flow from the highs and then dramatically change flow direction as it enters the lower surficial aquifer. The steady-state model has several limitations. The entire ground-water system is under transient hydraulic conditions, due mainly to seasonal and yearly changes in recharge and to withdrawal from irrigation wells. Yet this steady-state model is still considered to be an effective tool for understanding the ground-water-flow system u

Water-Resources Investigations Report↗
USGS2007

Hydrogeologic characterization of the Brazos River Alluvium Aquifer, Bosque County to Fort Bend County, Texas

Introduction The Brazos River alluvium aquifer underlies the Brazos River in Texas from Bosque County to Fort Bend County. The aquifer, one of 21 minor aquifers in the State, supplies water for irrigation, domestic, stock, and commercial use. The Brazos River alluvium aquifer likely will become more important in the future as demand for water increases statewide. A thorough understanding of the hydrogeology of the alluvium aquifer will be the foundation for future studies in the area. During October 2006-April 2007, the U.S. Geological Survey, in cooperation with the Texas Water Development Board, conducted a study to delineate the altitude of the top, altitude of the base, and thickness of the Brazos River alluvium aquifer, and to compile and summarize available hydraulic property (specific capacity, transmissivity, and hydraulic conductivity) data. A digital elevation model was used as the altitude of the top of the aquifer. The altitude of the base of the aquifer was generated using data from wells. The study area encompasses the Brazos River alluvium aquifer in parts of Bosque, Hill, McLennan, Falls, Robertson, Milam, Brazos, Burleson, Grimes, Washington, Waller, Austin, and Fort Bend Counties and a 1.5-mile-wide lateral buffer adjacent to the aquifer. The results of this study will be used by the Texas Water Development Board for input into a ground-water availability model.

Scientific Investigations Map↗
USGS2007

Application of surface geophysical methods, with emphasis on magnetic resonance soundings, to characterize the hydrostratigraphy of the Brazos River alluvium aquifer, College Station, Texas, July 2006: A pilot study

The U.S. Geological Survey, in cooperation with the Texas Water Development Board, used surface geophysical methods at the Texas A&M University Brazos River Hydrologic Field Research Site near College Station, Texas, in a pilot study, to characterize the hydrostratigraphic properties of the Brazos River alluvium aquifer and determine the effectiveness of the methods to aid in generating an improved ground-water availability model. Three non-invasive surface geophysical methods were used to characterize the electrical stratigraphy and hydraulic properties and to interpret the hydrostratigraphy of the Brazos River alluvium aquifer. Two methods, time-domain electromagnetic (TDEM) soundings and two-dimensional direct-current (2D–DC) resistivity imaging, were used to define the lateral and vertical extent of the Ships clay, the alluvium of the Brazos River alluvium aquifer, and the underlying Yegua Formation. Magnetic resonance sounding (MRS), a recently developed geophysical method, was used to derive estimates of the hydrologic properties including percentage water content and hydraulic conductivity. Results from the geophysics study demonstrated the usefulness of combined TDEM, 2D–DC resistivity, and MRS methods to reduce the need for additional boreholes in areas with data gaps and to provide more accurate information for ground-water availability models. Stratigraphically, the principal finding of this study is the relation between electrical resistivity and the depth and thickness of the subsurface hydrostratigraphic units at the site. TDEM data defined a three-layer electrical stratigraphy corresponding to a conductor-resistor-conductor that represents the hydrostratigraphic units—the Ships clay, the alluvium of the Brazos River alluvium aquifer, and the Yegua Formation. Sharp electrical boundaries occur at about 4 to 6 and 20 to 22 meters below land surface based on the TDEM data and define the geometry of the more resistive Brazos River alluvium aquifer. Variations in resistivity in the alluvium aquifer range from 10 to more than 175 ohm-meters possibly are caused by lateral changes in grain size. Resistivity increases from east to west along a profile away from the Brazos River, which signifies an increase in grain size within the alluvium aquifer and therefore a more productive zone with more abundant water in the aquifer. MRS data can help delineate the subsurface hydrostratigraphy and identify the geometric boundaries of the hydrostratigraphic units by identifying changes in the free water content, transmissivity, and hydraulic conductivity. MRS data indicate that most productive zones of the alluvium aquifer occur between 12 and 25 meters below land surface in the western part of the study area where the hydraulic conductivity can be as high as 250 meters per day. Hydrostratigraphically, individual hydraulic conductivity values derived from MRS were consistent with those from aquifer tests conducted in 1996 in the study area. Average hydraulic conductivity values from the aquifer tests range from about 61 to 80 meters per day, whereas the MRS-derived hydraulic conductivity values range from about 27 to 97 meters per day. Interpreting an interpolated profile of the hydraulic conductivity values and individual values derived from MRS can help describe the hydrostratigraphic framework of an area and constrain ground-water models for better accuracy.

Texas↗
USGS2005

Surface-geophysical investigation of a formerly used defense site, Machiasport, Maine, February 2003

Surface-geophysical surveys were conducted in February 2003 at a formerly used defense site in Maine, where residual chlorinated solvents are affecting off-site domestic water-supply wells. The U.S. Geological Survey and Argonne National Laboratory used surface-geophysical methods, including ground-penetrating radar and seismic-refraction tomography, to characterize the lithology and structure of the bedrock at the site and to identify highly fractured areas that may provide pathways for ground-water flow and contaminant transport. Multifrequency electromagnetic and inductive terrain-conductivity methods also were evaluated, but these techniques were adversely affected by a nearby naval computer and telecommunications station. Interpretation of the data from ground-penetrating radar indicates that depth to the weathered bedrock surface is approximately 0.5 to 3 meters. Reflections from within the bedrock are visible throughout all ground-penetrating radar profiles, and zones of scattered electromagnetic energy may correlate to zones of highly fractured bedrock. Interpretation of the data from seismic-refraction tomography inversion indicates that zones of relatively low seismic velocity and topographic lows may correlate with fractured and water-producing intervals within the bedrock. Integrated interpretation of the results from ground-penetrating radar and seismic-refraction tomography was used to locate boreholes along the surface-geophysical profiles. An integrated analysis of information obtained from the surface- and borehole-geophysical surveys and test drilling will be used by the U.S. Army Corps of Engineers to develop a conceptual model of ground-water flow and solute transport at the site.

Maine↗
USGS2001

Simulation of the interaction of karstic lakes Magnolia and Brooklyn with the upper Floridan Aquifer, southwestern Clay County, Florida

The stage of Lake Brooklyn, in southwestern Clay County, Florida, has varied over a range of 27 feet since measurements by the U.S. Geological Survey began in July 1957. The large stage changes have been attributed to the relation between highly transient surface-water inflow to the lake and subsurface conduits of karstic origin that permit a high rate of leakage from the lake to the Upper Floridan aquifer. After the most recent and severe stage decline (1990-1994), the U.S. Geological Survey began a study that entailed the use of numerical ground-water flow models to simulate the interaction of the lake with the Upper Floridan aquifer and the large fluctuations of stage that were a part of that process. A package (set of computer programs) designed to represent lake/aquifer interaction in the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model (MODFLOW-96) and the Three-Dimensional Method-of-Characteristics Solute-Transport Model (MOC3D) simulators was prepared as part of this study, and a demonstration of its capability was a primary objective of the study. (Although the official names are Brooklyn Lake and Magnolia Lake (Florida Geographic Names), in this report the local names, Lake Brooklyn and Lake Magnolia, are used.) In the simulator of lake/aquifer interaction used in this investigation, the stage of each lake in a simulation is updated in successive time steps by a budget process that takes into account ground-water seepage, precipitation upon and evaporation from the lake surface, stream inflows and outflows, overland runoff inflows, and augmentation or depletion by artificial means. The simulator was given the capability to simulate both the division of a lake into separate pools as lake stage falls and the coalescence of several pools into a single lake as the stage rises. This representational capability was required to simulate Lake Brooklyn, which can divide into as many as 10 separate pools at sufficiently low stage. In the first of two calibrated models, recharge to the water table, specified as a monthly rate, was set equal to 40 percent of the monthly rainfall rate. The specified rate of inflow to the uppermost stream segment was set equal to outflows from Lake Lowry estimated from lake stage and the 1994-97 rating table. Leakage to the intermediate and Upper Floridan aquifers was assumed to occur from the surficial aquifer system through the confining layers directly beneath deeper parts of the lake bottom. A leakance coefficient value of 0.001 feet per day per foot of thickness was used beneath Lake Magnolia, and a value of 0.005 feet per day per foot of thickness was used beneath most of Lake Brooklyn. With these values, the conductance through the confining layers beneath Lake Brooklyn was about 19 times that beneath Lake Magnolia. The simulated stages of Lake Brooklyn matched the measured stages reasonably well in the early (1957-72) and later (1990-98) parts of the simulation time period, but the match was unsatisfactory in an intermediate time period (1973-89). To resolve this discrepancy, the hypothesis was proposed that undocumented losses of water from Alligator Creek upstream from Lake Brooklyn or from the lake itself occurred between 1973 and 1989 when there was sufficient streamflow. The resulting simulation of lake stages matched the measured lake stages accurately during the entire simulation time period. The model was then revised to incorporate the assumption that only 20 percent of precipitation recharged the water table (the second calibrated model). Recalibration of the model required that leakance values for the confining units under deeper parts of the lakes also be reduced by nearly 50 percent. The stages simulated with the new parameter assumptions, but retaining the assumption of surface-water losses, were an excellent match of the measured values. The stage of Lake Magnolia was also simulated accurately. The results of sensitivity analyses show that simulated streamflow between Lakes Magnolia and Brooklyn tends to be water-budget controlled, and is not appreciably affected by the specified outflow altitude or channel characteristics of the receiving stream. To match heads measured in observation wells of the surficial aquifer network, the assigned hydraulic conductivity values were zoned, and ranged from a minimum of 4 feet per day to a maximum of 400 feet per day in the first calibrated model. These values were reduced by about 50 percent in the second calibrated model. Differences between observation wells were noted in the abruptness of changes of measured head values, and in the relation of the timing of peak measured heads and simulated peak heads. These differences seemed to be correlated with the depth of the water table below land surface. Spatially uniform values of transmissivity were specified for the intermediate (10,000 feet squared per day) and Upper Floridan (100,000 feet squared per day) aquifers. Simulated heads in the Upper Floridan aquifer layer follow the trend of the heads measured in a long-term observation well with data beginning in 1960. This result suggests that the observed head decline could be explained entirely in terms of the stage decline in Lake Brooklyn and may not indicate a regional trend.

Florida↗
USGS2006

User's guide to the Variably Saturated Flow (VSF) process to MODFLOW

A new process for simulating three-dimensional (3-D) variably saturated flow (VSF) using Richards' equation has been added to the 3-D modular finite-difference ground-water model MODFLOW. Five new packages are presented here as part of the VSF Process--the Richards' Equation Flow (REF1) Package, the Seepage Face (SPF1) Package, the Surface Ponding (PND1) Package, the Surface Evaporation (SEV1) Package, and the Root Zone Evapotranspiration (RZE1) Package. Additionally, a new Adaptive Time-Stepping (ATS1) Package is presented for use by both the Ground-Water Flow (GWF) Process and VSF. The VSF Process allows simulation of flow in unsaturated media above the ground-water zone and facilitates modeling of ground-water/surface-water interactions. Model performance is evaluated by comparison to an analytical solution for one-dimensional (1-D) constant-head infiltration (Dirichlet boundary condition), field experimental data for a 1-D constant-head infiltration, laboratory experimental data for two-dimensional (2-D) constant-flux infiltration (Neumann boundary condition), laboratory experimental data for 2-D transient drainage through a seepage face, and numerical model results (VS2DT) of a 2-D flow-path simulation using realistic surface boundary conditions. A hypothetical 3-D example case also is presented to demonstrate the new capability using periodic boundary conditions (for example, daily precipitation) and varied surface topography over a larger spatial scale (0.133 square kilometer). The new model capabilities retain the modular structure of the MODFLOW code and preserve MODFLOW's existing capabilities as well as compatibility with commercial pre-/post-processors. The overall success of the VSF Process in simulating mixed boundary conditions and variable soil types demonstrates its utility for future hydrologic investigations. This report presents a new flow package implementing the governing equations for variably saturated ground-water flow, four new boundary condition packages unique to unsaturated flow, the Adaptive Time-Stepping Package for use with both the GWF Process and the new VSF Process, detailed descriptions of the input and output files for each package, and six simulation examples verifying model performance.

Techniques and Methods↗
USGS1992

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 2: Derivation of finite-element equations and comparisons with analytical solutions

MODFE, a modular finite-element model for simulating steady- or unsteady-state, area1 or axisymmetric flow of ground water in a heterogeneous anisotropic aquifer is documented in a three-part series of reports. In this report, part 2, the finite-element equations are derived by minimizing a functional of the difference between the true and approximate hydraulic head, which produces equations that are equivalent to those obtained by either classical variational or Galerkin techniques. Spatial finite elements are triangular with linear basis functions, and temporal finite elements are one dimensional with linear basis functions. Physical processes that can be represented by the model include (1) confined flow, unconfined flow (using the Dupuit approximation), or a combination of both; (2) leakage through either rigid or elastic confining units; (3) specified recharge or discharge at points, along lines, or areally; (4) flow across specified-flow, specified-head, or head-dependent boundaries; (5) decrease of aquifer thickness to zero under extreme water-table decline and increase of aquifer thickness from zero as the water table rises; and (6) head-dependent fluxes from springs, drainage wells, leakage across riverbeds or confining units combined with aquifer dewatering, and evapotranspiration. The matrix equations produced by the finite-element method are solved by the direct symmetric-Doolittle method or the iterative modified incomplete-Cholesky conjugate-gradient method. The direct method can be efficient for small- to medium-sized problems (less than about 500 nodes), and the iterative method is generally more efficient for larger-sized problems. Comparison of finite-element solutions with analytical solutions for five example problems demonstrates that the finite-element model can yield accurate solutions to ground-water flow problems.

Techniques of Water-Resources Investigations↗