Search USGS⌕ Search

SEARCH · Search USGS

Results for “Hydrological Processes”

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

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

536 records · Page 8Linked to original sources
USGS2024

The 3D National Topography Model Call for Action—Part 1. The 3D Hydrography Program

The U.S. Geological Survey is initiating the 3D Hydrography Program (3DHP), the first systematic remapping of the Nation’s surface waters since the original 1:24,000-scale topographic mapping program was active from 1947 to 1992. Building on decades of experience maintaining the National Hydrography Dataset (NHD), the Watershed Boundary Dataset (WBD), and the NHDPlus High Resolution (NHDPlus HR), the 3DHP will completely refresh the Nation’s hydrography data and improve discovery and sharing of water-related data. The design of the 3DHP is based on the results of a study that estimated that the fully implemented program would have the potential to provide more than $1 billion in benefits to Federal, State, Tribal, Territorial, and local governments and to private and nonprofit organizations every year, in addition to myriad societal benefits. The 3DHP would directly support better decision making regarding water resources by providing more accurate, complete, and integrated information than is currently available. The 3DHP datasets will include a three-dimensional (3D) hydrography network generated from and integrated with elevation data from the 3D Elevation Program (3DEP) to better represent stream gradients and channel conditions, along with waterbodies, hydrologic units, hydrologically enhanced elevation and other surfaces, and more consistent and accurate attributes. The 3DHP datasets will inherit key attributes of the NHD, WBD, and NHDPlus HR, and they also will include new attributes and links to other data such as the U.S. Fish and Wildlife Service National Wetlands Inventory, groundwater data, and engineered hydrologic systems such as stormwater networks. The 3DHP will be designed to provide a set of open and interoperable web-based tools, maps, and data catalogs, creating a robust system for users to reference their information about water; the system elements are collectively referred to as the “infostructure.” The 3DHP and the infostructure can provide a foundational geospatial underpinning for the Internet of Water, a community-based effort to modernize tools and technologies to share water data. As proposed, the 3DHP would begin providing products and services to the public in 2024.

Circular↗
USGS2026

Emergent headwater types for watershed monitoring, research, and planning

Headwaters—that is headwater streams and the watersheds they drain—strongly influence water quality, ecosystem services, and hydrologic connectivity, yet remain poorly mapped, monitored, and understood at large spatial scales. Using available spatial data to classify headwater “types” could therefore enhance the utility of sparse monitoring data to advance research and watershed planning—and help update literature-based perceptual models of headwaters as steep, forested, and minimally disturbed watersheds. Here, we analyze headwaters across the contiguous United States to quantify how their attributes (climate, topography, land cover, soils) vary regionally and differ from their downstream watersheds. Our analysis (a) identified seven dominant headwater types with distinct combinations of biophysical and climatic attributes and (b) revealed systematic contrasts with downstream watersheds, with strongest differences in attributes representing potential flowpath depth. Our findings offer a transferable framework for classifying, modeling, and managing these critical yet understudied components of the hydrologic network.

conterminous United States↗
USGS2025

Beavers in the Tualatin River Basin, northwestern Oregon

Growing interest in beaver-assisted restoration in the Tualatin River Basin of northwestern Oregon motivated a series of studies by the U.S. Geological Survey to assess the capacity of the stream network to support beaver dams and to evaluate the effects of beaver dams and ponds on urban streams. This multichapter volume describes the data collection from 2016–17 and the findings of these studies, which were done in partnership with Clean Water Services. Chapter A documents the locations of beaver dams in the Tualatin River Basin and how many beaver dams the stream network could support with existing and improved riparian vegetation. Beaver dam capacity was estimated by modifying existing tools to account for the low gradient of many streams in the Tualatin River Basin. Chapter B describes the effects of beaver dams and ponds on hydrologic and hydraulic responses of storm flows. Hydrologic and hydraulic responses for two urban stream reaches were compared with and without beaver dams and ponds and for a range of streamflow conditions using two-dimensional hydraulic models. Chapter C characterizes the effects of beaver dams and ponds on the transport and deposition of suspended sediment. Continuous turbidity, discrete suspended-sediment samples, and streamflow measurements collected during storms and base-flow periods were used to assess: (1) suspended-sediment loads upstream and downstream from two beaver-affected reaches, and (2) seasonal and longitudinal turbidity patterns. Chapter D describes the effects of beaver dams and ponds on longitudinal, spatial, and seasonal water-quality patterns. Continuous and synoptic water-quality data were collected along urban stream reaches, and net ecosystem production was calculated for two beaver-affected reaches. The findings of these studies illustrate that the effects of beaver dams and ponds on hydrology, hydraulics, suspended-sediment transport and deposition, and water quality are dependent on the characteristics of a stream reach (for example, channel gradient, groundwater exchange, and riparian vegetation) and the characteristics of beaver dams and ponds along that reach. This information can be used to consider the implications of beaver-assisted restoration in the Tualatin River Basin and the effects of beaver dams and ponds in urban streams.

Oregon↗
USGS2025

Hydrogeologic framework and conceptual model of the Red River alluvial aquifer east of Lake Texoma, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. At present (2025), the Oklahoma Water Resources Board has not established a maximum annual yield for the Red River alluvial aquifer east of Lake Texoma. To support the evaluation and determination of a maximum annual yield, a hydrogeologic framework and conceptual groundwater-flow model were developed to assess groundwater availability in the Red River alluvial aquifer east of Lake Texoma. The scope of this hydrologic investigation is the alluvium and terrace containing the Red River alluvial aquifer in Oklahoma between Lake Texoma, the Texas State line, and the Arkansas State line, an extent referred to in this report as “the eastern part of the Red River alluvial aquifer.” Parts of the alluvium and terrace extent in Arkansas and Texas are included in some analyses to address hydrologic influences from outside the aquifer’s boundaries in Oklahoma. The eastern part of the Red River alluvial aquifer in southeastern Oklahoma consists of approximately 401,280 acres of Quaternary alluvium and terrace deposits associated with the Red River and its major tributaries. Mean annual recharge to the aquifer for the 1980–2022 study period was estimated to be 8.62 inches per year, or 17.98 percent of the mean annual precipitation over the same period (47.94 inches). This mean annual recharge rate is equivalent to an inflow of approximately 288,250 acre-feet per year for the eastern part of the Red River alluvial aquifer. Recharge estimated using the Soil-Water-Balance code accounts for 98.7 percent of the conceptual-model inflows to the eastern part of the Red River alluvial aquifer. Saturated-zone evapotranspiration accounts for 11.9 percent and net streambed seepage accounts for 87.4 percent of the outflows in the conceptual model.

Arkansas, Oklahoma, Texas↗
USGS2025

Disturbance is the primary determinant of food chain length when the top predator is constant

Food chain length (FCL) is a primary determinant of food web structure and is hypothesized to be influenced by habitat size, productivity, and disturbance. Understanding the environmental characteristics that determine food chain length can assist in understanding how food webs may be impacted due to changes in habitats and environmental characteristics. This study examines the impact of hydrologic disturbance on stream food webs when the top predator is constant. We analyzed FCL in less disturbed groundwater flashy streams and more disturbed runoff flashy streams using stable isotopes. Despite no difference in species richness or fish density, food chains in more disturbed streams had a lower FCL compared to food chains in more stable streams. Assemblage analysis showed that flow regime and drainage area significantly impacted individual species abundances. The more disturbed runoff flashy streams had higher proportions of primary consumer fish, such as the algivorous Campostoma sp. (Stonerollers), which likely drives the reduced FCL. Drainage area and land cover had non-significant relationships with FCL. Shifting community structure due to hydrologic variability likely leads to differences in diet of Micropterus dolomieu (Smallmouth Bass), and thus a difference in FCL.

Arkansas, Missouri, Oklahoma↗
USGS1996

Technique for estimating magnitude and frequency of peak flows in Maryland

A convenient and reliable technique for estimating flood magnitudes is required for effective flood-plain management and for the efficient design of bridges, culverts, embankments, and flood-protection structures. Methods are presented for estimating peak-flow magnitudes of selected frequencies, ranging from 2 to 500 years, for all nontidal drainage basins in Maryland. The methods were developed by generalized least-squares regression techniques using data from 219 gaged basins in and near Maryland. The State is divided into five hydrologic regions: the Appalachian Plateaus and Allegheny Ridges region, the Blue Ridge and Great Valley region, the Piedmont region, the Western Coastal Plain region, and the Eastern Coastal Plain region. These regions correspond to the physiographic provinces of the State, with the exceptions that (1) the Coastal Plain Province is divided into two hydrologic regions, and (2) there is no distinct hydrologic region corresponding to the Valley and Ridge Province as it is divided into its constituent Allegheny Ridges and Great Valley subdivisions. Sets of equations for calculating peak discharges based on physical basin characteristics are provided for each of the regions. Based on the peak-flow equations, methods for estimating peak flows are presented for ungaged and gaged streams in Maryland. The methods and equations are supported by generalized least-squares analysis of basin and flood-frequency characteristics data from 219 drainage basins in and near Maryland. estimates for each of the five regions are calculated using combinations of the fol-lowing basin characteristics: drainage area, forest cover, basin relief, carbonate rock coverage, storage, and runoff-curve number. Drainage area contributes to the estimate in all five study regions. Carbonate rock coverage is used only in the Blue Ridge and Great Valley region. Storage and runoff-curve number are used solely in the Eastern Coastal Plain region. All other basin characteristics are used in two or more regions. Standard errors of estimate for the regression equations range from 19 to 31 percent in the Appalachian Plateaus and Allegheny Ridges region, 34 to 47 percent in the Blue Ridge and Great Valley region, 33 to 48 percent in the Piedmont region, 45 to 64 percent in the Western Coastal Plain region, and 36 to 42 percent in the Eastern Coastal Plain region.

Maryland↗
USGS2024

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗
USGS2024

Will there be water? Climate change, housing needs, and future water demand in California

Climate change in California is expected to alter future water availability, impacting water supplies needed to support future housing growth and agriculture demand. In groundwater-dependent regions like California's Central Coast, new land-use related water demand and decreasing recharge is already stressing depleted groundwater basins. We developed a spatially explicit state-and-transition simulation model that integrates climate, land-use change, water demand, and groundwater gain-loss to examine the impact of future climate and land use change on groundwater balance and water demand in five counties along the Central Coast from 2010 to 2060. The model incorporated downscaled groundwater recharge projections based on a Warm/Wet and a Hot/Dry climate future from a spatially explicit hydrological process-based model. Two urbanization projections from a parcel-based, regional urban growth model representing 1) recent historical and 2) state-mandated housing growth projections were used as alternative spatial targets for future urban growth. Agricultural projections were based on recent historical trends from remote sensing data. Annual projected changes in groundwater balance were calculated as the difference between land-use related water demand, based on historical estimates, and climate-driven recharge plus agriculture return flows. Results indicate that future changes in climate-driven groundwater recharge, coupled with cumulative increases in agricultural water demand, result in overall declines in future groundwater balance, with a Hot/Dry future resulting in cumulative groundwater decline in all but Santa Cruz County. Cumulative declines by 2060 are especially prominent in San Luis Obispo (−2.9 to −5.1 Bm 3 ) and Monterey counties (−6.5 to −8.7 Bm 3 ), despite limited changes in agricultural water demand over the model period. These two counties show declining groundwater reserves in a Warm/Wet future as well, while San Benito and Santa Barbara County barely reach equilibrium. These results suggest future groundwater supplies may not be able to keep pace with regional demand and declining climate-driven recharge, resulting in a potential reduction in water security in the region. However, our county-scale projections showed new housing and associated water demand does not conflict with California's groundwater sustainability goals. Rather, future climate coupled with increasing agricultural groundwater demand may reduce water security in some counties, potentially limiting available groundwater supplies for new housing.

California↗
USGS2025

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

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

conterminous United States↗
USGS2025

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,↗
USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗
USGS2024

Three-dimensional geologic framework model of the Rio San Jose groundwater basin and adjacent areas, New Mexico

As part of a U.S. Geological Survey study in cooperation with the Bureau of Reclamation and the Pueblo of Acoma, New Mexico, and the Pueblo of Laguna, New Mexico, a digital three-dimensional geologic framework model was constructed for the Rio San Jose and its surface-water drainage basin in west-central New Mexico. This three-dimensional model defines the altitude, thickness, and extent of 18 geologic units for use in a regional numerical hydrologic model. The model included an undifferentiated Proterozoic basement layer, 13 consolidated Paleozoic and Mesozoic rock units, and 4 Cenozoic units. Model input data were compiled from published cross sections, well data, structure contour maps, selected geophysical data, and data derived from geologic maps and structural features in the study area. These data were used to construct faulted surfaces that represent the upper and lower subsurface geologic unit boundaries. The digital three-dimensional geologic framework model combines faults, the altitude of the tops of each geologic unit, and boundary lines depicting the subsurface extent of each geologic unit. The digital three-dimensional geologic model described in this report and the corresponding data release represents the generalized geometry of the subsurface geologic units; it reproduces with reasonable accuracy the input geologic data and is consistent with previously published subsurface conceptualizations of the region. The geologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

New Mexico↗
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↗
USGS2025

Beyond the mangroves: A global synthesis of tidal forested wetland types, drivers and future information opportunities

There is increasing awareness of the global diversity of tidal forested wetlands (TFWs) and their significance in the provision of ecosystem services. These ecosystems, including mangrove forests, tidal freshwater forested wetlands, supratidal forests and transitional forests together span multiple climatic zones, geomorphic settings, and inundation and salinity regimes. We utilise case studies across five continents to demonstrate the state of knowledge among TFWs. Intertidal mangroves are the best-defined of the TFWs thanks to decades of research on their geomorphology, hydrology and ecology across their broad distribution. Non-mangrove forest settings, however, demonstrate more diverse hydrological, biochemical and vegetation conditions. In many cases, non-mangrove forests are situated at upper intertidal or supratidal elevations, where surface waters and groundwater are subject to interactions between tides freshwater inputs. Salinity datasets show variations ranging from tidal freshwater forested wetlands and ‘low-salinity mangroves’ to mesohaline or marine salinities, often with high temporal variability. While the floristic composition of non-mangrove forests vary among biogeographic regions, locally dominant TFW species are commonly distributed beyond the tidal niche into non-tidal wetland and upland forests. This presents challenges for traditional remote sensing approaches to ecosystem mapping, which are mostly lacking for non-mangrove forests. Geomorphic approaches and developments in machine learning offer opportunities to address this.

Preprint↗
USGS2026

Continuous monitoring and temporal variability of fluorescence of dissolved organic matter in karst groundwater of the Edwards aquifer, south-central Texas, 2019–24

Widespread urbanization on the Edwards aquifer recharge zone in south-central Texas has prompted the development of new approaches and tools for evaluating the current and future status of water quality in the San Antonio and Barton Springs segments of the Edwards aquifer. The U.S. Geological Survey, in cooperation with the San Antonio Water System and the City of Austin, applied continuous monitoring of fluorescence of dissolved organic matter (fDOM) to characterize the sources and transport of dissolved organic matter and associated constituents in the Edwards aquifer. Continuous fDOM time-series data were adjusted using empirically derived temperature, turbidity, and inner-filter-effect correction factors prior to analysis and interpretation. Fully corrected fDOM time-series data from the study sites were evaluated for short- and long-term temporal variability in fDOM in the context of the hydrologic and climatic conditions that occurred during the 5-year study period to understand vulnerability of the aquifer to potential contaminants in recharge. Study results indicate that fDOM is an effective proxy for dissolved organic carbon that, in turn, is indicative of the influx of recent surface water and associated contaminants, and thus, aquifer vulnerability to those contaminants. Although other continuously monitored water-quality parameters provide insights into changes in water quality in response to varying hydrologic conditions, the quantitative relation between fDOM and dissolved organic carbon directly indicates the timing and magnitude of pulses of organic constituents derived from the land surface, whereas other water-quality parameters do not.

Texas↗
USGS2026

Widespread terrestrial ecosystem disruption at the onset of the Paleocene–Eocene Thermal Maximum

The Paleocene–Eocene Thermal Maximum (PETM, ~56 Mya) interval was marked by massive 13 C-depleted carbon emissions into the ocean/atmosphere system, manifested as a negative carbon isotope excursion (CIE) in sedimentary components, and ~5 °C global average warming. Episodes of hydrological perturbations and soil-erosion have been widely documented for the PETM but their link with vegetation- and carbon cycle changes remain poorly constrained. Here, we present organic microfossil evidence showing a strong increase in fern-dominated pioneer vegetation that replaced coniferous forests on the margin of the Norwegian Sea during the first millennia of the CIE. With the present stratigraphic constraints, the “fern spike” occurred simultaneously in terrestrial settings along the North Sea, Arctic Ocean, the US east coast and in southern Australia, indicating that pioneer vegetation persisted for several millennia following a partial collapse of previously stable terrestrial ecosystems. Both the ferns and influx of microcharcoal imply recurrent physical disturbance, including soil destabilization and erosion, potentially linked to droughts, wildfires, and strong hydrological forcing resulting from extreme climate change. Together with evidence for reworked clay minerals and ancient organic matter (kerogen), these findings show that highly disturbed terrestrial ecosystems were widespread across mid- and high-latitude regions globally. Carbon cycle model simulations suggest that a substantial loss of standing and buried biomass, along with oxidation of soil organic matter, acted as important positive feedbacks during the onset of the CIE. Additionally, enhanced kerogen weathering likely contributed as another major positive feedback throughout both the onset and main phase of the CIE.

Proceedings of the National Academy of Sciences↗