Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Urban Ecosystems”

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.

550 records · Page 31Linked to original sources
USGS2024

Evaluation of extinction risk for stream fishes within an urban riverscape using population viability analysis

1. The Santa Ana River in the Los Angeles region of California demonstrates common habitat degradation symptoms that are characteristic of the urban stream syndrome. These impacts have altered the Santa Ana River community structure, with few species as impacted as the native Santa Ana sucker (sucker; Pantosteus santaanae ). 2. Consequently, a recovery plan developed for sucker identified the need for a population viability analysis (PVA) to assess sucker extirpation risk. However, PVAs can be data-intensive and are subject to several sources of bias when standardized protocols are absent. 3. More than 20 years of sucker and arroyo chub (chub; Gila orcuttii ) surveys using different methods were compiled to build an integrated hierarchical multi-population PVA to estimate trends in abundance and extirpation probability of these native fishes from the Santa Ana River. 4. PVA modelling indicated similar patterns in sucker and chub abundance along the Santa Ana River, with the highest abundance of both species in the upper regions of the river during the early 2000s and downstream in recent years (2018–2022). Extirpation risk was estimated to be greatest near wastewater treatment facilities, where native fish abundance estimates have been zero since 2018. Extirpation risk was lower downstream of the wastewater treatment facilities for both species, although extinction risk was higher for sucker than chub throughout the river. 5. As the model evolves and more data are collected, the PVA could be used to assess the effects of various management actions, such as non-native predator removals and native fish re-introductions, on sucker and chub persistence.

California↗
USGS2010

Two mechanisms of aquatic and terrestrial habitat change along an Alaskan Arctic coastline

Arctic habitats at the interface between land and sea are particularly vulnerable to climate change. The northern Teshekpuk Lake Special Area (N-TLSA), a coastal plain ecosystem along the Beaufort Sea in northern Alaska, provides habitat for migratory waterbirds, caribou, and potentially, denning polar bears. The 60-km coastline of N-TLSA is experiencing increasing rates of coastline erosion and storm surge flooding far inland resulting in lake drainage and conversion of freshwater lakes to estuaries. These physical mechanisms are affecting upland tundra as well. To better understand how these processes are affecting habitat, we analyzed long-term observational records coupled with recent short-term monitoring. Nearly the entire coastline has accelerating rates of erosion ranging from 6 m/year from 1955 to 1979 and most recently peaking at 17 m/year from 2007 to 2009, yet an intensive monitoring site along a higher bluff (3–6 masl) suggested high interannual variability. The frequency and magnitude of storm events appears to be increasing along this coastline and these patterns correspond to a greater number of lake tapping and flooding events since 2000. For the entire N-TLSA, we estimate that 6% of the landscape consists of salt-burned tundra, while 41% is prone to storm surge flooding. This offset may indicate the relative frequency of low-magnitude flood events along the coastal fringe. Monitoring of coastline lakes confirms that moderate westerly storms create extensive flooding, while easterly storms have negligible effects on lakes and low-lying tundra. This study of two interacting physical mechanisms, coastal erosion and storm surge flooding, provides an important example of the complexities and data needs for predicting habitat change and biological responses along Arctic land–ocean interfaces.

Alaska↗
USGS2009

Sediment losses and gains across a gradient of livestock grazing and plant invasion in a cool, semi-arid grassland, Colorado Plateau, USA

Large sediment fluxes can have significant impacts on ecosystems. We measured incoming and outgoing sediment across a gradient of soil disturbance (livestock grazing, plowing) and annual plant invasion for 9 years. Our sites included two currently ungrazed sites: one never grazed by livestock and dominated by perennial grasses/well-developed biocrusts and one not grazed since 1974 and dominated by annual weeds with little biocrusts. We used two currently grazed sites: one dominated by annual weeds and the other dominated by perennial plants, both with little biocrusts. Precipitation was highly variable, with years of average, above-average, and extremely low precipitation. During years with average and above-average precipitation, the disturbed sites consistently produced 2.8 times more sediment than the currently undisturbed sites. The never grazed site always produced the least sediment of all the sites. During the drought years, we observed a 5600-fold increase in sediment production from the most disturbed site (dominated by annual grasses, plowed about 50 years previously and currently grazed by livestock) relative to the never grazed site dominated by perennial grasses and well-developed biocrusts, indicating a non-linear, synergistic response to increasing disturbance types and levels. Comparing sediment losses among the sites, biocrusts were most important in predicting site stability, followed by perennial plant cover. Incoming sediment was similar among the sites, and while inputs were up to 9-fold higher at the most heavily disturbed site during drought years compared to average years, the change during the drought conditions was small relative to the large change seen in the sediment outputs. ?? 2009 Elsevier B.V. All rights reserved.

Aeolian Research↗
USGS2007

The Yellowstone hotspot, Greater Yellowstone ecosystem, and human geography

Active geologic processes associated with the Yellowstone hotspot are fundamental in shaping the landscapes of the greater Yellowstone ecosystem (GYE), a high volcanic plateau flanked by a crescent of still higher mountainous terrain. The processes associated with the Yellowstone hotspot are volcanism, faulting, and uplift and are observed in the geology at the surface. We attribute the driving forces responsible for the northeastward progression of these processes to a thermal plume rising through the Earth’s mantle into the base of the southwest-moving North American plate. This progression began 16 million years ago (Ma) near the Nevada-Oregon border and arrived at Yellowstone about 2 Ma. Before arrival of the hotspot, an older landscape existed, particularly mountains created during the Laramide orogeny about 70–50 Ma and volcanic terrain formed by Absaroka andesitic volcanism mostly between 50–45 Ma. These landscapes were more muted than the present, hotspot-modified landscape because the Laramide-age mountains had worn down and an erosion surface of low relief had developed on the Absaroka volcanic terrain. The Yellowstone Plateau was built by hotspot volcanism of rhyolitic lavas and caldera-forming rhyolite tuffs (ignimbrites). Streams eroding back into the edges of this plateau have created scenic waterfalls and canyons such as the Grand Canyon of the Yellowstone and Lewis Canyon. Rhyolite is poor in plant nutrients and forms sandy, well-drained soils that support the monotonous, fire-adapted lodgepole pine forests of the Yellowstone Plateau. Non-rhyolitic rocks surround this plateau and sustain more varied vegetation, including spruce, fir, and whitebark pine forests broken by grassy meadows. Heat from the hotspot rises upward and drives Yellowstone’s famed geysers, hot springs, and mudpots. These thermal waters are home to specialized, primitive ecosystems, rich in algae and bacteria. The rock alteration associated with hydrothermal systems creates the bright colors of Yellowstone’s Grand Canyon. Basin-and-range-style faulting has accompanied migration of the hotspot to Yellowstone and formed the linear mountains and valleys that occur north and south of the hotspot track, which is the present-day eastern Snake River Plain. High rates of basin-and-range faulting occurred adjacent to the migrating Yellowstone hotspot, creating distinctive landscapes within the GYE such as the Teton Range/Jackson Hole, with characteristic rugged, forested ranges and adjacent flat-floored grassy valleys. The difference in altitude between the mountains and valleys provides a topographic gradient in which vegetation maturation advances with altitude; animal-migration patterns also follow this trend. The valleys provide natural meadows, agricultural land, town sites, and corridors for roads. Uplift of the GYE by as much as 1 km (3,000 ft) during the last 5 million years has resulted in ongoing erosion of deep, steep-walled valleys. Many prominent ecological characteristics of Yellowstone derive from this hotspot-induced uplift, including the moderate- to high- altitude terrain and associated cool temperatures and deep snowfall. Modern and Pleistocene climate and associated vegetation patterns strongly relate to the topography created by the hotspot and its track along the eastern Snake River Plain. Winter air masses from the moist northern Pacific Ocean traverse the topographic low of the Snake River Plain to where orographic rise onto the Yellowstone Plateau and adjacent mountains produces deep snow. A winter precipitation shadow forms on the lee (eastern) sides of the GYE. During Pleistocene glacial times, this moisture conduit provided by the hotspot-track-produced ice-age glaciers that covered the core of the present GYE. These glaciers sculpted bedrock and produced glacial moraines that are both forested and unforested, sand and gravel of ice-marginal streams and outwash gravels that are commonly covered with sagebrush-grassland, and silty lake sediments that are commonly covered by lush grassland such as Hayden Valley. The effects of the Yellowstone hotspot also profoundly shaped the human history in the GYE. Uplift associated with the hotspot elevates the GYE to form the Continental Divide, and streams drain radially outward like spokes from a hub. Inhabitants of the GYE 12,000–10,000 years ago, as well as more recent inhabitants, followed the seasonal green-up of plants and migrating animals up into the mountain areas. During European immigration, people settled around Yellowstone in the lower parts of the drainages and established roads, irrigation systems, and cultural associations. The core Yellowstone highland is too harsh for agriculture and inhospitable to people in the winter. Beyond this core, urban and rural communities exist in valleys and are separated by upland areas. The partitioning inhibits any physical connection of communities, which in turn complicates pursuit of common interests across the whole GYE. Settlements thus geographically isolated evolved as diverse, independent communities

Idaho, Montana, Utah, Wyoming↗
USGS2011

Causes and implications of the correlation between forest productivity and tree mortality rates

At global and regional scales, tree mortality rates are positively correlated with forest net primary productivity (NPP). Yet causes of the correlation are unknown, in spite of potentially profound implications for our understanding of environmental controls of forest structure and dynamics and, more generally, our understanding of broad-scale environmental controls of population dynamics and ecosystem processes. Here we seek to shed light on the causes of geographic patterns in tree mortality rates, and we consider some implications of the positive correlation between mortality rates and NPP. To reach these ends, we present seven hypotheses potentially explaining the correlation, develop an approach to help distinguish among the hypotheses, and apply the approach in a case study comparing a tropical and temperate forest. Based on our case study and literature synthesis, we conclude that no single mechanism controls geographic patterns of tree mortality rates. At least four different mechanisms may be at play, with the dominant mechanisms depending on whether the underlying productivity gradients are caused by climate or soil fertility. Two of the mechanisms are consequences of environmental selection for certain combinations of life-history traits, reflecting trade-offs between growth and defense (along edaphic productivity gradients) and between reproduction and persistence (as manifested in the adult tree stature continuum along climatic and edaphic gradients). The remaining two mechanisms are consequences of environmental influences on the nature and strength of ecological interactions: competition (along edaphic gradients) and pressure from plant enemies (along climatic gradients). For only one of these four mechanisms, competition, can high mortality rates be considered to be a relatively direct consequence of high NPP. The remaining mechanisms force us to adopt a different view of causality, in which tree growth rates and probability of mortality can vary with at least a degree of independence along productivity gradients. In many cases, rather than being a direct cause of high mortality rates, NPP may remain high in spite of high mortality rates. The independent influence of plant enemies and other factors helps explain why forest biomass can show little correlation, or even negative correlation, with forest NPP.

Ecological Monographs↗
USGS2021

Simulation of groundwater budgets and travel times for watersheds on the north shore of Long Island Sound, with implications for nitrogen-transport studies

Aquatic systems in and around the Long Island Sound (LIS) provide a variety of ecological and economic benefits, but in some areas of the LIS, aquatic ecosystems have become degraded by excess nitrogen. A substantial fraction of the nitrogen inputs to the LIS are transported through the groundwater-flow system. Because groundwater travel times in surficial aquifers can exceed 100 years, multiyear lags are introduced between inputs at the water table in recharge areas and discharge to inland or coastal receiving waters. The U.S. Geological Survey, in cooperation with the Connecticut Department of Energy and Environmental Protection and the U.S. Environmental Protection Agency’s Long Island Sound Study, developed a steady-state groundwater model of the watersheds draining from the northern shore of the LIS for the purpose of calculating groundwater budgets and travel times to coastal waters. The model was developed by using the MODFLOW–NWT software and existing spatial data on aquifers, river networks, land-surface altitudes, land cover, groundwater recharge, and water use. Coastal waters were delineated on the basis of the National Wetland Inventory; all non-coastal waters were collectively termed “inland waters.” A coarse-resolution model was calibrated by using the PEST++ software, long-term records of water levels in 65 wells, stream altitudes from 477 streams, base-flow records for 14 streamgages that are relatively unaffected by withdrawals, and error metrics based on incorrectly simulated flooding and incorrectly simulated dry streams. The calibrated values were used in a fine-resolution model in which the mean absolute residuals were 4.5 meters for groundwater levels, 1.3 meters for stream altitudes, and 7,200 cubic meters per day (2.9 cubic feet per second) for base flow. About 89 percent of the terrestrial cells were correctly simulated with the water table below land surface, and nearly 90 percent of the cells representing streams were correctly simulated as having the water table above the stream bottom. Together, these metrics suggest that this model is robust for simulating regional-scale groundwater patterns. Simulated groundwater budgets were compiled for the entire study area, for each HUC12 (Hydrologic Unit Code no. 12) watershed and its adjacent coastal waters, if applicable, within the study area, and for 14 coastal-embayment watersheds. Most groundwater (90.6 percent of inflows) discharged to inland waters, with smaller fractions to coastal waters (7.0 percent) and well withdrawals (2.4 percent). When computed for HUC12 watersheds with coastal discharge, the portions of groundwater discharging to coastal waters ranged from 0.02 to 66 percent of groundwater outflows, with a median of 13 percent. Within priority-embayment watersheds, the portions of groundwater discharging to coastal waters ranged from 2 to 56 percent, with a median of 15 percent. Groundwater travel times also were simulated for the entire study area, for each HUC12 watershed and its adjacent coastal waters, if applicable, within the study area and for 14 priority coastal embayments. Within the entire study area, the median groundwater travel time was 1.9 years, with an interquartile range of 0.1 to 5.9 years. Sensitivity analysis of groundwater travel times within a subbasin in the study area indicates that the travel times are a function of the grid resolution, with coarser grids resulting in shorter median travel times. Travel times for groundwater discharging to coastal waters were similar to travel times for groundwater discharging to inland waters, with a median of 1.9 years. Median travel times for the HUC12 watersheds ranged from 0.9 to 53.5 years, with a median of 1.8 years. Among HUC12 watersheds that include coastal areas, travel times for groundwater discharging to coastal waters ranged from less than 1 to 61.6 years, with a median of 2.8 years. The HUC12 watersheds with the longest simulated travel times were in the urban area near New York City where the model performance is less accurate. Median travel times for groundwater discharging to coastal waters within the priority-embayment watersheds ranged from less than 1 to 18.6 years, with a median of 2.3 years. A more focused analysis was conducted for the Niantic River watershed to demonstrate the applicability of the regional model to local-scale nitrogen-transport analyses by using nitrogen-input and -attenuation rates from literature sources. Nitrogen inputs were estimated by using land-cover-based loading factors, and attenuation was estimated by using attenuation factors based on geologic zones and soil properties. Based on this analysis, groundwater transports an estimated 22,000 kilograms of nitrogen per year (2.9 kilograms of nitrogen per hectare per year) to streams, rivers, and coastal waters within the Niantic River watershed. Approximately 36 percent of discharging nitrogen is from atmospheric-deposition sources, 38 percent is from fertilizers, and 26 percent is from septic systems. Most of the groundwater-transported nitrogen (88 percent) discharges first to streams and rivers, with only 12 percent discharging directly to coastal waters. Travel times for groundwater-transported nitrogen ranged from less than 1 day to more than 100 years, with a median of 1.6 years.

Connecticut, New York, Rhode Island↗
USGS2015

California State Waters Map Series — Offshore of San Francisco, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of San Francisco map area is centered on the City of San Francisco and the Golden Gate channel, a waterway that connects the Pacific Ocean to the San Francisco Bay between the Marin Headlands and San Francisco Peninsula. The San Francisco Bay Area is the second-largest urban area on the U.S. West Coast with a combined population of over seven million. The bay supports several major cargo ports and the Port of San Francisco’s Fisherman’s Wharf is a major center for Northern California’s commercial and sport fishing fleets. The coastal part of the map area predominantly consists of high bluffs and vertical sea cliffs shaped by uplift and erosion of the Marin Headlands and San Francisco Peninsula east of the San Andreas and San Gregorio Fault Zones. The seafloor in the map area extends from the shoreline and western end of the Golden Gate channel to water depths of about 30 to 50 m, except for the San Andreas graben area, where water depths reach 75 m. Sea-level rise, tidal currents, and tectonics have shaped bathymetry in the map area. During the Last Glacial Maximum, Sea level was about 125 m lower than present day and the shoreline was more than 45 km west of San Francisco near the Farallon Islands. At that time, the map area was part of a large alluvial plain connected to a drainage basin that included much of California’s Central Valley. A river system flowed westward through the narrows of the Golden Gate channel and an alluvial valley bounded to the north and south by bedrock highlands, including the present-day Pacifica-Pescadero and Bolinas shelves. Rising seas entered the Golden Gate about 11,000 to 10,000 years ago and subsequent marine flooding led to progressive growth of the San Francisco Bay. Strong tidal currents, accelerating through the relatively narrow Golden Gate, have scoured the bedrock channel to a depth of 113 m. East and west of the channel, tidal currents decelerate and form large fields of sand waves. Offshore of the Marin Headlands, eastward transfer of right-lateral fault slip in a complex of faults northwest of the map area has caused extension and the formation of a sediment basin called the San Andreas graben on the continental shelf. The accommodation space created by extension on the shelf and the proximity to sediment transported to the ocean through San Francisco Bay results in a sand-dominated offshore shelf environment. Seafloor habitats in the Offshore of San Francisco map area comprise significant sand-dominated sediment habitat with sand wave and ripple bedforms indicative of high wave and current energy. North of the Golden Gate, biological productivity resulting from coastal upwelling supports populations of Sooty Shearwater, Western Gull, Common Murre, Cassin’s Auklet, and many other less populous bird species. In addition, an observable recovery of Humpback and Blue Whales has occurred in the area; both species are dependent on coastal upwelling to provide nutrients. For the first time in 65 years, Pacific Harbor Porpoise returned to San Francisco Bay in 2009. On the coast north of the Golden Gate, the large extent of exposed inner shelf bedrock supports large forests of “bull kelp,” which is well adapted for high wave-energy environments. Common fish species found in the kelp beds and rocky reefs include painted greenling, kelp greenling, lingcod, and several varieties of rockfish. Circulation over the continental shelf in the Offshore of San Francisco map area is dominated by the southward-flowing California Current, an eastern limb of the North Pacific Gyre that flows from Oregon to Baja California. At its midpoint offshore of central California, the California Current transports subarctic surface waters southeastward, about 150 to 1,300 km from shore. Seasonal northwesterly winds that are, in part, responsible for the California Current, generate coastal upwelling. Ocean temperatures offshore of central California have increased over the past 50 years, driving an ecosystem shift from the productive subarctic regime towards a depopulated subtropical environment.

California↗
USGS2018

Preliminary synthesis and assessment of environmental flows in the middle Verde River watershed, Arizona

A 3-year study was undertaken to evaluate the suitability of the available modeling tools for characterizing environmental flows in the middle Verde River watershed of central Arizona, describe riparian vegetation throughout the watershed, and estimate sediment mobilization in the river. Existing data on fish and macroinvertebrates were analyzed in relation to basin characteristics, flow regimes, and microhabitat, and a pilot study was conducted that sampled fish and macroinvertebrates and the microhabitats in which they were found. The sampling for the pilot study took place at five different locations in the middle Verde River watershed. This report presents the results of this 3-year study. The Northern Arizona Groundwater Flow Model (NARGFM) was found to be capable of predicting long-term changes caused by alteration of regional recharge (such as may result from climate variability) and groundwater pumping in gaining, losing, and dry reaches of the major streams in the middle Verde River watershed. Over the period 1910 to 2006, the model simulated an increase in dry reaches, a small increase in reaches losing discharge to the groundwater aquifer, and a concurrent decrease in reaches gaining discharge from groundwater. Although evaluations of the suitability of using the NARGFM and Basin Characteristic Model to characterize various streamflow intervals showed that smallerscale basin monthly runoff could be estimated adequately at locations of interest, monthly stream-flow estimates were found unsatisfactory for determining environmental flows. Orthoimagery and Moderate Resolution Imaging Spectroradiometer data were used to quantify stream and riparian vegetation properties related to biotic habitat. The relative abundance of riparian vegetation varied along the main channel of the Verde River. As would be expected, more upland plant species and fewer lowland species were found in the upper-middle section compared to the lower-middle section, and vice-versa. Vegetation changes within the upper-middle and lower-middle reaches are related to differences in climate and hydrology. In general, the riparian vegetation of the middle Verde River watershed is that of a healthy ecosystem’s mixed age, mixed patch structure, likely a result of the mostly unaltered disturbance regime. The frequency of in-river hydrogeomorphic features (pool, riffle, run) varied along the middle Verde River channel. There was a greater abundance of riffle habitat in the upper-middle reach; the lower-middle reach included more pool habitat. The Oak Creek tributary was more homogenous in geomorphic stream habitat composition than West Clear Creek, where runs dominated the upper reaches and pools dominated many of the lower reaches. On the basis of the period of record and discharges recorded at 15-minute intervals, five flows were found to reach the gravel-transport threshold. Sediment mobilization computed with flows averaged over daily time steps yielded just three flows that reached the gravel-transport threshold, and monthly averaged flows yielded none. In the middle Verde River watershed, 15-minute data should be used when possible to evaluate sediment transport in the river system. Data from more than 300 fish surveys conducted from 1992 to 2011 were analyzed using two schemes, one that divided the river into five reaches based on basin characteristics, and a second that divided the river into five reaches based on degree of flow alteration (specifically, diversions). Fish community metrics and assemblage data were used to analyze patterns of species composition and abundance in the two approaches. Overall, native and non-native species were regularly interacting and probably competing for similar resources. Fish abundances were also analyzed in response to floods and other flow metrics. Although the data are limited, native fish abundances increased more rapidly than non-native fish abundances in response to large floods. The basin-characteristic reach analysis showed native fish in greater abundance in the upper-middle reaches of the Verde River watershed and generally decreasing with downstream distance. The median relative abundance of native fish decreased by 50 percent from reach 1 to reach 5. Using the reach scheme based on degree of flow alteration, nondiverted reaches were found to have a greater abundance of native fish than diverted reaches. In heavily diverted reaches, non-native species outnumbered native species. Fish metrics and stream-flow metrics for the 30, 90, and 365-day periods before collection were computed and the results analyzed statistically. Only abundance of all fish species was associated with the 30-day flow metrics. The 90-day flow metrics were generally positively associated with fish metrics, whereas the 365-day flow metrics had more negative correlations. In particular, significant relations were found between fish metrics and the magnitude and frequency of high flows, including maximum monthly flow, median annual number of high-flow events, and median annual maximum streamflow. Native sucker (Catostomidae) populations tended to decrease in periods of extended base flow, and fish in the non-native sunfish family (Centrarchidae) decreased in periods of flashy, high magnitude flows. A pilot study surveyed fish at five locations in the upper part of the middle Verde River watershed as a means to measure microhabitat availability and quantify native and non-native fish use of that available microhabitat. Results indicated that native and non-native species exhibit some clear differences in microhabitat use. Although at least some native and non-native fish were found in each velocity, depth, and substrate category, preferential microhabitat use was common. On a percentage basis, non-native species had a strong preference for slow-moving and deeper water with silt and sand substrate, with a secondary preference for faster moving and very shallow water and a coarse gravel substrate. Native species showed a general preference for somewhat faster, moderate depth water over coarse gravel and had no clear secondary preference. Macroinvertebrate-variables index period, high-flow year, and collection location (upper-middle Verde River, lowermiddle Verde River, or Verde River tributaries) were found to be important explanatory variables in differentiating among community metrics. Overall richness (number of unique taxa), Shannon’s diversity index, and the percent of the most dominant taxa were all highly correlated, but their response to each macroinvertebrate variable was different. The percentage of mayfly (order Ephemeroptera) taxa was significantly higher in Oak Creek and the upper-middle and lower-middle Verde River reaches, locations which have higher flows and more urbanization than other reaches. When community metrics were related to hydrologic metrics, caddisfly (order Trichoptera) populations appeared to increase and mayfly populations to decrease in response to less flashy and more stable streamflows. Conversely, caddisfly populations appeared to decrease and mayfly populations to increase in response to greater flow variability. Six locations along the Verde River were sampled for macroinvertebrates as part of a pilot study associated with this report—(1) below Granite Creek, (2) near Campbell Ranch, (3) at the U.S. Geological Survey Paulden gage, (4) at the Perkinsville Bridge, (5) at the USGS Clarkdale gage, and (6) near the Reitz Ranch property. A nonmetric multidimensional scaling ordination of macroinvertebrate assemblages showed that the Verde River below Granite Creek site was different from the five other sites and that the Perkinsville Bridge and near Reitz Ranch samples had similar community structure. The near Campbell Ranch and Paulden gage locations had similar microhabitat characteristics, with the exception of riparian cover, yet the assemblage structure was very different. The different community composition at Verde River below Granite Creek was likely due to it having the smallest substrate sizes, lowest velocities, shallowest depths, and most riparian cover of the six sites.

Arizona↗
USGS2013

Method to support Total Maximum Daily Load development using hydrologic alteration as a surrogate to address aquatic life impairment in New Jersey streams

More than 300 ambient monitoring sites in New Jersey have been identified by the New Jersey Department of Environmental Protection (NJDEP) in its integrated water-quality monitoring and assessment report (that is, the 305(b) Report on general water quality and 303(d) List of waters that do not support their designated uses) as being impaired with respect to aquatic life; however, no unambiguous stressors (for example, nutrients or bacteria) have been identified. Because of the indeterminate nature of the broad range of possible impairments, surrogate measures that more holistically encapsulate the full suite of potential environmental stressors need to be developed. Streamflow alteration resulting from anthropogenic changes in the landscape is one such surrogate. For example, increases in impervious surface cover (ISC) commonly cause increases in surface runoff, which can result in “flashy” hydrology and other changes in the stream corridor that are associated with streamflow alteration. The NJDEP has indicated that methodologies to support a hydrologically based Total Maximum Daily Load (hydro-TMDL) need to be developed in order to identify hydrologic targets that represent a minimal percent deviation from a baseline condition (“minimally altered”) as a surrogate measure to meet criteria in support of designated uses. The primary objective of this study was to develop an applicable hydro-TMDL approach to address aquatic-life impairments associated with hydrologic alteration for New Jersey streams. The U.S. Geological Survey, in cooperation with the NJDEP, identified 51 non- to moderately impaired gaged streamflow sites in the Raritan River Basin for evaluation. Quantile regression (QR) analysis was used to compare flow and precipitation records and identify baseline hydrographs at 37 of these sites. At sites without an appropriately long period of record (POR) or where a baseline hydrograph could not be identified with QR, a rainfall-runoff model was used to develop simulated baseline hydrographs. The hydro-TMDL approach provided an opportunity to evaluate proportional differences in flow attributes between observed and baseline hydrographs and to develop complementary flow-ecology response relations at a subset of Raritan River Basin sites where available flow and ecological information overlapped. The New Jersey Stream Classification Tool (NJSCT) was used to determine the stream class of all 51 study sites by using either an observed or a simulated baseline hydrograph. Two New Jersey stream classes (A and C) were evaluated to help characterize the unique hydrology of the Raritan River Basin. In general, class C streams (1.99–40.7 square miles) had smaller drainage areas than class A streams (0.7–785 square miles). Many of the non-impaired and moderately impaired class A and C streams in the Raritan River Basin were found to have significant hydrologic alteration as indicated by numerous flow values that fell outside the established 25th-to-75th- and the more conservative 40th-to-60th-percentile boundaries. However, percent deviations for the class C streams (defined as moderately stable streams with moderately high base-flow contributions) were, in general, much larger than those for the class A streams (defined as semiflashy streams characterized by moderately low base flow). The greater deviations for class C streams in the hydro-TMDL assessments likely resulted from comparisons that were based solely on simulated baseline hydrographs, which were developed without considering any anthropogenic influences in the basin. In contrast, comparisons for many of the class A streams were made by using an observed baseline, which already includes an implicit level of ISC and other human influences on the landscape. By using the hydro-TMDL approach, numerous flow deviations were identified that were indicative of streams that are highly regulated by reservoirs or dams, streams that are affected by increasing amounts of surface runoff resulting from ISC, and streams that are affected by water abstraction (that is, groundwater or surface-water withdrawals used for agricultural and human supply). Eight of the reservoir- and (or) dam-affected sites showed flow deviations that are indicative of flow-managed systems. For example, indices that account for the timing and magnitude of high and low flows were often found to fall outside the 25th-to-75th-percentile range. In general, at regulated class C streams, annual summer low flows are arriving later and tend to be lower, and high flows are arriving earlier with higher magnitudes of longer duration. At class A streams, high and low flows are arriving later with an overall increase in discharge with respect to the prereservoir baseline conditions. The drainage basins of eight of the study sites had large values of ISC (>10 percent), most likely as a result of expanding urban development. In general, the magnitude and frequency of high flows at class A and C sites with high ISC are increasing and were commonly found to fall outside the 25th-to-75th-percentile range. Additionally, magnitudes of low flows are becoming lower and, although the timing of high flows was highly variable, low-flow events appeared to be arriving earlier than would be expected under normal low-flow conditions. Three of the study sites appeared to be affected by hydrologic changes associated with water abstraction. At these sites, the timing of flows appeared to be altered. For example, low flows tended to arrive earlier and high flows arrived later at two of the three sites. Additionally, the magnitude and duration of low flows were commonly less than the 25th-percentile value and the duration of high flows appeared to increase. A reduced set of hydrologic and ecological variables was used to develop univariate and multivariate flow-ecology response models for the aquatic-invertebrate assemblage. Many hydrologic variables accounting for the duration, magnitude, frequency, and timing of flows were significantly correlated with ecological response. Multiple linear regression (MLR) models were developed to provide a more holistic evaluation of the combined effects of hydrologic alteration and to identify models with two or three hydrologic variables that account for a significant proportion of the variability in invertebrate-assemblage condition as represented by assemblage metric scores. MLR models, derived on the basis of hydrologic attributes, accounted for 35 to 75 percent of the variability in assemblage condition. The hydro-TMDL method developed herein for non- to moderately impaired Raritan River Basin streams utilizes a “surrogate” approach in place of the traditional “pollutant of concern” approach commonly used for TMDL development. Managers can use the results obtained by using the hydro-TMDL method to offset the effects of impervious-surface runoff and altered streamflow and to implement measures designed to achieve the necessary load reductions for the “pollutant of concern” (that is, percentage deviations of stream-class-specific flow-index values outside the established 25th-to-75th-percentile range). In this case, such deviations could represent all or a subset of the altered flow indices that prevent the stream from meeting designated aquatic-life criteria. This hydro-TMDL uses a reference, or attainment stream approach for developing the TMDL endpoint. That is, either observed or simulated baseline hydrographs were selected as appropriate reference conditions on the basis of results of QR analysis and watershed modeling procedures, respectively. For any stream in the Raritan River Basin evaluated as part of this study, the hydro-TMDL can be expressed as the greatest amount of deviation in flow a stream can exhibit without violating the stream’s designated aquatic-life criteria. Use of this surrogate approach is appropriate because flows that fall outside the established percentile ranges are ultimately a function of many anthropogenic modifications of the landscape, including the amount of stormwater runoff generated from impervious surfaces within a given basin, the presence of manmade structures designed to retain or divert water, the magnitude of ground- and surface-water abstraction, and the presence of water-supply processes implemented to support human needs. In addition, the stream-type-specific flow indices used as the basis for the hydro-TMDL approach are useful for representing the hydrologic conditions of class A and C streams/basins because they incorporate the full spectrum of flow conditions (very low to very high) that occur in the stream system over a long period of time, as well as those flow properties that change as a result of seasonal variation. Ultimately, an estimate of the maximum percentage flow reduction that could be allowed will be needed to address the aquatic-life impairments in many of the study streams in the Raritan River Basin and will be necessary for identifying appropriate target flow conditions for hydro-TMDL implementation. As described in this report, a target flow value equal to the 25th- or 75th-percentile flow rate could be selected as the point useful for setting specific hydrologic targets. This selection, however, is a management decision that could vary depending on the designated use of the stream or other regulatory factors (for example, water-supply protection, trout production, antidegradation policies, or special protection designations). In New Jersey streams where no unambiguous stressors can be identified, State monitoring agencies, such as the NJDEP, could choose to require the implementation of a flow-based TMDL that not only supports designated uses, but meets the regulatory requirements under the Clean Water Act, and represents a balance between water supply intended to meet human needs and the conservation of ecosystem integrity.

New Jersey↗
USGS2020

Fate and behavior tools related to inland spill response—Workshop on the U.S. Geological Survey’s role in Federal science support

Executive Summary There is a growing body of tools available for science support for determining the fate and behavior of industrial and agricultural chemicals that are rapidly injected (“spilled”) into aquatic environments. A 2-day roundtable-style workshop was held by the U.S. Geological Survey (USGS) in Middleton, Wisconsin, in December 2017 to describe and explore existing Federal science support for spill fate and behavior tools used for inland spills, ongoing and new fate and behavior studies, and science gaps in planning and response tools as part of the USGS Midcontinent Region’s efforts to include spill response as part of its strategic plans. A total of 28 attendees representing a variety of Federal, State, and regional entities presented on programs and tools used in various aspects of spill response. Most programs and tools discussed were for spills in riverine environments but tools and applications for spills in lakes, on land surfaces, in urban storm sewer networks, and groundwater also were discussed. A primary workshop focus was to facilitate communication and increase potential for future collaboration among agencies for inland spill science support. The role and need for more USGS science support within the inland spill community was discussed. Enhanced communication is needed within the USGS and the U.S. Department of the Interior science programs, as well as within and among other agencies that do emergency planning and response. A main conclusion of the workshop was that there are untapped resources of the USGS outlined in the agency’s science strategy that could strengthen science support for fate and behavior tools in inland areas, especially in the Upper Mississippi River, Ohio River, and Great Lakes Basins where large freshwater resources overlap with dense corridors of oil and hazardous substances, with transportation networks, and with large populations centers. Fate and behavior tools are being developed quickly for inland spill response by multiple Federal agencies in partnership with local and regional entities. Applicability of these tools ranges from planning and preparedness, to the early stages of spill response for protection of human life and property, and to the application of monitoring and models to assess the long-term consequences of spills. Key findings from the workshop, with an emphasis on potential further development of USGS science support, include the following: •The national and regional response to spills occurs within an established system that must be respected by all parties involved in spill response. The USGS’s role is to support spill responders who are physically working at a spill scene, deploying booms and using other efforts to contain and recover spilled materials. •The USGS has tools that have been used throughout spill response operations, from early response to recovery and restoration. Developing a more formal role for the USGS to participate in science support for inland spills on a consistent basis is a desired outcome. This will require the USGS to improve internal and external communication and would be best accomplished by assigning one or more coordinator positions within the agency to plan and oversee USGS spill-response efforts. More involvement of the USGS on National and Regional Response Teams, especially in the realm of the Science and Technology Subcommittees, will gofar in increasing external communication and integration of fate and behavior tools. •Rapid response to spills requires modeling and mapping of plumes and associated time-of-travel estimation for a range of stream sizes across the United States. Many existing models use USGS streamgage data and the USGS National Hydrography Dataset. Nearly all existing models would benefit from updated linkages to USGS StreamStats and its soon-to-be released time-of-travel estimates,real-time velocity, stream morphology, and slope data. Integrating USGS tools with those from other agencies could be done to better serve the larger spill response community. • A problem is that existing models to rapidly predict plume extent, as well as more followup/longer-term fate and transport models, can be unknown or unavailable to spill responders. Thus, creating and strengthening linkages among USGS scientists skilled at using these tools is needed to support spill response with the on-scene responders. • Research for inland spill fate and behavior done outside of an immediate spill response can assist with spill planning and preparedness by (1) revealing sites likely to experience spills in the future (high-risk sites) and (2) understanding how a spilled substance might behave under a range of environmental conditions. However, USGS research on this topic has been scarce and subject to funding availability. Examples include the 2010 Line 6B Spill release into the Kalamazoo River in Michigan, where the USGS provided science support for a variety of fate and behavior tools for stream and impoundment environments. A long-term research site in Bemidji, Minnesota, provides important insights into transformations and longevity of spilled oil in groundwater and groundwater-surface water interactions. • Linking stream models to other components of this inland environment, including groundwater, overland flow, and karst, is needed. Stream network data can be linked to underground conduits such as storm sewers and karst groundwater systems. Stream models can also be linked with geospatial data such as that contained in U.S. Environmental Protection Agency’s interactive mapping tools. • The USGS is uniquely qualified to collect water-quality data during spills in the United States because of its many geographically dispersed water science centers, its knowledge and preparedness for flood measurement and documentation, and its cadre of skilled water-quality employees. Rapid-deployment gages, used for floods, could also be used for spills if they included spill-specific sensors. Coordinated expertise at USGS water and environmental science centers can be used for monitoring spill effects and for assessing risk to water quality and ecological communities. • Scientists at the USGS have proven capable of providing science coordination and technical assistance within the Incident Command Structure at the request of the lead on-scene coordinator. This external coordination, as well as internal communication within USGS Water, Hazards, and Ecosystems Mission Areas, could be improved by establishing and naming a USGS spills coordinator. Scott Morlock, Jo Ellen Hinck, and Faith Fitzpatrick are currently (2017) serving in informal coordination roles in addition to their traditional duties.

Open-File Report↗