USGS Science⌕ Search

SEARCH · USGS Science

Results for “Integrated Environmental Assessment and Management”

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

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado↗
USGS2026

Structural evolution of iron coordination in proteins across Earth’s oxygenation history

Protein metal-binding sites support essential biological functions shaped by protein fold, subunit interactions, and cofactor chemistry. Because these sites encode both biochemical function and environmental constraint, they offer a route to connect protein evolution with changes in Earth’s surface environment through time. Of particular interest is iron (Fe), the most widely used metal in biology and a cofactor central to both anaerobic and aerobic metabolism. Here, we systematically compare the immediate chemical environments of functionally essential Fe-binding sites in three-dimensional protein structures to test whether Fe coordination spheres differ across oxygen contexts. Using a curated dataset of experimentally determined structures, we identify a clear shift in the local chemistry of Fe-binding environments from anaerobic to aerobic proteins. Aerobic Fe sites are significantly more hydrophilic than anaerobic ones, and amino-acid composition analyses show reduced cysteine use in aerobic Fe-binding neighborhoods. These patterns suggest that as Earth’s surface environments became more oxygenated, proteins retained Fe as a core redox metal while reconfiguring local coordination chemistry in ways less vulnerable to oxidative damage. More broadly, this study introduces and applies the Coordination Sphere Analysis and Comparison (CSAC) workflow, an open and archived Python workflow for extracting local metal-binding environments from structure datasets, providing a framework for linking metalloprotein structure to evolutionary and geobiological transitions across Earth history.

Discover Life↗
USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗
USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗
USGS2025

Characterization and validation of tidally calibrated strains from the Alto Tiberina Near Fault Observatory Strainmeter Array (TABOO-NFO-STAR)

Six horizontal borehole tensor strainmeters (TSM1-6) installed from Fall 2021 to Spring 2022 comprise the Alto Tiberina Near Fault Observatory Strainmeter Array (STAR), providing an unprecedented opportunity to investigate seismic and aseismic deformation from hazardous high- and low-angle normal faults in Italy. Prior to use in tectonic applications, they require in-situ calibration and correction for non-tectonic signals. We tidally calibrate the instruments, characterize the calibration uncertainty, and test the results against environmental and earthquake signals originating from local to teleseismic distances. The STAR sites demonstrably deviate from assumptions common to the standard manufacturer's calibrations, including negative areal coupling at TSM3-6. While the tidally calibrated strains have ~3-56% uncertainty, the calibrated dynamic strains show interstation precision and accuracy to nanostrain levels, and static coseismic offsets in the array footprint are within uncertainty. TSM3 records a complex series of strains that may arise from dynamically triggered near-borehole fracture slip and fluid flow that does not appear to affect its sensitivity to lower strain rate deformation. Future calibration improvement may be afforded with longer stable timeseries, particularly for TSM4. Overall, our analyses demonstrate expanded geodetic capability for detecting deformation in the Alto Tiberina Near Fault Observatory.

Alto Tiberina Near Fault Observatory↗
USGS2025

Engineering and risk research—A strategic vision for extending USGS earthquake science to risk assessment and mitigation

As a part of the U.S. Geological Survey (USGS) Natural Hazards Mission Area, the Engineering and Risk Project (ERP) of the Geologic Hazards Science Center conducts civil engineering research to extend natural hazards science to risk assessments. The ERP includes predominantly, but not exclusively, staff funded by the USGS Earthquake Hazards Program. A major focus of the ERP is to deliver actionable information to decisionmakers that helps promote the safety, security, and economic well-being of the Nation. The report discusses briefly the ERP’s evolution and scope, its organization structure, describes key products and deliverables, their efforts to engage with external partners, and its strategy and vision in conducting range of engineering and risk research. The research topics outlined in this report, as developed by the ERP, may be used to facilitate project planning activities during the next 10 years and serve as a longer-term planning resource for the USGS Earthquake Hazards Program and Geologic Hazards Science Center.

Circular↗
USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS2026

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate M W estimates for small earthquakes ( M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with M W > 3 . Using the relative amplitudes method to extend these M W magnitudes to many other events, we successfully provide relative moment magnitude ( M W,rel ) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated M W,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

New Mexico, Texas↗
USGS2026

Where will the cat cross the road? Comparing camera and GPS-based models for identifying wildlife corridors

Designing effective wildlife corridors is a critical conservation challenge in fragmented landscapes. GPS-based step selection functions strongly predict dispersal corridors and connectivity, but GPS collaring can be expensive and invasive. Camera-based occupancy models are widely used for connectivity analyses but may involve trade-offs in data resolution. Despite widespread use of both approaches, few studies have directly compared them using concurrent datasets. We developed a stacked single-species, single-season occupancy model and a Circuitscape connectivity surface for mountain lions (Puma concolor) on Washington’s Olympic Peninsula, USA, and compared them with a connectivity surface from an existing integrated step selection function. Both models predicted mountain lion GPS locations well, with binned Spearman rank correlations of 1 for Circuitscape and 0.96 for the step selection function, though step selection better identified habitat use by dispersers. Connectivity predictions were moderately correlated across the landscape ( r = 0.26), but agreement was strongest in human-dominated areas most critical for corridor planning. We conclude that GPS-based approaches are advantageous when data collection is feasible and the focus is on dispersal or fine-scale movement. However, camera-based approaches may be preferable for multi-species monitoring, large spatial and temporal scales, noninvasive sampling, when resources are limited, or when fine-scale or dispersal-specific inference is not required.

Washington↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS1997

Water-quality assessment of south-central Texas: Occurrence and distribution of volatile organic compounds in surface water and ground water, 1983-94, and implications for future monitoring

The study area of the South-Central Texas study unit of the National Water-Quality Assessment Program comprises the Edwards aquifer in the San Antonio region and its catchment area. The first phase of the assessment includes evaluation of existing water-quality data for surface water and ground water, including volatile organic compounds, to determine the scope of planned monitoring. Most analyses of volatile organic compounds in surface water are from the National Pollutant Discharge Elimination System sites in San Antonio, Texas. Nine volatile organic compounds were detected at the six sites. The three compounds with the most detections at National Pollutant Discharge Elimination System sites are 1,2,4-trimethylbenzene, toluene, and xylene. Analysis of volatile organic compounds in ground water was limited to Edwards aquifer wells. Twenty-eight volatile organic compounds were detected in samples from 89 wells. The five most commonly detected compounds in samples from wells, in descending order, are tetrachloroethene, trichloroethene, bromoform, chloroform, and dibromochloromethane. Detections of volatile organic compounds in surface water and ground water within the South-Central Texas study area are limited to site-specific sources associated with development; therefore, planned monitoring for possible detections of volatile organic compounds as part of the National Water-Quality Assessment Program will emphasize areas of expanding population and development. Monitoring of volatile organic compounds is planned at National Pollutant Discharge Elimination System sites, at basic fixed surface-water sites, and in the ground-water study-unit surveys.

Texas↗
USGS2023

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper↗
USGS2024

Simulated ground-motion records for the seismic assessment of monumental masonry structures

Earthquakes are natural disasters that can cause widespread devastation and loss of life. Simulated ground-motion records can be useful in regions with limited seismic stations or a history of damaging but infrequent earthquakes. This is especially true in areas with a high concentration of heritage masonry structures, which are especially susceptible to damage, as simulated records can be crucial in predicting their seismic response and protecting these buildings from seismic damage. Despite the importance of simulated earthquakes, few studies have investigated how effective they are compared to real earthquakes when assessing the structural response of heritage buildings. To address this knowledge gap, we employ two different simulation methods of the Mw 6.2 Faial earthquake, which occurred on July 9th, 1998, in the Azores, to replicate the recorded time-series at four available stations within an epicentral distance of 150 km. The study has two objectives: first, to validate the simulated records of the 1998 Faial earthquake using alternative stochastic ground-motion simulation approaches, and second, to determine how these approaches impact the seismic assessment of historic masonry structures. To accomplish these objectives, this study uses real and simulated ground-motion datasets to conduct non-linear response history analyses of the São Francisco Church, a monumental structure in Horta that sustained damage during the Faial earthquake. The results show that both simulation approaches yield structural responses similar to the observed records.

Conference Paper↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗
USGS2024

Assessing locations susceptible to shallow landslide initiation during prolonged intense rainfall in the Lares, Utuado, and Naranjito municipalities of Puerto Rico

Hurricane Maria induced about 70 000 landslides throughout Puerto Rico, USA, including thousands each in three municipalities situated in Puerto Rico's rugged Cordillera Central range. By combining a nonlinear soil-depth model, presumed wettest-case pore pressures, and quasi-three-dimensional (3D) slope-stability analysis, we developed a landslide susceptibility map that has very good performance and continuous susceptibility zones having smooth, buffered boundaries. Our landslide susceptibility map enables assessment of potential ground-failure locations and their use as landslide sources in a companion assessment of inundation and debris-flow runout. The quasi-3D factor of safety, F 3 , showed strong inverse correlation to landslide density (high density at low F 3 ). Area under the curve (AUC) of true positive rate (TPR) versus false positive rate (FPR) indicated success of F 3 in identifying head-scarp points (AUC = 0.84) and source-area polygons (0.85 ≤ AUC ≤ 0.88). The susceptibility zones enclose specific percentages of observed landslides. Thus, zone boundaries use successive F 3 levels for increasing TPR of landslide head-scarp points, with zones bounded by F 3 at TPR = 0.75, very high; F 3 at TPR = 0.90, high; and the remainder moderate to low. The very high susceptibility zone, with 118 landslides km −2 , covered 23 % of the three municipalities. The high zone (51 landslides km −2 ) covered another 10 %.

Puerto Rice↗
USGS2025

Laboratory assessment for recovery of porcine circovirus 2 and porcine reproductive and respiratory syndrome virus using two types of commercially available hollow-fiber ultrafilters

Groundwater near swine farms is an uninvestigated reservoir for porcine reproductive and respiratory syndrome virus (PRRSV) and porcine circoviruses (PCVs). Enteric microorganisms are often collected from groundwater via dead-end ultrafiltration, but recovery of PRRSV and PCV with this method has not been assessed. We recovered PRRSV2 and PCV2 by dead-end ultrafiltration followed by polyethylene glycol (PEG) precipitation, nucleic acid extraction, and reverse-transcription quantitative real-time PCR. We also compared 2 commercial hemodialysis ultrafilters (Asahi Kasei Rexeed-25A, Nipro Elisio-25H) and compared PRRSV2 recovery in these filters to other waterborne microorganisms. On average, 8 ± 1% of PRRSV2 was recovered by dead-end ultrafiltration and PEG precipitation, compared to 25 ± 6% for adenovirus 41. Full-process recovery of bacteria in the same filters was 5–15%; Cryptosporidium parvum recovery was 42 ± 12%. PCV2 was detected in 4 of 12 replicate filters, but low stock concentrations precluded quantitative recovery estimates. Elisio-25H ultrafilters performed similarly to Rexeed-25A filters for all organisms tested and is an effective replacement for the Rexeed-25A, which is no longer available in the United States. Our recovery of PRRSV2 and PCV2 by dead-end ultrafiltration in the laboratory suggests that PRRSV2 detection limits are as low as 3–50 genomic copies/L in sample volumes of 100–1,500 L. Based on quantitative microbial risk assessment, these concentrations are relevant to PRRSV2 infection rates in the U.S. swine herd.

Journal of Veterinary Diagnostic Investigation↗
USGS2026

PFAS remediation in a bioelectrochemical system inoculated with the west branch consortium (WBC-2)

Groundwater contamination by per - and polyfluoroalkyl substances (PFAS) poses a persistent environmental and public health concern. This study evaluates a two-chambered bioelectrochemical system (BES) inoculated with the West Branch Consortium (WBC-2) for PFAS remediation. Under an applied cathodic potential of −450 mV (versus Ag/AgCl), the BES with active WBC-2 achieved >99.0% perfluorooctanesulfonic acid (PFOS) removal within 21 days in deionized water with culture medium and > 98.9% removal of PFOS, perfluorooctanoic acid (PFOA), perfluorohexanoic acid (PFHxA), and perfluorohexanesulfonic acid (PFHxS) in contaminated groundwater after 102 days. Intermediate formation (e.g., PFOA, 6:2 fluorotelomer sulfonate (6:2 FTS), perfluoropropionic acid (PFPrA), perfluorobutanoic acid (PFBA)) and background-corrected fluoride release were consistent with PFOS transformation under anaerobic reducing conditions potentially involving defluorination. Following repeated PFOS spikes (100 μg/L on Days 0, 50, and 399), PFOA, PFPrA, and PFBA accumulated over 664 days. Despite being the dominant accumulated compound, PFOA accounted for <1.8% of the total spiked PFOS mass. Minimal PFOS transformation occurred in controls without active WBC-2, highlighting the importance of microbial metabolism. Biofilm analysis revealed dense colonization of rod-shaped bacteria on carbon fiber brushes. Enrichment of Bacillus , Agrobacterium , and other low-abundance taxa suggests selective adaptation to BES and PFAS conditions. These findings highlight BES driven by electrochemically stimulated microbial activity as a promising strategy for PFAS remediation.

Journal of Water Process Engineering↗
USGS2024

Estimating domestic self-supplied water use in Rhode Island, 2014–21

Water withdrawal from private groundwater wells is often unaccounted for in water planning studies, and water from private wells can be a source of exposure to environmental contaminants. The sizes of populations that depend on private wells for domestic water use and the amounts of water that are withdrawn from these wells are generally poorly represented in data collection efforts because of the challenges of locating, metering, or gathering withdrawal information from individual property owners. To address this problem, the U.S. Geological Survey, in cooperation with the Rhode Island Water Resources Board, estimated the volume of water withdrawn from domestic self-supply wells and the populations who use them for the State of Rhode Island at a 30-meter pixel spatial resolution and one-month temporal resolution between July 2014 and June 2021. The number of people reliant on domestic self-supply wells has increased in Rhode Island over the study period; however, the statewide estimate of total water withdrawal has not statistically increased. Withdrawals from private wells are largest in the inland areas of the western part of the State, and the towns of Scituate and Charlestown have the highest estimated withdrawals. Statewide monthly withdrawals ranged from 3.987 million gallons per day in March 2018 to 7.767 million gallons per day in September 2016. The median per capita domestic water use rate was 46.0 gallons per capita per day.

Rhode Island↗