Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pore size distribution”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗
DOE OSTI2024

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pore2Chip: All-in-one python tool for soil microstructure analysis and micromodel design

The Pore2Chip Python package is designed to create 2D micromodels using extracted data from 3D X-ray computed tomography (XCT) images. This package helps analyze soil structure and function, allowing for the investigation of hydro-biogeochemical processes that impact mineral extraction and reactivity, oxygen concentrations, and nutrient availability in disturbed or managed soils. Key metrics encompass pore size distributions, pore throat size distributions, and connectivity (pore coordination numbers). The final output is a 2D scalable SVG design representing a core or aggregate. Designs can be fabricated with methods such as laser etching, 3D printing, and photolithography.

lab-on-chip↗
DOE OSTI2026

Impact of Porous Transport Layer Morphology on the Performance of Proton Exchange Membrane Water Electrolyzers with Ultra-Low Iridium Loadings

Reducing Ir loadings in proton exchange membrane water electrolyzer anodes is critical for lowering capital expenses. Loading reduction could be achieved by improving the Ir activity via doping/alloying and/or the development of advanced microstructures. However, the anode porous transport layer (PTL) is a comparatively simple component whose properties also impact Ir utilization. Therefore, well-designed PTLs may also enable reduced Ir loadings. In this work, we survey eight PTLs from various manufacturers to observe their impact on cell performance at low (0.4 mg Ir cm -2 ) and ultralow (0.1 mg Ir cm -2 ) Ir loadings. The PTLs were characterized by their microstructural properties, including porosity, particle size distribution, and pore size distribution. Electrochemical cell performance was correlated to PTL morphology, and it was found that PTLs with lower porosities and smaller particle and pore radii enabled good performance even at ultralow Ir loadings. 1000-h durability testing indicated that using lower porosity PTLs can significantly improve durability behavior. A runaway voltage phenomenon was observed during durability testing of cells with ultralow Ir loadings, which was caused by increases in both anode and cathode overpotentials. Furthermore, we observed that the beginning of test performance of 0.1 mg Ir cm -2 cells correlates to the 1000-h degradation rates of 0.4 mg Ir cm -2 cells, suggesting that for the Ir catalyst used in this work, short-term testing at ultralow loadings can be used as an indicator of long-term degradation at higher loadings.

08 HYDROGEN↗
DOE OSTI2024

Thermally Insulating Transparent Barrier (THINNER) coatings on single pane windows

Conventional silica aerogel monoliths can provide remarkable thermal insulation but the presence of large pores (> 30 nm) tends to scatter visible light and render the material opaque or translucent instead of transparent. In addition, they are prone to cracking during synthesis and handling which makes them difficult to integrate in products and in particular in window solutions. This project developed two new solgel synthesis methods using silica precursors or preformed silica nanoparticles and ambient drying to produce mesoporous organo-silica monoliths. The monoliths were (i) thermally insulating, (ii) optically transparent, (iii) flexible, and (iv) hydrophobic. They feature porosity ranging from 50% to 90% with narrow pore size distribution with pore less than 20 nm resulting in excellent optical clarity (haze < 2%) and very low thermal conductivity (< 30 mW/mK). Interestingly, not only porosity but also pore size and mass fractal dimension were found to affect the thermal conductivity of the mesoporous silica. The superior transparency of the monoliths was shown to be attributed to dependent scattering among silica nanoparticles. Flexibility was achieved through trimethylchlorosilane surface modification. Furthermore, process scale-up and integration of the ambigel monoliths into window solutions using optically clear adhesives were also demonstrated in 6”x6” double-pane windows. The aerogel and the window solution were shown to be durable to accelerated aging under UV, moisture, and/or temperature gradient and thermal cycling. Overall, the window solution achieved the technical performance and the cost target of $10/sqft set for the SHIELD program. However, scaling to industry relevant scale (> 10’x10’) remains a challenge due to requirements on the drying process and fume hood size and to the propensity of large aerogel slab to crack during drying.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Insights into determining pore size properties of ultrafiltration membranes

The selectivity of porous membranes is often characterized using solute rejection tests, where membranes are challenged with dilute aqueous solutions of neutral solutes at operating conditions that minimize concentration polarization and fouling. In single solute tests, a membrane is challenged with one molecular weight (MW) solute at a time from low to high MW. Since single solute methods are time-intensive, mixed solute tests have become more common, where a mixture of several MW solutes challenges a membrane at once. However, the presence of large solutes in a mixture increases the rejection of smaller solutes. Furthermore, there are no universally accepted operating conditions or standard methods used by membrane manufacturers or researchers for the experiments, leading to difficulty in pore size and pore characteristic comparisons. In this paper, commercial ultrafiltration membranes were challenged with single and mixed solute polyethylene glycol (PEG) and dextran aqueous solutions. First, rejection values determined using total organic carbon (TOC) and high-performance liquid chromatography (HPLC) from single solute filtration experiments are compared. Differences in rejection curves obtained by the two techniques are attributed to solute polydispersity. Mixed solute filtration experiments with binary mixtures of solutes showcased co-solute interactions, which increase with both the size and weight percent of large solute in the mixture. Mixed solute filtration experiments at varying operating conditions (i.e., stir speed and flux) were conducted to determine operating conditions that mitigate co-solute interactions. Stir speed had a minimal effect on co-solute interactions. In contrast, low flux conditions can help minimize co-solute interactions, leading to pore size distributions that closely resemble results observed in single solute filtration using narrowly dispersed solutes. Additionally, at low flux conditions, the predicted membrane pore size distributions utilizing mixed solute experiments with PEG and dextran were similar.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A User-Friendly GUI Tool for Automated Microstructural Analysis of Fiber-Reinforced Composites and Porous Structures

Understanding and quantifying microstructural features such as fiber orientation and porosity is critical for predicting the mechanical behavior and performance of fiber-reinforced polymer composites. Traditional manual analysis is time-consuming, subjective, and unsuitable for high-throughput datasets. We present a graphical user interface (GUI) application that automates the analysis of microscopy images to extract key microstructural metrics, including fiber orientation tensors, fiber orientation distribution, porosity and pore size distribution. The app integrates multiple image segmentation techniques including global and local thresholding, clustering, and region-based approaches, offering flexibility for different types of image qualities and features. Users can load microstructural images, select regions of interest and segmentation techniques tailored to their image dataset. It also addresses a critical challenge in fiber orientation analysis: the ambiguities caused by touching, overlapping, or partially cut fibers. It supports autorun examples for standardized workflows, enabling reproducible analysis and facilitating training and benchmarking. This tool significantly reduces manual intervention, enhances consistency, and accelerates data generation for structure–property modeling, process optimization, and digital materials research. The tool is intended for use by materials scientists, engineers, and researchers engaged in composite characterization, quality control, and machine learning-based microstructural studies.

Chawla, Komal [ORNL] (ORCID:0000000190327565)↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗
DOE OSTI2025

Unravelling the origins of shale nanoporosity using small-angle neutron scattering (SANS)

Hydrocarbon production from tight rocks is constrained by slow diffusion within the shale matrix, limited by small pore sizes and low permeability. The nanopore proportion and size distribution significantly influence matrix permeability, a key property for optimizing hydrocarbon recovery and supporting hydrogen production while minimizing environmental impacts. Small-angle neutron scattering (SANS) has been an important tool for exploring the characteristics and structure of shale nanopores. This study used SANS to analyze nanoporosity and pore size distribution (<100 nm) in tight rocks with varying compositions to determine the influence of rock heterogeneity on SANS measurements. Results showed that nanoporosity correlates with clay content, with the highest clay-rich shale (52.48 wt% clay) exhibiting 8.8 % nanoporosity. SANS also revealed more nanopores than traditional nitrogen adsorption measurements, affirming its ability to reflect bulk mineralogy and upholding the relevance of experimental findings using this technique to optimize field operational approaches.

03 NATURAL GAS↗
DOE OSTI2024

Fabrication of model ultrafiltration membranes with uniform, high aspect ratio pores

Here, in this manuscript, we report the facile fabrication of large-area model membranes with highly uniform and high aspect ratio pores with diameters <20 nm. These membranes are useful for fundamental investigations of separation by size exclusion in the ultrafiltration regime, where species to be separated from solution have dimensions of 1–100 nm. Such investigations require membranes with narrow pores and high aspect ratios such that the Hagen–Poiseuille equation is followed, enabling well-known models such as the hindered transport model to be evaluated and other affecting factors to be ignored. We demonstrate that the sub-20 nm pores in the membrane are of sufficiently high aspect ratio such that water flux through the membrane is consistent with the Hagen–Poiseuille equation. The fabrication relies on self-assembling block copolymers to form uniform, densely packed patterns with sub-20 nm resolution, sequential infiltration synthesis to convert the block copolymer in situ into a mask with adequate contrast to etch pores with an aspect ratio >5, and low-resolution photolithography to transfer the pattern over a large area into a silicon nitride membrane. Model membranes with narrow pore-size distribution fabricated in this way provide the means to investigate parameters that impact size-selective ultrafiltration separations such as the relationships between solute or particle size and pore size, their distributions, and rejection profiles, and, therefore, test the validity or limits of separation models.

Chen, Wen↗
DOE OSTI2025

3D pore structure reconstruction of fluff region in top of EBR-II irradiated fuel by X-ray micro-computed tomography

Metallic fuels comprised of U-Zr binary and U-Pu-Zr ternary alloys are a promising fuel candidate for sodium-cooled fast reactors and possess many advantages over ceramic fuels such as higher uranium density and improved thermal conductivity. Fluff is a highly porous structure forming at the top of metallic fuel slugs during irradiation and was extensively observed in past sodium-cooled fast reactor (SFR) experiments such as Idaho National Laboratory’s Experimental Breeder Reactor II (EBR-II). Previous investigations have shown that fluff is rich in fissile atoms and hence could have neutronics implications, but its formation mechanisms are not understood. In this work, two cube-shaped lift outs were taken from the fluff and bulk fuel regions of a high burnup U-19Pu-10Zr composition EBR-II fuel slug and characterized by micro computed tomography (Micro-CT) and energy dispersive x-ray spectroscopy (EDS). This was used to characterize differences in microscale porosity and composition between fluff and bulk fuel, and to support or oppose postulated fluff formation mechanisms. Surface EDS showed minimal differences between the elemental compositions of the fluff and fuel cubes. However, three-dimensional (3D) pore analysis from the tomography scans showed significant differences between the fluff and fuel cube pore size distributions at the microscale level. The high porosity and irregular pore distribution seen in macroscale fluff was also present within the micro-scale fluff cube, having a much higher porosity than the fuel cube which was concentrated in a few large pore networks. In conclusion, these results indicate support for a creep driven mechanism of fluff formation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

SAN-Based Block Polymers as a Platform for Manufacturing Strong Isoporous Membranes

Ultrafiltration (UF) membranes are ubiquitous in water purification and bioprocessing. However, co-designing their mechanical and transport properties remains challenging because of the broad pore size distributions at the surface and within the bulk that result from nonsolvent-induced phase separation (NIPS) – their typical manufacturing process. These distributions influence the hydrodynamic resistance to water flow and the stress concentrations around the pores. Developing advanced UF membranes requires innovative molecular designs that offer control over the surface and bulk pores, as well as the mechanical properties of the load-bearing, polymer. Here, we introduce a platform for designing UF membranes by leveraging solution self-assembly of block polymers and chain architectures with pendant polar groups. The block polymers consist of a poly(styrene-co-acrylonitrile) hydrophobic block, which is known for its strength, and a poly(4-vinyl pyridine) hydrophilic block, which drives solution self-assembly. We focus on a series of block polymers with constant molecular weight, M n ≈ 115 kDa, SAN fraction, 75 wt.%, and varying acrylonitrile content, 0 to 40 mol%, to demonstrate that: (i) RAFT dispersion copolymerization of acrylonitrile and styrene provides a facile route to synthesize strong block polymers, (ii) incorporation of acrylonitrile into the hydrophobic block enhances membrane strength by facilitating chain entanglements and dipole-dipole interactions, and (iii) acrylonitrile alters the balance between membrane permeance and rejection, even when the membranes feature similar surface and bulk pores. Overall, our results provide insights into the molecular design of UF membranes with enhanced mechanical and separation properties, contributing to the development of materials for water and energy technologies.

deformation↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2025

Miscanthus × giganteus increases soil maximum water holding capacity compared to maize

Soil ecosystem services, like the ability to store water, have been depleted after a century of conventional, annual cropping, and perennial crops offer a solution to this and other agricultural environmental issues. We assessed the impact of Miscanthus × giganteus (miscanthus), a perennial biomass crop, on soil water holding capacity and structure compared to continuous maize (Zea mays L.) at two sites in Iowa. After three growing seasons, we measured the following: (1) maximum water holding capacity (MWHC) with and without soil structure, and (2) total porosity and pore size distribution (PSD) via micro-computed tomography (microCT). Miscanthus increased MWHC by 14.7% across both sites relative to maize (p = 0.002), and we attributed this to structural changes due to the lack of a crop effect when measured on structureless soils. No significant changes were detected in soil organic matter, texture, total porosity, or PSD that could explain the increase in MWHC under miscanthus. Our findings suggest that the increases in MWHC are primarily due to structural changes rather than increases in soil organic matter or porosity (at least porosity detectable by microCT). This study highlights miscanthus' potential to enhance soil water storage and underscores the need for further investigation to clarify the mechanisms through which this biomass crop influences soil structural properties.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Alternating salt and freshwater floods of coastal soils impact soil structure, hydraulic properties, and oxygen dynamics

Coastal soils are increasingly impacted by hydrologic intensification in the form of rising sea level and flooding from storm surges and precipitation. Alternating saltwater (SW) and freshwater (FW) exposure has the potential to disperse colloids, which can lead to disintegration of soil structure, clogging of pore spaces, and reduction in ecologically and biogeochemically important functions like infiltration or gas exchange. To investigate how hydrologic intensification affects soil structure and oxygen dynamics, we conducted a series of laboratory-based flood simulations. Intact soil cores of the A and B horizons (22 total) from the toe slope of an upland coastal forest along the western shore of Chesapeake Bay. We subjected the cores to 24 h of saturation with either FW or alternating brackish SW and FW, with a 24-h draining event in between flood events. Oxygen diffusion into soil during draining was reduced by up to 30% for soils flooded alternating SW-FW compared to soils flooded only with FW. We attributed the reduced oxygen diffusion to clogging of smaller pores by colloids, with colloid redistribution observed as an increase in specific surface area down the core profile of up to 59%. After three SW-FW floods (six floods total), there were significant changes in pore size distribution, significant redistribution of colloids, and the A horizon became sodic. We concluded that a small number of SW flooding events can induce a measurable change in soil physical properties that directly impacts the biogeochemical dynamics.

Rod, Kenton A. [Pacific Northwest National Laborat↗
DOE OSTI2024

Wetting mechanism and alteration of nano-sized shale pores: Insights from contrast variation small angle neutron scattering

Wettability of tight shale is crucial for fluid flow and mass transport process in energy geosciences. However, understanding the interfacial chemistry and wetting mechanisms at sub-nano-pore scales remains a formidable challenge. Here, in this study, the Contrast Variation technique of Small Angle Neutron Scattering (CV-SANS) is employed to investigate shale’s interfacial chemistry using reagents that possess a range of different polarities, including water, n-decane, toluene, and dimethyl methanamide. Through five different experimental strategies, we have demonstrated a successful modification of shale wettability, ranging from enhancement, weakening, to reversal. Delving into the mechanisms, we illustrated the crucial role of pre-existing liquid films in these changes, where the uniquely co-existing polar and non-polar functional groups in dimethyl methanamide acted as a conduit for interfacial chemistry adjustments. Furthermore, a solvent immersion led to matrix dilation as well as liberation of residual oil-occupied pores, resulting in altered pore size distributions, with hydrogen bonding playing a significant role in the polar groups. Interestingly, despite shale exhibiting a stronger affinity for oil over water, hydrophilic solvents induced more substantial dilation than lipophilic ones. Collectively, this work elucidates the dynamic change of interfacial chemistry via the configuration of polarity using chemical reagents, and the CV-SANS technique underscores its invaluable utilities in decoding the interfacial wettability traits in nanopore space of shale.

58 GEOSCIENCES↗