Geology ReportsSearch

SEARCH · Geology Reports

Results for “Cells”

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 415 records · Page 23Linked to original sources

Potential effects of existing and proposed groundwater withdrawals on water levels and natural groundwater discharge in Snake Valley, Juab and Millard Counties, Utah, White Pine County, Nevada, and surrounding areas in Utah and Nevada

Applications have been filed for several water-right changes and new water rights, with total withdrawals of about 1,800 acre-feet per year, in Snake Valley near Eskdale and Partoun, Utah. The Bureau of Land Management has identified 11 sites where the Bureau of Land Management holds water rights and 7 other springs of interest that could be affected by these proposed groundwater withdrawals. This report presents a hydrogeologic analysis of areas within Snake Valley to assess the potential effects on Bureau of Land Management water rights and other springs of interest resulting from existing and proposed groundwater withdrawals. A previously developed numerical groundwater-flow model was used to quantify potential groundwater drawdown and the capture, or groundwater withdrawals that results in depletion, of natural discharge resulting from existing and proposed groundwater withdrawals within Snake Valley. Existing groundwater withdrawals were simulated for a 50-year period prior to adding the newly proposed withdrawals to bring the model from pre-development conditions to the start of 2014. After this initial 50-year period, existing withdrawals, additional proposed withdrawals, and consequent effects were simulated for periods of 5, 10, 25, 50, and 100 years. Downward trends in water levels measured in wells indicate that the existing groundwater withdrawals in Snake Valley are affecting water levels. The numerical model simulated similar downward trends in water levels. The largest simulated drawdowns caused by existing groundwater withdrawals ranged between 10 and 26 feet and were near the centers of the agricultural areas by Callao, Eskdale, Baker, Garrison, and along the Utah-Nevada state line in southern Snake Valley. The largest simulated water-level declines were at the Bureau of Land Management water-rights sites near Eskdale, Utah, where simulated drawdowns ranged between 2 and 8 feet at the start of 2014. These results were consistent with, but lower than, observations from several wells monitored by the U.S. Geological Survey that indicated water-level declines of 6 to 18 feet near the Eskdale area since the mid-1970s and 1980s. The model cells where the simulated capture of natural groundwater discharge resulting from the existing withdrawals was greatest were those containing Kane Spring, Caine Spring, and Unnamed Spring 5, where existing groundwater withdrawals capture 13 to 29 percent of the total simulated natural discharge in these cells. Simulated drawdown and simulated capture of natural groundwater discharge resulting from the proposed withdrawals started in as few as 5 years at seven of the sites. After 100 years, four sites showed simulated drawdowns ranging between 1 and 2 feet; eight sites showed simulated drawdowns ranging between 0.1 and 0.9 feet; and five sites showed no simulated drawdown resulting from the proposed withdrawals. The largest amounts of simulated capture of natural groundwater discharge resulting from the proposed withdrawals after 100 years were in the model cells containing Coyote Spring, Kane Spring, and Caine Spring, which had capture amounts ranging between 5.5 and 9.1 percent of the total simulated natural discharge in these cells.

Nevada, Utah

Gas chemistry and isotope data for volcano monitoring at the Lassen Volcanic Center, Lassen Volcanic National Park

This report presents chemical and isotopic compositions of volcanic gases collected from thermal areas within Lassen Volcanic National Park in northern California from 1974 through 2019. As the southernmost volcano in the Cascade Range and designated a very-high-threat volcano by the U.S. Geological Survey, the Lassen Volcanic Center (LVC) requires consistent monitoring to assess potential volcanic hazards. In 2014, the California Volcano Observatory established a gas geochemical monitoring program at LVC to provide baseline data to evaluate future changes.</p> Results demonstrate consistent spatial patterns in bulk gas chemistry that support a two-circulation-cell hydrothermal model previously established for LVC. Gas samples from circulation cell 1 thermal areas have higher helium isotope ratios (6.59–7.50 times the air value) than those from circulation cell 2 (5.86–6.52 times the air value), indicating a stronger magmatic signature. The Sulphur Works and Pilot Pinnacle thermal areas within circulation cell 1 consistently emit gases with the highest magmatic helium contents, suggesting gas at these areas best represents conditions in the underlying volcanic system. A slight decrease in helium isotope values since 1974 may indicate progressive dilution of magmatic helium-3 (3He) by radiogenic helium-4 (4He) in the absence of recent magma intrusion. Carbon isotope compositions of carbon dioxide across all thermal areas are relatively uniform (−9.7–−7.3 per mil), falling within the range observed at other Cascade Range volcanoes. Based on gas geochemical characteristics and site accessibility, the Sulphur Works and Pilot Pinnacle thermal areas represent optimal targets for continued monitoring of the LVC magmatic-hydrothermal system. This study includes the most comprehensive helium isotope dataset collected at LVC currently available and establishes critical baseline data for future volcanic monitoring efforts.</p>

Open-File Report

Osmotic potential and projected drought tolerance of four phreatophytic shrub species in Owens Valley, California

A large part of the water used by plant communities growing on the floor of Owens Valley, California, is derived from a shallow unconfined aquifer. Fluctuations in the water table caused by groundwater withdrawal may result in periods when this water supply is not accessible to plants. The capacity of the plants to adapt to these periods of water loss depend on the availability of water stored in the soil and on physiological characteristics related to the ability of the plants to resist dehydration and wilting. Osmotic adjustment occurred in four phreatophytic shrub species at sites near bishop, California, where the water table had been lowered by a system of pump-equipped wells installed in the vicinity of vegetation transects. The pressure-volume techniques was used to determine osmotic potential and cell-wall elasticity between March 1985 and September 1986 for Atriplex torreyi, Chrysothamnus nauseosus , Sarcobatus vermiculatus, and Artemisia tridentata. Although not usually classified as a phreatophyte, Artemisia tridentata, where it grows on the valley floor, is apparently dependent on the depth to the water table. During late summer, osmotic potentials were 0.37 to 0.41 megapascal lower in plants growing on the site where the water table had been lowered compared to an adjacent site where the water table remained at its natural levels. Measurements of soil matric potential at the two sites indicated that osmotic adjustment occurred in response to stress caused by lowering the water table. A theoretical lower limit of osmotic adjustment was determined by comparing initial cell osmotic potentials with initial xylem water potentials. These experimentally derived limits indicated that A. torreyi and S. vermiculatus may maintain leaf cell turgor at significantly lower cell water potentials (about -4.5 megapascals) than C. nauseosus or A. tridentata (about -2.5 megapascals) and allows them to function in dryer soil environments. (Author 's abstract)

California

Documentation of computer programs to compute and display pathlines using results from the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model

A particle tracking post-processing package was developed to compute three-dimensional path lines based on output from steady-state simulations obtained with the U.S. Geological Survey modular 3-dimensional finite difference groundwater flow model. The package consists of two FORTRAN 77 computer programs: (1) MODPATH, which calculates pathlines, and (2) MODPATH-PLOT, which presents results graphically. MODPATH uses a semi-analytical particle tracking scheme. The method is based on the assumption that each directional velocity component varies linearly within a grid cell in its own coordinate direction. This assumption allows an analytical expression to be obtained describing the flow path within a grid cell. Given the initial position of a particle anywhere in a cell, the coordinates of any other point along its path line within the cell, and the time of travel between them, can be computed directly. Data is input to MODPATH and MODPATH-PLOT through a combination of files and interactive dialogue. Examples of how to use MODPATH and MODPATH-PLOT are provided for a sample problem. Listings of the computer codes and detailed descriptions of input data format and program options are also presented. (Author 's abstract)

Open-File Report

Simulation of ground-water flow in coastal Georgia and adjacent parts of South Carolina and Florida-predevelopment, 1980, and 2000

A digital model was developed to simulate steady-state ground-water flow in a 42,155-square-mile area of coastal Georgia and adjacent parts of South Carolina and Florida. The model was developed to (1) understand and refine the conceptual model of regional ground-water flow, (2) serve as a framework for the development of digital subregional ground-water flow and solute-transport models, and (3) serve as a tool for future evaluations of hypothetical pumping scenarios used to facilitate water management in the coastal area. Single-density ground-water flow was simulated using the U.S. Geological Survey finite-difference code MODFLOW-2000 for mean-annual conditions during predevelopment (pre?1900) and the years 1980 and 2000. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units. A combination of boundary conditions was applied, including a general-head boundary condition on the top active cells of the model and a time-variable fixed-head boundary condition along part of the southern lateral boundary. Simulated heads for 1980 and 2000 conditions indicate a good match to observed values, based on a plus-or-minus 10-foot (ft) calibration target and calibration statistics. The root-mean square of residual water levels for the Upper Floridan aquifer was 13.0 ft for the 1980 calibration and 9.94 ft for the 2000 calibration. Some spatial patterns of residuals were indicated for the 1980 and 2000 simulations, and are likely a result of model-grid cell size and insufficiently detailed hydraulic-property and pumpage data in some areas. Simulated potentiometric surfaces for predevelopment, 1980, and 2000 conditions all show major flow system features that are indicated by estimated peotentiometric maps. During 1980?2000, simulated water levels at the centers of pumping at Savannah and Brunswick rose more than 20 ft and 8 ft, respectively, in response to decreased pumping. Simulated drawdown exceeded 10 ft in the Upper Floridan aquifer across much of the western half of the model area, with drawdown exceeding 20 ft along parts of the western, northern, and southern boundaries where irrigation pumping increased during this period. From predevelopment to 2000 conditions, the simulated water budget showed an increase in inflow from, and decrease in outflow to, the general-head boundaries, and a reversal from net seaward flow to net landward flow across the coastline. Simulated changes in recharge and discharge distribution from predevelopment to 2000 conditions showed an increase in extent and magnitude of net recharge cells in the northern part of the model area, and a decrease in discharge or change to recharge in cells containing major streams and beneath major pumping centers. The model is relatively sensitive to pumping and the controlling head at the fixed-head boundary and less sensitive to the distribution of aquifer properties in general. Model limitations include: (1) its spatial scale and discretization, (2) the extent to which data are available to physically define the flow system, (3) the type of boundary conditions and controlling parameters used, (4) uncertainty in the distribution of pumping, and (5) uncertainty in field-scale hydraulic properties. The model could be improved with more accurate estimates of ground-water pumpage and better characterization of recharge and discharge.

Florida, Georgia, South Carolina

Development, Testing, and Sensitivity and Uncertainty Analyses of a Transport and Reaction Simulation Engine (TaRSE) for Spatially Distributed Modeling of Phosphorus in South Florida Peat Marsh Wetlands

Alterations to the predevelopment delivery of water and nutrients into the Everglades of southern Florida have been occurring for nearly a century. Major regional drainage projects, large-scale agricultural development, and changes to the hydrology of the Kissimmee River-Lake Okeechobee watershed have resulted in substantial phosphorus transport increases by surface waters. Excess phosphorus has accumulated in the soils of northern Everglades marshes to levels that have impaired the natural resources of the region. Regulations now limit the amount of phosphorous that enters the Everglades through an extensive network of water-control structures. This study involved the development and application of water-quality modeling components that may be applied to existing hydrologic models of southern Florida to evaluate the effects of different management scenarios. The result of this work is a spatially distributed water-quality model for phosphorus transport and cycling in wetlands. The model solves the advection-dispersion equation on an unstructured triangular mesh and incorporates a wide range of user-selectable mechanisms for phosphorus uptake and release parameters. In general, the phosphorus model contains transfers between stores; examples of stores that can be included are soil, water column (solutes), pore water, macrophytes, suspended solids (plankton), and biofilm. Examples of transfers are growth, senescence, settling, diffusion, and so forth, described with first order, second order, and Monod types of transformations. Local water depths and velocities are determined from an existing two-dimensional, overland-flow hydrologic model. The South Florida Water Management District Regional Simulation Model was used in this study. The model is applied to three case studies: intact cores of wetland soils with water, outdoor mesocosoms, and a large constructed wetland; namely, Cell 4 of Stormwater Treatment Area 1 West (STA-1W Cell 4). Different levels of complexity in the phosphorus cycling mechanisms were simulated in these case studies using different combinations of phosphorus reaction equations. Changes in water column phosphorus concentrations observed under the controlled conditions of laboratory incubations, and mesocosm studies were reproduced with model simulations. Short-term phosphorus flux rates and changes in phosphorus storages were within the range of values reported in the literature, whereas unknown rate constants were used to calibrate the model output. In STA-1W Cell 4, the dominant mechanism for phosphorus flow and transport is overland flow. Over many life cycles of the biological components, however, soils accrue and become enriched in phosphorus. Inflow total phosphorus concentrations and flow rates for the period between 1995 and 2000 were used to simulate Cell 4 phosphorus removal, outflow concentrations, and soil phosphorus enrichment over time. This full-scale application of the model successfully incorporated parameter values derived from the literature and short-term experiments, and reproduced the observed long-term outflow phosphorus concentrations and increased soil phosphorus storage within the system. A global sensitivity and uncertainty analysis of the model was performed using modern techniques such as a qualitative screening tool (Morris method) and the quantitative, variance-based, Fourier Amplitude Sensitivity Test (FAST) method. These techniques allowed an in-depth exploration of the effect of model complexity and flow velocity on model outputs. Three increasingly complex levels of possible application to southern Florida were studied corresponding to a simple soil pore-water and surface-water system (level 1), the addition of plankton (level 2), and of macrophytes (level 3). In the analysis for each complexity level, three surface-water velocities were considered that each correspond to residence times for the selected area (1-kilometer long) of 2, 10, and 20

Scientific Investigations Report

Three-dimensional numerical model of ground-water flow in northern Utah Valley, Utah County, Utah

A three-dimensional, finite-difference, numerical model was developed to simulate ground-water flow in northern Utah Valley, Utah. The model includes expanded areal boundaries as compared to a previous ground-water flow model of the valley and incorporates more than 20 years of additional hydrologic data. The model boundary was generally expanded to include the bedrock in the surrounding mountain block as far as the surface-water divide. New wells have been drilled in basin-fill deposits near the consolidated-rock boundary. Simulating the hydrologic conditions within the bedrock allows for improved simulation of the effect of withdrawal from these wells. The inclusion of bedrock also allowed for the use of a recharge model that provided an alternative method for spatially distributing areal recharge over the mountains. The model was calibrated to steady- and transient-state conditions. The steady-state simulation was developed and calibrated by using hydrologic data that represented average conditions for 1947. The transient-state simulation was developed and calibrated by using hydrologic data collected from 1947 to 2004. Areally, the model grid is 79 rows by 70 columns, with variable cell size. Cells throughout most of the model domain represent 0.3 mile on each side. The largest cells are rectangular with dimensions of about 0.3 by 0.6 mile. The largest cells represent the mountain block on the eastern edge of the model domain where the least hydrologic data are available. Vertically, the aquifer system is divided into 4 layers which incorporate 11 hydrogeologic units. The model simulates recharge to the ground-water flow system as (1) infiltration of precipitation over the mountain block, (2) infiltration of precipitation over the valley floor, (3) infiltration of unconsumed irrigation water from fields, lawns, and gardens, (4) seepage from streams and canals, and (5) subsurface inflow from Cedar Valley. Discharge of ground water is simulated by the model to (1) flowing and pumping wells, (2) drains and springs, (3) evapotranspiration, (4) Utah Lake, (5) the Jordan River and mountain streams, and (6) Salt Lake Valley by subsurface outflow through the Jordan Narrows. During steady-state calibration, variables were adjusted within probable ranges to minimize differences between model-computed and measured water levels as well as between model-computed and independently estimated flows that include: recharge by seepage from individual streams and canals, discharge by seepage to individual streams and the Jordan River, discharge to Utah Lake, discharge to drains and springs, discharge by evapotranspiration, and subsurface flows into and out of northern Utah Valley from Cedar Valley and to Salt Lake Valley, respectively. The transient-state simulation was calibrated to measured water levels and water-level changes with consideration given to annual changes in the flows listed above.

Utah

Effects of Potential Changes in Groundwater Withdrawals from the Sparta Aquifer on Water-Level Altitudes in Jefferson County, Arkansas

A groundwater-flow model of the Sparta aquifer was used to evaluate changes in water-level altitudes associated with the withdrawal of groundwater at varying rates from a well field near Pine Bluff, Arkansas, in Jefferson County. Water-level altitudes at three different model cell locations from five different scenarios for varying withdrawal rates from the well field were compared for the period 1998 to 2048. The three model cells used for the comparison were located (1) near the center of the well field, (2) near the center of the city of Pine Bluff (about 5 miles west of the center of the well field), and (3) about 15 miles north of the well field. Pumping rates at the well field were varied from 7.2 million gallons per day to 27 million gallons per day for the five scenarios analyzed, and water-level hydrographs were constructed for each scenario for each of the three model cell locations. Water-level altitudes near the center of the well field changed the most of the three model cell locations analyzed. Water-level altitudes were approximately 90 feet higher for the 7.2 million gallon per day scenario in 2048 compared to the baseline scenario of 25.4 million gallons per day. Whereas, water-level altitudes at the same location were 9 feet lower for the 27 million gallon per day scenario in 2048 compared to the baseline scenario.

Scientific Investigations Report

Approaches to highly parameterized inversion: Pilot-point theory, guidelines, and research directions

Pilot points have been used in geophysics and hydrogeology for at least 30 years as a means to bridge the gap between estimating a parameter value in every cell of a model and subdividing models into a small number of homogeneous zones. Pilot points serve as surrogate parameters at which values are estimated in the inverse-modeling process, and their values are interpolated onto the modeling domain in such a way that heterogeneity can be represented at a much lower computational cost than trying to estimate parameters in every cell of a model. Although the use of pilot points is increasingly common, there are few works documenting the mathematical implications of their use and even fewer sources of guidelines for their implementation in hydrogeologic modeling studies. This report describes the mathematics of pilot-point use, provides guidelines for their use in the parameter-estimation software suite (PEST), and outlines several research directions. Two key attributes for pilot-point definitions are highlighted. First, the difference between the information contained in the every-cell parameter field and the surrogate parameter field created using pilot points should be in the realm of parameters which are not informed by the observed data (the null space). Second, the interpolation scheme for projecting pilot-point values onto model cells ideally should be orthogonal. These attributes are informed by the mathematics and have important ramifications for both the guidelines and suggestions for future research.

Scientific Investigations Report

Hydrology and water quality of a dune-and-swale wetland adjacent to the Grand Calumet River, Indiana, 2019–22

Adverse ecological and water-quality effects associated with industrial land-use changes are common for littoral wetlands connected to river mouth ecosystems in the Grand Calumet River-Indiana Harbor Canal Area of Concern. These effects can be exacerbated by recent high Lake Michigan water levels that are problematic for wetland restoration. Wetlands in the adjacent Clark and Pine Nature Preserve and Pine Station Nature Preserve are intended to mitigate wetland destruction in the area of concern by restoring residual dune-and-swale wetlands and preserving habitat for endangered and threatened plant species. Physical hydrology and water-quality monitoring of restored wetland cells at the preserves were initiated during 2019 to evaluate changes after wetland restoration efforts in 2015 and near record-low water levels in early 2013. Lake Michigan water levels rose steadily between late 2013 and 2018 to record-high water levels in 2019 and 2020. In this report, precipitation, evapotranspiration, and groundwater and surface-water levels are analyzed to better understand wetland inundation controls and flow directions in restored northern dune-and-swale wetland settings relative to the Grand Calumet River. Continuous specific conductance data and discrete water-quality samples were collected and analyzed to provide a synoptic view of water quality for the restored wetlands. High Lake Michigan water levels affected Grand Calumet River stage and shallow groundwater elevations in the study area after the onset of peak lake levels in June 2019, that persisted through summer 2020, before finally receding in September 2020. Grand Calumet River stage peaked soon after lake levels in July 2019, whereas groundwater elevations in the study area peaked in October 2019. Specific conductance values in closed-basin wetland cells in the western and central parts of the nature preserves indicated a dilution trend and contrasted those of interconnected wetland cells along an eastern corridor, where alterations to wetland cells were more pronounced. Monitoring results indicate that varying seasonal wetland inundation trends with low stands in autumn have returned after high water table conditions owing to high water levels on Lake Michigan. Wetland water balance results during the study period indicated that the wetland ecosystem partially moderated flooding during high lake levels through summer evapotranspiration.

Indiana

Mineralogical applications of electron diffraction. 1. Theory and techniques

The small wavelengths used in electron-diffraction experiments and the thinness of the crystals necessary for the transmission of the electron beam combine to require a somewhat different diffraction geometry for the interpretation of electron-diffraction patterns than is used in the interpretation of X-ray diffraction patterns. This geometry, based on the reciprocal lattice concept and geometrical construction of Ewald, needed for the interpretation. of transmission electron-diffraction single-crystal patterns is here reviewed. Transmission electron-diffraction single-crystal patterns of two monoclinic substances, colemanite [CaB 3 O 3 (OH) 3 •H 2 O] and potassium chlorate (KC10 3 ), are examined and the .theory necessary for their interpretation is given in detail. The study of these patterns furnishes a basis for the interpretation of single-crystal patterns of materials belonging to any crystal system. It is shown that useful unit-cell data, accurate to a few tenths of a percent, can be obtained from the patterns of colemanite and KClO 3 . A method of evaluating unit-cell data from measurements of such single-crystal patterns is given. The transmission electron-diffraction powder pattern obtained from an oriented aggregate of thin crystals gives the same unit-cell data as are given by the electron-diffraction single-crystal pattern obtained from one crystal of the aggregate., A graphical method is given for precisely evaluating unit-cell constants from measurements of such a powder pattern.

Trace Elements Investigations

Computational considerations for collecting and using data in the equidistant cylindrical map projection and the bounds of sampling geographic data at progressively higher resolution

The Equidistant Cylindrical Map projection is popular with digital modelers and others for storing and processing worldwide data sets because of the simple association of latitude and longitude to cell values or pixels in the resulting grid. This projection does not accurately display area, and the diminished geographic area represented by cells at high latitudes is not often carefully considered. A simple mathematical analysis quantifies the discrepancy in area sampled by cells at different latitudes. The presence of this discrepancy indicates that the use of this projection can induce bias in data sets when both sampling and reporting data. It is demonstrated that as the resolution requirements of input data for models increase, the necessity of providing data to accurately describe smaller cells, particularly at high latitude, will be a challenge.

Techniques and Methods

Documentation for the “XT3D” option in the Node Property Flow (NPF) Package of MODFLOW 6

This report describes the “XT3D” option in the Node Property Flow (NPF) Package of MODFLOW 6. The XT3D option extends the capabilities of MODFLOW by enabling simulation of fully three-dimensional anisotropy on regular or irregular grids in a way that properly takes into account the full, three-dimensional conductivity tensor. It can also improve the accuracy of groundwater-flow simulations in cases in which the model grid violates certain geometric requirements. Three example problems demonstrate the use of the XT3D option to simulate groundwater flow on irregular grids and through three-dimensional porous media with anisotropic hydraulic conductivity. Conceptually, the XT3D method of estimating flow between two MODFLOW 6 model cells can be viewed in terms of three main mathematical steps: construction of head-gradient estimates by interpolation; construction of fluid-flux estimates by application of the full, three-dimensional form of Darcy’s Law, in which the conductivity tensor can be heterogeneous and anisotropic; and construction of the flow expression by enforcement of continuity of flow across the cell interface. The resulting XT3D flow expression, which relates the flow across the cell interface to the values of heads computed at neighboring nodes, is the sum of terms in which conductance-like coefficients multiply head differences, as in the conductance-based flow expression the NPF Package uses by default. However, the XT3D flow expression contains terms that involve “neighbors of neighbors” of the two cells for which the flow is being calculated. These additional terms have no analog in the conductance-based formulation. When assembled into matrix form, the XT3D formulation results in a larger stencil than the conductance-based formulation; that is, each row of the coefficient matrix generally contains more nonzero elements. The “RHS” suboption can be used to avoid expanding the stencil by placing the additional terms on the right-hand side of the matrix equation and evaluating them at the previous iteration or time step. The XT3D option can be an alternative to the Ghost-Node Correction (GNC) Package. However, the XT3D formulation is typically more computationally intensive than the conductance-based formulation the NPF Package uses by default, either with or without ghost nodes. Before deciding whether to use the GNC Package or XT3D option for production runs, the user should consider whether the conductance-based formulation alone can provide acceptable accuracy for the particular problem being solved.

Techniques and Methods

Testing the sensitivity of pumpage to increases in surficial aquifer system heads in the Cypress Creek well-field area, west-central Florida — An optimization technique

Tampa Bay depends on ground water for most of the water supply. Numerous wetlands and lakes in Pasco County have been impacted by the high demand for ground water. Central Pasco County, particularly the area within the Cypress Creek well field, has been greatly affected. Probable causes for the decline in surface-water levels are well-field pumpage and a decade-long drought. Efforts are underway to increase surface-water levels by developing alternative sources of water supply, thus reducing the quantity of well-field pumpage. Numerical ground-water flow simulations coupled with an optimization routine were used in a series of simulations to test the sensitivity of optimal pumpage to desired increases in surficial aquifer system heads in the Cypress Creek well field. The ground-water system was simulated using the central northern Tampa Bay ground-water flow model. Pumping solutions for 1987 equilibrium conditions and for a transient 6-month timeframe were determined for five test cases, each reflecting a range of desired target recovery heads at different head control sites in the surficial aquifer system. Results are presented in the form of curves relating average head recovery to total optimal pumpage. Pumping solutions are sensitive to the location of head control sites formulated in the optimization problem and as expected, total optimal pumpage decreased when desired target head increased. The distribution of optimal pumpage for individual production wells also was significantly affected by the location of head control sites. A pumping advantage was gained for test-case formulations where hydraulic heads were maximized in cells near the production wells, in cells within the steady-state pumping center cone of depression, and in cells within the area of the well field where confining-unit leakance is the highest. More water was pumped and the ratio of head recovery per unit decrease in optimal pumpage was more than double for test cases where hydraulic heads are maximized in cells located at or near the production wells. Additionally, the ratio of head recovery per unit decrease in pumpage was about three times more for the area where confining-unit leakance is the highest than for other leakance zone areas of the well field. For many head control sites, optimal heads corresponding to optimal pumpage deviated from the desired target recovery heads. Overall, pumping solutions were constrained by the limiting recovery values, initial head conditions, and by upper boundary conditions of the ground-water flow model.

Florida

Water quality of Lake Arlington on Village Creek, north-central Texas; 1973 to 1981

Water in Lake Arlington on Village Creek in north-central Texas had volume-weighted average concentrations of less than 240 milligrams per liter of dissolved solids, less than 30 milligrams per liter of dissolved chloride, and less than 40 milligrams per liter of dissolved sulfate between January 29, 1973, and August 20, 1981. The water was moderately hard (hardness greater than 60 but less than 120 milligrams per liter as calcium carbonate). The concentrations of each of these constituents were usually greatest during winter, especially during the first few years of the study, but decreased significantly as discharges of municipal wastes were systematically and progressively diverted to treatment facilities downstream from the reservoir. Thermal stratification in Lake Arlington usually begins during March or April and persists until October. Thermal stratification has resulted in significant seasonal and areal variations in the concentration of dissolved oxygen, dissolved iron, dissolved manganese, total inorganic nitrogen, and total phosphorus. Oxygen utilized in the decay of organic matter is not replenished during periods of summer stagnation, and water below depths of 30 to 40 feet usually contains less than 2 milligrams per liter of dissolved oxygen during stagnation. Even though heated effluent from an electrical generating plant is returned to the reservoir at site GC and causes an elevation of water temperature, average dissolved oxygen levels at this site are not significantly different from levels at other sites. During summer stagnation, reducing conditions result in the dissolution of iron and manganese from bottom deposits at deep sites in the reservoir. At site AQ, a deep site near Arlington Dam, dissolved iron concentrations in water near the bottom during summer stagnation ranged from less than 10 to 1,100 micrograms per liter and averaged about 640 micrograms per liter. Dissolved manganese concentrations near the bottom at site AC during summer staynation ranged from 20 to 2,700 micrograms per liter and averaged about 1,500 micrograms per liter. The concentrations of dissolved iron and dissolved manganese in water throughout the reservoir during winter circulation and in water near the reservoir surface during summer stagnation averaged less than 50 micrograms per liter. Seasonal temperature and dissolved oxygen cycles resulted in the recycling of dissolved iron and dissolved manganese between the water and bottom sediments. However, no significant accumulation of these constituents within the reservoir was detected during the study. The concentrations of total inorganic nitrogen and total phosphorus are greatest during summer stagnation in water near the bottom at deep sites. At site AC during the summer, the concentrations of total inorganic nitrogen in the hypolimnion averaged about 0.9 milligram per liter, and the concentration of total phosphorus near the bottom averaged about 0.2 milligram per liter. The concentrations of total inorganic nitrogen in the epilimnion at site AC averaged about 0.1 milligram per liter; the concentrations of total phosphorus averaged less than 0.1 milligram per liter. The densities and composition of algal populations varied seasonally. At site AC, total algae counts ranged from 200 to 240,000 cells per mi Hi liter and averaged about 50,000 cells per milliliter. At site FC, total algae counts ranged from 1,000 to 290,000 cells per milliliter and averaged about 56,000 cells per milliliter. Algal densities were greatest during the summer with blue-green algae being the predominant phyla.

Texas

A computer program for converting rectangular coordinates to latitude-longitude coordinates

A computer program was developed for converting the coordinates of any rectangular grid on a map to coordinates on a grid that is parallel to lines of equal latitude and longitude. Using this program in conjunction with groundwater flow models, the user can extract data and results from models with varying grid orientations and place these data into grid structure that is oriented parallel to lines of equal latitude and longitude. All cells in the rectangular grid must have equal dimensions, and all cells in the latitude-longitude grid measure one minute by one minute. This program is applicable if the map used shows lines of equal latitude as arcs and lines of equal longitude as straight lines and assumes that the Earth 's surface can be approximated as a sphere. The program user enters the row number , column number, and latitude and longitude of the midpoint of the cell for three test cells on the rectangular grid. The latitude and longitude of boundaries of the rectangular grid also are entered. By solving sets of simultaneous linear equations, the program calculates coefficients that are used for making the conversion. As an option in the program, the user may build a groundwater model file based on a grid that is parallel to lines of equal latitude and longitude. The program reads a data file based on the rectangular coordinates and automatically forms the new data file. (USGS)

Water-Resources Investigations Report

Water-quality conditions at selected landfills in Mecklenburg County, North Carolina, 1986-92

Water-quality conditions at five municipal landfills in Mecklenburg County, North Carolina, were studied during 1986-92. Analytical results of water samples from monitoring wells and streams at and near the landfills were used to evaluate effects of leachate on surface and ground water. Ground-water levels at monitoring wells were used to determine directions of ground-water flow at the landfills. Data from previous studies were used for analysis of temporal trends in selected water-quality properties and chemical constituents. Effects of leachate, such as large biochemical- and chemical-oxygen demands, generally were evident in small streams originating within the landfills, whereas effects of leachate generally were not evident in most of the larger streams. In larger streams, surface-water quality upstream and downstream from most of the landfills was similar. However, the chemical quality of water in Irwin Creek appears to have been affected by the Statesville Road landfill. Concentrations of several constituents indicative of leachate were larger in samples collected from Irwin Creek downstream from the Statesville Road landfill than in samples collected from Irwin Creek upstream from the landfill. The effect of leachate on ground-water quality generally was largest in water from wells adjacent to waste-disposal cells. Concentrations of most constituents considered indicative of leachate generally were smaller with increasing distance from waste-disposal cells. Water samples from offsite wells generally indicated no effect or very small effects of leachate. Action levels designated by the Mecklenburg County Engineering Department and maximum contaminant levels established by the U.S. Environmental Protection Agency were exceeded in some samples from the landfills. Ground-water samples exceeded action levels and maximum contaminant levels more commonly than surface-water samples. Iron and manganese were the constituents that most commonly exceeded action levels in water samples from the landfills. Synthetic organic compounds were detected more commonly and in larger concentrations in ground-water samples than in surface-water samples. Concentrations of synthetic organic compounds detected in water samples from monitoring sites at the landfills generally were much less than maximum contaminant levels. However, concentrations of some chlorinated organic compounds exceeded maximum contaminant levels in samples from several monitoring wells at the Harrisburg Road and York Road landfills. Trend analysis indicated statistically significant temporal changes in concentrations of selected water-quality constituents and properties at some of the monitoring sites. Trends detected for the Holbrooks Road and Statesville Road landfills generally indicated an improvement in water quality and a decrease in effects of leachate at most monitoring sites at these landfills from 1979 to 1992. Water-quality trends detected for monitoring sites at the Harrisburg Road and York Road landfills, the largest landfills in the study, differed in magnitude and direction. Upward trends generally were detected for sites near recently closed waste-disposal cells, whereas downward trends generally were detected for sites near older waste-disposal cells. Temporal trends in water quality generally reflected changes in degradation processes associated with the aging of landfill wastes.

North Carolina

Ground-water hydrology, historical water use, and simulated ground-water flow in Cretaceous-age Coastal Plain aquifers near Charleston and Florence, South Carolina

A quasi-three-dimensional, transient, digital, ground-water flow model representing the Coastal Plain aquifers of South Carolina, has been constructed to assist in defining the ground- water-flow system of Cretaceous aquifers near Charleston and Florence, S.C. Both cities are near the centers of large (greater than 150 feet) potentiometric declines in the Middendorf aquifer. In 1989, the diameter of the depressions was approximately 40 miles at Charleston and 15 miles at Florence. The potentiometric decline occurred between predevelopment (1926) and 1982 near Florence, and between predevelopment (1879) and 1989 near Charleston. The city of Charleston does not withdraw water from these aquifers; however, some of the small communities in the area use these aquifers for a potable water supply. The model simulates flow in and between four aquifer systems. The model has a variable-cell-size grid, and spans the Coastal Plain from the Savannah River in the southwest to the Cape Fear Arch in the northeast, and from the Fall Line in the northwest to approximately 30 miles offshore to the southeast. Model-grid cell size is 1 by 1 mile in a 48 by 48 mile area centered in Charleston, and in a 36 by 48 mile area centered in Florence. The model cell size gradually increases to a maximum of 4 by 4 miles outside the two study areas. The entire grid consists of 115 by 127 cells and covers an area of 39,936 square miles. The model was calibrated to historical water-level data. The calibration relied on three techniques: (1) matching simulated and observed potentiometric map surfaces, (2) statistical comparison of observed and simulated heads, and (3) comparison of observed and simulated well hydrographs. Systematic changes in model parameters showed that simulated heads are most sensitive to changes in aquifer transmissivity. Eight predictive ground-water-use scenarios were simulated for the Mount Pleasant area, which presently (1993) uses the Middendorf aquifer as a sole-source of potable water. These simulations use various combinations of spatial distribution, and injection of treated wastewater effluent for existing and future Middendorf aquifer wells.

South Carolina