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

Assigning boundary conditions to the Southern Inland and Coastal Systems (SICS) model using results from the South Florida Water Management Model (SFWMM)

The Comprehensive Everglades Restoration Plan (CERP) requires the testing and evaluation of different water-management scenarios for southern Florida. As part of CERP, the South Florida Water Management District is using its regional hydrologic model, the South Florida Water Management Model (SFWMM), to evaluate different hydrologic scenarios. The SFWMM was designed specifically for the inland freshwater areas in southern Florida, and extends only slightly into Florida Bay. Thus, the U.S. Geological Survey developed the Southern Inland and Coastal Systems (SICS) model, which is an integrated surface-water and ground-water model designed to simulate flows, stages, and salinities in the southern Everglades and Florida Bay. Modifications to the SICS boundary conditions allow the local-scale SICS model to be linked to the regional-scale SFWMM. The linked model will be used to quantify the effects of restoration alternatives on flows, stages, and salinities in the SICS area. This report describes the procedure for linking the SICS model with the SFWMM. The linkage is shown to work by comparing the results of a linked 5-year simulation with the results from a simulation in which the model boundaries are assigned using field data. The surface-water module of the SICS model is driven by areal influences and lateral boundaries. The areal influences (wind, rainfall, and evapotranspiration) remain the same when the SICS model is modified to link to the SFWMM. Four types of lateral boundaries (discharge, water level, no flow, and salinity) are used in the SICS model. Two of three discharge boundaries (at Taylor Slough Bridge and C-111 Canal) in the current SICS model domain are converted to water-level boundaries to increase accuracy. The only change to the third discharge boundary (at Levee 31W) is that the flow data are derived from SFWMM model output instead of using measured field data flows. Three water-level boundaries are modified only by receiving their data from SFWMM model output data. Additionally, two marine water-level boundaries remain the same because the SFWMM does not include Florida Bay and, therefore, this model cannot provide input data for these boundaries. The SICS no-flow boundaries remain intact because no additional data, provided by the SFWMM, suggest that any significant flow occurs along these boundaries. The Florida Bay salinity boundary is not modified because the SFWMM does not contain any salinity data that can be used to modify the model. The ground-water module of the SICS model contains a general-head boundary and a no-flow boundary. The general-head boundary, which extends along the edges of the wetland part of the SICS model domain, is modified by acquiring stage values from SFWMM cells that correspond in location to the SICS model cells. Values from the SFWMM cells are bilinearly interpolated and assigned to the appropriate SICS general-head boundary cells in all layers of the ground-water model. The ground-water no-flow boundary in Florida Bay is unaltered because the SFWMM does not include this area. A 5-year simulation was developed to test the linkage of the SICS model with the SFWMM. Results from the linked model are similar to those obtained from the original SICS model in which boundaries are assigned using field data. The simulated discharges at the coastal creeks along Florida Bay are about 5 percent lower than the field data simulation; water levels in the wetlands are about 4 percent lower, and salinities at the various coastal creeks are slightly higher.

Florida↗
USGS2007

Conceptual design of the Everglades Depth Estimation Network (EDEN) grid

The Everglades Depth Estimation Network (EDEN) offers a consistent and documented dataset that can be used to guide large-scale field operations, to integrate hydrologic and ecological responses, and to support biological and ecological assessments that measure ecosystem responses to the Comprehensive Everglades Restoration Plan (Telis, 2006). Ground elevation data for the greater Everglades and the digital ground elevation models derived from them form the foundation for all EDEN water depth and associated ecologic/hydrologic modeling (Jones, 2004, Jones and Price, 2007). To use EDEN water depth and duration information most effectively, it is important to be able to view and manipulate information on elevation data quality and other land cover and habitat characteristics across the Everglades region. These requirements led to the development of the geographic data layer described in this techniques and methods report. Relying on extensive experience in GIS data development, distribution, and analysis, a great deal of forethought went into the design of the geographic data layer used to index elevation and other surface characteristics for the Greater Everglades region. To allow for simplicity of design and use, the EDEN area was broken into a large number of equal-sized rectangles ('Cells') that in total are referred to here as the 'grid'. Some characteristics of this grid, such as the size of its cells, its origin, the area of Florida it is designed to represent, and individual grid cell identifiers, could not be changed once the grid database was developed. Therefore, these characteristics were selected to design as robust a grid as possible and to ensure the grid's long-term utility. It is desirable to include all pertinent information known about elevation and elevation data collection as grid attributes. Also, it is very important to allow for efficient grid post-processing, sub-setting, analysis, and distribution. This document details the conceptual design of the EDEN grid spatial parameters and cell attribute-table content.

Florida↗
USGS2008

TRIGRS - A Fortran Program for Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Analysis, Version 2.0

The Transient Rainfall Infiltration and Grid-Based Regional Slope-Stability Model (TRIGRS) is a Fortran program designed for modeling the timing and distribution of shallow, rainfall-induced landslides. The program computes transient pore-pressure changes, and attendant changes in the factor of safety, due to rainfall infiltration. The program models rainfall infiltration, resulting from storms that have durations ranging from hours to a few days, using analytical solutions for partial differential equations that represent one-dimensional, vertical flow in isotropic, homogeneous materials for either saturated or unsaturated conditions. Use of step-function series allows the program to represent variable rainfall input, and a simple runoff routing model allows the user to divert excess water from impervious areas onto more permeable downslope areas. The TRIGRS program uses a simple infinite-slope model to compute factor of safety on a cell-by-cell basis. An approximate formula for effective stress in unsaturated materials aids computation of the factor of safety in unsaturated soils. Horizontal heterogeneity is accounted for by allowing material properties, rainfall, and other input values to vary from cell to cell. This command-line program is used in conjunction with geographic information system (GIS) software to prepare input grids and visualize model results.

Open-File Report↗
USGS2010

Microbial and geochemical investigations of dissolved organic carbon and microbial ecology of native waters from the Biscayne and Upper Floridan Aquifers

Groundwater resources in the United States are under ever-increasing demands for potable, irrigation, and recreational uses. Additionally, aquifer systems are being used or targeted for use as storage areas for treated surface waters and (or) groundwaters via injection (for example, aquifer storage and recovery). To date, the influence that the nutrients, including carbon, in the injected water have on native microbial communities and the biogeochemistry in the subsurface zones used for storage of the injectate has not been determined. In this report, we describe a series of experiments that establishes a baseline dataset for the quantity and quality of organic and inorganic carbon and nutrients in the Biscayne Aquifer (BA) and Upper Floridan Aquifer (UFA) in south Florida. The most significant differences between the BA (26 meters below surface) and UFA (366 meters below surface) are the average specific conductance (0.552 and 6.12 microsiemens per centimeter, respectively), dissolved oxygen (1.6 and 0 milligrams per liter, respectively), and oxidation-reduction potential (40.3 and -358 millivolts, respectively). The dissolved organic carbon from the BA is characterized by carbon originating from terrestrial sources and microbial activities, while the UFA has a distinctive microbial signature. Acetate and lactate are the dominant carbon constituents in both aquifers. Additionally, components of the dissolved organic carbon from the UFA have a total trihalomethane-formation potential that is approximately threefold greater than the maximum contaminat level of 80 micrograms per liter established by the U.S. Environmental Protection Agency. The average native bacterial abundances in the aquifers are similar with 4.69x10^4 cells per milliliter in the BA and 1.33x10^4 cells per milliliter in the UFA. The average bacteriophage abundances are also similar with 1.15x10^5 virus-like particles in the BA and 1.92x10^5 virus-like particles in the UFA. Interestingly, ciliated protozoa are present in both aquifers. The average abundance of ciliates in the BA (2.97x10^3 ciliates per milliliter) is approximately twentyfold greater than abundances in the UFA (1.39x10^2 ciliates per milliliter). Collectively, these data indicate that microbial processes are the dominant contributor to the cycling of carbon and inorganic carbon in the BA and may be the only carbon cycling process in the UFA, as this aquifer has not had a terrestrial influx of carbon for more than 15,000 years. The rates of carbon, in the form of acetate, utilization by the native microbial communities are significantly different between the two aquifers. Based on data from 14C-acetate-utilization experiments, the microbial communities in the BA turn over the native acetate in 2.5 years, whereas communities in the UFA turn over native acetate in 6.8 years. These data support the hypothesis derived from the microbial-abundance data, in that the carbon for bacterial maintainence and growth is recycled from bacterial biomass released during cell lysis, especially in the UFA. An in situ diffusion chamber was designed to retain bacterial cells within the chamber while allowing native water constituents to move through the chamber. A series of 1-week deployments of chambers filled with fluorescent beads, inactivated native bacteria and laboratory grown and viable bacteria into the UFA, permitted by the State of Florida Environmental Protection Agency, was successfully completed. This was the first time this type of deployment into an aquifer system that is used for potable water supply has been permitted within the United States. This technology will allow, for the first time, in situ studies on the survival of microbial indicators of fecal pollution and true pathogens in groundwater systems.

Florida↗
USGS2012

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 ) and to consult with other Federal and State agencies to locate the pertinent geological data needed for the assessment. The geologic sequestration of CO 2 is one possible way to mitigate its effects on climate change. The methodology used for the national CO 2 assessment (Open-File Report 2010-1127; http://pubs.usgs.gov/of/2010/1127/) is based on previous USGS probabilistic oil and gas assessment methodologies. The methodology is non-economic and intended to be used at regional to subbasinal scales. The operational unit of the assessment is a storage assessment unit (SAU), composed of a porous storage formation with fluid flow and an overlying sealing unit with low permeability. Assessments are conducted at the SAU level and are aggregated to basinal and regional results. This report identifies and contains geologic descriptions of SAUs in separate packages of sedimentary rocks within the assessed basin 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, 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.

Open-File Report↗
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↗
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↗