SEARCH · USGS Science

Results for “Journal of Geophysical Research: Solid Earth”

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,663 records · Page 75Linked to original sources

USGS2026

Mineral commodity summaries 2026

Introduction Each mineral commodity chapter of the 2026 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2025 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the 2026 MCS are welcomed.

Mineral Commodity Summaries↗

USGS2025

Relationship between peak and cumulative ground motions from 49 Mw 3-6 earthquakes in the San Francisco Bay Area

We examine the relationship between peak ground velocity (PGV) and cumulative absolute displacement (CAD) for 49 M w 3 – 6 earthquakes in the San Francisco Bay Area (SFBA) and gain insight into the spatiotemporal partitioning of seismic energy in ground motion records with respect to source, path, and site effects. PGV and CAD are positively correlated, but there can be large deviations from the average trend. For example, ground motion records with either very long duration resonance or short pulse-like motions will have higher or lower CAD, respectively, but could have very similar PGV. We perform principal component analysis (PCA) on PGV-CAD for >7000 records in the SFBA with the goal of investigating what influences positive or negative anomalies in cumulative motions. PCA rotates the PGV-CAD datapoints into two principal components, where the one with the larger variance, which we call the “primary intensity” represents mostly the distance-dependence of ground motion amplitudes. The other principal component, which we call the “excess motion”, represents the deviation from cumulative motions that would be typical for a ground motion record with a given PGV. The excess motion will be positive in the case of records with long duration ringing and will be negative for short duration pulse-like ground motions. We find that excess motion is generally positive at sites in sedimentary basins and in soft sediments around the SF Bay. Excess motion is generally negative in the very near field, as well as at sites on hard bedrock. We discuss the findings here in terms of implications for seismic hazard applications and other wave propagation phenomena.

California↗

USGS2024

Iodine-129 in the eastern Snake River Plain aquifer at and near the Idaho National Laboratory, Idaho, 2021–22

Between the 1950s and 1980s, wastewater generated at the Idaho National Laboratory contained Iodine-129 ( 129 I); this wastewater was discharged directly into the eastern Snake River Plain (ESRP) aquifer through a deep disposal well, unlined infiltration ponds, or leaked from distribution systems below industrial facilities. During 2021–22, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy and the Idaho Department of Environmental Quality Idaho National Laboratory Oversight Program, collected groundwater samples from 64 monitoring wells in the ESRP aquifer, 6 of which are part of a multilevel monitoring system, to determine the concentration of 129 I in the groundwater. These samples were analyzed by accelerator mass spectrometry as part of a long-term ongoing study to track trends and occurrences of this carcinogenic, long-lived radionuclide in the environment. Concentrations ranged from slightly above the locally determined background concentration of 5.4×10 −6 picocuries per liter, to just below the U.S. Environmental Protection Agency’s maximum contaminant level of 1 picocurie per liter. Discharge of wastewater containing 129 I has been discontinued to the aquifer, and long-term trends from a subset (n=15) of sampled wells show decreasing 129 I concentrations over the last three decades. Concentrations of 129 I in groundwater from monitoring wells near facilities at the Idaho National Laboratory are affected by episodic recharge from an ephemeral surface-water source and by the fracture-flow dominated hydrologic regime in the ESRP aquifer. The spatially focused sampling effort has also identified a low-level 129 I plume that affects long-term water quality near and downgradient from the Advanced Test Reactor Complex in the southwestern part of the facility that had not been clearly defined in previous sampling efforts, although the definition of the plume is somewhat limited by available data.

Idaho↗

USGS2025

Lithium from magma to mine in an early Yellowstone hotspot caldera

Renewable energy technologies rely on the extraction of metals not historically in high demand, such as lithium (Li), for which ore deposit models are incompletely understood. One of the world’s largest Li deposits is hosted in lake sediments of the 16.4 Ma McDermitt caldera, which formed during the early stages of Yellowstone hotspot volcanism in the western United States. Eruptive and posteruptive mobility of Li are major challenges in elucidating deposit formation. Melt inclusions preserved in quartz crystals provide a means to assess pre-eruptive magmatic Li contents. Concentrations of Li determined by ion microprobe for melt inclusions in a McDermitt rhyolite lava are 400−1350 ppm, compared to 20−70 ppm Li in matrix rhyolite glasses. Synthesis with melt inclusion data for eight additional calderas demonstrates a recurrence of Li-rich rhyolitic magmas (200−2000 ppm Li) in the western part of the Yellowstone hotspot track. However, unlike the multicyclic caldera complexes with overlapping fault networks that may have compromised Li retention, the McDermitt caldera remained a closed hydrologic system throughout its evolution. Modeling indicates 100 km 3 of resurgent magma could yield 25−150 Mt Li in a magmatic fluid and supports accumulation of Li-rich magmatic fluid in a closed intracaldera lake, followed by evaporative concentration and sequestration of Li within clay minerals to generate the McDermitt deposit.

Idaho, Nevada, Oregon, Wyoming↗

USGS2024

Simulated mean monthly groundwater-transported nitrogen loads in watersheds on the north shore of Long Island Sound, 1993–2022

Elevated nitrogen loads are pervasive in the Long Island Sound, an estuary that receives freshwater and nutrients from both surface-water and groundwater discharge. Surface-water nitrogen loads to the Long Island Sound are relatively well characterized, but less is known about groundwater-transported nitrogen loads. Prior work on the northern shore of Long Island Sound (Connecticut and areas of New York and Rhode Island) suggested that groundwater travel times are relatively short (median less than 2 years) and that decade-long nutrient legacies are not widespread. Because the travel times are short, groundwater flow and nutrient loads likely vary substantially between months. In the current study, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency’s Long Island Sound Study and the Connecticut Department of Energy and Environmental Protection, developed a set of models to better characterize spatial and temporal patterns of groundwater-transported nitrogen loading from atmospheric deposition, septic systems, and fertilizers within the study area. The models provide an estimate, with uncertainty, of groundwater-transported nitrogen loads in the study area, filling a key gap in the nitrogen budget for Long Island Sound. The models also highlight the spatial and temporal variation in nitrogen loading throughout the study area. The modeling workflow involved four models. (1) A soil-water-balance model was developed by using the Soil-Water-Balance software to simulate groundwater recharge across the study area for water years 2005 through 2022. The simulated mean monthly recharge from the soil-water-balance model was used as input into a groundwater-flow model. (2) The groundwater-flow model was developed by using the MODFLOW 6 software and data for water years 1993 through 2022 and simulates average monthly hydrologic conditions. The groundwater-flow model was calibrated by using the Iterative Ensemble Smoother method within the PEST++ software. The Iterative Ensemble Smoother method generates an ensemble of sets of parameter values, with each set producing reasonable simulated hydrologic parameter values. (3) An ensemble of MODPATH particle-tracking simulations were run to generate particle flow paths and travel times, with each simulation using a different set of the flow model parameters. (4) A nitrogen load model uses the MODPATH simulation outputs to track nitrogen from the land surface through multiple attenuation zones until it discharges into fresh or saline surface water. As with the groundwater-flow model, the nitrogen model simulated average monthly groundwater-transported nitrogen loads for water years 1993 through 2022. One novel aspect of the nitrogen load model is that the nitrogen attenuation parameters were calibrated to observed nitrogen loads. Across the ensemble of simulated nitrogen loads, the median study-area-wide monthly simulated nitrogen loads from the aquifer to Long Island Sound throughout the year ranged from 900 to 18,600 kilograms of nitrogen per day, with a median load of 5,100 kilograms of nitrogen per day. The simulated loads were based on average monthly conditions for water years 1993 through 2022. Loads were highest during the winter and early spring and lowest during the late summer. However, simulated travel times for groundwater and nitrogen loads discharged to Long Island Sound during summer were longer than travel times for groundwater and loads discharged during the winter, indicating that, on average, groundwater discharged during summer traveled along different, and longer, flow paths, than groundwater discharged during winter. This indicates that summer loads would respond more slowly to changes in nitrogen inputs at the water table than winter loads. Over the entire study area, approximately 15 percent of the simulated load is from atmospheric deposition sources, 30 to 40 percent is from fertilizer, and 50 to 60 percent is from septic systems. The final analysis of the study involved simulating the change in groundwater-transported nitrogen load in response to upgrading septic systems or reducing fertilizing inputs to areas of turf grass. Both management interventions reduced the groundwater-transported nitrogen load, and reductions were greater in areas with greater loads from septic systems or turf-grass fertilizers. The delay between management actions and substantial reductions in groundwater-transported nitrogen loads varied seasonally; loads during the late summer months remained elevated longer than the winter loads.

Connecticut, Rhode Island↗

USGS2000

Water-quality conditions and relation to drainage-basin characteristics in the Scituate Reservoir Basin, Rhode Island, 1982-95

The Scituate Reservoir Basin covers about 94 square miles in north central Rhode Island and supplies more than 60 percent of the State of Rhode Island's drinking water. The basin includes the Scituate Reservoir Basin and six smaller tributary reservoirs with a combined capacity of about 40 billion gallons. Most of the basin is forested and undeveloped. However, because of its proximity to the Providence, Rhode Island, metropolitan area, the basin is subject to increasing development pressure and there is concern that this may lead to the degradation of the water supply. Selected water-quality constituent concentrations, loads, and trends in the Scituate Reservoir Basin, Rhode Island, were investigated locate parts of the basin likely responsible for exporting disproportionately large amounts of water-quality constituents to streams, rivers, and tributary reservoirs, and to determine whether water quality in the basin has been changing with time. Water-quality data collected between 1982 and 1995 by the Providence Water Supply Board PWSB) in 34 subbasins of the Scituate Reservoir Basin were analyzed. Subbasin loads and yields of total coliform bacteria, chloride, nitrate, iron, and manganese, estimated from constituent concentrations and estimated mean daily discharge records for the 1995 water year, were used to determine which subbasins contributed disproportionately large amounts of these constituents. Measurements of pH, color, turbidity, and concentrations of total coliform bacteria, sodium, alkalinity, chloride, nitrate, orthophosphate, iron, and manganese made between 1982 and 1995 by the PWSB were evaluated for trends. To determine the potential effects of human-induced changes in drainage- basin characteristics on water quality in the basin, relations between drainage-basin characteristics and concentrations of selected water-quality constituents also were investigated. Median values for pH, turbidity, total coliform bacteria, sodium, alkalinity, chloride, nitrate, and iron were largest in subbasins with predominately residential land use. Median instantaneous loads reflected drainage-basin size. However, loads normalized by drainage area (median instantaneous yields) also were largest in residential areas where point and non-point sources are likely, and in areas of poorly drained soils. Significant trends in water-quality constituents from 1982 to 1995 in the Scituate Reservoir Basin indicate that the quality of the water resources in the basin may be slowly changing. Scituate Reservoir subbasins with large amounts of residential land use showed increasing trends in alkalinity and chloride. In contrast, subbasins distributed throughout the drainage basin showed increasing trends in pH, color, nitrate, and iron concentrations, indicating that these characteristics and constituents may be affected more by atmospheric deposition. Although changing, water-quality constituent concentrations in the Scituate Reservoir Basin only occasionally exceeded Rhode Island and USEPA water-quality guidelines and standards. Result of correlation analysis between pH, color, turbidity, and concentrations total coliform bacteria, sodium, alkalinity, chloride, nitrate, orthophosphate, iron, and manganese and land use, geology, wetlands, slope, soil drainability, and roads indicated that the percentage of wetlands, roads, and slope appear have the greatest effect on water-quality in the Scituate Reservoir Basin. The percentage of urban, residential, and commercial land use also are important, but to a lesser degree than wetlands, roads, and slope. Finally, geology appears to have the least effect on water quality compared to other drainage-basin characteristics investigated.

Rhode Island↗

USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗

USGS2024

Estimating groundwater level records using MOVE.1 and computing monthly percentiles from estimated groundwater records in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, performed record extensions on groundwater levels at select wells using the Maintenance of Variance Extension type 1 (MOVE.1) method. The groundwater levels estimated from these record extensions were used to compute monthly percentiles to improve future determinations of a groundwater index. In Massachusetts, 27 of 29 short-record study wells with continuous groundwater levels between 0.8 and 8.1 years were suitable for record extensions; 37 long-record index wells were used to extend the groundwater level records at the study wells. The index well selected to pair with a study well was chosen based on Pearson correlation coefficient values; cross-correlation between the two wells; geologic and topographic similarity; and smallest distance spanning the wells. Each study well and its corresponding index well have 1 or more years of concurrent, overlapping data; a Pearson correlation coefficient that exceeded a threshold value of 0.8; and a similar aquifer type and hydrologic characteristics. Of the 29 study wells, 2 showed poor correlations with all index wells and were not considered for record extensions. Performance metrics used to assess the accuracy of the MOVE.1 models indicated that most models provided reasonable estimates of groundwater levels. Root mean square error values ranged from 0.097 to 2.292 feet, with a median of 0.536 foot. Nash-Sutcliffe efficiency coefficient values ranged from 0.623 to 0.996, with a median value of 0.759. Generally, study wells in close geographical proximity to their index well resulted in stronger model performance. The average length of groundwater level records was extended by 14.1 years to a new average of 18.1 years. The estimated groundwater level records from the MOVE.1 models resulted in an increase in the range of highest and lowest groundwater levels at 23 of 27 wells. The increase in range of groundwater levels was between 0.08 to 7.95 feet. Monthly percentiles for State drought indices were computed from the estimated MOVE.1 records and observed records through December 31, 2021. Percentiles computed from estimated records show an average groundwater level about 1.0 foot lower than observed data at the 2d percentile and 0.1 foot lower at the 30th percentile.

Massachusetts↗

USGS2024

Seismic techniques and suggested instrumentation to monitor volcanoes

Introduction Changes in the pressure or location of magma can stress or break surrounding rocks and trigger flow of nearby waters and gases, causing seismic signals, such as discrete earthquakes and tremor. These phenomena are types of seismic unrest that commonly precede eruption and can be used to forecast volcanic activity. Mass movements at the surface, including avalanches, debris flows, and lahars, may also generate seismic signals that are specifically addressed in chapter H , this volume (Thelen and others, 2024). Our focus in this chapter is to determine the levels of instrumentation recommended to produce high-quality, well-constrained seismic observations important for early warning of impending eruptions, detecting changes in ongoing eruptions, and characterizing other hazardous volcanic events. There are emerging techniques and new types of instrumentation, such as distributed acoustic sensing or rotational seismometers, that we do not consider here. These types of instrumentation show promise for monitoring but still require maturation before being considered more generally in volcano monitoring. Most of the capabilities mentioned below are universal for all types of volcanic systems, although some are best applied to stratovolcanoes with an apical single vent. In some settings, such as calderas or shield volcanoes, we must broaden coverage to include multiple possible storage regions or vent locations. As an example, Thelen (2014) discretized the long rift zones of shield volcanoes in Hawaiʻi as a set of evenly spaced “vents.” In this construct, each vent comes with recommendations, and several thousand network configurations were simulated to assess the effect on network quality levels and to determine the most efficient network design. The same process could be applied in a caldera setting or a volcanic field, where an evenly spaced grid of potential vents is considered. Localized recommendations for each unique system are beyond the scope of this report and left up to local experts to assess based on the conditions, restrictions, and requirements of each volcano.

Scientific Investigations Report↗

USGS2025

Tracing fluid evolution during batholith and deposit formation: In-situ and bulk rock Mo isotopes from the Questa porphyry Mo deposit, New Mexico, USA

The Oligocene Questa porphyry Mo deposit, located in New Mexico, USA, is a well-characterized system that formed over 400 kyr, punctuating the 8 Myr formation of the Questa-Latir batholith. Detailed sampling of the intrusions and minerals associated with magmatichydrothermal mineralization, as well as bookending barren intrusions, allows for examination of the Mo isotope record of the deposit and its host batholith. Through this established chronologic framework, such data can yield insights into the metallogenesis of deposit formation and long-term magma genesis evolution for upper crustal magmatism. The Mo isotopic compositions of barren intrusions yield a younging trend towards lighter δ 98 Mo. However, intrusions associated with Mo deposit formation yield heavier δ 98 Mo relative to barren magmatism, by up to 2.0‰. This indicates that they were affected by an isotopically heavy and Mo-rich fluid. The in-situ data from molybdenite covers the range observed in mineralized bulk rocks. However, individual samples consistently yield Mo isotopic variations of 0.4‰ to 0.5‰. The consistency of the δ 98 Mo variation in in-situ analyses of molybdenite within hand samples and individual mineral grains is likely the result of MoS 2 crystallization in a closed system and reflects an almost complete precipitation of Mo.

New Mexico↗

USGS2026

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li 2 O, respectively, for global Li pegmatite deposits ( n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs 2 O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb 2 O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Ore Geology Reviews↗

USGS1995

Streamflow, ground-water recharge and discharge, and characteristics of surficial deposits in Buzzards Bay Basin, southeastern Massachusetts

Streamflow measurements at 14 low-flow partial-record stations and two discontinued streamflow-gaging stations were related to concurrent streamflows at six long-term streamflow-gaging stations to estimate streamflows at selected flow durations from 50 to 99 percent for water years 1967-91. At low flows, stream discharge per square mile generally increased with increasing percentage of subbasin underlain by stratified-drift deposits. At the 70-, 90-, and 99-percent flow durations, subbasins underlain primarily by stratified-drift deposits had stream discharges per square mile at least 2, 4, and 8 times greater, respectively, than subbasins underlain primarily by till and bedrock deposits. Streamflow measured three times at the Paskamanset River showed that streamflows downstream of municipal pumped wells increased at a lower rate per unit drainage area than streamflows upstream of the pumped wells. When pumpage from wells near the river was added to measured streamflow at the most downstream station, the sum was similar to the potential streamflow at the station if its streamflow per unit area was the same as that upstream from the pumped wells. Mean ground-water recharge and discharge rates were computed from continuous records of daily mean discharge during water years 1967-91 for six streamflow-gaging stations in southeastern Massachusetts and Rhode Island. Estimates of mean ground-water recharge were 19.7 to 22.6 inches per year for stations with drainage areas primarily underlain by till and bedrock deposits, and 23.8 to 25.2 inches per year for stations with drainage areas primarily underlain by stratified-drift deposits. During drought years, annual ground-water recharge can be less than 50 percent of the mean ground-water recharge rate for water years 1967-91. Estimates of ground-water discharge rates generally were 1 to 4 inches per year less than estimates of ground-water recharge rates. Ground-water discharge-duration curves were calculated from hydrograph separation of streamflow records at the six continuous streamflow-gaging stations for water years 1967-91. When streamflows were less than the median streamflow (50-percent duration), ground-water discharges averaged 85.7 percent of the total streamflow at the six streamflow-gaging stations. Ground-water discharges were estimated at the 14 low-flow partial-record stations and two discontinued streamflow-gaging stations by multiplying streamflows at the selected flow durations from 50 to 99 percent for water years 1967-91 by 85.7 percent. Data on surficial deposits in the basin were obtained by drilling 10 wells, by conducting 12 seismic-refraction surveys, and by collecting well log data for 313 wells from consulting firms and the Massachusetts Department of Environmental Management. These data indicated that a few minor corrections could be made to existing saturated thickness maps for stratified-drift deposits in Buzzards Bay Basin. Few well logs were detailed enough to estimate transmissivities; thus existing transmissivity maps were not updated for stratified-drift deposits in the basin.

Massachusetts↗

USGS2000

Stormflow chemistry in the Santa Ana River below Prado Dam and at the diversion downstream from Imperial Highway, southern California, 1995-98

The Santa Ana River drains about 2,670 square miles of the densely populated coastal area of southern California, near Los Angeles. Almost all the flow in the river, more than 200,000 acre-feet annually, is diverted into ponds where it infiltrates and recharges underlying aquifers. About 2 million people are dependent on these aquifers for water supply. Stormflow in the Santa Ana River is considered a source of "high-quality" water suitable for use as a source of ground-water recharge. To test this assumption, stormflow samples were collected at two locations--below Prado Dam and at the diversion point downstream from Imperial Highway--for 12 winter storms between 1995 and 1998. Nitrate concentrations decreased during stormflow from a median concentration of 7.8 milligrams per liter in base flow to concentrations less than 1 milligram per liter in some large storms. Concentrations of chemically reduced forms of nitrogen (nitrite, ammonia, and organic nitrogen) increased during stormflow and are the predominant forms of nitrogen in large stormflows. Dissolved organic carbon (DOC) concentrations increased from a median concentration of 4.6 milligrams per liter in base flow to more than 20 milligrams per liter in some stormflows. Concentrations of DOC were especially high during the first storm of the rainy season, and large increases in DOC concentrations were measured even as a result of small early season storms that did not cause large increases in streamflow. DOC present during early season stormflow had less ultraviolet absorbance at 254 nanometers (UV 254 ) per unit of carbon than did DOC from late season stormflows. DOC in water held in storage behind Prado Dam had the highest UV 254 absorbance per unit of carbon. Maximum pesticide concentrations in stormflow did not exceed U.S. Environmental Protection Agency Maximum Contaminant Levels. Most pesticide concentrations were less than 1 microgram per liter and less than the detection limits obtained using standard drinking water analyses. Increases in concentrations of pesticides such as diazinon, malathion, and chlorpyrifos in stormflow result from runoff from urban areas downstream from Prado Dam. In general, large late season stormflows have the most pesticide detections of all stormflows sampled. Concentrations of methyl tert-butyl ether (MTBE), a gasoline additive, during base flow were as high as 0.9 microgram per liter and concentrations decreased during stormflow. Like pesticides, the concentrations did not exceed the U.S. Environmental Protection Agency Maximum Contaminant Levels for MTBE.

California↗

USGS2025

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter↗

USGS2025

Critical minerals in mine waste

Introduction Critical minerals are commodities with vulnerable supply chains that play a vital role in supporting the United States’ economy, national defense and security, emerging technologies, and energy independence. The prosperity of our Nation depends on generating a resilient supply of domestic critical minerals; mine waste may be an untapped source of these commodities. Mine waste from centuries of legacy mining persist on the landscape and may contain critical minerals and other valuable commodities previously deemed uneconomic to recover. At modern mines, the financial viability of recovering byproduct critical minerals, which are not the primary target, may be marginal and can ultimately destine them to mine waste. Further, mine waste can be a liability for the mining company or, at legacy mines, the taxpayer because of its effect on the landscape. The U.S. Geological Survey (USGS) has several initiatives to evaluate critical mineral resources in various types of waste. This factsheet highlights studies of mine waste carried out by USGS scientists at the Geology, Energy & Minerals Science Center in collaboration with other science centers funded through the USGS Mineral Resources Program. Recovery of critical minerals from mine waste can aid in remediation efforts and increase domestic supply of vital mineral resources.

Fact Sheet↗

USGS2026

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗

USGS2025

Water use across the conterminous United States, water years 2010–20

Withdrawals of water for human use are fundamental to the evaluation of the Nation’s water availability. This chapter provides an analysis of public supply, crop irrigation, and thermoelectric power water use for the conterminous United States (CONUS) during water years 2010–20. These three categories account for about 90 percent of water withdrawals in the Nation. The values presented here are based on modeling approaches that estimate water use at temporal (monthly) and spatial scales (12-digit hydrologic unit code—small watersheds sized 50–100 square kilometers) compatible for integration into a broader national assessment of water availability. Models also provide an understanding of factors that influence water use. An estimated 244,817 million gallons per day (Mgal/d; 28,677 million cubic meters per month [Mm 3 /mo]) were withdrawn on average within the CONUS during water years 2010–20 from fresh water and saline water for crop irrigation, public supply, and thermoelectric power, with shares of 43, 14.5, and 42.5 percent for each of these categories, respectively. In the same period, estimated withdrawals and consumptive use (1) for public supply were 35,400 and 4,219 Mgal/d (4,081 and 486 Mm 3 /mo), respectively; (2) for crop irrigation were 105,497 and 75,698 Mgal/d (12,147 and 8,716 Mm 3 /mo), respectively; and (3) for thermoelectric power from fresh water were 82,656 and 2,904 Mgal/d (9,952 and 345 Mm 3 /mo), respectively. Withdrawals for these categories of water use are highly spatially variable, with western States dominated by crop irrigation and eastern States dominated by thermoelectric-power water use. Public supply accounts for the largest percentage of water use in several heavily populated northeastern States. Reliance on groundwater compared to surface water depends on the availability of water sources and the type of water use. For public supply, withdrawals from groundwater are greater than withdrawals from surface water in the Western aggregated hydrologic regions, whereas the balance shifts to more surface water for the rest of the CONUS. In all aggregated hydrologic regions, the predominant source of water for crop irrigation is groundwater. Most thermoelectric power facilities in the eastern half of the CONUS use surface water from freshwater and saline sources; most thermoelectric power facilities in the western half of the CONUS use groundwater.

conterminous United States↗

USGS2025

How do hydrological variability and human activities control the spatiotemporal changes of riverine nitrogen export in the Upper Mississippi River Basin?

Excessive nitrogen export from agricultural watersheds remains a critical water quality challenge, with the Upper Mississippi River Basin (UMRB) significantly contributing to downstream eutrophication and hypoxia in the Gulf. This study investigates the spatiotemporal dynamics of riverine nitrate plus nitrite (NO 3 – + NO 2 – -N) export across the UMRB at high spatial resolution (12-digit Hydrologic Unit Codes or HUC12 subwatershed scale) during 2001–2020 and quantifies the effects of anthropogenic activities and hydrological variability on riverine NO 3 – + NO 2 – -N export changes in the region between 2001–2005 and 2016–2020. Our results revealed hotspots of substantial increases in NO 3 – + NO 2 – -N yields across the UMRB, with distinct regional patterns in driving factors. Over the entire UMRB, NO 3 – + NO 2 – -N yields increased by 9.7 kg/ha/yr on average from 2001–2005 to 2016–2020, with anthropogenic activities contributing 4.8 kg/ha/yr and hydrological variability contributing 4.9 kg/ha/yr. The northern and western UMRB had combined influences from both anthropogenic activities and hydrological variability, while the east-central regions had predominantly hydrologically driven changes. Agricultural sources, including fertilizer, manure, and biological nitrogen fixation, collectively contributed over 80% of NO 3 – + NO 2 – -N loading throughout the basin. This framework for disentangling human and hydrological impacts provides critical insights for developing effective and targeted watershed management strategies to reduce nutrient losses and improve water quality.

Upper Mississippi River Basin↗