Geology Reports⌕ Search

Geology topics

Richard B. Winston

Publications and source records attributed to Richard B. Winston.

17 recordsLinked to original sources

Revision of ModelMuse to support the use of PEST software with MODFLOW and SUTRA models

Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.

Techniques and Methods↗

Getting Started with MODFLFOW

Numerical modeling of groundwater flow systems was once accessible only to modeling specialists in the hydrogeological community. Software such as MODFLOW—the most frequently used groundwater modeling program in the world—and associated graphical user interfaces (GUIs) have made modeling possible for most groundwater scientists. This book provides the bridge from understanding to implementing models by introducing the basics of MODFLOW version 6 and providing readers who have a working knowledge of groundwater flow with a guide through construction of their first groundwater model.

Book↗

Simulation of water-table response to sea-level rise and change in recharge, Sandy Hook unit, Gateway National Recreation Area, New Jersey

The Sandy Hook Unit, Gateway National Recreation Area (hereafter Sandy Hook) in New Jersey is a 10-kilometer-long spit visited by thousands of people each year who take advantage of the historical and natural resources and recreational opportunities. The historical and natural resources are threatened by global climate change, including sea-level rise (SLR), changes in precipitation and groundwater recharge, and changes in the frequency and severity of coastal storms. Fresh groundwater resources are important to the ecosystems of Sandy Hook. The Bayside Holly Forest, one of only two known old-growth American holly ( Ilex opaca ) maritime forests, is particularly vulnerable to global climate change because of the proximity of the water table to land surface in low-lying areas and the potential for saltwater intrusion and inundation. The shallow groundwater-flow system on Sandy Hook is dominated by recharge from precipitation, fresh groundwater discharge to evapotranspiration (ET), discharge to surface seeps, and submarine groundwater discharge (groundwater discharging directly to the ocean). A three-dimensional groundwater-flow model that simulates the shallow groundwater-flow system and interaction with surrounding saltwater boundaries was constructed to simulate multi-density groundwater flow, treating the freshwater/saltwater transition zone as a sharp interface that represents the half-seawater surface. Groundwater-flow simulations completed for this study include a Baseline scenario, three SLR scenarios (0.2, 0.4, and 0.6 meter [m]), two Recharge scenarios—a 10-percent Increased Recharge scenario and a 10-percent Decreased Recharge scenario—and a scenario with 0.6 m of SLR and 10-percent increase in recharge. The Recharge scenarios indicate the system is not sensitive to a 10-percent increase or decrease in recharge from the Baseline scenario. In the SLR scenarios, SLR causes the water table to rise, resulting in increased fresh groundwater discharge to ET and seeps, and reduced submarine discharge compared to the Baseline scenario. The increased discharge to ET and seeps causes the magnitude of water-table rise to be less than that of SLR, which in turn causes the thickness of the freshwater lens to thin, reducing the depth to the half-seawater surface. Water-table rise associated with SLR diminishes the thickness of the unsaturated zone; comparing the Baseline and the 0.6-m SLR scenarios, the area where the simulated water table is above land surface increases by 50.6 hectares, from about 0.9 to 7.4 percent of the land area of Sandy Hook. Areas where the simulated water table is above land surface are likely to be emergent wetlands and contain freshwater if they are tens of meters or more from the shoreline. The steady-state simulations indicate that the percentage of land where the half-seawater surface is less than 9 m below the water table increases from about 2.5 percent (20 hectares) to about 9 percent (74 hectares) with 0.6 m of SLR. In low-lying areas close to the Sandy Hook Bay shoreline, the half-seawater surface is simulated to be as much as 20 m closer to the water table with SLR of 0.6 m. Transient salinization, if any, of shallow groundwater from increased frequency or severity of storm-driven inundation is not included in the analysis. Natural resources on Sandy Hook, particularly the Bayside Holly Forest, may be adversely affected by the rising water table associated with SLR. Freshwater emergent wetlands may increase in area at the expense of other ecosystem assemblages occurring in or on the edges of low-lying enclosed depressions. Cultural resources close to the water table, such as existing basements of structures, may be adversely affected.

New Jersey↗

ModelMuse Version 4: A graphical user interface for MODFLOW 6

ModelMuse, a graphical user interface for groundwater-modeling software, was modified to support MODFLOW 6. ModelMuse works with two types of spatial discretization in MODFLOW 6: structured grids (DIS) and discretization by vertices (DISV). Quadtree refinement is used to generate a DISV model from a structured-grid model. The locations and weights for ghost nodes used to improve DISV model accuracy are computed automatically by ModelMuse using a new algorithm. ModelMuse does not support other types of DISV grids and unstructured grids. ModelMuse supports options in MODFLOW 6 that designate individual cells as confined or convertible and remove inactive cells associated with discontinuous layers, thereby reducing the computational burden. ModelMuse can specify fully three-dimensional (3D), spatially variable anisotropy in hydraulic conductivity. Although MODFLOW 6 does not support the parameters supported by MODFLOW–2005, ModelMuse provides backward compatibility by allowing ModelMuse parameters to specify scale-factor variables in MODFLOW 6 time-series files within packages that support time-series files. ModelMuse can automatically convert the data for many of the packages from other MODFLOW models to the new data for these packages in MODFLOW 6. Some packages, such as the Streamflow-Routing (SFR) package, changed significantly enough that only a partial conversion is possible. Head and flow observations in older models are also converted to observation locations in the MODFLOW 6 Observation utility. ModelMuse accommodates the ability of MODFLOW 6 to store specific discharge components by allowing the user to visualize the components of a specific discharge on the model grid. ModelMuse supports the versions of MODPATH and ZONEBUDGET compatible with MODFLOW 6.

Scientific Investigations Report↗

ModelArchiver—A program for facilitating the creation of groundwater model archives

ModelArchiver is a program designed to facilitate the creation of groundwater model archives that meet the requirements of the U.S. Geological Survey (USGS) policy (Office of Groundwater Technical Memorandum 2016.02, https://water.usgs.gov/admin/memo/GW/gw2016.02.pdf , https://water.usgs.gov/ogw/policy/gw-model/ ). ModelArchiver version 1.0 leads the user step-by-step through the process of creating a USGS groundwater model archive. The user specifies the contents of each of the subdirectories within the archive and provides descriptions of the archive contents. Descriptions of some files can be specified automatically using file extensions. Descriptions also can be specified individually. Those descriptions are added to a readme.txt file provided by the user. ModelArchiver moves the content of the archive to the archive folder and compresses some folders into .zip files. As part of the archive, the modeler must create a metadata file describing the archive. The program has a built-in metadata editor and provides links to websites that can aid in creation of the metadata. The built-in metadata editor is also available as a stand-alone program named FgdcMetaEditor version 1.0, which also is described in this report. ModelArchiver updates the metadata file provided by the user with descriptions of the files in the archive. An optional archive list file generated automatically by ModelMuse can streamline the creation of archives by identifying input files, output files, model programs, and ancillary files for inclusion in the archive.

Open-File Report↗

Volume-weighted particle-tracking method for solute-transport modeling; Implementation in MODFLOW–GWT

In the traditional method of characteristics for groundwater solute-transport models, advective transport is represented by moving particles that track concentration. This approach can lead to global mass-balance problems because in models of aquifers having complex boundary conditions and heterogeneous properties, particles can originate in cells having different pore volumes and (or) be introduced (or removed) at cells representing fluid sources (or sinks) of varying strengths. Use of volume-weighted particles means that each particle tracks solute mass. In source or sink cells, the changes in particle weights will match the volume of water added or removed through external fluxes. This enables the new method to conserve mass in source or sink cells as well as globally. This approach also leads to potential efficiencies by allowing the number of particles per cell to vary spatially—using more particles where concentration gradients are high and fewer where gradients are low. The approach also eliminates the need for the model user to have to distinguish between “weak” and “strong” fluid source (or sink) cells. The new model determines whether solute mass added by fluid sources in a cell should be represented by (1) new particles having weights representing appropriate fractions of the volume of water added by the source, or (2) distributing the solute mass added over all particles already in the source cell. The first option is more appropriate for the condition of a strong source; the latter option is more appropriate for a weak source. At sinks, decisions whether or not to remove a particle are replaced by a reduction in particle weight in proportion to the volume of water removed. A number of test cases demonstrate that the new method works well and conserves mass. The method is incorporated into a new version of the U.S. Geological Survey’s MODFLOW–GWT solute-transport model.

Techniques and Methods↗

Visualization of groundwater withdrawals

Generating an informative display of groundwater withdrawals can sometimes be difficult because the symbols for closely spaced wells can overlap. An alternative method for displaying groundwater withdrawals is to generate a “footprint” of the withdrawals. WellFootprint version 1.0 implements the Footprint algorithm with two optional variations that can speed up the footprint calculation. ModelMuse has been modified in order to generate the input for WellFootprint and to read and graphically display the output from WellFootprint.

Open-File Report↗

Performance assessments of a novel well design for reducing exposure to bedrock‐derived arsenic

Arsenic in groundwater is a serious problem in New England, particularly for domestic well owners drawing water from bedrock aquifers. The overlying glacial aquifer generally has waters with low arsenic concentrations but is less used because of frequent loss of well water during dry periods and the vulnerability to surface‐sourced bacterial contamination. An alternative, novel design for shallow wells in glacial aquifers is intended to draw water primarily from unconsolidated glacial deposits, while being resistant to drought conditions and surface contamination. Its use could greatly reduce exposure to arsenic through drinking water for domestic use. Hypothetical numerical models were used to investigate the potential hydraulic performance of the new well design in reducing arsenic exposure. The aquifer system was divided into two parts, an upper section representing the glacial sediments and a lower section representing the bedrock. The location of the well, recharge conditions, and hydraulic properties were systematically varied in a series of simulations and the potential for arsenic contamination was quantified by analyzing groundwater flow paths to the well. The greatest risk of arsenic contamination occurred when the hydraulic conductivity of the bedrock aquifer was high, or where there was upward flow from the bedrock aquifer because of the position of the well in the flow system.

Groundwater↗

The role of uplift and erosion in the persistence of saline groundwater in the shallow subsurface

In many regions of the world, the shallow (<300 m) subsurface is replenished with meteoric recharge within a few centuries or millennia, but in some regions saline groundwater persists despite abundant rainfall. Analyses of the flushing rate of shallow groundwater usually consider the permeability and recharge rate and a static landscape. The influence of landscape evolution can become important over millions of years, however. Here we present numerical simulations of fluid flow and transport in the top 1 km of a sedimentary foreland basin dominated by aquitards, where the rate of uplift and erosion (20 m Ma −1 ) balances that of meteoric flushing. Paleozoic age saline groundwater and brine persist at shallow depths that might otherwise have contained potable water. Similar hydrogeologic conditions, and uplift and erosion rates, likely exist in many other regions of the world, where a moving landscape has probably never been considered as an important contributor to groundwater quality.

Geophysical Research Letters↗

Modifications made to ModelMuse to add support for the Saturated-Unsaturated Transport model (SUTRA)

This report (1) describes modifications to ModelMuse,as described in U.S. Geological Survey (USGS) Techniques and Methods (TM) 6–A29 (Winston, 2009), to add support for the Saturated-Unsaturated Transport model (SUTRA) (Voss and Provost, 2002; version of September 22, 2010) and (2) supplements USGS TM 6–A29. Modifications include changes to the main ModelMuse window where the model is designed, addition of methods for generating a finite-element mesh suitable for SUTRA, defining how some functions shouldapply when using a finite-element mesh rather than a finite-difference grid (as originally programmed in ModelMuse), and applying spatial interpolation to angles. In addition, the report describes ways of handling objects on the front view of the model and displaying data. A tabulation contains a summary of the new or modified dialog boxes.

Techniques and Methods↗

ListingAnalyst: A program for analyzing the main output file from MODFLOW

ListingAnalyst is a Windows® program for viewing the main output file from MODFLOW-2005, MODFLOW-NWT, or MODFLOW-LGR. It organizes and displays large files quickly without using excessive memory. The sections and subsections of the file are displayed in a tree-view control, which allows the user to navigate quickly to desired locations in the files. ListingAnalyst gathers error and warning messages scattered throughout the main output file and displays them all together in an error and a warning tab. A grid view displays tables in a readable format and allows the user to copy the table into a spreadsheet. The user can also search the file for terms of interest.

Ground Water↗

ModelMuse - A Graphical User Interface for MODFLOW-2005 and PHAST

ModelMuse is a graphical user interface (GUI) for the U.S. Geological Survey (USGS) models MODFLOW-2005 and PHAST. This software package provides a GUI for creating the flow and transport input file for PHAST and the input files for MODFLOW-2005. In ModelMuse, the spatial data for the model is independent of the grid, and the temporal data is independent of the stress periods. Being able to input these data independently allows the user to redefine the spatial and temporal discretization at will. This report describes the basic concepts required to work with ModelMuse. These basic concepts include the model grid, data sets, formulas, objects, the method used to assign values to data sets, and model features. The ModelMuse main window has a top, front, and side view of the model that can be used for editing the model, and a 3-D view of the model that can be used to display properties of the model. ModelMuse has tools to generate and edit the model grid. It also has a variety of interpolation methods and geographic functions that can be used to help define the spatial variability of the model. ModelMuse can be used to execute both MODFLOW-2005 and PHAST and can also display the results of MODFLOW-2005 models. An example of using ModelMuse with MODFLOW-2005 is included in this report. Several additional examples are described in the help system for ModelMuse, which can be accessed from the Help menu.

Techniques and Methods↗

BAT3 Analyzer: Real-time data display and interpretation software for the multifunction bedrock-aquifer transportable testing tool (BAT3)

The BAT3 Analyzer provides real-time display and interpretation of fluid pressure responses and flow rates measured during geochemical sampling, hydraulic testing, or tracer testing conducted with the Multifunction Bedrock-Aquifer Transportable Testing Tool (BAT3) (Shapiro, 2007). Real-time display of the data collected with the Multifunction BAT3 allows the user to ensure that the downhole apparatus is operating properly, and that test procedures can be modified to correct for unanticipated hydraulic responses during testing. The BAT3 Analyzer can apply calibrations to the pressure transducer and flow meter data to display physically meaningful values. Plots of the time-varying data can be formatted for a specified time interval, and either saved to files, or printed. Libraries of calibrations for the pressure transducers and flow meters can be created, updated and reloaded to facilitate the rapid set up of the software to display data collected during testing with the Multifunction BAT3. The BAT3 Analyzer also has the functionality to estimate calibrations for pressure transducers and flow meters using data collected with the Multifunction BAT3 in conjunction with corroborating check measurements. During testing with the Multifunction BAT3, and also after testing has been completed, hydraulic properties of the test interval can be estimated by comparing fluid pressure responses with model results; a variety of hydrogeologic conceptual models of the formation are available for interpreting fluid-withdrawal, fluid-injection, and slug tests.

Open-File Report↗

GoPhast: A graphical user interface for PHAST

GoPhast is a graphical user interface (GUI) for the USGS model PHAST. PHAST simulates multicomponent, reactive solute transport in three-dimensional, saturated, ground-water flow systems. PHAST can model both equilibrium and kinetic geochemical reactions. PHAST is derived from HST3D (flow and transport) and PHREEQC (geochemical calculations). The flow and transport calculations are restricted to constant fluid density and constant temperature. The complexity of the input required by PHAST makes manual construction of its input files tedious and error-prone. GoPhast streamlines the creation of the input file and helps reduce errors. GoPhast allows the user to define the spatial input for the PHAST flow and transport data file by drawing points, lines, or polygons on top, front, and side views of the model domain. These objects can have up to two associated formulas that define their extent perpendicular to the view plane, allowing the objects to be three-dimensional. Formulas are also used to specify the values of spatial data (data sets) both globally and for individual objects. Objects can be used to specify the values of data sets independent of the spatial and temporal discretization of the model. Thus, the grid and simulation periods for the model can be changed without respecifying spatial data pertaining to the hydrogeologic framework and boundary conditions. This report describes the operation of GoPhast and demonstrates its use with examples. GoPhast runs on Windows 2000, Windows XP, and Linux operating systems.

Techniques and Methods↗

SutraGUI, a graphical-user interface for SUTRA, a model for ground-water flow with solute or energy transport

This report describes SutraGUI, a flexible graphical user-interface (GUI) that supports two-dimensional (2D) and three-dimensional (3D) simulation with the U.S. Geological Survey (USGS) SUTRA ground-water-flow and transport model (Voss and Provost, 2002). SutraGUI allows the user to create SUTRA ground-water models graphically. SutraGUI provides all of the graphical functionality required for setting up and running SUTRA simulations that range from basic to sophisticated, but it is also possible for advanced users to apply programmable features within Argus ONE to meet the unique demands of particular ground-water modeling projects. SutraGUI is a public-domain computer program designed to run with the proprietary Argus ONE? package, which provides 2D Geographic Information System (GIS) and meshing support. For 3D simulation, GIS and meshing support is provided by programming contained within SutraGUI. When preparing a 3D SUTRA model, the model and all of its features are viewed within Argus 1 in 2D projection. For 2D models, SutraGUI is only slightly changed in functionality from the previous 2D-only version (Voss and others, 1997) and it provides visualization of simulation results. In 3D, only model preparation is supported by SutraGUI, and 3D simulation results may be viewed in SutraPlot (Souza, 1999) or Model Viewer (Hsieh and Winston, 2002). A comprehensive online Help system is included in SutraGUI. For 3D SUTRA models, the 3D model domain is conceptualized as bounded on the top and bottom by 2D surfaces. The 3D domain may also contain internal surfaces extending across the model that divide the domain into tabular units, which can represent hydrogeologic strata or other features intended by the user. These surfaces can be non-planar and non-horizontal. The 3D mesh is defined by one or more 2D meshes at different elevations that coincide with these surfaces. If the nodes in the 3D mesh are vertically aligned, only a single 2D mesh is needed. For nonaligned meshes, two or more 2D meshes of similar connectivity are used. Between each set of 2D meshes (and model surfaces), the vertical space in the 3D mesh is evenly divided into a user-specified number of layers of finite elements. Boundary conditions may be specified for 3D models in SutraGUI using a variety of geometric shapes that may be located freely within the 3D model domain. These shapes include points, lines, sheets, and solids. These are represented by 2D contours (within the vertically-projected Argus ONE view) with user-defined elevations. In addition, boundary conditions may be specified for 3D models as points, lines, and areas that are located exactly within the surfaces that define the model top and the bottoms of the tabular units. Aquifer properties may be specified separately for each tabular unit. If the aquifer properties vary vertically within a unit, SutraGUI provides the Sutra_Z function that can be used to specify such variation.

Open-File Report↗

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

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

Open-File Report↗