Search USGS⌕ Search

SEARCH · Search USGS

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 577 records · Page 32Linked to original sources
USGS2025

Climate change impacts on plant communities in the sagebrush region—A science synthesis to inform Bureau of Land Management resource management

This report synthesizes current (2024) science-based knowledge related to the impacts of climate change on big sagebrush vegetation in Western North America. This effort was conducted through the U.S. Geological Survey working with the Bureau of Land Management as part of multiple science syntheses to aid management agencies developing environmental impacts assessments in response to human-related or caused events. This report reviews the potential impacts climate change may have on sagebrush vegetation and related management decisions. The body of the synthesis introduces the diverse impacts of climate change across the region by first focusing directly on what climate change may entail in terms of altered temperature and precipitation patterns. The report then discusses how these changes could likely affect individual plant species based on experimental results and scale the impacts to species distributions and community composition. The synthesis section ends by surveying efforts to model potential future changes in habitat. The report goes on to link the synthesis conclusions to individual land uses or land management decisions, such as forage resources, restoration or fuel treatment. Finally, the report provides a section that discusses the pros and cons of available datasets that model the potential future of vegetation in the sagebrush region.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2018

Simulation of potential groundwater recharge for the glacial aquifer system east of the Rocky Mountains, 1980–2011, using the Soil-Water-Balance Model

An understanding of the spatial and temporal extent of groundwater recharge is critical for many types of hydrologic assessments involving water quality, contaminant transport, ecosystem health, and sustainable use of groundwater. Annual potential groundwater recharge was simulated at a 1-kilometer resolution with the Soil-Water-Balance (SWB) model for the glacial aquifer system east of the Rocky Mountains, from central Montana east to Maine, for calendar years 1980–2011. The SWB model used high resolution meteorological, land cover, and soil hydrology datasets that are nationally consistent and publicly available. The SWB model computed daily potential groundwater recharge as precipitation in excess of interception, runoff, evapotranspiration, and soil-water storage capacity. Daily potential recharge values within each year of the simulation were summed to produce annual potential recharge rates. Potential recharge as described in this report is water that infiltrates vertically below the plant rooting zone and is assumed to reach the water table. The calibrated SWB model in this report is called the glacial SWB model. Model calibration assumed that the area contributing to groundwater discharge equaled the surface watershed. The model was calibrated to stream base flows from 39 watersheds throughout the model domain that had hydrologic conditions appropriate for hydrograph separation. Base flows were calculated from daily streamflow records with the HYSEP local minimum hydrograph separation method The glacial SWB model reproduced the mean annual base-flow calibration targets well; the Nash-Sutcliffe efficiency coefficient was 0.94, and the root mean squared error was 1.28 inches per year. The glacial SWB model provides insight into the spatial and temporal variability in potential annual recharge across the glacial aquifer system. About 20 percent of the active model area had an average potential recharge rate of less than 1 inch per year. Total precipitation, total recharge, and recharge as a percentage of precipitation increased from west to east. A substantial amount of the recharge water (39 percent) entering the glacial aquifer system travels through developed (urbanized) and agricultural landscapes, which are known to cause water-quality impairments. Regional climatic events, such as the 1988 to 1989 drought, are apparent in the potential recharge time series. Potential recharge generally increased across the glacial aquifer system between 2001 and 2011. A comparison of the potential recharge from the glacial SWB model to previous broad-scale recharge estimates reveals several important considerations for future SWB modeling applications. Shifts in the overall distribution of potential recharge between separate models can be explained by methods used to generate base-flow calibration target datasets. Spatial patterns in potential recharge simulated by SWB models are strongly dependent on the data and assumptions used to assign model cells to hydrologic soil groups. A review of several SWB models used to estimate groundwater recharge (and not surface runoff) revealed that model results are most sensitive to input climatic data, followed by surface runoff (curve number) and root-zone depth parameters.

Scientific Investigations Report↗
USGS2000

Documentation of a computer program to simulate lake-aquifer interaction using the MODFLOW ground water flow model and the MOC3D solute-transport model

Heads and flow patterns in surficial aquifers can be strongly influenced by the presence of stationary surface-water bodies (lakes) that are in direct contact, vertically and laterally, with the aquifer. Conversely, lake stages can be significantly affected by the volume of water that seeps through the lakebed that separates the lake from the aquifer. For these reasons, a set of computer subroutines called the Lake Package (LAK3) was developed to represent lake/aquifer interaction in numerical simulations using the U.S. Geological Survey three-dimensional, finite-difference, modular ground-water flow model MODFLOW and the U.S. Geological Survey three-dimensional method-of-characteristics solute-transport model MOC3D. In the Lake Package described in this report, a lake is represented as a volume of space within the model grid which consists of inactive cells extending downward from the upper surface of the grid. Active model grid cells bordering this space, representing the adjacent aquifer, exchange water with the lake at a rate determined by the relative heads and by conductances that are based on grid cell dimensions, hydraulic conductivities of the aquifer material, and user-specified leakance distributions that represent the resistance to flow through the material of the lakebed. Parts of the lake may become "dry" as upper layers of the model are dewatered, with a concomitant reduction in lake surface area, and may subsequently rewet when aquifer heads rise. An empirical approximation has been encoded to simulate the rewetting of a lake that becomes completely dry. The variations of lake stages are determined by independent water budgets computed for each lake in the model grid. This lake budget process makes the package a simulator of the response of lake stage to hydraulic stresses applied to the aquifer. Implementation of a lake water budget requires input of parameters including those representing the rate of lake atmospheric recharge and evaporation, overland runoff, and the rate of any direct withdrawal from, or augmentation of, the lake volume. The lake/aquifer interaction may be simulated in both transient and steady-state flow conditions, and the user may specify that lake stages be computed explicitly, semi-implicitly, or fully-implicitly in transient simulations. The lakes, and all sources of water entering the lakes, may have solute concentrations associated with them for use in solute-transport simulations using MOC3D. The Stream Package of MODFLOW-2000 and MOC3D represents stream connections to lakes, either as inflows or outflows. Because lakes with irregular bathymetry can exist as separate pools of water at lower stages, that coalesce to become a single body of water at higher stages, logic was added to the Lake Package to allow the representation of this process as a user option. If this option is selected, a system of linked pools (sublakes) is identified in each time step and stages are equalized based on current relative sublake surface areas.

Water-Resources Investigations Report↗
USGS2005

Hydrogeologic setting and conceptual hydrologic model of the Spring Creek Basin, Centre County, Pennsylvania, June 2005

The Spring Creek Basin, Centre County, Pa., is experiencing some of the most rapid growth and development within the Commonwealth. This trend has resulted in land-use changes and increased water use, which will affect the quantity and quality of stormwater runoff, surface water, ground water, and aquatic resources within the basin. The U.S. Geological Survey (USGS), in cooperation with the ClearWater Conservancy (CWC), Spring Creek Watershed Community (SCWC), and Spring Creek Watershed Commission (SCWCm), has developed a Watershed Plan (Plan) to assist decision makers in water-resources planning. One element of the Plan is to provide a summary of the basin characteristics and a conceptual model that incorporates the hydrogeologic characteristics of the basin. The report presents hydrogeologic data for the basin and presents a conceptual model that can be used as the basis for simulating surface-water and ground-water flow within the basin. Basin characteristics; sources of data referenced in this text; physical characteristics such as climate, physiography, topography, and land use; hydrogeologic characteristics; and water-quality characteristics are discussed. A conceptual model is a simplified description of the physical components and interaction of the surface- and ground-water systems. The purpose for constructing a conceptual model is to simplify the problem and to organize the available data so that the system can be analyzed accurately. Simplification is necessary, because a complete accounting of a system, such as Spring Creek, is not possible. The data and the conceptual model could be used in development of a fully coupled numerical model that dynamically links surface water, ground water, and land-use changes. The model could be used by decision makers to manage water resources within the basin and as a prototype that is transferable to other watersheds.

Pennsylvania↗
USGS2013

An integrated multi-criteria scenario evaluation web tool for participatory land-use planning in urbanized areas: The Ecosystem Portfolio Model

Land-use land-cover change is one of the most important and direct drivers of changes in ecosystem functions and services. Given the complexity of the decision-making, there is a need for Internet-based decision support systems with scenario evaluation capabilities to help planners, resource managers and communities visualize, compare and consider trade-offs among the many values at stake in land use planning. This article presents details on an Ecosystem Portfolio Model (EPM) prototype that integrates ecological, socio-economic information and associated values of relevance to decision-makers and stakeholders. The EPM uses a multi-criteria scenario evaluation framework, Geographic Information Systems (GIS) analysis and spatially-explicit land-use/land-cover change-sensitive models to characterize changes in important land-cover related ecosystem values related to ecosystem services and functions, land parcel prices, and community quality-of-life (QoL) metrics. Parameters in the underlying models can be modified through the interface, allowing users in a facilitated group setting to explore simultaneously issues of scientific uncertainty and divergence in the preferences of stakeholders. One application of the South Florida EPM prototype reported in this article shows the modeled changes (which are significant) in aggregate ecological value, landscape patterns and fragmentation, biodiversity potential and ecological restoration potential for current land uses compared to the 2050 land-use scenario. Ongoing refinements to EPM, and future work especially in regard to modifiable sea level rise scenarios are also discussed.

Environmental Modelling and Software↗
USGS1986

Hydrologic characteristics of soils in parts of Arkansas, Colorado, Kansas, Missouri, Nebraska, New Mexico, Oklahoma, South Dakota, and Texas

Certain physical characteristics of soils, including permeability, available water capacity, thickness, and topographic position, have a definite effect on the hydrology of an area. They control the rate at which precipitation infiltrates or is transmitted through the soil, and thus they have a significant role in determining the rates both of ground-water recharge and surface runoff. In studies of ground-water hydrology, it commonly is useful to differentiate soils areally according to their physical characteristics and to assign values indicating the hydrologic responses of soils in the different areas. In 1980 the U.S. Geologic Surey began a study termed the “Central Midwest Region Aquifer System Analysis,” Jorgensen and Signor, 1981). The regional aquifer system is defined as including all aquifers of Cretaceous age or older in all or parts of nine States included in this report. For part of this analysis, use of a soil-moisture model is planned to obtain estimates of actual evapotranspiraion and recharge to the regional aquifer system. Data for the hydrologic characteristics of the soils in the study area are required as input to this soil-moisture model.

Arkansas, Colorado, Kansas, Missouri, Nebraska, Ne↗
USGS2009

Volcanogenic Massive Sulfide Deposits of the World - Database and Grade and Tonnage Models

Grade and tonnage models are useful in quantitative mineral-resource assessments. The models and database presented in this report are an update of earlier publications about volcanogenic massive sulfide (VMS) deposits. These VMS deposits include what were formerly classified as kuroko, Cyprus, and Besshi deposits. The update was necessary because of new information about some deposits, changes in information in some deposits, such as grades, tonnages, or ages, revised locations of some deposits, and reclassification of subtypes. In this report we have added new VMS deposits and removed a few incorrectly classified deposits. This global compilation of VMS deposits contains 1,090 deposits; however, it was not our intent to include every known deposit in the world. The data was recently used for mineral-deposit density models (Mosier and others, 2007; Singer, 2008). In this paper, 867 deposits were used to construct revised grade and tonnage models. Our new models are based on a reclassification of deposits based on host lithologies: Felsic, Bimodal-Mafic, and Mafic volcanogenic massive sulfide deposits. Mineral-deposit models are important in exploration planning and quantitative resource assessments for two reasons: (1) grades and tonnages among deposit types vary significantly, and (2) deposits of different types occur in distinct geologic settings that can be identified from geologic maps. Mineral-deposit models combine the diverse geoscience information on geology, mineral occurrences, geophysics, and geochemistry used in resource assessments and mineral exploration. Globally based deposit models allow recognition of important features and demonstrate how common different features are. Well-designed deposit models allow geologists to deduce possible mineral-deposit types in a given geologic environment and economists to determine the possible economic viability of these resources. Thus, mineral-deposit models play a central role in presenting geoscience information in a useful form to policy makers. The foundation of mineral-deposit models is information about known deposits. The purpose of this publication is to present the latest geologic information and newly developed grade and tonnage models for VMS deposits in digital form. This publication contains computer files with information on VMS deposits from around the world. It also presents new grade and tonnage models for three subtypes of VMS deposits and a text file allowing locations of all deposits to be plotted in geographic information system (GIS) programs. The data are presented in FileMaker Pro and text files to make the information available to a wider audience. The value of this information and any derived analyses depends critically on the consistent manner of data gathering. For this reason, we first discuss the rules used in this compilation. Next, we provide new grade and tonnage models and analysis of the information in the file. Finally, the fields of the data file are explained. Appendix A gives the summary statistics for the new grade-tonnage models and Appendix B displays the country codes used in the database.

Open-File Report↗
USGS2010

Simulated effects of groundwater pumping and artificial recharge on surface-water resources and riparian vegetation in the Verde Valley sub-basin, Central Arizona

In the Verde Valley sub-basin, groundwater use has increased in recent decades. Residents and stakeholders in the area have established several groups to help in planning for sustainability of water and other resources of the area. One of the issues of concern is the effect of groundwater pumping in the sub-basin on surface water and on groundwater-dependent riparian vegetation. The Northern Arizona Regional Groundwater-Flow Model by Pool and others (in press) is the most comprehensive and up-to-date tool available to understand the effects of groundwater pumping in the sub-basin. Using a procedure by Leake and others (2008), this model was modified and used to calculate effects of groundwater pumping on surface-water flow and evapotranspiration for areas in the sub-basin. This report presents results for the upper two model layers for pumping durations of 10 and 50 years. Results are in the form of maps that indicate the fraction of the well pumping rate that can be accounted for as the combined effect of reduced surface-water flow and evapotranspiration. In general, the highest and most rapid responses to pumping were computed to occur near surface-water features simulated in the modified model, but results are not uniform along these features. The results are intended to indicate general patterns of model-computed response over large areas. For site-specific projects, improved results may require detailed studies of the local hydrologic conditions and a refinement of the modified model in the area of interest.

Scientific Investigations Report↗
USGS1993

Observations and modeling of seismic background noise

The preparation of this report had two purposes. One was to present a catalog of seismic background noise spectra obtained from a worldwide network of seismograph stations. The other purpose was to refine and document models of seismic background noise that have been in use for several years. The second objective was, in fact, the principal reason that this study was initiated and influenced the procedures used in collecting and processing the data. With a single exception, all of the data used in this study were extracted from the digital data archive at the U.S. Geological Survey's Albuquerque Seismological Laboratory (ASL). This archive dates from 1972 when ASL first began deploying digital seismograph systems and collecting and distributing digital data under the sponsorship of the Defense Advanced Research Projects Agency (DARPA). There have been many changes and additions to the global seismograph networks during the past twenty years, but perhaps none as significant as the current deployment of very broadband seismographs by the U.S. Geological Survey (USGS) and the University of California San Diego (UCSD) under the scientific direction of the IRIS consortium. The new data acquisition systems have extended the bandwidth and resolution of seismic recording, and they utilize high-density recording media that permit the continuous recording of broadband data. The data improvements and continuous recording greatly benefit and simplify surveys of seismic background noise. Although there are many other sources of digital data, the ASL archive data were used almost exclusively because of accessibility and because the data systems and their calibration are well documented for the most part. Fortunately, the ASL archive contains high-quality data from other stations in addition to those deployed by the USGS. Included are data from UCSD IRIS/IDA stations, the Regional Seismic Test Network (RSTN) deployed by Sandia National Laboratories (SNL), and the TERRAscope network deployed by the California Institute of Technology in cooperation with other institutions. A map showing the approximate locations of the stations used in this study is provided in Figure 1. One might hope for a better distribution of stations in the southern hemisphere, especially Africa and South America, in order to look for regional variations in seismic noise (apart from the major differences between continental, coastal and island sites). Unfortunately, anyone looking for subtle regional variations in seismic noise is probably going to be disappointed by the spectral data presented in this report because much of the station data appear to be dominated by local disturbances caused by instrumental, environmental, cultural, or surf noise. Better instruments and better instrument siting, or a well-funded field program, will be needed before a global isoseismal noise map can be produced. However, by assembling a composite of background noise from a large network of stations, many of the local station variables are masked, and it is possible to create generalized spectral plots of Earth noise for hypothetical quiet and noisy station sites.

Open-File Report↗
USGS2015

Web based visualization of large climate data sets

We have implemented the USGS National Climate Change Viewer (NCCV), which is an easy-to-use web application that displays future projections from global climate models over the United States at the state, county and watershed scales. We incorporate the NASA NEX-DCP30 statistically downscaled temperature and precipitation for 30 global climate models being used in the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC), and hydrologic variables we simulated using a simple water-balance model. Our application summarizes very large, complex data sets at scales relevant to resource managers and citizens and makes climate-change projection information accessible to users of varying skill levels. Tens of terabytes of high-resolution climate and water-balance data are distilled to compact binary format summary files that are used in the application. To alleviate slow response times under high loads, we developed a map caching technique that reduces the time it takes to generate maps by several orders of magnitude. The reduced access time scales to >500 concurrent users. We provide code examples that demonstrate key aspects of data processing, data exporting/importing and the caching technique used in the NCCV.

Environmental Modelling and Software↗
USGS1980

Research on interactive genetic-geological models to evaluate favourability for undiscovered uranium resources

Current methods of evaluating favourability for undiscovered uranium resources are unduly subjective, quite possibly inconsistent and, as a consequence, of questionable reliability. This research is aimed at reducing the subjectivity and increasing the reliability by designing an improved method that depends largely on geological data and their statistical frequency of occurrence. This progress report outlines a genetic approach to modelling the geological factors that controlled uranium mineralization in order to evaluate the favourability for the occurrence of undiscovered uranium deposits of the type modeled. A genetic model is constructed from all the factors that describe the processes, in chronological sequence, that formed uranium deposits thought to have a common origin. The field and laboratory evidence for the processes constitute a geologic occurrence base that parallels the chronological sequence of events. The genetic model and the geologic-occurrence base are portrayed as two columns of an interactive matrix called the "genetic-geologic model". For each column, eight chronological stages are used to describe the overall formation of the uranium deposits. These stages consist of (1) precursor processes; (2) host-rock formation; (3) preparation of host-rock; (4) uranium-source development; (5) transport of uranium; (6) primary uranium deposition; (7) post-deposition modification; and (8) preservation. To apply the genetic-geological model to evaluate favourability, a question is posed that determines the presence or absence of each attribute listed under the geologic-occurrence base. By building a logic circuit of the attributes according to either their essential or non-essential nature, the resultant match between a well-documented control area and the test area may be determined. The degree of match is a measure of favourability for uranium occurrence as hypothesized in the genetic model. This process of geological decision analysis results in a series of favourability maps that can be combined into a final composite favourability map.

Conference Paper↗
USGS1996

Progress report on daily flow-routing simulation for the Carson River, California and Nevada

A physically based flow-routing model using Hydrological Simulation Program-FORTRAN (HSPF) was constructed for modeling streamflow in the Carson River at daily time intervals as part of the Truckee-Carson Program of the U.S. Geological Survey (USGS). Daily streamflow data for water years 1978-92 for the mainstem river, tributaries, and irrigation ditches from the East Fork Carson River near Markleeville and West Fork Carson River at Woodfords down to the mainstem Carson River at Fort Churchill upstream from Lahontan Reservoir were obtained from several agencies and were compiled into a comprehensive data base. No previous physically based flow-routing model of the Carson River has incorporated multi-agency streamflow data into a single data base and simulated flow at a daily time interval. Where streamflow data were unavailable or incomplete, hydrologic techniques were used to estimate some flows. For modeling purposes, the Carson River was divided into six segments, which correspond to those used in the Alpine Decree that governs water rights along the river. Hydraulic characteristics were defined for 48 individual stream reaches based on cross-sectional survey data obtained from field surveys and previous studies. Simulation results from the model were compared with available observed and estimated streamflow data. Model testing demonstrated that hydraulic characteristics of the Carson River are adequately represented in the models for a range of flow regimes. Differences between simulated and observed streamflow result mostly from inadequate data characterizing inflow and outflow from the river. Because irrigation return flows are largely unknown, irrigation return flow percentages were used as a calibration parameter to minimize differences between observed and simulated streamflows. Observed and simulated streamflow were compared for daily periods for the full modeled length of the Carson River and for two major subreaches modeled with more detailed input data. Hydrographs and statistics presented in this report describe these differences. A sensitivity analysis of four estimated components of the hydrologic system evaluated which components were significant in the model. Estimated ungaged tributary streamflow is not a significant component of the model during low runoff, but is significant during high runoff. The sensitivity analysis indicates that changes in the estimated irrigation diversion and estimated return flow creates a noticeable change in the statistics. The modeling for this study is preliminary. Results of the model are constrained by current availability and accuracy of observed hydrologic data. Several inflows and outflows of the Carson River are not described by time-series data and therefore are not represented in the model.

Open-File Report↗
USGS2025

Random forest regression models for estimating low-streamflow statistics at ungaged locations in New York, excluding Long Island

This report is currently being updated with re‑calibrated random‑forest models for StreamStats implementation. The basin‑characteristic inputs were modified to align with StreamStats requirements. Abstract Models to estimate low-streamflow statistics at ungaged locations in New York, excluding Long Island and including hydrologically connected basins from bordering States, were developed for the first time by the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation. A total of 224 basin characteristics were developed for 213 unaltered streamgages (locations where the human effects on streamflow were limited), across the following categories: basin geometry, climate, land cover, soils, surficial geology, and other characteristics. The basins with unaltered streamgages were evaluated for potential redundancy, and streamgages in close proximity and with similar drainage areas were flagged and removed from the testing and cross-validation datasets to prevent data leaking from the training dataset to the testing dataset. Random forest regression models were created by using basin characteristics as predictor variables and by developing a workflow to train, tune, and test the model. Models were developed to estimate the ungaged lowest annual 7-day and 30-day average streamflow that occurs (on average) once every 10 years (7Q10 and 30Q10). The top four basin characteristics used for the 7Q10 and 30Q10 models were drainage area, total stream length, perimeter of the basin, and length of the longest flow path. Results for the 7Q10 and 30Q10 models had coefficients of determination ( R 2 ) of 0.796 and 0.853, respectively. The output model results were bias-corrected for ungaged locations across New York and are available within the interactive StreamStats tool.

New York↗
USGS2008

Effects of Land-Use Changes and Ground-Water Withdrawals on Stream Base Flow, Pocono Creek Watershed, Monroe County, Pennsylvania

The Pocono Creek watershed drains 46.5 square miles in eastern Monroe County, Pa. Between 2000 and 2020, the population of Monroe County is expected to increase by 70 percent, which will result in substantial changes in land-use patterns. An evaluation of the effect of reduced recharge from land-use changes and additional ground-water withdrawals on stream base flow was done by the U.S. Geological Survey (USGS) in cooperation with the U.S. Environmental Protection Agency (USEPA) and the Delaware River Basin Commission as part of the USEPA?s Framework for Sustainable Watershed Management Initiative. Two models were used. A Soil and Water Assessment Tool (SWAT) model developed by the USEPA provided areal recharge values for 2000 land use and projected full buildout land use. The USGS MODFLOW-2000 ground-water-flow model was used to estimate the effect of reduced recharge from changes in land use and additional ground-water withdrawals on stream base flow. This report describes the ground-water-flow-model simulations. The Pocono Creek watershed is underlain by sedimentary rock of Devonian age, which is overlain by a veneer of glacial deposits. All water-supply wells are cased into and derive water from the bedrock. In the ground-water-flow model, the surficial geologic units were grouped into six categories: (1) moraine deposits, (2) stratified drift, (3) lake deposits, (4) outwash, (5) swamp deposits, and (6) undifferentiated deposits. The unconsolidated surficial deposits are not used as a source of water. The ground-water and surface-water systems are well connected in the Pocono Creek watershed. Base flow measured on October 13, 2004, at 27 sites for model calibration showed that streams gained water between all sites measured except in the lower reach of Pocono Creek. The ground-water-flow model included the entire Pocono Creek watershed. Horizontally, the modeled area was divided into a 53 by 155 cell grid with 6,060 active cells. Vertically, the modeled area was discretized into four layers. Layers 1 and 2 represented the unconsolidated surficial deposits where they are present and bedrock where the surficial deposits are absent. Layer 3 represented shallow bedrock and was 200 ft (feet) thick. Layer 4 represented deep bedrock and was 300 ft thick. A total of 873 cells representing streams were assigned to layer 1. Recharge rates for model calibration were provided by the USEPA SWAT model for 2000 land-use conditions. Recharge rates for 2000 for the 29 subwatersheds in the SWAT model ranged from 6.11 to 22.66 inches per year. Because the ground-water-flow model was calibrated to base-flow data collected on October 13, 2004, the 2000 recharge rates were multiplied by 1.18 so the volume of recharge was equal to the volume of streamflow measured at the mouth of Pocono Creek. During model calibration, adjustments were made to aquifer hydraulic conductivity and streambed conductance. Simulated base flows and hydraulic heads were compared to measured base flows and hydraulic heads using the root mean squared error (RMSE) between measured and simulated values. The RMSE of the calibrated model for base flow was 4.7 cubic feet per second for 27 locations, and the RMSE for hydraulic heads for 15 locations was 35 ft. The USEPA SWAT model was used to provide areal recharge values for 2000 and full buildout land-use conditions. The change in recharge ranged from an increase of 37.8 percent to a decrease of 60.8 percent. The ground-water-flow model was used to simulate base flow for 2000 and full buildout land-use conditions using steady-state simulations. The decrease in simulated base flow ranged from 3.8 to 63 percent at the streamflow-measurement sites. Simulated base flow at streamflow-gaging station Pocono Creek above Wigwam Run near Stroudsburg, Pa. (01441495), decreased 25 percent. This is in general agreement with the SWAT model, which estimated a 30.6-percent loss in base flow at the streamflow-gaging station.

Scientific Investigations Report↗
USGS2004

A theory for modeling ground-water flow in heterogeneous media

Construction of a ground-water model for a field area is not a straightforward process. Data are virtually never complete or detailed enough to allow substitution into the model equations and direct computation of the results of interest. Formal model calibration through optimization, statistical, and geostatistical methods is being applied to an increasing extent to deal with this problem and provide for quantitative evaluation and uncertainty analysis of the model. However, these approaches are hampered by two pervasive problems: 1) nonlinearity of the solution of the model equations with respect to some of the model (or hydrogeologic) input variables (termed in this report system characteristics) and 2) detailed and generally unknown spatial variability (heterogeneity) of some of the system characteristics such as log hydraulic conductivity, specific storage, recharge and discharge, and boundary conditions. A theory is developed in this report to address these problems. The theory allows construction and analysis of a ground-water model of flow (and, by extension, transport) in heterogeneous media using a small number of lumped or smoothed system characteristics (termed parameters). The theory fully addresses both nonlinearity and heterogeneity in such a way that the parameters are not assumed to be effective values. The ground-water flow system is assumed to be adequately characterized by a set of spatially and temporally distributed discrete values, ?, of the system characteristics. This set contains both small-scale variability that cannot be described in a model and large-scale variability that can. The spatial and temporal variability in ? are accounted for by imagining ? to be generated by a stochastic process wherein ? is normally distributed, although normality is not essential. Because ? has too large a dimension to be estimated using the data normally available, for modeling purposes ? is replaced by a smoothed or lumped approximation y?. (where y is a spatial and temporal interpolation matrix). Set y?. has the same form as the expected value of ?, y 'line' ? , where 'line' ? is the set of drift parameters of the stochastic process; ?. is a best-fit vector to ?. A model function f(?), such as a computed hydraulic head or flux, is assumed to accurately represent an actual field quantity, but the same function written using y?., f(y?.), contains error from lumping or smoothing of ? using y?.. Thus, the replacement of ? by y?. yields nonzero mean model errors of the form E(f(?)-f(y?.)) throughout the model and covariances between model errors at points throughout the model. These nonzero means and covariances are evaluated through third and fifth-order accuracy, respectively, using Taylor series expansions. They can have a significant effect on construction and interpretation of a model that is calibrated by estimating ?.. Vector ?.. is estimated as 'hat' ? using weighted nonlinear least squares techniques to fit a set of model functions f(y'hat' ?) to a. corresponding set of observations of f(?), Y. These observations are assumed to be corrupted by zero-mean, normally distributed observation errors, although, as for ?, normality is not essential. An analytical approximation of the nonlinear least squares solution is obtained using Taylor series expansions and perturbation techniques that assume model and observation errors to be small. This solution is used to evaluate biases and other results to second-order accuracy in the errors. The correct weight matrix to use in the analysis is shown to be the inverse of the second-moment matrix E(Y-f(y?.))(Y-f(y?.))', but the weight matrix is assumed to be arbitrary in most developments. The best diagonal approximation is the inverse of the matrix of diagonal elements of E(Y-f(y?.))(Y-f(y?.))', and a method of estimating this diagonal matrix when it is unknown is developed using a special objective function to compute 'hat' ?. When considered to be an estimate of f

Professional Paper↗
USGS2010

Estimation of groundwater use for a groundwater-flow model of the Lake Michigan Basin and adjacent areas, 1864-2005

The U.S. Geological Survey, at the request of Congress, is assessing the availability and use of the Nation's water resources to help characterize how much water is available now, how water availability is changing, and how much water can be expected to be available in the future. The Great Lakes Basin Pilot project of the U.S. Geological Survey national assessment of water availability and use focused on the Great Lakes Basin and included detailed studies of the processes governing water availability in the Great Lakes Basin. One of these studies included the development of a groundwater-flow model of the Lake Michigan Basin. This report describes the compilation and estimation of the groundwater withdrawals in those areas in Wisconsin, Michigan, Indiana, and Illinois that were needed for the Lake Michigan Basin study groundwater-flow model. These data were aggregated for 12 model time intervals spanning 1864 to 2005 and were summarized by model area, model subregion, category of water use, aquifer system, aquifer type, and hydrogeologic unit model layer. The types and availability of information on groundwater withdrawals vary considerably among states because water-use programs often differ in the types of data collected and in the methods and frequency of data collection. As a consequence, the methods used to estimate and verify the data also vary. Additionally, because of the different sources of data and different terminologies applied for the purposes of this report, the water-use data published in this report may differ from water-use data presented in other reports. These data represent only a partial estimate of groundwater use in each state because estimates were compiled only for areas in Wisconsin, Michigan, Indiana, and Illinois within the Lake Michigan Basin model area. Groundwater-withdrawal data were compiled for both nearfield and farfield model areas in Wisconsin and Illinois, whereas these data were compiled primarily for the nearfield model area in Michigan and Indiana. Overall water use for the selected areas in Wisconsin, Michigan, Indiana, and Illinois was less during early time intervals than during more recent intervals, with large increases beginning around the 1960s. Total estimated groundwater withdrawals for model input range from 18.01 million gallons per day (Mgal/d) for interval 1 (1864-1900) to 1,280.25 Mgal/d for interval 12 (2001-5). Withdrawals for the public-supply category make up the majority of the withdrawals in each of the four states. In Wisconsin and Michigan, the second largest withdrawals are for the irrigation category; in Indiana and Illinois, industrial withdrawals account for the second largest withdrawal amounts. The smallest withdrawals are for miscellaneous uses in Wisconsin and irrigation uses in Indiana and Illinois. Estimated groundwater withdrawals in the Southern Lower Peninsula of Michigan, Northeastern Illinois, and the farfield model area are generally larger than in the other model subregions. Withdrawals in Michigan and Indiana are predominantly from the Quaternary aquifer system, whereas withdrawals in Illinois are predominantly from the Cambrian-Ordovician aquifer systems. Withdrawals in Wisconsin are about equal from the Quaternary and Cambrian-Ordovician aquifer systems. Estimated groundwater withdrawals in Michigan and Indiana are predominantly from the unconfined unconsolidated aquifer type. Withdrawals in Illinois are largely from the deep confined bedrock aquifer type, although they decreased considerably in more recent time intervals. Wisconsin withdrawals are about equal from unconfined unconsolidated and deep confined bedrock aquifer types. Groundwater-withdrawal estimates in Wisconsin were compiled for the 47 easternmost counties within the boundary of the Lake Michigan Basin model, of which 32 counties, though not entirely contained, are at least partly within the Lake Michigan Basin. Overall, 6,457 withdrawal locations were estima

Scientific Investigations Report↗
USGS2015

A model for evaluating stream temperature response to climate change in Wisconsin

Expected climatic changes in air temperature and precipitation patterns across the State of Wisconsin may alter future stream temperature and flow regimes. As a consequence of flow and temperature changes, the composition and distribution of fish species assemblages are expected to change. In an effort to gain a better understanding of how climatic changes may affect stream temperature, an approach was developed to predict and project daily summertime stream temperature under current and future climate conditions for 94,341 stream kilometers across Wisconsin. The approach uses a combination of static landscape characteristics and dynamic time-series climatic variables as input for an Artificial Neural Network (ANN) Model integrated with a Soil-Water-Balance (SWB) Model. Future climate scenarios are based on output from downscaled General Circulation Models (GCMs). The SWB model provided a means to estimate the temporal variability in groundwater recharge and provided a mechanism to evaluate the effect of changing air temperature and precipitation on groundwater recharge and soil moisture. The Integrated Soil-Water-Balance and Artificial Neural Network version 1 (SWB-ANNv1) Model was used to simulate daily summertime stream temperature under current (1990–2008) climate and explained 76 percent of the variation in the daily mean based on validation at 67 independent sites. Results were summarized as July mean water temperature, and individual stream segments were classified by thermal class (cold, cold transition, warm transition, and warm) for comparison of current (1990–2008) with future climate conditions. Integrating the SWB Model with the ANN Model provided a mechanism by which downscaled global or regional climate model results could be used to estimate the potential effects of climate change on future stream temperature on a daily time step. To address future climate scenarios, statistically downscaled air temperature and precipitation projections from 10 GCMs and 2 time periods were used with the SWB-ANNv1 Model to project future stream temperature. Projections of future stream temperatures at mid- (2046–65) and late- (2081–2100) 21st century showed the July mean water temperature increasing for all stream segments with about 80 percent of stream kilometers increasing by 1 to 2 degrees Celsius (°C) by mid-century and about 99 percent increasing by 1 to 3 °C by late-century. Projected changes in stream temperatures also affected changes in thermal classes with a loss in the total amount of cold-water, cold-transition, and warm-transition thermal habitat and a gain in warm-water and very warm thermal habitat for both mid- and late-21st century time periods. The greatest losses occurred for cold-water streams and the greatest gains for warm-water streams, with a contraction of cold-water streams in the Driftless Area of western and southern Wisconsin and an expansion of warm-water streams across northern Wisconsin. Results of this study suggest that such changes will affect the composition of fish assemblages, with a loss of suitable habitat for cold-water fishes and gain in suitable habitat for warm-water fishes. In the end, these projected changes in thermal habitat attributable to climate may result in a net loss of fisheries, because many warm-water species may be unable to colonize habitats formerly occupied by cold-water species because of other habitat limitations (e.g., stream size, gradient). Although projected stream temperatures may vary greatly, depending on the emissions scenario and models used, the results presented in this report represent one possibility. The relative change in stream temperature can provide useful information for planning for potential climate impacts to aquatic ecosystems. Model results can be used to help identify vulnerabilities of streams to climate change, guide stream surveys and thermal classifications, prioritize the allocation of scarce financial resources, identify approaches to climate adaptation to best protect and enhance resiliency in stream thermal habitat, and provide information to make quantitative assessments of statewide stream resources.

Wisconsin↗
USGS2008

Audiomagnetotelluric data and preliminary two-dimensional models from Spring, Dry Lake, and Delamar Valleys, Nevada

This report presents audiomagnetotelluric (AMT) data along fourteen profiles in Spring, Delamar, and Dry Lake Valleys, and the corresponding preliminary two-dimensional (2-D) inverse models. The AMT method is a valuable tool for estimating the electrical resistivity of the Earth over depth ranges from a few meters to less than one kilometer, and it is important for revealing subsurface structure and stratigraphy within the Basin and Range province of eastern Nevada, which can be used to define the geohydrologic framework of the region. We collected AMT data by using the Geometrics StrataGem EH4 system. Profiles were 0.7 - 3.2 km in length with station spacing of 50-400 m. Data were recorded in a coordinate system parallel to and perpendicular to the regional geologic-strike direction with Z positive down. We show AMT station locations, sounding curves of apparent resistivity, phase, and coherency, and 2-D models of subsurface resistivity along the profiles. The 2-D inverse models are computed from the transverse electric (TE), transverse magnetic (TM), and TE+TM mode data by using a conjugate gradient, finite-difference method. Preliminary interpretation of the 2-D models defines the structural framework of the basins and the resistivity contrasts between alluvial basin-fill, volcanic units, and carbonate basement rocks.

Nevada↗