Multiple stressors alter greenhouse gas concentrations in streams through local and distal processes

Streams are significant contributors of greenhouse gases (GHG) to the atmosphere, and the increasing number of stressors degrading freshwaters may exacerbate this process, posing a threat to climatic stability. However, it is unclear whether the influence of multiple stressors on GHG concentrations in streams results from increases of in-situ metabolism (i.e., local processes) or from changes in upstream and terrestrial GHG production (i.e., distal processes). Here, we hypothesize that the mechanisms controlling multiple stressor effects vary between carbon dioxide (CO2) and methane (CH4), with the latter being more influenced by changes in local stream metabolism, and the former mainly responding to distal processes. To test this hypothesis, we measured stream metabolism and the concentrations of CO2 (pCO2) and CH4 (pCH4) in 50 stream sites that encompass gradients of nutrient enrichment, oxygen depletion, thermal stress, riparian degradation and discharge. Our results indicate that these stressors had additive effects on stream metabolism and GHG concentrations, with stressor interactions explaining limited variance. Nutrient enrichment was associated with higher stream heterotrophy and pCO2, whereas pCH4 increased with oxygen depletion and water temperature. Discharge was positively linked to primary production, respiration and heterotrophy but correlated negatively with pCO2. Our models indicate that CO2-equivalent concentrations can more than double in streams that experience high nutrient enrichment and oxygen depletion, compared to those with oligotrophic and oxic conditions. Structural equation models revealed that the effects of nutrient enrichment and discharge on pCO2 were related to distal processes rather than local metabolism. In contrast, pCH4 responses to nutrient enrichment, discharge and temperature were related to both local metabolism and distal processes. Collectively, our study illustrates potential climatic feedbacks resulting from freshwater degradation and provides insight into the processes mediating stressor impacts on the production of GHG in streams.

Non-interactive effects drive multiple stressor impacts on the taxonomic and functional diversity of atlantic stream macroinvertebrates

Freshwaters are considered among the most endangered ecosystems globally due to multiple stressors, which coincide in time and space. These local stressors typically result from land-use intensification or hydroclimatic alterations, among others. Despite recent advances on multiple stressor effects, current knowledge is still limited to manipulative approaches minimizing biological and abiotic variability. Thus, the assessment of multiple stressor effects in real-world ecosystems is required. Using an extensive survey of 50 stream reaches across North Portugal, we evaluated taxonomic and functional macroinvertebrate responses to multiple stressors, including marked gradients of nutrient enrichment, flow reduction, riparian vegetation structure, thermal stress and dissolved oxygen depletion. We analyzed multiple stressor effects on two taxonomic (taxon richness, Shannon-diversity) and two trait-based diversity indices (functional richness, functional dispersion), as well as changes in trait composition. We found that multiple stressors had additive effects on all diversity metrics, with nutrient enrichment identified as the most important stressor in three out of four metrics, followed by dissolved oxygen depletion and thermal stress. Taxon richness, Shannon-diversity and functional richness responded similarly, whereas functional dispersion was driven by changes in flow velocity and thermal stress. Functional trait composition changed along a major stress gradient determined by nutrient enrichment and oxygen depletion, which was positively correlated with organisms possessing fast-living strategies, aerial respiration, adult phases, and gathering-collector feeding habits. Overall, our results reinforce the need to consider complementary facets of biodiversity to better identify assembly processes in response to multiple stressors. Our data suggest that stressor interactions may be less frequent in real-word streams than predicted by manipulative experiments, which can facilitate mitigation strategies. By combining an extensive field survey with an integrative consideration of multiple biodiversity facets, our study provides new insights that can help to better assess and manage rivers in a global change context.

Multiple drying aspects shape dissolved organic matter composition in intermittent streams

Water availability is a fundamental driver of biogeochemical processing in highly dynamic ecosystems such as intermittent rivers and ephemeral streams (IRES), which are recognized as the most common fluvial ecosystem globally. Because of their global extent, IRES have a remarkable contribution to organic matter processing, which is expected to intensify as climate change and water extraction expand IRES extension. Nevertheless, the effect of the complexity of the drying process on river biogeochemistry remains unclear. This study investigated how drying aspects affect the dissolved organic carbon (DOC) concentration and composition in 35 streams along a wide flow-intermittence gradient in the NE Iberian Peninsula. To do that, four drying aspects: annual drying duration, annual frequency, duration of the last drying event, and time since the last drying event were characterized. Results showed that DOC concentration and the contribution of humic-like compounds were positively associated with intensifying drying conditions. In addition, protein-like compounds decreased over the drying gradient. More specifically, changes in DOC concentration were driven mainly by annual drying duration, whereas annual drying frequency and the duration of the last drying event jointly explained dissolved organic matter composition. These results suggest that the quantity and composition of dissolved organic matter in streams respond differently to the temporal aspects of the drying process. Our study can help to better anticipate changes in organic matter in the context of climate change.

biomonitoR: an R package for managing ecological data and calculating biomonitoring indices

The monitoring of biological indicators is required to assess the impacts of environmental policies, compare ecosystems and guide management and conservation actions. However, the growing availability of ecological data has not been accompanied by concomitant processing tools able to facilitate data handling and analysis. Multiple common challenges limit the usefulness of biomonitoring information across ecosystems and biological groups. Biomonitoring data analysis is currently constrained by time-consuming steps for data preparation and a data processing environment with limited integration in terms of software, biological groups, and protocols. We introduce biomonitoR, a package for the R programming language that addresses technical challenges for the management of ecological data and metrics calculation. biomonitoR implements most of the biological indices currently used or proposed in different fields of ecology and water resource management. Its combination of customizable functions aims to support a transferable and comprehensive biomonitoring workflow in a user-friendly environment. biomonitoR represents a versatile toolbox with five main assets: (i) it checks taxonomic information against reference datasets allowing for customization of trait and sensitivity scores; (ii) it supports heterogeneous taxonomic resolution allowing computations at multiple taxonomic levels; (iii) it calculates multiple biological indices, including metrics for both broad and stressor-specific ecological assessments; (iv) it enables user-friendly data visualization, helping both decision-making processes and data interpretation; and (v) it allows working with an interactive web application straight from R. Overall, biomonitoR can benefit the wide biomonitoring community, including environmental private consultants, ecologists and natural resource managers.

Towards women-inclusive ecology: Representation, behavior, and perception of women at an international conference

Conferences are ideal platforms for studying gender gaps in science because they are important cultural events that reflect barriers to women in academia. Here, we explored women’s participation in ecology conferences by analyzing female representation, behavior, and personal experience at the 1st Meeting of the Iberian Society of Ecology (SIBECOL). The conference had 722 attendees, 576 contributions, and 27 scientific sessions. The gender of attendees and presenters was balanced (48/52% women/men), yet only 29% of the contributions had a woman as last author. Moreover, men presented most of the keynote talks (67%) and convened most of the sessions. Our results also showed that only 32% of the questions were asked by women, yet the number of questions raised by women increased when the speaker or the convener was a woman. Finally, the post-conference survey revealed that attendees had a good experience and did not perceive the event as a threatening context for women. Yet, differences in the responses between genders suggest that women tended to have a worse experience than their male counterparts. Although our results showed clear gender biases, most of the participants of the conference failed to detect it. Overall, we highlight the challenge of increasing women’s scientific leadership, visibility and interaction in scientific conferences and we suggest several recommendations for creating inclusive meetings, thereby promoting equal opportunities for all participants.

Functional traits driving species role in the structure of terrestrial vertebrate scavenger networks

Species assemblages often have a non-random nested organization, which in vertebrate scavenger (carrion-consuming) assemblages is thought to be driven by facilitation in competitive environments. However, not all scavenger species play the same role in maintaining assemblage structure, as some species are obligate scavengers (i.e., vultures) and others are facultative, scavenging opportunistically. We used a database with 177 vertebrate scavenger species from 53 assemblages in 22 countries across five continents to identify which functional traits of scavenger species are key to maintaining the scavenging network structure. We used network analyses to relate ten traits hypothesized to affect assemblage structure with the "role" of each species in the scavenging assemblage in which it appeared. We characterized the role of a species in terms of both the proportion of monitored carcasses on which that species scavenged, or scavenging breadth (i.e., the species "normalized degree"), and the role of that species in the nested structure of the assemblage (i.e., the species "paired nested degree"), therefore identifying possible facilitative interactions among species. We found that species with high olfactory acuity, social foragers, and obligate scavengers had the widest scavenging breadth. We also found that social foragers had a large paired nested degree in scavenger assemblages, probably because their presence is easier to detect by other species to signal carcass occurrence. Our study highlights differences in the functional roles of scavenger species and can be used to identify key species for targeted conservation to maintain the ecological function of scavenger assemblages.

Remote sensing depicts riparian vegetation responses to water stress in a humid Atlantic region

Riparian areas in the Cantabrian Atlantic ecoregion (northwest Portugal) play a key role in soil formation and conservation, regulation of nutrient and water cycle, creation of landscape aesthetic value and the preservation of biodiversity. The maintenance of their ecological integrity is crucial given the ever increase in multiple anthropogenic (water demand and agriculture) and climatic pressures (droughts and extreme events). We developed a transferable remote sensing approach, taking advantage of the latest freely available technologies (Sentinel-2 and Copernicus Land products), to detect intra-annual and inter-annual changes in riparian vegetation productivity at the river basin scale related to water stress. This study has used the normalized difference vegetation index (NDVI) to investigate riparian vegetation productivity dynamics on three different vegetation types (coniferous, broadleaved and grassland) over the past 5 years (2015-2019). Our results indicated that inter-annual seasonality differed between drier (2017) and wetter (2016) years. We found that intra-annual dynamics of NDVI were influenced by the longitudinal river zonation. Our model ranked first (r²m = 0.73) showed that the productivity of riparian vegetation during the dry season was positively influenced by annual rainfall and by the type of riparian vegetation. The emergent long lags between climatic variation and riparian plant productivity provides opportunities to forecast early warnings of climatically-driven impacts. In addition, the different average productivity levels among vegetation types should be considered when assessing climatic impacts on riparian vegetation. Future applications of Sentinel 2 products could seek to distinguish riparian areas that are likely to be more vulnerable to changes in the annual water balance from those that are more resistant under longer-term changes in climate.

Making waves. Bridging theory and practice towards multiple stressor management in freshwater ecosystems

Despite advances in conceptual understanding, single-stressor abatement approaches remain common in the management of fresh waters, even though they can produce unexpected ecological responses when multiple stressors interact. Here we identify limitations restricting the development of multiple-stressor management strategies and address these, bridging theory and practice, within a novel empirical framework. Those critical limitations include that (i) monitoring schemes fall short of accounting for theory on relationships between multiple-stressor interactions and ecological responses, (ii) current empirical modelling approaches neglect the prevalence and intensity of multiple-stressor interactions, and (iii) mechanisms of stressor interactions are often poorly understood. We offer practical recommendations for the use of empirical models and experiments to predict the effects of freshwater degradation in response to changes in multiple stressors, demonstrating this approach in a case study. Drawing on our framework, we offer practical recommendations to support the development of effective management strategies in three general multiple-stressor scenarios.

DISPERSE, a trait database to assess the dispersal potential of European aquatic macroinvertebrates

Dispersal is an essential process in population and community dynamics, but is difficult to measure in the field. In freshwater ecosystems, information on biological traits related to organisms' morphology, life history and behaviour provides useful dispersal proxies, but information remains scattered or unpublished for many taxa. We compiled information on multiple dispersal-related biological traits of European aquatic macroinvertebrates in a unique resource, the DISPERSE database. DISPERSE includes nine dispersal-related traits subdivided into 39 trait categories for 480 taxa, including Annelida, Mollusca, Platyhelminthes, and Arthropoda such as Crustacea and Insecta, generally at the genus level. Information within DISPERSE can be used to address fundamental research questions in metapopulation ecology, metacommunity ecology, macroecology and evolutionary ecology. Information on dispersal proxies can be applied to improve predictions of ecological responses to global change, and to inform improvements to biomonitoring, conservation and management strategies. The diverse sources used in DISPERSE complement existing trait databases by providing new information on dispersal traits, most of which would not otherwise be accessible to the scientific community.

The interruption of longitudinal hydrological connectivity causes delayed responses in dissolved organic matter

Hydrology is the main driver of dissolved organic matter (DOM) dynamics in intermittent rivers and ephemeral streams. However, it is still unclear how the timing and the spatial variation in flow connectivity affect the dynamics of DOM and inorganic solutes. This study focuses on the impact of flow cessation on the temporal and spatial heterogeneity of DOM quantity and quality along an intermittent stream. We monitored a headwater intermittent stream at high spatial and temporal frequencies during a summer drying episode and analysed dissolved organic carbon (DOC) and its spectroscopic properties, inorganic solutes and dissolved CO₂. The drying period determined the disruption of the fluvial continuum with a recession of stream continuum at a rate of ~60 m/d and the gradual formation of a patched system of isolated pools of different sizes. Our results showed that the period of time that had elapsed since isolated pool formation (CI-days) was an essential factor for understanding how drying shaped the biogeochemistry of the fluvial system. Overall, drying caused a high DOC concentration and an increase in the humic-like fluorescence signal. Additionally, solutes showed contrasting responses to hydrological disconnection. Electrical conductivity, for instance, is a clear "sentinel" of the fragmentation process because it starts to increase before the hydrological disruption occurs. In contrast, DOC, most spectroscopic DOM descriptors and CO₂ showed delayed responses of approximately 5-21 days after the formation of isolated pools. Furthermore, the spatial location and volume of each isolated pool seemed to exert a significant impact on most variables. In contrast, the temperature did not follow a clear pattern. These findings indicate that the fragmentation of longitudinal hydrological connectivity does not induce a single biogeochemical response but rather stimulates a set of solute-specific responses that generates a complex biogeochemical mosaic in a single fluvial unit.

Rethinking megafauna

Concern for megafauna is increasing among scientists and non-scientists. Many studies have emphasized that megafauna play prominent ecological roles and provide important ecosystem services to humanity. But, what precisely are 'megafauna'? Here, we critically assess the concept of megafauna and propose a goal-oriented framework for megafaunal research. First, we review definitions of megafauna and analyse associated terminology in the scientific literature. Second, we conduct a survey among ecologists and palaeontologists to assess the species traits used to identify and define megafauna. Our review indicates that definitions are highly dependent on the study ecosystem and research question, and primarily rely on ad hoc size-related criteria. Our survey suggests that body size is crucial, but not necessarily sufficient, for addressing the different applications of the term megafauna. Thus, after discussing the pros and cons of existing definitions, we propose an additional approach by defining two function-oriented megafaunal concepts: 'keystone megafauna' and 'functional megafauna', with its variant 'apex megafauna'. Assessing megafauna from a functional perspective could challenge the perception that there may not be a unifying definition of megafauna that can be applied to all eco-evolutionary narratives. In addition, using functional definitions of megafauna could be especially conducive to cross-disciplinary understanding and cooperation, improvement of conservation policy and practice, and strengthening of public perception. As megafaunal research advances, we encourage scientists to unambiguously define how they use the term 'megafauna' and to present the logic underpinning their definition.

Evaluating anthropogenic impacts on naturally stressed ecosystems: Revisiting river classifications and biomonitoring metrics along salinity gradients

Naturally stressed ecosystems hold a unique fraction of biodiversity. However, they have been largely ignored in biomonitoring and conservation programmes, such as the EU Water Framework Directive, while global change pressures are threatening their singular values. Here we present a framework to classify and evaluate the ecological quality of naturally stressed rivers along a water salinity gradient. We gathered datasets, including aquatic macroinvertebrate assemblages and environmental information, for 243 river locations across the western Mediterranean to: a) gauge the role of natural stressors (salinity) in driving aquatic community richness and composition; b) make river classifications by encompassing the wide range of environmental and biological variation exhibited by Mediterranean rivers; c) provide effective biomonitoring metrics of ecological quality for saline rivers. Our results showed that water salinity played a pivotal role in explaining the community richness and compositional changes in rivers, even when considering other key and commonly used descriptors, such as elevation, climate or lithology. Both environmental and biologically-based classifications included seven river types: three types of freshwater perennial rivers, one freshwater intermittent river type and three new saline river types. These new saline types were not included in previous classifications. Their validation by independent datasets showed that the saline and freshwater river types represented differentiable macroinvertebrate assemblages at family and species levels. Biomonitoring metrics based on the abundance of indicator taxa of each saline river type provided a much better assessment of the ecological quality of saline rivers than other widely used biological metrics and indices. Here we demonstrate that considering natural stressors, such as water salinity, is essential to design effective and accurate biomonitoring programmes for rivers and to preserve their unique biodiversity.

Insect communities in saline waters consist of realized but not fundamental niche specialists

Considering how organisms adapt to stress is essential if we are to anticipate biological responses to global change in ecosystems. Communities in stressful environments can potentially be assembled by specialists (i.e. species that only occur in a limited range of environmental conditions) and/or generalist species with wider environmental tolerances. We review the existing literature on the salinity tolerance of aquatic insects previously identified as saline specialists because they were exclusively found in saline habitats, and explore if these saline realized niche specialists are also specialists in their fundamental niches or on the contrary are fundamental niche generalist species confined to the highest salinities they can tolerate. The results suggest that species inhabiting saline waters are generalists in their fundamental niches, with a predominant pattern of high survival in freshwater-low salinity conditions, where their fitness tends to be similar or even higher than in saline waters. Additionally, their performance in freshwater tends to be similar to related strictly freshwater species, so no apparent trade-off of generalization is shown. These results are discussed in the framework of the ecological and evolutionary processes driving community assembly across the osmotic stress gradient, and their potential implications for predicting impacts from saline dilution and freshwater salinization.This article is part of the theme issue 'Salt in freshwaters: causes, ecological consequences and future prospects'.

Do all roads lead to Rome? Exploring community trajectories in response to anthropogenic salinization and dilution of rivers

Abiotic stress shapes how communities assemble and support ecological functions. However, it remains unclear whether artificially increasing or decreasing stress levels would lead to communities assembling predictably along a single axis of variation or along multiple context-dependent trajectories of change. In response to stress intensity alterations, we hypothesize that a single trajectory of change occurs when trait-based assembly prevails, while multiple trajectories of change arise when dispersal-related processes modify colonization and trait-filtering dynamics. Here, we tested these hypotheses using aquatic macroinvertebrates from rivers exposed to gradients of natural salinity and artificially diluted or salinized ion contents. Our results showed that trait-filtering was important in driving community assembly in natural and diluted rivers, while dispersal-related processes seemed to play a relevant role in response to salinization. Salinized rivers showed novel communities with different trait composition, while natural and diluted communities exhibited similar taxonomic and trait compositional patterns along the conductivity gradient. Our findings suggest that the artificial modification of chemical stressors can result in different biological communities, depending on the direction of the change (salinization or dilution), with trait-filtering, and organism dispersal and colonization dynamics having differential roles in community assembly. The approach presented here provides both empirical and conceptual insights that can help in anticipating the ecological effects of global change, especially for those stressors with both natural and anthropogenic origins.This article is part of the theme issue 'Salt in freshwaters: causes, ecological consequences and future prospects'.

Effects of salinity changes on aquatic organisms in a multiple stressor context

Under global change, the ion concentration of aquatic ecosystems is changing worldwide. Many freshwater ecosystems are being salinized by anthropogenic salt inputs, whereas many naturally saline ones are being diluted by agricultural drainages. This occurs concomitantly with changes in other stressors, which can result in additive, antagonistic or synergistic effects on organisms. We reviewed experimental studies that manipulated salinity and other abiotic stressors, on inland and transitional aquatic habitats, to (i) synthesize their main effects on organisms' performance, (ii) quantify the frequency of joint effect types across studies and (iii) determine the overall individual and joint effects and their variation among salinity-stressor pairs and organism groups using meta-analyses. Additive effects were slightly more frequent (54%) than non-additive ones (46%) across all the studies (n = 105 responses). However, antagonistic effects were dominant for the stressor pair salinity and toxicants (44%, n = 43), transitional habitats (48%, n = 31) and vertebrates (71%, n = 21). Meta-analyses showed detrimental additive joint effects of salinity and other stressors on organism performance and a greater individual impact of salinity than the other stressors. These results were consistent across stressor pairs and organism types. These findings suggest that strategies to mitigate multiple stressor impacts on aquatic ecosystems should prioritize restoring natural salinity concentrations.This article is part of the theme issue 'Salt in freshwaters: causes, ecological consequences and future prospects'.

Evaluating riparian solutions to multiple stressor problems in river ecosystems - A conceptual study

Rivers are among the most sensitive of all ecosystems to the effects of global change, but options to prevent, mitigate or restore ecosystem damage are still inadequately understood. Riparian buffers are widely advocated as a cost-effective option to manage impacts, but empirical evidence is yet to identify ideal riparian features (e.g. width, length and density) which enhance ecological integrity and protect ecosystem services in the face of catchment-scale stressors. Here, we use an extensive literature review to synthesise evidence on riparian buffer and catchment management effects on instream environmental conditions (e.g. nutrients, fine sediments, organic matter), river organisms and ecosystem functions. We offer a conceptual model of the mechanisms through which catchment or riparian management might impact streams either positively or negatively. The model distinguishes scale-independent benefits (shade, thermal damping, organic matter and large wood inputs) that arise from riparian buffer management at any scale from scale-dependent benefits (nutrient or fine sediment retention) that reflect stressor conditions at broader (sub-catchment to catchment) scales. The latter require concerted management efforts over equally large domains of scale (e.g. riparian buffers combined with nutrient restrictions). The evidence of the relationships between riparian configuration (width, length, zonation, density) and scale-independent benefits is consistent, suggesting a high certainty of the effects. In contrast, scale-dependent effects as well as the biological responses to riparian management are more uncertain, suggesting that ongoing diffuse pollution (nutrients, sediments), but also sources of variability (e.g. hydrology, climate) at broader scales may interfere with the effects of local riparian management. Without concerted management across relevant scales, full biological recovery of damaged lotic ecosystems is unlikely. There is, nevertheless, sufficient evidence that the benefits of riparian buffers outweigh potential adverse effects, in particular if located in the upstream part of the stream network. This supports the use of riparian restoration as a no-regrets management option to improve and sustain lotic ecosystem functioning and biodiversity.

Multiple stressor effects on biological quality elements in the Ebro River: Present diagnosis and predicted responses

Multiple abiotic stressors affect the ecological status of water bodies. The status of waterbodies in the Ebro catchment (NE Spain) is evaluated using the biological quality elements (BQEs) of diatoms, invertebrates and macrophytes. The multi-stressor influence on the three BQEs was evaluated using the monitoring dataset available from the catchment water authority. Nutrient concentrations, especially total phosphorus (TP), affected most of the analyzed BQEs, while changes in mean discharge, water temperature, or river morphology did not show significant influences. Linear statistical models were used to evaluate the change of water bodies' ecological status under different combinations of future socioeconomic and climate scenarios. Changes in land use, rainfall, water temperature, mean discharge, TP and nitrate concentrations were modeled according to the future scenarios. These revealed an evolution of the abiotic stressors that could lead to a general decrease in the ecosystem quality of water bodies within the Ebro catchment. This deterioration was especially evidenced on the diatoms and invertebrate biological indices, mainly because of the foreseen increase in TP concentrations. Water bodies located in the headwaters were seen as the most sensitive to future changes.

Modelling the effects of climate and land-use change on the hydrochemistry and ecology of the River Wye (Wales)

Interactions between climate change and land use change might have substantial effects on aquatic ecosystems, but are still poorly understood. Using the Welsh River Wye as a case study, we linked models of water quality (Integrated Catchment - INCA) and climate (GFDL - Geophysical Fluid Dynamics Laboratory and IPSL - Institut Pierre Simon Laplace) under greenhouse gas scenarios (RCP4.5 and RCP8.5) to drive a bespoke ecosystem model that simulated the responses of aquatic organisms. The potential effects of economic and social development were also investigated using scenarios from the EU MARS project (Managing Aquatic Ecosystems and Water Resources under Multiple Stress). Longitudinal position along the river mediated response to increasing anthropogenic pressures. Upland locations appeared particularly sensitive to nutrient enrichment or potential re-acidification compared to lowland environments which are already eutrophic. These results can guide attempts to mitigate future impacts and reiterate the need for sensitive land management in upland, temperate environments which are likely to become increasingly important to water supply and biodiversity conservation as the effects of climate change intensify.

Analysing the impact of multiple stressors in aquatic biomonitoring data: A 'cookbook' with applications in R

Multiple stressors threaten biodiversity and ecosystem integrity, imposing new challenges to ecosystem management and restoration. Ecosystem managers are required to address and mitigate the impact of multiple stressors, yet the knowledge required to disentangle multiple-stressor effects is still incomplete. Experimental studies have advanced the understanding of single and combined stressor effects, but there is a lack of a robust analytical framework, to address the impact of multiple stressors based on monitoring data. Since 2000, the monitoring of Europe's waters has resulted in a vast amount of biological and environmental (stressor) data of about 120,000 water bodies. For many reasons, this data is rarely exploited in the multiple-stressor context, probably because of its rather heterogeneous nature: stressors vary and are mixed with broad-scale proxies of environmental stress (e.g. land cover), missing values and zero-inflated data limit the application of statistical methods and biological indicators are often aggregated (e.g. taxon richness) and do not respond stressor-specific. Here, we present a 'cookbook' to analyse the biological response to multiple stressors using data from biomonitoring schemes. Our 'cookbook' includes guidance for the analytical process and the interpretation of results. The 'cookbook' is accompanied by scripts, which allow the user to run a stepwise analysis based on his/her own data in R, an open-source language and environment for statistical computing and graphics. Using simulated and real data, we show that the recommended procedure is capable of identifying stressor hierarchy (importance) and interaction in large datasets. We recommend a minimum number of 150 independent observations and a minimum stressor gradient length of 75% (of the most relevant stressor's gradient in nature), to be able to reliably rank the stressor's importance, detect relevant interactions and estimate their standardised effect size. We conclude with a brief discussion of the advantages and limitations of this protocol.

Functional redundancy as a tool for bioassessment: A test using riparian vegetation

There is an urgent need to track how natural systems are responding to global change in order to better guide management efforts. Traditionally, taxonomically based metrics have been used as indicators of ecosystem integrity and conservation status. However, functional approaches offer promising advantages that can improve bioassessment performance. In this study, we aim to test the applicability of functional redundancy (FR), a functional feature related to the stability, resistance and resilience of ecosystems, as a tool for bioassessment, looking at woody riparian communities in particular. We used linear mixed-effect models to investigate the response of FR and other traditional biomonitoring indices to natural (drought duration) and anthropogenic stress gradients (flow regulation and agriculture) in a Mediterranean basin. Such indices include species richness, a taxonomic index, and the Riparian Quality Index, which is an index of ecological status. Then, we explored the ability of FR and the other indices to discriminate between different intensities of human alteration. FR showed higher explanatory capacity in response to multiple stressors, although we found significant negative relationships between all the biological indices (taxonomic, functional and ecological quality) and stress gradients. In addition, FR was the most accurate index to discriminate among different categories of human alteration in both perennial and intermittent river reaches, which allowed us to set threshold values to identify undisturbed (reference condition), moderately disturbed and highly disturbed reaches in the two types of river. Using these thresholds and the best-fitting model, we generated a map of human impact on the functional redundancy of riparian communities for all the stretches of the river network. Our results demonstrate that FR presents clear advantages over traditional methods, which suggests that it should be part of the biomonitoring toolbox used for environmental management so as to obtain better predictions of ecosystem response to environmental changes.

Similarity in the difference: changes in community functional features along natural and anthropogenic stress gradients

The effect of stressors on biodiversity can vary in relation to the degree to which biological communities have adapted over evolutionary time. We compared the responses of functional features of stream insect communities along chronic stress gradients with contrasting time persistence. Water salinity and land use intensification were used as examples of natural (long-term persistent) and anthropogenic (short-term persistent) stressors, respectively. A new trait-based approach was applied to quantify functional diversity components and functional redundancy within the same multidimensional space, using metrics at the taxon and community levels. We found similar functional responses along natural and anthropogenic stress gradients. In both cases, the mean taxon functional richness and functional similarity between taxa increased with stress, whereas community functional richness and functional redundancy decreased. Despite the differences in evolutionary persistence, both chronic stressors act as strong nonrandom environmental filters, producing convergent functional responses. These results can improve our ability to predict functional effects of novel stressors at ecoloiical and evolutionary scales.

Effects of repeated salt pulses on ecosystem structure and functions in a stream mesocosm

Rivers and streams affected by mining activities often receive short-term sharp salinity increases due to water-soluble stockpiled materials being washed into receiving water bodies. We conducted a mesocosm study to explore the response of structural (diatom and stream invertebrate communities) and functional descriptors (chlorophyll a concentration, fungal biomass and leaf decomposition) to repeated short salinity pulses (3h of duration, with nominal electrical conductivities of 5, 10 and 15 mS cm(-1)), mimicking the exposure pattern occurring at salt-mine affected rivers. The experiment was conducted in 12 artificial flow-through stream systems over 16 days. The effect of the salt pulses on the ecosystem structure and functioning did not fully match most of our initial hypotheses, with the community response being weaker than predicted. The diatom community was, however, dominated by salt-tolerant species throughout the experiment, showing no consistent response to the treatment. The invertebrate response was associated with statistically significant changes in community structure (i.e. abundance of the different taxa) but no statistically significant changes in taxa richness. The salt pulses affected some functional descriptors of the ecosystem: fungal biomass exhibited a unimodal response to treatment magnitude, algal growth (i.e. chl a biomass) was hampered with increasing conductivity and leaf decomposition was significantly reduced in the high treatment.