USGS Science⌕ Search

SEARCH · USGS Science

Results for “The ISME Journal”

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,663 records · Page 60Linked to original sources
USGS1977

Summary of results from electromagnetic and galvanic soundings on Kilauea Iki lava lake, Hawaii

The purpose of this report is to summarize various electrical sounding studies made on the Kilauea Iki lava lake and to present some of the preliminary interpretations that resulted from the measurements. The following discussion is not intended as a complete interpretation, but is intended to summarize the major points made in an invited oral presentation at the American Geophysical Union fall meeting held in San Francisco, Calif. on December 6, 1976. This paper was part of a group of papers dealing with various geophysical studies recently made on the Kilauea Iki lava lake. A general plan map of Kilauea Iki lava lake is shown in figure 1. Electromagnetic (EM) loop-loop soundings, D.C. resistivity soundings, and VLF apparent resistivity measurements were made on the lave lake surface. All of these surveys were made to determine the geoelectric structure of the lava lake and to determine parameters of the molten basaltic lens that exist in the partially solidified lava lake.

Hawaii↗
USGS2024

Utilization of stochastic ground motion simulations for scenario-based performance assessment of geo-structures

Probabilistic seismic performance assessments of engineered structures can be highly sensitive to the seismic input excitation and its variability. In the present study, the scenario-based performance assessment recommended by Federal Emergency Management Agency (FEMA) P-58 guidelines is adopted to estimate seismic fragility of concrete dams for various seismic hazard scenarios. Due to the scarcity of recorded ground motions and thereby their poor representation of uncertainties, stochastic ground motion simulation methods are utilized to obtain the required input excitations. Moreover, to understand the uncertainty in ground motion simulation models, two broadband stochastic simulation models are used to generate input excitations representing six seismic hazard scenarios defined by earthquake magnitude, source-to-site distance, and soil conditions. Optimal intensity measure parameters for each scenario are identified using a systematic procedure that considers criteria such as efficiency, practicality, proficiency, sufficiency, and hazard compatibility. Fragility curves and surfaces are derived using the cloud analysis technique, taking into account various damage measures and limit state functions. The study finds that the derived fragility curves are particularly sensitive to the selection of earthquake scenarios, the choice of records, and the methods used to calculate fragility curves, with less sensitivity observed to different engineering demand parameters. Given this sensitivity, particularly to ground motion selection, the study highlights the necessity of incorporating both model-to-model variability (epistemic uncertainty) and record-to-record variability (aleatory uncertainty), alongside the established material and modeling uncertainties, in the probabilistic seismic assessment.

Reliability Engineering and System Safety (RESS)↗
USGS2024

Geophysical modeling of a possible blind geothermal system near Battle Mountain, NV

The northeastern portion of the Reese River basin in north-central Nevada is the focus of detailed geophysical and geological studies as part of the INGENIOUS project, which aims to identify new, commercially viable hidden geothermal systems in the Great Basin region of the western U.S. This location, herein referred to as Argenta Rise, occupies a broad (~15km wide) left-step between major range-front fault systems along the northwestern edge of the Shoshone Range and Argenta Rim, with numerous ENE-striking intra-basin faults presumably accommodating sinistral-normal oblique slip across the step-over. Four discrete regions have been identified within the study area that have favorable structural settings for hosting a blind hydrothermal system. However, with no definitive or extensive surface manifestations of an active hydrothermal system (e.g., geysers, steam vents, sinter, etc.), detailed geophysical studies are necessary to resolve subsurface geology and structure, and identify zones of enhanced structural complexity that may promote hydrothermal fluid flow. Hence, we collected high-resolution gravity, MT, and rock property data (density, magnetic susceptibility), and analyzed the recently acquired GeoDAWN aeromagnetic data to characterize potential geothermal resources in this region. Using the new geophysical datasets, we jointly modeled gravity and magnetic data along a series of intersecting 2D profiles that integrated information from recent, local-scale fault mapping. Rock property measurements performed on outcrops and hand samples throughout the study area constrained the models. The MT data were used to construct a 3D resistivity model that highlights the location of inferred alteration and fluids in the subsurface. Combined MT and potential field results reveal which structures may be most important for controlling hydrothermal fluid migration, as well as which geologic units may host hydrothermal fluids. Our gravity derived depth to basement surface coincides well with the base of shallow conductive anomalies, suggesting hydrothermal fluids may be confined to basin fill sediments and volcanics. This work supports our development of 3D geophysical and geologic models that are focused along the western flank of the northern Shoshone Range and aids the process of selecting sites for temperature gradient drilling.

Nevada↗
USGS2025

Estimating the probability of export restrictions to inform mineral criticality

As demand for advanced technologies rises, mineral commodities will increase in geopolitical importance. To assess risks associated with mineral commodity supply chain disruptions, governmental agencies and others have developed "criticality" assessments, with criticality described using the economic impact and probability of supply chain disruptions. In previous work, subjective supply risk indicators were developed to approximate this probability, typically combining several factors such as supply diversity and political stability of trading partners, where indicator weightings can substantially impact results. This work explicitly quantifies trade barrier probability using an ensemble of several machine learning classifiers, with probability estimates informed by exogenous variables such as prior trade barrier implementation and global export dominance. Major differences in the high-probability countries and commodities are observed across models, but the ensemble method highlights Indonesia, China, Tanzania, and the United States as particularly high risk. This approach enables a direct, quantitative, objective approach to assessing trade barrier probability, enhancing risk identification and prioritization for policymakers.

SSRN↗
USGS2026

A geographic dataset of rocky reefs habitat areas of particular concern for the United States West Coast

The United States National Marine Fisheries Service determines “essential fish habitat (EFH)” for federally managed species in coordination with regional fishery management councils, considers adverse effects to those habitats, and provides information to further habitat conservation and enhancement. Identifying discrete subsets of EFH as “habitat areas of particular concern (HAPC)” can help focus conservation, management, and research efforts. In 2006, the Pacific Fishery Management Council designated rocky reefs along the United States (U.S.) West Coast as HAPCs for groundfishes because of their ecological significance, sensitivity to human impacts, and relative rarity. To better understand where rocky reefs occur, we (1) located rocky reef areas that were not included in the 2006 rocky reef dataset, and (2) incorporated best available data into a refined geographic dataset that enables visualization. Our update shows that rocky reefs are distributed throughout the U.S. West Coast continental margin, are patchier than previously known, and comprise 8% of the extent of all data inputs. This updated dataset will inform resource management decisions in coastal and marine environments.

California, Oregon, Washington↗
USGS2025

Over, under, and through: Hydrologic connectivity and the future of coastal landscape salinization

Seawater intrusion (SWI) affects coastal landscapes worldwide. Here we describe the hydrologic pathways through which SWI occurs - over land via storm surge or tidal flooding, under land via groundwater transport, and through watersheds via natural and artificial surface water channels—and how human modifications to those pathways alter patterns of SWI. We present an approach to advance understanding of spatiotemporal patterns of salinization that integrates these hydrologic pathways, their interactions, and how humans modify them. We use examples across the East Coast of the United States that exemplify mechanisms of salinization that have been reported around the planet to illustrate how hydrologic connectivity and human modifications alter patterns of SWI. Finally, we suggest a path for advancing SWI science that includes (a) deploying standardized and well-distributed sensor networks at local to global scales that intentionally track SWI fronts, (b) employing remote sensing and geospatial imaging techniques targeted at integrating above and belowground patterns of SWI, and (c) continuing to develop data analysis and model-data fusion techniques to measure the extent, understand the effects, and predict the future of coastal salinization.

Water Resources Research↗
USGS2026

Tectonic implications of transitional melting regimes from petrological, geochronological, and compositional characterization of the ophiolitic Seventymile terrane, Alaska, USA

New geochemical, U-Pb geochronology, and Sr-Nd-Hf isotope data provide evidence for the tectonic evolution of the Seventymile terrane in interior Alaska, USA. Ultramafic and mafic rocks of the Seventymile terrane are thought to represent components of a dismembered ophiolite and provide unique constraints on regional terrane evolution and accretion. The Seventymile ophiolite represents fragments of the Devonian to Permian Slide Mountain Ocean (SMO) that separated allochthonous and parautochthonous continental fragments of western North America. It now occurs as multiple thrust sheets containing Permian mafic and ultramafic rocks overlying and/or possibly imbricated with amphibolite-facies supracrustal rocks of allochthonous Yukon-Tanana terrane and parautochthonous North America. Seventymile klippen contain variably serpentinized peridotite, primarily harzburgite, low-grade meta-mafic rocks, and minor oceanic sedimentary rocks (argillite, chert, limestone, and metasandstone). Mafic rocks include gabbro to diabase, typically as dikes, veinlets, or rare massive stocks intruding peridotite. Mafic rocks also include greenstones of the Seventymile assemblage in klippen structurally underlying, and in shear zone contact with, ultramafic klippen. New trace element and radiogenic isotope data from mafic magmatic rocks associated with the Seventymile ultramafic bodies show evidence for a weakly subduction-modified mantle source, like the mantle source of normal mid-ocean-ridge basalt (N-MORB) or back-arc basin basalt (BABB). Seventymile assemblage greenstones are more heterogeneous. They range from N-MORB to enriched mid-ocean-ridge basalt (E-MORB) and ocean-island basalt (OIB), with a subset of samples indicative of continental arc affinity. These geochemistry results indicate that distinct tectonic environments are represented by at least two, and possibly three, lithological and structural units comprising the Seventymile terrane. Hf-Nd isotope systematics are consistent with a depleted MORB mantle (DMM)−like component that overlaps with Pacific MORB. Primary zircon is rare, but new in situ U-Pb data for gabbro and greenstone indicate ca. 274−272 Ma peak zircon and titanite crystallization. Scattered younger zircons define a ca. 255 Ma zircon peak and correspond to secondary crystallization associated with baddeleyite reaction of high-Si fluids during low-grade metamorphism. If Seventymile suites are contemporaneous, obduction associated with the closure of the SMO resulted in the stacking of ophiolitic packages representing distinct tectonomagmatic settings across the transition from pericontinental, to epicontinental, to distal ocean back-arc. Intrusions hosted in klippe of ultramafic rocks, plus the least subduction-modified greenstones underlying them, geologically and compositionally resemble Slide Mountain rocks of the Campbell Range formation in eastern Yukon and may provide a new piercing point across the Tintina fault.

Alaska, British Columbia↗
USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗
USGS2025

Stories in stone: The geoheritage significance of U.S. National Park Service stratotypes

The U.S. National Park Service stratotype inventory has systematically documented hundreds of published stratotypes across the country that represent a quintessential component of America’s geoheritage and possess significant scientific, educational, cultural, historic, and aesthetic values. As valuable geologic reference standards, stratotypes record intervals of Earth history relating to our nation’s unique geologic evolution, biological progression, cultural traditions, historical narratives, and stunning landforms. However, the general lack of awareness and knowledge regarding the scope, significance, distribution, and scale of stratotypes can result in these resources being overlooked. Consequently, stratotypes may inadvertently experience accidental alteration, damage, or destruction through both natural and anthropogenic processes. As with many American geoheritage sites, stratotypes lack specific recognition through legislation, policies, or other management strategies to ensure their preservation and protection. Through education and outreach, we hope to inform, increase awareness, and engage the public about these internationally and nationally significant geologic reference sites, unite various geoheritage communities, and promote their preservation for future generations.

Geoheritage↗
USGS2025

Timing of hydrothermal alteration and Au-Sb-W mineralization, Stibnite-Yellow Pine district, Idaho

The Stibnite-Yellow Pine district of central Idaho was mined from the early 1900s until the 1990s, extracting gold, antimony, tungsten, and mercury from veins and disseminated and replacement ores in mountainous terrain along the headwaters of the Salmon River. Mining during the two World Wars supplied critical antimony and tungsten to the war efforts. Recent exploration has delineated mineral resources of over 187 metric tons Au, 274 metric tons Ag, and 93,000 metric tons Sb. Mineralization is hosted in Cretaceous Idaho batholith granitic rocks and a sequence of Neoproterozoic to Paleozoic metasedimentary strata of carbonate and siliciclastic compositions. Historical studies outlined some of the complex paragenesis but debated the absolute age of mineralization. New petrographic and geochronologic work documents a sequence of five hydrothermal events in the Stibnite-Yellow Pine district. Event 1 is related to Cretaceous magmatic and hydrothermal activity and includes events ranging in age from 86 to 75 Ma, including sparse quartz-molybdenite veins dated at 86 Ma. Disseminated gold mineralization of event 2 is associated with sericitic alteration and sulfidation of igneous biotite and replacement of plagioclase by potassium feldspar, largely in granodiorite. Gold is present in zoned arsenian pyrite in both disseminated ores and in crosscutting carbonate-quartz veins containing pyrite and arsenopyrite. The large Yellow Pine deposit, localized at a dilatant bend in the Meadow Creek fault, hosts such disseminated and vein gold. Event 2 is interpreted as the major gold-forming event; 40 Ar/ 39 Ar ages of sericite and potassium feldspar alteration range, respectively, from 70 to 59 and 66 to 56 Ma. The long span is interpreted to reflect the age of gold mineralization and local overprinting by event 3. A narrower range from 66 to 61 Ma is interpreted to date the peak of gold mineralization and alteration. Event 3, tungsten mineralization with scheelite, is texturally later than event 2 gold and localized along the Meadow Creek structure. Event 3 scheelite has been dated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb methods at 57 Ma. Event 4, best developed in the West End area, includes gold-silver–bearing quartz-carbonate-pyrite veins and breccias with epithermal textures and potassium feldspar alteration envelopes. Adularia from event 4 yields 40 Ar/ 39 Ar plateau ages of 52 to 51 Ma. Event 5 antimony and mercury mineralization consists of stibnite veins and breccia cements at the Yellow Pine and Hangar Flat deposits as well as cinnabar veins and replacements at the peripheral Fern and Hermes deposits; it is constrained by an LA-ICP-MS U-Pb date on scheelite (ca. 47 Ma) intergrown with stibnite. Minor propylitic and argillic alteration is evident in 47 Ma igneous dikes, which do not contain economic mineralization. The Au-Sb-W ores in the Stibnite-Yellow Pine mining district formed over an extended time period from about 70 to 45 Ma in multiple pulses that were localized along the Meadow Creek fault zone. Each event corresponds to episodes of magmatism and/or hydrothermal activity in the region. Insignificant event 1 skarn and molybdenum mineralization is similar in age to the Thompson Creek porphyry molybdenum deposit in central Idaho. Event 2 gold mineralization occurred during a magmatic gap in central Idaho but was synchronous with magmatism in the Bitterroot lobe further north; event 2 is similar in age to orogenic gold-arsenic mineralization at the Beartrack mine in eastern Idaho. Event 3 scheelite mineralization coincides with tungsten mineralization at the Quartz Creek deposit, late magmatism in the Bitterroot lobe, and rapid exhumation of the Atlanta lobe of the Idaho batholith. Event 4 gold mineralization is coincident with the onset of regional Challis magmatism and extension. Event 5 antimony and mercury mineralization is time-equivalent to epithermal gold mineralization in the nearby Thunder Mountain volcanic field and the peak of Challis magmatism.

Idaho↗
USGS2026

Sackung at Bald Eagle ridge, central Colorado: An updated interpretation of ridge-spreading movement, structures, and mechanisms from 50 years (1975–2025) of U.S. Geological Survey research

Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.

Colorado↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗
USGS2025

Local, regional, and distal recordings of seismic unrest at Tau Island volcano, American Samoa

A seismic swarm near Taʻū Island, a volcanic island in eastern American Samoa, occurred from July to October 2022. The earliest unrest was noted as felt shaking reports in late July, and instrumentation varied in the beginning of the sequence as the U.S. Geological Survey (USGS) Hawaiian Volcano Observatory responded by installing temporary and then permanent seismometers to monitor the activity. This network variability made it difficult to characterize the earliest seismicity and contextualize the entire sequence to discriminate between an underlying tectonic or volcanic source. Here, we present results analyzing hydroacoustic detections from an International Monitoring System hydrophone array near Wake Island, 4500 km northwest of Taʻū Island volcano. Using least-squares beamforming analysis, we create a catalog of T-wave detections from the direction of Taʻū Island to track the earthquakes, some of which were located by the USGS National EarthquakeInformation Center. Both the rate and hydroacoustic pressures, which we interpret as a proxy for earthquake size, gradually increased from late July to August, peaking on August 19 (rate) and August 24 (size), before decreasing to background in late September. Minutes-long bursts of tremor were also contemporaneously recorded as local network data became avail-able on August 20. Tremor activity continued throughout the rest of August, peaking on August 25, before ending in earlySeptember. These tremor bursts were band-limited to ~ 1–5 Hz and recorded as S waves at a regional station on the island of Upolu in Samoa, 250 km to the west of Taʻū Island. Our results do not constrain the tremor locations, but comparisons of earthquake and tremor reduced displacements recorded locally and regionally suggest a deeper tremor source. We interpret the increase in earthquake size and rate, together with the occurrence, characteristics, and relative depth of the tremor to be the result of magmatic activity beneath Taʻū Island volcano.

American Samoa, Taʻū Island↗
USGS2025

Paleoproterozoic vein graphite mineralization caused by decarbonation in the Ruby Range, Montana, USA

Hydrothermal graphite veins are a possible source for modern battery materials and require better understanding of their carbon source(s) and absolute timing to develop mapable criteria for exploration models. We present new observations of graphite vein and alteration paragenesis and U-Pb LA-ICP-MS titanite age data from the Ruby prospect, Montana, USA, that constrain mineralization timing and source. The graphite veins cut high-temperature metamorphic rocks of the lower Christensen Range suite and are associated with intense diopside (Di0.69Hd0.27Jo0.04) alteration of marble. The oldest titanite ages in fresh marble and partially altered calc-silicate gneiss are ca. 2,500 – 2,450 Ma, show elevated REE values, and Eu/Eu* <1.5, consistent with growth during early regional metamorphism. Titanite in diopside-graphite alteration, interpreted as part of the hydrothermal vein-forming episode, cluster around 1,750 Ma, are characterized by lower REE values, and show Eu/Eu* >2; suggesting a low fO 2 fluid generated from metamorphism during the Big Sky orogeny. Our paragenetic observations and titanite ages indicate graphite vein formation via skarnoid decarbonation reactions in marble late in regional orogenesis. Granulite-facies, carbonate-bearing, supracrustal rock terranes are thus favorable for hydrothermal graphite vein deposits.

Montana↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS2024

U.S. Geological Survey Karst Interest Group proceedings, Nashville, Tennessee, October 22-24, 2024

Karst hydrogeologic systems represent challenging and unique conditions to scientists studying groundwater flow and contaminant transport. Karst terrains are characterized by distinct and beautiful landscapes, caverns, and springs, and many of the exceptional karst areas are designated as national or state parks. The range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than any other aquifer type. The U.S. Geological Survey (USGS) Karst Interest Group (KIG), formed in 2000, is a loosely knit, grassroots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst aquifers. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. To accomplish its mission, the KIG has organized a series of workshops. To date (2024), nine KIG workshops, including the workshop documented in this report, have been held. The abstracts and extended abstracts provide a snapshot in time of past and current karst related studies. The USGS Water Availability and Use Science Program funded the workshop and proceedings. The planning committee for the ninth workshop includes Thomas D. Byl (USGS and Tennessee State University), Allan K. Clark (USGS), Laura M. DeMott (USGS), Eve L. Kuniansky (USGS, Emeritus), Benjamin V. Miller (USGS), and Lawrence E. Spangler (USGS, Emeritus). The workshop proceedings are edited by Eve L. Kuniansky and Lawrence E. Spangler. The field trip guide was produced by Benjamin V. Miller and Brian Ham (Tennessee Department of Environment and Conservation) and included in the proceedings from the KIG’s 2021 virtual workshop to be used on the optional field trip held on Thursday, October 24, 2024.

Open-File Report↗
USGS2025

Properties of solid bitumen formed during hydrous, anhydrous, and brine pyrolysis of oil shale: Implications for solid bitumen reflectance in source-rock reservoirs

Solid bitumen reflectance (BR o in %) is widely used as a thermal maturity proxy in source-rock reservoirs, yet solid bitumen texture and physical attributes may be affected by varying environmental constraints during its formation, e.g., water saturation, mineral catalysis, or salinity. Herein we investigated the development of solid bitumen properties during artificial maturation using three diverse (lacustrine to marine) oil shale samples containing abundant amorphous H-rich sedimentary organic matter (bituminite). The oil shales were treated via pyrolysis (320 °C, 72 h) using hydrous, anhydrous, and brine conditions, causing the development of a newly formed solid bitumen in the experiment residues from the thermal conversion of bituminite. The properties of the newly formed solid bitumen then were evaluated via geochemical screening tests, optical and electron microscopy, and infrared spectroscopy. Experimental residues also were treated via solvent extraction, allowing characterization of the effects of extraction on solid bitumen. Results showed that bituminite with higher resin and asphaltene components evolved to a solid bitumen with higher reflectance (as a percentage of the original value) when exposed to the same heat treatment, relative to bituminite with higher saturate and aromatic components. Aromatization of solid bitumen also was generally more pronounced in the presence of deionized water relative to anhydrous conditions, supporting prior observations. These results suggest the compositions of primary sedimentary organic matter, and the local concentration of water, affect the development of solid bitumen during thermal advance, potentially explaining the origin of multiple solid bitumen populations with varying reflectance in source-rock reservoir samples or in a single microscope field. Experiments using brine were inconclusive regarding enhanced/decreased aromatization, which could be attributed to the salinities of the brines used in the tests. Extraction of residues caused a consistent increase in BR o values, suggesting that migrating or expelled oils could cause an increase in BR o via natural solvation and absent an increase in temperature. This work provides insights into the development and evolution of BR o in source-rock reservoirs as a function of the composition of the original bituminite and changing environmental conditions, with potentially broad application for petroleum prospecting and resource estimation.

Colorado↗