Seasonal and longitudinal water quality dynamics in three effluent-dependent rivers in Arizona.

Effluent-fed streams, which receive inputs from wastewater treatment plants, are becoming increasingly common across the globe as urbanization intensifies. In semi-arid and arid regions, where many natural streams have dried up due to over extraction of water, many streams rely completely on treated effluent to sustain baseflow during dry seasons. These systems are often thought of as 'second-class' or highly disturbed stream ecosystems, but they have the potential to serve as refuges for native aquatic biota if water quality is high, especially in areas where few natural habitats remain. In this study, we investigated seasonal and longitudinal water quality dynamics at multiple sites across six reaches of three effluent-dependent rivers in Arizona (USA) with the objective (1) to quantify changes in effluent water quality due to distance traveled and season/climate and (2) to qualify whether water quality conditions in these systems are sufficient to support native aquatic species. Study reaches ranged in length from 3 to 31 km and in geographic setting from low desert to montane conifer forest. We observed the lowest water quality conditions (e.g., elevated temperature and low dissolved oxygen) during the summer in low desert reaches, and significantly greater natural remediation of water quality in longer vs. shorter reaches for several factors, including temperature, dissolved oxygen and ammonia. Nearly all sites met or exceeded water quality conditions needed to support robust assemblages of native species across multiple seasons. However, our results also indicated that temperature (max 34.2 °C), oxygen levels (min 2.7 mg/L) and ammonia concentrations (max 5.36 mg/L N) may occasionally be stressful for sensitive taxa at sites closest to effluent outfalls. Water quality conditions may be a concern during the summer. Overall, effluent-dependent streams have the capacity to serve as refuges for native biota in Arizona, and they may become the only aquatic habitat available in many urbanizing arid and semi-arid regions.

A global perspective on the functional responses of stream communities to flow intermittence.

The current erosion of biodiversity is a major concern that threatens the ecological integrity of ecosystems and the ecosystem services they provide. Due to global change, an increasing proportion of river networks are drying and changes from perennial to non-perennial flow regimes represent dramatic ecological shifts with potentially irreversible alterations of community and ecosystem dynamics. However, there is minimal understanding of how biological communities respond functionally to drying. Here, we highlight the taxonomic and functional responses of aquatic macroinvertebrate communities to flow intermittence across river networks from three continents, to test predictions from underlying trait-based conceptual theory. We found a significant breakpoint in the relationship between taxonomic and functional richness, indicating higher functional redundancy at sites with flow intermittence higher than 28%. Multiple strands of evidence, including patterns of alpha and beta diversity and functional group membership, indicated that functional redundancy did not compensate for biodiversity loss associated with increasing intermittence, contrary to received wisdom. A specific set of functional trait modalities, including small body size, short life span and high fecundity, were selected with increasing flow intermittence. These results demonstrate the functional responses of river communities to drying and suggest that on-going biodiversity reduction due to global change in drying river networks is threatening their functional integrity. These results indicate that such patterns might be common in these ecosystems, even where drying is considered a predictable disturbance. This highlights the need for the conservation of natural drying regimes of intermittent rivers to secure their ecological integrity.

Impacts of baseflow and flooding on microplastic pollution in an effluent-dependent arid land river in the USA.

Effluent discharge from wastewater treatment plants can be a substantial source of microplastics in receiving water bodies including rivers. Despite growing concern about microplastic pollution in freshwater habitats, the literature has not yet addressed effluent-dependent rivers, which derive 100% of their baseflow from effluent. The objective of this study was to document and explore trends in microplastic pollution within the effluent-dependent lower Santa Cruz River near Tucson, Arizona (USA). We examined microplastic concentrations in the water column and benthic sediment and microplastic consumption by mosquitofish (Gambusia affinis) at 10 sites along a ~40 km stretch of the lower Santa Cruz River across two time periods: baseflow (effluent only) and post-flood (effluent immediately following urban runoff). In total, across both sampling periods, we detected microplastics in 95% of water column samples, 99% of sediment samples, and 6% of mosquitofish stomachs. Flow status (baseflow vs post-flood) was the only significant predictor of microplastic presence and concentrations in our models. Microplastic fragment concentrations in the water column were higher post-flood, microplastic fiber concentrations in benthic sediment were lower post-flood, and mosquitofish were more likely to have consumed microplastics post-flood than during baseflow. The additional microplastics detected after flooding was likely due to a combination of allochthonous material entering the channel via runoff and bed scour that exhumed microplastics previously buried in the riverbed. Effluent-dependent urban streams are becoming increasingly common; more work is needed to identify microplastic pollution baselines and trends in effluent rivers worldwide.

River ecosystem conceptual models and non-perennial rivers: A critical review.

Conceptual models underpin river ecosystem research. However, current models focus on continuously flowing rivers and few explicitly address characteristics such as flow cessation and drying. The applicability of existing conceptual models to nonperennial rivers that cease to flow (intermittent rivers and ephemeral streams, IRES) has not been evaluated. We reviewed 18 models, finding that they collectively describe main drivers of biogeochemical and ecological patterns and processes longitudinally (upstream-downstream), laterally (channel-riparian-floodplain), vertically (surface water-groundwater), and temporally across local and landscape scales. However, perennial rivers are longitudinally continuous while IRES are longitudinally discontinuous. Whereas perennial rivers have bidirectional lateral connections between aquatic and terrestrial ecosystems, in IRES, this connection is unidirectional for much of the time, from terrestrial-to-aquatic only. Vertical connectivity between surface and subsurface water occurs bidirectionally and is temporally consistent in perennial rivers. However, in IRES, this exchange is temporally variable, and can become unidirectional during drying or rewetting phases. Finally, drying adds another dimension of flow variation to be considered across temporal and spatial scales in IRES, much as flooding is considered as a temporally and spatially dynamic process in perennial rivers. Here, we focus on ways in which existing models could be modified to accommodate drying as a fundamental process that can alter these patterns and processes across spatial and temporal dimensions in streams. This perspective is needed to support river science and management in our era of rapid global change, including increasing duration, frequency, and occurrence of drying.

What's in a Name? Patterns, Trends, and Suggestions for Defining Non-Perennial Rivers and Streams.

Rivers that cease to flow are globally prevalent. Although many epithets have been used for these rivers, a consensus on terminology has not yet been reached. Doing so would facilitate a marked increase in interdisciplinary interest as well as critical need for clear regulations. Here we reviewed literature from Web of Science database searches of 12 epithets to learn (Objective 1-O1) if epithet topics are consistent across Web of Science categories using latent Dirichlet allocation topic modeling. We also analyzed publication rates and topics over time to (O2) assess changes in epithet use. We compiled literature definitions to (O3) identify how epithets have been delineated and, lastly, suggest universal terms and definitions. We found a lack of consensus in epithet use between and among various fields. We also found that epithet usage has changed over time, as research focus has shifted from description to modeling. We conclude that multiple epithets are redundant. We offer specific definitions for three epithets (non-perennial, intermittent, and ephemeral) to guide consensus on epithet use. Limiting the number of epithets used in non-perennial river research can facilitate more effective communication among research fields and provide clear guidelines for writing regulatory documents.

If you build it, they will come: rapid colonization by dragonflies in a new effluent-dependent river reach.

BACKGROUND: Aquatic ecosystems are greatly altered by urban development, including the complete loss of natural habitat due to water diversions or channel burial. However, novel freshwater habitats also are created in cities, such as effluent-dependent streams that rely on treated wastewater for flow. It is unclear how diverse these novel ecosystems are, or how quickly aquatic species are able to colonize them. In this study, we (1) quantify odonate (Insecta, Odonata) colonization of a novel effluent-dependent river reach, (2) examine how drying events affect odonates in these novel habitats, and (3) explore whether effluent-dependent streams can support diverse odonate assemblages. METHODS: We conducted monthly odonate surveys at three sites along the Santa Cruz River (Tucson, AZ, USA) between June 2019 and May 2020. One site was in a long-established effluent-dependent reach (flowing since the 1970s) that served as a reference site and two sites were in a newly-established reach that began flowing on June 24, 2019 (it was previously dry). We compared odonate species richness, assemblage composition, and colonization patterns across these reaches, and examined how these factors responded to flow cessation events in the new reach. RESULTS: Seven odonate species were observed at the study sites in the new reach within hours of flow initiation, and species rapidly continued to arrive thereafter. Within 3 months, species richness and assemblage composition of adult odonates were indistinguishable in the new and reference reaches. However, drying events resulted in short-term and chronic reductions in species richness at one of the sites. Across all three sites, we found over 50 odonate species, which represent nearly 40% of species known from the state of Arizona. DISCUSSION: Odonates were surprisingly diverse in the effluent-dependent Santa Cruz River and rapidly colonized a newly established reach. Richness levels remained high at study sites that did not experience drying events. These results suggest that consistent discharge of high-quality effluent into dry streambeds can be an important tool for promoting urban biodiversity. However, it remains to be seen how quickly and effectively less vagile taxa (e.g., mayflies, caddisflies) can colonize novel reaches. Effluent-dependent urban streams will always be highly managed systems, but collaboration between ecologists and urban planners could help to maximize aquatic biodiversity while still achieving goals of public safety and urban development.

Zero or not? Causes and consequences of zero-flow stream gage readings.

Streamflow observations can be used to understand, predict, and contextualize hydrologic, ecological, and biogeochemical processes and conditions in streams. Stream gages are point measurements along rivers where streamflow is measured, and are often used to infer upstream watershed-scale processes. When stream gages read zero, this may indicate that the stream has fully dried; however, zero-flow readings can also be caused by a wide range of other factors. Our ability to identify whether or not a zero-flow gage reading indicates a dry fluvial system has far reaching environmental implications. Incorrect identification and interpretation by the data user can lead to hydrologic, ecological, and/or biogeochemical predictions from models and analyses. Here, we describe several causes of zero-flow gage readings: frozen surface water, flow reversals, instrument error, and natural or human-driven upstream source losses or bypass flow. For these examples, we discuss the implications of zero-flow interpretations. We also highlight additional methodss for determining flow presence, including direct observations, statistical methods, and hydrologic models, which can be applied to interpret causes of zero-flow gage readings and implications for reach- and watershed-scale dynamics. Such efforts are necessary to improve our ability to understand and predict surface flow activation, cessation, and connectivity across river networks. Developing this integrated understanding of the wide range of possible meanings of zero-flows will only attain greater importance in a more variable and changing hydrologic climate.

Simulating rewetting events in intermittent rivers and ephemeral streams: A global analysis of leached nutrients and organic matter.

Climate change and human pressures are changing the global distribution and the extent of intermittent rivers and ephemeral streams (IRES), which comprise half of the global river network area. IRES are characterized by periods of flow cessation, during which channel substrates accumulate and undergo physico-chemical changes (preconditioning), and periods of flow resumption, when these substrates are rewetted and release pulses of dissolved nutrients and organic matter (OM). However, there are no estimates of the amounts and quality of leached substances, nor is there information on the underlying environmental constraints operating at the global scale. We experimentally simulated, under standard laboratory conditions, rewetting of leaves, riverbed sediments, and epilithic biofilms collected during the dry phase across 205 IRES from five major climate zones. We determined the amounts and qualitative characteristics of the leached nutrients and OM, and estimated their areal fluxes from riverbeds. In addition, we evaluated the variance in leachate characteristics in relation to selected environmental variables and substrate characteristics. We found that sediments, due to their large quantities within riverbeds, contribute most to the overall flux of dissolved substances during rewetting events (56%-98%), and that flux rates distinctly differ among climate zones. Dissolved organic carbon, phenolics, and nitrate contributed most to the areal fluxes. The largest amounts of leached substances were found in the continental climate zone, coinciding with the lowest potential bioavailability of the leached OM. The opposite pattern was found in the arid zone. Environmental variables expected to be modified under climate change (i.e. potential evapotranspiration, aridity, dry period duration, land use) were correlated with the amount of leached substances, with the strongest relationship found for sediments. These results show that the role of IRES should be accounted for in global biogeochemical cycles, especially because prevalence of IRES will increase due to increasing severity of drying events.

Seasonality and predictability shape temporal species diversity.

Temporal environmental fluctuations, such as seasonality, exert strong controls on biodiversity. While the effects of seasonality are well known, the predictability of fluctuations across years may influence seasonality in ways that are less well understood. The ability of a habitat to support unique, non-nested assemblages of species at different times of the year should depend on both seasonality (occurrence of events at specific periods of the year) and predictability (the reliability of event recurrence) of characteristic ecological conditions. Drawing on tools from wavelet analysis and information theory, we developed a framework for quantifying both seasonality and predictability of habitats, and applied this using global long-term rainfall data. Our analysis predicted that temporal beta diversity should be maximized in highly predictable and highly seasonal climates, and that low degrees of seasonality, predictability, or both would lower diversity in characteristic ways. Using stream invertebrate communities as a case study, we demonstrated that temporal species diversity, as exhibited by community turnover, was determined by a balance between temporal environmental variability (seasonality) and the reliability of this variability (predictability). Communities in highly seasonal mediterranean environments exhibited strong oscillations in community structure, with turnover from one unique community type to another across seasons, whereas communities in aseasonal New Zealand environments fluctuated randomly. Understanding the influence of seasonal and other temporal scales of environmental oscillations on diversity is not complete without a clear understanding of their predictability, and our framework provides tools for examining these trends at a variety of temporal scales, seasonal and beyond. Given the uncertainty of future climates, seasonality and predictability are critical considerations for both basic science and management of ecosystems (e.g., dam operations, bioassessment) spanning gradients of climatic variability.

Linking multidimensional functional diversity to quantitative methods: a graphical hypothesis--evaluation framework.

Functional trait analysis is an appealing approach to study differences among biological communities because traits determine species' responses to the environment and their impacts on ecosystem functioning. Despite a rapidly expanding quantitative literature, it remains challenging to conceptualize concurrent changes in multiple trait dimensions ("trait space") and select quantitative functional diversity methods to test hypotheses prior to analysis. To address this need, we present a widely applicable framework for visualizing ecological phenomena in trait space to guide the selection, application, and interpretation of quantitative functional diversity methods. We describe five hypotheses that represent general patterns of responses to disturbance in functional community ecology and then apply a formal decision process to determine appropriate quantitative methods to test ecological hypotheses. As a part of this process, we devise a new statistical approach to test for functional turnover among communities. Our combination of hypotheses and metrics can be applied broadly to address ecological questions across a range of systems and study designs. We illustrate the framework with a case study of disturbance in freshwater communities. This hypothesis-driven approach will increase the rigor and transparency of applied functional trait studies.

Dispersal ability and habitat requirements determine landscape-level genetic patterns in desert aquatic insects.

Species occupying the same geographic range can exhibit remarkably different population structures across the landscape, ranging from highly diversified to panmictic. Given limitations on collecting population-level data for large numbers of species, ecologists seek to identify proximate organismal traits-such as dispersal ability, habitat preference and life history-that are strong predictors of realized population structure. We examined how dispersal ability and habitat structure affect the regional balance of gene flow and genetic drift within three aquatic insects that represent the range of dispersal abilities and habitat requirements observed in desert stream insect communities. For each species, we tested for linear relationships between genetic distances and geographic distances using Euclidean and landscape-based metrics of resistance. We found that the moderate-disperser Mesocapnia arizonensis (Plecoptera: Capniidae) has a strong isolation-by-distance pattern, suggesting migration-drift equilibrium. By contrast, population structure in the flightless Abedus herberti (Hemiptera: Belostomatidae) is influenced by genetic drift, while gene flow is the dominant force in the strong-flying Boreonectes aequinoctialis (Coleoptera: Dytiscidae). The best-fitting landscape model for M. arizonensis was based on Euclidean distance. Analyses also identified a strong spatial scale-dependence, where landscape genetic methods only performed well for species that were intermediate in dispersal ability. Our results highlight the fact that when either gene flow or genetic drift dominates in shaping population structure, no detectable relationship between genetic and geographic distances is expected at certain spatial scales. This study provides insight into how gene flow and drift interact at the regional scale for these insects as well as the organisms that share similar habitats and dispersal abilities.

Demographic stability metrics for conservation prioritization of isolated populations.

Systems of geographically isolated habitat patches house species that occur naturally as small, disjunct populations. Many of these species are of conservation concern, particularly under the interacting influences of isolation and rapid global change. One potential conservation strategy is to prioritize the populations most likely to persist through change and act as sources for future recolonization of less stable localities. We propose an approach to classify long-term population stability (and, presumably, future persistence potential) with composite demographic metrics derived from standard population-genetic data. Stability metrics can be related to simple habitat measures for a straightforward method of classifying localities to inform conservation management. We tested these ideas in a system of isolated desert headwater streams with mitochondrial sequence data from 16 populations of a flightless aquatic insect. Populations exhibited a wide range of stability scores, which were significantly predicted by dry-season aquatic habitat size. This preliminary test suggests strong potential for our proposed method of classifying isolated populations according to persistence potential. The approach is complementary to existing methods for prioritizing local habitats according to diversity patterns and should be tested further in other systems and with additional loci to inform composite demographic stability scores.

Evolution of aquatic insect behaviours across a gradient of disturbance predictability.

Natural disturbance regimes--cycles of fire, flood, drought or other events--range from highly predictable (disturbances occur regularly in time or in concert with a proximate cue) to highly unpredictable. While theory predicts how populations should evolve under different degrees of disturbance predictability, there is little empirical evidence of how this occurs in nature. Here, we demonstrate local adaptation in populations of an aquatic insect occupying sites along a natural gradient of disturbance predictability, where predictability was defined as the ability of a proximate cue (rainfall) to signal a disturbance (flash flood). In controlled behavioural experiments, populations from predictable environments responded to rainfall events by quickly exiting the water and moving sufficiently far from the stream to escape flash floods. By contrast, populations from less predictable environments had longer response times and lower response rates, reflecting the uncertainty inherent to these environments. Analysis with signal detection theory showed that for 13 out of 15 populations, observed response times were an optimal compromise between the competing risks of abandoning versus remaining in the stream, mediated by the rainfall-flood correlation of the local environment. Our study provides the first demonstration that populations can evolve in response to differences in disturbance predictability, and provides evidence that populations can adapt to among-stream differences in flow regime.