SEARCH · USGS Science

Results for “Hydrology 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.

At least 1,747 records · Page 97Linked to original sources

USGS2010

Use of regression‐based models to map sensitivity of aquatic resources to atmospheric deposition in Yosemite National Park, USA

An abundance of exposed bedrock, sparse soil and vegetation, and fast hydrologic flushing rates make aquatic ecosystems in Yosemite National Park susceptible to nutrient enrichment and episodic acidification due to atmospheric deposition of nitrogen (N) and sulfur (S). In this study, multiple linear regression (MLR) models were created to estimate fall‐season nitrate and acid neutralizing capacity (ANC) in surface water in Yosemite wilderness. Input data included estimated winter N deposition, fall‐season surface‐water chemistry measurements at 52 sites, and basin characteristics derived from geographic information system layers of topography, geology, and vegetation. The MLR models accounted for 84% and 70% of the variance in surface‐water nitrate and ANC, respectively. Explanatory variables (and the sign of their coefficients) for nitrate included elevation (positive) and the abundance of neoglacial and talus deposits (positive), unvegetated terrain (positive), alluvium (negative), and riparian (negative) areas in the basins. Explanatory variables for ANC included basin area (positive) and the abundance of metamorphic rocks (positive), unvegetated terrain (negative), water (negative), and winter N deposition (negative) in the basins. The MLR equations were applied to 1407 stream reaches delineated in the National Hydrography Data Set for Yosemite, and maps of predicted surface‐water nitrate and ANC concentrations were created. Predicted surface‐water nitrate concentrations were highest in small, high‐elevation cirques, and concentrations declined downstream. Predicted ANC concentrations showed the opposite pattern, except in high‐elevation areas underlain by metamorphic rocks along the Sierran Crest, which had relatively high predicted ANC (>200 μ eq L −1 ). Maps were created to show where basin characteristics predispose aquatic resources to nutrient enrichment and acidification effects from N and S deposition. The maps can be used to help guide development of water‐quality programs designed to monitor and protect natural resources in national parks.

Yosemite Nataional Park↗

USGS2012

Permeability of rock samples from Cajon Pass, California

Hydrologic rock properties are important to the debate about the state of stress and heat flow regime of the San Andreas fault. In particular, these properties are pertinent to whether frictional heat generated by high shear stresses can be convected away by circulating ground water, thus accounting for the absence of a measurable heat flow anomaly across the fault. To help resolve this debate, we have conducted laboratory permeability measurements on intact core samples from the Cajon Pass drill hole extracted from depths between 2100 and 3500 m. These rocks were all crystalline in nature, including granodiorites, tonalites, monzogranites and gneisses. Confining and pore pressures matched the in situ pressures for each depth, and distilled water was used as the permeating medium. Permeabilities ranged from 10 −22 to 10 −19 m 2 for effective pressures between 36 and 56 MPa. In general, the permeability values decreased with depth in a manner consistent with earlier studies of rocks between 500 and 2100 m in the drill hole. Petrographic observations indicate that repeated episodes of crack healing and sealing are the mechanisms responsible for the extremely low permeability values. Secondary minerals that seal the microfractures include abundant laumontite, calcite, quartz, chlorite, and other phyllosilicates. The permeability results suggest that water circulation in the vicinity of the drill hole must be restricted to major fracture zones, with minimal exchange of water in more intact regions. This finding is consistent with geochemical evidence of little mixing of the pore waters sampled from different sections of the borehole. The results also suggest that massive water circulation through the bulk of the rock is unlikely as a mechanism for obscuring the heat flow anomaly that would be expected if shear stresses along the San Andreas fault are high.

Geophysical Research Letters↗

USGS2021

Snowpack signals in North American tree rings

Climate change has contributed to recent declines in mountain snowpack and earlier runoff, which in turn has intensified hydrological droughts in western North America. Climate model projections suggest that continued and severe snowpack reductions are expected over the 21st century, with profound consequences for ecosystems and human welfare. Yet the current understanding of trends and variability in mountain snowpack is limited by the relatively short and strongly temperature forced observational record. Motivated by the urgent need to better understand snowpack dynamics in a long-term, spatially coherent framework, here we examine snow-growth relationships in western North American tree-ring chronologies. We present an extensive network of snow-sensitive proxy data to support high space/time resolution paleosnow reconstruction, quantify and interpret the type and spatial density of snow related signals in tree-ring records, and examine the potential for regional bias in the tree-ring based reconstruction of different snow drought types (dry versus warm). Our results indicate three distinct snow-growth relationships in tree-ring chronologies: moisture-limited snow proxies that include a spring temperature signal, moisture-limited snow proxies lacking a spring temperature signal, and energy-limited snow proxies. Each proxy type is based on distinct physiological tree-growth mechanisms related to topographic and climatic site conditions, and provides unique information on mountain snowpack dynamics that can be capitalized upon within a statistical reconstruction framework. This work provides a platform and foundational background required for the accelerated production of high-quality annually-resolved snowpack reconstructions from regional to high (<12 km) spatial scales in western North America, and by extension, will support an improved understanding of the vulnerability of snowmelt-derived water resources to natural variability and future climate warming.

Environmental Research Letters↗

USGS2012

Seismicity and shear strain in the southern Great Basin of Nevada and California

This study examines the relationship between the distribution of small earthquakes (M L ≤4.3) and mechanisms of strain accumulation and relaxation in an area with long repeat times between large events, the Southern Great Basin Seismic Network (SGBSN) region. The Great Basin is a unique continental extensional province characterized by normal and strike-slip faulting, high heat flow, crust of thin to normal thickness, and high elevations. The SGBSN is operated to provide data to address suitability issues pertaining to Yucca Mountain, Nevada which is being evaluated as a potential site for a national mined geologic nuclear waste repository. Suitability issues include estimation of the probability of occurrence of future damaging earthquakes, the characterization of the mechanisms that drive hydrologic flow, and the identification of fractures (faults) that might act as flow conduits or barriers. This study attempts to explain the distribution of small earthquakes in terms of spatial variations in the shear strain field; where strain concentrates there should be a greater number of small earthquakes. Strain field models are constructed under the assumption that long term fault behavior perturbs an otherwise uniform strain field. These strain field models are then interpreted with regard to the regional tectonics and site suitability issues. Modeling results provide one possible explanation of why earthquake clusters cover regions much larger than the surface projections of any of mapped major faults; clusters in a wide band along and extending northeast of the northern half of the Furnace Creek fault may correspond to elevated shear strains along the fault and a broad cluster in the Pahranagat Shear Zone may be associated with shear strain arising from a distribution of smaller localized faults. The relatively large number of small earthquakes in the southern and eastern portions of the Nevada Test Site is consistent with the strain field models. A minimum in shear strain at Yucca Mountain is predicted by all models consistent with an almost total lack of earthquakes observed there. The region to the west of the Death Valley/Furnace Creek fault system, the portion of the study area with the most active deformation but few small earthquakes, is an area of low shear strain. A possible reason for this is that the fault configuration in the area is optimal for accommodating regional deformation via large earthquakes or creep. While there is also a relative lack of earthquakes at Yucca Mountain, this may be indicative of a lack of accumulating strain energy and thus, a lower potential for a large earthquake.

Journal of Geophysical Research↗

USGS2023

Understanding the spatiotemporal distribution of snow refugia in the rain-snow transition zone of north-central Idaho

Knowledge of snow cover distribution and disappearance dates over a wide range of scales is imperative for understanding hydrological dynamics and for habitat management of wildlife species that rely on snow cover. Identification of snow refugia, or places with relatively late snow disappearance dates (SDDs) compared to surrounding areas, is especially important as climate change alters snow cover timing and duration. The purpose of this study was to increase understanding of snow refugia in complex terrain spanning the rain-snow transition zone at fine spatial and temporal scales. To accomplish this objective, we used remote cameras to provide relatively high temporal and spatial resolution measurements on snowpack conditions. We built linear models to relate SDDs at the monitoring sites to topoclimatic and canopy cover metrics. One model to quantify SDDs included elevation, aspect, and an interaction between canopy cover and cold-air pooling potential. High-elevation, north-facing sites in cold-air pools (CAPs) had the latest SDDs, but isolated lower-elevation points also exhibited relatively late potential SDDs. Importantly, canopy cover had a much stronger effect on SDDs in CAPs than in non-CAPs, indicating that best practices in forest management for snow refugia could vary across microtopography. A second model that included in situ hydroclimate observations (December – February (DJF) temperature and March 1 snow depth) indicated that March 1 snow depth had little impact on SDD at the coldest winter temperatures, and that DJF temperatures had a stronger effect on SDD at lower snow depths, implying that the relative importance of snowfall and temperature could vary across hydroclimatic contexts in their impact on snow refugia. This new understanding of factors influencing snow refugia can guide forest management actions to increase snow retention and inform management of snow-dependent wildlife species in complex terrain.

Environmental Research Letters↗

USGS2022

Using ensemble data assimilation to estimate transient hydrologic exchange flow under highly dynamic flow conditions

Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high-frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady-state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub-daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation-based approach to estimate the sub-daily flux under highly dynamic flow conditions using multi-depth temperature observations at a 5-min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi-depth temperature data at 1- or 2-hr time intervals assuming one-dimensional flow and heat transport governing equation. By comparing estimated fluxes with model-generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation-based flux estimation method was able to capture the vertical sub-daily fluxes using multi-depth high-resolution temperature data alone, even in the presence of multi-dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1-D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

Washington↗

USGS2024

Controls on stable methane isotope signatures in northern peatlands and potential shifts in signatures under permafrost thaw scenarios

Northern peatlands are a globally significant source of methane (CH 4 ), and emissions are projected to increase due to warming and permafrost loss. Understanding the microbial mechanisms behind patterns in CH 4 production in these systems will be key to predicting annual emissions changes, with stable carbon isotopes (δ 13 C-CH 4 ) being a powerful tool for characterizing these drivers. Given that δ 13 C signatures of CH 4 are used in top-down atmospheric inversion models to partition sources, our ability to model CH 4 production pathways and associated δ 13 C-CH 4 signatures in peatland types impacted by a changing climate is critical. We sought to characterize the role of environmental conditions, including both hydrologic and vegetation patterns associated with permafrost thaw, on δ 13 C-CH 4 signatures from a diverse set of high-latitude peatlands. We measured porewater and emitted CH 4 stable isotopes, pH, and vegetation composition from five boreal-Arctic peatlands. Porewater δ 13 C-CH 4 was strongly associated with peatland type, with δ 13 C enriched values obtained from more minerotrophic fens (-61.2 ± 9.1‰) compared to permafrost-free bogs (-74.1 ± 9.4‰) and raised permafrost bogs (-81.6 ± 11.5‰). Variation in porewater δ 13 C-CH 4 was best explained by sedge cover, CH 4 concentration, and the interactive effect of peatland type and pH ( r 2 = 0.50, p < 0.001). Emitted δ 13 C-CH 4 varied greatly but was positively correlated with porewater δ 13 C-CH 4 , suggesting that porewater data can be used to predict changing emissions signatures from these systems. We calculated a weighted mean mixed atmospheric CH 4 signature for northern peatlands of -65.3 ± 7‰ and show that this signature is more sensitive to landscape drying (4 to 10 % depletion in δ 13 C) than wetting (1.5 to 5% enrichment in δ 13 C) under permafrost thaw scenarios. Our results suggest northern peatland δ 13 C-CH 4 signatures are likely to shift in the future which has important implications for source partitioning in atmospheric inversion models.

Journal of Geophysical Research: Biogeosciences↗

USGS2008

Effects of residence time on summer nitrate uptake in Mississippi River flow-regulated backwaters

Nitrate uptake may be improved in regulated floodplain rivers by increasing hydrological connectivity to backwaters. We examined summer nitrate uptake in a series of morphologically similar backwaters on the Upper Mississippi River receiving flow-regulated nitrate loads via gated culverts. Flows into individual backwaters were held constant over a summer period but varied in the summers of 2003 and 2004 to provide a range of hydraulic loads and residence times (??). The objectives were to determine optimum loading and ?? for maximum summer uptake. Higher flow adjustment led to increased loading but lower ?? and contact time for uptake. For highest flows, ?? was less than 1 day resulting in lower uptake rates (Unet, < 300 mg m-2 day-1), low uptake efficiency (U% < 20%) and a long uptake length (Snet > 4000 m). For low flows, ?? was greater than 5 days and U% approached 100%, but Unet was 200 mg m-2 day-1. Snet was < half the length of the backwaters under these conditions indicating that most of the load was assimilated in the upper reaches, leading to limited delivery to lower portions. Unet was maximal (384-629 mg m-2 day-1) for intermediate flows and ?? ranging between 1 and 1.5 days. Longer Snet (2000-4000 m) and lower U% (20-40%) reflected limitation of uptake in upper reaches by contact time, leading to transport to lower reaches for additional uptake. Uptake by ???10 000 ha of reconnected backwaters along the Upper Mississippi River (13% of the total backwater surface area) at a Unet of ???630 mg m-2 day-1 would be the equivalent of ???40% of the summer nitrate load (155 mg day-1) discharged from Lock and Dam 4. These results indicate that backwater nitrate uptake can play an important role in reducing nitrate loading to the Gulf of Mexico. Copyright ?? 2008 John Wiley & Sons, Ltd.

River Research and Applications↗

USGS2005

Effects of urban development in the Puget Lowland, Washington, on interannual streamflow patterns: Consequences for channel form and streambed disturbance

Recovery and protection of streams in urban areas depend on a comprehensive understanding of how human activities affect stream ecosystems. The hydrologic effects of urban development and the consequences for stream channel form and streambed stability were examined in 16 streams in the Puget Lowland, Washington, using three streamflow metrics that integrate storm‐scale effects of urban development over annual to decadal timescales: the fraction of time that streamflow exceeds the mean streamflow (T Qmean ), the coefficient of variation of annual maximum streamflow (CV AMF ), and the fraction of time that streamflow exceeds the 0.5‐year flood (T 0.5 ). Urban streams had low interannual variability in annual maximum streamflow and brief duration of frequent high flows, as indicated by significant correlations between road density and both CV AMF and T 0.5 . The broader distribution of streamflow indicated by T Qmean may be affected by urban development, but differences in T Qmean between streams are also likely a result of other physiographic factors. The increase in the magnitude of frequent high flows due to urban development but not their cumulative duration has important consequences for channel form and bed stability in gravel bed streams because geomorphic equilibrium depends on moderate duration streamflow (e.g., exceeded 10% of the time). Streams with low values of T Qmean and T 0.5 are narrower than expected from hydraulic geometry. Dimensionless boundary shear stress (t*) for the 0.5‐year flood was inversely related to T 0.5 among the streams, indicating frequent and extensive bed disturbance in streams with low values of T 0.5 . Although stream channels expand and the size of bed material increases in response to urban streamflow patterns, these adjustments may be insufficient to reestablish the disturbance regime in urban streams because of the differential increase in the magnitude of frequent high flows causing disturbance relative to any changes in longer duration, moderate flows that establish a stable channel.

Water Resources Research↗

USGS2000

Influence of net freshwater supply on salinity in Florida Bay

An annual water budget for Florida Bay, the large, seasonally hypersaline estuary in the Everglades National Park, was constructed using physically based models and long‐term (31 years) data on salinity, hydrology, and climate. Effects of seasonal and interannual variations of the net freshwater supply (runoff plus rainfall minus evaporation) on salinity variation within the bay were also examined. Particular attention was paid to the effects of runoff, which are the focus of ambitious plans to restore and conserve the Florida Bay ecosystem. From 1965 to 1995 the annual runoff from the Everglades into the bay was less than one tenth of the annual direct rainfall onto the bay, while estimated annual evaporation slightly exceeded annual rainfall. The average net freshwater supply to the bay over a year was thus approximately zero, and interannual variations in salinity appeared to be affected primarily by interannual fluctuations in rainfall. At the annual scale, runoff apparently had little effect on the bay as a whole during this period. On a seasonal basis, variations in rainfall, evaporation, and runoff were not in phase, and the net freshwater supply to the bay varied between positive and negative values, contributing to a strong seasonal pattern in salinity, especially in regions of the bay relatively isolated from exchanges with the Gulf of Mexico and Atlantic Ocean. Changes in runoff could have a greater effect on salinity in the bay if the seasonal patterns of rainfall and evaporation and the timing of the runoff are considered. One model was also used to simulate spatial and temporal patterns of salinity responses expected to result from changes in net freshwater supply. Simulations in which runoff was increased by a factor of 2 (but with no change in spatial pattern) indicated that increased runoff will lower salinity values in eastern Florida Bay, increase the variability of salinity in the South Region, but have little effect on salinity in the Central and West Regions.

Florida↗

USGS2006

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of 'interactive' ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d -1. 3H/3He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d-1) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.

Journal of Hydrology↗

USGS2011

Monitoring a large volume CO2 injection: Year two results from SECARB project at Denbury’s Cranfield, Mississippi, USA

The Southeast Regional Carbon Sequestration Partnership (SECARB) early project in western Mississippi has been testing monitoring tools and approaches to document storage efficiency and storage permanence under conditions of CO 2 EOR as well as downdip injection into brine. Denbury Onshore LLC is host for the study and has brought a depleted oil and gas reservoir, Cranfield Field, under CO 2 flood. Injection was started in July 2008 and has now achieved injection rates greater than 1.2 million tons/year though 23 wells, with cumulative mass injected as of August, 2010 of 2.2 million metric tons. Injection is into coarse grained fluvial deposits of the Cretaceous lower Tuscaloosa Formation in a gentle anticline at depths of 3300 m. A team of researchers from 10 institutions has collected data from five study areas, each with a different goal and different spatial and temporal scale. The Phase 2 study began at the start of injection and has been using pressure and temperature as a tool for assessing permanence mostly in the oil productive interval. Real-time read-out shows high sensitivity to distant changes in injection rate and confirms the geologic model of reservoir compartmentalization. Above-zone pressure monitoring &sim;120 m above the injection interval is used to test the sensitivity of this approach for documentation of integrity of the confining system in an area of numerous well completions as pressure increase is induced in the reservoir by more than 70 bar. Monitoring of the High Volume Injection Test (HiVIT) area includes repeat measurements of aqueous geochemistry in the injection zone. Rock-water- CO 2 interactions in the reservoir as CO 2 dissolves are minimized by mineral “armoring” by abundant chlorite cement in high permeability reservoir sandstone. Geochemical monitoring of confined freshwater aquifers at depths of 70–100 m is underway. Groundwater analysis focuses on assessment of the sensitivity of this method to detect leakage above background variability. A repeat seismic survey of the HiVIT is planned for late 2010 to assess saturation change especially in downdip brine-only areas. A study focused on feasibility of monitoring the shallow subsurface to separate leakage from normal complex surface fluxes is underway at an monitoring array installed in October 2009 to assess the interactions of recharge, soil gas, and shallow groundwater aquifers. Recent well re-entry and tracer injection will provide further information to interpret observed elevated deep-sourced methane. The Detailed Area Study (DAS) is collecting dense time-lapse data from closely-spaced three well array of an injector and two observation wells. The observation wells were completed with fiberglass casing to facilitate electrical resistance tomography (ERT) measurements, and a diverse array of instrumentation was both cemented behind casing and suspended on tubing. Injection started at the DAS December 1, 2009. We have measured pulsed neutron and resistivity via wireline, downhole and above-zone pressure, distributed temperature, and fluid chemistry including introduced pulses of perfluorocarbons, noble gases, and SF6 as tracers. Between wells, time-lapse crosswell seismic and electrical resistance tomography (ERT) are used to measure saturation change. The goals are to measure changes as fluids evolve from single phase (brine) to two phase (CO 2 –brine) in order to document linkages between pressure and sweep efficiency. A time-lapse VSP survey bridges the vertical resolution and areal coverage between cross-well and surface seismic. The repeat surveys for many tools are scheduled for September, 2010. Reservoir characterization based on cores, historic and new wireline log data, production history, hydrologic tests, fluid analysis, and a three-D seismic survey have been used in multiple numerical models to predict reservoir response in order to design effective monitoring strategies and optimize deployment. History matching of observed response to predicted response is used to interpret results and improve confidence in conceptual models and numerical approaches. Probabilistic methods have been used to assess the significant uncertainties resulting from reservoir heterogeneity.

Mississippi↗

USGS2021

A machine learning approach to modeling streamflow with sparse data in ungaged watersheds on the Wyoming Range, Wyoming, 2012–17

Scant availability of streamflow data can impede the utility of streamflow as a variable in ecological models of aquatic and terrestrial species, especially when studying small streams in watersheds that lack streamgages. Streamflow data at fine resolution and broad extent were needed by collaborators for ecological research on small streams in several ungaged watersheds of southwestern Wyoming, where streamflow data are sparse. To improve the utility of sparse streamflow data to ecological research in ungaged watersheds, we developed a machine learning approach in R for modeling spatially and temporally continuous monthly streamflow from 2012 through 2017 in three semiarid montane-steppe watersheds (with drainage areas of 26–55 square miles and mean elevations of 8,031–8,455 feet) on the Wyoming Range in the upper Green River Basin. A machine learning streamflow (MLFLOW) model was calibrated and validated with 971 discrete streamflow observations and 24 static and dynamic predictor variables derived from geospatial and time series data on climatic, physiographic, and anthropogenic characteristics affecting streamflow. The predictor variables were temporally and spatially conditioned to amplify the relation of predictor variables to monthly streamflow. The MLFLOW model had satisfactory agreement between observed and predicted streamflow (coefficient of determination [ R 2 ]=0.80, Nash-Sutcliffe efficiency [NSE]=0.79, NSE with log-transformed data [logNSE]=0.82, and percent bias [PBIAS]=0.7 percent). NSE and logNSE indicated the MLFLOW model performed equally well for high and low flows, and PBIAS indicated the MLFLOW model did not overpredict or underpredict monthly streamflow. Streamflow predictions seemed to well represent the annual hydrograph within the study area during the study period. The most important variables (statistically important in the MLFLOW model) for explaining monthly streamflow were temporally and spatially conditioned dynamic climatic variables, mostly precipitation and snow water equivalent. Importance of the static and dynamic variables did not differ substantially among the three watersheds but differed considerably among the 6 years. Monthly streamflow increased with increasing precipitation, snow water equivalent, and drainage area but decreased with increasing forest cover, elevation, evapotranspiration, and temperature. The MLFLOW model was most sensitive to selection of dynamic climatic variables. Unconditioned dynamic climatic variables alone explained 54 percent of the variance ( R 2 =0.54) in monthly streamflow, whereas adding static physiographic and anthropogenic variables only explained 12 percent more of the variance ( R 2 =0.66). Also, spatial conditioning of all variables together with temporal conditioning of dynamic variables increased the variance explained in the MLFLOW model by another 14 percent ( R 2 =0.80). The MLFLOW model also had greater sensitivity to temporal than to spatial differences in the data. For the MLFLOW model trained with observations from all watersheds and years or for models trained with observations from all except one watershed or 1 year left out sequentially, performance was better in testing on observations from each watershed than from each year separately. Also, performance was better for models fitted to fewer sites than to fewer months of observations. The greatest utility of the modeling approach is the ease of use and the speed of processing input data, running the model, and interpreting the model output, whereas the greatest limitation is the need for spatially and temporally representative streamflow observations to drive the model. Although familiarity with R is necessary, only a working knowledge of hydrology (for selecting appropriate predictor variables and evaluating the quality of streamflow observations) and a rudimentary understanding of machine learning models are needed. Therefore, this modeling approach is practicable for other scientists who work with water but who are not hydrologists.

Wyoming↗

USGS2010

Simulating the volatilization of solvents in unsaturated soils during laboratory and field infiltration experiments

This paper describes laboratory and field experiments which were conducted to study the dynamics of trichloroethylene (TCE) as it volatilized from contaminated groundwater and diffused in the presence of infiltrating water through the unsaturated soil zone to the land surface. The field experiments were conducted at the Picatinny Arsenal, which is part of the United States Geological Survey Toxic Substances Hydrology Program. In both laboratory and field settings the gas and water phase concentrations of TCE were not in equilibrium during infiltration. Gas-water mass transfer rate constants were calibrated to the experimental data using a model in which the water phase was treated as two phases: a mobile water phase and an immobile water phase. The mass transfer limitations of a volatile organic compound between the gas and liquid phases were described explicitly in the model. In the laboratory experiment the porous medium was nonsorbing, and water infiltration rates ranged from 0.076 to 0.28 cm h −1 . In the field experiment the water infiltration rate was 0.34 cm h −1 , and sorption onto the soil matrix was significant. The laboratory-calibrated gas-water mass transfer rate constant is 3.3×10 −4 h −1 for an infiltration rate of 0.076 cm h −1 and 1.4×10 −3 h −1 for an infiltration rate of 0.28 cm h −1 . The overall mass transfer rate coefficients, incorporating the contribution of mass transfer between mobile and immobile water phases and the variation of interfacial area with moisture content, range from 3×10 −4 h −1 to 1×10 −2 h −1 . A power law model relates the gas-water mass transfer rate constant to the infiltration rate and the fraction of the water phase which is mobile. It was found that the results from the laboratory experiments could not be extrapolated to the field. In order to simulate the field experiment the very slow desorption of TCE from the soil matrix was incorporated into the mathematical model. When desorption from the soil matrix was added to the model, the calibrated gas-water mass transfer rate constant is 2 orders of magnitude lower than that predicted using the power law model developed for the nonsorbing laboratory soil.

Water Resources Research↗

USGS2024

Physics to fish—Understanding the factors that create and sustain native fish habitat in the San Francisco Estuary

Executive Summary The Bureau of Reclamation (Reclamation) operates the Central Valley Project (CVP), one of the nation’s largest water projects. Reclamation has an ongoing need to improve the scientific basis for adaptive management of the CVP and, by extension, joint operations with California’s State Water Project. The U.S. Geological Survey (USGS) works cooperatively with the Bureau of Reclamation to provide scientific support for the management of Reclamation’s CVP project. Major habitat restoration efforts and a new water-diversion point are planned to benefit delta smelt ( Hypomesus transpacificus ) and other species of concern while ensuring the reliability of water supply. In addition, various flow actions and management activities have been identified as possible methods to increase populations of delta smelt and salmonid ( Oncorhynchus spp.) runs of concern. The overarching goal of this cooperative project was to provide Reclamation with the scientific information needed to evaluate the efficacy of ongoing and future adaptive management actions and to improve the scientific basis for more flexible CVP operations that would achieve water-supply reliability and fish protection. The research and monitoring described in this report comprises the period 2015–19 and focuses on management issues related to native fish species of concern, especially delta smelt. Conserving the delta smelt population while providing a reliable water supply is a primary management and policy issue in California. Our approach for this cooperative project is based on the “physics to fish” concept, the idea that high-quality habitat is generated and sustained by the interaction between physical processes and the landscape. These interactions create a template for chemical and biological processes that can change across a variety of spatial and temporal scales. Following this concept, this project (hereafter referred to as “the physics to fish project”) included monitoring and studies of water flows, sediments, water quality, and invertebrate and fish dynamics across a range of spatial and temporal scales and in regions relevant to resource managers tasked with managing water supplies and ecosystem health in the San Francisco Estuary. The intent of this approach was to document the habitat conditions, important processes, and interactions among them that create high-quality habitat for native fishes so that the likely effects of future management actions (for example, habitat restoration) can be objectively assessed at the local (site-specific), regional (within subregions of the estuary), and landscape (across the entire estuary and beyond) scales. Hydrodynamically, the upper estuary (landward of Carquinez Strait) is characterized by a fixed volume of tidally exchanged water (for example, tidal prism) that interacts with the existing channel network and bathymetry to create regions with differing hydrodynamics. Our results indicate that careful study of construction or reoperation of existing infrastructure to perform management actions can help (1) improve the accuracy of hydrodynamic models; (2) further understanding of ecological effects; and (3) enhance abilities to predict ecological outcomes. At the local scale, we developed a new concept called the Lagrangian to Eulerian (LE) ratio that can be used as a tool for understanding the importance of various hydrodynamic processes in specific channels or channel networks and for forecasting transport dynamics. Channels with LE ratios<1 in a channel network or in a dead-end slough are hydrodynamically able to develop an exchange zone between two parcels of water that may have different chemical and physical properties. In a dead-end channel, there is a landward region with long residence time (no-exchange zone) and a seaward region with short residence time (high-exchange zone) that are well mixed with seaward waters. At the transition (exchange zone) between the high and no-exchange regions, a gradient will form in water-quality constituents that differ in concentration between the landward and seaward waters. Turbidity affects fish habitat and has declined through time in the San Francisco Estuary. Average turbidity across the Sacramento–San Joaquin Delta (hereafter referred to as “the Delta”) is dependent on annual hydrology. In dry years, the region around Cache Slough (known regionally as the “Cache Slough Complex”) in the northern Delta is generally more turbid than Suisun Bay and the lower Sacramento River. When the Yolo By-Pass (known regionally as “Yolo Bypass”), a large flood bypass that runs parallel to the Sacramento River in the northern Delta, is not flooding and river flows are lower, sediment is usually transported into the Cache Slough Complex because flood tides dominate ebb tides, resulting in transport of suspended sediment from seaward areas of the upper estuary into the Cache Slough Complex. These hydrodynamic conditions also favor the formation of turbidity maximums (TMs) in the Cache Slough Complex. The TMs are areas of higher suspended-sediment concentration, providing higher-turbidity habitat favored by some fishes, including delta smelt, and they can also concentrate other constituents, including phytoplankton and organic carbon that can be important in food webs. Pelagic primary production by phytoplankton is the basis for Delta food webs supporting pelagic fishes such as delta smelt; however, phytoplankton abundance in the Delta has declined during recent decades. We examined how nutrients, hydrodynamics, and other factors affect phytoplankton blooms. Based on our results, we developed three new concepts of phytoplankton bloom formation in the Delta, each associated with a distinct set of hydrologic conditions. First, productivity cascades highlighted how local processes can contribute to phytoplankton blooms observed at the regional scale. Second, we observed phytoplankton blooms in the upper San Francisco Estuary that were associated with transport out of Yolo By-Pass (transport blooms). Third, we also documented a series of phytoplankton blooms that were in the confluence area at the landward edge of Suisun Bay. The conditions leading to creation of confluence phytoplankton blooms are not yet understood, but the confluence region connects the Cache Slough Complex with Suisun Marsh. Therefore, blooms in this area have the potential to spread to large areas of the Delta. At the landscape scale, the distribution of the invasive clams ( Potamocorbula amurensis and Corbicula fluminea , hereafter referred to as “ Corbicula ”) is driven by salinity. At smaller spatial scales, the distribution of either species is sensitive to multiple factors affecting survival and reproduction, complicating efforts to predict distribution and abundance without considering local-scale conditions across the area of interest. In the Cache Slough Complex, the area landward of the exchange zone in regions with LE ratio<1 were characterized by low abundances of Corbicula probably because recruits from seaward areas are not transported past the exchange zone and because there are no landward tributaries with adult Corbicula to provide an upstream source of recruits. Corbicula biomass was highest near or downstream from the exchange zone consistent with Corbicula grazing on phytoplankton produced in the exchange zone or transported from the no-exchange zone. The severity of Corbicula grazing could be reduced by manipulating the hydrodynamic characteristics of waterways; however, the beneficial and harmful effects on the organisms meant to benefit from increased phytoplankton production, including zooplankton and fish species of concern, should be thoroughly examined before manipulating hydrodynamic characteristics. The distribution of fishes at the landscape scale is generally driven by the position of the salinity field in the estuary. The physics to fish project compared distributions of fishes at Ryer Island, a tidal wetland in Suisun Bay and a region of variable salinity, with fish distributions at the Cache Slough Complex, a freshwater region. At Ryer Island, there was an absence of freshwater invasive species and an abundance of native species, such as Sacramento splittail ( Pogonichthys macrolepidotus ), tule perch ( Hysterocarpus traskii ), and Sacramento pikeminnow ( Ptychocheilus grandis ). The native species were almost exclusively captured in wetland and nearshore shallow-water habitat regardless of water-quality conditions. In the Cache Slough Complex, our regional scale objective was to elucidate how hydrodynamic-physical habitat interactions drive fish-community structure. Our studies showed that dendritic channel systems were better able to support native species, while intertidal habitats supported those species best able to exploit the transient character of the habitat. Habitats upstream from the exchange zone were especially important in supporting high numbers of native fishes relative to within or downstream from the exchange zone. Many of the native species were associated with tidal marsh in the no-exchange zone. More pelagic-oriented, mobile species, such as Striped Bass ( Morone saxatilis ), threadfin shad ( Dorosoma petenense ), and Sacramento pikeminnow, were more affected by water-quality conditions, such as turbidity. The physics to fish concept developed in this project provides a framework for designing individual projects and for considering the cumulative effects of multiple projects in a region, using the LE ratio as a guiding metric. The physics to fish concept may also provide a suitable framework for coordinating management actions. Tidal wetlands can function in several ways in the hydrodynamic framework. Relatively small tidal wetlands with short channel networks and with LE ratios>1 are not able to maintain a landward no-exchange zone or an exchange zone. This likely means that any contributions to pelagic food webs would be limited to resources derived from wetland vegetation, which can include dissolved and particulate organic matter (detritus) and populations of consumers that can increase in abundance based on those resources. The fate of the contributed production from these channels depends on the characteristics of the receiving waters seaward of the tidal wetland. If these channels join a large system such as Suisun Bay, then any contribution is likely to be rapidly dispersed in the larger volume; however, the channel junction might provide a focal point for consumers, such as fishes, to congregate and feed on material leaving the wetland on ebb tides before it is dispersed in the larger volume. Fishes might also access these resources by entering the wetland. The physics to fish project has established a foundation and several new concepts for understanding how habitat restoration can benefit native fish populations at the local and regional levels. Many of the ideas regarding habitat restoration and channel modifications outlined in this report could help guide management actions that could improve conditions for native fishes at little or no water cost beyond water already dedicated to other management actions. A complete list of products originating from this work is provided in appendix 1.

California↗

USGS2009

Ground and surface temperature variability for remote sensing of soil moisture in a heterogeneous landscape

At the Little River Watershed (LRW) heterogeneous landscape near Tifton Georgia US an in situ network of stations operated by the US Department of Agriculture-Agriculture Research Service-Southeast Watershed Research Lab (USDA-ARS-SEWRL) was established in 2003 for the long term study of climatic and soil biophysical processes. To develop an accurate interpolation of the in situ readings that can be used to produce distributed representations of soil moisture (SM) and energy balances at the landscape scale for remote sensing studies, we studied (1) the temporal and spatial variations of ground temperature (GT) and infra red temperature (IRT) within 30 by 30 m plots around selected network stations; (2) the relationship between the readings from the eight 30 by 30 m plots and the point reading of the network stations for the variables SM, GT and IRT; and (3) the spatial and temporal variation of GT and IRT within agriculture landuses: grass, orchard, peanuts, cotton and bare soil in the surrounding landscape. The results showed high correlations between the station readings and the adjacent 30 by 30 m plot average value for SM; high seasonal independent variation in the GT and IRT behavior among the eight 30 by 30 m plots; and site specific, in-field homogeneity in each 30 by 30 m plot. We found statistical differences in the GT and IRT between the different landuses as well as high correlations between GT and IRT regardless of the landuse. Greater standard deviations for IRT than for GT (in the range of 2-4) were found within the 30 by 30 m, suggesting that when a single point reading for this variable is selected for the validation of either remote sensing data or water-energy models, errors may occur. The results confirmed that in this landscape homogeneous 30 by 30 m plots can be used as landscape spatial units for soil moisture and ground temperature studies. Under this landscape conditions small plots can account for local expressions of environmental processes, decreasing the errors and uncertainties in remote sensing estimates caused by landscape heterogeneity.

Journal of Hydrology↗

USGS2024

Hydrologic study of green infrastructure in poorly drained urbanized soils at RecoveryPark, Detroit, Michigan, 2014–21

Uncontrolled stormwater runoff volume is a legacy stressor on sewer-system capacity that is further compromised by the effects of aging infrastructure. Green stormwater infrastructure (GSI) has been used in a variety of designs and configurations (for example, bioretention) with the goal of increasing evapotranspiration and infiltration in the local water cycle. In practice, GSIs have variable effectiveness in reducing runoff volume. An urban residential site near Detroit, Michigan, called RecoveryPark was monitored for 8 years before and after GSI construction to evaluate how effectively the GSI reduced volumes of stormwater flowing to Detroit’s Water Resource Recovery Facility through combined sewer systems. In addition to the GSI, the study site included an urban farm where salad crops were grown in hoop houses. The monitoring approach was to characterize the urban water cycle through high-frequency measurements of inflows and outflows. Datasets included meteorological data, soils and sediment characteristics, groundwater levels, flows within the combined sewer system, and soils and water chemistry with specific focus on the disposition of road salt. Although land cover within the RecoveryPark sewershed was high-density residential in the 1950s, the sewershed included only one residence within the 8.74-acre sewershed during this study. Measurements of annual precipitation at the site exceeded long-term annual averages by more than 10 inches during 3 of the 8 years of study. Potential evapotranspiration was often greater than the measured precipitation that averaged 28–34 inches per year. As compared to underlying clay-rich sediments, soils data indicated relatively permeable sediments near land surface with estimated hydraulic conductivity of 0.75 inches per hour; however, these values decreased with increasing depth. Groundwater-level data revealed increases in groundwater storage as indicated by increases in seasonal groundwater levels and development of a groundwater mound adjacent to the GSI. These increases in groundwater levels were directly adjacent to swales designed to infiltrate stormwater and only became evident after installing the GSI. Flows within the combined sewer system included rainwater, septic effluent, groundwater infiltration, leakage from water-supply lines, and release of water stored in abandoned foundations. Dry-weather flows (no rain fell within the prior 3 days) averaged 7–10 gallons per minute, which were much greater than flows estimated by septic outflow alone. A set of estimated water budgets were compiled, and results showed large discrepancies in unaccounted flows. To further examine these discrepancies, dye-tracing within the combined sewer system helped examine the sources of water by relating flow volumes to drainage area. For one of the monitoring sites within the combined sewer system along the southeast side of the study area, flows estimated by dye concentrations were more than 10 percent greater than those measured by standard methods. Through peak-flow-regression analysis, a minimum of 2.4 million gallons of water per year were infiltrated or lost to evapotranspiration because of GSI construction. After site modifications were made by excavating gravel drains to improve drainage characteristics, estimated stormwater volumes within the combined sewer system returned to near preconstruction levels. The GSI was effectively bypassed to address slow infiltration rates and standing water; the bypass all but eliminated the potential benefits of volume reduction. Late in the project, a water-quality study was added to examine the transport of road salt and associated chloride within the GSI and the combined sewer system. Continuous specific conductance was used as a surrogate for chloride concentrations to estimate that 2,790 pounds of dissolved chloride passed through the sewershed during the winter months of late 2020 through early 2021. These data were collected after GSI modification, therefore most, if not all, of the chloride was transported directly to Detroit’s Water Resource Recovery Facility via the combined sewer system. Mixing diagrams using chloride and bromide concentrations of road salt, potable water, rainwater, groundwater, and water from the combined sewer system confirmed that water within the combined sewer system is a mix of these sources. The poor condition of the combined sewer system pipes and resulting unaccounted inflows added to the challenge of accurately monitoring and identifying sources and sinks of water within the RecoveryPark sewershed. Our research results suggest that—along with clear and quantifiable objectives—the catchment and site conditions should be well-characterized before determining the GSI design. In addition, the work presented in this report provides implications and lessons learned for effectiveness and future studies of GSI in urban settings. These efforts can be improved through increased communication between stakeholders, use of high-quality soils in GSI that have suitable hydraulic characteristics, redundant data-collection networks for critical data streams, and focusing meteorological-data collection within the GSI to obtain relevant evapotranspiration data.

Michigan↗

USGS1999

Long-term experimental manipulation of winter snow regime and summer temperature in arctic and alpine tundra

Three 60 m long, 2·8 m high snowfences have been erected to study long-term effects of changing winter snow conditions on arctic and alpine tundra. This paper describes the experimental design and short-term effects. Open-top fiberglass warming chambers are placed along the experimental snow gradients and in controls areas outside the fences; each warming plot is paired with an unwarmed plot. The purpose of the experiment is to examine short- and long-term changes to the integrated physical-biological systems under simultaneous changes of winter snow regime and summer temperature, as part of the Long-Term Ecological Research network and the International Tundra Experiment. The sites were at Niwot Ridge, Colorado, a temperate high altitude site in the Colorado Rockies, and Toolik Lake, Alaska, a high-latitude site. Initial results indicate that although experimental designs are essentially identical at the arctic and alpine sites, experimental effects are different. The drift at Niwot Ridge lasts much longer than do the Toolik Lake drifts, so that the Niwot Ridge fence affects both summer and winter conditions, whereas the Toolik Lake fence affects primarily winter conditions. The temperature experiment also differs in effect between the sites. Although the average temperature increase at the two sites is similar (daily increase 1·5 °C at Toolik and 1·9 °C at Niwot Ridge), at Toolik Lake there is only minor diurnal variation, whereas at Niwot Ridge the daytime increases are extreme on sunny days (as much as 7–10 °C), and minimum nighttime temperatures in the chambers are often slightly cooler than ambient (by about 1 °C). The experimental drifts resulted in wintertime increases in temperature and CO 2 flux. Temperatures under the deep drifts were much more consistent and warmer than in control areas, and at Niwot Ridge remained very close to 0 °C all winter. These increased temperatures were likely responsible for observed increases in system carbon loss. Initial changes to the aboveground biotic system included an increase in growth in response to both snow and warming, despite a reduced growing season. This is expected to be a transient response that will eventually be replaced by reduced growth. At least one species, Kobresia myosuroides , had almost completely died at Niwot Ridge three years after fence construction, whereas other species were increasing. We expect in both the short- and long-term to see the strongest effects of snow at the Niwot Ridge site, and stronger effects of temperature at Toolik Lake.

Hydrological Processes↗