Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Data Science in Science”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,642 records · Page 80Linked to original sources
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗
USGS1986

Regionalization of low-flow characteristics of Tennessee streams

Procedures for estimating 3-day 2-year, 3-day 10-year, 3-day 20-year, and 7-day 10-year low flows at ungaged stream sites in Tennessee are based on surface geology and drainage area size. One set of equations applies to west Tennessee streams, and another set applies to central and east Tennessee streams. The equations do not apply to streams where flow is significantly altered by activities of man. Standard errors of estimate of equations for west Tennessee are 24 to 32% and for central and east Tennessee 31 to 35%. Streamflow recession indexes, in days/log cycle, are used to account for effects of geology of the drainage basin on low flow of streams. The indexes in Tennessee range from 32 days/log cycle for clay and shale to 350 days/log cycle for gravel and sand, indicating different aquifer characteristics of the geologic units that sustain streamflows during periods of no surface runoff. Streamflow recession rate depends primarily on transmissivity and storage characteristics of the aquifers, and the average distance from stream channels to basin divides. Geology and drainage basin size are the most significant variables affecting low flow in Tennessee streams according to regression analyses. (Author 's abstract)

Tennessee↗
USGS1998

Methods and guidelines for effective model calibration; with application to UCODE, a computer code for universal inverse modeling, and MODFLOWP, a computer code for inverse modeling with MODFLOW

This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.

Water-Resources Investigations Report↗
USGS2026

Overstorey mortality promotes juvenile piñon pine growth during favourable weather at cooler, wetter sites

Hotter droughts have resulted in widespread tree die-off events globally, frequently leading to regeneration failure. Dry forest recovery often depends on the growth and survival of extant juvenile trees. However, it is unclear how microenvironmental changes following overstorey tree die-off affect juvenile trees, particularly in dryland systems where tree recruitment is typically limited by water availability and heat stress. We simulated an overstorey tree die-off event by girdling trees in piñon-juniper woodlands across the south-western United States. We sampled juvenile piñon pine growth from live and dead overstorey treatments across six study sites spanning a regional latitudinal gradient and local elevational gradients. We examined how juvenile branch and needle growth differed between live and dead overstorey treatments, and whether responses varied with weather conditions and juvenile tree size following overstorey mortality. We found greater juvenile branch and needle growth under dead compared with live overstorey trees for 2 years following overstorey mortality at mid- and high-elevation sites which are typically cooler and wetter than the other sites. These observed growth releases were contingent on favourable post-mortality weather conditions. Higher growth under dead overstorey occurred at sites experiencing near-average climatic water deficits compared with sites experiencing above-average climatic water deficits. Growth at all sites increased from the first to second year after overstorey mortality. Across sites, growth was unrelated to juvenile tree size. Synthesis . Our results underscore differentiation in juvenile responses to overstorey tree die-off driven by local site conditions and weather across the range of Pinus edulis . Overstorey mortality resulted in consistently higher juvenile growth only at climatically favourable sites and during favourable weather, while unmeasured microsite differences could help account for variation observed at the hottest and driest site. Results from less climatically favourable sites suggest that overstorey trees neither directly limit nor facilitate juvenile growth, though further study over longer timeframes is needed to resolve the pace and magnitude of potential recovery or decline. Overall, juvenile vigour may be promoted following overstorey mortality only in a narrow spatial (site) and temporal (weather) environmental context, suggesting additional vulnerabilities for piñon populations under more arid conditions.

Arizona, Colorado↗
USGS2023

Stratigraphic notes—Volume 1, 2022

This is the first volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 . This first volume ("Stratigraphic notes—Volume 1, 2022") includes articles that provide guidance for those who wish to submit papers to “Stratigraphic Notes,” as well as information on how to make your manuscripts compliant for geologic names reviews and how to organize your paper’s content to facilitate inclusion of new or revised names in Geolex. This volume also includes some specific guidance on conducting geologic names reviews of geologic and hydrogeologic reports.

Professional Paper↗
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗
USGS2017

Field-trip guide to Mount Hood, Oregon, highlighting eruptive history and hazards

This guidebook describes stops of interest for a geological field trip around Mount Hood volcano. It was developed for the 2017 International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Scientific Assembly in Portland, Oregon. The intent of this guidebook and accompanying contributions is to provide an overview of Mount Hood, including its chief geologic processes, magmatic system, eruptive history, local tectonics, and hazards, by visiting a variety of readily accessible localities. We also describe coeval, largely monogenetic, volcanoes in the region. Accompanying the field-trip guidebook are separately authored contributions that discuss in detail the Mount Hood magmatic system and its products and behavior (Kent and Koleszar, this volume); Mount Hood earthquakes and their relation to regional tectonics and the volcanic system (Thelen and Moran, this volume); and young surface faults cutting the broader Mount Hood area whose extent has come to light after acquisition of regional light detection and ranging coverage (Madin and others, this volume). The trip makes an approximately 175-mile (280-kilometer) clockwise loop around Mount Hood, starting and ending in Portland. The route heads east on Interstate 84 through the Columbia River Gorge National Scenic Area. The guidebook points out only a few conspicuous features of note in the gorge, but many other guides to the gorge are available. The route continues south on the Mount Hood National Scenic Byway on Oregon Route 35 following Hood River, and returns to Portland on U.S. Highway 26 following Sandy River. The route traverses rocks as old as the early Miocene Eagle Creek Formation and overlying Columbia River Basalt Group of middle Miocene age, but chiefly lava flows and clastic products of arc volcanism of late Miocene to Holocene age.

Oregon↗
USGS2025

Oligocene–Miocene development and evolution of the south Dome Rock Mountains basin, lower Colorado River corridor, Arizona, USA

Sedimentary basins in the Colorado River extensional corridor record large-magnitude Basin and Range extension and younger dextral shear deformation in the evolving Pacific−North America plate boundary. The south Dome Rock Mountains basin is located in west-central Arizona (USA), where the history of basin evolution, style of deformation, and timing of the transition between extension and dextral shear are not well constrained. We integrate new geologic mapping of the south Dome Rock Mountains basin with zircon U-Pb geochronology of six marker beds to characterize the timing of basin evolution and the slip history of the south Dome Rock Mountains normal fault. Structural analysis defines structures consistent with regional extension and younger dextral shear deformation. We use sedimentological and petrological analysis to interpret the depositional environments of three basin sequences. A lower basin sequence consists of fluvial strata deposited in an internally drained intermontane basin from ca. 35 Ma to 24.4 Ma. A sequence of volcanic rocks was emplaced in the basin between 24.4 Ma and 23.3 Ma, before the onset of local extensional faulting. An upper basin sequence of coarse conglomerate and sedimentary breccia was deposited on the hanging wall of a half graben structure after 23.3 Ma and through ca. 12.8 Ma, synchronous with 3−7 km of dip slip on the west-dipping, listric south Dome Rock Mountains fault. The basin was subsequently deformed by dextral shear after 12.8 Ma and before 4.8 Ma. These results document how distributed deformation related to the evolving Pacific−North America plate boundary occurred >100 km from the primary plate boundary.

Arizona, California↗
USGS2025

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone↗
USGS2025

Microbial ecology of permafrost soils: Populations, processes, and perspectives

Permafrost microbial research has flourished in the past decades, due in part to improvements in sampling and molecular techniques, but also the increased focus on the permafrost greenhouse gas feedback to climate change and other ecological processes in high latitude and alpine permafrost soils. Permafrost microorganisms are adapted to these extreme environments and remain active at low temperatures and when resources are limited. They are also an important component of global elemental cycles as they regulate organic matter turnover and greenhouse gas production, particularly as permafrost thaws. Here we review the permafrost microbiology literature coupled with an exploration of its historical aspects, with a particular focus on a new understanding advanced by molecular biology techniques. We further identify knowledge gaps and ways forward to improve our understanding of microbial contributions to ecosystem biogeochemistry of permafrost-affected systems.

Permafrost and Periglacial Processes↗
USGS2010

Divisions of geologic time—Major chronostratigraphic and geochronologic units

Effective communication in the geosciences requires consistent uses of stratigraphic nomenclature, especially divisions of geologic time. A geologic time scale is composed of standard stratigraphic divisions based on rock sequences and is calibrated in years. Over the years, the development of new dating methods and the refinement of previous methods have stimulated revisions to geologic time scales. Advances in stratigraphy and geochronology require that any time scale be periodically updated. Therefore, Divisions of Geologic Time, which shows the major chronostratigraphic (position) and geochronologic (time) units, is intended to be a dynamic resource that will be modified to include accepted changes of unit names and boundary age estimates. This fact sheet is a modification of USGS Fact Sheet 2007-3015 by the U.S. Geological Survey Geologic Names Committee.

Fact Sheet↗
USGS2024

Seismicity and anisotropic imaging reveal an active detachment beneath the northern Alaska Range foothills

North of the Denali Fault, the collision between the Yakutat block with North America is accommodated by a fold-thrust belt giving rise to the northern Alaska Range foothills. At the western end, the Kantishna Hills anticline hosts prominent microseismicity and surface deformation, interpreted as active folding of the Kantishna Hills anticline above a midcrustal detachment. We test for this detachment by using anisotropy-aware receiver functions to image fabric contrasts within the crust in context with seismicity. Seismic stations near the crest of the Kantishna Hills anticline and near its southern flank show a single strong contrast in dipping fabric at depths of 12–13 km near microseismicity clustering depths, consistent with a detachment plane beneath the fold. A minimum b -value at 10–13 km depth is consistent with seismicity on the detachment, compatible with the imaged anisotropic contrast, while off-fault seismicity is shallower and deeper with smaller magnitudes. South-dipping imbricate thrusts in schist characterize the northern Alaska Range foothills structure. This supports our interpretation of the observed anisotropy as reflecting SSW–SSE-dipping foliation above a detachment at ∼10–13 km depth that exploits existing crustal weaknesses along subtle fabric contrasts observed in the seismically quiescent region north of the actively deforming belt.

Alaska↗
USGS2025

Road salt collection and redistribution at an urban rain garden on sandy soil, Gary, Indiana

Rain gardens installed as green infrastructure to divert storm runoff from entering combined sewers also collect dissolved constituents and particulates. An urban rain garden in northwestern Indiana, USA, was continuously monitored from November 2019 to May 2021 to evaluate the fate of dissolved constituents entering the rain garden in runoff. Physical and chemical properties of soils in the rain garden were also monitored, along with underlying groundwater. Linear regression models relating specific conductance to chloride concentration indicated that the 0.0371-ha (3998 square feet) rain garden collected approximately 1490 kg (3285 pounds) of road salt from the surrounding 0.2228 ha (24,500 square feet) of impervious surfaces. Soils and groundwater were seasonally affected by road salt application but carryover from year to year was not indicated. Rain garden soil permeability (5.20 × 10 −5 to 9.72 × 10 −5 m/s) remained unchanged during the study period and soil organic carbon generally increased under native vegetation. The results suggest that a rain garden built on sandy soil can divert substantial quantities of runoff and dissolved constituents from combined sewers; however, chloride is transported to sub-infrastructure groundwater that eventually discharges to adjacent waterways with concentrations lower than those observed in runoff.

Indiana↗
USGS2024

Developing, testing, and communicating earthquake forecasts: Current practices and future directions

While deterministically predicting the time and location of earthquakes remains impossible, earthquake forecasting models can provide estimates of the probabilities of earthquakes occurring within some region over time. To enable informed decision-making of civil protection, governmental agencies, or the public, Operational Earthquake Forecasting (OEF) systems aim to provide authoritative earthquake forecasts based on current earthquake activity in near-real time. Establishing OEF systems involves several nontrivial choices. This review captures the current state of OEF worldwide and analyzes expert recommendations on the development, testing, and communication of earthquake forecasts. An introductory summary of OEF-related research is followed by a description of OEF systems in Italy, New Zealand, and the United States. Combined, these two parts provide an informative and transparent snapshot of today's OEF landscape. In Section 4, we analyze the results of an expert elicitation that was conducted to seek guidance for the establishment of OEF systems. The elicitation identifies consensus and dissent on OEF issues among a non-representative group of 20 international earthquake forecasting experts. While the experts agree that communication products should be developed in collaboration with the forecast user groups, they disagree on whether forecasting models and testing methods should be user-dependent. No recommendations of strict model requirements could be elicited, but benchmark comparisons, prospective testing, reproducibility, and transparency are encouraged. Section 5 gives an outlook on the future of OEF. Besides covering recent research on earthquake forecasting model development and testing, upcoming OEF initiatives are described in the context of the expert elicitation findings.

Reviews of Geophysics↗
USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗
USGS2024

Cross section N–N' through the Valley and Ridge province of the southern Appalachian basin, from Greene County, west-central Alabama, to Bibb County, central Alabama

Introduction Geologic cross section N–N′ is the sixth in a series of geologic cross sections constructed by the U.S. Geological Survey to document and improve understanding of the geologic framework and petroleum systems of the Appalachian basin. Cross section N–N′ provides a regional view of the structural and stratigraphic framework of the Appalachian basin in the Valley and Ridge province in western and central Alabama; it spans approximately 69 miles (mi) (111 kilometers [km]). This geologic cross section is a companion to geologic cross sections E–E′ , D–D′ , C–C′ , I–I′ , and A–A′ that are located approximately 350 to 550 mi (563 to 885 km) to the northeast. Cross section N–N' complements earlier geologic cross sections through the Alabama part of the Appalachian basin. Although some of the other cross sections show more structural and stratigraphic detail, they are of more limited extent geographically and stratigraphically. Cross section N–N′ contains information that is useful for evaluating energy resources in the Appalachian basin. Although the Appalachian basin petroleum systems are not shown on the cross section, many of their key elements (such as source rocks, reservoir rocks, seals, and traps) can be inferred from lithologic units, unconformities, and geologic structures shown on the cross section. Other aspects of petroleum systems (such as the timing of petroleum generation and petroleum migration pathways) may be evaluated by burial history, thermal history, and fluid flow models based on what is shown on the cross section. In addition, cross section N–N′ may be used as a reconnaissance tool to identify plausible geologic structures and strata for the subsurface storage of liquid waste or for the sequestration of carbon dioxide.

Alabama↗
USGS2024

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon↗
USGS2026

Deep subsurface organic-rich shale supports abundant, diverse, and novel fungi

As Earth’s principal reservoir of organic carbon and microbial biomass, the deep subsurface hosts microorganisms capable of mobilizing this once-sequestered carbon. Contrary to standard assumptions of eukaryotic scarcity, this study documents abundant fungal communities, ranging from 4.2 × 10 3 to 6.8 × 10 3 fungal cells mL −1 , across a methane-producing organic-rich shale 247–556 meters below the surface. Although fungal:bacterial cell ratios ranged from 1:7028 to 1:713, application of biomass conversion factors developed for oceanic systems yielded a median fungal:bacterial biomass ratio of 1:4.7. 16S rRNA gene amplicons revealed bacterial and archaeal communities mirroring those found in well-characterized extremophilic, carbon-degrading environments, while sequencing of 18S rRNA gene and ITS rRNA spacer amplicons collectively identified a eukaryotic hotspot with 689 fungal OTUs across six phyla. The dominant fungal classes, Agaricomycetes and Dothideomycetes, are well-established degraders of recalcitrant carbon compounds at the surface, suggesting they may similarly contribute to organic matter degradation and ecosystem maintenance in the subsurface. Cultivation and isolation efforts yielded 205 fungal strains, including 13 candidate novel taxa, underscoring the deep subsurface as an underexplored eukaryotic habitat. Stable carbon isotopes indicate methane is predominantly generated via microbial conversion of the fossil carbon, while water isotopes suggest in situ geochemical conditions have been relatively stable since the Late Pleistocene, with subglacial recharge as a plausible mechanism for microbial introduction. Collectively, these findings suggest that fungi are underrecognized contributors to organic matter transformation and functional diversity in the deep biosphere, revealing a critical gap in our understanding of deep subsurface ecosystem processes.

Indiana, Michigan, Ohio↗