A Multicompartment Approach - Diatoms, Macrophytes, Benthic Macroinvertebrates and Fish - To Assess the Impact of Toxic Industrial Releases on a Small French River

The River Luzou flows through a sandy substrate in the South West of France. According to the results of two assessment surveys, the Water Agency appraised that this river may not achieve the good ecological status by 2015 as required by the Water Framework Directive (2000/60/EC). This ecosystem is impacted by industrial effluents (organic matter, metals and aromatic compounds). In order to assess and characterize the impact, this study aimed to combine a set of taxonomic and non-taxonomic metrics for diatoms, macrophytes, macroinvertebrates and fish along the up- to downstream gradient of the river. Diversity metrics, biological indices, biological and ecological traits were determined for the four biological quality elements (BQE). Various quantitative metrics (biomass estimates) were also calculated for diatom communities. The results were compared to physicochemical analysis. Biological measurements were more informative than physicochemical analysis, in the context of the study. Biological responses indicated both the contamination of water and its intensity. Diversity metrics and biological indices strongly decreased with pollution for all BQE but diatoms. Convergent trait selection with pollution was observed among BQE: reproduction, colonization strategies, or trophic regime were clearly modified at impaired sites. Taxon size and relation to the substrate diverged among biological compartments. Multiple anthropogenic pollution calls for alternate assessment methods of rivers' health. Our study exemplifies the fact that, in the case of complex contaminations, biological indicators can be more informative for environmental risk, than a wide screening of contaminants by chemical analysis alone. The combination of diverse biological compartments provided a refined diagnostic about the nature (general mode of action) and intensity of the contamination.


Introduction
Rivers are a key component for the development of civilizations and, therefore, human-induced impacts on these ecosystems are diverse and always increasing. Since 2000 in Europe, the Water Framework Directive [1] requires water bodies to achieve good ecological status by 2015, and consider biology as the central element of the assessment. This Directive recommends using four biological quality elements (BQEs) -diatoms, macrophytes, benthic macroinvertebrates and fish -to assess the rivers' ecological status.
The River Luzou flows through a sandy substrate in the South West of France (Landes). From the results of two assessment surveys (2004 and 2006), the Water Agency appraised that this river may not achieve the good ecological status by 2015 [2]. The ecosystem is impacted by the effluent of an industrial plant producing rubber: the pollution is very diverse (mainly organic matter, metals and aromatic compounds including aniline), variable in composition over the year, and released in pulses (51 releases in 2010).
Nevertheless, even when combined, indices often lead to excessive reduction of environmental information [17], so the use of non-taxonomic measures such as biological and ecological traits provides new perspectives for biological assessment methods [18][19]. Traits represent qualitative and quantitative information related to the biology of organisms and their relationship with the environment. Traits describe taxon ecological preferenda, life cycle, morphology, physiology or behaviour [12,[20][21][22][23][24][25].
In this context, we hypothesized that the assessment of the ecological status of the River Luzou according to the recommendations of the WFD, e.g. biological indices calculation and punctual physicochemical analysis of the water, would be insufficient to highlight the impact suffered by the studied ecosystem. In other words, we assumed that the combination of relevant taxonomic and non-taxonomic metrics for the four BQEs cited along the up-to downstream gradient of the River Luzou, is a consistent way to assess the impacts of the industrial releases on the ecosystem. We tested our hypothesis by sampling benthic macroinvertebrates, diatoms, macrophytes, fish and water at three sites on an up-to downstream gradient of the river, in different seasons, in 2009 and 2010. We compared the relevance of the different types of results, e.g. classical biological indices, taxonomic and non-taxonomic metrics, punctual physicochemical parameter values, to accurately monitor the changes in major ecosystem components.

Study area and sampling sites
The River Luzou (Landes, South West France) is a small river (28 km long) flowing through a sandy substrate, and characterized by low pH, low conductivity and low nutrient concentrations (Fig. 1). The industrial plant, situated about 18 km from the river source, introduces highly toxic pollution and in addition increases the river water temperature through use of water for cooling and tank washing processes. The wastewater is stored in a settling pool and overflow is released into the River Luzou directly downstream of the plant. Three sampling sites publicly accessible were selected along the study reach: a site called ''Up'' situated upstream of the factory (Lambert 93 coordinates: X385686, Y6312556) was considered as the reference station (i.e. unimpacted by the toxic effluents); sites ''Down 1'' (X387414, Y6312572) and ''Down 2'' (X389923, Y6310416) were situated about 500 m and 4500 m downstream of the factory, respectively. Benthic macroinvertebrates, diatoms, macrophytes, fish and water were sampled at the three sites once or several times, in different seasons, in 2009 and 2010 (see below for details). The water, macrophyte, macroinvertebrate, and phytobenthos sampling procedures were approved by the Adour-Garonne Water Agency.
The substances monitored were chosen because they are known to be present in the releases, as stated in the self-monitoring data from the plant itself. Parameters showing values below the detection limit for all upstream and downstream sites were removed from the analysis. For those occasionally below, half the value of the detection limit was used [26].
Additionally, maximal theoretical concentration of aniline (a target component of the toxic wastes released by the industrial site) was evaluated using industrial monitoring data from year 2010 and the QMNA5 value (monthly low flow value that may not occur more than once every 5 years).

Fish sampling
Ethics statement: All fish were properly collected and handled in an ethical manner, with all required permissions from the Adour-Garonne Water Agency. No other permissions were required for completion of this research, and this study does not include endangered or otherwise protected species.
Fish surveys were carried out by experienced fisheries staff of the departmental fishing federation AAPPMA64 (accredited by ECCEL Environnement) and all sampling procedures complied with the French and European Union legislation on animal welfare. Electric fishing was carried out at the minimum power settings needed to incapacitate the fish and thus no adverse impact was expected. The fish were handled with great care. This includes electrofishing and manipulations (counting, weighing and measurements) where fish were maintained in river water. After measurements, all fish were returned alive into the river.
All sites were electrofished (pulsed direct current waveform), during low flow period (September 2009) and according to the NF EN 14011 standard [27]. Electrofishing was authorized by the Adour Garonne Water Agency and performed by operators accredited by ECCEL Environnement. Fish were sorted and stored in a large basin, and then counted, measured and weighed, then released alive in the water. In the case of very numerous individuals, fish were counted and weighed in homogenous sets. The sampling area was systematically reported.

Macrophyte sampling
All sampling was carried out according the standard NF T90-395 [5] during vegetation periods (September 2009, July 2010). The areas covered by macrophyte beds and by each taxon were evaluated for each site. Taxa difficult to identify in situ were collected, packed and transported to the laboratory for determination.

Benthic macroinvertebrate sampling
Macroinvertebrates were collected at each of the 4 seasons (Autumn: September 2009, Winter: November 2009, Spring: April 2010 and Summer: July 2010). Macroinvertebrates were collected with a Surber sampler (mesh size 500 mm, sampling area 0.05 m 2 ), the device required by the French standardized protocol XP T90-333 [28]. Micro-samples were taken at each site from among twelve mesohabitats defined as visually distinct units within the stream, considering apparent physical uniformity (sensu Armitage et al. 1995), and described by a combination of substrate types and current velocities. Mesohabitats were sampled in a hierarchical order according to the IBGN standard [6] to maximize the taxonomic richness of the faunal assemblage at the site scale, after gathering sample units. Macroinvertebrates were sorted and identified to the family level following standard XP T90-388 [29], except for some groups identified at a higher taxonomic level (i.e. Oligochaeta, Bryozoa, Nematoda, Hydracarina). The twelve micro-samples were pooled to finally constitute a single sample, used to calculate the IBGN index [6]. Quantitative measurements. Artificial substrates were immersed in the river to quantify dry weight, ash-free dry mass of biofilms, chlorophyll a concentrations, and the number of live and dead diatom cells. Six glass slides (total surface area reaching 300 cm 2 ) were placed in a rack at the three sites [30]. After 15 days immersion (t15) and 30 days immersion (t30), three slides were removed from the water and scraped into a standard volume of mineral water, to obtain three replicates per sampling date then separated into aliquots.

Phytobenthos sampling
A 20-mL aliquot was used to determine the dry weight (DW) and ash-free dry mass (AFDM) of the biofilm, expressed as mg cm 22 , according to European standard NF EN 872 [31]. Ten millilitres of the suspension were filtered through a Whatman GF/ C filter, then extracted with acetone for 24 h before spectrophotometric analyses. Chlorophyll a concentrations were calculated following Lorenzen [32]. A 5-mL aliquot was preserved with 0.5 mL of formalin solution for diatom cell density enumeration [live and dead, 33] and taxonomic identification.
Diatom community characterization. Diatom samples were collected from natural surfaces (pebbles or macrophytes, at t30) and artificial substrates (at t15 and t30), according to a standardized method NF T 90-354 [7]. For each slide, 400 valves were determined to the lowest taxonomical level possible. Diatom species were identified at 10006 magnification (Leitz DMRB light microscope), mainly according to Krammer and Lange-Bertalot [34] and Lange-Bertalot [35], by examining permanent slides of cleaned diatom frustules, digested in boiling H 2 O 2 (30%) and HCl (35%) and mounted in a high refractive index medium (Naphrax, Northern Biological Supplies Ltd., UK; RI = 1.74). A total of 234 diatom taxa were identified. The 107 taxa with abundances higher than five individuals (considering all samples) were used to describe community structure.

Biofilm-related quantitative metrics
Seasonal variations in biomass: dry weight, ash-free dry mass, chlorophyll a concentrations, and the number of live and dead diatoms were reduced by normalization using the mean values calculated on the reference site (Up) at each season.

Data analysis
Major differences between sampling sites considering physicochemical parameters were investigated using Principal Component Analysis (PCA). Data were normalized and redundant variables identified through Spearman pairwise correlation tests. Among redundant variables, only one was kept. Concentrations of aniline were taken as supplementary variable, assuming that the upstream site (Up) was free of aniline.
Concerning biological data, diatom and macroinvertebrate community structures were described using PCA based on species relative abundances. Functional traits-related information was obtained from the tables Taxa 6 Samples and Taxa 6 Traits for each compartment. Taxa 6 Samples tables were expressed in abundances, except for macrophytes (percentage cover) and fish (biomass). In order to produce functional profiles for each trait (relative distribution of the information among the categories), the following process was applied: i) for a given site and for each category of traits, taxon scores were weighted by abundances (cover or biomass); ii) the sums of the weighted scores were then expressed as a relative abundance distribution (within a trait), giving the site trait profile. For diatoms, this analysis was based on the 107 dominant taxa.
Only diatoms and benthic macroinvertebrates enabled statistical analyses (four sampling dates). Taxonomic, indicial and functional metrics were compared between sites by Kruskal-Wallis tests. For macrophytes and fish (respectively two and one replicates), only simple visual comparisons were possible.
All analyses were performed using the R software, version 3.0.2 [42], with packages ade4 [43] for descriptive analyses of data and agricolae [44] for Kruskal-Wallis tests.

Analysis of environmental data
Alkalinity and HCO 3 2 giving redundant ecological information (Spearman test = 1), HCO 3 2 was removed from the dataset. The following parameters, even though they were known to be released in the effluents, were systematically below the detection limit and were thus removed from the analysis: copper, lead, manganese, nickel, cadmium, mercury, cyanide, arsenic, benzene and nitrobenzene, toluene and 2-nitrotoluene. For the same reason, pesticides were also removed.
The environmental parameters finally kept for analysis and their contributions to axes 1 and 2 are listed in Table 1.
Axes 1 and 2 account for 70% of the total inertia (Fig. 2a). Axis 1 discriminates upstream from downstream conditions but does not clearly discriminate Down 1 from Down 2, while both Axes 1 and 2 separates winter samples from the others. Downstream sites are characterized by higher concentrations of Zn, Nkjeldahl, SO 4 22 , NH 4 + , BOD5, higher conductivity, alkalinity and temperature (Fig. 2b).
Mean calculated concentration of aniline downstream was 0.44 mg L 21 . Aniline, represented a posteriori, shows a cos 2 equal to 0.42 on axis 1 and to 0.28 on axis 2.

Analysis of biological data
Upstream communities were typical of the Landes ecoregion. Minnow (Phoxinus phoxinus), gudgeon (Gobio gobio), stone loach (Barbatula barbatula), brook lamprey (Lampetra planeri) and eel (Anguilla anguilla) were present in the River Luzou. The benthic macroinvertebrate fauna (thirty taxa) is dominated by Chironomidae and Gammaridae. Acidophilous and neutrophilous diatom taxa from the genera Eunotia and Brachysira were abundant, associated with numerous Karayevia oblongella, Tabellaria flocculosa and Peronia fibula. The macrophytic community was mainly composed of phanerogams characteristic of oligotrophic and acidic waters: Myriophyllum alterniflorum and Potamogeton polygonifolium. Table 2 gathers information about differences between upstream and downstream conditions, according to Kruskal-Wallis tests. Taxonomic metrics based on fish, macrophytes and benthic macroinvertebrates showed similar responses to toxic pollution: a strong decrease in species richness, diversity and equitability was observed at station Down 1. Fish response was very marked, as only two minnow specimens were found at station Down 1 (considered as a null biomass), and a still low biomass was harvested at station Down 2. Macrophyte community structure was strongly shifted from a diversity of algae, bryophytes (e.g. Fontinalis antipyretica) and spermaphytes (in particular, Callitriche platycarpa), to a Sphaerotilus sp.-dominated community at Down 1, causing the dramatic decrease in taxonomic metrics observed. At Down 2, the diversity of taxonomic groups increased, and high percentages of algae (such as Cladophora sp.) were found. Taxonomic metrics for benthic macroinvertebrates were also significantly lower at stations Down 1 and Down 2, reflecting that macroinvertebrate species composition was mainly influenced by the up-to downstream gradient (Fig. 3a), more than by season (cold vs. warm waters). No significant differences between downstream and upstream sites were globally observed for diatoms, but the flora showed large seasonal variations (Fig. 3b).

Taxonomical metrics and indices
Poor IPR, IBMR and IBGN scores classified station Down 1 in a ''bad'' ecological status, and station Down 2 in a ''medium'' to ''poor'' status ( Table 2). The high abundance (up to 60%) at downstream sites of Achnanthidium minutissimum, a species considered as oligo-to mesotrophic by IBD and IPS, maintained good index scores whereas the SPEARorganic index differentiated upstream and downstream sites.
Diatom communities sampled on glass slides 30 days after immersion presented a strong decrease in live cell density and in chlorophyll a concentrations downstream. At station Down 1, dry weight and ash-free dry mass increased for communities sampled on glass slides 15 days after immersion, but diatom mortality was significantly higher (decrease in the ratio of live to dead cells).

Functional metrics
Diatoms. Even though Kruskall Wallis test results were not significant, some functional traits from communities sampled on natural substrates (pebbles) showed trends between upstream and downstream sites (Table 3). Pioneer species were more abundant at station Down 2 while ''high profile'' and ''motile'' species Table 1. Average (6standard errors, n = 4) values and contributions of environmental parameters to axis 1 and axis 2 of the PCA performed on physicochemical parameters (see Figure 2).  (Table 4). Taxa with tegumentary respiration, burrowers, interstitial or permanently attached were also favoured downstream. Taxa without any form of resistance decreased in favour of taxa able to produce cocoons. Active dispersal (aerial and aquatic) decreased while aquatic passive dispersion increased (station Down 1). Feeding behaviours were also modified downstream whereas absorbers and deposit feeders, eating detritus, fine sediments and microorganisms dominated. As a result shredders, filter-feeders and scrapers (only at station Down 1) eating plants (alive or dead) or living microinvertebrates decreased. Concerning ecological traits, polysaprobic and eutrophic taxa increased downstream while xenosaprobic, oligotrophic and mesotrophic taxa decreased.
Macrophytes. The different floristic group proportions were clearly modified downstream (Table 5). Station Down 1 presented more than 99% of heterotrophic forms, hence phanerogramrelated traits were not reliable for its description. At station Down 2 the phanerogams did not recover the upstream reference status, allowing the installation of filamentous algae. At this site, taxa were preferentially annual, with asexual reproduction mainly by fragmentation and sexual reproduction based on higher numbers of seeds, with entomophilous dispersion. These taxa showed larger emergent leaves, related to a higher morphological index (large size).
Fish. Functional metrics were not applicable in station Down 1, where only one species (2 individuals) was found. At station Down 2 the biomass and abundance of invertivorous species decreased in favour of omnivores (Table 6). Pelagic species became dominant over benthic ones.

Discussion
The importance of biology for the assessment of ecological status PCA performed on environmental data discriminated upstream from downstream conditions, but sites Down 1 and Down 2 were not well distinguished, in contrast to results obtained with biological data (especially for macro-invertebrates). A high seasonal variation can also be noted, with different dilution conditions according to whether the flow was low or high.
Moreover, several toxicants known to be released in the plant effluents were not detected by physicochemical analysis. Thus,     Table 3. Diatom-related traits: calculation (mean and standard deviation) and differences between sites.   such analysis does not seem to be particularly reliable in the case of the River Luzou: to characterize intermittent pollution, measurements should be performed exactly during the releases, or by high resolution analyses. The problem of toxicity related to cocktail effects and degradation compounds can also hardly be tackled in this way. In this context, biological related metrics are potentially more informative than chemical analyses. Seasonal variations in community structure (Fig. 3) were also observed, and the different metrics used allowed the nature, and intensity, of the pollution present in the River Luzou to be highlighted.

Temporal scales of biological responses
The seasonal variations observed in the PCA performed on environmental data were correlated to changes in diatom responses over time. Although quantitative measurements indicated decreasing diatom biomass from up-to downstream whatever the sampling season, the number of cells settled was up to 15 fold higher in warmer conditions compared to winter (data not shown). The seasonality in water contamination was highlighted by, e.g., higher percentages of teratologies (up to 3% on natural substrates) and of the species Achnanthidium minutissimum downstream in Autumn, indicating toxic pollution [45][46].
Diatoms have fast growth rates (from hours to days), and thus respond very quickly to variations in their environment. BQE with longer life span reflect more averaged water quality on different time scales. Seasonal patterns were thus less pronounced for macrophytes and macroinvertebrates that integrated global quality over the year. Ultimately, fish responses were expected to reveal environmental conditions on the longer term (years).

Taxonomic metrics reveal pollution intensity
Richness, diversity and equitability indices are classically used to evaluate ecological status of water and their decline is indicative of a disturbed environment [47]. Except for diatoms, these metrics clearly decreased for all the biological compartments studied, at station Down 1, which is in accordance with the literature [16,[48][49]. They recovered to variable extents at Down 2, indicating weaker biological impact likely due to lower toxicant availability (by dilution and/or toxicant degradation). For diatoms, the phenomenon observed contrasted with the literature, rather reporting a decrease of these taxonomic metrics in toxic conditions [45,50]. For this BQE in our study, richness, diversity and equitability indices were not relevant to highlight the toxic impact of pollution (see below).

Biological indices rather indicate the nature of pollution
Biological indices except IBD and IPS drastically decreased downstream, due to the presence of numerous tolerant taxa for IBMR and IBGN (Sphaerotilus sp., Oligochaeta and Chironomidae), or to the lack of fish populations for IPR. SPEARorganic index being negatively correlated to toxicants (icides, surfactants, petrochemicals) [11], its decrease downstream was also consistent. In contrast IPS and IBD scores remained good due to the high proportion of Achnanthidium minutissimum (up to 60%). However, the occurrence of abnormal forms, mainly affecting A. minutissimum, clearly characterized downstream conditions. Above 1%, this rate of occurrence is considered to reflect the impact of toxic pollutants on diatom communities [45][46]. Moreover, recent works also suggested that, due to its pioneering character, A. minutissimum was indicative of toxic pollution [45]. Therefore, the use of diatom-based biological indices to highlight toxic pollution is not recommended, but a careful analysis of community composition (as used for index calculations) can also provide information regarding the nature of the contamination.
Additionally, quantitative metrics (decrease of chlorophyll a concentrations and diatom density downstream) were consistent with Morin et al. [45], and with Wang et al. [51] who demonstrated that a derivative of aniline could inhibit adhesion in certain diatom species. The percentage of live cells tended to decrease downstream, as already reported by Stevenson and Bahls [52] and Gillet et al. [53] who examined whether the percentage of live diatoms in periphyton communities could be used as a metric of human disturbance in streams and rivers.
Periphytic biomass (dry weight and ash-free dry mass) tended to increase downstream, whereas diatom cell numbers and chlorophyll a did not. This unexpected result can perhaps be related to the periphyton becoming more heterotrophic, which would be consistent with the massive development of Sphaerotilus sp. observed at Down 1 in the macrophytic community. One could also hypothesize that rivers from Landes ecoregion, characterized by naturally acidic waters and nutrient depletion, represent a particular ecosystem where artificial eutrophication (concomitant to industrial release) can enhance the periphyton growth and richness. Other recent works [54] reported reduced impacts on biomass of metal toxicity in acid-adapted biofilms. These phenomena may explain the typical quantitative response of phytobenthos towards toxic pollution.  Survival strategies under high pollution Functional metrics illustrated how the ecosystem of the River Luzou changed under toxic pollution. According to the results, many modifications concerning the four BQE were convergent, like reproduction and colonization strategies, or the trophic regime. Southwood [55] wrote that those physiological adaptations were typically found in impaired sites, as they induce tolerance to harsh conditions. First, when escaping pollution, as fish attempt to do, was not possible, resistant taxa became dominant (Sphaerotilus sp., Oligochaeta and Chironomidae), with production of cocoons as extreme resistance forms for macroinvertebrates.
Early colonizers with ruderal strategies [56] were also favoured, like the diatom Achnanthidium minutissimum which reached an abundance of up to 60% downstream. Morin et al. [57] already observed this high abundance of A. minutissimum under toxic conditions, implying a thinning of the biofilm by loss of highprofile taxa (according to Passy [22], high profile guild reaches a maximum in nutrient-rich sites and in conditions of low flow disturbance). Achnanthidium minutissimum was probably less disfavoured thanks to its small size reducing the exposure time to toxicants [58], or its adnate posture. Macrophyte communities evolved towards the predominance of annual taxa producing a great number of reproductive organs per year and per individual (e.g. Berula erecta, Scirpus fluitans or Sparganium erectum), ensuring a rapid spread across the river. Macroinvertebrates' high colonization ability was provided by the combination of a greater number of reproductive cycles and dispersion by drift, those two strategies being already reported in the literature [59]. Drift is an important way of species dissemination and recolonisation of river systems by lotic macroinvertebrates [60], and is known to increase with chemical disturbance [61].
Finally, taxon sizes showed opposite trends. Diatom taxa with lower biovolumes were favoured, whereas bigger macroinvertebrate taxa, or macrophytes with large emergent leaves, increased downstream. This, however, resulted from the same strategy of organisms facing toxicants. Indeed, a decrease of diatom cell sizes are generally observed in environments exposed to toxic pollutions: reduction of size expresses increased vegetative multiplication, and reduced sexual reproduction (not measured here, but concordant with traits of other compartments: higher reproduction rates and asexual reproduction), and selection of smaller species (reducing uptake of toxicants into the cell).

Toxic pollution also drives indirect changes across the trophic web
Our results highlighted changes in feeding habits, which represent an important aspect of the community trophic structure  Table 6. Fish-related traits: calculation (mean and standard deviation) and differences between sites. modifications. Summarizing, the reference ecosystem (Up) showed balanced communities composed of species typical of the Landes ecoregion [62]. The different compartments are driven by complex interactions. Primary producers provide food resources to primary consumers (like macroinvertebrates and fish feeding on phytobenthos) but also refuge and habitat and/or egg-laying substrates (especially macrophytes). The modifications observed in these communities (disappearance or sharp decline of phanerogams, biofilm thinning) represent a real impact on invertebrates and may explain the changes observed downstream. In particular, in such streams like the river Luzou where there is a poor diversification of abiotic substrates, macrophytes not only provide a food source, but also a shelter for invertebrates [63]. Macroinvertebrate structure and biomass are also modified by the presence or absence of predators. With high pollution (Down 1), biomass of primary producers decreased to be replaced by filamentous heterotrophs, affecting the subsequent components in a cascade. Aside from some potential direct toxicity (according to the toxicants mode of action), macroinvertebrates were probably driven by the resource, selecting absorbers and deposit feeders feeding on detritus and microorganisms in accordance with Archaimbault et al. [60] and Schultheis et al. [64]. Moreover, the decline of macroinvertebrate communities may have forced invertivorous macroinvertebrates and fish to leave the site. Further downstream (Down 2), primary producers tended to diversify, as well as primary consumers and predators. This partial recovery thus reflected both direct (reduction of toxic pressure, slight increase in nutrient availability) and indirect (return to a more balanced ecosystem) improvement of the ecosystem.

Conclusions
In conclusion, our study revealed that in a context of multiple contaminants, of pulse inputs, of complex cocktails, and/or of release of unknown substances, a combination of different biological measurements, from different aquatic communities, can be much more informative than punctual physicochemical analysis and single biotic indices.
Under such strong and diverse anthropogenic pressure, a multicompartment approach allows the integrated observation of community trajectories towards adaptation, accounting for the complex biotic relationships in aquatic ecosystems. Colonization strategies, reproduction and trophic regimes seem to be key indicators of this adaptation.
A further step would be to continue this multi-compartment survey during restoration programmes, to identify the behaviour, and/or biological elements that tend to recover more rapidly. Data from the literature suggest that mobile organisms would have the greatest ability to recolonize sites after water quality improvement. The time necessary to reach complete recovery of the ecosystem when pollution ceases is also likely to be very variable among biological compartments and is an important component to be determined in the context of the implementation of the Water Framework Directive.