USGS Science⌕ Search

SEARCH · USGS Science

Results for “Energy and Fuels”

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,656 records · Page 70Linked to original sources
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2025

A method to obtain remotely sensed grain size distributions from nonplanar granular deposits

Constraining the grain size distribution of granular deposits with complex surfaces is difficult with existing approaches. Field and laboratory techniques are time consuming and limited by the maximum grain size that laboratories can accommodate. In this study, we present a new method to identify the coarse fraction of the grain size distribution at a debris-flow fan deposit surveyed with terrestrial laser scanning (TLS) in Glenwood Canyon, Colorado, USA. This method is a novel grain segmentation algorithm developed for application to point cloud data of deposits with complex surfaces and angular grains ranging in size from centimeters to a meter. This approach combines an existing random forest machine learning method with a novel iterative clustering algorithm. We compared the grain size distribution from our algorithm with a Wolman pebble count conducted in the field, and found a root mean squared error of less than 2 cm from the 5th to 95th percentile of the grain size distribution of grains ranging from cobble to boulder sized (6.3–78 cm in our application). Finally, we compared our new algorithm with an existing open-source grain segregation algorithm, and our method outperformed the selected alternative when applied to the debris-flow deposit point cloud.

Colorado↗
USGS2024

Coastal Science Navigator companion guide—Discover the U.S. Geological Survey coastal science products you need

The Coastal Science Navigator is an online gateway to a wide variety of U.S. Geological Survey (USGS) coastal change hazards-related information, data, and tools relevant to stakeholders’ scientific and decision-making needs. The products within the Coastal Science Navigator provide data related to past, present, and future threats to our coastlines. The filter search allows users to see all available products and identify relevant options by time scale, geographic scope, coastal hazard theme, and other filters. The guided search suggests products based on users’ answers to a short series of questions. A comprehensive summary is available for each product. The idea for the Coastal Science Navigator arose in 2020 in response to stakeholder feedback identifying the need for a central source for USGS coastal science information. It was published in July 2023 and initially included 55 products. Regular updates are planned to integrate other existing and new products. This guide introduces some of the many coastal change hazards-related products available through the USGS. In it, we showcase the products included in the Coastal Science Navigator’s initial publication in July 2023. While it is not representative of all the information, tools, and data available, we hope it serves as a compelling snapshot of what the USGS has to offer and encourages you to explore the Coastal Science Navigator to discover more of the products you need. To navigate this guide, the products have been organized by the time scale they are best suited for—past, present, or future—although many products cover multiple time scales. An additional section features software, one of the many product types available as filters within the Coastal Science Navigator. Other products include downloadable data, websites, and geonarratives (web pages that combine text, images, and interactive maps into narratives you can scroll through). Featured geographic scopes are also highlighted within this guide, detailing some of the many regions in which the USGS conducts research and illustrating another way to filter products within the Coastal Science Navigator.

Circular↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗
USGS2025

Did they feel it? Legacy maroseismic data illuminates an engimatic 20th century earthquake

The challenges and the importance of preserving legacy instrumental records of earthquakes are now well-recognized (e.g., Richards & Hellweg, 2020, https://doi.org/10.1785/0220200053 ). Seismologists may not be aware of parallel challenges and opportunities with legacy macroseismic data for earthquakes in the United States. For much of the 20th century, macroseismic data were collected by a series of U.S. government agencies using a standard questionnaire distributed on postcards. Published summaries of postcards provide macroseismic data akin to modern Did You Feel It? questionnaire responses. In this paper we focus on the M 6.5 Fickle Hill, California earthquake, on 21 December 1954 (Hellweg et al., 2025) as a proof-of-concept, illustrating the potential of what we dub Did They Feel It? (DTFI) data to improve our understanding of significant 20th century U.S. earthquakes for which instrumental data are sparse. Legacy macroseismic data interpreted following modern conventions can potentially constrain traditional ShakeMaps at a level of detail and accuracy that in some respects rival maps for modern earthquakes. The updated ShakeMap for the 1954 Fickle Hill earthquake, also drawing from recently published media and first-person accounts, supports the location, depth, and stress drop value estimated from available instrumental data (Hellweg et al., 2025).

California↗
USGS2025

Widespread anhydrite saturation in Laramide-age arc magmas of southwestern USA

Anhydrite is considered a rare mineral phase in magmas, with only ∼33 documented occurrences worldwide. However, anhydrite readily decomposes in the near-surface environment, making it difficult to recognize its former presence in rocks collected at or near Earth’s surface. In such samples, only small anhydrite inclusions fully shielded within other minerals can have survived. During a recent field trip to the southwestern USA, we sampled 17 Laramide-age (ca. 40−80 Ma) magma systems, most of which are associated with porphyry copper deposits. A systematic search for anhydrite inclusions preserved within apatite, amphibole, plagioclase, and quartz phenocrysts in ∼100 rock samples by optical microscopy and Raman spectroscopy revealed that each of these 17 magma systems was at least temporarily anhydrite-saturated. Also, most previously identified magmatic anhydrite-bearing intrusions are associated with porphyry copper deposits, and both intrusive and volcanic rocks containing magmatic anhydrite show high Sr/Y ratios. These observations suggest that anhydrite saturation and porphyry copper formation are linked via magma fractionation at high pressure. Compared to average arc magmas, anhydrite-bearing magmas are unusually oxidized and sulfur-rich and seem to also be unusually water-rich. Hence, our preferred interpretation is that magma generation and/or fractionation at high pressure promotes the formation of superhydrous and oxidized magmas, which in turn promotes high sulfur contents and ultimately the precipitation of anhydrite. The high mineralization potential of these magmas does not need to result from their high sulfur content but could be caused by other properties of high-pressure magmas.

Arizona, New Mexico↗
USGS2026

Black rails in Baja California

Summarizing records from 1905 to 2025, we document the presence and distribution of the California black rail ( Laterallus jamaicensis coturniculus ) in northwestern México, particularly Baja California and the Colorado River Delta. This subspecies has a restricted distribution in the western U.S. and northwestern México. We found few records over the past century: sparse detections in the early 1900s and an increase in reported records after 2000. Over 120 years, we compiled 35 documented records of California black rail occurrence in Baja California with a minimum reported count of 56 rails. Our call-broadcast survey methods reconfirmed previously occupied sites and documented the species at a previously unreported site. Surveys in the Baja California portion of the Colorado River Delta yielded 20 detections across nine survey events between 2000 and 2022. In April 2003, standardized surveys detected 12 California black rails: 7 in Bahía de San Quintín, 3 in Arroyo San Telmo, and 2 in Arroyo El Rosario. Our recent survey work (2022–2025) recorded maximum monthly counts of 31 detections at Bahía de San Quintín and 15 detections at Estero de Punta Banda, the latter representing a previously undocumented location for the subspecies. These surveys also yielded the first photographic and acoustic documentation of California black rail in México. The contribution of spatial and temporal data provided by citizen science platforms increased our understanding of the subspecies’ distribution. Enforceable protection of this subspecies’ habitat is needed to safeguard its persistence.

Baja California↗
USGS2024

Environmental monitoring of groundwater, surface water, and soil at the Ammonium Perchlorate Rocket Motor Destruction Facility at the Letterkenny Army Depot, Chambersburg, Pennsylvania, 2021

Letterkenny Army Depot in Chambersburg, Pennsylvania, built an Ammonium Perchlorate Rocket Motor Destruction (ARMD) Facility in 2016 to centralize rocket motor destruction and contain all waste during the destruction process. The U.S. Geological Survey has collected environmental samples from groundwater, surface water, and soils at ARMD since 2016. During 2021, samples were collected from four groundwater wells in September, one surface-water site in October, and five soil sites in November near the facility. Samples were analyzed for nutrients, trace metals, major ions, total volatile organic compounds, and perchlorate. Perchlorate was not detected in any 2021 samples. Groundwater results showed no constituents exceeded any U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL). Dissolved arsenic (As) was detected in one well above the reporting detection level (RDL) of 3 micrograms per liter (μg/L) at 5.4 μg/L but below its MCL of 10 μg/L. Dissolved iron (Fe) was the only inorganic constituent measured above an EPA secondary maximum contaminant level (SMCL). All groundwater samples collected in 2021 exceeded the Fe SMCL of 300 μg/L, with concentrations ranging from 390 μg/L to 3,500 μg/L. Surface-water data collected during 2021 showed no measured constituents in the surface-water sample that exceeded any EPA MCL or SMCL. Soil samples collected from 2016 through 2021 showed all concentrations of As exceeded the EPA soil screening levels of 3 milligrams per kilogram (mg/kg) but did not exceed the Pennsylvania medium-specific concentrations for As of 61 mg/kg. Arsenic concentrations in 2021 ranged from 9.1 mg/kg to 12.9 mg/kg. The 2021 results for the ARMD Facility indicate no increases in concentrations of reported compounds compared to data from 2016 to 2020. The contained burn treatment facility for demilitarization of rocket motors during 2021 appears to have operated without elevating concentrations of target compounds compared to previous years.

Pennsylvania↗
USGS2024

U.S.-Mexico Borderland & vegetation community map

People on both sides of the United States-Mexico border need a high-resolution, binational vegetation community map that spans the entire United States-Mexico borderlands. Traditionally, mapping efforts in this region were impeded by complex logistics related to the international border, differing national needs and plans, and resource allocations and priorities. To address this need, scientists from the U.S. Geological Survey (USGS) Southwest Biological Science Center partnered with the Sonoran Joint Venture, the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program, data engineers from the Department of Biosystems Engineering at the University of Arizona, and collaborators from the Wildlands Network, the Borderlands Program to produce the first prototype land cover map within the overlapping Mojave Desert, Sonoran Desert, and the North American Bird Conservation Initiative’s Bird Conservation Region 33 (BCR33) using Landsat satellite data . BCR33 is an area of high biodiversity, providing habitat for bird species of concern and other wildlife. The land cover map supports FWS recovery plan efforts related to conservation planning activities for many species, including Yellow-billed Cuckoo ( Coccyzus americanus ), Cactus Ferruginous Pygmy-Owl ( Glaucidium brasilianum cactorum ), Southwestern Willow Flycatcher ( Empidonax traillii extimus ), Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis ), Bendire’s thrasher ( Toxostoma bendirei ), LeConte’s thrasher ( Toxostoma lecontei ), Masked Bobwhite ( Colinus virginianus ridgwayi ), jaguar ( Panthera onca ), and endangered plants such as Bartram’s stonecrop ( Graptopetalum bartramii ) and the Pima pineapple cactus ( Coryphantha robustispina ssp. robustispina ). In 2024, a Phase-II map for the full BCR33 region was completed, increasing the understanding of the binational nature of natural communities. The published map and associated paper can be found here .

Borderland↗
USGS2025

Irrigated agriculture influences selenium levels in an endangered marsh bird

Selenium bioaccumulation in aquatic food webs poses risks to wildlife, particularly in wetlands receiving irrigation runoff. The Salton Sea, California’s largest lake, is primarily sustained by agricultural drainage. This drainage creates wetland habitat along the lakeshore that many bird species depend on, including the federally endangered Yuma Ridgway’s rail ( Rallus obsoletus yumanensis ). However, these marshes may pose an ecological trap – attracting rails despite high selenium exposure. We captured rails during the 2020–2023 breeding seasons and compared rail selenium levels within three types of marshes (fed with irrigation runoff, Colorado River water, or groundwater). We collected blood, breast feathers, and head feathers of rails in all three water sources for selenium comparisons. We tagged adult rails with GPS transmitters to locate nests and foraging locations where we collected eggshells, unhatched eggs, and prey. We assessed selenium exposure by collecting multiple prey species commonly eaten by rails in all three water sources. Selenium concentrations varied among sampling locations. Selenium concentrations in most sample types were predominately influenced by water source and marsh inflow velocity (sometimes in combination with marsh size). Distance to inflow, however, did not influence selenium concentrations in any sample type. Selenium concentrations were highest in agricultural-fed marshes compared to river-fed and spring-fed marshes. Increased marsh inflow velocities resulted in lower selenium concentrations. Given the risk of an ecological trap, our results suggest that supplementing wetlands with Colorado River water could mitigate selenium bioaccumulation in Yuma Ridgway’s rails.

California↗
USGS2023

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2024

Timing and geometry of the Chemehuevi Formation reveal a late Pleistocene sediment pulse into the Lower Colorado River

The Chemehuevi Formation is a distinctive 50−150-m-thick wedge-shaped Pleistocene sedimentary unit deposited by the Colorado River. It lines the perimeters of the river’s floodplains and bedrock canyons for more than 600 km between the mouth of the Grand Canyon and the delta region in the Gulf of California. The formation is composed of a basal tan to light-yellowish-brown and pale-orange mud-dominated facies overlain and interbedded by a light-yellow-brown sand-dominated facies. The unit is one of two extensively exposed aggradational packages in the Lower Colorado River corridor, in addition to a series of other smaller alluvial terrace deposits. The Chemehuevi Formation appears to represent the response of a fully integrated Colorado River system to a significant perturbation, in contrast to the Bullhead Alluvium, which is likely a unique result of Pliocene river integration. The aggradation of the Chemehuevi Formation in the Lower Colorado River corridor may be similarly due to a unique event in the Colorado River system, or it may instead be a well-preserved sedimentary sequence recording typical behavior of the Colorado River below the Grand Canyon in the late Pleistocene. As such, multiple causal mechanisms have been proposed, but no study to date has conclusively explained the Chemehuevi Formation. To help resolve its timing, duration, and origin, we applied post-infrared infrared stimulated luminescence, carbonate U-Th series, and zircon sensitive high-resolution ion microprobe U-Th series geochronology to determine the ages of key exposures of the unit over a wide spatial area. These new data demonstrate that the Chemehuevi Formation was deposited ca. 110−90 ka. The depositional ages collectively overlap, suggesting that deposition occurred rapidly relative to the resolution of the geochronometers. The new depositional timing coincides with a shift from glacial to interglacial conditions after the marine isotope stage 5-6 transition. This observation is consistent with a climate-induced sediment pulse as a causal mechanism, yet correlations with similar deposits in the Colorado River headwaters or in neighboring catchments appear elusive. Potentially, climate transitions between glacial and interglacial periods induced a sediment pulse from hillslopes of the Colorado River system that resulted in the Chemehuevi Formation. An alternative or additional explanation is that the Chemehuevi Formation represents release of lava dam−impounded sediment in the Grand Canyon. The surface geometry of the Chemehuevi Formation projects upstream to the approximate location of lava dams, and the largest possible lava dam impoundment (the Upper Prospect dam) is comparable in volume to the formation. The lava dam hypothesis appears to be a possible explanation for the Chemehuevi Formation. However, tying deposition to a specific lava dam or series of lava dams remains challenging due to discrepancies in timing and volume. The combined effects of a series of lava dams may have led to the Chemehuevi Formation, as the last Pleistocene lava dam eruption coincides with the onset of deposition. Alternatively, the formation may result from the combined effects of both regional climate transitions and the lava dams that created a transient reservoir to compound a climate transition−driven sediment pulse. The geochronologic data presented here do not allow us to distinguish between the lava dam or climate transition hypotheses but will need to be reconciled with any future proposed depositional model.

Arizona, California, Nevada↗
USGS2024

Using stable oxygen isotope dual-inlet isotope-ratio mass spectrometry to elucidate uranium transport and mixed 230Th/U calcite formation ages at the seminal Devils Hole, Nevada, natural laboratory

Rationale Vein calcite in Devils Hole has been precipitating continuously in oxygen-isotope equilibrium at a constant temperature for over 500 000 years, providing an unmatched δ 18 O paleoclimate time series. A substantial issue is that coeval calcite (based on matching δ 18 O values) has uranium-series ages differing by 12 000 years. Methods An unparalleled high-accuracy δ 18 O chronology series from continuously submerged calcite was used to correct the published uranium-series ages of non-continuously formed calcite in two cores, cyclically exposed by water-table decline during glacial–interglacial transitions. This method relies on the premise that the δ 18 O values of coevally precipitated calcite are identical, allowing matching calcite δ 18 O values to establish formation ages. Results Exposed calcite can have apparent ages that are 12 000 years too young due to unrecognized uranium mobility and resulting mixed ages identified in over 50 mixed uranium-series ages from previous studies. Secondary uranium in fluids, sourced from the formation or dissolution of porous carbonate deposits (folia) with high uranium-238 ( 238 U) concentrations, has migrated up to 10 mm into vein calcite. Conclusions The continuously submerged Devils Hole δ 18 O chronology is not explained by orbital forcing. Rather, this chronology represents a regional climate record in the southern Great Basin of sea-surface-temperature (SST) variations off California, variations that preceded the last and penultimate deglaciations by 5000 to approximately 10 000 years. Temporal discrepancies between the continuously submerged Devils Hole chronology and other regional δ 18 O records (e.g., the Leviathan chronology) can be explained by unrecognized cryptic, pernicious uranium mobility, leading to model estimations that may be thousands of years younger than actual ages. Consequently, paleo-moisture availability, water-table, and groundwater recharge models based on these mixed uranium-series ages are too young by as much as 12 000 years. The potential for post-formation uranium addition in subaerial cores and speleothems underscores the need for caution in uranium-series dating, highlighting δ 18 O time-series comparisons as a method for identifying mixed ages.

Nevada↗
USGS2026

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗
USGS2025

lasertram: A Python library for time resolved analysis of laser ablation inductively coupled plasma mass spectrometry data

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) data has a wide variety of uses in the geosciences for in-situ chemical analysis of complex natural materials. Improvements to instrument capabilities and operating software have drastically reduced the time required to generate large volumes of data relative to previous methodologies. Raw data from LA-ICP-MS, however, is in counts per unit time (typically counts per second), not elemental concentrations and converting these count ratesto concentrations requires additional processing. For complex materials where the ablated volume may contain a range of material compositions, a moderate amount of user input is also required if appropriate concentrations are to be accurately calculated. In geologic materials such as glasses and minerals that potentially have numerous heterogeneities (e.g., microlites or other inclusions) within them, this is typically determiningwhether the total ablation signal should be filtered to remove these heterogeneities. This necessitates that the LA-ICP-MS data processing pipeline is one that is not automated, but is also designed to enable rapid and efficient processing of large volumes of data. Here we introduce , a Python library for the time resolved analysis of LA-ICP-MS data. We outline its mathematical theory, code structure, and provide an example of how it can be used to provide the time resolved analysis necessitated by LA-ICP-MS data of complex geologic materials. Throughout the pipeline we show how metadata and data are incrementally added to the objects created such that virtually any aspect of an experiment may be interrogated and its quality assessed. We also show, that when combined with other Python libraries for building graphical user interfaces, it can be utilized outside of a pure scripting environment. can be found at https://doi.org/10.5066/P1DZUR3Z

Applied Computing and Geosciences↗