USGS ScienceSearch

SEARCH · USGS Science

Results for “Field Guide Series”

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.

56 records · Page 4Linked to original sources

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

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