Geology Reports⌕ Search

Geology topics

Michael N. Fienen

Publications and source records attributed to Michael N. Fienen.

At least 37 records · Page 2Linked to original sources
USGS2024

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education↗
USGS2023

Effects of auto-adaptive localization on a model calibration using ensemble methods

Simulations of the natural systems for environmental decision-making typically benefit from a highly parameterized approach (Hunt et al. 2007; Doherty and Hunt 2010), which enhances the flow of information contained in state observations to the parameters and improves application to decision support. However, parameter estimation (PE) with highly parameterized environmental models using traditional approaches (e.g., Doherty and Hunt 2010) is computationally intensive. Attempts at addressing the computational burden include improved computing approaches (e.g., Schreüder 2009; Hunt et al. 2010) and advances in algorithmic approaches (e.g., Tonkin and Doherty 2005; Welter et al. 2012, 2015). Recently, the iterative ensemble smoother (IES) approach (Chen and Oliver 2013; White 2018; White et al. 2020a) has greatly improved the efficiency of the PE calibration process compared to previous algorithms while concurrently providing nonlinear estimates of uncertainty (Hunt et al. 2021).

Groundwater↗
USGS2022

Modflow-setup: Robust automation of groundwater model construction

In an age of both big data and increasing strain on water resources, sound management decisions often rely on numerical models. Numerical models provide a physics-based framework for assimilating and making sense of information that by itself only provides a limited description of the hydrologic system. Often, numerical models are the best option for quantifying even intuitively obvious connections between human activities and water resource impacts. However, despite many recent advances in model data assimilation and uncertainty quantification, the process of constructing numerical models remains laborious, expensive, and opaque, often precluding their use in decision making. Modflow-setup aims to provide rapid and consistent construction of MODFLOW groundwater models through robust and repeatable automation. Common model construction tasks are distilled in an open-source, online code base that is tested and extensible through collaborative version control. Input to Modflow-setup consists of a single configuration file that summarizes the workflow for building a model, including source data, construction options, and output packages. Source data providing model structure and parameter information including shapefiles, rasters, NetCDF files, tables, and other (geolocated) sources to MODFLOW models are read in and mapped to the model discretization, using Flopy and other general open-source scientific Python libraries. In a few minutes, an external array-based MODFLOW model amenable to parameter estimation and uncertainty quantification is produced. This paper describes the core functionality of Modflow-setup, including a worked example of a MODFLOW 6 model for evaluating pumping impacts to a lake in central Wisconsin, United States.

Wisconsin↗
USGS2022

Simulation of regional groundwater flow and groundwater/lake interactions in the Central Sands, Wisconsin

A multiscale, multiprocess modeling approach was applied to the Wisconsin Central Sands region in central Wisconsin to quantify the connections between the groundwater system, land use, and lake levels in three seepage lakes in Waushara County, Wisconsin: Long and Plainfield (The Plainfield Tunnel Channel Lakes), and Pleasant Lakes. A regional groundwater-flow model, the Newton Raphson formulation of the U.S. Geological Survey modular finite-difference flow model groundwater-modeling package (MODFLOW-NWT), centered on the lakes, was used to extend regional surface-water boundaries to provide boundary conditions for two focused inset models, in the hydrologic simulation modeling package (MODFLOW 6), at higher resolution around the lakes. Land use and groundwater use were simulated at a regional scale using the Soil Water Balance model, which provided recharge and water-use boundary conditions for the MODFLOW models. Agricultural irrigation is the primary groundwater use in the area. Land and groundwater-use scenarios representing no irrigation, current (2018) irrigation, and potential future irrigation were simulated with the groundwater-flow model and the lake levels over a 38-year representative climate period.

Wisconsin↗
USGS2022

Areas contributing recharge to priority wells in valley-fill aquifers in the Neversink River and Rondout Creek drainage basins, New York

In southeastern New York, the villages of Ellenville, Wurtsboro, Woodridge, the hamlet of Mountain Dale, and surrounding communities in the Neversink River and Rondout Creek drainage basins rely on wells that pump groundwater from valley-fill glacial aquifers for public water supply. Glacial aquifers are vulnerable to contamination because they are highly permeable and have a shallow depth to water table. To protect the quality of these water resources, water managers need accurate information about the areas that contribute recharge to production wells that pump from these aquifers. The New York State Department of Environmental Conservation and the New York State Department of Health designated eight priority wells in this region for which water supply protection is of primary concern. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 with the general objectives of (1) improving understanding of regional groundwater-flow system, (2) delineating areas contributing recharge to eight priority production wells, and (3) quantifying the uncertainty of these contributing areas in a probabilistic way that can be used to inform decision-making related to priority well source-water protection. To complete these objectives, a MODFLOW 6 groundwater model was created encompassing the eight priority wells and the surrounding flow system, which includes parts of the Neversink River and Rondout Creek Basins in Sullivan County and Ulster County, New York. The model was built using Python tools (such as flopy , modflow-setup , and sfrmaker ) that facilitate transparent and repeatable model development using existing datasets. The model parameters were estimated with a stepwise approach using an iterative ensemble smoother implementation of the Parameter ESTimation software PEST++ (version 5.0.0). We evaluated initial “best guess” parameter bounds with a prior Monte Carlo analysis. Results of the first prior Monte Carlo analysis were used to make informed adjustments to model parameter bounds (typically resulting in expanded bounds), and a second prior Monte Carlo analysis was run to identify improved ranges for model parameters during history matching. The history matching effort produced an ensemble of parameter values for the groundwater-flow model that spans the range of values within prior uncertainty bounds. The ensemble is informed by the historical observation data, within a reasonable range of uncertainty on those observations. This history-matched ensemble was used in a particle tracking Monte Carlo analysis to delineate the areas contributing recharge to priority wells. The groundwater-flow and particle tracking (MODPATH7) models were run once for each ensemble member. Deterministic contributing areas computed for each ensemble member were aggregated to produce maps showing the probability that a location contributes recharge to priority wells. Finally, the particle tracking Monte Carlo analysis was repeated for six pumping scenarios, representing a wide range of possible pumping levels, to incorporate uncertainty in future pumping rates related to population growth or other management decisions. Increasing pumping rates generally led to larger contributing recharge areas and larger areas of high probability that a location contributes recharge to priority wells. These maps show the overall uncertainty of the areas contributing recharge to priority wells in the study area and provide a tool for risk-based decision making for protection of well source water.

New York↗
USGS2022

A scalable model-independent iterative data assimilation tool for sequential and batch estimation of high dimensional model parameters and states

Ensemble-based data assimilation (DA) methods have displayed strong potential to improve model state and parameter estimation across several disciplines due to their computational efficiency, scalability, and ability to estimate uncertainty in the dynamic states and the parameters. However, a barrier to adoption of ensemble DA methods remains. Namely, there is currently a lack of available tools that enable efficient and scalable DA in a non-intrusive fashion and that support implementation flexibility. This paper presents an open-source software tool (PESTPP-DA) that implements a range of data assimilation methods—Ensemble Kalman filter, Ensemble Kalman Smoother and Ensemble Smoother—using the widely known PEST model-interface protocols, to interact with any model. Two iterative solutions can be used for nonlinear and/or non-Gaussian assimilation problems. To demonstrate the broad range of PESTPP-DA applications, two synthetic case studies are presented: (1) the Lorenz model and (2) a groundwater pumping test in the presence of a non-Gaussian hydraulic conductivity field.

Environmental Modelling & Software↗
USGS2022

A model-independent tool for evolutionary constrained multi-objective optimization under uncertainty

An open-source tool has been developed to facilitate constrained single- and multi-objective optimization under uncertainty (CMOU) analyses. The tool uses the well-known PEST interface protocols to communicate with the underlying forward simulation, making it non-intrusive. The tool contains a built-in parallel run manager to make use of heterogeneous and distributed computing resources. Several popular and well-known evolutionary algorithms are implemented and can be combined with a range of approaches to represent uncertainty in model-derived constraint/objective values. These attributes serve to address the current barrier to adopt advanced CMOU analyses for a wide range of decision-support problems across the environmental modeling spectrum. We demonstrate the capabilities of the CMOU tool on a well-known analytical benchmark problem that we augmented to include uncertainty, as well as on a synthetic density-dependent coastal groundwater management benchmark problem. Both demonstrations highlight the importance of explicitly accounting for uncertainty to convey risk and reliability in pareto-optimal design.

Environmental Modelling & Software↗
USGS2022

Assessing spatial transferability of a random forest metamodel for predicting drainage fraction

Fully distributed hydrological models are widely used in groundwater management, but model speed and data requirements impede their use for decision support purposes. Metamodels provide a simpler and faster model which emulates the underlying complex model using machine learning techniques. However, metamodel predictions beyond the ranges, in space and/or time, of training data are highly uncertain, and thus it is important to assess the predictive model performance to ranges outside the training data, i.e., model transferability . We present a novel methodology for evaluating model transferability to areas not contained in the training data set, based on various metrics that quantify the differences in covariate distributions between training and testing data. The transferability method can be employed as a screening tool to assess the suitability of a metamodel for spatial prediction beyond its training domain. We evaluated this transferability approach on a Random Forest metamodel of a 1000 km 2 fully distributed coupled groundwater model for predicting drainage fraction, the partitioning of infiltrating water between drains and groundwater. We conducted spatial cross-validation on 9 holdout sub-basins to assess metamodel transferability beyond sampling locations and compared this estimate with a random split-sample validation test. Using mappable covariates only, the metamodel showed high performance (R 2 = 0.79) tested on a 20% randomly sampled holdout. Conversely, metamodel performance significantly decreased for the 9 spatial holdouts (R 2 ranging from 0.13 to 0.61). We document that the proposed transferability metric correlates with metamodel predictive performance, and demonstrate its use to assess model transferability to datasets outside the training data spatial domain.

Journal of Hydrology↗
USGS2021

Groundwater/surface-water interactions in the Partridge River Basin and evaluation of hypothetical future mine pits, Minnesota

The Partridge River Basin (PRB) covers 156 square miles in northeastern Minnesota with headwaters in the Mesabi Iron Range. The basin is characterized by extensive wetlands, lakes, and streams in poorly drained and often thin glacial material overlying Proterozoic bedrock. To better understand the interaction between these extensive surface water features and the groundwater system, a three-dimensional, steady-state, groundwater-flow model of the PRB was developed by the U.S. Geological Survey in cooperation with the Great Lakes Indian Fish & Wildlife Commission using the finite-difference computer code MODFLOW-NWT. The model simulates steady-state base flow in streams and groundwater interactions using the streamflow routing (SFR2) package. Existing mining features including tailings basins, stockpiles, pumped mine pits, and flooded mine pits were simulated using either high hydraulic conductivity zones or the drain (DRN) package. The unsaturated zone flow (UZF) package was used to better represent the groundwater system in areas with a high water table and for wetlands often associated with such areas. UZF typically is used to represent unsaturated zone processes but also can simulate the rejection of recharge and groundwater discharge to the land surface when the water table is near land surface. The steady-state model used data from the 2011 to 2013 period when 2011 high-resolution land surface (light detecting and ranging [lidar]) data were available that reflected land-surface and water elevations from mining activity in the basin. The parameter-estimation software suite PEST_HP was used to obtain a best fit of the modeled to measured groundwater levels, streamflow, pit inflow rates, and mapped peat deposits. The PEST calibration used the target residuals from two models with the same model parameters and targets from two separate periods: (1) a 1995–2015 calibration model, which provided a larger number of calibration targets, and (2) a 2011–2013 mining conditions model, which included calibration targets that reflected conditions consistent with the modeled mine-workings topography. Calibration of the PRB model resulted in ranges of glacial horizontal hydraulic conductivity parameters that generally agreed with literature values and other models of the region. Horizontal hydraulic conductivity of the bedrock was higher in the upper bedrock layers where numerous and continuous fractures have been observed and lower in the deeper bedrock layers. Average basin-wide calibrated infiltration was 5.3 inches per year. An average of 4.6 inches per year of infiltration crosses the water table and becomes recharge and 0.7 inch per year is rejected by UZF due to saturated conditions at the land surface. Simulated groundwater runoff (the sum of rejected recharge and groundwater seepage to the land surface) can either be routed to streams or removed from the model as evapotranspiration. The calibrated model indicates relatively shallow groundwater-flow paths dominating and approximately 50 percent of the stream base flow coming from groundwater runoff. The 2011–2013 mining conditions model was then used to develop five model scenarios simulating the response of the groundwater and surface-water system to potential hydrologic stress. The purpose of these mine pit scenarios is to present a possible workflow to quantify a model’s uncertainty for a given model forecast and serve as a possible guide for initial data collection that may improve a future model’s ability to make such a forecast. The scenarios included one scenario with the currently existing Peter Mitchell pit at final buildout and flooded to an elevation of 1,500 feet, and four scenarios with a hypothetical, new mine pit plus the flooded Peter Mitchell at final buildout. The five model scenarios were used to forecast streamflow at six locations in the PRB, pit inflow rates for the new mine pits and the flooded Peter Mitchell pit, and the average depth to water in 12 wetlands. A linear uncertainty analysis was performed using information from the PEST calibration and tools in the PyEMU python package to assess model uncertainty propagation to the model forecasts. Streamflows generally were reduced with future mining and the greatest streamflow reductions occurred from the flooded Peter Mitchell Pit, probably due to its large size. Average depth to groundwater in wetlands was most affected the closer the wetland was to a new mine pit. Linear uncertainty methods were also used to evaluate data worth, which is the ability for potential new groundwater elevation observations to reduce the uncertainty in scenario forecasts. Data worth was performed for a grid of new hydraulic head observations. Overall, areas with nonnegligible data worth generally corresponded to wetland areas with no groundwater seepage to land surface from UZF. These model behaviors indicated that the land-surface boundary condition simulated by the UZF package was pinning the groundwater elevations to the land surface in areas with groundwater seepage (33 percent of the 2011–2013 base conditions model) such that the sensitivity to new observations in these areas was minimal. Therefore, representing wetlands as boundary conditions minimized the usefulness of data worth calculations because wetland areas were present over a large part of the model domain. Probabilistic capture zones were estimated for each of the mines in the model scenarios. A capture zone represents the area contributing recharge to a model feature, like a well or a mine pit, and can be calculated by forward tracking particles from the water table. By using Monte Carlo techniques, it is possible to generate estimated capture zones that include the probability of recharge capture given the uncertainty present in the model. Monte Carlo techniques use randomly generated model parameter sets sampled from a plausible parameter range to create many possible realizations. The resulting capture zone arrays were calculated by tallying the total number of realizations in which a particle from a model cell was captured by the feature. Probabilities from the Monte Carlo runs ranged from 1 (captured in 100 percent of the runs) near the pits to 0 (captured in 0 percent of the runs) at the edges of the capture zone. Capture zones were not always spatially continuous; for example, the capture zone for the proposed mine pits south of the flooded Peter Mitchell pit was discontinuous with capture surrounding the proposed mine pit and north of the flooded Peter Mitchell pit. This northern section represents deeper groundwater flow paths that originate in the topographic high, move under the flooded pit, and discharge into the proposed pit. This pattern of capture indicates the possibility of some deeper flow through the upper fractured bedrock when the shallow groundwater flow system is modified. These results underscore that future site-specific applications of the base condition model require the input of site-specific data and recalibration to focus on the site of interest.

Minnesota↗
USGS2021

Risk-based wellhead protection decision support: A repeatable workflow approach

Environmental water management often benefits from a risk-based approach where information on the area of interest is characterized, assembled, and incorporated into a decision model considering uncertainty. This includes prior information from literature, field measurements, professional interpretation, and data assimilation resulting in a decision tool with a posterior uncertainty assessment accounting for prior understanding and what is learned through model development and data assimilation. Model construction and data assimilation are time consuming and prone to errors, which motivates a repeatable workflow where revisions resulting from new interpretations or discovery of errors can be addressed and the analyses repeated efficiently and rigorously. In this work, motivated by the real world application of delineating risk-based (probabilistic) sources of water to supply wells in a humid temperate climate, a scripted workflow was generated for groundwater model construction, data assimilation, particle-tracking and post-processing. The workflow leverages existing datasets describing hydrogeology, hydrography, water use, recharge, and lateral boundaries. These specific data are available in the United States but the tools can be applied to similar datasets worldwide. The workflow builds the model, performs ensemble-based history matching, and uses a posterior Monte Carlo approach to provide probabilistic capture zones describing source water to wells in a risk-based framework. The water managers can then select areas of varying levels of protection based on their tolerance for risk of potential wrongness of the underlying models. All the tools in this workflow are open-source and free, which facilitates testing of this repeatable and transparent approach to other environmental problems.

Groundwater↗
USGS2021

SFRmaker and Linesink-Maker: Rapid construction of streamflow routing networks from hydrography data

Groundwater models have evolved to encompass more aspects of the water cycle, but the incorporation of realistic boundary conditions representing surface water remains time-consuming and error-prone. We present two Python packages that robustly automate this process using readily available hydrography data as the primary input. SFRmaker creates input for the MODFLOW SFR package, while Linesink-maker creates linesink string input for the GFLOW analytic element program. These programs can reduce weeks or even months of manual effort to a few minutes of execution time, and carry the added advantages of reduced potential for error, improved reproducibility and facilitation of step-wise modeling through reduced dependency on a particular conceptual model or discretization. Two real-world examples at the county to multi-state scales are presented.

Groundwater↗
USGS2021

Assessment of groundwater trends near Crex Meadows, Wisconsin

Crex Meadows Wildlife Area (Crex) is a 30,000-acre property in Burnett County, Wisconsin. Crex is managed by the Wisconsin Department of Natural Resources (WDNR) with the goal of providing public recreation opportunities while also protecting the quality of native ecological communities and species on the property. The WDNR’s management strategy includes controlling water levels at flowages in Crex using a system of dikes, water control structures, ditches, and a diversion pump. For the past several decades there has been concern among nearby landowners that the water manage-ment strategy at Crex may be contributing to groundwater flooding in adjacent, privately held properties. This issue has been particularly contentious during periods when regional groundwater elevations are already high. This study was conducted in response to those concerns. For the study, a network of 12 monitoring wells was installed in and to the west of Crex. Groundwater elevations were recorded in the wells before, during, and after water-level changes in the western Crex flowages to assess if groundwater elevations to the west of Crex are detectably affected by the flowage water levels. This study successfully collected groundwater elevations in 11 study wells during a 3-month period in 2019 when water elevations in the Dike 6 flowage and Erickson flowage were lowered and then raised. The data logger at a 12th location failed and no data were recorded. The groundwater elevation trends in these study wells were compared with groundwater elevation trends at a regional U.S. Geological Survey well to provide information for determining if changing the flowage elevations had a noticeable response in the study wells west of Crex Meadows. This analysis was done by (1) evaluating study well groundwater elevation trends compared to the regional well, (2) using a scatter plot of study well and regional well data during raising and lowering periods, (3) assessing horizontal hydraulic gradient data during the study period, and (4) assessing the cumulative departure from the mean groundwater elevation for each well. Overall, regional groundwater elevations had a down-ward trend before and during the flowage lowering period and then had an upward trend during the flowage raising period. This pattern was observed in the regional well and in all the study wells adjacent to and several miles from the flowages. The similarity in patterns indicates that precipitation and regional groundwater flow conditions were the dominant drivers of the system during the study period. The scatter plot and cumulative departure from the mean analysis showed that in addition to regional trends, wells 1, 6, and 7 were likely affected by the changes in the flowage water levels. Overall, at least on the timescale of this study, water management at Crex likely did not have detectable effects on wells outside the Crex property. Wells installed on the Crex property including the wells in the lakebeds of the flowages (wells 1 and 7) and possibly well 6 east of the flowages showed what seems to be minor affects due to water management at Crex.

Wisconsin↗
USGS2021

Extending the capture map concept to estimate discrete and risk-based streamflow depletion potential

A popular and contemporary use of numerical groundwater models is to estimate the discrete relation between groundwater extraction and surface-water/groundwater exchange. Previously, the concept of a “capture map” has been put forward as a means to effectively summarize this relation for decision-making consumption. While capture maps have enjoyed success in the environmental simulation industry, they are deterministic, ignoring uncertainty in the underlying model. Furthermore, capture maps are not typically calculated in a manner that facilitates analysis of varying combinations of extraction locations and/or reaches. That is, they are typically constructed with focus on a single reach or group of reaches. The former of these limitations is important for conveying risk to decision makers and stakeholders, while the latter is important for decision-making support related to surface-water management, where future foci may include reaches that were not the focus of the original capture analysis. Herein, we use the concept of a response matrix to generalize the theory of the capture-map approach to estimate spatially discrete streamflow depletion potential. We also use first-order, second-moment uncertainty estimation techniques with the concept of “risk shifting” to place capture maps and streamflow depletion potential in a stochastic, risk-based framework. Our approach is demonstrated for an integrated groundwater/surface-water model of the lower San Antonio River, Texas, USA.

Groundwater↗
USGS2021

Multi‐constrained catchment scale optimization of groundwater abstraction using linear programming

Due to increasing water demands globally, freshwater ecosystems are under constant pressure. Groundwater resources, as the main source of accessible freshwater, are crucially important for irrigation worldwide. Over‐abstraction of groundwater leads to declines in groundwater levels; consequently, the groundwater inflow to streams decreases. The reduction in base flow and alteration of the stream flow regime can potentially have an adverse impact on groundwater‐dependent ecosystems. A spatially distributed, coupled groundwater‐surface water model can simulate the impacts of groundwater abstraction on aquatic ecosystems. A constrained optimization algorithm and a simulation model in combination can provide an objective tool for the water practitioner to evaluate the interplay between economic benefits of groundwater abstractions and requirements to environmental flow. In this study, a holistic catchment‐scale groundwater abstraction optimization framework has been developed that allows for a spatially explicit optimization of groundwater abstraction, while fulfilling a pre‐defined maximum allowed reduction of stream flow (base flow (Q95) or median flow (Q50)) as constraint criteria for 1484 stream locations across the catchment. A balanced K‐Means clustering method was implemented to reduce the computational burden of the optimization. The model parameters and observation uncertainties calculated based on Bayesian linear theory allow for a risk assessment on the optimized groundwater abstraction values. The results from different optimization scenarios indicated that using the linear programming optimization algorithm in conjunction with integrated models provides valuable information for guiding the water practitioners in designing an effective groundwater abstraction plan with the consideration of environmental flow criteria important for the ecological status of the entire system.

Groundwater↗