SEARCH · USGS Science

Results for “Geological Society of America Bulletin”

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

SourceResearch connectionFrom yearThrough year

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

1,666 records · Page 38Linked to original sources

USGS2024

The addition of 144Nd atomic mass to routine ICP-MS analysis as a Quick Screening Tool for Approximating Rare Earth Elements (Q-STAR) in natural waters

Rare earth elements (REEs) are a class of critical minerals, all of which can have supply chain vulnerability that impacts economic security. These elements are widely measured in environmental matrices via inductively coupled plasma mass spectrometry (ICP-MS); however, successful quantification can require time-consuming, sample-specific optimization. While a sample-by-sample approach is appropriate for targeted quantification studies, this approach is not suitable for mineral exploration efforts where rapidly screening thousands of samples for the presence of REEs is desired. Here, we demonstrated the use of a Quick Screening Tool for Approximating REEs (Q-STAR) to detect REEs in surface water and groundwater matrices, collected as part of existing environmental studies. A mass-to-charge ratio of 144 ( m / z = 144) was added to an ICP-MS method to screen for REEs in filtered water samples submitted for metals analyses to the U.S. Geological Survey (USGS) National Water Quality Laboratory. We detected the presence of REEs above a reference threshold of 1200 counts per second in 18 % of pre-selected 6626 samples. Using this screened dataset, we mapped estimated dissolved REE concentrations across the United States in relation to ecoregions and underlying geology. Data are constrained to where sample collection took place but nevertheless show estimated aqueous dissolved REE concentrations on a geographic scale that has not yet been studied. To validate Q-STAR, REEs were measured in a USGS standard reference sample, a subset of 88 archived filtered water samples, and in fresh filtered surface water samples. Our targeted analyses demonstrated a strong linear relationship between Q-STAR predicted and measured values in all archived samples for Nd (r 2 = 0.94), and light REEs (LREEs) such as lanthanum (La) (r 2 = 0.93), praseodymium (Pr) (r 2 = 0.94) and samarium (Sm) (r 2 = 0.94). Using Q-STAR screen values, nine field sites were identified and surface water samples recollected to confirm the continued presence of Nd and LREEs. Q-STAR can be used to screen an unlimited number of water samples for the presence of REEs prior to time-intensive and costly quantitative analyses and to generate large REE datasets for further investigation.

Journal of Geochemical Exploration↗

USGS2025

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa↗

USGS2024

Analysis of factors affecting plume remediation in a sole-source aquifer system, southeastern Nassau County, New York

Several plumes of dissolved, chlorinated solvents, including trichloroethylene, have been identified in a sole-source aquifer near the former Northrop Grumman Bethpage Facility and Naval Weapons Industrial Reserve Plant sites in southeastern Nassau County, New York. Past investigations have documented that the groundwater contamination originated from this industrial area and now extends to the south, in the direction of groundwater flow. The intermixed plumes are commonly referred to as the “Navy Grumman groundwater plume.” Detailed groundwater-flow modeling was needed for the New York State Department of Environmental Conservation (NYSDEC) to evaluate design options necessary for the construction, operation, optimization, maintenance, and monitoring of a groundwater extraction and treatment cleanup plan selected in a December 2019 Amended Record of Decision by the NYSDEC to comprehensively address these plumes. Consequently, the NYSDEC began a cooperative study with the U.S. Geological Survey in 2020 to better understand the local hydrogeologic framework using two independent approaches to characterize aquifer heterogeneity and update an existing regional groundwater-flow model to provide transient boundary conditions for new inset groundwater-flow models of the plume area. We developed these detailed inset models for the two independent aquifer characterizations using history-matching techniques coupled with a novel approach to risk-based management optimization of the remedial design. We also used the updated regional model to assess this optimized groundwater extraction and treatment design for potential saltwater intrusion. The ensembles of parameters resulting from history matching provided a platform with which to evaluate capture by water-supply and remedial wells using particle-tracking techniques. Using the ensemble to select a risk stance, we performed multiobjective optimization to identify various configurations of remedial pumping that are consistent with external constraints and that favor potentially competing objectives. Multiple solutions provide tradeoffs that NYSDEC can consider. In general, pumping redistribution may help to prevent further contamination migration downgradient. These and other study results are intended to support decisions for the remedial design focused on the local area encompassing the full extent of the Navy Grumman groundwater plume.

New York↗

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↗

USGS2000

Use of a ground-penetrating radar system to detect pre- and post-flood scour at selected bridge sites in New Hampshire, 1996-98

Ground-penetrating radar was used to measure the depth and extent of existing and infilled scour holes and previous scour surfaces at seven bridges in New Hampshire from April 1996 to November 1998. Ground-penetrating-radar survey techniques initially were used by the U.S. Geological Survey to study streambed scour at 30 bridges. Sixteen of the 30 bridges were re-surveyed where floods exceeded a 2-year recurrence interval. A 300-megahertz signal was used in the ground-penetrating radar system that penetrated through depths as great as 20 feet of water and as great as 32 feet of streambed materials. Existing scour-hole dimensions, infilled thickness, previous scour surfaces, and streambed materials were detected using ground-penetrating radar. Depths to riprap materials and pier footings were identified and verified with bridge plans. Post data-collection-processing techniques were applied to assist in the interpretation of the data, and the processed data were displayed and printed as line plots. Processing included distance normalization, migration, and filtering but processing was kept to a minimum and some interference from multiple reflections was left in the record. Of the 16 post-flood bridges, 22 ground-penetrating-radar cross sections at 7 bridges were compared and presented in this report. Existing scour holes were detected during 1996 (pre-flood) data collection in nine cross sections where scour depths ranged from 1 to 3 feet. New scour holes were detected during 1998 (post-flood) data collection in four cross sections where scour depths were as great as 4 feet deep. Infilled scour holes were detected in seven cross sections, where depths of infilling ranged from less than 1 to 4 feet. Depth of infilling by means of steel rod and hammer was difficult to verify in the field because of cobble and boulder streambeds or deep water. Previous scour surfaces in streambed materials were identified in 15 cross sections and the depths to these surfaces ranged from 1 to 10 feet below the streambed. Riprap materials or pier footings were identified in all cross sections. Calculated record depths generally agree with bridge plans. Pier footings were exposed at two bridges and steel pile was exposed at one bridge. Exposures were verified by field observations.

New Hampshire↗

USGS1986

Water-quality appraisal of NASQAN stations below impoundments, eastern Tennessee

The National Stream Quality Accounting Network (NASQAN) is a network of stations at which systematic and continuing water-quality data are collected. Major objectives of this U.S. Geological Survey program are (1) to depict areal variability of streamflow and water-quality conditions nationwide on a year-by-year basis and (2) to detect long-term changes in streamflow and stream quality. Several NASQAN stations in East Tennessee are downstream from impoundments which have a significant effect on water quality. NASQAN data obtained from the Tennessee River below Watts Bar Dam and the Clinch River below Melton Hill Dam were compared to water-quality data from the basins upstream. The comparison indicates that NASQAN data obtained below impoundment may not be adequate to describe a composite picture of water quality in the accounting unit. Detention time of storage in the impoundments is believed to moderate the range of constituent values observed at the NASQAN stations. Data obtained upstream and downstream from Watts Bar Dam indicate that the water sampled at the NASQAN station comes from stratified layers of the impoundment and is not representative of an integrated sample of water from the impoundment. Values of total recoverable iron suggest that, because of adsorption to sediments in impoundments, some constituents are not accurately described by sampling below impoundments. Relations between water-quality constituents and flow at stations on the Clinch River and Tennessee River are not well defined due to regulation. Direct load computations for many constituents were therefore not possible, which diminished the utility of data from these NASQAN stations to account for quantity versus quality of the water. Load computations were only possible for ionic constituents through use of a continuous specific-conductance record as an intermediary. Compensation for the effects of discharge prior to application of the Seasonal Kendall test for trends could not be done and identification of trends in water-quality constituents caused by some process (source) change was not possible. Some water-quality trends indicated by data from the Clinch and Tennessee Rivers might reflect the decreasing trend in discharge during the 1972-82 water years. Thus the stations below Watts Bar Dam and below Melton Hill Dam do not adequately meet the NASQAN objective to detect and assess long-term changes in stream quality.

Tennessee↗

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

Surficial geology and Quaternary tectonics of the Madison Valley and fault zone, Madison, Gallatin, and Beaverhead Counties, southwest Montana

The north-northwest-striking Madison fault is approximately 95 kilometers in length, lying at the confluence of the northeastern Basin and Range province and the Yellowstone tectonic parabola. The fault zone consists primarily of west-dipping normal faults that have east-dipping antithetic faults, which create the Madison Valley graben and several northeast-trending intrabasin faults. The Madison fault and associated sections discussed herein refer to the main west-dipping, range-bounding fault along the eastern side of the valley. Detailed geologic mapping (1:12,000 scale) of the entire fault zone and fault scarp profiling (total of 102 profiles) of the Madison fault reveal greater late Quaternary paleoseismic activity towards the south, including at least three paleoevents along the southern part of the fault that postdate Pinedale glaciation. Early to middle Holocene alluvial fans have vertical surface offsets that average between 2.0 and 3.0 meters and define the characteristic single-event surface offset. Pinedale lateral moraines have vertical surface offsets as great as 12.0 meters. Late Pleistocene to Holocene multiple-event fault scarps show little evidence of beveling, suggesting short seismic recurrence intervals and potential late Pleistocene and Holocene temporal clustering. Long-term average tectonic activity rates indicate slip rates ranging from 0.18–0.6 millimeters per year. Based on a comparison of fault-scarp height versus maximum slope angle of known regression lines developed from other paleoseismic investigations, the most recent event ranges from 5–1 ka. The northern section of the fault zone is defined by multiple normal faults, which detached the hanging walls of Laramide thrust faults within the Paleozoic and Mesozoic strata. This resulted in the partitioning of extension along multiple preexisting structures and less displacement along individual normal fault strands. Structural controls on lateral propagation of individual paleoevents involve the position of lateral ramps along preexisting Laramide contractional faults. This resulted in greater displacement within the larger basement-cored structures along the southern section of the fault zone, where extension is accommodated by one inferred principal basement-involved normal fault. Inferred east-northeast trending, intrabasin, normal faults within the southern half of the fault zone have no late Pleistocene displacement.

Idaho, Montana, Wyoming↗

USGS2026

Final project peport for “Mapping riparian vegetation response to climate change on the San Carlos Apache Reservation and Upper Gila River watershed to inform restoration priorities: 1935 to present – Phase 2: Focus on tamarisk vegetation”

Riparian ecosystems play a critical role in supporting wildlife habitats, maintaining water quality, and sustaining ecological resilience in arid landscapes. In the Upper Gila River level-4 Hydrologic Unit Code (HUC-4; Identification Number – 1504) watershed of Arizona and New Mexico, and across areas of the San Carlos Apache Tribe of the San Carlos Apache Reservation (hereafter, Tribe/Tribal – entity; or Reservation - place), riparian ecosystems have been substantially altered by the widespread expansion of tamarisk ( Tamarix spp.), an invasive and non-native riparian species. Tamarisk has high water use and increased flammability, making it a growing concern as droughts intensify across the southwestern United States (U.S.). Recent research indicates that tamarisk is increasingly stressed under prolonged drought conditions, which can elevate wildfire risk and further degrade riparian habitat. In response, the U.S. Geological Survey (USGS) Western Geographic Science Center, in collaboration with the Tribe, developed remote sensing–based tools to map riparian vegetation composition and monitor vegetation condition over time. These tools enable accurate identification of tamarisk extent, detection of vegetation stress, and comparison with native species such as cottonwood and willow. This information directly supports restoration and management actions, including targeted tamarisk removal, protection of endangered species habitat, and prioritization of areas for native vegetation recovery. The tools also provide insight into how riparian vegetation responds to changing climate and hydrologic conditions, strengthening long-term planning for riparian forest management. Findings from this work demonstrate that riparian vegetation responses vary across river systems, indicating differences in plant composition and hydro-climatic conditions. Tamarisk on both the Gila and San Carlos Rivers generally exhibit greater declines in greenness in response to higher temperature as well as lower precipitation and river flow, while cottonwood ( Populus fremontii ) and willow (Salix spp.) are particularly sensitive to changes in discharge along the San Carlos River and respond more directly to both high and low flows, with more moderate responses to precipitation and temperature. By providing actionable, science-based information, this research supports climate adaptation planning, wildfire risk reduction, and improved riparian ecosystem health, with benefits that extend beyond the Reservation to regional water and wildlife conservation efforts.

Arizona, New Mexico↗

USGS2025

Simulated hydrologic responses to proposed wastewater-returnflow scenarios in Falmouth, Massachusetts

The Cape Cod aquifer is the sole source of drinking water for communities on Cape Cod, Massachusetts, including the Town of Falmouth, where the aquifer is currently threatened by contamination from septic-system-derived nitrogen. To address this problem, the Town is proposing to sewer areas of Falmouth, treat the wastewater at the Town’s Main Wastewater Treatment Facility (a nitrogen removing/tertiary treatment facility), and discharge the treated wastewater to an ocean outfall pipe in Nantucket Sound. The U.S. Geological Survey, in cooperation with the Town of Falmouth, updated a three-dimensional steady-state groundwater flow model to represent current (defined as 2019–23) average hydrologic conditions and to simulate the long-term average freshwater hydrologic response to two wastewater-return-flow scenarios. Scenario 1 involves the sewering of all properties south of Route 28 in Falmouth, which approximates the Town’s possible sewer expansion over the next 20–30 years. Scenario 2 involves sewering of all properties in Falmouth to demonstrate the maximum potential effect of sewering on the aquifer. Overall, the simulated hydrologic response of water-table altitudes and streamflow in both scenarios was relatively small compared to fluctuations from natural recharge. In scenario 1, the water-table altitude decreased by about 0.1 feet south of Route 28, where the conversion to municipal sewers removed wastewater-return flow from onsite septic systems. The water-table altitude decreased by about 0.1–0.2 feet over a larger area in Falmouth under town-wide sewering in scenario 2. The greatest decrease in water-table altitude in both scenarios occurred near the Main Wastewater Treatment Facility, with a decrease of about 1.1 feet in scenario 1 and about 1.3 feet in scenario 2. Simulated decreases in streamflow also were estimated for six selected streams in Falmouth and Mashpee. In both scenarios, the largest simulated decreases in streamflow were at the Coonamessett River, which is the closest stream to the Main Wastewater Treatment Facility. In scenario 1, the average annual decrease in flow at the Coonamessett River was 0.1 cubic feet per second, a 1.1 percent decrease from current (2019–23) conditions. In scenario 2, streamflow at the Coonamessett River decreased by 0.6 cubic feet per second, a 5.4 percent decrease from current (2019–23) conditions.

Massachusetts↗

USGS2024

Thick- and thin-skinned contractional styles and the tectonic evolution of the northern Sangre de Cristo Mountains, Colorado, USA

The Sangre de Cristo Mountains of southern Colorado and northern New Mexico, USA, contain an unusual combination of thick- and thin-skinned contractional structures involving both basement and cover rocks in the Laramide Rocky Mountain foreland. These structures are truncated by down-faulted extensional basins to the east and west. Together with synorogenic sediments, these structures preserve a record of the rise of the Ancestral Rocky Mountains, the Laramide orogeny, and Rio Grande rifting. Laramide structures within the mountains provide clues to processes that link the three events and to necessary conditions for thin-skinned and thick-skinned contractional structures to form together in continental interiors. To examine the full variety of structural styles, a portion of the northern Sangre de Cristo fold- and-thrust belt in Colorado was described and interpreted using geologic maps and structural cross-sections. Stratigraphic relations of the Ancestral Rocky Mountain highlands and basin fill were reconstructed from existing maps. These relations allow identification of faults inherited from the Ancestral Rocky Mountains, differentiation of thrust sheets, and in some cases, estimation of the magnitude of displacement. To examine relations between Laramide thrusts and Rio Grande rifting, kinematic data were collected from a thrust fault adjacent to rift faults. Three thrust fault styles were recognized: thin-skinned basement, thin-skinned cover rocks, and thick-skinned basement. Thin-skinned thrusts arising from a hinterland beneath the present San Luis Valley carried sheets of Proterozoic basement rocks northeast over a Laramide foreland. These basement thrusts are interpreted to be faults of the Ancestral Rocky Mountains that reactivated during the Laramide orogeny. The Laramide foreland consists of thin-skinned thrusts and folds in sedimentary cover rocks as young as 49 Ma. Both thin-skinned thrusts in basement and cover rocks are bounded by thick-skinned basement thrusts that moved intermittently throughout the Laramide orogeny. We infer that thin-skinned thrusts form in continental interiors where deformation is focused in weak strata of thick basin fill and in fluid-reaction weakened preexisting faults in basement rocks. Both conditions are met in the Sangre de Cristo Mountains. Basement thrusts adjacent to the San Luis Valley contain evidence of plastic contractional microstructures overprinted by extensional microstructures that may record the transition from Laramide contraction to Rio Grande extension of the crust.

Colorado↗

USGS2024

Investigation of land cover within wetland complexes at Dixie Meadows, Churchill County, Nevada, from October 2015 to January 2022

The U.S. Geological Survey investigated land cover at subannual time steps within six wetland areas in Dixie Valley, Churchill County, Nevada, from October 2015 to January 2022. As requested by the U.S. Fish and Wildlife Service, we used aerial photography and satellite remote sensing data to map surface water and other land cover types within the wetland complexes. We identified five land cover classes using the green normalized difference vegetation index (gNDVI) and its inverse relationship to the normalized difference water index (NDWI) within three U.S. Department of Agriculture National Agriculture Imagery Program aerial images (acquired in 2015, 2017, and 2019) and 110 European Space Agency Sentinel-2 satellite images (acquired 2015–2022). The relative wetness of soil conditions within each land cover class is estimated by comparison to previously published observations of relative conductivity measured by 79 field-based sensors within the wetlands from 2019 to 2021. We mapped the areal coverage of the five land cover classes for approximately 385 acres (1,559,000 square meters [m²]) comprising six individual wetland complexes as well as a larger 1,298- acre (5,254,000-m 2 ) area of interest inclusive of the wetland complexes and adjacent landscape. Land cover of open water (Class 5) primarily within ponds at one of the wetland complexes comprised 8,333 m 2 , on average, of the wetland complexes. Land cover of mixed shallow surface water, saturated soil, and vegetation (Class 4) comprised 111,723 m 2 on average of the wetland complexes. Land cover of dense green vegetation canopy cover (Class 3) that often (46 percent of observations) had underlying surface water or saturated soil conditions comprised 592,522 m 2 on average of the wetland complexes. The remaining areas of the wetland complexes not mapped as these three land cover types (Classes 2 and 1) had sparse vegetation or bare soil cover and commonly (greater than or equal to 67 percent of observations) had dry soil conditions. The investigation of land cover detailed in this report could inform future efforts to map land cover more precisely via higher resolution remote sensing or ground-based surveying or could be incorporated with other environmental monitoring data to characterize habitat and hydrology of the wetland complexes at Dixie Meadows.

Nevada↗

USGS2024

StreamStats—A quarter century of delivering web-based geospatial and hydrologic information to the public, and lessons learned

StreamStats is a U.S. Geological Survey (USGS) web application that provides streamflow statistics, such as the 1-percent annual exceedance probability peak flow, the mean flow, and the 7-day, 10-year low flow, to the public through a map-based user interface. These statistics are used in many ways, such as in the design of roads, bridges, and other structures; in delineation of floodplains for land-use zoning and setting of insurance rates; for regulatory purposes, such as the permitting of wastewater discharges; and for hydrologic and climate change studies. StreamStats was first developed for Massachusetts and released in 2001. The application provided users with the ability to obtain streamflow statistics computed from data collected at USGS streamgages and to obtain estimates of streamflow statistics for user-selected ungaged sites. Massachusetts StreamStats used geographic information system software and digital mapping to compute drainage-basin characteristics, which were then used in statistical models to estimate streamflow statistics for the user-selected sites. The statistical models were in the form of equations that were developed through a process known as regression analysis. StreamStats was the first known web application with the ability to do interactive geoprocessing. The utility of Massachusetts StreamStats was instantly apparent, leading the USGS to develop a version of StreamStats that could be implemented nationally. USGS State offices normally were required to develop custom regression equations and prepare local digital mapping data needed for implementing StreamStats for their States. Funding needed to complete this work usually was provided through cooperative agreements between the USGS and State agencies. In 2004, Idaho became the first to be released in the national version of StreamStats. By 2023, 44 States were fully implemented and six were undergoing implementation. StreamStats has undergone many modifications over the years to keep up with changes to the underlying software and to add functionality. Customized functionality and separate linked StreamStats applications were developed for several States. Meeting the high demand for additions and improvements to StreamStats while also adhering to budgetary constraints has, at times, been challenging. The StreamStats development team has identified numerous additional improvements that could be made to provide better performance and more functionality. The lessons learned from the experience of building and operating StreamStats for nearly 25 years could be relevant to others interested in pursuing efforts of a similar scale.

Circular↗

USGS2025

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗

USGS2026

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho↗

USGS1999

Water-quality assessment of south-central Texas — Descriptions and comparisons of nutrients, pesticides, and volatile organic compounds at three intensive fixed sites, 1996-98

Water-quality samples were collected during April 1996-April 1998 at three intensive fixed sites in the San Antonio region of the South-Central Texas study unit as part of the U.S. Geological Survey National Water-Quality Assessment Program. The sampling strategy for the intensive fixed-site assessment is centered on obtaining information about the occurrence and seasonal patterns of selected constituents including nutrients, pesticides, and volatile organic compounds. The three sites selected to determine the effects of agriculture and urbanization on surface-water quality in the study unit are Medina River at LaCoste (agriculture indicator site), Salado Creek (lower station) at San Antonio (urban indicator site), and San Antonio River near Elmendorf (integrator site). Concentrations of two nutrients, dissolved nitrite plus nitrate nitrogen and total phosphorus, were largest at the integrator site, which is downstream of municipal wastewater treatment plants. Nitrite plus nitrate nitrogen concentrations at this site often exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) for drinking water. All total phosphorus concentrations at the site exceeded the EPA recommended maximum concentration for streams not discharging directly into reservoirs. Nitrite plus nitrate nitrogen concentrations at the integrator site tended to be smaller, and total phosphorus concentrations at the urban site tended to be larger in samples collected during stormflow than during base flow. The most detections and largest concentrations of three pesticides (atrazine, diazinon, and prometon) were in samples collected at the urban site. Some pesticide concentrations at the agriculture site showed a seasonal pattern of increasing concentrations during spring, the peak application season. Four pesticides (atrazine, deethylatrazine, diazinon, and prometon) were detected in at least 38 percent of samples collected at all three sites. The concentrations of all detected pesticides that have an MCL were less than the MCL at the three sites. More volatile organic compounds (VOC) were detected at the urban indicator site than at the agriculture indicator site, mostly likely because more sources are located in urbanized areas. The most VOCs detected and the largest concentrations of two VOCs (chloroform and tetrahydrofuran) were in samples from the integrator site. More VOCs were detected in samples collected at the integrator site during stormflow than during base flow. The concentrations of all detected VOCs that have an MCL were less than the MCL at the three sites.

Texas↗

USGS2026

Color-Shaded Relief Bathymetry of Lake Chelan, Washington

In 2023, the U.S. Geological Survey (USGS) completed a 3-week field operation in Lake Chelan, Washington, the third deepest lake in the conterminous United States, to acquire sediment cores, high-resolution bathymetry data, and chirp sub-bottom profiles. The high-resolution bathymetry displayed in this map publication was collected using a Norbit Winghead i67-curved array multibeam system; raw data were cleaned and gridded (3 by 3 meters [m]) in the Teledyne CARIS software. In some locations where the slopes of Lake Chelan are near-vertical, multibeam returns were sparse, resulting in artifacts resembling holes that appear as grey dots within the color-shaded bathymetric relief. Although the 2023 bathymetric mapping of Lake Chelan resulted in similar morphometric characteristics to previous hydrographic campaigns (table 1; Kendra and Singleton, 1987; Pelletier and others, 1989), these new data provide a much more detailed understanding of the two physiographic basins that comprise Lake Chelan: the deeper (as deep as 445.8 m [fig. 1]; previously reported as 453 m) and narrower (1–2 kilometers [km] wide) Lucerne Basin to the north, and the shallower (as deep as 127.8 m [fig. 1]), and wider (2–3 km wide) Wapato Basin to the south. With this detailed imaging of Lake Chelan, we find pervasive morphological evidence of slope failure and mass-transport deposits (MTDs) in at least 89 locations within the lake and 63 subaqueous deltas (fig. 2). The deeper Lucerne Basin is dominated by steep slopes and numerous subaqueous deltas fed by steep fluvial catchments, whereas the Wapato Basin contains numerous landslide scarps, scars, and debris aprons along the more gently sloping lake margins, and abundant relict kettle structures from the Pleistocene deglaciation. Lake Chelan’s geomorphic, geophysical, and sediment data contain important records of glacial retreat and advance, regional volcanic eruptions, and earthquakes that generate ground motions severe enough to induce landslides recorded in the sediment as MTDs. The epicenter of the 1872 Washington State Earthquake (M 6.5–7.2; Bakun and others, 2002; Brocher and others, 2017; Sherrod and others, 2021) is along the Spencer Canyon Fault, near the southern end of Lake Chelan. Shaking from this earthquake was strong enough to affect local aquifers, induce liquefaction of unconsolidated sediments, produce large seiches (Brocher and others, 2018), and likely cause landslides within Lake Chelan.

Washington↗

USGS2025

Collaborative drought science planning in the Colorado River Basin

The U.S. Geological Survey (USGS) is using collaborative, interdisciplinary planning to develop data and tools needed to optimize the management of water resources and land use by resource management agencies during an ongoing, multidecadal drought in the Colorado River Basin. The USGS Actionable and Strategic Integrated Science and Technology team works to build relationships with resource management agencies and other stakeholders who can benefit from the use of USGS data and products. In 2023, the Actionable and Strategic Integrated Science and Technology team hosted a series of collaborative workshops to bring together representatives of resource management agencies and other stakeholders (any person or entity with interests in a resource or location) with USGS program managers, scientists, and multidisciplinary subject matter experts to codevelop concepts for interdisciplinary drought science and technology projects to address pressing needs related to drought in the Colorado River Basin. Workshop participants identified current and recent scientific data that could be shared through a centralized online data portal. Workshop participants also identified drought science and technology needs and developed project concepts to address those science needs. Participants categorized project concepts based on their potential to develop short-, mid-, and long-term drought science data and tools, provide for the spatial or temporal expansion of ongoing USGS science projects, and address high-priority science needs. Participants developed nine project concepts: (1) understanding shifting ecohydrologic baselines, (2) San Juan River Basin synthesis, (3) incorporating dynamic land cover into hydrologic models, (4) aridification compared to drought, (5) surface water-groundwater interactions, (6) cascading effects of drought on dust, (7) cascading effects of drought on water availability, (8) cascading effects of drought on socioeconomic factors, and (9) the value of water in the Colorado River Basin. This report provides an overview of the 2023 Codesign Workshop Series, synthesized outcomes from workshop materials and discussions, and science project concepts that emerged from the collaborative meetings that will continue to be refined into science project proposals through codevelopment processes. This report also highlights lessons learned and next steps needed to receive feedback and testing of the USGS Science Collaboration Portal, continue collaboration to develop detailed specifics and steps for short-term wins, develop interdisciplinary project proposals, and implement science planning and studies.

Arizona, California, Colorado, Nevada, New Mexico,↗