Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geothermal”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Optimization of fluorescence and surface adsorption of citric acid/ethanolamine carbon nanoparticles for subsurface tracers

The tunable fluorescence and biocompatibility of fluorescent carbon nanoparticles (FCNs) makes them appealing tracers in geothermal, environmental, and clinical applications. Here, we address FCNs synthesized from reaction and pyrolysis of citric acid and ethanolamine. We examine the fluorescence and the adsorption of these FCNs as critical parameters that influence their use in applications. FCN fluorescence was observed to develop during pyrolysis and measured at pyrolysis temperatures ranging from 190 to 250 °C. Transmission electron microscopy shows that the FCN average diameter of approximately 13 nm is unaffected by pyrolysis temperature. However, pyrolysis temperature strongly affects fluorescence, which peaks at 210 °C. Here, the surface hydrophobicity of the FCNs increases with pyrolysis temperature as measured by retention on reverse phase chromatography. Together, these results suggest that particle surface functionalities provide both fluorescence and hydrophilicity, and that these surface functionalities are destroyed with excessive pyrolysis. To examine their surface adsorption, quartz crystal microbalance measurements demonstrate that pyrolysis temperature increases FCN adsorption on mineral surfaces. Surfaces with greater hydrophobicity show greater FCN adsorption. Adsorption is unaffected by the presence of mono- or divalent ions in solution. This suggests adsorption by hydrophobic interaction for the Fe 2 O 3 , SiO 2 , and AlSiO surfaces studied.

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

Influence of extreme temperature conditions on CO 2 direct air capture using amino-acid solutions

Geological features play a pivotal role in determining the feasibility of deploying CO₂ direct air capture (DAC) technologies, primarily because they influence the availability of cost-effective energy sources, such as natural gas and geothermal energy, and also due to the potential for CO₂ sequestration. Many regions face challenges due to variable weather conditions including seasonal temperature fluctuations, high or low humidity, and sub-ambient temperatures. These extremes can reduce DAC performance or even lead to catastrophic events. Aqueous solvents considered for DAC systems are particularly vulnerable to seasonal variations in colder climates, where the solvent may underperform or freeze. It is therefore essential to investigate the CO₂ capture efficiency of aqueous solvents across a broad range of environmental temperatures, spanning sub-zero to hot conditions (>30 °C). In this study, DAC operation is examined using a high-flux solvent–air crossflow contactor under two major weather scenarios: (i) cold conditions below 0 °C and (ii) hot conditions above 30 °C. A parametric study is conducted to investigate the contactor performance regarding CO₂ removal efficiency, uptake capacity, and reaction kinetics versus temperature when the air velocity through the contactor exceeds 1 m/s. The efficacy of the contactor is systematically investigated using various anti-freeze amino-acid solvent formulations. A mass-transfer mechanistic model is developed to assess the process performance over a wide temperature range and propose scalable design guidelines. Machine learning is also employed to identify key parameters affecting the CO₂ capture efficiency. It is shown that air velocity and temperature are the primary factors influencing CO₂ uptake. Based on performance data obtained under subfreezing temperatures, a technoeconomic analysis is conducted to evaluate the feasibility of using aqueous solvents in seasonal cold regions. In conclusion, the findings of this study provide valuable insights into siting considerations for deploying solvent-based DAC, thereby contributing to the advancement of sustainable carbon removal solutions.

Air–liquid contactor↗
DOE OSTI2022

Graphene nanoplatelets reinforced cement as a solution to leaky wellbores reinforcing weak points in hydrated Portland cement with graphene nanoparticles improves mechanical and chemical durability of wellbore cements

In this work, to improve the performance of wellbore cement under subsurface environments, graphene nanoplatelets (GNPs) were added in various percentages to Class-H cement slurry. Microstructural characterization of cement slurries cured at 90 °C and 95% RH indicates that GNP modifies the microstructure of hydrated cement by reinforcing pore spaces. As a result, the mechanical properties, such as Young's modulus and axial peak stress, obtained from high temperature triaxial compression tests are significantly improved based on different percentages of added GNPs. Furthermore, the hydrated 1 × 2 inch GNP- Portland cement cores tested under simulated deep wellbore conditions of high-temperature and high pressure (HTHP), appear to have a ductile-like behavior, when compared to a typical brittle nature of Portland cement pastes. From our observations and published data on graphene resistance to fracture, GNP addition to cement is evidently enhancing the flexibility of cement. These improvements would reduce the risks associated with wellbore cement deterioration and a consequent leakage in fossil fuel production, geothermal energy production, CO 2 storage as well as long-term sealing materials in plugging and abandonment of all wellbores at the end of their service life.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The MOOSE fluid properties module

The Fluid Properties module within the Multiphysics Object-Oriented Simulation Environment (MOOSE) is used to compute fluid properties for numerous applications, ranging from nuclear reactor thermal hydraulics to geothermal energy. Those applications drove the development of the module to enable numerous different fluid equations of states, property lookups with primitive and conserved flow variable to cater to pressure and density-driven solvers, and an object-oriented design facilitating expansion and maintenance. Each fluid property is implemented in its own class but inherits capabilities such as automatic differentiation, automated out-of-bounds handling or variable conversion capabilities. Here, this paper presents the module, its design, its user and developer interface, its content in terms of fluids and properties, and several of its applications showing its major role in the MOOSE simulation ecosystem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Numerical modeling and parametric study of a dual purpose underground thermal battery

We report integrating thermal energy storage with building energy systems can enable flexible building electric demands at buildings to help mitigate the mismatch between electricity supply and demand. A novel building heating and cooling system that integrates a dual-source heat pump with hybrid thermal storage named dual-purpose underground thermal battery (DPUTB) has been developed for reshaping building electric demands. The proposed DPUTB integrated geothermal heat pump system is an original innovation that enables Grid-interactive Energy Efficient buildings. This paper focuses on the study of the novel DPUTB. The DPUTB works as both a thermal storage tank (an inner tank) and a ground heat exchanger (an outer tank separated from the inner tank by the insulation material). High fidelity and computationally effective models are needed to predict the performance of the novel DPUTB. This study has developed a simplified dynamic model for the DPUTB according to heat transfer and energy conservation principles and validated it by using experimental data obtained from testing a small-scale DPUTB apparatus. A parametric study was conducted to identify a design that can achieve the target thermal storage performance of load shift and energy efficiency. The parametric study results show that the inner tank shell thermal conductivity and the phase change material's melting point are the two most influencing factors on the performance of the DPUTB. One single full-size DPUTB with the identified design could provide 1-ton cooling (3.51 kW) with the supply temperature lower than 11 °C for 4 hours in summer after being fully charged in 8 hours. The inner tank filled with phase change material is for cooling thermal storage as a latent tank in the design. However, its capacity can be as high as 60 MJ as a sensible water tank for heating storage in winter. In the future study, the DPUTB model will be incorporated into the dual-source heat pump system for evaluating the overall system performance of demand side management in the long term.

25 ENERGY STORAGE↗
DOE OSTI2021

Coupling dynamic in situ X-ray micro-imaging and indentation: A novel approach to evaluate micromechanics applied to oil shale

Oil and gas shales are a system where understanding the mechanical properties at the microscale is of paramount importance, e.g. to better understand the behavior of proppant-shale contacts and their role in the evolution of propped fractures in unconventional reservoirs. Here we show for the first time an experiment coupling indentation testing with in situ X-ray imaging in a Green River shale sample. A full compliance curve has been measured with the sample in water, allowing to visualize the indentation of the sample in function of axial load, in a purpose-built system for combined in situ indentation and X-ray imaging. A series of 3D datasets were used for a digital volume correlation study to obtain local strain fields. This analysis has been complemented with the analysis of cracks. Finally, geomechanical modeling has been carried out to replicate and generalize the observed behavior in the shale. This study validated this experimental approach, providing a breakthrough in understanding micro-mechanics in shales, and demonstrates how this class of experiments can be important for studies involving the prediction of the evolution of propped fractures in shale reservoirs, with possible applications in a much larger number of application fields (geothermal, materials science, etc.)

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗
DOE OSTI2023

Mechanisms of dissolution from gibbsite step edges elucidated by ab initio molecular dynamics with enhanced sampling

Surface reactions underpin our collective understanding of mineral dissolution. This impacts a range of predictive geochemical models (for weathering, the fate and transport of metals) as well as industrial utilization of minerals. Yet atomistic details are rarely known due to the complex mineral/fluid interfacial environment. There is significant need to understand the mechanistic details of dissolution reactions and how they depend on surface morphology and solution conditions. This work utilizes surface pit models in conjunction with changes to solution composition that mimic pH to explore surface detachment of aluminate from gibbsite, which is a primary source of Al in soils and within the industrial processing of aluminum. Ab initio molecular dynamics simulations with enhanced sampling has been used to explore the detailed process of the detachment, the results of which indicate two potential pathways that are differentiated based upon the extent of water hydration. Here, the heights of the energy barriers depend upon the local morphology which influence the number of bridges (quasi-)simultaneously broken (1 or 2) or the Al-O coordination of the neighboring aluminum atoms (5 or 6) at the armchair edge. pH effects are significant, with a nearly 50% reduction in barrier height under alkaline conditions that are relevant to geothermal fluids and Al extraction from minerals.

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

Evaluation of accuracy and convergence of numerical coupling approaches for poroelasticity benchmark problems

Accurate modeling of subsurface flow and transport processes is vital as the prevalence of subsurface activities such as carbon sequestration, geothermal recovery, and nuclear waste disposal increases. Computational modeling of these problems leverages poroelasticity theory, which describes coupled fluid flow and mechanical deformation. Although fully coupled monolithic schemes are accurate for coupled problems, they can demand significant computational resources for large problems. In this work, a fixed stress scheme is implemented into the Sandia Sierra Multiphysics toolkit. Two implementation methods, along with the fully coupled method, are verified with one-dimensional (1D) Terzaghi, 2D Mandel, and 3D Cryer sphere benchmark problems. The impact of a range of material parameters and convergence tolerances on numerical accuracy and efficiency was evaluated. Overall the fixed stress schemes achieved acceptable numerical accuracy and efficiency compared to the fully coupled scheme. However, the accuracy of the fixed stress scheme tends to decrease with low permeable cases, requiring the finer tolerance to achieve a desired numerical accuracy. For the fully coupled scheme, high numerical accuracy was observed in most of cases except a low permeability case where an order of magnitude finer tolerance was required for accurate results. Finally, a two-layer Terzaghi problem and an injection–production well system were used to demonstrate the applicability of findings from the benchmark problems for more realistic conditions over a range of permeability. Simulation results suggest that the fixed stress scheme provides accurate solutions for all cases considered with the proper adjustment of the tolerance. Furthermore, this work clearly demonstrates the robustness of the fixed stress scheme for coupled poroelastic problems, while a cautious selection of numerical tolerance may be required under certain conditions with low permeable materials.

58 GEOSCIENCES↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2023

Prediction and demonstration of periodic tensile cracking in rate-dependent porous cement

Periodic dilation bands are shown to occur when cylinder-shaped saturated cement specimens are subjected to rapid decompression, with the spacing between the bands scaling as a power law of the decompression time. This behavior is predicted by theory, which is specified to this scenario based on prior work that was more general and applied to compaction. Specifically, the formation of the bands is shown to coincide with periodicity that naturally arises in the effective stress when the material follows a Terzaghi-type consolidation law, although in reverse for dilation and with material deformation described by a rate-dependent viscoplastic law. When strain rate is proportional to effective stress to a power that is greater than one, periodic regions of tensile effective stress arise through the poromechanical fluid-solid coupling. The theoretically-predicted exponent of the resulting power-law relationship between dilation band spacing and unloading time successfully brackets the experimental results, for which it is observed that the exponent is slightly stronger than a square-root relationship (0.28 to 0.67) at testing temperatures of 20 °C and 90 °C. Further, these results demonstrate that rapid-depressurization leads to periodic fracturing that could, on the one hand, be detrimental to the isolation provided by cement used to seal wellbores in the petroleum industry. On the other hand, the dilation bands could also be favorable to production if generated in low-permeability reservoir rocks such as shales that are targeted for petroleum production or granites that are targeted for geothermal energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2023

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗
DOE OSTI2022

Large-scale paleo water-table rise in a deep desert aquifer recorded by dissolved noble gases

Illuminating past hydrological and climatological conditions in arid regions may provide insights into future trends in groundwater availability. The main objective of this study is to explore paleorecharge processes in the deep regional Nubian Sandstone Aquifer (NSA), which stretches below the hyperarid deserts of the Sinai Peninsula (Egypt) and the Negev (Israel), using dissolved noble gas data. Extremely high amounts of dissolved excess air (Ne concentrations up to 4 × solubility equilibrium) were observed in the ancient ( 81 Kr-depleted) groundwater that was recharged during past pluvial epochs in the basin. Here, the observed unique excess air signal and the clear spatiotemporal structure of the noble gas compositions in the aquifer are hypothesized to reflect two major characteristics of the groundwater system: (a) large-scale, long-term rises in the water table that facilitate the entrapment and dissolution of air bubbles under increased hydrostatic pressure, and (b) the geological settings in the southern recharge area in Sinai, including an almost horizontal position of the layers and intercalation of low-permeability rock formations within the aquifer, which inefficiently reject air bubbles from groundwater, thereby allowing for substantial entrapment of bubbles that dissolve as the water table rises. Enhanced dissolution of entrapped air was also observed in other paleo-groundwater-containing regional sandstone aquifers across North Africa, which, together with the recent findings from the NSA, suggest that largescale water table fluctuations have likely extensively occurred during past pluvial periods over these (and possibly other) arid regions. A preliminary assessment of noble gas recharge temperatures (NGTs) indicates an apparent decoupling of surface and water table temperatures in the case of deep aquifers with a thick (hundreds of m) unsaturated zone. This observed decoupling calls for a re-evaluation of previously obtained NGT records and a need for future work to consider the modification of NGTs due to soil air fractionation and geothermal heating of the underlying deep unsaturated zones.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hydraulic fracture mapping using wavelet-based fusion of wave transmission and emission measurements

Fracture characterization is essential to hydrocarbon and geothermal exploration and production as well as carbon geo-sequestration. Active and passive measurements have been extensively applied to map fractures across various scales. Active and passive methods have inherent advantages and limitations that complement each other. Here, in this study, to improve the mapping of embedded fractures and the surrounding geomechanically altered regions, we integrate active shear-wave transmission measurement with passive acoustic-emission measurement collected during a lab-scale hydraulic-fracturing experiment. The proposed data fusion leverages the advantages of the two modalities of measurements, while minimizing their individual limitations. Polarized shear-wave transmission (active) measurement was collected before and after the hydraulic fracturing, and the acoustic emission (passive) measurement was collected during hydraulic fracturing. Two sets of two-dimensional maps of fracture and fracture-induced damage in axial, median and frontal planes were obtained by separately processing active and passive measurements. These 2D maps are then fed to wavelet-based image fusion technique to integrate the two sources of information for the reliable mapping/imaging of the embedded fractures and the surrounding geomechanically altered regions.

58 GEOSCIENCES↗
DOE OSTI2021

Recent seismicity at Ceboruco Volcano (Mexico)

Magma and related hydrothermal fluid movement, and their interaction with solid rock, in active volcanic regions, generate a wide variety of seismic waves whose characterization can mitigate the risk of a potential eruption. Located in the western region of the Trans-Mexican Volcanic Belt, Ceboruco Volcano, whose last eruptive period was 1870–1875, is considered to be one of the most hazardous volcanoes in Mexico. We have conducted a detailed study of the seismicity in the surroundings of Ceboruco's volcanic edifice to assess the current state of this volcano. A dense temporary seismic network with 25 seismic stations in an area of 16 km × 16 km was deployed between November 2016 and July 2017, as part of the P-24 project of the CeMIEGeo consortium; this effort has allowed the detection of 81 earthquakes concentrated beneath the crater with depths between 4 and 8 km. In this study, we observe that the recorded seismicity occurs in swarms, and we specifically identify four sequences that we characterize in detail via the first focal mechanisms available for this volcano. Our results suggest a change in the local seismicity distribution compared to earlier observations, which reported seismic activity near the volcano edifice associated with fluid migration along zones of weakness related to the extensional stresses of the Tepic-Zacoalco rift. The changes in seismic patterns and obtained focal mechanisms are consistent with observed fluid effects at many geothermal sites worldwide, but also could suggest resumption of activity at this currently dormant volcano.

58 GEOSCIENCES↗
DOE OSTI2025

Ecology of methyl-coenzyme M reductase encoding Thermoproteota

The recent demonstration that members of at least three classes of archaea affiliated with the Thermoproteota superphylum are involved in the production of the climate-active gas methane has sparked discussions about how well we understand the diversity of methanogens. Here, in this work, we show that members of all three of these lineages, as well as several other, yet uncultured and physiologically uncharacterized groups within the Thermoproteota that encode the key enzyme of anaerobic methane cycling, methyl-coenzyme M reductase (MCR), are widely distributed in anoxic ecosystems. We postulate that the taxonomic, metabolic, and ecological diversity of methanogenic and MCR-encoding Thermoproteota are poorly understood, and that the contribution of methylotrophic and thermoproteotal methanogenesis to methane production is largely unknown. We hypothesize that thermoproteotal methanogens could contribute, potentially substantially, to methane emissions in many anoxic environments that harbor methylated precursors, including wetlands, sediments, peat, rice paddies, wastewater sludge, and geothermal systems. We highlight the necessity to experimentally test the (eco)physiology of these widely distributed archaea using both culture-dependent ( in vitro ) and culture-independent ( in situ ) approaches to assess their potential contribution to methane emissions. Last, we stress the importance of remaining agnostic about the physiology of MCR-encoding Thermoproteota in the absence of experimental data because most of these archaea also carry the genetic potential to grow non-methanogenically.

MCR↗