Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Open Source Education”

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,244 recordsLinked to original sources
USGS2024

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education↗
USGS2025

ShakeAlert®—Communication, education, outreach and technical engagement strategic vision

Executive Summary In 2006, the U.S. Geological Survey (USGS) began directly supporting ShakeAlert ® research and in 2012 the ShakeAlert demonstration system began testing ( Given and others, 2018 ). The ShakeAlert earthquake early warning (EEW) system is a partnership between the U.S. Geological Survey (USGS) and the three West Coast States (Washington, Oregon, and California) served by the ShakeAlert System, which is part of the larger Advanced National Seismic System (ANSS). With more than 143 million people exposed to potentially damaging shaking in the United States ( Jaiswal and others, 2015 ), earthquakes are a national hazard. Most of our Nation’s earthquake risk is concentrated in the highly populated areas on the active plate tectonic boundaries on the West Coast of the conterminous United States. ShakeAlert is the first public alert system in the United States to provide rapid mass notification of earthquake detection, potentially offering seconds of warning before strong shaking arrives. A few seconds may not seem like much time, but the information in ShakeAlert Messages can be used to trigger automated actions that can prevent injury or death, reduce immediate damage, and speed recovery from earthquakes. The information product issued by the ShakeAlert system is called a ShakeAlert Message and is one of the information products and tools of the ANSS. The ShakeAlert System includes the USGS component, plus the pathways by which ShakeAlert-powered products and (or) services are delivered to end users. Alerts can be delivered to cell phones or be used to trigger automated systems to protect equipment, facilities, and infrastructure, such as slowing or stopping a train. ShakeAlert-powered automated actions can include fire house doors that can be opened to prevent jamming, heavy equipment (for example, trains, elevators, and cranes) that can be automatically stopped or parked in safe positions, and pipeline valves that can be closed to prevent surges and spills. A few seconds of warning also may be sufficient for people to take protective actions, such as drop, cover, and hold on or modified protective actions for a broad range of populations. Advance training may increase the benefit of a speedy response to an alert. Outreach and education about EEW may raise awareness of the overall earthquake threat and how people can best react when they receive an alert or feel shaking. ShakeAlert communication, education, outreach, and technical engagement (CEO&TE) efforts are highly collaborative and essential for the success of the ShakeAlert System. This strategic vision informs how the vast ShakeAlert CEO&TE Community operates and works together. The CEO&TE Community delineates a strategic framework that is intended to set the path for a long-term, sustainable approach to CEO&TE through three focus areas and five priorities. Focus Areas Technical engagement.— The objective of this focus area is to expand and broaden ShakeAlert technical engagement and implementation. External engagement.— This focus area targets engagement in long-term communication, education, and outreach planning, implementation, and evaluation. Internal engagement.— The purpose of this focus area is to build and sustain a robust infrastructure to optimize collaboration, information sharing, and project planning among the community of stakeholders charged with implementing ShakeAlert CEO&TE. Priorities The enumeration of the five priorities listed below does not suggest priority ranking. Public safety, preparedness, and resilience; Technical implementation and engagement; Consistent messaging and communication; Integration with other Federal and State earthquake hazards products; and Educational resources development and dissemination. This strategic vision is a tangible outcome of collaboration among many stakeholders beginning in July 2016. Since then, the work of the ShakeAlert CEO&TE Community has grown into an international effort. The USGS has developed, tested, and implemented a broad spectrum of communication, education, and outreach tools and resources—all of which recognize that seconds matter when it comes to safety and mitigating harm from earthquake hazards. The CEO&TE social science research effort has provided invaluable insights into the ShakeAlert System’s human interface. USGS-licensed technical partners develop, test, and implement real-world applications using ShakeAlert Messages. The success of ShakeAlert CEO&TE efforts is predicated on robust collaboration across numerous agencies, organizations, and groups. As such, this strategic vision outlines a “partnership model” that delineates roles and responsibilities to ensure alignment with focus areas and priorities. The partnership model includes the CEO&TE lead agency (USGS); its principal partners (State agencies and university partners); its implementation partners (for example, technical partners who build systems to deliver ShakeAlert-powered products and (or) services [focus area one]), earthquake education partners who work to increase public preparedness for seismic events (focus area two); and other organizations that work together to enhance the adoption and effectiveness of the ShakeAlert System. These partners collaborate and convene through a variety of working groups and forums, which are also described in this strategic vision and align with focus area three (internal engagement). The CEO&TE Community collaboratively developed its operating principles and a consensus-based, decision-making strategic framework to guide its collective work. Performance metrics are used to continually measure success. Ultimately, the USGS and ShakeAlert CEO&TE Community are advancing the ShakeAlert System that as of the publication of this strategic vision to “provide earthquake early warning for all” serves more than 50 million people.

Circular↗
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

Computational electromagnetic geophysics for groundwater system studies: A review on established practices and recent advances

Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies. As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.

Journal of Hydrology↗
USGS2026

Critical minerals memory match game

An educational information packet about the 2025 List of Critical Minerals, which includes a memory match game about select critical minerals and how they are used.

General Information Product↗
USGS2025

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

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

Geoheritage↗
USGS2025

Preliminary field report of landslide hazards following Hurricane Helene

Executive Summary This report reflects our knowledge regarding the widespread landslide activity associated with Hurricane Helene observed during the U.S. Geological Survey’s (USGS) mission assignment to North Carolina in October 2024. The material in this report was originally prepared for the Federal Emergency Management Agency under mission assignment DR-4827-NC. The data and commentary in this report are reflective of a report provided to the Federal Emergency Management Agency (FEMA) on October 18, 2024, as well as information provided in briefings at the Buncombe County Emergency Operations Center. The report has been modified for public dissemination. This assessment was based on systematic visual examination and mapping of landslide locations from aerial and satellite imagery, visual and photographic observations from low-level helicopter overflights and conversations with local landslide experts from the North Carolina Geological Survey and Appalachian Landslide Consultants PLLC, and more than 50 years of combined landslide hazard professional experience of the mission-assigned field team. No systematic field investigations were done by the USGS. While responding to the event, the USGS did not identify any landslides that posed an immediate major threat to recovery personnel in parts of nine counties in North Carolina (Avery, Buncombe, Henderson, McDowell, Mitchell, Polk, Rutherford, Watauga, and Yancey); however, threats from renewed landslide activity may remain heightened in localized areas for months or even years. Known areas of the most abundant landslide occurrence include Bat Cave, Lake Lure, Chimney Rock, Swannanoa, Black Mountain, Fairview, steep areas in Asheville, and the Blue Ridge Parkway. The USGS shared detailed locations of known landslides with the Emergency Operations Centers. The thousands of landslide scars on hillsides and landslide deposits on flatter ground may present some threat to recovery activities. Soil and rocks will continue to erode from newly exposed landslide scars and may pose a threat to people and infrastructure who are immediately nearby. In general, the steeper and taller the landslide scar, the greater the potential threat. This threat is heightened during periods of rainfall and increases with the duration and intensity of rainstorms. Very heavy rainfall, or repeated rainfall events during short periods, could also initiate new landslides on steep slopes. Excavation of landslide deposits, particularly excavation of those deposits directly adjacent to steep slopes, may also pose a threat to nearby people and equipment. An interagency collaborative mapping effort led by the USGS that informed this assessment identified 1,155 landslide locations by the October 2024 briefings, but that number increased to 2,217 in a final reviewed version of the locations published in January 2025. Locations were mapped from satellite imagery, fixed-wing and helicopter surveys, media and social media, and field reports in the 3 weeks following the passage of the remnants of Hurricane Helene. USGS products outlined in this report are publicly available and include geotagged photographs from aerial reconnaissance, hazard models, an interactive view of mapped landslide locations, and landslide safety and education resources.

North Carolina, South Carolina, Tennessee, Virgini↗
USGS2025

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska↗
USGS2026

Open-source gravity reduction workflows for geothermal resource assessment

Potential-field geophysical data such as gravity can enhance understanding of geothermal resources at all stages of the resource life cycle, including assessment, exploration, development, and monitoring, and at multiple scales, from the reservoir scale to regional scale. However, to make gravity data useful for geothermal resource characterization, several processing steps are required to isolate the effects of density variations in the Earth’s crust to enable the identification of structural features associated with geothermal resources. Although this process is well-established, standard computational implementations for processing gravity data that are FAIR (Findable, Accessible, Interoperable, and Reproduceable) are still lacking. This paper details ongoing efforts at the U.S. Geological Survey (USGS) to develop a standard set of open-source Python tools for gravity data reduction that align with the FAIR principles. This workflow makes use of existing open-source tools for geophysical data processing with the goal of maximizing opportunities for rapid improvements, interoperability, and adaptability to other types of geophysical data.

Conference Paper↗
USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗
USGS2026

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

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

Seismological Research Letters↗
USGS2025

Applying U-Th disequilbrium for dating siliceous sinters

Continental hydrothermal systems are critical avenues for the crustal transport of heat and mass captured for geothermal energy and mineral exploration. Thus, understanding their temporal evolution and longevity is important for resource characterization. Deposits of microlaminated siliceous sinter, common surface expressions of high temperature reservoirs (> 170 °C), have the potential to trace hydrothermal histories. Geothermal reservoirs are often located within uranium-bearing silicic volcanic rock where subsurface fluid-rock interactions extract U into hydrothermal fluids. U incorporated in the surface sinter deposit has the potential to provide a U— Th disequilibrium dating option. We focus on samples from El Tatio geyser field in the Altiplano of northern Chile, the largest geothermal system in the Andes. Our resulting 230 Th/U ages, along with the water and deposit elemental compositions, suggest concentrations of U and Th vary predictably along the sinter apron. While distal facies containing the highest U concentrations (> 50 μg/g) are least affected by detrital Th corrections, they can display suspected open-system behavior. In contrast, more medial facies, where bacterial mats and other porous textures are co mmonly concentrated, have only trace amounts of U (< 0.1 μg/g), which leads to unreliable or geologically improbable dates. Proximal facies tend to date most consistently. By comparing existing 14 C ages with 230 Th/U results, 230 Th/U ages tend to be younger than the 14 C ages, supporting the presence of a 14 C-dead carbon influence. New data confirm that the onset of geothermal activity at El Tatio goes back to the late Pleistocene.

Altiplano, El Tatio geyser field↗
USGS2025

Statewide cumulative human health risk assessment of inorganics-contaminated groundwater wells, Montana, USA

Across the United States, rural residents rely on unregulated and generally unmonitored private wells for drinking water, which may pose serious health risks due to unrecognized contaminants. We assessed the nature, degree, and spatial distribution of cumulative health risks from inorganic contaminants in groundwater. Our analysis included nearly 84,000 data points from 6500+ wells, across 51 of Montana's 98 watersheds, using a public groundwater database. We compared a drinking water screening level cumulative risk assessment (CRA) for inorganics based on the U.S. Environmental Protection Agency (EPA) protective health thresholds (Maximum Contaminant Level Goals, Health Advisories [MCLG-HAs]) to a CRA based on EPA public supply enforceable Maximum Contaminant Levels (MCLs). Based on median concentrations of 19 inorganics (antimony, arsenic, barium, beryllium, boron, cadmium, chromium, copper, fluoride, manganese, molybdenum, nickel, nitrate, lead, selenium, strontium, thallium, uranium, zinc), 75% of watersheds had MCLG-HA-based cumulative risk values > 1.0; arsenic and uranium contributed the most risk, followed by strontium, fluoride, manganese and boron. Hence, this screening level (Tier I) CRA indicated widespread potential for unrecognized human health risk to private well users from inorganic contaminants considering both carcinogenic and non-carcinogenic risks. Sensitivity analysis showed that benchmarks applied (MCLG-HAs versus MCLs) exerted the largest control on results. Our findings identify priority regions for Tier 2 risk assessments to elucidate local sources and distributions of geogenic versus anthropomorphic contaminants. Our study is the first statewide assessment of cumulative health risk from groundwater that we are aware of, and results support increased statewide drinking water education and testing to reduce human health risks from contaminated private well water.

Montana↗
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

Distinguishing natural from mining-related metal sources by including streambank groundwater data in a stream mass loading study

Distinguishing stream metal loading caused by mine features from that caused by natural background sources remains challenging, yet this distinction is essential for making effective remedial decisions at many legacy mine sites. We combine a stream tracer injection and synoptic sampling study with data from shallow near-stream groundwater wells to estimate left-bank versus right-bank metal loading contributions at the 100-m spatial scale. The study was performed in the East Mancos River, a mountain headwater stream in Colorado, USA. The dominant source of elevated stream metal concentrations could be either groundwater infiltration through right-bank Doyle Mine waste piles or natural acid rock drainage from hydrothermally altered bedrock located mainly on the left bank. For the five metals of concern (Cu, Al, Zn, Cd, and Mn), we find that 15 % of the load contributed by diffuse groundwater inputs in the section potentially influenced by Doyle mine waste originates from the right bank. This right-bank potential mine contribution equates to only 3 % of the total watershed instream load for these metals. Furthermore, apparent 3 H/ 3 He groundwater ages in segments contributing most of the right-bank metal loading are sufficiently old (9–12 yr) to suggest that infiltration through the waste piles, located only 140–180 m from the stream, is unlikely. Estimated potential Doyle mine loading contributions can therefore be considered maximum values. Study results thus indicate that Doyle mine waste piles are a minor source of metal loading under low-flow conditions, and streambank groundwater data can provide valuable additional information in stream mass loading studies.

Colorado↗
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↗
USGS2026

Simulated ground motion dataset in the Azores Plateau, Portugal, on rock and soil sites

Building on a previously developed bedrock dataset, this study extends the Azores Plateau ground motion simulations to include soil-amplified records and introduces a comprehensive validation framework. Soil amplification is modeled using one-dimensional soil profiles. A stochastic source-based approach is employed to generate the dataset, incorporating randomization of input-model parameters to account for the aleatory uncertainty in seismic activity. The accuracy of the dataset is verified through a comprehensive validation framework, showing that the randomization effectively captures variance and inter-period correlation observed in records. This work provides a robust dataset for advancing seismic hazard and risk assessment in the Azores Plateau.

central and eastern Azores islands↗
USGS2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah↗