USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of the Mechanics and Physics of Solids”

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

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2026

An improved empirical model for predicting postfire debris-flow volume in the western United States

Reliable estimates of debris-flow volume can be used to help predict the magnitude of debris-flow hazards following wildfire in the western United States. In this study, we compiled and used a database of 227 postfire debris-flow volumes that were collected across the western United States to develop a multiple linear regression model for predicting postfire debris-flow volume. We explored 36 predictor variables related to rainfall, terrain, and fire characteristics, and selected the model with the combination of variables that yielded the most accurate predictions of debris-flow volume. We evaluated model performance against the entire volume database, as well as against four subsets of volume data from southern California, the Intermountain West, the Southwest, and regions with limited volume data, such as northern California and Washington. We also compared model performance against 3 existing postfire debris-flow volume models that were developed for use in southern California, the Intermountain West, and the Southwest. We demonstrate that the new volume model performs as well as the regional models in the regions for which they were developed and outperforms existing models when applied to volumes from data-limited regions in the western United States. These results indicate that the debris-flow volume model introduced in this study can be used to improve postfire hazard assessments across the western United States, especially outside of southern California.

Arizona, California, Colorado, New Mexico, Utah, W↗
USGS2026

3D semantic mapping of surface geological features

Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.

Computers & Geosciences↗
USGS2026

Uranium—Deposits, production and resources, market dynamics, and supply chain risks

Introduction Interest in nuclear power for the generation of electricity has risen with the increase in the need for more diverse baseload power, enhanced energy security, and the development of new technologies, such as small modular reactors (SMRs), which could provide power for remote areas, industrial applications, and artificial intelligence (AI) data centers. In 2024, the U.S. Department of Energy received $2.7 billion in congressional funding to bolster the domestic uranium production and nuclear fuel supply chain and address reliance on imports from foreign suppliers. In 2025, the U.S. Government issued several Executive and Secretary’s orders aimed at revitalizing the U.S. nuclear sector. If SMRs are to be as widely deployed in the United States and worldwide as envisioned, demand for uranium (nuclear reactor fuel) will likely increase. After the Fukushima nuclear accident in 2011, the market spot price of uranium began a decline, followed by a decrease in U.S. and global uranium exploration and mine development expenditures that led to a uranium supply deficit until 2020, when prices started to recover, prompting a resurgence in uranium exploration and development. In January of 2024, the uranium spot price rose to a 17-year high $106 (U.S. dollars) per pound of U3O8 (triuranium oxide, commonly known as “yellowcake”), which is expected to increase uranium exploration, mine development, and uranium production domestically and worldwide.

Fact Sheet↗
USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗
USGS2024

The 3D Elevation Program—Supporting Ohio's economy

Introduction High-quality elevation data are proving to be a resource of great economic value in dealing with many important issues in Ohio. Current and accurate high-resolution elevation data support flood risk management, water quantity and quality assessment, precision farming, conservation planning, impervious-surface modeling, forest and other natural resources management, abandoned mine and geologic hazard assessment, karst mapping, and siting of wellhead pads for horizontal drilling. These data also support coastal zone management, traffic safety and preliminary engineering site-selection studies for transportation infrastructure, solar potential and other renewable energy planning, aviation safety, and identification of features of interest or concern such as archaeological sites and orphan oil and gas wells. Critical applications that meet the State’s management needs depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features.

Ohio↗
USGS2024

U.S. Geological Survey Karst Interest Group proceedings, Nashville, Tennessee, October 22-24, 2024

Karst hydrogeologic systems represent challenging and unique conditions to scientists studying groundwater flow and contaminant transport. Karst terrains are characterized by distinct and beautiful landscapes, caverns, and springs, and many of the exceptional karst areas are designated as national or state parks. The range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than any other aquifer type. The U.S. Geological Survey (USGS) Karst Interest Group (KIG), formed in 2000, is a loosely knit, grassroots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst aquifers. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. To accomplish its mission, the KIG has organized a series of workshops. To date (2024), nine KIG workshops, including the workshop documented in this report, have been held. The abstracts and extended abstracts provide a snapshot in time of past and current karst related studies. The USGS Water Availability and Use Science Program funded the workshop and proceedings. The planning committee for the ninth workshop includes Thomas D. Byl (USGS and Tennessee State University), Allan K. Clark (USGS), Laura M. DeMott (USGS), Eve L. Kuniansky (USGS, Emeritus), Benjamin V. Miller (USGS), and Lawrence E. Spangler (USGS, Emeritus). The workshop proceedings are edited by Eve L. Kuniansky and Lawrence E. Spangler. The field trip guide was produced by Benjamin V. Miller and Brian Ham (Tennessee Department of Environment and Conservation) and included in the proceedings from the KIG’s 2021 virtual workshop to be used on the optional field trip held on Thursday, October 24, 2024.

Open-File Report↗
USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2023

Connecting flood-related fluvial erosion and deposition with vulnerable downstream road-stream crossings

Fluvial erosion is increasingly responsible for infrastructure and building damages associated with floods as the intensity of extreme rainfalls hit rural and urban rivers in a variety of climate settings across the United States. Extreme floods in 2016 and 2018 caused widespread culvert blockages and road failures, including extensive damage along steep tributaries and ravines in the Marengo River, Wisconsin, watershed during 2016 and 2018. A study conducted by the U.S. Geological Survey (USGS), Wisconsin Wetlands Association (WWA), Ashland County, and the Northwest Wisconsin Regional Planning Commission (NWRPC) investigated the special concern of fluvial erosion hazards (FEHs) associated with gullying, streamside landslides, and the loss of wetland storage in headwaters. In 2019, a pilot study was begun to map and classify ephemeral and perennial streams and wetlands in terms of their sensitivity to FEHs. This study combined data from field-based rapid geomorphic assessments (RGAs) coupled with a stream network-wide geographic information system (GIS) approach for mapping stream segments, referred to as fluvial process zones (FPZ), sensitive to erosion, deposition, and channel change. The GIS approach used nationally available 10-meter (m) resolution topology and an extended stream network to map FPZs based on Strahler stream order, stream power, channel slope, presence of adjacent steep valley sides and headwater flats, and adjacent landform setting. Bankfull channel widths derived from RGA-based hydraulic geometry curves combined with drainage areas, an estimate of bankfull flow, and channel slope were used to calculate specific stream power for the FPZs. Lastly, the FPZs were characterized by their location within three major landform settings that affect erosion potential. The resulting vulnerability maps provided a screening framework to identify FPZs that are sensitive to incision, gullying and mass wasting along steep headwater ephemeral channels, as well as downstream perennial channels that have the potential for valley-side landslides, coarse sediment deposition, and channel change. Lastly, each FPZ was characterized in terms of hydrologic alteration associated with ditching. The vulnerability mapping products and rankings of sensitivity of FPZs will ultimately be used by Ashland County and their collaborators to prioritize natural flood management projects that mitigate FEHs, restore hydrology, and reconnect channels with adjacent wetlands and floodplains.

Wisconsin↗
USGS2025

Bayesian mapping of regionally grouped, sparse, univariate earth science data

Some earth science data are naturally grouped by region, and it is often desirable to map these data by region. However, if there are only a few samples within each region, then the map should be smoothed in an appropriate way to mitigate the problems that arise from having only a few samples. A smoothing algorithm based on a Bayesian hierarchical model is developed and presented in this report. This algorithm has several features that make it especially suitable for mapping earth science data: it can account for measurements that are censored, it can process multiple datasets with different measurement errors and different censoring thresholds, and it can calculate the uncertainty in any statistic that is mapped. The algorithm is demonstrated by mapping gold concentrations that are measured in streambed sediments in the Taylor Mountains quadrangle in southwestern Alaska.

Alaska↗
USGS2025

Modeling interconnected minerals markets with multicommodity supply curves: Examining the copper-cobalt-nickel system

Demand for many of the metals used in the energy transition is expected to grow rapidly. Many of these are by-products, often considered critical because their production responds weakly to prices and is instead tied to the economics of the host mineral. We present a model of prices and production for jointly produced commodities that accounts for interconnectivity between host and by-product markets at the mine level. We demonstrate this method using the copper–cobalt–nickel system, in which approximately 99% of cobalt is a by-product of copper or nickel mining. Our results show that the model more accurately captures the economic benefits of diversified mine outputs than previous approaches. Furthermore, changes in demand drivers for any two commodities produce non-linear effects on production and price. We challenge the prior best-practice assumption that cobalt cannot impact the copper or nickel markets. Recognizing the importance of both copper and cobalt for future electrification, we emphasize that incentivizing the copper industry to reduce cobalt supply risks could inadvertently undermine copper supply.

Nature Communications↗
USGS1990

Large springs in the Valley and Ridge Province in Tennessee

Approximately 2,700 miscellaneous discharge measurements for 171 large springs in 28 counties of East Tennessee, predominantly within the Valley and Ridge physiographic province, were analyzed statistically and results tabulated to summarize data useful to the Appalachian Valleys-Piedmont Regional Aquifer System Analysis study. The number of measurements at each spring ranged from 1 to 65. Information from each spring includes station numbers, latitude, longitude, spring name and where data are sufficient, an estimate of the mean discharge and the discharge exceeded by 75%, 50%, and 25% of the miscellaneous measurements. Data are referenced to locations on a 1:1,000 ,000-scale map. The highest mean spring discharge was 32.2 cu ft/sec (14,5000 gallons/min).

Tennessee↗
USGS2024

Managing basin-scale carbon sequestration: A tragedy of the commons approach

The Tragedy of the Commons is a well studied problem in the literature of ecology, economics, and environmental policy which illustrates the deleterious consequences of managing common pool resources when individual and social incentives are misaligned. In this work, we apply a simple model of carbon sequestration in a deep saline aquifer by two neighboring geologic CO 2 storage (GCS) operators to begin investigating if a Tragedy of the Commons framework applies to GCS. Specifically, we consider the pressure space as a “commons” because the injection by each firm at its own well increases the downhole injection pressure at both wells. We assume that a firm will decrease its injection rate if the downhole pressure at its well exceeds a predefined maximum (i.e., exceeds the “pressure limit”). With this assumption in place, we find that the same injection flowrates are optimal for both wells, regardless of whether they are owned by the same firm or competing firms. This suggests that GCS may not be best represented by a pure Tragedy of the Commons framework under our initial assumptions. However, there could be economic incentives or contractual obligations that may result in either or both GCS operators being unwilling to reduce their injection rates. Thus, we conclude the conference paper with a discussion of future extensions of our approach that may demonstrate closer alignment with the Tragedy of the Commons, including explicit definitions of pore-space rights, firm uncertainty regarding the parameters of the Theis equation, and the potential role of unitization.

Conference Paper↗
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

Evidence for organic and igneous carbon sources in the Mesoproterozoic Adirondack Highlands graphite district, New York

This study examined the origin of graphite mineralization in the Mesoproterozoic Adirondack Highlands. Flake graphite was mined from the historic mining district around Ticonderoga, New York, USA, from the 1850s to the 1920s. Most of the production and reserves from the district are disseminated graphite from feldspathic quartz schists (e.g., Dixon-American Graphite Co. mine, Faxon property, Flake Graphite Company). Protolith sediments of the schists were deposited in the Trans-Adirondack backarc basin at ca. 1.25 Ga and were metamorphosed to high grade at 1.18–1.15 Ga during the accretionary Shawinigan orogeny and at 1.09–1.02 Ga during the Ottawan phase of the Grenvillian orogeny. Major elements and carbon isotopes of graphite from schists (δ 13 C = −28.1‰ to −5.4‰ relative to Vienna Peedee belemnite) are consistent with derivation from organic carbon in protolith sandstones that devolatilized during metamorphism. Some small mines and prospects in the district mined graphitic calcite marbles and clinopyroxene ± scapolite skarns associated with a pegmatite suite dated by U-Pb in zircon to 1.04–1.01 Ga (e.g., Lead Hill and Crown Point deposits), and rare graphite veins are observed at several deposits. Hydrothermal flake graphite from pegmatites, skarns, and veins are distinct from graphite in schist-hosted deposits and have δ 13 C values that cluster around −7‰ to −5‰, consistent with deposition from igneous CO 2 . Covariation of carbon and oxygen stable isotopes and U-Pb zircon geochronology are not compatible with precipitation from fluids derived from country rock metapelites or marbles. Some graphite districts elsewhere show compelling evidence for remobilization and redeposition of sedimentary-derived carbon as hydrothermal graphite, sometimes as thick veins producing world-class, high-grade deposits. In the Adirondacks, there is evidence for two generations of graphite mineralization with distinct carbon sources: syngenetic graphite with an organic origin in Shawinigan quartz schists and epigenetic graphite with an igneous origin in Ottawan skarns and pegmatites.

New York↗
USGS2025

Groundwater age estimates for the Mississippi River Valley alluvial aquifer based on tracer data collected during 2018–20

This study characterized groundwater age across the Mississippi River Valley alluvial aquifer (MRVA). Groundwater samples from 69 MRVA wells and 19 wells in Tertiary units of the Mississippi embayment aquifer system (MEAS) were analyzed for sulfur hexafluoride (SF 6 ), tritium ( 3 H), helium (He), and (or) carbon-14 of dissolved inorganic carbon ( 14 C). The age distributions of 89 samples were estimated by fitting lumped parameter models to processed tracer concentrations with the U.S. Geological Survey software TracerLPM. Mean ages of MRVA groundwater samples ranged from 12 to 22,000 years, with a median of 140 years. Mean ages of MEAS groundwater samples ranged from 230 to 52,000 or more years, with a median of 13,500 years. The spatial distribution of MRVA groundwater ages was found to be influenced by depth, inflow of groundwater from deeper units, and soil saturated hydraulic conductivity. In parts of the MRVA, the spatial distribution of MRVA groundwater ages was found to be influenced by annual recharge and (or) annual groundwater pumpage.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗
USGS2026

Teach me how to pycap: A high-capacity well decision support tool using analytical solutions in Python

Regulatory agencies in humid temperate environments rely on timely evaluations of streamflow depletion and drawdown to protect aquatic ecosystems and existing water users. Numerical models offer detailed insights, but their complexity and time demands often preclude their practical use in rapid decision-making. We present pycap-dss, an open-source Python package that implements a suite of analytical solutions for estimating streamflow depletion and drawdown. The tool supports superposition of multiple wells and time-varying pumping, enabling cumulative impact assessments in situations with multiple wells and streams. The software is modular and extensible, allowing users to interchange solutions or add new analytical methods. A YAML-based configuration supports batch processing of multiple wells, and an optional AnalysisProject class facilitates integration with regulatory workflows. Rigorous unit and regression testing ensures computational reliability, and continuous integration supports ongoing development. We demonstrate deterministic examples of drawdown where multiple solutions are readily compared and streamflow depletion with multiple wells in the Central Sands region of Wisconsin. We also show the value of Monte Carlo analyses of streamflow depletion in the same Central Sands example, leveraging computational efficiency to evaluate the uncertainty of individual and cumulative streamflow depletion calculations from over 200 high-capacity wells.

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