Next Article in Journal
The Status of Arsenic Pollution in the Greek and Cyprus Environment: An Overview
Next Article in Special Issue
Assessing the Ecological Status of European Rivers and Lakes Using Benthic Invertebrate Communities: A Practical Catalogue of Metrics and Methods
Previous Article in Journal
Planning the Urban Waterfront Transformation, from Infrastructures to Public Space Design in a Sea-Level Rise Scenario: The European Union Prize for Contemporary Architecture Case
Previous Article in Special Issue
Selection of Macroinvertebrate Indices and Metrics for Assessing Sediment Quality in the St. Lawrence River (QC, Canada)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

How to Assess the Ecological Status of Highly Humic Lakes? Development of a New Method Based on Benthic Invertebrates

1
Laboratory of Hydrobiology, Institute of Biology, University of Latvia, Jelgavas 1, LV-1004 Riga, Latvia
2
Statistics Unit, Riga Stradiņš University, Baložu 14, LV-1048 Riga, Latvia
3
Institute of Food Safety, Animal Health and Environment BIOR, Lejupes 3, LV-1076 Riga, Latvia
4
European Commission, Joint Research Centre, 21027 Ispra, Italy
*
Authors to whom correspondence should be addressed.
Water 2021, 13(2), 223; https://doi.org/10.3390/w13020223
Submission received: 7 December 2020 / Revised: 13 January 2021 / Accepted: 14 January 2021 / Published: 18 January 2021
(This article belongs to the Special Issue Freshwater Macroinvertebrates: Main Gaps and Future Trends)

Abstract

:
Highly humic lakes are typical for the boreal zone. These unique ecosystems are characterised as relatively undisturbed habitats with brown water, high acidity, low nutrient content and lack of macrophytes. Current lake assessment methods are not appropriate for ecological assessment of highly humic lakes because of their unique properties and differing human pressures acting on these ecosystems. This study proposes a new approach suitable for the ecological status assessment of highly humic lakes impacted by hydrological modifications. Altogether, 52 macroinvertebrate samples from 15 raised bog lakes were used to develop the method. The studied lakes are located in the raised bogs at the central and eastern parts of Latvia. Altered water level was found as the main threat to the humic lake habitats since no other pressures were established. A multimetric index based on macroinvertebrate abundance, littoral and profundal preferences, Coleoptera taxa richness and the Biological Monitoring Working Party (BMWP) Score is suggested as the most suitable tool to assess the ecological quality of the highly humic lakes.

1. Introduction

Humic lakes, also known as brown-water lakes, are typical for the boreal zone, located between 50° to 70° N latitude. These ecosystems are characterized by dark water colour, low water transparency, and low pH caused by the high concentration of dissolved organic matter (DOM), mostly originating from the catchment area and consisting of refractory humic substances [1,2]. Over the last decades, an impressive body of evidence has been accumulated which suggests that DOM is a major modulator of the structure and function of lake ecosystems, affecting numerous features such as light regime, thermal stratification, nutrient availability, primary production, and microbial metabolism [3,4,5,6]. Numerous studies have shown that the biological communities differ considerably from those of clear-water lakes (phytoplankton: [7]; macrophytes: [8]; periphyton: [9]; zooplankton: [10]; fish fauna: [11]). Moreover, their response to human stressors might differ too [12,13,14], asking for monitoring and assessment approaches different to those used for clear-water lakes [15,16].
Multiple human pressures, such as nutrient enrichment, hydrological and morphological alterations, invasion of non-native species and climate change, affect humic lakes [17,18]. Some of these pressures are similar to those impacting clear water lakes, but some are different, e.g., bog lakes are affected by artificial peatland drainage and peat extraction associated with habitat degradation, erosion, increased leaching of nutrients and dissolved organic carbon [19,20]. Similarly, peatland lakes are impacted by forestry practices (afforestation, fertilization, and clear-cutting) which have been shown to increase catchment loadings of nutrients, sediments and dissolved organic carbon (DOC) to these ecosystems [21,22,23]. These impacts can have a profound effect on the water quality, lake habitats, associated biological assemblages and conservation value of these lakes [22,24,25].
In Europe, the Water Framework Directive (WFD) [26] establishes a framework for the protection of inland and coastal waters. According to the WFD, lakes have to be classified into five status classes (high, good, moderate, poor and bad) based on four biological quality elements (BQEs)—phytoplankton, benthic flora, benthic invertebrates, and fish fauna. In addition, physico-chemical elements (e.g., nutrient conditions) and hydromorphological elements (e.g., flow conditions) are used to support ecological classification. The European Union (EU) member states have to identify degraded water bodies (i.e., less than good status) and to establish programmes of measures for each river basin district to reduce significant anthropogenic pressures and achieve good water status.
A large number of lake assessment methods have been intercalibrated and included in the European countries’ monitoring toolkits [27]. In the recent years, all of these methods have been intercalibrated (i.e., compared and harmonized) among the EU member states [28,29]. Lake assessment methods include both primary producers, e.g., phytoplankton (e.g., [30,31]), macrophytes (e.g., [32,33]), phytobenthos (e.g., [34]), and heterotrophs such as benthic invertebrates and fish fauna (e.g., [33,35,36]).
However, two problems are still overlooked. At first, the majority of assessment systems target nutrient enrichment, while other key pressures are largely neglected. This is especially true regarding hydromorphological pressures, which affect a considerable number of lakes across Europe [37]. Only few assessment systems tackle the ecological effects of these pressures (shore degradation: [38]; water level fluctuations: [39]) and only two of these systems have been intercalibrated among member states [36].
Second, despite the well-known differences among clear and brown-water lakes [1], the current lake assessment systems are adopted mostly for clear-water lakes. Recently, several studies have raised the issue that assessment systems might not be appropriate for humic lake assessment [13,17,40]. Hence, there is an imperative need for the development of appropriate assessment methods targeting humic lakes.
Studies on macroinvertebrates in highly humic lakes are mainly focused either on biodiversity [41,42] or specific taxonomic groups, e.g., chironomids [43] or Coleoptera [44]. Raised bog water bodies are also known as habitats for rare and protected macroinvertebrate species [41].
Mires and bogs cover 4.9% of the territory of Latvia [45] ranking Latvia number 9 by the total area of peatlands among all the European countries [46]. Bog lakes are listed as protected habitats within the EU Habitats Directive emphasizing their high conservation value [47]. However, these lakes have been impacted by a range of anthropogenic activities, most importantly anthropogenic drainage and peat harvesting which can lead to water level fluctuations, loss of biodiversity and degradation of the lake ecological status [48,49]. Nevertheless, the effects of these impacts are poorly understood and there are no assessment tools in place to assess the ecological condition of humic lakes. According to the River Basin Management Plans (RBMPs) of Latvia, current methods reflect bog lakes at poor or bad status, though the anthropogenic pressures are irrelevant [50]. As a result, there is a need to develop new methods for the ecological assessment of the highly humic lakes.
The objectives of this study are (i) to explore littoral benthic invertebrate community response to hydrological alterations in highly humic lakes; (ii) to develop a biotic index for assessing hydrological alterations in these aquatic ecosystems.

2. Materials and Methods

2.1. Study Sites

Altogether 15 highly humic lakes were studied at seven raised bogs comprising the national monitoring data from the Latvian Environmental, Geology and Meteorology Centre and studies from the Institute of Biology, University of Latvia. The lakes were divided into two groups: lakes with altered water level due to drainage and lakes with natural or restored hydrological regime (see Figure 1).
All of the studied lakes are located in the territories included in the Natura 2000 network. In the central part of Latvia, Cena Mire and Melnais Lake Mire are represented by the samples from five small bog lakes each. In the Cena Mire, lakes are located in pristine areas, while in the Melnais Lake Mire—close to the peat excavation fields, thus representing a hydrologically disturbed state. In the northern part of Latvia, sampling was conducted in Lake Ramatas Lielezers and Lake Soku, both being with altered water level due to outflowing drainage ditches. The water level was receded at these lakes with visible open peat outcrops at the shoreline (see Figure 2).
Eastern part of Latvia is represented by three lakes (Pieslaista, Deguma and Orlovas) that have an unaltered or restored hydrological regime. The water level at these lakes was natural and not affected due to the drainage ditches or peat excavation (see Figure 3).

2.2. Physical and Chemical Parameters

Waterbodies of Melnais Lake Mire and Cena Mire were sampled in May 2015 and analyses were conducted at the Laboratory of Soils of the University of Latvia. Soku, Ramatas Lielezers, Deguma, Pieslaista and Orlovas lakes were sampled four times in different seasons in 2017 and analysis was conducted at the Laboratory of Latvian Environment, Geology and Meteorology Centre. In 2015, pH and electric conductivity (EC) were measured in-situ by using a portable pH tester (HI 98127, HANNA instruments, Sarmeola di Rubano, Italy) and conductivity tester (The Original Dist HI 98300, HANNA instruments, Sarmeola di Rubano, Italy). In 2017, these parameters were measured in-situ by using a portable probe (HQ40d, Hach Companies, Loveland, CO, USA). Total phosphorus (TP) was detected by ascorbic acid method after the digestion using potassium persulfate. Total nitrogen (TN) samples were digested by potassium persulfate, then nitrates were reduced to nitrites in a Cd column and analysed spectrophotometrically. Water colour was analysed spectrophotometrically using the Pt/Co scale [51].

2.3. Macroinvertebrate Sampling and Sample Processing

In larger lakes, sampling was conducted in a 100 m long, representative shoreline section, while in smaller lakes sampling was conducted around all of the shoreline. Benthic macroinvertebrate samples were taken by hand net (frame size 0.25 × 0.25 m, mesh of 0.5 mm) using the sweeping technique. A hand net was placed on the bottom parallel to the shore, if the depth was less than 0.8 m or under the vegetation overhang at the same depth and moved upwards over the vegetation stands to the surface. Five replicates of 0.5 m sweeps were chosen in proportion to habitat types, e.g., of stands of Menyanthes trifoliata, Sphagnum cuspidatum, Batrachospermum turfosum, Carex spp., bare littoral, etc. At the small lakes of Melnais Lake Mire and Cena Mire, samples were taken in May 2015. Sampling at Soku, Ramatas Lielezers, Deguma, Pieslaista and Orlovas lakes was conducted in May and October 2017. Additional samples were taken at lakes Soku, Ramatas Lielezers and Deguma in May 2019.
Sampled material was washed through a sieve with a mesh size of 0.5 mm at field. All replicates were placed in the polyethylene containers, labelled and preserved in 70% ethyl alcohol (final concentration).
Preserved samples were washed at the laboratory; all specimens were picked out from the vegetation, detritus and peat particles. Macroinvertebrates were identified to the smallest achievable taxonomical (species, genera) level, excluding Oligochaeta and juvenile Hydrachnidia. Specimens of Diptera were identified to the family level.

2.4. Data Analysis

2.4.1. Selection of Metrics

According to the WFD, the ecological quality assessment indices are required multimetric consisting of composition, abundance, sensitive/tolerant taxa and diversity metrics [26]. Macroinvertebrate metrics were calculated using ASTERICS 4.0.4. software (Wageningen Software Labs, Wageningen, The Netherlands) [52]. Numerically unsuitable metrics and majority of the metrics specific for the lotic habitats were excluded from further analysis. In total, 139 indices describing 52 samples from water bodies in open raised bogs were tested for a selection of multimetric index according to requirements of the EU WFD [26]. We generally followed the procedure described by Hering et al. [53], beginning with the reduction of dimensionality by the evaluation of each metric value distribution between the altered and the natural water bodies.

2.4.2. Sensitivity to Stressor

Only descriptors correlating to the stressor gradient can be used in the development of a multimetric index. We used Mann–Whitney–Wilcoxon U tests (each metric as dependent in groups of stressor, thereinafter U test) together with boxplots and simple binary logistic regression (groups of stressor as dependant: altered = 0, natural = 1, thereinafter binary logistic regression (BLR)) to select only statistically significant metrics for further processing. Metric values prior to BLR were centred by mean and scaled by standard deviation to ease the convergence, while unscaled values were used in the boxplots and U test.
With type I statistical error rate 0.05, the U test returned 66 statistically significant indices and BLR-44 (Table S1). All the indices found significant by BLR were significant also by U test. As a reason for this step is the reduction of dimensionality with selection of potentially most important indices, we did not account for possible false discoveries due to the multiple testing. We included metrics (except life index due to specific relation to lotic environment) found significant by both methods in further investigation of their suitability.

2.4.3. Numerical Suitability

The multimetric index must consist of metrics that tend to describe large gradients with possibly low skew and preferably without outliers [53]. These numeric properties can be the best evaluated graphically if values are from the same or a similar scale. We used the violin (density) plots with a point overlay (with stressor group indicated by a shape) to select one to six metrics per a metric type. Before plotting, each metric was scaled by its observed maximum value. Graphs used in evaluation are provided in Figure 4 with 14 metrics included in the further investigation marked with the dark background.

2.4.4. Ecological Relevance

The group of experts authoring this paper evaluated each previously selected metric to exclude possibly non-explanatory correlations and those of ecologically low meaning. Additionally, metrics created for ecoregions other than boreal were excluded resulting in ten metrics for further use.

2.4.5. Correlations

Metrics with Spearman`s correlation coefficient |≥0.8| are considered redundant and only one of them must be used [53]. Additionally, we tried to avoid possible multicollinearity by selecting metrics with even lower (|≤0.6|) values [54] to obtain at least one metric per one WFD criteria. Correlation coefficients of the selected metrics are available in Table S2.
At this point, we decided to keep following metrics:
  • For composition/abundance: abundance [ind/m2];
  • For richness/diversity: Coleoptera (taxa);
  • For sensitivity/tolerance: BMWP (Biological Monitoring Working Party) Score;
  • For functional metrics: (%) littoral, (%) littoral + profundal, (%) profundal.
As the selected functional metrics are highly correlated both statistically and ecologically, we decided to compare three models with inclusion of every previously mentioned metric from above and functional metrics as follows:
  • (%) littoral;
  • (%) littoral + profundal;
  • (%) profundal.
As in every principal model where metric types have the same weights, we used the mean value as the result for the multimetric index.

2.4.6. Scaling

To ensure that the multimetric index is limited between 0 and 1, and to avoid potential influence of some extremely high- or low-quality sites anchoring is suggested [53]. We used the same approach for metrics decreasing with increasing impairment, but corrected the approach for the metrics increasing with increasing impairment to:
Value = 1—(Metric result—Lower Anchor)/(Upper Anchor—Lower Anchor). We corrected values >1 to 1 and negative values to 0.
To ensure good spread of sites within multimetric index, we compared several anchor values (Table S3):
  • 5th percentile and 95th percentile;
  • 10th percentile and 90th percentile;
  • 10th percentile and 80th percentile;
  • Each of the previous with prespecified values for “Coleoptera (taxa)” as 0 for lower and 4 for upper.

2.4.7. Quality Classification

We used the quality classes in accordance with WFD demands, following suggestion of Hering et al. [53]:
  • Reference ≥0.8;
  • Good ≥0.6 < 0.8;
  • Moderate ≥0.4 < 0.6;
  • Poor ≥0.2 < 0.4;
  • Bad <0.2.
We consider an index to be the best, if natural lakes are concentrated at the reference and good quality classes, while altered lakes are at bad- and poor-quality classes with some mixture present in a class of moderate quality.
We used the software R 4.0.3 (The R Foundation for Statistical Computing, Vienna, Austria) for data analysis [55]. Data processing and visualisations were performed within the tidyverse ecosystem [56].

3. Results

3.1. Characterisation of Chemical and Environmental Variables

The studied lakes are poly-humic lakes as indicated by high water colour values (114–666 mg Pt/L) and low pH values (3.35–6.09). Electric Conductivity (EC) is in the range of 21–65 μS/cm, concentrations of TN are 0.43–1.68 mg/L and TP are 0.017–0.061 mg/L. The highest conductivity, water color and total nitrogen values were observed in waterbodies of the Melnais Lake Mire (see Table 1).

3.2. Benthic Invertebrate Taxa

A list of benthic invertebrate taxa found in highly humic lakes is presented in the Supplemental Material (Table S4). Altogether, 18,808 individuals from 106 macroinvertebrate taxa are recorded at the studied bog lakes, of which the orders Coleoptera, Odonata and Trichoptera have the highest species richness. Mayflies (Ephemeroptera) are represented by four species, of which Leptobphlebia vespertina is the most widespread, missing only in the water bodies of the Cena Mire. In the humic lakes of Melnais Lake Mire and Cena Mire, larvae of dragonflies Leucorrhinia albifrons and L. pectoralis are recorded. These two species are included in the Bern Convention [57] and in the Habitats Directive [47]. The highest abundance of macroinvertebrates is recorded from waterbodies of Melnais Lake Mire and Cena Mire, ranging from 813 to 2014 individuals per sample, while the benthic invertebrate abundance at larger lakes ranges from 51 to 1202 individuals.
The macroinvertebrate orders Diptera and Coleoptera are the most abundant taxa at the lakes of the Cena Mire, while chironomids dominate in lakes of the Melnais Lake Mire. Taxonomic structure at lakes Orlovas, Deguma, Pieslaista, Ramatas Lielezers and Soku vary due to repeated sampling in different years and seasons, generally with Diptera, Coleoptera and Ephemeroptera as the dominant taxa. Molluscs (Gastropoda and Bivalvia) are completely absent at the studied lakes (see Figure 5).

3.3. Metrics

The comparison of macroinvertebrate abundance (ind/m2) between samples of bog lakes with altered water level and those with natural water level show significant differences (p < 0.001) (Figure 6). The macroinvertebrate abundance in lakes with natural level varies from 109 to 1202 individuals, while in lakes with altered level from 51 to 502. Similar significant differences are found between the number of taxa varying from 9 to 22 taxa in natural lakes and 3 to 18 taxa in altered lakes (p < 0.001). Additionally, our results show higher values of BMWP Score at natural lakes rather than the altered ones (p = 0.002). The BMWP values in natural lakes vary from 33 to 74, whereas from 10 to 67 in bog lakes with altered water level. Taxa richness of Coleoptera and ETCO (Ephemeroptera, Trichoptera, Coleoptera, Odonata) show the same significant differences between the studied lakes (p < 0.001). Natural bog lakes are represented by 1 to 5 Coleoptera species, while in the altered lakes 0 to 3 species are found. The number of ETCO varies from 5 to 14 species in natural lakes while 1 to 9 species in bog lakes with altered water level, respectively. Nevertheless, altered lakes are represented by taxa preferring littoral and profundal habitats (p < 0.001) and higher percentage of Diptera (p < 0.001). Number of taxa of Trichoptera (p = 0.023), Odonata (p = 0.009) and Heteroptera (p = 0.007) is higher at humic lakes with natural water level.

3.4. Multimetric Index

We found the model 2 (Figure 7) ((%) littoral + profundal and abundance (ind/m2) and Coleoptera (taxa) and BMWP Score) with P10 and P80 anchoring and predefined Coleoptera (taxa) anchoring to be the best classifier for the highly humic lakes.
The distribution of values of the highly humic lake multimetric index are shown in Figure 8.

4. Discussion

The total number of taxa found in the studied highly humic lakes is relatively low and similar to those of other studies [40,41]. Desrochers and van Duinen [58] noted that the chemical constraints of humic lakes almost entirely exclude several macroinvertebrate taxa, e.g., lumbricid worms, isopods, and snails and the low nutrient availability may limit the presence of species with high nutrient preferences. In addition, acidity is known as a limiting factor for many macroinvertebrate groups, e.g., freshwater snails [59], mussels [60] and mayflies [61]. In general, macroinvertebrate species specialized on different mire habitats are obviously able to realize their life cycle at water pH of 4.0–5.0 [62]. The macroinvertebarte taxonomical composition and lack of molluscs indicate the peculiarities of the studied highly humic lakes.
The invertebrate assemblage in altered lakes showed a reduced taxonomic richness, especially in Trichoptera, Coleoptera, Heteroptera and Odonata, an increased proportion of Diptera and an overall lower numerical abundance.
The best-performing metrics were total abundance (for composition/abundance category), number of Coleoptera taxa (for richness/diversity category), BMWP score (for sensitivity/tolerance category), and % littoral and profundal (for functional metric category).

4.1. Use of Benthic Invertebrates in Lake Ecological Assessment

In lakes, phytoplankton and macrophytes are the most widely used communities for ecological assessment. However, in the last decades numerous systems using benthic invertebrates have been developed, following the requirements of the WFD [36]. These systems differ by habitat sampled (mostly littoral, but some sampled the profundal), putative pressure assessed and metrics used.
According to the concept of multimetric approach [53,63] different types of metrics should be included into the assessment system (composition/abundance metrics, richness/diversity metrics; sensitivity/tolerance metrics; functional metrics).

4.1.1. Sensitivity/Tolerance Metrics

Almost all lake benthic invertebrate assessment systems include some kind of sensitivity index, most common are Average Score per Taxon (ASPT) index [64], Acid Water Indicator Community (AWIC) index [65], Benthic Quality Index [66] and Fauna index [67].
Originally, ASPT and BMWP indices have been developed for river water quality in Britain [68], however they have proved suitable for lake assessment. For instance, Šidagytė et al. [69] have demonstrated relationships between ASPT and total phosphorus, biochemical oxygen demand (BOD) and tropho-morphoindex in lakes of Lithuania, while Mavromati et al. [70]—relationships between ASPT and the total phosphorus and shoreline modification in lakes for Greece. Similarly, relationships between ASPT and eutrophication and other pressures have been demonstrated in lakes of Denmark [71]. Currently, ASPT index is used in the lake assessment in Denmark, Greece, Estonia, Latvia, Lithuania, Poland, and Sweden [36] while BMWP in Hungary and Bulgaria [72]. Additionally, in our study, we found that ASPT and BMWP differentiate between natural and altered humic lakes, though we chose to include BMWP (p = 0.002) rather than ASPT (p = 0.45).

4.1.2. Richness/Diversity Metrics

Richness diversity metrics are the widely represented metric category in the lake assessment: almost all the countries use some of these metrics: total taxa richness, Shannon diversity, Margalef diversity or other. Many studies have revealed relationships between total taxa richness/diversity metrics and different human pressures: morphological degradation [73], water-level fluctuation [74], total phosphorus [70] and integrated pressure index [71]. However, in many cases, a pooled taxa number of stressor-sensitive macroinvertebrate orders can be more informative: for instance, CEP taxa richness (Coleoptera, Ephemeroptera, Plecoptera) has shown relatively strong relationships with a range of eutrophication indicators and hydromorphological index in Lithuania, so it was included in the Lithuanian assessment system [69]. Similarly, number of EPTCBO (Ephemeroptera, Plecoptera, Trichoptera, Coleoptera, Bivalvia, Odonata) taxa is used in Denmark, based on demonstrated relationships with eutrophication pressure [71]. In contrast, we found a significant difference in the total taxa richness, number of genera, number of Odonata taxa and the number of Coleoptera taxa (p < 0.001) between humic lakes with natural and altered water level. We selected Coleoptera taxa richness to include into the multimetric index. Additionally, Šidagytė et al. [69] demonstrated strong relationship between Coleoptera taxa richness and different pressure descriptors, both eutrophication and hydromorphological alterations.

4.1.3. Composition/Abundance Metrics

Composition/abundance metrics are widely used in lake assessment systems, mostly expressed as proportion (relative abundance) of specific taxa. Thus, Lithuanian and Danish assessment systems include % COP (Coleoptera, Odonata, Plecoptera) but Greek littoral assessment system—% of Diptera, and German assessment system—% Odonata. Several studies demonstrated between proportion of different taxa % Gastropoda, % Odonata and anthropogenic pressures [69,73]. We have found that % Coleoptera, % Odonata and % Diptera differ strikingly between natural and altered lakes; however, abundance of total community (expressed as ind/m2) was selected as a core metric due to the strong inter-correlations between the proportion and the diversity metrics [53].

4.1.4. Functional Metrics

An alternative to species identity-based methods is the use of functional metrics based on species traits. This approach is recommended by many studies [75,76,77]. However, only few lake assessment systems include functional metrics, e.g., % abundance of feeding type collectors, is included in Austrian and German lake assessment systems, and % abundance of habitat type lithal in the German alpine lake assessment system. Percentage of feeding type predators is used in Sweden to assess acidification [78]. In our study, several functional metrics showed the difference between natural and altered lakes: % littoral, % profundal and % littoral + profundal (p < 0.001). Percentage of organisms with littoral and profundal preference was selected to be included in the core metrics. We assume that this group consists of generalist organisms that might indicate the altered water level.

4.2. Assessment of Hydrological Modifications

Traditionally, lake assessment has focused on eutrophication using primary producers—phytoplankton, macrophytes and, recently phytobenthos communities [79]. However, hydrological alterations, e.g., regulation of lake water level for power production and flood control is among the major anthropogenic disturbances in boreal aquatic ecosystems [80]. Several studies have demonstrated strong effects of water level fluctuations on lake biota, mainly benthic invertebrates, fish fauna and macrophytes [39,80,81]. For instance, Aroviita and Hämäläinen [74] showed a marked decrease in species richness of benthic macroinvertebrate with increasing intensity of water-level regulation, especially for Ephemeroptera, Trichoptera, Coleoptera or Megaloptera. Similarly, changes in species composition were reported in regulated lakes of Ireland, i.e., decrease of Crustaceans, increase of Chironomids, Oligochaetes, and invasive amphipods [82] and in regulated lakes in Italy—increase in mobile and/or feeding opportunistic taxa and decrease in sessile and/or herbivorous taxa [83]. Additionally, Brauns et al. [84] described the decrease in Coleoptera, Odonata, Trichoptera and functional groups of piercers, predators, shredders and xylophagous as the potential effects of lake water level fluctuation in lakes of Germany.
Furthermore, the decrease in benthic invertebrate biomass and/or abundance in the littoral area of regulated lakes has been reported by several studies [85,86,87]. However, other studies did not find any significant effect on numerical abundance of invertebrates [82], probably because water level regulation exerts stronger effect on the biomass of invertebrates than on the numerical abundance affecting larger taxa more [74].
So far, only two lake ecological assessment systems addressed hydromorphological alterations: Slovenian lake assessment system [38] and German assessment system for alpine lakes [36]. However, several lake assessments include multiple pressures including hydromorphological alterations (e.g., [69]). Similar to other studies, we found a marked change in composition and abundance of benthic invertebrates, which further can be used in the development of assessment systems, specifically addressing effect of hydrological alterations in highly humic lakes.

4.3. Assessment of Highly Humic Lakes

Humic lakes constitute a considerable portion of lakes in the boreal zone [1]. However, most assessment systems are developed for clear-water lakes characterized by neutral to alkaline pH, low level of DOM, and water transparency depending on the number of phytoplankton. Humic lakes differ substantially from these clear-water systems and therefore might need different assessment methods [13,17,40]. For instance, Benthic Quality Index (BQI)—a widely used metric in the assessment of lake status in Sweden and Finland [87] was deemed inadequate of assessment of the humic lakes as oxygen depletion and dominance by the tolerant species Chironomus anthracinus and C. plumosus are natural phenomena and not an effect of human impacts [17]. Further, several phytoplankton, macrophyte and fish metrics:
(i)
Classified reference lakes as impacted;
(ii)
Did not differentiate between reference and impacted lakes [17].
This can be explained by the fact that humic lake communities are more tolerant to the environmental fluctuations [88] and are less taxa-rich and diverse comparing to clear-water communities [8]. Further, light limitation due to high level of humic substances plays an important role in these lakes, so several widely used lake assessment parameters as macrophyte colonization depth cannot be used in these lakes [14].
In Latvia, lake assessment system is based on number of taxa, number of EPTCBO taxa, Shannon-Wiener diversity index, ASPT index and acidity index [89]. However, this assessment system was not appropriate for assessment of humic lakes, as also near-natural lakes were classified as impacted according to this system [50]. This problem was encountered also in lakes of Finland [17]. Unsuitability of certain metrics for highly humic lakes can be solved by developing an assessment system targeted to specific human pressure and lake types, as shown by this study. Whether this multimetric index is applicable to humic lakes in other regions needs to be tested in future studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/13/2/223/s1, Table S1: The results of the binary logistic regression (BLR). Z and p values of the chosen metrics between natural and altered highly humic lakes, Table S2: Correlation coefficients of the selected metrics for the highly humic lake multimetric index, Table S3: Anchor values of the macroinvertebrate metrics for a potential multimetric index for the highly humic lake quality assessment, Table S4: Macroinvertebrate taxa found in the studied highly humic lakes.

Author Contributions

Conceptualization, D.O., S.P., I.K.; funding acquisition, D.O.; field works, D.O., A.S., J.J.; sample analysis, D.O., A.S.; data analysis, A.A., D.O., A.S., I.K., J.J.; Writing of original draft, D.O., S.P.; Review and editing, I.K., J.J., A.A., A.S. All authors have read and agreed to the published version of the manuscript.

Funding

The data were collected in the frame of ESF supported project PuReST, No. 1DP/1.1.1.2.0/13/APIA/VIAA/044 and a project supported by The Administration of Latvian Environmental Protection Fund, No. 1-08/202/2018.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the data acquisition from two different projects and the National Surface Water Monitoring program, all requiring different data availability policy.

Acknowledgments

Authors want to express their gratitude to the Latvian Environment, Geology and Meteorology Centre.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keskitalo, J.; Eloranta, P. (Eds.) Limnology of Humic Waters; Backhuys Publishers: Leiden, The Netherlands, 1999; p. 284. [Google Scholar]
  2. Wetzel, R.G. Limnology: Lake and River Ecosystems, 3rd ed.; Academic Press: San Diego, CA, USA, 2001; p. 1006. [Google Scholar]
  3. Salonen, K.; Arvola, L.; Tulonen, T.; Hammar, T.; Metsälä, T.R.; Kankaala, P.; Münster, U. Planktonic food chains of a highly humic lake. I. A mesocosm experiment during the spring primary production maximum. Hydrobiology 1992, 229, 125–142. [Google Scholar] [CrossRef]
  4. Jones, R.I. The influence of humic substances on lacustrine planktonic food-chains. Hydrobiology 1992, 229, 73–91. [Google Scholar] [CrossRef]
  5. Solomon, C.T.; Jones, S.E.; Weidel, B.C.; Buffam, I.; Fork, M.L.; Karlsson, J.; Larsen, S.; Lennon, J.T.; Read, J.S.; Sadro, S.; et al. Ecosystem consequences of changing inputs of terrestrial dissolved organic matter to lakes: Current knowledge and future challenges. Ecosystem 2015, 18, 376–389. [Google Scholar] [CrossRef]
  6. Steinberg, C.E.W.; Kamara, S.; Prokhotskaya, V.Y.; Manusadzianas, L.; Karasyova, T.A.; Timofeyev, M.A.; Jie, Z.; Paul, A.; Meinelt, T.; Farjalla, V.F.; et al. Dissolved humic substances—Ecological driving forces from the individual to the ecosystem level? Freshw. Biol. 2006, 51, 1189–1210. [Google Scholar] [CrossRef]
  7. Lepistö, L.; Holopainen, A.L.; Vuoristo, H. Type-specific and indicator taxa of phytoplankton as a quality criterion for assessing the ecological status of Finnish boreal lakes. Limnology 2004, 34, 236–248. [Google Scholar] [CrossRef] [Green Version]
  8. Toivonen, H.; Huttunen, P. Aquatic macrophytes and ecological gradients in 57 small lakes in southern Finland. Aquat. Bot. 1995, 51, 197–221. [Google Scholar] [CrossRef]
  9. Vesterinen, J.; Devlin, S.P.; Syväranta, J.; Jones, R.I. Influence of littoral periphyton on whole-lake metabolism relates to littoral vegetation in humic lakes. Ecology 2017, 98, 3074–3085. [Google Scholar] [CrossRef] [Green Version]
  10. Sarvala, J.; Kankaala, P.; Zingel, P.; Arvola, L. Zooplankton. In Limnology of Humic Waters; Keskitalo, J., Eloranta, P., Eds.; Backhuys Publishers: Leiden, The Netherlands, 1999; pp. 173–191. [Google Scholar]
  11. Rask, M.; Viljanen, M.; Sarval, J. Humic lakes as fish habitats. In Limnology of Humic Waters; Keskitalo, J., Eloranta, P., Eds.; Backhuys Publishers: Leiden, The Netherlands, 1999; pp. 209–224. [Google Scholar]
  12. Ptacnik, R.; Solimini, A.G.; Brettum, P. Performance of a new phytoplankton composition metric along a eutrophication gradient in Nordic lakes. Hydrobiology 2009, 633, 75–82. [Google Scholar] [CrossRef]
  13. Rask, M.; Vuori, K.-M.; Hämäläinen, H.; Järvinen, M.; Hellsten, S.; Mykrä, H.; Arvola, L.; Ruuhijärvi, J.; Jyväsjärvi, J.; Kolari, I.; et al. Ecological classification of large lakes in Finland: Comparison of classification approaches using multiple quality elements. Hydrobiology 2011, 660, 37–47. [Google Scholar] [CrossRef]
  14. Søndergaard, M.; Phillips, G.; Hellsten, S.; Kolada, A.; Ecke, F.; Mäemets, H.; Mjelde, M.; Azzella, M.M.; Oggioni, A. Maximum growing depth of submerged macrophytes in European lakes. Hydrobiology 2013, 704, 165–177. [Google Scholar] [CrossRef]
  15. Lepistö, L.; Holopainen, A.L.; Vuoristo, H.; Rekolainen, S. Phytoplankton assemblages as a criterion in the ecological classification of lakes in Finland. Boreal Environ. Res. 2006, 11, 35–44. [Google Scholar]
  16. Birk, S.; Ecke, F. The potential of remote sensing in ecological status assessment of coloured lakes using aquatic plants. Ecol. Indic. 2014, 46, 398–406. [Google Scholar] [CrossRef]
  17. Alahuhta, J.; Vuori, K.M.; Hellsten, S.; Järvinen, M.; Olin, M.; Rask, M.; Palomäki, A. Defining the ecological status of small forest lakes using multiple biological quality elements and palaeolimnological analysis. Fund. Appl. Limnol. 2009, 175, 203–216. [Google Scholar] [CrossRef]
  18. Baars, J.R.; Murray, D.A.; Hannigan, E.; Kelly-Quinn, M. Macroinvertebrate assemblages of small upland peatland lakes in Ireland. Biol. Environ. 2014, 114, 233–248. [Google Scholar] [CrossRef]
  19. Armstrong, A.; Holden, J.; Kay, P.; Francis, B.; Foulger, M.; Gledhill, S.; McDonald, A.T.; Walker, A. The impact of peatland drain-blocking on dissolved organic carbon loss and discolouration of water; results from a national survey. J. Hydrol. 2010, 381, 112–120. [Google Scholar] [CrossRef]
  20. Beadle, J.M.; Brown, L.E.; Holden, J. Biodiversity and ecosystem functioning in natural bog pools and those created by rewetting schemes. Wiley Interdiscip. Rev. Water 2015, 2, 65–84. [Google Scholar] [CrossRef]
  21. Lepistö, L.; Saura, M. Effects of forest fertilization on phytoplankton in a boreal brown-water lake. Boreal Environ. Res. 1998, 3, 33–43. [Google Scholar]
  22. Drinan, T.J.; Graham, C.T.; O’Halloran, J.; Harrison, S.S.C. The impact of conifer plantation forestry on the Chydoridae (Cladocera) communities of peatland lakes. Hydrobiology 2013, 700, 203–219. [Google Scholar] [CrossRef] [Green Version]
  23. Deininger, A.; Jonsson, A.; Karlsson, J.; Bergström, A.K. Pelagic food webs of humic lakes show low short-term response to forest harvesting. Ecol. Appl. 2018, 29, 1–12. [Google Scholar] [CrossRef]
  24. Drinan, T.J.; Foster, G.N.; Nelson, B.H.; O’Halloran, J.; Harrison, S.S.C. Macroinvertebrate assemblages of peatland lakes: Assessment of conservation value with respect to anthropogenic land-cover change. Biol. Conserv. 2013, 158, 175–187. [Google Scholar] [CrossRef]
  25. Graham, C.T.; Drinan, T.J.; Harrison, S.S.; O’Halloran, J. Relationship between plantation forest and brown trout growth, energetics and population structure in peatland lakes in western Ireland. For. Ecol. Manag. 2014, 321, 71–80. [Google Scholar] [CrossRef]
  26. European Council Directive. European Commission 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. OJEC 2000, 327, 1–72. [Google Scholar]
  27. Poikane, S.; Salas Herrero, F.; Kelly, M.G.; Borja, A.; Birk, S.; van de Bund, W. European aquatic ecological assessment methods: A critical review of their sensitivity to key pressures. Sci. Total Environ. 2020, 740, 140075. [Google Scholar] [CrossRef] [PubMed]
  28. Poikane, S.; Zampoukas, N.; Borja, A.; Davies, S.P.; van de Bund, W.; Birk, S. Intercalibration of aquatic ecological assessment methods in the European Union: Lessons learned and way forward. Environ. Sci. Policy 2014, 44, 237–246. [Google Scholar] [CrossRef]
  29. Poikane, S.; Birk, S.; Böhmer, J.; Carvalho, L.; de Hoyos, C.; Gassner, H.; Hellsten, S.; Kelly, M.; Solheim, A.L.; Olin, M.; et al. A hitchhiker’s guide to European lake ecological assessment and intercalibration. Ecol. Indic. 2015, 52, 533–544. [Google Scholar] [CrossRef]
  30. Carvalho, L.; Poikane, S.; Solheim, A.L.; Phillips, G.; Borics, G.; Catalan, J.; De Hoyos, C.; Drakare, S.; Dudley, B.J.; Järvinen, M.; et al. Strength and uncertainty of phytoplankton metrics for assessing eutrophication impacts in lakes. Hydrobiology 2013, 704, 127–140. [Google Scholar] [CrossRef] [Green Version]
  31. Phillips, G.; Lyche-Solheim, A.; Skjelbred, B.; Mischke, U.; Drakare, S.; Free, G.; Järvinen, M.; de Hoyos, C.; Morabito, G.; Poikane, S.; et al. A phytoplankton trophic index to assess the status of lakes for the Water Framework Directive. Hydrobiology 2013, 704, 75–95. [Google Scholar] [CrossRef] [Green Version]
  32. Penning, W.E.; Dudley, B.; Mjelde, M.; Hellsten, S.; Hanganu, J.; Kolada, A.; van den Berg, M.; Poikane, S.; Phillips, G.; Willby, N.; et al. Using aquatic macrophyte community indices to define the ecological status of European lakes. Aquat. Ecol. 2008, 42, 253–264. [Google Scholar] [CrossRef]
  33. Poikane, S.; Portielje, R.; Denys, L.; Elferts, D.; Kelly, M.; Kolada, A.; Mäemets, H.; Phillips, G.; Søndergaard, M.; Willby, N.; et al. Macrophyte assessment in European lakes: Diverse approaches but convergent views of ‘good’ ecological status. Ecol. Indic. 2018, 94, 185–197. [Google Scholar] [CrossRef]
  34. Kelly, M.; Urbanic, G.; Acs, E.; Bennion, H.; Bertrin, V.; Burgess, A.; Denys, L.; Gottschalk, S.; Kahlert, M.; Karjalainen, S.M.; et al. Comparing aspirations: Intercalibration of ecological status concepts across European lakes for littoral diatoms. Hydrobiology 2014, 734, 125–141. [Google Scholar] [CrossRef] [Green Version]
  35. Lyche-Solheim, A.; Feld, C.K.; Birk, S.; Phillips, G.; Carvalho, L.; Morabito, G.; Mischke, U.; Willby, N.; Søndergaard, M.; Hellsten, S.; et al. Ecological status assessment of European lakes: A comparison of metrics for phytoplankton, macrophytes, benthic invertebrates and fish. Hydrobiology 2013, 704, 57–74. [Google Scholar] [CrossRef] [Green Version]
  36. Poikane, S.; Johnson, R.K.; Sandin, L.; Schartau, A.K.; Solimini, A.G.; Urbanič, G.; Arbačiauskas, K.; Aroviita, J.; Gabriels, W.; Miler, O.; et al. Benthic macroinvertebrates in lake ecological assessment: A review of methods, intercalibration and practical recommendations. Sci. Total Environ. 2016, 543, 123–134. [Google Scholar] [CrossRef] [PubMed]
  37. EEA. European Waters—Assessment of Status and Pressures; EEA Report; European Environment Agency: København, Denmark, 2018; pp. 1–90. [Google Scholar] [CrossRef]
  38. Urbanič, G. A Littoral Fauna Index for assessing the impact of lakeshore alterations in Alpine lakes. Ecohydrology 2014, 7, 703–716. [Google Scholar] [CrossRef]
  39. Mjelde, M.; Hellsten, S.; Ecke, F. A water level drawdown index for aquatic macrophytes in Nordic lakes. Hydrobiology 2013, 704, 141–151. [Google Scholar] [CrossRef]
  40. Olin, M.; Rask, M.; Ruuhijärvi, J.; Tammi, J. Development and evaluation of the Finnish fish-based lake classification method. Hydrobiology 2013, 713, 149–166. [Google Scholar] [CrossRef]
  41. Van Duinen, G.A.; Brock, A.M.; Kuper, J.T.; Leuven, R.S.E.W.; Peeters, T.M.J.; Roelofs, J.G.M.; van der Velde, G.; Verberk, W.C.E.P.; Esselink, H. Do restoration measures rehabilitate fauna diversity in raised bogs? A comparative study on aquatic macroinvertebrates. Wetl. Ecol. Manag. 2003, 11, 447–459. [Google Scholar] [CrossRef]
  42. Hannigan, E.; Kelly-Quinn, M. Composition and structure of macroinvertebrate communities in contrasting open-water habitats in Irish peatlands: Implications for biodiversity conservation. Hydrobiology 2012, 692, 19–28. [Google Scholar] [CrossRef]
  43. Walker, I.R.; Fernando, C.H.; Paterson, C.G. Associations of Chironomidae (Diptera) of shallow, acid, humic lakes and bog pools in Atlantic Canada, and a comparison with an earlier paleoecological investigation. Hydrobiology 1985, 120, 11–22. [Google Scholar] [CrossRef]
  44. Pakulnicka, J.; Zawal, A. Model of disharmonic succession of dystrophic lakes based on aquatic beetle fauna (Coleoptera). Mar. Freshw. Res. 2019, 70, 195–211. [Google Scholar] [CrossRef]
  45. Pakalne, M.; Kalnina, L. Mire ecosystems in Latvia. Stapfia 2005, 85, 147–174. [Google Scholar]
  46. Montanarella, L.; Jones, R.J.A.; Hiederer, R. The distribution of peatland in Europe. Mires Peat 2006, 1, 1–10. [Google Scholar]
  47. European Council Directive. Council Directive 92/43/EEC on the conservation of natural habitats and of wild fauna and flora. OJEC 1992, 206, 7–42. [Google Scholar]
  48. Pakalne, M.; Kalniņa, L. Mires in Latvia. Suo 2000, 51, 213–226. [Google Scholar]
  49. Klavins, M.; Kokorite, I.; Springe, G.; Skuja, A.; Parele, E.; Rodinov, V.; Druvietis, I.; Strāķe, S.; Urtans, A. Water quality in cutaway peatland lakes in Seda mire, Latvia. Ecohydrol. Hydrobiol. 2010, 10, 61–70. [Google Scholar] [CrossRef]
  50. Latvian Environment, Geology and Meteorology Centre. River Basin Management Plans 2016–2021. Available online: https://videscentrs.lvgmc.lv/lapas/udens-apsaimniekosana-un-pludu-parvaldiba (accessed on 7 December 2020).
  51. American Public Health Association; American Water Works Association; Water Environment Federation. Standard Methods for the Examination of Water and Wastewater, 21th ed.; American Public Health Association: Washington, DC, USA, 2005; p. 1496. [Google Scholar]
  52. AQEM European Stream Assessment Program. English Manual, Version 2.3. Available online: http://www.aqem.de/mains/products.php (accessed on 7 December 2020).
  53. Hering, D.; Feld, C.K.; Moog, O.; Ofenbock, T. Cook book for the development of a Multimetric Index for biological condition of aquatic ecosystems: Experiences from the European AQEM and STAR projects and related initiatives. Hydrobiology 2006, 566, 311–324. [Google Scholar] [CrossRef]
  54. Sokal, R.R.; Rohlf, F.J. Biometry: The Principles and Practice of Statistics in Biological Research, 3rd ed.; W.H. Freeman and Co: New York, NY, USA, 1995; p. 880. [Google Scholar]
  55. R Core Team. R: A Language and Environment for Statistical Computing. 2020. Available online: http://www.r-project.org/index.html (accessed on 3 November 2020).
  56. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  57. Convention on the Conservation of European Wildlife and Natural Habitats (Bern, Switzerland). 1979. Available online: https://rm.coe.int/1680078aff (accessed on 12 January 2021).
  58. Desrochers, A.; van Duinen, G.J. Peatland fauna. In Boreal Peatland Ecosystems, Ecological Studies 188; Wieder, R.K., Vitt, D.H., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 67–100. [Google Scholar]
  59. Ewald, M.L.; Feminella, J.W.; Lenertz, K.K.; Henry, R.P. Acute physiological responses of the freshwater snail Elimia flava (Mollusca: Pleuroceridae) to environmental pH and calcium. Comp. Biochem. Physiol. Part C: Toxicol. Pharmacol. 2009, 150, 237–245. [Google Scholar] [CrossRef]
  60. Mäkelä, T.P.; Oikari, A.O.J. The effects of low water pH on the ionic balance in freshwater mussel Anodonta anatina L. Ann. Zool. Fennici. 1992, 29, 169–175. [Google Scholar]
  61. Schartau, A.K.; Moe, S.J.; Sandin, L.; McFarland, B.; Raddum, G.G. Macroinvertebrate indicators of lake acidification: Analysis of monitoring data from UK, Norway and Sweden. Aquat. Ecol. 2008, 42, 293–305. [Google Scholar] [CrossRef]
  62. Soldán, T.; Bojková, J.; Vrba, J.; Bitušík, P.; Chvojka, P.; Papáček, M.; Peltanová, J.; Sychra, J.; Tátosová, J. Aquatic insects of the Bohemian Forest glacial lakes: Diversity, long-term changes, and influence of acidification. Silva Gabreta 2012, 18, 123–283. [Google Scholar]
  63. Karr, J.R.; Chu, E.W. Biological Monitoring and Assessment: Using Multimetric Indexes Effectively; EPA 235-R97-001; US Environmental Protection Agency: Washington, DC, USA, 1997; p. 149.
  64. Armitage, P.D.; Moss, D.; Wright, J.F.; Furse, M.T. The performance of a new biological water quality score system based on macroinvertebrates over a wide range of unpolluted running-water sites. Water Res. 1983, 17, 333–347. [Google Scholar] [CrossRef]
  65. Davy-Bowker, J.; Murphy, J.F.; Rutt, G.P.; Steel, J.E.C.; Furse, M.T. The development and testing of a macroinvertebrate biotic index for detecting the impact of acidity on streams. Arch. Hydrobiol. 2005, 163, 383–403. [Google Scholar] [CrossRef]
  66. Wiederholm, T. Use of benthos in lake monitoring. J. Water Pollut. Con. Fed. 1980, 52, 537–547. [Google Scholar]
  67. Miler, O.; Porst, G.; McGoff, E.; Pilotto, F.; Donohue, L.; Jurca, T.; Solimini, A.G.; Sandin, L.; Irvin, K.; Aroviita, J.; et al. Morphological alterations of lake shores in Europe: A multimetric ecological assessment approach using benthic macroinvertebrates. Ecol. Indic. 2013, 34, 398–410. [Google Scholar] [CrossRef]
  68. Hawkes, H.A. Origin and development of the biological monitoring working party score system. Water Res. 1998, 32, 964–968. [Google Scholar] [CrossRef]
  69. Šidagytė, E.; Višinskienė, G.; Arbačiauskas, K. Macroinvertebrate metrics and their integration for assessing the ecological status and biocontamination of Lithuanian lakes. Limnologica 2013, 43, 308–318. [Google Scholar] [CrossRef]
  70. Mavromati, E.; Kemitzoglou, D.; Tsiaoussi, V.; Lazaridou, M. Report on the Development of the National Method for the Assessment of Ecological Status of Natural Lakes in Greece, with the Use of Littoral Benthic Invertebrates; Greek Biotope/Wetland Centre and Special Secretariat for the Natural Environment and Waters, Ministry of Environment and Energy: Thermi, Greece, 2020. [Google Scholar]
  71. Wiberg-Larsen, P.; Rasmussen, J.J. Revised Danish Macroinvertebrate Index for Lakes—A Method to Assess Ecological Quality. Aarhus University, DCE—Danish Centre for Environment and Energy. 2020. Available online: http://dce2.au.dk/pub/SR373.pdf (accessed on 17 January 2020).
  72. Böhmer, J.; Chiriac, G.; Varbiro, G.; Wolfram, G.; Poikane, S. Intercalibrating the National Classifications of Ecological Status for Eastern Continental Lakes: Biological Quality Element: Benthic Invertebrates; EUR 29342 EN; Publications Office of the European Union: Luxembourg, 2018; p. 54. [Google Scholar] [CrossRef]
  73. Urbanič, G.; Petkovska, V.; Pavlin, M. The relationship between littoral benthic invertebrates and lakeshore modification pressure in two alpine lakes. Fundam. Appl. Limnol. 2012, 180, 157–173. [Google Scholar] [CrossRef]
  74. Aroviita, J.; Hämäläinen, H. The impact of water-level regulation on littoral macroinvertebrate assemblages in boreal lakes. Hydrobiology 2008, 613, 45–56. [Google Scholar] [CrossRef]
  75. Friberg, N.; Bonada, N.; Bradley, D.C.; Dunbar, M.J.; Edwards, F.K.; Grey, J.; Hayes, R.B.; Hildrew, A.G.; Lamouroux, N.; Trimmer, M.; et al. Biomonitoring of human impacts in freshwater ecosystems. The good, the bad and the ugly. Adv. Ecol. Res. 2011, 44, 1–68. [Google Scholar] [CrossRef]
  76. Reyjol, Y.; Argillier, C.; Bonne, W.; Borja, A.; Buijse, A.D.; Cardoso, A.C.; Daufresne, M.; Kernan, M.; Ferreira, M.T.; Poikane, S.; et al. Assessing the ecological status in the context of the European Water Framework Directive: Where do we go now? Sci. Total Environ. 2014, 497, 332–344. [Google Scholar] [CrossRef] [PubMed]
  77. Brucet, S.; Poikane, S.; Lyche-Solheim, A.; Birk, S. Biological assessment methods for European lakes: Ecological rationale and human impacts. Freshw. Biol. 2013, 58, 1106–1115. [Google Scholar] [CrossRef]
  78. Johnson, R.K.; Goedkoop, W. Bedömningsgrunder för Bottenfauna i Sjöar Och Vattendrag: Användarmanual och Bakgrundsdokument; Sveriges lantbruksuniversitet: Uppsala, Sweden, 2007; p. 84. [Google Scholar]
  79. Kelly, M.G.; Birk, S.; Willby, N.J.; Denys, L.; Drakare, S.; Kahlert, M.; Karjalainen, S.M.; Marchetto, A.; Pitt, J.A.; Urbanič, G.; et al. Redundancy in the ecological assessment of lakes: Are phytoplankton, macrophytes and phytobenthos all necessary? Sci. Total Environ. 2016, 568, 594–602. [Google Scholar] [CrossRef] [PubMed]
  80. Sutela, T.; Vehanen, T.; Rask, M. Assessment of the ecological status of regulated lakes: Stressor-specific metrics from littoral fish assemblages. Hydrobiology 2011, 675, 55–64. [Google Scholar] [CrossRef]
  81. Sutela, T.; Aroviita, J.; Keto, A. Assessing ecological status of regulated lakes with littoral macrophyte, macroinvertebrate and fish assemblages. Ecol. Indic. 2013, 24, 185–192. [Google Scholar] [CrossRef]
  82. Evtimova, V.V.; Donohue, I. Water-level fluctuations regulate the structure and functioning of natural lakes. Freshw. Biol. 2016, 61, 251–264. [Google Scholar] [CrossRef]
  83. Mastrantuono, L.; Pilotto, F.; Rossopinti, A.; Bazzanti, M.; Solimini, G. Response of littoral macroinvertebrates to morphological disturbances in Mediterranean lakes: The case of Lake Piediluco (Central Italy). Fund. Appl. Limnol. 2015, 186, 297–310. [Google Scholar] [CrossRef]
  84. Brauns, M.; Garcia, X.F.; Pusch, M. Potential effects of water level fluctuations on littoral invertebrates in lowland lakes. Hydrobiology 2008, 613, 5–12. [Google Scholar] [CrossRef]
  85. Palomäki, R. Response by macrozoobenthos biomass to water level regulation in some Finnish lake littoral zones. Hydrobiology 1994, 286, 17–26. [Google Scholar] [CrossRef]
  86. Tikkanen, P.; Kantola, L.; Niva, T.; Hellsten, S.; Alasaarela, E. Ecological Aspects of Lake Regulation in Northern Finland. Part 3. Macrozoobenthos and Feeding of Fish; VTT Research Notes 987; Technical Research Centre of Finland: Espoo, Finland, 1989; (in Finnish with English). [Google Scholar]
  87. Jyväsjärvi, J.; Nyblom, J.; Hämäläinen, H. Palaeolimnological validation of estimated reference values for a lake profundal macroinvertebrate metric (Benthic Quality Index). J. Paleolimnol. 2010, 44, 253–264. [Google Scholar] [CrossRef]
  88. Tammi, J.; Lappalainen, A.; Mannio, J.; Rask, M.; Vuorenmaa, J. Effects of eutrophication on fish and fisheries in Finnish lakes: A survey based on random sampling. Fish. Manag. Ecol. 1999, 6, 173–186. [Google Scholar] [CrossRef]
  89. Skuja, A.; Ozoliņš, D. Fitting New Method—Latvian Lake Macroinvertebrate Multimetric Index (LLMMI) to Results of Central—Baltic Geographical Intercalibration Group (CB—GIG) Lake Benthic Macroinvertebrate Intercalibration; Report; Institute of Biology, University of Latvia: Salaspils, Latvia, 2016; p. 12. [Google Scholar]
Figure 1. The macroinvertebrate sampling sites and distribution of fens in Latvia. Black triangles: highly humic lakes with altered water level; black dots: highly humic lakes with natural or restored water level.
Figure 1. The macroinvertebrate sampling sites and distribution of fens in Latvia. Black triangles: highly humic lakes with altered water level; black dots: highly humic lakes with natural or restored water level.
Water 13 00223 g001
Figure 2. Lake Soku with receded water level and visible open peat outcrops at the shoreline in May 2019.
Figure 2. Lake Soku with receded water level and visible open peat outcrops at the shoreline in May 2019.
Water 13 00223 g002
Figure 3. Lake Deguma with natural water level and mire vegetation at the shoreline in May 2019.
Figure 3. Lake Deguma with natural water level and mire vegetation at the shoreline in May 2019.
Water 13 00223 g003
Figure 4. The evaluation graphs of (a) composition/abundance; (b) richness/diversity; (c) sensitive/tolerance and (d) functional metric groups of highly humic lakes.
Figure 4. The evaluation graphs of (a) composition/abundance; (b) richness/diversity; (c) sensitive/tolerance and (d) functional metric groups of highly humic lakes.
Water 13 00223 g004
Figure 5. The dominance structure of macroinvertebrate assemblages in highly humic lakes of Latvia.
Figure 5. The dominance structure of macroinvertebrate assemblages in highly humic lakes of Latvia.
Water 13 00223 g005
Figure 6. Differences of macroinvertebrate (a) abundance, (b) taxa richness, (c) BMWP (Biological Monitoring Working Party) Score, (d) Coleoptera taxa, (e) ETCO (Ephemeroptera, Trichoptera, Coleoptera, Odonata), (f) littoral + profundal preference, (g) % Diptera, (h) Trichoptera taxa, (i) Odonata taxa and (j) Heteroptera taxa between bog lakes with altered and natural level.
Figure 6. Differences of macroinvertebrate (a) abundance, (b) taxa richness, (c) BMWP (Biological Monitoring Working Party) Score, (d) Coleoptera taxa, (e) ETCO (Ephemeroptera, Trichoptera, Coleoptera, Odonata), (f) littoral + profundal preference, (g) % Diptera, (h) Trichoptera taxa, (i) Odonata taxa and (j) Heteroptera taxa between bog lakes with altered and natural level.
Water 13 00223 g006aWater 13 00223 g006b
Figure 7. The model comparison of the humic lake multimetric indices at the lakes with natural and altered water level. The model caption consists of model_Coleoptera (c = yes; a = no) percentile anchors.
Figure 7. The model comparison of the humic lake multimetric indices at the lakes with natural and altered water level. The model caption consists of model_Coleoptera (c = yes; a = no) percentile anchors.
Water 13 00223 g007
Figure 8. The distribution of values of the newly developed highly humic lake multimetric index at the lakes with natural and altered water level.
Figure 8. The distribution of values of the newly developed highly humic lake multimetric index at the lakes with natural and altered water level.
Water 13 00223 g008
Table 1. Mean annual chemical and environmental parameters at the studied highly humic lakes (Cond—electric conductivity, TN—total nitrogen, TP—total phosphorus).
Table 1. Mean annual chemical and environmental parameters at the studied highly humic lakes (Cond—electric conductivity, TN—total nitrogen, TP—total phosphorus).
LakesAltered Water Level (+)YearpHCond μS/cmColour mg Pt/LTN mg/LTP mg/L
Cenas Mire Lake 1 20153.93281240.950.017
Cenas Mire Lake 2 20154.44291440.990.020
Cenas Mire Lake 3 20154.49261140.900.019
Cenas Mire Lake 4 20153.75321890.830.021
Cenas Mire Lake 5 20153.46433040.920.020
Melnais Lake Mire Lake 1+20153.68653931.310.019
Melnais Lake Mire Lake 2+20153.53494021.250.022
Melnais Lake Mire Lake 3+20153.58443651.170.022
Melnais Lake Mire Lake 4+20153.43455051.360.017
Melnais Lake Mire Lake 5+20153.35486661.680.028
Deguma Lake 20175.09312220.950.032
Orlovas Lake 20175.42212050.870.037
Pieslaista Lake 20175.02362380.650.061
Ramatas Lielezers Lake+20176.09231340.650.032
Soku Lake+20175.88251300.430.035
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ozoliņš, D.; Skuja, A.; Jēkabsone, J.; Kokorite, I.; Avotins, A.; Poikane, S. How to Assess the Ecological Status of Highly Humic Lakes? Development of a New Method Based on Benthic Invertebrates. Water 2021, 13, 223. https://doi.org/10.3390/w13020223

AMA Style

Ozoliņš D, Skuja A, Jēkabsone J, Kokorite I, Avotins A, Poikane S. How to Assess the Ecological Status of Highly Humic Lakes? Development of a New Method Based on Benthic Invertebrates. Water. 2021; 13(2):223. https://doi.org/10.3390/w13020223

Chicago/Turabian Style

Ozoliņš, Dāvis, Agnija Skuja, Jolanta Jēkabsone, Ilga Kokorite, Andris Avotins, and Sandra Poikane. 2021. "How to Assess the Ecological Status of Highly Humic Lakes? Development of a New Method Based on Benthic Invertebrates" Water 13, no. 2: 223. https://doi.org/10.3390/w13020223

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