Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface data models”

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 397 records · Page 22

Play Fairway analysis of geothermal resources across the State of Hawai‘i: 4. Updates with new groundwater chemistry, subsurface stress analysis, and focused geophysical surveys

This paper is the fourth in a series on a play fairway analysis of geothermal resources across the State of Hawai‘i. Here we describe recent exploration activities that include groundwater sampling in ten locations statewide, as well as geophysical surveys on Lana‘i, across the SW rift of Haleakala Volcano (Maui), and surrounding Mauna Kea Volcano (Hawai‘i Island). We derive a first-order method for computing topographic stress using Green’s functions, finding that topography can induce appreciable crustal stress. We develop a new method for incorporating depth information about resistivity, density, and topographic stresses into our previously published equations to calculate resource probability and confidence. We incorporate newly collected and modeled data into our calculations to update statewide maps of probability and confidence. Lana‘i Island and southeast Mauna Kea are identified as the top targets for exploratory drilling. The east rift of Haleakala Volcano, the southern-most region of Mauna Loa’s SW rift (Hawai‘i Island), and central Kauai are identified as targets for geophysical surveying.

15 GEOTHERMAL ENERGY↗

Deep Learning for Simultaneous Inference of Hydraulic and Transport Properties

Abstract Identification of a heterogeneous conductivity field and reconstruction of a contaminant release history are key aspects of subsurface remediation. These two goals are achieved by combining model predictions with sparse and noisy hydraulic head and concentration measurements. Solution of this inverse problem is notoriously difficult due to, in part, high dimensionality of the parameter space and high computational cost of repeated forward solves. We use a convolutional adversarial autoencoder (CAAE) to parameterize a heterogeneous non‐Gaussian conductivity field via a low‐dimensional latent representation. A three‐dimensional dense convolutional encoder‐decoder (DenseED) network serves as a forward surrogate of the flow and transport model. The CAAE‐DenseED surrogate is fed into the ensemble smoother with multiple data assimilation (ESMDA) algorithm to sample from the Bayesian posterior distribution of the unknown parameters, forming a CAAE‐DenseED‐ESMDA inversion framework. The resulting CAAE‐DenseED‐ESMDA inversion strategy is used to identify a three‐dimensional contaminant source and conductivity field. A comparison of the inversion results from CAAE‐ESMDA with physical flow and transport simulator and from CAAE‐DenseED‐ESMDA shows that the latter yields accurate reconstruction results at the fraction of the computational cost of the former.

Zhou, Zitong↗

Relaxation behavior in low-frequency complex conductivity of sands caused by bacterial growth and biofilm formation by Shewanella oneidensis under a high-salinity condition

Complex electrical conductivity is increasingly used to monitor subsurface processes associated with microbial activities because microbial cells mostly have surface charges and thus electrical double layers. Although highly saline environments are frequently encountered in coastal and marine sediments, there are limited data available on the complex conductivity associated with microbial activities under a high-salinity condition. Therefore, we have developed the spectral responses of complex conductivity of sand associated with bacterial growth and biofilm formation under a highly saline condition of approximately 1% salinity and approximately 2 S/m pore water conductivity with an emphasis on relaxation behavior. A column test is performed, in which the model bacteria Shewanella oneidensis MR-1 are stimulated for cell growth and biofilm formation in a sand pack, whereas the complex conductivity is monitored from 0.01 Hz to 10 kHz. The test results indicate that the real conductivity increases in the early stage due to the microbial metabolites and the increased surface conduction with cell growth but soon begin to decrease because of the reduction of charge passages due to bioclogging. However, the imaginary conductivity significantly increases with time, and clear bell-shaped relaxation behaviors are observed with the peak frequency of 0.1–1 Hz, associated with the double-layer polarization of cells and electrically conductive pili and biofilms. The Cole-Cole relaxation model appears to capture such relaxation behaviors well, and the modeling results indicate gradual increases in normalized chargeability and decreases in relaxation time during bacterial growth and biofilm formation in the highly saline condition. Comparison with previous literature confirms that the high-salinity condition further increases the normalized chargeability, whereas it suppresses the phase shift and thus the imaginary conductivity. Our results suggest that the complex conductivity can effectively capture microbial biomass formation in sands under a highly saline condition.

58 GEOSCIENCES↗

Permafrost thermal conditions are sensitive to shifts in snow timing

Abstract Changes in snow precipitation at high latitudes can significantly affect permafrost thermal conditions and thaw depth, potentially exposing more carbon-laden soil to microbial decomposition. A fully coupled process-based surface/subsurface thermal hydrology model with surface energy balance is used to analyze the impact of intra-annual variability in snow on permafrost thermal regime and the active layer thickness. In the four numerical scenarios considered, simulations were forced by the same meteorological data, except the snow precipitation, which was systematically altered to change timing of snowfall. The scenarios represent subtle shifts in snow timing, but the snow onset/melt days, the end of winter snowpack depth, and total annual snow precipitation are unchanged among scenarios. The simulations show that small shifts in the timing of snow accumulation can have significant effects on subsurface thermal conditions leading to active layer deepening and even talik formation when snowfall arrives earlier in the winter. The shifts in snow timing have a stronger impact on wetter regions, especially soil underneath small ponds, as compared to drained regions. This study highlights the importance of understanding potential changes in winter precipitation patterns for reliable projections of active-layer thickness in a changing Arctic climate.

54 ENVIRONMENTAL SCIENCES↗

Geothermal play fairway analysis, part 1: Example from the Snake River Plain, Idaho

The Snake River Plain (SRP) volcanic province overlies the track of the Yellowstone hotspot, a thermal anomaly that extends deep into the mantle. Most of the area is underlain by a basaltic volcanic province that overlies a mid-crustal intrusive complex, which in turn provides the long-term heat flux needed to sustain geothermal systems. Previous studies have identified several known geothermal resource areas within the SRP. For the geothermal study presented herein, our goals were to: (1) adapt the methodology of Play Fairway Analysis (PFA) for geothermal exploration to create a formal basis for its application to geothermal systems, (2) assemble relevant data for the SRP from publicly available and private sources, and (3) build a geothermal PFA model for the SRP and identify the most promising plays, using GIS-based software tools that are standard in the petroleum industry. The study focused on identifying three critical resource parameters for exploitable hydrothermal systems in the SRP: heat source, reservoir and recharge permeability, and cap or seal. Data included in the compilation for heat source were heat flow, distribution and ages of volcanic vents, groundwater temperatures, thermal springs and wells, helium isotope anomalies, and reservoir temperatures estimated using geothermometry. Reservoir and recharge permeability was inferred from the analysis of stress orientations and magnitudes, post-Miocene faults, and subsurface structural lineaments based on magnetics and gravity data. Data for cap or seal included the distribution of impermeable lake sediments and clay-seal associated with hydrothermal alteration below the regional aquifer. These data were used to compile Common Risk Segment maps for heat, permeability, and seal, which were combined to create a Composite Common Risk Segment map for all southern Idaho that reflects the risk associated with geothermal resource exploration and identifies favorable resource tracks. Our regional assessment indicated that undiscovered geothermal resources may be located in several areas of the SRP. Two of these areas, the western SRP and Camas Prairie, were selected for more detailed assessment, during which heat, permeability, and seal were evaluated using newly collected field data and smaller grid parameters to refine the location of potential resources. These higher resolution assessments illustrate the flexibility of our approach over a range of scales.

54 ENVIRONMENTAL SCIENCES↗

Hillslope-Channel Transitions and the Role of Water Tracks in a Changing Permafrost Landscape

The Arctic is experiencing rapid climate change, and the effect on hydrologic processes and resulting geomorphic changes to hillslopes and channels is unclear because we lack quantitative models and theory for rapid changes resulting from thawing permafrost. Here, the presence of permafrost modulates water flow and the stability of soil-mantled slopes, implying that there should be a signature of permafrost processes, including warming-driven disturbance, in channel network extent. To inform understanding of hillslope-channel dynamics under changing climates, we examined soil-mantled hillslopes within a ~300 km 2 area of the Seward Peninsula, western Alaska, where discontinuous permafrost is particularly susceptible to thaw and rapid landscape change. In this study, we pair high-resolution topographic and satellite data to multi-annual observations of InSAR-derived surface displacement over a 5-year period to quantify spatial variations in topographic change across an upland landscape. We find that neither the basin slope nor the presence of knickzones controls the magnitude of recent surface displacements within the study basin, as may be expected under conceptual models of temperate hillslope evolution. Rather, the highest displacement magnitudes tended to occur at the broad hillslope-channel transition zone. In this study area, this zone is occupied by water tracks, which are zero-order ecogeomorphic features that concentrate surface and subsurface flow paths. Our results suggest that water tracks, which appear to occupy hillslope positions between saturation and incision thresholds, are vulnerable to warming-induced subsidence and incision. We hypothesize that gullying within water tracks will outpace infilling by hillslope processes, resulting in the growth of the channel network under future warming.

58 GEOSCIENCES↗

Technoeconomic Analysis of Infrastructure Buildout Scenarios

The success of CCUS deployment in the southeast will depend on the optimized pipeline network from CO 2 point sources to geologic sinks for CO 2 storage. An optimal pipeline framework will reduce both environmental impacts and pipeline building cost. However, finding an optimal pipeline/transport infrastructure design is a non-trivial task and requires simultaneous usage of large volumes of data, computational resources, and a state-of-the-art simulator. Such a transport infrastructure model must consider point sources for CO 2 capture and associated volumes, sinks for CO 2 storage, and transportation from source to sink via pipeline networks. These considerations must be addressed simultaneously and systemically in an optimization process, in which a defined objective function (i.e., capital, variable CO 2 capture, transport, and storage cost function) is required to be minimized with the consideration of practical constraints (e.g., CO 2 flow through pipelines is less than the maximum capacity). Although optimization has been widely used in subsurface resources production and CO 2 storage, its application in large scale CCUS infrastructure design is rarely reported in the literature. SimCCS, developed by Los Alamos National Laboratory, is an open-source CCUS pipeline infrastructure design toolset that facilitates the optimization of pipeline infrastructure networks.

99 GENERAL AND MISCELLANEOUS↗

Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS) - Deliverable Task 7: Risk Assessment

The “Optimizing Alabama’s CO 2 Storage in Shelby County: Project OASIS” CarbonSAFE Phase II Project seeks to build on regional data sets that demonstrate that the subsurface within Shelby County, Alabama has the potential to store commercial volumes of CO 2 safely, permanently, and economically. The primary target CO 2 Storage Complex is the deep Ketona Dolomite located within a 140 square mile area of the Valley and Ridge Region of Alabama. This deep saline reservoir is beneath a confining system encompassing at least 6,500 ft of shales and other low permeability sediments. Project OASIS drilled a deep stratigraphic test well to confirm the geological properties of the confining system and storage reservoir(s) within the Storage Complex. The geological data was incorporated into numerical models to establish the areal extent of the CO 2 plume and help design the storage site and its monitoring system. The project is managed by the Southern States Energy Board (SSEB), an interstate compact organization consisting of governors and state legislative leaders from sixteen southern states, Puerto Rico, and the U.S. Virgin Islands, as well as an appointee by the President of the United States. The organizational compact provides it with access to state government organizations and legislatures. The Board also maintains an Associate Members program comprised of energy resource companies, utilities, trade groups, academic R&D science and technology experts and energy consultants. Further, the SSEB staff is experienced in managing and coordinating complex energy and environmental programs, from research programs to full-scale design and demonstrations of new and innovative technologies. Project OASIS is a public-private partnership of six entities with multiple principal investigators (PIs). SSEB’s Lead PI and Co-PI are responsible for all aspects of project performance in accordance with the DOE-NETL Cooperative Agreement. SSEB has issued subgrants to Advanced Resources International, Inc., Alabama A&M University, Auburn University, Crescent Resource Innovation, and Oklahoma State University. Advanced Resources International, Inc., issued subgrants to Baker Hughes and Loudon Technical Services for field services

20 FOSSIL-FUELED POWER PLANTS↗

Detrital U-Pb zircon and 40 Ar/ 39 Ar muscovite geochronology from Middle Pennsylvanian strata in the Anadarko Basin, Texas Panhandle, USA

The Late Mississippian-Permian Anadarko Basin formed in Texas and Oklahoma, USA as the result of inversion of Neoproterozic and Cambrian rift structures. Subsidence was driven by flexural loading of the Amarillo-Wichita Uplift, and this uplift may represent the easternmost element of the Ancestral Rocky Mountains system. The northwestern part of this basin has generally been interpreted to have been filled by sediment derived from the Ancestral Front Range Uplift, ~475 km to the northwest during the early stages of basin filling. We test this model using U-Pb detrital zircon and 40 Ar/ 39 Ar detrital muscovite results from three subsurface samples of the Morrow B sandstone in the northwestern part of the Anadarko Basin. We provide a new maximum depositional age of 310.9 ± 4.9 Ma that indicates the age of the Morrow B to be no older than late Atokan to early Desmoinesian Age (Moscovian), ~10 Myr younger than previously interpreted. In contrast to some previous interpretations, we propose that the most likely source for the sediment in the Morrow B is the Amarillo Uplift to the south. Detrital zircon and detrital muscovite data have age peaks at 900–1300 Ma, 1370 Ma and 1600–1800 Ma corresponding to derivation from Grenville, Granite-Rhyolite and Yavapai-Mazatzal basement provinces, respectively. A dominant detrital zircon peak at 1370 Ma suggests that Mesoproterozoic granites in the Amarillo Uplift were exposed by Middle Pennsylvanian time, and detritus eroded from the Amarillo Uplift dominated the lowstand sediment budget of the Texas Panhandle during this time; small volumes of sediment were likely sourced from the Ancestral Front Range. In conclusion, this study presents the first detrital geochronology data from the subsurface Anadarko Basin and the first detrital muscovite data from late Paleozoic southwestern Laurentia. The results presented here highlight the interpretive power of combined detrital zircon and muscovite datasets.

58 GEOSCIENCES↗

3D Geologic Framework Modelling of the Los Alamos National Laboratory Site and Pajarito Plateau: Integrating a realistic 3D fault network and modelling subsurface relationships in a sparsely sampled and complex geologic region

The subsurface geology beneath the Pajarito Plateau is critical to understanding the seismic hazard of the Pajarito Fault System, yet our understanding of this geology is relatively poor. While previous 3D geologic framework models of the area have been created for the purposes of understanding hydrogeologic flow, they are inadequate for the purposes of understanding the Pajarito Fault System. The specific challenges of using oil and gas software for this purpose include: (1) the geologic complexities resulting from volcanism and tectonism; (2) a need for a high level of stratigraphic detail over a large area; (3) a near complete lack of seismic data; and (4) sparse wellbore data. Presented here is a workflow that handles these challenges of adapting commercially available software used by the oil and gas industries to this seismic hazard problem.

58 GEOSCIENCES↗

Big Data Meets Geothermal Exploration (CRADA Final Report)

As part of the Cyclotron Road program, Zanskar Geothermal & Minerals, Inc. investigated the application of micro-earthquake and ambient noise seismology methods to imaging and characterizing the structural characteristics and hydrothermal flux of subsurface faults. Significant advances in what could be resolved were enabled by two major developments in seismology: 1) the availability of large-n arrays of low-cost seismometers, and 2) the availability of increased computational power and semi-automated data reduction algorithms. In tandem, these advances may improve the signal-to-noise ratio and spatial precision of the data collected and enable higher-resolution characterization of subsurface fracture systems and their spatio-temporal evolution. These tools supported efforts to reduce dry-hole risk and to improve wellfield productivity for geothermal resource development. In particular, two applications of these advances were evaluated: 1) fracture-seismic imaging, which was used to detect ambient emissions from fluid-filled fractures, and 2) reservoir tomography, which used information about travel paths, source locations, and source parameters of micro-earthquakes to identify areas of enhanced permeability. Integration of these methods provided guidance for siting wells and served as prior constraints for reservoir models, informing forecasts of power potential and production and injection strategies aimed at minimizing temperature decline and improving overall resource productivity.

15 GEOTHERMAL ENERGY↗

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

60 APPLIED LIFE SCIENCES↗

Model Package Report: Geoframework Model of the Hanford Site 100 Area

The purpose of the 100 Area Geoframework Model (GFM) is to provide a reasonable, consistent, and defensible three dimensional representation of the hydrostratigraphic units below the River Corridor at the Hanford Site to support contaminant fate and transport models. The GFM is a three dimensional representation of the subsurface geologic structure. From this, three dimensional geologic model-exported results, in the form of points or surfaces, are used as inputs to populate and assemble the various numerical model architectures. The objective of this report is to define the process used to produce a hydrostratigraphic model for the hydrostratigraphic units beneath the Hanford Site 100 Area. The GFM may support several other CH2M HILL Plateau Remediation Company project needs and objectives, including providing geologic information to support remedial investigations and actions and as a tool to present River Corridor geology and contaminant extents. The GFM was constructed based on information through 2019 available in the Integrated Data Management System, Hanford Environmental Information System, and the Hanford Site geologic contacts (GeoContacts) datasets. Revisions to the 100 Area GFM may occur to incorporate new data and information. Each version will be maintained in configuration control using date stamps for identifying supporting databases, figures, and interpretations. Supporting data include the final three dimensional geoframework surfaces, two dimensional structure and isopach maps, and all the geologic contact inputs and interpreted geologic log data. These final products provide a supporting set of geologic information that together allow the creation of the GFM. This information is managed, updated, and maintained via CH2M HILL Plateau Remediation Company under the Environmental Model Management Archive.

58 GEOSCIENCES↗

Regime shifts in Arctic terrestrial hydrology manifested from impacts of climate warming

Abstract. Anthropogenic warming in the Arctic is causing hydrological cycle intensification and permafrost thaw, with implications for flows of water, carbon, and energy from terrestrial biomes to coastal zones. To better understand the likely impacts of these changes, we used a hydrology model driven by meteorological data from atmospheric reanalysis and two global climate models for the period 1980–2100. The hydrology model accounts for soil freeze–thaw processes and was applied across the pan-Arctic drainage basin. The simulations point to greater changes over northernmost areas of the basin underlain by permafrost and to the western Arctic. An acceleration of simulated river discharge over the recent past is commensurate with trends drawn from observations and reported in other studies. Between early-century (2000–2019) and late-century (2080–2099) periods, the model simulations indicate an increase in annual total runoff of 17 %–25 %, while the proportion of runoff emanating from subsurface pathways is projected to increase by 13 %–30 %, with the largest changes noted in summer and autumn and across areas with permafrost. Most notably, runoff contributions to river discharge shift to northern parts of the Arctic Basin that contain greater amounts of soil carbon. Each season sees an increase in subsurface runoff; spring is the only season where surface runoff dominates the rise in total runoff, and summer experiences a decline in total runoff despite an increase in the subsurface component. The greater changes that are seen in areas where permafrost exists support the notion that increased soil thaw is shifting hydrological contributions to more subsurface flow. The manifestations of warming, hydrological cycle intensification, and permafrost thaw will impact Arctic terrestrial and coastal environments through altered river flows and the materials they transport.

54 ENVIRONMENTAL SCIENCES↗

A computational pipeline to generate a synthetic dataset of metal ion sorption to oxides for AI/ML exploration

The charged mineral/electrolyte interfaces are ubiquitous in the surface and subsurface–including the surroundings of the geological disposal sites for radioactive waste. Therefore, understanding how ions interact with charged surfaces is critically important for predicting radionuclide mobility in the case of waste leakage. At present, the Surface Complexation Models (SCMs) are the most successful thermodynamic frameworks to describe ion retention by mineral surfaces. SCMs are interfacial speciation models that account for the effect of the electric field generated by charged surfaces on sorption equilibria. These models have been successfully used to analyze and interpret a broad range of experimental observations including potentiometric and electrokinetic titrations or spectroscopy. Unfortunately, many of the current procedures to solve and fit SCM to experimental data are not optimal, which leads to a non-transferable or non-unique description of interfacial electrostatics and consequently of the strength and extent of ion retention by mineral surfaces. Recent developments in Artificial Intelligence (AI) offer a new avenue to replace SCM solvers and fitting algorithms with trained AI surrogates. Unfortunately, there is a lack of a standardized dataset covering a wide range of SCM parameter values available for AI exploration and training–a gap filled by this study. Here, we described the computational pipeline to generate synthetic SCM data and discussed approaches to transform this dataset into AI-learnable input. First, we used this pipeline to generate a synthetic dataset of electrostatic properties for a broad range of the prototypical oxide/electrolyte interfaces. The next step is to extend this dataset to include complex radionuclide sorption and complexation, and finally, to provide trained AI architectures able to infer SCMs parameter values rapidly from experimental data. Here, we illustrated the AI-surrogate development using the ensemble learning algorithms, such as Random Forest and Gradient Boosting. These surrogate models allow a rapid prediction of the SCM model parameters, do not rely on an initial guess, and guarantee convergence in all cases.

Li, Chunhui↗

Photoelectrochemically Self Improving Si/GaN Photocathode: Figure 4a Raw Data

Calculated formation energy of (0001 ̅) c-plane GaN surfaces. The density-functional theory (DFT) calculations were performed with the generalized gradient approximation, using the plane-wave pseudopotential quantum-ESPRESSO package. Ultrasoft pseudopotentials were used to describe the interaction between valence electrons and ionic cores, where Ga 3d electrons were explicitly treated as valence electrons. The wave functions and electronic density were expanded in a plane-wave basis set truncated at a cutoff energy of 30 Ry and 240 Ry, respectively. The surface was modeled using a repeated slab geometry using a (2×2) unit cell with a thickness of seven GaN bilayers, where the Ga-terminated bottom side was passivated by pseudo-hydrogen with charge of 1.25e. In addition, the (1010) surface was modeled using 8 GaN atomic layers with a lateral dimension of 9.57 Å x 10.38 Å. A vacuum width of at least 15.0 Å was introduced between consecutive slabs, and a k-point mesh of 5x5x1 and 2x2x1 was used to sample the surface Brillouin zone for the polar and non-polar surfaces, respectively. Calculated formation energy of (0001 ̅) c-plane GaN surfaces with different oxygen configuration within the anion and cation-rich limits: ideal GaN surfaces of c-plane and m-plane; configurations with 100% replacement of nitrogen by oxygen on the surface c-plane and m-plane; 100% replacement of nitrogen by oxygen in the subsurface bilayer c-plane and m-plane; and configurations with 50% replacement of nitrogen of the surface and subsurface bilayers by oxygens c-plane and m-plane. In both cases, the ideal surface is chosen as the reference.

photocathode↗