Next Article in Journal
Epiphytic Diatom-Based Biomonitoring in Mediterranean Ponds: Traditional Microscopy versus Metabarcoding Approaches
Next Article in Special Issue
Swin Transformer for Complex Coastal Wetland Classification Using the Integration of Sentinel-1 and Sentinel-2 Imagery
Previous Article in Journal
Special Issue: Environmental Fate of Contaminants in the Aquatic Environment
Previous Article in Special Issue
How does Urban Pollution Influence Macroinvertebrate Traits in Forested Riverine Systems?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Insecticides and Drought as a Fatal Combination for a Stream Macroinvertebrate Assemblage in a Catchment Area Exploited by Large-Scale Agriculture

1
South Bohemian Research Center of Aquaculture and Biodiversity of Hydrocenoses, Faculty of Fisheries and Protection of Waters, University of South Bohemia in České Budějovice, Zátiší 728/II, 389 25 Vodňany, Czech Republic
2
Povodí Labe, Víta Nejedlého 951/8, Slezské Předměstí, 500 03 Hradec Králové, Czech Republic
3
Institute of Environmental and Chemical Engineering, Faculty of Chemical Technology, University of Pardubice, Studentská 573, 532 10 Pardubice, Czech Republic
*
Authors to whom correspondence should be addressed.
Water 2021, 13(10), 1352; https://doi.org/10.3390/w13101352
Submission received: 12 March 2021 / Revised: 4 May 2021 / Accepted: 5 May 2021 / Published: 13 May 2021
(This article belongs to the Special Issue Freshwater Communities in Human-Altered Ecosystems)

Abstract

:
This case study documents responses in a headwater macroinvertebrate assemblage to insecticide pollution and hydrological drought. In 2014, the Doubravka brook (Czech Republic) was damaged by a large overflow of a mixture of chlorpyrifos (CPS) and cypermethrin (CP). In 2016–2017, this brook was then affected by severe drought that sometimes led to an almost complete absence of surface water. We found significant relationships between the strength of both these disturbances and the deeper taxonomic levels of both the overall macroinvertebrate assemblage (classes) and the arthropod assemblage alone (orders and dipteran families), as well as the functional feeding groups (FFGs). The CPS-CP contamination was mostly negatively correlated to arthropod and non-arthropod taxa and was positively correlated only with FFG collector-gatherers; on the other hand, the drought was negatively correlated to Simuliidae, Ephemeroptera, Trichoptera, and the FFG of grazer-scrapers and passive filterers. Drought conditions correlated most positively with Isopoda, Ostracoda, Heteroptera, adult Coleoptera, and predator and active filterer FFGs. The chosen eco-indicators (SPEARpesticides, SPEARrefuge, BMWP, and EPT) used as support information reveal the poor ecological status of the whole assemblage, including the control site, the cause of which is most likely to be the exploitation of the adjacent catchment area by large-scale agriculture. This type of agricultural exploitation will undoubtedly affect macroinvertebrate assemblages as a result of agrochemical and soil inputs during run-off events and will also exacerbate the effect of droughts when precipitation levels drop.

1. Introduction

Macroinvertebrate assemblages inhabiting headwaters in a cultural landscape have to cope with exposure to anthropogenic activities and their implications [1,2,3,4,5]. The vast majority of published studies generally associate human activities with negative impacts on stream invertebrate biota and their ecological functioning; that is, they provoke a fall in diversity, abundance, biomass, and organic matter processing [6,7,8,9,10,11,12]. Nevertheless, human-triggered effects can result in greater biomass of specific aquatic organisms [13,14,15] or, occasionally, in higher species richness [16,17], and some ostensibly heavily affected sites may sometimes even act as refuges for endangered species [18,19,20].
Over the last few years, much of Europe has suffered from severe droughts in combination with unusually high temperatures [21,22], and climatological prognoses suggest that these events will become more frequent in the future [23,24]. Combined with intensive land use and its effects on functioning water regimes, climate change could cause huge ecological and economic damage [25,26,27,28]. Affected headwaters will dry up more often and more quickly, thereby multiplying the impact of pollution by agricultural or urban wastewater [29,30,31]. Aside from the plethora of pollutants that occur in surface waters, pesticides (e.g., insecticides, herbicides, and fungicides) are of concern due to the negative effect they have, in particular, on aquatic biota and, in general, on whole ecosystem functioning [32,33,34,35]. The main reason is their massive use worldwide and their typically high specific toxicity for non-target aquatic organisms [5,36,37,38]. Pesticides are an integral part of conventional agriculture and enter the water at much higher concentrations after short-term runoff events than those that are usually detected by standard monitoring sampling methods [39,40]. Therefore, pesticides or their residues may have acute or chronic effects on aquatic organisms [41,42]. Furthermore, toxicants with great adsorption to soil or organic matter may persist in sediments; hence, they can threat sediment dwelling organism for a prolonged time [10,43]. The presence of pesticides in aquatic ecosystems is often accompanied by a wide range of anthropogenic influences, such as habitat degradation, drainage [44,45], artificial siltation and sedimentation [46,47], increased nutrient input, and the occurrence of more frequent extreme hydrological events (droughts and floods) at interconnected sites [44,48]. These factors often exacerbate the negative effects of pollutants on aquatic biota, generally due to the losses they induce in many kinds of resources in an ecological sense (e.g., microhabitats, refuges, food, and the strength of interactions) [49,50,51].
Invertebrate populations should be able to withstand specific disturbances whose impact will presumptively be reversible. Nevertheless, certain regularly recurring processes in combination with other biotic or abiotic conditions will undoubtedly damage exposed assemblages by eliminating the most susceptible species (e.g., Margaritifera margaritifera or Prosopistoma pennigerum) or even entire deeper taxonomic groups (e.g., heptageniid mayflies, perlid stoneflies, and goerid caddisflies), thereby eventually negating ecological traits [4,49,52,53,54,55,56,57]. For this reason, laboratory and field experiments such as those defined by Diamond [58] that manipulate the effects of xenobiotics under various experimental settings (simulation of diverse possible natural scenarios) are necessary [11]. However, experiments carried out in natural scenarios (natural experiments) are also essential as they can verify obtained knowledge in real environments [59].
The aim of this study was thus to assess the effects of serious insecticide contamination on a headwater macroinvertebrate assemblage and observe the ability of species to recolonize affected stretches of the damaged water course. We used a dataset that was not created purely for research purposes; rather, it was used by the state environmental protection authorities as part of their legal investigation into the accident. Thus, it lacks some of the information needed for a thorough assessment of certain aspects of the macroinvertebrate assemblage. Our study reveals the response of this macroinvertebrate assemblage to an extreme event, a situation that is rarely observed, probably because—unlike chronic effects—these accidents only occur rarely. In terms of climate change and evolving strategies in agriculture, our results reflect a situation that will increasingly threaten organisms in watercourses in agriculture landscapes, as well as all other lifeforms that are dependent on these damaged water bodies.

2. Materials and Methods

2.1. Locality Description

The Doubravka brook is situated in central Bohemia (Czech Republic, Central Europe, spring to confluence: N 49.7775686, E 15.5793436–N 49.8619361, E 15.4977614; Figure 1). This small, third-order stream (the Stahler order used) is 13.8 km long, has a catchment area of 21.6 km2 and average discharge of up to 9 L s−1. It rises at an altitude of 460 m a.s.l. and flows into a fifth-order stream at 251 m a.s.l. The riverbed of the brook is natural and it is surrounded by typical alluvial vegetation, together with part of a cultural coniferous forest (in all, an approximately 150–600 m wide bio-corridor). This sector of the brook has been preserved due to the steep terrain. Most of its catchment area is occupied by uniform blocks of fields (each with an average area of 30 ha) where rapeseed, corn, and cereals are cultivated. The fields are drained by underground water collectors connected to straight channels, as well as to artificial and/or originally temporary or permanent tributaries. The brook ecosystem harbors microhabitats with shallow riffles with macro-, meso-, micro-lithal, and pools anchored by the root systems of alder trees (Alnus glutinosa, A. incana). In these microhabitats detritus, particulate organic matter (POM), smaller-sized pieces of gravel and xylal accumulate. However, in recent years this locality has been affected by severe droughts leading to zero surface discharge, a phenomenon that was often observed during the study period, mainly during extremely hot summers. The riffles often disappeared and the bottom of the isolated pools became almost completely covered by silt and organic matter.

2.2. Insecticide Contamination

In March 2014, the brook was accidentally contaminated by a commercial insecticide product containing organophosphate chlorpyrifos (CPS) and pyrethroid cypermethrin (CP) (at a ratio of 10:1, respectively). It flowed down a straight drainage channel receiving wastewater and runoff-water from a local agricultural company that had illegally stored unused insecticide. The concentration of CPS in the sediment dry weight (d.w.) was 13 mg kg−1 at the confluence of the brook and the drainage channel (Site 1, Figure 1) five days after the accident. Mass mortality of over 10,000 individuals of the critically endangered noble crayfish Astacus astacus, macrozoobenthos, and fish (brown trout Salmo trutta and stone loach Barbatula barbatula) was observed for over 6 km downstream from Site 1. A. astacus and B. barbatula samples were analysed for CPS (77 and 31,000 μg of CPS kg−1 of body weight, respectively). CPS sediment concentrations decreased rapidly over the following months (Table 1). The bottom of the contaminated stretch of the brook also became abnormally overgrown by filamentous algae and periphyton during the first months after the accident.

2.3. Sediment Sampling

Sediment was sampled according to the standardized accredited methods (ČSN EN ISO 5667-1, ČSN EN ISO 5667-3, ČSN ISO 5667-12, ČSN ISO 5667-14, and ČSN ISO 5667-15) by a telescopic sampler with a wide-neck stainless steel container. A layer of the bottom sediment was scrapped using the edge of the container. Regarding the character of the riverbed containing a small amount of fine sediment, sediment was sampled in several places where it was possible within one sampling site. These individual sub-samples were homogenized into one mixed sample, put inside polyethylene sampling bottles, and stored in a polystyrene thermo box filled with ice during transport prior to the laboratory analysis.

2.4. Analysis for Pesticides

2.4.1. Reagents and Materials

The analytical standards of chlorpyrifos (chlorpyrifos-ethyl, CAS 2921-88-2) and deuterated chlorpyrifos-D10 (CAS 285138-81-0) (Dr. Ehrenstorfer; Augsburg, Germany) were of 99.49% and 99.1% purity, respectively. Methanol and mobile phase additive ammonium acetate were LC-MS purity (Merck; Darmstadt, Germany). An SG Water Ultra Clear TWF UV plus TM with TOC-Monitoring water system from SG Water (Hamburg, Germany) was used throughout the study to obtain the LC-MS grade water.

2.4.2. Sample Processing and Analyses

The ultrasonic assisted extraction with methanol solvent was applied for the preparation of sediment samples as it was described by Ferenčík and Schovánková [60]. The samples were dried at 20 °C using a lyophilizer Christ Alpha 1–4 (Osterode, Germany), and sifted through a 2 mm Retsch sieve (stainless steel). The 10 mL aliquot of each sample was homogenized by grinding using a mixer mill MM 200 from Retsch (Haan, Germany) and zirconium oxide 25 mL mixing jars with balls. 5 µL of the internal standard solution (chlorpyrifos-D10, 100 pg mL−1) were spiked into 0.5 g homogenized subsamples, then extracted using 5 mL of methanol in 50 mL glass Erlenmeyer flasks in ultrasonic bath Binder Electronic (Berlin, Germany) for 60 min at 40 °C (maximum ultrasonic power). The extracts were centrifuged by centrifuge Jouan B4I (Thermo Electron Industries, Chateau-Gontier, France) at 2400 g. The 100 µL aliquot of extract was diluted with 900 µL of 100 mmol L−1 ammonium acetate in a 2 mL sample vial used for LC-MS/MS measurement. A smaller aliquot of a sample or higher extraction volume was selected for the more contaminated samples.

2.4.3. LC-MS/MS Conditions

Based on published literature [61], development and optimization of the LC-MS/MS method was done. The chromatographic separations were conducted using a Waters Acquity UPLC HSS T3 column (100.0 × 2.1 mm, 1.8 µm particle size) and guard column 2.1 × 5.0 mm with the same chemistry, thermostated at 40 °C. 5% methanol and 5 mmol L−1 ammonium acetate in water (A) and 5 mmol L−1 ammonium acetate in methanol (B) were used as mobile phases. Gradient elution program was set up from 0.1% B to 99.9% B, the method duration was 18.5 min. The injection volume was 250 µL.
The LC-MS/MS analysis was carried out using triple quadrupole Waters Premier XE (Manchester, UK) connected with ultra-high-performance liquid chromatograph Waters Acquity UPLC® (Milford, MA, USA). The electrospray ionization source heated at 120 °C in the positive ionization mode (3.5 kV) was used. For acquisition and processing of data, the Waters MassLynx software was used.
Chlorpyrifos SRM transitions (precursor > product) were as follows: quantifying 349.8 > 197.7 (cone voltage—25 V and collision energy—19 V), qualifying 349.8 > 96.7 (cone voltage—25 V and collision energy—34 V). In addition, internal standard chlorpyrifos-D10 transitions were as follows: quantifying 359.8 > 198.7 (cone voltage—25 V and collision energy—19 V).

2.4.4. Validation of the Analytical Procedure

Six points’ calibration curve was linear in the range of 10–500 ng mL−1. Internal standard and matrix matching standard methods were used for quantification of chlorpyrifos. Recoveries of chlorpyrifos in sediments were calculated by spiking the matrix at two concentration levels—50 and 500 ng g−1—in triplicates. The range of recoveries of the spiked samples was 78–121%. The repeatability of the method was determined as relative standard deviation (RSD) of repeating analysis of spiked samples and was lower than 20%. The method detection and quantitation limits of chlorpyrifos was expressed based on S/N ratio of 3 and 10, respectively [61]. The limit of quantitation (LOQ) was calculated to 20 µg kg−1 (d.w.) with a measurement uncertainty of 30% covering the whole sample preparation and measurement procedure assessed using control chart (2 s) constructed at a concentration level 50 µg kg−1.

2.5. Macrozoobenthos Sampling and Analysis

2.5.1. Sampling and Determination

Monitoring took place in 2014–2017. Three sites along the brook were chosen (Figure 1). Site C (control) was not contaminated by the insecticides. Sampling of macrozoobenthos was performed using kick sampler (25 × 25 cm net frame dimensions, 500 µm mesh) by 3 min-long sampling multihabitat method, following rules of European Water Framework Directive (WFD) [62]. Briefly, this method is based on an a priori evaluation of presence and spatial proportions of various microhabitats at the chosen sampling site. Subsequently, an interval taken from the total time (3 min) is allotted for each sampling at particular chosen micro-habitats regarding to their spatial proportionality compared with the other ones. The sub-samples of micro-habitats were always pooled into only one sample, processed using a round steel sieve (40 cm diameter, 500 µm mesh) and preserved in 70% technical ethanol. Organisms were determined to the highest possible taxonomic level; their frequencies in the samples were recorded (Table S6). Before determination, single portions of material were individually separated from smaller particles using a small plastic sieve (500 µm mesh) and tap water.
The first sampling was performed in June 2014 (approximately 2.5 months after the accident). The others in 2015–2017 were restricted between April and September (Table 1). The sampling was initiated by the state monitoring authorities investigating the accident. Therefore, the experimental design lacked replicates of macrozoobenthos samples at individual Sites and timepoints and relied instead on a multihabitat sampling approach.

2.5.2. Estimation of Taxa Traits

A dataset containing the ecological and plasticity traits of all detected taxa was created (Table S7). The categories of these ecological traits are represented by the following five subsets: Functional feeding groups (Collector-gatherers, Shredders, Xylophages, Active filterers, Passive filterers, Predators, and Grazer-scrapers,) Macro-habitat preference (Limnophiles, Rheophiles, and Torrenticoles), Meso-habitat preference (Epibenthos, Endobenthos, and Hyponeuston), Micro-habitat preference (Xylal—tree trunks, branches, and roots; Akal—fine to medium-sized gravel; POM—deposits of fine particulate organic matter; Pelal—mud and sludge; Detritus—deposits of coarse particulate organic matter;, Psammal—sand; Lithal—coarse gravel, cobbles and blocks; Phytal—algae, mosses, and plants) and Reproductive ability (<0.5 generation per year, 0.5 generation per year, 1 generation per year, 2 generations per year, and >2 generations per year). The plasticity trait used was Body length of adult or final larval instar (<4 mm, 4–7 mm, 8–12 mm, 13–18 mm, 19–25 mm, and >26 mm). The classification was empirically encoded by the proportion of a particular type of trait for detected taxa within each subset category.

2.5.3. Calculating Eco-Indicator Parameters

The base parameters (total richness, total abundance, Shannon–Wiener index) commonly used to describe a macroinvertebrate assemblage status were generated in R studio software using the BiodiversityR package [63] for all samples. Three types of ‘Species At Risk’ (SPEAR) classifications (coded by dummy variables) were used to reveal the proportions and sums of the abundances of the most susceptible species, as well as and their richness, as follows: (i) ‘SPEARpesticides (old)’ were classified taking into account certain lethal concentration levels and species-specific ecological traits that could be critical for the organism exposed to the pesticides [40]; (ii) the ‘SPEARrefuge’ concept represents the adjustment of the original SPEARpesticides list with a category of sensitive species with a strong potential for successfully recolonising contaminated stretches of streams and brooks [49]; and (iii) ‘SPEARpesticides’ are a subset of the ‘SPEARpesticides (old)’ that do not belong to the ‘SPEARrefuge’. SPEAR indices were calculated according to the following adjusted equation:
SPEAR   index = i = 1 n log ( 4 x i + 1 ) × y i = 1 n log ( 4 x i + 1 )
where n is the total number of taxa, xi is the abundance of taxon i, and y is a binary code representing taxon i thus: 1 if it belongs to the SPEAR, otherwise 0. The log (4x + 1) transformation of the abundances aims to decrease the influence of populations with mass development [49].
EPT (Ephemeroptera–Plecoptera–Trichoptera) richness, abundance, and relative abundance (%) were calculated because many species groups belonging to these three orders are often very sensitive to several types of disturbances including insecticide pollution. EPT-derived parameters are frequently used as a measure of the quality of the environment and/or the seriousness of a disturbance effect by comparing cases. The EPT abundance usually correlates with the SPEAR abundance.
The original biological monitoring working party (BMWP) score and its average score per taxon (ASPT), which are based on family presence rather than on abundance, was also derived [64]. The aim of the BMWP score is to indicate the status of organic pollution. Families sensitive to organics are usually more susceptible to insecticides.

2.6. Classification of Hydrological Status

The hydrological status semi-quantitative scoring system for implementation in the ordination analysis was designed following Boulton and Lake [53] (Table 2).The prior hydrological status at each sampling point was based on our observations (Figures S3–S5), precipitation and temperature data from previous months for the particular area, and discharge data in the river Doubrava—a second-order river receiving water of the Doubravka brook (data taken from the website of the Czech hydrometeorological institute, Prague, Czech Republic, Table S4 and Figure S2). In 2015–2017, abnormally high temperatures in this area were accompanied by a lack of precipitation. The scores assigned to each sample are shown in Table S5.

2.7. Statistical Analysis

The dataset containing the abundance of individual taxa as response variables and two environmental predictor variables Chlorpyrifos-Cypermethrin contamination status (CPS-CP) and drought status (Drought), and one covariate (Transect), were analysed using the multivariate Canonical correspondence analysis (CCA) in CANOCO 5 software [65]. Similarly, the chosen deeper taxonomic levels for all detected taxa and arthropods (representing the most diverse and, simultaneously, the most sensitive groups to insecticides) were used as a response in the redundancy analysis (RDA)—classes and orders (dipterans divided into families), respectively. The variables CPS-CP and taxa abundance were log-transformed; the response variables were centered but not standardized. Additionally, the ecological and biological traits of the invertebrates (see Section 2.4.2. Estimation of species traits) were implemented as standardized composition-weighted trait averages (using the abundance of taxa in the samples) into independent RDA against the CPS-CP and Drought variables. The significance of the relationships between these environmental variables and individual compositions of taxa or traits was tested using particular canonical axes and their eigenvalues with a Monte Carlo permutation test (number of permutations: 1999). The cases were always shifted into two permutation blocks, defined by the Transect covariate, to maintain autocorrelation. Principal component analyses (PCA) were used for the trait compositions, which did not reveal any significant relationship with the explanatory variables.

3. Results

3.1. Recolonisation of Poisoned Stretches by Sensitive Insect Species

Eighty-three days after the accident, the sample from non-polluted Site C contained 24 insect taxa and one oligochaete taxon. At polluted Site 1, the species richness was reduced to eight taxa (two oligochaetes and six insects) (Table 3). Site 1 was recolonized by two species classified as SPEARpesticides (Ecdyonurus sp. and Rhyacophila nubila) and three species considered as SPEARrefuge (Ephemera danica, Halesus digitatus, and Hydropsyche instabilis) after the poisoning. The two remaining sensitive insect taxa belonged to the Chrironomidae family, represented by an early Chironomidae gen. sp. instar and Brillia bifida.
At Site 2, 10 Insecta taxa, 3 Oligochaeta taxa and 2 Gastropoda taxa were detected. Two Diptera families (Chironomidae: Chironomidae gen. sp., Chironomus riparius gr., Corynoneura sp., Micropsectra sp., Paratrichocladius rufiventris, Thienemannimyia sp. and Tvetenia verralli and Simuliidae: Simulium vernum) dominated the insect assemblage. The remaining two detected insect taxa, both caddisflies classified as SPEARpesticides (Rhyacophyla sp.) and SPEARrefuge (Hydropsyche sp.), were detected in very low numbers (one and eight individuals, respectively).
The SPEARpesticides and SPEARrefuge detected at Site C but not in the poisoned stretches were Ephemerella mucronata, Habrophlebia lauta, Serratella ignita, Leuctra sp., Nemoura sp., Hydropsyche siltalai, Polycentropus flavomaculatus and Potamophylax luctuosus. The species not at risk (SPEnotARpesticides) that were not found to colonise at least one of the two contaminated Sites were Microtendipes chloris gr., Orthocladius sp., Tanypodinae gen. sp., Tvetenia verralli, Ceratopogonidae gen. sp., Dicranota sp., Tipula maxima and Elmis sp. lv.

3.2. Relationship between Disturbances and Macroinvertebrate Taxonomic Composition

A Monte-Carlo permutation test performed within independent RDA revealed the significant power of both variables (adjusted p < 0.05) for predicting the taxonomic compositions of the whole macroinvertebrate assemblage, which could be classified into seven deeper taxonomic levels (mostly classes) (Figure 2A), and the arthropod assemblage classified into 19 deeper taxonomic groups (mostly orders) (Figure 2B). Sediment contamination by insecticides (CPS-CP) negatively correlated with all the deeper taxonomic levels (in terms of their abundance) other than with Platyhelminthes (Polycelis nigra) (Figure 2A). Platyhelminthes, Crustacea, Bivalvia, Gastropoda and Insecta were the most reliably fitted groups. Hence, sensitivity to contamination, from the most to the least affected taxonomic groups, is as follows: Bivalvia (only Pisidium casertanum and P. personatum) and Gastropoda < Insecta and Crustacea (Ostracoda and Isopoda) < Platyhelminthes (Polycelis nigra) (no effect). The remaining two annelid taxa (Hirudinea and Oligochaeta) did not reach the same quality of fit; nevertheless, the Oligochaeta group does seem to be sensitive to CPS-CP contamination. The increasing intensity of drying up positively correlated with densities of all given taxa, but mostly with the Bivalvia species, Crustacea, Hirudinea, and Oligochaeta.
Predominantly lotic Heteroptera and adult Coleoptera were the commonest insect taxa in assemblages found at an advanced stage of drying up (Figure 2B). Conversely, the family Simuliidae was highly sensitive to drought; the EPT group was negatively correlated with both factors. This group showed relatively higher resilience to contamination compared to the other arthropods; however, it was almost completely lacking in samples taken under the worst hydrological conditions (Table 3).
The test performed using the CCA to ordinate the whole assemblage into the highest possible taxonomic level revealed the marginal effect of both predictors (adjusted p = 0.055 for the effect of both disturbances).

3.3. Effect of Disturbances on Ecological and Biological Species Traits

Of the ecological and plasticity traits described in Materials and Methods (see Section 2.4.2 Estimation of Species Traits), only the Functional feeding groups for all detected taxa (Figure 3) showed a significant relationship with the first ordination axis and both predictors (adjusted p < 0.05) within the partial RDA. The CPS-CP contamination negatively correlated with the abundance of predators and active filterers, both of which, however, prevailed over the other feeding strategies during drought conditions (Figure 3). This trend was caused by increased densities in, especially, chironomid, ostracod and bivalvian assemblages in isolated pools: namely, the facultative active filterers from the Tanytarsini tribe (Chironominae) such as Micropsetra apposita and obligatory active filterers represented by Ostracoda and Bivalvia (Pisidium sp.). Numerous predators or facultative predators from the Tanypodinae subfamily including Macropelopia nebulosa, Procladius (Holotanypus) sp., and Thienemannimyia sp. were followed by the nectonic heteropteran Notonecta sp. and the adults of the coleopterans Platambus maculatus and Agabus didymus during drying up conditions. On the other hand, contamination correlated positively with collector-gatherers for a short time period; nevertheless, this strategy was suppressed during the hydrological droughts, as occurred with passive filtration and grazing-scraping. The shredders and xylophages were the least involved in material processing during both peaks of the perturbations (Table S3).
The remaining compositions of species traits extrapolated from taxa abundance did not significantly correspond with the variables. Nevertheless, relationships between disturbances and compositions of species traits (listed in the Section 2.4.2 Estimation of Species Traits) are shown in the diagrams resulting from the PCA (Figure S1).

4. Discussion

4.1. Response of Macroinvertebrate Taxa to the Disturbances

Macroinvertebrate species assemblages can substantially change after acute poisoning with insecticides and during prolonged periods of drought [53,66,67]. The significant relationship plotted in Figure 2A suggests that no groups of organisms prefer the CPS-CP contamination, and that only the Platyhelminthes (Polycelis nigra) tolerate it. Contrary to our prediction (Table S1) that the stream macroinvertebrate assemblage would shift from being dominated by arthropods to being dominated by non-arthropod species, the insect and crustacean groups, in general, did not seem to be the most affected groups [67]. The most likely cause of this effect is that the downstream part of the brook was colonized by species from the upper non-affected part. In addition, the insect species with the ability to drift downstream or with winged imagines (ephemeropterans, plecopterans, and chironomids) are often perfectly adapted for colonizing new habitats, and so insects may in fact be very resilient to insecticides [34]. These adaptations may be much more advantageous than a resistance to a prolonged or severe disturbance [68,69]. However, if species populations are too weak or non-existing in adjacent brooks or catchments, they can become temporarily or permanently extinct [70].
Hydrological conditions—and especially droughts—heavily influenced the macroinvertebrate assemblage. The abundance of all groups correlated more or less positively with an increase in droughts (Figure 2A). Bivalves (P. casertanum and P. personatum), gastropods, oligochaetes, crustaceans (Asellus aquaticus and Ostracoda), as well as leeches, all showed greater correlation than insects (bivalves and gastropods may have increased their relative density due to a decrease in available space during droughts). Furthermore, the drought period was accompanied by increased sedimentation of organic matter caused by very low water discharge. A further meaningful biological factor that could be partially responsible for both the redundant accumulation of organic matter and assemblage shift is the total extinction of A. astacus in this brook. Before the accident, noble crayfish were found at a density of around one individual per m2 at Sites 1 and 2, and its population did not recover at any point during this study (no individuals caught in traps) [71]. These bottom-dwellers shred and release leaf litter from under obstacles on riverbed into the current, dig burrows in riverbanks, and also feed on macrozoobenthos, presumably slow-movable macroinvertebrates [72].
The Crustacea group represented by the isopod A. aquaticus and the class Ostracoda (Figure 2A) appeared at all Sites during the dry period in 2016–2017, as did the Bivalvia. The relationship between CPS-CP and this group was similar to the relationship between CPS-CP and the insect group, which demonstrates the high specific toxicity of these compounds to arthropods in general [73,74]. Nevertheless, the drought had a generally positive effect on these two crustaceans (Figure 2A,B). The relationship between insecticide pollution and the abundance of A. aquaticus within ecosystems that are often affected by toxicants as well as organic pollution remains unclear [40]. Although A. aquaticus is very sensitive to acute exposure to insecticide, it can tolerate lower concentrations of various toxic compounds, including insecticides [75,76,77]. It is also able to avoid desiccation [78]; furthermore, Extence [79] reported a significant increase in the abundance of A. aquaticus in a lowland river during a drought period as a result of greater accumulation of organic matter. The response of A. aquaticus to these two types of disturbances can be summed as intolerance to acute poisoning but a general preference for an ecosystem disturbed by partial drying up and the increased deposition of organic material that can serve as a food resource.
The relationships between arthropod groups and disturbances shown in Figure 2B suggests that very few taxa were positively correlated with the CPS-CP contamination; however, the individual response to greater hydrological drought was more complex. Members of the family Simuliidae were the most negatively correlated with greater drought (and had a positive relationship with CPS-CP contamination), a logical finding given that, due to their feeding adaptations, they are dependent on flowing water current [80]. We believe that the positive relationship with increasing CPS-CP contamination is due to their ability to recolonize contaminated stretches of brook if favourable hydrological conditions are present (i.e., flowing water) [81]. The general intolerance of many species from the EPT group to pesticides, organic pollution, and a lack of water (factors that are typical in water courses affected by large-scale agriculture) has been well documented and is discussed in more detail below (Section 4.2) [82,83,84].
Taxonomic groups ordinated near the Chironomidae family were represented by the larvae of other dipteran species and coleopteran larvae and adults. This group of taxa seems to be very sensitive to acute pollution; however, it colonised the most contaminated Sites immediately after the mayflies and caddisflies and, unlike them, were tolerant to drought. This also occurred in the families such as the Dytiscidae that colonising shrinking pools. Chironomids are the most complex family and individual species have very plastic ecological characteristics [85,86,87,88]. Their high susceptibility to CPS-CP and other insecticides under laboratory conditions does not correspond to the generally high resilience of many Chironomid species to various types of disturbances [89,90]. However, the preferences of chironomid larvae for certain hydrological conditions varies between species within subfamilies or even within a single genus. For instance, the subfamily Orthocladiinae decreased in abundance and richness during the dry period, whilst the Tanypodinae and Tanytarsini (Micropsectra sp.) dominated in the samples during drought conditions in 2016 (Table S2). The high densities of Micropsectra sp. could be a product of the stress caused by the drought, and numerous tanypodins would either feed on them or be similarly forced to move into refuges [84,91]. It is worth mentioning the fact that stream macroinvertebrates often concentrate in small wet areas during periods of drought, and that they can increase or decrease their densities depending on taxon-specific life histories [54]. This system may be much more at risk to different types of pollution and other interferences, and higher densities and weakening may favor the spread of disease within assemblages [92].

4.2. Response of Species-Trait Assemblages to Disturbances

Taxonomic and functional feeding group (FFG) assemblages often have intimately interrelated characteristics that change along environmental gradients and may react in different ways to disturbances [54]. Less time is usually required for the recovery of trophic structures than for the recovery of taxonomic assemblage structures [93]. Although shredders are considered to be the most susceptible group to insecticides, and despite the fact that they are often also sensitive to droughts, the approximal correlation between shredders and both drought and CPS-CP was close to zero [53,67]. Almost identical trends were estimated for xylophages. The relative abundance of shredders and xylophages was low at all Sites (average proportion to other FFGs: 6.96% and 1.39%, respectively), whilst collector-gatherers dominated in most cases (average proportion to other FFGs: 50.01%; Table S3). This could have been caused by the time of sampling (April–September) because shredder abundance commonly increases in autumn when more leaf litter is available [94]. However, if the whole assemblage were chronically influenced by insecticides at all Sites, an arthropod-shredder population, represented especially by caddisfly larvae (no gammarids were detected during the study), would be suppressed. Specifically, the shredder FFG was represented by the caddisfly H. digitatus, which prefers permanent waters, at the contaminated Sites after the accident, whilst the isopod A. aquaticus, which is considerably more tolerant to drought conditions, was the most abundant representative of the facultative shredding strategy during the period when the brook was drying up [78,95]. The exclusion of shredders from headwaters is thought to trigger ecological changes in overlapping assemblages due to the crucial role they play in the smooth cycling of nutrients by processing coarse particulate organic matter into smaller transportable particles (fecal pellets and orts). The headwater streams lacking shredders are colonized by saprotrophic microorganisms (bacteria, fungi, and protozoa), followed by microphagous collector-gatherers and filterer macroinvertebrates. This type of assemblage is not very resistant to the highly turbulent water flow found in headwaters after heavy precipitation [25,67,96,97].
The positive relationship between drying up and predator abundance has been described by Boulton and Lake [53], who named this phenomenon as ‘predator soup’. It occurs as organism density increases in shrinking pools. Herbst et al. [98] reported an expansion of micro-predators (e.g., Tanypodinae, and Ceratopogonidae) that was partially the case in our study. The possible proliferation of an active filterer FFG assemblage in temporary standing water (e.g., Culicidae) was also suggested [53]. Nevertheless, a correct interpretation is more difficult because part of the active filtering organisms, e.g., Tanytarsini chironomids, would be forced to move into residual pools with negative consequences for their survival [84,91]. Passive filterers rely on water currents; thus, this feeding strategy may be useless when stream flow ceases. Nevertheless, within a species this strategy is usually combined with other options for getting food. Similarly, the grazer-scrapers, being the second most dominant FFG after collector-gatherers, were negatively correlated with the drought variable. Despite the fact that this strategy could co-exist with others in a single species, many studies describe it as existing on its own without any other feeding mechanism, e.g., in certain heptagenid nymphs, gastropods, etc. [99]. Herbst et al. [98] reported a decrease in the relative abundance of grazers during a severe drought. According to our results, the abundance of grazers correlated negatively with the Drought variable, while its correlation with CPS-CP was close to zero.

4.3. Parameters of a Macroinvertebrate Assemblage at the Studied Localities

Stressor-specific eco-indicators are a useful tool for monitoring ecosystems because they reveal ecological effects and assess how single and combined stressors affect ecosystem structure and function [40,50,100]. The parameters calculated for each sampling Site (Table 3), which showed an obvious shift in a macroinvertebrate assemblage between localities, were not statistically supported. Although the multihabitat-sampling method used to assess the macrozoobenthos status is robust [62], this approach is often inadequate for a standard research analysis. However, in this particular study, the dataset was of use due to the very strong effect of the studied ecological disturbances. It is hard to interpret differences in total abundances from the data; nevertheless, fivefold lower macroinvertebrate abundance 2.5 months after poisoning in 2014 at Site 1 compared to Site C does suggest that macroinvertebrate density did not recover. Cuffney et al. [67] propose that approximately four months are needed for macroinvertebrate density to recover after an acute poisoning event (even though the biomass will remain significantly lower in the first year before recovering in the second [101]). However, in these authors’ study the recolonization process from an upstream section did not occur, unlike in our study. The richness, BMWP score and ASPT index estimated for Site C were regularly higher than at Sites 1 and 2, which reveals the negative impact at both Sites during the monitoring period. The BMWP approach, which considers only the richness of specific family (not abundance), indicates the presence of a higher level of saprobity at these Sites [64].
Insecticide pollution chiefly affects the most susceptible stream insect assemblages [67]. The abundance (Abu) and richness (N) of the indicator species (SPEAR and EPT) decreased at the contaminated Sites; however, the indices calculated from the proportion of the sum of the sensitive species abundances vs. the sum of all species abundance (even if log-transformed) are intended for use in cases of chronical stress, not for events of massive spills of toxicants leading to the total extermination of all species [40]. Therefore, the conclusions resulting from these indices could be misinterpreted. For example, when comparing Site C with Site 1 in 2014, both samples had almost the same value for the SPEARpesticide index (old), despite the huge contamination that had occurred there!
In terms of assemblage parameters (richness, H’, SPEAR, EPT-derived values; see Table 3), site C was the least disturbed environment—despite its similar or even worse hydrological situation—when compared to Sites 1 and 2 in 2016 and 2017. Nevertheless, the values for the SPEARpesticides index calculated for each sample, including these from Site C, reflect its ‘poor’ or ‘bad’ ecological status (according to the EU Water Framework Directive) [49,102]. These findings could be explained by the lack of any well-conserved refuge area. Sufficiently long stretches of uncontaminated headwaters situated in riparian forests are thought to act as significant functional refuges enhancing stream hydro-morphology and enabling the recolonization of a damaged downstream section [103]. However, pollution by an insecticide or any other kind of pesticides can often occur even in economically exploited forests [104].
Many authors consider that ‘edge-of-field’ runoff after heavy rain represents the main vector for pesticides entering water ecosystems [9,40,105,106]. According to Stehle and Schulz [52], the strength of their effect can be modulated by application patterns, geographical and meteorological conditions, the physicochemical properties of the insecticide (e.g., stronger sorption properties can reduce the impact) and its intrinsic toxicity. The basin of the Doubravka brook including its source is exploited agriculturally. The excessive size of local fields (up to 100 ha) drained by surface and sometimes by sub-surface systems, combined with steep slopes, undoubtedly increases the risk of pollution by agricultural chemicals, of extreme hydrological situations and of soil erosion [44,45].
Immediately after the accident (2014), total invertebrate abundance and richness were substantially lower near the contamination source. The first pioneers to colonize Site 1 were the SPEARpesticides (old) and EPT species, although there were also a number of SPEARrefuge. Despite their high sensitivity to pesticides, these species are frequently detected in contaminated stretches near uncontaminated stream sections regardless of the pesticide pressure. The occurrence of SPEARrefuge in contaminated or disturbed streams can be caused by dispersal-based resilience and the ability of these organisms to disperse from uncontaminated to contaminated stream sections (via drift and adult dispersal) [49]. Given the spatial proximity between uncontaminated Site C situated upstream from Site 1 (Figure 1) and the time lapse between the accident and the sampling time (2.5 months), the drift of these insect species is the most likely explanation of their occurrence there. In addition, the proportion of stream organisms drifting downstream from stretches damaged by repeated insecticide pollution can be much higher than in undamaged stretches [67]. The situation at all sampling Sites (low values of SPEARpesticide indices, agricultural exploitation of the basin, and the period of pesticide application between the accident and the first sampling) suggests chronic pollution by pesticides or another disturbance. Nevertheless, the lethality of short-term pollution and the level of macroinvertebrate drift from Site C is unknown; thus, recolonization via drift downstream from this Site could either be accelerated by a lower-effect disturbance, especially at the level of sensitive species richness, or inhibited by a higher-effect disturbance at this level and, especially, at a level of sensitive species abundance. It is worth mentioning that fresh generations of chironomids and Oligochaeta (the Naidinae subfamily) were detected at Site 1, although both groups are considered to be SPEnotARpesticides, and they could even have been reproducing at this Site after the accident. Since the bottom of contaminated stretches of brook became covered by large amounts of filamentous algae and periphyton during the first few months after the accident, these organisms could have prospered in this micro-habitat despite the persistence of the contamination in sediments (Figure S1B).
The distance between sampling Site C and sampling Site 2 (approximately 6 km) could represent the threshold distance needed for fast recolonisation by SPEARrefuge species [49]. This agrees with our findings because the detected abundance and richness of SPEARpesticides and SPEARrefuge species were very low here compared to Site 1, although Site 2 was not damaged as much by the gradual flow of contaminated sediment from the drainage channel. On the other hand, the sensitive insect species (SPEnotARpesticides) represented by Diptera (Chironomidae and Simuliidae) might have recolonized this stretch given their reproductive potential. Since the detected dipteran larvae were mostly bi- or multi-voltine (and in a case of later instars, above all epibenthic), survival in the hyporheic zone is less likely. In the case of caddisflies, migration from the upstream section and/or from a possible small tributary refuge situated between Sites 1 and 2 cannot be excluded.
In 2016 and 2017, the almost total absence of SPEARpesticides and EPT species at Sites 1 and 2 (Table 3) raises several possible questions about another pesticide contamination event (or another disturbance) during this period. This event was partially revealed by the increased concentration of cypermethrin recorded in sediments sampled in the drainage upstream from the confluence with the Doubravka brook. Despite this, a higher concentration of cypermethrin was not detected at the observed Sites. During this period, the area was affected by severe hydrological drought. Given that the drainage channel during the time period supplemented Sites 1 and 2 with muddy, pesticide-enriched water, this situation is a good example of the combined effects of drought, greater sedimentation and pesticide pollution (’ramp’, ‘press’, and ‘pulse’ types of disturbance, respectively, according to Lake [107]), which led to partial absences of the EPT species, especially at Site 1 (Table 3). Despite the fact that locality C was affected by the same or even greater lack of water than these two localities (Table S5), the EPT groups resisted here. Nevertheless, the redundant sedimentation caused by the lack of water current probably triggered the high development of chironomids and oligochaetes.

5. Conclusions

This case study attempts to describe how the macroinvertebrate assemblage responds to acute poisoning by insecticides (Chlorpyrifos and Cypermethrin) and severe hydrological droughts in a headwater stream that is chronically affected by agrochemicals and large-scale agriculture exploitation in its catchment area. In the case of the accidental insecticide contamination, the assemblage reacted with a drop in richness and abundance. The contaminated Sites were colonized by resilient taxa and the lotic insecticide-susceptible taxa were able to spread from non-contaminated refugial stretches situated upstream from the source of the contamination. In the case of hydrological droughts, lotic organisms declined and the assemblage shifted to favor taxa that prefer higher levels of organic pollution and lentic hypo-neustonic (nektonic) predatory insects.
A non-stable hydrological regime is a serious issue that many once-permanent aquatic ecosystems have to confront. Damage caused by agricultural pollution will be more serious for sensitive taxa if headwater ecosystems turn into semi-temporary systems. Nevertheless, the consequences of droughts are of greater concern to ecologists and administrative bodies. We are convinced that the rationale management of headwater basins is key in ensuring the proper ecological status of whole river networks, which are essential for human society in many aspects.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/w13101352/s1, Table S1: Presumed negative effect of chlorpyrifos and cypermethrin on chosen invertebrate taxonomic groups; Table S2: Abundance and richness of chironomid subfamilies and tribes in individual samples; Table S3: Percentage proportion of functional feeding groups in terms of abundance in individual samples; Figure S1: PCA ordination diagrams showing relationships within the composition of given taxa traits and between these traits and particular explanatory variables; Table S4: Meteorological data of the Vysočina Region (Czech Republic), containing the catchment area of the Doubravka brook, for individual years of the invertebrate sampling; Figure S2: The hydrologic situation during the sampling campaign 2014–2017 in the Doubrava river (Spačice, Czech Republic); Figure S3: Photographic documentation of different hydrologic situations through the year 2016 at the Site C (control = unaffected, Transect 1); Figure S4: Photographic documentation of different hydrologic situations through the year 2016 at the Site 1 (affected, Transect 1); Figure S5. Photographic documentation of different hydrologic situations through the year 2016 at the Site 2 (affected, Transect 2); Table S5: The hydrological status score assigned to each sample during the sampling period 2014–2017; Table S6: Detected taxa and their semi-quantitative abundances in all samples; Table S7: Detected taxa and their ecological and plasticity traits.

Author Contributions

Conceptualization, J.Š., A.K., and M.B.; methodology, J.Š., M.L., and M.F.; formal analysis, J.Š. and M.L.; writing—original draft preparation, M.L. and M.B.; writing—review and editing, A.K., M.B., and J.Š.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Agriculture of the Czech Republic, project no. QK1910282.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the supplementary material.

Acknowledgments

We would like to thank Lukáš Veselý for valuable advices and Miloš Buřič and Jan Kubec for their contribution during field works.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Měkotová, J.; Šarapatka, B.; Štěrba, O.; Harper, D. Restoration of a river landscape: Biotopes as a basis for quantification of species diversity and evaluation of landscape quality. Ecohydrol. Hydrobiol. 2006, 6, 43–51. [Google Scholar] [CrossRef]
  2. Bedient, P.B.; Huber, W.C.; Vieux, B.E. Hydrology and Floodplain Analysis, 4th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2008; pp. 721–722. [Google Scholar]
  3. Poole, G.C. Fluvial landscape ecology: Addressing uniqueness within the river discontinuum. Freshw. Biol. 2002, 47, 641–660. [Google Scholar] [CrossRef] [Green Version]
  4. Ward, J. Riverine landscapes: Biodiversity patterns, disturbance regimes, and aquatic conservation. Biol. Conserv. 1998, 83, 269–278. [Google Scholar] [CrossRef]
  5. Schäfer, R.B.; von der Ohe, P.C.; Rasmussen, J.; Kefford, B.J.; Beketov, M.A.; Schulz, R.; Liess, M. Thresholds for the effects of pesticides on invertebrate communities and leaf breakdown in stream ecosystems. Environ. Sci. Tech. 2012, 46, 5134–5142. [Google Scholar] [CrossRef] [Green Version]
  6. Worthington, T.; Shaw, P.; Daffern, J.; Langford, T. The effects of a thermal discharge on the macroinvertebrate community of a large British river: Implications for climate change. Hydrobiologia 2015, 753, 81–95. [Google Scholar] [CrossRef] [Green Version]
  7. Kennedy, T.L.; Turner, T.F. River channelization reduces nutrient flow and macroinvertebrate diversity at the aquatic terrestrial transition zone. Ecosphere 2011, 2, 1–13. [Google Scholar] [CrossRef]
  8. Johnson, R.C.; JIN, H.S.; Carreiro, M.M.; Jack, J.D. Macroinvertebrate community structure, secondary production and trophic-level dynamics in urban streams affected by non-point-source pollution. Freshw. Biol. 2013, 58, 843–857. [Google Scholar] [CrossRef]
  9. Martin, S.; Bertaux, A.; Le Ber, F.; Maillard, E.; Imfeld, G. Seasonal changes of macroinvertebrate communities in a stormwater wetland collecting pesticide runoff from a vineyard catchment (Alsace, France). AECT 2012, 62, 29–41. [Google Scholar] [CrossRef]
  10. Maund, S.J.; Hamer, M.J.; Lane, M.C.; Farrelly, E.; Rapley, J.H.; Goggin, U.M.; Gentle, W.E. Partitioning, bioavailability, and toxicity of the pyrethroid insecticide cypermethrin in sediments. Environ. Toxicol. Chem. 2002, 21, 9–15. [Google Scholar] [CrossRef]
  11. Bláha, M.; Grabicova, K.; Shaliutina, O.; Kubec, J.; Randák, T.; Zlabek, V.; Buřič, M.; Veselý, L. Foraging behaviour of top predators mediated by pollution of psychoactive pharmaceuticals and effects on ecosystem stability. Sci. Total Environ. 2019, 662, 655–661. [Google Scholar] [CrossRef]
  12. Veselý, L.; Boukal, D.S.; Buřič, M.; Kozák, P.; Kouba, A.; Sentis, A. Effects of prey density, temperature and predator diversity on nonconsumptive predator-driven mortality in a freshwater food web. Sci. Rep. 2017, 7, 1–9. [Google Scholar] [CrossRef]
  13. Racey, P.R.; Swift, S.M.; Rydell, J.; Brodie, L. Bats and insects over two Scottish rivers with contrasting nitrate status. Anim. Conserv. 1998, 1, 195–202. [Google Scholar] [CrossRef]
  14. Benítez-Mora, A.; Camargo, J.A. Ecological responses of aquatic macrophytes and benthic macroinvertebrates to dams in the Henares River Basin (Central Spain). Hydrobiologia 2014, 728, 167–178. [Google Scholar] [CrossRef]
  15. Anderson, C.B.; Rosemond, A.D. Ecosystem engineering by invasive exotic beavers reduces in-stream diversity and enhances ecosystem function in Cape Horn, Chile. Oecologia 2007, 154, 141–153. [Google Scholar] [CrossRef] [PubMed]
  16. Iñiguez-Armijos, C.; Hampel, H.; Breuer, L. Land-use effects on structural and functional composition of benthic and leaf-associated macroinvertebrates in four Andean streams. Aquat. Ecol. 2018, 52, 77–92. [Google Scholar] [CrossRef] [Green Version]
  17. Hawkins, C.P.; Mykrä, H.; Oksanen, J.; Vander Laan, J.J. Environmental disturbance can increase beta diversity of stream macroinvertebrate assemblages. Glob. Ecol. Biogeogr. 2015, 24, 483–494. [Google Scholar] [CrossRef]
  18. Chester, E.; Robson, B. Anthropogenic refuges for freshwater biodiversity: Their ecological characteristics and management. Biol. Conserv. 2013, 166, 64–75. [Google Scholar] [CrossRef]
  19. Oertli, B. Freshwater biodiversity conservation: The role of artificial ponds in the 21st century. Aquat. Conserv. 2018, 28, 264–269. [Google Scholar] [CrossRef]
  20. Harabiš, F.; Dolný, A. Odonates need natural disturbances: How human-induced dynamics affect the diversity of dragonfly assemblages. Freshw. Sci. 2015, 34, 1050–1057. [Google Scholar] [CrossRef]
  21. Schuldt, B.; Buras, A.; Arend, M.; Vitasse, Y.; Beierkuhnlein, C.; Damm, A.; Gharun, M.; Grams, T.E.; Hauck, M.; Hajek, P. A first assessment of the impact of the extreme 2018 summer drought on Central European forests. BAAE 2020, 45, 86–103. [Google Scholar] [CrossRef]
  22. Hänsel, S.; Ustrnul, Z.; Łupikasza, E.; Skalak, P. Assessing seasonal drought variations and trends over Central Europe. Adv. Water Resour. 2019, 127, 53–75. [Google Scholar] [CrossRef]
  23. Beniston, M. The 2003 heat wave in Europe: A shape of things to come? An analysis based on Swiss climatological data and model simulations. Geophys. Res. Lett. 2004, 31, 1–4. [Google Scholar] [CrossRef] [Green Version]
  24. Karger, D.N.; Schmatz, D.; Dettling, G.; Zimmermann, N.E. High resolution monthly precipitation and temperature timeseries for the period 2006–2100. Sci. Data 2020, 7, 1–10. [Google Scholar] [CrossRef] [PubMed]
  25. Štěrba, O.; Měkotová, J.; Benář, V.; Šarapatka, B.; Rychnovská, M.; Kubíček, F.; Řehořek, V. River Landscape and Its Ecosystems, 1st ed.; In Czech Republic; Univerzita Palackého: Olomouc, Czech Republic, 2008; pp. 295–335. [Google Scholar]
  26. Müller, C.; Cramer, W.; Hare, W.L.; Lotze-Campen, H. Climate change risks for African agriculture. Proc. Natl. Acad. Sci. USA 2011, 108, 4313–4315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bhattachan, A.; Emanuel, R.E.; Ardon, M.; Bernhardt, E.S.; Anderson, S.M.; Stillwagon, M.G.; Ury, E.A.; Bendor, T.K.; Wright, J.P. Evaluating the effects of land-use change and future climate change on vulnerability of coastal landscapes to saltwater intrusion. Elementa Sci. Anthrop. 2018, 6, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Rulík, M.; White, S. The Role of Water in the Landscape. In Assessment and Protection of Water Resources in the Czech Republic; Zelenakova, M., Fialová, J., Negm, A.M., Eds.; Springer Water: Cham, Switzerland, 2020; pp. 71–90. [Google Scholar]
  29. Benateau, S.; Gaudard, A.; Stamm, C.; Altermatt, F. Climate Change and Freshwater Ecosystems: Impacts on Water Quality and Ecological Status; Hydro-CH2018Project; Federal Office for the Environment (FOEN): Bern, Switzerland, 2019; pp. 1–110. [Google Scholar] [CrossRef]
  30. Piggott, J.J.; Townsend, C.R.; Matthaei, C.D. Climate warming and agricultural stressors interact to determine stream macroinvertebrate community dynamics. Glob. Chang. Biol. 2015, 21, 1887–1906. [Google Scholar] [CrossRef] [PubMed]
  31. Munz, N.A.; Burdon, F.J.; De Zwart, D.; Junghans, M.; Melo, L.; Reyes, M.; Schönenberger, U.; Singer, H.P.; Spycher, B.; Hollender, J. Pesticides drive risk of micropollutants in wastewater-impacted streams during low flow conditions. Water Res. 2017, 110, 366–377. [Google Scholar] [CrossRef]
  32. Matsumura, F. Toxicology of Insecticides; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; pp. 355–402. [Google Scholar]
  33. Carvalho, F.P. Pesticides, environment, and food safety. Food Energy Secur. 2017, 6, 48–60. [Google Scholar] [CrossRef]
  34. López, Ó.; Fernández-Bolaños, J.G.; Gil, M.V. Classical Insecticides: Past, Present, and Future. In Green Trends in Insect Control; López, Ó., Fernández-Bolaños, J.G., Eds.; Royal Society of Chemistry: Cambridge, UK, 2011; pp. 53–93. [Google Scholar]
  35. Nowell, L.H. Pesticides in Stream Sediment and Aquatic Biota: Distribution, Trends, and Governing Factors; CRC Press: Boca Raton, FL, USA, 1999; pp. 47–463. [Google Scholar]
  36. Beketov, M.; Foit, K.; Schäfer, R.; Schriever, C.; Sacchi, A.; Capri, E.; Biggs, J.; Wells, C.; Liess, M. SPEAR indicates pesticide effects in streams–comparative use of species-and family-level biomonitoring data. Environ. Pollut. 2009, 157, 1841–1848. [Google Scholar] [CrossRef]
  37. Muenze, R.; Orlinskiy, P.; Gunold, R.; Paschke, A.; Kaske, O.; Beketov, M.A.; Hundt, M.; Bauer, C.; Schüürmann, G.; Moeder, M. Pesticide impact on aquatic invertebrates identified with Chemcatcher® passive samplers and the SPEARpesticides index. Sci. Total Environ. 2015, 537, 69–80. [Google Scholar] [CrossRef]
  38. Belgers, J.D.M.; Van Lieverloo, R.J.; Van der Pas, L.J.; Van den Brink, P.J. Effects of the herbicide 2, 4-D on the growth of nine aquatic macrophytes. Aquat. Bot. 2007, 86, 260–268. [Google Scholar] [CrossRef]
  39. Stehle, S.; Knäbel, A.; Schulz, R. Probabilistic risk assessment of insecticide concentrations in agricultural surface waters: A critical appraisal. Environ. Monit. Assess. 2013, 185, 6295–6310. [Google Scholar] [CrossRef] [PubMed]
  40. Liess, M.; Ohe, P.C.V.D. Analyzing effects of pesticides on invertebrate communities in streams. Environ. Toxicol. Chem. 2005, 24, 954–965. [Google Scholar] [CrossRef] [PubMed]
  41. Pereira, J.L.; Antunes, S.C.; Castro, B.B.; Marques, C.R.; Gonçalves, A.M.; Gonçalves, F.; Pereira, R. Toxicity evaluation of three pesticides on non-target aquatic and soil organisms: Commercial formulation versus active ingredient. Ecotoxicology 2009, 18, 455–463. [Google Scholar] [CrossRef]
  42. Pestana, J.L.; Loureiro, S.; Baird, D.J.; Soares, A.M. Fear and loathing in the benthos: Responses of aquatic insect larvae to the pesticide imidacloprid in the presence of chemical signals of predation risk. Aquat. Toxicol. 2009, 93, 138–149. [Google Scholar] [CrossRef]
  43. Ankley, G.T.; Call, D.J.; Cox, J.S.; Kahl, M.D.; Hoke, R.A.; Kosian, P.A. Organic carbon partitioning as a basis for predicting the toxicity of chlorpyrifos in sediments. Environ. Toxicol. Chem. 1994, 13, 621–626. [Google Scholar] [CrossRef]
  44. Blann, K.L.; Anderson, J.L.; Sands, G.R.; Vondracek, B. Effects of agricultural drainage on aquatic ecosystems: A review. Crit. Rev. Environ. Sci. Tech. 2009, 39, 909–1001. [Google Scholar] [CrossRef]
  45. Brown, C.D.; Van Beinum, W. Pesticide transport via sub-surface drains in Europe. Environ. Pollut. 2009, 157, 3314–3324. [Google Scholar] [CrossRef]
  46. Matthaei, C.D.; Piggott, J.J.; Townsend, C.R. Multiple stressors in agricultural streams: Interactions among sediment addition, nutrient enrichment and water abstraction. J. Appl. Ecol. 2010, 47, 639–649. [Google Scholar] [CrossRef]
  47. Townsend, C.R.; Scarsbrook, M.R.; Dolédec, S. The intermediate disturbance hypothesis, refugia, and biodiversity in streams. Limnol. Oceanogr. 1997, 42, 938–949. [Google Scholar] [CrossRef] [Green Version]
  48. Ebert, S.; Hulea, O.; Strobel, D. Floodplain restoration along the lower Danube: A climate change adaptation case study. Clim. Dev. 2009, 1, 212–219. [Google Scholar] [CrossRef]
  49. Knillmann, S.; Orlinskiy, P.; Kaske, O.; Foit, K.; Liess, M. Indication of pesticide effects and recolonization in streams. Sci. Total Environ. 2018, 630, 1619–1627. [Google Scholar] [CrossRef] [PubMed]
  50. von der Ohe, P.C.; Goedkoop, W. Distinguishing the effects of habitat degradation and pesticide stress on benthic invertebrates using stressor-specific metrics. Sci. Total Environ. 2013, 444, 480–490. [Google Scholar] [CrossRef] [PubMed]
  51. Arenas-Sánchez, A.; Rico, A.; Vighi, M. Effects of water scarcity and chemical pollution in aquatic ecosystems: State of the art. Sci. Total Environ. 2016, 572, 390–403. [Google Scholar] [CrossRef]
  52. Stehle, S.; Schulz, R. Agricultural insecticides threaten surface waters at the global scale. Proc. Natl. Acad. Sci. USA 2015, 112, 5750–5755. [Google Scholar] [CrossRef] [Green Version]
  53. Boulton, A.J.; Lake, P.S. Effects of Drought on Stream Insects and its Ecological Consequences. In Aquatic Insects: Challenges to Populations: Proceedings of the Royal Entomological Society’s 24th Symposium; Lancaster, J., Briers, R.A., Eds.; CABI: Wallingford, UK, 2008; pp. 81–102. [Google Scholar]
  54. Boulton, A.J.; Peterson, C.G.; Grimm, N.B.; Fisher, S.G. Stability of an aquatic macroinvertebrate community in a multiyear hydrologic disturbance regime. Ecology 1992, 73, 2192–2207. [Google Scholar] [CrossRef]
  55. Pastorino, P.; Zaccaroni, A.; Doretto, A.; Falasco, E.; Silvi, M.; Dondo, A.; Elia, A.C.; Prearo, M.; Bona, F. Functional Feeding Groups of Aquatic Insects Influence Trace Element Accumulation: Findings for Filterers, Scrapers and Predators from the Po Basin. Biology 2020, 9, 288. [Google Scholar] [CrossRef]
  56. Geist, J. Strategies for the conservation of endangered freshwater pearl mussels (Margaritifera margaritifera L.): A synthesis of conservation genetics and ecology. Hydrobiologia 2010, 644, 69–88. [Google Scholar] [CrossRef] [Green Version]
  57. Schletterer, M.; Füreder, L. The family Prosopistomatidae (Ephemeroptera): A review on its ecology and distribution, with particular emphasis on the European species Prosopistoma pennigerum Müller, 1785. Aquat. Insects 2009, 31, 603–620. [Google Scholar] [CrossRef]
  58. Diamond, J. Overview: Laboratory experiments, field experiments, and natural experiments. Nature 1983, 304, 586–587. [Google Scholar] [CrossRef]
  59. Ligeiro, R.; Hughes, R.M.; Kaufmann, P.R.; Heino, J.; Melo, A.S.; Callisto, M. Choice of field and laboratory methods affects the detection of anthropogenic disturbances using stream macroinvertebrate assemblages. Ecol. Indic. 2020, 115, 106382. [Google Scholar] [CrossRef]
  60. Ferenčík, M.; Schovánková, J. Occurrence of Polar Contaminants (Pesticides and Pharmaceuticals) in Suspended Matter in Povodí Labe Basin (in Czech). In Proceedings of the VII. Conference “Sedimenty vodných tokov a nádrží“ (in Slovak), Bratislava, Slovak Republic, 22–23 May 2013; Hucko, P., Tölgyessy, P., Eds.; Slovenská vodohospodárska spoločnosť člen ZSVTS pri Výskumnom ústave vodného hospodárstva: Bratislava, Slovak Republic, 2013; pp. 125–132. [Google Scholar]
  61. Rocha, A.A.; Monteiro, S.H.; Andrade, G.C.; Vilca, F.Z.; Tornisielo, V.L. Monitoring of pesticide residues in surface and subsurface waters, sediments, and fish in center-pivot irrigation areas. J. Braz. Chem. Soc. 2015, 26, 2269–2278. [Google Scholar] [CrossRef]
  62. Kokeš, J.; Němejcová, D. Metodika Odběru a Zpracování vzorků Makrozoobentosu Tekoucích vod Metodou PERLA; VÚV TGM Praha: Prague, Czech Republic, 2006; pp. 1–10. (In Czech) [Google Scholar]
  63. Kindt, R. Package ‘BiodiversityR’. Package for Community Ecology and Suitability Analysis. Version 2.11–12. 2019. Available online: https://CRAN.R-project.org/package=BiodiversityR (accessed on 11 March 2021).
  64. Paisley, M.; Trigg, D.; Walley, W. Revision of the biological monitoring working party (BMWP) score system: Derivation of present-only and abundance-related scores from field data. River Res. Appl. 2014, 30, 887–904. [Google Scholar] [CrossRef] [Green Version]
  65. Lepš, J.; Šmilauer, P. Multivariate Analysis of Ecological Data Using CANOCO; Cambridge University Press: Cambridge, UK, 2003; pp. 60–75. [Google Scholar]
  66. Friberg, N.; Lindstrøm, M.; Kronvang, B.; Larsen, S.E. Macroinvertebrate/sediment relationships along a pesticide gradient in Danish streams. Hydrobiologia 2003, 493, 103–110. [Google Scholar] [CrossRef]
  67. Cuffney, T.F.; Wallace, J.B.; Webster, J.R. Pesticide manipulation of a headwater stream: Invertebrate responses and their significance for ecosystem processes. Freshw. Invertebr. Biol. 1984, 3, 153–171. [Google Scholar] [CrossRef]
  68. Fritz, K.M.; Dodds, W.K. Resistance and resilience of macroinvertebrate assemblages to drying and flood in a tallgrass prairie stream system. Hydrobiologia 2004, 527, 99–112. [Google Scholar] [CrossRef]
  69. Stanley, E.H.; Buschman, D.L.; Boulton, A.J.; Grimm, N.B.; Fisher, S. Invertebrate resistance and resilience to intermittency in a desert stream. Amer. Midl. Nat. 1994, 288–300. [Google Scholar] [CrossRef] [Green Version]
  70. Clarke, A.; Mac Nally, R.; Bond, N.; Lake, P.S. Macroinvertebrate diversity in headwater streams: A review. Freshw. Biol. 2008, 53, 1707–1721. [Google Scholar] [CrossRef]
  71. Buřič, M.; University of South Bohemia, České Budějovice, Czech Republic; Kouba, A.; University of South Bohemia, České Budějovice, Czech Republic. Personal communication, 2018.
  72. Reynolds, J.; Souty-Grosset, C.; Richardson, A. Ecological roles of crayfish in freshwater and terrestrial habitats. Freshwater Crayfish 2013, 19, 197–218. [Google Scholar] [CrossRef]
  73. Watts, M. Ecotoxicity. In Chlorpyrifos as a Possible Global POP; Pesticide Action Network North America: Oakland, CA, USA, 2012; pp. 1–34. Available online: http://www.panna.org/sites/default/files/Chlorpyrifos_as_POP_final_2013_0.pdf (accessed on 12 May 2021).
  74. European Commission. Review report for the active substance alpha-cypermethrin. In Health and Consumer Protection Directorate-General; European Commission: Brussels, Belgium, 2004; pp. 1–70. Available online: https://epdfs.co.uk/doc/304998/european-commission-alpha-cypermethrin-sanco/4335/2000- (accessed on 12 May 2021).
  75. Rubach, M.N.; Ashauer, R.; Maund, S.J.; Baird, D.J.; Van den Brink, P.J. Toxicokinetic variation in 15 freshwater arthropod species exposed to the insecticide chlorpyrifos. Environ. Toxicol. Chem. 2010, 29, 2225–2234. [Google Scholar] [CrossRef]
  76. Plahuta, M.; Tišler, T.; Toman, M.J.; Pintar, A. Toxic and endocrine disrupting effects of wastewater treatment plant influents and effluents on a freshwater isopod Asellus aquaticus (Isopoda, Crustacea). Chemosphere 2017, 174, 342–353. [Google Scholar] [CrossRef]
  77. Spänhoff, B.; Bischof, R.; Böhme, A.; Lorenz, S.; Neumeister, K.; Nöthlich, A.; Küsel, K. Assessing the impact of effluents from a modern wastewater treatment plant on breakdown of coarse particulate organic matter and benthic macroinvertebrates in a lowland river. Water Air Soil Pollut. 2007, 180, 119–129. [Google Scholar] [CrossRef]
  78. Iversen, T.M.; Wiberg-Larsen, P.; Hansen, S.B.; Hansen, F.S. The effect of partial and total drought on the macroinvertebrate communities of three small Danish streams. Hydrobiologia 1978, 60, 235–242. [Google Scholar] [CrossRef]
  79. Extence, C. The effect of drought on benthic invertebrate communities in a lowland river. Hydrobiologia 1981, 83, 217–224. [Google Scholar] [CrossRef]
  80. Finelli, C.M.; Hart, D.D.; Merz, R. Stream insects as passive suspension feeders: Effects of velocity and food concentration on feeding performance. Oecologia 2002, 131, 145–153. [Google Scholar] [CrossRef]
  81. Currie, D.C.; Craig, D.A. Feeding strategies of larval black flies. In Black Flies: Ecology, Population Management, and Annotated World List; Kim, K.C., Merritt, R.W., Eds.; Pennsylvania State University Press: University Park, PA, USA, 1987; pp. 155–170. [Google Scholar]
  82. Horak, C.N.; Assef, Y.A.; Grech, M.G.; Miserendino, M.L. Agricultural practices alter function and structure of macroinvertebrate communities in Patagonian piedmont streams. Hydrobiologia 2020, 847, 3659–3676. [Google Scholar] [CrossRef]
  83. Lenat, D.R.; Crawford, J.K. Effects of land use on water quality and aquatic biota of three North Carolina Piedmont streams. Hydrobiologia 1994, 294, 185–199. [Google Scholar] [CrossRef]
  84. Aspin, T.W.; Hart, K.; Khamis, K.; Milner, A.M.; O’Callaghan, M.J.; Trimmer, M.; Wang, Z.; Williams, G.M.; Woodward, G.; Ledger, M.E. Drought intensification alters the composition, body size, and trophic structure of invertebrate assemblages in a stream mesocosm experiment. Freshw. Biol. 2019, 64, 750–760. [Google Scholar] [CrossRef]
  85. Vallenduuk, H.J.; Pillot, H.K.M. Chironomidae Larvae, Vol. 1: Tanypodinae: General Ecology and Tanypodinae; KNNV Publishing: Zeist, The Netherlands, 2007; pp. 74–132. [Google Scholar]
  86. Pillot, H.K.M. Chironomidae Larvae, Vol. 2: Chironomini: Biology and Ecology of the Chironomini; KNNV Publishing: Zeist, The Netherlands, 2009; pp. 6–247. [Google Scholar]
  87. Pillot, H.K.M. Chironomidae Larvae, Vol. 3: Orthocladiinae: Biology and Ecology of the Aquatic Orthocladiinae; KNNV Publishing: Zeist, The Netherlands, 2014; Volume 3, pp. 7–287. [Google Scholar]
  88. Tokeshi, M. Life cycles and population dynamics. In The Chironomidae: The Biology and Ecology of Non-Biting Midges; Armitage, P.D., Cranston, P.S., Pinder, L.C.V., Eds.; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1995; pp. 225–268. [Google Scholar]
  89. Barron, M.G.; Woodburn, K.B. Ecotoxicology of chlorpyrifos. Rev. Environ. Contam. Toxicol. 1995, 1–93. [Google Scholar] [CrossRef]
  90. Mehler, W.T.; Du, J.; Lydy, M.J.; You, J. Joint toxicity of a pyrethroid insecticide, cypermethrin, and a heavy metal, lead, to the benthic invertebrate Chironomus dilutus. Environ. Toxicol. Chem. 2011, 30, 2838–2845. [Google Scholar] [CrossRef]
  91. Lancaster, J.; Ledger, M.E. Population-level responses of stream macroinvertebrates to drying can be density-independent or density-dependent. Freshw. Biol. 2015, 60, 2559–2570. [Google Scholar] [CrossRef]
  92. Marcogliese, D. The impact of climate change on the parasites and infectious diseases of aquatic animals. Rev. Sci. Tech. 2008, 27, 467–484. [Google Scholar] [CrossRef]
  93. Hutchens, J.J., Jr.; Chung, K.; Wallace, J.B. Temporal variability of stream macroinvertebrate abundance and biomass following pesticide disturbance. J. N. Am. Benthol. Soc. 1998, 17, 518–534. [Google Scholar] [CrossRef]
  94. Johnson, R.C.; Carreiro, M.M.; Jin, H.-S.; Jack, J.D. Within-year temporal variation and life-cycle seasonality affect stream macroinvertebrate community structure and biotic metrics. Ecol. Indic. 2012, 13, 206–214. [Google Scholar] [CrossRef]
  95. Graf, W.; Murphy, J.; Dahl, J.; Zamora-Munoz, C.; Lopez-Rodriguez, M.J. Distribution and Ecological Preferences of European Freshwater Organisms. Volume 1. Trichoptera; Pensoft Publishers: Sofia, Bulgaria, 2008; p. 284. [Google Scholar]
  96. Vannote, R.L.; Minshall, G.W.; Cummins, K.W.; Sedell, J.R.; Cushing, C.E. The river continuum concept. Can. J. Fish. Aquat. Sci. 1980, 37, 130–137. [Google Scholar] [CrossRef]
  97. Wallace, J.B.; Eggert, S.L.; Meyer, J.L.; Webster, J.R. Multiple trophic levels of a forest stream linked to terrestrial litter inputs. Science 1997, 277, 102–104. [Google Scholar] [CrossRef] [Green Version]
  98. Herbst, D.B.; Cooper, S.D.; Medhurst, R.B.; Wiseman, S.W.; Hunsaker, C.T. Drought ecohydrology alters the structure and function of benthic invertebrate communities in mountain streams. Freshw. Biol. 2019, 64, 886–902. [Google Scholar] [CrossRef]
  99. Zilli, F.L.; Montalto, L.; Marchese, M.R. Benthic invertebrate assemblages and functional feeding groups in the Paraná River floodplain (Argentina). Limnologica 2008, 38, 159–171. [Google Scholar] [CrossRef] [Green Version]
  100. Adámek, Z.; Helešic, J.; Maršálek, B.; Rulík, M. Applied Hydrobiology; University South Bohemia in České Budějovice, Faculty of Fisheries and Protection of Waters: Vodňany, Czech Republic, 2014; pp. 205–244. [Google Scholar]
  101. Wallace, J.B.; Vogel, D.S.; Cuffney, T. Recovery of a headwater stream from an insecticide-induced community disturbance. J. N. Am. Benthol. Soc. 1986, 5, 115–126. [Google Scholar] [CrossRef]
  102. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 Establishing a Framework for Community Action in the Field of Water Policy. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=celex%3A32000L0060 (accessed on 23 February 2021).
  103. Orlinskiy, P.; Münze, R.; Beketov, M.; Gunold, R.; Paschke, A.; Knillmann, S.; Liess, M. Forested headwaters mitigate pesticide effects on macroinvertebrate communities in streams: Mechanisms and quantification. Sci. Total Environ. 2015, 524, 115–123. [Google Scholar] [CrossRef]
  104. Michael, J.L. Pesticides used in forestry and their impacts on water quality. South. Weed Sci. Soc. 2000, 53, 81–91. [Google Scholar]
  105. Schulz, R.; Liess, M. A field study of the effects of agriculturally derived insecticide input on stream macroinvertebrate dynamics. Aquat. Toxicol. 1999, 46, 155–176. [Google Scholar] [CrossRef]
  106. Hunt, L.; Bonetto, C.; Marrochi, N.; Scalise, A.; Fanelli, S.; Liess, M.; Lydy, M.J.; Chiu, M.-C.; Resh, V.H. Species at Risk (SPEAR) index indicates effects of insecticides on stream invertebrate communities in soy production regions of the Argentine Pampas. Sci. Total Environ. 2017, 580, 699–709. [Google Scholar] [CrossRef] [Green Version]
  107. Lake, P.S. Disturbance, patchiness, and diversity in streams. J. N. Am. Benthol. Soc. 2000, 19, 573–592. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of the study area: (A) Locations of the sampling sites along the longitudinal profile of the Doubravka brook (in blue; flow direction indicated by the blue arrow). The contaminated drainage channel is highlighted in red. The sampling sites are shown as yellow circles: C and Site 1 = Transect 1; Site 2 = Transect 2. (B) Study area shown as the red rectangle in the Czech Republic.
Figure 1. Map of the study area: (A) Locations of the sampling sites along the longitudinal profile of the Doubravka brook (in blue; flow direction indicated by the blue arrow). The contaminated drainage channel is highlighted in red. The sampling sites are shown as yellow circles: C and Site 1 = Transect 1; Site 2 = Transect 2. (B) Study area shown as the red rectangle in the Czech Republic.
Water 13 01352 g001
Figure 2. RDA ordination diagrams showing the relationships between disturbances (CPS-CP = chlorpyrifos-cypermethrin contamination derived from the chlorpyrifos sediment contamination, and Drought = drought status; see Materials and Methods) as solid black arrows, with taxa represented by grey arrowed lines. The power of the both disturbances to predict taxonomic compositions was significant (adjusted p < 0.05). The length of arrows corresponds to the quality of fit. (A) Relationships between disturbances and the composition of deeper taxonomic levels extrapolated from the abundances of all detected taxa (CPS-CP: explained variation 26.9%, pseudo-F = 3.3; Drought: explained variation 22.9%, pseudo-F = 3.7, respectively; adjusted variation explained by both predictors = 37.28%). Additional information for the labels marked by the asterisks: Platyhelminthes is represented by a single species—Polycelis nigra; Crustacea is represented by Ostracoda (not determined into a deeper taxonomic level) and Asellus aquaticus; Bivalvia is represented by only two species—Pisidium casertanum and P. personatum. (B) Relationships between disturbances and the composition of deeper taxonomic levels extrapolated from the abundances of arthropods. (CPS-CP: explained variation 24.7%, pseudo-F = 3.0; Drought: explained variation 27.0%, pseudo-F = 4.5, respectively; adjusted variation explained by both predictors = 39.65%). The asterisk marking Ostracoda indicates that this group was not determined into a deeper taxonomic level.
Figure 2. RDA ordination diagrams showing the relationships between disturbances (CPS-CP = chlorpyrifos-cypermethrin contamination derived from the chlorpyrifos sediment contamination, and Drought = drought status; see Materials and Methods) as solid black arrows, with taxa represented by grey arrowed lines. The power of the both disturbances to predict taxonomic compositions was significant (adjusted p < 0.05). The length of arrows corresponds to the quality of fit. (A) Relationships between disturbances and the composition of deeper taxonomic levels extrapolated from the abundances of all detected taxa (CPS-CP: explained variation 26.9%, pseudo-F = 3.3; Drought: explained variation 22.9%, pseudo-F = 3.7, respectively; adjusted variation explained by both predictors = 37.28%). Additional information for the labels marked by the asterisks: Platyhelminthes is represented by a single species—Polycelis nigra; Crustacea is represented by Ostracoda (not determined into a deeper taxonomic level) and Asellus aquaticus; Bivalvia is represented by only two species—Pisidium casertanum and P. personatum. (B) Relationships between disturbances and the composition of deeper taxonomic levels extrapolated from the abundances of arthropods. (CPS-CP: explained variation 24.7%, pseudo-F = 3.0; Drought: explained variation 27.0%, pseudo-F = 4.5, respectively; adjusted variation explained by both predictors = 39.65%). The asterisk marking Ostracoda indicates that this group was not determined into a deeper taxonomic level.
Water 13 01352 g002
Figure 3. Ordination diagram showing the relationships between the average functional feeding groups (FFGs) represented by solid grey arrows and disturbances (CPS-CP = chlorpyrifos-cypermethrin contamination estimated from the chlorpyrifos sediment contamination and Drought = drought status; see Materials and Methods) represented by solid black arrows. Both disturbances have significant (p < 0.05) power to predict the FFG composition. Arrow length corresponds to the quality of fit. The FFG was extrapolated from abundances of all detected organisms (CPS-CP: explained variation 25.3%, pseudo-F = 3.1; Drought: explained variation 22.3%, pseudo-F = 3.4, respectively; adjusted variation explained by both predictors = 34.50%).
Figure 3. Ordination diagram showing the relationships between the average functional feeding groups (FFGs) represented by solid grey arrows and disturbances (CPS-CP = chlorpyrifos-cypermethrin contamination estimated from the chlorpyrifos sediment contamination and Drought = drought status; see Materials and Methods) represented by solid black arrows. Both disturbances have significant (p < 0.05) power to predict the FFG composition. Arrow length corresponds to the quality of fit. The FFG was extrapolated from abundances of all detected organisms (CPS-CP: explained variation 25.3%, pseudo-F = 3.1; Drought: explained variation 22.3%, pseudo-F = 3.4, respectively; adjusted variation explained by both predictors = 34.50%).
Water 13 01352 g003
Table 1. Chlorpyrifos content in the sediments analysed in the samples from Sites C (control = unaffected, Transect 1), 1 (affected, Transect 1) and 2 (affected, Transect 2) during the sampling period 2014–2017. The asterisks indicate macroinvertebrate sampling. LOQ = 20 µg kg−1 (d.w.).
Table 1. Chlorpyrifos content in the sediments analysed in the samples from Sites C (control = unaffected, Transect 1), 1 (affected, Transect 1) and 2 (affected, Transect 2) during the sampling period 2014–2017. The asterisks indicate macroinvertebrate sampling. LOQ = 20 µg kg−1 (d.w.).
Sampling DateSite C
(µg kg−1 d.w.)
Site 1
(µg kg−1 d.w.)
Site 2
(µg kg−1 d.w.)
Time Passed after Accident (Days)
4th IV 2014<LOQ13,0003505
20th VI 2014 *<LOQ1285183
17th IV 2015 *<LOQ87<LOQ384
26th VII 2016 *<LOQ33<LOQ850
19th IX 2017 *<LOQ25<LOQ1476
Table 2. Hydrological status score.
Table 2. Hydrological status score.
Description:Score:
Fast-flowing1
Loss of fast-flowing habitats2
Loss of lateral connectivity to stream-edge habitats3
Loss of longitudinal connectivity, flow stops, isolated pools form4
Pools shrink water quality deteriorates5
Total absence of water6
Table 3. Summary of macroinvertebrate community parameters calculated for each sample within the sampling period (2014–2017). Abbreviations: Shannon–Wiener index (H’), number of species (N), Abundance (Abu). For meaning of parameters, see Section 2.4.3.
Table 3. Summary of macroinvertebrate community parameters calculated for each sample within the sampling period (2014–2017). Abbreviations: Shannon–Wiener index (H’), number of species (N), Abundance (Abu). For meaning of parameters, see Section 2.4.3.
2014201520162017
ParameterC12C12C12C12
Abundance11422282752644966447228422313722733087761432
Richness26915392623342221432326
H‘2.401.441.513.102.061.692.052.441.532.992.172.77
SPEARpesticides index0.230.150.020.110.130.090.200.060.050.1600.08
SPEARrefuge index0.250.290.050.140.150.020.06000.120.080.10
SPEARpesticides index (old)0.480.440.070.250.290.110.260.060.050.270.080.18
N SPEARpesticides 621543931903
N SPEARrefuge731651200523
N SPEARpesticides (old)1352119411311426
Abu SPEARpesticides3068140302694968136024
Abu SPEARrefuge1424465836240002801672
Abu SPEARpesticides (old)4485279866281349684161696
N EPT1452121047011124
Abu EPT44852714894281000683961680
EPT%42.7322.810.2522.989.730.633.5200.9411.972.065.59
Original BMWP score1004226107774210249381416572
ASPT index (±SEM)6.67 ± 0.776.00 ± 1.353.71 ± 0.786.29 ± 0.785.50 ± 0.794.20 ± 0.635.10 ± 0.624.08 ± 0.533.80 ± 0.795.88 ± 0.634.64 ± 0.704.80 ± 0.76
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Let, M.; Špaček, J.; Ferenčík, M.; Kouba, A.; Bláha, M. Insecticides and Drought as a Fatal Combination for a Stream Macroinvertebrate Assemblage in a Catchment Area Exploited by Large-Scale Agriculture. Water 2021, 13, 1352. https://doi.org/10.3390/w13101352

AMA Style

Let M, Špaček J, Ferenčík M, Kouba A, Bláha M. Insecticides and Drought as a Fatal Combination for a Stream Macroinvertebrate Assemblage in a Catchment Area Exploited by Large-Scale Agriculture. Water. 2021; 13(10):1352. https://doi.org/10.3390/w13101352

Chicago/Turabian Style

Let, Marek, Jan Špaček, Martin Ferenčík, Antonín Kouba, and Martin Bláha. 2021. "Insecticides and Drought as a Fatal Combination for a Stream Macroinvertebrate Assemblage in a Catchment Area Exploited by Large-Scale Agriculture" Water 13, no. 10: 1352. https://doi.org/10.3390/w13101352

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop