Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Water Storage Trends in High Mountain Asia

Changes in terrestrial water storage (TWS) in High Mountain Asia (HMA) could have major societal impacts, as the region’s large reservoirs of glaciers, snow, and groundwater provide a freshwater source to more than one billion people. We seek to quantify and close the budget of secular changes in TWS over the span of the GRACE satellite mission (2003–2016). To assess the TWS trend budget we consider a new high-resolution mass trend product determined directly from GRACE L1B data, glacier mass balance derived from Digital Elevation Models (DEMs), groundwater variability determined from confined and unconfined well observations, and terrestrial water budget estimates from a suite of land surface model simulations with the NASA Land Information System (LIS). This effort is successful at closing the aggregated TWS trend budget over the entire HMA region, the glaciated portion of HMA, and the Indus and Ganges basins, where the full-region trends are primarily due to the glacier mass balance and groundwater signals. Additionally, we investigate the closure of TWS trends at individual 1-arc-degree mascons (area ≈12,000 km2); a significant improvement in spatial resolution over previous analyses of GRACE-derived trends. This mascon-level analysis reveals locations where the TWS trends are well-explained by the independent datasets, as well as regions where they are not; identifying specific geographic areas where additional data and model improvements are needed. The accurate characterization of total TWS trends and its components presented here is critical to understanding the complex dynamics of the region, and is a necessary step toward projecting future water mass changes in HMA.

Bryant D Loomis↗
DOE OSTI2022

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Role of Atmospheric Rivers on Groundwater: Lessons Learned From an Extreme Wet Year

Abstract In the coastal regions of the western United States, atmospheric rivers (ARs) are associated with the largest precipitation generating storms and contribute up to half of annual precipitation, but the impact of ARs on the integrated hydrologic cycle, specifically on groundwater storage and hydrodynamics, is largely unknown. To better explore the hydrologic behavior of AR versus non‐AR event precipitation, we present a novel combination of two water tracking methods (one in the atmosphere and one in the subsurface) to explicitly track the full lifecycle of water parcels generated by ARs. Simulations of northern California's Cosumnes River watershed during the record wet 2017 water year are performed via the coupling of a high‐resolution regional climate model and a land surface‐groundwater model accounting for lateral groundwater flow. Despite ARs contributing more precipitation than non‐AR storms, we find less AR water is preferentially stored in aquifers by year end. Fractionally, ARs result in 300% less snow derived groundwater‐recharged compared to non‐AR precipitation. Rain‐on‐snow (RoS) plays an important role in AR‐driven discharge, where over 50% of total discharge from ARs snow is from RoS events. Finally, despite record‐breaking annual precipitation, simulated groundwater depletion occurs by year end due to estimates of groundwater pumping activities. The results from these simulations serve as a partial analogue of future hydrologic conditions where ARs are expected to intensify and provide a greater fraction of annual precipitation due to climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

A Preliminary Relationship Between the Depth of Martian Gullies and the Abundance of Hydrogen on Near-Surface Mars

Recent compelling evidence has been presented to suggest the presence of near-surface water or water ice on Mars. The Mars Orbiter Camera (MOC) has photographed relatively young fluvial features in the form of gullies which have been attributed to groundwater seepage. The Gamma Ray Spectrometer (GRS) from the Mars Odyssey spacecraft has detected large amounts of hydrogen in the Martian soil, inferred to come from water or water ice within the upper meter of the surface. We explore the model of groundwater seepage as the mechanism of gully formation as opposed to other mechanisms. We investigate the abundance of hydrogen to the depth at which gully-like features form. We see a positive correlation, which implies the presence of a groundwater system for Mars.

Frey, E. L.↗
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Relating land surface, water table, and weathering fronts with a conceptual valve model for headwater catchments

Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the earth-air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re-orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes toward an “ideal hill”. In this steady state, reaction-front valves partition water volumes into shallow and deep flowpaths. Furthermore, these flows dissolve low- and high-solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Groundwater and Intruder Radionuclide Screening

Solid Waste Management (SWM) operates the E-Area Low-Level Waste Facility (ELLWF) where low level on-site and off-site solid waste streams are buried. The facility has been in operation since late 1994 and is currently projected to remain in operation until 2065. This facility can accommodate a broad range of waste forms resulting from the six different types of disposal unit options (i.e., varying degrees of engineered barriers → trenches to concrete vaults). This facility is currently operating under a Performance Assessment (PA) issued back in 2008, along with several subsequent supporting Special Analyses (SA). The Savannah River National Laboratory (SRNL) developed the prior PAs and SAs and has been tasked to update the facility’s PA. For operating the E-Area facility, a Waste Inventory Tracking System (e.g., WITS) is actively employed by waste generators where every radionuclide entering the facility (to be buried in one of its many disposal units) must be either directly or indirectly accounted for. Since there is a large number of radionuclides in existence (>3,000), the International Commission on Radiological Protection (ICRP-107) has provided guidance on the subset of radionuclides requiring further assessment in landfills such as the E-Area LLWF. The ICRP-107 publication provides critical radiological information on 1,252 radionuclides of 97 elements. This database is the critical starting point for developing a consistent inventory limit system. The explicit measurement and tracking of all radionuclides are not necessary when process knowledge, burial history, and radiological aspects are factored into conservative groundwater and intruder screening processes. Across the DOE complex these screening processes have been historically performed using the methodology suggested by the National Council on Radiation Protection and Measurements (NCRP) as presented in their original report (NCRP 1984) and then refined in a later report (NCRP 1996). In the recommended screening models employed within this report the traditional NCRP models are updated to better handle progeny and better reflect the known characteristics of E-Area. This improved screening process is referred to as the “NCRP-like” method. Various upgrades to the traditional NCRP methodology have been used by others (e.g., NRC funded effort by Kennedy and Strenge 1992). The more detailed models by Kennedy and Strenge (1992) were considered in this effort and are compared to the results from the more traditional NCRP-like models. The groundwater and inadvertent intruder screening analyses presented in this report start with this 1,252 radionuclide list and reduces it down to more manageable lists that are applicable to the various disposal unit types contained within E-Area. In order to reduce this starting list, some level of exposure risk must be considered acceptable. Historically, a dose (or concentration level) has been compared with a screening criterion set to 1% of a performance measure (e.g., a beta-gamma dose not to exceed the 4 mrem/yr beta-gamma performance measure x 0.01 = 0.04 mrem/yr). Thus, if a radionuclide produced a bounding or screening-level dose (or concentrations) less than the screening criterion, it could be safely removed from further consideration. Several tiers of screening and bounding level analyses have been considered in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Groundwater and Intruder Radionuclide Screening

Solid Waste Management (SWM) operates the E-Area Low-Level Waste Facility (ELLWF) where low level on-site and off-site solid waste streams are buried. The facility has been in operation since late 1994 and is currently projected to remain in operation until 2065. This facility can accommodate a broad range of waste forms resulting from the six different types of disposal unit options (i.e., varying degrees of engineered barriers → trenches to concrete vaults). This facility is currently operating under a Performance Assessment (PA) issued back in 2008, along with several subsequent supporting Special Analyses (SA). The Savannah River National Laboratory (SRNL) developed the prior PAs and SAs and has been tasked to update the facility’s PA. For operating the E-Area facility, a Waste Inventory Tracking System (e.g., WITS) is actively employed by waste generators where every radionuclide entering the facility (to be buried in one of its many disposal units1) must be either directly or indirectly accounted for. Since there is a large number of radionuclides in existence (>3,000), the International Commission on Radiological Protection (ICRP-107) has provided guidance on the subset of radionuclides requiring further assessment in landfills such as the E-Area LLWF. The ICRP-107 publication provides critical radiological information on 1,252 radionuclides of 97 elements. This database is the critical starting point for developing a consistent inventory limit system. The explicit measurement and tracking of all 1,252 ICRP-107 radionuclides can be reduced when process knowledge, burial history, and radiological aspects are factored into conservative groundwater and intruder screening processes. Across the DOE complex these screening processes have been historically performed using the methodology suggested by the National Council on Radiation Protection and Measurements (NCRP) as presented in their original report (NCRP 1984) and then refined in a later report (NCRP 1996). In the recommended screening models employed within this report the traditional NCRP models are updated to better handle progeny and better reflect the known characteristics of E-Area. This improved screening process is referred to as the “NCRP-like” method. Various upgrades to the traditional NCRP methodology have been used by others (e.g., NRC funded effort by Kennedy and Strenge 1992). The more detailed models by Kennedy and Strenge (1992) were considered in this effort and are compared to the results from the more traditional NCRP-like models. The groundwater and inadvertent intruder screening analyses presented in this report start with this 1,252 radionuclide list and reduces it down to more manageable lists that are applicable to the various disposal unit types contained within E-Area. In order to reduce this starting list, some level of exposure risk must be considered acceptable. Historically, a dose (or concentration level) has been compared with a screening criterion set to 1% of a performance measure (e.g., a betagamma dose not to exceed the 4 mrem/yr beta-gamma performance measure x 0.01 = 0.04 mrem/yr). Thus, if a radionuclide produced a bounding or screening-level dose (or concentrations) less than the screening criterion, it could be safely removed from further consideration. Several tiers of screening and bounding level analyses have been considered in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2023

South Repeater Building (Swmu 121) Pilot Study Work Plan, Kennedy Space Center, Florida

This document is a Pilot Study Work Plan for a groundwater remediation system at the South Repeater Building (N6-1118) located at the John F. Kennedy Space Center (KSC), Florida. Aqueous fire-fighting foam was applied near the building in the late 1990s during a wildfire. Per- and polyfluoroalkyl substances (PFAS) are components of the foam. Sampling completed at and near the South Repeater Building identified PFAS in groundwater, surface water and soil. KSC is planning to design and install a groundwater remediation system to address the potential for off-Center migration of PFAS compounds in groundwater. This Pilot Study Work Plan outlines subsurface hydrogeologic data that will be collected to use in the construction of a groundwater flow model that will guide the design of a full-scale remediation system. This document was prepared by AECOM Technical Services, Inc., for NASA under Indefinite Delivery Indefinite Quantity Contract 80KSC019D0010, Task Order 80KSC019F0096.

PFAS↗
DOE OSTI2021

Unclassified Source Term and Radionuclide Data for the Groundwater Flow and Contaminant Transport Model for Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada, Revision 1

This report documents the evaluation of information and data available on the unclassified source term and radionuclide (RN) contamination for Central and Western Pahute Mesa (PM): Corrective Action Units (CAUs) 101 and 102. The Central and Western PM CAUs constitute two areas of the Nevada National Security Site (NNSS) used for underground nuclear testing (Figure 1-1). Eighty-five underground nuclear tests were conducted at the PM CAUs between 1965 and 1992 (NNSA/NFO, 2015b). Underground nuclear tests deposit radioactive elements in the subsurface, resulting in groundwater contamination in the vicinity of the underground test locations. At the PM CAUs, the testing resulted in more than 2.68 × 10 7 curies (Ci) of radioactivity (decay-corrected to September 30, 2012) remaining underground (Finnegan et al., 2016). Of the 85 underground nuclear tests conducted at the PM CAUs, three were shallow tests detonated within 200 meters (m) of the surface and referred to as the “near-surface tests”; the remaining 82 were detonated deeper than 200 m below ground surface (bgs) and are referred to as the “deep underground tests” in this document.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗
DOE OSTI2024

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2020

Groundwater Drought: Environmental Controls and Monitoring

Arising from a long-term deficit of precipitation, groundwater drought often lags meteorological drought by months to years, depending on climate and subsurface hydrogeological conditions. Due to the paucity of in situ groundwater data in many parts of the world, global monitoring of groundwater storage variations and drought may be best achieved using satellite observations and/or groundwater time series simulated by hydrological models. The Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow On (GRACE-FO) missions have greatly benefitted the modeling and monitoring of groundwater storage changes and drought at the global scale. In this chapter, we first review environmental controls on the temporal variability of groundwater using in situ data. We then describe an approach that infuses GRACE and GRACE-FO observations into a land surface model for assessing groundwater storage changes and drought globally. We also discuss characteristics of simulated groundwater drought and the limitations of current groundwater drought monitoring approaches.

Bailing Li↗
DOE OSTI2024

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Inverse Analysis with Variational Autoencoders: A Comparison of Shallow and Deep Networks

Inverse problems are applied to determine unknown properties by matching observational data with a physical model that often takes many parameters as input. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented involving regularization with a technique known as a variational autoencoder (VAE), which is trained to map a high-dimensional parameter space with a complex structure to a low-dimensional latent space with a simple structure. We apply this approach to unconditioned realizations of the parameters (heterogeneous hydraulic fields) for a hydrogeological inverse problem. Two types of hydraulic conductivity fields are used to evaluate the characterization for the different levels of heterogeneity complexity of the physical inputs. This approach keeps the computational cost of generating the training data low. The reason is unconditioned realizations neither rely on the observational data used to perform the inverse analysis nor require any groundwater flow model forward runs. In addition, this approach applies regularization on a low-dimensional latent space from the VAE and increases optimization efficiency through automatic differentiation. Furthermore, two different neural network (NN) structures are tested for their utility in using the VAE for inverse analysis. The performance of a deep, convolutional neural network strongly depends on the dimensionality of the latent space, which requires tuning. In contrast, a shallow, dense neural network provides consistently accurate characterization without tuning. Furthermore, our approach evaluates the advantages of the shallow, dense neural network over the deep, convolutional one and enables future application to a wide range of inverse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Calibration and Uncertainty Estimation Using the Ensemble Kalman Filter with a Large Subsurface Flow and Transport Model - 20321

At routinely monitored groundwater contamination sites, periodically measured environmental conditions such as groundwater levels and contaminant concentrations are used to inform and confirm a conceptual site model (CSM) and guide the development and calibration of a numerical groundwater flow and transport model. The calibration of groundwater flow and transport models after each measurement (sampling) event can illuminate deficiencies in a CSM, identify areas where additional monitoring is warranted, and predict the behavior of the system to guide decision making. However, manual and automated (e.g. PEST) model calibration tools can be time-consuming and computationally expensive to implement after each sampling event. Perhaps as a result, such calibration tools generally utilize all available monitoring data simultaneously rather than sequentially assimilating monitoring data one sampling event at a time as the results from sampling become available. A more real-time data assimilation approach may reduce parameter uncertainty, quantify the value of additional monitoring data, and produce a usable model more quickly and with less effort. To mitigate the potential time-consuming aspects of manual and widely applied automated calibration techniques, a data assimilation algorithm called the ensemble Kalman filter (EnKF) was evaluated as a relatively efficient method of model calibration and uncertainty assessment via the sequential integration of monitoring data into a model. The EnKF was able to successfully and efficiently assimilate monitoring and modeling data to calibrate a complex flow and transport model at a real-world site with significant subsurface heterogeneity, uncertainty, and 12 years of monitoring data (over 4,000 individual measurements of groundwater levels and over 2,500 measurements of contaminant concentrations). Starting with an uncalibrated model data from annual sampling events were sequentially assimilated, and the resultant predication errors and estimated parameter uncertainties were tracked. After all monitoring data were assimilated, both flow and transport residuals at the end of the EnKF process were comparable to those produced via a concurrent PEST calibration effort but required fewer model simulations. Both uncertainty and prediction errors decreased over time. In a real-time application, the adequacy of the model could be assessed after each sampling event. The benefits of such a real-time approach to utilizing monitoring data include reduced costs (in the form of model updates or site characterization efforts), early flagging of possible errors in the CSM, and a reduced risk of overfitting and corresponding increased confidence in model predictions. This tool may be particularly useful compared to other calibration techniques (e.g. manual, PEST) when model runtimes are long, calibration parameters are many, or parameter uncertainty is large. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗