Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Computers and Geology”

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 1,045 records · Page 58Linked to original sources
USGS2008

Documentation of Computer Program INFIL3.0 - A Distributed-Parameter Watershed Model to Estimate Net Infiltration Below the Root Zone

This report documents the computer program INFIL3.0, which is a grid-based, distributed-parameter, deterministic water-balance watershed model that calculates the temporal and spatial distribution of daily net infiltration of water across the lower boundary of the root zone. The bottom of the root zone is the estimated maximum depth below ground surface affected by evapotranspiration. In many field applications, net infiltration below the bottom of the root zone can be assumed to equal net recharge to an underlying water-table aquifer. The daily water balance simulated by INFIL3.0 includes precipitation as either rain or snow; snowfall accumulation, sublimation, and snowmelt; infiltration into the root zone; evapotranspiration from the root zone; drainage and water-content redistribution within the root-zone profile; surface-water runoff from, and run-on to, adjacent grid cells; and net infiltration across the bottom of the root zone. The water-balance model uses daily climate records of precipitation and air temperature and a spatially distributed representation of drainage-basin characteristics defined by topography, geology, soils, and vegetation to simulate daily net infiltration at all locations, including stream channels with intermittent streamflow in response to runoff from rain and snowmelt. The model does not simulate streamflow originating as ground-water discharge. Drainage-basin characteristics are represented in the model by a set of spatially distributed input variables uniquely assigned to each grid cell of a model grid. The report provides a description of the conceptual model of net infiltration on which the INFIL3.0 computer code is based and a detailed discussion of the methods by which INFIL3.0 simulates the net-infiltration process. The report also includes instructions for preparing input files necessary for an INFIL3.0 simulation, a description of the output files that are created as part of an INFIL3.0 simulation, and a sample problem that illustrates application of the code to a field setting. Brief descriptions of the main program routine and of each of the modules and subroutines of the INFIL3.0 code, as well as definitions of the variables used in each subroutine, are provided in an appendix.

Scientific Investigations Report↗
USGS1982

Obtaining maps and data from the U.S. Geological Survey

The U.S. Geological Survey produces a variety of resource information for the United States. This includes many data bases of particular interest to planners such as land use and terrain information prepared by the National Mapping Division, water quantity and quality data collected by Water Resources Division, and coal resource information gathered by the Geologic Division. These data are stored in various forms, and information on their availability can be obtained from appropriate offices in the U.S. Geological Survey as well as from USGS Circular 777. These data have been used for the management, development, and monitoring of our Nation's resources by Federal, State, and local agencies.

Computers, Environment and Urban Systems↗
USGS2016

Uncertainty in V s30 -based site response

Methods that account for site response range in complexity from simple linear categorical adjustment factors to sophisticated nonlinear constitutive models. Seismic‐hazard analysis usually relies on ground‐motion prediction equations (GMPEs); within this framework site response is modeled statistically with simplified site parameters that include the time‐averaged shear‐wave velocity to 30 m ( V S 30 ) and basin depth parameters. Because V S 30 is not known in most locations, it must be interpolated or inferred through secondary information such as geology or topography. In this article, we analyze a subset of stations for which V S 30 has been measured to address effects of V S 30 proxies on the uncertainty in the ground motions as modeled by GMPEs. The stations we analyze also include multiple recordings, which allow us to compute the repeatable site effects (or empirical amplification factors [EAFs]) from the ground motions. Although all methods exhibit similar bias, the proxy methods only reduce the ground‐motion standard deviations at long periods when compared to GMPEs without a site term, whereas measured V S 30 values reduce the standard deviations at all periods. The standard deviation of the ground motions are much lower when the EAFs are used, indicating that future refinements of the site term in GMPEs have the potential to substantially reduce the overall uncertainty in the prediction of ground motions by GMPEs.

Bulletin of the Seismological Society of America↗
USGS2012

Real-time seismic monitoring of instrumented hospital buildings

In collaboration with the Department of Veterans Affairs (VA), the U.S. Geological Survey's National Strong Motion Project has recently installed sophisticated seismic monitoring systems to monitor the structural health of two hospital buildings at the Memphis VA Medical Center in Tennessee. The monitoring systems in the Bed Tower and Spinal Cord Injury buildings combine sensing technologies with an on-site computer to capture and analyze seismic performance of buildings in near-real time.

Tennessee↗
USGS2012

Advanced earthquake monitoring system for U.S. Department of Veterans Affairs medical buildings–Instrumentation

In collaboration with the U.S. Department of Veterans Affairs (VA), the National Strong Motion Project (NSMP; http://nsmp.wr.usgs.gov/ ) of the U.S. Geological Survey has been installing sophisticated seismic systems that will monitor the structural integrity of 28 VA hospital buildings located in seismically active regions of the conterminous United States, Alaska, and Puerto Rico during earthquake shaking. These advanced monitoring systems, which combine the use of sensitive accelerometers and real-time computer calculations, are designed to determine the structural health of each hospital building rapidly after an event, helping the VA to ensure the safety of patients and staff. This report presents the instrumentation component of this project by providing details of each hospital building, including a summary of its structural, geotechnical, and seismic hazard information, as well as instrumentation objectives and design. The structural-health monitoring component of the project, including data retrieval and processing, damage detection and localization, automated alerting system, and finally data dissemination, will be presented in a separate report.

Open-File Report↗
USGS2017

Compilation of streamflow statistics calculated from daily mean streamflow data collected during water years 1901–2015 for selected U.S. Geological Survey streamgages

Streamflow statistics are needed by decision makers for many planning, management, and design activities. The U.S. Geological Survey (USGS) StreamStats Web application provides convenient access to streamflow statistics for many streamgages by accessing the underlying StreamStatsDB database. In 2016, non-interpretive streamflow statistics were compiled for streamgages located throughout the Nation and stored in StreamStatsDB for use with StreamStats and other applications. Two previously published USGS computer programs that were designed to help calculate streamflow statistics were updated to better support StreamStats as part of this effort. These programs are named “GNWISQ” (Get National Water Information System Streamflow (Q) files), updated to version 1.1.1, and “QSTATS” (Streamflow (Q) Statistics), updated to version 1.1.2. Statistics for 20,438 streamgages that had 1 or more complete years of record during water years 1901 through 2015 were calculated from daily mean streamflow data; 19,415 of these streamgages were within the conterminous United States. About 89 percent of the 20,438 streamgages had 3 or more years of record, and about 65 percent had 10 or more years of record. Drainage areas of the 20,438 streamgages ranged from 0.01 to 1,144,500 square miles. The magnitude of annual average streamflow yields (streamflow per square mile) for these streamgages varied by almost six orders of magnitude, from 0.000029 to 34 cubic feet per second per square mile. About 64 percent of these streamgages did not have any zero-flow days during their available period of record. The 18,122 streamgages with 3 or more years of record were included in the StreamStatsDB compilation so they would be available via the StreamStats interface for user-selected streamgages. All the statistics are available in a USGS ScienceBase data release.

Open-File Report↗
USGS1973

Atlas of Mars, Oxia Palus quadrangle, MC-11

Th is mosaic is made with computer enhanced Mariner 9 pictures, rectified and scaled to a Mercator projection. Computer processing was performed by the Image Processing Laboratory (IPL) of the Jet Propulsion Laboratory (JPL). Processing included contrast enhancement and high-pass filtering which accentuates topographic features while subduing albedo variations. Photographic reproductions of the images were generated from the IPL digital magnetic tapes and mosaicked by the U.S. Geological Survey. Relative horizontal placement of pictures is controlled by spacecraft tracking data, with adjustment to fit the Mariner 9 control net of August, 1972 (Davies, M.E., 1972, Mariner 9 control net of Mars: August, 1972: Rand Corp., R-1 122- JPL) where possible. Individual analytical adjustment of pictures was not performed. Triangles are shown at the correct latitude and longitude of the control points identified and numbered by Davies (ibid.). Broken crosses are used to identify the images of the control points. Standard deviation of image placement discrepancy is less than 5 km. The projection is based on a theoretical dynamical ellipticity of 1/192. The scale is based on an equatorial radius of 3393.4 km.

Open-File Report↗
USGS1973

Atlas of Mars, Amenthes quadrangle, MC-14

This mosaic is made with computer enhanced Mariner 9 pictures, rectified and scaled to a Mercator projection. Computer processing was performed by the Image Processing Laboratory (IPL) of the Jet Propulsion Laboratory (JPL). Processing included contrast enhancement and high-pass filtering which accentuates topographic features while subduing albedo variations. Photographic reproductions of the images were generated from the IPL digital magnetic tapes and mosaicked by the U.S. Geological Survey. Relative horizontal placement of pictures is controlled by spacecraft tracking data, with adjustment to fit the Mariner 9 control net of August, 1972 (Davies, M.E., 1972, Mariner 9 control net of Mars: August, 1972: Rand Corp., R-1122-JPL) where possible. Individual analytical adjustment of pictures was not performed. Triangles are shown at the correct latitude and longitude of the control points identified and numbered by Davies (ibid.). Broken crosses are used to identify the images of the control points. Standard deviation of image placement discrepancy is less than 5 km. The-projection is based on a theoretical dynamical ellipticity of 1/192. The scale is based on an equatorial radius of 3393.4 km.

Open-File Report↗
USGS1973

Atlas of Mars, Phoenicis Lacus quadrangle, MC-17

This mosaic is made with computer enhanced Mariner 9 pictures, rectified and scaled to a Mercator projection. Computer processing was performed by the Image Processing Laboratory (IPL) of the Jet Propulsion Laboratory (JPL). Processing included contrast enhancement and high-pass filtering which accentuates topographic features while subduing albedo variations. Photographic reproductions of the images were generated from the IPL digital magnetic tapes and mosaicked by the U.S. Geological Survey. Horizontal controls were derived from spacecraft tracking data, which are based on astronomical definitions of the spin axis and system of longitudes in use during the Mariner 9 mission. The projection is based on a theoretical dynamical ellipticity of 1/192. The scale is based on an equatorial radius of 3393.4 km. Standard deviation in image placement is less than 5 km at map scale.

Open-File Report↗
USGS1976

Field verification of method for distributing flow through multiple-bridge openings

Field data collected at three sites in Louisiana having multiple-bridge openings were used to check a method of distribution of flow through two or more bridges. Actual peak discharges in the various main and overflow bridges ranged from 878 to 11,055 cubic feet per second (24.9 to 313 m 2 /s). The error of the computed discharge to the actual discharge ranged between -20 and +6 percent, with the exception of one opening that had an error of +101 percent.

Louisiana↗
USGS2024

Challenges in developing vertical hazard for seismic analysis of concrete dams

The seismic response of concrete dams depends on both the horizontal (H) and vertical (V) components of ground motion (GM), and excluding the V components when conducting response history analyses (RHAs) may underestimate the seismic fragility. Although V components of GM time series could be selected to be consistent with the hazard, hazard curves for the V component (or for short, V hazard curves) are rarely available for evaluating such hazard consistency. One of the challenges in developing V hazard curves in the U.S. Geological Survey (USGS) National Seismic Hazard Model (NSHM) is the lack of ground motion models (GMMs) that directly estimate intensity measures (IMs) for the V component of GM. In this study, we compare alternative approaches of developing V hazard curves. Using a case study in California, we compare V hazard curves computed from GMMs for the V component of GM against those from applying GMMs for the V/H ratios and those from selecting multicomponent time series that are consistent with the H hazard. We discuss and invite feedback on issues such as availability of GMMs, applicability of GMMs, and choice of inputs for ground motion selection and modification (GMSM).

Conference Paper↗
USGS1996

Evaluation and modification of five techniques for estimating stormwater runoff for watersheds in west-central Florida

Several traditional techniques have been used for estimating stormwater runoff from ungaged watersheds. Applying these techniques to watersheds in west-central Florida requires that some of the empirical relations be extrapolated beyond tested ranges. As a result, there is uncertainty as to the accuracy of these estimates. Sixty-six storms occurring in 15 west-central Florida watersheds were initially modeled using the Rational Method, the U.S. Geological Survey Regional Regression Equations, the Natural Resources Conservation Service TR-20 model, the U.S. Army Corps of Engineers HEC-1 model, and the U.S. Environmental Protection Agency Storm Water Management Model. The techniques were applied according to the guidelines specified in the user manuals or standard engineering textbooks as though no field data were available and the selection of input parameters was not influenced by observed data. Computed estimates were compared with observed runoff to evaluate the accuracy of the techniques. One watershed was eliminated from further evaluation when it was determined that the area contributing runoff to the stream varies with the amount and intensity of rainfall. Therefore, further evaluation and modification of the input parameters were made for only 62 storms in 14 watersheds. Runoff ranged from 1.4 to 99.3 percent of rainfall. The average runoff for all watersheds included in this study was about 36 percent of rainfall. The average runoff for the urban, natural, and mixed land use watersheds was about 41, 27, and 29 percent of rainfall, respectively. Initial estimates of peak discharge using the Rational Method produced average watershed errors that ranged from an underestimation of 50.4 percent to an overestimation of 767 percent. The coefficient of runoff ranged from 0.20 to 0.60. Calibration of the technique produced average errors that ranged from an underestimation of 3.3 percent to an over estimation of 1.5 percent. The average calibrated coefficient of runoff for each watershed ranged from 0.02 to 0.72. The average values of the coefficient of runoff necessary to calibrate the urban, natural, and mixed land use watersheds were 0.39, 0.16, and 0.08, respectively. The U.S. Geological Survey regional regression equations for determining peak discharge produced errors that ranged from an underestimation of 87.3 percent to an overestimation of 1,140 percent. The regression equations for determining runoff volume produced errors that ranged from an underestimation of 95.6 percent to an overestimation of 324 percent. Regression equations developed from data used for this study produced errors that ranged between an underestimation of 82.8 percent and an overestimation of 328 percent for peak discharge and from an underestimation of 71.2 percent to an overestimation of 241 percent for runoff volume. Use of the equations developed for west-central Florida streams produced average errors for each type of watershed that were lower than errors associated with use of the U.S. Geological Survey regional equations. Initial estimates of peak discharges and runoff volumes using the Natural Resources Conservation Service TR-20 model, produced average errors of 44.6 and 42.7 percent, respectively, for all the watersheds. Curve numbers and times of concentration were adjusted to match estimated and observed peak discharges and runoff volumes. The average change in the curve number for all the watersheds was a decrease of 2.8 percent. The average change in the time of concentration was an increase of 59.2 percent. The shape of the input dimensionless unit hydrograph also had to be adjusted to match the shape and peak time of the estimated and observed flood hydrographs. Peak rate factors for the modified input dimensionless unit hydrographs ranged from 162 to 454. The mean errors for peak discharges and runoff volumes were reduced to 18.9 and 19.5 percent, respectively, using the average calibrated input parameters for each watershed. Initial estimates of peak discharges and runoff volumes using the U.S. Army Corp of Engineers Hydrologic Engineering Center-1 model, produced average errors of 105 and 26.8 percent respectively, for all the watersheds. Curve numbers and lag times were adjusted to match estimated and observed peak discharges and runoff volumes. The average change in the curve number for all the watersheds was a decrease of 2.5 percent. The average change in the lag time was an increase of 169 percent. The mean errors for peak discharges and runoff volumes were reduced to 5.8 and 1.4 percent, respectively, using the average calibrated input parameters for each watershed. The observed and estimated peak discharges and runoff volumes could be matched by adjusting curve numbers and lag time using the U.S. Army Corp of Engineers Hydrologic Engineering Center-1 model; however, the shape of the estimated flood hydrograph and timing of the peak could not be matched. The input dimensionless unit hydrograph must also be changed to increase the accuracy of the Hydrologic Engineering Center-1 model for watersheds in west-central Florida. The source code has to be modified and recompiled to enter different dimensionless unit hydrographs into the HEC-1 program. During application of the U.S. Environmental Protection Agency Storm Water Management Model, two separate infiltration methods were evaluated. Initial estimates of peak discharges and runoff volumes produced mean errors of 46.5 and 6.8 percent, respectively, for all watersheds using the Green-Ampt infiltration method, and 48.8 and 9.5 percent, respectively, using the Horton infiltration method. The mean errors were reduced to 18 and 0.3 percent for the Green-Ampt method and 20.9 and 7.2 percent for the Horton method using the average calibrated input parameters for each watershed. Estimates of peak discharges and runoff volumes were initially made for watersheds in west-central Florida using recommended procedures, then compared to observed peak discharges and runoff volumes. Subsequently, the procedures were modified to increase accuracy for this area. The same methods used during the study could be used in other parts of the world to evaluate the accuracy of standard methods for estimating stormwater runoff.

Florida↗
USGS2019

Xenon hydrate as an analogue of methane hydrate in geologic systems out of thermodynamic equilibrium

Methane hydrate occurs naturally under pressure and temperature conditions that are not straightforward to replicate experimentally. Xenon has emerged as an attractive laboratory alternative to methane for studying hydrate formation and dissociation in multiphase systems, given that it forms hydrates under milder conditions. However, building reliable analogies between the two hydrates requires systematic comparisons, which are currently lacking. We address this gap by developing a theoretical and computational model of gas hydrates under equilibrium and nonequilibrium conditions. We first compare equilibrium phase behaviors of the Xe·H 2 O and CH 4 ·H 2 O systems by calculating their isobaric phase diagram, and then study the nonequilibrium kinetics of interfacial hydrate growth using a phase field model. Our results show that Xe·H 2 O is a good experimental analog to CH 4 ·H 2 O, but there are key differences to consider. In particular, the aqueous solubility of xenon is altered by the presence of hydrate, similar to what is observed for methane; but xenon is consistently less soluble than methane. Xenon hydrate has a wider nonstoichiometry region, which could lead to a thicker hydrate layer at the gas‐liquid interface when grown under similar kinetic forcing conditions. For both systems, our numerical calculations reveal that hydrate nonstoichiometry coupled with hydrate formation dynamics leads to a compositional gradient across the hydrate layer, where the stoichiometric ratio increases from the gas‐facing side to the liquid‐facing side. Our analysis suggests that accurate composition measurements could be used to infer the kinetic history of hydrate formation in natural settings where gas is abundant.

Geochemistry, Geophysics, Geosystems↗
USGS2004

Development of a geodatabase and conceptual model of the hydrogeologic units beneath air force plant 4 and Naval Air Station-Joint Reserve Base Carswell Field, Fort Worth, Texas

Air Force Plant 4 and adjacent Naval Air Station-Joint Reserve Base Carswell Field at Fort Worth, Texas, constitute a government-owned, contractor-operated facility that has been in operation since 1942. Contaminants from AFP4, primarily volatile organic compounds and metals, have entered the ground-water-flow system through leakage from waste-disposal sites and from manufacturing processes. The U.S. Geological Survey developed a comprehensive geodatabase of temporal and spatial environmental information associated with the hydrogeologic units (alluvial aquifer, Goodland-Walnut confining unit, and Paluxy aquifer) beneath the facility and a three-dimensional conceptual model of the hydrogeologic units integrally linked to the geodatabase. The geodatabase design uses a thematic layer approach to create layers of feature data using a geographic information system. The various features are separated into relational tables in the geodatabase on the basis of how they interact and correspond to one another. Using the geodatabase, geographic data at the site are manipulated to produce maps, allow interactive queries, and perform spatial analyses. The conceptual model for the study area comprises computer-generated, three-dimensional block diagrams of the hydrogeologic units. The conceptual model provides a platform for visualization of hydrogeologic-unit sections and surfaces and for subsurface environmental analyses. The conceptual model is based on three structural surfaces and two thickness configurations of the study area. The three structural surfaces depict the altitudes of the tops of the three hydrogeologic units. The two thickness configurations are those of the alluvial aquifer and the Goodland-Walnut confining unit. The surface of the alluvial aquifer was created using a U.S. Geological Survey 10-meter digital elevation model. The 2,130 point altitudes of the top of the Goodland-Walnut unit were compiled from lithologic logs from existing wells, available soil-boring logs, and previous studies. Data from 120 wells, primarily from existing reports, were used to create a map of the approximate altitude of the Paluxy aquifer.

Texas↗
USGS2013

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗
USGS1992

Sampling design for spatially distributed hydrogeologic and environmental processes

A methodology for the design of sampling networks over space is proposed. The methodology is based on spatial random field representations of nonhomogeneous natural processes, and on optimal spatial estimation techniques. One of the most important results of random field theory for physical sciences is its rationalization of correlations in spatial variability of natural processes. This correlation is extremely important both for interpreting spatially distributed observations and for predictive performance. The extent of site sampling and the types of data to be collected will depend on the relationship of subsurface variability to predictive uncertainty. While hypothesis formulation and initial identification of spatial variability characteristics are based on scientific understanding (such as knowledge of the physics of the underlying phenomena, geological interpretations, intuition and experience), the support offered by field data is statistically modelled. This model is not limited by the geometric nature of sampling and covers a wide range in subsurface uncertainties. A factorization scheme of the sampling error variance is derived, which possesses certain atttactive properties allowing significant savings in computations. By means of this scheme, a practical sampling design procedure providing suitable indices of the sampling error variance is established. These indices can be used by way of multiobjective decision criteria to obtain the best sampling strategy. Neither the actual implementation of the in-situ sampling nor the solution of the large spatial estimation systems of equations are necessary. The required values of the accuracy parameters involved in the network design are derived using reference charts (readily available for various combinations of data configurations and spatial variability parameters) and certain simple yet accurate analytical formulas. Insight is gained by applying the proposed sampling procedure to realistic examples related to sampling problems in two dimensions.

Advances in Water Resources↗
USGS2010

An efficient deterministic-probabilistic approach to modeling regional groundwater flow: 2. Application to Owens Valley, California

The applicability of a deterministic-probabilistic model for predicting water tables in southern Owens Valley, California, is evaluated. The model is based on a two-layer deterministic model that is cascaded with a two-point probability model. To reduce the potentially large number of uncertain variables in the deterministic model, lumping of uncertain variables was evaluated by sensitivity analysis to reduce the total number of uncertain variables to three variables: hydraulic conductivity, storage coefficient or specific yield, and source-sink function. Results demonstrate that lumping of uncertain parameters reduces computational effort while providing sufficient precision for the case studied. Simulated spatial coefficients of variation for water table temporal position in most of the basin is small, which suggests that deterministic models can predict water tables in these areas with good precision. However, in several important areas where pumping occurs or the geology is complex, the simulated spatial coefficients of variation are over estimated by the two-point probability method.

California↗
USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗