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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47Linked to original sources

USGS2024

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗

USGS2026

Phytoplankton responses to experimental nitrogen and phosphorus loading in the eutrophic and colored Caloosahatchee River, Florida

The Caloosahatchee River, located in southwest Florida, is a eutrophic and colored river that flows from Lake Okeechobee westward into its estuary and the Gulf of America. Cyanobacterial harmful algal blooms (HABs) are a documented problem along this freshwater-to-marine waterway where nutrient enrichment has been identified as a key factor in bloom occurrence but has not been experimentally tested in the river. This study is the first to test the effects of inorganic nutrient loading on phytoplankton assemblages in the Caloosahatchee River and the effects of different nutrient sources on phytoplankton dynamics at different times of the year. Three independent, in situ experiments were conducted to test the effects of daily, incrementally increased ammonium, nitrate, and phosphate loading on phytoplankton at different times of the year (summer, fall, winter). Over the 72-hour enclosure period, phytoplankton abundance metrics (cell concentration, chlorophyll- a , and phycocyanin), dissolved oxygen, and pH increased, and fluorescent dissolved organic matter and turbidity decreased in all treatments and controls. Increased phytoplankton abundance metrics relative to controls were observed after 72 hours of exposure to elevated ammonium and nitrate in summer and only ammonium in winter, suggesting periodic nitrogen limitation; however, no treatment effects on phytoplankton assemblage structure in terms of resemblance and diversity metrics were found. Increases in total cell concentrations were driven by elevated growth rates of already dominant taxa but not sufficiently to form a visible bloom. Cyanobacteria consistently dominated the phytoplankton, particularly Aphanocapsa and Merismopedia , whereas the common HAB-forming Microcystis maintained consistently low abundance. This study provides new information on the ecology of phytoplankton assemblages in the Caloosahatchee River and could be used by water resources managers to evaluate strategies for controlling cyanobacterial HABs in the river.

Florida↗

USGS1956

Duttonite, a new quadrivalent vanadium oxide from the Peanut mine, Montrose County, Colorado

Duttonite, a new quadrivalent vanadium oxide from the Peanut mine, Montrose County, Colo., has the formula VO(OH) 2 . The mineral occurs as crusts and coatings of pale-brown transparent platy crystals, as one of the first oxidation products of montroseite ore. It is associated with melanovanadite and abundant crystals of hexagonal native selenium. Duttonite is biaxial positive, 2V is about 60°, dispersion is r < v , moderate; X = a , pale pinkish brown; Y = c , pale yellow-brown; Z = b , pale brown; α = 1.810 ± 0.003, β = 1.900 ± 0.003, γ > 2.01. The hardness is about 2.5; the calculated specific gravity is 3.24. The chemical analysis shows, in percent: V 2 o 3 2.6, V 2 O 4 75.3, FeO 0.4, H 2 O 18.1, insoluble 4.2, total 100.6. Duttonite is monoclinic, a o = 8.80 ± 0.02A, b o - 3.95 ± 0.01A, c o - 5.96 ± 0.02A, β = 90°40 1 ± 5 1 . The space group is I 2/ c , ( C 6 2 h ); the cell contents are 4[VO(OH) 2 ]. The crystals are strongly pseudo-orthorhombic, and the structure departs only slightly from the space group Imcm . Duttonite is named for Captain Clarence Edward Dutton (1841-1912). A detailed study of the geology, geochemistry, and mineralogy of the vanadium-uranium ore at the Peanut mine, Montrose County, Colo., was begun early in 1954 by Carl H. Roach of the U. S. Geological Survey. A number of rare and new minerals were found in the ore and the study of these samples was undertaken by Mary E. Thompson. Duttonite is the first new vanadium mineral to be described from the Peanut mine. It is named for Captain Clarence Edward Dutton (1841-1912), who was one of the first geologists to work in the Colorado Plateau region and who was a member of the U. s. Geological Survey from 1879-91. We are indebted to the following members of the Geological Surbey: K. E. Valentine for spectrographic analyses of duttonite, and M. E. Mrose and H. T. Evans, Jr., for measurement of the unit cell constants. This work is part of a program being conducted by the U. S. Geological Survey on behalf of the Division of Raw Materials of the U. S. Atomic Energy Commission.

Colorado↗

USGS2006

User guide for the farm process (FMP1) for the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model, MODFLOW-2000

There is a need to estimate dynamically integrated supply-and-demand components of irrigated agriculture as part of the simulation of surface-water and ground-water flow. To meet this need, a computer program called the Farm Process (FMP1) was developed for the U.S. Geological Survey three-dimensional finite-difference modular ground-water flow model, MODFLOW- 2000 (MF2K). The FMP1 allows MF2K users to simulate conjunctive use of surface- and ground water for irrigated agriculture for historical and future simulations, water-rights issues and operational decisions, nondrought and drought scenarios. By dynamically integrating farm delivery requirement, surface- and ground-water delivery, as well as irrigation-return flow, the FMP1 allows for the estimation of supplemental well pumpage. While farm delivery requirement and irrigation return flow are simulated by the FMP1, the surface-water delivery to the farm can be simulated optionally by coupling the FMP1 with the Streamflow Routing Package (SFR1) and the farm well pumping can be simulated optionally by coupling the FMP1 to the Multi-Node Well (MNW) Package. In addition, semi-routed deliveries can be specified that are associated with points of diversion in the SFR1 stream network. Nonrouted surface-water deliveries can be specified independently of any stream network. The FMP1 maintains a dual mass balance of a farm budget and as part of the ground-water budget. Irrigation demand, supply, and return flow are in part subject to head-dependent sources and sinks such as evapotranspiration from ground water and leakage between the conveyance system and the aquifer. Farm well discharge and farm net recharge are source/sink terms in the FMP1, which depend on transpiration uptake from ground water and other head dependent consumptive use components. For heads rising above the bottom of the root zone, the actual transpiration is taken to vary proportionally with the depth of the active root zone, which can be restricted by anoxia or wilting. Depths corresponding to anoxia- or wilting-related pressure heads within the root zone are found using analytical solutions of a vertical pseudo steady-state pressure- head distribution over the depth of the total root zone (Consumptive Use Concept 1). Alternatively, a simpler, conceptual model is available, which defines how consumptive use (CU) components vary with changing head (CU Concept 2). Subtracting the ground water and precipitation transpiration components from the total transpiration yields a transpiratory irrigation requirement for each cell. The total farm delivery requirement (TFDR) then is determined as cumulative transpiratory and evaporative irrigation requirements of all farm cells and increased sufficiently to compensate for inefficient use from irrigation with respect to plant consumption. The TFDR subsequently is satisfied with surface- and ground-water delivery, respectively constrained by allotments, water rights, or maximum capacities. Five economic and noneconomic drought response policies can be applied optionally, if the potential supply of surface water and ground water is insufficient to meet the crop demand: acreage-optimization with or without a water conservation pool, deficit irrigation with or without water-stacking, and zero policy.

Techniques and Methods↗

USGS2011

MODFLOW-NWT, a Newton formulation for MODFLOW-2005

This report documents a Newton formulation of MODFLOW-2005, called MODFLOW-NWT. MODFLOW-NWT is a standalone program that is intended for solving problems involving drying and rewetting nonlinearities of the unconfined groundwater-flow equation. MODFLOW-NWT must be used with the Upstream-Weighting (UPW) Package for calculating intercell conductances in a different manner than is done in the Block-Centered Flow (BCF), Layer Property Flow (LPF), or Hydrogeologic-Unit Flow (HUF; Anderman and Hill, 2000) Packages. The UPW Package treats nonlinearities of cell drying and rewetting by use of a continuous function of groundwater head, rather than the discrete approach of drying and rewetting that is used by the BCF, LPF, and HUF Packages. This further enables application of the Newton formulation for unconfined groundwater-flow problems because conductance derivatives required by the Newton method are smooth over the full range of head for a model cell. The NWT linearization approach generates an asymmetric matrix, which is different from the standard MODFLOW formulation that generates a symmetric matrix. Because all linear solvers presently available for use with MODFLOW-2005 solve only symmetric matrices, MODFLOW-NWT includes two previously developed asymmetric matrix-solver options. The matrix-solver options include a generalized-minimum-residual (GMRES) Solver and an Orthomin / stabilized conjugate-gradient (CGSTAB) Solver. The GMRES Solver is documented in a previously published report, such that only a brief description and input instructions are provided in this report. However, the CGSTAB Solver (called X MD) is documented in this report. Flow-property input for the UPW Package is designed based on the LPF Package and material-property input is identical to that for the LPF Package except that the rewetting and vertical-conductance correction options of the LPF Package are not available with the UPW Package. Input files constructed for the LPF Package can be used with slight modification as input for the UPW Package. This report presents the theory and methods used by MODFLOW-NWT, including the UPW Package. Additionally, this report provides comparisons of the new methodology to analytical solutions of groundwater flow and to standard MODFLOW-2005 results by use of an unconfined aquifer MODFLOW example problem. The standard MODFLOW-2005 simulation uses the LPF Package with the wet/dry option active. A new example problem also is presented to demonstrate MODFLOW-NWT's ability to provide a solution for a difficult unconfined groundwater-flow problem.

Techniques and Methods↗

USGS2012

Documentation of the Surface-Water Routing (SWR1) Process for modeling surface-water flow with the U.S. Geological Survey Modular Ground-Water Model (MODFLOW-2005)

A flexible Surface-Water Routing (SWR1) Process that solves the continuity equation for one-dimensional and two-dimensional surface-water flow routing has been developed for the U.S. Geological Survey three-dimensional groundwater model, MODFLOW-2005. Simple level- and tilted-pool reservoir routing and a diffusive-wave approximation of the Saint-Venant equations have been implemented. Both methods can be implemented in the same model and the solution method can be simplified to represent constant-stage elements that are functionally equivalent to the standard MODFLOW River or Drain Package boundary conditions. A generic approach has been used to represent surface-water features (reaches) and allows implementation of a variety of geometric forms. One-dimensional geometric forms include rectangular, trapezoidal, and irregular cross section reaches to simulate one-dimensional surface-water features, such as canals and streams. Two-dimensional geometric forms include reaches defined using specified stage-volume-area-perimeter (SVAP) tables and reaches covering entire finite-difference grid cells to simulate two-dimensional surface-water features, such as wetlands and lakes. Specified SVAP tables can be used to represent reaches that are smaller than the finite-difference grid cell (for example, isolated lakes), or reaches that cannot be represented accurately using the defined top of the model. Specified lateral flows (which can represent point and distributed flows) and stage-dependent rainfall and evaporation can be applied to each reach. The SWR1 Process can be used with the MODFLOW Unsaturated Zone Flow (UZF1) Package to permit dynamic simulation of runoff from the land surface to specified reaches. Surface-water/groundwater interactions in the SWR1 Process are mathematically defined to be a function of the difference between simulated stages and groundwater levels, and the specific form of the reach conductance equation used in each reach. Conductance can be specified directly or calculated as a function of the simulated wetted perimeter and defined reach bed hydraulic properties, or as a weighted combination of both reach bed hydraulic properties and horizontal hydraulic conductivity. Each reach can be explicitly coupled to a single specific groundwater-model layer or coupled to multiple groundwater-model layers based on the reach geometry and groundwater-model layer elevations in the row and column containing the reach. Surface-water flow between reservoirs is simulated using control structures. Surface-water flow between reaches, simulated by the diffusive-wave approximation, can also be simulated using control structures. A variety of control structures have been included in the SWR1 Process and include (1) excess-volume structures, (2) uncontrolled-discharge structures, (3) pumps, (4) defined stage-discharge relations, (5) culverts, (6) fixed- or movable-crest weirs, and (7) fixed or operable gated spillways. Multiple control structures can be implemented in individual reaches and are treated as composite flow structures. Solution of the continuity equation at the reach-group scale (a single reach or a user-defined collection of individual reaches) is achieved using exact Newton methods with direct solution methods or exact and inexact Newton methods with Krylov sub-space methods. Newton methods have been used in the SWR1 Process because of their ability to solve nonlinear problems. Multiple SWR1 time steps can be simulated for each MODFLOW time step, and a simple adaptive time-step algorithm, based on user-specified rainfall, stage, flow, or convergence constraints, has been implemented to better resolve surface-water response. A simple linear- or sigmoid-depth scaling approach also has been implemented to account for increased bed roughness at small surface-water depths and to increase numerical stability. A line-search algorithm also has been included to improve the quality of the Newton-step upgrade vector, if possible. The SWR1 Process has been benchmarked against one- and two-dimensional numerical solutions from existing one- and two-dimensional numerical codes that solve the dynamic-wave approximation of the Saint-Venant equations. Two-dimensional solutions test the ability of the SWR1 Process to simulate the response of a surface-water system to (1) steady flow conditions for an inclined surface (solution of Manning's equation), and (2) transient inflow and rainfall for an inclined surface. The one-dimensional solution tests the ability of the SWR1 Process to simulate a looped network with multiple upstream inflows and several control structures. The SWR1 Process also has been compared to a level-pool reservoir solution. A synthetic test problem was developed to evaluate a number of different SWR1 solution options and simulate surface-water/groundwater interaction. The solution approach used in the SWR1 Process may not be applicable for all surface-water/groundwater problems. The SWR1 Process is best suited for modeling long-term changes (days to years) in surface-water and groundwater flow. Use of the SWR1 Process is not recommended for modeling the transient exchange of water between streams and aquifers when local and convective acceleration and other secondary effects (for example, wind and Coriolis forces) are substantial. Dam break evaluations and two-dimensional evaluations of spatially extensive domains are examples where acceleration terms and secondary effects would be significant, respectively.

Techniques and Methods↗

USGS2012

User guide for MODPATH version 6—A particle-tracking model for MODFLOW

MODPATH is a particle-tracking post-processing model that computes three-dimensional flow paths using output from groundwater flow simulations based on MODFLOW, the U.S. Geological Survey (USGS) finite-difference groundwater flow model. This report documents MODPATH version 6. Previous versions were documented in USGS Open-File Reports 89-381 and 94-464. The program uses a semianalytical particle-tracking scheme that allows an analytical expression of a particle's flow path to be obtained within each finite-difference grid cell. A particle's path is computed by tracking the particle from one cell to the next until it reaches a boundary, an internal sink/source, or satisfies another termination criterion. Data input to MODPATH consists of a combination of MODFLOW input data files, MODFLOW head and flow output files, and other input files specific to MODPATH. Output from MODPATH consists of several output files, including a number of particle coordinate output files intended to serve as input data for other programs that process, analyze, and display the results in various ways. MODPATH is written in FORTRAN and can be compiled by any FORTRAN compiler that fully supports FORTRAN-2003 or by most commercially available FORTRAN-95 compilers that support the major FORTRAN-2003 language extensions.

Techniques and Methods↗

USGS2014

Two graphical user interfaces for managing and analyzing MODFLOW groundwater-model scenarios

Scenario Manager and Scenario Analyzer are graphical user interfaces that facilitate the use of calibrated, MODFLOW-based groundwater models for investigating possible responses to proposed stresses on a groundwater system. Scenario Manager allows a user, starting with a calibrated model, to design and run model scenarios by adding or modifying stresses simulated by the model. Scenario Analyzer facilitates the process of extracting data from model output and preparing such display elements as maps, charts, and tables. Both programs are designed for users who are familiar with the science on which groundwater modeling is based but who may not have a groundwater modeler’s expertise in building and calibrating a groundwater model from start to finish. With Scenario Manager, the user can manipulate model input to simulate withdrawal or injection wells, time-variant specified hydraulic heads, recharge, and such surface-water features as rivers and canals. Input for stresses to be simulated comes from user-provided geographic information system files and time-series data files. A Scenario Manager project can contain multiple scenarios and is self-documenting. Scenario Analyzer can be used to analyze output from any MODFLOW-based model; it is not limited to use with scenarios generated by Scenario Manager. Model-simulated values of hydraulic head, drawdown, solute concentration, and cell-by-cell flow rates can be presented in display elements. Map data can be represented as lines of equal value (contours) or as a gradated color fill. Charts and tables display time-series data obtained from output generated by a transient-state model run or from user-provided text files of time-series data. A display element can be based entirely on output of a single model run, or, to facilitate comparison of results of multiple scenarios, an element can be based on output from multiple model runs. Scenario Analyzer can export display elements and supporting metadata as a Portable Document Format file.

Techniques and Methods↗

USGS2016

MT3D-USGS version 1: A U.S. Geological Survey release of MT3DMS updated with new and expanded transport capabilities for use with MODFLOW

MT3D-USGS, a U.S. Geological Survey updated release of the groundwater solute transport code MT3DMS, includes new transport modeling capabilities to accommodate flow terms calculated by MODFLOW packages that were previously unsupported by MT3DMS and to provide greater flexibility in the simulation of solute transport and reactive solute transport. Unsaturated-zone transport and transport within streams and lakes, including solute exchange with connected groundwater, are among the new capabilities included in the MT3D-USGS code. MT3D-USGS also includes the capability to route a solute through dry cells that may occur in the Newton-Raphson formulation of MODFLOW (that is, MODFLOW-NWT). New chemical reaction Package options include the ability to simulate inter-species reactions and parent-daughter chain reactions. A new pump-and-treat recirculation package enables the simulation of dynamic recirculation with or without treatment for combinations of wells that are represented in the flow model, mimicking the above-ground treatment of extracted water. A reformulation of the treatment of transient mass storage improves conservation of mass and yields solutions for better agreement with analytical benchmarks. Several additional features of MT3D-USGS are (1) the separate specification of the partitioning coefficient ( K d ) within mobile and immobile domains; (2) the capability to assign prescribed concentrations to the top-most active layer; (3) the change in mass storage owing to the change in water volume now appears as its own budget item in the global mass balance summary; (4) the ability to ignore cross-dispersion terms; (5) the definition of Hydrocarbon Spill-Source Package (HSS) mass loading zones using regular and irregular polygons, in addition to the currently supported circular zones; and (6) the ability to specify an absolute minimum thickness rather than the default percent minimum thickness in dry-cell circumstances. Benchmark problems that implement the new features and packages test the accuracy of new code through comparison to analytical benchmarks, as well as to solutions from other published codes. The input file structure for MT3D-USGS adheres to MT3DMS conventions for backward compatibility: the new capabilities and packages described herein are readily invoked by adding three-letter package name acronyms to the name file or by setting input flags as needed. Memory is managed in MT3D-USGS using FORTRAN modules in order to simplify code development and expansion.

Techniques and Methods↗

USGS2016

T-COMP — A suite of programs for extracting transmissivity from MODFLOW models

Simulated transmissivities are constrained poorly by assigning permissible ranges of hydraulic conductivities from aquifer-test results to hydrogeologic units in groundwater-flow models. These wide ranges are derived from interpretations of many aquifer tests that are categorized by hydrogeologic unit. Uncertainty is added where contributing thicknesses differ between field estimates and numerical models. Wide ranges of hydraulic conductivities and discordant thicknesses result in simulated transmissivities that frequently are much greater than aquifer-test results. Multiple orders of magnitude differences frequently occur between simulated and observed transmissivities where observed transmissivities are less than 1,000 feet squared per day. Transmissivity observations from individual aquifer tests can constrain model calibration as head and flow observations do. This approach is superior to diluting aquifer-test results into generalized ranges of hydraulic conductivities. Observed and simulated transmissivities can be compared directly with T-COMP, a suite of three FORTRAN programs. Transmissivity observations require that simulated hydraulic conductivities and thicknesses in the volume investigated by an aquifer test be extracted and integrated into a simulated transmissivity. Transmissivities of MODFLOW model cells are sampled within the volume affected by an aquifer test as defined by a well-specific, radial-flow model of each aquifer test. Sampled transmissivities of model cells are averaged within a layer and summed across layers. Accuracy of the approach was tested with hypothetical, multiple-aquifer models where specified transmissivities ranged between 250 and 20,000 feet squared per day. More than 90 percent of simulated transmissivities were within a factor of 2 of specified transmissivities.

Nevada↗

USGS2021

Probabilistic methodology for the assessment of original and recoverable coal resources, illustrated with an application to a coal bed in the Fort Union Formation, Wyoming

Executive Summary The U.S. Geological Survey (USGS) has been using its Circular 891 for evaluating uncertainty in coal resource assessments for more than 35 years. Calculated cell tonnages are assigned to four qualitative reliability classes depending exclusively on distance to the nearest drill hole. The main appeal of this methodology, simplicity, is also its main drawback. Reliability may depend so marginally on distance to the nearest drill hole that, over time, it has become evident that Circular 891 is inadequate for modeling reliability and is limited by other shortcomings. The present publication describes the use of geostatistics as an approach allowing a more satisfactory performance than that which is achieved following Circular 891. Geostatistics takes advantage of partly random and partly organized fluctuations in attributes such as coal thickness, coal density, and elevation of the top of a coal bed, borrowing concepts and tools that have been standard features in statistics and risk analysis for decades. Considering that readers interested in this study may not have the background to go directly into the details of the methodology, we start by explaining geostatistical concepts and modeling techniques. The remainder of the publication is devoted to formulating the assessment methodology, applying it to data from the Fillmore Ranch coal bed in the Fort Union Formation in Wyoming, and explaining the computer software applied for performing calculations and displays. The assessment methodology has been designed to report three different forms of resources: coal in place, coal mineable by surface mining methods, and coal mineable by underground mining methods. These three types of resources are reported graphically by displaying both the magnitude and the reliability of total coal resources and resources at the cell scale. In the case of the Fillmore Ranch coal bed example, there is a 90-percent probability that the resources in place are 9.687 ± 0.383 billion short tons (bst), while the coal available for underground mining is 2.279 ± 0.160 bst, and that available for surface mining is only 0.240 ± 0.025 bst because of the steep dip to the west away from the outcrop. These magnitudes are derived from numerical probability distributions not following any specific form.

Wyoming↗

USGS2000

Assessment of freshwater mussels in the Allegheny River at Foxburg, Pennsylvania, 1998

The upper reaches of the Allegheny River are a high-quality resource that supports populations of a number of endangered species. Two endangered species of freshwater mussel, the northern riffleshell, Epioblasma torulosa rangiana , and clubshell, Pleurobema clava , are present in this river reach. Prior to a bridge-replacement project at the Allegheny River at Foxburg, Pa. (river mile 86.2), a mussel survey was conducted to determine if either of the protected endangered species are present and to assess the quality of the mussel bed at this site. Shore and near-shore searches followed by SCUBA diving transects determined that a mussel bed is present at this location. The freshwater-mussel community consists of at least 14 species and includes both endangered species. Substrate type, stream velocity, and channel morphology combined with results from the transect searches documented that the mussel bed is largely limited to the east side of the river. A two-stage sampling design was employed to estimate the overall mussel abundance within the bed. Twelve 4-m 2 (square meter) cells were sampled with four 0.25-m 2 quadrants per cell. The mean population density within the mussel bed is 8.4 mussel per square meter. The overall mussel population in the survey area is 225,567 individuals (95-percent confidence interval, 135,973 to 374,195). River-bed scour and acid mine drainage appear to limit mussel distribution in the center and left section of the channel. Overall, a healthy and diverse mussel bed exists at this site.

Water-Resources Investigations Report↗

USGS2000

Characterization and simulation of ground-water flow in the Kansas River Valley at Fort Riley, Kansas, 1990-98

Hydrologic data and a ground-water flow model were used to characterize ground-water flow in the Kansas River alluvial aquifer at Fort Riley in northeast Kansas. The ground-water flow model was developed as a tool to project ground-water flow and potential contaminant-transport paths in the alluvial aquifer on the basis of past hydrologic conditions. The model also was used to estimate historical and hypothetical ground-water flow paths with respect to a private- and several public-supply wells. The ground-water flow model area extends from the Smoky Hill and Republican Rivers downstream to about 2.5 miles downstream from the city of Ogden. The Kansas River Valley has low relief and, except for the area within the Fort Riley Military Reservation, is used primarily for crop production. Sedimentary deposits in the Kansas River Valley, formed after the ancestral Kansas River eroded into bedrock, primarily are alluvial sediment deposited by the river during Quaternary time. The alluvial sediment consists of as much as about 75 feet of poorly sorted, coarse-to-fine sand, silt, and clay, 55 feet of which can be saturated with ground water. The alluvial aquifer is unconfined and is bounded on the sides and bottom by Permian-age shale and limestone bedrock. Hydrologic data indicate that ground water in the Kansas River Valley generally flows in a downstream direction, but flow direction can be quite variable near the Kansas River due to changes in river stage. Ground-water-level changes caused by infiltration of precipitation are difficult to detect because they are masked by larger changes caused by fluctuation in Kansas River stage. Ratios of strontium isotopes Sr87 and Sr86 in water collected from wells in the Camp Funston Area indicate that the ground water along the northern valley wall originates, in part, from upland areas north of the river valley. Water from Threemile Creek, which flows out of the uplands north of the river valley, had Sr87:Sr86 ratios similar to those in ground water from wells in the northern Camp Funston Area. In addition, comparison of observed water levels from wells CF90-06, CF97-101, and CF97-401 in the Camp Funston Area and ground-water levels simulated for these wells using floodwave-response analysis indicates that ground-water inflow from bedrock is a hydraulic stress that, in addition to the changing stage in the Kansas River, acts on the aquifer. This hydraulic stress seems to be located near the northern valley wall because the effect of this stress is greater for well CF97-101, which is the well closest to the valley wall. Ground-water flow was simulated using a modular, three-dimensional, finite-difference ground-water flow model (MODFLOW). Particle tracking, used to visualize ground-water flow paths in the alluvial aquifer, was accomplished using MODPATH. Forward-in-time particle tracking indicated that, in general, particles released near the Kansas River followed much more variable paths than particles released near the valley wall. Although particle tracking does not simulate solute transport, this increased path variability indicates that, near the river, ground-water contaminants could follow many possible paths towards the river, whereas more distant from the river, ground-water contaminants likely would follow a narrower corridor. Particle tracks in the Camp Funston Area indicate that, for the 1990-98 simulation period, contaminants from the ground-water study sites in the Camp Funston Area would be unlikely to move into the vicinity of Ogden's supply wells. Backward-in-time particle tracking indicated that the flow-path and recharge areas for model cells corresponding to Ogden's supply wells lie near the northern valley wall and extend into the northern Camp Funston Area. The flow-path and recharge areas for model cells corresponding to Morris County Rural Water District wells lie within Clarks Creek Valley and probably extend outside the model area. Three hypothetical simulations, i

Water-Resources Investigations Report↗

USGS2000

Quality-assurance design applied to an assessment of agricultural pesticides in ground water from carbonate bedrock aquifers in the Great Valley of eastern Pennsylvania

Assessments to determine whether agricultural pesticides are present in ground water are performed by the Commonwealth of Pennsylvania under the aquifer monitoring provisions of the State Pesticides and Ground Water Strategy. Pennsylvania's Department of Agriculture conducts the monitoring and collects samples; the Department of Environmental Protection (PaDEP) Laboratory analyzes the samples to measure pesticide concentration. To evaluate the quality of the measurements of pesticide concentration for a groundwater assessment, a quality-assurance design was developed and applied to a selected assessment area in Pennsylvania. This report describes the quality-assurance design, describes how and where the design was applied, describes procedures used to collect and analyze samples and to evaluate the results, and summarizes the quality assurance results along with the assessment results. The design was applied in an agricultural area of the Delaware River Basin in Berks, Lebanon, Lehigh, and Northampton Counties to evaluate the bias and variability in laboratory results for pesticides. The design—with random spatial and temporal components—included four data-quality objectives for bias and variability. The spatial design was primary and represented an area comprising 30 sampling cells. A quality-assurance sampling frequency of 20 percent of cells was selected to ensure a sample number of five or more for analysis. Quality-control samples included blanks, spikes, and replicates of laboratory water and spikes, replicates, and 2-lab splits of groundwater. Two analytical laboratories, the PaDEP Laboratory and a U.S. Geological Survey Laboratory, were part of the design. Bias and variability were evaluated by use of data collected from October 1997 through January 1998 for alachlor, atrazine, cyanazine, metolachlor, simazine, pendimethalin, metribuzin, and chlorpyrifos. Results of analyses of field blanks indicate that collection, processing, transport, and laboratory analysis procedures did not contaminate the samples; there were no false-positive results. Pesticides were detected in water when pesticides were spiked into (added to) samples. There were no false negatives for the eight pesticides in all spiked samples. Negative bias was characteristic of analytical results for the eight pesticides, and bias was generally in excess of 10 percent from the ‘true’ or expected concentration (34 of 39 analyses, or 87 percent of the ground-water results) for pesticide concentrations ranging from 0.31 to 0.51 mg/L (micrograms per liter). The magnitude of the negative bias for the eight pesticides, with the exception of cyanazine, would result in reported concentrations commonly 75-80 percent of the expected concentration in the water sample. The bias for cyanazine was negative and within 10 percent of the expected concentration. A comparison of spiked pesticide-concentration recoveries in laboratory water and ground water indicated no effect of the ground-water matrix, and matrix interference was not a source of the negative bias. Results for the laboratory-water spikes submitted in triplicate showed large variability for recoveries of atrazine, cyanazine, and pendimethalin. The relative standard deviation (RSD) was used as a measure of method variability over the course of the study for laboratory waters at a concentration of 0.4 mg/L. An RSD of about 11 percent (or about ?0.05 mg/L)characterizes the method results for alachlor, chlorpyrifos, metolachlor, metribuzin, and simazine. Atrazine and pendimethalin have RSD values of about 17 and 23 percent, respectively. Cyanazine showed the largest RSD at nearly 51 percent. The pesticides with low variability in laboratory-water spikes also had low variability in ground water. The assessment results showed that atrazinewas the most commonly detected pesticide in ground water in the assessment area. Atrazine was detected in water from 22 of the 28 wells sampled, and recovery results for atrazine were some of the worst (largest negative bias). Concentrations of the eight pesticides in ground water from wells were generally less than 0.3 µg/L. Only six individual measurements of the concentrations in water from six of the wells were at or above 0.3 µg/L, ﬁve for atrazine and one for metolachlor. There were eight additional detections of metolachlor and simazine at concentrations less than 0.1 µg/L. No well water contained more than one pesticide at concentra-tions at or above 0.3 µg/L. Evidence exists, how-ever, for a pattern of co-occurrence of metolachlor and simazine at low concentrations with higher concentrations of atrazine. Large variability in replicate samples and negative bias for pesticide recovery from spiked samples indicate the need to use data for pesticide recovery in the interpretation of measured pesti-cide concentrations in ground water. Data from samples spiked with known amounts of pesticides were a critical component of a quality-assurance design for the monitoring component of the Pesti-cides and Ground Water Strategy. Trigger concentrations, the concentrations that require action under the Pesticides and Ground Water Strategy, should be considered maximums for action. This consideration is needed because of the magnitude of negative bias.

Water-Resources Investigations Report↗

USGS2003

Analytical and numerical simulation of the steady-state hydrologic effects of mining aggregate in hypothetical sand-and-gravel and fractured crystalline-rock aquifers

Analytical solutions and numerical models were used to predict the extent of steady-state drawdown caused by mining of aggregate below the water table in hypothetical sand-and-gravel and fractured crystalline-rock aquifers representative of hydrogeologic settings in the Front Range area of Colorado. Analytical solutions were used to predict the extent of drawdown under a wide range of hydrologic and mining conditions that assume aquifer homogeneity, isotropy, and infinite extent. Numerical ground-water flow models were used to estimate the extent of drawdown under conditions that consider heterogeneity, anisotropy, and hydrologic boundaries and to simulate complex or unusual conditions not readily simulated using analytical solutions. Analytical simulations indicated that the drawdown radius (or distance) of influence increased as horizontal hydraulic conductivity of the aquifer, mine penetration of the water table, and mine radius increased; radius of influence decreased as aquifer recharge increased. Sensitivity analysis of analytical simulations under intermediate conditions in sand-and-gravel and fractured crystalline-rock aquifers indicated that the drawdown radius of influence was most sensitive to mine penetration of the water table and least sensitive to mine radius. Radius of influence was equally sensitive to changes in horizontal hydraulic conductivity and recharge. Numerical simulations of pits in sand-and- gravel aquifers indicated that the area of influence in a vertically anisotropic sand-and-gravel aquifer of medium size was nearly identical to that in an isotropic aquifer of the same size. Simulated area of influence increased as aquifer size increased and aquifer boundaries were farther away from the pit, and simulated drawdown was greater near the pit when aquifer boundaries were close to the pit. Pits simulated as lined with slurry walls caused mounding to occur upgradient from the pits and drawdown to occur downgradient from the pits. Pits simulated as refilled with water and undergoing evaporative losses had little hydro- logic effect on the aquifer. Numerical sensitivity analyses for simulations of pits in sand-and-gravel aquifers indicated that simulated head was most sensitive to horizontal hydraulic conductivity and the hydraulic conductance of general-head boundaries in the models. Simulated head was less sensitive to riverbed conductance and recharge and relatively insensitive to vertical hydraulic conductivity. Numerical simulations of quarries in fractured crystalline-rock aquifers indicated that the area of influence in a horizontally anisotropic aquifer was elongated in the direction of higher horizontal hydraulic conductivity and shortened in the direction of lower horizontal hydraulic conductivity compared to area of influence in a homogeneous, isotropic aquifer. Area of influence was larger in an aquifer with ground-water flow in deep, low-permeability fractures than in a homogeneous, isotropic aquifer. Area of influence was larger for a quarry intersected by a hydraulically conductive fault zone and smaller for a quarry intersected by a low-conductivity fault zone. Numerical sensitivity analyses for simulations of quarries in fractured crystalline-rock aquifers indicated simulated head was most sensitive to variations in recharge and horizontal hydraulic conductivity, had little sensitivity to vertical hydraulic conductivity and drain cells used to simulate valleys, and was relatively insensitive to drain cells used to simulate the quarry.

Water-Resources Investigations Report↗

USGS2018

Evaluation of the Source and Transport of High Nitrate Concentrations in Ground Water, Warren Subbasin, California

Ground water historically has been the sole source of water supply for the Town of Yucca Valley in the Warren subbasin of the Morongo ground-water basin, California. An imbalance between ground-water recharge and pumpage caused ground-water levels in the subbasin to decline by as much as 300 feet from the late 1940s through 1994. In response, the local water district, Hi-Desert Water District, instituted an artificial recharge program in February 1995 using imported surface water to replenish the ground water. The artificial recharge program resulted in water-level recoveries of as much as 250 feet in the vicinity of the recharge ponds between February 1995 and December 2001; however, nitrate concentrations in some wells also increased from a background concentration of 10 milligrams per liter to more than the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 44 milligrams per liter (10 milligrams per liter as nitrogen). The objectives of this study were to: (1) evaluate the sources of the high-nitrate concentrations that occurred after the start of the artificial-recharge program, (2) develop a ground-water flow and solute-transport model to better understand the source and transport of nitrates in the aquifer system, and (3) utilize the calibrated models to evaluate the possible effect of a proposed conjunctive-use project. These objectives were accomplished by collecting water-level and water-quality data for the subbasin and assessing changes that have occurred since artificial recharge began. Collected data were used to calibrate the ground-water flow and solute-transport models. Data collected for this study indicate that the areal extent of the water-bearing deposits is much smaller (about 5.5 square miles versus 19 square miles) than that of the subbasin. These water-bearing deposits are referred to in this report as the Warren ground-water basin. Faults separate the ground-water basin into five hydrogeologic units: the west, the midwest, the mideast, the east and the northeast hydrogeologic units. Water-quality analyses indicate that septage from septic tanks is the primary source of the high-nitrate concentrations measured in the Warren ground-water basin. Water-quality and stable-isotope data, collected after the start of the artificial recharge program, indicate that mixing occurs between imported water and native ground water, with the highest recorded nitrate concentrations in the midwest and the mideast hydrogeologic units. In general, the timing of the increase in measured nitrate concentrations in the midwest hydrogeologic unit is directly related to the distance of the monitoring well from a recharge site, indicating that the increase in nitrate concentrations is related to the artificial recharge program. Nitrate-to-chloride and nitrogen-isotope data indicate that septage is the source of the measured increase in nitrate concentrations in the midwest and the mideast hydrogeologic units. Samples from four wells in the Warren ground-water basin were analyzed for caffeine and selected human pharmaceutical products; these analyses suggest that septage is reaching the water table. There are two possible conceptual models that explain how high-nitrate septage reaches the water table: (1) the continued downward migration of septage through the unsaturated zone to the water table and (2) rising water levels, a result of the artificial recharge program, entraining septage in the unsaturated zone. The observations that nitrate concentrations increase in ground-water samples from wells soon after the start of the artificial recharge program in 1995 and that the largest increase in nitrate concentrations occur in the midwest and mideast hydrogeologic units where the largest increase in water levels occur indicate the validity of the second conceptual model (rising water levels). The potential nitrate concentration resulting from a water-level rise in the midwest and mideast hydrogeologic units was estimated using a simple mixing-cell model. The estimated value is within the range of concentrations measured in samples from wells, further indicating the validity of the second conceptual model. A ground-water flow model and a solute-transport model were developed for the Warren ground-water basin for the period 1956-2001. MODFLOW-96 was used for the ground-water flow model and MOC3D was used for the solute-transport model. The model cell size is about 500 feet by 500 feet and the models were discretized vertically into three layers. The models were calibrated using a trial-and-error approach using water-level and nitrate-concentration data collected between 1956-2001. In order to better match the measured data, low fault hydraulic characteristic values were required, thereby compartmentalizing the ground-water basin. In addition, it was necessary to parameterize the specific yield distribution for the top model layer where unconfined ground-water conditions occur into three homogeneous zones. Separate sets of specific- yield values were needed to simulate the drawdown and subsequent water-level recovery. In addition, the calibrated natural recharge was about 83 acre-feet per year. The entrainment of unsaturated-zone septage was simulated as recharge having an associated nitrate concentration. The volume of recharge was a function of the measured water-level rise between 1994-98 and the moisture content of the unsaturated zone. The nitrate concentration of the recharge water was a weighted function of the assumed nitrate concentration in the infiltrating water associated with the overlying land use. The simulated hydraulic head and nitrate concentration results were in good agreement with the measured data indicating that the mechanism for the increase in nitrate concentrations was rising water levels entraining high-nitrate septage in the unsaturated-zone. The calibrated models were used to simulate the possible effects of a planned conjunctive-use project in the western part of the ground-water basin. The simulated project included the addition of a new recharge pond and a new extraction well. In addition, recharge at two existing recharge ponds was increased and three existing production wells were pumped, treated in a nitrate-removal facility, and used for water supply. The simulated hydraulic heads increased in the west, the mideast, and parts of the east hydrogeologic units; however, the simulated hydraulic heads decreased in the midwest and northeast hydrogeologic units. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter (10 milligrams per liter as nitrogen) in parts of the west as a result of the increase in simulated hydraulic head. The simulated nitrate concentrations decreased in part of the midwest hydrogeologic unit as a result of the artificial recharge and pumping from the nitrate-removal wells. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter in part of the mideast and parts of the east hydrogeologic units beneath commercial land-use areas.

California↗

USGS2003

Recalibration of a ground-water flow model of the Mississippi River Valley alluvial aquifer of northeastern Arkansas, 1918-1998, with simulations of water levels caused by projected ground-water withdrawals through 2049

A digital model of the Mississippi River Valley alluvial aquifer in eastern Arkansas was used to simulate ground-water flow for the period from 1918 to 2049. The model results were used to evaluate effects on water levels caused by demand for ground water from the alluvial aquifer, which has increased steadily for the last 40 years. The model results showed that water currently (1998) is being withdrawn from the aquifer at rates greater than what can be sustained for the long term. The saturated thickness of the alluvial aquifer has been reduced in some areas resulting in dry wells, degraded water quality, decreased water availability, increased pumping costs, and lower well yields. The model simulated the aquifer from a line just north of the Arkansas-Missouri border to south of the Arkansas River and on the east from the Mississippi River westward to the less permeable geologic units of Paleozoic age. The model consists of 2 layers, a grid of 184 rows by 156 columns, and comprises 14,118 active cells each measuring 1 mile on a side. It simulates time periods from 1918 to 1998 along with further time periods to 2049 testing different pumping scenarios. Model flux boundary conditions were specified for rivers, general head boundaries along parts of the western side of the model and parts of Crowleys Ridge, and a specified head boundary across the aquifer further north in Missouri. Model calibration was conducted for observed water levels for the years 1972, 1982, 1992, and 1998. The average absolute residual was 4.69 feet and the root-mean square error was 6.04 feet for the hydraulic head observations for 1998. Hydraulic-conductivity values obtained during the calibration process were 230 feet per day for the upper layer and ranged from 230 to 730 feet per day for the lower layer with the maximum mean for the combined aquifer of 480 feet per day. Specific yield values were 0.30 throughout the model and specific storage values were 0.000001 inverse-feet throughout the model. Areally specified recharge rates ranged from 0 to about 30 inches and total recharge increased from 1972 to 1998 by a factor of about four. Water levels caused by projected ground-water withdrawals were simulated using the calibrated model. Simulations represented a period of 50 years into the future in three scenarios with either unchanged pumpage, pumpage increased by historic trends, or pumpage increased by historic trends except in two areas of the Grand Prairie. If pumping remains at 1997 rates, this produces extreme water-level declines (areas where model cells have gone dry or where the water level in the aquifer is equal to or less than the original saturated thickness, assuming confined conditions in the aquifer everywhere in the formation in predevelopment times) in the aquifer in two areas of the aquifer (one in the Grand Prairie area between the Arkansas and White Rivers and the other west of Crowleys Ridge along the Cache River) with about 400 square miles going dry. Increasing the pumping rates to that which would be projected using historic data led to increased extreme water-level declines in both areas with about 1,300 square miles going dry. Declines in both scenarios generally occurred most rapidly between 2009 and 2019. Reducing the pumping rates to 90 percent of that used for projected historic rates in areas between the Arkansas and White Rivers relating to two diversion projects of the U.S. Army Corps of Engineers and other agencies did little to decrease the extreme water-level declines. However, these pumpage reductions are small (amounting to about 16 percent of the reductions that could result from implementation of these diversion projects).

Water-Resources Investigations Report↗

USGS2003

Biogeochemical processes that produce dissolved organic matter from wheat straw

The chemical reactions that lead to the formation of dissolved organic matter (DOM) in natural waters are poorly understood. Studies on the formation of DOM generally are complicated because almost all DOM isolates have been derived from mixtures of plant species composed of a wide variety of different types of precursor compounds for DOM formation. This report describes a study of DOM derived mainly from bales of wheat straw that had been left in a field for several years. During this period of time, black water from the decomposing wheat straw accumulated in pools in the field. The nuclear magnetic resonance and infrared spectra of the black water DOM indicate that it is composed almost entirely of lignin and carbohydrate polymeric units. Analysis by high-performance size-exclusion chromatography with multi-angle laser-light scattering detection indicates that the number average molecular weight of the DOM is 124,000 daltons. The results presented in this report indicate that the black water DOM is composed of hemicellulose chains cross-linked to lignin oligomers. These types of structures have been shown to exist in the hemicellulose matrix of plant cell walls. The cross-linked lignin-hemicellulose complexes apparently were released from partially degraded wheat-straw cell walls with little alteration. In solution in the black water, these lignin-hemicellulose polymers fold into compact globular particles in which the nonpolar parts of the polymer form the interiors of the particles and the polar groups are on the exterior surfaces of the particles. The tightly folded, compact conformation of these particles probably renders them relatively resistant to microbial degradation. This should be especially the case for the aromatic lignin structures that will be buried in the interiors of the particles.

Water-Resources Investigations Report↗