USGS Science⌕ Search

SEARCH · USGS Science

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 487 records · Page 27Linked to original sources
USGS2012

Geologic framework for the national assessment of carbon dioxide storage resources: Bighorn Basin, Wyoming and Montana: Chapter A in Geologic framework for the national assessment of carbon dioxide storage resources

The 2007 Energy Independence and Security Act ( Public Law 110–140 ) directs the U.S. Geological Survey (USGS) to conduct a national assessment of potential geologic storage resources for carbon dioxide (CO 2 ). The methodology used for the national CO 2 assessment follows that of previous USGS work. The methodology is non-economic and intended to be used at regional to subbasinal scales. This report identifies and contains geologic descriptions of twelve storage assessment units (SAUs) in six separate packages of sedimentary rocks within the Bighorn Basin of Wyoming and Montana and focuses on the particular characteristics, specified in the methodology, that influence the potential CO 2 storage resource in those SAUs. Specific descriptions of the SAU boundaries as well as their sealing and reservoir units are included. Properties for each SAU such as depth to top, gross thickness, net porous thickness, porosity, permeability, groundwater quality, and structural reservoir traps are provided to illustrate geologic factors critical to the assessment. Although assessment results are not contained in this report, the geologic information included here will be employed, as specified in the methodology of earlier work, to calculate a statistical Monte Carlo-based distribution of potential storage space in the various SAUs. Figures in this report show SAU boundaries and cell maps of well penetrations through the sealing unit into the top of the storage formation. Wells sharing the same well borehole are treated as a single penetration. Cell maps show the number of penetrating wells within one square mile and are derived from interpretations of incompletely attributed well data, a digital compilation that is known not to include all drilling. The USGS does not expect to know the location of all wells and cannot guarantee the amount of drilling through specific formations in any given cell shown on cell maps.

Wyoming, Montana↗
USGS2012

Geologic framework for the national assessment of carbon dioxide storage resources: Hanna, Laramie, and Shirley Basins, Wyoming: Chapter C in Geologic framework for the national assessment of carbon dioxide storage resources

The 2007 Energy Independence and Security Act (Public Law 110-140) directs the U.S. Geological Survey (USGS) to conduct a national assessment of potential geologic storage resources for carbon dioxide (CO 2 ). The methodology used for the national CO 2 assessment is non-economic and intended to be used at regional to subbasinal scales. This report identifies and contains geologic descriptions of twelve storage assessment units (SAUs) in six separate packages of sedimentary rock within the Hanna, Laramie, and Shirley Basins of Wyoming. It focuses on the particular characteristics, specified in the methodology, that influence the potential CO 2 storage resource in those SAUs. Specific descriptions of SAU boundaries as well as their sealing and reservoir units are included. Properties for each SAU, such as depth to top, gross thickness, net porous thickness, porosity, permeability, groundwater quality, and structural reservoir traps are provided to illustrate geologic factors critical to the assessment. Although assessment results are not contained in this report, the geologic information included herein will be employed, as specified in the methodology, to calculate a statistical Monte Carlo-based distribution of potential storage space in the various SAUs. Figures in this report show SAU boundaries and cell maps of well penetrations through the sealing unit into the top of the storage formation. Cell maps show the number of penetrating wells within one square mile and are derived from interpretations of incompletely attributed well data in a digital compilation that is known not to include all drilling. The USGS does not expect to know the location of all wells and cannot guarantee the amount of drilling through specific formations in any given cell shown on cell maps.

Wyoming↗
USGS2013

Geologic framework for the national assessment of carbon dioxide storage resources: Columbia Basin of Oregon, Washington, and Idaho, and the Western Oregon-Washington basins

The 2007 Energy Independence and Security Act (Public Law 110–140) directs the U.S. Geological Survey (USGS) to conduct a national assessment of potential geologic storage resources for carbon dioxide (CO 2 ). The methodology used by the USGS for the national CO 2 assessment follows that of previous USGS work. The methodology is non-economic and intended to be used at regional to subbasinal scales. This report identifies and contains geologic descriptions of three storage assessment units (SAUs) in Eocene and Oligocene sedimentary rocks within the Columbia, Puget, Willapa, Astoria, Nehalem, and Willamette Basins of Oregon, Washington, and Idaho, and focuses on the characteristics, specified in the methodology, that influence the potential CO 2 storage resource in those SAUs. Specific descriptions of the SAU boundaries as well as their sealing and reservoir units are included. Properties for each SAU, such as depth to top, gross thickness, porosity, permeability, groundwater quality, and structural reservoir traps, are provided to illustrate geologic factors critical to the assessment. The designated sealing unit in the Columbia Basin is tentatively chosen to be the ubiquitous and thick Miocene Columbia River Basalt Group. As a result of uncertainties regarding the seal integrity of the Columbia River Basalt Group, the SAUs were not quantitatively assessed. Figures in this report show SAU boundaries and cell maps of well penetrations through sealing units into the top of the storage formations. The cell maps show the number of penetrating wells within one square mile and are derived from interpretations of incompletely attributed well data, a digital compilation that is known not to include all drilling. The USGS does not expect to know the location of all wells and cannot guarantee the amount of drilling through specific formations in any given cell shown on the cell maps.

Idaho, Oregon, Washington↗
USGS2013

Geologic framework for the national assessment of carbon dioxide storage resources: Arkoma Basin, Kansas Basins, and Midcontinent Rift Basin study areas

2007 Energy Independence and Security Act (Public Law 110–140) directs the U.S. Geological Survey (USGS) to conduct a national assessment of potential geologic storage resources for carbon dioxide (CO 2 ). The methodology used by the USGS for the national CO 2 assessment follows that of previous USGS work. This methodology is non-economic and intended to be used at regional to subbasinal scales. This report identifies and contains geologic descriptions of three storage assessment units (SAUs) in Upper Cambrian to Mississippian sedimentary rocks within the Arkoma Basin study area, and two SAUs in Upper Cambrian to Mississippian sedimentary rocks within the Kansas Basins study area. The Arkoma Basin and Kansas Basins are adjacent with very similar geologic units; although the Kansas Basins area is larger, the Arkoma Basin is more structurally complex. The report focuses on the characteristics, specified in the methodology, that influence the potential CO 2 storage resource in the SAUs. Specific descriptions of the SAU boundaries as well as their sealing and reservoir units are included. Properties for each SAU, such as depth to top, gross thickness, porosity, permeability, groundwater quality, and structural reservoir traps, are usually provided to illustrate geologic factors critical to the assessment. Although assessment results are not contained in this report, the geologic information herein was employed, as specified in the USGS methodology, to calculate a probabilistic distribution of potential storage resources in each SAU. The Midcontinent Rift Basin study area was not assessed, because no suitable storage formations meeting our size, depth, reservoir quality, and regional seal guidelines were found. Figures in this report show study area boundaries along with the SAU boundaries and cell maps of well penetrations through sealing units into the top of the storage formations. The cell maps show the number of penetrating wells within one-square mile and are derived from interpretations of incompletely attributed well data and from a digital compilation that is known not to include all drilling. The USGS does not expect to know the location of all wells and cannot guarantee the amount of drilling through specific formations in any given cell shown on the cell maps.

Arkansas;Louisiana;Oklahoma;Texas↗
USGS2014

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↗
USGS2025

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↗
USGS1988

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↗
USGS1989

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↗
USGS1992

A generalized finite-difference formulation for the U.S. Geological Survey modular three-dimensional finite-difference ground-water flow model

The U.S. Geological Survey's Modular Ground-Water Flow Model assumes that model nodes are in the center of cells and that transmissivity is constant within a cell. Based on these assumptions, the model calculates coefficients, called conductance, that are multiplied by head difference to determine flow between cells. Although these are common assumptions in finite-difference models, other assumptions are possible. A new option to the model program reads conductance as input data rather than calculating it. This optional lows the user to calculate conductance outside of the model. The user thus has the flexibility to define conductance using any desired assumptions. For a water-table condition, horizontal conductance must change as water level varies. To handle this situation, the new option reads conductance divided by thickness (CDT) as input data. The model calculates saturated thickness and multiplies it by CDT to obtain conductance. Thus, the user is still free from the assumptions of centered nodes and constant transmissivity in cells. The model option is written in FORTRAN77 and is fully compatible with the existing model. This report documents the new model option; it includes a description of the concepts, detailed input instructions, and a listing of the code.

Open-File Report↗
USGS1996

Conversion of the Twin Cities metropolitan area numerical ground-water-flow model from the Trescott-Larson computer code to the McDonald-Harbaugh computer code

A numerical model of ground-water flow in the Twin Cities metropolitan area was converted from the TrescottLarson computer code to the McDonald-Harbaugh computer code to facilitate current and future use of the model using up-to-date computer software and hardware. Differences exist between the two computer codes in how headdependent source-sink functions, including leaky rivers, springs and seepage faces (drains), and head-dependent flux boundaries, are formulated. Because of differences in the formulation and calculation of the conductance terms, the conductance values for the head-dependent source-sink functions from the Trescott-Larson Twin Cities model were multiplied by the area of the appropriate model cell for conversion to the McDonald-Harbaugh Twin Cities model. Leaky rivers were simulated using the river package in the McDonald-Harbaugh Twin Cities model, springs and seepage faces were simulated using the drain package, and head-dependent flux boundaries were simulated using the general-head boundary package. Hydraulic heads and flows in the aquifer system calculated by the McDonald-Harbaugh Twin Cities model were compared to those calculated by the Trescott-Larson Twin Cities model to verify that the results calculated by the two models are similar. Mean differences in calculated hydraulic heads for the drift, drift and St. Peter aquifer, and drift and Prairie du Chien-Jordan aquifer model layers ranged from 1.1 to 1.6 feet. Mean differences in calculated hydraulic heads for the drift and Ironton-Galesville aquifer, drift and Eau Claire confining unit, and drift and Mount Simon-Hinckley aquifer model layers were 0.4 feet or less. Large differences in calculated hydraulic heads at a few cells probably were due to differences in the solution methods used in the McDonald-Harbaugh and Trescott-Larson computer codes to calculate hydraulic heads at and near model layer boundaries. Differences in calculated flow rates for the McDonald-Harbaugh and Trescott-Larson Twin Cities models were less than 0.3 cubic feet per second for recharge and the head-dependent source-sink functions. Differences in calculated flow in and flow out of constant-head cells for each model layer for the two models ranged from about 2 cubic feet per second for the Mount Simon-Hinckley aquifer model layer to about 45 cubic feet per second for the Prairie du ChienJordan aquifer model layer. The differences between the net flow rates at constant-head cells calculated by the two models for each model layer were much smaller, ranging from 0.01 cubic feet per second for the Ironton-Galesville and Mount Simon-Hinckley aquifer model layers to 8.39 cubic feet per second for the St. Peter aquifer model layer. Differences in the calculated ground-water budgets for the two models were 5.3 percent for the total sources and 4.4 percent for the total sinks.

Minnesota↗
USGS2005

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↗
USGS2008

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↗
USGS2009

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↗
USGS2009

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↗
USGS2011

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↗
USGS2024

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↗
USGS1958

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↗
USGS2011

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↗