SEARCH · USGS Science

Results for “PLoS ONE”

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

SourceResearch connectionFrom yearThrough year

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

1,186 records · Page 62Linked to original sources

USGS2025

Submarine groundwater discharge creates cold‐water refugia that can mitigate exposure of heat stress in nearshore corals

Coral reef mortality around the world is accelerating due to human activities and rising sea temperatures that cause bleaching, which is expected to become more frequent. Our ability to predict which corals will be most resilient, however, remains limited due to insufficient information characterizing nearshore temperature and habitat conditions. In this study, we examine how submarine groundwater discharge (SGD) reduces nearshore water temperatures and exposure of corals to heat stress, complementing the understanding that SGD can adversely affect coral when it contains elevated nutrient concentrations. Data from fixed nearshore sensors and vertical depth profiles along ~100 km of the western shoreline of the Island of Hawai’i from 2003 to 2014 demonstrate that submarine groundwater discharge (SGD) can reduce nearshore water temperatures by 1 °C–5°C and create estuarine-like conditions with salinities as low as 20 PSU, where the prevalent coral species, Pocillopora meandrina, Porites lobata , and Montipora capitata , thrive. Time-series temperature records reveal that exposure to high ambient ocean temperatures, which are known to initiate bleaching events, are reduced up to 5%–46% of the time. Coral health surveys indicated coral bleaching in response to moderately high annual temperatures in 2010 and 2011, with more colonies affected farther from cold, SGD-fed waters. Synthesis of these results, along with coral response data following the more extreme marine heat wave of 2014–2015, demonstrates lower coral loss and greater coral recovery near groundwater seeps, particularly those with higher flux and influence on reducing nearshore water temperatures. Our results demonstrate that SGD may therefore provide a beneficial ecosystem service and enhance coral reef resilience, particularly where human-related nutrient additions to groundwater can be mitigated. The implications of our findings are relevant across tropical coasts where groundwater inputs can be substantial, such as the Caribbean and Indo-Pacific, and contribute to improving our understanding of coral sensitivity to gradients in temperature and nutrient stress. Improved management of groundwater resources could thus be vital to local–regional strategies for mitigating future heat stress.

Hawaii↗

USGS2026

Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA

Methane is a potent greenhouse gas that plays an important role in atmospheric chemistry and global warming. The current global methane budget has large uncertainties, and a better understanding of the budget would help to guide strategies for reducing anthropogenic emissions to fight climate change. Natural geologic methane emissions are a particularly poorly constrained source, with top-down estimates from 14 C in ice cores suggesting much lower geologic emissions than bottom-up scaling of direct flux measurements. Our study aims to contribute to resolving this discrepancy through improved bottom-up characterization of geologic methane seepage in the San Juan Basin in southwestern Colorado and northwestern New Mexico, USA. We performed 983 new flux chamber measurements in this basin during summer 2022 and winter 2023 field campaigns. Our results, in combination with prior measurements, suggest that natural seepage in the San Juan Basin only occurs on or near the Fruitland coal outcrop. Specifically, our new measurements confirm previous measurements of seepage along the northwestern exposure of the Fruitland outcrop in Colorado (a known hydrodynamic overpressure region) and for the first time, identified seepage locations along the southernmost Fruitland outcrop exposure in New Mexico, in association with a coal cleat and a fault. Overall, seepage along the Fruitland coal outcrop is heterogeneously distributed, with both positive and negative (interpreted as microbial soil sink) methane fluxes. Features that are hypothesized to be predictive of seepage (e.g., faults) were not associated with positive methane fluxes in areas outside of the Fruitland outcrop. Our best estimate for total geologic methane seepage in the San Juan Basin from spatial interpolation and statistical upscaling is approximately 0.14 Tg CH 4 /yr, with a range from 0.029 to 0.48 Tg CH 4 /yr. This best-estimate value is lower than a previous bottom-up estimate from a gridded seepage inventory, but higher than a previous top-down estimate.

Colorado, New Mexico↗

USGS2025

Land cover change within wetland complexes at Dixie Meadows, Churchill County, Nevada: 2015 – 2023

Dixie Meadows, Nevada, is a system of geothermal springs and seeps that feed a complex of marshes and wetland meadows that are located within lands managed by the Bureau of Land Management (BLM) and the Department of Defense (DOD). A previous U.S. Geological Survey report documented variability in satellite imagery-based land cover classifications for seven wetland complexes at near monthly time intervals between October 2015 and January 2022. This report presents additional data, extending analysis to November 2023. Land cover classifications between October 2015 and November 2023 demonstrated an association between vegetation cover characteristics and surface moisture, with Class 1 having dry, bare soil or sparse upland vegetation, Class 2 having moist, bare soil or sparse to small vegetation, Class 3 having dense green vegetation with potentially saturated soil conditions, Class 4 having a mix of shallow surface water, saturated soil, and dense green vegetation, and Class 5 having open surface water. Most of the wetland complexes occur close to spring outflows primarily within land managed by the DOD, though portions are also within BLM lands. The intervening and surrounding landscape outside of the wetland complexes assessed in this study are managed by the BLM. As a result, Class 1 land covers had the largest areal coverage for BLM managed lands. Classes 2 and 3 land covers were primarily mapped inside the wetland complexes and thus had the largest area coverage within DOD managed lands. Class 4 was almost exclusively mapped within the wetland complexes and thus was largely contained within DOD managed lands. Class 5 (open water) was exclusively mapped in and adjacent to a single wetland complex with catchment ponds on land managed by the BLM. The distribution of these land cover classes over the study period was seasonally and annually variable. Land cover areas of Classes 1 and 2 were larger during the spring months. Conversely, land cover areas of Classes 3 and 4 tended to be greatest during the summer or fall. These patterns might be influenced by differences in seasonal water sources and phenology.

Nevada↗

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↗

USGS1998

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Cecil Field Naval Air Station, Jacksonville, Florida

As part of the Installation Restoration Program, Cecil Field Naval Air Station, Jacksonville, Florida, is considering remedialaction alternatives to control the possible movement of contaminants from sites that may discharge to the surface. This requires a quantifiable understanding of ground-water flow through the surficial aquifer system and how the system will respond to any future stresses. The geologic units of interest in the study area consist of sediments of Holocene to Miocene age that extend from land surface to the base of the Hawthorn Group. The hydrogeology within the study area was determined from gamma-ray and geologists' logs. Ground-water flow through the surficial aquifer system was simulated with a seven-layer, finite-difference model that extended vertically from the water table to the top of the Upper Floridan aquifer. Results from the calibrated model were based on a long-term recharge rate of 6 inches per year, which fell in the range of 4 to 10 inches per year, estimated using stream hydrograph separation methods. More than 80 percent of ground-water flow circulates within the surficial-sand aquifer, which indicates that most contaminant movement also can be expected to move through the surficial-sand aquifer alone. The surficial-sand aquifer is the uppermost unit of the surficial aquifer system. Particle-tracking results showed that the distances of most flow paths were 1,500 feet or less from a given site to its discharge point. For an assumed effective porosity of 20 percent, typical traveltimes are 40 years or less. At all of the sites investigated, particles released 10 feet below the water table had shorter traveltimes than those released 40 feet below the water table. Traveltimes from contaminated sites to their point of discharge ranged from 2 to 300 years. The contributing areas of the domestic supply wells are not very extensive. The shortest traveltimes for particles to reach the domestic supply wells from their respective contributing areas ranged from 70 to 200 years.

Florida↗

USGS2025

Relating systematic molecular and textural properties of graptolite pyrolyzed via gold tube hydrous pyrolysis: Implications for thermal proxies in lower Paleozoic marine shales

A series of gold tube pyrolysis experiments (72 h, 300–550 °C, 50 MPa) conducted on a graptolite-rich lower Paleozoic marine shale generated pyrolysis residues for a comprehensive evaluation of the molecular and structural variability of three types of graptolite periderm. Organic petrology, Raman spectroscopy, and field emission scanning electron microscopy (FE-SEM) with energy dispersive spectroscopy (EDS) were combined to evaluate the thermal evolution process. The three types of graptolite periderm, namely granular, non-granular, and nodular graptolite, were analyzed by Raman spectroscopy wherein point measurements were obtained after the maceral was identified and the location verified by organic petrology. Distinct thermal evolution pathways among non-granular, granular, and nodular graptolite periderms were recorded. The evolution patterns of the Raman parameters, particularly D1 and G bands, highlight the differences in geochemical composition of the graptolite periderm types and the alteration of molecular structure with increasing thermal maturity. Raman parameters D1 (position of the D1 peak), G-FWHM (full width at half maximum of the G peak), and ratios D1-FWHM/G-FWHM (full width at half maximum of the D1 peak ratioed to G-FWHM) and A D1 /A G (ratio of D1 and G peak intensities) showed effectiveness in assessing thermal maturity. Bireflectance with increasing gold tube pyrolysis temperature followed a hierarchy: non-granular > granular > nodular, reflecting different molecular alignment intensities. Qualitative FE-SEM evaluation showed that fine-grained mineral inclusions (primarily Fe-sulfide as determined via EDS) were associated with the graptolite populations, with granular graptolite containing greater amounts of coarser-grained (e.g., ∼300–1400 nm) mineral inclusions relative to non-granular and nodular graptolite, which contain finer-grained (e.g., ∼100–200 nm) inclusions difficult to resolve with optical microscopy. These findings are investigated to highlight the mechanisms that drive organic matter evolution within graptolite during thermal maturation, as well as to explore some of the limitations of using spectroscopic parameters as thermal maturity proxies.

International Journal of Coal Geology↗

USGS2025

Decomposing the Tea Bag Index and finding slower organic matter loss rates at higher elevations and deeper soil horizons in a minerogenic salt marsh

Environmental gradients can affect organic matter decay within and across wetlands and contribute to spatial heterogeneity in soil carbon stocks. We tested the sensitivity of decay rates to tidal flooding and soil depth in a minerogenic salt marsh using the Tea Bag Index (TBI). Tea bags were buried at 10 and 50 cm depths across an elevation gradient in a subtropical Spartina alterniflora marsh in Georgia (USA). Plant and animal communities and soil properties were characterized once, while replicate tea bags and porewaters were collected several times over 1 year. TBI decay rates were faster than prior litterbag studies in the same marsh, largely due to rapid green tea loss. Rooibos tea decay rates were more comparable to natural marsh litter, potentially suggesting that is more useful as a standardized organic matter proxy than green tea. Decay was slowest at higher marsh elevations and not consistently related to other biotic (e.g., plants, crab burrows) or abiotic factors (e.g., porewater chemistry), indicating that local hydrology strongly affected organic matter loss rates. TBI rates were 32 %–118 % faster in the 10 cm horizon than at 50 cm. Rates were fastest in the first 3 months and slowed 54 %–60 % at both depths between 3 and 6 months. Rates slowed further between 6 and 12 months, but this was more muted at 10 cm (17 %) compared to 50 cm (50 %). Slower rates at depth and with time were unlikely due to the TBI stabilization factor, which was similar across depths and decreased from 6 to 12 months. Slower decay at 50 cm demonstrates that rates were constrained by environmental conditions in the deeper horizon rather than the composition of this highly standardized litter. Overall, these patterns suggest that hydrological setting, which affects oxidant introduction and reactant removal and is often overlooked in marsh decomposition studies, may be a particularly important control on organic matter loss in the short term (3–12 months).

Georgia↗

USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗

USGS2024

Plan to coordinate post-earthquake investigations supported by the National Earthquake Hazards Reduction Program (NEHRP)

Introduction This report presents a plan supported by the National Earthquake Hazards Reduction Program (NEHRP) to coordinate domestic and international post-earthquake investigations (herein called “the Plan”). Post-earthquake scientific and engineering investigations are undertaken to capture critical information to understand the causes and impacts of the event, lessons from which can substantially improve the Nation’s resilience after future earthquakes. NEHRP is the Federal Government’s coordinated nationwide program to reduce risks to life and property from earthquakes. The Plan describes the activation and coordination of the four designated NEHRP Agencies in the Federal Government: Federal Emergency Management Agency, National Institute of Standards and Technology, U.S. National Science Foundation, and U.S. Geological Survey (USGS). The Plan also describes coordination between NEHRP Agencies and other organizations that may participate in pre-event and post-earthquake investigations, including non-NEHRP Federal agencies; State, regional, local, Tribal, and territorial agencies; domestic nongovernmental organizations; academic institutions and affiliated organizations; private companies; foreign governmental agencies and nongovernmental organizations; and international organizations. The Plan delineates the coordination of NEHRP post-earthquake scientific and engineering investigations to document the direct, indirect, and cascading physical and societal impacts from fault rupture and ground shaking hazards and from secondary hazards such as landslides, liquefaction, and tsunamis. In addition, the Plan identifies pre-event activities necessary to ensure that post-earthquake investigations are executed effectively. The USGS is the lead NEHRP Agency for activating and coordinating NEHRP post-earthquake investigations and for implementing this Plan. The USGS also leads coordination of the NEHRP Agencies in completing the pre-event activities identified in the Plan. The Plan has new information and supersedes USGS Circular 1242, “The Plan to Coordinate NEHRP Post-Earthquake Investigations,” which was published in 2003. This second-generation Plan was developed with the assistance of the Applied Technology Council (ATC) of Redwood City, California, under USGS contract 140G0121P0309, ATC-155. A 17-member Project Review Panel provided guidance on plan development, and input was solicited from subject matter experts representing key stakeholder groups and from participants at a public workshop.

Circular↗

USGS2026

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin↗

USGS2025

Determination of per- and polyfluoroalkyl substances in water by direct injection of matrix-modified centrifuge supernatant and liquid chromatography/tandem mass spectrometry with isotope dilution

A direct-injection liquid chromatography/tandem mass spectrometry method was developed to determine 34 per- and polyfluoroalkyl substances (PFAS), including selected branched isomers, in centrifuge supernatant of matrix-modified (amended with approximately 50 percent methanol) water samples. The method has been validated in reagent water, surface water, groundwater, and wastewater effluent. Other water types (for example, drinking water, untreated wastewater, and landfill leachate) have been analyzed by the method but not systematically validated. Recovery of isotope-dilution standards, added to each sample, may be used to assess method performance in nonvalidated matrices on a sample-by-sample basis. Using this method, PFAS concentrations were determined in the range of 2–2,000 nanograms per liter in water samples. This range can be extended by diluting concentrated samples. At circumneutral pH, most compounds are present in the environment in their ionized form, and data are reported as such (for example, perfluorooctanoic acid is referred to as “perfluorooctanoate” [PFOA], perfluorooctane sulfonic acid is referred to as “perfluorooctane sulfonate” [PFOS]). Sample preparation procedures were designed without the use of filtration and with minimum sample handling steps to mitigate procedural losses of target compounds due to sorption to surfaces. Further, isotope-dilution quantification allowed for the correction of bias that may result from procedural losses, matrix-induced signal suppression or enhancement, and other factors. Validation experiments to characterize bias and variability, method detection level, and holding time were done in four distinct water matrices—reagent water, surface water, treated wastewater effluent, and groundwater—at multiple concentration levels. Mean PFAS recoveries met data quality objectives of bias and variability studies in all four validation matrices except for two compounds with low and variable recovery in the reagent water matrix only. Isotope-dilution standards, treated as surrogate compounds, were analyzed in more than 1,500 customer-submitted environmental samples with aggregate recovery of 102.5±6.5 percent (mean±standard deviation). Maximum holding times for all target compounds in the four validation matrices were 28 days for refrigerated samples and 90 days for frozen samples.

Techniques and Methods↗

USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗

USGS2024

Distribution of ancient carbon in groundwater and soil gas from degradation of petroleum near the Red Hill Bulk Fuel Storage Facility, O‘ahu, Hawai‘i

The groundwater below the Red Hill Bulk Fuel Storage Facility (the facility) in Oʻahu, Hawaiʻi, contains fuel compounds from past spills. This study used carbon-14 analyses to distinguish fuel-derived carbon from background carbon, along with other biodegradation indicators, to address two goals: (1) determine the extent and migration direction of groundwater affected by residual fuel below the facility and (2) determine if residual fuel locations in the subsurface could be identified by analyzing soil gas at the surface above the facility. Groundwater from 19 wells was sampled between September 2022 and April 2023. Nonvolatile dissolved organic carbon (NVDOC) from a well presumed to be unaffected by past spills contained 38 percent ancient carbon indicating a natural source of ancient carbon in the subsurface. The NVDOC concentrations and ancient carbon percentages indicate fuel biodegradation products are likely present on the north and south of Red Hill with the greatest effects at well RHMW02 near the 2014 spill site. The NVDOC concentrations are almost three times higher than diesel range organic (DRO) concentrations in groundwater from the same sites. Major ion data indicate that iron reduction is an important biodegradation process. Soil probe samples and soil carbon traps were used to determine the carbon-14 content of soil carbon dioxide. Ancient carbon from fuel biodegradation was not detected at any soil probe or carbon trap site in contrast to a 2017 study which reported ancient carbon detections. A reanalysis of the 2017 results using a range of local values for background carbon-14 indicates that ancient carbon from fuel biodegradation was probably only detected in lower tunnel exhaust system samples and not in any soil carbon trap samples. Measurements of carbon dioxide efflux with a dynamic closed chamber were highly variable. The soil gas results indicate that soil gas measurements at land surface were not useful for detecting residual fuel at the facility.

Hawaii↗

USGS2025

Ungulate migrations of the Western United States, volume 5

Many ungulates migrate between distinct summer and winter ranges to take advantage of spatially and temporally variable food sources and avoid threats such as predators and deep snow. In 2018, the U.S. Department of the Interior established Secretarial Order 3362, which provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations by enhancing habitat quality for ungulates across the Western United States. In response to the order, the U.S. Geological Survey (USGS) created the Corridor Mapping Team, which is a collaboration among 11 State agencies, regional and Federal partners, and an expanding number of Tribal wildlife agencies. Together, the Corridor Mapping Team maps ungulate migrations throughout the Western United States and publishes them in the USGS “Ungulate Migrations of the Western United States” report series. This report details migrations and seasonal ranges from 36 additional herds and includes 2 herd updates detailed in previous reports. The Corridor Mapping Team has mapped the migrations and seasonal ranges of 218 unique herds for the report series, including this report. The report series serves as a map-based inventory of the ungulate migrations across the Western United States for biologists, managers, policymakers, and conservation practitioners. Building on the previous report volumes in the series, volume 5 additionally describes some of the local and national initiatives that are incorporating the products, tools, and information from this growing USGS report series.

Navajo Nation↗

USGS2026

Magmatic source of the opening phase of Kīlauea’s 2018 Lower East Rift Zone eruption

The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea.

Hawaii↗

USGS2026

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper↗

USGS2024

Layered intrusions in the Precambrian: Observations and perspectives

Layered intrusions are plutonic bodies of cumulates that form by the crystallization of mantle-derived melts. These intrusions are characterized by igneous layering distinguishable by shifts in mineralogy, texture, or composition. Layered intrusions have been fundamental to our understanding of igneous petrology; however, it is their status as important repositories of critical metals – such as platinum-group elements, chromium, and vanadium – that has predominantly driven associated research in recent decades. Many layered intrusions were emplaced during the Precambrian, predominantly at the margins of ancient cratons during intervals of supercontinent accretion and destruction. It appears that large, layered intrusions require rigid crust to ensure their preservation, and their geometry and layering is primarily controlled by the nature of melt emplacement. Layered intrusions are best investigated by integrating observations from various length-scales. At the macroscale, intrusion geometries can be discerned, and their presence understood in the context of the regional geology. At the mesoscale, the layering of an intrusion may be characterized, intrusion-host rock contact relationships studied, and the nature of stratiform mineral occurrences described. At the microscale, the mineralogy and texture of cumulate rocks and any mineralization are elucidated, particularly when novel microtextural and mineral chemical datasets are integrated. For example, here we demonstrate how mesoscale observations and microscale datasets can be combined to understand the petrogenesis of the perplexing snowball oiks outcrop located in the Upper Banded Series of the Stillwater Complex. Our data suggest that the orthopyroxene oikocrysts did not form in their present location, but rather formed in a dynamic magma chamber where crystals were transported either by convective currents or within crystal-rich slurries. Critical metals may be transported to the level of a nascent intrusion as dissolved components in the melt. Alternatively, ore minerals are entrained from elsewhere in a plumbing system, potentially facilitated by volatile-rich phases. There are many ore-forming processes propounded by researchers to occur at the level of emplacement; however, each must address the arrival of the ore mineral, its concentration of metals, and its accumulation into orebodies. In this contribution, several of these processes are described as well as our perspectives on the future of layered intrusion research.

Precambrian Research↗

USGS2026

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗