SEARCH · USGS Science

Results for “Modeling Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15Linked to original sources

USGS1981

Digital-computer model of the principal ground-water reservoir in Beryl-Enterprise area, Escalante Desert, Utah

The computer model presented in this report was used to simulate the principal ground-water reservoir in the Beryl-Enterprise area, Escalante Desert, Beaver, Iron, and Washington Counties, Utah (Mower, 1981). The details of the formulation of the model, testing of its validity, and the results of predictions are discussed in the cited report. This report was prepared as part of a cooperative program with the Utah Department of Natural Resources, Division of Water Rights, to investigate the water resources of the State. It is an addendum to the principal interpretive report, and it is presented in order to make the model available to anyone desiring to use it for additional predictions. The main program used was the finite-difference model for aquifer simulation in two dimensions documented by Trescott, Pinder, and Larson, (1976). Minor modifications were made to adapt the program to the principal ground-water reservoir in the Beryl-Enterprise area. All the modifications are listed at the top of table 1, and were related to parameter input and output, thus none of the computational subroutines were affected. The parameter arrays (table 1) and map of the area with a grid overlay (pi. 1) are given on following pages. The model simulates an aquifer- under water-table conditions, mostly composed of unconsoliuated basin-fill deposits. The boundaries of the modeled area (pi. 1) generally coincide with the boundaries of the saturated basin fill. However, in the southwest-central part of the model, permeable consolidated rock is included; and that part of the northern boundary between the Black and Wah Wah Mountains is an arbitrary boundary in basin fill between the Beryl-Enterprise area and the Milford area that lies to the northeast. The ignimbrite at Table Butte also was included in the active part of the model. The model includes simulation of discharge by evapotranspiration from phreatophytes. The areal recharge array was used to simulate recharge entering the modeled area at its boundaries and from stream infiltration in the southern corner near Enterprise. In addition, this array included discharge by wells operated during the period simulated as being under steady-state conditions (virtually 1937), and discharging wells simulating flow of water northeast to the Milford area. These wells also were included in the transient-state simulation (1937-77), although any changes in this discharge were modeled using the pumpage array (Group IV, table 1). The wells simulating outflow to the Milford area are shown on plate 1, but the wells pumping in 1937 are not shown unless they also were pumped during 1937-77. The pumpage array was used to simulate: (1) Discharge from wells, (2) discharge after 1977 from a mine in the southwest-central part of the model and recharge resulting form the mine discharge (pi. 1), and (3) changes in discharge in wells operated during the steady-state period. Recharge from irrigation was simulated by reducing pumpage from nodes where irrigation occurs. Discharge from all wells was reduced by 5 percent by multiplying all pumpage by 0.95 in the computer program. North of Newcastle, in T. 35 S., R. 15 W., pumpage was reduced by 35 percent because surface materials are very permeable.

Utah↗

USGS1998

Hydrogeology and simulation of ground-water flow in the Paluxy aquifer in the vicinity of Landfills 1 and 3, U.S. Air Force Plant 4, Fort Worth, Texas

Ground-water contamination of the surficial terrace alluvial aquifer has occurred at U.S. Air Force Plant 4, a government-owned, contractor-operated facility, northwest of Fort Worth, Texas. A poorly constructed monitoring well, P–22M, open to the underlying middle zone of the Paluxy aquifer was installed at landfill 3, October 1987, allowing leakage of contaminated ground water to reach the Paluxy aquifer. This well was plugged and abandoned in November 1995. Additionally, volatile organic compounds have been detected in fractures in the Goodland-Walnut confining unit, the hydrogeologic unit separating the terrace alluvial aquifer from the underlying Paluxy aquifer, beneath the western part of landfill 1. Volatile organic compounds in concentrations near the analytical detection limit were detected in the upper Paluxy prior to the drilling of well P–22M. The ground-water-flow simulation model described in this report was developed to examine the best logistically feasible location to install recovery wells to capture the low concentration (less than 100 micrograms per liter) trichloroethylene plume beneath landfills 1 and 3 (west Paluxy plume). Once the recovery wells were installed (1996), the simulation model was recalibrated with new data. This report documents the capture area of the installed recovery wells. Four geologic units are pertinent to this site-specific model. From oldest to youngest, these are the Glen Rose Formation, Paluxy Formation, Walnut Formation, and Goodland Limestone. The Glen Rose Formation is relatively impermeable in the study area and forms the confining unit underlying the Paluxy Formation. The Paluxy Formation forms the Paluxy aquifer, which is a public drinking water supply for the City of White Settlement. The Walnut Formation and Goodland Limestone form the Goodland-Walnut confining unit overlying the Paluxy aquifer. Near landfill 3, gamma-ray logs indicate three distinct zones of the Paluxy Formation; upper, middle, and lower. The formation is about 170-feet thick near landfill 3, and each zone is about 57-feet thick. Two steady-state simulations using the computer program MODFLOW were analyzed using the particle-tracking computer program, MODPATH. One simulation is the calibration simulation using Paluxy aquifer water-level data for May 1993. The second simulation includes the installed recovery wells. A variably spaced grid was designed for the model. The smallest grid cells, 25 by 25 feet, are in the vicinity of landfills 1 and 3. The largest cells, 4,864.5 by 1,441.5 feet, are at the northwestern corner of the model grid near the Parker-Tarrant County line. The modeling was accomplished with three layers representing the upper, middle, and lower zones of the Paluxy aquifer. Particles, which represent contaminant molecules moving in solution with the ground water, were tracked from well P–22M and an area below landfill 1, at the top of the upper zone of the Paluxy aquifer, for 9 years (forward tracking). The forward tracking estimates where contaminants might move by advection from 1987 to 1996. Analysis of backward tracking from the new recovery wells indicates that the simulated contributing area to the recovery wells intercepts the contaminant plume, minimizing off-site migration of the west Paluxy plume. To determine the effectiveness of the recovery wells, monitoring wells southeast of Building 14 have been installed (1996–97) for sampling.

Texas↗

USGS1992

Modifications of the U.S. Geological Survey modular, finite-difference, ground-water flow model to read and write geographic information system files

This report documents modifications to the U.S. Geological Survey modular, three-dimensional, finite-difference, ground-water flow model, commonly called MODFLOW, so that it can read and write files used by a geographic information system (GIS). The modified model program is called MODFLOWARC. Simulation programs such as MODFLOW generally require large amounts of input data and produce large amounts of output data. Viewing data graphically, generating head contours, and creating or editing model data arrays such as hydraulic conductivity are examples of tasks that currently are performed either by the use of independent software packages or by tedious manual editing, manipulating, and transferring data. Programs such as GIS programs are commonly used to facilitate preparation of the model input data and analyze model output data; however, auxiliary programs are frequently required to translate data between programs. Data translations are required when different programs use different data formats. Thus, the user might use GIS techniques to create model input data, run a translation program to convert input data into a format compatible with the ground-water flow model, run the model, run a translation program to convert the model output into the correct format for GIS, and use GIS to display and analyze this output. MODFLOWARC, avoids the two translation steps and transfers data directly to and from the ground-water-flow model. This report documents the design and use of MODFLOWARC and includes instructions for data input/output of the Basic, Block-centered flow, River, Recharge, Well, Drain, Evapotranspiration, General-head boundary, and Streamflow-routing packages. The modification to MODFLOW and the Streamflow-Routing package was minimized. Flow charts and computer-program code describe the modifications to the original computer codes for each of these packages. Appendix A contains a discussion on the operation of MODFLOWARC using a sample problem.

Open-File Report↗

USGS2002

User's guide to Model Viewer, a program for three-dimensional visualization of ground-water model results

Model Viewer is a computer program that displays the results of three-dimensional groundwater models. Scalar data (such as hydraulic head or solute concentration) may be displayed as a solid or a set of isosurfaces, using a red-to-blue color spectrum to represent a range of scalar values. Vector data (such as velocity or specific discharge) are represented by lines oriented to the vector direction and scaled to the vector magnitude. Model Viewer can also display pathlines, cells or nodes that represent model features such as streams and wells, and auxiliary graphic objects such as grid lines and coordinate axes. Users may crop the model grid in different orientations to examine the interior structure of the data. For transient simulations, Model Viewer can animate the time evolution of the simulated quantities. The current version (1.0) of Model Viewer runs on Microsoft Windows 95, 98, NT and 2000 operating systems, and supports the following models: MODFLOW-2000, MODFLOW-2000 with the Ground-Water Transport Process, MODFLOW-96, MOC3D (Version 3.5), MODPATH, MT3DMS, and SUTRA (Version 2D3D.1). Model Viewer is designed to directly read input and output files from these models, thus minimizing the need for additional postprocessing. This report provides an overview of Model Viewer. Complete instructions on how to use the software are provided in the on-line help pages.

Open-File Report↗

USGS2014

Stochastic modeling of a lava-flow aquifer system

This report describes preliminary three-dimensional geostatistical modeling of a lava-flow aquifer system using a multiple-point geostatistical model. The purpose of this study is to provide a proof-of-concept for this modeling approach. An example of the method is demonstrated using a subset of borehole geologic data and aquifer test data from a portion of the Calico Hills Formation, a lava-flow aquifer system that partially underlies Pahute Mesa, Nevada. Groundwater movement in this aquifer system is assumed to be controlled by the spatial distribution of two geologic units—rhyolite lava flows and zeolitized tuffs. The configuration of subsurface lava flows and tuffs is largely unknown because of limited data. The spatial configuration of the lava flows and tuffs is modeled by using a multiple-point geostatistical simulation algorithm that generates a large number of alternative realizations, each honoring the available geologic data and drawn from a geologic conceptual model of the lava-flow aquifer system as represented by a training image. In order to demonstrate how results from the geostatistical model could be analyzed in terms of available hydrologic data, a numerical simulation of part of an aquifer test was applied to the realizations of the geostatistical model.

Nevada↗

USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗

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↗

USGS1997

Vulnerability of ground water to atrazine leaching in Kent County, Michigan

A steady-state model of pesticide leaching through the unsaturated zone was used with readily available hydrologic, lithologic, and pesticide characteristics to estimate the vulnerability of the near-surface aquifer to atrazine contamination from non-point sources in Kent County, Michigan. The modelcomputed fraction of atrazine remaining at the water table, RM , was used as the vulnerability criterion; time of travel to the water table also was computed. Model results indicate that the average fraction of atrazine remaining at the water table was 0.039 percent; the fraction ranged from 0 to 3.6 percent. Time of travel of atrazine from the soil surface to the water table averaged 17.7 years and ranged from 2.2 to 118 years. Three maps were generated to present three views of the same atrazine vulnerability characteristics using different metrics (nonlinear transformations of the computed fractions remaining). The metrics were chosen because of the highly (right) skewed distribution of computed fractions. The first metric, rm = RM λ (where λ was 0.0625), depicts a relatively uniform distribution of vulnerability across the county with localized areas of high and low vulnerability visible. The second metric, rm λ-0.5 , depicts about one-half the county at low vulnerability with discontinuous patterns of high vulnerability evident. In the third metric, rm λ-1.0 ( RM ), more than 95 percent of the county appears to have low vulnerability; small, distinct areas of high vulnerability are present. Aquifer vulnerability estimates in the RM metric were used with a steady-state, uniform atrazine application rate to compute a potential concentration of atrazine in leachate reaching the water table. The average estimated potential atrazine concentration in leachate at the water table was 0.16 μg/L (micrograms per liter) in the model area; estimated potential concentrations ranged from 0 to 26 μg/L. About 2 percent of the model area had estimated potential atrazine concentrations in leachate at the water table that exceeded the USEPA (U.S. Environmental Protection Agency) maximum contaminant level of 3 μg/L. Uncertainty analyses were used to assess effects of parameter uncertainty and spatial interpolation error on the variability of the estimated fractions of atrazine remaining at the water table. Results of Monte Carlo simulations indicate that parameter uncertainty is associated with a standard error of 0.0875 in the computed fractions (in the rm metric). Results of kriging analysis indicate that errors in spatial interpolation are associated with a standard error of 0.146 (in the rm metric). Thus, uncertainty in fractions remaining is primarily associated with spatial interpolation error, which can be reduced by increasing the density of points where the leaching model is applied. A sensitivity analysis indicated which of 13 hydrologic, lithologic, and pesticide characteristics were influential in determining fractions of atrazine remaining at the water table. Results indicate that fractions remaining are most sensitive to the unit changes in pesticide half life and in organic-carbon content in soils and unweathered rocks, and least sensitive to infiltration rates. The leaching model applied in this report provides an estimate of the vulnerability of the near-surface aquifer in Kent County to contamination by atrazine. The vulnerability estimate is related to water-quality criteria developed by the USEPA to help assess potential risks from atrazine to the near-surface aquifer. However, atrazine accounts for only 28 percent of the herbicide use in the county; additional potential for contamination exists from other pesticides and pesticide metabolites. Therefore, additional work is needed to develop a comprehensive understanding of the relative risks associated with specific pesticides. The modeling approach described in this report provides a technique for estimating relative vulnerabilities to specific pesticides and for helping to assess potential risks.

Michigan↗

USGS2015

Development of a HEC-RAS temperature model for the North Santiam River, northwestern Oregon

A one-dimensional, unsteady streamflow and temperature model (HEC-RAS) of the North Santiam and Santiam Rivers was developed by the U.S. Geological Survey to be used in conjunction with previously developed two-dimensional hydrodynamic water-quality models (CE-QUAL-W2) of Detroit and Big Cliff Lakes upstream of the study area. In conjunction with the output from the previously developed models, the HEC-RAS model can simulate streamflows and temperatures within acceptable limits (mean error [bias] near zero; typical streamflow errors less than 5 percent; typical water temperature errors less than 1.0 °C) for the length of the North Santiam River downstream of Big Cliff Dam under a series of potential future conditions in which dam structures and/or dam operations are modified to improve temperature conditions for threatened and endangered fish. Although a two-dimensional (longitudinal, vertical) CE-QUAL-W2 model for the North Santiam and Santiam Rivers downstream of Big Cliff Dam exists, that model proved unstable under highly variable flow conditions. The one-dimensional HEC-RAS model documented in this report can better simulate cross-sectional-averaged stream temperatures under a wide range of flow conditions. The model was calibrated using 2011 streamflow and temperature data. Measured data were used as boundary conditions when possible, although several lateral inflows and their associated water temperatures, including the South Santiam River, were estimated using statistical models. Streamflow results showed high accuracy during low-flow periods, but predictions were biased low during large storm events when unmodeled ephemeral tributaries contributed to the actual streamflow. Temperature results showed low annual bias against measured data at two locations on the North Santiam River and one location on the Santiam River. Mean absolute errors using 2011 hourly data ranged from 0.4 to 0.7 °C. Model results were checked against 2012 data and showed a positive bias at the Santiam River station (+0.6 ˚C). Annual mean absolute errors using 2012 hourly data ranged from 0.4 to 0.8 °C. Much of the error in temperature predictions resulted from the model’s inability to accurately simulate the full range of diurnal fluctuations during the warmest months. Future iterations of the model could be improved by the collection and inclusion of additional streamflow and temperature data, especially near the mouth of the South Santiam River. Presently, the model is able to predict hourly and daily water temperatures under a wide variety of conditions with a typical error of 0.8 and 0.7 °C, respectively.

Oregon↗

USGS2021

American and Sacramento Rivers, California, erodibility measurements and model

Executive Summary A previous report by the authors described sediment sampling and drilling by the U.S. Geological Survey (USGS) beside the American and Sacramento Rivers near Sacramento, California, in support of a U.S. Army Corps of Engineers project focused on regional flood control. The drilling was performed to define lithology, extract samples for laboratory testing, and perform borehole erosion tests (BETs). The U.S. Department of Agriculture (USDA) performed jet erodibility tests (JETs) near each drilling site, and a team from Texas A&M University performed laboratory tests with an erosion function apparatus (EFA). Collectively, the effort was intended to reveal spatial variations in sediment erodibility and provide data for use in a model to simulate morphological response to a major flood. The data collected by the USGS are available in a public data release. This report, developed in cooperation with the U.S. Army Corps of Engineers, provides comparisons of the three types of measurements of the erodibility of riverbed sediments. The BET is performed in the field and reveals erodibility of sediments below the bed surface. The JET is likewise performed in the field but reveals only erodibility of exposed sediments. The EFA test is done in the laboratory and was performed on soils extracted from different depths beneath the bed surface, in many cases reconstituted for laboratory testing. Tests were performed at nominally similar locations but differed by meters to tens of meters in horizontal locations. The comparison was undertaken to investigate differences among results obtained by the individual measurement approaches and to elucidate pros and cons of each method. The critical shear stress to initiate erosion and the rate of change of erosion rate per unit increase of excess shear stress, sometimes referred to as the erosion coefficient, served as the primary basis for comparison. The three test methods in some cases resulted in order of magnitude differences in estimates of these parameters. Some differences could be attributed to variances in site location or result from testing surface sediment versus a deeper layer, but systematic differences are also evident in the results. The tests performed in the laboratory using the EFA resulted in much lower values of critical shear stress and much higher values of the erosion coefficient compared to the JET tests performed by the USDA team on surface sediments. Critical shear stress was poorly resolved in the BET results because of the limited number of results per site, but the erosion coefficients derived from BET results were systematically lower than those obtained using the EFA. A new, simplified approach is also proposed to estimate the increase in channel cross-sectional area during a large flood, given data describing the initial river cross section, riverbed erodibility parameters, and peak flood discharge and duration. The model runs until the cross section erodes to an equilibrium condition or the flood ends. Output describes the area of the cross section at the end of the simulation and the time required to reach equilibrium if it was reached within the simulated period. The model assumes unique, constant values for both the critical shear stress and the erosion coefficient and represents the fluid mechanics in a simplified way, making it of limited value for quantitative predictions. It does, however, provide an indication of which cross sections are most likely to undergo the greatest change in the design event and can be used to investigate sensitivity of erosion predictions to variability in sediment erodibility measurements.

California↗

USGS2021

The basin characterization model—A regional water balance software package

This report documents the computer software package, Basin Characterization Model, version 8 (BCMv8)—a monthly, gridded, regional water-balance model—and provides detailed operational instructions and example applications. After several years of many applications and uses of a previous version, CA-BCM, published in 2014, the BCMv8 was refined to improve the accuracy of the water-balance components, particularly the recharge estimate, which is the most difficult to accurately assess. The improvement of the various water-balance components targeted the actual evapotranspiration component, which, in turn, reduced the uncertainty of the recharge estimate. The improvement of this component was enabled by the availability of a national, gridded actual-evapotranspiration product from the U.S. Geological Survey that was unique in its scope to combine remotely sensed spatial variability and ground-based long-term water-balance constraints. This dataset provided the ability to assess monthly actual evapotranspiration for 62 vegetation types and to perform regional calibration in watersheds throughout California with the objective of closing the water balance using improved estimates for each component. The refinements, including vegetation-specific evapotranspiration, enabled the development of applications that could explore various aspects of landscape disturbance, such as wildfire, forest management, or urbanization. The improvements to BCMv8 also provided the ability to assess long-term sustainability of water resources under a variety of management applications or future climate projections.

California↗

USGS2016

Groundwater-flow model for the Wood River Valley aquifer system, south-central Idaho

A three-dimensional numerical model of groundwater flow was developed for the Wood River Valley (WRV) aquifer system, Idaho, to evaluate groundwater and surface-water availability at the regional scale. This mountain valley is located in Blaine County and has a drainage area of about 2,300 square kilometers (888 square miles). The model described in this report can serve as a tool for water-rights administration and water-resource management and planning. The model was completed with support from the Idaho Department of Water Resources, and is part of an ongoing U.S. Geological Survey effort to characterize the groundwater resources of the WRV. A highly reproducible approach was taken for constructing the WRV groundwater-flow model. The collection of datasets, source code, and processing instructions used to construct and analyze the model was distributed as an R statistical-computing and graphics package. Flow in the WRV aquifer was simulated using the MODFLOW-USG groundwater flow model. The transient flow model simulates groundwater flow between 1995 and 2010. The model uses a 100-meter (328-feet) uniform grid spacing with 54,922 active model cells distributed over three model layers. A confining unit in the south-central part of the Bellevue fan necessitated the use of a multi-layer model. Specified-flow boundaries were used to simulate the groundwater inflows from each of the major tributary basins (also known as tributary basin underflow) and the areal recharge of precipitation and applied irrigation. Head‑dependent flow boundaries were used to simulate the stream-aquifer flow exchange in river reaches and the groundwater discharge at the outlet boundaries of Stanton Crossing and Silver Creek. The model was calibrated by adjusting aquifer hydraulic properties to match simulated and measured water levels and stream-aquifer flow exchange, using the parameter-estimation program PEST. The model reasonably simulated the measured water-table elevation, orientation, and gradients. Stream-aquifer flow exchange along river reaches also was reasonably simulated by the model. Inflow into the WRV aquifer system originates from three sources (from largest to smallest): Streamflow loss from the Big Wood River and Silver Creek, Areal recharge of precipitation and applied irrigation, and Tributary basin underflow. Outflow from the WRV aquifer system originates from five sources (from largest to smallest): Aquifer discharge into the Big Wood River and Silver Creek, Production-well pumping, Discharge across the water table into the vadose zone, Subsurface outflow beneath Silver Creek near Picabo, and Subsurface outflow beneath the Big Wood River near Stanton Crossing. Temporal changes in aquifer storage are most affected by areal recharge and groundwater pumping, and also contribute to changes in streamflow gains.

Idaho↗

USGS2016

Review and synthesis: Changing permafrost in a warming world and feedbacks to the Earth System

The permafrost component of the cryosphere is changing dramatically, but the permafrost region is not well monitored and the consequences of change are not well understood. Changing permafrost interacts with ecosystems and climate on various spatial and temporal scales. The feedbacks resulting from these interactions range from local impacts on topography, hydrology, and biology to complex influences on global scale biogeochemical cycling. This review contributes to this focus issue by synthesizing its 28 multidisciplinary studies which provide field evidence, remote sensing observations, and modeling results on various scales. We synthesize study results from a diverse range of permafrost landscapes and ecosystems by reporting key observations and modeling outcomes for permafrost thaw dynamics, identifying feedbacks between permafrost and ecosystem processes, and highlighting biogeochemical feedbacks from permafrost thaw. We complete our synthesis by discussing the progress made, stressing remaining challenges and knowledge gaps, and providing an outlook on future needs and research opportunities in the study of permafrost–ecosystem–climate interactions.

Environmental Research Letters↗

USGS2010

Building hierarchical models of avian distributions for the State of Georgia

To predict the distributions of breeding birds in the state of Georgia, USA, we built hierarchical models consisting of 4 levels of nested mapping units of decreasing area: 90,000 ha, 3,600 ha, 144 ha, and 5.76 ha. We used the Partners in Flight database of point counts to generate presence and absence data at locations across the state of Georgia for 9 avian species: Acadian flycatcher (Empidonax virescens), brownheaded nuthatch (Sitta pusilla), Carolina wren (Thryothorus ludovicianus), indigo bunting (Passerina cyanea), northern cardinal (Cardinalis cardinalis), prairie warbler (Dendroica discolor), yellow-billed cuckoo (Coccyxus americanus), white-eyed vireo (Vireo griseus), and wood thrush (Hylocichla mustelina). At each location, we estimated hierarchical-level-specific habitat measurements using the Georgia GAP Analysis18 class land cover and other Geographic Information System sources. We created candidate, species-specific occupancy models based on previously reported relationships, and fit these using Markov chain Monte Carlo procedures implemented in OpenBugs. We then created a confidence model set for each species based on Akaike's Information Criterion. We found hierarchical habitat relationships for all species. Three-fold cross-validation estimates of model accuracy indicated an average overall correct classification rate of 60.5%. Comparisons with existing Georgia GAP Analysis models indicated that our models were more accurate overall. Our results provide guidance to wildlife scientists and managers seeking predict avian occurrence as a function of local and landscape-level habitat attributes.

Journal of Wildlife Management↗

USGS2008

SIM_ADJUST -- A computer code that adjusts simulated equivalents for observations or predictions

This report documents the SIM_ADJUST computer code. SIM_ADJUST surmounts an obstacle that is sometimes encountered when using universal model analysis computer codes such as UCODE_2005 (Poeter and others, 2005), PEST (Doherty, 2004), and OSTRICH (Matott, 2005; Fredrick and others (2007). These codes often read simulated equivalents from a list in a file produced by a process model such as MODFLOW that represents a system of interest. At times values needed by the universal code are missing or assigned default values because the process model could not produce a useful solution. SIM_ADJUST can be used to (1) read a file that lists expected observation or prediction names and possible alternatives for the simulated values; (2) read a file produced by a process model that contains space or tab delimited columns, including a column of simulated values and a column of related observation or prediction names; (3) identify observations or predictions that have been omitted or assigned a default value by the process model; and (4) produce an adjusted file that contains a column of simulated values and a column of associated observation or prediction names. The user may provide alternatives that are constant values or that are alternative simulated values. The user may also provide a sequence of alternatives. For example, the heads from a series of cells may be specified to ensure that a meaningful value is available to compare with an observation located in a cell that may become dry. SIM_ADJUST is constructed using modules from the JUPITER API, and is intended for use on any computer operating system. SIM_ADJUST consists of algorithms programmed in Fortran90, which efficiently performs numerical calculations.

Ground Water Modeling Investigation Report↗

USGS2011

A digital model for planning water management at Benton Lake National Wildlife Refuge, west-central Montana

Benton Lake National Wildlife Refuge is an important area for waterfowl production and migratory stopover in west-central Montana. Eight wetland units covering about 5,600 acres are the essential features of the refuge. Water availability for the wetland units can be uncertain owing to the large natural variations in precipitation and runoff and the high cost of pumping supplemental water. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service, has developed a digital model for planning water management. The model can simulate strategies for water transfers among the eight wetland units and account for variability in runoff and pumped water. This report describes this digital model, which uses a water-accounting spreadsheet to track inputs and outputs to each of the wetland units of Benton Lake National Wildlife Refuge. Inputs to the model include (1) monthly values for precipitation, pumped water, runoff, and evaporation; (2) water-level/capacity data for each wetland unit; and (3) the pan-evaporation coefficient. Outputs include monthly water volume and flooded surface area for each unit for as many as 5 consecutive years. The digital model was calibrated by comparing simulated and historical measured water volumes for specific test years.

Montana↗

USGS2022

Hydrologic and geochemical characterization of the Petaluma River watershed, Sonoma County, California

Executive Summary The objectives of the study are to (1) develop an updated assessment of the hydrogeology and geochemistry of the Petaluma valley watershed (PVW) and (2) develop an integrated hydrologic model for the PVW. The purpose of this report is to describe the conceptual model of the hydrologic, hydrogeologic, and water-quality characteristics of the PVW and a numerical groundwater-flow model of PVW.

California↗

USGS2007

An initial-abstraction, constant-loss model for unit hydrograph modeling for applicable watersheds in Texas

Estimation of representative hydrographs from design storms, which are known as design hydrographs, provides for cost-effective, riskmitigated design of drainage structures such as bridges, culverts, roadways, and other infrastructure. During 2001?07, the U.S. Geological Survey (USGS), in cooperation with the Texas Department of Transportation, investigated runoff hydrographs, design storms, unit hydrographs,and watershed-loss models to enhance design hydrograph estimation in Texas. Design hydrographs ideally should mimic the general volume, peak, and shape of observed runoff hydrographs. Design hydrographs commonly are estimated in part by unit hydrographs. A unit hydrograph is defined as the runoff hydrograph that results from a unit pulse of excess rainfall uniformly distributed over the watershed at a constant rate for a specific duration. A time-distributed, watershed-loss model is required for modeling by unit hydrographs. This report develops a specific time-distributed, watershed-loss model known as an initial-abstraction, constant-loss model. For this watershed-loss model, a watershed is conceptualized to have the capacity to store or abstract an absolute depth of rainfall at and near the beginning of a storm. Depths of total rainfall less than this initial abstraction do not produce runoff. The watershed also is conceptualized to have the capacity to remove rainfall at a constant rate (loss) after the initial abstraction is satisfied. Additional rainfall inputs after the initial abstraction is satisfied contribute to runoff if the rainfall rate (intensity) is larger than the constant loss. The initial abstraction, constant-loss model thus is a two-parameter model. The initial-abstraction, constant-loss model is investigated through detailed computational and statistical analysis of observed rainfall and runoff data for 92 USGS streamflow-gaging stations (watersheds) in Texas with contributing drainage areas from 0.26 to 166 square miles. The analysis is limited to a previously described, watershed-specific, gamma distribution model of the unit hydrograph. In particular, the initial-abstraction, constant-loss model is tuned to the gamma distribution model of the unit hydrograph. A complex computational analysis of observed rainfall and runoff for the 92 watersheds was done to determine, by storm, optimal values of initial abstraction and constant loss. Optimal parameter values for a given storm were defined as those values that produced a modeled runoff hydrograph with volume equal to the observed runoff hydrograph and also minimized the residual sum of squares of the two hydrographs. Subsequently, the means of the optimal parameters were computed on a watershed-specific basis. These means for each watershed are considered the most representative, are tabulated, and are used in further statistical analyses. Statistical analyses of watershed-specific, initial abstraction and constant loss include documentation of the distribution of each parameter using the generalized lambda distribution. The analyses show that watershed development has substantial influence on initial abstraction and limited influence on constant loss. The means and medians of the 92 watershed-specific parameters are tabulated with respect to watershed development; although they have considerable uncertainty, these parameters can be used for parameter prediction for ungaged watersheds. The statistical analyses of watershed-specific, initial abstraction and constant loss also include development of predictive procedures for estimation of each parameter for ungaged watersheds. Both regression equations and regression trees for estimation of initial abstraction and constant loss are provided. The watershed characteristics included in the regression analyses are (1) main-channel length, (2) a binary factor representing watershed development, (3) a binary factor representing watersheds with an abundance of rocky and thin-soiled terrain, and (4) curve numb

Scientific Investigations Report↗