SEARCH · USGS Science

Results for “Bulletin, Southern California Academy of Sciences”

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

SourceResearch connectionFrom yearThrough year

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

1,666 records · Page 77Linked to original sources

USGS2026

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗

USGS2026

A practical decision tool for marine bird mortality assessments

Given the rise in anthropogenic, environmental, and disease events contributing to marine bird mortality, there is a critical need to improve the rigor of mortality assessments. Deficits in data collection and mortality estimation can hinder a manager’s ability to document the scale of events and assess population level impacts. Therefore, to inform decisions required during activities, such as conservation status assessments or harvest management, organizations may choose to incorporate mortality assessments into response plans. Resources, capacity, and assets to assess mortality vary across jurisdictions (federal, state, Indigenous, local, etc.), and clear guidance to support mortality assessments is often unavailable or not clearly addressed. Here, we present a decision support tool to help managers identify and evaluate survey options to assess bird mortality in a diverse array of scenarios. The objective of the decision tool is to improve data collection and availability, which will increase the ability to estimate mortality robustly, given situation-specific attributes and constraints. This decision tool is designed to guide the response when a mortality event is initially encountered and offers suggestions for assessment and reporting procedures in the absence of other guidance or to complement existing protocols. The decision tool is also meant to inform decision making for response determination and resource allocation. The tool facilitates examination of options for further assessment and monitoring, which users determine by examining questions pertaining to species prioritization, determination of mortality minimum spatial extent, and the potential magnitude of impacts on affected species. Finally, identification of appropriate survey methods that address imperfect detection when a complete census is not possible are determined by exploring location, spatial and temporal extent, and the type of species affected. Ultimately, this decision tool aims to facilitate and improve the standardization of mortality assessments, equipping managers with a practical resource to navigate the decision-making process for marine bird mortality estimation.

Ornithological Applications↗

USGS2025

Chemical and isotopic fractionation during melt inclusion formation

Melt inclusions are used to study the origin and evolution of magmas. The extent to which they represent equilibrium melt compositions, however, critically hinges on the ratio of crystal growth rate to diffusion rate in melts. If the rate of crystal growth is limited by the supply of nutrients and the dissipation of unwanted components, the trapped melt will be depleted in compatible elements and enriched in incompatible elements. Despite widespread recognition of the potential for melt inclusions to trap disequilibrium compositions, their identification remains a challenge. Here, we demonstrate how stable isotopes provide a solution to this problem. Melt inclusions in basaltic tephra from Kīlauea volcano, Hawai‘i, display Mg isotope fractionation up to 0.49‰ ( 26 Mg/ 24 Mg) relative to bulk rock. These observations indicate that kinetic processes, such as chemical diffusion, are at play, as equilibrium mineral-melt isotope fractionation at high temperatures is an order of magnitude smaller. We propose that the heavy Mg isotopic compositions of the melt inclusions were generated by the faster diffusion and thus preferential incorporation of 24 Mg by the growing olivine, leaving the slower diffusant ( 26 Mg) behind. The incompatible elements (e.g., Ca, Ti, and K) are predicted to display similarly large fractionations when boundary layer effects are significant. These findings show that stable isotopes can be a useful tool in identifying disequilibrium compositions in melt inclusions.

Geology↗

USGS2026

Coastal barrier resilience and resistance: Analysis and metrics for characterizing coastal state

Barrier islands are shaped by a variety of short- and long-term environmental processes such as storms and relative sea-level rise. These islands, found along the estuarine-marine interface, provide ecosystem services including storm surge and wave attenuation, erosion protection to inland marshes, habitat for fish and wildlife, and recreation. Natural resource managers require actionable information on how barrier island resilience and resistance changes over time to understand how an island’s current state relates to past conditions and to inform restoration prioritization and implementation. The U.S. Geological Survey and The Water Institute collaborated on a study to develop indicators of resilience and resistance for barrier islands in Louisiana. Here, resilience captures island persistence on yearly to decadal time scales, and resistance captures persistence on event time scales of days to weeks. The indicators fall in two categories: Tier 1 Screening Metrics, that can be readily calculated from available data, are easily interpretable as an evaluation of barrier condition, and provide a high-level snapshot of overall barrier resilience and resistance; and Tier 2 Analysis Metrics, which are detailed metrics that required specialized analysis and interpretation and are more applicable to answering specific questions managers may have about barrier state. The research team derived Tier 1 resilience indicators from subaerial land and vegetation cover calculated from publicly available maps and products based on satellite imagery. By benchmarking the total land and vegetation extent against their respective historical maxima, this metric provides a snapshot of an island’s current state in the context of its long-term trajectory. The research team developed Tier 1 resistance indicators based on subaerial island configuration and water level recurrence as a proxy for evaluating island resistance to storms, which are the primary driver of short-term change. These Tier 1 metrics can be analyzed over time to provide a high-level assessment of how an island’s resistance decreases because of elevation loss or sea-level rise or increases due to restoration or natural recovery. The research team developed Tier 2 resilience and resistance indicators and associated analyses to provide detailed information for specific time periods or applications (e.g., wildlife management). These metrics include habitat coverage from high-resolution maps, which show composition changes over time to capture the evolving resilience of specific habitat types; high tide flooding analysis, which evaluate island area relative to specified flooding thresholds to characterize resistance in the short-term or, if analyzed over time, indicate changes in resilience; and hypsometric curve analysis, which allows managers to evaluate island area changes above their own elevation benchmarks of interest and similarly characterize resistance in the short-term or indicate changes in resilience if assessed over time. The research team calculated Tier 1 metrics of the barrier islands and headlands along the coast of Louisiana for the period of 1984 through 2021 and Tier 2 metrics for select times during that period depending on data available and quality. The results were captured in a report card for each barrier, which also includes an overview of the metrics and their interpretation; a restoration and storm history; and Tier 1 and Tier 2 metric analysis, including benchmarking against coastwide and regional values as well as to an island’s pre-restoration trajectory. These report cards provide a readily digestible synthesis of barrier condition and trajectory that coastal managers can use to support restoration prioritization and other decisions.

Louisiana↗

USGS2024

Developments in African industrial minerals for renewable energy

Introduction Africa is emerging as a leading source for minerals used in the manufacture of batteries for electric vehicles and in other renewable energy applications. New graphite, lithium, and rare-earth mines have or could be opened in African countries from 2017 through 2026. Estimates of production capacities for graphite, lithium, and rare-earth mines for 2023 and beyond are based upon supply-side assumptions, such as announced plans for new capacity construction and bankable feasibility studies, as well as projected trends that could affect current producing facilities in 2023 and planned new facilities projected to come online by 2026. Forward-looking information, including estimates of future production capacities, graphite flake distributions, and timing of the start of operations, are subject to risk factors and uncertainties that could cause actual events or results to differ significantly from expected outcomes. Projects listed in this report are presented as an indication of industry plans and are not a U.S. Geological Survey (USGS) prediction of what will take place. Only projects with planned startup dates are included in this report; ther graphite, lithium, and rare-earth projects in Africa without startup dates were known to be in various stages of development but are not included in this fact sheet.

Fact Sheet↗

USGS2024

Earthquake relocations delineate discrete a fault network and deformation corridor throughout Southeast Alaska and Southwest Yukon

Deformation in southeastern Alaska and southwest Yukon is governed by the subduction and translation of the Pacific-Yakutat plates relative to the North American plate in the St. Elias region. Despite notable historical seismicity and major regional faults, studies of the region between the Fairweather and Denali faults are complicated by glacial coverage and the remote setting. In the last decade, significant improvements have been made to the density of regional broadband seismometer networks. We relocate more than 5,000 earthquakes between 2010 and 2021 in the region of southeastern Alaska and southwestern Yukon utilizing these improved seismic networks. With reductions in catalog uncertainty, particularly in depth, we quantify the thickness of the seismogenic layer in the crust throughout the region and locate seismicity on a shallow network of upper-crustal faults. Relocated earthquakes, combined with an updated focal-mechanism catalog, permit estimating and classifying motion of active faults. This includes mapping the Totschunda-Fairweather “Connector” fault, which plays an important role in explaining regional deformation, and identifying new faults like the Kathleen Lake fault. We draw similarities between our seismic observations and simplified conceptual models of regional tectonics, which describe a dominant transpressional regime and localized slip partitioning. Our results support a hypothesis where current deformation is taking place on a well-defined and evolved network of shallow faults in the corridor between the Totschunda-Fairweather “Connector” and Denali faults.

Alaska, Yukon↗

USGS2025

Forecast, monitor, adapt: A multi-agency strategy to protect people from postfire debris flows

In 2020, a wildfire burned across Glenwood Canyon in Colorado, USA. A history of postfire debris flows in the region and a hazard assessment for the burn area indicated that potentially life-threatening debris flows could be triggered by rainfall within months of a wildfire. As a result, four government agencies evaluated strategies to help mitigate hazards, including the loss of human life, that may be associated with debris-flow events. After the fire, 26 large debris flows occurred in the summer of 2021 and three sediment-laden flows occurred in the summer of 2023, but there were no major injuries or fatalities reported. We found that integrating hazard assessment/ forecasting, monitoring, and adaptation scenarios was a successful strategy for reducing postfire debris-flow risks to human life (including injuries and fatalities). Weather forecasts and estimates of debris-flow triggering rainfall thresholds, likelihood, and volume were used to anticipate the timing, location, and magnitude of debris-flow events. Rainfall monitoring and detailed recordkeeping of storms that triggered debris flows were used to validate and update debris-flow warning thresholds that varied with time following the wildfire. Although the governmental agencies working in this burn area had distinct and differing agency mandates, they were able to integrate information to reduce the risk of debris-flow events to human life.

Colorado↗

USGS2026

Competition between invasive round goby and native slimy sculpin in a coldwater laboratory setting

Slimy sculpin ( Cottus cognatus ) abundance has substantially declined in many regions of the Laurentian Great Lakes coincident with the arrival and range expansions of invasive round goby ( Neogobius melanostomus ). Previous studies have shown that round goby can outcompete native mottled sculpin ( C. bairdii ) for shelter and negatively affect native spoonhead sculpin ( C. ricei ) growth, but direct behavioral interactions of slimy sculpin and round goby have not been observed. We conducted a laboratory experiment in which a pair of slimy sculpin were placed in an aquarium that contained an artificial spawning shelter and one round goby (treatment, n = 10) or no other fish (control, n = 10) for at least 30 days at temperatures reflective of their potential offshore overlap. We measured slimy sculpin egg deposition, aggressive interactions, artificial shelter occupancy, growth, and survival. Slimy sculpin were less likely to occupy shelters and experienced increased intraspecific competition when a round goby was present. In contrast, round goby rarely chased slimy sculpin, slimy sculpin were more likely to chase round goby, and slimy sculpin growth, egg deposition, and survival were not significantly different between groups. Differences in density between treatments and the relative size of fish could confound our observations. However, the dominant slimy sculpin occupied the shelter at similar rates between experimental groups if the round goby was smaller. Our results indicate that larger round goby overlapping with slimy sculpin could affect slimy sculpin shelter and space use, but the Great Lakes environment may create more stress for slimy sculpin than we observed.

Lake Michigan, Lake Superior↗

USGS2000

Suspended-sediment characteristics in the Housatonic River Basin, western Massachusetts and parts of eastern New York and northwestern Connecticut, 1994-96

Suspended-sediment concentrations, discharges, loads, and yields were determined for eight subbasins in the Housatonic River Basin in western Massachusetts, eastern New York, and northwestern Connecticut from April 1994 through March 1996. Suspended-sediment samples were collected at three continuous-record sediment stations and at four partial-record sediment stations. Suspended-sediment concentrations in samples collected during the period of study ranged from less than 0.5 to 3,400 milligrams per liter, and concurrent streamflows ranged from 0.03 to 126 cubic feet per second per square mile at the seven stations. Median suspended-sediment concentrations in samples collected at each station ranged from 7 to 61 milligrams per liter. Median streamflows during suspended-sediment sampling ranged from 1.86 to 5.88 cubic feet per second per square mile. Instantaneous suspended-sediment yields ranged from less than 0.005 to 185 tons per day per square mile, and medians ranged from 0.03 to 1.12 tons per day per square mile at the seven stations. Total suspended-sediment loads (mass) from April 1994 through March 1996 at the continuous-record sediment stations were 11,603 tons at Housatonic River near Great Barrington, 7,929 tons at Green River, and 54,347 tons at Housatonic River near Ashley Falls. Suspended-sediment load during January 1996 at the Green River station accounted for about 54 percent of the total suspended-sediment load for the Green River during the 2-year study. Suspended-sediment load on January 19 and 20, 1996, at the Green River station accounted for about 50 percent of the January 1996 suspended-sediment load, or about 27 percent of the total suspended-sediment load during the 2-year study. This large suspended-sediment transport was the result of rainfall and snowmelt on January 19 and 20, 1996--the equivalent of a 5- to 6-inch rain storm in the Green River subbasin. Total suspended-sediment loads during the 2-year study at the partial-record sediment stations were 3,052 tons at Williams River, 1,758 tons at Ironworks Brook, and 17,927 tons at Konkapot River. Suspended-sediment yields from April 1994 through March 1996 at the continuous-record sediment stations were 21 (tons/yr)/mi 2 at Housatonic River near Great Barrington, 78 (tons/yr)/mi 2 at Green River, and 58 (tons/yr)/mi 2 at Housatonic River near Ashley Falls. Suspended-sediment yields during the 2- year study at the partial-record sediment stations were 35 (tons/yr)/mi 2 at Williams River, 78 (tons/yr)/mi 2 at Ironworks Brook, and 147 (tons/yr)/mi 2 at Konkapot River. Suspended-sediment yields were estimated for two subbasins in the Housatonic River Basin--Schenob Brook at Sheffield, a partial-record sediment station, and the area adjacent to the Housatonic River between Great Barrington and Ashley Falls. The estimate of suspended-sediment yield at Schenob Brook at Sheffield of 82 (tons/yr)/mi 2 is comparable to the yield determined for the Green River and Ironworks Brook. The estimate of suspended-sediment yield for the area adjacent to the Housatonic River between Great Barrington and Ashley Falls was 395 (tons/yr)/mi 2 . This estimated suspended-sediment yield was 2.7 to 18.8 times greater than that estimated for any of the other subbasins. Several basin and land-use characteristics thought to affect suspended-sediment transport in the subbasins were compared to the suspended-sediment yields. The characteristics that seemed to affect suspended-sediment discharge were dams, which contributes to decreased yields; Hadley, Limerick, Linlithgo, Saco, and Winooski silt loam soils (high erodibility soils), which contributes to increased yields; stratified-drift deposits, which contributes to increased yields; and agricultural and open land, which contributes to increased yields. The effect of stratified-drift deposits on suspended-sediment discharge is thought to be greater when those deposits are of glaciolacustrine (generally clay, silt, and fine sand), rather than glaciofluvial (clay, silt, sand, gravel, and cobbles) origin. The silt loam soils, glaciolacustrine deposits, and agricultural and open land are interrelated, inasmuch as the silt loam soils generally are associated with glaciolacustrine deposits and agricultural activities.

Connecticut, Massachusetts, New York↗

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↗

USGS1989

Techniques for simulating flood hydrographs and estimating flood volumes for ungaged basins in east and west Tennessee

A dimensionless hydrograph developed for a variety of basin conditions in Georgia was tested for its applicability to streams in East and West Tennessee by comparing it to a similar dimensionless hydrograph developed for streams in East and West Tennessee. Hydrographs of observed discharge at 83 streams in East Tennessee and 38 in West Tennessee were used in the study. Statistical analyses were performed by comparing simulated (or computed) hydrographs, derived by application of the Georgia dimensionless hydrograph, and dimensionless hydrographs developed from Tennessee data, with the observed hydrographs at 50 and 75% of their peak-flow widths. Results of the tests indicate that the Georgia dimensionless hydrography is virtually the same as the one developed for streams in East Tennessee, but that it is different from the dimensionless hydrograph developed for streams in West Tennessee. Because of the extensive testing of the Georgia dimensionless hydrograph, it was determined to be applicable for East Tennessee, whereas the dimensionless hydrograph developed from data on streams in West Tennessee was determined to be applicable in West Tennessee. As part of the dimensionless hydrograph development, an average lagtime in hours for each study basin, and the volume in inches of flood runoff for each flood event were computed. By use of multiple-regression analysis, equations were developed that relate basin lagtime to drainage area size, basin length, and percent impervious area. Similarly, flood volumes were related to drainage area size, peak discharge, and basin lagtime. These equations, along with the appropriate dimensionless hydrograph, can be used to estimate a typical (average) flood hydrograph and volume for recurrence-intervals up to 100 years at any ungaged site draining less than 50 sq mi in East and West Tennessee.

Tennessee↗

USGS2026

Preface to the focus section on intraplate earthquakes

More than a half century after plate tectonics provided an overarching framework to explain earthquakes along active plate boundaries, numerous theories have been proposed to explain where, why, and how often earthquakes occur well away from active plate boundaries, but a paradigm remains elusive. Even the classification of earthquakes away from active plate boundaries as "intraplate" raises issues, with potentially important distinctions between Stable Continental Regions and more actively deforming regions including passive margins and failed rifts. Some of the largest known intraplate earthquakes themselves remain enigmatic, having occurred before the modern instrumental era. Hazard assessments are often data-limited: low fault-slip rates relative to landscape modification rates result in poor discoverability of fault sources, challenging the characterization of source zones and earthquake recurrence; the completeness and homogenization of instrumental earthquake catalogs using uncertain magnitude conversions can lead to uncertainties in earthquake recurrence; and, limited strong-motion observations for large-magnitude events at near-source distances leads to uncertainties in the selection and development of ground-motion models for seismic hazard studies. Data from recent intraplate earthquakes around the world—from the moment magnitude M 7.7 2001 Bhuj, India, earthquake 25 years ago to the 2024 M 4.8 Tewksbury, New Jersey earthquake—have yielded both new insights and new questions. The papers in this special focus discuss many of the long-standing challenges involved with intraplate earthquake investigations and provide a snapshot of the state of the art with current research to advance our understanding.

Seismological Research Letters↗

USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗

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↗

USGS2025

Earthquake recurrence estimates for northern Caribbean faults from combinatorial optimization

We use combinatorial optimization to find the optimal spatial distribution of random samples of earthquakes (≥6.5) that minimize the misfit in target slip rates for all faults in the northeast Caribbean, and we derive magnitude-frequency relationships with uncertainties for these faults. Slip rates for many faults are derived from GPS block models, not direct measurements, because of their underwater locations. Predicted recurrence rates for eastern Hispaniola and Puerto Rico Trench faults are 220-450 yr for M7 and 3-5 kyr for M8, with maximum feasible magnitude of M8.2. The most frequent earthquakes with magnitudes ≥7.0 are predicted on the large upper plate strike-slip faults, Enriquillo (EF) and Septentrional Fault, commensurate with the historical record. Calais et al. (2023) suggested that shortening in western Hispaniola is accommodated on the offshore Jérémie and onshore Malpasse faults north and south of EF, instead of on terrestrial faults in western Hispaniola and EF. Because of our system-modeling approach, such a configuration predicts less frequent earthquakes on EF and on western Hispaniola and Muertos convergent zones. Recurrence times of a few 100s yr for M6.7 earthquakes is predicted on the submerged faults in Mona Passage, and infrequent M>7 earthquakes are predicted on the Virgin Islands faults.

Lesser Antilles↗

USGS2025

RioM-1: A new calcite reference material for U-Pb LA-ICP-MS geochronology

Determining absolute ages of carbonate diagenesis, faulting, fossil formation, speleothem growth, carbonate-hosted hydrocarbon deposits, vein mineralisation and hydrothermal alteration has become increasingly accessible through LA-ICP-MS U-Pb dating of calcite, complementing traditional isotope dilution methods still applicable to certain materials (e.g., speleothems via micro-drilling). However, well-calibrated reference materials for LA-ICP-MS calcite geochronology remain scarce. Here, we characterise the Rio Maior calcite, designated ‘RioM-1’, as a potential reference material for U-Pb dating by LA-ICP-MS. Fragments (0.1 to 1 cm 3 ) from a single scalenohedral crystal were analysed by ID-TIMS (U-Pb), LA-ICP-MS (U-Pb and 87 Sr/ 86 Sr), and SIMS (O isotopes). RioM-1 displays high U mass fraction and low, though variable, proportions of common Pb. Combined ID-TIMS analyses from two independent laboratories yielded a Tera-Wasserburg lower intercept date of 63.93 ± 0.11 Ma (2 s , MSWD = 1.3, n = 16). LA-ICP-MS U-Pb data from three independent laboratories are concordant with the TIMS age, producing a pooled date of 64.10 ± 0.12/1.2 Ma (2 s , MSWD = 7.7, n = 708) and an initial 207 Pb/ 206 Pb ratio of 0.85 ± 0.01. Other calcite reference materials, when normalised to RioM-1, yielded dates within uncertainty of their published values. SIMS measurements returned a mean 18 O/ 16 O of 0.002044450 ± 181 (1 s ) and δ 18 O SMOW of 19.57 ± 0.92‰ (1 s ), while LA-ICP-MS strontium isotope measurement yielded a mean 87 Sr/ 86 Sr of 0.708177 ± 9 (2 s ).

Geostandards and Geoanalytical Research↗

USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗

USGS2025

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗