Search NASA⌕ Search

SEARCH · Search NASA

Results for “geochemical interactions”

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 73 records · Page 4

Noctis Landing: A Proposed Landing Site/Exploration Zone for Human Missions to the Surface of Mars

The proposed Noctis Landing Site/Exploration Zone (LS/EZ) is shown in Figure 1. Our preliminary study suggests that the proposed site meets all key Science and Resources (incl. Civil Engineering) requirements. The site is of significant interest, as the EZ not only offers a large number and wide range of regions of interest (ROIs) for short-term exploration, it is also located strategically at the crossroads between Tharsis and Valles Marineris, which are key for long-term exploration. The proposed site contains Regions of Interest (ROIs) that meet the following Science requirements: -­‐ Access to (1) deposits with a high preservation potential for evidence of past habitability and fossil biosignatures and (2) sites that are promising for present habitability. The site presents a wide variety of ROIs qith likely aqueous features and deposits, including sinous channels and valleys, slope gullies, lobate debris aprons, impact craters with lobate ejecta flows, and "bathtub ring" deposits. Neutron spectrometry also suggests hydrogen is present within the topmost 0.3 m or so of 4 to 10 wt% WEH (Water Equivalent Hydrogen). -­‐ Noachian and/or Hesperian rocks in a stratigraphic context that have a high likelihood of containing trapped atmospheric gases. Collapsed canyon rim material with preserved stratigraphy is abundantly present and accessible. -­‐ Exposures of at least two crustal units that have regional or global extents, that are suitable for radiometric dating, and that have relative ages that sample a significant range of martian geological time. Canyons floors in Ius Chasma, Tithonium Chasma, and plateau tops on Tharsis and in Sinai Planum offer access to distinct crustal units of regional extent. -­‐ Access to outcrops with linked morphological and/or geochemical signatures indicative of aqueous or groundwater/ mineral interactions. Iron and sulfur-bearing deposits on canyon floors in Noctis Labyrinthus, and in Ius Chasma (IC) and Tithonium Chasma (TC) offer many such outcrop options. -­‐ Identifiable stratigraphic contacts and cross-cutting relationships from which relative ages can be determined. In place and collapsed canyon walls in NL, TC, and IC offer such opportunities. -­‐ Other types of ROIs include access points to surrounding plateau top areas for longer term regional exploration. A key attribute of the proposed Noctic Landing site is its strategic location to allow the shortest possible surface excusions to Tharsis and Valles Marineris (VM). VM is the feature and region on Mars that exposes the longest record of Mars' geology and evolution through time. Tharsis is the region of Mars that has experienced the longest and most extensive volcanic history, and might still be volcanically active. Some of the youngest lava flows on Mars have been identified on the western flanks of the Tharsis Bulge, i.e., within driving range of future longrange (500 - 1000 km) pressurized rover traverses. The proposed site also contains ROIs that offer the following Resources (incl. Civil Engineering) characteristics: -­‐ Access to raw material that exhibits the potential to (1) be used as feedstock for water-generating in situ resource utilization (ISRU) processes and (2) yield significant quantities (greater than 100 MT) of water. The raw material is likely in the form of hydrated minerals, and possibly ice/regolith mix. The top of the raw material deposit is at the surface. -­‐ Access to a region where infrastructure can be emplaced or constructed. This region is less than 5 km from the LS and contains flat, stable terrain. The region exhibits evidence for an abundant source of loose regolith. Several deep pits in the area combined with the availability of sand suggests that some natural terrain features can be adapted for construction purposes. -­‐ Access to raw material that exhibits the potential to be used as metal feedstock for ISRU and construction purposes. Iron and sulfur-rich mineral surface deposits have been identified in CRISM data in many locations in this area. Noctis Landing is the lowest-altitude location on Mars that straddles both the Tharsis region (above average geothermal gradients) and Valles Marineris (minimal crustal thickness from surface (valley floor) to a subsurface liquid water table. Noctis Landing has the potential for being an ideal site for eventual deep drilling on Mars to access deep subsurface liquid water and potentially encountering extant life. Available data remains insufficient to fully qualify the Noctis Landing site. Additional remote sensing data (visible, Near and Mid-IR, and radar) and surface reconnaissance via a high-mobility robotic rover are recommended. In particular, it will be important to assess the trafficability of the site, and its potential for yielding water and metals as a resource. Access to plateau tops from the Noctis Landing site on the canyon floor should be demonstrated. Future exploration of the site would also be enhanced significantly by the availability of robotic (tele-operatable) surveying and sample-collecting drones. Testing of the use of such collaborative science and exploration technologies should be conducted at terrestrial sites such as the Haughton-Mars Project site on Devon Island, High Arctic, among others. Note: Noctis Landing is not an official Mars nomenclature name for this location. Because the area of the proposed LS/EZ had no name, and because it is close to Noctis Labyrinthus to the West while being distinct from it, the provisional name Noctis Landing is proposed. Noctis means night in Latin.

Lee, Pascal↗

Nanoscale Interfacial Smoothing and Dissolution during Unconventional Reservoir Stimulation Implications for Hydrocarbon Mobilization and Transport

Hydraulic fracturing of low-permeability rocks significantly enhances hydrocarbon production from unconventional reservoirs. However, fluid transport through low-permeability rocks, and the influence of geochemical transformations on pore networks, are poorly constrained. Mineral reactivity during interactions with injected water may alter the physical nature of the rock, which may affect hydrocarbon mobility. To assess alteration to the rock, we previously conducted a hydrothermal experiment that reacted cubed rock samples (1 cm3) with synthetic hydraulic fracturing fluid to simulate physicochemical reactivity during hydraulic fracturing. We analyzed unreacted and reacted rocks by small-angle neutron scattering and high-pressure mercury intrusion to determine how nanoscale (~1-1400 nm) pore networks of unconventional reservoir rocks are influenced by reaction with hydraulic fracturing injectates. We show that fluid-rock interactions have a two-fold influence on hydrocarbon recovery, promoting both hydrocarbon mobilization and transport. Pore-matrix interfaces smooth via removal of clay mineral surface asperities, reducing available surface area for hydrocarbon adsorption by ca. 12-75%. Additionally, hydraulic fracturing fluid-induced dissolution creates new pores with diameters ranging in size from 800-1400 nm, increasing the permeability of the rocks by a factor of 5-10. These two consequences of mineral dissolution act in concert to release hydrocarbons from the host rock and facilitate transport through the rock during unconventional reservoir production.

Herz-Thyhsen, Ryan↗

The Mineralogical Record of Oxygen Fugacity Variation and Alteration in Northwest Africa 8159: Evidence for Interaction Between a Mantle Derived Martian Basalt and a Crustal Component(s)

A prominent geochemical feature of basaltic magmatism on Mars is the large range in initial Sr isotopic ratios (approx. 0.702 - 0.724) and initial epsilon-Nd values (approx. -10 to greater than +50). Within this range, the shergottites fall into three discreet subgroups. These subgroups have distinct bulk rock REE patterns, mineral chemistries (i.e. phosphate REE patterns, Ni, Co, V in olivine), oxygen fugacity of crystallization, and stable isotopes, such as O. In contrast, nakhlites and chassignites have depleted epsilon-Nd values (greater than or equal to +15), have REE patterns that are light REE enriched, and appear to have crystallized near the FMQ buffer. The characteristics of these various martian basalts have been linked to different reservoirs in the martian crust and mantle, and their interactions during the petrogenesis of these magmas. These observations pose interesting interpretive challenges to our understanding of the conditions of the martian mantle (e.g. oxygen fugacity) and the interaction of mantle derived magmas with the martian crust and surface. Martian meteorite NWA 8159 is a unique fine-grained augite basalt derived from a highly depleted mantle source as reflected in its initial epsilon-Nd value, contains a pronounced light REE depleted pattern, and crystallized presumably under very oxidizing conditions. Although considerably older than both shergottites and nahklites, it has been petrogenetically linked to both styles of martian magmatism. These unique characteristics of NWA 8159 may provide an additional perspective for deciphering the petrogenesis of martian basalts and the nature of the crust of Mars.

Shearer, Charles K.↗

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Borehole Disposal: Update on the Science, Technology, and Potentially Suitable Radioactive Wastes - 20170

A recent Special Issue of the Energies Journal on Deep Borehole Disposal of Nuclear Waste has delivered a timely update on the science and technology of borehole disposal and the types of radioactive wastes it could potentially accommodate. The Special Issue papers discuss (i) circumstances under which a national waste management programme might wish to consider deep borehole disposal [1], (ii) a status report of deep borehole disposal options in Germany [2], (iii) the use of deep horizontal drill-holes in sedimentary, metamorphic or igneous rocks for disposal of spent nuclear fuel and high-level radioactive waste [3], (iv) the analysis of corrosion performance of engineered barrier systems [4], (v) a review of the potential cementing systems suitable for deep borehole disposal [5], (vi) the thermal evolution around heat-generating waste for a wide range of material properties and disposal configurations [6], (vii) a geochemical analysis of deep brines focussed on fluid-rock interactions [7], (viii) post-closure performance assessment calculations for deep borehole disposal of Cs/Sr capsules [8], and (ix) an example safety case for deep borehole disposal of nuclear wastes [9]. This paper provides an overview of the main findings from these publications and highlights some of the technological and scientific challenges that need to be overcome when developing a deep borehole disposal concept. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Geochemical Phosphorus Sequestration in Tundra Soils Impedes Delivery of Bioavailable Phosphorus to the Kuparuk River, Alaska, USA: Implications for the Broader Arctic Region

Long-term river monitoring of the Kuparuk River (North Slope, Alaska, USA) confirms significant increases in solutes that are indicative of active layer thickening due to thawing permafrost. However, there is no evidence of an increase in total dissolved phosphorus (TDP) or soluble reactive phosphorus (SRP), the nutrient that limits primary production in this and similar rivers in the region. Here, we show that Mehlich-3 extractable iron (Fe) and aluminum (Al) in active layer soils impart high P geochemical sorption capacities across a range of landscape features that we would expect to promote lateral movement of water and solutes to headwater streams in our study watershed. Reanalysis of a recently published pan-Arctic soils database that includes active layer and permafrost soil samples suggests that this high P sorption capacity could be common in other parts of the Arctic region. We conclude that soil minerals enhance P retention on hillslopes and propose pedogenic secondary Fe and Al minerals may continue to retain P in these soils and limit biological productivity in the adjacent river even as active layer thickening increases potential P mobility in the watershed. We suggest that similar interactions may occur in other areas of the Arctic where comparable geochemical conditions prevail.

Sutor, Frederick W. [Univ. of Vermont, Burlington,↗

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗

Project Report: Undergraduate Student Research Program (USRP)

To better understand geochemical processes occurring on Europa's seafloor, we investigated the effects of varying Fe?content in the seafloor rock and varying temperature. Iron is important in such geochemical processes as the production of methane through serpentinization (e.g. Allen and Seyfried, 2003) and can be a nutrient for microbes (Russell and Hall, 2006; Park and Kim, 2001). It can also offer clues as to the state of differentiation of Europa's core/mantle. If Europa is fully differentiated and contains an iron core, we would expect there to be little iron in the mantle and ocean floor whereas a homogeneous Europa would have iron evenly dispersed throughout the ocean floor. Furthermore, the composition of the ocean is a result of water?rock interactions at the seafloor. This project investigated the effects of temperature on geochemical processes, comparing high temperature (> 250oC) hydrothermal vents (Kelley et al., 2001) to lower temperature (~20oC) cold seeps (e.g. Orphan et al., 2002).

Europa↗

Ultramafic Terranes and Associated Springs as Analogs for Mars and Early Earth

Putative extinct or extant Martian organisms, like their terrestrial counterparts, must adopt metabolic strategies based on the environments in which they live. In order for organisms to derive metabolic energy from the natural environment (Martian or terrestrial), a state of thermodynamic disequilibrium must exist. The most widespread environment of chemical disequilibrium on present-day Earth results from the interaction of mafic rocks of the ocean crust with liquid water. Such environments were even more pervasive and important on the Archean Earth due to increased geothermal heat flow and the absence of widespread continental crust formation. The composition of the lower crust and upper mantle of the Earth is essentially the-same as that of Mars, and the early histories of these two planets are similar. It follows that a knowledge of the mineralogy, water-rock chemistry and microbial ecology of Earth's oceanic crust could be of great value in devising a search strategy for evidence of past or present life on Mars. In some tectonic regimes, cross-sections of lower oceanic crust and upper mantle are exposed on land as so-called "ophiolite suites." Such is the case in the state of California (USA) as a result of its location adjacent to active plate margins. These mafic and ultramafic rocks contain numerous springs that offer an easily accessible field laboratory for studying water/rock interactions and the microbial communities that are supported by the resulting geochemical energy. A preliminary screen of Archaean biodiversity was conducted in a cold spring located in a presently serpentinizing ultramafic terrane. PCR and phylogenetic analysis of partial 16s rRNA, sequences were performed on water and sediment samples. Archaea of recent phylogenetic origin were detected with sequences nearly identical to those of organisms living in ultra-high pH lakes of Africa.

Blake, David↗

Convergence of micro-geochemistry and micro-geomechanics towards understanding proppant shale rock interaction: A Caney shale case study in southern Oklahoma, USA

As a direct outcome of economic development coupled with an increase in population, global energy demand will continue to rise in the coming decades. Although renewable energy sources are increasingly investigated for optimal production, the immediate needs require focus on energy sources that are currently available and reliable, with a minimal environmental impact; the efficient exploration and production of unconventional hydrocarbon resources is bridging the energy needs and energy aspirations, during the current energy transition period. The main challenges are related to the accurate quantification of the critical rock properties that influence production, their heterogeneity and the multiscale driven physico-chemical nature of rock–fluid interactions. A key feature of shale reservoirs is their low permeability due to dominating nanoporosity of the clay-rich matrix. As a means of producing these reservoirs in a cost-effective manner, a prerequisite is creation of hydraulic fracture networks capable of the highest level of continued conductivity. Fracturing fluid chemical design, formation brine geochemical composition, and rock mineralogy all contribute to swelling-induced conductivity damage. The Caney Shale is an organic-rich, often calcareous mudrock. Many studies have examined the impact that clay has on different kinds of shale productivity but there is currently no data reported on the Caney Shale in relation to horizontal drilling; all reported data on the Caney Shale is on vertical wells which are shallow, compared to an emerging play that is at double the depth. Here, in this work we develop geochemical–geomechanical integration of rock properties at micro-and nanoscales that can provide insights into the potential proppant embedment and its mitigation. The novel methodology amalgamates the following: computed X-ray tomography, scanning electron microscopy, energy dispersive spectroscopy, micro-indentation, and Raman spectroscopy techniques. Here, our results show that due to the multiscale heterogeneity in the Caney Shale, these geochemical and structural properties translate into a variation in mechanical properties that will impact interaction between the proppant and the host shale rock.

03 NATURAL GAS↗

The CORSAGE Programme: Continuous Orbital Remote Sensing of Archipelagic Geochemical Effects

Current and pending oceanographic remote sensing technology allows the conceptualization of a programme designed to investigate ocean island interactions that could induce short-term nearshore fluxes of particulate organic carbon and biogenic calcium carbonate from pelagic island archipelagoes. These events will influence the geochemistry of adjacent waters, particularly the marine carbon system. Justification and design are provided for a study that would combine oceanographic satellite remote sensing (visible and infrared radiometry, altimetry and scatterometry) with shore-based facilities. A programme incorporating the methodology outlined here would seek to identify the mechanisms that cause such events, assess their geochemical significance, and provide both analytical and predictive capabilities for observations on greater temporal and spatial scales.

Acker, J. G.↗

Applied genomics for industrially relevant microalgal strain development & optimization: characterizing genotype-phenotype relationships towards scalable environment-enhancing energy systems

Multiple planetary boundaries considered a safe operating space for humanity have been exceeded in recent years, with twice as many boundaries transgressed in 2023 as in 2009. Bioenergy systems are unique in that they can interact with and improve many of the transgressed boundaries directly, including multiple geochemical cycles, water and land use, and climate change. Among bioenergy systems, microalgae-based environment enhancing energy (E 2 -energy) are promising bioenergy systems for drop-in biofuels, valuable materials, chemicals and therapeutics, while making deep emission cuts and remediating wastewater, all without competing for agricultural resources.

09 BIOMASS FUELS↗

MEVTV Workshop on Nature and Composition of Surface Units on Mars

Topics addressed include: SNC meteorites and their potential for providing information about the geochemical evolution of Mars; remote sensing; photogeological inferences of Martian surface compositions; and interactions of the surface with volatiles in either the surface or the atmosphere.

Zimbelman, J. R.↗

Integrated experimental studies of pore structure and fluid uptake in the Bossier Shale in eastern Texas, USA

Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.

58 GEOSCIENCES↗

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗