Next Article in Journal
Hydro-Energy Suitability of Rivers Regarding Their Hydrological and Hydrogeological Characteristics
Next Article in Special Issue
Impact of EU Environmental Policy Implementation on the Quality and Status of Greek Rivers
Previous Article in Journal
Grading Evaluation of the Structural Connectivity of River System Networks Based on Ecological Functions, and a Case Study of the Baiyangdian Wetland, China
Previous Article in Special Issue
Suggested Sampling Methodology for Lake Benthic Macroinvertebrates under the Requirements of the European Water Framework Directive
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

“One Out–All Out” Principle in the Water Framework Directive 2000—A New Approach with Fuzzy Method on an Example of Greek Lakes

1
Department of Civil Engineering, School of Engineering, Democritus University of Thrace, 67100 Xanthi, Greece
2
Department of Zoology, School of Biology, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
3
Greek Biotope/Wetland Centre (EKBY), The Goulandris Natural History Museum, 57001 Thermi, Greece
*
Authors to whom correspondence should be addressed.
Water 2021, 13(13), 1776; https://doi.org/10.3390/w13131776
Submission received: 12 May 2021 / Revised: 23 June 2021 / Accepted: 24 June 2021 / Published: 27 June 2021

Abstract

:
The “One Out–All Out” (OOAO) principle imposed by the WFD selects the worst ecological status assessed by different biological quality elements (BQEs). Since it is a precautionary rule that can lead to problems of underestimation of the overall status, its amendment has been a matter of debate for WFD 20+. The use of fuzzy methods that express the functional relationships between variables in ecology and management has been gaining more ground recently. Here is attempted the inclusion of a fuzzy regression among the frequently monitored BQE (phytoplankton) and the outcome of OOAO application in six Greek lakes. The latter was determined by the comparison of four BQE indices in order to assess the extent to which BQEs might underpin the optimal/actual qualitative classification of a waterbody. This approach encompasses the uncertainty and the possibility to broaden the acceptable final EQR based on the character and status of each lake. We concluded that the fuzzy OOAO is an approach that seems to allow a better understanding of the WFD implementation and case-specific evaluation, including the uncertainty in classification as an asset. Moreover, it offers a deeper understanding through self-learning processes based on the existing datasets.

1. Introduction

The Water Framework Directive (WFD) [1], which was formally adopted in 2000, requires that all European Member States (MSs) have to assess the ecological status or potential of their inland, transitional, and coastal waterbodies. A key feature of the WFD is the use of the catchment-scale as the management unit, thereby allowing for the integrated effects of all pressures on waterbodies to be considered [2]. In particular, it introduces an innovative approach to manage and protect aquatic ecosystems in a holistic way rather than focusing only on specific aspects of water quality. Classification of ecological status or potential is based on different biological quality elements (BQEs) representing main ecosystem components, such as phytoplankton, benthic macroinvertebrates, other aquatic flora (macrophytes, phytobenthos), and fish [1].
Up to now, WFD complemented by its daughter and sister directives (as the Environmental Quality Standards Directive [3], the Groundwater Directive [4], the Flood Directive [5]), has achieved the set-up of a management framework across all European Union (EU) waterbodies, including lakes, mitigating their deterioration. However, at the same time, around 57% of rivers and 44% of lakes in EU have failed to achieve good ecological status or potential [6], while according to the WFD Fitness check [7], no substantial progress in waterbodies’ overall status has been made between the first and the second river basin management cycles. This means that waterbodies classified at less than good ecological status or potential must be restored, and further management responses should be applied to improve the overall system’s health. These responses are presented under the term “Program of Measures, PoMs” in the River Basin Management Plans (RBMPs) and target the pressures and their drivers contributing to ecosystem dysfunction [8]. The overall system’s classification is based on the “One Out–All Out” (OOAO) principle, meaning that the worst quality of any of the BQEs used in the assessment regulates the overall ecological status of a waterbody.
The majority of PoMs require significant cost while, often, the lack of financial resources for implementing a long list of measures can be a factor for not achieving better results. Thus, the diagnosis of the efficacy of each BQE in reflecting the actual status, which in turn guides the classification of each waterbody, is of high importance. Misclassification could lead either to over-precautionary results and imposition of restoration costs disproportionate or neglect of the need for strict measures to target pressures. It is acknowledged that it is more difficult to make progress visible in the WFD goals due to the OOAO principle [7]. Indeed, the OOAO principle was originally set up as a precautionary principle. The ecological basis for the debate on the usefulness of the OOAO principle has been clearly demonstrated by several researchers [9,10,11]. They state that the OOAO principle is well suited when different stressors are responsible for the degradation of the individual BQEs and further, that the tendency for misclassification increases with the BQEs uncertainty. Waterbodies are either affected by predominant pressures (e.g., organic pollution or eutrophication) or by multiple pressures, which can interact in additive, synergistic, or antagonistic ways (e.g., nutrient enrichment, hydromorphological degradation, toxic substances, overfishing) [12]. At the same time, the ecological knowledge of the response of different organism groups to pressures varies across Europe [13]. In addition, new pressures are appearing (e.g., microplastic pollution, pharmaceuticals, light and noise, freshwater salinization) [14].
The uncertainty of the BQEs and their significance for the ecological status classification have been addressed in many studies [15,16,17,18,19,20]. Currently, in relevant EU working groups, there are discussions on the use of a new set of indicators, which will be based on BQEs, and can serve as a supplementary guide to overall status or potential. For instance, in assessing the progress of ecological status, since the beginning of monitoring and the implementation of PoMs, the comparison of OOAO with other integrative methods in coastal and transitional waterbodies showed a low (18%) disagreement for the coastal and a much higher (58%) for the transitional [21]. This can be explained by the fact that one or more BQEs appear to have a low or a moderate correlation with various stressors, thus often describing similar inconsistencies [22]. Thus, the establishment of alternative combination rules is suggested [23], taking also into consideration the type I (i.e., a waterbody is below good status even if the waterbody in reality has good status) and type II (i.e., reducing the likelihood that a waterbody is classified as good status when, in reality, it is below good status) errors, as highlighted by Hering et al. [13].
Based on the above, Borja and Rodriguez [21] have suggested that MSs should avoid adopting this simplistic principle, although it may present an effective indication in the ecological status assessment. Thus, amendments on the application of OOAO or alternative approaches should be further considered. In a recent paper reporting the results of a survey on the suitable options of how the OOAO principle should be implemented, Carvalho et al. [11] show that the most popular responses opted for the boosting of selected BQEs and the emphasis on progress reporting. Therefore, the linkage of pressures with ecological status sets the basis for establishing management targets and restoration measures [24]. According to Carvalho et al. [11], the only state adopting management decisions dependent on the strongest evidence of BQE is the United Kingdom.
So far, the fundamental target of WFD implementation has been the development of assessment methods for all BQEs required in the WFD. For EU lakes, a wide variety of assessment methods have been developed based on different BQEs that respond to specific pressures [15,25,26] or multiple pressures [27,28]. In Greece, until recently, there were gaps in the assessment of the ecological quality of lakes with WFD-compatible indices. Five national indices [29,30,31,32,33], applicable in the broader Mediterranean area, were developed for all BQEs reflecting the local peculiarities. It cannot be generally stated which BQEs underestimate or overestimate the actual status. It is empirically observed that the qualitative assessment depends more on the lake type, its special traits, and the type and intensity of pressures that is subject to.
Fuzzy logic in ecology has been gaining ground during the last two decades. The natural complexity can be reflected in fuzzy ecological modeling and scenarios simulation [34], while such “unconventional” data analyses can be applied for sustainability evaluation [35]. These approaches allow for overcoming the usual obstacle of precise quantification processes [36]. Moreover, fuzzy logic can assist in management planning and decision making by enhancing objectivity [37]. The fuzzy regression tool is useful in expressing functional relationships between variables, especially when the available dataset is insufficient [38]. Kitsikoudis et al. [39] employed fuzzy regression and set a fuzzy band to produce the lower and upper acceptable limits (left and right boundaries) of critical Shields stress. Thus, the ambiguity of selecting a threshold can be avoided, and a smoother transition to the actual state can be provided [40,41].
In this study, we first assessed the strictness of the OOAO application based on different ecological assessment approaches and BQEs (phytoplankton, macrophytes, and benthic macroinvertebrates). Then, we examined a method utilizing fuzzy logic to express the functional relationships between two variables. fuzzy regression was applied for phytoplankton (as the main BQE for lakes and more frequently monitored) and OOAO that derived from all available BQE indices. This approach took place to assess the extent to which different BQEs can underpin the optimal/actual qualitative classification of a waterbody.

2. Materials and Methods

2.1. Study Area

As a case study, we used the BQEs from six natural lakes in northern and central Greece (Figure 1). These lakes are included in the National Water Monitoring Network, while one of them is transboundary (Doirani), shared with the Republic of North Macedonia. According to the Hellenic typology [29], these lakes belong to two types: (a) warm monomictic, deep, natural lakes with mean depth >9 m (type GR-DNL, 3 lakes) and (b) polymictic, shallow lakes with mean depth 3–9 m (type GR-SNL, 3 lakes). The studied lakes are characterized by a variability in their limnological and trophic state attributes (Table 1).

2.2. Biological Assessment Methods

In the present study, each BQE was assessed by one index, except for benthic macroinvertebrates, as follows:
Phytoplankton: Hellenic Phytoplankton Assessment System (HeLPhy), which is composed of four metrics, i.e., chlorophyll-a, total biovolume, modified Nygraard index, and biovolume of cyanobacteria and responds to eutrophication [29].
Macrophytes: Hellenic Lake Macrophytes (HeLM), which consists of two metrics, i.e., trophic index and maximum depth of colonization and responds to eutrophication [30].
Benthic macroinvertebrates: Greek Lake Benthic Invertebrate Index (GLBiI), which is composed of three metrics, i.e., the number of taxa, the Simpson’s diversity index in the profundal and sublittoral zones, and the relative contribution of Chironomidae (%) family in the profundal zone [32]; and Hellenic Lake Littoral Benthic Invertebrate Assessment System (HeLLBI), which is composed of three metrics, i.e., relative abundance of Odonata (% of abundance classes), average score per taxon, and Simpson Diversity Index [33]. The GLBiI responds to eutrophication, while HeLLBI to hydromorphological alterations and eutrophication.
Specifically, phytoplankton samplings took place during the warm period, in the pelagic zone, covering the euphotic zone water column (2.5 × Secchi disk depth) using a Nansen type sampler. Benthic macroinvertebrates from the profundal/sublittoral zones were sampled in spring and autumn, with the use of a grab [32]. Littoral benthic macroinvertebrates were sampled in spring using the three-minute kick/sweep method with standard hand net (500 μm mesh size) [33]. Regarding macrophytes, sampling took place during the vegetative period, and the belt transect-mapping method was applied [30]. Finally, the dataset used comprised 46 lake years with phytoplankton EQRs, 11 with macrophyte, 6 with profundal/sublittoral zoobenthos, and 7 with littoral zoobenthos EQRs.
For each assessment method, as requested by the WFD, the Ecological Quality Ratio (EQR) was calculated as each index value was divided by the reference condition values, ranging from 0 (bad) to 1 (high). The reference values for GLBiI were estimated by the hindcasting procedure for each lake and do not correspond to actual reference sites’ values, while for the rest of the indices, the reference values are presented in Appendix A (Table A1). The class boundaries were evenly spaced for all BQEs (i.e., high > 0.8, good: (0.8–0.6), moderate: (0.6–0.4), poor: (0.4–0.2), bad ≤ 0.2), and thus, the comparison among them was feasible. EQR values for each BQE were available at different time intervals (Table 2). The values for the application of the OOAO principle were estimated for the two river basin management cycles (2012–2015 and 2016–2019). Analyses were applied using different combinations for integrating multiple BQEs at waterbody level. These combinations were:
(a)
OOAO: the lowest EQR of BQEs was attributed for the whole waterbody;
(b)
average: the arithmetic average of the EQRs for all BQEs was calculated and rounded to the nearest class; and
(c)
median: the median of EQRs for all BQEs was calculated and rounded to the nearest class.
It is worth noting that we performed all the statistical analyses based on the above-mentioned BQEs since in Greece, for these monitoring periods, data for lakes are limited only to those. Moreover, we did not use diatoms since the national index based on diatoms is under development and not available yet.

2.3. The Proposed Fuzzy Regression Model

The data of the fuzzy regression can be either fuzzy or crisp. Usually, the data are rather crisp numbers, and thus, the uncertainty arises from the adopted fuzzy model, that is, from the fuzzy coefficients. In this work, we deal only with crisp data. Fuzzy regression analysis gives a fuzzy functional relationship between the dependent and independent variables [42]. In contrast to the statistical regression, fuzzy regression analysis has no error term, while the uncertainty is incorporated in the model by using fuzzy numbers [39,43].
From a computational point of view, according to the Tanaka [44] approach, the problem of fuzzy linear regression is finally formulated as a constrained optimization problem. In the case that symmetrical fuzzy triangular numbers are used, the problem is transformed into a linear programming problem [45,46]. The main fuzzy concept that is used in fuzzy linear regression is the fuzzy number. The general description of fuzzy numbers can be found in Klir and Yuan [47]. Usually, fuzzy triangular symmetrical numbers are used. Since in the examined problems there is a lack of data, the authors selected the simplest membership function, which is the fuzzy symmetrical triangular number, to avoid any overtraining. Hence, here, A ˜ means a fuzzy triangular symmetrical number (Figure 1), which has the membership function (μA) presented below:
μ A a = 1 a 𝘢 i c i ,   i f   𝘢 i c i a 𝘢 i + c i 0
Many times, the fuzzy triangular symmetrical numbers are denoted as A ˜ = 𝘢 , c . The membership function expresses the degree according to which a member of the real axis (in general, the general set) belongs to the examined fuzzy number. The fuzzy numbers can be described as “approximately (or almost) ai” (e.g., almost one). The term ai represents the central value and the term c the semi-width (Figure 2).
The fuzzy linear regression model proposed by Tanaka [44] and Tanaka et al. [48] has the following form in this study:
Y ˜ j = A ˜ 0 + A ˜ 1 x j   with   j = 1 , ,   m  
where xj is the independent variable (here, the phytoplankton index), m is the number of data, and Y ˜ j is the fuzzy predicted value of the dependent variable (OOAO principle) considering the jth data. In other words, the fuzzy linear regression uses fuzzy numbers as coefficients instead of crisp numbers.
In addition, as it will be explained below, the fuzzy estimation of the dependent variable Y ˜ j must contain the observed data yj (Figure 3). The observed data are crisp numbers (OOAO principle), whilst the fuzzy regression provides a fuzzy estimation of the OOAO.
For instance, if the regression concludes to
Y ˜ j = 0.0457 , 0.0457 + 1 , 0 x j ,
this means that the main parameter that influences the OOAO is the phytoplankton index (xj). However, an uncertainty arises from the constant term (−0.0457, 0.0457). In contrast, if the model concludes to:
Y ˜ j = 0.4237 ,   0.0137 + 0 , 0 x j
this means that the phytoplankton index has no influence in the evaluation of the OOAO. Then the outcome, that is, the OOAO (fuzzy) estimation, will be a fuzzy tube parallel to xx’ axis, and the estimation will be between   a 0 c 0 , a 0 + c 0 , = 0.41 , 0.4374 . The fuzzy constant term expresses the determinant influence of the other indices.
The h-cut set of the fuzzy number A (with 0 < h ≤ 1) is defined as follows (Figure 3):
A h = α μ A ( a ) h ,   α
where is the real number axis.
Notice that the h-cut set is a crisp set determined from the fuzzy set according to a selected value of the membership function and, alternatively, a fuzzy set can be practically derived from a significant number of h-cut sets. In case of h = 0, the above definition of Equation (5) can be modified without the equality in order to describe the zero-cut [39].
The constraints express the concept of inclusion of the observed data (which are crisp numbers) within the produced fuzzy band. In general, the inclusion of a fuzzy set A into the fuzzy set B with the associated degree is defined as follows:
A h B h
In our case, the set A is a crisp set whilst the set B is a fuzzy set. Therefore, for each point of data, the observed dependent variable yj (OOAO) must be included into the produced h-cut of the fuzzy number Y ˜ j (blue line in Figure 4).
According to the presented methodology, the selection of the used h-level is dependent on the quality of data. Another interesting point is that even if the initial definition of the inclusion property incorporates all the h-cuts, the definition of inclusion that is used in Equation (6) is based on a unique (pre-selected) level h instead of all the values.
In this study, since the data are crisp (for each individual data), the set A is only a crisp value (a point of data that must be included in the produced fuzzy band), and the fuzzy set B is a fuzzy triangular number. Hence, the inclusion constraints in our problem become [40,49]:
a 1 x j + a 0 ( 1 h ) c 1 x j + c 0 y j a 1 x j + a 0 + ( 1 h ) c 1 x j + c 0 , j = 1 , , m
where a1, a0 are the centers of the coefficients that correspond to the independent variable (phytoplankton index) and the constant term correspondingly. The terms c1, c0 indicate the semi-widths of the coefficients that correspond to the independent variable and the constant term correspondingly. More analytically, the term a 1 x j + a 0 ( 1 h ) c 1 x j + c 0 expresses the left boundary of the h-cut of the fuzzy number Y ˜ j , whilst the term a 1 x j + a 0 + ( 1 h ) c 1 x j + c 0 expresses the right boundary of the h-cut of the fuzzy number Y ˜ j .
In case that there is not enough data, as in our case, the minimization of the achieved fuzzy band is proposed. Indeed, it is obvious that if the produced fuzzy band has a large magnitude, then it will contain all data. However, no functional fuzzy relation would be produced.
Even if, initially, Tanaka [45] proposed the minimization of the semi-width of the fuzzy coefficients, finally, he proposed the sum of the semi-widths (fuzzy spreads) for the produced dependent variable for all the data:
j = 1 m y R - y L = 2 j = 1 m c 1 x j + c 0 = 2 m c 0 + j = 1 m c 1 x j
where y R ,   y L are the right-hand boundary and the left-hand boundary, respectively, of the fuzzy set, which are the boundaries of the zero-cut.
Therefore, in condition of fuzzy triangular numbers as coefficients and by using the mentioned objective function, the problem of fuzzy linear regression is concluded to a linear programming problem [39,43]. The linear programming can be solved by using a plenty of commercial packages. In addition, the linear programming has very good mathematical properties, for instance, a local minimum solution is global and hence, in contrast with the no linear optimization, in this case, the LP provides a global optimum solution.
In the examined case (by using fuzzy triangular symmetrical numbers), the equation concludes to:
J = m i n m c 0 + j = 1 m c 1 x i j
subject to restrictions:
y j a 1 x j + a 0 ( 1 h ) c 1 x j + c 0
y j a 1 x j + a 0 + ( 1 h ) c 1 x j + c 0
c1, c0, ≥ 0, where j = 0, 1,…, m
The above mathematical functions, a 1 x j + a 0 ( 1 h ) c 1 x j + c 0 ,   a 1 x j + a 0 + ( 1 h ) c 1 x j + c 0 present the lower and the upper boundaries, respectively, of the corresponding h-cut of Y ˜ j . As y j , the jth observed data is meant, considering the dependent variable, which in this application is a crisp number.
It is well known that the solution for any other cut-level, h’ can be obtained from the optimal h-cut-level solution as follows [50,51]:
A h = a h ,   1 h 1 h c h
the term h means the initial h-cut and the term h’ means the new level h’. From Equation (10), it is obvious that the new level h’ affects only the semi-widths. Hence, in case that there is not enough data, after the solution that is achieved from h = 0, the optimal solution for any other level h’ can be achieved. The total fuzzy band will be:
J c h = 1 1 h J c 0
The term c h , c 0 are the matrices which represent the semi-width of the solution with h’ and for h = 0, respectively. For instance, the new matrix of the semi-width for h = 0.5 will be equal to:
c h = 0.5 = 2 c h = 0

2.4. Tested Scenarios and Basic Interpretation

The testing that we included in our research was split in two distinct scenarios:
Scenario a: the fuzzy regression was applied between phytoplankton, as the most sensitive and more frequently monitored BQE, and the OOAO values. The values for OOAO were generated using the unique yearly EQR values of phytoplankton, but the same values were kept for the indices of macrophytes and benthic macroinvertebrates for all years within each monitoring period (2012–2015 and 2016–2019). That is, the EQR values of these indices were extended beyond their actual monitored year (Table 2). According to the WFD, the monitoring frequency of macrophytes and benthic macroinvertebrates for operational monitoring is once every three years.
Scenario b: the fuzzy regression took place between phytoplankton yearly EQR values and the OOAO values generated using the yearly EQR values of phytoplankton and only the values of the indices of macrophytes and benthic invertebrates corresponding to the year they were actually monitored (Table 2).
The presented fuzzy regression model is applied between the yearly values of the independent (phytoplankton index) variable and the quality expressed by the dependent variable (OOAO). Therefore, when the coefficients of the phytoplankton index have low values, that is, small central value (ai) and semi-width (ci), then it is obvious that the other quality indices affect the OOAO. On the other hand, if the coefficients of the phytoplankton index have a high value, which is large central value and significant semi-width, then the phytoplankton index significantly affects the OOAO. To work properly, there has to be originally a wide range of EQR values from all indices included in the regression. When a lake has quality values only within an EQR class range, it “cripples” the lake, not allowing it to go beyond these values. When both EQR values in both axes are identical, the Phytoplankton index “commands” the OOAO, and in that case, the observation value is placed closest to the central fuzzy value. The boundaries designate the maximum acceptable value for the OOAO principle acknowledging every value in the dataset. That means that the final acceptable assessment value can be within the range set by the central value plus or minus the semi-width (ci, the second coefficient). This coefficient is the uncertainty created by the existing dataset at hand each time. When a boundary value falls over an existing observed value, this means that another BQE designated the range of acceptable values. It can be easily understood that the higher the x (phytoplankton) coefficient value, the more important is phytoplankton for the OOAO principle. Similarly, the larger the second coefficient is, the wider the acceptable ranges. It should also be noted that for every new added value in the dataset (i.e., values for 2020, 2021, etc.), the regression can produce a slightly differentiated equation, meaning new semi-widths and thus new acceptable EQR ranges.

3. Results

3.1. BQEs Comparison and Status Classification Approaches

The ecological status assessed for the case-studied lakes differ up to two quality scales (Table 3). The macrophyte-based index appears to be the most lenient, while the index based on macroinvertebrates of the sublittoral and profundal zones appears as the strictest. On the other hand, the index based on the macroinvertebrates of the littoral zone remains more “austere” than phytoplankton, except in the cases of Vegoritis and Lysimachia (Table 3). The results for both monitoring periods using the OOAO principle end in 16.7% of the cases in poor and good ecological quality and 66.7% in moderate, while when using the median and average approach, the results were similar, classifying the cases at 50% moderate and 50% good ecological quality. As for the final classification, the OOAO principle, selecting the worst case, showed the poorest ecological status, downgrading up to two scales the status assessment, whereas the median and average indices’ EQRs were more kind in total, reflecting the quality distance among BQEs mentioned above.

3.2. Implementation of the Proposed Methodology

3.2.1. Fuzzy Regression Scenario A

In the cases of the Vegoritis and Volvi lakes, the estimated OOAO values form a tube parallel to the xx’ axis (Figure 5a,b). This means that the phytoplankton index has no significant influence on the water quality classification for the examined lakes. To be precise, phytoplankton has a negligible coefficient (0.0025) for Lake Vegoritis, while for Lake Volvi, phytoplankton has no effect in shaping the OOAO values (Figure 5a,b). In this case, it is obvious that all the other indices rather than the phytoplankton index influence the OOAO. In contrast, the OOAO equation for Lake Yliki is proportional to the phytoplankton index, while the semi-width indirectly indicates the influence of the other indices (Figure 5c) encompassing the generated uncertainty (0.3123), which is also dependent on phytoplankton EQR.
Accordingly, Figure 6 depicts the fuzzy regressions for shallow natural lakes. A simple line is produced for Lake Doirani (Figure 6a) since the phytoplankton index has the lowest value among the indices, and hence, it is identical with the OOAO values, demonstrating no uncertainty. In Ozeros and Lysimachia lakes, the produced fuzzy curves represent a tube almost parallel to the xx’ axis (Figure 6b,c). In Lake Lysimachia (Figure 6b), a fuzzy coefficient (0.4038) is produced from the EQR values generated from all other indices than the phytoplankton one. Hence, the final ecological status would be within the range of 0.4038 ± 0.0138. The OOAO values for Lake Ozeros are slightly dependent on phytoplankton (0.0038) (Figure 6c).

3.2.2. Fuzzy Regression Scenario B

It seems that performing a fuzzy regression with phytoplankton EQR values (one for each year of the monitoring period) and two values for all other indices (one per four years) shapes the OOAO in favor of phytoplankton due to reduced yearly EQR comparisons. So, regardless of the lake type (shallow or deep), the independent variable (phytoplankton index) becomes a driving factor for the regression (Figure 7 and Figure 8). In all lakes, the OOAO values are almost equivalent to that of the phytoplankton index, except in the case of Lake Ozeros (Figure 8). In that case, the phytoplankton coefficient is also high (0.7865), and the uncertainty is phytoplankton related, in contradiction with Lysimachia, Yliki, and Volvi lakes, where the uncertainty is an independent term.

4. Discussion

Assigning nature to boxes is difficult, and therefore, it is a challenge to fulfil the requirements of the WFD to define an overall ecological status, considering the ecosystem as a whole and using multiple BQEs in different waterbody types. The OOAO principle, where the status is based on the lowest score of any of the BQEs, is a precautionary rule that can lead to the strictest classification of waterbodies’ ecological status [52]. In this study, we attempted to consider combinations for multiple BQEs to assess ecological status in two different lake types. Moreover, we performed a fuzzy method to overcome the obstacle of the blind application of the OOAO approach, taking into account the mandatory monitoring periods according to the WFD.
As for the mere EQR values comparison, the macrophyte-based index appears as the less strict index, with an exception in the case of Lysimachia, possibly because it focuses on the euphotic littoral zone. It is considered that the abundance and the composition of plant species very often have a quite wide ecological scale [19,20]. Additionally, which BQE reflects better the lake quality is a matter of the lake itself, since all differ according to their morphometry and hydrology, and they are usually subjected to different pressures. Even within a lake, multiple BQEs can respond differently in various lakes’ parts, seasons, and time scales [53,54]. According to previous studies [21,52,55] regarding the different approaches to evaluate the overall status, the OOAO principle was the most conservative, downgrading the ecological status of lakes. The OOAO principle is also considered as an inadequate tool for addressing progress or deterioration in status [56], as it masks the real changes that WFD implementation can achieve, especially in cases where PoMs have been implemented with delay. The suggested fuzzy method in all scenarios retains the strictness of the OOAO principle but in a way that the main BQE (i.e., phytoplankton) can be boosted and the character of each lake showcased. This approach comes in line with the most voted choices in a query posed to all experts and attendees during EU WFD 2017 e-conference [11]. Moreover, the fuzzy regression can keep learning the waterbody as the dataset expands. In Scenario A, using repeatedly the same value for the indices of macrophytes and benthic invertebrates for four consecutive years, is a bias against the main BQE (i.e., phytoplankton) for lakes status classification as the regression frames any picture that this BQE could give. This was depicted by the fuzzy regression that generated tubes parallel to xx’ axis and the equations that lessened the importance of phytoplankton. The only exception was Lake Yliki because EQR values were based on macroinvertebrate and macrophyte indices, which were not available during the second monitoring period (2016–2019), hence affecting the OOAO and consequently the fuzzy regression. This bias is imposed by the WFD recommendations requiring different monitoring frequencies for each BQE responding to different life spans and community variability. Regarding Scenario B (using only the value when the BQEs were monitored), one would claim that it is biased in favor of phytoplankton, but we deemed that this approach allows the main BQE (i.e., phytoplankton) to express the yearly quality variation, respecting the OOAO principle when other BQEs were also monitored. In any of the scenarios, this fuzzy method enables the experts to assess if the most often monitored BQE (i.e., phytoplankton) plays a crucial role in defining the OOAO values. This enables the assessment of the waterbody with the uncertainty expressed by the boundaries, which work as a safety net. Additionally, it guides the experts to check which of the other BQEs characterize better the qualitative status of the lake.
As for the progress reporting of individual BQEs, according to Carvahlo et al. [11], this requires a larger (lakes and monitoring years) dataset to apply a statistically solid fuzzy regression. However, to evaluate the progress as expressed by indices is an ambitious task, given the result is practically often a circumstance of the relationship among the nature of administrative governance, the nature of PoMs taken [57], and the minimum of WFD obligations (along with the expected lag from the waterbody). Thus, the framework for overall ecological status classification should in fact include a form of elasticity. It is expected that different governance approaches for the same nature of measures could end in diverse results due to design and implementation progress. Tsakiris [58] highlighted that the first RBMPs showed a significant progress on the fulfilment of environmental objectives, while the second reporting period did not support the same rate in progress [59,60]. Following the OOAO principle as safer, only 20% of the EU waterbodies showed improvement in their status [61] despite costly measures [62]. Instead, there was a call for more stringent nutrient thresholds [63,64], signifying the need for new approaches for the sustainability of WFD 20+. Moreover, Greece and especially its lake basins host mainly agroecosystems, where diffuse pollution cannot easily be tackled, and green infrastructure measures need more time to pay off in quality status terms [8,57]. This issue is aggravated under the Mediterranean climate prism and some specific hydromorphological characteristics (i.e., shore alterations, high water retention time). Moreover, a six-year period for RBMP cycles renewal can be considered as adequate for reporting progress in ecological point of view, but such a period can hinder trends and the effectiveness of policy measures [56] or include the possible benefits from action plans in a management point of view. The next step that goes beyond the OOAO principle, and is nowadays under consultation [56], is the use of new quality indicators which are focused on implementation shifts (i.e., pressure reduction, PoMs application) [56]. Novel reporting protocols can perhaps improve the digitalization and administrative simplification but add extra burden on MSs to follow the changes. It is unclear, though, if this change can positively impact the OOAO principle in any means.
In this study, we suggest the development of an overall assessment that will be based on the precautionary OOAO principle through the fuzzy OOAO approach, which we consider that it enables a better implementation of the WFD and allows for a case-specific evaluation, including the uncertainty in classification as an asset. We argue that such an approach offers a deeper system understanding through self-learning processes based on the existing datasets, while it also recognizes the different character of the lakes and adds a special gravity to the most sensitive and most frequently monitored BQE. However, we believe that all BQEs should be monitored, as freshwater ecosystems are highly complex systems [65], and an efficient management should depend on multiple components [66]. Additionally, all BQEs samples should be taken in a timing able to support a more proper application of OOAO (as it was treated in Scenario B) under the same environmental conditions (i.e., meteorological, hydrological) to establish more stable stressor-response relationships.

5. Conclusions

  • The inclusion of a fuzzy regression among the frequently monitored BQE (phytoplankton) and the outcome of OOAO (determined by the comparison of four BQE indices) application in lakes encompasses the uncertainty and the possibility to broaden the acceptable final EQR based on the character and status of each lake;
  • The fuzzy OOAO is an approach that seems to allow a better understanding of the WFD implementation and case-specific evaluation, including the uncertainty in classification as an asset;
  • It offers a deeper understanding through self-learning processes based on the existing datasets;
  • As for the progress reporting of individual BQEs, this requires a more complete dataset to apply a statistically solid fuzzy regression.

Author Contributions

Conceptualization, I.K. and D.L.; methodology, M.S., D.L., and C.N.; software, M.S.; formal analysis, M.S., D.L., and C.N.; resources, V.T.; data curation, D.L. and C.N.; writing—original draft preparation, D.L., C.N., M.S., and I.K.; writing—review and editing, D.L., C.N., M.S., I.K., M.L., D.B., and V.T.; visualization, M.S. and C.N.; supervision, I.K. and M.L.; project administration, D.B. All authors have read and agreed to the published version of the manuscript.

Funding

Part of the data used in this research come from Act MIS 5001204, financed by the European Union Cohesion Fund (Partnership Agreement 2014–2020) and from Acts MIS 371010, 371138, 371140, 371145, financed by the European Regional Development Fund (National Strategic Reference Framework 2007–2013). Publication costs were covered by No 98561 private grant managed by Special Account for Research Funds of Aristotle University of Thessaloniki for consulting services towards the Special Secretariat for Waters-Greece.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data that support reported results for phytoplankton, macrophytes, and littoral benthos were used under license from The Goulandris Natural History Museum, Greek Biotope/Wetland Centre (EKBY). They are available upon reasonable request from V.T. and with permission of The Goulandris Natural History Museum, Greek Biotope/Wetland Centre (EKBY). Regarding the data for sublittoral and profundal benthos, they are available upon reasonable request from C.N.

Acknowledgments

Part of the present study was conducted in the frame of the Greek National Water Monitoring Network, according to the JMD 140384/2011, implemented by The Goulandris Natural History Museum, Greek Biotope/Wetland Centre (EKBY). The network is supervised by the General Directorate for Waters of the Ministry of Environment and Energy. EKBY’s personnel conducted samplings and contributed to sample analysis. Part of phytoplankton analysis was commissioned to AUTH.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Reference values of phytoplankton, macrophytes, and littoral invertebrates indices.
Table A1. Reference values of phytoplankton, macrophytes, and littoral invertebrates indices.
IndexMetricsLake Type GR-DNLLake Type GR-SNL
HeLPhyTotal Phytoplankton Biovolume (mm3 L−1)1.290.74
Cyanobacteria Biovolume (mm3 L−1)0.010.01
modNygaard Index1.031.11
Chlorophyll a (μg L−1)1.563.59
HeLMTIHelm7.147.14
Cmax (m)12.26.1
HeLLBIASPT5.475.47
Odonata (% AC)16.6716.67
Simpson0.800.80

References

  1. European Community. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 establishing a framework for com-munity action in the field of water policy. Eur. Dir. 2000, L327, 1–72. [Google Scholar]
  2. Moe, S.; Couture, R.-M.; Haande, S.; Lyche Solheim, A.; Jackson-Blake, L. Predicting Lake Quality for the Next Generation: Impacts of Catchment Management and Climatic Factors in a Probabilistic Model Framework. Water 2019, 11, 1767. [Google Scholar] [CrossRef] [Green Version]
  3. European Community. Directive 2008/105/EC of the European Parliament and of the Council of 16 December 2008 on environmental quality standards in the field of water policy, amending and subsequently repealing Council Directives 82/176/EEC, 83/513/EEC, 84/156/EEC, 84/491/EEC, 86/280/EEC and amending Directive 2000/60/EC of the European Parliament and of the Council. Eur. Dir. 2008, L348, 84–97. [Google Scholar]
  4. European Community. Directive 2006/118/EC of the European Parliament and of the Council of 12 December 2006 on the protection of groundwater against pollution and deterioration. Eur. Dir. 2006, L372, 19–31. [Google Scholar]
  5. European Community. Directive 2007/60/EC of the European Parliament and of the Council of 23 October 2007 on the assessment and management of flood risks (Text with EEA relevance). Eur. Dir. 2007, L288, 27–34. [Google Scholar]
  6. Kristensen, P.; Whalley, C.; Zal, F.N.N.; Christiansen, T. European Waters Assessment of Status and Pressures 2018; European Environment Agency: Copenhagen, Denmark, 2018; p. 85. [Google Scholar]
  7. European Commission. Fitness check of the Water Framework Directive, Groundwater Directive, Environmental Quality Standards Directive and Floods Directive; European Commission: Brussels, Belgium, 2019. [Google Scholar]
  8. Kagalou, I.; Latinopoulos, D. Filling the Gap between Ecosystem Services Concept and River Basin Management Plans: The Case of Greece in WFD 20+. Sustainability 2020, 12, 7710. [Google Scholar] [CrossRef]
  9. Prato, S.; La Valle, P.; De Luca, E.; Lattanzi, L.; Migliore, G.; Morgana, J.G.; Munari, C.; Nicoletti, L.; Izzo, G.; Mistri, M. The “one-out, all-out” principle entails the risk of imposing unnecessary restoration costs: A study case in two Mediterranean coastal lakes. Mar. Pollut. Bull. 2014, 80, 30–40. [Google Scholar] [CrossRef]
  10. Moe, S.J.; Lyche Solheim, A.; Soszka, H.; Gołub, M.; Hutorowicz, A.; Kolada, A.; Picińska-Fałtynowicz, J.; Białokoz, W. Integrated assessment of ecological status and misclassification of lakes: The role of uncertainty and index combination rules. Ecol. Indic. 2015, 48, 605–615. [Google Scholar] [CrossRef] [Green Version]
  11. Carvalho, L.; Mackay, E.B.; Cardoso, A.C.; Baattrup-Pedersen, A.; Birk, S.; Blackstock, K.L.; Borics, G.; Borja, A.; Feld, C.K.; Ferreira, M.T.; et al. Protecting and restoring Europe’s waters: An analysis of the future development needs of the Water Framework Directive. Sci. Total Environ. 2019, 658, 1228–1238. [Google Scholar] [CrossRef]
  12. Birk, S.; Chapman, D.; Carvalho, L.; Spears, B.M.; Andersen, H.E.; Argillier, C.; Auer, S.; Baattrup-Pedersen, A.; Banin, L.; Beklioglu, M.; et al. Impacts of multiple stressors on freshwater biota across spatial scales and ecosystems. Nat. Ecol. Evol. 2020, 4, 1060–1068. [Google Scholar] [CrossRef]
  13. Hering, D.; Borja, A.; Carstensen, J.; Carvalho, L.; Elliott, M.; Feld, C.K.; Heiskanen, A.S.; Johnson, R.K.; Moe, J.; Pont, D.; et al. The European Water Framework Directive at the age of 10: A critical review of the achievements with recommendations for the future. Sci. Total Environ. 2010, 408, 4007–4019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Reid, A.J.; Carlson, A.K.; Creed, I.F.; Eliason, E.J.; Gell, P.A.; Johnson, P.T.J.; Kidd, K.A.; MacCormack, T.J.; Olden, J.D.; Ormerod, S.J.; et al. Emerging threats and persistent conservation challenges for freshwater biodiversity. Biol. Rev. 2019, 94, 849–873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Carvalho, L.; Poikane, S.; Lyche Solheim, A.; 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. Hydrobiologia 2013, 704, 127–140. [Google Scholar] [CrossRef] [Green Version]
  16. Clarke, R.T.; Hering, D. Errors and uncertainty in bioassessment methods—major results and conclusions from the STAR project and their application using STARBUGS. Hydrobiologia 2006, 566, 433–439. [Google Scholar] [CrossRef]
  17. Kelly, M.; Bennion, H.; Burgess, A.; Ellis, J.; Juggins, S.; Guthrie, R.; Jamieson, J.; Adriaenssens, V.; Yallop, M. Uncertainty in ecological status assessments of lakes and rivers using diatoms. Hydrobiologia 2009, 633, 5–15. [Google Scholar] [CrossRef]
  18. Thackeray, S.J.; Nõges, P.; Dunbar, M.J.; Dudley, B.J.; Skjelbred, B.; Morabito, G.; Carvalho, L.; Phillips, G.; Mischke, U.; Catalan, J.; et al. Quantifying uncertainties in biologically-based water quality assessment: A pan-European analysis of lake phytoplankton community metrics. Ecol. Indic. 2013, 29, 34–47. [Google Scholar] [CrossRef] [Green Version]
  19. Kolada, A. The use of aquatic vegetation in lake assessment: Testing the sensitivity of macrophyte metrics to anthropogenic pressures and water quality. Hydrobiologia 2010, 656, 133–147. [Google Scholar] [CrossRef]
  20. Kuczyńska-Kippen, N. The use of bdelloids in reference to rotifer biocoenotic indices as an indicator of the ecological state of small field water bodies: The effect of macrophytes, shading and trophic state of water. Ecol. Indic. 2018, 89, 576–583. [Google Scholar] [CrossRef]
  21. Borja, A.; Rodriguez, J.G. Problems associated with the ‘one-out, all-out’ principle, when using multiple ecosystem components in assessing the ecological status of marine waters. Mar. Pollut. Bull. 2010, 60, 1143–1146. [Google Scholar] [CrossRef]
  22. Moss, B. The Water Framework Directive: Total environment or political compromise? Sci. Total Environ. 2008, 400, 32–41. [Google Scholar] [CrossRef]
  23. Borja, A.; Ranasinghe, A.; Weisberg, S.B. Assessing ecological integrity in marine waters, using multiple indices and ecosystem components: Challenges for the future. Mar. Pollut. Bull. 2009, 59, 1–4. [Google Scholar] [CrossRef]
  24. Hering, D.; Johnson, R.K.; Kramm, S.; Schmutz, S.; Szoszkiewicz, K.; Verdonschot, P.F.M. Assessment of European streams with diatoms, macrophytes, macroinvertebrates and fish: A comparative metric-based analysis of organism response to stress. Freshw. Biol. 2006, 51, 1757–1785. [Google Scholar] [CrossRef]
  25. McFarland, B.; Carse, F.; Sandin, L. Littoral macroinvertebrates as indicators of lake acidification within the UK. Aquat. Conserv. Mar. Freshw. Ecosyst. 2010, 20, S105–S116. [Google Scholar] [CrossRef]
  26. Urbanič, G. A Littoral Fauna Index for assessing the impact of lakeshore alterations in Alpine lakes. Ecohydrology 2014, 7, 703–716. [Google Scholar] [CrossRef]
  27. Poikane, S.; Johnson, R.K.; Sandin, L.; Schartau, A.K.; Solimini, A.G.; Urbanic, G.; Arbaciauskas, 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]
  28. Poikane, S.; Ritterbusch, D.; Argillier, C.; Bialokoz, W.; Blabolil, P.; Breine, J.; Jaarsma, N.G.; Krause, T.; Kubecka, J.; Lauridsen, T.L.; et al. Response of fish communities to multiple pressures: Development of a total anthropogenic pressure intensity index. Sci. Total Environ. 2017, 586, 502–511. [Google Scholar] [CrossRef]
  29. Tsiaousi, V.; Zervas, D.; Tsiripidis, I. Report on the Development of the National Method for the Assessment of the Ecological Status of Natural Lakes in Greece, Using the Biological Quality Element “Phytoplankton”, 1st ed.; Greek Biotope/Wetland Centre and Special Secretariat for Waters, Ministry of Environment: Thermi, Greece, 2017; p. 17.
  30. Zervas, D.; Tsiaoussi, V.; Tsiripidis, I. HeLM: A macrophyte-based method for monitoring and assessment of Greek lakes. Environ. Monit. Assess. 2018, 190, 326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Petriki, O.; Lazaridou, M.; Bobori, D.C. A fish-based index for the assessment of the ecological quality of temperate lakes. Ecol. Indic. 2017, 78, 556–565. [Google Scholar] [CrossRef]
  32. Ntislidou, C.; Lazaridou, M.; Tsiaoussi, V.; Bobori, D.C. A new multimetric macroinvertebrate index for the ecological assessment of Mediterranean lakes. Ecol. Indic. 2018, 93, 1020–1033. [Google Scholar] [CrossRef]
  33. Mavromati, E.K.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; p. 11.
  34. Salski, A. Ecological applications of fuzzy logic. In Ecological Informatics; Springer: Berlin/Heidelberg, Germany, 2006; pp. 3–14. [Google Scholar]
  35. Prato, T. A fuzzy logic approach for evaluating ecosystem sustainability. Ecol. Model. 2005, 187, 361–368. [Google Scholar] [CrossRef]
  36. Pouw, F.A.; Kwiatkowska, M. An overview of fuzzy-logic based approaches to ecology: Addressing uncertainty. In Proceedings of the 2013 Joint IFSA World Congress and NAFIPS Annual Meeting (IFSA/NAFIPS), Edmonton, AB, Canada, 24–28 June 2013; pp. 540–545. [Google Scholar]
  37. Spiliotis, M.; Panagiotou, L.; Kagalou, I.; Latinopoulos, D. A Fuzzified Multicriteria Outranking Method for Water Framework Directive Implementation in a Heavily Modified Urban Lake (Pamvotis, Greece). Water Resour. Manag. 2020, 34, 4491–4510. [Google Scholar] [CrossRef]
  38. Ganoulis, J. Engineering Risk Analysis of Water Pollution Probabilities and Fuzzy Sets; Verlagsgesellschaft: Weinheim, Germany, 1994; p. 306. [Google Scholar]
  39. Kitsikoudis, V.; Spiliotis, M.; Hrissanthou, V. Fuzzy Regression Analysis for Sediment Incipient Motion under Turbulent Flow Conditions. Environ. Process. 2016, 3, 663–679. [Google Scholar] [CrossRef]
  40. Spiliotis, M.; Hrissanthou, V. Fuzzy and crisp regression analysis between sediment transport rates and stream discharge in the case of two basins in northeastern Greece. Conv. Fuzzy Regres. Theory Eng. Appl. 2018, 1, 49. [Google Scholar]
  41. Spiliotis, M.; Angelidis, P.; Papadopoulos, B. Assessment of annual hydrological drought based on fuzzy estimators. In Proceedings of the 4th IAHR Europe Congress, “Sustainable hydraulics in the era of global change”, Liege, Belgium, 27–29 July 2016. [Google Scholar]
  42. Papadopoulos, B.K.; Sirpi, M.A. Similarities in fuzzy regression models. J. Optim. Theory Appl. 1999, 102, 373–383. [Google Scholar] [CrossRef]
  43. Spiliotis, M.; Bellos, C. Flooding risk assessment in mountain rivers. In Proceedings of the Water resources management in a changing World: Challenges and opportunities, Instabul, Turkey, 10–13 June 2015. [Google Scholar]
  44. Tanaka, H. Fuzzy data analysis by bossibilistic linear models. Fuzzy Sets Syst. 1987, 24, 262–375. [Google Scholar] [CrossRef]
  45. Tanaka, H.; Uejima, S.; Asai, K. Linear regression analysis with fuzzy model. IEEE Trans. Syst. Man. Cybern. 1982, 12, 903–907. [Google Scholar]
  46. Tsakiris, G.; Tigkas, D.; Spiliotis, M. Assessment of interconnection between two adjacent watersheds using deterministic and fuzzy approaches. Eur. Water 2006, 15, 15–22. [Google Scholar]
  47. Klir, G.; Yuan, B. Fuzzy Sets and Fuzzy Logic: Theory and Applications; Prentice Ыall: Hoboken, NJ, USA, 1995; p. 574. [Google Scholar]
  48. Tanaka, H.; Hayashi, I.; Watada, J. Possibilistic linear regression analysis for fuzzy data. Eur. J. Oper. Res. 1989, 40, 389–396. [Google Scholar] [CrossRef]
  49. Tzimopoulos, C.; Papadopoulos, K.; Papadopoulos, B. Fuzzy regression with applications in hydrology. Int. J. Eng. Innov. Technol. 2016, 5, 69–75. [Google Scholar]
  50. Papadopoulos, B.K.; Sirpi, M.A. Similarities and distances in fuzzy regression modeling. Soft Comput. 2004, 8, 556–561. [Google Scholar] [CrossRef]
  51. Profillidis, V.A.; Papadopoulos, B.K.; Botzoris, G.N. Similarities in fuzzy regression models and application on transportation. Fuzzy Econ. Rev. 1999, 4, 4. [Google Scholar] [CrossRef]
  52. Caroni, R.; Bund, W.; Clarke, R.T.; Johnson, R.K. Combination of multiple biological quality elements into waterbody assessment of surface waters. Hydrobiologia 2013, 704, 437–451. [Google Scholar] [CrossRef]
  53. 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. Hydrobiologia 2013, 704, 57–74. [Google Scholar] [CrossRef] [Green Version]
  54. Poikane, S.; Portielje, R.; Denys, L.; Elferts, D.; Kelly, M.; Kolada, A.; Maemets, H.; Phillips, G.; Sondergaard, 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] [PubMed]
  55. Alahuhta, J.; Vuori, K.-M.H.; Olin, M.; Rask, M.; Palomäki, A. Defining the ecological status of small forest lakes using multiple biological quality elements and palaeolimnological analysis. Fundam. Appl. Limnol. 2009, 175, 203–216. [Google Scholar] [CrossRef]
  56. European Commission. Draft Technical Report on Water Quality Indicators. Available online: https://circabc.europa.eu/ui/group/9ab5926d-bed4-4322-9aa7-9964bbe8312d/library/0d521a5a-2f8e-49b6-a011-204ceb144960/details (accessed on 25 March 2021).
  57. Wiering, M.; Boezeman, D.; Crabbé, A. The Water Framework Directive and Agricultural Diffuse Pollution: Fighting a Running Battle? Water 2020, 12, 1447. [Google Scholar] [CrossRef]
  58. Tsakiris, G. The Status of the European Waters in 2015: A Review. Environ. Process. 2015, 2, 543–557. [Google Scholar] [CrossRef] [Green Version]
  59. Lyche Solheim, A.; Globevnik, L.; Austnes, K.; Kristensen, P.; Moe, S.J.; Persson, J.; Phillips, G.; Poikane, S.; Van de Bund, W.; Birk, S. A new broad typology for rivers and lakes in Europe: Development and application for large-scale environmental assessments. Sci. Total Environ. 2019, 697, 134043. [Google Scholar] [CrossRef]
  60. European Commission. Commission Staff Working Document. Report on the Implementation of the Water Framework Directive (2000/60/EC) and the Floods Directive (2007/60/EC)—Second River Basin Management Plans and First Flood Risk Management Plans; Report from the Commission to the European Parliament and the Council; European Commission: Brussels, Belgium, 2019. [Google Scholar]
  61. Zacharias, I.; Liakou, P.; Biliani, I. A Review of the Status of Surface European Waters Twenty Years after WFD Introduction. Environ. Process. 2020, 7, 1023–1039. [Google Scholar] [CrossRef]
  62. Voulvoulis, N.; Arpon, K.D.; Giakoumis, T. The EU Water Framework Directive: From great expectations to problems with implementation. Sci. Total Environ. 2017, 575, 358–366. [Google Scholar] [CrossRef] [Green Version]
  63. Poikane, S.; Várbíró, G.; Kelly, M.G.; Birk, S.; Phillips, G. Estimating river nutrient concentrations consistent with good ecological condition: More stringent nutrient thresholds needed. Ecol. Indic. 2021, 121, 107017. [Google Scholar] [CrossRef]
  64. Poikane, S.; Kelly, M.G.; Salas Herrero, F.; Pitt, J.A.; Jarvie, H.P.; Claussen, U.; Leujak, W.; Lyche Solheim, A.; Teixeira, H.; Phillips, G. Nutrient criteria for surface waters under the European Water Framework Directive: Current state-of-the-art, challenges and future outlook. Sci. Total Environ. 2019, 695, 133888. [Google Scholar] [CrossRef] [PubMed]
  65. Moss, B. Climate change, nutrient pollution and the bargain of Dr Faustus. Freshw. Biol. 2010, 55, 175–187. [Google Scholar] [CrossRef]
  66. Kelly, M.G.; Birk, S.; Willby, N.J.; Denys, L.; Drakare, S.; Kahlert, M.; Karjalainen, S.M.; Marchetto, A.; Pitt, J.A.; Urbanic, 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]
Figure 1. Map of the studied lakes.
Figure 1. Map of the studied lakes.
Water 13 01776 g001
Figure 2. Fuzzy triangular symmetrical number “almost a”.
Figure 2. Fuzzy triangular symmetrical number “almost a”.
Water 13 01776 g002
Figure 3. Triangular fuzzy number, its h-cut (is a crisp set), and zero-cut (α, α+).
Figure 3. Triangular fuzzy number, its h-cut (is a crisp set), and zero-cut (α, α+).
Water 13 01776 g003
Figure 4. The concept of inclusion in the case of fuzzy linear regression.
Figure 4. The concept of inclusion in the case of fuzzy linear regression.
Water 13 01776 g004
Figure 5. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario A for three deep natural lakes: (a) Vegoritis, (b) Volvi, and (c) Yliki.
Figure 5. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario A for three deep natural lakes: (a) Vegoritis, (b) Volvi, and (c) Yliki.
Water 13 01776 g005
Figure 6. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario A for three shallow natural lakes: (a) Doirani, (b) Lysimachia, and (c) Ozeros.
Figure 6. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario A for three shallow natural lakes: (a) Doirani, (b) Lysimachia, and (c) Ozeros.
Water 13 01776 g006
Figure 7. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario B for three deep natural lakes: (a) Vegoritis, (b) Volvi, and (c) Yliki.
Figure 7. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario B for three deep natural lakes: (a) Vegoritis, (b) Volvi, and (c) Yliki.
Water 13 01776 g007
Figure 8. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario B for three shallow natural lakes: (a) Doirani, (b) Lysimachia, and (c) Ozeros.
Figure 8. Graphical representation of the fuzzy regression for the OOAO (“One Out–All Out”) principle (dimensionless) with respect to phytoplankton index (dimensionless) in case of Scenario B for three shallow natural lakes: (a) Doirani, (b) Lysimachia, and (c) Ozeros.
Water 13 01776 g008aWater 13 01776 g008b
Table 1. Limnological characteristics of the studied lakes. OL, oligotrophic; MT, mesotrophic; ET, eutrophic.
Table 1. Limnological characteristics of the studied lakes. OL, oligotrophic; MT, mesotrophic; ET, eutrophic.
LakeAltitude (m.a.s.l.)Mean Depth (m) aMaximum Depth (m) aLake Area (km2) aTrophic Status
Doirani1424.55.532.4 *ET
Lysimachia163.57.713.0ET
Ozeros223.86.110.4ET
Vegoritis51026.152.447.4MT-ET
Volvi3712.527.372.9ET
Yliki8020.938.521.6OL-MT
a Data available from the national monitoring program implemented by the Greek Biotope-Wetland Centre (EKBY), * approximately 44% of the lake area being within the territory of Greece.
Table 2. Frequency of monitoring for each biological quality element. P, phytoplankton; M, macrophytes; B-L, littoral zoobenthos; B-SP, sublittoral/profundal zoobenthos.
Table 2. Frequency of monitoring for each biological quality element. P, phytoplankton; M, macrophytes; B-L, littoral zoobenthos; B-SP, sublittoral/profundal zoobenthos.
Lake1st Monitoring Period2nd Monitoring Period
20122013201420152016201720182019
DoiraniPP, MP, B-SPPP, MP P
LysimachiaPPP, M, B-SPP, B-LPP, MP, B-LP
OzerosPPP, M, B-SPP, B-LPP, B-LPP
VegoritisPP, MP, B-SPPP, MPP, B-LP
VolviPP, MP, B-SPPP, MP, B-L P
YlikiPPP, M, B-SPP, B-LPPPP
Table 3. Ecological quality of six Greek lakes, based on multi-metric indices for different biological quality elements. The ecological status assessed by the “One Out–All Out” (OOAO), median, and average approaches. Quality color scale and Ecological Quality Ratio values are presented. P, phytoplankton [29]; M, macrophytes [30]; B-SP, sublittoral/profundal zoobenthos [32]; B-L, littoral zoobenthos [33].
Table 3. Ecological quality of six Greek lakes, based on multi-metric indices for different biological quality elements. The ecological status assessed by the “One Out–All Out” (OOAO), median, and average approaches. Quality color scale and Ecological Quality Ratio values are presented. P, phytoplankton [29]; M, macrophytes [30]; B-SP, sublittoral/profundal zoobenthos [32]; B-L, littoral zoobenthos [33].
Lake TypeLakeMonitoring PeriodPMB-SPB-LOOAOMedianAverage
DeepVegoritis1st (2012–2015)0.660.750.54 0.540.660.65
2nd (2016–2019)0.640.62 0.690.620.640.65
Volvi1st (2012–2015)0.450.700.41 0.410.450.52
2nd (2016–2019)0.460.71 0.440.440.460.53
Yliki1st (2012–2015)0.770.690.340.480.340.590.57
2nd (2016–2019)0.75 0.750.750.75
ShallowDoirani1st (2012–2015)0.560.770.69 0.560.690.68
2nd (2016–2019)0.570.77 0.570.670.67
Lysimachia1st (2012–2015)0.590.590.390.510.390.530.51
2nd (2016–2019)0.530.42 0.630.420.530.52
Ozeros1st (2012–2015)0.710.450.530.490.450.510.55
2nd (2016–2019)0.720.62 0.520.520.620.62
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Latinopoulos, D.; Spiliotis, M.; Ntislidou, C.; Kagalou, I.; Bobori, D.; Tsiaoussi, V.; Lazaridou, M. “One Out–All Out” Principle in the Water Framework Directive 2000—A New Approach with Fuzzy Method on an Example of Greek Lakes. Water 2021, 13, 1776. https://doi.org/10.3390/w13131776

AMA Style

Latinopoulos D, Spiliotis M, Ntislidou C, Kagalou I, Bobori D, Tsiaoussi V, Lazaridou M. “One Out–All Out” Principle in the Water Framework Directive 2000—A New Approach with Fuzzy Method on an Example of Greek Lakes. Water. 2021; 13(13):1776. https://doi.org/10.3390/w13131776

Chicago/Turabian Style

Latinopoulos, Dionissis, Mike Spiliotis, Chrysoula Ntislidou, Ifigenia Kagalou, Dimitra Bobori, Vasiliki Tsiaoussi, and Maria Lazaridou. 2021. "“One Out–All Out” Principle in the Water Framework Directive 2000—A New Approach with Fuzzy Method on an Example of Greek Lakes" Water 13, no. 13: 1776. https://doi.org/10.3390/w13131776

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