Next Article in Journal
Recombinant Clostridium difficile Toxin Fragments as Carrier Protein for PSII Surface Polysaccharide Preserve Their Neutralizing Activity
Next Article in Special Issue
The Use of Recently Developed Histochemical Markers for Localizing Neurotoxicant Induced Regional Brain Pathologies
Previous Article in Journal
Response to John G. Brock-Utne. Comment on Hurley, J.C. Towards Clinical Application of Anti-endotoxin Antibodies; A Re-Appraisal of the Disconnect
Previous Article in Special Issue
Influence of Two Depuration Periods on the Activity and Transcription of Antioxidant Enzymes in Tilapia Exposed to Repeated Doses of Cylindrospermopsin under Laboratory Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trophic State and Toxic Cyanobacteria Density in Optimization Modeling of Multi-Reservoir Water Resource Systems

1
Hydraulic Sector, Department of Civil an Environmental Engineering and Architecture, University of Cagliari, via Marengo 2, Cagliari 09123, Italy
2
Sardinian Water System Management Authority, Regional Government, via Mameli 88, Cagliari 09123, Italy
*
Author to whom correspondence should be addressed.
Submission received: 8 January 2014 / Revised: 23 March 2014 / Accepted: 24 March 2014 / Published: 22 April 2014
(This article belongs to the Collection Marine and Freshwater Toxins)

Abstract

:
The definition of a synthetic index for classifying the quality of water bodies is a key aspect in integrated planning and management of water resource systems. In previous works [1,2], a water system optimization modeling approach that requires a single quality index for stored water in reservoirs has been applied to a complex multi-reservoir system. Considering the same modeling field, this paper presents an improved quality index estimated both on the basis of the overall trophic state of the water body and on the basis of the density values of the most potentially toxic Cyanobacteria. The implementation of the index into the optimization model makes it possible to reproduce the conditions limiting water use due to excessive nutrient enrichment in the water body and to the health hazard linked to toxic blooms. The analysis of an extended limnological database (1996–2012) in four reservoirs of the Flumendosa-Campidano system (Sardinia, Italy) provides useful insights into the strengths and limitations of the proposed synthetic index.

1. Introduction

Integrated Water Resources Management (IWRM) requires the analysis of quantitative, as well as qualitative aspects in a common strategy. Over the past few decades, a number of simulation and optimization models have been developed to manage multi-reservoir systems [3]. Mathematical optimization procedures can be constructed to efficiently solve an optimization model that adequately approximates the real water system management requirements, and a subsequent simulation phase can greatly narrow down the search for the optimum [4,5,6,7]. When considering the tools needed for adequately approaching mathematical optimization to real problems, water quality characterization requirements are critical components in constructing the optimization model. Optimization models defining water quality policy can be developed at different levels of complexity [8], from conceptual models to physically-based models, the latter often suffering from limited computational feasibility and great computational complexity.
In the Mediterranean region, eutrophication is one of the most serious problems affecting the quality of water stored in surface reservoirs. The increase in nutrients leads to the greater productivity of the water system, which may result in an excessive increase in algal biomass or other primary producers, such as macrophytes. Excessive algal biomass can seriously affect water quality, especially if it creates anaerobic conditions. A trophic state is a multidimensional concept, involving different aspects; nevertheless, single criterion alternatives to multi-parameter indices have also been developed. Several researchers have proposed relationships between algal biomass, nutrient concentration and nutrient load. Carlson [9] suggested that algal biomass is the key indicator for trophic state classification, and chlorophyll-a (Chl-a) is the best for estimating algal biomass in most reservoirs. The classification of the quality of water held in reservoirs based on the trophic state index (TSI) proposed by Carlson [9] can be inserted into an optimization model to make preliminary assessment of quality related to overall algal biomass concentrations in the period of maximum stratification. TSI can be easily used to insert management constraints on water use possibilities, to predict short-term trophic evolution in different climate scenarios and to check the relationships between parameters [2]. This approach can be adopted in simplified models for completely mixed lakes without stratification when considering complex multi-reservoir systems. Nevertheless, the use of TSI alone does not provide data concerning algal composition nor does it allow us to establish whether the algal blooms will be made up of a specific type of algae. Some families of microscopic phytoplanktonic algae produce toxins during water body eutrophication phenomena. Particularly in artificial reservoirs, toxin-producing phytoplankton belongs almost exclusively to Cyanobacteria. Toxic blooms of Cyanobacteria are a growing problem throughout the world. Monitoring of the toxicity of these blooms showed that almost half are indeed toxic. The most frequently found toxins, microcystins, are new, oncogenic-risk substances that cannot be eliminated by standard water purification processes. Special filters or costly treatments are needed to prevent them from entering the supply network. In such cases, limitations on resource use based solely on TSI values may not be sufficient. This article aims to present a combined TSI and algal classification system that could be used to provide a single parameter for classifying reservoir water quality in an optimization model framework treating complex water systems.
The paper is structured as follows: Section 2 briefly presents the different approaches from the literature on the identification of interrelations between different parameters to characterize the trophic state of lakes and reservoir; Section 3 focuses on the phenomena of harmful algal blooms (HAB); Section 4 describes the limnological survey and the evaluation of water quality in the period of 1996–2012 in four main reservoirs in the Flumendosa-Campidano water supply system (Sardinia, Italy); the water quality evaluation (QE) is based on the validation of Carlson’s approach [9] and the analysis of the density of algal species specifically considering the class of Cyanobacteria; finally, a synthetic quality index is proposed and evaluated to be inserted in the Water Resources Graphical Interface Decision Support System (WARGI-DSS) [10] in a predictive screening-level model. Some conclusions are drawn to synthesize the major results of the overall analysis and its limitations, and the need for further research in the field of water quality optimization modeling is noted.

2. Reservoir Trophic State Characterization

The trophic status in reservoirs is a multidimensional concept involving a wide range of data types (morphological, physical, chemical and biological) that are most notably related to human activity in the basins. To estimate future trends of the eutrophication rate under changing conditions, several models and statistical analyses use biomass, an indirect measure of which is given by Chl-a concentration. Several researchers [9,11,12,13] have found a log-linear relationship between Chl-a and nutrients (nitrogen and phosphorus). The input-output models developed by Dillon and Rigler [14] and by Vollenweider [15] can be used to predict phytoplankton biomass (expressed as Chl-a). Smith [13] developed a theoretical framework that reduces the error of chlorophyll prediction in the Dillon-Rigler model by taking into account the influence of the total phosphorus (TP)/total nitrogen (TN) weight ratio. Smith [13] showed that TN concentration influences chlorophyll concentration, even in lakes where phosphorus alone had previously been presumed to be a limiting factor. In these simple input-output models, other factors that can also influence the trophic status, such as biological interaction, internal nutrient loading and physical conditions, are not assessed. Although these models are only applicable in discrete periods of the year under the very restrictive assumption that a reservoir is a mixed system, these approaches are the most common in the literature.
When analyzing complex phenomena, understanding how these factors are related and inserting them into a model using analytical relations requires great effort. Brylinsky and Mann [16] analyzed the interactions between a large number of variables using a combination of hierarchical models and statistical analyses. Hierarchical models provided a framework for the identification of the relationships between these variables that were tested using correlation and factor analysis. The variables were grouped into those influencing primary production through nutrient availability (depending on the nutrient input, characteristics and area of the drainage basin, the nutrient dilution, precipitation, evaporation, surface area and mean depth and the nutrient distribution, wind, surface area and temperature range of water) and through energy availability (depending on the total incident radiation, reflection and turbidity). They used simple or multiple regression analysis for the comprehension of the type and degree of the relationship between morphological, chemical, physical and biological variables. The results of regression analysis showed that phytoplankton biomass and Chl-a are closely correlated with each other and with phytoplankton production. Multiple regression using biological indices to estimate phytoplankton production also showed that Chl-a and TP explain 45% of phytoplankton production variance.
The various attempts to analyze these complex phenomena in a simplified manner enable a description of the trophic condition of the reservoir using a mathematical modeling approach for complex water supply system management optimization.
To define the water quality index in complex multi-reservoir systems, we use the TSI proposed by Carlson [9], which in recent years seems to have gained general acceptance from the limnological community for characterizing a reservoir’s trophic state. TSI is evaluated using Chl-a, TP and Secchi disk (SD) transparency measurements. The TSI scale shown in Table 1 ranges from zero (ultra-oligotrophic) to 100 (hyper-eutrophic) [17]. High and/or increasing trophic status values indicate an increase in eutrophic conditions (higher biomass).
In accordance with recent legislation in European countries, and specifically in Italy [18], the environmental state is based on four macro descriptors (Chl-a; transparency; total phosphorous; hypolimnion oxygen) and summarized in a QE index of five possible numerical values:
  • QE = 1 → excellent;
  • QE = 2 → good;
  • QE = 3 → acceptable;
  • QE = 4 → poor;
  • QE = 5 → bad.
Combining QE and TSI, the proposed water quality classification for reservoirs is presented in Table 1.
Table 1. Relations between trophic state index (TSI) values, lake attributes and quality evaluation (QE) values.
Table 1. Relations between trophic state index (TSI) values, lake attributes and quality evaluation (QE) values.
QETSIAttributes
1<40Oligotrophy
240–50Mesotrophy
350–70Eutrophy
470–80Hypereutrophy
5>80Over-hypereutrophy

3. Problems Associated with the Presence of Phytoplankton in Water Stored in Reservoirs

In internal low hydro-dynamism water bodies (such as regulation reservoirs with a significant ratio between storage capacity and mean annual inflow), the microscopic organisms making up phytoplankton play a leading role in the primary production of living organic substance or biomass.
While phytoplankton plays a key role in lake ecosystems, its excessive proliferation becomes a serious problem for water stored in reservoirs used for potable and/or recreational supply. The immediate consequences of this excess growth are numerous, ranging from the simple abundance of suspended particles (phytoplankton itself, zooplankton, bacteria, fungi and detritus) to increased concentrations of ammonia, nitrites, hydrogen sulfide, methane, ethane and humic acids, to bad flavor and smell in fish and water, due to the presence of particular algae, and the possible development of toxic algae.
At times, the phenomenon is so obvious, that the naked eye can see the mass of microscopic algae, which produce blooms giving a particular coloring to the water body. The term “algal blooms” indicates a situation in which 80%–90% of the mass of microscopic algae consists of one or two species. In particular, Cyanobacteria may be considered as blooming when their cell number exceeds one million per liter. This coating of microscopic algae covers the surface of the water and decreases its transparency. This, in turn, prevents the penetration of sunlight, something that, coupled with the thermal stratification typical of reservoirs in the Mediterranean area, inevitably engenders conditions of anoxia and, hence, the above-mentioned consequences.
A further alarming aspect associated with eutrophication is the fact that algal species responsible for blooms belong to the taxonomic group of Cyanobacteria, present in surface waters all round the world. Cyanobacteria produce a wide range of toxins, which, according to their effects, are classified as hepatotoxins, neurotoxins, skin irritants and other toxins (Table 2) [19].
Table 2. Toxins and producer organisms.
Table 2. Toxins and producer organisms.
NameProducer
NeurotoxinsAnatoxin-aAnabaena, Aphanizomenon
Homo-anatoxin-aOscillatoria (Planktothrix)
Anatoxin-a(s)Anabaena, Oscillatoria (Planktothrix)
Paralytic shellfish poisons (saxitoxins)Anabaena, Aphanizomenon
HepatotoxinsCylindrospermopsinAphanizomenon, Cylindrospermopsis
MicrocystinsAnabaena, Aphanocapsa, Microcystis, Oscillatoria (Planktothrix)
NodularinsNodularia
Contact irritant-dermal toxinsAplysiatoxinSchizothrix
Cyanotoxins are produced and contained in the cells during their growth and released into the water mainly during cell lysis (e.g., for senescence or death) rather than by continuous excretion [20]. Thus, a water body populated by toxic species will inevitably become enriched with their toxins, which can affect human beings or animals using the water. The literature reports frequent cases of deadly poisoning of many animal species, both wild and domestic, following consumption of water with potentially toxic Cyanobacteria blooms, the poisoning being due to neurotoxins or hepatotoxins.
Neurotoxins are produced by species, such as Dolichospermum flos-aquae, Dolichospermum spiroides, Anabaena circinalis and Aphanizomenon flos-aquae, and genera, such as Oscillatoria, Gomphosphaeria and Trichodesmium. These organisms can, at the same time, produce hepatotoxins, but the latter act more slowly; therefore, the clinical picture is dominated by neurotoxins.
Hepatotoxins are more frequent than neurotoxins and are produced by Cyanobacteria belonging to a number of genera, the most common being Microcystis, Anabaena, Oscillatoria (Planktothrix), Nostoc, Nodularia, Aphanizomenon [21,22] and Cylindrospermopsis. Among hepatotoxins, the most abundant are microcystins, 50 types of which are known. Their effects on human beings and animals may be summarized as:
  • Acute hepatotoxicosis, due to direct ingestion;
  • Carcinogenic effect, if ingested in sub-acute doses over a long period of time (hepatic cancer);
  • Acute pneumonias, if breathed in.
In Italy, at the present time, microcystins each year affect the reservoirs of seven regions out of twenty. While the World Health Organization has addressed the problem of Cyanobacteria blooms, neither the European Union (EU) nor Italian legislation have established concentration limits for cyanotoxins in water for human consumption and/or recreational use. On the assumption of a three-stage water purification system [2]. In this study, we adopted a threshold value of Cyanobacteria density for potable use (and for cautionary purposes, also for irrigation use) of 100 × 106 cell/L.

4. Limnological Survey and Water QE in the Flumendosa-Campidano Reservoirs

The Flumendosa-Campidano water supply system extends over southeastern Sardinia (Italy), reaching to the center of the Island. The hydrology is typically Mediterranean, with the alternation of several droughty years with years marked by intense rainfall. Besides being the most extensive on the island, this water system is also the most complex, since it interconnects with other systems, and it is a multi-reservoir and multi-use system. Its pivot is a series of reservoirs, linked in a cascading sequence, whence depart pressure pipelines for residential use and open channels supplying irrigation water to the Campidano plain. An extended description of the main characteristics of the system and the connection scheme made using the graphical interface of WARGI-DSS are reported in Sechi and Sulis [10].
Since the early 1990s, the Sardinian Water System Management Authority (Ente Acque della Sardegna, ENAS) has conducted an intensive monitoring program to identify the water quality status in the most important reservoirs in the system, particularly the following four reservoirs, whose main characteristics are shown in Table 3: (1) Flumendosa at Nuraghe Arrubiu (latitude = 39°42′38′′ N; longitude = 9°17′25′′ E); (2) Mulargia at Monte su Rei (latitude = 39°37′30′′ N; longitude = 9°14′23′′ E); (3) Is Barrocus (latitude = 39°45′10′′ N; longitude = 9°5′32′′ E); and (4) Cixerri at Genna Is Abis (latitude = 39°16′53′′ N; longitude = 8°52′56′′ E).
Table 3. The main characteristics of the reservoirs in the Flumendosa-Campidano system.
Table 3. The main characteristics of the reservoirs in the Flumendosa-Campidano system.
Reservoir featuresCixerriMulargiaFlumendosaIs Barrocus
Total catchment basin (km2)426.001183.161004.5193.00
Reservoir surface at maximum level (km2)4.9012.409.006.30
Elevation at maximum regulation level (m)39.00258.00267.00413.00
Volume at maximum regulation level (×106 m3)23.90320.70292.9011.96
Measurements of the main chemical and biological parameters have been performed by the ENAS since the appearance of algae blooms in reservoirs in the early 1990s. These data seem to suggest that eutrophication and high mineral nutrient load in inflowing water are indeed the major water quality problems of the system. In general terms, only the two larger reservoirs in the system, Flumendosa and Mulargia, have a trophic level variable between oligotrophy and mesotrophy [23]. The other two lakes, Cixerri and Is Barrocus, are in a state of more or less advanced eutrophication, at times hovering between eutrophy and hypertrophy [24]. ENAS measured sampling parameters selected to represent the trophic conditions of the reservoirs, particularly Chl-a, TP and SD transparency. Measurements were taken at least once a month from January 1996 to December 2012, at different depths in all reservoirs.
Following several studies [11,12,13,14], we used Chl-a as a simple and useful estimator of phytoplankton biomass and found, for the reservoirs in Table 3, a relationship between the average concentration of Chl-a in the summer by using measures of TP during the spring (Figure 1):
log10[Chl] = 1.065log10 TP − 0.758
Figure 1. Summer average chlorophyll-a (Chl-a) concentration vs. total phosphorus (TP) concentration at the spring overturn in four reservoirs.
Figure 1. Summer average chlorophyll-a (Chl-a) concentration vs. total phosphorus (TP) concentration at the spring overturn in four reservoirs.
Toxins 06 01366 g001
We found that Chl-a in the four reservoirs of the system was quite accurately related to phosphorus concentration (R = 0.83). This is a significant factor in water quality management, because it makes it possible to predict the trophic state of reservoirs once the Chl-a concentration is known.
A log-log transformation of the SD transparency and Chl-a concentration data (Figure 2) was used to define the resulting log-linear combination (R = 0.64):
log10[SD] = 0.551 − 0.328 log10[Chl]
where the SD transparency is in meters and the Chl-a concentration is in milligrams per cubic meter taken near the surface. One possible explanation for the exponential form of Equation (2) is given by Steele [25].
We defined the TSI equations specifically suited for the Flumendosa-Campidano water system based on the trophic state definition of the Carlson’s approach [9]:
TSI = 10 × (6 − log2 SD)
Combining Equation (3) with the regression of Chl-a against TP (Equation (1)) and of SD against Chl-a (Equation (2)), site-specific equations for TSI in the Flumendosa-Campidano water system were assessed. Specifically, TSI can be computed from SD transparency, Chl-a concentration and TP using the equations in the forms:
TSI(SD) = 60 − 14.41ln SD
TSI(Chl) = 2.05 ln Chl + 52
TSI(TP) = 2.19 ln TP + 48.47
Figure 2. Secchi disk (SD) transparency vs. Chl-a concentration near the surface in four reservoirs.
Figure 2. Secchi disk (SD) transparency vs. Chl-a concentration near the surface in four reservoirs.
Toxins 06 01366 g002
As previously explained, the TSI based on Chl, TP and SD can be obtained for the classification of reservoirs. Using the previous equations (Equations (4)−(6)), TSIs have been calculated from the observed data, and Table 4 reports the statistical values.
Table 4. Statistical summaries of the TSI data evaluated from the samples.
Table 4. Statistical summaries of the TSI data evaluated from the samples.
LakeWater quality parameterNumber of samplesMax.Min.MeanStandard deviation
FlumendosaTSI (Chl)1463058454.4
TSI (TP)1455533443.2
TSI (SD)1466526416.8
CixerriTSI (Chl)2296936565.8
TSI (TP)2496446552.6
TSI (SD)2298738657.6
Is BarrocusTSI (Chl)2496537514.8
TSI (TP)2504359503.1
TSI (SD)2527036546.8
MulargiaTSI (Chl)2396236484.5
TSI (TP)2375637463.1
TSI (SD)2377330479.2
As regards phytoplankton composition, the microscopic observations were carried out on fresh and Lugol fixed samples in 10-mL sedimentation chambers, using an inverted microscope (Axiovert 100, Carl Zeiss Microscopy, LLC, Thornwood, NY, USA). The cell density was determined according to Utermöhl [26] and Lund et al. [27]. Species were determined according to Anagnostidis and Komárek [28], Komárek and Anagnostidis [29,30], Huber-Pestalozzi [31], Germain [32], and Cronberg and Annadotter [33]. Species in the four reservoirs mostly belong to five classes: Cyanobacteria, Chlorophyceae, Bacillariophyceae (Diatoms), Cryptophyceae and Conjugatophyceae. Furthermore, Dinophyceae and Chrysophyceae were recorded, but not analyzed, in this paper, because of their low densities. To assess the contribution of each class to total density, each reservoir density value, collected between 1996 and 2012, was processed as follows:
  • Mean annual density values (cells/liter) were calculated for each class (absolute and percent value);
  • Minimum and maximum values reached by each class in the period under study (1996–2012) were recorded;
  • The most representative species in each class were identified.
Figure 3, Figure 4, Figure 5 and Figure 6 show the data processing output. The graphic depiction enables us to gain an idea of how the density of each class changed over time in each reservoir.
Figure 3. The densities (percentage of cells per liter) of (a) classes and (b) species in Cixerri reservoir.
Figure 3. The densities (percentage of cells per liter) of (a) classes and (b) species in Cixerri reservoir.
Toxins 06 01366 g003
In the Cixerri reservoir (Figure 3), the dominance of Cyanobacteria has been particularly massive in the period between 1996 and 2002, with a maximum peak density value in 2000 equal to 99%, corresponding to a density of 736.051 (×103 cell/L). The Aphanizomenon genus (Nostocales order) presented dominance with the species A. flos-aquae, A. gracile and A. issatschenkoi. Some of these blooms have passed one million cells per liter [34]. Negative toxicity assays (enzyme-linked immunosorbent assay (ELISA) kit for saxitoxin) were obtained for these species. The interesting presence should be noticed, in September 2012, of the tropical genus, Cylindrospermopsis, which was tested using the ELISA kit for cilyndrospermopsin for the potential toxicity of some of its strains. All the assays were negative. Within the Dolichospermum genus, D. flos-aquae have dominated more frequently than D. spiroides and D. planctonicum. Within Oscillatoriales, the Planktothrix genus showed a variable concentration with high values in the early years of the survey and in the last two years. Regarding the genus, Microcystis, M. aeruginosa presented several flowering episodes with high density values.
Figure 4. The densities (percentage of cells per liter) of (a) classes and (b) species in Is Barrocus reservoir.
Figure 4. The densities (percentage of cells per liter) of (a) classes and (b) species in Is Barrocus reservoir.
Toxins 06 01366 g004
With regard to the Is Barrocus reservoir (Figure 4), the long-term dynamic shows that the most abundant class over the period of 1996–2002 was that of Cyanobacteria with a peak in the year of 2009 equal to 64.749 (×103 cell/L) composed by Aphanocapsa sp. and Anathece sp. at 93%. Furthermore, Aphanocapsa and Anathece gave a strong contribution in terms of density over the entire period, while within the Nostocales, the Dolichospermum with the species D. flos-aquae showed a significant density over the period of 1996–2004 and regressed in the subsequent years, almost tending to disappear. Another Nostocales that significantly contributed to the composition over the period of 1999–2002, is represented by Anabaena aphanizomenoides. Microcystis with M. aeruginosa showed a similar trend with some blooms until 2006, while disappearing in the following years.
Figure 5. The densities (percentage of cells per liter) of (a) classes and (b) species in Flumendosa reservoir.
Figure 5. The densities (percentage of cells per liter) of (a) classes and (b) species in Flumendosa reservoir.
Toxins 06 01366 g005
In Flumendosa reservoir (Figure 5), the dominance was of the Bacillariophyceae and Chlorophyceae classes during the early investigated period. Within the class of Chlorophyceae, the species, Chlorella sp. and Coelastrum reticulatum, show the highest average density from 1996 to 1999. We observed a massive dominance of the class of Cyanobacteria from 2000 to 2012, with a significant increase in terms of the average density, especially in the last two years. The Cyanobacteria with the highest density values was the Planktothrix agardhii-rubescens group, which was very abundant in the years of 2000–2004, with a peak in 2002 of 35.230 (×103 cell/L) corresponding to 99% of the total density class. Planktothrix agardhii-rubescens has been tested over the years for their potential toxicity. Furthermore, this species has been responsible for frequent and abundant flowering events and showed an increase in the density values from 2001 to 2012. The Microcystis genus was responsible for the onset of blooms in the reservoir (Figure 7), even if with a lower impact than the Planktothrix genus. In the species composition, the Dolichospermum genus significantly disappeared. Only one bloom occurred in 2003, and sporadically low density values were recorded in the subsequent years.
A dominance of Chlorophyceae with Coelastrum reticulatum and Pediastrum simplex was noted in the Mulargia reservoir from 1996 to 1999 (Figure 6) with a peak density of 2.618 (×103 cell/L) in 1999. In the long-term dynamic of algal classes from 2000 to 2012, there is a clear dominance of Cyanobacteria and a good correspondence with the species dynamic. The Cyanobacteria species with the highest occurrence was of the Planktothrix agardhii-rubescens group, which became more and more frequent over the years with the occurrence of algal blooms (Figure 8) and which presented a peak density in 2002 with 16.7960 (×103 cell/L) with a dominance of 99%. These species were observed until 2012. Another species of Cyanobacteria contributing significantly to the density of the entire class was Microcystis sp. with an important bloom in 2001 recording a density of 115.390 (×103 cell/L) with a dominance of 80%.
Figure 6. The densities (percentage of cells per liter) of (a) classes and (b) species in Mulargia reservoir.
Figure 6. The densities (percentage of cells per liter) of (a) classes and (b) species in Mulargia reservoir.
Toxins 06 01366 g006
Figure 7. Bloom of Microcystis spp. in Flumendosa reservoir.
Figure 7. Bloom of Microcystis spp. in Flumendosa reservoir.
Toxins 06 01366 g007
Figure 8. Bloom of Planktothrix agardhii-rubescens group in Mulargia reservoir.
Figure 8. Bloom of Planktothrix agardhii-rubescens group in Mulargia reservoir.
Toxins 06 01366 g008
Data processed up to this point indicate that in these four reservoirs, by far the most dominant class in the period under consideration was that of Cyanobacteria, responsible for extensive and regular blooms in the system’s reservoirs. Furthermore, the Chlorophyceae, mainly with Chlorella, Coelastrum, Oocystis, Pediastrum and Scenedesmus, showed significant density values. The third class is represented by Bacillariophyceae with Cyclotella, Aulacoseira, Fragilaria and Asterionella, the most frequently found. During algal bloom episodes, their Cyanobacteria concentration can be as high as several million cells per liter, thus creating a significant restraint on the use of the water resource. If we add to this the fact that all three genera (Planktothrix, Aphanizomenon and Microcystis) are known as being toxin-producing, we can easily see how important it is to monitor the quality of water bodies, also from the viewpoint of their phytoplankton composition.
Table 5 shows the results of the algal toxicity test through the water column (0–20 m) carried out in Mulargia and Flumendosa reservoirs during four events of Planktothrix agardhii-rubescens blooms in the period of 2010–2012. The results of immunoassays have highlighted the presence of toxic strains of microcystins in both reservoirs. Recently, the “agardhii-rubescens group” has received greater attention as being responsible for frequent toxic blooms. The table reports the mean density values of the water column of P. agardhii-rubescens, the corresponding date, the intracellular and extracellular concentrations of the microcystins produced, as well as their total value. In fact, the toxins were extracted following the procedure described by Meriluoto and Eriksson [35], which allowed one to separate the two fractions, intracellular and extracellular, which were subsequently treated and assayed with an enzyme immunoassay ELISA. Results show that the highest levels of toxins in Flumendosa were found inside the cells of Planktothrix agardhii-rubescens (1.35 µg/L) with a maximum value in dissolved form equal to 0.37 at 2.5 m in depth during the spring of 2011. In Mulargia, the species presented the highest values of toxins inside the cells, with values almost always >2 µg/L and reaching values >5 µg/L at 10 m in depth, those maximum values exceeding the range (0–5 ppb) of the enzyme immunoassay test (ELISA microtiter plate microcystins-ADDA) validity. The extracellular concentrations of microcystins were always lower than 1 µg/L, which is the threshold value for total microcystin-LR (free plus cell-bound) in drinking water established by the World Health Organization (WHO) guideline [36].
Table 5. Concentrations of microcystins produced by Planktothrix agardhii-rubescens in Flumendosa and Mulargia reservoirs at different depths from four sampling dates.
Table 5. Concentrations of microcystins produced by Planktothrix agardhii-rubescens in Flumendosa and Mulargia reservoirs at different depths from four sampling dates.
ReservoirsDateDepth (m)Density (×106 cells L−1)MCYSTs (in)* (µg/L)MCYSTs (ex)** (µg/L)MCYSTs (tot)*** (µg/L)
FlumendosaMarch 20111.0851.30.41.6
2.5841.40.41.5
5.0761.10.21.3
10.01191.00.21.2
MulargiaMarch 20100.0552.90.02.9
1.0592.70.02.7
2.5592.00.02.0
5.0401.60.01.6
7.5252.50.02.5
10.0322.10.02.1
15.0261.70.01.7
20.0222.30.02.3
March 20122.5441.30.31.6
7.5481.30.31.6
10.052>5.00.4>5.0
15.0341.60.31.9
May 20127.5422.70.23.0
10.0105>5.00.4>5.0
15.0811.70.32.0
20.0761.20.31.4
* MCYSTs (in) = Intracellular microcystins concentration; ** MCYSTs (ex) = Extracellular microcystins concentration; and ***MCYSTs (tot) = Total microcystins concentration.

5. Quality Characterization in Management Optimization Modeling

To provide a preliminary evaluation of water system management performances, it is necessary to consider various synthetic quantity and quality datasets that define different hydrological and demand scenarios and related management rules to be optimized. This approach has been partially performed using the WARGI-DSS; here, a more extended time horizon than the historical dataset of Sulis et al. [2] is presented to define a comprehensive predictive modeling approach.
A higher level of uncertainty is to be expected in system forecasting using a synthetic quality index dataset (scenario) in the WARGI-DSS optimization of system management. Dealing with synthetic scenarios, several types of relationships have been considered in the regression analysis with the aim of modeling the dependency of TSI from other data characterizing reservoirs. Indeed, in order to determine the regression equation, a simple linear least-squares curve was determined between hydrological parameters and TSI values at each reservoir. More complicated approaches do not produce significantly better results, at least for the reservoirs under consideration. Therefore, the equation here, used to model the dependence of stored water TSIs on stored volumes in a reservoir, is a simple multiple linear regression that incorporates month-explanatory TSI variation to allow the curve’s slope to change depending on the month of the year.
As stated above, monthly TSI values in the period between 1996 and 2012, were used to calibrate the regression model in the following form:
TSIt = aTSI(t−1) + bVt + cT + d
where:
TSLt 
evaluated variable at current period t;
TSL(t−1) 
evaluated variable at the previous time period;
T 
transformed time date;
Vt 
stored volume in the reservoir at current period t;
a,b,c,d 
regression coefficients.
The number of data pairs, correlation coefficient R and residual standard error are listed in Table 6.
Table 6. Regression results for TSI evaluation in Flumendosa, Mulargia, Is Barrocus and Cixerri reservoirs.
Table 6. Regression results for TSI evaluation in Flumendosa, Mulargia, Is Barrocus and Cixerri reservoirs.
ParameterFlumendosaMulargiaIs BarrocusCixerri
Coefficient R0.520.560.590.66
Standard Error4.253.923.714.22
Number of samplings145186197242
A large number of chemical, physical and biological factors and their interactions determine the Cyanobacteria density in a reservoir [37,38]. It is widely accepted that TP and TN are the most important factors enhancing the growth of Cyanobacteria in situ. In general, nutrient concentrations influence the growth of Cyanobacteria more than the TN:TP ratio. Several observations do not support the significance of the TN:TP ratio [39,40]. The phosphorus-Chl-a models of Dillon and Rigler [14] and Vollenweider [15], already used to implicate phosphorus as an important nutrient in phytoplankton biomass, have also been useful for predicting the concentration of Cyanobacteria or microcystin in reservoirs. Using regression analysis, the best fits for both Chl-a and Cyanobacteria suggest that similar patterns of change in Chl-a and Cyanobacteria occur with increasing TP concentration [41].
Preliminary analyses have been made to verify the interaction over time between Chl-a and Cyanobacteria in the main reservoirs of the Flumendosa-Campidano system and to evaluate the possibility of forecasting Cyanobacteria density. As stated above, monthly TSI values in the period between 1996 and 2005, were used to calibrate the regression model in the following form:
ln[D(cyano)t] = a∙ln[D(cyano)t−1] + bVt + cT + d
where:
D(cyano)t
evaluated variable at current period t;
D(cyano)t−1
evaluated variable at the previous time period;
T 
transformed time date;
Tt 
stored volume in the reservoir at current period t;
a,b,c,d 
regression coefficients.
The number of data pairs, correlation coefficient R and residual standard error are listed in Table 7.
Table 7. Regression results for Cyanobacteria density evaluation in Flumendosa, Mulargia, Is Barrocus and Cixerri reservoirs.
Table 7. Regression results for Cyanobacteria density evaluation in Flumendosa, Mulargia, Is Barrocus and Cixerri reservoirs.
ParameterFlumendosaMulargiaIs BarrocusCixerri
Coefficient R0.480.620.430.58
Standard Error1.321.341.201.52
Number of samplings144132171178
Following Sulis et al. [2], a combined approach evaluates the QE index in a reservoir, j, at time t, using both TSI and D(cyano):
Toxins 06 01366 i001
Considering the threshold of 100 × 106 cell/L, in the case of lower D(cyano), the QE index is based only on TSI as reported in Table 1; otherwise, the QE value is assumed to equal five (water unsafe for human consumption) no matter the TSI value. Based on the value of the state parameters of the water system (stored volume V in each reservoir, j, at time t) under a synthetic hydrological scenario able to reproduce historical drought periods [42], Equations (9)–(11) give the time series of QE values to be inserted in the WARGI-DSS in a predictive screening-level model. This integrated approach can provide a preliminary estimation of the impact of different planning alternatives on some quantity and quality issues at the water system scale. Furthermore, it provides an estimate of how the uncertainty of quality data and parameter values can affect different water uses.

6. Conclusions

The definition of minimum quality requirements with regard to possible uses is an essential element in the management of complex, multi-user and multi-resource water supply systems, since it limits the actual availability of the resource for use. This paper presents an approach to define a QE index to be inserted in a mathematical optimization model for complex water system management. Details on using the water quality index in optimization modeling have been presented in Sulis et al. [2]. In this paper, a more extended database is considered to properly define interactions between the most important parameters related to trophic status of lakes and reservoirs, and the relationships from the literature are validated in the Flumendosa-Campidano water system. With regards to this system, the paper presents log-linear relationships between Chl-a, TP and SD to define new TSI equations. Furthermore, interactions over time between Chl-a, Cyanobacteria and stored volume were analyzed to evaluate the possibility of forecasting Chl-a and Cyanobacteria density. Further efforts will include the validation of the proposed equation to a larger number of reservoirs in the Island of Sardinia. The analysis shows that the Cyanobacteria was the dominant phytoplankton class in the four reservoirs, and the Planktothrix genus was the most abundant in three of them. The results highlight the necessity of paying greater attention to these toxic species in the IWRM of lakes and reservoirs affected by these flowering events. Since the proposed approach is a semi-quantitative immunoassay, a more in-depth and analytical approach could be of value in order to link flowering events with the toxic nature of the responsible species.

Acknowledgments

The authors gratefully acknowledge the Sardinia Regional Government for the financial support (Programma Operativo Regionale del Fondo Sociale Europeo (P.O.R. Sardinia F.S.E.) Operational Programme of the Autonomous Region of Sardinia, European Social Fund 2007–2013; Axis IV Human Resources, Objective I.3, Line of Activity I.3.1).

Author Contributions

Andrea Sulis has carried out all the statistical analysis reported in the paper and the data analysis reported in Paragraph 5; Andrea Sulis and Giovanni M. Sechi wrote Paragraphs 1–4 and the part of Paragraph 5 related to the TSI estimation. Also Andrea Sulis has coordinated all authors’ efforts and contributions to the paper.
Oriana Soru has carried out the analysis of phytoplankton density and Elisa test; Paola Buscarinu and Oriana Soru have processed the data, elaborate tables and graphs and wrote the part of the paper related to the cyanobacteria.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sechi, G.M.; Sulis, A. Multi-reservoir system optimization using chlorophyll-a trophic indexes. Water Resour. Manag. 2007, 21, 849–860. [Google Scholar] [CrossRef]
  2. Sulis, A.; Buscarinu, P.; Sechi, G.M. Using reservoir trophic-state indexes in optimisation modelling of water-resource systems. Environ. Model. Softw. 2011, 40, 241–225. [Google Scholar]
  3. Loucks, D.P.; van Beek, E. Water Resources Systems Planning and Management: An Introduction to Methods, Models and Applications; The United Nations Organization for Education, Science and Culture (UNESCO) Press: Paris, France, 2005. [Google Scholar]
  4. Loucks, D.P.; Stedinger, J.R.; Haith, D.A. Water Resource Systems Planning and Analysis; Prentice Hall: Englewood Cliffs, NJ, USA, 1981. [Google Scholar]
  5. Yeh, W.-G. Reservoir management and operations models: A state-of-the-art review. Water Resour. Res. 1985, 25, 1797–1818. [Google Scholar] [CrossRef]
  6. Simonovic, S.P. Tools for water management: One view of the future. Water Int. 2000, 25, 76–88. [Google Scholar] [CrossRef]
  7. Labadie, J.W. Optimal operation of multireservoir systems: State-of-the-art review. J. Water Resour. Plan. Manag. 2004, 130, 93–111. [Google Scholar] [CrossRef]
  8. Schlüter, M.; Savitsky, A.G.; McKinney, D.C.; Lieth, H. Optimizing long-term water allocation in the Amudarya River delta: A water management model for ecological impact assessment. Environ. Model. Softw. 2005, 20, 529–545. [Google Scholar] [CrossRef]
  9. Carlson, R.E. A trophic state index for lakes. Limnol. Oceanogr. 1977, 22, 361–369. [Google Scholar] [CrossRef]
  10. Sechi, G.M.; Sulis, A. Water system management through a mixed optimization-simulation approach. J. Water Resour. Plan. Manag. 2009, 135, 160–170. [Google Scholar] [CrossRef]
  11. Sakamoto, M. Primary production by phytoplankton community in some Japanese lakes and its dependence on lake depth. Arch. Hydrobiol. 1966, 62, 1–28. [Google Scholar]
  12. Schindler, D.W. Factors regulating phytoplankton production and standing crop in the world’s freshwater. Limnol. Oceanogr. 1978, 23, 478–486. [Google Scholar] [CrossRef]
  13. Smith, V.H. The nitrogen and phosphorus dependence of algal biomass in lakes: An empirical and theoretical analysis. Limnol. Oceanogr. 1982, 27, 1101–1112. [Google Scholar] [CrossRef]
  14. Dillon, P.J.; Rigler, F.H. The phosphorus-chlorophyll relationship in lakes. Limnol. Oceanogr. 1974, 19, 767–773. [Google Scholar] [CrossRef]
  15. Vollenweider, R.A. Input-output models with special reference to the phosphorus loading concept in limnology. Schweiz. Z. Hydrol. 1975, 37, 53–84. [Google Scholar]
  16. Brylinsky, M.; Mann, K.H. An analysis of factors governing productivity in lakes and reservoirs. Limnol. Oceanogr. 1973, 18, 1–14. [Google Scholar] [CrossRef]
  17. Carlson, R.E.; Simpson, J. A Coordinator’s Guide to Volunteer Lake Monitoring Methods; North American Lake Management Society: Madison, WI, USA, 1966. [Google Scholar]
  18. Norme in Materia Ambientale; Decree Law 152; Gazzetta Ufficiale: Rome, Italy, 2006.
  19. Guidelines for Drinking Water Quality: Addendum to Volume 2: Health Criteria and Other Supporting Information; World Health Organization (WHO): Geneva, Switzerland, 1998.
  20. Sivonen, K.; Jones, G. Cyanobacterial Toxins. In Toxic Cyanobacteria in Water: A Guide to Public Health Significance, Monitoring and Management; Chorus, I., Bartram, J., Eds.; E & FN Spon: London, UK, 1999; pp. 41–111. [Google Scholar]
  21. Harada, K.; Kondo, F.; Lawton, L. Laboratory Analysis of Cyanotoxins. In Toxic Cyanobacteria in Water; Chorus, I., Bartram, J., Eds.; E & FN Spon: London, UK, 1999; pp. 369–399. [Google Scholar]
  22. Carmichael, W.W. Health Effects of Toxin Producting Cyanobacteria: “The Cyanohabs”. Freshwater Harmful Algal Blooms: Health Risk and Control Management; Istituto Superiore di Sanità: Roma, Italy, 2000; pp. 70–80. [Google Scholar]
  23. Sechi, N. La Struttura algale e lo stato trofico dei laghi della Sardegna. In Rapporti ISTISAN, Workshop, Aspetti sanitari della problematica dei cianobatteri nelle acque superficiali italiane; Istituto Superiore di Sanità: Roma, Italy, 2000; pp. 101–110. (in Italian) [Google Scholar]
  24. Sechi, N.; Buscarinu, P.; Lugliè, A. Limnologia del nuovo invaso Cixerri nei primi due anni di esistenza. In Proceedings of the A.I.O.L.: Atti del 10. Congresso dell'Associazione Italiana di Oceanologia e Limnologia, Genova, Italy, 4–6 November 1992; pp. 201–208. (in Italian).
  25. Steele, J.H. Environmental control of phytosynthesis in the sea. Limnol. Oceanogr. 1962, 7, 137–150. [Google Scholar] [CrossRef]
  26. Utermöhl, H. Neue wege in der quantitativen Erfassung des Planktons. Verh. Int. Ver. Limnol. 1931, 5, 567–569. (in German). [Google Scholar]
  27. Lund, J.W.G.; Kipling, C.; le Cren, E. The inverted microscope method of estimating algal numbers and the statistical basis of estimations by counting. Hydrobiology 1958, 11, 143–170. [Google Scholar] [CrossRef]
  28. Anagnostidis, K.; Komárek, J. Modern approach to the classification system of cyanophytes. 3-Oscillatoriales. Algol. Stud. 1990, 50–53, 327–472. [Google Scholar]
  29. Komárek, J.; Anagnostidis, K. Modern approach to the classification system of Cyanophytes 4-Nostocales. Arch. Hydrobiol. 1989, 56, 247–345. [Google Scholar]
  30. Komárek, J.; Anagnostidis, K. Süßwasserflora von Mitteleuropa: Cyanoprokaryota; Spektrum Akademischer Verlag: Heidelberg, Germany, 2005. (in German) [Google Scholar]
  31. Huber-Pestalozzi, G. Das Phytoplankton des Süßwassers; Systematic und Biologie: Stuttgard, Germany; pp. 1938–1983. (in German)
  32. Germain, H. Flore des Diatomées: Eaux Douces et Saumâtres; Societe nouvelle des Editions Boubee: Paris, France, 1981. (in French) [Google Scholar]
  33. Cronberg, G.; Annadotter, H. Manual on Aquatic Cyanobacteria; International Society for the Study of Harmful Algae (ISSHA): Copenhagen, Denmark, 2006. [Google Scholar]
  34. Sechi, N.; Buscarinu, P.; Pilo, E.; Lugliè, A. Il fitoplancton del lago Cixerri nei primi cinque anni di esistenza. In Proceedings of the A.I.O.L.: Atti del 12. Congresso dell'Associazione Italiana di Oceanologia e Limnologia, Genova, Italy, 18–21 September 1996; Volume 2, pp. 41–51. (in Italian).
  35. Meriluoto, J.A.O.; Eriksson, J.E. Rapid analysis of peptide toxins in cyanobacteria. J. Chromatogr. 1988, 438, 93–99. [Google Scholar] [CrossRef]
  36. Guidelines for Drinking Water Quality; World Health Organization: Geneva, Switzerland, 2011.
  37. Pick, F.R.; Lean, D.R. The role of macronutrients (C, N, P) in controlling cyanobacterial dominance in temperate lakes. N. Z. J. Mar. Freshw. Res. 1987, 21, 425–434. [Google Scholar] [CrossRef]
  38. Dokulil, M.T.; Teubner, K. Cyanobacterial dominance in lakes. Hydrobiologia 2000, 438, 1–12. [Google Scholar] [CrossRef]
  39. Xie, L.; Xie, P.; Li, S.; Tang, H.; Liu, H. The low TN:TP ratio, a cause or a result of Microcystis blooms? Water Res. 2003, 37, 2073–2080. [Google Scholar] [CrossRef]
  40. Vaitomaa, J. The Effects of Environmental Factors on Biomass and Microcystin Production by the Freshwater Cyanobacterial Genera Microcystis and Anabaena. Ph.D. Thesis, University of Helsinki, Helsinki, Finland, 27 June 2006. [Google Scholar]
  41. Haney, J.F.; Ikawa, M. A Survey of 50 NH Lakes for Microcystins (MCs); Final Report; University of New Hampshire: Durham, NH, USA, 2000. [Google Scholar]
  42. Arena, C.; Cannarozzo, M.; Mazzola, M.R. Multi-year drought frequency analysis at multiple sites by operational hydrology—A Comparison of methods. Phys. Chem. Earth 2006, 31, 1146–1163. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Sulis, A.; Buscarinu, P.; Soru, O.; Sechi, G.M. Trophic State and Toxic Cyanobacteria Density in Optimization Modeling of Multi-Reservoir Water Resource Systems. Toxins 2014, 6, 1366-1384. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins6041366

AMA Style

Sulis A, Buscarinu P, Soru O, Sechi GM. Trophic State and Toxic Cyanobacteria Density in Optimization Modeling of Multi-Reservoir Water Resource Systems. Toxins. 2014; 6(4):1366-1384. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins6041366

Chicago/Turabian Style

Sulis, Andrea, Paola Buscarinu, Oriana Soru, and Giovanni M. Sechi. 2014. "Trophic State and Toxic Cyanobacteria Density in Optimization Modeling of Multi-Reservoir Water Resource Systems" Toxins 6, no. 4: 1366-1384. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins6041366

Article Metrics

Back to TopTop