Metamodels to Bridge the Gap Between Modeling and Decision Support

Insights from process-based models are a mainstay of many groundwater investigations; however, long runtimes often preclude their use in the decision-making process. Screening-level predictions are often needed in areas lacking time or funding for rigorous process-based modeling. The U.S. Geological Survey (USGS) Groundwater Resources and National Water Quality Assessment Programs are addressing these issues by evaluating the “metamodel” to bridge these gaps. A metamodel is a statistical model founded on a computationally expensive model. Although faster, the question remains: Can a statistical model provide similar insights to a numerical model with faster results

Metamodels to Bridge the Gap Between Modeling and Decision Support

Insights from process-based models are a mainstay of many groundwater investigations; however, long runtimes often preclude their use in the decision-making process. Screening-level predictions are often needed in areas lacking time or funding for rigorous process-based modeling. The U.S. Geological Survey (USGS) Groundwater Resources and National Water Quality Assessment Programs are addressing these issues by evaluating the “metamodel” to bridge these gaps. A metamodel is a statistical model founded on a computationally expensive model. Although faster, the question remains: Can a statistical model provide similar insights to a numerical model with faster results

Obtaining Parsimonious Hydraulic Conductivity Fields Using Head and Transport Observations: A Bayesian Geostatistical Parameter Estimation Approach

Flow path delineation is a valuable tool for interpreting the subsurface hydrogeochemical environment. Different types of data, such as groundwater flow and transport, inform different aspects of hydrogeologic parameter values (hydraulic conductivity in this case) which, in turn, determine flow paths. This work combines flow and transport information to estimate a unified set of hydrogeologic parameters using the Bayesian geostatistical inverse approach. Parameter flexibility is allowed by using a highly parameterized approach with the level of complexity informed by the data. Despite the effort to adhere to the ideal of minimal a priori structure imposed on the problem, extreme contrasts in parameters can result in the need to censor correlation across hydrostratigraphic bounding surfaces. These partitions segregate parameters into facies associations. With an iterative approach in which partitions are based on inspection of initial estimates, flow path interpretation is progressively refined through the inclusion of more types of data. Head observations, stable oxygen isotopes (18O/16O ratios), and tritium are all used to progressively refine flow path delineation on an isthmus between two lakes in the Trout Lake watershed, northern Wisconsin, United States. Despite allowing significant parameter freedom by estimating many distributed parameter values, a smooth field is obtained

Bridging groundwater models and decision support with a Bayesian network

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Water Resources Research 49 (2013): 6459–6473, doi:10.1002/wrcr.20496.Resource managers need to make decisions to plan for future environmental conditions, particularly sea level rise, in the face of substantial uncertainty. Many interacting processes factor in to the decisions they face. Advances in process models and the quantification of uncertainty have made models a valuable tool for this purpose. Long-simulation runtimes and, often, numerical instability make linking process models impractical in many cases. A method for emulating the important connections between model input and forecasts, while propagating uncertainty, has the potential to provide a bridge between complicated numerical process models and the efficiency and stability needed for decision making. We explore this using a Bayesian network (BN) to emulate a groundwater flow model. We expand on previous approaches to validating a BN by calculating forecasting skill using cross validation of a groundwater model of Assateague Island in Virginia and Maryland, USA. This BN emulation was shown to capture the important groundwater-flow characteristics and uncertainty of the groundwater system because of its connection to island morphology and sea level. Forecast power metrics associated with the validation of multiple alternative BN designs guided the selection of an optimal level of BN complexity. Assateague island is an ideal test case for exploring a forecasting tool based on current conditions because the unique hydrogeomorphological variability of the island includes a range of settings indicative of past, current, and future conditions. The resulting BN is a valuable tool for exploring the response of groundwater conditions to sea level rise in decision support.This work was funded by the USGS Climate and Land Use Mission Area, Research and Development Program and the USGS Natural Hazards Mission Area, Coastal and Marine Geology Program

An Interactive Bayesian Geostatistical Inverse Protocol for Hydraulic Tomography

Hydraulic tomography is a powerful technique for characterizing heterogeneous hydrogeologic parameters. An explicit trade-off between characterization based on measurement misfit and subjective characterization using prior information is presented. We apply a Bayesian geostatistical inverse approach that is well suited to accommodate a flexible model with the level of complexity driven by the data and explicitly considering uncertainty. Prior information is incorporated through the selection of a parameter covariance model characterizing continuity and providing stability. Often, discontinuities in the parameter field, typically caused by geologic contacts between contrasting lithologic units, necessitate subdivision into zones across which there is no correlation among hydraulic parameters. We propose an interactive protocol in which zonation candidates are implied from the data and are evaluated using cross validation and expert knowledge. Uncertainty introduced by limited knowledge of dynamic regional conditions is mitigated by using drawdown rather than native head values. An adjoint state formulation of MODFLOW-2000 is used to calculate sensitivities which are used both for the solution to the inverse problem and to guide protocol decisions. The protocol is tested using synthetic two-dimensional steady state examples in which the wells are located at the edge of the region of interest

Effects of sea-level rise on barrier island groundwater system dynamics – ecohydrological implications

This paper is not subject to U.S. copyright. The definitive version was published in Ecohydrology 7 (2014): 1064–1071, doi:10.1002/eco.1442.We used a numerical model to investigate how a barrier island groundwater system responds to increases of up to 60 cm in sea level. We found that a sea-level rise of 20 cm leads to substantial changes in the depth of the water table and the extent and depth of saltwater intrusion, which are key determinants in the establishment, distribution and succession of vegetation assemblages and habitat suitability in barrier islands ecosystems. In our simulations, increases in water-table height in areas with a shallow depth to water (or thin vadose zone) resulted in extensive groundwater inundation of land surface and a thinning of the underlying freshwater lens. We demonstrated the interdependence of the groundwater response to island morphology by evaluating changes at three sites. This interdependence can have a profound effect on ecosystem composition in these fragile coastal landscapes under long-term changing climatic conditions.We used a numerical model to investigate how a barrier island groundwater system responds to increases of up to 60 cm in sea level. We found that a sea-level rise of 20 cm leads to substantial changes in the depth of the water table and the extent and depth of saltwater intrusion, which are key determinants in the establishment, distribution and succession of vegetation assemblages and habitat suitability in barrier islands ecosystems. In our simulations, increases in water-table height in areas with a shallow depth to water (or thin vadose zone) resulted in extensive groundwater inundation of land surface and a thinning of the underlying freshwater lens. We demonstrated the interdependence of the groundwater response to island morphology by evaluating changes at three sites. This interdependence can have a profound effect on ecosystem composition in these fragile coastal landscapes under long-term changing climatic conditions. Published 2013. This article is a U.S. Government work and is in the public domain in the USA

Understanding the DayCent model: Calibration, sensitivity, and identifiability through inverse modeling

AbstractThe ability of biogeochemical ecosystem models to represent agro-ecosystems depends on their correct integration with field observations. We report simultaneous calibration of 67 DayCent model parameters using multiple observation types through inverse modeling using the PEST parameter estimation software. Parameter estimation reduced the total sum of weighted squared residuals by 56% and improved model fit to crop productivity, soil carbon, volumetric soil water content, soil temperature, N2O, and soil NO3− compared to the default simulation. Inverse modeling substantially reduced predictive model error relative to the default model for all model predictions, except for soil NO3− and NH4+. Post-processing analyses provided insights into parameter–observation relationships based on parameter correlations, sensitivity and identifiability. Inverse modeling tools are shown to be a powerful way to systematize and accelerate the process of biogeochemical model interrogation, improving our understanding of model function and the underlying ecosystem biogeochemical processes that they represent

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

A multi-resolution, non-parametric, Bayesian framework for identification of spatially-varying model parameters

This paper proposes a hierarchical, multi-resolution framework for the identification of model parameters and their spatially variability from noisy measurements of the response or output. Such parameters are frequently encountered in PDE-based models and correspond to quantities such as density or pressure fields, elasto-plastic moduli and internal variables in solid mechanics, conductivity fields in heat diffusion problems, permeability fields in fluid flow through porous media etc. The proposed model has all the advantages of traditional Bayesian formulations such as the ability to produce measures of confidence for the inferences made and providing not only predictive estimates but also quantitative measures of the predictive uncertainty. In contrast to existing approaches it utilizes a parsimonious, non-parametric formulation that favors sparse representations and whose complexity can be determined from the data. The proposed framework in non-intrusive and makes use of a sequence of forward solvers operating at various resolutions. As a result, inexpensive, coarse solvers are used to identify the most salient features of the unknown field(s) which are subsequently enriched by invoking solvers operating at finer resolutions. This leads to significant computational savings particularly in problems involving computationally demanding forward models but also improvements in accuracy. It is based on a novel, adaptive scheme based on Sequential Monte Carlo sampling which is embarrassingly parallelizable and circumvents issues with slow mixing encountered in Markov Chain Monte Carlo schemes

A Gibbs sampler for inequality‐constrained geostatistical interpolation and inverse modeling

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94819/1/wrcr11621.pd