SEARCH · USGS Science

Results for “ES&T Water”

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.

801 records · Page 11Linked to original sources

USGS2026

Cambrian and Ordovician stratigraphy, conodont biostratigraphy, and microfacies analysis to support 1:24,000-scale geologic mapping of the southern Lake Champlain valley, New York and Vermont

Introduction Geologic mapping in the southern Lake Champlain valley of New York and Vermont (fig. 1) has required evaluation of stratigraphic nomenclature used since the early 1900s. The paleogeography of the Cambrian and Ordovician adjacent to the Adirondack Highlands, the high-relief mountains of Proterozoic igneous and metamorphic rocks of northern New York, has impacted the lithostratigraphic succession, facies changes, and several hiatuses that occur within the package of clastic and carbonate rocks that demonstrates overall deepening consistent from the southern to northern Appalachians. The stratigraphic nomenclature used for the geologic mapping of the Cambrian and Ordovician strata of the study area is derived from evaluation of names put forth by previous workers from areas of northern New York and western Vermont. Descriptions of type localities or type sections and other criteria set forth by the North American Stratigraphic Code (North American Commission on Stratigraphic Nomenclature [NACSN], 2021) were consulted. The North American Stratigraphic Code states that lithostratigraphic units, such as the fundamental unit, the formation, be defined by lithic character and be mappable (NACSN, 2021, Article 24). However, some early workers (Walcott, 1912; Kay, 1937; Fisher and Hanson, 1951; Oxley and Kay, 1959) have defined stratigraphic units based on fossil content, which may or may not be facies dependent, and hence may not be mappable as lithostratigraphic units. Although preservation of older names takes priority in stratigraphic nomenclature (NACSN, 2021, Article 7), some do not have adequately defined stratotypes or are located a distance away with different lithologies than that in the study area. Biostratigraphy can be a helpful tool for understanding correlations and facies changes. Herein, conodont biostratigraphy is utilized to date and correlate separate lithologic units. Further, microfacies analyses of units through inspection of petrographic thin sections provide paleogeographic information that helps in understanding the genesis of geologic units.

New York, Vermont↗

USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗

USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗

USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗

USGS2025

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut↗

USGS2026

Ten new insights in climate science 2025

Interdisciplinary understanding is vital for delivering sound climate policy advice. However, navigating the ever-growing and increasingly diverse scholarly literature on climate change is challenging for any individual researcher. This annual synthesis highlights and explains recent advances across a variety of fields of climate change research. This year, the 10 insights focus on: (1) the record-warmth of 2023/2024 and the elevated Earth energy imbalance; (2) acceleration of ocean warming and intensifying marine heatwaves; (3) northern land carbon sinks under strain; (4) reinforcing feedback between biodiversity loss and climate change; (5) accelerated depletion of groundwater; (6) global dengue incidence; (7) global income losses and labour productivity declines; (8) strategic scaling of CDR; (9) integrity challenges in carbon credit markets and emerging responses; and (10) effective policy mixes for emissions reductions. The insights have been written to be accessible to researchers from different fields, serving as entry-points to specific topics, as well as providing an overview of the evolving landscape of climate change research. In the final section, the insights are used to develop overarching policy-relevant messages. This paper provides the basis for a science-policy report that was shared with all Party delegations ahead of COP30 in Belém, Brazil.

Global Sustainability↗

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↗

USGS2025

Exploration for blind geothermal systems in the eastern Great Basin of Utah: An update on the “Lund North” INGENIOUS detailed study site

Existing geothermal production in Utah is commonly collocated with surficial expressions of geothermal heat including active hot springs and hot spring deposits. However, geothermal potential across the Great Basin region is thought to be much higher for hidden or blind geothermal systems. Accordingly, exploration techniques that can locate geothermal resources that lack surface thermal features could support future development of these systems. The goal of the INGENIOUS project is to reduce exploration risk and discover new, economically viable hidden geothermal systems in the Great Basin region. This paper summarizes the efforts and preliminary results of blind geothermal resource prospecting in Utah as part of the larger INGENIOUS project. The Lund North site in the Basin and Range province of southwest Utah was designated as the fourth detailed study site for the project. The site was initially identified from the preliminary Play Fairway Analysis (PFA) geothermal favorability map, which is based on various data elements, and the identification of favorable structural settings across Utah. This site contains existing legacy data including a thermal gradient borehole with a heat-flow value exceeding 220 mW/m 2 . New data collected at this site include terrestrial gravity, magnetotellurics, transient electromagnetics, fluid geochemistry, geologic mapping, UAV-based lidar, paleomagnetism, aeromagnetic surveys, and a shallow temperature survey. Preliminary data and model interpretations corroborate the presence of a large, east-dipping, northeast-southwest-trending normal fault on the edge of a typical Basin and Range graben. Quaternary fault mapping refined by high-resolution lidar surveys suggests a complex step over fault geometry. Utilizing multi-disciplinary datasets for further site assessment at Lund North will facilitate the development of a local-scale PFA and geothermal conceptual model. This new data will be used to site successive thermal gradient drilling to verify heat flow with temperatures at depth at the Lund North site and further investigate the existence of a local geothermal system.

Utah↗

USGS2024

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington↗

USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Mauna Loa volcano erupts crystal-poor material at its summit and more crystal-rich material on its rift zones. Some of the more olivine-rich lava flows contain xenoliths with diverse mineralogy, including cumulate harzburgites with high-Mg# orthopyroxenes and high-Fo olivines (both > 84). Previous experimental work and thermodynamic modelling has proposed that high-Mg# orthopyroxenes only crystallize from Mauna Loa melts at high pressures (> 6 kbar, > 20 km), leading to suggestions that there is a region of sub-Moho magma storage at Mauna Loa in addition to the geophysically imaged magma reservoir at 2–5 km depth below the summit. We use melt and fluid inclusion barometry combined with thermodynamic models to further investigate this suggestion. Fluid inclusion data from harzburgites and dunitic xenoliths yield storage depths remarkably similar to those found in non-xenolithic crystals from lavas and tephras, with a clear peak at ~ 2–3 km (below the summit). Depths from melt inclusions in these xenoliths overlap with fluid inclusion pressures, ruling out the possibility of fluid inclusion re-equilibration during a period of stalling in a shallower reservoir. We examine five different thermodynamic models and find that the minimum pressure of olivine-orthopyroxene co-saturation varies by ~ 4 kbar (~ 12 km). These models also fail to predict that orthopyroxene is stable in ~ 15–80% of compositionally relevant experimental charges which grew orthopyroxene. Overall, this shows that phase stability modelling is an unreliable method of determining magma storage depth at Mauna Loa. We suggest that model discrepancies reflect a lack of experimental constraints on orthopyroxene stability at > 1200 ℃ and 0.01–5 kbar. Based on the presence of large oikocrystic orthopyroxenes completely enclosing rounded olivine chadacrysts, we suggest that these harzburgitic xenoliths formed through the reaction of intruding melts with olivine mush piles within the Mauna Loa edifice at ~ 3 km depth below the summit, with no need for a deeper storage reservoir. The predominance of pre-eruptive shallow storage means that there is more chance of detecting reservoir destabilization with geophysical monitoring techniques compared to a scenario where melts are supplied from sub-Moho reservoirs.

Hawaii↗

USGS2026

Unprecedented burning in tropical peatlands during the 20th century compared to the previous two millennia

Tropical peatland wildfire incidence has risen in recent decades, driven by drainage for land use and intensified by severe droughts with global climate change. These disturbances have altered vegetation structure, disrupted ecosystem functioning, and increased carbon emissions, particularly in Southeast Asia. However, the long-term history and characteristics of wildfires in tropical peatlands remain largely unknown. Here, we compiled fifty-eight macro-charcoal records from peatlands across the tropics, ranging from lowland forested to montane peatlands, to assess millennia-scale changes and controlling factors of tropical peatland burning. We divided the datasets into four main sub-regions: Neotropical, Afrotropical, Indomalayan and Australasian ecoregions to explore regional variability. Tropical peatlands had high burning levels between 0 and 850 ce , followed by a relatively low and stable period until a marked increase during the 20th century. The general trend in tropical peatland burning follows changes in global temperature, and climate variables that control the length and severity of drought events have a notable influence on peat burning before 1900 ce . During the 20th century, regional differences were observed, with declining fire trends in the Neotropical and Afrotropical regions and increasing fire trends in the Indomalayan and Australasian regions. This difference is likely attributable to human activities, and such intervention is also evident in palm swamps and hardwood swamps under similar wet, weakly seasonal climates. With the increase in anthropogenic pressures on peatlands and greater climate variability, future wildfires in peatlands are likely to become more frequent and widespread across all tropical ecoregions. Conservation and sustainable land-use practices could be used to mitigate and control peatland burning and protect these carbon-rich sinks.

Global Change Biology↗

USGS2026

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Ungulate migrations of the Western United States, volume 6

This report, volume 6 in the “Ungulate Migrations of the Western United States” report series, showcases the migrations of 23 ungulate herds in the Western United States. The report series is produced by the Corridor Mapping Team (CMT). Led by the U.S. Geological Survey, the CMT is a collaboration among 11 State agencies, as well as regional and Federal partners, and an expanding number of Tribal wildlife agencies. The CMT was initiated in response to the U.S. Department of the Interior Secretarial Order 3362, which was signed in 2018 and provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations throughout the Western United States. Including this volume, the report series has detailed the migrations of 237 unique ungulate herds throughout the Western United States and continues to serve as a valuable resource to guide local and regional management, policy, and on-the-ground work necessary to maintain intact and functional ungulate migrations. This report highlights several guiding principles of the CMT that facilitate collaboration among the diverse set of partners and contribute to the program’s continued successes. Notably, raw global positioning system data are not shared among participating agencies and the U.S. Geological Survey, delineating migration corridors and seasonal ranges relies on empirical data, the CMT provides flexible approaches to participating State and Tribal partners, and regular CMT meetings create a framework for open communication among agency partners that supports transboundary mapping of migrations. The 237 ungulate migrations that have been included in the report series are an expanding inventory, which can help maintain ungulate migrations in perpetuity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2026

The value of electronic tagging and tracking studies for understanding fish–hypoxia interactions

Objective Hypoxia (i.e., low dissolved oxygen availability) is a natural phenomenon but can also be induced by human activities (e.g., nutrient enrichment from runoff). Given that dissolved oxygen is essential for aquatic life, periods of hypoxia tend to have negative consequences (e.g., sublethal disturbances, mortality) for most fishes. Extensive laboratory research has documented hypoxia thresholds and physiological and behavioral outcomes for a variety of freshwater and marine fishes and there is also an extensive body of fieldwork assessing population-level responses (e.g., survival, distribution). Yet, studying how individual fish respond to hypoxia in the wild has proved challenging; electronic tagging and tracking tools (e.g., biotelemetry, biologging) have made it easier to study individual responses to hypoxia and complement other tools like hydroacoustics that tend to focus on population-level responses. Methods We review what has been learned from contemporary studies that employ electronic tagging and tracking tools to understand how fish respond to hypoxia in the wild. Topics explored include identifying and validating hypoxia thresholds, habitat compression, connectivity, mortality, physiological and bioenergetic consequences, and extreme weather conditions. We also consider what we have learned about evaluating various management plans for hypoxia and reflect on how electronic tags have also been used in aquaculture systems and in hybrid studies that combine laboratory and field research. We highlight fishes in Lake Erie as a research narrative to demonstrate how electronic tagging and tracking have markedly improved our understanding of a longstanding hypoxia issue. Results Our synthesis revealed that electronic tagging and tracking have provided critical information on how fish respond to hypoxia in the field, revealing complex trade-offs and compensatory mechanisms as well as cryptic hypoxia-induced mortality. Beyond just illuminating space use and mortality, tags equipped with various sensors are revealing how fish deal with hypoxia in real time in terms of physiology, bioenergetics, and behavior. Conclusions As electronic tagging and tracking methods experience further innovation and are increasingly applied to understand the effects of hypoxia on fish, we expect more unanticipated findings about the effects of hypoxia on fish in all aquatic ecosystems, which will strengthen our ability to manage and mitigate hypoxia. Combining tools and approaches (e.g., lab and field) is perhaps the best way to generate comprehensive understanding.

Transactions of the American Fisheries Society↗

USGS2025

The Hardscrabble Creek complex: A newly discovered, mostly buried, Mesoproterozoic mafic-ultramafic pluton in the Wet Mountains, Colorado, USA

The origin of prolific ca. 1.4 Ga ferroan magmatism between the southwestern USA and eastern Canada is enigmatic and has been explained by various models, including extensional, mantle plume, and convergent plate-margin models. Rare mafic plutons are associated with the ferroan plutons, which may help constrain their mantle source and tectonic setting. In the southwestern USA, only two such mafic plutons are known to exist. We present the first evidence for a third, mostly buried, potentially layered, mafic-ultramafic Mesoproterozoic pluton, informally referred to as the Hardscrabble Creek complex, in the central Wet Mountains of Colorado, USA. Recent geophysical data show an elliptical magnetic and gravity high spatially coincident with local gabbroic outcrops. New field and petrographic analyses of these exposed rocks reveal that they consist of ultramafic to mafic cumulates, including orthopyroxenite, olivine norite, norite, and anorthosite. High-precision U-Pb dating of zircon from orthopyroxenite and norite yield weighted mean 206 Pb/ 238 U dates of 1352.36 ± 1.60 Ma and 1352.37 ± 1.71 Ma, respectively. These dates indicate that the complex formed over a narrow timeframe, after the adjacent 1362 ± 7 Ma ferroan San Isabel Granite, and during the waning stages of the regional ca. 1.4 Ga ferroan magmatism. Whole-rock geochemistry and Nd-Sr-Pb isotope compositions of samples from the Hardscrabble Creek complex are similar to those of the San Isabel Granite, suggesting that they were derived from the same or a similar mantle source. The mineral chemistry of the samples is comparable to Proterozoic massif-type anorthosites and related mafic intrusions, indicating that the Hardscrabble Creek complex and San Isabel Granite together represent a rare anorthosite-mangerite-charnockite-granite (AMCG) suite in the southwestern USA. The Hardscrabble Creek complex is unique because it formed ~80 m.y. after the other few mafic plutons in the southwestern USA, and it contains an ultramafic section that is absent from these plutons and rare to the AMCG suite in general. A combination of arc-like whole-rock geochemistry, chondrite uniform reservoir-like Nd-Sr-Pb isotopes, and ocean island basalt (OIB)-like zircon trace element chemistry suggests that the complex was derived from a partial melt of OIB-like mantle and interacted with metasomatically enriched lithospheric mantle. The enriched lithospheric mantle signature, combined with the long ~160 m.y. duration of magmatism in the region, is consistent with a period of protracted convergent tectonism.

Colorado↗

USGS2026

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10 6 m 3 of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H 2 O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Hawaii↗

USGS2023

Overview of the Cenozoic geology of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia

The Harrat Rahat volcanic field, located in the west-central part of the Kingdom of Saudi Arabia, is one of the larger Cenozoic harrats among the more than 17 harrats situated upon the Arabia Plate. The map plate contained herein shows, at a scale of 1:100,000, the mapped volcanic geology of northern Harrat Rahat, which consists of the northernmost one-fifth of Harrat Rahat. Northern Harrat Rahat has an area of about 3,340 square kilometers (km 2 ), of which 2,567 km 2 is covered by Harrat Rahat volcanic rocks and deposits, and it encompasses more than 900 exposed vents (that is, craters, cryptodomes, maars, and scoria cones), 289 of which are isolated by younger volcanic deposits and have not been correlated with the 234 volcanic rock units distinguished by geologic mapping. Northern Harrat Rahat is of special interest owing to its proximity to the city of Al Madinah al Munawwarah, which sits within, and is continuing to expand southward over, the north end of the volcanic field. Al Madinah is home to an expanding population, currently at more than 2 million residents, together with the intermittent addition of approximately 3 million pilgrims during Hajj and Umrah (religious visitations). The center of Al Madinah is less than 8 km from lava flows of the only confirmed historically documented eruption, which occurred in 1256 C.E. (654 A.H.). Earlier prehistoric lava flows also encroached into the area of the present-day city limits, as demonstrated by volcanic rocks exposed widely throughout the city in roadcuts, parks, and excavations for new buildings, although no evidence has been found of any other than the 1256 C.E. lava having reached that area during times of human habitation. Geologic mapping was undertaken by the U.S. Geological Survey in collaboration with the Saudi Geological Survey. The features of primary interest within the map area are scoria cones, lava flows, lava domes, craters, and pyroclastic deposits from the Quaternary, which have compositions of basalt, hawaiite, mugearite, benmoreite, and trachyte. The geologic mapping was published by Downs and others (2019) on a single sheet at 1:75,000 scale and two detailed sheets at 1:25,000 scale, accompanied by detailed explanations of the geology. This report presents the geology on a single sheet at 1:100,000 scale and provides condensed geologic explanations for the convenience of readers of this volume. Some minor errors of nomenclature and lava source regions that do not change fundamental interpretations have been corrected herein.

Professional Paper↗