Community Data Mining Approach for Surface Complexation Database Development.

This paper presents a comprehensive data-to-model workflow, including a findable, accessible, interoperable, reusable (FAIR) community sorption database (newly developed LLNL Surface Complexation/Ion Exchange (L-SCIE) database) along with a data fitting workflow to efficiently optimize surface complexation reaction constants with multiple surface complexation model (SCM) constructs. This workflow serves as a universal framework to mine, compile, and analyze large numbers of published sorption data as well as to estimate reaction constants for parameterizing reactive transport models. The framework includes (1) data digitization from published papers, (2) data unification including unit conversions, and (3) data-model integration and reaction constant estimation using geochemical software PHREEQC coupled with the universal parameter estimation code PEST. We demonstrate our approach using an analysis of U(VI) sorption to quartz based on a first L-SCIE implementation, concluding that a multisite SCM construct with carbonate surface species yielded the best fit to community data. Surface complexation reaction constants extracted from this approach captured all available sorption data available in the literature and provided insight into previously published reaction constants and surface complexation model constructs. The L-SCIE sorption database presented herein allows for automating this approach across a wide range of metals and minerals and implementing novel machine learning approaches to reactive transport in the future.

PyLEnM: A Machine Learning Framework for Long-Term Groundwater Contamination Monitoring Strategies.

In this study, we have developed a comprehensive machine learning (ML) framework for long-term groundwater contamination monitoring as the Python package PyLEnM (Python for Long-term Environmental Monitoring). PyLEnM aims to establish the seamless data-to-ML pipeline with various utility functions, such as quality assurance and quality control (QA/QC), coincident/colocated data identification, the automated ingestion and processing of publicly available spatial data layers, and novel data summarization/visualization. The key ML innovations include (1) time series/multianalyte clustering to find the well groups that have similar groundwater dynamics and to inform spatial interpolation and well optimization, (2) the automated model selection and parameter tuning, comparing multiple regression models for spatial interpolation, (3) the proxy-based spatial interpolation method by including spatial data layers or in situ measurable variables as predictors for contaminant concentrations and groundwater levels, and (4) the new well optimization algorithm to identify the most effective subset of wells for maintaining the spatial interpolation ability for long-term monitoring. We demonstrate our methodology using the monitoring data at the Savannah River Site F-Area. Through this open-source PyLEnM package, we aim to improve the transparency of data analytics at contaminated sites, empowering concerned citizens as well as improving public relations.

Improving Long-term Monitoring of Contaminated Groundwater at Sites where Attenuation-based Remedies are Deployed.

This study presents an effective approach to tackle the challenge of long-term monitoring of contaminated groundwater sites where remediation leaves residual contamination in the subsurface. Traditional long-term monitoring of contaminated groundwater sites focuses on measuring contaminant concentrations and is applicable to sites where contaminant mass is removed or degraded to a level below the regulatory standard. The traditional approach is less effective at sites where risk from metals or radionuclides continues to exist in the subsurface after remedial goals are achieved. We propose a long-term monitoring strategy for this type of waste site that focuses on measuring the hydrological and geochemical parameters that control attenuation or remobilization of contaminants while de-emphasizing contaminant-concentration measurements. We demonstrate how this approach would be more effective than traditional long-term monitoring, using a site in South Carolina, USA, where groundwater is contaminated by several radionuclides. A comprehensive enhanced attenuation remedy has been implemented at the site to minimize discharge of contamination to surface water. The immobilization of contaminants occurs in three locations by manipulation of hydrological and geochemical parameters, as well as by natural attenuation processes. Deployment of our proposed long-term monitoring strategy will combine subsurface and surface measurements using spectroscopic tools, geophysical tools, and sensors to monitor the parameters controlling contaminant attenuation. The advantage of this approach is that it will detect the potential for contaminant remobilization from engineered and natural attenuation zones, allowing potential adverse changes to be mitigated before contaminant attenuation is reversed.

In Situ Monitoring of Groundwater Contamination Using the Kalman Filter.

This study presents a Kalman filter-based framework to establish a real-time in situ monitoring system for groundwater contamination based on in situ measurable water quality variables, such as specific conductance (SC) and pH. First, this framework uses principal component analysis (PCA) to identify correlations between the contaminant concentrations of interest and in situ measurable variables. It then applies the Kalman filter to estimate contaminant concentrations continuously and in real-time by coupling data-driven concentration-decay models with the previously identified data correlations. We demonstrate our approach with historical groundwater data from the Savannah River Site F-Area: We use SC and pH data to estimate tritium and uranium concentrations over time. Results show that the developed method can estimate these contaminant concentrations based on in situ measurable variables. The estimates remain reliable with less frequent or no direct measurements of the contaminant concentrations, while capturing the dynamics of short- and long-term contaminant concentration changes. In addition, we show that data mining, such as PCA, is useful to understand correlations in groundwater data and to design long-term monitoring systems. The developed in situ monitoring methodology is expected to improve long-term groundwater monitoring by continuously confirming the contaminant plume's stability and by providing an early warning system for unexpected changes in the plume's migration.

Water Table Dynamics and Biogeochemical Cycling in a Shallow, Variably-Saturated Floodplain.

Three-dimensional variably saturated flow and multicomponent biogeochemical reactive transport modeling, based on published and newly generated data, is used to better understand the interplay of hydrology, geochemistry, and biology controlling the cycling of carbon, nitrogen, oxygen, iron, sulfur, and uranium in a shallow floodplain. In this system, aerobic respiration generally maintains anoxic groundwater below an oxic vadose zone until seasonal snowmelt-driven water table peaking transports dissolved oxygen (DO) and nitrate from the vadose zone into the alluvial aquifer. The response to this perturbation is localized due to distinct physico-biogeochemical environments and relatively long time scales for transport through the floodplain aquifer and vadose zone. Naturally reduced zones (NRZs) containing sediments higher in organic matter, iron sulfides, and non-crystalline U(IV) rapidly consume DO and nitrate to maintain anoxic conditions, yielding Fe(II) from FeS oxidative dissolution, nitrite from denitrification, and U(VI) from nitrite-promoted U(IV) oxidation. Redox cycling is a key factor for sustaining the observed aquifer behaviors despite continuous oxygen influx and the annual hydrologically induced oxidation event. Depth-dependent activity of fermenters, aerobes, nitrate reducers, sulfate reducers, and chemolithoautotrophs (e.g., oxidizing Fe(II), S compounds, and ammonium) is linked to the presence of DO, which has higher concentrations near the water table.

Extrapolating active layer thickness measurements across Arctic polygonal terrain using LiDAR and NDVI data sets.

Landscape attributes that vary with microtopography, such as active layer thickness (ALT), are labor intensive and difficult to document effectively through in situ methods at kilometer spatial extents, thus rendering remotely sensed methods desirable. Spatially explicit estimates of ALT can provide critically needed data for parameterization, initialization, and evaluation of Arctic terrestrial models. In this work, we demonstrate a new approach using high-resolution remotely sensed data for estimating centimeter-scale ALT in a 5 km2 area of ice-wedge polygon terrain in Barrow, Alaska. We use a simple regression-based, machine learning data-fusion algorithm that uses topographic and spectral metrics derived from multisensor data (LiDAR and WorldView-2) to estimate ALT (2 m spatial resolution) across the study area. Comparison of the ALT estimates with ground-based measurements, indicates the accuracy (r2 = 0.76, RMSE ±4.4 cm) of the approach. While it is generally accepted that broad climatic variability associated with increasing air temperature will govern the regional averages of ALT, consistent with prior studies, our findings using high-resolution LiDAR and WorldView-2 data, show that smaller-scale variability in ALT is controlled by local eco-hydro-geomorphic factors. This work demonstrates a path forward for mapping ALT at high spatial resolution and across sufficiently large regions for improved understanding and predictions of coupled dynamics among permafrost, hydrology, and land-surface processes from readily available remote sensing data.

Resilient design in nuclear energy: Critical lessons from a cross-disciplinary analysis of the Fukushima Dai-ichi nuclear accident.

This paper presents a multidisciplinary analysis of the Fukushima Dai-ichi Nuclear Power Plant accident. Along with the latest observations and simulation studies, we synthesize the time-series and event progressions during the accident across multiple disciplines, including in-plant physics and engineering systems, operators' actions, emergency responses, meteorology, radionuclide release and transport, land contamination, and health impacts. We identify three key factors that exacerbated the consequences of the accident: (1) the failure of Unit 2 containment venting, (2) the insufficient integration of radiation measurements and meteorology data in the evacuation strategy, and (3) the limited risk assessment and emergency preparedness. We conclude with new research and development directions to improve the resilience of nuclear energy systems and communities, including (1) meteorology-informed proactive venting, (2) machine learning-enabled adaptive evacuation zones, and (3) comprehensive risk-informed emergency planning while leveraging the experience from responses to other disasters.

Estimating geographic variation of infection fatality ratios during epidemics.

Objectives: We aim to estimate geographic variability in total numbers of infections and infection fatality ratios (IFR; the number of deaths caused by an infection per 1,000 infected people) when the availability and quality of data on disease burden are limited during an epidemic. Methods: We develop a noncentral hypergeometric framework that accounts for differential probabilities of positive tests and reflects the fact that symptomatic people are more likely to seek testing. We demonstrate the robustness, accuracy, and precision of this framework, and apply it to the United States (U.S.) COVID-19 pandemic to estimate county-level SARS-CoV-2 IFRs. Results: The estimators for the numbers of infections and IFRs showed high accuracy and precision; for instance, when applied to simulated validation data sets, across counties, Pearson correlation coefficients between estimator means and true values were 0.996 and 0.928, respectively, and they showed strong robustness to model misspecification. Applying the county-level estimators to the real, unsimulated COVID-19 data spanning April 1, 2020 to September 30, 2020 from across the U.S., we found that IFRs varied from 0 to 44.69, with a standard deviation of 3.55 and a median of 2.14. Conclusions: The proposed estimation framework can be used to identify geographic variation in IFRs across settings.

Nuclear waste Educator's workshop: What and how do we teach about nuclear waste?

A workshop was held at the Massachusetts Institute of Technology (MIT) on July 25th and 26th, 2022. The objective was to develop a blueprint for educating next-generation engineers and scientists about nuclear waste management and disposal, which requires knowledge from diverse disciplines, including nuclear, chemical, civil, environmental, and geological science and engineering. The 49 participants included university professors, researchers, industry experts, and government officials from different areas. First, we have developed a list of key fundamental knowledge on waste management and disposal across the nuclear fuel cycle. In addition, we discussed strategies on how to teach students with diverse backgrounds through innovative teaching strategies as well as how to attract students into this area. Through the workshop, we identified the critical needs to (1) develop community resources for nuclear waste education; (2) synthesize historical perspectives, including past contamination and the management of general hazardous waste; (3) emphasize a complete life-cycle perspective, including proper waste management as the key component for energy sustainability; (4) teach students how to communicate about the key facts and risks to technical and non-technical audiences; and (5) accelerate the use of the state-of-art-technologies to attract and retain a young workforce. Furthermore, we aim to build a diverse, inclusive community that supports students in developing their own narratives about nuclear waste, particularly in recognizing that antagonistic views have been important to improving safety and protecting public health and the environment.

Spatial and temporal prediction of radiation dose rates near Fukushima Daiichi Nuclear Power Plant.

In this paper, we have developed a methodology to estimate the spatiotemporal distribution of radiation air dose rates around the Fukushima Daiichi Nuclear Power Plant (FDNPP). In our exploratory data analysis, we found that (1) the temporal evolution of dose rates is composed of a log-linear decay trend and fluctuations of air dose rates that are spatially correlated among adjacent monitoring posts; and (2) the slope of the log-linear environmental decay trend can be represented as a function of the apparent initial dose rates, coordinate position, land-use type, and soil type. From these observations, we first estimated the log-linear decay trend at each location based on these predictors, using the random forest method. We then developed a modified Kalman filter coupled with a Gaussian process model to estimate the dose-rate time series at a given location and time. We applied this method to the Fukushima evacuation zone (as of March 2017), which included 17 monitoring post locations (with monitoring datasets collected between 2014 and 2018) and generated a time series of dose-rate maps. Our results show that this approach allows us to produce accurate spatial and temporal predictions of radiation dose-rate maps using limited spatiotemporal measurements.

Surface parameters and bedrock properties covary across a mountainous watershed: Insights from machine learning and geophysics.

Bedrock property quantification is critical for predicting the hydrological response of watersheds to climate disturbances. Estimating bedrock hydraulic properties over watershed scales is inherently difficult, particularly in fracture-dominated regions. Our analysis tests the covariability of above- and belowground features on a watershed scale, by linking borehole geophysical data, near-surface geophysics, and remote sensing data. We use machine learning to quantify the relationships between bedrock geophysical/hydrological properties and geomorphological/vegetation indices and show that machine learning relationships can estimate most of their covariability. Although we can predict the electrical resistivity variation across the watershed, regions of lower variability in the input parameters are shown to provide better estimates, indicating a limitation of commonly applied geomorphological models. Our results emphasize that such an integrated approach can be used to derive detailed bedrock characteristics, allowing for identification of small-scale variations across an entire watershed that may be critical to assess the impact of disturbances on hydrological systems.

Influence of soil heterogeneity on soybean plant development and crop yield evaluated using time-series of UAV and ground-based geophysical imagery.

Understanding the interactions among agricultural processes, soil, and plants is necessary for optimizing crop yield and productivity. This study focuses on developing effective monitoring and analysis methodologies that estimate key soil and plant properties. These methodologies include data acquisition and processing approaches that use unmanned aerial vehicles (UAVs) and surface geophysical techniques. In particular, we applied these approaches to a soybean farm in Arkansas to characterize the soil-plant coupled spatial and temporal heterogeneity, as well as to identify key environmental factors that influence plant growth and yield. UAV-based multitemporal acquisition of high-resolution RGB (red-green-blue) imagery and direct measurements were used to monitor plant height and photosynthetic activity. We present an algorithm that efficiently exploits the high-resolution UAV images to estimate plant spatial abundance and plant vigor throughout the growing season. Such plant characterization is extremely important for the identification of anomalous areas, providing easily interpretable information that can be used to guide near-real-time farming decisions. Additionally, high-resolution multitemporal surface geophysical measurements of apparent soil electrical conductivity were used to estimate the spatial heterogeneity of soil texture. By integrating the multiscale multitype soil and plant datasets, we identified the spatiotemporal co-variance between soil properties and plant development and yield. Our novel approach for early season monitoring of plant spatial abundance identified areas of low productivity controlled by soil clay content, while temporal analysis of geophysical data showed the impact of soil moisture and irrigation practice (controlled by topography) on plant dynamics. Our study demonstrates the effective coupling of UAV data products with geophysical data to extract critical information for farm management.

Meanders as a scaling motif for understanding of floodplain soil microbiome and biogeochemical potential at the watershed scale.

BACKGROUND: Biogeochemical exports from watersheds are modulated by the activity of microorganisms that function over micron scales. Here, we tested the hypothesis that meander-bound regions share a core microbiome and exhibit patterns of metabolic potential that broadly predict biogeochemical processes in floodplain soils along a river corridor. RESULTS: We intensively sampled the microbiomes of floodplain soils located in the upper, middle, and lower reaches of the East River, Colorado. Despite the very high microbial diversity and complexity of the soils, we reconstructed 248 quality draft genomes representative of subspecies. Approximately one third of these bacterial subspecies was detected across all three locations at similar abundance levels, and ~ 15% of species were detected in two consecutive years. Within the meander-bound floodplains, we did not detect systematic patterns of gene abundance based on sampling position relative to the river. However, across meanders, we identified a core floodplain microbiome that is enriched in capacities for aerobic respiration, aerobic CO oxidation, and thiosulfate oxidation with the formation of elemental sulfur. Given this, we conducted a transcriptomic analysis of the middle floodplain. In contrast to predictions made based on the prominence of gene inventories, the most highly transcribed genes were relatively rare amoCAB and nxrAB (for nitrification) genes, followed by genes involved in methanol and formate oxidation, and nitrogen and CO2 fixation. Within all three meanders, low soil organic carbon correlated with high activity of genes involved in methanol, formate, sulfide, hydrogen, and ammonia oxidation, nitrite oxidoreduction, and nitrate and nitrite reduction. Overall, the results emphasize the importance of sulfur, one-carbon and nitrogen compound metabolism in soils of the riparian corridor. CONCLUSIONS: The disparity between the scale of a microbial cell and the scale of a watershed currently limits the development of genomically informed predictive models describing watershed biogeochemical function. Meander-bound floodplains appear to serve as scaling motifs that predict aggregate capacities for biogeochemical transformations, providing a foundation for incorporating riparian soil microbiomes in watershed models. Widely represented genetic capacities did not predict in situ activity at one time point, but rather they define a reservoir of biogeochemical potential available as conditions change. Video abstract.

Optimizing long-term monitoring of radiation air-dose rates after the Fukushima Daiichi Nuclear Power Plant.

Radiation air dose rates near the Fukushima Daiichi Nuclear Power Plant (FDNPP) have been steadily decreasing over the past eight years since the release of radioactive elements in March 2011. Currently, the radiation monitoring program is expected to transition to long-term monitoring after most of the remediation activities are completed. The main long-term monitoring objectives are to (1) confirm the continuing reduction of contaminant and hazard levels, (2) provide assurance for the public, (3) accumulate the basic datasets for scientific knowledge and future preparation, and (4) detect changes or anomalies in contaminant mobility (if they occur), or any unexpected processes or events. In this work, we have developed a methodology for optimizing the monitoring locations of radiation air dose-rate monitoring. Our approach consists of three steps in order to determine monitoring locations in a systematic manner: (1) prioritizing the critical locations, such as schools or regulatory requirement locations, (2) diversifying locations that cover the key environmental controls that are known to influence contaminant mobility and distributions, and (3) capturing the heterogeneity of radiation air-dose rates across the domain. For the second step, we use a Gaussian mixture model to identify the representative locations among multiple environmental variables, such as elevation and land-cover types. For the third step, we use a Gaussian process model to capture and estimate the heterogeneity of air-dose rates across the domain. Employing an integrated dose-rate map derived from Bayesian geostatistical methods as a reference map, we distribute the monitoring locations in such a way as to capture the heterogeneity of the reference map. Our results have shown that this approach allows us to select monitoring locations in a systematic manner such that the heterogeneity of air dose rates is captured by the minimal number of monitoring locations.

Characterizing regional-scale temporal evolution of air dose rates after the Fukushima Daiichi Nuclear Power Plant accident.

In this study, we quantify the temporal changes of air dose rates in the regional scale around the Fukushima Dai-ichi Nuclear Power Plant in Japan, and predict the spatial distribution of air dose rates in the future. We first apply the Bayesian geostatistical method developed by Wainwright et al. (2017) to integrate multiscale datasets including ground-based walk and car surveys, and airborne surveys, all of which have different scales, resolutions, spatial coverage, and accuracy. This method is based on geostatistics to represent spatial heterogeneous structures, and also on Bayesian hierarchical models to integrate multiscale, multi-type datasets in a consistent manner. We apply this method to the datasets from three years: 2014 to 2016. The temporal changes among the three integrated maps enables us to characterize the spatiotemporal dynamics of radiation air dose rates. The data-driven ecological decay model is then coupled with the integrated map to predict future dose rates. Results show that the air dose rates are decreasing consistently across the region. While slower in the forested region, the decrease is particularly significant in the town area. The decontamination has contributed to significant reduction of air dose rates. By 2026, the air dose rates will continue to decrease, and the area above 3.8 µSv/h will be almost fully contained within the non-residential forested zone.

Evaluating temporal controls on greenhouse gas (GHG) fluxes in an Arctic tundra environment: An entropy-based approach.

There is significant spatial and temporal variability associated with greenhouse gas (GHG) fluxes in high-latitude Arctic tundra environments. The objectives of this study are to investigate temporal variability in CO2 and CH4 fluxes at Barrow, AK and to determine the factors causing this variability using a novel entropy-based classification scheme. In particular, we analyzed which geomorphic, soil, vegetation and climatic properties most explained the variability in GHG fluxes (opaque chamber measurements) during the growing season over three successive years. Results indicate that multi-year variability in CO2 fluxes was primarily associated with soil temperature variability as well as vegetation dynamics during the early and late growing season. Temporal variability in CH4 fluxes was primarily associated with changes in vegetation during the growing season and its interactions with primary controls like seasonal thaw. Polygonal ground features, which are common to Arctic regions, also demonstrated significant multi-year variability in GHG fluxes. Our results can be used to prioritize field sampling strategies, with an emphasis on measurements collected at locations and times that explain the most variability in GHG fluxes. For example, we found that sampling primary environmental controls at the centers of high centered polygons in the month of September (when freeze-back period begins) can provide significant constraints on GHG flux variability - a requirement for accurately predicting future changes to GHG fluxes. Overall, entropy results document the impact of changing environmental conditions (e.g., warming, growing season length) on GHG fluxes, thus providing clues concerning the manner in which ecosystem properties may be shifted regionally in a future climate.

Climate change impact on residual contaminants under sustainable remediation.

This study investigates the potential impact of climate change on residual contaminants in vadose zones and groundwater. We assume that the effect of climate changes can be represented by perturbations in the natural recharge through the aquifer system. We perform numerical modeling of unsaturated/saturated flow and transport and consider different performance metrics: contaminant concentrations at observation wells and contaminant export at the site's boundary. We evaluate the effect of increasing and decreasing recharge as well as the impact of potential failure of surface capping structures employed to immobilize vadose zone contaminants. Our approach is demonstrated in a real case study by simulating transport of non-reactive radioactive tritium at the U.S. Department of Energy's Savannah River Site. Results show that recharge changes significantly affect well concentrations: after an initial slight dilution we identify a significant concentration increase at different observation wells some years after the recharge increase and/or the cap failure, as a consequence of contaminants' mobilization. This effect is generally emphasized and occurs earlier as the recharge increases. Under decreased aquifers' recharge the concentration could slightly increase for some years, due to a decrease of dilution, depending on the magnitude of the negative recharge shift. We identify trigger levels of recharge above which the concentration/export breakthrough curves and the time of exceedance of the Maximum Contaminant Level for tritium are remarkably affected. Moreover, we observe that the contaminant export at the control plane, identified as the risk pathway to the downgradient population, may only be minimally affected by shifts in the natural recharge regime, except for some extreme cases. We conclude that more frequent sampling and in-situ monitoring near the source zone should be adopted to better explain concentrations' anomalies under changing climatic conditions. Moreover, the maintenance of the cap is critical not only to sequester residual contaminants in the vadose zone, but also to reduce the uncertainty associated with future precipitation changes. Finally, realistic flow and transport simulations achieved through proper calibration processes, rather than conservative modeling, should be adopted to identify non-trivial trade-offs which enable better allocation of resources towards reducing uncertainty in decision making.

Characterizing regional-scale temporal evolution of air dose rates after the Fukushima Daiichi Nuclear Power Plant accident.

In this study, we quantify the temporal changes of air dose rates in the regional scale around the Fukushima Dai-ichi Nuclear Power Plant in Japan, and predict the spatial distribution of air dose rates in the future. We first apply the Bayesian geostatistical method developed by Wainwright et al. (2017) to integrate multiscale datasets including ground-based walk and car surveys, and airborne surveys, all of which have different scales, resolutions, spatial coverage, and accuracy. This method is based on geostatistics to represent spatial heterogeneous structures, and also on Bayesian hierarchical models to integrate multiscale, multi-type datasets in a consistent manner. We apply this method to the datasets from three years: 2014 to 2016. The temporal changes among the three integrated maps enables us to characterize the spatiotemporal dynamics of radiation air dose rates. The data-driven ecological decay model is then coupled with the integrated map to predict future dose rates. Results show that the air dose rates are decreasing consistently across the region. While slower in the forested region, the decrease is particularly significant in the town area. The decontamination has contributed to significant reduction of air dose rates. By 2026, the air dose rates will continue to decrease, and the area above 3.8 µSv/h will be almost fully contained within the non-residential forested zone.

Comparison of Electrostatic and Non-Electrostatic Models for U(VI) Sorption on Aquifer Sediments.

A non-electrostatic generalized composite surface complexation model (SCM) was developed for U(VI) sorption on contaminated F-Area sediments from the U.S. Department of Energy Savannah River Site, South Carolina. The objective of this study was to test if a simpler, semi-empirical, non-electrostatic U(VI) sorption model (NEM) could achieve the same predictive performance as a SCM with electrostatic correction terms in describing U(VI) plume evolution and long-term mobility. One-dimensional reactive transport simulations considering key hydrodynamic processes, Al and Fe minerals, as well as H+ and U surface complexation, with and without electrostatic correction terms, were conducted. The NEM was first calibrated with laboratory batch H+ and U(VI) sorption data on F-Area sediments, and then the surface area of the NEM was adjusted to match field observations of dissolved U(VI). Modeling results indicate that the calibrated NEM was able to perform as well as the previously developed electrostatic model in predicting the long-term evolution of H+ and U(VI) at the site, given the variability of field-site data. The electrostatic and NEM models yield somewhat different results for the time period when basin discharge was active; however, it is not clear which modeling approach may be better to model this early time period because groundwater quality data during this period were not available. A key finding of this study is that the applicability of NEM (and thus robustness of its predictions) to the field system evolves with time and is strongly dependent on the pH range that was used to develop the model.