USGS Science⌕ Search

SEARCH · USGS Science

Results for “Input”

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

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

At least 1,243 records · Page 69Linked to original sources

Three-dimensional hydrogeologic framework model of the Russian River watershed, California

The Russian River watershed is in northern Sonoma County and southern Mendocino County, California, in the northern part of the California Coast Ranges. The Russian River serves as a supply for agricultural irrigation and for municipal, domestic, and commercial uses. Through a cooperative agreement with the California State Water Resources Control Board and Sonoma County Water Agency, the U.S. Geological Survey has completed studies to better understand the hydrogeologic system and develop numerical hydrologic modeling tools to evaluate and aid in managing groundwater resources. This report focuses on the development of a digital three-dimensional hydrogeologic framework model of the Russian River watershed for use in groundwater resource assessment and numerical models. The digital three-dimensional hydrogeologic framework model of the Russian River watershed portrays the altitude, thickness, and extent of five hydrogeologic units. These five hydrogeologic units include (1) a basement unit, (2) the Sonoma Volcanics, (3) a consolidated sedimentary rock unit, (4) an unconsolidated sediment unit, and (5) channel alluvium. Model input data were compiled from published geologic maps, interpreted well data, and a model of the top of basement derived from gravity data. These data were used to construct surfaces that represent the upper and lower subsurface boundaries of each hydrogeologic unit. Top surfaces were created for the five hydrogeologic units and then stacked in three dimensions to create a solid-volume digital model. The digital three-dimensional hydrogeologic framework model described in this report and the corresponding data represent the generalized geometry of the subsurface geologic units; the model reproduces the input geologic data with reasonable accuracy and is consistent with previously published subsurface conceptualizations of the region. The model indicates the overall geometry of the basement within the watershed and the spatial extent, altitude, and thickness of the basin-filling units. The hydrogeologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

California↗

Evaluating drought risk of the Red River of the North Basin using historical and stochastic streamflow upstream from Emerson, Manitoba

Drought and its effect on streamflow are important to understand because of the potential to adversely affect water supply, agricultural production, and ecological conditions. The Red River of the North Basin in north-central United States and central Canada is susceptible to dry conditions. During an extended drought, streamflow conditions in the Red River of the North may become inadequate to support existing water supply needs in the basin for agriculture, industry, human use, and aquatic life. To understand potential future low-streamflow conditions in the Red River of the North Basin, the U.S. Geological Survey, in cooperation with the International Joint Commission, North Dakota Department of Water Resources, Red River Joint Water Resource District, and Red River Watershed Management Board, developed a water-balance model of the Red River of the North Basin upstream from Emerson, Manitoba, Canada, and coupled the model with stochastic weather inputs to simulate possible future low-streamflow conditions. Historical changes in low-streamflow conditions were characterized across the Red River of the North Basin using multiple change-point analysis for 12 streamgages. Across these stations, significant change-point years in 1943 and 1994 marked increases in the magnitude of low-streamflow conditions. During 1920–2015, conversion of primary land (not affected by human use) to agricultural and secondary land was followed by a conversion from smalls grains to corn and soybeans as the dominant crop type. From land-use analysis, 1940–2000 was determined to have relatively stable land use and therefore was used as the calibration period for the water-balance model. A deterministic water-balance model was developed for the Red River of the North Basin upstream from Emerson, Manitoba. The water-balance model was calibrated with data from 37 U.S. Geological Survey streamgages for 1940–2000 and verified using data for 2001–15. The calibrated water-balance model simulated streamflow distributions that mirrored the seasonal patterns of the observed mean monthly streamflow and the standard deviation of the monthly streamflow data, especially during the fall and winter months when streamflow was lowest. For the verification period, during the low-streamflow months of December through January, the difference between simulated and observed data was similar to the calibration comparison and successfully reproduced seasonal trends in the distribution of streamflow, even when using weather data that were outside the calibration period. To determine the future risk of low-streamflow conditions in the Red River of the North Basin, a block-bootstrap method was used to generate multiple possible future climates. These stochastically generated weather time series were then input to a water-balance model to simulate a distribution of possible streamflows. Three sets of experiments were performed, with each experiment containing a set of scenarios. The first set of experiments from the stochastic streamflow model were designed to investigate how changes in reservoir management would affect the distribution of low streamflow. Relative to scenario 1 (present-day [2023] reservoir operation), scenario 2 (no reservoir operation) shifted the low-streamflow frequency curves downward, reducing the annual minimum monthly streamflow for the Emerson subbasin. Subbasins were defined by the contributing area upstream from a selected streamgage station. Relative to scenario 1, scenario 3 (regulated streamflow with an increased reservoir capacity of 10 percent) shifted the low-streamflow frequency curves upward for the Emerson subbasin. The magnitude of this upward shift, caused by increased reservoir capacity, was lower than the magnitude of the shift caused by the absence of the reservoirs, which indicates that the streamflow was most affected when the reservoirs were first constructed. The second set of experiments from the stochastic streamflow model included two scenarios that were performed to better understand how the Red River of the North Basin responds to long periods of low or high precipitation. The results indicate that the model consistently overestimated streamflow, but the relative change between a wet and dry climate state of simulated streamflow distribution reasonably matched the relative change of historical streamflow. Across the subbasins, the model was most accurate for low-streamflow conditions associated with nonexceedance probabilities between 20 and 40 percent. The third set of experiments from the stochastic streamflow model were done to investigate low-streamflow response across the basin to several drought events. Low-end streamflow was reduced when the basin was exposed to a drought, and the magnitude of the reduction increased with longer or more intense droughts. Compared to the low-intensity drought scenarios, the range of percent reductions (as indicated by the interquartile range) was larger for the high-intensity drought scenarios for all subbasins, and the subbasins of Grand Forks and Emerson had a smaller range of reductions compared to the other three subbasins. The larger drainage area—combined with the large contribution of the Red Lake River and several other Minnesota tributaries that generally experience wetter climate conditions—upstream from the Emerson and Grand Forks subbasins may contribute to the smaller range in reductions under the high intensity scenarios. Comparison of the percent reduction in low-end streamflow among subbasins also indicated that the effects of drought duration and intensity could be cumulative. Combining factors of time and intensity produced a larger reduction in streamflow than when each effect was isolated. The array of drought scenarios can be used to determine how a subbasin would respond to multiple possible future conditions. Based on climate predictions, the drought scenario that best matches a future anticipated drought scenario can be used to estimate a low streamflow response for a given subbasin.

Manitoba, Minnesota, North Dakota↗

MODFLOW-2005, the U.S. Geological Survey modular ground-water model - documentation of shared node local grid refinement (LGR) and the boundary flow and head (BFH) package

This report documents the addition of shared node Local Grid Refinement (LGR) to MODFLOW-2005, the U.S. Geological Survey modular, transient, three-dimensional, finite-difference ground-water flow model. LGR provides the capability to simulate ground-water flow using one block-shaped higher-resolution local grid (a child model) within a coarser-grid parent model. LGR accomplishes this by iteratively coupling two separate MODFLOW-2005 models such that heads and fluxes are balanced across the shared interfacing boundary. LGR can be used in two-and three-dimensional, steady-state and transient simulations and for simulations of confined and unconfined ground-water systems. Traditional one-way coupled telescopic mesh refinement (TMR) methods can have large, often undetected, inconsistencies in heads and fluxes across the interface between two model grids. The iteratively coupled shared-node method of LGR provides a more rigorous coupling in which the solution accuracy is controlled by convergence criteria defined by the user. In realistic problems, this can result in substantially more accurate solutions and require an increase in computer processing time. The rigorous coupling enables sensitivity analysis, parameter estimation, and uncertainty analysis that reflects conditions in both model grids. This report describes the method used by LGR, evaluates LGR accuracy and performance for two- and three-dimensional test cases, provides input instructions, and lists selected input and output files for an example problem. It also presents the Boundary Flow and Head (BFH) Package, which allows the child and parent models to be simulated independently using the boundary conditions obtained through the iterative process of LGR.

Techniques and Methods↗

Revised multi-node well (MNW2) package for MODFLOW ground-water flow model

Wells that are open to multiple aquifers can provide preferential pathways to flow and solute transport that short-circuit normal fluid flowlines. Representing these features in a regional flow model can produce a more realistic and reliable simulation model. This report describes modifications to the Multi-Node Well (MNW) Package of the U.S. Geological Survey (USGS) three-dimensional ground-water flow model (MODFLOW). The modifications build on a previous version and add several new features, processes, and input and output options. The input structure of the revised MNW (MNW2) is more well-centered than the original verion of MNW (MNW1) and allows the user to easily define hydraulic characteristics of each multi-node well. MNW2 also allows calculations of additional head changes due to partial penetration effects, flow into a borehole through a seepage face, changes in well discharge related to changes in lift for a given pump, and intraborehole flows with a pump intake located at any specified depth within the well. MNW2 also offers an improved capability to simulate nonvertical wells. A new output option allows selected multi-node wells to be designated as 'observation wells' for which changes in selected variables with time will be written to separate output files to facilitate postprocessing. MNW2 is compatible with the MODFLOW-2000 and MODFLOW-2005 versions of MODFLOW and with the version of MODFLOW that includes the Ground-Water Transport process (MODFLOW-GWT).

Techniques and Methods↗

MODFLOW–LGR—Documentation of ghost node local grid refinement (LGR2) for multiple areas and the boundary flow and head (BFH2) package

This report documents the addition of ghost node Local Grid Refinement (LGR2) to MODFLOW-2005, the U.S. Geological Survey modular, transient, three-dimensional, finite-difference groundwater flow model. LGR2 provides the capability to simulate groundwater flow using multiple block-shaped higher-resolution local grids (a child model) within a coarser-grid parent model. LGR2 accomplishes this by iteratively coupling separate MODFLOW-2005 models such that heads and fluxes are balanced across the grid-refinement interface boundary. LGR2 can be used in two-and three-dimensional, steady-state and transient simulations and for simulations of confined and unconfined groundwater systems. Traditional one-way coupled telescopic mesh refinement methods can have large, often undetected, inconsistencies in heads and fluxes across the interface between two model grids. The iteratively coupled ghost-node method of LGR2 provides a more rigorous coupling in which the solution accuracy is controlled by convergence criteria defined by the user. In realistic problems, this can result in substantially more accurate solutions and require an increase in computer processing time. The rigorous coupling enables sensitivity analysis, parameter estimation, and uncertainty analysis that reflects conditions in both model grids. This report describes the method used by LGR2, evaluates accuracy and performance for two-and three-dimensional test cases, provides input instructions, and lists selected input and output files for an example problem. It also presents the Boundary Flow and Head (BFH2) Package, which allows the child and parent models to be simulated independently using the boundary conditions obtained through the iterative process of LGR2.

Techniques and Methods↗

Two graphical user interfaces for managing and analyzing MODFLOW groundwater-model scenarios

Scenario Manager and Scenario Analyzer are graphical user interfaces that facilitate the use of calibrated, MODFLOW-based groundwater models for investigating possible responses to proposed stresses on a groundwater system. Scenario Manager allows a user, starting with a calibrated model, to design and run model scenarios by adding or modifying stresses simulated by the model. Scenario Analyzer facilitates the process of extracting data from model output and preparing such display elements as maps, charts, and tables. Both programs are designed for users who are familiar with the science on which groundwater modeling is based but who may not have a groundwater modeler’s expertise in building and calibrating a groundwater model from start to finish. With Scenario Manager, the user can manipulate model input to simulate withdrawal or injection wells, time-variant specified hydraulic heads, recharge, and such surface-water features as rivers and canals. Input for stresses to be simulated comes from user-provided geographic information system files and time-series data files. A Scenario Manager project can contain multiple scenarios and is self-documenting. Scenario Analyzer can be used to analyze output from any MODFLOW-based model; it is not limited to use with scenarios generated by Scenario Manager. Model-simulated values of hydraulic head, drawdown, solute concentration, and cell-by-cell flow rates can be presented in display elements. Map data can be represented as lines of equal value (contours) or as a gradated color fill. Charts and tables display time-series data obtained from output generated by a transient-state model run or from user-provided text files of time-series data. A display element can be based entirely on output of a single model run, or, to facilitate comparison of results of multiple scenarios, an element can be based on output from multiple model runs. Scenario Analyzer can export display elements and supporting metadata as a Portable Document Format file.

Techniques and Methods↗

MT3D-USGS version 1: A U.S. Geological Survey release of MT3DMS updated with new and expanded transport capabilities for use with MODFLOW

MT3D-USGS, a U.S. Geological Survey updated release of the groundwater solute transport code MT3DMS, includes new transport modeling capabilities to accommodate flow terms calculated by MODFLOW packages that were previously unsupported by MT3DMS and to provide greater flexibility in the simulation of solute transport and reactive solute transport. Unsaturated-zone transport and transport within streams and lakes, including solute exchange with connected groundwater, are among the new capabilities included in the MT3D-USGS code. MT3D-USGS also includes the capability to route a solute through dry cells that may occur in the Newton-Raphson formulation of MODFLOW (that is, MODFLOW-NWT). New chemical reaction Package options include the ability to simulate inter-species reactions and parent-daughter chain reactions. A new pump-and-treat recirculation package enables the simulation of dynamic recirculation with or without treatment for combinations of wells that are represented in the flow model, mimicking the above-ground treatment of extracted water. A reformulation of the treatment of transient mass storage improves conservation of mass and yields solutions for better agreement with analytical benchmarks. Several additional features of MT3D-USGS are (1) the separate specification of the partitioning coefficient ( K d ) within mobile and immobile domains; (2) the capability to assign prescribed concentrations to the top-most active layer; (3) the change in mass storage owing to the change in water volume now appears as its own budget item in the global mass balance summary; (4) the ability to ignore cross-dispersion terms; (5) the definition of Hydrocarbon Spill-Source Package (HSS) mass loading zones using regular and irregular polygons, in addition to the currently supported circular zones; and (6) the ability to specify an absolute minimum thickness rather than the default percent minimum thickness in dry-cell circumstances. Benchmark problems that implement the new features and packages test the accuracy of new code through comparison to analytical benchmarks, as well as to solutions from other published codes. The input file structure for MT3D-USGS adheres to MT3DMS conventions for backward compatibility: the new capabilities and packages described herein are readily invoked by adding three-letter package name acronyms to the name file or by setting input flags as needed. Memory is managed in MT3D-USGS using FORTRAN modules in order to simplify code development and expansion.

Techniques and Methods↗

One-Water Hydrologic Flow Model: A MODFLOW based conjunctive-use simulation software

The U.S. Geological Survey’s (USGS) Modular Ground-Water Flow Model (MODFLOW-2005) is a computer program that simulates groundwater flow by using finite differences. The MODFLOW-2005 framework uses a modular design that allows for the easy development and incorporation of new features called processes and packages that work with or modify inputs to the groundwater-flow equation. A process solves a flow equation or set of equations. For example, the central part of MODFLOW is the groundwater-flow process that solves the groundwater-flow equation; the surface-water routing process is an additional process that solves the surface-water flow equation. Packages are code related to the groundwater-flow process. For example, the subsidence package modifies the groundwater-flow process by including aquifer compaction effects on flow. With the development of new packages and processes, the MODFLOW-2005 base framework diverged into multiple independent versions designed for specific simulation needs. This divergence limited each independent MODFLOW release to its specific purpose, so that there was no longer a single, comprehensive, general-purpose hydraulic-simulation framework. The MODFLOW One-Water Hydrologic Flow Model (MF-OWHM, also informally known as OneWater) is an integrated hydrologic flow model that combines multiple MODFLOW-2005 variants in one cohesive simulation software; changes were made to enable multiple capabilities in one code. This fusion of the MODFLOW-2005 versions resulted in a simulation software that can be used to address and analyze a wide class of conjunctive-use, water-management, water-food-security, and climate-crop-water scenarios. As a second core version of MODFLOW-2005, MF-OWHM maintains backward compatibility with existing MODFLOW-2005 versions, with features that include the following: Process-based simulation. Saturated groundwater flow (three-dimensional). Surface-water flow (one- and two-dimensional). Stream and river flow. Lake and reservoir storage. Landscape simulation and irrigated agriculture. Land-use and crop simulation. Root uptake of groundwater. Actual evapotranspiration. Estimated irrigation demand. Reservoir operations. Aquifer compaction and subsidence by vertical model-grid deformation. Seawater intrusion by a sharp-interface assumption. Karst-aquifer and fractured-bedrock flow. Turbulent and laminar-pipe network flow. Unsaturated groundwater flow (one-dimensional). Internal linkages among the processes that couple hydraulic head, flow, and deformation. Redesigned code for faster simulation, increased user-input options, easier model updates, and more robust error reporting than in previous models. MF-OWHM is a MODFLOW-2005 based integrated hydrologic model that can simulate and analyze varying environmental conditions to allow for the evaluation of management options from many components of human and natural water movement through a physically based, supply and demand framework. The term “integrated,” in the context of this report, refers to the tight coupling of groundwater flow, surface-water flow, landscape processes, aquifer compaction and subsidence, reservoir operations, and conduit (karst) flow. Another benefit of this integrated hydrologic model is that models developed to run by MODFLOW-2005, MODFLOW-NWT, MODFLOW-CFP, or MODFLOW-FMP can also be simulated with MF-OWHM. At the time of this report’s publication, MF-OWHM version 2 (MF-OWHM2) does not include a direct internal simulation of snowmelt, advanced mountainous watershed rainfall-runoff simulation, detailed shallow soil-moisture accounting, or atmospheric moisture content. Atmospheric moisture may be accounted for indirectly by, optionally, specifying a pan-evaporation rate, reference evapotranspiration, and precipitation. These features are not included to ensure that simulation runtime remains short enough to enable the use of automated methods of calibrating model parameters to field observations, which typically require many simulation model runs. The MF-OWHM approach is to include as much detail as possible to simulate hydrological processes, providing the simulation runtimes remain reasonable enough to allow for robust parameter estimation and model calibration. To represent both natural and human-influenced flow, MF-OWHM integrates physically based flow processes derived from MODFLOW-2005 in a supply and demand framework. From this integration, the physically based movement of groundwater, surface water, imported water, and precipitation serve as supply to meet consumptive demands associated with irrigated and non-irrigated agriculture, natural vegetation, and urban water uses. Water consumption is determined by balancing the available water supply with water demand, leading to the concept of a demand-driven, supply-constrained simulation. The MF-OWHM Supply-and-Demand Framework is especially useful for the analysis of agricultural water use, where there are often few data available to describe changes in land-use through time, such as crop type and distribution, and the associated changes in groundwater pumpage. This framework attempts to satisfy each land-use water demand with available water supplies—that is, groundwater uptake, precipitation, and irrigation. An option provided in MF-OWHM2 is to automatically increase groundwater pumping for irrigation, which often is unknown, by the calculated residual between demand and the other available sources of supply. From large- to small-scale applications, the physically based supply and demand framework provides key capabilities for simulating and analyzing historical, current, and future conjunctive-use of surface water and groundwater. To achieve the physically based supply and demand framework, the MODFLOW-2005 standard of no inter-package and -process communication was relaxed for MF-OWHM2. Traditional MODFLOW simulation models required that all packages and processes interact through the groundwater-flow equation or by removing the water flow from the simulation domain. For example, the MODFLOW-2005 representation of a groundwater well extracts water from the groundwater-flow equation (by subtraction) and removes it from the simulation domain. This feature is available in the MF-OWHM framework, but options have been added to allow the specification of a use or destination of pumped groundwater within the model domain, for example, it can be used for irrigation, managed aquifer recharge, or return-flow to streams.

Techniques and Methods↗

PHAST--a program for simulating ground-water flow, solute transport, and multicomponent geochemical reactions

The computer program PHAST simulates multi-component, reactive solute transport in three-dimensional saturated ground-water flow systems. PHAST is a versatile ground-water flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. PHAST is applicable to the study of natural and contaminated ground-water systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock-water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, density-dependent flow, or waters with high ionic strengths. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux, and leaky conditions, as well as the special cases of rivers and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, gases, surface complexation sites, ion exchange sites, and solid solutions; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, gases, exchangers, surfaces, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and post-processing programs; or in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST. Operator splitting of the flow, transport, and geochemical equations is used to separate the three processes into three sequential calculations. No iterations between transport and reaction calculations are implemented. A three-dimensional Cartesian coordinate system and finite-difference techniques are used for the spatial and temporal discretization of the flow and transport equations. The non-linear chemical equilibrium equations are solved by a Newton-Raphson method, and the kinetic reaction equations are solved by a Runge-Kutta or an implicit method for integrating ordinary differential equations. The PHAST simulator may require large amounts of memory and long Central Processing Unit (CPU) times. To reduce the long CPU times, a parallel version of PHAST has been developed that runs on a multiprocessor computer or on a collection of computers that are networked. The parallel version requires Message Passing Interface, which is currently (2004) freely available. The parallel version is effective in reducing simulation times. This report documents the use of the PHAST simulator, including running the simulator, preparing the input files, selecting the output files, and visualizing the results. It also presents four examples that verify the numerical method and demonstrate the capabilities of the simulator. PHAST requires three input files. Only the flow and transport file is described in detail in this report. The other two files, the chemistry data file and the database file, are identical to PHREEQC files and the detailed description of these files is found in the PHREEQC documentation.

Techniques and Methods↗

Section 3. The SPARROW Surface Water-Quality Model—Theory, application and user documentation

SPARROW (SPAtially Referenced Regressions On Watershed attributes) is a watershed modeling technique for relating water-quality measurements made at a network of monitoring stations to attributes of the watersheds containing the stations. The core of the model consists of a nonlinear regression equation describing the non-conservative transport of contaminants from point and diffuse sources on land to rivers and through the stream and river network. The model predicts contaminant flux, concentration, and yield in streams and has been used to evaluate alternative hypotheses about the important contaminant sources and watershed properties that control transport over large spatial scales. This report provides documentation for the SPARROW modeling technique and computer software to guide users in constructing and applying basic SPARROW models. The documentation gives details of the SPARROW software, including the input data and installation requirements, and guidance in the specification, calibration, and application of basic SPARROW models, as well as descriptions of the model output and its interpretation. The documentation is intended for both researchers and water-resource managers with interest in using the results of existing models and developing and applying new SPARROW models. The documentation of the model is presented in two parts. Part 1 provides a theoretical and practical introduction to SPARROW modeling techniques, which includes a discussion of the objectives, conceptual attributes, and model infrastructure of SPARROW. Part 1 also includes background on the commonly used model specifications and the methods for estimating and evaluating parameters, evaluating model fit, and generating water-quality predictions and measures of uncertainty. Part 2 provides a user's guide to SPARROW, which includes a discussion of the software architecture and details of the model input requirements and output files, graphs, and maps. The text documentation and computer software are available on the Web at http://usgs.er.gov/sparrow/sparrow-mod.html .

Techniques and Methods↗

Water, Energy, and Biogeochemical Model (WEBMOD), user’s manual, version 1

The Water, Energy, and Biogeochemical Model (WEBMOD) uses the framework of the U.S. Geological Survey (USGS) Modular Modeling System to simulate fluxes of water and solutes through watersheds. WEBMOD divides watersheds into model response units (MRU) where fluxes and reactions are simulated for the following eight hillslope reservoir types: canopy; snowpack; ponding on impervious surfaces; O-horizon; two reservoirs in the unsaturated zone, which represent preferential flow and matrix flow; and two reservoirs in the saturated zone, which also represent preferential flow and matrix flow. The reservoir representing ponding on impervious surfaces, currently not functional (2016), will be implemented once the model is applied to urban areas. MRUs discharge to one or more stream reservoirs that flow to the outlet of the watershed. Hydrologic fluxes in the watershed are simulated by modules derived from the USGS Precipitation Runoff Modeling System; the National Weather Service Hydro-17 snow model; and a topography-driven hydrologic model (TOPMODEL). Modifications to the standard TOPMODEL include the addition of heterogeneous vertical infiltration rates; irrigation; lateral and vertical preferential flows through the unsaturated zone; pipe flow draining the saturated zone; gains and losses to regional aquifer systems; and the option to simulate baseflow discharge by using an exponential, parabolic, or linear decrease in transmissivity. PHREEQC, an aqueous geochemical model, is incorporated to simulate chemical reactions as waters evaporate, mix, and react within the various reservoirs of the model. The reactions that can be specified for a reservoir include equilibrium reactions among water; minerals; surfaces; exchangers; and kinetic reactions such as kinetic mineral dissolution or precipitation, biologically mediated reactions, and radioactive decay. WEBMOD also simulates variations in the concentrations of the stable isotopes deuterium and oxygen-18 as a result of varying inputs, mixing, and evaporation. This manual describes the WEBMOD input and output files, along with the algorithms and procedures used to simulate the hydrology and water quality in a watershed. Examples are presented that demonstrate hydrologic processes, weathering reactions, and isotopic evolution in an alpine watershed and the effect of irrigation on water flows and salinity in an intensively farmed agricultural area.

Techniques and Methods↗

PDEPTH—A computer program for the geophysical interpretation of magnetic and gravity profiles through Fourier filtering, source-depth analysis, and forward modeling

PDEPTH is an interactive, graphical computer program used to construct interpreted geological source models for observed potential-field geophysical profile data. The current version of PDEPTH has been adapted to the Windows platform from an earlier DOS-based version. The input total-field magnetic anomaly and vertical gravity anomaly profiles can be filtered to produce derivative products such as reduced-to-pole magnetic profiles, pseudogravity profiles, pseudomagnetic profiles, and upward-or-downward-continued profiles. A variety of source-location methods can be applied to the original and filtered profiles to estimate (and display on a cross section) the locations and physical properties of contacts, sheet edges, horizontal line sources, point sources, and interface surfaces. Two-and-a-half-dimensional source bodies having polygonal cross sections can be constructed using a mouse and keyboard. These bodies can then be adjusted until the calculated gravity and magnetic fields of the source bodies are close to the observed profiles. Auxiliary information such as the topographic surface, bathymetric surface, seismic basement, and geologic contact locations can be displayed on the cross section using optional input files. Test data files, used to demonstrate the source location methods in the report, and several utility programs are included.

Techniques and Methods↗

Analytical solutions for one-, two-, and three-dimensional solute transport in ground-water systems with uniform flow

Analytical solutions to the advective-dispersive solute-transport equation are useful in predicting the fate of solutes in ground water. Analytical solutions compiled from available literature or derived by the author are presented for a variety of boundary condition types and solute-source configurations in one-, two-, and three-dimensional systems having uniform ground-water flow. A set of user-oriented computer programs was created to evaluate these solutions and to display the results in tabular and computer-graphics format. These programs incorporate many features that enhance their accuracy, ease of use, and versatility. Documentation for the programs describes their operation and required input data, and presents the results of sample problems. Derivations of selected solutions, source codes for the computer programs, and samples of program input and output also are included.

Techniques of Water-Resources Investigations↗

User's guide to SEAWAT; a computer program for simulation of three-dimensional variable-density ground-water flow

The SEAWAT program was developed to simulate three-dimensional, variable-density, transient ground-water flow in porous media. The source code for SEAWAT was developed by combining MODFLOW and MT3DMS into a single program that solves the coupled flow and solute-transport equations. The SEAWAT code follows a modular structure, and thus, new capabilities can be added with only minor modifications to the main program. SEAWAT reads and writes standard MODFLOW and MT3DMS data sets, although some extra input may be required for some SEAWAT simulations. This means that many of the existing pre- and post-processors can be used to create input data sets and analyze simulation results. Users familiar with MODFLOW and MT3DMS should have little difficulty applying SEAWAT to problems of variable-density ground-water flow. MODFLOW was modified to solve the variable-density flow equation by reformulating the matrix equations in terms of fluid mass rather than fluid volume and by including the appropriate density terms. Fluid density is assumed to be solely a function of the concentration of dissolved constituents; the effects of temperature on fluid density are not considered. Temporally and spatially varying salt concentrations are simulated in SEAWAT using routines from the MT3DMS program. SEAWAT uses either an explicit or implicit procedure to couple the ground-water flow equation with the solute-transport equation. With the explicit procedure, the flow equation is solved first for each timestep, and the resulting advective velocity field is then used in the solution to the solute-transport equation. This procedure for alternately solving the flow and transport equations is repeated until the stress periods and simulation are complete. With the implicit procedure for coupling, the flow and transport equations are solved multiple times for the same timestep until the maximum difference in fluid density between consecutive iterations is less than a user-specified tolerance. The SEAWAT code was tested by simulating five benchmark problems involving variable-density ground-water flow. These problems include two box problems, the Henry problem, Elder problem, and HYDROCOIN problem. The purpose of the box problems is to verify that fluid velocities are properly calculated by SEAWAT. For each of the box problems, SEAWAT calculates the appropriate velocity distribution. SEAWAT also accurately simulates the Henry problem, and SEAWAT results compare well with those of SUTRA. The Elder problem is a complex flow system in which fluid flow is driven solely by density variations. Results from SEAWAT, for six different times, compare well with results from Elder's original solution and results from SUTRA. The HYDROCOIN problem consists of fresh ground water flowing over a salt dome. Simulated contours of salinity compare well for SEAWAT and MOCDENSE.

Techniques of Water-Resources Investigations↗

Finite difference model for aquifer simulation in two dimensions with results of numerical experiments

The model will simulate ground-water flow in an artesian aquifer, a water-table aquifer, or a combined artesian and water-table aquifer. The aquifer may be heterogeneous and anisotropic and have irregular boundaries. The source term in the flow equation may include well discharge, constant recharge, leakage from confining beds in which the effects of storage are considered, and evapotranspiration as a linear function of depth to water. The theoretical development includes presentation of the appropriate flow equations and derivation of the finite-difference approximations (written for a variable grid). The documentation emphasizes the numerical techniques that can be used for solving the simultaneous equations and describes the results of numerical experiments using these techniques. Of the three numerical techniques available in the model, the strongly implicit procedure, in general, requires less computer time and has fewer numerical difficulties than do the iterative alternating direction implicit procedure and line successive overrelaxation (which includes a two-dimensional correction procedure to accelerate convergence). The documentation includes a flow chart, program listing, an example simulation, and sections on designing an aquifer model and requirements for data input. It illustrates how model results can be presented on the line printer and pen plotters with a program that utilizes the graphical display software available from the Geological Survey Computer Center Division. In addition the model includes options for reading input data from a disk and writing intermediate results on a disk.

Techniques of Water-Resources Investigations↗

A model for simulation of flow in singular and interconnected channels

A one-dimensional numerical model is presented for simulating the unsteady flow in singular riverine or estuarine reaches and in networks of reaches composed of interconnected channels. The model is both general and flexible in that it can be used to simulate a wide range of flow conditions for various channel configurations. The channel geometry of the network to be modeled should be sufficiently simple so as to lend itself to characterization in one spatial dimension. The flow must be substantially homogenous in density, and hydrostatic pressure must prevail everywhere in the network channels. The slope of each channel bottom ought to be mild and reasonably constant over its length so that the flow remains subcritical. The model accommodates tributary inflows and diversions and includes the effects of wind shear on the water surface as a forcing function in the flow equations. Water-surface elevations and flow discharges are computed at channel junctions, as well as at specified intermediate locations within the network channels. The one-dimensional branch-network flow model uses a four-point, implicit, finite-difference approximation of the unsteady-flow equations. The flow equations are linearized over a time step, and branch transformations are formulated that describe the relationship between the unknowns at the end points of the channels. The resultant matrix of branch-transformation equations and required boundary-condition equations is solved by Gaussian elimination using maximum pivot strategy. Five example applications of the flow model are illustrated. The applications cover such diverse conditions as a singular upland river reach in which unsteady flow results from hydropower regulations, coastal rivers composed of sequentially connected reaches subject to unsteady tide-driven flow, and a multiply connected network of channels whose flow is principally governed by wind tides and seiches in adjoining lakes. The report includes a listing of the FORTRAN IV computer program and a description of the input data requirements. Model supporting programs for the processing and input of initial and boundary-value data are identified, various model output formats are illustrated, and instructions are given to permit the production of graphical output using the line printer, electromechanical pen plotters, cathode-ray-tube display units, or microfilm recorders.

Techniques of Water-Resources Investigations↗

Effects of a Cattail Wetland on Water Quality of Irondequoit Creek near Rochester, New York

A 6-year (1990-96) study of the Ellison Park wetland, a 423-acre, predominantly cattail ( Typha glauca ) marsh in Monroe County, N.Y., was conducted to document the effect that this wetland has on the water quality of Irondequoit Creek, which flows through it. Irondequoit Creek drains 151 square miles of mostly urban and suburban land and is the main tributary to Irondequoit Bay on Lake Ontario. The wetland was a sink for total phosphorus and total suspended solids (28 and 47 percent removal efficiencies, respectively, over the 6-year study period). Sedimentation and vegetative filtration appear to be the primary mechanisms for the decrease in loads of these constituents. Total nitrogen loads were decreased slightly by the wetland; removal efficiencies for ammonia-plus-organic nitrogen and nitrate-plus-nitrite were 6 and 3 percent, respectively. The proportions of total phosphorus and total nitrogen constituents were altered by the wetland. Orthophosphate and ammonia nitrogen were generated within the wetland and represented 12 percent of the total phosphorus output load and 1.8 percent of total nitrogen output load, respectively. Conservative chemicals, such as chloride and sulfate, were littleaffected by the wetland. Concentrations of zinc, lead, and cadmium showed statistically significant decreases, which are attributed to sedimentation and filtration of sediment and organic matter to which these elements adsorb. Sediment samples from open-water depositional areas in the wetland contained high concentrations of (1) trace metals, including barium, manganese, strontium, zinc (each of which exceeded 200 parts per million), as well as chromium, copper, lead, and vanadium, and (2) some polycyclic aromatic hydrocarbons. Persistent organochlorine pesticides, such as chlordane, dieldrin, DDT and its degradation products (DDD and DDE), and polychlorinated biphenyls (PCB's), also were detected, but concentrations of these compounds were within the ranges often found in depositional environments in highly urbanized areas. Cattail shoots attained a maximum height of 350 centimeters, a density of more than 30 shoots per square meter, and total biomass of more than 5,600 grams per square meter (46 percent of which was in above-ground tissues during the growing season). Nitrogen and potassium were three times more abundant in above-ground tissues (2.4 and 1.5 percent by dry weight, respectively) than in below-ground tissues (0.8 and 0.5 percent, respectively). Concentrations of phosphorus, molybdenum, and manganese in above-ground tissues were similar to those in below-ground tissues, but the concentrations of all other constituents were considerably higher in below-ground tissues. Concentrations of several elements exceeded those typically found in natural wetlands; these included manganese (417 ppm, parts per million) and sodium (3,600 ppm) in above-ground tissues, and aluminum (1,540 ppm), iron (15,400 ppm), manganese (433 ppm), and sodium (10,000 ppm) in below-ground tissues. Large quantities of nutrients are assimilated by wetland vegetation during the growing season, but neither tissue production nor microbial metabolic processes appeared to play a significant role in the observed patterns of surface-water chemical input-to-output relations on a seasonal basis. Presumably, internal cycling of nutrients sequestered in the sediments and detritus, combined with a summer increase in microbially mediated chemical transformations, obscured the effects of vegetative assimilation during the summer on surface-water chemical loads. Additionally, the natural confinement of most flows within the banks of Irondequoit Creek, which resulted in passage of stormwater through the wetland with little dispersion or detention in the cattail and backwater areas, diminished the capability of the wetland to improve water quality. Additional factors that probably affected the chemical-removal efficiency of the wetland included chemical inflow loading rates, storage and release mechanisms of the sediments (sedimentation, adsorption, filtration, precipitaton, dissolution, and resuspension), and accretion and burial of organic matter. Measurements of chlorophyll_ a concentrations, and calculations of potential phosphorus concentrations, since the 1970’s indicate an improvement in the trophic state of Irondequoit Bay. Estimated average annual loads (1990-96) of selected constituents entering Irondequoit Bay indicate that, since 1980, the loads of all major forms of nitrogen have decreased, chloride loads have increased, and sulfate loads have changed little. Inputs of total phosphorus and suspended solids to the wetland have increased since 1980, possibly as a result of increased erosion by stormflows from an increasingly developed watershed. The wetland decreases the loads of these constituents, but the trends of these loads entering Irondequoit Bay cannot be reliably defined because the removal efficiencies during the two earlier study periods (1980–81 and 1984–88) are known.

Water-Resources Investigations Report↗

Factors affecting nutrient trends in major rivers of the Chesapeake Bay Watershed

Trends in nutrient loads and flow-adjusted concentrations in the major rivers entering Chesapeake Bay were computed on the basis of water-quality data collected between 1985 and 1998 at 29 monitoring stations in the Susquehanna, Potomac, James, Rappahannock, York, Patuxent, and Choptank River Basins. Two computer models-the Chesapeake Bay Watershed Model (WSM) and the U.S. Geological Survey's 'Spatially Referenced Regressions on Watershed attributes' (SPARROW) Model-were used to help explain the major factors affecting the trends. Results from WSM simulations provided information on temporal changes in contributions from major nutrient sources, and results from SPARROW model simulations provided spatial detail on the distribution of nutrient yields in these basins. Additional data on nutrient sources, basin characteristics, implementation of management practices, and ground-water inputs to surface water were analyzed to help explain the trends. The major factors affecting the trends were changes in nutrient sources and natural variations in streamflow. The dominant source of nitrogen and phosphorus from 1985 to 1998 in six of the seven tributary basins to Chesapeake Bay was determined to be agriculture. Because of the predominance of agricultural inputs, changes in agricultural nutrient sources such as manure and fertilizer, combined with decreases in agricultural acreage and implementation of best management practices (BMPs), had the greatest impact on the trends in flow-adjusted nutrient concentrations. Urban acreage and population, however, were noted to be increasing throughout the Chesapeake Bay Watershed, and as a result, delivered loads of nutrients from urban areas increased during the study period. Overall, agricultural nutrient management, in combination with load decreases from point sources due to facility upgrades and the phosphate detergent ban, led to downward trends in flow-adjusted nutrient concentrations atmany of the monitoring stations in the watershed. The loads of nutrients, however, were not reduced significantly at most of the monitoring stations. This is due primarily to higher streamflow in the latter years of the monitoring period, which led to higher loading in those years. Results of this study indicate a need for more detailed information on BMP effectiveness under a full range of hydrologic conditions and in different areas of the watershed; an internally consistent fertilizer data set; greater consideration of the effects of watershed processes on nutrient transport; a refinement of current modeling efforts; and an expansion of the non-tidal monitoring network in the Chesapeake Bay Watershed.

Delaware, Maryland, New York, Pennsylvania, Virgin↗