Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Modeling Report”

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 271 records · Page 15Linked to original sources
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↗
USGS2001

River and Reservoir Operations Model, Truckee River basin, California and Nevada, 1998

The demand for all uses of water in the Truckee River Basin, California and Nevada, commonly is greater than can be supplied. Storage reservoirs in the system have a maximum effective total capacity equivalent to less than two years of average river flows, so longer-term droughts can result in substantial water-supply shortages for irrigation and municipal users and may stress fish and wildlife ecosystems. Title II of Public Law (P.L.) 101-618, the Truckee?Carson?Pyramid Lake Water Rights Settlement Act of 1990, provides a foundation for negotiating and developing operating criteria, known as the Truckee River Operating Agreement (TROA), to balance interstate and interbasin allocation of water rights among the many interests competing for water from the Truckee River. In addition to TROA, the Truckee River Water Quality Settlement Agreement (WQSA), signed in 1996, provides for acquisition of water rights to resolve water-quality problems during low flows along the Truckee River in Nevada. Efficient execution of many of the planning, management, or environmental assessment requirements of TROA and WQSA will require detailed water-resources data coupled with sound analytical tools. Analytical modeling tools constructed and evaluated with such data could help assess effects of alternative operational scenarios related to reservoir and river operations, water-rights transfers, and changes in irrigation practices. The Truckee?Carson Program of the U.S. Geological Survey, to support U.S. Department of the Interior implementation of P.L. 101-618, is developing a modeling system to support efficient water-resources planning, management, and allocation. The daily operations model documented herein is a part of the modeling system that includes a database management program, a graphical user interface program, and a program with modules that simulate river/reservoir operations and a variety of hydrologic processes. The operations module is capable of simulating lake/ reservoir and river operations including diversion of Truckee River water to the Truckee Canal for transport to the Carson River Basin. In addition to the operations and streamflow-routing modules, the modeling system is structured to allow integration of other modules, such as water-quality and precipitation-runoff modules. The USGS Truckee River Basin operations model was designed to provide simulations that allow comparison of the effects of alternative management practices or allocations on streamflow or reservoir storages in the Truckee River Basin over long periods of time. Because the model was not intended to reproduce historical streamflow or reservoir storage values, a traditional calibration that includes statistical comparisons of observed and simulated values would be problematic with this model and database. This report describes a chronology and background of decrees, agreements, and laws that affect Truckee River operational practices; the construction of the Truckee River daily operations model; the simulation of Truckee River Basin operations, both current and proposed under the draft TROA and WQSA; and suggested model improvements and limitations. The daily operations model uses Hydrological Simulation Program?FORTRAN (HSPF) to simulate flow-routing and reservoir and river operations. The operations model simulates reservoir and river operations that govern streamflow in the Truckee River from Lake Tahoe to Pyramid Lake, including diversions through the Truckee Canal to Lahontan Reservoir in the Carson River Basin. A general overview is provided of daily operations and their simulation. Supplemental information that documents the extremely complex operating rules simulated by the model is available.

Water-Resources Investigations Report↗
USGS2022

Methods for evaluating Gap Analysis Project habitat distribution maps with species occurrence data

The National Gap Analysis Project created species habitat distribution models for all terrestrial vertebrates in the United States to support conservation assessments and explore patterns of species richness. Those models link species to specific habitats throughout the range of each species. For most vertebrates, there are not enough occurrence data to drive inductive, range-wide species habitat distribution models at high spatial and thematic resolution. However, it is possible to use occurrence data for model evaluation. The combination of citizen science, formal species survey work, and digitized specimen archives are making millions of observations available to the scientific community. Our challenge is to combine the mostly unstructured data into metrics that help us characterize and understand patterns of biodiversity. In this work, we propose two model-evaluation metrics. The first, a buffer proportion assessment, is based on the proportion of habitat in the range relative to the mean proportion of habitat around each of the species’ occurrence records. The second is a measure of the sensitivity (proportion of true presence) to buffer distances around occurrence records. The buffer proportion is a modification of model prevalence versus point prevalence metric, whereby comparison to a null model allows us to determine if the model performs better or worse than random. In this report, we describe the workflow used to compile and filter the species occurrence records from online resources (for example, the Global Biodiversity Information Facility) and show results for a single species, Desmognathus quadramaculatus (black-bellied salamander). For the salamander, 222 occurrence points met our criteria for inclusion in the evaluation. We found the model performed better than random with a buffer proportion index of 1.745, indicating about 5 times as much habitat was found adjacent to known occurrence records than would be expected from randomly located sites throughout the range. Sensitivity increased with larger buffer distances and leveled off to around 0.7 between 1,000- and 2,000-meter buffer distances, indicating the model is likely best suited for scales exceeding 1,000 meters. We plan to report the buffer proportion assessment and sensitivity metrics along with the full species model reports to increase understanding of the model’s performance and to use the metrics to help prioritize revisions to the models.

Alabama, Georgia, North Carolina, Tennessee, Virgi↗