SEARCH · USGS Science

Results for “Water Information Series Report”

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 829 records · Page 46Linked to original sources

USGS1980

Data-base system for northern Midwest regional aquifer-system analysis

The U.S. Geological Survey is conducting a study of the Cambrian and Ordovician aquifer system of the northern Midwest as part of a national series of Regional Aquifer-Systems Analysis (RASA). An integral part of this study will be a simulation of the ground-water flow regime using the Geological Survey's three-dimensional finite-difference model. The first step in the modeling effort is the design and development of a systematic set of processes to facilitate the collection, evaluation, manipulation, and use of large quantities of information. A computerized data-base system to accomplish these goals has been completed for the northern Midwest RASA. The input to the data-base system consists of either point values or contour maps of the hydrogeologic data required by the model. Digital samples of contoured surfaces are obtained by machine digitization. Uniformly spaced model-node values are then computed from the discrete data by two-dimensional interpolation. The interpolator uses the regional (long wavelength) characteristics of the data to compute a surface with minimum total curvature. Local (short wavelength) data are merged with the regional surface to produce the model-node values at the desired spacing. Nonuniformly spaced nodal values may be obtained by fitting two-dimensional polynomials (bicubic splines) to surfaces formed by uniformly spaced points and solving for the surface values at the node locations. The uniformly or nonuniformly spaced node values constitute the output of the data-base system and in turn are the input arrays of data values for the ground-water flow model. Management of the data base is facilitated by forming a data file for each model input parameter and for other hydrogeologic data used in the study. The data files are subdivided into elements consisting of individual model-layer arrays. In this form the data base can be readily accessed and edited. In addition to the interpolation programs, software components of the system include programs: (1) to transform geographic coordinates on a Lambert conformal conic projection to cartesian coordinates and vice versa, (2) to machine contour gridded data, (3) to expedite manipulation and editing of data, and (4) to perform various interarray computations. The data-base system has the following attributes: (1) manual handling of data is minimized and machine handling of data is maximized, (2) given unequally spaced point data over the extent of a study area, model input arrays for any reasonable uniform node spacing can be rapidly computed, (3) accuracy of computed node values are generally compatible with accuracy and spatial distribution of point data, (4) a relatively large class of nonuniformly spaced node configurations can be computed, (5) data within data-base files can be easily accessed and readily edited, and (6) the occurrence of data processing errors at various stages of data-base generation is monitored by machine contouring the computed grids. Functioning of the data-base system is illustrated by the sequence of steps required to produce a model-input array of transmissivity, given raw geologic and hydraulic conductivity data.

Water-Resources Investigations Report↗

USGS2015

Geohydrologic and water-quality characterization of a fractured-bedrock test hole in an area of Marcellus shale gas development, Tioga County, Pennsylvania

An integrated analysis of core, geophysical logs, gas isotopes, and specific-depth water-quality samples from the Cherry Flats test hole was used to characterize the stratigraphy, water-bearing zones, and groundwater quality at a site in southern Tioga County, Pennsylvania. The study was completed as a cooperative effort between the Pennsylvania Department of Natural Resources, Bureau of Topographic and Geologic Survey (BTGS) and the U.S. Geological Survey (USGS). The multi-disciplinary characterization of the test hole provided information to aid the bedrock mapping of the Cherry Flats 7.5-minute quadrangle by BTGS, and to help quantify the depth and character of fresh and saline groundwater in an area of shale-gas exploration. The Cherry Flats test hole was cored to a depth of 1,513 feet (ft) below land surface (bls) and cased to 189 ft through the collapsed mine workings of the former Arnot No. 2 underground coal mine. The test hole penetrated 128.0 ft of Allegheny Formation and 154.1 ft of Pottsville Formation of Pennsylvanian age, 564.8 ft of Huntley Mountain Formation of Mississippian and Devonian age, and 666.3 ft of Catskill Formation of Devonian age. Core recovery was nearly 100 percent, except where complete core loss occurred from a depth of 1,231.1 to 1,240.8 ft. Several coal beds and mined-out coal horizons were penetrated in the Allegheny and Pottsville Formations. The test hole penetrated the entire thickness of the Huntley Mountain Formation and was completed in the middle part of the Catskill Formation. Bedding features penetrated by the test hole were estimated to have a strike of 021 degrees and dip about 1.7 degrees to the southeast, consistent with the test-hole location on the north limb of the Blossburg syncline. Most fractures penetrated by the test hole were parallel to bedding, with steeply dipping fractures present but much less common. Fracture density, determined from optical televiewer, acoustic televiewer, and video logs, generally increased with depth from the base of casing to about 400 ft bls, then decreased with depth to the bottom of the hole except for an increase from 506 to 568 ft bls. Very few fractures were penetrated from 600 to 850 ft. The depths of fresh and saline water-bearing zones were identified in the test hole by geophysical-log analysis and, for inflow zones, verified by specific-depth groundwater sampling by the use of a wire-line point sampler. Under ambient conditions and during pumping of the test hole, fresh water entered the hole from fractures at 567 and 580.5 ft bls, within grayish-red siltstone and greenish-gray sandstone, respectively, and flowed upward and exited at fractures from 303 to 319.5 ft; a very minor amount exited into fractures within coal beds at 240.4 and 252 ft bls. Transmissivity, estimated from analysis of the specific-capacity data and flowmeter logs, was about 18 ft2/d for the fracture zones from 567 to 580.5 ft and 6.7 ft2/d for fracture zones from 240.4 to 252 ft bls. The analysis estimated the hydraulic head of the lower zone and that of the upper flow zone was 8 ft higher and 37 ft lower than the composite water level in the test hole, respectively. Water samples of the freshwater inflow from zones at 567 to 580.5 ft bls had a total dissolved solids concentration of 577 mg/L indicating that these zone are in the lower part of the active groundwater flow system. Below the freshwater-bearing zone at 580.5 ft, the flowmeter did not detect any vertical flow in the test hole, and the gradient of the temperature log approached the geothermal gradient, indicating little ambient fluid flow and minimal fracture transmissivity below this depth. However, small seeps of saline water having total dissolved solids concentrations of greater than about 6,200 mg/L at 945 and 946 ft bls, from dark-greenish-gray to greenish-gray silty beds, were delineated by a time series of specific conductancelogs and observed on the video log. A wat

Pennsylvania↗

USGS2009

Archive of digital boomer and CHIRP seismic reflection data collected during USGS field activity 08LCA03 in Lake Panasoffkee, Florida, May 2008

In May of 2008, the U.S. Geological Survey (USGS) conducted geophysical surveys in Lake Panasoffkee, located in central Florida, as part of the USGS Lakes and Coastal Aquifers (LCA) study. This report serves as an archive of unprocessed digital boomer and Compressed High Intensity Radar Pulse (CHIRP)* seismic reflection data, trackline maps, navigation files, Field Activity Collection System (FACS) logs, Geographic Information System (GIS) files, and formal Federal Geographic Data Committee (FGDC) metadata. Filtered and gained (a relative increase in signal amplitude) digital images of the seismic profiles and geospatially corrected interactive profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. *Due to poor data acquisition conditions associated with the lake bottom sediments, only two CHIRP tracklines were collected during this field activity. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG-Y format (Barry and others, 1975) and may be downloaded and processed with commercial or public domain software such as Seismic Unix (SU). Example SU processing scripts and USGS software for viewing the SEG-Y files (Zihlman, 1992) are provided. The USGS Florida Integrated Science Center (FISC) - St. Petersburg assigns a unique identifier to each cruise or field activity. For example, 08LCA03 tells us the data were collected in 2008 for the Lakes and Coastal Aquifers (LCA) study and the data were collected during the third field activity for that study in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity ID. The naming convention used for each seismic line is as follows: yye##a, where 'yy' are the last two digits of the year in which the data were collected, 'e' is a 1-letter abbreviation for the equipment type (for example, b for boomer and c for CHIRP), '##' is a 2-digit number representing a specific track, and 'a' is a letter representing the section of a line if recording was prematurely terminated or rerun for quality or acquisition problems. The boomer plate is an acoustic energy source that consists of capacitors charged to a high voltage and discharged through a transducer in the water. The transducer is towed on a sled floating on the water surface and, when discharged, emits a short acoustic pulse, or shot, which propagates through the water, sediment column, or rock beneath. The acoustic energy is reflected at density boundaries (such as the seafloor, sediment, or rock layers beneath the seafloor), detected by the receiver, and recorded by a PC-based seismic acquisition system. This process is repeated at timed intervals (for example, 0.5 s) and recorded for specific intervals of time (for example, 100 ms). In this way, a two-dimensional (2-D) vertical profile of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the boomer acquisition geometry. The EdgeTech SB-424 CHIRP system used for this survey has a vertical resolution of 4 - 8 cm, a penetration depth that is usually less than 2 m beneath the seafloor, and uses a signal of continuously varying frequency. The towfish is a sound source and receiver, which is typically towed 2 - 5 m above the seafloor. The acoustic energy is reflected at density boundaries (such as the seafloor or sediment layers beneath the seafloor), detected by a receiver, and recorded by a PC-based seismic acquisition system. This process is repeated at timed intervals (for example, 0.125 s) and recorded for specific intervals of time (for example, 50 ms); the resulting profile is a two-dimensional vertical image of the shallow geologic structure beneath the ship track. Figure 2 displays the acquisition geometry for the CHIRP system. Refer to table 1 for a summary of acquisition parameters and table 2 for trackline statistics.

Florida↗

USGS2014

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

USGS2024

Annotated bibliography of scientific research relevant to oil and gas reclamation best management practices in the western United States, published from 1969 through 2020

Integrating recent scientific knowledge into management decisions supports effective natural resource management and can lead to better resource outcomes. However, finding and accessing scientific knowledge can be time consuming and costly. To assist in this process, the U.S. Geological Survey has created a series of annotated bibliographies on topics of management concern for lands in the western United States (U.S.). Oil and gas development on public lands is a long-standing and substantial component of local and regional economies and has expanded in recent decades, particularly on public lands in the western U.S. This development is associated with extensive networks of pipelines, roads, and processing facilities, across which reclamation is Federally mandated following initial well pad development (“interim” reclamation) and once resource extraction is complete (“final” reclamation). Reclamation is critical for recovering ecological services to energy-affected lands, including vegetation productivity, wildlife habitat, water and air quality, and soil stability (for example, resistance to wind and water erosion). However, reclamation of oil and gas affected lands in the western U.S. has proved challenging due to an array of regulatory and environmental factors, such as minimally developed soils, short growing seasons, herbivory, high winds, invasive species, rugged terrain, and in particular, arid climates associated with low total precipitation, high evapotranspiration rates, and highly variable precipitation patterns. We compiled and summarized journal articles, government reports, technical reports, proceedings, and theses and dissertations relevant to oil and gas reclamation. We constrained our search to products published on or before December 31, 2020 but did not limit our search by a starting date; the earliest product resulting from this effort was published in March 1969. Second, we manually scanned the last 15 years (2005-2020) of tables of contents in journals, bibliographies, and proceedings of which we were aware would contain articles highly relevant to this bibliography. We carried out the search for these products through multiple means: (1) performing a structured search of two reference databases, (2) examining articles published since 2005 in highly relevant scientific journals and conference proceedings, and (3) reviewing additional material suggested by authors of products identified in steps 1 and 2. Our search was intentionally broad in order to identify as much relevant work as possible, much of which is professionally applied and tested within the industry of oil and gas reclamation, but which remains unpublished in scientific journals. We refined the initial list of products by removing: (1) duplicates, (2) products not written in English, (3) products that were not relevant to the arid ecosystems of western North America, (4) products that were not released as research, data products, or review articles in journals or as formal scientific reports, and (5) products with data which were not relevant to reclamation of oil and gas-affected lands, or for which the study did not present new data, findings, or syntheses relevant to reclamation of oil and gas-affected lands. We summarized each product using a consistent structure (background, objectives, methods, location, findings, and implications) and assigned standardized management topics to each. Management topics are intended to aid online searching within the bibliography and are described in more detail in the Methods Section of this report; they include what type of disturbance the product addresses (well pads, mining, pipelines), what aspect of oil and gas reclamation they pertain to (practices, standards, monitoring), what type of data are present in the product (for instance soil or vegetation recovery data), and an indication if the product were from a source other than a published, peer-reviewed outlet (such as dissertations or unpublished professional reports – these are identified as grey literature). The review process for this annotated bibliography included an initial internal colleague review of each summary, requesting input on each summary from an author of the original product, and a formal peer-review. Our initial searches resulted in 3,197 total products, of which 290 met our criteria for inclusion. “Reclamation Practices” is by far the management topic most addressed, followed by “Reclamation Monitoring,” for example, products assessing what and how monitoring methods are used to track and measure reclamation outcome. This document may be accessed at https://doi.org/10.3133/ofr20231068 or from the U.S. Geological Survey Publication Warehouse ( https://pubs.usgs.gov/ ). The 1-page product summaries herein will also be used to create a bibliography at https://apps.usgs.gov/science-for-resource-managers that includes links to each original product, where available, and in which subject matter will be searchable by topic, location, and year. The studies compiled and summarized here may inform planning and management actions that seek to reclaim landscapes across the western U.S. which have been affected by oil and gas development.

Open-File Report↗

USGS2010

Aquatic assessment of the Ely Copper Mine Superfund site, Vershire, Vermont

The Ely Mine, which operated from 1821 to 1905, and its area of downstream impact constitute the Ely Copper Mine Superfund site. The site was placed on the National Priorities List in 2001. The mine comprises underground workings, foundations from historical structures, several waste-rock piles, roast beds associated with the smelting operation, and slag piles resulting from the smelting. The mine site is drained by Ely Brook, which includes several tributaries, one of which drains a series of six ponds. Ely Brook empties into Schoolhouse Brook, which flows 3.3 kilometers and joins the Ompompanoosuc River. The aquatic ecosystem at the site was assessed using a variety of approaches that investigated surface-water quality, sediment quality, and various ecological indicators of stream-ecosystem health. The degradation of surface-water quality is dominated by copper with localized effects caused by iron, aluminum, cadmium, and zinc. Chronic water-quality criteria for copper are exceeded in the surface water of four of the six ponds on the Ely Brook tributary, and all of Ely Brook and Schoolhouse Brook, and of the Ompompanoosuc River downstream of the confluence with Schoolhouse Brook. Comparison of hardness-based and Biotic Ligand Model-based water-quality criteria for copper yields similar results with respect to extent of impairment. However, the Biotic Ligand Model criteria are mostly lower than the hardness-based criteria and thus suggest a greater degree of impairment, particularly in the Ely Brook watershed, where dissolved organic carbon concentrations and pH values are lower. Surface-water toxicity testing correlates strongly with the extent of impact. Likewise, riffle-habitat benthic invertebrate richness and abundance data support these results through the stream environment. Similarly, the index of biotic integrity for the fish community in Schoolhouse Brook and the Ompompanoosuc River document degraded habitats throughout Schoolhouse Brook from Ely Brook down to the Ompompanoosuc River. The sediment environment shows similar extents of impairment also dominated by copper, although localized degradation due to chromium, nickel, lead, and zinc was documented on the basis of probable effects concentrations. In contrast, equilibrium-partitioning sediment benchmarks indicate no toxic effects would be expected in sediments at the reference sites, and uncertain toxic effects throughout Ely Brook and Schoolhouse Brook, except for the reference sites and site EB-600M. The results for site EB-600M indicate predicted toxic effects. Acute toxicity testing of in situ pore waters using Hyalella azteca indicates severe impacts in Ely Brook reaching 100 percent lethality at site EB-90M. Acute toxicity testing of in situ pore waters using Chironomus dilutus shows similar, but not as severe, toxicity. Neither set of in situ pore-water toxicity tests showed significant impairment in Schoolhouse Brook or the Ompompanoosuc River. Chronic sediment toxicity testing using Hyalella azteca indicated significant toxicity in Ely Brook, except at site EB-90M, and in Schoolhouse Brook. The low toxicity of EB-90M may be a reflection of the low lability of copper in that sediment as indicated by a low proportion of extractable copper (1.1 percent). Depositional-targeted habitat invertebrate richness and abundance data support these conclusions for the entire watershed, as do the index of biotic integrity data from the fish community. The information was used to develop an overall assessment of the impact on the aquatic system that appears to be a result of the acid rock drainage at the Ely Mine. More than 700 meters of Ely Brook, including two of the six ponds, were found to be severely impacted, on the basis of water-quality data and biological assessments. The reference location was of good quality based on the water quality and biological assessment. More than 3,125 meters of Schoolhouse Brook are also severely impacted, on the basis of water-quality data and biological assessments. The biological community begins to recover near the confluence with the Ompompanoosuc River. The evidence is less conclusive regarding the Ompompanoosuc River. The sediment data suggest that the sediments could be a source of toxicity in Ely Brook and Schoolhouse Brook. The surface-water assessment is consistent with the outcome of a surface-water toxicity testing program performed by the U.S. Environmental Protection Agency for Ely Brook and Schoolhouse Brook and a surface-water toxicity testing program and in situ amphibian testing program for the ponds.

Vermont↗

USGS2016

Decision analysis to support development of the Glen Canyon Dam long-term experimental and management plan

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, National Park Service, and Argonne National Laboratory, completed a decision analysis to use in the evaluation of alternatives in the Environmental Impact Statement concerning the long-term management of water releases from Glen Canyon Dam and associated management activities. Two primary decision analysis methods, multicriteria decision analysis and the expected value of information, were used to evaluate the alternative strategies against the resource goals and to evaluate the influence of uncertainty. A total of 18 performance metrics associated with 8 out of 12 resource goals (fundamental objectives) were developed by the Bureau of Reclamation and National Park Service in partnership with subject-matter teams composed of Federal, State, tribal, and private experts. A total of 19 long-term strategies associated with 7 alternatives were developed by the Bureau of Reclamation, National Park Service, Argonne National Laboratory, U.S. Geological Survey, and Cooperating Agencies. The 19 long-term strategies were evaluated against the 18 performance metrics using a series of coupled simulation models, taking into account the effects of several important sources of uncertainty. A total of 27 Federal, State, tribal, and nongovernmental agencies were invited by the Assistant Secretary of Interior to participate in a swing-weighting exercise to understand the range of perspectives about how to place relative value on the resource goals and performance metrics; 14 of the 27 chose to participate. The results of the swing-weighting exercise were combined with the evaluation of the alternatives to complete a multicriteria decision analysis. The effects of uncertainty on the ranking of long-term strategies were evaluated through calculation of the value of information. The alternatives and their long-term strategies differed across performance metrics, producing unavoidable tradeoffs; thus, there was no long-term strategy that was dominated by another across all performance metrics. When the performance of each alternative was weighted across performance metrics, three alternatives (B, D, and G) were top-ranked depending on the set of weights proposed: Alternative B was favored by those stakeholders that placed a high value on hydropower; Alternative G was favored by those stakeholders that placed a high value on the restoration of natural processes, like beachbuilding and natural vegetation; and Alternative D was favored by the remaining stakeholders. Surprisingly, these rankings were not sensitive to the critical uncertainties that were evaluated; that is, the choice of a preferred long-term strategy was sensitive to the value-based judgment about how to place relative weight on the resource goals but was not sensitive to the uncertainties in the system dynamics that were evaluated in this analysis. The one area of uncertainty that did slightly affect the ranking of alternatives was the long-term pattern of hydrological input; because of this sensitivity, some attention to the possible effects of climate change is warranted. The results of the decision analysis are meant to serve as only one of many sources of information that can be used to evaluate the alternatives proposed in the Environmental Impact Statement. These results only focus on those resource goals for which quantitative performance metrics could be formulated and evaluated; there are other important aspects of the resource goals that also need to be considered. Not all the stakeholders who were invited to participate in the decision analysis chose to do so; thus, the Bureau of Reclamation, National Park Service, and U.S. Department of Interior may want to consider other input.

Arizona↗

USGS2022

Estimating stream temperature in the Willamette River Basin, northwestern Oregon—A regression-based approach

The alteration of thermal regimes, including increased temperatures and shifts in seasonality, is a key challenge to the health and survival of federally protected cold-water salmonids in streams of the Willamette River basin in northwestern Oregon. To better support threatened fish species, the U.S. Army Corps of Engineers (USACE) and other water managers seek to improve the thermal regime in the Willamette River and key tributaries downstream of USACE dams by utilizing strategically timed flow releases from USACE dams. To inform flow management decisions, regression relations were developed for 12 Willamette River basin locations below USACE dams relating stream temperature with streamflow and air temperature utilizing publicly available datasets spanning 2000–18. The resulting relations provide simple tools to investigate stream temperature responses to changes in streamflow and climatic conditions in the Willamette River system. Regression relations on the Willamette River and key tributaries show that, at locations sufficiently distant from the direct temperature influence of upstream dam releases, air temperature and streamflow are reasonable proxies to predict the 7-day average of the daily mean (7dADMean) and 7-day average of the daily maximum (7dADMax) water temperature with errors generally ≤1 degrees Celsius (°C). To account for seasonal variations in the relation between air temperature, streamflow, and stream temperature, a transition-smoothed, seasonal regression approach was used. Stream temperature is inversely correlated with streamflow in all seasons except “winter” (January–March), when it is relatively independent. Stream temperature is positively correlated with air temperature in all seasons, but the slope decreases at very low or very high air temperatures. Generally, fit is best for seasonal models “winter” (January–March), “spring” (April–May), “summer” (June–August), and “early autumn” (September–October). Error in “autumn” (November–December) is larger, probably due to variation in the onset timing of winter storms. Simulated results from a climatological analysis of predicted stream temperature suggest that, excluding extremes and accounting for some seasonal variability, the 7dADMean and 7dADMax stream temperature sensitivity to air temperature and streamflow varies by location on the river. To investigate the potential range of stream temperature variability based on historical air temperature and streamflow conditions, stream temperature predictions were calculated using synthetic time series comprised of daily temperature values representing the 0.10, 0.33, 0.50, 0.67, and 0.90 quantile of air temperature and streamflow from 1954 (the year meaningful streamflow augmentation began) to 2018. Results show that from a “very hot” (0.90 quantile) and “very dry” (0.10 quantile) year to a “very cool” (0.10 quantile) and “very wet” (0.90; all quantiles from 1954 to 2018) year, the stream temperature sensitivity to air temperature and streamflow is about 3 °C at Harrisburg (river mile 161.0) and increases to about 5 °C at Keizer (river mile 82.2). While the number of days exceeding regulatory criteria are fewer in cooler, wetter years than in warmer, dryer years, the models suggest that the Willamette River will likely continue to exceed the State of Oregon maximum water-temperature criterion of 18 °C for sustained periods from late spring to early autumn and that the flow management practices evaluated in this study, while effective at influencing stream temperature, likely cannot prevent many or all such exceedances. As modeled for 2018, a representative very hot year with normal to below-normal streamflow, stream temperature sensitivity to changes in streamflow of ±100 to ±1000 cubic feet per second produced mean monthly temperature changes from 0.0 to 1.4 °C at Keizer, Albany, and Harrisburg during summer. For a specified change in flow, temperature sensitivity is greater at upstream locations where streamflow is less than that at downstream locations because the change in streamflow is a greater percentage of total streamflow at upstream locations. Similarly, temperature response to a set change in flow is greater in the summer and early autumn low-flow season than in spring when flows are higher. The regression models developed in this study thus indicate that flow management is likely to have a greater effect on stream temperature at upstream locations (such as Harrisburg or Albany) and during the low-flow season than at downstream locations (such as Keizer) or during periods of higher streamflow.

Oregon↗

USGS2018

Analysis of groundwater response to tidal fluctuations, Operable Unit 2, Area 8, Naval Base Kitsap, Keyport, Washington

Operable Unit 2, Area 8, at Naval Base Kitsap, Keyport is the site of a former chrome-plating facility that released metals (primarily chromium and cadmium), chlorinated volatile organic compounds, and petroleum compounds into the local environment. To ensure long-term protectiveness, as stipulated in the Fourth Five-Year Review for the site, Naval Facilities Engineering Command Northwest collaborated with the U.S. Environmental Protection Agency, the Washington State Department of Ecology, and the Suquamish Tribe, to collect data to monitor the contamination left in place and to ensure the site does not pose a risk to human health or the environment. To support these efforts, refined information was needed on the interaction of fresh groundwater with seawater in response to the up-to 13-ft tidal fluctuations at this nearshore site adjacent to Port Orchard Bay. The information was analyzed to meet the primary objective of this investigation, which was to determine the optimal time during the semi-diurnal and the neap-spring tidal cycles to sample groundwater for freshwater contaminants in Area 8 monitoring wells. Groundwater levels and specific conductance in five monitoring wells, along with marine water-levels (tidal levels) in Port Orchard Bay, were monitored every 15 minutes during a 3-week duration to determine how nearshore groundwater responds to tidal forcing. Time series data were collected from October 24, 2017, to November 16, 2017, a period that included neap and spring tides. Vertical profiles of specific conductance were also measured once in the screened interval of each well prior to instrument deployment to determine if a freshwater/saltwater interface was present in the well during that particular time. The vertical profiles of specific conductance were measured only one time during an ebbing tide at approximately the top, middle, and bottom of the saturated thickness within the screened interval of each well. The landward-most well, MW8-8, was completely freshwater, while one of the most seaward wells, MW8-9, was completely saline. A distinct saltwater interface was measured in the three other shallow wells (MW8-11, MW8-12, and MW8-14), with the topmost groundwater occurring fresh underlain by higher conductivity water. Lag times between minimum spring-tide level and minimum groundwater levels in wells ranged from about 2 to 4.5 hours in the less-than 20-ft deep wells screened across the water table, and was about 7 hours for the single 48-ft deep well screened below the water table. Those lag times were surprisingly long considering the wells are all located within 200-ft of the shoreline and the local geology is largely coarse-grained glacial outwash deposits. Various manmade subsurface features, such as slurry walls and backfilled excavations, likely influence and confuse the connectivity between seawater and groundwater. The specific-conductance time-series data showed clear evidence of substantial saltwater intrusion into the screened intervals of most shallow wells. Unexpectedly, the intrusion was associated with the neap part of the tidal cycle around November 13–16, when relatively low barometric pressure and high southerly winds led to the highest high and low tides measured during the monitoring period. The data consistently indicated that the groundwater had the lowest specific conductance (was least mixed with seawater) during the prior neap tides around October 30, the same period when the shallow groundwater levels were lowest. Although the specific conductance response is somewhat different between wells, the data do suggest that it is the heights of the actual high-high and low-low tides, regardless of whether or not they occur during the neap or spring part of the cycle, that allows seawater intrusion into the nearshore aquifer at Area 8. With all the data taken into consideration, the optimal time for sampling the shallow monitoring wells at Area 8 would be centered on a 2–5-hour period following the predicted low-low tide during neap tide, with due consideration of local atmospheric pressure and wind conditions that have the potential to generate tides that can be substantially higher than those predicted from lunar-solar tidal forces. The optimal time for sampling the deeper monitoring wells at Area 8 would be during the 6–8-hour period following a predicted low-low tide, also during the neap tide part of the tidal cycle. The specific time window to sample each well following a low tide can be found in table 5. Those periods are when groundwater in the wells is most fresh and least diluted by seawater intrusion. In addition to timing, consideration should be given to collecting undisturbed samples from the top of the screened interval (or top of the water table if below the top of the interval) to best characterize contaminant concentrations in freshwater. A downhole conductivity probe could be used to identify the saltwater interface, above which would be the ideal depth for sampling.

Washington↗

USGS1914

The Transportation of Debris by Running Water

Scope.-The finer debris transported by a stream is borne in suspension. The coarser is swept along the channel bed. The suspended load is readily sampled and estimated, and much is known as to its quantity. The bed load is inaccessible and we are without definite information as to its amount. The primary purpose of the investigation was to learn the laws which control the movement of bed load, and especially to determine how the quantity of load is related to the stream's slope and discharge and to the degree of comminution of the debris. Method.-To this end a laboratory was equipped at Berkeley, Cal., and experiments were performed in which each of the three conditions mentioned was separately varied and the resulting variations of load were observed and measured. Sand and gravel were sorted by sieves into grades of uniform size. Determinate discharges were used. In each experiment a specific load was fed to a stream of specific width and discharge, and measurement was made of the slope to which the stream automatically adjusted its bed so as to enable the current to transport the load. The slope factor.-For each combination of discharge, width, and grade of debris there is a slope, called competent slope, which limits transportation. With lower slopes there is no load, or the stream has no capacity for load. With higher slopes capacity exists; and increase of slope gives increase of capacity. The value of capacity is approximately proportional to a power of the excess of slope above competent slope. If S equal the stream's slope and sigma equal competent slope, then the stream's capacity varies as (S - sigma)n. This is not a deductive, but an empiric law. The exponent n has not a fixed value, but an indefinite series of values depending on conditions. Its range of values in the experience of the laboratory is from 0.93 to 2.37, the values being greater as the discharges are smaller or the debris is coarser. The discharge factor.-For each combination of width, slope, and grade of debris there is a competent discharge, k. Calling the stream's discharge Q, the stream's capacity varies as (Q - k)o. The observed range of values for o is from 0.81 to 1.24, the values being greater as the slopes are smaller or the debris is coarser. Under like conditions o is less than n; or, in other words, capacity is less sensitive to change3 of discharge than to changes of slope. The fineness factor.-For each combination of width, slope, and discharge there is a limiting fineness of debris below which no transportation takes place. Calling fineness (or degree of comminution) F and competent fineness o, the stream's capacity varies with (F - o)p. The observed range of values for p is from 0.50 to 0.62, the values being greater as slopes and discharges are smaller. Capacity is less sensitive to changes in fineness of debris than to changes in discharge or slope. The form factor.-Most of the experiments were with straight channels. A few with crooked channels yielded nearly the same estimates of capacity. The ratio of depth to width is a more important factor. For any combination of slope, discharge, and fineness it is possible to reduce capacity to zero by making the stream very wide and shallow or very narrow and deep. Between these extremes is a particular ratio of depth to width, p, corresponding to a maximum capacity. The values of p range, under laboratory conditions, from 0.5 to 0.04, being greater as slope, discharge, and fineness are less. Velocity.-The velocity which determines capacity for bed load is that near the stream's bed, but attempts to measure bed velocity were not successful. Mean velocity was measured instead. To make a definite comparison between capacity and mean velocity it is necessary to postulate constancy in some accessory condition. If slope be the constant, in which case velocity changes with discharge, capacity varies on the average with the 3.2 power of velocity. If discharge be the constant, in w

Professional Paper↗

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↗

USGS2010

Glaciers of Asia

This chapter is the ninth to be released in U.S. Geological Survey Professional Paper 1386, Satellite Image Atlas of Glaciers of the World, a series of 11 chapters. In each of the geographic area chapters, remotely sensed images, primarily from the Landsat 1, 2, and 3 series of spacecraft, are used to analyze the specific glacierized region of our planet under consideration and to monitor glacier changes. Landsat images, acquired primarily during the middle to late 1970s and early 1980s, were used by an international team of glaciologists and other scientists to study various geographic regions and (or) to discuss related glaciological topics. In each glacierized geographic region, the present areal distribution of glaciers is compared, wherever possible, with historical information about their past extent. The atlas provides an accurate regional inventory of the areal extent of glacier ice on our planet during the 1970s as part of a growing international scientific effort to measure global environmental change on the Earth’s surface. The chapter is divided into seven geographic parts and one topical part: Glaciers of the Former Soviet Union (F–1), Glaciers of China (F–2), Glaciers of Afghanistan (F–3), Glaciers of Pakistan (F–4), Glaciers of India (F–5), Glaciers of Nepal (F–6), Glaciers of Bhutan (F–7), and the Paleoenvironmental Record Preserved in Middle-Latitude, High-Mountain Glaciers (F–8). Each geographic section describes the glacier extent during the 1970s and 1980s, the benchmark time period (1972–1981) of this volume, but has been updated to include more recent information. Glaciers of the Former Soviet Union are located in the Russian Arctic and various mountain ranges of Russia and the Republics of Georgia, Kyrgyzstan, Tajikistan, and Kazakstun. The Glacier Inventory of the USSR and the World Atlas of Ice and Snow Resources recorded a total of 28,881 glaciers covering an area of 78,938 square kilometers (km 2 ). China includes many of the mountain-glacier systems of the world including the Himalaya, Karakorum, Tien Shan and Altay mountain ranges. The glaciers are widely scattered and cover an area of about 59,425 km 2 . The mountain glaciers may be classified as maritime, subcontinental or extreme continental. In Afghanistan, more than 3,000 small glaciers occur in the Hindu Kush and Pamir mountains. Most glaciers occur on north-facing slopes shaded by mountain peaks and on east and southeast slopes that are shaded by monsoon clouds. The glaciers provide vital water resources to the region and cover an area of about 2,700 km 2 . Glaciers of northern Pakistan are some of the largest and longest mid-latitude glaciers on Earth. They are located in the Hindu Kush, Himalaya, and Karakoram mountains and cover an area of about 15,000 km 2 . Glaciers here are important for their role in providing water resources and their hazard potential. The glaciers in India are located in the Himalaya and cover about 8,500 km 2 . The Himalaya contains one of the largest reservoirs of snow and ice outside the polar regions. The glaciers are a major source of fresh water and supply meltwater to all the rivers in northern India, thereby affecting the quality of life of millions of people. In Nepal, the glaciers are located in the Himalaya as individual glaciers; the glacierized area covers about 5,324 km 2 . The region is the highest mountainous region on Earth and includes the Mt. Everest region. Glaciers in the Bhutan Himalaya have a total area of about 1,317 km 2 . Many recent glacier studies are focused on glacier lakes that have the potential of generating dangerous glacier lake outburst floods. Research on the glaciers of the middle-latitude, high-mountain glaciers of Asia has also focused on the information contained in the ice cores from the glaciers. This information helps in the reconstruction of paleoclimatic records, and the computer modeling of global climate change.

Professional Paper↗

USGS2004

Hydrogeologic framework of the North Fork and surrounding areas, Long Island, New York

Ground water on the North Fork of Long Island is the sole source of drinking water, but the supply is vulnerable to saltwater intrusion and upconing in response to heavy pumping. Information on the area's hydrogeologic framework is needed to analyze the effects of pumping and drought on ground-water levels and the position of the freshwater-saltwater interface. This will enable water-resource managers and water-supply purveyors to evaluate a wide range of water-supply scenarios to safely meet water-use demands. The extent and thickness of hydrogeologic units and position of the freshwater-saltwater interface were interpreted from previous work and from exploratory drilling during this study. The fresh ground-water reservoir on the North Fork consists of four principal freshwater flow systems (referred to as Long Island mainland, Cutchogue, Greenport, and Orient) within a sequence of unconsolidated Pleistocene and Late Cretaceous deposits. A thick glacial-lake-clay unit appears to truncate underlying deposits in three buried valleys beneath the northern shore of the North Fork. Similar glacial-lake deposits beneath eastern and east-central Long Island Sound previously were inferred to be younger than the surficial glacial deposits exposed along the northern shore of Long Island. Close similarities in thickness and upper-surface altitude between the glacial-lake-clay unit on the North Fork and the glacial-lake deposits in Long Island Sound indicate, however, that the two are correlated at least along the North Fork shore. The Matawan Group and Magothy Formation, undifferentiated, is the uppermost Cretaceous unit on the North Fork and constitutes the Magothy aquifer. The upper surface of this unit contains a series of prominent erosional features that can be traced beneath Long Island Sound and the North Fork. Northwest-trending buried ridges extend several miles offshore from areas southeast of Rocky Point and Horton Point. A promontory in the irregular, north-facing cuesta slope extends offshore from an area southwest of Mattituck Creek and James Creek. Buried valleys that trend generally southeastward beneath Long Island Sound extend onshore northeast of Hashamomuck Pond and east of Goldsmith Inlet. An undifferentiated Pleistocene confining layer, the lower confining unit, consists of apparently contiguous units of glacial-lake, marine, and nonmarine clay. This unit is more than 200 feet thick in buried valleys filled with glacial-lake clay along the northern shore, but elsewhere on the North Fork, it is generally less than 50 feet thick and presumably represents an erosional remnant of marine clay. Its upper surface is generally 75 feet or more below sea level where it overlies buried valleys, and is generally 100 feet or less below sea level in areas where marine clay has been identified. A younger unit of glacial-lake deposits, the upper confining unit, is a local confining layer and underlies a sequence of late Pleistocene moraine and outwash deposits. This unit is thickest (more than 45 feet thick) beneath two lowland areas--near Mattituck Creek and James Creek, and near Hashamomuck Pond--but pinches out close to the northern and southern shores and is locally absent in inland areas of the North Fork. Its upper-surface altitude generally rises to near sea level toward the southern shore. Freshwater in the Orient flow system is limited to the upper glacial aquifer above the top of the lower confining unit. The upper confining unit substantially impedes the downward flow of freshwater in inland parts of the Greenport flow system. Deep freshwater within the lower confining unit in the east-central part of the Cutchogue flow system probably is residual from an interval of lower sea level. The upper confining unit is absent or only a few feet thick in the west-central part of the Cutchogue flow system and does not substantially impede the downward flow of freshwater, but the lower confining unit probably impedes the downward flow of freshwater within a southeast-trending buried valley in this area.

New York↗

USGS2020

Smallmouth buffalo (Ictiobus bubalus) growth across a 1200km human use and ecological disturbance gradient in the Upper Mississippi River System

Smallmouth buffalo (Ictiobus bubalus) is a common and widely distributed large-bodied species of the family Catostomidae. It inhabits large rivers and reservoirs of the eastern continental United States (east of the continental Divide) and is most abundant and common in the large rivers of the Midwest and Central Plains, though it does occur as far north and east as the Hudson Bay drainage and as far south and west as Arizona (Edwards and Twoney 1982). Historically, smallmouth buffalo were an important component of commercial fisheries on both the Mississippi and Illinois Rivers. However, following the introduction of common carp (Cyprinus carpio) in the mid-1800s (Carlander 1954), the construction of a system of navigation dams on Upper Mississippi and Illinois River in the 1930s (USGS 1999), and water quality/pollution issues through the 1980s (Weiner 2010), the role of smallmouth buffalo in the overall UMRS fish community and commercial fishery has generally diminished relative to historical standards. Still, smallmouth buffalo remains an important and valued component of the UMRS commercial fishery. The study area is represented by three study reaches on the Illinois River and three study reaches on the Upper Mississippi River (Figure 1). Collectively, these study reaches represent nearly 1200 river km and exist across strong and pronounced ecological and disturbance gradients. For example, habitat composition, water quality, commercial navigation intensity, aquatic plant prominence, and the number and abundance of nonnative fish species vary strongly across the study domain, with northern Mississippi River reaches exhibiting less navigation traffic, better water quality, markedly greater aquatic plant prominence, more diverse habitat composition, and comparably much smaller numbers of nonnative species than the lower Mississippi River study reach and those on the Illinois River (USGS 1999; Johnson and Hagerty [eds] 2008; Irons et al. 2009). Long term monitoring efforts conducted under the auspices of the Upper Mississippi River Restoration program over the past 27 years have provided tremendous insights into shifts and changes of the overall UMRS fish community (Ickes et al. 2005; Garvey et al. 2010; Schramm and Ickes 2016). However, these monitoring efforts observe only the most basic aspects of the UMRS fish community (i.e., catch, length, weight, distribution, and occurrence). To gain a greater understanding of forces driving community level shifts and changes, more directed study is needed on the functional attributes of fish populations (i.e., growth, mortality, recruitment). Collectively, these functional attributes of populations are termed population dynamics and/or vital rates. It is important to note, the population dynamics of fishes in large rivers is generally poorly understood, especially for non-game species (Ickes 2018). The prevailing view is that abiotic factors largely govern inter-annual population dynamics, typically based upon rather short-term observations and correlations with assorted abiotic river attributes that vary on a seasonal or annual basis (for example, Risotto and Turner 1985). However, the role that longer-term abiotic factors play in regulating population abundance, or that biotic factors internal to the population (e.g., spawner-recruit dynamics, growth dynamics) or external to the population (e.g., predator-prey dynamics, sympatric competitors, disease) remain poorly understood. Achieving a greater understanding of these dynamics is important for stock, game, and invasive species management. In 2017, as part of a larger study designed to gain vital population rate information for smallmouth buffalo in the Upper Mississippi and Illinois Rivers (“Smallmouth Buffalo population demographics of the Upper Mississippi River System”; UMRR LTRM 2018SOW project items 2018MMBF1-2018MMBF6) annual growth patterns in smallmouth buffalo were determined and evaluated. This was accomplished by measuring growth histories recorded in annual growth increments on hard bony parts (here otoliths), a method known generically as biochronology, and somewhat analogous to dendrochronology practiced by foresters. These methods allow one to generate time-series of annual growth histories that depend upon age, year class (i.e., cohort), and annual environmental conditions experienced by the population over time (Weisberg, 1993). Biochronology methods were used to develop a 36-year time series of smallmouth buffalo growth in the Upper Mississippi and Illinois Rivers across a 1200 km ecological and human use disturbance gradient. Annual growth intervals were identified and measured from otoliths to determine fish age and growth history. A mixed model that parses the growth increment into age and year effects was fit to these data. Given the pronounced ecological and disturbance gradients inherent to the UMRS and the study domain, an a priori expectation of differing patterns in growth is accepted as a null hypothesis to test. The goal of this study was to model smallmouth buffalo growth as a function of the age of the fish and the growth year in which the growth was gained. The primary modeling objective was to parse growth observed on each annulus into a portion attributable to the age of the fish and the portion attributable to the year in which the growth was gained. In effect, this modeling approach removes the somewhat trivial age effects on growth so that a non-confounded growth year effect can be gained. Results attributable to growth year provide a time series of growth information that is of the same duration as the oldest fish observed and solely reflects environmental influences on growth. These model responses can then be investigated relative to environmental covariate time-series suspected of influencing growth of smallmouth buffalo in the Upper Mississippi and Illinois Rivers (e.g., temperature, discharge, population density, population mortality, forage availability, sympatric competition, habitat composition, navigation intensity, nonnative fish prominence, etc.). Thus, the primary scientific objective was to investigate if and how smallmouth buffalo growth varies in accordance with innate ecological and disturbance gradients across the study domain.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

USGS1988

The Manti, Utah, landslide

PART A: The Manti landslide is in Manti Canyon on the west side of the Wasatch Plateau in central Utah. In early June 1974, coincident with the melting of a snowpack, a rock slump/debris flow occurred on the south rim of Manti Canyon. Part of the slumped material mixed with meltwater and mobilized into a series of debris flows that traveled down the slope a distance of as much as 1.2 km. Most of the flows were deposited either at the base of the steep rocks of the canyon rim or at the site of an old, silted reservoir. A small part of the debris flow deposit stopped on the head of the very large, relatively inactive Manti landslide. The upper part of the landslide began moving as cracks propagated downslope. A little more than a year later, August 1975, movement extended the full length of the old landslide, and about 19 million m 3 of debris about 3 km long and as much as 800 m wide threatened to block the canyon. The upper part of the landslide apparently had moved small amounts between 1939 and 1974. This part of the landslide, identifiable on pre-1974 aerial photographs, consisted of well-defined linears on the landslide flanks and two large internal toe bulges about 2 km downslope from the head. The abrupt reactivation in 1974 proceeded quickly after the debris flows had provided a surcharge in the head and crown area. Movement propagated downslope at 4-5 m/h for the first few days following reactivation. During 1974, the reactivation probably encompassed all the parts of the landslide that had moved small amounts between 1939 and 1974. Movement nearly or completely stopped during the winter of 1974-75, but began again in the spring of 1975. The landslide enlarged from the flanks of the internal toe bulges to Manti Creek at a rate of 2-3 m/h. Movement stopped again during the winter of 1975-76 and began again in the spring of 1976. Thereafter, the displacements have been small compared to earlier. The displacement rates for the landslide were variable depending on where and when they were measured. At the waterline crossing, about 500 m downslope from the head, the maximum rate was about 1 mid, and the peak rate occurred within the first 30 days following reactivation. Almost a year later, during the spring of 1975, the rates there were 0.1-0.3 mid. By the time movement had propagated 2.5 km downslope to Manti Creek, more than 40 m of displacement had occurred at the waterline. Cracks were first noted at Manti Creek on June 21,1975, but movement was initially very slow there. The maximum rate of sliding of about 3.1 mid occurred during the period September 1-19,1975, and the movement decreased and stopped during the winter of 1975-76. At the time the lower part of the landslide was moving rapidly, the rate farther upslope was small. The result was that the landslide changed from being in compression, which was caused by loading from above, to being in extension, which was caused by the lower part moving faster than the upper part. One of the results of the extension was that the landslide pulled apart on a steeper segment of the slope and exposed the failure surface. Since the fall of 1975, the landslide has been separated into two independent parts. There has been no movement in the lower part, whereas movements in the upper part have continued at a rate of a few meters per year. Although the landslide is 6.5 km from the nearest permanent dwelling, it cost nearly $2 million in actual damages and in preparation for a major disaster that did not occur. The waterline for the city of Manti was replaced and a well was drilled to provide an emergency water supply. PART B: The Manti landslide abuts, at its toe, another large landslide called the North slide. In 1975, public officials were concerned that continued movement of the Manti landslide would trigger reactivation of the North slide. In response to this threat, four borings were placed in the North slide to obtain samples for testing and information

Utah↗

USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗

USGS1993

Field observations, preliminary model analysis, and aquifer thermal efficiency

In May 1980, the University of Minnesota began a project to evaluate the feasibility of storing heated (150 degrees Celsius (°C) water in the deep (180 to 240 meters (m)) Franconia-Ironton-Galesville aquifer and later recovering it for space heating. The Aquifer Thermal-Energy Storage (ATES) system was a doublet-well design in which the injection and withdrawal wells were spaced approximately 250 m apart. High-temperature water from the university's steam-generation facilities supplied heat for injection. Water was pumped from one of the wells through a heat exchanger, where heat was added or removed. Water then was injected back into the aquifer through the other well. The experimental plan for testing the ATES system consisted of a series of short-term hot-water injection, storage, and withdrawal cycles. Each cycle was 24 days long, and each injection, storage, and withdrawal step of the cycle was 8 days. The Franconia-Ironton-Galesville aquifer is a consolidated sandstone, approximately 60 m thick, the top of which is approximately 180 m below the land surface. It is confined above by the St. Lawrence Formation a dolomitic sandstone 8 m thick and below by the Eau Claire Formation a shale 30 m thick. Initial hydraulic testing with inflatable packers indicated that the aquifer has four hydraulic zones with distinctly different values of relative horizontal hydraulic conductivity. The thickness of each zone was determined by correlating data from geophysical logs, core samples, and the inflatable-packer tests. A comprehensive network for data collection, storage, and analysis was designed to monitor temperature and pressure changes during the ATES test cycles. A total of 22 pressure transducers and 56 thermocouples monitored pressures and temperatures in the aquifer and in the upper and lower confining units. Temperature and pressure measurements were collected in observation well nests at distances of approximately 7 and 14 m from the production wells. All pressure and temperature data were transmitted through buried cables to a central data logger, where the measurements were viewed independently or stored on computer magnetic tape for later analysis. Interactive computer programs were available to display data stored on magnetic tapes as individual measurements or as plots of pressure and temperature versus time. Analyses of step-drawdown and constant-discharge aquifer tests indicate that the Franconia-Ironton-Galesville aquifer is anisotropic in the horizontal plane. Major and minor transmissivities are 101.5 and 44.6 m 2 /d (square meters per day), respectively, for the Ironton and Galesville Sandstones and 40.0 and 24.0 m 2 /d, respectively, for the upper part of the Franconia Formation. The average transmissivity of the entire Franconia-Ironton-Galesville aquifer is about 98 m 2 /d. Effective porosity ranges from 0.25 to 0.31, and the average storage coefficient is 4.5x10 -5 . Two computer models were constructed to simulate the movement of ground water and heat. The first was a nonisothermal, isotropic, single-phase, radial, ground-water flow and thermal-energy-transport model that was constructed to examine the sensitivity of model results to various hydraulic and thermal properties. The model also was used to study the potential for buoyancy flow within the aquifer and the effect of various cyclic injection and withdrawal schemes on the relative thermal efficiency of the aquifer. The second model was a threedimensional ground-water flow and thermal-energy-transport model that was constructed to incorporate the anisotropy of the aquifer. In the first model, the sensitivity analysis assumed 8 days of injection of 150°C water at 18.9 liters per second (L/s), 8 days of storage, and 8 days of withdrawal of hot water at 18.9 L/s. The analysis indicates that, for practical ranges of hydraulic and thermal properties, the ratio of horizontal to vertical hydraulic conductivity is the least important property and thermal dispersivity is the most important property used to compute temperature and aquifer thermal efficiency. Buoyancy flow was examined for several values of hydraulic conductivity and ratios of horizontal to vertical hydraulic conductivities. For the assumed base values of hydraulic and thermal properties, buoyancy flow was negligible. The greatest simulated buoyancy flow resulted from simulations in which horizontal hydraulic conductivity was increased to ten times the base value and in which the vertical hydraulic conductivity was set equal to the horizontal hydraulic conductivity. The effects of various injection and withdrawal rates and durations on computed values of aquifer relative-thermal efficiency and final well-bore temperature were studied for five 1-year hypothetical test cycles of injection and withdrawal. The least efficient scheme was 8 months injection of 150°C water at 18.9 L/s and 4 months of withdrawal of hot water at 18.9 L/s. The most efficient scheme was obtained with 6 months of injection of 150°C water at 18.9 L/s and 6 months of withdrawal of hot water at 37.8 L/s. The hypothetical simulations indicate that the calibrated model of the doublet-well system would be a valuable tool for use by the university in selecting a highly efficient system operation. In the second model, analytical solutions of anisotropic hydraulic flow around the doublet-well system were obtained to provide fluid-flux boundary conditions around the heat-injection well in three dimensions. This information simplified simulation of the doublet-well system because only the heat injection well needed to be simulated. This second model was calibrated with data from an 8-day ambienttemperature injection test at 18.9 L/s. Boundary-flux conditions were examined for nonisothermal conditions by simulating 8 days of injection of 150°C water at 18.9 L/s. Results of simulations using both models indicate that the fluxboundary conditions are adequate for simulations of short-term heatinjection testing.

Professional Paper↗

USGS2026

Small creatures, big challenge: A review of current knowledge and assessment of research needs to inform water withdrawal management for the Duck River's freshwater mussel fauna

The Duck River is a remarkable national resource, boasting outstanding biodiversity (Ahlstedt et al. 2017). It contains 151 fish species, over 60 mussel species, and 22 snail species, many of which are federally listed (Etnier and Starnes 1993; Johnson et al. 2013; Ahlstedt et al. 2017; Womble and Rosenberger 2025a). The Duck River watershed is at a pivotal time in its management and ecological history due to significant growth in human population and economic activity. The consequent increase in surface water demand presents challenges for sustainable water resource management. The Tennessee Department of Environment and Conservation (TDEC) has permitted or proposed to permit eight water utilities to withdraw over 77 million gallons per day (mgd) from the Duck River, representing > 40% increase from previously permitted amounts (USFWS 2025). On November 20, 2024, Tennessee Governor Bill Lee signed Executive Order No 108 to ensure long-term stability of the Duck River watershed through balancing economic growth, water resource management, and environmental conservation. Some primary components of this plan are the creation of a Comprehensive Regional Drought Management Plan, the development of a Habitat Conservation Plan, and the establishment of the Duck River Watershed Planning Partnership (DRWPP). Although the planning partnership is not focused on a sole aspect of the Duck River, but the watershed as a whole, the mussel assemblages of the Duck River have been a particular focus due to their importance in providing ecosystem services, sensitivity to water quality, conservation status, and sensitivity to periodic drought conditions in this system. With multiple state and federal agencies involved in managing the Duck River’s unique water resources, a focused review on how future water management could affect the freshwater mussel communities in the Duck River will contribute to this process. Our objectives of this review are to: 1. Review the effects of changing flow regimes on mussels, with emphasis on the effect of droughts/low flows, highlighting key processes/factors affecting mussels. 2. Highlight critical data gaps that could allow future, evidence-based water management decisions in this system. 3. Propose research topics to address critical data gaps and enable evidence-based water management decisions that consider the needs of freshwater mussels. Freshwater mussels are a critical component of riverine ecosystems, providing ecosystem services including water filtration, nutrient cycling, and habitat provision (Vaughn et al. 2004; Haag 2012) that benefit human communities. Environmental flows, also referred to as ecological or e-flows, refer to the quantity, timing, and quality of water flows necessary to maintain the ecological health of rivers, wetlands, and other aquatic ecosystems (Acreman and Dunbar 2004; Acreman 2016). For freshwater mussels, environmental flows include maintaining a flow regime and associated water quality regimes that are favorable for both juveniles and adults while also seasonally supporting host fish and the survival of settling glochidia (i.e., larval mussels) in suitable environments (Gates et al. 2015). Freshwater mussels are adapted to the natural flow regime of their native streams, including seasonal variability in flows. As sedentary filter feeders, mussels require adequate water flow to provide oxygen and food (Allen and Vaughn 2009); however, they also require periods of low flow to increase access to fish hosts and to ensure juveniles are not displaced (Gates et al 2015). Mussels are sensitive to water quality conditions, including pollutants, sedimentation, and artificial alterations of temperature and pH (Strayer and Malcom 2012). Mussel community characteristics, including presence, abundance, and species richness, generally increase with stream size (e.g., stream order); however, responses to changes in flow conditions are often species-specific, and the magnitude and timing of extreme flow conditions have more predictive power for mussel community response than general changes in flow (Lopez and Vaughn 2021). Complex hydraulic variables, such as relative shear stress and Reynolds number, can be strong predictors of mussel community responses because they combine simple hydraulic variables (e.g., velocity, depth) with stream channel characteristics, such as substrate composition and dynamics. Droughts directly and indirectly affect mussels by influencing their physiology, behavior, reproduction, and ecosystem roles. Mass mussel die-offs during droughts and low water conditions are widely documented (Vaughn et al. 2015; Sousa et al. 2018; Dubose et al. 2019). Droughts can also lead to a shift in mussel community composition, with drought-tolerant species becoming more dominant while more sensitive species experience declines (Galbraith et al, 2010; Lopez et al. 2022; Vaughn and Atkinson 2025). Droughts affect mussels through direct mortality from stranding as well as physiological stress from altered water temperature and quality, reduced habitat and food availability and quality, increased predation and competition, disruption of reproductive cycles, and increased susceptibility to disease and pathogens. Although the effects of droughts, water withdrawals, and low flow conditions on mussels are an increasing area of study, there is a general paucity of data-driven consensus regarding optimal mitigation approaches and strategies (Cushway et al. 2024). Long-term management decisions regarding water withdrawals, particularly during extended or multi-year drought conditions, may consider factors beyond direct mortality such as the sublethal effects of drought and low flow conditions on freshwater mussels for a comprehensive understanding of the effects of prolonged low flows on freshwater mussels and aquatic communities.

Tennessee↗