Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗
DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗
DOE OSTI2024

Electrical Resistivity Changes During Heating Experiments Unravel Heterogeneous Thermal‐Hydrological‐Mechanical Processes in Salt Formations

Abstract Rock salt is considered a suitable medium for the permanent disposal of heat‐generating radioactive waste due to its isolation properties. However, excavation damage and heating induce complex and heterogeneous thermal‐hydrological‐mechanical (THM) processes across different zones. Quantifying this heterogeneity is crucial for accurate long‐term performance assessment models, but traditional methods lack the necessary resolution. This study employs 4D electrical resistivity tomography (ERT) monitoring during controlled heating experiments in a salt formation to unravel the spatiotemporal dynamics of THM processes. Advanced time‐lapse inversion and clustering analysis quantify subsurface properties and map the heterogeneity of THM dynamics. The ERT results can estimate subsurface properties and delineate the damaged and intact zones, enabling appropriate parameterization and representation of processes for long‐term modeling. This approach may be used in further improving the predictive models and ensuring the safe long‐term disposal of radioactive waste in rock salt.

58 GEOSCIENCES↗
DOE OSTI2024

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗
DOE OSTI2024

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗
DOE OSTI2025

Surface Energy Balance Responses to Radiative Forcing in the Central Arctic From MOSAiC and Models

The Arctic surface energy budget (SEB) couples the atmosphere with the sea ice, making it useful for both studying surface processes as well as evaluating models. Improved understanding of atmosphere-ice interactions is required to improve models, requiring year-round observations to address seasonally dependent biases. This work uses novel observations from the MOSAiC expedition to quantify the responses of surface fluxes to radiative forcing over sea ice throughout a complete annual cycle. We identify two primary regimes of flux response: an ice growth regime in winter and an ice melt regime in summer. In the growth regime, changes in radiative forcing impact upwelling longwave, sensible heat, and subsurface heat fluxes, whereas in the melt regime changes in radiative forcing primarily alter the amount of melt and subsurface transmission because the surface temperature is fixed. These observed responses of surface fluxes to radiative forcing are used to evaluate seven weather forecast models during the ice growth regime. In most models, the responses of surface fluxes to radiative forcing do not match observations. Many models also have biased downwelling longwave. One model (the Coupled Arctic Forecast System; CAFS) adequately captures both the mean radiative forcing and the flux responses in winter. CAFS is further evaluated against observations spanning the full MOSAiC year, demonstrating sufficient agreement to provide a more generalized understanding of these SEB process relationships across the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Stormwater Storage and Retention Within an Urban Prairie Wetland Complex

Climate change is expected to increase the frequency and severity of flooding in the Great Lakes region. In many cities, flood-control infrastructure is insufficient to protect against future climate conditions. Consequently, there is increasing focus on stormwater storage provided by urban greenspace, such as wetlands and prairies, but the ecohydrological behavior of these ecosystems is not well understood when they are embedded within cities. To improve understanding of hydrological connectivity between urban areas and natural greenspaces, we deployed a sensor network in Gensburg Markham Prairie (GMP), a large intact prairie-wetland complex in south suburban Chicago. We used the resulting high-frequency time-series data to assess surface-subsurface hydrologic dynamics between upland and low-lying wetland areas, interactions between the prairie and surrounding environment, and stormwater storage provided by the prairie. Rapid infiltration within the prairie during and after storm events provides subsurface flow that stores considerable water, flattens storm hydrographs, and increases the wetland hydroperiod. Much of the stormwater input to GMP derives from the surrounding cityscape. Consequently, storage within the prairie-wetland system reduces and slows stormwater discharge to downstream urban communities. For a typical 5-year 24-hr storm with 10.9 cm of rain, GMP stores 77,100 m 3 , 64% greater than the estimated direct rainfall volume onto the prairie, yielding 30,000 m 3 of offsite stormwater storage. This improved understanding of ecohydrological dynamics in urban prairies and wetlands informs the design and implementation of green infrastructure to meet growing needs for stormwater management.

Rivera, Vivien Anne [Northwestern University, Evan↗
DOE OSTI2026

Anisotropic Melt Inclusions as a Confounding Signal for Ice‐Penetrating Radar Observations

Ice-penetrating radar is a powerful geophysical tool for understanding the subsurfaces of Earth, Mars, and icy moons. Radar reflectivity, attenuation, and birefringence are used to infer subsurface hydrology, englacial temperature, water content, and crystal orientation fabric. However, conventional radar sounding analyses either ignore melt or use classical mixing models which assume spherical melt inclusions, obscuring anisotropic contributions of melt. Here, we use geometric mixing models to calculate the reflectivity, attenuation, and birefringence of temperate ice containing anisotropic melt. We find that anisotropic melt can introduce significant deviations in radar measurements. For instance, melt anisotropy may impact reflectivity-based estimates of water content by up to 30% volume fraction, while attenuation-based estimates may vary by up to 43%. Critically, a melt volume fraction of just $f$ ~ 10 -5 can reproduce birefringence signals previously attributed to ice fabric. We conclude that melt anisotropy may significantly impact analyses of englacial hydrology, subglacial hydrology, and ice fabric.

58 GEOSCIENCES↗
DOE OSTI2026

Constraining the Hydration of Clay Minerals and Abundances of Amorphous Phases in Gale Crater, Mars

Both water and organic matter are required for the development and persistence of life. Phyllosilicates (clay minerals) have high surface areas that easily sorb water and organic matter. The Curiosity rover has investigated several hundred meters of stratigraphy in Gale crater, including where clays were detected from orbit. Previous results have suggested that subsurface hydration is greatest in units with the most abundant clays, suggesting that these minerals may be hydrated. Organics have also been found throughout Gale crater. Smectites are the most common and abundant phyllosilicates in Gale crater samples and can expand and sorb water and organics in interlayer sites. The most common organic sorption processes on Earth typically involve water or hydroxyl, so hydrated phyllosilicates are good candidates for organic preservation. Using newly derived subsurface hydration results with previously published mineralogy and geochemistry, we derived modeled constraints on the abundances of hydrated amorphous phases, “excess” water, and “excess” cations. These “excess” phases are not accounted for by published crystalline phase abundances or by amorphous phases constrained here. We found correlations between smectites and both “excess” water and “excess” cation abundances, indicating that smectites in Gale crater are hydrated and that cation bridging could be a mechanism for sorption of organics. Our results also show the persistence of amorphous sulfates, opal-A, and volcanic or impact glass, which indicate low water-rock interactions. Increased abundances of sulfates and glass in stratigraphically higher samples may indicate lower water availability and environmental aridification during the time these units were being deposited.

58 GEOSCIENCES↗
DOE OSTI2025

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗
DOE OSTI2026

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mechanochemically responsive polymer enables shockwave visualization

Abstract Understanding the physical and chemical response of materials to impulsive deformation is crucial for applications ranging from soft robotic locomotion to space exploration to seismology. However, investigating material properties at extreme strain rates remains challenging due to temporal and spatial resolution limitations. Combining high-strain-rate testing with mechanochemistry encodes the molecular-level deformation within the material itself, thus enabling the direct quantification of the material response. Here, we demonstrate a mechanophore-functionalized block copolymer that self-reports energy dissipation mechanisms, such as bond rupture and acoustic wave dissipation, in response to high-strain-rate impacts. A microprojectile accelerated towards the polymer permanently deforms the material at a shallow depth. At intersonic velocities, the polymer reports significant subsurface energy absorption due to shockwave attenuation, a mechanism traditionally considered negligible compared to plasticity and not well explored in polymers. The acoustic wave velocity of the material is directly recovered from the mechanochemically-activated subsurface volume recorded in the material, which is validated by simulations, theory, and acoustic measurements. This integration of mechanochemistry with microballistic testing enables characterization of high-strain-rate mechanical properties and elucidates important insights applicable to nanomaterials, particle-reinforced composites, and biocompatible polymers.

Science & Technology - Other Topics↗
DOE OSTI2025

Complex carbonate phases drive geologic CO2 mineralization

Abstract Geologic carbon sequestration in mafic and ultramafic reservoirs is a scalable strategy for carbon dioxide removal, offering permanent storage via mineralization as stable carbonates. However, there is limited information on the structure and composition of key mineralization endpoints during sequestration. Here, we unravel the atomic structure, composition, and nanoscale morphology of carbonates recovered from a field-scale demonstration of CO 2 mineralization in basalt. Using transmission electron microscopy, we mapped mineralogical variations from the initial to later stages of subsurface carbonate growth and identified a previously unknown cation-ordered ankerite phase that exerts a primary control over carbonation processes. This study has provided a new understanding of subsurface carbonation pathways which will impact the parameterization of predictive geochemical models for future sequestration efforts in basalt formations.

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

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2024

Efficient and generalizable nested Fourier-DeepONet for three-dimensional geological carbon sequestration

Geological carbon sequestration (GCS) involves injecting CO2 into subsurface geological formationsfor permanent storage. Numerical simulations could guide decisions in GCS projects by predictingCO 2 migration pathways and the pressure distribution in storage formation. However, these simula-tions are often computationally expensive due to highly coupled physics and large spatial-temporalsimulation domains. Surrogate modelling with data-driven machine learning has become a promis-ing alternative to accelerate physics-based simulations. Among these, the Fourier neural operator(FNO) has been applied to three-dimensional synthetic subsurface models. Despite its good accuracyin simulating CO 2 plume migration, it requires large computational resources in training and alsolacks generalizability. Here, to further improve performance, we have developed a nested Fourier-DeepONet by combining the expressiveness of the FNO with the modularity of a deep operatornetwork (DeepONet). This new framework is twice as efficient as a nested FNO for training and has atleast 80% lower GPU memory requirement due to its flexibility to treat temporal coordinates sepa-rately. These performance improvements are achieved without compromising prediction accuracy.In addition, the generalization and extrapolation ability of nested Fourier-DeepONet beyond thetraining range has been thoroughly evaluated. Nested Fourier-DeepONet outperformed the nestedFNO for extrapolation in time with more than 50% reduced error. It also exhibited good extrapolationaccuracy beyond the training range in terms of reservoir properties, number of wells, and injectionrate.

Lee, Jonathan E. [Department of Chemical and Envir↗
DOE OSTI2024

Canopy cover and microtopography control precipitation-enhanced thaw of ecosystem-protected permafrost

Abstract Northern high-latitudes are projected to get warmer and wetter, which will affect rates of permafrost thaw and mechanisms by which thaw occurs. To better understand the impact of rain, as well as other factors such as snow depth, canopy cover, and microtopography, we instrumented a degrading permafrost plateau in south-central Alaska with high-resolution soil temperature sensors. The site contains ecosystem-protected permafrost, which persists in unfavorable climates due to favorable ecologic conditions. Our study (2020–2022) captured three of the snowiest years and three of the four wettest years since the site was first studied in 2015. Average thaw rates along an across-site transect increased nine-fold from 6 ± 5 cm yr −1 (2015–2020) to 56 ± 12 cm yr −1 (2020–2022). This thaw was not uniform. Hummock locations, residing on topographic high points with relatively dense canopy, experienced only 8 ± 9 cm yr −1 of thaw, on average. Hollows, topographic low points with low canopy cover, and transition locations, which had canopy cover and elevation between hummocks and hollows, thawed 44 ± 6 cm yr −1 and 39 ± 13 cm yr −1 , respectively. Mechanisms of thaw differed between these locations. Hollows had high warm-season soil moisture, which increased thermal conductivity, and deep cold-season snow coverage, which insulated soil. Transition locations thawed primarily due to thermal energy transported through subsurface taliks during individual rain events. Most increases in depth to permafrost occurred below the ∼45 cm thickness seasonally frozen layer, and therefore, expanded existing site taliks. Results highlight the importance of canopy cover and microtopography in controlling soil thermal inputs, the ability of subsurface runoff from individual rain events to trigger warming and thaw, and the acceleration of thaw caused by consecutive wet and snowy years. As northern high-latitudes become warmer and wetter, and weather events become more extreme, the importance of these controls on soil warming and thaw is likely to increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Methane emission hotspots in a boreal forest-fen mosaic potentially linked to deep taliks

Permafrost thaw is transforming boreal forests into mosaics of wetlands and drier uplands. Topographic controls on hydrological and ecological conditions impact methane (CH 4 ) fluxes, contributing to uncertainty in local and regional CH 4 budgets and underlying drivers. The objective of this study was to explore CH 4 fluxes and their drivers in a transitioning boreal forest-fen ecosystem (Goldstream Valley, Alaska, USA). This landscape is characterized by thawing discontinuous permafrost and heterogeneous mosaics of fens, collapse-scar channels, and small mounds of permafrost soils. From a survey in July 2021, observed chamber CH4 fluxes included fen areas with intermediate to very high emissions (29.8–635.3 mg CH 4 m −2 d −1 ), clustered locations with CH 4 uptake (−2.11 to −0.7 mg CH 4 m −2 d −1 ), and three anomalous emission hotspots (342.4–772.4 mg CH 4 m −2 d −1 ) that were located near samples with lower emissions. Some surface and near-surface variables partially explained the spatial variation in CH 4 flux. Log-transformed CH 4 flux had a positive linear relationship with soil moisture at 20 cm depth ( R 2 = 0.31, p -value < 1e-5) and negative linear relationships with microtopography ( R 2 = 0.13, p -value < 0.006) and slope ( R 2 = 0.28, p -value < 2e-5). Methane emissions generally occurred in flat, wet, graminoid-dominated fens, whereas CH 4 uptake occurred on permafrost mounds dominated by feather mosses and woody vegetation. However, the CH 4 hotspots occurred on drier, slightly sloped locations with low or undetectable near-surface methanogen abundance, suggesting that CH 4 was produced in deeper soils. When the hotspot samples were omitted, log-transformed CH 4 flux had a positive linear relationship with near-surface methanogen abundance ( R 2 = 0.29, p -value = 0.0023), and stronger linear relationships with soil moisture, slope, and soil macronutrient concentrations. Our findings suggest that some CH 4 emission hotspots could arise from CH 4 in deep taliks. The inference that methanogenesis occurs in deep taliks was strengthened by the identification of intrapermafrost taliks across the study area using low-frequency geophysical induction. This study assesses surface spatial heterogeneity in the context of subsurface permafrost conditions and highlights the complexity of CH 4 flux patterns in transitioning forest-wetland ecosystems. To better inform regional CH 4 budgets, further research is needed to understand the spatial distribution of terrestrial CH 4 hotspots and to resolve their surface, near-surface, and subsurface drivers.

boreal↗