Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface diffusion”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Elucidating the Initial Steps in α-Uranium Hydriding Using First-Principles Calculations

Hydrogen embrittlement of uranium, which arises due to the formation of a structurally weak pyrophoric hydride, poses a major safety risk in material applications. Previous experiments have shown that hydriding begins on the top or near the surface (i.e., subsurface) of α-uranium. However, the fundamental molecular-level mechanism of this process remains unknown. Here, in this work, starting from pristine α-U bulk and surfaces, we present a systematic investigation of possible mechanisms for the formation of metal hydride. Specifically, we address this problem by examining the individual steps of hydrogen embrittlement, including surface adsorption, subsurface absorption, and the interlayer diffusion of atomic hydrogen. Furthermore, by examining these processes across different facets, we highlight the importance of both (1) hydrogen monolayer coverage and (2) applied tensile strain on hydriding kinetics. Taken together, by studying previously overlooked phenomena, this study provides foundational insights into the initial steps of this overall complex process. We anticipate that this work will guide near-term future development of multiscale kinetic models for uranium hydriding and subsequently identify potential strategies to mitigate this undesired process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic and static corrosion of chromium-alumina refractory in simulated nuclear waste glass

Monofrax® K-3, a chromium–alumina refractory, is widely used in nuclear waste glass melters due to its high durability. However, aggressive vitrification conditions still lead to corrosion that limits melter lifespan. Here, this study investigates two key corrosion modes, subsurface and melt-line corrosion, under static and dynamic conditions, using a custom setup enabling precise control of melt velocity and temperature. Experimental results, interpreted using diffusion-limited dissolution models, show that (i) subsurface corrosion increases by approximately 10% per mm/s increase in melt velocity, while melt-line corrosion remains virtually unaffected, (ii) corrosion rate dependence on temperature is inversely proportional to melt viscosity, and (iii) glass composition affects both corrosion modes similarly, based on Cr 2 O 3 , Al 2 O 3 , and alkali contents. The experimental data also allowed the determination of diffusion coefficients of major dissolving species, which were found to be consistent between the melt-line and subsurface corrosion models, indicating their suitability for future CFD studies.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

Importance of the Antarctic Slope Current in the Southern Ocean Response to Ice Sheet Melt and Wind Stress Change

Abstract We use two coupled climate models, GFDL‐CM4 and GFDL‐ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre‐industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near‐surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large‐scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean‐state representation and meltwater‐driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.

Beadling, R. L.↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Lithium-Ion Battery Materials as Tunable, “Redox Non-Innocent” Catalyst Supports

The development of general strategies for the electronic tuning of a catalyst’s active site is an ongoing challenge in heterogeneous catalysis. To this end, herein, we describe the application of Li-ion battery cathode and anode materials as redox non-innocent catalyst supports that can be continuously modulated as a function of lithium intercalation. A zero-valent nickel complex was oxidatively grafted onto the surface of lithium manganese oxide (Li x Mn 2 O 4 ) to yield isolated Ni2+ occupying the vacant interstitial octahedral site in the Li diffusion channel on the surface and subsurface of the spinel structure (Ni/Li x Mn 2 O 4 ). Additionally, the activity of Ni/Li x Mn 2 O 4 for olefin hydrogenation, as a representative probe reaction, was found to increase monotonically as a function of support reductive lithiation. Simulation of Ni/Li x Mn 2 O 4 reveals the dramatic impact of surface redox states on the viability of the homolytic oxidative addition mechanism for H 2 activation. Catalyst control through support lithiation was extended to an organotantalum complex on Li x TiO 2 , demonstrating the generality of this phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Observation of molecular hydrogen produced from bridging hydroxyls on anatase TiO2(101)

Anatase TiO2 is used extensively in a wide range of catalytic and photocatalytic processes and is a promising catalyst for hydrogen production. Here, we show that molecular hydrogen was produced from bridging hydroxyls, HOb, on the (101) surface of single crystal anatase (a-TiO2(101)). Electron bombardment at 30 K produced bridging oxygen vacancies in the surface. Deuterated bridging hydroxyls, DOb, were subsequently formed via dissociation of adsorbed D2O and confirmed by infrared reflection-absorption spectroscopy. During temperature programmed desorption (TPD) spectroscopy, D2 desorption was observed at 520 K. Density functional theory calculations show that both H2 and H2O production from HOb are endothermic at 0 K on a-TiO2(101), but H2 (H2O) desorption is entropically driven above 230 K (800 K). The calculated activation barrier for H2 desorption is 1.40 eV, which is similar to the desorption energy obtained from analysis of the D2 TPD spectra. The H2 desorption likely proceeds in two steps: H atom diffusion on the surface and then recombination. The calculations indicate that diffusion of the hydrogen atoms into the subsurface could also be an important process. We thank Drs. Z. Dohnalek and R. Rousseau for stimulating discussions. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, & Biosciences and performed in EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle. We also acknowledge the WPI computational facilities for providing resources to perform the DFT studies.

Deskins, N A.↗
DOE OSTI2024

Directly resolving surface vs. lattice self-diffusion in iron at the nanoscale using in situ atom probe capabilities

Surface self-diffusion studies on metals under elevated reaction conditions are limited, as it is inherently challenging to unambiguously follow atomic transport across highly-reactive surfaces. Here, quantitative and mechanistic insight into thermally induced atomic transport processes in bcc α-iron at the sub-nanometer level was achieved using isotopic tracer techniques coupled with in situ atom probe tomography (APT) capabilities. Specifically, using a reactor directly connected to the APT, needle-shaped specimens fabricated from epitaxial thin films with an embedded 57 Fe tracer layer were annealed in Ar at 500 °C and 350 °C for 1 hour. Furthermore, the tracer was positioned at various depths in the APT specimen by field evaporation, enabling targeted and simultaneous analysis of lattice and surface diffusion. 57 Fe concentration profiles reveal lattice self-diffusion occurs at 500 °C on the order of ~7 – 9 monolayers, while lattice diffusion is not resolvable at 350 °C. Considerable surface transport was, however, observed at both conditions, where atomic transport over the specimen surface led to the formation of a thin (≤1 nm), isotopically-intermixed layer at the surface. Further, the observed isotopic redistributions at 500 °C were convoluted by additional processes occurring in the subsurface, such as atomic intermixing in correlation with lattice diffusion. However, surface diffusion was determined to be the primary transport process at 350 °C and was thereby quantified. Ultimately, these results demonstrate the significance of surface self-diffusion as a short circuit pathway. More broadly, this approach has the potential to provide detailed insight into (self-)diffusion mechanisms across various materials while targeting site-specific reactions under elevated reaction conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Soil gas probes for monitoring trace gas messengers of microbial activity

Abstract Soil microbes vigorously produce and consume gases that reflect active soil biogeochemical processes. Soil gas measurements are therefore a powerful tool to monitor microbial activity. Yet, the majority of soil gases lack non-disruptive subsurface measurement methods at spatiotemporal scales relevant to microbial processes and soil structure. To address this need, we developed a soil gas sampling system that uses novel diffusive soil probes and sample transfer approaches for high-resolution sampling from discrete subsurface regions. Probe sampling requires transferring soil gas samples to above-ground gas analyzers where concentrations and isotopologues are measured. Obtaining representative soil gas samples has historically required balancing disruption to soil gas composition with measurement frequency and analyzer volume demand. These considerations have limited attempts to quantify trace gas spatial concentration gradients and heterogeneity at scales relevant to the soil microbiome. Here, we describe our new flexible diffusive probe sampling system integrated with a modified, reduced volume trace gas analyzer and demonstrate its application for subsurface monitoring of biogeochemical cycling of nitrous oxide (N 2 O) and its site-specific isotopologues, methane, carbon dioxide, and nitric oxide in controlled soil columns. The sampling system observed reproducible responses of soil gas concentrations to manipulations of soil nutrients and redox state, providing a new window into the microbial response to these key environmental forcings. Using site-specific N 2 O isotopologues as indicators of microbial processes, we constrain the dynamics of in situ microbial activity. Unlocking trace gas messengers of microbial activity will complement -omics approaches, challenge subsurface models, and improve understanding of soil heterogeneity to disentangle interactive processes in the subsurface biome.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Utah FORGE: Geochemical Data for Cold Groundwaters and Produced Geothermal Fluids

Geochemical data for cold groundwaters and produced geothermal fluids around the Utah FORGE site. The data is compiled into four tables in the attached Excel File. Table 1 is a compilation of compositions (anions, cations, weak acids, oxygen, hydrogen, and carbon isotopes) for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 2 is a compilation of noble gas (He, Ne, Ar) and He and Ne isotopic compositions for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 3 provides values for calculated advective and diffusive fluxes of helium. Table 4 provides values of calculated subsurface stored heat between the Opal Mound fault and the Utah FORGE site, which are related to volumes of recently solidified magmatic heat sources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Time-Resolved Neutron Imaging for Hydrogen Uptake in Subsurface Lithologies

Geologic hydrogen production and underground storage are increasingly important for meeting rising energy demands while providing clean-combustion advantages. However, hydrogen’s high diffusivity and propensity for leakage through porous media necessitate direct evaluation of its transport behavior in subsurface materials. Whereas X-ray microcomputed tomography (μCT) studies often employ contrast agents or surrogate gases, this study leverages neutron transmission radiography/CT to observe hydrogen migration in situ. This work represents the first demonstration of real-time neutron radiography of hydrogen migration in reservoir and caprock lithologies. Cylindrical cores of Indiana limestone, Amherst Gray sandstone, and Tumey shale were subjected to constant-pressure hydrogen charging and scanned in real time using high-resolution neutron radiography. Results indicate immediate hydrogen infiltration in sandstone and limestone, with homogeneous distribution detected throughout their pore structure. In contrast, hydrogen remained largely absent from fine-grained shale under the same pressure, except in an apparently localized fracture zone, where neutron signatures confirmed the presence of hydrogen. Subsequent neutron CT of the sandstone sample, using image subtraction against an uncharged reference, corroborated hydrogen distribution patterns. Even under lowpressure, single-phase conditions, distinct neutron imaging signatures of hydrogen were achieved. These preliminary findings underscore the potential of neutron imaging for advancing subsurface hydrogen migration research.

08 HYDROGEN↗
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Electron Trapping and Ion Leaching at the Li-Modified Quartz–Water Interface

Electrons and ions at metal oxide-water interfaces have a critical role in many phenomena and applications, thus making their properties of considerable interest. We here investigate the behavior of a lithium impurity at the water interface with Li-doped SiO 2 using hybrid density functional based ab initio molecular dynamics simulations. We find that the excess electron donated by the lithium dopant localizes on a surface silicon atom and is partially solvated at the aqueous interface. While the position of the excess electron does not change during our ~60 ps simulation, the lithium cation diffuses from its initial position in the SiO 2 subsurface toward the interface and eventually leaches out of the oxide surface and becomes solvated by four water molecules, forming an aqueous Li+- electron complex at the interface. The degree of interaction between the localized electron and the lithium ion controls not only the energy level of the excess electron but also the extent of structural distortion and hydrogen-bonding at the interface. The results of our study can be relevant for understanding the initial stages of glass corrosion, where alkali ions are leached from the surface region by interactions with water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Phenethylammonium Functionalization Enhances Near-Surface Carrier Diffusion in Hybrid Perovskites

Understanding semiconductor surface properties and manipulating them chemically are critical for improving their performance in optoelectronic devices. Hybrid halide perovskites have emerged as an exciting class of highly efficient solar materials; however, their device performance could be limited by undesirable surface properties that impede carrier transport and induce recombination. Here we show that surface functionalization of methylammonium lead iodide (MAPbI 3 ) perovskite with phenethylammonium iodide (PEAI), a commonly employed spacer cation in two-dimensional perovskites, can enhance carrier diffusion in the near-surface regions and reduce defect density by more than one order of magnitude. Using transient transmission and reflection microscopy, we selectively imaged the transport of the carriers near the (001) surface and the bulk carriers in single-crystal MAPbI 3 microplates. The surface functionalization increases the diffusion coefficient of the carriers in a 40-nm subsurface region from ~ 0.6 cm 2 s -1 to ~ 1 cm 2 s -1 , similar to the value for bulk carriers. These results suggest the PEA ligands are effective in suppressing surface defect and phonon scattering and shed light on the mechanisms for enhancing photophysical properties and improving solar cell efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Iodine effective diffusion coefficients through volcanic rock: Influence of iodine speciation and rock geochemistry

Accurate prediction of the subsurface transport of iodine species is important for the assessment of long-term nuclear waste repository performance, as well as monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty, given that radioiodine decays into radioxenon. However, the transport of iodine through intact geologic media is not well understood, compromising our ability to assess risk associated with radioiodine migration. The current study's goal is to quantify the matrix diffusion of iodine species through saturated volcanic rock, with particular attention paid to the redox environment and potential speciation changes. Here, diffusion experiments were run for iodide through lithophysae-rich lava, lithophysae-poor lava, and welded tuff, whereas iodate diffusion was studied through welded tuff. Iodine transport was compared with a conservative tracer, HDO, and effective diffusion coefficients were calculated. Likely due to a combination of size and anion exclusion effects, iodine species diffused more slowly than the conservative tracer through all rock types tested. Furthermore, oxidation of iodide to iodate was observed in the lithophysae-poor lava, affecting transport. Results provide much needed data for subsurface transport models that predict radioiodine migration from underground sources, and indicate the pressing need for geochemical and redox interactions to be incorporated into these models.

58 GEOSCIENCES↗
DOE OSTI2023

Machine Learning Predictions of Simulated Self-Diffusion Coefficients for Bulk and Confined Pure Liquids

Diffusion properties of bulk fluids have been predicted using empirical expressions and machine learning (ML) models, suggesting that predictions of diffusion also should be possible for fluids in confined environments. The ability to quickly and accurately predict diffusion in porous materials would enable new discoveries and spur development in relevant technologies such as separations, catalysis, batteries, and subsurface applications. Here in this work, we apply artificial neural network (ANN) models to predict the simulated self-diffusion coefficients of real liquids in both bulk and pore environments. The training data sets were generated from molecular dynamics (MD) simulations of Lennard-Jones particles representing a diverse set of 14 molecules ranging from ammonia to dodecane over a range of liquid pressures and temperatures. Planar, cylindrical, and hexagonal pore models consisted of walls composed of carbon atoms. Our simple model for these liquids was primarily used to generate ANN training data, but the simulated self-diffusion coefficients of bulk liquids show excellent agreement with experimental diffusion coefficients. ANN models based on simple descriptors accurately reproduced the MD diffusion data for both bulk and confined liquids, including the trend of increased mobility in large pores relative to the corresponding bulk liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Probing the Surface Chemistry of Lithium Nitridation

Chemical synthesis of Li 3 N through lithium nitridation has potential to advance rechargeable battery and nitrogen fixation technology. However, studies of the conditions for forming Li 3 N on the lithium surface via nitrogen gas exposure report contradictory findings, such as the spontaneous reaction of Li with pure N 2 , the impossibility of forming Li 3 N through pure Li and N 2 interaction, the requirement of trace H 2 O to catalyze the reaction, and evidence to the contrary. In this study, ambient pressure X-ray photoelectron spectroscopy (APXPS) was applied to evaluate the in situ chemical evolution of the lithium metal surface under nitrogen gas up to 800 mTorr. At pressures ≤10 mTorr, no Li 3 N was detected. At higher pressures, surface Li 3 N rapidly reacts with trace CO 2 . Additionally, because metallic lithium is readily oxidized by trace gases, the atomic nitrogen concentration of the lithium surface remains below 2%. When nitridation follows oxidation by O 2 gas, CO 2 gas, or H 2 O vapor, surface Li 3 N formation is inhibited. These results suggest that nitrogen gas can diffuse through the oxidized lithium metal surface to react with subsurface metallic lithium.

Binding energy↗
DOE OSTI2024

Subsurface Redox Interactions Regulate Ebullitive Methane Flux in Heterogeneous Mississippi River Deltaic Wetland

Abstract As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive‐transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient‐driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

03 NATURAL GAS↗