USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Enviornmental Science”

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,605 records · Page 78Linked to original sources
USGS2024

Deep lithospheric controls on surface deformation and seismicity around the East Anatolian Fault Zone and A3 Triple Junction

The East Anatolian Fault Zone (EAFZ) is a plate-bounding strike-slip fault capable of hosting large earthquakes, as demonstrated by the extremely damaging February 2023 M w 7.8 and M w 7.7 mainshocks of the Kahramanmaraş earthquake sequence. Deformation related to this boundary, part of the Anatolia-Arabia-Africa (A 3 ) Triple Junction, is diffuse, as was shown by part of this earthquake sequence occurring on a northern splay of the EAFZ (the Sürgü-Çardak Fault Zone; SCFZ). Controls on surface deformation are commonly linked to stress in the brittle upper crust, but the complex deformation and seismicity patterns in this region may also reflect deeper processes, such as variations in the location and extent of the strong Arabian Plate lithospheric mantle. Seismic tomography indicates that the Arabian Plate underthrusts Anatolia as far north as the SCFZ and extends as far west as the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced M w >6 earthquakes. By investigating the relationship between deformation since the inception of the EAFZ (ca. 5 Ma), seismic structure, and seismicity, we infer that the SCFZ will become the future SE boundary of the Anatolian Plate as part of the evolving A 3 Triple Junction.

GSA Today↗
USGS2025

Groundwater salinity: Applying the specific conductance and water type proxy

Groundwater is increasingly needed for water supplies but may have limited utility in some locations because of its salinity. Salinity, often expressed as total dissolved solid (TDS), is frequently estimated using specific conductance (SC) measurements. However, the commonly used proxy (0.65 multiplied by SC to indicate TDS, common in many handheld meters) can result in inaccurate TDS estimates. First, the TDS–SC relationship is not linear over the entire concentration range of groundwater. Furthermore, the TDS (and salinity)–SC relationships vary substantially depending on the major-ion composition. Here we develop a proxy method utilizing SC and major-ion water type to estimate TDS and salinity specifically for groundwaters. Compared to most surface waters, groundwater tends to have a wider range of salinity (fresh to highly saline) and higher concentrations of bedrock-derived solutes such as carbonate ions, silica, and many other ions. The dataset used to develop the proxies includes water chemistry data from 149,059 discrete groundwater samples. The groundwater proxies, which employ nonlinear log–log relations, utilize five water types (HCO 3 , Cl, Ca-Mg-SO 4 , Na-K-SO 4 , and mixed waters), are accurate (median percent difference between TDS and salinity determined using the proxy compared to discrete measurements was <±0.8%) over a wide range of SC (up to 200 mS/cm), rapid, cost-effective, and can be measured on-site.

Groundwater↗
USGS2024

Exploring the dynamic interactions between the Southern San Andreas Fault and a normal fault under the Salton Sea

We investigate the dynamic interactions between the Southern San Andreas Fault (SSAF) and a proximal normal fault (NF) beneath the Salton Sea in southern California. The NF, positioned near the SSAF terminus at Bombay Beach, exhibits 11–15 displacement events across 14 stratigraphic sequences, with a range of 0.2–1.4 m of vertical offset since ∼2–3 ka. Notably, four of these events may align temporally with SSAF earthquakes, raising questions about the possible interplay between the two faults. Utilizing dynamic rupture models, we analyze the coseismic interactions between the SSAF and NF, addressing under what conditions the SSAF induces slip on the NF. Our findings reveal that a suite of SSAF ruptures, particularly those propagating from north to south, can trigger slip on the normal fault and replicate observed vertical offsets. If the SSAF extends beneath the Salton Sea, earthquakes originating south of the NF intersection are less likely to trigger normal fault slip, although we cannot exclude this possibility. Some SSAF ruptures do not trigger discernible slip on the NF, rendering such events undetectable in the stratigraphic record. Our research contributes toward discussions regarding the seismic hazard in southern California, shedding light on the interplay between the SSAF and NF.

California↗
USGS2024

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone↗
USGS2025

Grammar to graph—An approach for semantic transformation of annotations to triples

Data annotation is the process of labeling data to show the outcome that a related data model should predict. In this study, annotation data were transformed into semantic graph triples, mainly for use with the Resource Description Framework (RDF), a type of entity-relationship-attribute data model for graph databases. The transformation of annotation data to semantic graph triples provides complex linguistic meaning with data handling advantages such as reduced data storage needs, improved logical specification of relations between objects, and reusable classes and properties that support logic and inference. A grammar-based framework in graph form supports user questions and queries. The words defining approximately 334 topographic feature types compiled by the U.S. Geological Survey were tokenized as units of analysis and grouped by part of speech. Their dependency relations were identified for this study using natural language processing libraries. Dependency concepts are used as structured semantic relations among part-of-speech classes. Tokens, units equivalent to words, form instances of classes and were quantified within a tabular output format using PostgreSQL data storage software. Table data were logically aligned as triples following a mapping file and stored with an ontology file using Ontop virtual triplestore software. A grammar ontology schema for the data was synchronized to match queries whose results validated the graph’s structure. The text analysis produced 8 part-of-speech classes of content words for object representations and 4 classes of function words for operational applications. Dependency relations formed 27 ontology properties for topographic subgraph structures. Token occurrences shaped overall ontology salience and formed a lexicon of syntactic terms for subgraph objects and properties. The schema ontology of class and property population shapes formed the lexicon of English terms. SPARQL Protocol and RDF Query Language (SPARQL) was used with the lexicon to conform data to RDF guidelines. This study confirms the hypothesis that although linguistic logic varies from description logic, its approximation applies to ontology design. Property and query use case patterns extracted from the analysis support queries concerning complex topographic relations and patterns normally embedded within text definitions. The method used in this study could be applied to text forms in other domains, such as survey notes.

Scientific Investigations Report↗
USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗
USGS2024

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

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

Nevada↗
USGS2026

The Sedimentary Geochemistry and Paleoenvironments Project Phase 2 data release: An open data resource for the study of Earth's environmental history

Geochemical data from sedimentary rocks are the primary source of information regarding Earth's surface evolution through time, including its air and water envelopes and interactions with life and deep Earth processes. The Sedimentary Geochemistry and Paleoenvironments Project (SGP) is a scientific consortium centered around open data and community-driven development of cyberinfrastructure tools and resources for sedimentary geochemistry and Earth history. Here we describe the SGP Phase 2 data release, which focused on incorporating Paleoproterozoic and Mesoproterozoic (2500–1000 million years ago) data and better accommodating carbonate data. This data release was built through the involvement of >200 researchers worldwide in academia, government, and industry, and provides the largest available public data resource for our user community in the academic fields of geochemistry, sedimentology, tectonics, paleontology, Earth history, and paleoclimate, as well as the petroleum and minerals industries. The dataset now encompasses 126,006 samples and 4,132,371 geochemical analyses. In addition to direct entry by SGP Team Members, we have ingested and incorporated datasets from the Geoscience Australia OZCHEM database, the Alberta Geological Survey, and the Deep-Time Marine Sedimentary Element Database (DM-SED) compilation. This paper details sampling in the Phase 2 dataset with respect to age, geography, lithology, and other geological characteristics, documents access via our search website and API, discusses possible issues and/or biases in the dataset that could impact analyses, describes plans for governance and stewardship of data from Indigenous lands, and serves as the citable reference paper for the data release.

Chemical Geology↗
USGS2025

Framework for assessing intranational risk of critical mineral sectors: Case study of India

Recent geopolitical actions regarding the supply of critical minerals have highlighted both the threats posed by countries dominating a critical mineral supply chain, and the potential mitigations served by alternative supplier countries. However, no systematic quantification exists for the intranational factors that determine such threats or mitigations. This study outlines a methodological framework for quantifying the intranational risk profile of critical mineral sectors and compares them with the sector's international potential for leverage. The four intranational factors of concentration, synchronization, private sector involvement, and regulatory obstacles were analyzed across multiple scales in order to determine an intranational risk score that quantifies the resilience or fragility of the critical mineral sector in question. Network analysis was also separately conducted to identify the most important entities operating across all critical mineral sectors. Finally, this study assessed the critical mineral sector's global share of reserves and production to determine an international risk score that quantifies the country's critical mineral sector's leverage potential. A case study of India showed that two critical mineral sectors with high leverage potential, chromite and barite, were respectively intranationally resilient and fragile. Of the remaining sectors with minimal leverage, five sectors such as aluminum were intranationally resilient and may serve as viable mitigation options for future critical mineral supply. Three sectors such as titanium were intranationally fragile, but recent regulatory developments may both decrease intranational risk and increase international risk, potentially making them major critical mineral sectors in the future.

Resources Policy↗
USGS2025

U.S. Geological Survey Integrated Water Availability Assessment—2010–20

This professional paper is a multichapter report that assesses water availability in the United States for water years 2010–20. This work was conducted as part of the fulfillment of the mandates of Subtitle F of the Omnibus Public Land Management Act of 2009 (Public Law 111-11), also known as the SECURE Water Act. As such, this work examines the spatial and temporal distribution of water quantity and quality in surface water and groundwater, as related to human and ecosystem needs and as affected by human and natural influences. Chapter A introduces the National Integrated Water Availability Assessment and provides important background and definitions for how the report characterizes water availability and its components. Chapter A also presents the key findings of Chapters B–F and thus acts as a summary of the entire report. Chapter B is a national assessment of water supply, which is the quantity of water supplied through climatic inputs. Chapter C is a national assessment of water quality, which is the chemical and physical characteristics of water. Chapter D assesses water use including withdrawals and consumptive use in the conterminous United States. Chapter E presents an analysis of factors affecting future water availability under changing climate conditions. The National Integrated Water Availability Assessment culminates with Chapter F, which is an integrated assessment of water availability that considers the amount and quality of water coupled with the suitability of that water for specific uses. Together, these six chapters constitute the National Integrated Water Availability Assessment for water years 2010–20.

Professional Paper↗
USGS2023

An inventory of three-dimensional geologic models—U.S. Geological Survey, 2004–22

A database of spatial footprints and characteristics of three-dimensional geological models that were constructed by the U.S. Geological Survey between 2004 and 2022 was compiled as part of ongoing development of subsurface geologic information by the USGS National Cooperative Geologic Mapping Program. This initial inventory resulted in the compilation of 38 three-dimensional geological models that vary widely in their spatial extent, the type and purpose of the model, the number of subsurface units characterized by the model, and the software platforms used to create the model. This Data Report provides the scientific rationale and explanation of the contents of a companion USGS digital data release of spatial data and attributes associated with each three-dimensional model.

Data Report↗
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↗
USGS2024

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

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

Reliability Engineering and System Safety (RESS)↗
USGS2024

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota↗
USGS2026

Streamflow, water quality, and constituent loads and yields, Scituate Reservoir drainage area, Rhode Island, water year 2023

As part of a long-term cooperative program to monitor water quality within the Scituate Reservoir drainage area, the U.S. Geological Survey, in cooperation with Providence Water (formerly the Providence Water Supply Board), collected streamflow and water-quality data in tributaries to the Scituate Reservoir, Rhode Island. Streamflow and concentrations of chloride and sodium estimated from records of specific conductance for 16 tributaries were used to calculate loads of chloride and sodium during water year 2023 (October 1, 2022, through September 30, 2023). Water-quality samples were collected by Providence Water at 38 sampling stations on tributaries to the Scituate Reservoir during water year 2023. These water-quality data are summarized by using values of central tendency and are used, in combination with measured (or estimated) streamflows, to calculate loads and yields of selected water-quality constituents for water year 2023. Annual mean streamflows for monitoring stations in this study ranged from 0.47 to 34.4 cubic feet per second during water year 2023. At the 16 continuous-record streamgages, tributaries transported about 3,000 metric tons of chloride and 1,900 metric tons of sodium to the Scituate Reservoir; annual chloride yields for the tributaries ranged from 19 to 200 metric tons per square mile, and annual sodium yields ranged from 13 to 120 metric tons per square mile. At the stations where water-quality samples were collected by Providence Water, the medians of the median daily loads were 65,000 million colony forming units per day for coliform bacteria, 870 million colony forming units per day for Escherichia coli , 200 kilograms per day for chloride, 8.2 grams per day as nitrogen for nitrite, 520 grams per day as nitrogen for nitrate, and 300 grams per day as orthophosphate for phosphate. The medians of the median yields were 25,000 million colony forming units per day per square mile for coliform bacteria, 600 million colony forming units per day per square mile for Escherichia coli , 65 kilograms per day per square mile for chloride, 4.6 grams per day per square mile as nitrogen for nitrite, 180 grams per day per square mile as nitrogen for nitrate, and 180 grams per day per square mile as orthophosphate for phosphate.

Rhode Island↗
USGS2025

Re-Os geochronology of molybdenite by LA-ICP-MS/MS

Re-Os geochronology is a powerful tool to directly date sulphide mineralisation using phases such as molybdenite, pyrite and chalcopyrite. We present a snapshot of a method in development utilizing laser ablation inductively coupled-plasma ‘triple quad’ mass spectrometry (LA-ICP-MS/MS) for Re-Os geochronology of molybdenite. While traditional LA-(MC)-ICP-MS measurements are limited by the >90% isobaric interference correction of parent radioactive 187 Re on radiogenic daughter 187 Os, the MS/MS technology minimizes this correction through careful instrument tuning and use of a (CH 4 + H 2 + He) reaction gas mixture to mass shift 187 Os away from 187 Re. We demonstrate the potential of the MS/MS technique by presenting robust Re-Os dates for three different molybdenite samples ranging in age from ~28 to 1520 Ma.

Conference Paper↗
USGS2025

The U.S. Geological Survey National Water Quality Network—Groundwater, 2024

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, 8 of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,089 wells in the network. All wells are sampled for physical properties, nutrients, major ions, trace elements, per- and polyfluoroalkyl substances (PFAS), and a subset of wells are sampled for pesticides, volatile organic compounds, radionuclides, and microbiological contaminants.

Contiguous United States↗
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↗