Geology topics

Michael N. Fienen

Publications and source records attributed to Michael N. Fienen.

At least 55 records · Page 3Linked to original sources

USGS2021

Three-dimensional distribution of residence time metrics in the glaciated United States using metamodels trained on general numerical models

Residence time distribution (RTD) is a critically important characteristic of groundwater flow systems; however, it cannot be measured directly. RTD can be inferred from tracer data with analytical models (few parameters) or with numerical models (many parameters). The second approach permits more variation in system properties but is used less frequently than the first because large‐scale numerical models can be resource intensive. Using a novel automated approach, a set of 115 inexpensive general simulation models (GSMs) was used to create RTD metrics (fraction of young groundwater, defined as < 65 years old; mean travel time of young fraction; median travel time of old fraction; and mean path length). GSMs captured the general trends in measured tritium concentrations in 431 wells. Boosted Regression Tree metamodels were trained to predict these RTD metrics using available wall‐to‐wall hydrogeographic digital sets as explanatory features. The metamodels produced a three‐dimensional distribution of predictions throughout the glacial system that generally matched with the numerical model RTD metrics. In addition to the expected importance of aquifer thickness and recharge rate in predicting RTD metrics, two new data sets, Multi‐Order Hydrologic Position (MOHP) and hydrogeologic terrane were important predictors. These variables by themselves produced metamodels with Nash‐Sutcliffe efficiency close to the full metamodel. Metamodel predictions showed that the volume of young groundwater stored in the glaciated U.S. is about 6,000 km 3 , or about 0.5% of globally stored young groundwater.

Water Resources Research↗

USGS2021

Towards improved environmental modeling outcomes: Enabling low-cost access to high-dimensional, geostatistical-based decision-support analyses

Computer models of environmental systems routinely inform decision making for water resource management. In this context, quantifying uncertainty in the important simulated outputs, and reducing uncertainty through assimilating historic system-state observations, is as important as the numerical model. However, implementing high-dimensional and stochastic workflows are challenging, often requiring that practitioners have theoretical and practical understanding of several advanced topics. Worse, implementing these important analyses can take substantial time and effort. This additional effort is often cited as justification for postponing, or even forgoing, these analyses. Herein, we present scripting tools to facilitate the efficient and repeatable construction of high-dimensional, geostatistical-based PEST interfaces, including uncertainty analyses. As demonstrated, these tools can be applied with minimal effort to a model with varied temporal and spatial discretization. Ultimately, these tools can enable low-cost access to valuable decision-support analyses earlier and more frequently during the modeling workflow.

Environmental Modelling & Software↗

USGS2021

Appendix C: Central sands lakes study technical report: Modeling documentation

This report provides the necessary documentation of the numerical models developed for the Central Sands Lake study in central Wisconsin and will be included as a technical appendix in the report to the Wisconsin State Legislature by the Wisconsin Department of Natural Resources (WDNR) in response to 2017 Wisconsin Act 10. This legislation directed WDNR to determine whether existing and potential groundwater withdrawals are causing or are likely to cause significant reduction of mean seasonal water levels at Pleasant Lake, Long Lake, and Plainfield Lake (s. 281.34(7m)(2)(b), Wis. Stats.) in Waushara County, Wisconsin. To evaluate the potential hydrologic connection between groundwater withdrawals and the nearby study lakes, hydrologic models were created that focused on the lakes of interest and yet were large enough to cover a broad enough region to extend to the major hydrologic boundaries of the natural flow system. The areas near the lakes require finer-scale grid discretization (or spacing) to better represent the lakes and streams in the model, but also need to cover a large enough area to include the groundwater withdrawal locations that have the potential to cause reduction in water levels in the lakes. To accomplish these goals, three groundwater models were created: a regional model extending to major hydrologic boundaries; and two inset models, inheriting boundaries from the regional model but focused near the lakes. Each of the inset models, in turn, included a detailed area close to the lakes surrounded by an area at the same spatial scale as the regional model (Figure 1). To support WDNR in evaluating the connection between groundwater withdrawals and lake levels, a representative time period was required over which to compare land use with and without irrigated agriculture and for WDNR to evaluate potential lake stage and flux changes related to irrigated agriculture. WDNR chose the climate period of 1981-2018 to be representative of a typical period and provided two land use scenarios—one with no irrigated agriculture and one with assumed crop rotations similar to current conditions—to simulate with groundwater models to, then, compare lake responses with. As a result, simulations over this climate record are not intended to recreate the history of 1981-2018 because land use changed over that time. These runs are, instead, intended to provide a basis on which to compare land use with and without irrigation-related groundwater withdrawals based on the current arrangement of land use and a varied climatic record. Groundwater withdrawals focused on irrigated-agriculture-related water use because greater than 95% of groundwater withdrawal in the two inset models around the study lakes is for irrigated agriculture water use. The period of 2012-2018 was used for parameter estimation (synonymously referred to as “history matching”) for the groundwater models. This time period was chosen because it includes the most complete water use records to simulate groundwater withdrawals. History matching was performed using groundwater elevations, lake stages, and streamflow observations over the 2012-2018 time period and processed observations derived from those raw data. Climatic data were incorporated into the model using a soil-water balance approach. A soil water balance model was constructed at the scale of the regional groundwater model to both calculate recharge based on land use and climate, and in the long-term climate-period runs, to estimate water use required by irrigated agriculture to apply as well boundary conditions in the groundwater model in the absence of reported water use values over that period.

Wisconsin↗

USGS2020

Approaches to highly parameterized inversion: PEST++ Version 5, a software suite for parameter estimation, uncertainty analysis, management optimization and sensitivity analysis

PEST++ Version 5 extends and enhances the functionality of the PEST++ Version 3 software suite, providing environmental modeling practitioners access to updated Version 3 tools as well as new tools to support decision making with environmental models. Version 5 of PEST++ includes tools for global sensitivity analysis (PESTPP-SEN); least-squares parameter estimation with integrated first-order, second-moment parameter and forecast uncertainty estimation (PESTPP-GLM); an iterative, localized ensemble smoother (PESTPP-IES); and a tool for management optimization under uncertainty (PESTPP-OPT). Additionally, all PEST++ Version 5 tools have a built-in fault-tolerant, multithreaded parallel run manager and are model independent, using the same protocol as the widely used PEST software suite. PEST++ Version 5 is consistent with PEST++ Version 3 conventions and design philosophy. The software’s emphasis continues to target efficient and optimized algorithms that have proven beneficial in decision-support settings and can accommodate large, highly parameterized problems. Expanded and new capabilities are now available to express uncertainty using Monte Carlo and analytical uncertainty approaches and allow evaluation of thousands to millions of parameters. New management optimization capabilities in Version 5 also allow environmental models to be used to answer management questions using multiple societal constraints in a risk-based framework. The PEST++ Version 5 software suite can be compiled for Microsoft Windows® and Unix-based operating systems such as Apple and Linux®; the source code is available with a Microsoft Visual Studio® 2019 solution; and CMake support for all three operating system is also provided. PEST++ Version 5 continues to build a foundation for an open-source framework capable of producing model-independent, robust, and efficient decision-support tools for large environmental models. The functionality of each of the PEST++ tools are demonstrated on a simple example problem. Implications of decisions used when using the PEST++ suite tools are also discussed.

Techniques and Methods↗

USGS2020

Simulation of groundwater flow in the regional aquifer system on Long Island, New York, for pumping and recharge conditions in 2005–15

A three-dimensional groundwater-flow model was developed for the aquifer system of Long Island, New York, to evaluate (1) responses of the hydrologic system to changes in natural and anthropogenic hydraulic stresses, (2) the subsurface distribution of groundwater age, and (3) the regional-scale distribution of groundwater travel times and the source of water to fresh surface waters and coastal receiving waters. The model also provides the groundwater flow components used to define model boundaries for possible inset models used for local-scale analyses. The three-dimensional, groundwater flow model developed for this investigation uses the numerical code MODFLOW–NWT to represent steady-state conditions for average groundwater pumping and aquifer recharge for 2005–15. The particle-tracking algorithm MODPATH, which simulates advective transport in the aquifer, was used to estimate groundwater age, delineate the areas at the water table that contribute recharge to coastal and freshwater bodies, and estimate total travel times of water from the water table to discharge locations. A three-dimensional, 1-meter (3.3-foot) topobathymetric model was used to determine land-surface altitudes for the island and seabed altitudes for the surrounding coastal waters. The mapped extents and surface altitudes of major geologic units were compiled and used to develop a three-dimensional hydrogeologic framework of the aquifer system, including aquifers and confining units. Lithologic data from deep boreholes and previous aquifer-test results were used to estimate the three-dimensional distribution of hydraulic conductivity in principal aquifers. Natural recharge from precipitation was estimated for 2005–15 using a modified Thornthwaite-Mather methodology as implemented in a soil-water balance model. Components of anthropogenic recharge—wastewater return flow, storm water inflow, and inflow from leaky infrastructure—also were estimated for 2005–15. Groundwater withdrawals for various sources, including public water supply, industrial, remediation, and agricultural, were compiled or estimated for the same period. These data were incorporated into the model development to represent the aquifer system geometry, boundaries, and initial hydraulic properties of the regional aquifers and confining units within the Long Island aquifer system. Average hydraulic conditions—water levels and streamflows—for 2005–15 were estimated using existing data from the U.S. Geological Survey National Water Information System database. Model inputs were adjusted to best match average hydrologic conditions using inverse methods as implemented in the parameter-estimating software PEST. The calibrated model was used to simulate average hydrologic conditions in the aquifer system for 2005–15. About 656 cubic feet per second of water was withdrawn on average annually for 2005–15 for water supply and an average of about 349 cubic feet per second of water recharged the aquifer annually from return flow and leaky infrastructure. Parts of New York City have drawdowns exceeding 25 feet, mostly because of urbanization and associated large decreases in recharge rates. Large areas in the western part of the island have drawdowns exceeding 10 feet, mostly from large groundwater withdrawals and sewering, which largely eliminates wastewater return flow. Water-table altitudes in eastern parts of the island increased by more than 2 feet in some areas as a result of wastewater return flow in unsewered areas and changes in land use. Changes in streamflows show a similar pattern as water-table altitudes. Streamflows generally decrease in western parts of the island where there are large drawdowns and increase in eastern parts of the island where water-table altitudes increase.

New York↗

USGS2020

Prioritizing river basins for intensive monitoring and assessment by the US Geological Survey

The US Geological Survey (USGS) is currently (2020) integrating its water science programs to better address the nation’s greatest water resource challenges now and into the future. This integration will rely, in part, on data from 10 or more intensively monitored river basins from across the USA. A team of USGS scientists was convened to develop a systematic, quantitative approach to prioritize candidate basins for this monitoring investment to ensure that, as a group, the 10 basins will support the assessment and forecasting objectives of the major USGS water science programs. Candidate basins were the level-4 hydrologic units (HUC04) with some of the smaller HUC04s being combined; median candidate-basin area is 46,600 km 2 . Candidate basins for the contiguous United States (CONUS) were grouped into 18 hydrologic regions. Ten geospatial variables representing land use, climate change, water use, water-balance components, streamflow alteration, fire risk, and ecosystem sensitivity were selected to rank candidate basins within each of the 18 hydrologic regions. The two highest ranking candidate basins in each of the 18 regions were identified as finalists for selection as “Integrated Water Science Basins”; final selection will consider input from a variety of stakeholders. The regional framework, with only one basin selected per region, ensures that as a group, the basins represent the range in major drivers of the hydrologic cycle. Ranking within each region, primarily based on anthropogenic stressors of water resources, ensures that settings representing important water-resource challenges for the nation will be studied.

contiguous United States↗

USGS2020

Simulation of the probabilistic plume extent for a potential replacement wastewater-infiltration lagoon, and probabilistic contributing areas for supply wells for the Town of Lac du Flambeau, Vilas County, Wisconsin

An existing two-dimensional, steady-state groundwater-flow model of the shallow groundwater-flow system of the Lac du Flambeau Reservation in Vilas County, Wisconsin, originally developed by the U.S. Geological Survey, was used to simulate the potential for wastewater from a proposed relocation of a wastewater lagoon to contaminate the Lac du Flambeau Band of Lake Superior Chippewa’s drinking-water-supply wells. This simulation was completed by the U.S. Geological Survey in cooperation with the Lac du Flambeau Band of Lake Superior Chippewa and Indian Health Service. The simulated scenarios consisted of removing wastewater infiltration from existing lagoons and re-applying that infiltration at the proposed location. Two analyses were performed for the scenarios. First, the probable extent of the plume discharging from the proposed infiltration lagoons was mapped with a Monte Carlo algorithm that used uncertainty identified during the calibration process to simulate thousands of possible outcomes. Second, the Monte Carlo method was again used to simulate a probabilistic contributing area for the Tribe’s nearby “Main Pumphouse” supply wells. The purpose of the simulations was to evaluate the potential for infiltrated wastewater to be captured by the public-supply wells. Most features of the previously developed model remained unchanged, including calibrated parameters such as hydraulic conductivity and recharge. Thus, the same covariance distributions that were generated during calibration of the regional model (Juckem and others, 2014) remained unchanged and were used to inform the Monte Carlo simulations for the scenario simulations described in this report. The reader is encouraged to read the full report by Juckem and others (available at https://doi.org/10.3133/sir20145020 ) for a detailed description of the model design and calibration, as well as a description of the Monte Carlo method, its limitations, and the original results. Results for these new scenarios indicate that the probabilistic plume extent for the proposed infiltration lagoons does not reach the Main Pumphouse wells using pumping rates and wastewater volumes estimated for 2010. Similarly, the contributing area for the Main Pumphouse wells does not capture water from within the proposed infiltration lagoon footprint. However, at higher pumping rates and wastewater volumes, as projected by the Tribe for about 2035, the contributing area for the Main Pumphouse wells do include particles that originated within the proposed lagoon footprint, albeit at low probabilities. That is, for a few of the thousands of simulations that represented a range of calibration-informed parameter covariances, some amount of infiltrated wastewater was captured by the Main Pumphouse wells under projected 2035 conditions.

Wisconsin↗

USGS2020

Towards reproducible environmental modeling for decision support: A worked example

A fully worked example of decision-support-scale uncertainty quantification (UQ) and parameter estimation (PE) is presented. The analyses are implemented for an existing groundwater flow model of the Edwards aquifer, Texas, USA, and are completed in a script-based workflow that strives to be transparent and reproducible. High-dimensional PE is used to history-match simulated outputs to corresponding state observations of spring flow and groundwater level. Then a hindcast of a historical drought is made. Using available state observations recorded during drought conditions, the combined UQ and PE analyses are shown to yield an ensemble of model results that bracket the observed hydrologic responses. All files and scripts used for the analyses are placed in the public domain to serve as a template for other practitioners who are interested in undertaking these types of analyses.

Texas↗

USGS2020

A non-intrusive approach for efficient stochastic emulation and optimization of model-based nitrate-loading management decision support

Use of physically-motivated numerical models like groundwater flow-and-transport models for probabilistic impact assessments and optimization under uncertainty (OUU) typically incurs such a computational burdensome that these tools cannot be used during decision making. The computational challenges associated with these models can be addressed through emulation. In the land-use/water-quality context, the linear relation between nitrate loading and surface-water/groundwater nitrate concentrations presents an opportunity for employing an efficient model emulator through the application of impulse-response matrices. When paired with first-order second-moment techniques, the emulation strategy gives rise to the “stochastic impulse-response emulator” (SIRE). SIRE is shown to facilitate non-intrusive, near-real time, and risk-based evaluation of nitrate-loading change scenarios, as well as nitrate-loading OUU subject to surface-water/groundwater concentration constraints in high decision variable and parameter dimensions. Two case studies are used to demonstrate SIRE in the nitrate-loading context.

Environmental Modelling and Software↗

USGS2020

Groundwater model simulations of stakeholder-identified scenarios in a high-conflict irrigated area

This study investigated collaborative groundwater‐flow modeling and scenario analysis in the Little Plover River basin, Wisconsin, USA where an unconfined aquifer supplies groundwater for agricultural irrigation, industrial processing, municipal water supply, and stream baseflow. We recruited stakeholders with diverse interests to identify, prioritize, and evaluate scenarios defined as management changes to the landscape. Using a groundwater flow model, we simulated the top 10 stakeholder‐ranked scenarios under historically informed dry, average, and wet weather conditions and evaluated the ability of scenarios to meet government‐defined stream flow performance measures. Results show that multiple changes to the landscape are necessary to maintain optimum stream flow, particularly during dry years. Yet, when landscape changes from three scenarios—transferring water from the local waste water treatment plant to basin headwaters, moving municipal wells further from the river and downstream, and converting 240 acre (97 ha) of irrigated land to unirrigated land—were simulated in combination, the probability of meeting or exceeding optimum flows rose to 75, 65, and 34% at upper, mid, and lower stream gages, respectively, in dry climate conditions. Discussions with stakeholders reveal that the collaborative model and scenario analysis process resulted in social learning that built upon the existing complex and dynamic institutional landscape. The approach provided a forum for solution‐based discussions, and the model served as an important mediation tool for the development and evaluation of community‐defined scenarios in a high conflict environment. Today, stakeholders continue to work collaboratively to overcome challenges and implement voluntary solutions in the basin.

Wisconsin↗

USGS2020

Integrated hydro-terrestrial modeling: Development of a national capability

Water is one of our most important natural resources and is essential to our national economy and security. Multiple federal government agencies have mission elements that address national needs related to water. Each water-related agency champions a unique science and/or operational mission focused on advancing a portion of the nation’s ability to meet our water-related challenges, often in close collaboration with scientists from the academic community. These diverse mission needs have engendered a rich and extensive base of water-related data and modeling capabilities. While useful for their intended purposes, these capabilities are not well integrated to address complex regional problems and overarching national problems. These major investments by several federal agencies and their scientific partners, however, lay the foundation for an integrated hydro-terrestrial modeling and data infrastructure that will enhance knowledge, understanding, prediction, and management of the nation’s diverse water challenges.

Report↗

USGS2019

A simple method for simulating groundwater interactions with fens to forecast development effects

Protection of fens–wetlands dependent on groundwater discharge–requires characterization of groundwater sources and stresses. Because instrumentation and numerical modeling of fens is labor intensive, easy-to-apply methods that model fen distribution and their vulnerability to development are desirable. Here we demonstrate that fen areas can be simulated using existing steady-state MODFLOW models when the unsaturated zone flow (UZF) package is included. In cells where the water table is near land surface, the UZF package calculates a head difference and scaled conductance at these “seepage drain” cells to generate average rates of vertical seepage to the land. This formulation, which represents an alternative to blanketing the MODFLOW domain with drains, requires very little input from the user because unsaturated flow-routing is inactive and results are primarily driven by easily obtained topographic information. Like the drain approach, it has the advantage that the distribution of seepage areas is not predetermined by the modeler, but rather emerges from simulated heads. Beyond the drain approach, it takes account of intracell land surface variation to explicitly quantify multiple surficial flows corresponding to infiltration, rejected recharge, recharge and land-surface seepage. Application of the method to a basin in southeastern Wisconsin demonstrates how it can be used as a decision-support tool to first, reproduce fen distribution and, second, forecast drawdown and reduced seepage at fens in response to shallow pumping.

Wisconsin↗

USGS2019

Cross-scale interactions dictate regional lake carbon flux and productivity response to future climate

Lakes support globally important food webs through algal productivity and contribute significantly to the global carbon cycle. However, predictions of how broad-scale lake carbon flux and productivity may respond to future climate are extremely limited. Here, we used an integrated modeling framework to project changes in lake-specific and regional primary productivity and carbon fluxes under 21st century climate for thousands of lakes. We observed high uncertainty in whether lakes collectively were to increase or decrease lake CO 2 emissions and carbon burial in our modeled region owing to divergence in projected regional water balance among climate models. Variation in projected air temperature influenced projected changes in lake primary productivity (but not CO 2 emissions or carbon burial) as warmer air temperatures decreased productivity through reduced lake water volume. Cross-scale interactions between regional drivers and local characteristics dictated the magnitude and direction of lake-specific carbon flux and productivity responses to future climate.

Geophysical Research Letters↗

USGS2019

Revisiting “An Exercise in Groundwater Model Calibration and Prediction” after 30 years: Insights and New Directions

In 1988, an important publication moved model calibration and forecasting beyond case studies and theoretical analysis. It reported on a somewhat idyllic graduate student modeling exercise where many of the system properties were known; the primary forecasts of interest were heads in pumping wells after a river was modified. The model was calibrated using manual trial‐and‐error approaches where a model's forecast quality was not related to how well it was calibrated. Here, we investigate whether tools widely available today obviate the shortcomings identified 30 years ago. A reconstructed version of the 1988 true model was tested using increasing parameter estimation sophistication. The parameter estimation demonstrated the inverse problem was non‐unique because only head data were available for calibration. When a flux observation was included, current parameter estimation approaches were able to overcome all calibration and forecast issues noted in 1988. The best forecasts were obtained from a highly parameterized model that used pilot points for hydraulic conductivity and was constrained with soft knowledge. Like the 1988 results, however, the best calibrated model did not produce the best forecasts due to parameter overfitting. Finally, a computationally frugal linear uncertainty analysis demonstrated that the single‐zone model was oversimplified, with only half of the forecasts falling within the calculated uncertainty bounds. Uncertainties from the highly parameterized models had all six forecasts within the calculated uncertainty. The current results outperformed those of the 1988 effort, demonstrating the value of quantitative parameter estimation and uncertainty analysis methods.

Groundwater↗

USGS2019

Growing pains of crowdsourced stream stage monitoring using mobile phones: The development of CrowdHydrology

Citizen science-based approaches to monitor the natural environment tend to be bimodal in maturity. Older and established programs such as the Audubon’s Christmas bird count and Community Collaborative Rain, Hail, and Snow Network (CoCoRaHS) have thousands of participants across decades of observations, while less mature citizen science projects have shorter lifespans often focused on local or regional observations with tens or hundreds of participants. For the latter, it can be difficult to transition into a more mature and sustainable citizen science-based research program. This paper focuses on this transition by evaluating CrowdHydrology (ca. 2010), a citizen science project that has transitioned from a regional to national network. It evaluates the data accuracy, citizen participation, and station popularity. The CrowdHydrology network asks citizens to send in text messages of water levels in streams and lakes, which has resulted in 16,294 observations submitted by over 8,000 unique participants at 120 unique locations. Using water level data and participation records from CrowdHydrology, we analyze the expansion and citizen participation from a regional to national citizen science network. We identify barriers to participation and evaluate why some citizen science observation stations are popular while others are not. We explore our chosen contributory program model for CrowdHydrology and the influence this model has had on long-term participation. Results demonstrate a highly variable rate of contributions of citizen scientists. This paper proposes hypotheses on why many of our observations are from one-time participants and why some monitoring stations are more popular than others. Finally, we address the future expansion of the CrowdHydrology network by evaluating successful monitoring locations and growing interest of watershed groups to expand the network of gauges.

Frontiers in Earth Science↗

USGS2019

Use of a Numerical Model to Simulate the Hydrologic System and Transport of Contaminants Near Joint Base Cape Cod, Western Cape Cod, Massachusetts

Historical training and operational activities at Joint Base Cape Cod (JBCC) on western Cape Cod, Massachusetts, have resulted in the release of contaminants into an underlying glacial aquifer that is the sole source of water to the surrounding communities. Remedial systems have been installed to contain and remove contamination from the aquifer. Groundwater withdrawals for public supply are expected to increase as the region continues to urbanize. Increases in water-supply withdrawals and wastewater return flow likely will affect the hydrologic system around JBCC and could affect the transport of any contamination that may remain in the aquifer following remediation of contamination from the JBCC. The U.S. Geological Survey, in cooperation with the Air Force Civil Engineer Center, developed a numerical, steady-state regional model of the Sagamore flow lens on western Cape Cod and evaluated the potential effects of future (2030) groundwater withdrawals on water levels, streamflows, hydraulic gradients, and advective transport near the JBCC. The aquifer consists generally of sandy sediments underlain by impermeable bedrock and is bounded laterally by a freshwater/saltwater interface. Data on the altitude of the bedrock surface, position of the freshwater/saltwater interface, lithology of the aquifer, spatial distribution of recharge, and hydrologic boundaries were incorporated into the three-dimensional, finite-difference groundwater flow model. Some inputs into the numerical model—aquifer properties, leakances, and recharge—are represented as parameters to facilitate estimation of optimal parameter values in an inverse calibration. A hybrid parameterization scheme, with both zones of piecewise constancy and pilot points, is used to represent hydraulic conductivity; other adjustable parameters include recharge, boundary leakance, and porosity. Data on water levels, the distribution of subsurface contamination, and groundwater ages were compiled, evaluated, and used to develop observations of long-term average hydraulic gradients and advective-transport patterns. These observations of steady-state hydrologic conditions were combined with the parameterized groundwater model in an inverse calibration to estimate model parameters that best fit the observations. Current (2010) and future (2030) conditions were simulated in the calibrated model to characterize the groundwater flow system and to determine potential effects of increased groundwater withdrawals on advective-transport patterns at the JBCC. Groundwater flow and advective transport are radially outward from a water-table divide in the northern part of the JBCC; flow diverges from the divide toward all points of the compass. Most groundwater flow and contaminant transport occur in shallow parts of the aquifer. On average, about one-half of the groundwater flux occurs in the shallowest 20 percent of the saturated thickness; shallow flow is even more predominant near streams and lakes. Projected (2030) increases in groundwater withdrawals decrease water levels by a maximum of about 1.2 feet in the northern part of the JBCC; drawdowns exceeding 1 foot generally are limited to areas near the largest increases in withdrawals, such as in the northern part of the JBCC, near Long Pond in Falmouth, and in eastern Barnstable. Streamflow decreases average about 6 percent; the largest decreases are in areas with the largest drawdowns. Changes in hydraulic-gradient directions at the water table exceed 1 degree in about 13 percent of the aquifer, generally near groundwater divides where gradient magnitudes are small and near large groundwater withdrawals. Predictions of advective transport from randomly selected locations at the water table are similar for current (2010) and future (2030) groundwater withdrawals. The results indicate that projected increases in groundwater withdrawals affect water levels and streamflows, but effects on hydraulic gradients and advective transport at the JBCC likely are small. Several underlying assumptions inherent in the model, including observations and weights used in the calibration, representation of local-scale heterogeneity, and simulation of the freshwater/saltwater interface, could affect model calibration and predictions; these assumptions were evaluated with alternative models and alternative inverse calibrations. Eight alternative calibrations were performed in which different, but reasonable, observations and weights were used. The preferred calibrated model had the best overall fit to the observations. Fine-grained silty sediments occur in many parts of the aquifer, and silt lenses can locally affect hydraulic gradients. A set of alternative models in which silts were represented with different correlation distances and hydraulic conductivities indicated that explicitly representing silt lenses could affect model calibration but that the implicit representation of local-scale heterogeneity may be sufficient at the regional scale to represent regional-scale hydraulic gradients. For the coastal boundary, two alternative models representing silty and sandy seabeds and their associated interface positions were developed to test the importance of the assumed coastal-boundary condition. The two alternative models resulted in different predictions of streamflow—streamflows increase with smaller (silty) seabed leakances. However, predictions of advective transport, particularly near the JBCC, generally were similar between the alternative and preferred calibrated models, indicating that the seabed leakance and associated interface position at the coastal boundary does not affect simulations of advective transport in inland parts of the aquifer.

Massachusetts↗

USGS2019

Groundwater flow model for Western Chippewa County–Including analysis of water resources related to industrial sand mining and irrigated agriculture

A groundwater flow model for western Chippewa County, Wisconsin, was developed by the Wisconsin Geological and Natural History Survey (WGNHS) and the U.S. Geological Survey (USGS) using the computer program MODFLOW. The model is the result of a five-year groundwater study commissioned by Chippewa County in 2012 to evaluate the effects of industrial sand mining and irrigated agriculture on the county’s water resources. The study incorporates existing data and newly acquired data from fieldwork conducted within the study area. The groundwater model may be useful for future investigations, such as evaluation of proposed high-capacity well sites, development of municipal wellhead protection plans, and studies that seek to further quantify surface water-groundwater relationships. The model conceptualizes the hydrostratigraphy of western Chippewa County as six stacked layers. Each layer is distinct, beginning with unlithified glacial material at the surface, and alternating between sandstones (that act as aquifers) and shale units (that serve as aquitards). The model is bounded below by Precambrian crystalline bedrock and its perimeter was derived from a regional-scale groundwater flow model. The MODFLOW model represented average conditions during 2011–2013 with “steady-state” assumptions, meaning that simulated water levels do not fluctuate seasonally or from year to year. Steady-state models simplify natural variability, making results of scenario simulations easier to interpret and compare while also maximizing effects of stressors because the simulated stress is always applied (not halted after a few months or years). Model calibration used the parameter estimation code (PEST), and calibration targets included heads (groundwater levels) and streamflows. Calibration focused on 2011–2013 because a large amount of head and streamflow data were available for that period. The MODFLOW model explicitly simulates all sources and sinks of water, including groundwater/surface-water interaction with streamflow routing. Model input included estimates of aquifer hydraulic conductivity and a spatial groundwater recharge distribution developed using a GIS-based soil-water-balance (SWB) model applied to the model area. Groundwater withdrawals were simulated for 269 high-capacity wells across the entire model domain, which includes western Chippewa County and adjacent portions of Dunn, Barron, and Rusk Counties. Collectively, these wells withdrew about 1.14 million gallons per year between 2011 and 2013. Once the model was calibrated, it was applied to two distinct scenarios of increased groundwater withdrawals: one evaluating hydrologic effects of more intensive industrial sand mining and the second evaluating the hydrologic effects of more intensive agricultural irrigation practices. Each scenario was developed with input from Chippewa County and a stakeholder group established expressly for this study. The scenarios were designed to represent reasonable future buildout conditions for both mining and irrigated agriculture. The mining scenario underscores the potential hydrologic effects related to changing land-use practices (i.e., hilltops and farmland becoming sand mines), while the irrigated agriculture scenario illustrates the potential hydrologic effects of intensifying existing land-use practices (i.e., installing new wells to irrigate farm fields). While each scenario evaluated distinctly different conditions, modeling results demonstrated the potential of both scenarios to lower the water table and reduce baseflows in headwater streams within the modeled area. In the case of irrigated agriculture, hydrologic effects were associated directly with groundwater withdrawals. By assuming that irrigation did not decrease, this steady-state simulation represented a sustained future effect. By contrast, hydrologic effects of industrial sand mining were the result of both groundwater withdrawals at mines and land-use changes that effectively reduced recharge to groundwater over distinct phases of active mining. This scenario included a post-mining phase, during which groundwater withdrawals stopped and mined areas were reclaimed to undeveloped prairie grass cover. If reclamation to undeveloped prairie indeed occurs as simulated, long-term increases in the water table and stream baseflows are possible. In this sense, the scenario representing build out of irrigated agriculture led to long-term baseflow declines while the future buildout of industrial sand mining led to declines that dissipated following mine reclamation to undisturbed prairie. Future investigations in similar hydrogeologic settings may find the following insights gleaned from this study useful: ❚❚ The characterization of hydrogeologic properties, delineation of hydrogeologic units, and calibration of groundwater flow models benefited from incorporation of accurate well construction reports, high-quality borehole geophysical logs, and streamflow gaging data. ❚❚ Infiltration testing performed in active mining areas provided evidence that reducing the degree and extent of compaction and enhancing areas designed to retain and infiltrate stormwater runoff could potentially reduce runoff and increase groundwater recharge. ❚❚ Similarly, reclaiming mined areas to prairie grasses would be expected to reduce runoff and increase groundwater recharge by reducing compaction and improving soil structure and vegetation that can slow runoff and enhance infiltration.

Wisconsin↗

USGS2018

Wrangling distributed computing for high-throughput environmental science: An introduction to HTCondor

Biologists and environmental scientists now routinely solve computational problems that were unimaginable a generation ago. Examples include processing geospatial data, analyzing -omics data, and running large-scale simulations. Conventional desktop computing cannot handle these tasks when they are large, and high-performance computing is not always available nor the most appropriate solution for all computationally intense problems. High-throughput computing (HTC) is one method for handling computationally intense research. In contrast to high-performance computing, which uses a single "supercomputer," HTC can distribute tasks over many computers (e.g., idle desktop computers, dedicated servers, or cloud-based resources). HTC facilities exist at many academic and government institutes and are relatively easy to create from commodity hardware. Additionally, consortia such as Open Science Grid facilitate HTC, and commercial entities sell cloud-based solutions for researchers who lack HTC at their institution. We provide an introduction to HTC for biologists and environmental scientists. Our examples from biology and the environmental sciences use HTCondor, an open source HTC system.

PLOS Computational Biology↗