Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Modeling Report”

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

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

At least 91 records · Page 5Linked to original sources
USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗
USGS2024

Modeling the impacts of Glen Canyon Dam operations on Colorado River resources

At the time of this report, the Bureau of Reclamation (Reclamation) is writing two supplemental Environmental Impact Statements (sEIS ) and a new Environmental Impact Statement (EIS) that will analyze the effects of changing water flow out of Glen Canyon Dam (GCD) (U.S. Department of Interior, 2024). These actions have the potential to affect downstream resources, including threatened and endangered species, in the Grand Canyon. This report covers modeling support provided on the two sEIS by the USGS Grand Canyon Monitoring and Research Center (GCMRC; U.S. Geological Survey, Southwest Biological Science Center). The first sEIS (Interim Guidelines sEIS) modifies Reclamation’s 2007 Colorado River Interim Guidelines for Lower Basin Shortages and Coordinated Operations for Lake Powell and Lake Mead (USBR, 2007) that determines annual water releases from GCD based on inflow to Lake Powell, the power generating requirements of GCD, and relative lake levels of Lake Powell and Lake Mead. Drought conditions have lowered the level of Lake Powell and may require GCD to release less water than was analyzed in the Interim EIS (7-million-acre feet/year). The effect of less water released, as well as lower lake levels and associated water quality concerns, on downstream resources was not analyzed in the 2007 Interim Guidelines EIS. Reclamation requested GCMRC support to provide models predicting the effects to resources of water releases lower than 7M acre feet/year, including models predicting effects to threatened and endangered species for use in a Biological Assessment.

Report↗
USGS2004

Guidelines for Evaluating Ground-Water Flow Models

Ground-water flow modeling is an important tool frequently used in studies of ground-water systems. Reviewers and users of these studies have a need to evaluate the accuracy or reasonableness of the ground-water flow model. This report provides some guidelines and discussion on how to evaluate complex ground-water flow models used in the investigation of ground-water systems. A consistent thread throughout these guidelines is that the objectives of the study must be specified to allow the adequacy of the model to be evaluated.

Scientific Investigations Report↗
USGS2020

An enhanced hydrologic stream network based on the NHDPlus medium resolution dataset

The National Hydrography Dataset Plus, Version 2.1 (NHDPlusV2.1) is an attribute-rich digital stream network for the conterminous United States, serving as a foundational infrastructure for reporting hydrologic information at both regional and national scales. SPAtially Referenced Regressions On Watershed attributes (SPARROW) is a process-based statistical model that relies on a digital hydrologic network like NHDPlusV2.1 to establish spatial relations between quantities of monitored contaminant loads and contaminant sources, accounting for the physical characteristics along flow paths affecting contaminant transport. The U.S. Geological Survey National Water Quality Assessment project adopted and modified the medium-resolution NHDPlusV2.1 network for use as the primary framework supporting SPARROW modeling. This report describes the enhancements made to improve the routing capabilities and the value-added attributes of NHDPlusV2.1 to support modeling and other hydrologic analyses. These enhancements include corrections to inconsistencies in network/routing information, filling in missing attribute values of associated characteristics, accounting of water use affecting flow, new variables useful for interpreting network data, revised flowline attributes such as slope and flow, and incorporation of ancillary spatial data into the network. The resulting dataset containing the enhancements to the network is named E2NHDPlusV2_US. Although the enhancements described in the report were developed for use in SPARROW modeling, the enhancements are expected to be useful for a wide variety of hydrologic studies within the United States.

Conterminous United States↗
USGS2000

The derivation of World Magnetic Model 2000

This report contains a detailed summary of the data used, the modelling techniques employed and the results obtained in the production of the World Magnetic Model 2000. This model is designed for use in air and sea navigation systems and is valid till 2005.0. The derivation of World Magnetic Model 2000 has been the joint responsibility of the British Geological Survey (BGS) and the United States Geological Survey (USGS). The World Magnetic Model is the standard model in UK Ministry of Defence and US Department of Defense, the North Atlantic Treaty Organization (NATO), and the World Hydrographic Office (WHO) navigation and attitude/heading referencing systems. It is also used widely in civilian navigation systems.

British Geological Survey Technical Report↗
USGS2014

Modeling a historical mountain pine beetle outbreak using Landsat MSS and multiple lines of evidence

Mountain pine beetles are significant forest disturbance agents, capable of inducing widespread mortality in coniferous forests in western North America. Various remote sensing approaches have assessed the impacts of beetle outbreaks over the last two decades. However, few studies have addressed the impacts of historical mountain pine beetle outbreaks, including the 1970s event that impacted Glacier National Park. The lack of spatially explicit data on this disturbance represents both a major data gap and a critical research challenge in that wildfire has removed some of the evidence from the landscape. We utilized multiple lines of evidence to model forest canopy mortality as a proxy for outbreak severity. We incorporate historical aerial and landscape photos, aerial detection survey data, a nine-year collection of satellite imagery and abiotic data. This study presents a remote sensing based framework to (1) relate measurements of canopy mortality from fine-scale aerial photography to coarse-scale multispectral imagery and (2) classify the severity of mountain pine beetle affected areas using a temporal sequence of Landsat data and other landscape variables. We sampled canopy mortality in 261 plots from aerial photos and found that insect effects on mortality were evident in changes to the Normalized Difference Vegetation Index (NDVI) over time. We tested multiple spectral indices and found that a combination of NDVI and the green band resulted in the strongest model. We report a two-step process where we utilize a generalized least squares model to account for the large-scale variability in the data and a binary regression tree to describe the small-scale variability. The final model had a root mean square error estimate of 9.8% canopy mortality, a mean absolute error of 7.6% and an R 2 of 0.82. The results demonstrate that a model of percent canopy mortality as a continuous variable can be developed to identify a gradient of mountain pine beetle severity on the landscape.

Montana↗
USGS2023

Updated estimates of water budget components for the Mississippi Embayment Region using a soil-water-balance model, 2000–2020

A Soil-Water-Balance (SWB) model for the Mississippi embayment region in Arkansas, Tennessee, Mississippi, and Louisiana was constructed and calibrated to gain insight into potential recharge patterns for the Mississippi River Valley alluvial aquifer, which has had substantial drawdown under intense pumping stress over the last several decades. An analysis of the net infiltration term from the SWB model combined with newly gathered airborne electromagnetic geophysical data on the surficial sediments in a calibrated modular three-dimensional finite-difference (MODFLOW 6) groundwater flow model of one area in the alluvial plain found that the distribution of net infiltration was significantly different from the recharge that gets to the water table through the complicated silt and clay stratigraphy of the unsaturated zone. The net infiltration of water through the rooting zone as simulated by SWB ranges from 5.7 to 12.3 inches per year in the alluvial plain part of the model domain, and is fairly evenly distributed within local areas. Recharge to the underlying aquifer is less and is much more focused in particular zones where the connectivity through the upper layers of the unsaturated zone above the water table is greater, indicating possible horizontal flow and perched water table conditions in the unsaturated zone. Runoff and net infiltration together account for 32 percent of the incoming precipitation overall and somewhat higher percentages in the alluvial plain area on an annual basis. These terms are much higher in the fall and winter than in the summer. Actual evapotranspiration accounts for between 62 and 72 percent on average of the annual precipitation but dominates all other terms in the summer months. Without irrigation, summertime net infiltration and runoff would be near zero in the crop-dominated alluvial plain area. The SWB model reproduced reported irrigation rates for corn, soybeans, rice, and cotton on an annual basis fairly well. The SWB model for the Mississippi embayment region was calibrated using more than 15,000 observations representing four parts of the calculated water budget: actual evapotranspiration, surface runoff, net infiltration, and irrigation. Using a Monte Carlo approach to determine the uncertainty in the model results stemming from the uncertainty in the model parameters used in the calibration, the uncertainty in the annual actual evapotranspiration values was around 5 percent, whereas the uncertainty in the irrigation, net infiltration, and runoff was around 20 percent.

Mississippi Embayment Region↗
USGS2014

Relations of water-quality constituent concentrations to surrogate measurements in the lower Platte River corridor, Nebraska, 2007 through 2011

The lower Platte River, Nebraska, provides drinking water, irrigation water, and in-stream flows for recreation, wildlife habitat, and vital habitats for several threatened and endangered species. The U.S. Geological Survey (USGS), in cooperation with the Lower Platte River Corridor Alliance (LPRCA) developed site-specific regression models for water-quality constituents at four sites (Shell Creek near Columbus, Nebraska [USGS site 06795500]; Elkhorn River at Waterloo, Nebr. [USGS site 06800500]; Salt Creek near Ashland, Nebr. [USGS site 06805000]; and Platte River at Louisville, Nebr. [USGS site 06805500]) in the lower Platte River corridor. The models were developed by relating continuously monitored water-quality properties (surrogate measurements) to discrete water-quality samples. These models enable existing web-based software to provide near-real-time estimates of stream-specific constituent concentrations to support natural resources management decisions. Since 2007, USGS, in cooperation with the LPRCA, has continuously monitored four water-quality properties seasonally within the lower Platte River corridor: specific conductance, water temperature, dissolved oxygen, and turbidity. During 2007 through 2011, the USGS and the Nebraska Department of Environmental Quality collected and analyzed discrete water-quality samples for nutrients, major ions, pesticides, suspended sediment, and bacteria. These datasets were used to develop the regression models. This report documents the collection of these various water-quality datasets and the development of the site-specific regression models. Regression models were developed for all four monitored sites. Constituent models for Shell Creek included nitrate plus nitrite, total phosphorus, orthophosphate, atrazine, acetochlor, suspended sediment, and Escherichia coli (E. coli) bacteria. Regression models that were developed for the Elkhorn River included nitrate plus nitrite, total Kjeldahl nitrogen, total phosphorus, orthophosphate, chloride, atrazine, acetochlor, suspended sediment, and E. coli. Models developed for Salt Creek included nitrate plus nitrite, total Kjeldahl nitrogen, suspended sediment, and E. coli. Lastly, models developed for the Platte River site included total Kjeldahl nitrogen, total phosphorus, sodium, metolachlor, atrazine, acetochlor, suspended sediment, and E. coli.

Nebraska↗
USGS2016

Visual Basic, Excel-based fish population modeling tool - The pallid sturgeon example

The model presented in this report is a spreadsheet-based model using Visual Basic for Applications within Microsoft Excel ( http://dx.doi.org/10.5066/F7057D0Z ) prepared in cooperation with the U.S. Army Corps of Engineers and U.S. Fish and Wildlife Service. It uses the same model structure and, initially, parameters as used by Wildhaber and others (2015) for pallid sturgeon. The difference between the model structure used for this report and that used by Wildhaber and others (2015) is that variance is not partitioned. For the model of this report, all variance is applied at the iteration and time-step levels of the model. Wildhaber and others (2015) partition variance into parameter variance (uncertainty about the value of a parameter itself) applied at the iteration level and temporal variance (uncertainty caused by random environmental fluctuations with time) applied at the time-step level. They included implicit individual variance (uncertainty caused by differences between individuals) within the time-step level. The interface developed for the model of this report is designed to allow the user the flexibility to change population model structure and parameter values and uncertainty separately for every component of the model. This flexibility makes the modeling tool potentially applicable to any fish species; however, the flexibility inherent in this modeling tool makes it possible for the user to obtain spurious outputs. The value and reliability of the model outputs are only as good as the model inputs. Using this modeling tool with improper or inaccurate parameter values, or for species for which the structure of the model is inappropriate, could lead to untenable management decisions. By facilitating fish population modeling, this modeling tool allows the user to evaluate a range of management options and implications. The goal of this modeling tool is to be a user-friendly modeling tool for developing fish population models useful to natural resource managers to inform their decision-making processes; however, as with all population models, caution is needed, and a full understanding of the limitations of a model and the veracity of user-supplied parameters should always be considered when using such model output in the management of any species.

Open-File Report↗
USGS2015

2014 National Park visitor spending effects: economic contributions to local communities, states, and the nation

The National Park System covers more than 84 million acres and is comprised of more than 401 sites across the Nation. These lands managed by the National Park Service (NPS) serve as recreational destinations for visitors from across the Nation and around the world. On vacations or on day trips, NPS visitors spend time and money in the gateway communities surrounding NPS sites. Spending by NPS visitors generates and supports a considerable amount of economic activity within park gateway economies. The NPS has been measuring and reporting visitor spending and economic effects for the past 25 years. The 2012 analysis marked a major revision to the NPS visitor spending effects analyses, with the development of the Visitor Spending Effects model (VSE model) which replaced the previous Money Generation Model (see Cullinane Thomas et al. (2014) for a description of how the VSE model differs from the previous model). This report provides updated VSE estimates associated with 2014 NPS visitation. System-wide visitation estimates in 2014 increased by 7% (or 19.2 million visits) compared to 2013 (Ziesler, 2015). Visitation in 2014 rebounded from a 2013 decline that included a 16-day government shutdown and many park closures for repairs after Superstorm Sandy hit the Northeast in late 2012. The re-opening of the Washington Monument, some 21 months after it was rocked by an earthquake and repaired, also added to 2014 visitation numbers. Additionally, several national parks saw record-breaking visitation in 2014, including Joshua Tree, Rocky Mountain, Grand Teton and Glacier national parks. This report begins by presenting an overview of economic effects analyses, followed by details about the data and methods used for this analysis and 2014 model updates. Estimates of NPS visitor spending in 2014 and resulting economic effects at the local, state, regional, and national levels are then presented. The report concludes with a description of current data limitations. Park-level spending and economic effects estimates are included in the appendix. New this year, results from the Visitor Spending Effects report series are available online via an interactive tool. Users can explore current year visitor spending, jobs, labor income, value added, and output effects by sector for national, state, and local economies. This interactive tool is available via the NPS Social Science Program webpage at http://www.nature.nps.gov/socialscience/economics.cfm.

Natural Resource Report↗
USGS2009

Geologic Setting and Hydrogeologic Units of the Columbia Plateau Regional Aquifer System, Washington, Oregon, and Idaho

The Columbia Plateau Regional Aquifer System (CPRAS) covers approximately 44,000 square miles of northeastern Oregon, southeastern Washington, and western Idaho. The area supports a $6 billion per year agricultural industry, leading the Nation in production of apples and nine other commodities (State of Washington Office of Financial Management, 2007; U.S. Department of Agriculture, 2007). Groundwater availability in the aquifers of the area is a critical water-resource management issue because the water demand for agriculture, economic development, and ecological needs is high. The primary aquifers of the CPRAS are basalts of the Columbia River Basalt Group (CRBG) and overlying basin-fill sediments. Water-resources issues that have implications for future groundwater availability in the region include (1) widespread water-level declines associated with development of groundwater resources for irrigation and other uses, (2) reduction in base flow to rivers and associated effects on temperature and water quality, and (3) current and anticipated effects of global climate change on recharge, base flow, and ultimately, groundwater availability. As part of a National Groundwater Resources Program, the U.S. Geological Survey began a study of the CPRAS in 2007 with the broad goals of (1) characterizing the hydrologic status of the system, (2) identifying trends in groundwater storage and use, and (3) quantifying groundwater availability. The study approach includes documenting changes in the status of the system, quantifying the hydrologic budget for the system, updating the regional hydrogeologic framework, and developing a groundwater-flow simulation model for the system. The simulation model will be used to evaluate and test the conceptual model of the system and later to evaluate groundwater availability under alternative development and climate scenarios. The objectives of this study were to update the hydrogeologic framework for the CPRAS using the available geologic mapping and well information and to develop a digital, three-dimensional hydrogeologic model that could be used as the basis of a groundwater-flow model. This report describes the principal geologic and hydrogeologic units of the CPRAS and geologic map and well data that were compiled as part of the study. The report also describes simplified regional hydrogeologic sections and unit extent maps that were used to conceptualize the framework prior to development of the digital 3-dimensional framework model.

Idaho, Oregon, Washington↗
USGS2021

Statistical methods for simulating structural stormwater runoff best management practices (BMPs) with the Stochastic Empirical Loading and Dilution Model (SELDM)

This report documents statistics for simulating structural stormwater runoff best management practices (BMPs) with the Stochastic Empirical Loading and Dilution Model (SELDM). The U.S. Geological Survey developed SELDM and the statistics documented in this report in cooperation with the Federal Highway Administration to indicate the risk for stormwater flows, concentrations, and loads to exceed user-selected water-quality goals and the potential effectiveness of mitigation measures to reduce such risks. In SELDM, three treatment variables—hydrograph extension, volume reduction, and water-quality treatment—are simulated by using the trapezoidal distribution and the rank correlation with the associated runoff variables. This report describes methods for calculating the trapezoidal distribution statistics and rank correlation coefficients for these treatment variables and methods for estimating the minimum irreducible concentration (MIC), which is the lowest expected effluent concentration from a BMP site or a category of BMPs. These statistics are different from the statistics commonly used to characterize or compare BMPs; they are designed to provide a stochastic transfer function to approximate the quantity, duration, and quality of BMP effluent given the associated inflow values for a population of storm events. Analyses for this study were done with data extracted from a modified copy of the December 2019 version of the International Stormwater Best Management Practices Database. Statistics for volume reduction, hydrograph extension, and water-quality treatment were developed with selected data. The medians of the best-fit statistics for selected constituents were used to construct generalized cumulative distribution functions for the three treatment variables. For volume reduction and hydrograph extension, selection of a Spearman’s rank correlation coefficient (rho) value that is the average of the median and maximum values for the BMP category may help generate realistic simulation results in SELDM. The median rho value may be selected to help generate realistic simulation results for water-quality treatment variables. Water-quality treatment statistics, including trapezoidal ratios and MIC values, were developed for 51 runoff-quality constituents commonly measured in highway and urban runoff studies. Statistics were calculated for water-quality properties, sediment and solids, nutrients, major and trace inorganic elements, organic compounds, and biologic constituents. Analysis of MIC values provides information to guide professional judgement for selecting values for simulating water quality at sites of interest. The MIC is a lower bound for BMP discharge concentrations and will therefore replace simulated BMP discharge concentrations below the selected value. A new method for estimating MIC values, the lognormal variate of inflow concentrations, was developed in this report and these statistics were calculated for individual constituents and constituent categories. Inflow quality is correlated to MIC values for some constituents, but regional soil concentrations were not strongly correlated to MIC values.

Scientific Investigations Report↗
USGS2023

Groundwater residence times in glacial aquifers—A new general simulation-model approach compared to conventional inset models

Groundwater is important as a drinking-water source and for maintaining base flow in rivers, streams, and lakes. Groundwater quality can be predicted, in part, by its residence time in the subsurface, but the residence-time distribution cannot be measured directly and must be inferred from models. This report compares residence-time distributions from four areas where groundwater flow and travel time were simulated with conventional simulation-inset models (IMs) and with a new automated model-construction method called general simulation models (GSMs). The comparison provides an opportunity to explore controls on travel time and improve the methods used in the creation of GSMs. These models can be useful for three main-use cases: (1) rapid testing of relationships that govern groundwater flow and age, (2) generation of consistent examples for training a machine-learning metamodel, and (3) serving as a starting point for more detailed models. Comparison of the GSMs to IMs indicated a qualified pattern of agreement for residence-time distributions as indicated by the Nash-Sutcliffe efficiency and Spearman’s correlation coefficient. The agreement was best for the median values of the simulated residence times in young fractions of groundwater (defined as the fractions of groundwater in samples less than 65 years old) at the scale of the eight-digit hydrologic-unit code. Generally, the median values of the young fractions in the IMs were correlated with the median values from the GSMs. The relative trends across the four areas also were similar for the other residence-time metrics. The medians of residence-time metrics at finer scales show a fair degree of scatter. The GSM results compared most poorly for median travel times in the older fraction of groundwater (older than 65 years). The GSM approach is intended as a flexible framework for developing models that can be useful individually as screening tools or collectively to support projects in statistical learning. Although one set of GSM algorithms was presented here, the approach can accommodate many types of data and also different categories of prior information. Comparison of GSMs and IMs suggests ways in which the GSMs, while remaining easy to construct and calibrate, can be improved for estimating groundwater travel times. IMs do not yield exact travel times, and matching GSMs to IMs does not guarantee an improvement; however, IMs provide a convenient benchmark against which to explore relations between physical characteristics of watersheds and the distribution of travel times within them. This effort was undertaken as part of the National Water Quality Program of the U.S. Geological Survey to assist in determining the susceptibility of groundwater in glacial aquifers to a variety of natural and anthropogenic contaminants.

Illinois, Indiana, Michigan, Wisconsin↗
USGS2009

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↗
USGS2017

Coding conventions and principles for a National Land-Change Modeling Framework

This report establishes specific rules for writing computer source code for use with the National Land-Change Modeling Framework (NLCMF). These specific rules consist of conventions and principles for writing code primarily in the C and C++ programming languages. Collectively, these coding conventions and coding principles create an NLCMF programming style. In addition to detailed naming conventions, this report provides general coding conventions and principles intended to facilitate the development of high-performance software implemented with code that is extensible, flexible, and interoperable. Conventions for developing modular code are explained in general terms and also enabled and demonstrated through the appended templates for C++ base source-code and header files. The NLCMF limited-extern approach to module structure, code inclusion, and cross-module access to data is both explained in the text and then illustrated through the module templates. Advice on the use of global variables is provided.

Techniques and Methods↗
USGS2018

Preliminary geologic framework developed for a proposed environmental monitoring study of a deep, unconventional Marcellus Shale drill site, Washington County, Pennsylvania

Background In the fall of 2011, the U.S. Geological Survey (USGS) was afforded an opportunity to participate in an environmental monitoring study of the potential impacts of a deep, unconventional Marcellus Shale hydraulic fracturing site. The drill site of the prospective case study is the “Range Resources MCC Partners L.P. Units 1-5H” location (also referred to as the “RR–MCC” drill site), located in Washington County, southwestern Pennsylvania. Specifically, the USGS was approached to provide a geologic framework that would (1) provide geologic parameters for the proposed area of a localized groundwater circulation model, and (2) provide potential information for the siting of both shallow and deep groundwater monitoring wells located near the drill pad and the deviated drill legs. The lead organization of the prospective case study of the RR–MCC drill site was the Groundwater and Ecosystems Restoration Division (GWERD) of the U.S. Environmental Protection Agency. Aside from the USGS, additional partners/participants were to include the Department of Energy, the Pennsylvania Geological Survey, the Pennsylvania Department of Environmental Protection, and the developer Range Resources LLC. During the initial cooperative phase, GWERD, with input from the participating agencies, drafted a Quality Assurance Project Plan (QAPP) that proposed much of the objectives, tasks, sampling and analytical procedures, and documentation of results. Later in 2012, the proposed cooperative agreement between the aforementioned partners and the associated land owners for a monitoring program at the drill site was not executed. Therefore, the prospective case study of the RR–MCC site was terminated and no installation of groundwater monitoring wells nor the collection of nearby soil, stream sediment, and surface-water samples were made. Prior to the completion of the QAPP and termination of the perspective case study the geologic framework was rapidly conducted and nearly completed. This was done for three principal reasons. First, there was an immediate need to know the distribution of the relatively undisturbed surface to near-surface bedrock geology and unconsolidated materials for the collection of baseline surface data prior to drill site development (drill pad access road, drill pad leveling) and later during monitoring associated with well drilling, well development, and well production. Second, it was necessary to know the bedrock geology to support the siting of: (1) multiple shallow groundwater monitoring wells (possibly as many as four) surrounding and located immediately adjacent to the drill pad, and (2) deep groundwater monitoring wells (possibly two) located at distance from the drill pad with one possibly being sited along one of the deviated production drill legs. Lastly, the framework geology would provide the lateral extent, thickness, lithology, and expected discontinuities of geologic units (to be parsed or grouped as hydrostratigraphic units) and regional structure trends as inputs into the groundwater model. This report provides the methodology of geologic data accumulation and aggregation, and its integration into a geographic information system (GIS) based program. The GIS program will allow multiple data to be exported in various formats (shapefiles [.shp], database files [.dbf], and Keyhole Markup Language files [.KML]) for use in surface and subsurface geologic site characterization, for sampling strategies, and for inputs for groundwater modeling.

Pennsylvania↗
USGS2010

One-Dimensional Transport with Equilibrium Chemistry (OTEQ) - A reactive transport model for streams and rivers

OTEQ is a mathematical simulation model used to characterize the fate and transport of waterborne solutes in streams and rivers. The model is formed by coupling a solute transport model with a chemical equilibrium submodel. The solute transport model is based on OTIS, a model that considers the physical processes of advection, dispersion, lateral inflow, and transient storage. The equilibrium submodel is based on MINTEQ, a model that considers the speciation and complexation of aqueous species, acid-base reactions, precipitation/dissolution, and sorption. Within OTEQ, reactions in the water column may result in the formation of solid phases (precipitates and sorbed species) that are subject to downstream transport and settling processes. Solid phases on the streambed may also interact with the water column through dissolution and sorption/desorption reactions. Consideration of both mobile (waterborne) and immobile (streambed) solid phases requires a unique set of governing differential equations and solution techniques that are developed herein. The partial differential equations describing physical transport and the algebraic equations describing chemical equilibria are coupled using the sequential iteration approach. The model's ability to simulate pH, precipitation/dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between instream chemistry and hydrologic transport at the field scale. This report details the development and application of OTEQ. Sections of the report describe model theory, input/output specifications, model applications, and installation instructions. OTEQ may be obtained over the Internet at http://water.usgs.gov/software/OTEQ.

Techniques and Methods↗
USGS2014

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↗