Search NASA⌕ Search

SEARCH · Search NASA

Results for “shale gas”

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

Pressure-stable supported ionic liquid membranes using isoporous supports for evaluating pure- and mixed-gas light paraffin fractionation

Advances in horizontal drilling and hydraulic fracturing have spurred the growth of domestic U.S. energy production. Membranes, typically silicone rubbers, have found utility in shale gas treatment as fuel gas conditioning units to selectively remove C 2 + hydrocarbons at pressures up to 30 bar, producing clean CH 4 for gas engines. However, more selective materials could be beneficial for broader shale gas treatment applications, such as dew point control units. Supported ionic liquid membranes (SILMs) offer a potential opportunity for improving C 3 H 8 /CH 4 selectivity, but they lack pressure-stability. Here, we report C 3 H 8 /CH 4 selective and pressure-stable SILMs using isoporous supports. SILMs with supports that had minimal defects remained stable up to 15 bar of transmembrane pressure. This stability allowed for pure- and mixed-gas testing of the resulting membranes at elevated pressures. Furthermore, these tests revealed that low viscosity ILs (<100 cp) may display mixed-gas permeances and permselectivity nearly identical to pure-gas results. On the other hand, higher viscosity ILs may display increasing permeances and permselectivity with increasing C 3 H 8 fugacity, similar to rubbery polymers. Ultimately, the SILMs demonstrated relatively high pressure-stability due to the isoporous supports and competitive mixed-gas C 3 H 8 /CH 4 permselectivity compared to silicone rubber.

Rosenthal, Justin J. [University of Texas at Austi↗
DOE OSTI2025

Characterization of kerogen nanopores using 2D NMR relaxation and MD simulations

The characterization of kerogen nanopores is crucial for predicting the geostorage capacity and recoverability of natural gas in unconventional gas shale reservoirs. Towards this end, a powerful technique is presented which integrates 2D NMR T 1 -T 2 relaxation measurements with molecular dynamics (MD) simulations of hydrocarbons confined in the nanopores of kerogen. The integrated NMR-MD technique is demonstrated using T 1 -T 2 measurements of kerogen isolates and organic-rich chalks saturated with heptane, together with MD simulations of heptane completely dissolved in a realistic kerogen model. The NMR-MD results are used to extract the swelling ratio and nanopore size distribution of kerogen as a function of depth in the reservoir. The effects of organic nanoconfinement on the T 1 relaxation dispersion and T 2 residual dipolar coupling of heptane are investigated, as well as the effect of downhole effective stress on the kerogen nanopore size as a function of depth and compaction. Potential applications in partially depleted gas shale reservoirs are discussed, including CO 2 utilization/geostorage, geostorage of green H 2 , and integration of the NMR-MD technique with thermodynamic models for predicting the competitive sorption of gas mixtures in kerogen.

Compaction↗
DOE OSTI2025

Alaska Liquid Natural Gas Pipeline Front-End Engineering & Design (Final Technical Report)

The Alaska Gasline Development Corporation (AGDC) is Alaska’s natural gas infrastructure development corporation established in 2013. AGDC’s mission is to maximize the benefit of Alaska’s vast North Slope natural gas resources for Alaskans through the development of infrastructure necessary to move the gas into local and international markets. AGDC was identified for a Congressionally Directed Spending (CDS) project for funding in the Energy and Water Development and Related Agencies Appropriations Act, 2023 under the heading: “Congressionally Directed Energy Efficiency and Renewable Energy Projects.” The CDS included $\$$4,000,000 of direct funding, with required match funds, to move the project forward. Alaska’s North Slope holds America’s largest proven and conventional natural gas supply. The integrated Alaska LNG Project will deliver 3.5 billion cubic feet of natural gas per day from Alaska’s North Slope gas fields to Alaskans as well as to a marine terminal located at tidewater in Cook Inlet. Alaska LNG is an integrated gas infrastructure project with three major components: a gas treatment plant (GTP) located at Prudhoe Bay, an 807-mile (1,287 km) gas pipeline (Mainline Pipeline) to Southcentral Alaska with interconnections for in-state gas use, and a natural gas liquefaction facility (LNG Facility) in Nikiski, Alaska. The integrated Alaska LNG Project has several strategic advantages including proven gas resources, existing upstream infrastructure, an advantageous arctic climate for LNG production, proximity to LNG markets, a track record of reliability from a state that first began exporting LNG to Japan in 1969, and broad support from Alaskans. North Slope natural gas is a conventional resource and can be produced with minimal drilling at a fraction of the carbon dioxide emissions of shale gas from the Lower 48 states. Through the development of the Alaska LNG Project, Alaska can provide energy security to Alaskans and a stable source of LNG to the Asia-Pacific region for generations. The Alaska LNG Project has been progressed through Pre-Front-End Engineering Design (Pre-FEED) and has obtained all major federal and State of Alaska permits and authorizations to construct the project, including the Federal Energy Regulatory Commission (FERC) Order Granting Authorization Under Section 3 of the Natural Gas Act. On September 5, 2024, the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) awarded Project No. DE-FE0032307 to AGDC with the objective to progress the project to Front-End Engineering Design (FEED) entry for the Alaska LNG Project Phase 1 Pipeline. The award Start Date was made effective July 1, 2023, with a Period of Performance through June 30, 2025. On March 27, 2025, AGDC announced the execution of definitive commercial agreements with Glenfarne Alaska LNG, LLC, an affiliate of Glenfarne Group, LLC, (together as “Glenfarne”), to lead the development of the Alaska LNG Project and enter FEED for the Phase 1 Pipeline. Project activities are now funded and directed by this private sector partner who holds a 75% interest in 8 Star Alaska, LLC (8 Star). 8 Star holds the assets of the Alaska LNG Project. As planned, AGDC continues to hold 25% minority interest in 8 Star and will play a governance role moving forward with Alaska LNG. This definitive commercial agreement milestone led to the successful completion of AGDC’s Statement of Project Objectives (SOPO) for FEED entry and led to the completion of DOE Project No. DE-FE0032307. At conclusion of the SOPO, AGDC also reached the award’s maximum federal cost share of $\$$4,000,000. AGDC is, therefore, providing Final Technical Report to close out DOE Project No. DE-FE0032307.

02 PETROLEUM↗
DOE OSTI2025

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
DOE OSTI2025

A review on advances in oxidative coupling of methane (OCM) for industrial use and prospects of CO 2 –H 2 O splitting integration

The discovery of shale gas reserves has encouraged the development of direct methods for methane conversion into valuable chemicals, offering an alternative to indirect approaches that involve an energy-intensive and intermittent syngas production step, leading to high CO 2 emissions. Amongst the direct methods, the oxidative coupling of methane (OCM) is a potential pathway to reduce CO 2 emissions and can produce commodity chemicals such as ethylene, a chemical regarded as central to the petrochemical industry. Even though OCM has been studied for over four decades, the technology still has not found commercial application. Amongst the challenges regarding industrial deployment of OCM, the most significant one is the requirement of a high ethylene yield of 30 % which is currently reported to be around 20 %. Moreover, the highly exothermic nature of the process and controlling the carbon selectivity over oxides of carbon (COx) is the heart of the problem. Numerous researchers have presented promising results in terms of catalysts, reactor designs and feeding strategies for OCM. However, due to lack of inclusiveness in the results, none of the combination of catalysts, reactors and system optimizations has been able to bring about its industrial viability. The current paper presents an extensive review of the noteworthy attempts to achieve industrial targets for OCM. Moreover, a comprehensive criteria is presented which highlights the desired end state for the industrial deployment of OCM technology. Furthermore, the criteria is based on literature survey and a comparison with industrially deployed ethylene production plants i.e., naphtha or ethane steam cracker plants. Finally, a novel integration technology is presented which includes a combination of OCM and CO 2 -H 2 O splitting in a chemical looping reactor design to enable efficient energy utilization and minimal heat losses to the environment.

CO2 Splitting↗
DOE OSTI2025

Autocorrelation and Multifractal Detrended Fluctuation Analyses Reveal Superdiffusive Mass Transport in Solvent-Filled Nanoporous Media

Fluorescence fluctuation spectroscopy experiments were conducted to better understand the complex mass transport dynamics of organic molecules in liquid-filled nanoporous media. Anodic aluminum oxide (AAO) membranes incorporating 10 and 20 nm diameter cylindrical pores were employed as model materials. Nile red (NR) dye was used as a fluorescent tracer. The dye was dissolved separately in ethanol and toluene at a concentration of 20 nM and used to fill the membrane nanopores. Confocal fluorescence microscopy was employed to capture photon intensity time series data reflecting apparent diffusion of the dye within the pores. Autocorrelation of these data revealed that NR diffusion within the membranes occurred over a broad range of time scales. The autocorrelation decays were fit to a model for one-dimensional diffusion incorporating both fast and slow components having apparent diffusion coefficients, D f and D s , differing by a factor of ∼100. The fast mechanism was attributed to hindered bulk-like diffusion in the central pore cavity, while slow diffusion likely involved absorption of the dye to the pore surfaces. Unfortunately, important evidence of diffusion anomalies is lost in the broad autocorrelation decays obtained. The method of multifractal detrended fluctuation analysis (MF-DFA) was applied to the same data as a means to overcome this limitation. MF-DFA revealed that time series acquired from within the nanopores were multifractal and exhibited evidence of anomalous superdiffusion, likely resulting from the participation of a desorption-mediated diffusion mechanism. Monte Carlo simulations of time series modeling desorption-mediated diffusion in cylindrical nanopores provided support for this assignment. Here, the new knowledge gained affords an improved understanding of hydrocarbon dynamics within nanoporous oil and gas shales.

Diffusion↗
DOE OSTI2026

A Robust Methodology to Elucidate Kinetics of Room Temperature Electrochemical Propane Adsorption on Platinum

Electrocatalytic activation of alkanes can further decarbonize chemical manufacturing by leveraging affordable renewable electricity and readily available shale gas reserves in the United States. Earlier works have identified the unique role of Pt in adsorbing and activating alkanes, like propane, at room temperature in acidic, aqueous electrolytes, revealing spontaneous formation of deeply dehydrogenated propane-derived surface species with an intact C 3 - backbone. Although an adsorption mechanism was hypothesized, it has not been explicitly investigated to date, preventing the quantification of kinetic rate parameters. A robust methodology to investigate and benchmark propane adsorption kinetics on Pt is critical for the rational design of electrocatalysts that exhibit higher selectivity toward desired partially oxidized products. Herein, we analyze an oxidative current transience that appears during the adsorption of propane on Pt in aqueous electrochemical conditions and develop a methodology that elucidates the adsorption mechanism and enables quantification of rate parameters such as order dependences and apparent activation barriers. This method yields an expected first-order dependence with respect to propane concentration at low coverage and reveals a second-order dependence with respect to the concentration of surface active sites. Additionally, the apparent activation barrier for propane adsorption was calculated using an Arrhenius analysis of the current transience under temperature control. The experimentally measured activation barrier of 35 kJ mol –1 is in excellent agreement with the theoretical barrier calculated by density functional theory (DFT). The kinetic analysis was extended, via the use of transition state theory, to extract entropy and enthalpy of activation, yielding consistent results with the proposed two-step adsorption mechanism and DFT calculations. These results demonstrate reliable quantification of kinetic parameters for electrocatalytic activation of C–H bonds in alkanes that can be employed for rational catalyst development for a versatile range of electrocatalytic conditions.

alkane activation↗
DOE OSTI2026

Reacting CO 2 with Light Alkanes to Value-Added Products

Catalytic conversion of anthropogenic carbon dioxide (CO 2 ) into value-added products is a promising strategy to mitigate global carbon emissions. Concurrently, the shale gas revolution has provided an abundant supply of light alkanes (methane, ethane, propane, and butane), presenting a unique opportunity to employ these underutilized hydrocarbons as an effective, low-cost hydrogen source for CO 2 reduction. In this Perspective, we summarize past efforts, current state, and future opportunities for reacting CO 2 with light alkanes to generate a diverse range of value-added products. Compared with direct alkane conversion, the introduction of CO 2 fundamentally alters reaction thermodynamics and kinetics, enabling selective C–H and C–C bond activation while suppressing catalyst deactivation from coke formation. Building on decades of research in dry reforming and CO 2 -assisted dehydrogenation, recent advances in catalyst design have enabled CO 2 -assisted dehydrogenation processes that approach chemical equilibrium for the selective production of olefins and syngas. Importantly, advances in catalyst design and reactor engineering have further expanded the product scope beyond gas-phase (syngas and olefins) to include liquid-phase (oxygenates and aromatics), and solid-phase products (carbon nanomaterials). We highlight key catalyst design principles for controlling reaction pathways and discuss major challenges and opportunities in developing selective and versatile platforms for the simultaneous upgrading of CO 2 and light alkanes.

CO2↗
DOE OSTI2025

NETL Natural Gas Lifecycle Model 2020 Rev1

This is the excel-based life cycle model that contains all the parameters and Monte Carlo simulation capabilities to model the techno-regions contained in the report: Life Cycle Analysis of Natural Gas Extraction and Power Generation: U.S. 2020 Emissions Profile. This is an update to a previous version (https://edx.netl.doe.gov/dataset/netl-natural-gas-lifecycle-model) and corrects a modeling error in the post-processing natural gas composition.

Calcium Chloride↗
DOE OSTI2025

Predicting Dynamic-to-Static Correction Factor from Petrophysical Data and Chemostratigraphy using Unsupervised Machine Learning

Estimating static mechanical properties of stratigraphic layers is critical for optimizing subsurface engineering applications. To estimate dynamic-to-static correction factor F ds (static-to-dynamic Young’s modulus ratio) across the Caney shale interval in Oklahoma, USA, we integrated triaxial test measurements and petrophysical data, including well logs and X-ray fluorescence (XRF) using unsupervised machine learning (ML). We used a novel workflow that includes principal component analysis (PCA) to reduce data set dimensionality of well logs and XRF data sets—both separately and combined—creating three scenarios, and later applied inverse distance weighting (IDW) to derive F ds profiles for these scenarios. Furthermore, we applied K-means clustering on each scenario to predict depositional facies, and built a stiffness zonation profile through chemostratigraphic analysis of the terrigenous elements to validate the predicted F ds . The predicted F ds profile from each scenario using the PCA-IDW method was compared with the constant F ds approach from our previous study by calculating the root mean square error (RMSE). The combined data sets scenario yielded the lowest RMSE value of 0.113, while the RMSE values for the well logs and XRF scenarios were 0.131 and 0.129, respectively. In addition, the predicted F ds from the XRF scenario well-matched the stiffness zonation from the chemostratigraphic analysis that was built using the optimized K-means clustering of nine clusters for that scenario. These methods and findings offer a valuable tool for refining lithological classification and improving the F ds profile, potentially enhancing drilling and stimulation strategies for subsurface energy engineering applications.

clastic rock↗
DOE OSTI2025

Exploring Catalyst Compositions for Microwave-Assisted Methane Dehydroaromatization

The flaring of natural gas in U.S. shale regions remains a challenge for producers. One alternative to flaring is converting the wasted gas into valuable chemicals. Microwave-based processes offer a promising solution, potentially enabling the development of compact, modular systems for on-site production of chemicals, such as aromatics, from natural gas. This is due to the advantages of microwave heating, including efficient heating of compact volumes, accelerated reaction rates, and electrification of heating. However, developing effective catalysts for microwave-based processes is challenging, as conventional materials often require modification to be effectively heated by microwaves. This study provides an overview of a catalyst development project focused on a molybdenum-supported zeolite catalyst, optimized for the direct conversion of methane into aromatics under microwave irradiation. It details the synthesis, characterization, and the effects of promoters, as well as computational efforts undertaken to understand and enhance the catalyst's performance.

flare gas↗
DOE OSTI2025

Marcellus assessment shows continued expansion

Increasing uncertainty regarding the remaining volume of recoverable gas in Appalachian shale plays provided the original motivation for this study. Assessments of the Marcellus and Utica shale plays published by the US Geological Survey (USGS) in 2019 established resource volumes dramatically lower than ongoing well performance would indicate. Further, as production in the Marcellus increased and wells matured, operators observed that wells were producing more than the estimated gas-in-place (GIP).

03 NATURAL GAS↗
DOE OSTI2025

Kinetic and Mechanistic Understanding of Boron-Containing Catalysts for the Production of Light Olefins through Selective Oxidation

In recent years, an abundance of natural gas obtained from shale deposits has created opportunities for the US chemical industry to establish new and more efficient processes for the production of chemical intermediates and polymers. For example, the replacement of petroleum-derived naphtha as a to shale-derived ethane as feedstocks has made the production of ethylene/polyethylene less costly. However, this shift in feedstock has led to lower production of propylene from steam crackers, and new “on-purpose” light olefin production technologies have been implemented to bridge the gap between propylene supply and demand. The most important of these processes is direct propane dehydrogenation (PDH), but catalyst deactivation and high temperatures lead to high capital and operating costs. Oxidative dehydrogenation of propane (ODHP) is an efficient alternative to PDH, but high selectivity to CO2 from previous metal oxide catalyst prevents implementation. Hexagonal boron nitride (hBN) and other boron-containing materials have been recently developed as highly selective catalysts for both ODHP and oxidative cracking, surpassing previous systems. Much is still unknown in regard to the mechanism and the dynamic surface changes that occur under reaction conditions. This research project aims to combine in situ/operando spectroscopic investigation and kinetic studies to elucidate the interplay between the proposed gas-phase mechanism and the surface structures on the working catalyst. These insights will be leveraged to create a framework for optimizing reaction conditions and material synthesis for the next generation of boron-containing selective oxidation catalysts.

02 PETROLEUM↗
DOE OSTI2026

Effect of Mo precursors in Microwave-assisted Methane Dehydroaromatization over Mo/HZSM5 catalysts

Natural gas flaring occurs in remote shale regions due to limited pipeline takeaway capacity. The conversion of the associated natural gas into aromatics in modular microwave reactors is a viable alternative to monetize the wasted gas. Microwaves offer rapid, selective heating in compact reactor systems that can enable on-demand chemical production at the well-site. Mo-HZSM-5 catalysts are widely used for aromatic production, but the location and nature of the active sites are still under debate. This study focuses on the use of 6 different Mo precursors to elucidate insights into the Mo properties that are desirable for BTX production. The catalysts were characterized by different methods (XPS, TPR, Raman, etc) to determine differences in Mo catalytic properties that may affect performance and understand these differences through performance testing under microwave at 700C for methane dehydroaromatization. Metal precursors that enable a better distribution of Mo into the pores and over the surface lead to improved benzene yield, whereas those that limited Mo to primarily the surface suffer rapid deactivation and low benzene yields. Additionally, strong Lewis acidity that arises from the sodium containing precursor drastically shifts the product selectivity towards dehydrogenation, which produces more ethylene and carbon. This catalyst had the highest deactivation constant of all catalysts tested.

gas flaring reduction↗
DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- 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. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗
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↗