USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Research”

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

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

1,762 records · Page 98Linked to original sources

Landscape complexity and soil moisture variation in south Georgia, USA, for remote sensing applications

This research addressed the temporal and spatial variation of soil moisture (SM) in a heterogeneous landscape. The research objective was to investigate soil moisture variation in eight homogeneous 30 by 30 m plots, similar to the pixel size of a Landsat Thematic Mapper (TM) or Enhanced Thematic Mapper plus (ETM+) image. The plots were adjacent to eight stations of an in situ soil moisture network operated by the United States Department of Agriculture-Agriculture Research Service USDA-ARS in Tifton, GA. We also studied five adjacent agricultural fields to examine the effect of different landuses/land covers (LULC) (grass, orchard, peanuts, cotton and bare soil) on the temporal and spatial variation of soil moisture. Soil moisture field data were collected on eight occasions throughout 2005 and January 2006 to establish comparisons within and among eight homogeneous plots. Consistently throughout time, analysis of variance (ANOVA) showed high variation in the soil moisture behavior among the plots and high homogeneity in the soil moisture behavior within them. A precipitation analysis for the eight sampling dates throughout the year 2005 showed similar rainfall conditions for the eight study plots. Therefore, soil moisture variation among locations was explained by in situ local conditions. Temporal stability geostatistical analysis showed that soil moisture has high temporal stability within the small plots and that a single point reading can be used to monitor soil moisture status for the plot within a maximum 3% volume/volume (v/v) soil moisture variation. Similarly, t-statistic analysis showed that soil moisture status in the upper soil layer changes within 24 h. We found statistical differences in the soil moisture between the different LULC in the agricultural fields as well as statistical differences between these fields and the adjacent 30 by 30 m plots. From this analysis, it was demonstrated that spatial proximity is not enough to produce similar soil moisture, since t-test's among adjacent plots with different LULCs showed significant differences. These results confirm that a remote sensing approach that considers homogeneous LULC landscape fragments can be used to identify landscape units of similar soil moisture behavior under heterogeneous landscapes. In addition, the in situ USDA-ARS network will serve better in remote sensing studies in which sensors with fine spatial resolution are evaluated. This study is a first step towards identifying landscape units that can be monitored using the single point reading of the USDA-ARS stations network. ?? 2008 Elsevier B.V.

Journal of Hydrology↗

A landscape indicator approach to the identification and articulation of the consequences of land-cover change in the Mid-Atlantic Region, 1973-2001

Landscape indicators, derived from land-use and land-cover data, hydrology, nitrate deposition, and elevation data, were used by Jones and others (2001a) to calculate the ecological consequences of land-cover change. Nitrate loading and physical bird habitat were modeled from 1973 and 1992 land-cover and other spatial data for the Mid-Atlantic region. Utilizing the same methods, this study extends the analysis another decade with the use of the 2001 National Land Cover Dataset. Land-cover statistics and trends are calculated for three time periods: 1973-1992, 1992-2001 and 1973-2001. In addition, high-resolution aerial photographs (1 meter or better ground-sample distance) were acquired and analyzed for thirteen pairs of adjacent USGS 7.5 minute quadrangle maps in areas where distinct positive or negative changes to nitrogen loading and bird habitat were previously calculated. During the entire 30 year period, the data show that there was extensive loss of agriculture and forest area and a major increase in urban land-cover classes. However, the majority of the conversion of other classes to urban occurred during the 1992-2001 period. During the 1973-1992 period, there was only moderate increase in urban area, while there was an inverse relationship between agricultural change and forest change. In general, forest gain and agricultural loss was found in areas of improving landscape indicators, and forest loss and agricultural gain was found to occur in areas of declining indicators related to habitat and nitrogen loadings, which was generally confirmed by the aerial photographic analysis. In terms of the specific model results, bird habitat, which is mainly related to the extent of forest cover, declined overall with forest extent, but was also affected more in the decline of habitat quality. Nitrate loading, which is mainly related to agricultural land cover actually improved from 1992-2001, and in the overall study, mainly due to the conversion of agriculture to forests and urban. The high-resolution imagery analysis was significant in that it confirmed, at a very local level, the specific land-cover changes that were driving the landscape metrics and model results that were calculated from moderate resolution land-cover data and models. These were generally subtle changes in patch size of agriculture, forest, and urban areas, but had substantial effects on bird habitat and nitrogen loadings. This analysis of high-resolution imagery demonstrates and confirms the important ability of moderate-resolution land-cover data to capture significant landscape-level activity that is directly related to specific metrics of ecological significance. It also demonstrates consistent landscape-scale relationships between data derived from high-resolution, moderate-resolution and landscape-model sources. Finally, many of the areas of improvement and decline in bird habitat and nitrogen loadings appear to be potentially regional in nature and likely reflect some local trend in landscape activity. Although the use of ecoregions as sampling units has been criticized in recent years, these results show that basic changes in Level 1 land-cover categories, such as forest and agriculture, may still reflect ecoregional patterns and considerations at some scale of mapping and analysis. This is a potentially important area for future landscape-indicator research. This and other follow-on research opportunities are discussed.

Open-File Report↗

Comprehensive water quality of the Boulder Creek Watershed, Colorado, during high-flow and low-flow conditions, 2000

Executive Summary The Boulder Creek Watershed, Colorado, is 1160 square kilometers in area and ranges in elevation from 1480 to 4120 meters above sea level. Streamflow originates primarily as snowmelt near the Continental Divide, and thus discharge varies seasonally and annually (Chapter 1). Most of the water in Boulder Creek is diverted for domestic, agricultural, and industrial use. Some diverted water is returned to the creek as wastewater effluent and by ditch returns, and additional water enters as groundwater and by transbasin diversions. These diversions and returns lead to complex temporal and spatial variations in discharge. The variations in discharge, along with natural factors such as geology and climate, and anthropogenic factors such as wastewater treatment, agriculture, mining, and urbanization, can affect water chemistry. As with many watersheds in the American West, dependable water quality and sufficient water supply are issues facing local water managers and users. Detailed water-quality and sediment sampling allows the identification of sources and sinks of chemical constituents and an understanding of the processes at work in a river system. This study, the most comprehensive water-quality analysis performed for Boulder Creek to date, was a cooperative effort of the U.S. Geological Survey (USGS) and the city of Boulder. Geographic information systems and modeling programs were used to delineate watershed boundaries, land cover, and geology (Chapter 2). During high-flow (June 2000) and low-flow (October 2000) conditions, researchers evaluated 226 water-quality variables, including basic water-quality indicators (Chapter 3), major ions and trace elements (Chapter 4), wastewater-derived organic compounds (Chapter 5), and pesticides (Chapter 6). Discharge (Chapter 1) and bed-sediment particle size and mineralogy (Chapter 7) were also evaluated. This cooperative study was facilitated by the Boulder Area Sustainability Information Network (BASIN), which provides public access to environmental information about the Boulder Creek Watershed on a website, www.basin.org. In addition to the USGS and city of Boulder data, researchers at the Institute of Arctic and Alpine Research at the University of Colorado provided water chemistry data for the headwaters of North Boulder Creek, upstream of the reach of the USGS/city of Boulder sampling sites (Chapter 8). Snowmelt produces high flows in Boulder Creek in late spring to early summer (Chapter 1). Because precipitation falling in the headwaters is very dilute (specific conductance about 5 microsiemens per centimeter), most chemical constituents are present in lower concentrations during high flows (Chapters 3, 4, 5, 6, and 8). However, concentrations of some constituents, such as total suspended solids (Chapter 3) and organic carbon (Chapter 5), increase during the spring snowmelt flush. The upper basin, which consists of alpine, subalpine, montane, and foothills regions west of the mouth of Boulder Canyon, is underlain by Precambrian igneous and metamorphic rocks (Chapter 1). Major dissolved inorganic constituents in headwater sites were found to be enriched by factors of 10 to 20 relative to precipitation; this is consistent with minor weathering of the local crystalline bedrock (Chapter 4). Some anthropogenic input is observed in the headwaters; precipitation introduces nitrogen derived from fossil fuel combustion and agricultural activities (Chapter 8). The lower basin, which consists of the plains region east of the mouth of Boulder Canyon, is underlain by Mesozoic sedimentary rock and Quaternary alluvium, and has substantially more anthropogenic sources. Concentrations of most dissolved inorganic constituents increased in the lower basin. Differentiation between natural and anthropogenic sources of some dissolved constituents is difficult because both sources contribute to the water composition in this region. The increase of most major constituents (bicarbonate, calcium, chloride, magnesium, sodium, and sulfate) is consistent with weathering of the underlying sedimentary bedrock (Chapter 4). It is likely that anthropogenic loading of constituents in this reach occurs during storm events. Fecal coliform concentrations were variable and in some cases exceeded state standards, primarily during low-flow conditions (Chapter 3). Effluent from Boulder’s 75th Street Wastewater Treatment Plant (WWTP) has a substantial impact on the water chemistry of lower Boulder Creek. The WWTP increases the concentrations of nutrients such as nitrogen and phosphorus (Chapter 3), major ions and trace metals (Chapter 4), and organic carbon (Chapter 5) in Boulder Creek. The effluent contained a spike in gadolinium, a rare earth element that is ingested for magnetic resonance imaging as a contrasting agent and then excreted to the urban wastewater system. The effluent also contained trace organic compounds such as surfactants, pharmaceuticals, hormones (Chapter 5), and pesticides (Chapter 6), which also were detected at downstream Boulder Creek sites. Water chemistry of Boulder Creek downstream of the WWTP is largely controlled by the degree of dilution of the wastewater effluent, which varies depending on the baseflow of Boulder Creek, the volume of wastewater effluent, and depletion by agricultural diversions. Coal Creek, a tributary of Boulder Creek, contains wastewater effluent from four additional WWTPs, and increases the load of many constituents in Boulder Creek. In addition to the impact from wastewater effluent, lower Boulder Creek is affected by agricultural land use. Eleven of 84 analyzed pesticides were detected in Boulder Creek or its inflows, primarily in the eastern section of the watershed (Chapter 6). This collaborative study provides an in-depth evaluation of the hydrology, water chemistry, and sediment mineralogy of North Boulder Creek, Middle Boulder Creek, Boulder Creek, and major inflows. The detailed sampling and analysis in this report provide a baseline for future reference, as well as information on the effect of land use and geology on water chemistry.

Colorado↗

Applied Geochemistry Special Issue on Environmental geochemistry of modern mining

Environmental geochemistry is an integral part of the mine-life cycle, particularly for modern mining. The critical importance of environmental geochemistry begins with pre-mining baseline characterization and the assessment of environmental risks related to mining, continues through active mining especially in water and waste management practices, and culminates in mine closure. The enhanced significance of environmental geochemistry to modern mining has arisen from an increased knowledge of the impacts that historical and active mining can have on the environment, and from new regulations meant to guard against these impacts. New regulations are commonly motivated by advances in the scientific understanding of the environmental impacts of past mining. The impacts can be physical, chemical, and biological in nature. The physical challenges typically fall within the purview of engineers, whereas the chemical and biological challenges typically require a multidisciplinary array of expertise including geologists, geochemists, hydrologists, microbiologists, and biologists. The modern mine-permitting process throughout most of the world now requires that potential risks be assessed prior to the start of mining. The strategies for this risk assessment include a thorough characterization of pre-mining baseline conditions and the identification of risks specifically related to the manner in which the ore will be mined and processed, how water and waste products will be managed, and what the final configuration of the post-mining landscape will be. In the Fall 2010, the Society of Economic Geologists held a short course in conjunction with the annual meeting of the Geological Society of America in Denver, Colorado (USA) to examine the environmental geochemistry of modern mining. The intent was to focus on issues that are pertinent to current and future mines, as opposed to abandoned mines, which have been the focus of numerous previous short courses. The geochemical challenges of current and future mines share similarities with abandoned mines, but differences also exist. Mining and ore processing techniques have changed; the environmental footprint of waste materials has changed; environmental protection has become a more integral part of the mine planning process; and most historical mining was done with limited regard for the environment. The 17 papers in this special issue evolved from the Society of Economic Geologists’ short course. The relevant geochemical processes encompass the source, transport, and fate of contaminants related to the life cycle of a mine. Contaminants include metals and other inorganic species derived from geologic sources such as ore and solid mine waste, and substances brought to the site for ore processing, such as cyanide to leach gold. Factors, such as mine-waste mineralogy, hydrologic setting, mine-drainage chemistry, and microbial activity, that affect the hydrochemical risks from mining are reviewed by Nordstrom et al. In another paper, Nordstrom discusses baseline characterization at mine sites in a regulatory framework, and emphasizes the influence of mineral deposits in producing naturally elevated concentrations of many trace elements in surface water and groundwater. Surface water quality in mineralized watersheds is influenced by a number of processes that act on daily (diel) cycles and can produce dramatic variations in trace element concentrations as described by Gammons et al. Pre-mining baseline characterization studies should strive to capture the magnitude of these diel variations. Desbarats et al., using a case study of mine drainage from a gold mine, illustrate how elements that commonly occur as negatively charged species (anions) in solution, such as arsenic as arsenate, behave in an opposite fashion than most metals, which occur as positively charged species (cations). Significant improvement in the understanding of factors that influence the toxicity of metals to aquatic organisms in surface water has highlighted the importance of aqueous chemistry, particularly dissolved organic carbon, as described by Smith et al. Stream sediment contamination is another important pathway for affecting aquatic organisms, as reviewed by Besser et al. Understanding and predicting environmental consequences from mining begins with knowing the mineralogy and mineral reactivity of the ore, the wastes, and of secondary minerals formed later. Jamieson et al. review the importance of mineralogical studies in mine planning and remediation. A number of types of site-specific studies are needed to identify environmental risks related to individual mines. Lapakko reviews the general framework of mine waste characterization studies that are integral to the mine planning process. Hageman et al. present a comparative study of several static tests commonly used to characterize mine waste. The mining and ore processing practices employed at a specific mine site will vary on the basis of the commodities being targeted, the geology of the deposit, the geometry of the deposit, and the mining and ore processing methods used. Thus, these factors, in addition to the waste management practices used, can result in a variety of end-member mine waste features, each of which has its own set of challenges. Open pit mines and underground mines require waste rock to be removed to access ore. Waste rock presents unique problems because the rock is commonly mineralized at sub-economic grades and has not been processed to remove potentially problematic minerals, such as pyrite. Amos et al. examine the salient aspects of the geochemistry of waste rock. Mill tailings – the waste material after ore minerals have been removed – are a volumetrically important solid waste at many mine sites. Their fine grain size and the options for their management make their behavior in the environment distinct from that of waste rock. Lindsay et al. describe some of these differences through three case-study examples. Subaqueous disposal of tailings is another option described by Moncur et al. Cyanide leaching for gold extraction is a common method throughout the world. Johnson describes environmental aspects of cyanidation. Uranium mining presents unique environmental challenges, particularly since in-situ recovery has seen widespread use. Campbell et al. review the environmental geochemistry of uranium mining and current research on bioremediation. Ore concentrates from many types of metal mining undergo a pyrometallurgical technique known as smelting to extract the metal. Slag is the result of smelting, and it may be an environmental liability or a valuable byproduct, as described by Piatak et al. Finally, the open pits that result from surface mining commonly reach below the water table. At the end of mining, these pits may fill to form lakes that become part of the legacy of the mine. Castendyk et al., in two papers, review theoretical aspects of the environmental limnology of pit lakes. They also describe approaches that have been used to model pit lake water balance, wall-rock contributions to pit lake chemistry, pit lake water quality, and limnological processes, such as vertical mixing, through the use of three case studies.

Applied Geochemistry↗

The lacustrine carbon cycle as illuminated by the waters and sediments of two hydrologically distinct headwater lakes in north-central Minnesota, U.S.A

The accumulation rates of CaCO 3 and organic carbon (OC) in lake sediments are delicately balanced between production in the epilimnion and destruction in the hypolimnion. The cycling of these two forms of carbon makes a "carbon pump" that greatly affects the biogeochemical cycles of other elements. To further understand these biogeochemical dynamics, the lakes, streams, and wetlands of the Shingobee River headwater area of north-central Minnesota have been subjected to intensive hydrologic and biogeochemical studies. Williams Lake, situated close to the highest point in the regional flow system, is hydrologically closed, with no surface inlet or outlet, and ground water and precipitation as the only sources of water. Shingobee Lake, situated at the lowest point in the regional flow system, has the Shingobee River as an inlet and outlet. The surface waters of both lakes are oversaturated, and the bottom waters undersaturated, with respect to CaCO 3 during the summer. The small amount of CaCO 3 that is precipitated in the epilimnion of Williams Lake during the summer is dissolved in the undersaturated hypolimnion and sediments with the result that no CaCO 3 is incorporated into the profundal surface sediments. Because of the high phytoplankton productivity of Shingobee Lake, sufficient CaCO 3 is produced in the epilimnion that large amounts survive the corrosive hypolimnion and sediments, and an average of 46 wt. % accumulates in surface sediments. Another consequence of higher phytoplankton productivity in Shingobee Lake is that the hypolimnion becomes oxygen deficient within a month after overturn in both the spring and fall. Because of reducing conditions that develop in the hypolimnion of Shingobee Lake, high concentrations of dissolved Fe and Mn accumulate there during summer stratification. Precipitation of Fe and Mn oxyhydroxides during periods of fall and spring overturn results in high concentrations of Fe and Mn in surface sediments. In Williams Lake, high concentrations of Fe and Mn do not build up in the hypolimnion. The concentration of CaCO 3 is about 80 wt. % in lower Holocene sediments of both lakes. The lower Holocene sediments in both lakes also contain high concentrations of Fe and Mn, and the lower Holocene sediments of Shingobee are laminated. The waters of both lakes had identical values of δ 13 C and δ 18 O during the early Holocene, but the waters of Williams Lake "evolved" during the early Holocene, increasing about 10‰ in both δ 13 C and δ 18 O. Deposits of lacustrine marl occur as much as seven meters above the present elevation of Williams Lake, the highest of the two lakes. Taken together, these observations suggest that the lakes were once connected to form a larger lake called Lake Willobee with a hypolimnion that was anoxic, at least seasonally.

Minnesota↗

Surface monitoring of microseismicity at the Decatur, Illinois, CO2 sequestration demonstration site

Sequestration of CO2 into subsurface reservoirs can play an important role in limiting future emission of CO2 into the atmosphere (e.g., Benson and Cole, 2008). For geologic sequestration to become a viable option to reduce greenhouse gas emissions, large-volume injection of supercritical CO2 into deep sedimentary formations is required. These formations offer large pore volumes and good pore connectivity and are abundant (Bachu, 2003; U.S. Geological Survey Geologic Carbon Dioxide Storage Resources Assessment Team, 2013). However, hazards associated with injection of CO2 into deep formations require evaluation before widespread sequestration can be adopted safely (Zoback and Gorelick, 2012). One of these hazards is the potential to induce seismicity on pre-existing faults or fractures. If these faults or fractures are large and critically stressed, seismic events can occur with magnitudes large enough to pose a hazard to surface installations and, possibly more critical, the seal integrity of the cap rock. The Decatur, Illinois, carbon capture and storage (CCS) demonstration site is the first, and to date, only CCS project in the United States that injects a large volume of supercritical CO2 into a regionally extensive, undisturbed saline formation. The first phase of the Decatur CCS project was completed in November 2014 after injecting a million metric tons of supercritical CO2 over three years. This phase was led by the Illinois State Geological Survey (ISGS) and included seismic monitoring using deep borehole sensors, with a few sensors installed within the injection horizon. Although the deep borehole network provides a more comprehensive seismic catalog than is presented in this paper, these deep data are not publicly available. We contend that for monitoring induced microseismicity as a possible seismic hazard and to elucidate the general patterns of microseismicity, the U.S. Geological Survey (USGS) surface and shallow borehole network described below provides an adequate event detection threshold. The formation targeted for injection is the Mount Simon Sandstone, which is laterally extensive, has high porosity and permeability and has the potential to host future CCS projects due to its favorable hydrologic characteristics and proximity to industrial sources of CO2 (Birkholzer and Zhou, 2009). At Decatur, CO2, a byproduct of ethanol production at the Archer Daniels Midland (ADM) facility, is compressed to supercritical state and injected at 2.1 km depth into the 460 m thick Mount Simon Sandstone. This sandstone has varying properties, ranging from the lower, fine- to coarse-grained sandstone with high permeability and porosity, to the middle and upper Mount Simon, which consist of planar, cross-bedded layers of varied permeability and porosity (Leetaru and Freiburg, 2014). The changes in permeability and porosity within the Mount Simon Sandstone, due to depositional and diagenetic differences, create horizontal baffles, which inhibit vertical flow and restrict the injected CO2 to remain near the injection horizon (Bowen et al., 2011). The lowest portion of the Mount Simon Sandstone overlying the Precambrian rhyolite basement is the Pre-Mount Simon interval, generally < 15 m in thickness and composed of fine- to medium-grain size sandstone that is highly deformed (Leetaru and Freiburg, 2014). The basement rhyolite has a clayrich matrix and is fractured, with significant alterations within the fractures. The primary sealing cap rock is the Eau Claire Formation, a 100–150 m thick unit at a depth of roughly 1.69 km (Leetaru and Freiburg, 2014). The Maquoketa Shale Group and the New Albany Shale serve as secondary and tertiary seals at shallower depths of ∼820 and ∼650 m, respectively. The ISGS managed the Illinois Basin–Decatur Project (IBDP), a three-year project beginning in November 2011, during which carbon dioxide was injected at a rate of ∼1000 metric tons per day until November 2014 (Finley et al., 2011, 2013). ADM manages the Illinois Industrial CCS (ICCS) project, which will inject ∼3000 metric tons/day into a second injection well starting in the summer of 2015. The USGS began monitoring microseismicity with a 13- station seismic network at Decatur in July 2013 (Fig. 1). This network provides good detection capabilities and azimuthal (focal sphere) coverage for microseismicity with moment magnitudes (Mw) above about −0:5. Here, we report on 19 months of microseismicity monitoring at the Decatur CO2 sequestration site, which permits a detailed look at the evolution and character of injection-induced seismicity.

Illinois↗

Effect of subalpine canopy removal on snowpack, soil solution, and nutrient export, Fraser Experimental Forest, CO

Research on the effects of vegetation manipulation on snowpack, soil water, and streamwater chemistry and flux has been underway at the Fraser Experimental Forest (FEF), CO, since 1982. Greater than 95% of FEF snowmelt passes through watersheds as subsurface flow where soil processes significantly alter meltwater chemistry. To better understand the mechanisms accounting for annual variation in watershed streamwater ion concentration and flux with snowmelt, we studied subsurface water flow, its ion concentration, and flux in conterminous forested and clear cut plots. Repetitive patterns in subsurface flow and chemistry were apparent. Control plot subsurface flow chemistry had the highest ion concentrations in late winter and fall. When shallow subsurface flow occurred, its Ca 2+ , SO 4 2− , and HCO 3 − concentrations were lower and K + higher than deep flow. The percentage of Ca 2+ , NO 3 − , SO 4 2− , and HCO 3 − flux in shallow depths was less and K + slightly greater than the percentage of total flow. Canopy removal increased precipitation reaching the forest floor by about 40%, increased peak snowpack water equivalent (SWE)>35%, increased the average snowpack Ca 2+ , NO 3 − , and NH 4 + content, reduced the snowpack K + content, and increased the runoff four-fold. Clear cutting doubled the percentage of subsurface flow at shallow depths, and increased K + concentration in shallow subsurface flow and NO 3 − concentrations in both shallow and deep flow. The percentage change in total Ca 2+ , SO 4 2− , and HCO 3 − flux in shallow depths was less than the change in water flux, while that of K + and NO 3 − flux was greater. Relative to the control, in the clear cut the percentage of total Ca 2+ flux at shallow depths increased from 5 to 12%, SO 4 2− 5·4 to 12%, HCO 3 − from 5·6 to 8·7%, K + from 6 to 35%, and NO 3 − from 2·7 to 17%. The increases in Ca 2+ and SO 4 2− flux were proportional to the increase in water flux, the flux of HCO 3 − increased proportionally less than water flux, and NO 3 − and K + were proportionally greater than water flux. Increased subsurface flow accounted for most of the increase in non-limiting nutrient loss. For limiting nutrients, loss of plant uptake and increased shallow subsurface flow accounted for the greater loss. Seasonal ion concentration patterns in streamwater and subsurface flow were similar .

Hydrological Processes↗

Hydroecological factors governing surface water flow on a low-gradient floodplain

"Interrelationships between hydrology and aquatic ecosystems are better understood in streams and rivers compared to their surrounding floodplains. Our goal was to characterize the hydrology of the Everglades ridge and slough floodplain ecosystem, which is valued for the comparatively high biodiversity and connectivity of its parallel-drainage features but which has been degraded over the past century in response to flow reductions associated with flood control. We measured flow velocity, water depth, and wind velocity continuously for 3 years in an area of the Everglades with well-preserved parallel-drainage features (i.e., 200-m wide sloughs interspersed with slightly higher elevation and more densely vegetated ridges). Mean daily flow velocity averaged 0.32 cm s1 and ranged between 0.02 and 0.79 cm s1. Highest sustained velocities were associated with flow pulses caused by water releases from upstream hydraulic control structures that increased flow velocity by a factor of 2–3 on the floodplain for weeks at a time. The highest instantaneous measurements of flow velocity were associated with the passage of Hurricane Wilma in 2005 when the inverse barometric pressure effect increased flow velocity up to 5 cm s1 for several hours. Time-averaged flow velocities were 29% greater in sloughs compared to ridges because of marginally higher vegetative drag in ridges compared to sloughs, which contributed modestly (relative to greater water depth and flow duration in sloughs compared to ridges) to the predominant fraction (86%) of total discharge through the landscape occurring in sloughs. Univariate scaling relationships developed from theory of flow through vegetation, and our field data indicated that flow velocity increases with the square of water surface slope and the fourth power of stem diameter, decreases in direct proportion with increasing frontal area of vegetation, and is unrelated to water depth except for the influence that water depth has in controlling the submergence height of vegetation that varies vertically in its architectural characteristics. In the Everglades the result of interactions among controlling variables was that flow velocity was dominantly controlled by water surface slope variations responding to flow pulses more than spatial variation in vegetation characteristics or fluctuating water depth. Our findings indicate that floodplain managers could, in addition to managing water depth, manipulate the frequency and duration of inflow pulses to manage water surface slope, which would add further control over flow velocities, water residence times, sediment settling, biogeochemical transformations, and other processes that are important to floodplain function."

Water Resources Research↗

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗

Fate and pathways of injection-well effluent in the Florida Keys

Twenty-four wells (21 locations) were core drilled into the limestone beneath the Keys, reef tract, and outer reefs to determine if sewage effluents injected in Class V wells onshore are reaching offshore reef areas via underground flow. These wells were fitted with PVC casings and well screens and were sampled every three months for a period of one year. Analyses showed consistent hypersalinity in most wells and a marked increase in nitrogen (as ammonia) in offshore ground water. Other forms of nitrogen (NO 2 and NO 3 ) and phosphorous were not particularly elevated in offshore ground water but were above the levels found in surface marine water. The highest levels of nitrogen (NO 2 and NO 3 ) and phosphorous were in shallow onshore ground waters. Sources for the nutrients in the shallow onshore ground water consist of septic tanks and cesspools (@ 24,000 and 5,000 in the Florida Keys, respectively), agricultural fertilizers, and natural vegetation. Ammonia concentrations were low in shallow ground waters beneath the Florida Keys, probably because of oxidizing conditions. Tidal pumping is particularly active, especially nearshore. Hydraulic heads sufficient to elevate well water as much as 7 cm above sea level during falling tides were detected in all nearshore wells. During rising tides, the situation was reversed and water flowed into the wells. Tidal pumping implies considerable water movement both in and out of the upper few meters of limestone. Tidal pumping is a likely mechanism for mixing and transferring nutrient-rich ground water into the overlying marine waters. Although tidal pumping should cause rather complete mixing and dilution of any freshwater-based effluents entering the limestone via the more than 600 disposal wells in the Florida Keys, the ground waters in the 30- to 40-ft-depth range (9-12 m) nevertheless remained slightly hypersaline relative to sea water throughout the year. Fecal coliform and fecal streptococcal bacteria were associated with three Lower Keys offshore wells and two shallow onshore wells at Key Largo. On occasions, these bacteria were detected farther offshore, once in a well 4 miles off Key Largo. The bacterial analyses for Key Largo (both onshore and offshore) are supported by two independent bacteriological researchers using more sophisticated methods than the standard 100-ml membrane-filter method used in this study. Fecal bacteria can serve as tracers; thus, we conclude their presence is possible evidence for offshore transport of ground waters originating on Key Largo. Elevated nutrients (ammonia) and slightly elevated dissolved total phosphorous in offshore ground waters, however, cannot be tied to onshore sources with existing data. Rock analyses of material from our cores do not prove or disprove the hypothesis that limestone beneath the Keys or reef tract is serving as a sink for phosphorus or other nutrients. The data, however, do not rule out phosphorus uptake by limestone adjacent to disposal sources. For the purposes of this study, monitoring wells were not positioned sufficiently close to injection wells to determine if uptake of phosphorous is taking place. Ground waters were found to contain more dissolved solids than could be accounted for if hypersalinity resulted from simple evaporation of sea water. These data indicate that ground waters in the vicinity of our wells are dissolving solids from the rock rather than precipitating material within the rock framework; however, as mentioned above, our wells were not positioned sufficiently close to disposal wells to determine if localized uptake is occurring. Examination of rock cores from these wells revealed a general distribution of reef- and grainstone-facies belts. The Upper and Middle Keys are composed of a thin coral reef facies that extends only a few hundred feet seaward of the Keys. Reef facies give way to mudstone facies within a few yards of shore on the Florida Bay side of the Keys. On the seaward side of the Keys, beneath Hawk Channel and White Bank, the Pleistocene limestone is a mixed grainstone, packstone, and wackstone facies. Corals are rare or absent. The Pleistocene limestone beneath the outer reefs 4 to 5 miles offshore, however, consists of reef facies with the same coral fauna as that found on Key Largo. This pattern of two major reef-facies belts separated by a 2- to 4-mile-wide belt of grainstone facies may have as yet undetermined effects on groundwater circulation beneath the Florida reef tract. Grainstone is approximately an order of magnitude less permeable than the coralline Key Largo Limestone facies. The Q3 surface, a major subsurface unconformity thought to form an effective confining zone elsewhere in south Florida, was not detected in wells drilled more than 1 mile from shore. This unconformity, however, was detected in all wells drilled on or near the Keys. What was found to be a more effective and widespread confining layer is the Holocene sediment deposited on the Pleistocene limestone during the past 6,000 to 7,000 years. These relatively impermeable sediments are extensive, forming a belt up to 5 miles wide beginning about 0.5 mile offshore. Holocene sediments generally consist of low-permeability lime mud just above the Pleistocene surface, overlain by more permeable carbonate sands and reefs. Leakage of ground water by tidal pumping is not likely to occur through lime-mud-dominated areas such as Hawk Channel but is likely to occur through isolated porous and permeable Holocene reefs situated on Pleistocene limestone highs, and in places where Holocene sediment does not cover the limestone bedrock. Leakage is therefore limited to 1) a shallow-water 0.5-mile-wide nearshore belt of exposed Key Largo Limestone, 2) Holocene patch reefs, which grow on mud-free topographic rock highs, and 3) along the seaward side of the outermost reef in 35 to 65 ft (10-20 m) of water, where Holocene reef and sediment accumulations are thin or absent. This study did not address direct measurements of lateral groundwater movement or a hydrologic mechanism for transporting hypersaline ground water away from the Florida Keys. More recent work, however (Halley et al., 1994), shows that sea level in Florida Bay is higher than on the Atlantic side of the Keys more than 50% of the time. Higher sea level on the bay side of the Keys provides a potential for groundwater flow toward the Atlantic most of the time. Use of tracers (dyes or harmless bacteriological tracers) injected into the center of tightly spaced clusters of monitoring wells is a simple way to ascertain the net direction and rate of groundwater movement. Knowing the direction and rate of groundwater movement is needed for prediction and modeling efforts in the future

Florida↗

Characterization and evaluation of controls on post-fire streamflow response across western U.S. watersheds

This research investigates the impact of wildfires on watershed flow regimes, specifically focusing on evaluation of fire events within specified hydroclimatic regions in the western United States, and evaluating the impact of climate and geophysical variables on response. Eighty-two watersheds were identified with at least 10 years of continuous pre-fire daily streamflow records and 5 years of continuous post-fire daily flow records. Percent change in annual runoff ratio, low flows, high flows, peak flows, number of zero flow days, baseflow index, and Richards–Baker flashiness index were calculated for each watershed using pre- and post-fire periods. Independent variables were identified for each watershed and fire event, including topographic, vegetation, climate, burn severity, percent area burned, and soils data. Results show that low flows, high flows, and peak flows increase in the first 2 years following a wildfire and decrease over time. Relative response was used to scale response variables with the respective percent area of watershed burned in order to compare regional differences in watershed response. To account for variability in precipitation events, runoff ratio was used to compare runoff directly to PRISM precipitation estimates. To account for regional differences in climate patterns, watersheds were divided into nine regions, or clusters, through k-means clustering using climate data, and regression models were produced for watersheds grouped by total area burned. Watersheds in Cluster 9 (eastern California, western Nevada, Oregon) demonstrate a small negative response to observed flow regimes after fire. Cluster 8 watersheds (coastal California) display the greatest flow responses, typically within the first year following wildfire. Most other watersheds show a positive mean relative response. In addition, simple regression models show low correlation between percent watershed burned and streamflow response, implying that other watershed factors strongly influence response. Spearman correlation identified NDVI, aridity index, percent of a watershed's precipitation that falls as rain, and slope as being positively correlated with post-fire streamflow response. This metric also suggested a negative correlation between response and the soil erodibility factor, watershed area, and percent low burn severity. Regression models identified only moderate burn severity and watershed area as being consistently positively/negatively correlated, respectively, with response. The random forest model identified only slope and percent area burned as significant watershed parameters controlling response. Results will help inform post-fire runoff management decisions by helping to identify expected changes to flow regimes, as well as facilitate parameterization for model application in burned watersheds.

Hydrology and Earth System Sciences↗

Davis Pond freshwater prediversion biomonitoring study: freshwater fisheries and eagles

In January 2001, the construction of the Davis Pond freshwater diversion structure was completed by the U.S. Army Corps of Engineers. The diversion of freshwater from the Mississippi River is intended to mitigate saltwater intrusion from the Gulf of Mexico and to lessen the concomitant loss of wetland areas. In addition to the freshwater inflow, Barataria Bay basin would receive nutrients, increased flows of sediments, and water-borne and sediment-bound compounds. The purpose of this biomonitoring study was, therefore, to serve as a baseline for prediversion concentrations of selected contaminants in bald eagle ( Haliaeetus leucocephalus ) nestlings (hereafter referred to as eaglets), representative freshwater fish, and bivalves. Samples were collected from January through June 2001. Two similarly designed postdiversion studies, as described in the biological monitoring program, are planned. Active bald eagle nests targeted for sampling eaglet blood (n = 6) were generally located southwest and south of the diversion structure. The designated sites for aquatic animal sampling were at Lake Salvador, at Lake Cataouatche, at Bayou Couba, and along the Mississippi River. Aquatic animals representative of eagle prey were collected. Fish were from three different trophic levels and have varying feeding strategies and life histories. These included herbivorous striped mullet ( Mugil cephalus ), omnivorous blue catfish ( Ictalurus furcatus ), and carnivorous largemouth bass ( Micropterus salmoides ). Three individuals per species were collected at each of the four sampling sites. Freshwater Atlantic rangia clams ( Rangia cuneata ) were collected at the downstream marsh sites, and zebra mussels ( Dreissena spp.) were collected on the Mississippi River. The U.S. Geological Survey (USGS) Biomonitoring of Environmental Status and Trends (BEST) protocols served as guides for fish sampling and health assessments. Fish are useful for monitoring aquatic ecosystems because they accumulate pesticides and other contaminants. Biomarker data on individual fish, generated at the USGS National Wetlands Research Center (Lafayette, La.), included percent white blood cells in whole blood, spleen weight to body weight ratio, liver weight to body weight ratio, condition factor, splenic macrophage aggregates, and liver microsomal 7-ethoxyresorufin-o-deethylase (EROD) activity. Fish age was estimated by comparing total lengths with values from the same species in the Southeast United States as determined from the literature. Contaminant analyses were coordinated by the U.S. Fish and Wildlife Service (USFWS) Analytical Control Facility (Laurel, Md.), where residues of organochlorine (OC) pesticides, total polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons (AHs), and trace elements were determined. The organic contaminant data were generated at the Mississippi State University Chemical Lab (Mississippi State, Miss.), and the inorganic contaminant data were generated by the Texas A&M University Geochemical and Environmental Research Group (College Station, Tex.). Statistical tests were performed to assess relationships among contaminants, fish age, fish species, and collection sites. Trends in interspecific differences among fish in concentrations of contaminants were noted. Striped mullet (hereafter mullet) frequently displayed the highest chemical concentrations. Levels of contaminants were generally higher in samples obtained from the Mississippi River than in those collected from the diversion area and were higher in mussels and clams (hereafter bivalves) than in fish. Because the Mississippi River sampling site for mullet and largemouth bass was downriver of the structure and south of New Orleans and the catfish site was upriver, the downriver data may not be directly reflective of the results from the receiving waters at the Davis Pond structure. Compared to the Caernarvon freshwater prediversion study in 1990 that assessed possible influx of contaminants with the freshwater diversion, contaminant levels in fishes and bivalves in this study were generally lower, yet three nontoxic inorganic elements in Davis Pond fish samples exhibited ranges of concentrations that were more than two times higher than did those from Caernarvon. Levels in bivalves were different between diversions but about equal in the numbers of trace elements showing high levels per location. Contaminant values were compared to those listed in various literature and agency sources, both regional and national, including the National Contaminant Biomonitoring Program (NCBP), in which the 85th percentile and above represents what is considered to be an elevated contaminant concentration and cause for concern. Generally, bivalves were at the high end of their ranges for both organic and inorganic contaminants. In this study, OCs were detectable in 67 percent of fish from the Mississippi River site, ranging from 0.15 to 1.09 μg/g wet weight (ww) or fresh weight (fw), and in 11 percent of the fish from the marsh sites, ranging from 0.06 to 0.612 μg/g ww. Bivalves from the Mississippi River had OC levels of 0.096 μg/g ww, whereas none were detectable in bivalves at the marsh sites. In this study, p,p ’-dichlorodiphenyldichloroethylene ( p,p’ -DDE) (a biodegradation product of DDT [dichlorodiphenyl trichloroethane]) and total PCBs were the most frequently detected OCs and were primarily from the Mississippi River. For total OC content, using adjusted least squares means, some significant interactions were noted between fish species and sites. PAHs were detected in aquatic animals at all sites (range of 0.017–17.534 μg/g ww), as were AHs (range of 0.423–4.549 μg/g ww); the highest levels of PAHs and AHs were found in bivalves from the Mississippi River. When analysis of variance (α = 0.05) was performed with data from aquatic animals, there were only two significant relationships between PAHs, AHs, and OCs between species, site, and age or the interaction among these variables. There was an interaction between fish species and n-decane (an AH) in that mullet and largemouth bass had significantly higher levels than did catfish ( P = 0.0175). When general linear means were used to investigate associations of inorganic contaminants among fish species, site, and age or any interactions among these variables, no significant results were noted for arsenic, cadmium, lead, beryllium, boron, molybdenum, or nickel. The range of mercury in fish in this study was 0.04–0.14 μg/g ww (0.14– 0.48 μg/g dry weight [dw]), with the most elevated levels detected in predatory largemouth bass at the sampling point farthest downstream from the structure and within the marsh area. Mercury was positively correlated with fish age ( P = 0.0152), where levels were estimated to increase 0.0253 parts per million (ppm) dw per year. In the Mississippi River, catfish showed significantly higher levels of mercury than did mullet or largemouth bass ( P = 0.00167). Among fish species, mullet displayed the highest levels in fish of aluminum, barium, manganese, and iron, all considered to have low toxicity in hydrologic systems. An interaction between fish and site was seen with aluminum ( P = 0.0031), where concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with barium ( P = 0.0009), chromium ( P = <0.0001), manganese ( P = 0.0004), strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059). For iron ( P = 0.0.0001), mullet and largemouth bass at both the Mississippi River and Lake Salvador showed higher levels than did catfish, and these two species showed higher levels at two of the four sites. An interaction between fish and site was also seen with chromium ( P = <0.0001) in that concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059), metals for which deleterious effects have been demonstrated in other ecosystems. The NCBP program lists the 85th percentile for zinc at 34.2 μg/g fw (117.9 μg/g dw). In the Davis Pond prediversion biomonitoring study (hereafter the current study), one fish (MUL31RIVER, fish ID 8) showed values higher than that (125.4 μg/g dw or 37.54 μg/g ww), and the Mississippi River bivalve sample (MUSSRIVER) had a value of 140 μg/g dw (41.2 μg/g ww). In the current study, approximately 86 percent of the fish had measurable selenium levels, yet none reached the 85th percentile. The 85th percentile for selenium from the NCBP was 0.73 μg/g ww. Significantly higher levels of selenium were seen in mullet than in largemouth bass and catfish ( P = 0.0023). The NCBP 85th percentile for lead is 0.22 μg/g ww (0.76 μg/g dw). In the current study, the range of concentrations of lead was as much as 18.3 ppm dw (MUL31RIVER, fish ID 8), with the three most elevated values (range of 3.46–5.31 μg/g ww) coming from mullet from the Mississippi River. Biomarker data are measurable and directly reflect the condition of the animal, and measuring more than one biomarker in an individual increases confidence in health assessments. In the current study, biomarkers included macrophage aggregates (MAs), liver (hepatosomatic index [HSI]) and spleen (splenosomatic index [SSI]) weight to body weight ratios, percent white blood cells (WBCs) in whole blood, and condition factor. Few significant differences were noted with any of the biomarkers between sites, and there were no relationships between species and sites. For improved use of biomarker assessments, an increase in fish sample size would be useful for postdiversion sampling, as would comparisons of fish of the same sex and reproductive condition. During the current study, success for eagle nests in the diversion area and reference sites was similar as determined by numbers of nestlings fledged. When temperatures were below average during winter 2000, nests in both regions similarly failed. At each nest, the primary evidence of food items was small mammals. Eaglets (n = 6) generally appeared healthy, and whole blood concentrations of organic contaminants exceeded detection limits with three incidences of p,p’ -DDE (0.002–0.006 μg/L ww) and one incidence of oxychlordane (0.002 μg/L ww). The levels of p,p’ -DDE were well below those that have been inversely correlated with productivity and success rates of nesting bald eagles on a regional scale. The low values found in the whole blood samples for OC pesticides and PCBs were even lower when corrected for plasma volume. Aluminum values were 3.66 and 5.75 μg/L in two samples, zinc ranged from 5.21 to 6.77 μg/L ww in six samples, and silicon ranged from 1.7 to 4.6 μg/L in four samples. Selenium was detectable in each bird with the range at 0.332–0.566 μg/L ww, and strontium ranged from 0.0581 to 0.0975 μg/L ww. Mercury was detectable in blood samples from each bird and ranged from 0.0254 to 0.0845 μg/L ww, whereas lead was detectable in four samples and ranged from 0.0042 to 0.0136 μg/L ww. Although no detectable levels of total PCBs were found (also correlated with decreased reproductive productivity), 70 percent of the aquatic animals from the Mississippi River contained total PCBs (range 0.13–0.79 μg/L), whereas only about 7 percent of the aquatic animals sampled from the marsh area contained PCBs. Suggestions for postdiversion sampling include lowering the analytical detection limit for some metals, sampling aquatic animals over the course of a single season, obtaining a higher sample number of mature fish of one species (for example, blue catfish) within a range of total lengths for biomarker analyses, obtaining otoliths for estimating fish ages, assessing dioxins in eaglet blood, examining triazines in water, and obtaining all Mississippi River fish samples as close to the Davis Pond structure intake as possible. Because contaminants found in blood of eaglets reflect their prey species and because of the contaminant levels found in fish in the current study, eaglets may not be consuming primarily these species; therefore, obtaining juvenile nutria ( Myocastor coypus ) or turtle species for contaminant analyses might be considered, as well as collecting greater blood volume and using plasma to measure OCs and PCBs. Data obtained postdiversion will be compared with prediversion data to monitor changes.

Louisiana↗

Fishway Entrance Palisade

This technical report summarizes the work that was conducted by the University of Massachusetts Amherst and the United States Geological Survey (USGS), along with other project partners, on the Fishway Entrance Palisade (EP), a projected funded through the Department of Energy’s (DOE) funding opportunity titled ‘Innovative Solutions for Fish Passage at Hydropower Dams’ (DE‐FOA‐0001662). The period of performance ranged from September 1, 2018 through September 30, 2021. The EP is a novel fish passage engineering technology designed to provide more favorable entry conditions for fish and to reduce costs relative to conventional fishway auxiliary water systems (AWS). The EP project has four primary components. First, the Northeast United States Auxiliary Water Systems Database was created (Northeast Fishway Auxiliary Water Systems Database Section). The database, developed with material provided by the U.S. Fish and Wildlife Service, contains information on fishway type (e.g., lift, Denil, pool and weir) and Auxiliary Water System (AWS) details (e.g., water conveyance method, diffuser type) for 60 hydroelectric sites in the region. Findings indicate that nearly 4 out of every 10 fishway in the region is a fish lift and approximately 1 out of every 4 is a Denil ladder. The remainder are a mix of vertical slot fishways, pool and weirs, and Ice Harbor fishways. Furthermore, over half of all AWS systems use floor diffusers to discharge the auxiliary (or attraction) water into the entrance of a fishway, whereas only 14% use wall diffusers. Second, limited experiments on a conventional AWS with live, actively migrating fish were conducted at the USGS Easter Ecological Science Center (EESC) S.O. Conte Research Laboratory (Conventional Auxiliary Water System Experiments Section). This study determined how water velocity through a wall diffuser, without turning vanes or timber baffles to distribute the flow, affects the behavior and passage of adult American shad, a conservative surrogate species for migratory fish on the East Coast. Two gross diffuser velocity treatments were examined, 0.5 ft/s and 1.0 ft/s. These wall diffuser velocities represented current (0.5 ft/s) and past (1.0 ft/s) design criteria guidelines set forth by the USFWS North Atlantic-Appalachian Region (Rojas 2020; USFWS 2019). Six trials with a total of 151 American Shad were conducted in June of 2019 for the two treatments. No differences in American shad passage efficiency were discovered between the two treatments, while approximately 3 in every 4 attempts were successful at passing the diffuser. While these results may appear to indicate that the generally accepted gross wall diffuser velocity criteria for American shad of 0.5 ft/s could be safely increased to 1.0 ft/s, further analysis is warranted. Furthermore, it is unknown how other migratory and resident fish species that traverse these structures would be impacted by such a change. Studying the wall diffuser hydraulics led to an important AWS observation. Without turning vanes or timber baffles in this study, doubling the diffuser area was insufficient at producing the type of flow field change one may expect by halving the gross diffuser velocity. Instead, the flow fields throughout each treatments study area were similar, which led to similar results in shad performance. This not only highlights the importance of installing flow guidance devices like turning vanes, but also to the importance of properly maintaining them, which can be costly. Third, more expansive experiments on the novel EP were conducted in the spring of 2019 and 2021 (Fishway Entrance Palisade Experiments). The goal of this study was to determine how adult American shad responded to a variety of conditions at a full-scale EP. A total of six treatments were examined by changing the average auxiliary channel velocity between 1.0 and 5.0 ft/s in intervals of 1.0 ft/s and by inserting/removing an entrance gate at the opening of the fishway. Thirty trials with a total of 1,273 shad were conducted over the two years. In all treatments, at least ~7 out of every 10 fish successfully passed the EP diffuser and swam into the entrance channel within the 3.5-hour long trial, highlighting the general effectiveness of the novel AWS technology. In both study years, lower velocities through the EP diffuser led to increased shad performance, though performance peaked for the 2 ft/s velocity treatment. This treatment condition represents an approximate six-fold increase in gross diffuser velocity relative to conventional auxiliary water systems, which in turn presents opportunities for cost savings (e.g., reduction in diffuser size). Shad performance, in general, was worse in 2019 than in 2021, potentially due to the different run timing when our trials were conducted (2019 trials occurred near the end of the migration season, unlike in 2021). Treatments in 2019 had approximately a 20% reduction in entrance efficiency by the trial end, including a 16.7% drop for the 3 ft/s velocity treatment in 2019 relative to 2021 (the only carryover treatment between years). Lastly, adding an entrance gate caused a significant delay to entry. The time to 25% entry raised ~20 minutes from the near instantaneous 25% entry that was reported for the other treatments conducted in the same year (2021). Though by the end of the 3.5-hour trial, the overall entrance efficiency nearly matched those of the other 2021 treatments. The fourth and final component of the EP project was an economic analysis that focused on the cost of attraction and environmental flows (Modeling Power Generation Losses Due to Environmental and Fish Passage Attraction Flows at a Run-Of-River Hydroelectric Operation in the Northeast). The study assessed the economic impact of meeting environmental flow requirements at a representative hydroelectric facility and fish lift in the Northeast. An initial finding of the study was that there is a paucity of published data on the costs of meeting attraction and environmental flows. This is due, in part, to the proprietary nature of this data. To explore the costs associated with these flows, three types of environmental flows were assessed: upstream fishway attraction flows, downstream fishway attraction flows, and habitat maintenance minimum flows. A physics-based model was developed and calibrated with three years of hourly generation and flow data as inputs. Gage flow inputs were adjusted and used to calculate power generated. To address hydrologic variability, the model was executed to simulate 30 years of historical flows. Results indicate that both interannual and seasonal climatic factors impact the costs of meeting environmental flow requirements. Generation potential is most strongly curtailed during dry years in terms of maximizing the capacity factor (the percent of time a plant generates at capacity). Dry years, and especially dry summers, have the most significant costs associated with mitigation flows. Of the three types of flows, habitat flows are most costly in terms of power production, followed by upstream attraction flows. Downstream attraction flows are least costly. This finding is the likely result of differences in both flow rates and duration of the seasonal requirement for each flow. Overall, environmental flows represented a 2-12% loss in annual generation, but losses during a dry summer can reach over 20%.

Final Technical Report↗

Reorganization of vegetation, hydrology and soil carbon after permafrost degradation across heterogeneous boreal landscapes

The diversity of ecosystems across boreal landscapes, successional changes after disturbance and complicated permafrost histories, present enormous challenges for assessing how vegetation, water and soil carbon may respond to climate change in boreal regions. To address this complexity, we used a chronosequence approach to assess changes in vegetation composition, water storage and soil organic carbon (SOC) stocks along successional gradients within four landscapes: (1) rocky uplands on ice-poor hillside colluvium, (2) silty uplands on extremely ice-rich loess, (3) gravelly–sandy lowlands on ice-poor eolian sand and (4) peaty–silty lowlands on thick ice-rich peat deposits over reworked lowland loess. In rocky uplands, after fire permafrost thawed rapidly due to low ice contents, soils became well drained and SOC stocks decreased slightly. In silty uplands, after fire permafrost persisted, soils remained saturated and SOC decreased slightly. In gravelly–sandy lowlands where permafrost persisted in drier forest soils, loss of deeper permafrost around lakes has allowed recent widespread drainage of lakes that has exposed limnic material with high SOC to aerobic decomposition. In peaty–silty lowlands, 2–4 m of thaw settlement led to fragmented drainage patterns in isolated thermokarst bogs and flooding of soils, and surface soils accumulated new bog peat. We were not able to detect SOC changes in deeper soils, however, due to high variability. Complicated soil stratigraphy revealed that permafrost has repeatedly aggraded and degraded in all landscapes during the Holocene, although in silty uplands only the upper permafrost was affected. Overall, permafrost thaw has led to the reorganization of vegetation, water storage and flow paths, and patterns of SOC accumulation. However, changes have occurred over different timescales among landscapes: over decades in rocky uplands and gravelly–sandy lowlands in response to fire and lake drainage, over decades to centuries in peaty–silty lowlands with a legacy of complicated Holocene changes, and over centuries in silty uplands where ice-rich soil and ecological recovery protect permafrost.

Alaska↗

Swatara Creek basin of southeastern Pennsylvania: An evaluation of its hydrologic system

Local concentrations of population in the Swatara Creek basin of Pennsylvania find it necessary to store, transport, and treat water because local supplies are either deficient or have been contaminated by disposal of wastes in upstream areas. Water in the basin is available for the deficient areas and for dilution of the coal-mine drainage in the northern parts and the sewage wastes in the southern parts. Swatara Creek drains 576 square miles just east of Harrisburg, Pa., and is the largest tributary to the Susquehanna River from the north side below Harrisburg. It rises in the southern Pocono Mountains and flows southwestward across the Lebanon Plateau. On an average day Swatara Creek discharges more than 630 million gallons into the Susquehanna River at Middletown, Pa. In a year this amounts to about 23 inches of water over the entire basin and is the residual from an average annual precipitation of 45.5 inches. During an average year the flow in Swatara Creek from the upper third of the basin above Harper Tavern is always greater than 1,300 mgd (million gallons per day) for at least 15 days and is always greater than 25 mgd for at least 350 days. The daily streamflow from the basin averages 1.1 mgd per sq mi, but yields from different areas range from 0.97 to 1.22 mgd per sq mi. These variations are caused chiefly by differences in precipitation and land cover. The area of lowest yield is in the valleys west of Tremont, and the highest yields are in the Upper and Lower Little Swatara Creek subbasins. At high and medium stages the chemical character of the water in the streams is suitable for public and private supplies. At lower stages, defending on the areas and the amounts of contamination by coal-mine drainage and sewage pollution, the natural flow may require some treatment. At low stages the chemical characteristics of the natural flow not affected by man is almost identical with that of the ground water in the area drained by the stream. In general, the total dissolved solids range from about 25 to 400 parts per million and the hardness is as much as about 300 parts per million. The ground-water increment to the base flow of Swatara Creek averages about 240 mgd, or about 8.8 inches annually, for the basin. Generally, ground-water supplies in amounts of less than 0.5 mgd can be developed south of Blue Mountain. Supplies of several million gallons per day have been developed for industrial use from the permeable limestones in the south-central part of the basin. More intensive investigation in other parts of the basin would indicate areas where supplies of more than 0.5 mgd could be developed from properly spaced wells. The chemical character of water from wells depends largely on the host rock. In highly soluble rocks water contains large amount of dissolved solids; in more resistant rocks concentrations are lower. The chemical character of unpolluted ground water generally reflects the composition of the more readily soluble minerals in the local geologic environment. Areas contaminated by septic- tank effluent may have above normal amounts of nitrate and detergent products. Except where polluted, most ground water is suitable for public and industrial uses without extensive treatment. Sites for storage of surface water exist in the part of the basin lying in the valley and ridge area. As much as 30 to 40 percent of the annual flow could be impounded for release as low-flow augmentation for dilution of mine drainage and other wastes in the basin. Low sediment yields of supplying drainage areas would ensure a long life expectancy of reservoirs at these sites. Overbank flooding of the main stem of the Swatara Creek and its tributaries has occurred many times in the past. However, it has not been a hazard because urban development has not encroached on the flood plain. An inundation map of the August 1933 flood provides a basis that urban planners may use to avoid future damage. As water in the Swatara Creek moves downstream to the Susquehanna River, the flow is influenced consecutively by a large annual rainfall on the northern valley and ridge area, the wastes of surface and subsurface coal-mining activities, and less annual rainfall on the part of the basin lying in the Lebanon Plateau area; the flow is supplemented and further influenced by many tributaries and by the industrial and domestic wastes that are carried by these secondary streams. The annual precipitation ranges from 52 inches at the east edge and 49 inches at the west edge of the mountainous part of the basin to about 41 inches at the southwestern part at Middletown. The rainfall generally is adequate during the growing season to mature the crops. The mean annual temperature at Lebanon is about 52&deg;F, and the growing season is about 180 days. In this report the basin has been divided into eight hydrologic zones, leased on runoff, natural use of water, and chemical character of water. Four zones lie in the valley and ridge area, three lie in the Lebanon Plateau area, and one lies in the highland along the southeastern basin boundary. In each of the zones the hydrologic characteristics are virtually the same, but they may be completely different from those in adjacent zones. The boundaries of the zones generally coincide with boundaries between geologic formations, and the areas in each zone include rocks of similar influence on water. Streams in zone 4 at the northeast edge of the plateau have the highest average surface runoff from 1.2 to 1.1 mgd per sq mi whereas those in zone 2 at the northwest edge of the valley and ridge area have the lowest, about 1.0 mgd. Streams in zone 8, along the southeast edge of the basin, have the largest sustained low-flow yield, about 0.26 to 0.19 mgd per sq mi; those in zone 5 overlying the Martinsburg Shale east of Harrisburg have the smallest sustained low-flow yields, 0.03 to 0.01 mgd. Streams in the limestone area of zone 7 have the greatest range in low-flow yields in any one zone from 0.60 to 0 mgd per sq mi. Low-flow yields in zones 1 through 4 range from 0.13 to 0.03 mgd per sq mi. Surface flows from zones 1 and 2 are generally acidic and contain high concentrations of sulfate, iron, and total dissolved solids especially where contaminated with mine wastes. Surface flows from zones 3 and 4 are dilute, slightly alkaline, and suitable for public water supplies. Surface flows from zones 5, 6, and 7 are alkaline and contain moderate concentrations of dissolved solids with waters of highest hardness occurring in zone 7. Surface flows from zone 8 are dilute to moderately mineralized and are relatively high in silica concentration. Nitrate concentrations are high in surf Fee flows below sewage outfalls and in ground water contaminated by septic tank effluent and industrial wastes. Average annual sediment yields of 550 to 650 tons per square mile are characteristic of zones 1 and 2 where strip mining has destroyed the forest cover and coal culm is carried into the streams. From agricultural lands on the Martinsburg Shale in zones 5 and 6, annual sediment yields range from 300 to 350 tons per square mile; but from agricultural lands on the siliceous rocks in zone 8 and zones 3 and 4 in the valley and ridge area, the sediment yield ranges from 200 to 250 tons annually per square mile. Lowest annual sediment yields in the basin are in the forested areas of siliceous rocks in zones 2, 3, 4, and 5, and in the sinkhole topography of the limestones in zone 7 where the yield ranges from 30 to 35 tons and 50 to 60 tons per square mile, respectively. The amount of ground water that can be developed in the basin is dependent on the ability of the underlying rocks to yield water to wells. More than 300 gpm (gallons per minute) can be obtained from wells in alluvial materials in the valley bottoms and in some of the limestones where large solution channels and fractures are penetrated by the wells. From 50 to 300 gpm can be obtained from wells in loosely cemented sandstones and in fractured limestones. From 10 to 50 gpm can be developed from wells in the shales and harder sandstones. The most dense rocks will yield from 1 to 10 gpm from fractures and crevices. Most wells yield water from the upper 350 feet of the formation, for this part contains the most fractures or solution channels. Studies show that the velocity at which a contaminant will move downstream in the basin is related to the discharge of the stream at the time. At a stream discharge of about 400 mgd at Pine Grove, a contaminant in Swatara Creek would require about 40 hours to move from Pine Grove to Middletown. As a result of dispersion and dilution, the maximum concentration of the contaminant at Middletown would be less than 20 percent the concentration at Pine Grove under these conditions. An evaluation of the availability of water in the basin indicates that about I,239 mgd enters as precipitation, 630 mgd leaves as streamflow, 580 mgd is evaporated and transpired, and 56 mgd is diverted for use by man. Not all the diversions for man's use are lost to the basin, as about 27 mgd is returned as sewage for reuse. About one-fourth of the waste water is returned to the ground and the remainder to stream drainageways. Of that diverted by man, 11.6 mgd is used for public supply and 44.4 mgd for industrial and private supplies. Diversions of streamflow furnish 86 percent of the public supply and 27 percent of the industrial supply, and ground-water sources yield the remainder. Municipal and private sewage treatment plants are upgrading the waste water in many places, but no provisions are being made for treatment other than natural dilution and assimilation for the 15 mgd of coal-nine drainage in the northern part of the basin. Technology for economic treatment of mine water is not available at this time, although research in this field is being done. Urbanization eastward from Harrisburg and around Lebanon has increased the population density of the basin. Densities of 500 people per square mile and water use exceeding 2.0 mgd per sq mi can be expected in the future. By the year 2000 the population of the basin may increase 60 percent; and if the per capita rate of use increases 0.5 percent per year the domestic requirements for water will be about two times the present use, or 23 mgd. Similarly, if the present 1:4 ratio of domestic use to industrial use of water continues, at least 89 mgd will be needed for industry in the future. Although an increase to twice the present use of water can be foreseen, or 112 mgd, water for the dilution and assimilation of wastes from treatment systems are not included. Providing water for dilution of wastes from treatment plants has not been a problem, but in the future the amounts needed for this purpose will be greater as the population increases. As water becomes more valuable, treatment of sewage wastes to reduce the biochemical-oxygen-demand load by at least 80 to 90 percent will be necessary to conserve water for more productive uses. As much as 100 mgd may be needed for waste dilution in the basin by year 2000. The present trends in suburban and light industrial development will probably persist in the basin. Problems arising through changes in economic value of water, conflicts in use, and alternatives in development are typical of those confronting the manager of a water-resource system.

Pennsylvania↗

Selenium hazards in the Salton Sea environment—Summary of current knowledge to inform future wetland management

Quaternary marine and continental shales in the western United States are sources of selenium that can be loaded into the aquatic environment through mining, agricultural, and energy production processes. The mobilization of selenium from shales through agricultural irrigation has been recognized since the 1930s; however, discovery of deformities in birds and other wildlife using agricultural habitats during the 1980s spurred studies to determine the extent and effects of the contamination. Through these early studies, researchers determined that biota in the Salton Sea drainage basin was at risk from legacy selenium contamination in the Colorado River watershed. The Salton Sea and its surrounding managed and unmanaged wetlands provide vital inland habitat and trophic support for diverse assemblages of resident and migratory wildlife, and understanding regional selenium hazards for these trust species is a priority for many Federal and State agencies. The modern Salton Sea is a shallow, landlocked saline lake in Riverside and Imperial Counties (not shown) of California that is sustained by irrigation return and perennial river inflow. Changes in water transfer agreements under the 2003 Quantification Settlement Agreement (QSA) have resulted in reduced irrigation flow, declining lake levels, and the evolution of unmanaged wetlands in areas where drains and rivers no longer reach the Salton Sea. These wetlands provide additional habitat for some species of concern, but their potential to increase selenium hazards for trust species is largely unknown. From the 1980s to 2020, efforts to document selenium contamination and effects throughout the region have resulted in a considerable amount of selenium data from the Salton Sea and its surrounding drainage basin; however, no long-term (greater than 20 years), consistent sampling program has been established, and all data have been collected by different entities using a variety of protocols and analytical techniques. This lack of coordination has been previously documented in regional management plans and has led to difficulty in reliably assessing selenium hazards in the Salton Sea environment. This report provides a summary of the available disparate selenium information collected from water, sediment, and biota in the Salton Sea region since the 1980s and to identify data gaps that need to be filled to understand the potential effects of selenium on species of concern, including federally endangered desert pupfish ( Cyprinodon macularius ) and Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis; formerly Yuma Clapper Rail, Rallus longirostris yumanensis ). Available data from the Salton Sea drainage basin show that water from the Colorado River has the lowest selenium concentration of all surface water sources. All other surface water flowing into the Salton Sea has elevated selenium concentrations due to evaporation and evapotranspiration that occurs in agricultural fields and associated water delivery infrastructure or leaching of selenium from irrigated farmland soils. The Salton Sea has lower selenium concentrations because of various biogeochemical processes that recycle selenium into the sediment or volatilize it to the atmosphere; however, these mechanisms are not well defined, and it is not clear if selenium cycling will change in response to possible changes in the oxidation state of the Salton Sea bottom waters as water levels decline. Agricultural drains have the highest average selenium concentrations, but few drains have been sampled since changes in irrigation practices have occurred (due to the 2003 QSA). Groundwater selenium concentrations are variable; some wells south of the Salton Sea have selenium concentrations as high as 300 micrograms per liter (µg/L), whereas selenium concentrations are below detection in other wells. Groundwater and surface-water geothermal discharge zones around the margins of the Salton Sea and in unmanaged wetlands have not been studied in detail, and published selenium measurements are not available for these surface features. Selenium concentrations in the sediment of the Salton Sea drainage basin are highest in wetland particulate organic matter and the Salton Sea lakebed, indicating that removal of selenium from the water to the sediment has been a primary mechanism for keeping selenium concentrations low in the water column. Sediment selenium concentrations in wetlands are lower than in the Salton Sea but higher than inflowing drains and rivers, indicating the lentic wetland sites also may be important sinks for selenium because of biogeochemical processes. Sediment selenium data have not been collected in agricultural drains since changes in irrigation practices occurred (due to the 2003 QSA), and it is unknown if selenium sequestration from the water column has changed in these systems. We divided biological data into broad taxonomic categories, including primary producers, invertebrates, herpetofauna, mammals, fishes, and birds to facilitate evaluation of selenium concentrations and spatiotemporal trends observed in the Salton Sea. Overall, selenium concentrations were substantially greater in algae samples compared to all vascular plant samples combined. Median selenium concentrations in several invertebrate taxa (Chironomidae, Formicidae, Corixidae, Corbiculidae and Nereididae, and Decapoda) exceeded the maximum suggested dietary threshold of 3.0–4.0 micrograms per gram (µg/g) dry weight (dw) for predators consuming invertebrates in aquatic food webs. The greatest number of samples were collected from fish, and selenium distributions among species and locations showed that the range for most samples was lower than the U.S. Environmental Protection Agency selenium criterion for aquatic life (8.5 µg/g dw whole body, 11.3 µg/g dw fillets). The median selenium concentrations for whole body fish were below the selenium criterion in most locations, except for bairdiella ( Bairdiella icistia ) from the Salton Sea and irrigation drains, a few individual tilapia spp. (family Cichlidae, including genera Tilapia, Oreochromis , and their hybrids) from the river and river outlets, and several western mosquitofish ( Gambusia affinis ) and sailfin molly ( Poecilia latipinna ) from irrigation drain outlets. For avian samples combined among years and locations, median selenium concentrations in livers from all families except waders and Ibis (family Threskiornithidae) were higher than levels expected to cause selenium toxicosis (10–20 µg/g dw), and all median egg concentrations were above or near 6.0 μg/g dw, which is a conservative threshold value for reproductive impairment. Most knowledge gaps we identified for water, sediment, and biota were interrelated, and the use of integrated approaches to address knowledge gaps can provide greater insight into the drivers behind selenium hazards. Integrated water, sediment, and biota studies could help identify cost-effective management solutions that serve multiple purposes. A comprehensive analysis of the hydrology, biogeochemistry, and food-web processes in wetlands and other habitats can inform predictive models to identify drivers of selenium bioavailability, uptake from the environment and subsequent trophic transfer, ultimately forming the basis for experimental habitat management manipulations to minimize selenium hazards to wildlife. Furthermore, a comprehensive, long-term sampling and analytical laboratory plan would enable comparison of data among different entities that are sampling at the Salton Sea. Such efforts are well suited to help fill knowledge gaps that preclude understanding of selenium hazards and future management options for biota using Salton Sea habitats, including newly formed wetlands throughout the region. All data compiled for this report are available in two U.S. Geological Survey data releases: Groover and others (2022) for water and sediment samples and De La Cruz and others (2022) for biological samples. The data releases include all publicly available data for selenium concentrations in water, sediment, and biological samples collected in and around the Salton Sea, including the Coachella and Imperial Valleys. The data releases also include previously unpublished data.

California↗