SEARCH · USGS Science

Results for “Water Resources Research”

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

USGS2025

Flood-Inundation Maps of the Current and Jacks Fork Rivers including the Ozark National Scenic Riverways, Southeast Missouri, 2023

Digital flood-inundation maps for a 131.8-mile reach of the Current River and a 44.6-mile reach of the Jacks Fork River, in southeast Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the Ozark Foothills Regional Planning Commission and the South Central Ozark Council of Governments. The maps also encompass the 134 miles of the Current and Jacks Fork Rivers within the Ozark National Scenic Riverways, which is the first national park area to protect a river system. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (hereafter referred to as “stages”) at eight reference USGS streamgages—five on the Current River (USGS station numbers 07064440, 07064533, 07066510, 07067000, and 07068000) and three on the Jacks Fork River (USGS station numbers 07065200, 07065495, and 07066000). Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at four of these sites (USGS station numbers 07067000, 07068000, 07065495, and 07066000). Flood profiles were computed for seven of the eight map reaches by means of two-dimensional hydraulic models and the remaining reach by a one-dimensional hydraulic model. The models were calibrated by using stage-streamflow relations or streamflow measurements at the USGS streamgages and from high-flow stage measurements from water-level loggers distributed throughout the reaches. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “action stage” or near bankfull, to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root mean square error) to delineate the area flooded at each water level and the associated water depths. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that is critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Missouri↗

USGS2026

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri↗

USGS2025

Seamless Geologic Map Database for the Intermountain West, United States: A foundational dataset for mineral systems analysis

The Intermountain West has a complex geologic history, resulting in the formation of a diverse array of mineral deposits. Effective mineral exploration requires understanding the spatial and temporal relationships among geologic processes and events, a key focus of the mineral systems approach to exploration. This paper presents the Intermountain West Seamless Geologic Map Database, a unified dataset designed to support mineral exploration. Integrating geologic provinces, structural settings, and hydrothermal alteration, the database leverages the Seamless Integrated Geologic Mapping (SIGMa) extension to the USGS Geologic Map Schema (GeMS) to standardize geologic data from varied sources. SIGMa's hierarchical stratigraphic organization and feature-level metadata enhance data interoperability and reusability, enabling seamless query, analysis, and visualization of lithology, structural features, mineral deposits, geochronology, hydrothermal alteration. volcanic activity, and By providing a regionally consistent and dynamically evolving geologic map, this database provides a foundational framework for mineral exploration and geologic research. It also allows for an efficient workflow that expedites the publication of integrated geologic map databases.

Conference Paper↗

USGS2024

Formation of vertical columnar seismic structures and seafloor depressions by groundwater discharge in the drowned Miami Terrace platform and overlying deep-water carbonates, southeastern Florida

The presence of vertical cross-formational fluid migration passageways within sedimentary basins can profoundly impact aquifer and reservoir fluid-flow and their identification is fundamental to informing management of subsurface fluid resources (groundwater, oil, gas). In an onshore and offshore southeastern part of Florida, 2D/3D seismic-reflection and bathymetry data document ∼153 vertical columnar structures composed of reflection disruptions up to 790 m in the height and averaging 360 m in diameter, and ∼219 subcircular to circular seafloor depressions up to 1334 m wide. Our study focuses on these features found within the offshore shallow-marine carbonate Miami Terrace platform, which drowned approximately at the end of the middle Miocene, and within overlying Plio-Quaternary deep-water carbonate slope and drift deposits. Most columnar structures are rooted in stratiform aquifers of the Miami Terrace platform and associated with faults or fault intersections produced by Eocene and circa late Miocene tectonics. The columns commonly terminate within the platform or as subcircular depressions along an amalgamated karstic and drowning unconformity at the platform top. The columns typically stretch upwards from a zone of deep karst cavity collapse through the Miami Terrace platform with upward decreasing sag on internal reflections. Following drowning and Plio-Quaternary partial burial of the Miami Terrace platform by deep-water deposits, the subcircular depressions and faults along the platform top were points of origin for a second phase of column growth upward into the deep-water deposits. The continuation of deep platform cavity collapse and column evolution produced pockmarks along paleo-seafloors within the deep-water deposits and at the present-day sea floor. The Plio-Quaternary pockmarks formed at water depths too deep to suggest an origin related to meteoric karst above or near sea level, but rather their formation is suggested to be related to cyclic sea level falls that drove increased groundwater head and density gradients, and seafloor discharge of offshore freshened groundwater sourced from the underlying platform. Plausibly, mixing of freshened groundwater and seawater at the seafloor discharge sites drove dissolution of the host deep-water deposits, which together with erosion by groundwater venting and current scouring formed the pockmarks. Seaward of the Plio-Quaternary seafloor pockmarks, at the late-middle Miocene upper slope of the Miami Terrace platform and along the regional karst/drowning unconformity is a slope-parallel band of ∼189 densely distributed subcircular seafloor depressions with diameters up to 1334 m at water depths up to ∼660 m. It is plausible that along the upper slope, faults and fractures produced by gravity-driven slope instability and possibly tectonics formed a dense network of fluid passageways that promoted upward artesian freshened groundwater flow to sites of discharge where mixing with seawater generated limestone dissolution and the depressions. But tectonic uplift may have forced emersion and initial meteoric sinkhole formation circa late Miocene with later enhancement by freshened groundwater discharge and bottom current erosion.

Florida↗

USGS2024

Geologic map and structure sections along the southern part of the Bartlett Springs Fault Zone and adjacent areas from Cache Creek to Lake Berryessa, northern Coast Ranges, California

Introduction Located in the Coast Ranges of northern California, the Bartlett Springs Fault Zone is the easternmost fault in the San Andreas Fault system in northern California. The fault is a right-lateral, strike-slip structure considered capable of producing an earthquake of moment magnitude 7. The purpose of this mapping is to better characterize the geology and earthquake hazards associated with the southern part of the Bartlett Springs Fault Zone and to help identify any evidence of active uplift on the faults bounding the Coast Ranges. Although the area immediately surrounding the Bartlett Springs Fault Zone is sparsely populated, its southern segment presents a potential seismic hazard to northern California communities as far away as the San Francisco Bay region and Sacramento. There are also nearby water resources, mineral resources, and public lands used for public recreation. The Coast Ranges of northern California are a series of northwest-southeast-oriented mountain ranges and valleys located north of the San Francisco Bay region, between the Pacific Ocean to the west and the Sacramento Valley to the east. The region has rugged terrain, high mountain peaks that reach more than 2,400 meters above sea level, isolated and narrow valley bottoms on which most human settlements are located, and large drainage systems that tend to follow the northwest-southeast-oriented topographic grain. The physiographic character of the region is shaped by its bedrock geology, deformational history, and active faulting. The basement rocks of the northern Coast Ranges consist of the Franciscan Complex and the Great Valley complex, the latter of which consists of two informal units, the Coast Range ophiolite and the Great Valley sequence. The Franciscan Complex and the Great Valley complex are in structural contact along the Coast Range Fault, a regional-scale structure and fundamental crustal boundary. The Franciscan Complex and the Great Valley complex are superposed by active, northwest-southeast-striking strike-slip faults that are associated with seismicity swarms. These active strike-slip faults can produce moderate to large earthquakes that have moment magnitudes of 7–8. In places, these active structures bound large ranges and valleys, suggesting that much of the modern topographic expression is the result of active deformation processes. This report contains new 1:24,000-scale geologic mapping along the southern part of the Bartlett Springs Fault Zone between Clear Lake and Lake Berryessa. The map area spans 738 square kilometers in northern Napa County, southern Lake County, and parts of Yolo and Colusa Counties. The south and east borders of the map are 90 kilometers north of San Francisco and 70 kilometers west of Sacramento, respectively. The map area is within the Knoxville mining district, which has a history of mercury and gold mining dating back to the mid-19th century. The two main towns in the region, Lower Lake and Clearlake, California, are west-northwest of the map area. Approximately 71,000 people live in the cities and rural communities located within a 40-kilometer radius of the center of the map area. The bedrock geology, cross sections, and structural data presented herein are critical for evaluating the long-term evolution of the Bartlett Springs Fault Zone. This work will supplement studies on local seismic hazards, liquefaction potential, landslide hazards, earthquake geology, natural resources, groundwater resources, engineering geology, and tectonic history by providing the background information for site-specific investigations on these subjects.

California↗

USGS2025

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia↗

USGS2025

Scoping decision-maker needs and science availability to support regional natural capital accounting in the U.S. Colorado River Basin

Natural capital accounting has the potential to yield important policy insights at multiple scales, but there remains a disconnect between regional-scale natural capital accounts and their use for informing policy. In this paper, we propose a roadmap that could lead to the creation of policy-relevant regional accounts, with steps split across an initial scoping phase and a subsequent development phase. We demonstrate the scoping steps in action with an application to the Colorado River Basin (“Basin”), a large watershed in the southwestern United States (U.S.) that has faced aridification and substantial high-profile tradeoffs around the use of its water and other natural resources. Drawing on prior U.S. Geological Survey science co-production efforts, we conducted a series of eight discussion sessions with 41 scientists and science representatives whose work is relevant to Basin water, riparian and riverine ecosystems, upland ecosystems and energy and minerals. We summarise participants' thoughts on key topics and economic linkages, their insights and questions of interest and their recommendations on existing scientific data sources and gaps. We evaluate the suitability of the available data for construction of System of Environmental-Economic Accounting (SEEA) Central Framework and SEEA Ecosystem Accounting accounts, including those for land, water, forests, energy and minerals and ecosystems (covering extent, condition and ecosystem services). We present a series of lessons learned during the scoping phase, as well as lessons that could be relevant for future practitioners engaging in the development phase. The information can help guide the development of timely and relevant regional-scale environmental-economic accounts in the U.S. and beyond.

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

USGS2026

Satellite time series analysis to quantify changing climax ciénegas using a state and transition model approach

Ciénegas are rare wetlands in arid landscapes of the North American Southwest, historically providing critical ecological and hydrological functions but increasingly threatened by changing climate and land use pressures. This study quantifies changes in ciénega condition and floodplain dynamics using a state-and-transition model (STM) informed by expert knowledge and remote sensing. Key factors include woody plant encroachment, water availability, and soil aggradation. We mapped 31 ciénegas with high-resolution imagery and analyzed Landsat data (1985–2023) to assess vegetation health and moisture using the Normalized Difference Vegetation Index (NDVI) and Normalized Difference Infrared Index (NDII). Results show substantial interannual variability in phenology, water stress, and soil moisture, with regional drying and elevation strongly influencing ciénega resilience. We classified ciénegas into three functional states—healthy, desiccated, and dormant—and mapped their 2023 condition. Trend analyses indicate most ciénegas exhibit greening despite drought, though localized variability underscores the need for site-specific management. None are in a stable climax (reference) state; rather, they transition among states in response to external drivers. Increasing woody plant cover and surface drying, likely linked to declining regional water tables, favor deep-rooted species over wetland grasses—a pattern mirrored in adjacent control plots. Spatially explicit analysis revealed intra-ciénega variability often masked by aggregated data, highlighting the importance of high-resolution monitoring. Seasonal and long-term trends provide context for understanding ciénega dynamics, including degradation and restoration pathways. This study emphasizes the importance of groundwater conservation and demonstrates how remote sensing supports long-term monitoring. The STM framework offers a practical tool for adaptive management to sustain freshwater resources in arid environments.

Arizona, New Mexico↗

USGS2024

The advantages of electric-grade heat storage in long closed-loop wells

Though long (e.g., tens of kilometers) closed-loop geothermal wells are under consideration for the extraction of subsurface heat, these wells might also serve as an efficient energy storage mechanism for electricity generation. Using a semi-analytic model solution, the potential for electric-grade heat storage as a function of ambient temperature (e.g., corresponds to depth of loop), well length, well diameter, and flow rate is evaluated. The following simplified cases for comparison of standard closed loop operation with energy storage operation were considered: [1] constant flow rate and constant temperature (90 °C) injection to represent the standard closed-loop base-case; and [2] constant flow rate and annual cycle sinusoidal temperature (90-150 °C) to represent seasonal (summer) charging while solar resources are peak. Initial temperatures for all scenarios considered herein are a uniform 175 °C. Electric-grade heat is assumed to be delivered whenever temperatures at the extraction point exceed 90 °C. For calculation purposes only, if temperature falls below 100 °C, it is assumed that heat delivered is sub-economic, so no electricity would be produced. For all scenarios, temperatures at the extraction well asymptotically approach the flow-weighted average injection temperature, but energy storage scenarios exhibited a damped time-varying signal that diminishes in magnitude with length of the loop. The asymptotic approach depends on initial temperatures in the rock and the heat extraction rate (a function of well diameter and flowrate). This analysis demonstrates that shorter closed loops can produce more electricity over time than longer closed-loops previously proposed for electricity production over typical engineering design lifetimes (e.g., 30 years). Although only a high-temperature scenario is considered herein, rock that is initially below boiling temperature would not host a standard closed-loop resource, but injection of hot water seasonally would asymptotically heat this low-temperature system to temperatures capable of electricity production. In other words, regardless of initial temperatures, closed loops could be used to store electricity with no critical minerals in the geothermal battery.

Conference Paper↗

USGS2025

Estimating the probability of export restrictions to inform mineral criticality

To assess risks associated with advanced technologies’ supply chain disruptions, governmental agencies and others have developed mineral “criticality” assessments, with criticality described using the economic impact and probability of supply chain disruptions. Previous work developed subjective supply risk indicators to approximate this probability, typically combining several factors such as supply diversity and trading partners’ political stability, where indicator weightings can substantially impact results. This work explicitly quantifies export barrier probability using an ensemble of machine learning classifiers, with probability estimates informed by exogenous variables, including prior barrier implementation and global export dominance. Major differences in high-probability countries and commodities are observed across models, but the ensemble method highlights Indonesia, China, Tanzania, and the United States as particularly high risk. The Supplementary Data File provides export barrier probability estimates for each analyzed country-commodity pair, enabling a direct, quantitative, objective contribution to assessing mineral criticality, enhancing risk identification and prioritization for policymakers.

Resources, Conservation, and Recycling↗

USGS2025

Colored shaded-relief bathymetry and acoustic backscatter of Lake Sammamish, Washington

Evidence of strong earthquakes (such as underwater landslides and associated deposits) may be recorded within the lacustrine sediments of Pacific Northwest lakes. The floor of Lake Sammamish, Wash., an approximately 11 kilometer (6.8 mile) long, 2 kilometer (1.2 mile) wide, and 35 meter (114.8 feet) deep lake located in a populated region just east of Seattle, was mapped by the U.S. Geological Survey in November of 2021 to search for evidence of past earthquakes. Mapping was conducted using a SWATHplus-M 234-kHz interferometric side-scan sonar system was pole-mounted on the U.S. Geological Survey research vessel Parke Snavely, and the system collected full-coverage bathymetric and acoustic backscatter data, which were processed to 2-meter spatial resolution. Two maps were created, a colored shaded-relief bathymetric map showing lake floor morphology (sheet 1), and an acoustic-backscatter map showing backscatter intensities (sheet 2). The results may then be utilized together to investigate past earthquake activity.

Washington↗

USGS2024

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https://ca.water.usgs.gov/bay-delta/2018-delta-wide-mapping-surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California↗

USGS2024

Geomorphic map of the Umatilla River corridor, Oregon

This map portrays the distribution of landforms along the Umatilla River in northeastern Oregon and covers a corridor 127 kilometers long from the confluence of the Umatilla River with the Columbia River upstream to Meacham Creek. The map encompasses the valley bottom and extends about 1 kilometer up the adjoining hillslopes. Map data are intended to support water quality and fisheries enhancement efforts pursuant to the First Foods, a resource-management approach that focuses on traditionally gathered foods including water, fish, big game, roots, and berries and calls attention to the reciprocity between people and the foods upon which humans depend. The Umatilla River drains about 6,300 square kilometers on the northwest slope of the Blue Mountains in northeast Oregon. Most of the drainage basin is underlain by Miocene basalt flows of the Columbia River Basalt Group. Younger, weakly lithified, late Miocene and early Pliocene gravel deposits of local origin (for example, McKay Formation) are mapped in a few places. Upland surfaces are mantled with windborne silt (loess) correlative with deposits elsewhere known as the Palouse Formation. Surfaces below an elevation of about 340 meters were inundated repeatedly by large Pleistocene glacial outburst floods, most emanating from glacial Lake Missoula in western Montana. In backflooded areas such as the lower Umatilla River valley, Missoula floods deposited extensive slack-water silt. Areas mapped as open water, active channel and tie channel, flood basin, valley bottom, and modified land constitute the geomorphic floodplain: the area subject to occasional inundation by the Umatilla River. Deposits and landforms within the floodplain are inset into Missoula flood deposits and hence postdate the 20–15-kilo-annum Missoula floods. Some floodplain deposits are no more than a few centuries old, as indicated by substantial erosion and deposition during the Umatilla River flood of February 2020, the largest since systematic measurements began in October 1903. Deposits and landforms of the floodplain are transient features within the longer-term incision of the Umatilla River into mid-Miocene flood basalts and younger gravel of the McKay Formation.

Oregon↗

USGS2026

The saline groundwater legacy of a large buried coastal paleo-estuary

Elevated groundwater salinity in coastal regions threatens the beneficial use of fresh groundwater. Coastal groundwater management typically focuses on preventing intrusion from modern sources of seawater; however, past geological processes can also leave a legacy of saline groundwater now hidden in the subsurface. Here, multiple extensive airborne electromagnetic surveys provide detailed evidence of residual salinity from a paleo-estuary filling a late Pleistocene incised valley impacting more than 10,000 km 2 that is now hidden beneath coastal Louisiana’s deltaic plain. Our results show that the three-dimensional pattern of saline groundwater beneath Louisiana mimics that of near-surface aquifers surrounding the modern Delaware Bay estuary, fingerprinting the signature of the past drowning of a large, incised valley of the Mississippi River following post-glacial sea-level rise. These findings demonstrate a new framework for understanding legacy sources of saltwater critical for managing stressed water resources along global coastlines.

Louisiana↗

USGS2025

Spatiotemporal variability of algal biomass and nitrate in Owasco and Seneca Lakes in the Finger Lakes Region, New York, in 2019

Cyanobacterial harmful algal blooms (CyanoHABs) have become increasingly common, threatening the security of water resources globally. The U.S. Geological Survey conducted high-resolution nearshore mapping surveys using boat-mounted multiparameter sondes and nitrate sensors during the summer and fall of 2019 on Owasco Lake and Seneca Lake, two lakes with documented CyanoHABs in the Finger Lakes region of New York State. Discrete sensor measurements and water-quality samples were collected at fixed points along survey routes and continuous data were generated at open-water monitoring platforms. This investigation examined whether water-quality information from nearshore surveys was representative of open-water conditions and if nearshore surveys could be used to identify areas with localized nearshore CyanoHABs and potential sources of nutrients not captured by tributary sampling. In addition to comparisons across methods, nearshore concentrations of nitrate and chlorophyll were evaluated relative to tributary outlets, cyanobacterial abundance and biovolume at discrete locations, and the locations of near-surface CyanoHABs that were designated as “confirmed with high toxins” by the New York State Department of Environmental Conservation. Nitrate and chlorophyll concentrations were comparable across methods for each lake, although concentration ranges were typically higher for nearshore mapping datasets than for nearshore discrete datasets. Nearshore surveys indicated areas of nitrate enrichment that varied temporally in both lakes. Orthophosphate was not routinely detected. Across methods, median chlorophyll concentrations were higher for the summer survey than for the fall survey in Owasco Lake. Nearshore chlorophyll concentrations varied more temporally in Owasco Lake than in Seneca Lake. Phytoplankton and cyanobacterial abundance and biovolume were about five times higher in Owasco Lake than in Seneca Lake. Cyanobacteria dominated the phytoplankton community in most samples, and Microcystis comprised the bulk of the cyanobacterial biovolume. The most abundant potential cyanotoxin-producing (specifically microcystins) genera were Microcystis , Synechococcus , Aphanocapsa , and Pseudanabaena . The cyanobacterial community composition was comparable between open-water monitoring platforms and nearshore samples. Microcystins were detected in seven survey samples from Owasco Lake, in one survey sample from Seneca Lake, and in one sample each from the open-water monitoring platforms on Owasco and Seneca Lakes that were collected about 7 days before the fall surveys. Microcystin detections were not consistently associated with high cyanobacterial cell counts or cyanotoxin-producing genera. Results from nearshore surveys were comparable to open-water monitoring platforms and discrete nearshore observations in the absence of nearshore or open-water CyanoHABs in these systems during the study. Patterns of nearshore concentrations of nitrate and chlorophyll from nearshore surveys may aid in the identification of areas with localized nitrate loading and shifts in phytoplankton abundance and community composition.

New York↗

USGS2025

Extracting data from maps: Lessons learned from the artificial intelligence for critical mineral assessment competition

The U.S. Geological Survey (USGS), Defense Advanced Projects Research Agency (DARPA), NASA Jet Propulsion Laboratory (JPL), and MITRE ran a 12-week machine learning competition aimed at accelerating development of AI tools for critical mineral assessments. The Artificial Intelligence for Critical Mineral Assessment Competition solicited innovative solutions for two challenges: 1) automated georeferencing of historical maps, and 2) automated feature extraction from historical maps. Competitors used a new dataset of historical map images to train, validate, and evaluate their models. Automated georeferencing pipelines attained a median root-mean square error of 1.1 km. Prompt-based extraction (i.e., with user input) of polygons, polylines, and points from geologic maps yielded median F1-scores of 0.77, 0.56, 0.35, respectively. Geologic maps pose numerous challenges for AI workflows because they vary significantly. However, despite its short duration, the competition yielded promising results that have since spurred further innovation in this area and led to the development of new AI tools to semi-automate key, time-consuming parts of the assessment workflow.

Applied Computing and Geosciences↗

USGS2025

Gulf Coast Basin CORE-CM initiative final report

The Bureau of Economic Geology at the University of Texas at Austin (UT-BEG) is leading the Gulf Coast Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Initiative to assess the potential to produce critical minerals (CMs), including rare earth elements (REEs) from coal, coal ash, and produced water related to oil and gas production, and related materials (alumina processing waste [red mud], heavy mineral sands, graphite, and zeolite) within the Gulf Coast Basin. This project represents the first phase in a long-term program and provides reconnaissance data that will be foundational for future work by assessing resources and suggesting plans to be conducted in future work and expanding stakeholder engagement. The project includes several tasks designed to identify, characterize, and assess several necessary aspects for development of CMs and REEs in the Gulf Coast Basin.

Gulf Coast basin↗

USGS2024

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗