SANTOS, MARINS, AGUIAR, CHALAR, SILVA, and LACERDA: Hydrochemistry and Trophic State Change in a Large Reservoir in the Brazilian Northeast Region Under Intense Drought Conditions

Hydrochemistry and Trophic State Change in a Large Reservoir in the Brazilian Northeast Region Under Intense Drought Conditions

Abstract

The study shows changes on physical and chemical water parameters and of trophic state in a large reservoir in the Brazilian semiarid region following decreasing reservoir volume due to rainfall shortage during four consecutive years. The monitoring period, between November 2011 and May 2014, assessed approximately 50% water volume reduction and 10 meters’ decrease of reservoir water level that degraded water quality. Decrease in reservoir volume, strong evaporation and the permanent influence of anthropogenic activities, favored the concentration of salts and accumulation of nutrients and of increasing pH. Thermal stratification of the water column occurred when volume was maximum and lead to a significant reduction in dissolved oxygen in the hypolimnion (0.07 to 2.62 mg L–1). Diminishing volume resulted in mixing of the hypolimnion nutrient-rich and oxygen-poor waters in the entre water column and changed the initial oligotrophic condition to eutrophic. However, the temporal scale of the response of the reservoir’s trophic state differs in the different areas of the reservoir. Whereas deeper areas accumulating nutrients from aquaculture and agriculture progressively became mesotrophic and eventually eutrophic; shallower regions far from direct anthropogenic influences, changed their trophic sate much later, but rapidly turned into super-eutrophic conditions, probably due to more intense sediment resuspension and water mixing. Trophic State Index followed nutrient increase during most of the period. However, it also responded to an increase in chlorophyll a concentrations when the reservoir achieved its minimum volume, in particular in the shallower areas. The results suggest that this type of reservoir systems are vulnerable to eutrophication during extended drought periods and that a better assessment of the maximum support capacity for reservoir activities, particularly aquaculture, must be re-assessed taking into consideration worst case scenarios forecasted by global climate change.




INTRODUCTION

Reservoirs in semiarid regions, including those in Northeastern Brazil, witness unpredictable periods of rainfall when they receive extremely low watershed water inflow, exhibit high evaporation rates and a longer water residence time. These processes are the primary cause of water quality decrease due to increasing concentrations of dissolved salts (Freire et al., 2009). This scenario results in nutrient accumulation and concentration, leading to increasing algal density and the frequency of cyanobacteria blooms, thus rendering these systems much more vulnerable to eutrophication. The recent development of intensive fish cage aquaculture in many reservoirs in the Brazilian semiarid may contribute to a further deterioration of water quality, in particular during drought periods (Oliveira et al., 2015).

The trophic state of reservoirs in semiarid regions may vary according to reservoir volume (Braga et al., 2015), precipitation (Chaves et al., 2013), inputs of external loads of nutrients from surrounding soils (Lopes et al., 2014; Santos et al., 2014) and internal processes such as aquaculture (Bezerra et al., 2014). However, how these drivers interact to affect the trophic state of reservoirs are still poorly understood. Between 2011 and 2014, the northeastern semiarid region of Brazil witnessed a prolonged period of drought, resulting in a drastic reduction of the storage volume of artificial reservoirs. For example, the Castanhão reservoir, the largest multiuse reservoir in the region registered a drop from 88% to 27% of its storage capacity, rendering a unique opportunity to improve our understanding the relationship between volume and trophic state and the influence of anthropogenic pressures. The understanding of the reservoir’s response during this period is of high significance since several externalities act to enhance the negative consequences of it. This longer period of abnormal low rainfall years may become more frequent and unpredictable due to climate change (PBMC, 2013), with direct impacts on the level of reservoirs, inputs of nutrients from diffuse sources and its eutrophication (Moss et al., 2011; Dawadi and Ahmad, 2013; Umaña, 2014). On the other hand, water demand is increasing due to rapid regional development of irrigated agriculture and urbanization in recent years, in particular irrigate agriculture (Lacerda et al., 2008). Additionally, with the financial subsidies from the federal government to expand intensive fish-cage aquaculture to promote income and food safety for the region, the activity expanded by 20% per year, increasing environmental pressures within the reservoir proper (Oliveira et al., 2015). Therefore, sustainable use of the reservoirs services needs better modelling capability and scenarios construction in order to subsidize stockholder and decision makers. In this context, the present study aimed to spatially and temporally analyze and discuss the main physical and chemical water parameters of the Castanhão reservoir during a four-year period. Thereby contributing to the characterization of its trophic state in response to extended drought by applying the Trophic State Index approach to a comprehensive interpretation of the trophic variation in the reservoir.

METHODS

Study area

The study took place in the Castanhão reservoir (Latitude 5.50°S; Longitude 38.47°W) in the Middle Jaguaribe River watershed, which is located entirely within the semiarid region in the State of Ceará, NE Brazil (Fig. 1). The semiarid climate of Northeastern Brazil features peculiarities resulting from the behavior of its regulating weather systems marked by irregularities in rainfall across time and space, with annual rainfall means commonly ranging from 400 to 1000 mm and average of 756.5 mm during the past 80 years (FUNCEME, 2014). Rainfall occurs from January to June and is scarce from July to December. During the study period, the monthly precipitation varied from 0-181 mm. The years 2012, 2013 e 2014 were drier, with a mean annual rainfall of 302.3, 656.5 and 571.1 mm, respectively, significantly below the historical annual average. This extended drought period relates to a strong El Nino event.

The Castanhão reservoir flooded completely for the first time in 2004. The total storage capacity of the reservoir is 6.7 billion m3, and the normal operating capacity is 4.45 billion m3. The reservoir covers a flooded area of 325 km2 and is 48 km in length, with a depth exceeding 50 m in some areas (DNOCS, 2014). The categories of the World Commission on Dams (2000) classify the Castanhão as a large reservoir. The reservoir is a multiple use lake with major function as water storage for human and agriculture uses. It harbors the largest fish aquaculture facility in the state and is also use for recreational objectives, but not for energy production. The reservoir basin receives anthropogenic inputs of nutrients from different sources. Avelino (2015) developed a DPSIR analysis of the Castanhão reservoir and estimated, using established emission factors for these activities (Lacerda et al., 2008; Paula Filho et al., 2015) the total load of N and P from these sources. Annual direct emissions to the lake are mostly from fish farming, with an average annual production during the study period of about 18,000 tons of Nile Tilapia, and reached 519 t of N and 163 t of P). In a recent study, Molisani et al. (2015) calculated 12 tons of N being emitted per 6-month production cycle from a fish farm in the Castanhão reservoir that produces about 448 t of fish annually. Extrapolating this figure for the total average production in the entire reservoir one reaches about 964 tons of N. Those authors, however, did not indicate how emission factors were calculated and this number is probably an overestimation, but it is consistent with the annual emission of N estimated by Avelino (2015). Irrigated agriculture farms, covering 1254 ha and producing nearly 2 million tons mostly of fruit crops, are located around the lake, due to low-cost water availability. Runoff from these farms enters directly to the water body, totaling 198 t and 149 t of N and P, respectively. Outside the reservoir margins, agriculture is restricted to subsidence farming and contributes little to the total nutrient emissions to the reservoir. Urban wastes and waste waters from local villages, totaling about 7400 inhabitants, are partially treated. Sewage collection reaches about 70% of the total effluent in urban areas, whereas nearly 90% of solid wastes are adequately disposed (IBGE, 2010). Rural areas have no wastewater or sewage collection or treatment and no solid waste disposal systems (IPECE, 2011). However, population is sparse, reaching about 3200 inhabitants. As a result, annual emissions of N and P are relatively small, reaching 271 t of N and 73 t of P. Husbandry, mostly extensive cattle, although with a higher relative emission of about 1100 t for both nutrients, is located in the upper reaches of the basin and hardly reaches the reservoir, mostly due to semiarid conditions hampering an effective transport from soil runoff (COGERH, 2011; Avelino, 2015). In summary, fish farming is the largest direct source of anthropogenic derived N and P to the reservoir and contributes with about 52.5% and 42.3% of the total anthropogenic direct emissions of N and P, respectively.

Fig. 2 shows fluctuations in the volume of water stored in the reservoir during the studied period. These water level fluctuations are determined primarily by the dam system operation. However, during the sampling period, there was a drastic reduction in the stored volume due to the prolonged absence of rain in the reservoir’s basin, resulting in a volume decrease from about 85% of its total water storage capacity in November, 2011 to about 35% of its total volume in May, 2014.

Sampling

Sampling campaigns occurred in November 2011, March and August 2012, January and August 2013 and May 2014, covering 10 stations located throughout the reservoir (Fig. 1). The sampling dates were chosen according to the annual rainfall distribution. Campaigns in November, March and January, occurred in rainy seasons; whereas August and May in dry seasons. In this way 3 rain and 3 dry periods were covered.

In the first two campaigns, only four and six stations respectively were sampled due to logistical problems. Previous to water collection, physical and chemical variables were measured in situ in surface waters (0.5 to 1.0 m depth) as follows: dissolved oxygen (YSI 556 probe, YSI Inc., Yellow Springs; water temperature, turbidity and electrical conductivity (Compact-CTD model AST D687; JFE Advantech Co., Ltd., Nishinomiya); pH (Portable 826 pH-meter; Metrohm AG, Herisau); and transparency with a Secchi disk. In addition, CTD and dissolved oxygen depth profiles in each station described thermal structure of the water column in all campaigns.

Water samples collected from the subsurface (1.0 m) in Van Dorn bottles were analyzed for inorganic nutrients after filtering in the field lab through 47-mm-diameter AP40 glass fiber filters. Samples were immediately frozen for transport and later analysis. Unfiltered samples were used to determine the total phosphorus and total nitrogen concentrations. Variables were quantified in triplicate, with the final detection performed by visible spectrum spectrophotometry, according to: ammonia nitrogen (Koroleff, 1970), nitrate (Braga et al., 2015), nitrite (Bendschneider and Robinson, 1952), total nitrogen and total phosphorus (Valderrama, 1981) and soluble reactive phosphorus (Murphy and Rilley, 1962).

To quantify chlorophyll a concentrations, the samples were filtered immediately after sampling in a field lab through 47-mm-diameter AP40 glass fiber filters. Filters were stored frozen and transported to the main laboratory were quantification was obtained in a spectrophotometer according to the ISO 10260 (1992) protocol.

The trophic state of the reservoir was estimated using the trophic index (TSItsr) developed by Cunha et al. (2013). The index comprised the following equations that accounted for the concentrations of chlorophyll a (Chla; µg L–1) and total phosphorus (TP; µg L–1):

(eq. 1)
limno-2017-1-1433-e001.jpg
(eq. 2)
limno-2017-1-1433-e002.jpg
(eq. 3)
limno-2017-1-1433-e003.jpg

The defined limits for the trophic state classes were as follows: ultraoligotrophic: TSItsr ≤ 51.1; oligotrophic: TSItsr ≥ 51.2-53.1; mesotrophic: TSItsr ≥ 53.2-55.7; eutrophic: TSItsr ≥ 55.8-58.1; supereutrophic: TSItsr ≥ 58.2-59.0 and hypereutrophic: TSItsr≥ 59.1.

Principal Component Analysis (PCA), described the relationships between limnological variables and a cluster analysis to evaluate the longitudinal patterns in the reservoir and explore the similarities among sampling stations and define groups of associated ones. The analyses were performed using the STATISTICA 8.0 software package (StatSoft, Inc., Tulsa), assuming a significance level of α=0.05. For the cluster analysis, the data matrix was normalized in Z and Ward’s clustering method was used in combination with the Euclidean squared distances. The mean values of each parameter were subjected to means testing at a significance level of P≤0.05.

RESULTS

Tab. 1 shows mean, minimum and maximum values of the limnological parameters monitored in the reservoir during the sampling period and considering all sampling stations, whereas Supplementary Tab. 1 shows all individual results. Even considering that samples were collected in different periods of the day no significant change in temperature was observed between sampling and stations. However, thermal stratification of the water column during the period lead to a significant reduction in dissolved oxygen in the hypolimnion, reaching 0.07 to 2.62 mg L–1 (Ex.: Station 5) (Fig. 3 and Supplementary Fig. 1). Thermal stratification was particularly evident in the first two campaigns when the reservoir volume decreased very little, less than 5%. After that, thermal stratification disappeared, allowing wind induced mixing with the hypolimnion and a reduction in DO of the whole water column.

Average chlorophyll a concentrations differed significantly during the sampling period (P≤0.05) and increased progressively as the reservoir volume decreased (Tab. 1). Highest values occurred in May 2014 (12.4-37.9 µg L–1) in stations 1 and 2. High turbidity also occurred in these period and stations. Water transparency decreased 4-fold during the studied period. Secchi disk depth ranged from 4.0 m (station 9) in November 2011 to 1.0 m (station 1) in May 2014. Turbidity values remained below 4.2 NTU, with the highest values recorded in August 2013 (station 1) and in May 2014 (stations 1, 2, 3 and 4). PCA analysis strongly suggests volume decrease as the major factor influencing turbidity and the chlorophyll-a content. Average pH varied little throughout most of the sampling period from 7.5 to 8.2, and were not significantly different. However, significantly higher pH values were recorded in May 2014 (8.7), when reservoir volume was lowest, following the increased photosynthetic activity as well as conductivity. Average electrical conductivity at 25°C increased progressively as the reservoir volume decreased from the lowest values of 295 µS cm–1 in the beginning of the study in November 2011 and reached its maximum average value in May 2014, of 353 µS cm–1. The reservoir was well oxygenated at the surface (DO >5 mg L–1) throughout the study and DO concentrations did not differ significantly throughout the monitoring period (P>0.05).

Total phosphorus concentration were significantly higher (P<0.05) (>45 µg L–1) in May 2014, and lower and relatively similar (17 to 27 µg L–1) among the other sampling campaigns. Soluble reactive phosphorus concentrations (not shown in Tab. 1) were below the detection limit (<1 µg L–1) in most parts of the reservoir. Detectable concentrations occurred only in station 1 in January 2013 (2.5±0.8 µg L–1) and in station 5 in May 2014 (25.8±1.8 µg L–1). Nitrogen is the major nutrient component of the direct anthropogenic emissions to lake with a total annual emission of about 1,000 tons. Average total nitrogen concentrations did not differ significantly throughout the sampling period (397-598 µg L–1). In most sampling stations and campaigns, nitrate, ammonia nitrogen and nitrite concentrations were below the detection limit. The highest detected concentrations were 83.0±3.5 µg L–1 (station 9 in March 2012) for nitrate and 327.2±1.0 µg L–1 (station 5, in May 2014) for ammonia nitrogen. The availability of ammonia nitrogen in stations 5, 6, 7 and 8 in March 2012 might explain the algal blooms observed during that sampling campaign. Detectable nitrite concentrations (2.5±0.1 µg L–1) occurred only in January 2013 and May 2014.

Similarities among stations during the monitoring period were evaluated using a multivariate cluster analysis considering all variables presented in Tab. 1. Two situation appeared when data from each campaigns were successively added to the analysis. During most of the monitoring period up to August 2013, two groups of stations were clearly recognized (Fig. 4a). A Group 1 including the deeper stations closer to the dam (5, 6, 7, 8 and 9) with lower turbidity and N concentrations and a Group 2 including all shallower stations located along the fluvial axis of the reservoir (1, 2, 3, 4 and 10), differed by higher turbidity and N concentrations. When adding data from the last campaign, when reservoir level was at the lowest, the former Group 2 were subdivided, within two subgroups, (2A: stations 3, 4 and 10; 2B: stations 1 and 2); mostly differed by turbidity and Secchi disk depth (Fig. 4b). Group 1 included stations characterized by greater depths at the proximity to the dam. The area receives the largest effluents form fish farms and irrigated agriculture; and suggest an accumulation site for the nutrients drained from these activities, confirmed by the highest concentrations of P and N in stations of this group compared to Group 2, at least up to August 2013. In May 2014, Group 2B including stations 1 and 2, located at shallower depths and in regions further upstream of the dam, differed from the other by highest turbidity and lowest Secchi disk depth associated with the progressively shallower depths of the stations. In addition, Group 2B presented the highest chlorophyll a concentrations, which caused super-saturation of dissolved oxygen to occur (Tab. 2).

A principal component analysis of the data throughout the sampling periods illustrates the temporal changes in water quality of the reservoir through time (Fig. 5). Factor 1 (mostly chlorophyll a and total phosphorus) explain bout 40% of the variation in hydrochemistry, followed by factor 2 (total nitrogen) explaining about 14% of the variability (Tab. 3). The PCA clearly shows the impact of volume on chlorophyll a, total phosphorus and total nitrogen concentrations through time, which will then affect the reservoir’s trophic state index.

Fig. 6 shows the trophic state index variation in relation to the decreasing reservoir capacity. The trophic state index differed slightly between the spatial groups throughout the sampling period. The reservoir area encompassed by Group 1 remained oligotrophic from November 2011 to March 2012 and mesotrophic from August 2012 to May 2014, mostly due to a continuous increase in chlorophyll a. Phosphorus only showed a significant increase in May 2014 also contributing to the change in trophic state. In the area encompassed by Group 2A and Group 2B, waters remained oligotrophic throughout most of the monitored period, but with Group 2A, turning eutrophic and Group 2B, turning supereutrophic in May 2014, when the reservoir attained its minimum volume and a peak in chlorophyll a and phosphorus concentrations occurred.

DISCUSSION

The decrease in reservoir volume and the consequent lack of water renewal by rainfall, strong evaporation and the permanent influence of anthropogenic activities’ emissions, mostly from fish farming, irrigated agriculture and urban wastes, favored the concentration of salts, the accumulation of nutrients and increasing pH. Previous studies in lakes and reservoirs in semiarid regions showed high electrical conductivity that typically exceed 300 µS cm–1 and exhibit a negative correlation with water level (Bouvy et al., 1999; Eskinazi-Sant’Anna et al., 2007; Barbosa et al., 2012). In addition, wind forcing under decreasing volume produces rapid water mixing and sediment resuspension, changing water column properties and bringing up accumulated nutrients from the hypolimnion (Freire et al., 2009).

At the Castanhão reservoir decreasing volume is associated with increasing electrical conductivity, total phosphorus, turbidity, pH and to a lesser extent with total nitrogen, strengthen the interaction between reservoir volume and water chemistry. Thermal stratification of the water column led to oxycline formation along with a significant reduction in dissolved oxygen in the hypolimnion, but not at the surface. The increased concentration of organic detritus at the bottom of the water column and the subsequent decomposition of this material consume large amounts of oxygen and induce anoxia, as exemplified by station 5 (but see also Supplementary Tab. 1).

Chlorophyll a concentration is a good proxy of phytoplankton biomass and a valuable indicator of the trophic status of a water body. Salas and Martino (2001) used chlorophyll a concentrations to classify the trophic statuses of tropical lakes and suggested concentrations of 5 to 10 µg L–1 as indicative of a mesotrophic state. Average chlorophyll a concentrations in the Castanhão reservoir remained below 5 µg L–1 up to January 2013. Afterwards, concentrations in all stations exceeded 5 µg L–1, and exceeded 10 µg L–1, in May 2014 in stations 1, 2, 3, 4 e 9. Simultaneously, total phosphorus concentrations exceeded 35 µg L–1, within the range of values considered as the initial the values at which the deleterious environmental effects of eutrophication begin to appear (Dodds et al., 1998; Salas and Martino 2001).

The increasing nutrient concentrations observed after January 2013 might have resulted from the remobilization of accumulated nutrients in the hypolimnion and sediments consequent to the longer residence time of the reservoir. Previous observation on nutrient balance in the Castanhão reservoir showed a net retention of over 97% of the total fluvial input of N and P within the reservoir resulting in nutrient accumulation mostly in the hypolimnion (Molisani et al., 2013). With decreasing water depth, wind action can break the stability of the water column, thus favoring water column mixing processes and the resuspension of bottom material and the consequent increase in nutrient concentrations. These processes have been suggested to explain increasing nutrient concentrations in various reservoirs from the semiarid region during drought periods (Bouvy et al. 2003; Freire et al., 2009; Geraldes and George 2012).

Eventually, variations in the reservoir water level can result in changes in the physical and chemical structure of the system and consequently the dynamics of plankton communities, including zooplankton abundance (Geraldes and George, 2012), phytoplankton (Chalar, 2006), cyanobacteria blooms (Bouvy et al., 2003; Silva et al., 2013) and diversity (Costa et al., 2009). Naselli-Flores (2000) even suggested that water level fluctuations, influence phytoplankton composition more strongly than the nutrients levels. For example, a reduction in the reservoir volumes in Sicily during the spring because of the high demand from summer usage strongly affected the phytoplankton and nutrient dynamics, leading to eutrophication and the selection of cyanobacteria species that were harmful to human health (Naselli-Flores, 2003).

The trophic state index indicated a change in the trophic status of the Castanhão reservoir from oligotrophic to eutrophic during this extended drought period. In a short monitoring study between November 2006 and July 2007, Molisani et al. (2010) classified this reservoir as mesotrophic and detected anoxic events as well as the presence of cyanobacteria typical of mesotrophic/eutrophic environments, just after a two-year drought period. Chaves et al. (2013) also observed that the trophic state of a reservoir in the semiarid region relates to seasonal fluctuations in rainfall.

The observed deterioration of the TSI, however, was not evenly distributed in the entire reservoir. The three different groups identified and separated by their hydro chemical variables (Fig. 4) responded with different timing in developing eutrophic conditions. Group 1, including those deeper stations, receiving effluents from fish farming and irrigated agriculture, became mesotrophic and reached a eutrophic state earlier than the other groups. As water volume decreased further, Group 2A and 2B, characterized by much shallower water depth, became eutrophic. Further, when reservoir volume reached its minimum Group 2B, including the shallowest stations, became hypereutrophic. This suggest that a first process inducing changes in the trophic state was the concentration and accumulation of nutrients from anthropogenic activities, surrounding Group 1 region (Fig. 1) following the decrease in the reservoir volume. Afterwards, apart from continuing accumulation, and extra enrichment of nutrients was due to breaking thermal stratification, mixing of the water column (Groups 2A, 2B), and at least in the shallowest stations (Group 2B) sediment resuspension.

Molisani et al. (2015) reported, based on stable C and N distribution, current velocities and deposition rates of fecal pellets, that effluents from cage aquaculture (feed and feces) were not extensively reaching the water column and the bottom sediments of the cage area, on the contrary, effluents from farms are carried out to the deeper, central region of the reservoir. Hydrodynamic modeling of the reservoir (Molisani et al., 2015; Oliveira et al., 2015) showed surface wind-drive currents to move westwards from the center deeper part of the reservoir to shore, whereas bottom currents flush-out fish farm areas and incoming runoff eastwards to the deeper regions (Group 1 region), where stratification exists during high volume periods (e.g., in station 6 in November, 2011 - Supplementary Tab. 1), but is broken when the reservoir volume decreases. This contributes to the increase in nutrients in the hypolimnion of the reservoir. With the intense drought in recent years, this nutrient-rich hypolimnion waters mix with the entire column, increasing eutrophication, as showed by the TSI of the latter periods of the study.

Unfortunately, by comparing the emission of a single farm (67 t production) with the total N contribution from the watershed, Molisani et al. (2015) reached the conclusion that fish farming contributes relatively little to the total N inputs to the reservoir. When extrapolated to the entire fish production of the reservoir, however, this activity is the largest direct contributor of N and P to the reservoir. The hydrodynamics of the lake, cleanses the farm areas by exporting the majority of its effluents to deeper parts of the lake, where nutrients are accumulated under the thermocline establish when the reservoir is in its maximum volume, reaching over 54 m in depth in the central part. Under the intense drought verified during the study period, however, these accumulated nutrients are release to the water column following disruption of the thermocline due to shallower depths and triggering eutrophication.

The results obtained in this study raise concern on the sustainability of fish farming in semiarid reservoirs under a scenario of global climate change. The multiple use of reservoir waters in regions of hot and cold semiarid climate is presently constrained by global climate change (Jeppesen et al., 2015). Most global modeling showed these regions to be strongly affected by reduced precipitation and enhanced evaporation, as predicted for northeastern Brazil (Oyama and Nobre, 2003; Marengo et al., 2009). Actually, reduction in the annual rainfall has already been recorded for the region. Moncunill (2006) has estimated an average decrease of 5.6 mm.yr–1 of rainfall since the 1960’s over the Ceará State, whereas Godoy and Lacerda (2015) suggested that reduction of annual rainfall, in particular during the dry season, drastically reduced the Jaguaribe River flow over the middle Jaguaribe River Basin, where the Castanhão reservoir is located. Climate change is intensifying the symptoms of eutrophication in freshwaters, mostly affecting diffuse sources of nutrients to lakes and reservoirs (Jeppesen et al., 2015) increasing nutrient concentrations due to evaporation in reduced water volumes becomes an additional factor strengthening eutrophication. At the Castanhão reservoir, average decrease of rainfall will extend the duration of drought periods, such as the one covered by this study. This will result in volume decrease, breaking of the thermal stratification and triggering eutrophication. Falling reservoir levels apart from directly concentrating the nutrients already present in the water column, released by in-lake activities such as aquaculture, will expose sediments from margins and shallow areas receiving up-basin effluents from agriculture and husbandry and urban areas, to mineralization and increasing nutrient release to the lake.

CONCLUSIONS

The prolonged period of drought caused a dramatic decrease of the reservoir volume, which in turn induced changes in the physical and chemical structure of the system that resulted in eutrophication. The formation of spatially different compartments regarding the major physical and chemical parameters, however, suggests that the intensity of the eutrophication process in not uniformly distributed in the reservoir.

When normal rainfall periods keep elevated reservoir volumes, specific water current patterns wash out effluents from fish farms and irrigated agriculture transferring them to deeper layers, where nutrients accumulate in a relatively segregated hypolimnion. When the reservoir becomes shallower, breaking of the thermal stratification and sediment remobilization accelerates nutrient enrichment of the water column and triggers eutrophication.

The scenario of global climate change has to be taken into consideration to the decision-making processes aiming to ensure the sustainability of the multiple uses of the Castanhão reservoir and the preservation of high water quality for the noblest uses. Reducing rainfall and increasing frequency of extreme events, such as extended droughts and strong downpours, during the rainy seasons, of semiarid climate, will intensify the processes described for the Castanhão reservoir and will require re-dimensioning of fish aquaculture and irrigated agriculture, as well as better treatment of urban wastes.

ACKNOWLEDGMENTS

We thank the National Council for Scientific and Technological Development (CNPq-Brazil Proc. No. 573.601/2008-9) and the Ceará Foundation for Scientific and Technological Development (FUNCAP, Proc. Nos. 561.282/2010 and 120.100/2011) for financial support and grants to the authors.

REFERENCES

1 

IIF Avelino, 2015. [Avaliação ambiental do Reservatório Castanhão (CE) utilizando o método DPSIR].[Thesis in Portuguese]. Universidade Federal do Ceará, Fortaleza.

2 

JEL Barbosa, ESF Medeiros, J Brasil, RS Cordeiro, MCB Crispim, GHG Silva, 2012. Aquatic Systems in Semi-arid Brazil: limnology and management. Acta Limnol. Brasil. 24:103-118.

3 

K Bendschneider, RJ Robinson, 1952. A new spectrophotometric method for the determination of nitrite in seawater. J. Mar. Res. 11:87-96.

4 

LAV Bezerra, WD Paulino, DS Garcez, H Becker, JI Sánchez-Botero, 2014. Características limnológicas em um reservatório do semiárido brasileiro com cultivo intensivo de Orechromis niloticus (Linnaeus, 1758). Acta Limnol. Brasil. 26:47-59.

5 

M Bouvy, R Arfi, S Oliveira, M Marinho, B Beker, 1999. Dynamics of a toxic cyanobacterial bloom (Cylindrospermopsis raciborskii) in a shallow reservoir in the semi-arid region of Northeast Brazil. Aquat. Microb. Ecol. 20:285-297.

6 

M Bouvy, SM Nascimento, RJR Molica, A Ferreira, V Huszar, SMFO Azevedo, 2003. Limnological features in Tapacurá reservoir (northeast Brazil) during a severe drought. Hydrobiol. 493:115-130.

7 

GG Braga, V Becker, JNP Oliveira, Junior JR Mendonça, AFM Bezerra, LM Torres, AMF Galvão, A Mattos, 2015. Influence of extended drought on water quality in tropical reservoirs in a semiarid region. Acta Limnol. Brasil. 27:15-23.

8 

G Chalar, 2006. [Dinámica de la eutrofización a diferentes escalas temporales: embalse Salto Grande (Argentina-Uruguay)], p. 87-101. In: J.G. Tundisi, T. Matsumura-Tundisi, C.S. Galli (eds.), [Causas, consequências e tecnologias de gerenciamento e controle].[Book in Portuguese]. Instituto Internacional de Ecologia, São Paulo.

9 

FIB Chaves, PF Lima, RC Leitão, WD Paulino, ST Santaella, 2013. Influence of rainfall on the trophic status of a Brazilian semiarid reservoir. Acta Scientiarum 35:505-511.

10 

COGERH (Companhia de Gestão dos Recursos Hídricos), 2011. [Inventario Ambiental Açude Castanhão].[Book in Portuguese]. Geosolos Consultoria, Projetos e Serviços Ltda, Fortaleza: 142 pp.

11 

IAS Costa, SRS Cunha, RF Panosso, MF Araújo, JLS Melo, EM Eskinazi-Sant’anna, 2009. [Dinâmica de cianobactérias em reservatórios eutróficos do semiárido do Rio Grande do Norte].[Article in Portuguese]. Oecol. Brasil. 13:382-401.

12 

DGF Cunha, MC Calijuri, MC Lamparelli, 2013. A trophic state index for tropical/subtropical reservoirs (TSItsr). Ecol. Engin. 60:126-134.

13 

S Dawadi, S Ahmad, 2013. Evaluating the impact of demand-side management on water resources under changing climatic conditions and increasing population. J. Environ. Manag. 114:261-275.

14 

(Departamento Nacional de Obras Contra as Secas) DNOCS, 2014. [Açude Castanhão].[In Portuguese]. Accessed: 10 June 2014. Available from: http://www.dnocs.gov.br

15 

WK Dodds, JR Jones, EB Welch, 1998. Suggested classification of stream trophic state, distributions of temperate of stream trophic state: chlorophyll, total nitrogen, and phosphorus. Wat. Res. 32:1455-1462.

16 

EM Eskinazi-Sant’anna, R Menezes, IS Costa, RF Panosso, MF Araújo, JL Attayde, 2007. [Composição da comunidade zooplanctônica em reservatórios eutróficos do semiárido do Rio Grande do Norte].[Article in Portuguese]. Oecol. Brasil. 11:410-421.

17 

RHF Freire, MC Calijuri, ST Santaella, 2009. Longitudinal patterns and variations in water quality in a reservoir in the semiarid region of NE Brazil: responses to hydrological and climatic changes. Acta Limnol. Brasil. 21:251-262.

18 

(Fundação Cearense de Meteorologia e Recursos Hídricos) FUNCEME, 2014. [Chuvas: Séries Históricas].[In Portuguese]. Accessed in: October 2014. Available from: http://www.funceme.br/index.php/areas/tempo/download-de-series-historicas

19 

AM Geraldes, C George, 2012. Limnological variations of a deep reservoir in periods with distinct rainfall patters. Acta Limnol. Brasil. 24:417-426.

20 

MDP Godoy, LD Lacerda, 2014. River-island response to land-use change within the Jaguaribe River, Brazil. J. Coast. Res. 30:399-410.

21 

IBGE (Instituto Brasileiro de Geografia e Estatística), 2010. [Censo 2010].[In Portuguese]. Assessed in: June 2015. Available from: http://www.censo2010.ibge.gov.br/pt/

22 

IPECE (Instituto de Pesquisa e Estatística Econômica do Ceará), 2011. Assessed in June, 2015. Available from http://www2.ipece.ce.gov.br/publicacoes/ceara_em_numeros/2011/infra/index.htm

23 

ISO, 1992. Water quality measurement of biochemical parameters spectrophotometric determination of chlorophyll-a concentration. ISO Norm 10260. International Organization for Standardization, Geneva.

24 

E Jeppesen, S Brucet, L Naselli-Flores, E Papastergiadou, K Stefanidis, T Noges, JL Attayde, T Zohary, J Coppens, T Bucak, RF Menezes, FRS Freitas, M Kernan, M Søndergaard, M Beklioğlu, 2015. Ecological impacts of global warming and water abstraction on lakes and reservoirs due to changes in water level and related changes in salinity. Hydrobiologia 750:201-227.

25 

F Koroleff, 1970. Revised version of direct determination ammonia of natural waters as indophenol blue. Information on techniques and methods for sea water analysis. ICES Interlab. Rep. 3:19-22.

26 

LD Lacerda, MM Molisani, D Sena, LP Maia, 2008. Estimating the importance of natural and anthropogenic sources on N and P emission to estuaries along the Ceará state coast NE Brazil. Environ. Monit. Assess. 141:149-164.

27 

FB Lopes, EM Andrade, ACM Meireles, H Becker, AA Batista, 2014. Assessment of the water quality in a large reservoir in semiarid region of Brazil. R. Bras. Eng. Agríc. Ambiental 18:437-445.

28 

JA Marengo, R Jones, MLM Alvesa, MC Valverdea, 2009. Future change of temperature and precipitation extremes in South America as derived from the PRECIS regional climate modeling system. Inter. J. Climatol. 29: 2241-2255.

29 

MM Molisani, HS Barroso, H Becker, MOP Moreira, CAG Hijo, TM Monte, GH Vasconcellos, 2010. Trophic state, phytoplankton assemblages and limnological diagnosis of the Castanhão Reservoir, CE, Brazil. Acta Limnol. Brasil. 22:1-12.

30 

MM Molisani, H Becker, HS Barroso, CAG Hijo, TM Monte, GH Vasconcellos, LD Lacerda, 2013. The influence of Castanhão reservoir on nutriente and suspended matter transport during rainy season in the ephemeral Jaguaribe River (CE, Brazil). Brazil. J. Biol. 73:115-123.

31 

MM Molisani, TM Monte, GH Vasconcellos, HS Barroso, MOP Moreira, H Becker, CE Rezende, MAL Franco, EGG Farias, PB Camargo, 2015. Relative effects of nutrient emission from intensive cage aquaculture on the semiarid reservoir water quality. Environ. Monit. Assess. 187:707-721.

32 

DF Moncunill, 2006. The rainfall trend over Ceará and its implications, p. 315-323. Proc. 8th Int. Conf. Southern Hemisphere Meteorology and Oceanography, Foz do Iguaçu.

33 

B Moss, S Kosten, M Meerhoff, RW Battarbee, E Jeppesen, N Mazzeo, K Havens, G Lacerot, Z Liu, L De Meester, H Paerl, M Scheffer, 2011. Allied attack: climate change and eutrophication. Inland Waters 1:101-105.

34 

J Murphy, JP Riley, 1962. A modified single solution method for the determination of phosphate in natural waters. Anal. Chim. Acta 27:31-36.

35 

L Naselli-Flores, 2000. Phytoplankton assemblages in twenty-one Sicilian reservoirs: relationship between species composition and environmental factors. Hydrobiologia 424:1-11.

36 

L Naselli-Flores, 2003. Man-made lakes in Mediterranean semiarid climate: the strange case of Dr. Deep Lake and Mr. Shallow Lake. Hydrobiologia 506-509:13-21.

37 

KF Oliveira, LD Lacerda, TF Peres, MF Bezerra, FJS Dias, 2015. Emission factor and balance of mercury in fish farms in an artificial reservoir in NE-Brazil. Environ. Sci. Pollut. Res. 22:18278-18287.

38 

MD Oyama, CA Nobre, 2003. A new climate–vegetation equilibrium state for tropical South America. Geophys. Res. Lett. 30:2199-2203.

39 

Filho FJ Paula, RV Marins, LD Lacerda, 2015. Natural and anthropogenic emissions of N and P to the Parnaíba River Delta in NE Brazil. Estuar. Coast. Shelf Sci. 166:34-44.

40 

(Painel Brasileiro de Mudanças Climáticas) PBMC, 2013. [Base científica das mudanças climáticas. Contribuição do Grupo de Trabalho 1ao Primeiro de Avaliação Nacional do Painel Brasileiro de Mudanças Climáticas. Sumário Executivo do GT2].[in Portuguese]. Painel Brasileiro de Mudanças Climáticas, Rio de Janeiro, Brasil.

41 

HJ Salas, P Martino, 2001. Metodologías simplificadas para La evaluación de eutroficacion em lagos cálidos tropicales. Organización Panamericana de La Salud (OPS) – División de Salud y Ambiente. Centro Panamericano de Ingenieria Sanitaria y Ciências Del Ambiente (CEPIS).

42 

JCN Santos, EM Andrade, Neto JR Araújo, AM Meireles, HAQ Palácio, 2014. Land use and trophic state dynamics in a tropical semi-arid reservoir. Rev. Ciênc. Agr. 45:35-44.

43 

LP Silva, JA Costa BGB Santos, MB Bezerra, ATF Silva, RV Marins, 2013. Caracterização da comunidade fitoplanctônica no açude padre Cícero (Castanhão), CE. An. XIV Congr. Bras. Limnol. Associação Brasileira de Limnologia, Bonito, Mato Grosso do Sul, p:123.

44 

G Umaña, 2014. Ten years of limnological monitoring of a modified natural lake in the tropics: Cote Lake, Costa Rica. Rev. Biol. Trop. 62:567-578.

45 

JC Valderrama, 1981. The simultaneous analysis of total nitrogen and total phosphorus in natural waters. Mar Chem. 10:109-122.

46 

Commission on Dams World, 2000. Dams and development: a new framework for decision-making. Earthscam Publications, London: 404 pp.

Fig. 1.

Study area and location of sampling stations in the Castanhão reservoir, NE Brazil.

limno-2017-1-1433-g001.jpg
Fig. 2.

Fluctuation of the water volume of the Castanhão reservoir, NE Brazil between 2011 and 2014 (DNOCS, 2014).

limno-2017-1-1433-g002.jpg
Fig. 3.

Vertical profile of temperature and dissolved oxygen concentrations at station 5 during the monitoring period in the Castanhão reservoir, NE Brazil.

limno-2017-1-1433-g003.jpg
Fig. 4.

Cluster analysis of the sampling stations in the Castanhão reservoir. Including all the monitoring period and variables shown in Tab. 1. (a) From November 2011 to August 2013. (b) May 2014.

limno-2017-1-1433-g004.jpg
Fig. 5.

Principal component analyses (PCA) of the Castanhão hydrochemistry during the study period.

limno-2017-1-1433-g005.jpg
Fig. 6.

Trophic state index and its relationship with reservoir capacity fluctuations in the Castanhão reservoir in NE Brazil, along the monitoring period.

limno-2017-1-1433-g006.jpg
Tab. 1.

Mean, maximum and minimum values for the variables monitored in surface waters of the Castanhão reservoir, NE Brazil, during the monitoring period.

Variables Statistics November 2011
(n=6)
March 2012
(n=6)
August 2012
(n=10)
January 2013
(n=10)
August 2013
(n=10)
May 2014
(n=10)
Temperature (°C) Mean±SD 28.9±0.3a 30.2±0.6b 27.8±0.3c 28.8±0.3a 29.1±0.5a 30.4±1.0b
Min-Max 28.7-29.3 29.5-31.0 27.3-28.3 28.3-29.3 28.1-29.8 29.7-32.6
Secchi (m) Mean±SD 3.4±0.6ac 2.3±0.4b 3.4±0.3c 2.8±0.4ab 2.5±0.4b 2.3±0.8b
Min-Max 2.7-4.0 1.7-2.9 3.0-3.8 2.0-3.2 2.0-3.2 1.0-3.0
Turbidity (NTU) Mean±SD 1.3±0.1ab 1.1±0.1a 1.3±0.2ab 1.5±0.3ab 1.5±0.6ab 1.9±1.0b
Min-Max 1.1-1.5 1.0-1.3 1.1-1.6 1.1-2.2 1.1-2.9 1.2-4.2
Conductivity (µS cm–1) Mean±SD 295±2a 313±4ab 315±8b 343±7c 347±9c 353±20c
Min-Max 294-298 308-317 307-328 334-354 340-372 322-372
Dissolved oxygen (mg L–1) Mean±SD 7.0±0.1a 6.6±0.2a 7.0±0.3a 6.8±0.7a 6.6±0.4a 7.7±1.8a
Min-Max 6.9-7.1 6.4-6.8 6.6-7.5 6.1-8.4 6.0-7.4 6.5-11.7
pH Mean±SD 8.1±0.1abc 8.1±0.2abc 7.7±0.3ab 7.5±0.3a 8.2±0.4bc 8.7±0.7c
Min-Max 8.0-8.2 7.9-8.4 7.1-8.0 7.2-8.1 7.6-8.8 7.1-9.9
Chlorophyll a (µg L–1) Mean±SD 2.6±1.2a 3.7±1.2ab 4.0±1.2ab 4.1±1.3ab 5.2±2.5ab 14.8±11.4b
Min-Max 1.3-4.0 2.3-5.3 2.1-6.1 2.4-6.4 2.6-9.8 2.7-37.9
Tot-P (µg L–1) Mean±SD 27.2±8.9ab 17.2±7.1a 22.3±2.6a 30.1±10.7ab 20.6±8.7a 49.1±18.8b
Min-Max 22.2-40.5 9.5-28.0 17.8-25.3 15.8-48.5 7.6-33.7 23.6-77.2
Tot-N (µg L–1) Mean±SD 510±133a 397±167a 489±190a 405±170a 481±158a 598±312a
Min-Max 349-643 216-611 185-798 170-714 210-713 99-1.210

[i] Different letters are significant different mean values (HSD Tukey test, with P≤0.05).

Tab. 2.

Mean, maximum and minimum values of the studied variables in surface water observed in the Castanhão reservoir in the different groups depicted by the cluster analysis in May 2014.

Variables Statistics Group 1 (n=29) Group 2A (n=13) Group 2B (n=10)
Temperature (°C) Mean±SD 28.9±1.0a 28.6±0.7a 29.6±1.7a
Min-Max 27.4-31.0 27.3-29.9 28.1-32.6
Secchi (m) Mean±SD 3.0±0.6a 2.9±0.5ab 2.2±0.8b
Min-Max 1.7-4.0 2.0-3.6 1.0-3.4
Turbidity (NTU) Mean±SD 1.2±0.1a 1.5±0.3b 2.2±1.0b
Min-Max 1.0-1.5 1.2-2.2 1.5-4.2
Conductivity (µS cm–1) Mean±SD 315±16a 332±17a 336±12a
Min-Max 294-344 307-354 322-354
DO (mg L–1) Mean±SD 6.9±0.5a 6.8±0.3a 7.8±2.0a
Min-Max 6.3-8.4 6.1-7.1 6.2-11.7
pH Mean±SD 7.9±0.3a 7.6±0.3a 8.4±0.8a
Min-Max 7.2-8.4 7.1-8.0 7.4-9.9
Chlorophyll a (µg L–1) Mean±SD 4.1±1.4a 3.2±0.7a 10.6±14.7a
Min-Max 1.3-6.4 2.1-4.5 2.1-37.9
Tot-P (µg L–1) Mean±SD 25.4±10.6a 23.1±4.7a 31.5±21.3a
Min-Max 9.5-48.5 12.0-30.0 7.6-77.2
Tot-N (µg L–1) Mean±SD 492.0±162.2a 364.7±160.5a 512.5±176.7a
Min-Max 264.3-797.9 169.8-706.0 259.3-746.0

[i] Different letters are significant different mean values (P≤0.05).

Tab. 3.

Principal component analyses (PCA) of the Castanhão hydrochemistry. Underlined numbers are significant at P<0.05 (n=52).

Variable Factor 1 Factor 2
Temperature -0.674890 -0.072759
Turbidity -0.676510 -0.206673
Conductivity -0.595970 0.660054
Dissolved oxygen -0.086678 0.361447
pH -0.623495 -0.309359
Secchi 0.639051 -0.236606
Chlorophyll a -0.825686 -0.266973
Total phosphorus -0.724864 -0.240770
Total nitrogen -0.362774 -0.471254
Volume 0.765747 -0.485564
Expl. Var 3.998125 1.353106
Prp. totl 0.399813 0.135311
Abstract views:
1810

Views:
PDF
678
Supplementary
131
HTML
883

Article Metrics

Metrics Loading ...

Metrics powered by PLOS ALM


Copyright (c) 2016 Janaina A. Santos, Rozane V. Marins, José E. Aguiar, Guillermo Challar, Francisco A.T.F. Silva, Luiz D. Lacerda

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
 
© PAGEPress 2008-2017     -     PAGEPress is a registered trademark property of PAGEPress srl, Italy.     -     VAT: IT02125780185