Search USGS⌕ Search

Geology topics

Nicholas Corson-Dosch

Publications and source records attributed to Nicholas Corson-Dosch.

9 recordsLinked to original sources

Multi-objective optimization of a hydro-economic model in an over-allocated agricultural basin

Groundwater depletion for agricultural irrigation poses significant environmental and economic challenges. This study introduces a proof-of-concept that combines hydro-economic modeling, scenario-based modeling, and multi-objective optimization to manage pumping curtailment in an over-allocated basin in the western United States. Three optimization scenarios were evaluated, each offering different degrees of management flexibility. Results reveal that scenarios with finer spatial resolution achieved greater environmental benefits per unit profit loss. Additionally, strategies allowing fractional reductions in curtailed wells–rather than complete shutdowns based on water rights seniority–substantially improved efficiency, highlighting the value of increased decision-making flexibility. Although scenario testing can aid stakeholder engagement and strategy exploration, multi-objective optimization provides a systematic framework to quantify tradeoffs between competing objectives. This combined approach demonstrates promise for building consensus and supporting the design of sustainable water management strategies that balance agricultural livelihoods with ecosystem preservation.

Oregon↗

Assessing flood water infiltration and storage in a restored floodplain

In urban areas, floodplain restoration is gaining prominence as a strategy for restoring the natural functions of floodplain ecosystems and reducing flood risk. This has spurred research into potential interactions between floodwaters, the hyporheic zone, and the floodplain aquifer. An urban restored stream in Wisconsin, USA, was used as a case study to examine four methods to estimate floodplain infiltration and storage during overbank floods. We characterised flood-related infiltration over a 4-year period from 2018 through 2021 by simultaneously and continuously measuring groundwater levels and vertical temperature profiles with stream water levels linked to high-resolution flood inundation maps. High-resolution topographic data helped to quantify surface floodplain storage and the unsaturated soil volume relative to flood stage. Infiltration estimates from the simple methods align well with those from the more complex methods; however, the complex methods provide additional insights about the factors influencing infiltration. Results from all methods indicate that the volume of water that vertically infiltrates during floods is likely small relative to the total volume of the flood, with 0.08%–0.52% of flood water infiltrating into the floodplain, on average. Spatially variable vertical hydraulic gradients, driven by flood depth, groundwater level, and permeability, imply heterogeneous patterns of infiltration across the floodplain. Gradients favourable for infiltration typically occurred during the onset of flooding but, over the study period, were mostly (98% of the time) favourable for groundwater discharge to the channel (non-flood periods). These findings highlight the importance of considering surface-groundwater dynamics, floodplain soils, and unsaturated floodplain volume in defining the benefits of floodplain infiltration for flood attenuation.

Wisconsin↗

Analysis of factors affecting plume remediation in a sole-source aquifer system, southeastern Nassau County, New York

Several plumes of dissolved, chlorinated solvents, including trichloroethylene, have been identified in a sole-source aquifer near the former Northrop Grumman Bethpage Facility and Naval Weapons Industrial Reserve Plant sites in southeastern Nassau County, New York. Past investigations have documented that the groundwater contamination originated from this industrial area and now extends to the south, in the direction of groundwater flow. The intermixed plumes are commonly referred to as the “Navy Grumman groundwater plume.” Detailed groundwater-flow modeling was needed for the New York State Department of Environmental Conservation (NYSDEC) to evaluate design options necessary for the construction, operation, optimization, maintenance, and monitoring of a groundwater extraction and treatment cleanup plan selected in a December 2019 Amended Record of Decision by the NYSDEC to comprehensively address these plumes. Consequently, the NYSDEC began a cooperative study with the U.S. Geological Survey in 2020 to better understand the local hydrogeologic framework using two independent approaches to characterize aquifer heterogeneity and update an existing regional groundwater-flow model to provide transient boundary conditions for new inset groundwater-flow models of the plume area. We developed these detailed inset models for the two independent aquifer characterizations using history-matching techniques coupled with a novel approach to risk-based management optimization of the remedial design. We also used the updated regional model to assess this optimized groundwater extraction and treatment design for potential saltwater intrusion. The ensembles of parameters resulting from history matching provided a platform with which to evaluate capture by water-supply and remedial wells using particle-tracking techniques. Using the ensemble to select a risk stance, we performed multiobjective optimization to identify various configurations of remedial pumping that are consistent with external constraints and that favor potentially competing objectives. Multiple solutions provide tradeoffs that NYSDEC can consider. In general, pumping redistribution may help to prevent further contamination migration downgradient. These and other study results are intended to support decisions for the remedial design focused on the local area encompassing the full extent of the Navy Grumman groundwater plume.

New York↗

Comparing single and multiple objective constrained optimization algorithms for tuning a groundwater remediation system

Groundwater flow and particle tracking models are critical tools to simulate the natural system, contaminant fate and transport, and effects of remediation. Constrained optimization uses models to systematically explore the interplay between remedial design and contaminant fate, considering uncertainty. Sequential Linear Programming (SLP) provides a design alternative addressing a single goal (e.g. maximum hydraulic containment, maximum mass removal). Multi-objective algorithms like Nondominated Sorting Genetic Algorithm (NSGA-II) explore the tradeoffs among such objectives and more (e.g. cost, public-supply well contamination). We explore both approaches at a contaminated site in Long Island, New York USA. We compare the algorithms and ramifications on results. NSGA-II explores, at additional computational cost, explicit tradeoffs among multiple objectives, providing additional insights relative to SLP. The NGSA-II algorithm allows for graphical consideration of three objectives. SLP decision variables often settle at predetermined bounds. Bounds assignment thus differs from parameter estimation; bounds must be acceptable rather than safeguards.

New York↗

Design and calibration of a nitrate decision support tool for groundwater wells in Wisconsin, USA

This paper describes development of a nitrate decision support tool for groundwater wells (GW-NDST) that combines nitrate leaching and groundwater lag-times to compute well concentrations. The GW-NDST uses output from support models that simulate leached nitrate, groundwater age distributions, and nitrate reduction rates. The support models are linked through convolution to simulate nitrate transport to wells. Spatially distributed parameters were adjusted through calibration to 34,255 nitrate sample targets. Prediction uncertainty is illustrated via Monte Carlo realizations informed during calibration. Over 78% of target concentrations were within the simulated range of results from 450 realizations. An example forecasting scenario illustrates that a range of feasible outcomes exist and should be considered when interpreting forecasts for decision making. Uncertainty in forecasting is unavoidable; the intent of characterizing uncertainty in the GW-NDST is to facilitate decision making by increasing insight into the response of nitrate contamination to physical and chemical processes.

Wisconsin↗

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↗

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↗

Evaluating the dynamics of groundwater, lakebed transport, nutrient inflow and algal blooms in Upper Klamath Lake, Oregon, USA

Transport of nutrients to lakes can occur via surface-water inflow, atmospheric deposition, groundwater (GW) inflow and benthic processes. Identifying and quantifying within-lake nutrient sources and recycling processes is challenging. Prior studies in hypereutrophic Upper Klamath Lake, Oregon, USA, indicated that ~60% of the early summer phosphorus (P) load to the lake was internal and hypothesized to be lakebed sediment release. Dynamic nutrient transport processes were examined to better characterize the nutrient sources. One-dimensional heat transport models calibrated to observed lakebed temperatures and a cross-sectional GW flow model provided estimates of GW-inflow rates that were greatest in spring and decreased through summer. One-dimensional solute transport models calibrated to observed lakebed pore-water dissolved silica (Si) and dissolved phosphate-phosphorus (DP) concentrations indicated that nutrients were transported from the lakebed by advection, diffusion, and enhanced mixing by benthic organisms and waves, and that DP removal occurred near the lakebed interface. Estimated water, Si, DP and total-phosphorus (TP) budgets indicated that GW contributed 21% of lake water inflow and at least 26, 20 and 16% of total Si, DP and TP inflow, respectively, when conservatively assuming background GW nutrient concentrations. However, lakebed GW (LGW) is enriched in nutrients during flow through lakebed sediment and the estimated GW contribution increased to 29 (33), 49 (67) and 43% (61%) of total Si, DP and TP inflow, respectively, if 20% (50%) of GW inflow to the lake was assumed to have LGW concentrations. Net nutrient inflow to the lake was greatest in spring and coincident with the annual diatom bloom. Inflowing dissolved nutrients appear to be assimilated by diatoms during the spring and become available for the summer Aphanizomenon flos-aquae bloom when the diatoms senesce. Thus, nutrient-enriched GW inflow and nutrient recycling by successive algal blooms must be considered when evaluating internal nutrient loading to lakes.

Oregon↗

Benthic vertical hydraulic gradients in Upper Klamath Lake, Oregon, 2017

Groundwater piezometers and lake stilling wells were deployed as paired sets at 10 locations in Upper Klamath Lake in south-central Oregon from May to October 2017 to measure hydraulic heads in and beneath the lake. Continuous water-level data from piezometers and stilling wells were then used to calculate the vertical hydraulic gradient (VHG) across the sediment-water interface to determine the direction and relative magnitude of the movement of water between the lake and underlying sediments. Over the study period, heads in lake-bed sediments closely tracked lake levels, both decreasing from spring into autumn. Instantaneous VHG was highly dynamic at all sites and exhibited high-frequency (less than 1 day to less than 1 hour) variations in magnitude and direction. Instantaneous and weekly mean VHG values often exceeded, but were commonly within, the range of measurement uncertainty (VHG less than +0.009 foot per foot [ft/ft] and greater than -0.009 ft/ft). 63 percent of instantaneous VHG values and 66 percent of weekly mean VHG values were within this range. Study period mean VHG was within measurement uncertainty at seven of the nine sites that had continuous water-level data, but two littoral sites (LC03 and LS01) had positive (upward) values greater than measurement uncertainty and are likely locations of vertical groundwater seepage. Data collected in this study provide new information about the hydraulic conditions at the sediment-water interface in UKL and demonstrate that sediment-groundwater exchange in UKL is spatially and temporally heterogeneous.

Oregon↗