Correction: Generating H2 during the CO2 sequestration in basalt formations
In this article [1], the author’s name Carl Steefel was incorrectly written as Carl Stefeel. The original article has been corrected.
Huang, Lizhi↗
Engineering topics
Steefel, Carl
Publications and source records attributed to Steefel, Carl.
HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite
Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation
Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.
58 GEOSCIENCES↗
Aerobic respiration controls on shale weathering, Geochimica et Cosmochimica Acta, 2023: Dataset
This data package was generated in order to support the development of a deep-time weathering model and to assess the coupling between shale weathering and aerobic respiration in the paper “Aerobic respiration controls on shale weathering” by Stolze et al., Geochimica et Cosmochimica Acta (2023). The package contains two csv files providing the average CO2(g) concentration profiles [ppm] and mineral concentration profiles [wt%], respectively. The CO2(g) concentration profiles were measured in the vicinity of the monitoring well PLM2 between January 2018 and April 2019. The gas samples were collected in the unsaturated zone to a depth of 1.52 m. The mineral concentration profiles were determined by X-Ray diffraction (XRD). The XRD measurements were performed on sub-core samples collected in the monitoring well PLM3 down to a depth of 7.01 m. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.Update on 2024-05-28: Revised versions of the CSV data files (CO2_data_GCA_Stolze_et_al_2023.csv and XRD_data_GCA_Stolze_et_al_2023.csv) were made to apply ESS-DIVE's CSV reporting format guidelines. Updated versions of the File Level Metadata (v2_20240528_flmd.csv) and Data Dictionary (v2_20240528_dd.csv) files were updated to reflect the changes made to the CSV files.
54 ENVIRONMENTAL SCIENCES↗
Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling
Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.
54 ENVIRONMENTAL SCIENCES↗
Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model
Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.
58 GEOSCIENCES↗
Soil thickness map at two hillslopes near the pumphouse in the east river watershed, Colorado
The soil thickness maps were created by using a hybrid model-data approach. Field sampling and remote sensing data of the spatial distribution of two hillslopes in the Pump House area in the East River Watershed in the CO., the U.S. The data package includes the geospatial data of the soil thickness maps at two hillslopes near the pumphouse, and the associated remote sensing data, including lidar DEM and a shape file of the boundary of the study area. The data can be viewed in GIS software such as QGIS or ArcGIS desktop. The geospatial data can also be viewed in Python or Matlab. The data were generated for the purpose of modeling surface hydrology and near-surface chemistry. This work shows how to combine sampling data and a process-based model to predict one of the highest uncertainty in the land surface process, the soil thickness.
2D Geospatial Maps↗
ATS (Advanced Terrestrial Simulator) integrated hydrology and reactive transport model output in Copper Creek, Colorado
This dataset is generated using the ATS (Advanced Terrestrial Simulator) model at Copper Creek, Colorado, the largest catchment in the East River watershed. ATS is an integrated hydrology and reactive transport model to simulate the Concentration-Discharge (C-Q) relationship, and is used to quantify the geochemical export and understand how watersheds respond to climate disturbances. Key hydrogeochemical variables, including discharge, precipitation, evapotranspiration, and geochemical concentrations, were simulated and output at a daily frequency at the Copper Creek outlet. A group of hydrology and geochemical model variables over the spatial domain are output at a bi-monthly frequency. The model used 1-km resolution daily frequency Daymet meteorological forcing, and mineral composition from previous studies. The simulations are performed on NERSC. The datasets include the input files for running the ATS (Advanced Terrestrial Simulator) simulation, including the mesh file, input script and geochemical model files. The output files from the ATS (Advanced Terrestrial Simulator) are the daily time-series discharge and geochemical concentration simulated at the outlet of Copper Creek, and the bi-monthly spatial output over the simulation domain.
54 ENVIRONMENTAL SCIENCES↗
An Integrated Multiscale Modeling Framework for Unconventional Stimulation and Production (Final Report)
The production of oil and gas from unconventional reservoirs largely depends upon two main features operating at different scales: (1) the establishment of a reservoir scale stimulated fracture network that effectively communicates with the rock volume, enhancing permeability and transport to the wellbore and (2) the coupled multi-phase flow, chemical and mechanical processes affecting the migration of hydrocarbons from the low permeability country rock adjacent to the stimulated fracture network. To date, there has been no simulation framework that allows seamless and integrated prediction of these features across spatial scales extending from the pore structure of the reservoir rock to the volume of the reservoir. In addition, there has been a lack of suitable field measurements to test such models, as stimulation and production data are often proprietary and not freely available to national laboratories and academic institutions. New multi-scale simulation capabilities are needed that are validated against suitable field-based research experiments on hydraulic fracturing and shale production.
02 PETROLEUM↗
Hysteresis Patterns of Watershed Nitrogen Retention and Loss Over the Past 50 years in United States Hydrological Basins
Patterns of watershed nitrogen (N) retention and loss are shaped by how watershed biogeochemical processes retain, biogeochemically transform, and lose incoming atmospheric deposition of N. Loss patterns represented by concentration, discharge, and their associated stream exports are important indicators of integrated watershed N retention behaviors. We examined continental United States (CONUS) scale N deposition (e.g., wet and dry atmospheric deposition), vegetation trends, and stream trends as potential indicators of watershed N-saturation and retention conditions, and how watershed N retention and losses vary over space and time. By synthesizing changes and modalities in watershed nitrogen loss patterns based on stream data from 2200 U.S. watersheds over a 50 years record, our work revealed two patterns of watershed N-retention and loss. One was a hysteresis pattern that reflects the integrated influence of hydrology, atmospheric inputs, land-use, stream temperature, elevation, and vegetation. The other pattern was a one-way transition to a new state. We found that regions with increasing atmospheric deposition and increasing vegetation health/biomass patterns have the highest N-retention capacity, become increasingly N-saturated over time, and are associated with the strongest declines in stream N exports—a pattern, that is, consistent across all land cover categories. We provide a conceptual model, validated at an unprecedented scale across the CONUS that links instream nitrogen signals to upstream mechanistic landscape processes. Our work can aid in the future interpretation of in-stream concentrations of DOC and DIN as indicators of watershed N-retention status and integrators of watershed hydrobiogeochemical processes.
58 GEOSCIENCES↗
Knowledge-Guided Machine Learning (KGML) Platform to Predict Integrated Water Cycle and Associated extremes
Focal Area(s): Predictive modeling through the use of AI techniques and insight gleaned from complex data (both observed and simulated). Science Challenge: Although advanced predictive capabilities of the water cycle are critical to address environmental needs and develop sustainable solutions for energy demands, there is no robust framework, to say the least, that seamlessly integrates local to intermediate to global scales and a gamut of biogeophysical information to enhance understanding of the integrated water cycle and its associated extremes.
54 ENVIRONMENTAL SCIENCES↗
Learning from learning machines: improving the predictive power of energy-water-land nexus models with insights from complex measured and simulated data
Focal Area(s): Insights gleaned from complex data (observed and simulated) using AI; Predictive modeling through the use of AI techniques and AI-derived model components, including physics- and knowledge-informed models; Energy-water-land nexus and integrated energy systems – models of MultiSector dynamics. Science Challenge: Scientific communities are in need of tools for the computational integration of physics-based models, experimental data, and empirical/observational studies across a broad range of temporal and spatial scales to explore and meet EESSD grand challenges. We require AI algorithms for the discovery of process-drivers in the Earth-energy-human system and to link them with complex measured and simulated data. We aim to understand the energy-water-land nexus, particularly under extreme forcing scenarios and rare events, leveraging scale-aware AI process models, including probabilistic uncertainties, and benefiting from the emerging 5G-enabled landscape-scale sensing and edge computing capabilities. These models are needed to identify instabilities and tipping points that manifest extreme system behaviors with consequences for integrated energy systems and the environment. This work requires fundamental advances in uncertainty quantification, in particular, to identify and model unlikely but catastrophic outliers. Specifically, we need models that are interpretable to domain scientists and, essentially, explainable to public and private stake holders.
54 ENVIRONMENTAL SCIENCES↗
On Demand Machine Learning for Multi-Fidelity Biogeochemistry in River Basins Impacted by Climate Extremes
Predictive modeling of watershed to river basin scale biogeochemistry through the use of AI techniques and AI-derived model components and the use of AI to design a prediction system comprised of a hierarchy of multi-fidelity models, including AI driven model/component/parameterization selection.
58 GEOSCIENCES↗
Modeling the Impact of Riparian Hollows on River Corridor Nitrogen Exports
Recent studies in snowmelt-dominated catchments have documented changes in nitrogen (N) retention over time, such as declines in watershed exports of N, though there is a limited understanding of the controlling processes driving these trends. Working in the mountainous headwater East River Colorado watershed, we explored the effects of riparian hollows as N-cycling hotspots and as important small-scale controls on observed watershed trends. Using a modeling-based approach informed by remote sensing and in situ observations, we simulated the N-retention capacity of riparian hollows with seasonal and yearly hydrobiogeochemical perturbations imposed as drivers. We then implemented a scaling approach to quantify the relative contribution of riparian hollows to the total river corridor N budget. We found that riparian hollows primarily serve as N sinks, with N-transformation rates significantly limited by periods of enhanced groundwater upwelling and promoted at the onset of rainfall events. Given these observed hydrologic controls, we expect that the nitrate ( NO 3 - ) sink capacity of riparian hollows will increase in magnitude with future climatic perturbations, specifically the shift to more frequent rainfall events and fewer snowmelt events, as projected for many mountainous headwater catchments. Our current estimates suggest that while riparian hollows provision ~5–20% of NO 3 - to the river network, they functionally act as inhibitors to upland NO 3 - reaching the stream. Our work linking transient hydrological conditions to numerical biogeochemical simulations is an important step in assessing N-retaining features relative to the watershed N budget and better understanding the role of small-scale features within watersheds.
54 ENVIRONMENTAL SCIENCES↗
Process Discovery through Assimilation of Complex Biogeochemical Datasets
This white paper addresses two of three focal areas identified in the white paper call: 1) Biogeochemical data acquisition and assimilation enabled by machine learning and 3) Insight gleaned from complex data using AI. We focus on AI application to complex biogeochemistry (BGC) data (e.g. laboratory experimental data, field manipulation data, literature data), which is an untapped source of information for improving Earth System Predictability (ESP).
58 GEOSCIENCES↗
Scripts and raw/processed precipitation NRCS network data, Upper Colorado, 2008-2017
This data package consists of scripts and data that were used to generate all the results and figures for Mital et al., 2020 (https://doi.org/10.3389/frwa.2020.00020). The purpose of the study was to develop a new algorithm to impute (or gap-fill) missing daily precipitation data. Study area was the Upper Colorado Water Resources Region (UCWRR), over a time period of 2008-2017. In terms of data, the package consists of raw precipitation measurements from Natural Resources Conservation Service (NRCS) network. These raw measurements are compiled into a single csv file (NRCS_dates_2008_2017.csv). The corresponding metadata are compiled in Metadata_processed.csv. The package also consists of outputs generated during Baseline and Sequential imputation runs as documented in Mital et al., 2020. Scripts used to generate all the results and figures for Mital et al., 2020 are also uploaded. A readme file documenting the layout of the data archive is also uploaded, and a KML file shows the geographic range for the data.
54 ENVIRONMENTAL SCIENCES↗
Downscaled precipitation and mean air temperature datasets; East-Taylor subbasin; 2008-2019; daily temporal resolution; 400 m spatial resolution
This dataset provides gridded meteorological forcing data (specifically, daily precipitation and daily mean air temperature). The dataset has been generated by downscaling the Parameter-elevation Regressions on Independent Slopes Model (PRISM) dataset from a spatial resolution of 800 m to 400 m. The time period is 2008-2019, and the mapped area is East Taylor subbasin in Upper Colorado. The data are in the form of NetCDF files, arranged by year. Compared to the PRISM dataset, the dates of the downscaled dataset are shifted backwards by one (e.g., downscaled data for May 25 corresponds to PRISM data for May 26). This temporal shifting makes the meteorological forcing correspond more closely to the prescribed date. The NetCDF format is a standard raster format that can be read using any Geographic Information System (GIS) software; plenty of modules exist in popular scripting languages (such as Python, R and Matlab) that can also be used to read NetCDF files.East_Taylor_Tavg_PRISM400 corresponds to mean air temperature.East_Taylor_Precip_PRISM400.zip corresponds to precipitation.The proprietary PRISM data (800 m resolution) were purchased with funding from the Watershed Function Scientific Focus Area supported by U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under award no. DE-AC02-05CH11231.Dataset update on March 9, 2022:Uploaded new data files that are identical to the earlier files but are now in the NetCDF format, arranged by year. Earlier, the files were in the GeoTiff format, arranged by date.
54 ENVIRONMENTAL SCIENCES↗
Exascale applications: skin in the game
As noted in Wikipedia, skin in the game refers to having ‘incurred risk by being involved in achieving a goal’, where ‘ skin is a synecdoche for the person involved, and game is the metaphor for actions on the field of play under discussion’. For exascale applications under development in the US Department of Energy Exascale Computing Project, nothing could be more apt, with the skin being exascale applications and the game being delivering comprehensive science-based computational applications that effectively exploit exascale high-performance computing technologies to provide breakthrough modelling and simulation and data science solutions. These solutions will yield high-confidence insights and answers to the most critical problems and challenges for the USA in scientific discovery, national security, energy assurance, economic competitiveness and advanced healthcare. This article is part of a discussion meeting issue ‘Numerical algorithms for high-performance computational science’.
97 MATHEMATICS AND COMPUTING↗