Acessibilidade / Reportar erro

Potential of the retention capability of a Neotropical reservoir (São Paulo State, Brazil)

Potencial da capacidade de retenção de um reservatório Neotropical (Estado de São Paulo, Brasil)

Abstract:

Aim

Man-made reservoirs lead to several changes in their downstream rivers that depend on the hydraulic characteristics of the reservoirs. However, their multiple uses can also provide facilities that influence the ecosystem services that they provide. This study addressed the potential ability of a Neotropical reservoir to trap chemical species aiming to assess the role of these ecosystems to mitigate pollution.

Methods

Retention capability modeling was examined for a small subtropical reservoir with high hydraulic retention time (> 100 days). The temporal ranges of 9 physical and chemical water variables over a five-year period were used to calculate the mass balance and to determine the retentive capability (alpha parameter) of the Itupararanga Reservoir (São Paulo State, Brazil). To explain the long-term mass balance of these variables, it was assumed that the reservoir is a completely mixed system with a step input.

Results

The highest values of parametrized alpha (high retention capability) occurred in wet months (up to 500 mm) for all variables. High reaction rate constants (k) and low hydraulic flushing suggested that sink processes prevail over the export ones, mainly for total phosphorus. The rainfall pattern showed minor importance for trapping elements.

Conclusions

In the Neotropics, hydraulic characteristics of the ecosystem (e.g., low area:volume ratio) are a tool for pollution management in man-made reservoirs.

Keywords
ecosystem services; mass balance; eutrophication; water resources management; subtropical reservoir

Resumo:

Objetivo

Os reservatórios artificiais desencadeiam várias mudanças em seus rios a jusante que dependem das características hidráulicas dos reservatórios. No entanto, seus múltiplos usos também podem fornecer recursos que influenciam os serviços ecossistêmicos que eles fornecem. Este estudo abordou a capacidade potencial de um reservatório Neotropical em capturar espécies químicas, com o objetivo de avaliar o papel desses ecossistemas na mitigação da poluição.

Métodos

A modelagem da capacidade de retenção foi examinada para um pequeno reservatório subtropical com alto tempo de retenção hidráulica (> 100 dias). As variações temporais de 9 variáveis físicas e químicas da água ao longo de um período de cinco anos foram usados para calcular o balanço de massa e determinar a capacidade retentiva (parâmetro alfa) do reservatório Itupararanga (estado de São Paulo, Brasil). Para explicar o balanço de massa de longo prazo dessas variáveis, foi assumido que o reservatório seja um sistema completamente misturado com uma entrada em degrau.

Resultados

Os maiores valores de assimilação (i.e., alta capacidade de retenção) ocorreram nos meses úmidos (até 500 mm) para todas as variáveis. Altos coeficientes de reação (k) e baixa descarga hidráulica sugerem que os processos de sumidouro prevalecem sobre os de exportação, principalmente para fósforo total. O padrão de precipitação mostrou menor importância para o aprisionamento dos elementos.

Conclusões

Nos Neotrópicos, as características hidráulicas do ecossistema (e.g., baixa relação área: volume) são uma ferramenta para o gerenciamento da poluição em reservatórios artificiais.

Palavras-chave
serviços ecossistêmicos; balanço de massa; eutrofização; gestão de recursos hídricos; reservatório subtropical

1. Introduction

Rivers are ecosystems that provide a range of essential services (drinking water, nutrient cycling, water pollution management, fisheries and aquaculture, transport of dissolved and suspended materials, habitats for a range of natural species) that support life on Earth (Van Cappellen & Maavara, 2016Van Cappellen, P., & Maavara, T., 2016. Rivers in the anthropocene: global scale modifications of riverine nutrient fluxes by damming. Ecohydrol. Hydrobiol. 16(2), 106-111. http://dx.doi.org/10.1016/j.ecohyd.2016.04.001.
http://dx.doi.org/10.1016/j.ecohyd.2016....
). Damming of natural rivers is an anthropogenic activity that represents one of the main disturbances for rivers and for their riparian areas, leading to both local and global consequences (Van Cappellen & Maavara, 2016Van Cappellen, P., & Maavara, T., 2016. Rivers in the anthropocene: global scale modifications of riverine nutrient fluxes by damming. Ecohydrol. Hydrobiol. 16(2), 106-111. http://dx.doi.org/10.1016/j.ecohyd.2016.04.001.
http://dx.doi.org/10.1016/j.ecohyd.2016....
; Akbarzadeh et al., 2019Akbarzadeh, Z., Maavara, T., Slowinski, S., & Van Cappellen, P., 2019. Effects of damming on river nitrogen fluxes: a global analysis. Global Biogeochem. Cycles 33(11), 1339-1357. http://dx.doi.org/10.1029/2019GB006222.
http://dx.doi.org/10.1029/2019GB006222...
). Dams are mainly constructed for economic purposes, such as water irrigation, water supply, electricity generation and flood control. However, they can also act as a retention basin for pollutants and this has implications for eutrophication and/or for retention of biogenic elements that may cause concerns for water supply, navigation, biodiversity conservation and aesthetic issues (Bartoszek & Koszelnik, 2016Bartoszek, L., & Koszelnik, P., 2016. The qualitative and quantitative analysis of the coupled C, N, P and Si retention in complex of water reservoirs. Springerplus 5(1), 1157. PMid:27504255. http://dx.doi.org/10.1186/s40064-016-2836-7.
http://dx.doi.org/10.1186/s40064-016-283...
; Cunha-Santino et al., 2017Cunha-Santino, M.B., Fushita, A.T., & Bianchini Junior, I., 2017. A modeling approach for a cascade of reservoirs in the Juquiá-Guaçu River (Atlantic Forest, Brazil). Ecol. Modell. 356, 48-58. http://dx.doi.org/10.1016/j.ecolmodel.2017.04.008.
http://dx.doi.org/10.1016/j.ecolmodel.20...
; Qin et al., 2020Qin, L., Lei, P., Lei, Q., Liu, H., Li, X., Zhang, H., & Lindsey, S., 2020. Evaluating the effect of dam construction on the phosphorus fractions in sediments in a reservoir of drinking water source, China. Environ. Monit. Assess. 192(2), 99. PMid:31912244. http://dx.doi.org/10.1007/s10661-019-8053-4.
http://dx.doi.org/10.1007/s10661-019-805...
).

Hydraulic management (e.g., water level control) of reservoirs affects water quality due to the retention of elements derived from the drainage area and from point source pollution. Rainfall is essential for maintenance of the maximum water level and also influences the hydraulic retention time (HRT). HRT is an important variable that may control both sedimentation and biogeochemical reactions occurring in the water column. The retention capability of reservoirs is also affected by thermal stratification and hydrological conditions (Kerimoglu & Rinke, 2013Kerimoglu, O., & Rinke, K., 2013. Stratification dynamics in a shallow reservoir under different hydro-meteorological scenarios and operational strategies. Water Resour. Res. 49(11), 7518-7527. http://dx.doi.org/10.1002/2013WR013520.
http://dx.doi.org/10.1002/2013WR013520...
). For instance, the occurrence of anoxic conditions in the hypolimnion (Cole & Hannan, 1990Cole, T.M., & Hannan, H.H., 1990. Dissolved oxygen dynamics. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 71-107.) can also induce the release of soluble elements from the sediments (Søndergaard et al., 2003Søndergaard, M., Jensen, J.P., & Jeppesen, E., 2003. Role of sediment and internal loading of phosphorus in shallow lakes. Hydrobiologia 506-509(1-3), 135-145. http://dx.doi.org/10.1023/B:HYDR.0000008611.12704.dd.
http://dx.doi.org/10.1023/B:HYDR.0000008...
).

Rainfall water inputs modify the HRT, and consequently the downstream flow of elements, hence altering reservoir retention capability causing changes in downstream ecosystems due to exportation of elements. In addition, if water input and volume change, then water column circulation is likely to affect the thermal stratification of the reservoir, with resulting consequences for a range of chemical and biological processes within the reservoir (Armengol et al., 1999Armengol, J., Garcia, J.C., Comerma, M., Romero, M., Dolz, J., Roura, M., Han, B.H., Vidal, A., & Šimek, K., 1999. Longitudinal processes in canyon type reservoirs: the case of Sau (N.E. Spain). In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 313-345.; Teixeira et al., 2022Teixeira, A.P., Bortolini, J.C., & Carneiro, F.M., 2022. Taxonomic and functional spatial distribution model of phytoplankton in tropical cascading reservoirs. Rev. Bras. Bot. Braz. J. Bot. 45(2), 791-805. http://dx.doi.org/10.1007/s40415-022-00810-7.
http://dx.doi.org/10.1007/s40415-022-008...
). On the one hand, these changes can minimize the retention of elements, mimicking the original condition of the river (i.e., before reservoir construction). On the other hand, depending on water quality status, those conditions can increase the export of pollutants and toxic elements to downstream aquatic environments (Straškraba & Tundisi, 1999Straškraba, M., & Tundisi, J.G., 1999. Reservoir ecosystem functioning: theory and application. In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 565-583.).

Understanding the role of rainfall pattern on the retention capability of reservoirs is useful for water management purposes (Straškraba, 1999Straškraba, M., 1999. Retention time as a key variable of reservoir limnology. In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 385-410.). Most studies about river impoundment have focused on the loss of biodiversity, invasive species, landscape changes, water resources and social impacts (Hoeinghaus et al., 2009Hoeinghaus, D.J., Agostinho, A.A., Gomes, L.C., Pelicice, F.M., Okada, E.K., Latini, J.D., Kashiwaqui, E.A.L., & Winemiller, K.O., 2009. Effects of river impoundment on ecosystem services of large tropical rivers: embodied energy and market value of artisanal fisheries. Conserv. Biol. 23(5), 1222-1231. PMid:19459891. http://dx.doi.org/10.1111/j.1523-1739.2009.01248.x.
http://dx.doi.org/10.1111/j.1523-1739.20...
). Although several studies have examined the effects of dams on maintaining water quality and the amount of water they provide, relatively few studies have examined these effects in a quantitative, mathematically standardized way (e.g., Cunha-Santino et al., 2017Cunha-Santino, M.B., Fushita, A.T., & Bianchini Junior, I., 2017. A modeling approach for a cascade of reservoirs in the Juquiá-Guaçu River (Atlantic Forest, Brazil). Ecol. Modell. 356, 48-58. http://dx.doi.org/10.1016/j.ecolmodel.2017.04.008.
http://dx.doi.org/10.1016/j.ecolmodel.20...
; Bianchini Junior et al., 2019Bianchini Junior, I., Fushita, Â.T., & Cunha-Santino, M.B., 2019. Evaluating the retention capacity of a new subtropical run-of-river reservoir. Environ. Monit. Assess. 191(3), 161. PMid:30771013. http://dx.doi.org/10.1007/s10661-019-7295-5.
http://dx.doi.org/10.1007/s10661-019-729...
). This is of particular importance in tropical reservoirs, because in addition to sedimentation, material retention is closely linked to biological processes (e.g., autochthonous organic matter production), where high temperatures tend to support high values of chemical and biochemical reaction rates. On the contrary, the mineralization rates (also affected by temperature) tend to minimize the retention efficiency (Kennedy & Walker, 1990Kennedy, R.H., & Walker, W.W., 1990. Reservoir nutrient dynamics. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 109-131.; Davidson & Janssens, 2006Davidson, E.A., & Janssens, I.A., 2006. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440(7081), 165-173. PMid:16525463. http://dx.doi.org/10.1038/nature04514.
http://dx.doi.org/10.1038/nature04514...
).

In this study, we address the potential of the retention capability in a Neotropical reservoir in which thermal stratification is usually present during the rainy season. We hypothesize that periods of low precipitation (i.e., dry season: spring and winter) increase the HRT, and therefore, the retention of sediment particles and other materials, making the effect of temperature secondary, through indigenous production, on the retention of the elements. This has implications for internal loading of the reservoir, and consequently for the water supply downstream. We chose the Itupararanga Reservoir, located in subtropical climate (Brazil), to adjust a retention capability model because this system is small, with electricity generation roughly proportional to the amount of rainfall. However, the model can be applied for any ecosystem with steady-state flow and with HRT high enough to reach equilibrium.

2. Materials and Methods

2.1. Study area

Itupararanga Reservoir (23°36'43.35”S;47°23'49.87”W) located in the Middle Tietê River Basin in São Paulo State, Brazil is included in Water Resources Management Unit (UGRHI) 10, which has an area of 1 178 922.17 km2. According to the classification proposed by Straškraba (1999)Straškraba, M., 1999. Retention time as a key variable of reservoir limnology. In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 385-410., the reservoir is small (area ranging from 1 to 100 km2) and Class B (15 days <HRT <1 year). As stated by the total phosphorus results obtained in the period of this survey (2007-2012), this reservoir is predominantly mesotrophic (sensuVollenweider & Kerekes, 1982Vollenweider, R.A., & Kerekes, J., 1982. Eutrophication of waters: monitoring, assessment and control. Paris: Organization for Economic Cooperation and Development.); total phosphorus concentrations range from 20 ± 13 µg/L in lacustrine zone and 31 ± 22 µg/L in the riverine zone. The climate of the region is Cwa (rainy summer and dry winter) according to Köppen (1931)Köppen, W., 1931. Grundriss der klimakunde. Berlin: De Gruyter. http://dx.doi.org/10.1515/9783111667751.
http://dx.doi.org/10.1515/9783111667751...
. It is characterized by the occurrence of a very rainy and hot season, between October and March, when more than 70% of the volume of rainfall is concentrated. The dry season (between April ans September) is cold with little rain (Abreu & Tonello, 2017Abreu, M.C., & Tonello, K.C., 2017. Avaliação dos parâmetros hidrometeorológicos na bacia do rio Sorocaba/SP. Rev. Bras. Meteorol. 32(1), 99-109. http://dx.doi.org/10.1590/0102-778632120150164.
http://dx.doi.org/10.1590/0102-778632120...
). In the period 1991-2020, the average monthly rainfall in the region (State of São Paulo) is 120.6 mm (± 86.2 mm); on average, July and August are the months with the lowest rainfall (monthly accumulated averages: 26.8 and 23.7 mm, respectively) and December and January with the highest (228.3 and 286.6 mm, respectively). The average annual temperature is 20.8 ºC (± 2.0 ºC); on average, the lowest temperatures occur in June and July (17.8 and 17.7 ºC, respectively) and the highest in January and February (22.9 and 23.1 ºC, respectively); INMET (2022)Instituto Nacional de Meteorologia - INMET, 2022. Normais climatológicas do Brasil período 1991-2020. Retrieved in 2022, May 20, from https://portal.inmet.gov.br/normais.
https://portal.inmet.gov.br/normais...
.

The reservoir was built in 1914 by damming the Sorocaba River that is formed by the junction of the Sorocamirim, Sorocabuçu and Una rivers (Figure 1). The main use of the reservoir is electric power generation and water supply for about one million people. The reservoir is 26 km in length, with a drainage area of 936 km2, and a shoreline length of 192.9 km (Table 1). On the right margin of the reservoir, there is an Environmental Protection Area with the Brazilian Savannah and Atlantic Forest species. In the left margin, the land use is predominantly agricultural, covering about 48% of the basin area (Secchin, 2012Secchin, L.F., 2012. Caracterização ambiental e avaliação de cargas difusas da área de drenagem da represa Itupararranga, SP. [Master’s thesis in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/D.18.2012.tde-30082012-150106
https://doi.org/10.11606/D.18.2012.tde-3...
) (Figure 1). The presence of metals (Pb, Cu, Cr, Mn, Fe, Al, Zn) in the sediments of the Itupararanga Reservoir (usually with higher incidences and contents in riverine and transition regions) indicates a recent and weak pollution by agrochemicals and effluents (Rosa et al., 2015Rosa, A.H., Silva, A.A.M.J., Melo, C.A., Moschini-Carlos, V., Guandique, E.G., Fraceto, L.F., & Lourenço, R.W., 2015. Diagnóstico ambiental e avaliação do uso e ocupação do solo visando a sustentabilidade da represa de Itupararanga, importante área da bacia do médio Tietê. In: Pompêo, M., Moschini-Carlos, V., Nishimura, P.Y., Silva, S.C., & Doval, J.C.L., eds. Ecologia de reservatórios e interfaces. São Paulo: USP, 212-231. Retrieved in 2022, May 20, from http://www.livrosabertos.sibi.usp.br/portaldelivrosUSP/catalog/book/35
http://www.livrosabertos.sibi.usp.br/por...
).

Figure 1
Location of São Paulo State in Brazil, Itupararanga Reservoir with the main tributaries and sampling points arrowed (Inlet Point: Headwater; Outlet Point: Dam) and the land use of Itupararanga basin.
Table 1
Main characteristics of the Itupararanga Reservoir (CBA, 1993Companhia Brasileira de Alumínio - CBA, 1993. Relatório técnico administrativo das usinas hidroelétricas. São Paulo: Grupo Votarantim.; ANEEL, 2005Agência Nacional de Energia Elétrica - ANEEL, 2005. Atlas de energia elétrica (2. ed.). Brasília: ANEEL.; UFSCar, 2008Universidade Federal de São Carlos - UFSCar, 2008. Inventário limnológico para o monitoramento da qualidade da água e dos sedimentos do reservatório da UHE Itupararanga (Relatório Final). São Carlos: Fai-UFSCar.).

The adduction flows in the reservoir (generated mainly from the Sorocaba River) normally follow the rainfall regime. The effluent flows essentially result from: i) the maintenance of water for supply; ii) the operating rules of the hydroelectric power plant. Thus, normally, the volume of the reservoir tends to be greater in the spring-summer period. In turn, HRT temporal variations do not have a well-defined seasonality, as they result from the flow regime and water storage management (i.e., electricity generation, and water supply).

2.2. Sampling

Water samples were collected at 2 sampling points in the subsurface of the reservoir (z ≈ 10 cm): in the upstream region (headwater of the reservoir - inlet: 23.6245°S; 47.3364ºW) and in the dam region (limnetic region - outlet: 23.6153ºS; 47.3904ºW) (Figure. 1). The inlet point (located on the Sorocaba River) integrates the chemical and physical characteristics of the main tributaries of the reservoir; the tributaries (Una, Sorocabuçu and Sorocamirim rivers), which are primarily responsible for the tributary flow of the reservoir (over 90%). The choice of only two points (inlet and outlet) stems from the assumptions of calculating the mass balance of the reservoir (Bianchini Junior et al., 2019Bianchini Junior, I., Fushita, Â.T., & Cunha-Santino, M.B., 2019. Evaluating the retention capacity of a new subtropical run-of-river reservoir. Environ. Monit. Assess. 191(3), 161. PMid:30771013. http://dx.doi.org/10.1007/s10661-019-7295-5.
http://dx.doi.org/10.1007/s10661-019-729...
).

For mass balance calculations, the following variables were selected: electrical conductivity (EC; µS/cm), dissolved oxygen (DO; mg/L), water temperature (wT; ºC), pH, turbidity (NTU), total solids (TS; mg/L), chlorophyll-a (Chl-a; mg/L), total phosphorus (TP; µg/L), total Kjeldahl nitrogen (TKN; mg/L), and nitrate (N-NO3; mg/L). Vertical profiles of water temperature, pH, OD, EC and turbidity were determined at the sampling points with a multiparameter probe (YSI Model 556 MPS).

To analyze the TKN, nitrate, TP and TS, the samples were collected and kept refrigerated (4 ºC) during transport to the laboratory. The samples used for chlorophyll-a determinations were filtered at the time of collection, and the filters were kept cool and in the dark until the concentrations were determined.

TKN concentrations were determined by the titrimetric method; nitrate, TP and chlorophyll-a concentrations were measured by colorimetry and TS by gravimetry; these analyses were performed in accordance with the APHA (2005)American Public Health Association - APHA, 2005. Standard methods for examination of water and wastewater (25th ed.). Washington: APHA, AWWA, WEF.. Assessments were conducted quarterly during the period from September 2007 to July 2012 (UFSCar, 2008Universidade Federal de São Carlos - UFSCar, 2008. Inventário limnológico para o monitoramento da qualidade da água e dos sedimentos do reservatório da UHE Itupararanga (Relatório Final). São Carlos: Fai-UFSCar.; Bottino, 2011Bottino, F., 2011. Diversidade, biomassa e decomposição de macrófitas aquáticas no Reservatório Itupararanga - SP [Doctoral dissertation in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/T.18.2011.tde-08022012-104315
https://doi.org/10.11606/T.18.2011.tde-0...
; Casali, 2014Casali, S.P., 2014. A comunidade fitoplanctônica no reservatório de Itupararanga (Bacia do Rio Sorocaba, SP) [Doctoral dissertation in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/T.18.2014.tde-25092014-152955
https://doi.org/10.11606/T.18.2014.tde-2...
); the results of limnological surveys carried out by the environmental agency of the São Paulo State (CETESB) were also utilized (CETESB, 2010Companhia Ambiental do Estado de São Paulo - CETESB, 2010. Relatório de Águas Interiores. Retrieved in 2020, October 20, from https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios.
https://cetesb.sp.gov.br/aguas-interiore...
; 2011Companhia Ambiental do Estado de São Paulo - CETESB, 2011. Relatório de Águas Interiores. Retrieved in 2020, October 20, from https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios.
https://cetesb.sp.gov.br/aguas-interiore...
; 2012Companhia Ambiental do Estado de São Paulo - CETESB, 2012. Relatório de Águas Interiores. Retrieved in 2020, October 20, from https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios.
https://cetesb.sp.gov.br/aguas-interiore...
).

2.3. Mathematical modeling

To describe the mass balance of variables, it was assumed that the Itupararanga Reservoir is a completely mixed system with a step input (Equations 1 and 2; Chapra, 2008Chapra, S.C., 2008. Surface water-quality modeling. Long Grove: Waveland Press.). These assumptions were adopted because they make it possible to perform reservoir mass balances quickly and simply. The selected model (Equation 2) is widely disseminated as it helps decision makers in activities aimed at improving water quality, maintaining autochthonous biodiversity and ecosystem services (e.g., effluent control actions as a function of downstream water quality). According to seasonal vertical thermal profile simulation (Ryan & Harleman, 1971Ryan, P.J., & Harleman, D.R.F., 1971. Prediction of the annual cycle of temperature changes in stratified lake or reservoir: mathematical model and user’s manual. Cambridge, MA: MIT Department of Civil Engineering, Technical Report 137, 132 p.), this reservoir is classified as warm monomictic (sensuHutchinson & Löffler, 1956Hutchinson, G.E., & Löffler, H., 1956. The thermal classification of lakes. Proc. Natl. Acad. Sci. USA 42(2), 84-86. PMid:16589823. http://dx.doi.org/10.1073/pnas.42.2.84.
http://dx.doi.org/10.1073/pnas.42.2.84...
). However, it was possible to admit that it is a thoroughly mixed system because the main entrances and exits of the elements occur mostly primarily through the epilimnion and the metalimnion, since i) tributary water temperatures usually are closer to the epilimnion and hypolimnion temperatures, in this case, the inlet occurs over lake water (i.e., inlet type “overflow”); Ford (1990)Ford, D.E., 1990. Reservoir transport processes. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 15-41.; ii) the reservoir is a shallow environment (average depth: 12 m) and reasonably exposed to wind action (fetch ca. 4.5 km); iii) the water outlet is 15 m above the elevation of the base of the dam; iv) determinations of vertical thermal profiles carried out in the reservoir (this study) indicated that the occurrence of thermal stratification is usually restricted to the lacustrine region, close to the dam; v) usually hypolimnion occurs at high depths (below 8-10 m).

α = Q V + k (1)
C = W α V 1 e α t (2)

where: α = assimilation factor (d-1) (Σ sink); k = first order reaction rate constant (d-1);

QV = hydraulic flushing (HF) (d-1); Q = flow rate (average value on a monthly basis);

V = volume of the reservoir; C = steady state concentration (i.e., related to the dam area - outlet concentrations); W = the loading term (kg d-1) (i.e., the daily load of substance, referring to the reservoir inlets); t = the time required to reach the equilibrium concentration.

The study used three basic assumptions: (i) the reservoir can be represented by a zero-dimensional model (i.e., Continuous Stirred Tank Reactor; CSTR); (ii) the reservoir is in a steady state (for any month, the initial and final values are constant), necessarily generating a “step” loading function (Jørgensen & Fath, 2010Jørgensen, S.E., & Fath, B.D., 2010. Fundamentals of ecological modelling: application in environmental management and research. Amsterdam: Elsevier.), the step input is essentially an “on-off” function that has a jump discontinuity at t = 0 (Chapra, 2008Chapra, S.C., 2008. Surface water-quality modeling. Long Grove: Waveland Press.); (iii) the hydraulic retention time is long enough for the reservoir to reach equilibrium.

The selected model (Equations 1 and 2) assumes that the reservoir is in a stable phase (Cunha-Santino et al., 2017Cunha-Santino, M.B., Fushita, A.T., & Bianchini Junior, I., 2017. A modeling approach for a cascade of reservoirs in the Juquiá-Guaçu River (Atlantic Forest, Brazil). Ecol. Modell. 356, 48-58. http://dx.doi.org/10.1016/j.ecolmodel.2017.04.008.
http://dx.doi.org/10.1016/j.ecolmodel.20...
) and the limnological variable values downstream derive from the initial concentrations, except for the quantities of elements that were trapped (or released) in the reservoir (Teodoru & Wehrli, 2005Teodoru, C., & Wehrli, B., 2005. Retention of sediments and nutrients in the Iron Gate I Reservoir on the Danube River. Biogeochemistry 76(3), 539-565. http://dx.doi.org/10.1007/s10533-005-0230-6.
http://dx.doi.org/10.1007/s10533-005-023...
) during the HRT period. The assimilation factor (α) obtained from the parameterization indicate that there were losses (e.g., sedimentation, biological absorption, chemical reactions and adsorption) or gains of the elements in the reservoir compared to the initial concentrations (in this case, derived from the inputs of the Sorocaba River, upstream of the reservoir). If α is positive and higher than hydraulic flushing (HF), the value indicates that the element is retained in the reservoir; therefore the retention of the element is greater than the flushing. If α is equal to HF, the retention is null, the element is only carried by the flow of the reservoir. The negative value of alpha or α < HF indicate that the retention is null, and the reservoir is the source of element (e.g., internal loading process, lateral runoff, point source of the element) (Cunha-Santino et al., 2017Cunha-Santino, M.B., Fushita, A.T., & Bianchini Junior, I., 2017. A modeling approach for a cascade of reservoirs in the Juquiá-Guaçu River (Atlantic Forest, Brazil). Ecol. Modell. 356, 48-58. http://dx.doi.org/10.1016/j.ecolmodel.2017.04.008.
http://dx.doi.org/10.1016/j.ecolmodel.20...
).

The model parameterization was performed calculating the assimilation factor (i.e., EC, DO, pH, TP, TKN, N-NO3, TS, Turb, and Chl-a) in each month of the study. For this reason, the following were used in Equation 2: (i) the values of each input variable (W on a daily basis); (ii) the monthly average flow of the reservoir; (iii) the monthly average retention time; (iv) the average volume for the month. By substituting the value of alpha by an iterative method (Generalized Reduced Gradient algorithm; Fylstra et al., 1998Fylstra, D., Lasdon, L., Watson, J., & Waren, A., 1998. Design and use of the microsoft excel solver. Interfaces 28(5), 29-55. http://dx.doi.org/10.1287/inte.28.5.29.
http://dx.doi.org/10.1287/inte.28.5.29...
), the values of the calculated variables (C parameter; Equation 2) were the same as the values determined in situ (outlet). To choose the date of the output variable value, the hydraulic retention time criterion was used, i.e., the value of the variable on the day of entry into the head of the reservoir compared with the value for the same variable after the time of HRT, at the dam. For this purpose, the monthly output data of the variables were interpolated.

2.4. Statistical analysis

We drew up a box chart graphic using the software Origin 8 to examine the variability of retention factor among the different variables. The relationship between HRT and retention capability of the reservoir was assessed by performing simple linear regression for TP. Linear regression was performed using the Origin 8 software, which uses the least squares method. The choice of total phosphorus for this analysis was due to its importance for defining the trophic state of the reservoir (Zhou et al. 2016Zhou, B., Xu, Y., Vogt, R.D., Lu, X., Li, X., Deng, X., Yue, A., & Zhu, L., 2016. Effects of land use change on phosphorus levels in surface waters - a case study of a watershed strongly influenced by agriculture. Water Air Soil Pollut. 227(5), 160. http://dx.doi.org/10.1007/s11270-016-2855-6.
http://dx.doi.org/10.1007/s11270-016-285...
).

3. Results

The wet season occurred from January to March and the dry season was from April to June (Figure 2a). The inflow and outflow (Figure 2b) were high mainly between 2009 and 2010 when the accumulated rainfall was higher than 600 mm in the wet season increasing the reservoir volume and decreasing its HRT (Figure 2c).

Figure 2
Daily rainfall in the dam area of Itupararanga Reservoir (a). Hydraulic characteristics of Itupararanga Reservoir: Flow data (b), and volume and Hydraulic Retention Time (HRT) (c); (Bottino, 2011Bottino, F., 2011. Diversidade, biomassa e decomposição de macrófitas aquáticas no Reservatório Itupararanga - SP [Doctoral dissertation in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/T.18.2011.tde-08022012-104315
https://doi.org/10.11606/T.18.2011.tde-0...
; Casali, 2014Casali, S.P., 2014. A comunidade fitoplanctônica no reservatório de Itupararanga (Bacia do Rio Sorocaba, SP) [Doctoral dissertation in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/T.18.2014.tde-25092014-152955
https://doi.org/10.11606/T.18.2014.tde-2...
; Votorantim Energia weather station: 23°36'43.35”S; 47°23'49.87”W).

The maximum EC value (115 µS/cm) was in the dry season (August/08) and a sharp decrease occurred in the summer (37 µS/cm in March/11) (Table 2). The maximum EC retention was 32% (HRT = 182 d). In general, the DO concentrations were low (< 5.0 mg/L) mainly in the wet months during 2007-2011 (Table 2). The lowest DO concentrations (< 2.0 mg/L) were in 2009 (October/09 and December/09). The pH values (Table 2) were typically low in the dry season. However, October and December/09 showed the lowest values (≈ 5.9). The highest value was 8.6 in September/09. The highest pH decrease was 19% (HRT = 186 d). The wT variation corresponded to the seasonal variation (Table 2). The highest wT was 29 ºC (November/09) and the lowest was 16.3 ºC (August/08).

Table 2
Entry and exit times and respective values of electrical conductivity (EC), dissolved oxygen (DO), pH, water temperature (wT) and total phosphorus (TP) from the Itupararanga Reservoir.

The TS concentrations showed a low variation (Table 3). In the dry season of 2008, 2009 and 2010, the concentrations were high (80 to 104 mg/L). Turbidity (Table 3) was high in the wet months, and the maximum value (23 NTU) occurred in November/07. However, in August/08 and July/11 (dry season), the values increased (13 and 9.1 NTU, respectively). The maximum retention of TS and turbidity was 77% (HRT = 216 and 235 d, respectively).

Table 3
Entry and exit times and respective values of total Kjeldhal nitrogen (TKN), nitrate (nitrogen basis; N-NO3), total solids (TS), turbidity (Turb) and chlorophyll-a (Chl-a) from the Itupararanga Reservoir.

The highest nutrient concentrations occurred in periods of high precipitation ranging from 74 to 80 µg/L for TP and up to 1.0 mg/L for TKN and maximum of 0.5 mg/L for N-NO3 (Table 3). However, the wet season contributed to the increase of nutrient concentrations (from November to March). In September/09, the TKN concentration increased (1.13 mg/L). The Chl-a concentrations increased over time with the highest value (37 µg/L) in November/11 and the lowest values in November/09 (Table 3). The maximum retention was 74% for TP and 58.5% for TKN with HRT = 163 d and 126 d, respectively. For Chl-a, the maximum retention (49%) occurred when the HRT was lowest (99 d).

The alpha values showed the same pattern of temporal variation for EC, DO, pH, Turbidity and Chl-a (Figures 3 and 4). Regardless of seasonality, the highest values of parameterized alpha occurred in wet months (from 90 mm to 504 mm) for all variables. Nutrients (TP, NKT and N-NO3) showed the highest alpha values (3.82 d-1, 6.8 d-1, 3.67 d-1). In December/09 (monthly precipitation: 260 mm), the HF values were high for EC and pH (1.41% and 1.64% of alpha). When the precipitation exceeded 500 mm (January/10), the HF showed the highest values for DO, TS and Turbidity corresponding to 4.06%, 1.02% and 1.14% of alpha values, respectively. However, the highest HF percentage was 7.5% for Chl-a in November/09 when the highest monthly precipitation occurred (319 mm).

Figure 3
Water temperature changes in inlet and outlet sites of reservoir. Electric Conductivity, dissolved oxygen, pH, and Total Phosphorous (TP) retention parameters (α; black line and square), and hydraulic flows (HF; red line and circle) in terms of α percentage.
Figure 4
Total Kjeldhal nitrogen (TKN), nitrate (in N basis; N-NO3), total solids (TS), turbidity, and chlorophyll-a retention parameters (α; black line and square), and hydraulic flows (HF; red line and circle) in terms of α percentage.

The estimation of retention potential resulted in 273 values of alpha parameter for 9 water variables and there were no negative values (Figure 5a). Nutrients, mainly TP and TKN showed the highest variation in the study period. Turbidity and chlorophyll a concentration also showed high data dispersion, while EC, DO, and pH had low variability. Taking into account that α=HF+k (Equation 1), the comparison of the composition of the alpha values indicated that the reaction rate constants (k) were much higher than the hydraulic flushing values, pointing out that for the conditions of the Itupararanga Reservoir (small reservoir with high HRT), the sink processes prevailed over the exportations. According to the results (Figures 4 and 5), the percent contributions of HF to the alpha values were low (less than 2%) for CE, pH, TS and turbidity. For DO and nutrients, the HF contribution was higher than 2% and for chlorophyll a, it was 7.53%. Overall, the TP retention capability of this reservoir did not show a substantial rise with HRT increase (Figure 5b). However, increasing flow makes HF assume more relevance (Figure 5c) in the retention capability due to the decrease of the reaction rate constant (k) to offset the assimilation factor (α).

Figure 5
Alfa values for different water variables (a). Relationships between hydraulic retention time (HRT)of the Itupararanga Reservoir, and the retention capability (α) for total phosphorus (b), hydraulic flows (HF) (c).

4. Discussion

Damming a river change both water flow and the transport of materials, by altering settling or mineralization/immobilization patterns in the system (Kumwimba et al., 2022Kumwimba, M.N., Bao, L., Jie, Z., Li, X., Huang, J., Wang, W., Li, X., Su, J., Muyembe, D.K., Guide, A., & Dzakpasu, M., 2022. Nutrients retention of a series of small dam-impacted urban rivers in northern China. J. Environ. Chem. Eng. 10(3), 107967. http://dx.doi.org/10.1016/j.jece.2022.107967.
http://dx.doi.org/10.1016/j.jece.2022.10...
). The potential of retention and recycling by lentic ecosystems manifests itself across different spatial scales (Cheng & Basu, 2017Cheng, F.Y., & Basu, N.B., 2017. Biogeochemical hotspots: role of small water bodies in landscape nutrient processing. Water Resour. Res. 53(6), 5038-5056. http://dx.doi.org/10.1002/2016WR020102.
http://dx.doi.org/10.1002/2016WR020102...
). Our results demonstrated that substances sink from inlet to outlet, despite seasonality, the high HRT of the system (in general > 120 days during the study period), a large drainage area, and a low variation of water level. The maintenance of similarity in the inflow and outflow even during low rainfall periods supports high HRT, which is the key factor controlling the mass balance in man-made reservoirs.

Retention time is the main factor related to the trapping of elements in reservoirs as it underpins the settling and biological processes that drive the nutrient loss, toxicity of sediments and the increase in phytoplankton (Kawara et al., 1998Kawara, O., Yura, E., Fujii, S., & Matsumoto, T., 1998. A study on the role of hydraulic retention time in eutrophication of the Asahi River Dam reservoir. Water Sci. Technol. 37(2), 245-252. http://dx.doi.org/10.2166/wst.1998.0146.
http://dx.doi.org/10.2166/wst.1998.0146...
; Kõiv et al., 2011Kõiv, T., Nõges, T., & Laas, A., 2011. Phosphorus retention as a function of external loading, hydraulic turnover time, area and relative depth in 54 lakes and reservoirs. Hydrobiologia 660(1), 105-115. http://dx.doi.org/10.1007/s10750-010-0411-8.
http://dx.doi.org/10.1007/s10750-010-041...
; Hansen et al., 2016Hansen, E., Chan, K.-S., Jones, C.S., & Schilling, K., 2016. Assessing the relative importance of nitrogen-retention processes in a large reservoir using time-series modeling. J. Agric. Biol. Environ. Stat. 21(1), 152-169. http://dx.doi.org/10.1007/s13253-015-0218-1.
http://dx.doi.org/10.1007/s13253-015-021...
). Many studies report the role of the reservoirs as retention basins that contribute to decreases of up to 70% of sediment flux (Dai et al., 2009Dai, A., Qian, T., Trenberth, K.E., & Milliman, J.D., 2009. Changes in continental freshwater discharge from 1948 to 2004. J. Clim. 22(10), 2773-2792. http://dx.doi.org/10.1175/2008JCLI2592.1.
http://dx.doi.org/10.1175/2008JCLI2592.1...
), nearly 90% of suspended solids (Némery et al., 2016Némery, J., Gratiot, N., Doan, P.T.K., Duvert, C., Alvarado-Villanueva, R., & Duwig, C., 2016. Carbon, nitrogen, phosphorus, and sediment sources and retention in a small eutrophic tropical reservoir. Aquat. Sci. 78(1), 171-189. http://dx.doi.org/10.1007/s00027-015-0416-5.
http://dx.doi.org/10.1007/s00027-015-041...
), approximately 12% of the global phosphorus loading (Maavara et al., 2015Maavara, T., Parsons, C.T., Ridenour, C., Stojanovic, S., Dürr, H.H., Powley, H.R., & Van Cappellen, P., 2015. Global phosphorus retention by river damming. Proc. Natl. Acad. Sci. USA 112(51), 15603-15608. PMid:26644553. http://dx.doi.org/10.1073/pnas.1511797112.
http://dx.doi.org/10.1073/pnas.151179711...
) and about 47% of nitrate (Cheng & Basu, 2017Cheng, F.Y., & Basu, N.B., 2017. Biogeochemical hotspots: role of small water bodies in landscape nutrient processing. Water Resour. Res. 53(6), 5038-5056. http://dx.doi.org/10.1002/2016WR020102.
http://dx.doi.org/10.1002/2016WR020102...
). Thus, HRT regulates the boundary of the role of a reservoir as a sink or source of compounds controlling the main ecological processes occurring within the system (Pacheco et al., 2015Pacheco, F.S., Soares, M.C.S., Assireu, A.T., Curtarelli, M.P., Roland, F., Abril, G., Stech, J.L., Alvalá, P.C., & Ometto, J.P., 2015. The effects of river inflow and retention time on the spatial heterogeneity of chlorophyll and water-air CO2 fluxes in a tropical hydropower reservoir. Biogeosciences 12(1), 147-162. http://dx.doi.org/10.5194/bg-12-147-2015.
http://dx.doi.org/10.5194/bg-12-147-2015...
).

Modeling of the biogeochemistry characteristics of reservoirs, including the mass balances, is a valuable tool for water management, especially in the tropics where eutrophication is a major concern (Bartoszek & Koszelnik, 2016Bartoszek, L., & Koszelnik, P., 2016. The qualitative and quantitative analysis of the coupled C, N, P and Si retention in complex of water reservoirs. Springerplus 5(1), 1157. PMid:27504255. http://dx.doi.org/10.1186/s40064-016-2836-7.
http://dx.doi.org/10.1186/s40064-016-283...
; Némery et al., 2016Némery, J., Gratiot, N., Doan, P.T.K., Duvert, C., Alvarado-Villanueva, R., & Duwig, C., 2016. Carbon, nitrogen, phosphorus, and sediment sources and retention in a small eutrophic tropical reservoir. Aquat. Sci. 78(1), 171-189. http://dx.doi.org/10.1007/s00027-015-0416-5.
http://dx.doi.org/10.1007/s00027-015-041...
). By using 9 water variables, we observed in general decreased outlet concentrations in relation to the inlet conditions despite seasonality effects, except for DO. This reduction is an expected pattern for long-established reservoirs as the low variation of hydraulic characteristics (i.e. water level, drainage flow) allows the internal processing of the elements (Gonzaga et al., 2007Gonzaga, M., Cunha-Santino, M.B., & Bianchini Junior, I., 2007. Methodological test of efficiency of heterotrophic potential. Acta Sci. Biol. Sci. 29(2), 197-201. https://doi.org/10.4025/actascibiolsci.v29i2.526.
https://doi.org/10.4025/actascibiolsci.v...
; Araújo et al., 2011Araújo, F.G., Azevedo, M.C.C., & Ferreira, M.N.L., 2011. Seasonal changes and spatial variation in the water quality of a eutrophic tropical reservoir determined by the inflowing river. Lake Reserv. Manage. 27(4), 343-354. http://dx.doi.org/10.1080/07438141.2011.627753.
http://dx.doi.org/10.1080/07438141.2011....
).

The high values of assimilation factor (α) indicated a high capacity of retention which is related to the high values of the reaction coefficient (k) that, in turn, depends on the HRT. High HRT decreases the HF supporting the reactions of mineralization or immobilization over time, hence contributing to the internal processing of elements, mainly in the hypolimnion, or in the sediments that may be a reactive fraction in many tropical reservoirs. In small reservoirs (low area:volume ratio), the reaction coefficient is inversely proportional to the hydraulic residence time. Inverse relationships were also observed between k and the surface area and depth of the water bodies, alluding to the strong size control on the reactivity of a water body (Cheng & Basu, 2017Cheng, F.Y., & Basu, N.B., 2017. Biogeochemical hotspots: role of small water bodies in landscape nutrient processing. Water Resour. Res. 53(6), 5038-5056. http://dx.doi.org/10.1002/2016WR020102.
http://dx.doi.org/10.1002/2016WR020102...
). As reservoir management flow tends to follow the pattern of precipitation, the hydrological conditions may be important to predict the retention. Therefore, changes in the rainfall regime influence the flushing of elements downstream. However, for the Itupararanga Reservoir, this relationship was not clear considering the increase in the precipitation in the dry season (from 235 to 661 mm) that generated HRT similar to the wet period but showed minor importance to the exportation of pollutants.

The model parameterization indicated a weak positive relationship of HRT and the assimilation factor, pointing out that, under current hydraulic conditions, the trapping of elements is prevalent. However, the HRT reduction (owing to the higher intensities and duration of rainfall) can contribute to the rise of internal loading (by changing the thermal stratification regime), which may increase the eutrophication due to the increase in water input derived from a cultured drainage area. It is widely known that dams affect the sediment transport and water transparency, and our results showed a similar pattern for retention of solids and turbidity (> 70%). Total phosphorus also showed retention up to 70%. Phosphorus is the main nutrient causing concerns for freshwater management and our results suggest that the outflow is a good predictor of retention capability of TP. For relatively constant volume, decreasing or increasing the outflow favors and disfavors, respectively, the assimilation factor. However, for intermediate values of outflow (≈ 9 - 20 m3 s-1), the retention capability is high because HRT remains high.

Phosphorus retention in reservoirs increases the reactivity of sediments in these ecosystems contributing to P internal loading in the long term, responsible for algal and cyanobacterial blooms and eutrophication (Harper, 1982Harper, D., 1982. Eutrophication of freshwaters. London: Chapman & Hall.; Cooke et al., 1993Cooke, G.D., Welch, E.B., Peterson, S.A., & Newroth, P.R., 1993. Restoration and management of lakes and reservoirs. Boca Raton: Lewis Publishers.; Horne & Goldman, 1994Horne, A.J., & Goldman, C.R., 1994. Limnology. New York: McGraw-Hill.; Vo et al., 2014Vo, N.X.Q., Doan, T.V., & Kang, H., 2014. Impoundments increase potential for phosphorus retention and remobilization in an urban stream. Environ. Eng. Res. 19(2), 175-184. http://dx.doi.org/10.4491/eer.2014.19.2.175.
http://dx.doi.org/10.4491/eer.2014.19.2....
; Qin et al., 2020Qin, L., Lei, P., Lei, Q., Liu, H., Li, X., Zhang, H., & Lindsey, S., 2020. Evaluating the effect of dam construction on the phosphorus fractions in sediments in a reservoir of drinking water source, China. Environ. Monit. Assess. 192(2), 99. PMid:31912244. http://dx.doi.org/10.1007/s10661-019-8053-4.
http://dx.doi.org/10.1007/s10661-019-805...
). The model parameterization indicated that across the spatial scale from the headwater to the dam, the phosphorus trapping on the sediment increases due to the rising of oxidant conditions and neutral pH. These results have important implications for water quality and eutrophication control. However, when we consider this reservoir, there is an increase in cyanobacterial blooms (https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios/) even with low phosphorus concentration in the water column. Therefore, phosphorus trapping alone did not prevent phytoplankton growth, which is possibly influenced by transport of phytoplankton from tributary rivers, and nitrogen limitation (Beghelli et al., 2016Beghelli, F.G.S., Frascareli, D., Pompêo, M.L.M., & Moschini-Carlos, V., 2016. Trophic state evolution over 15 years in a tropical reservoir with low nitrogen concentrations and cyanobacteria predominance. Water Air Soil Pollut. 227(3), 95. http://dx.doi.org/10.1007/s11270-016-2795-1.
http://dx.doi.org/10.1007/s11270-016-279...
).

Thermal stratification plays an important role in retention capability of reservoirs as part of the elements are trapped in the hypolimnion while the epilimnion remains element-poor due to uptake and metabolism by organisms, especially algae and cyanobacteria, in that part of the water column (Némery et al., 2016Némery, J., Gratiot, N., Doan, P.T.K., Duvert, C., Alvarado-Villanueva, R., & Duwig, C., 2016. Carbon, nitrogen, phosphorus, and sediment sources and retention in a small eutrophic tropical reservoir. Aquat. Sci. 78(1), 171-189. http://dx.doi.org/10.1007/s00027-015-0416-5.
http://dx.doi.org/10.1007/s00027-015-041...
). This is especially important in warm monomictic systems (Hutchinson & Löffler, 1956Hutchinson, G.E., & Löffler, H., 1956. The thermal classification of lakes. Proc. Natl. Acad. Sci. USA 42(2), 84-86. PMid:16589823. http://dx.doi.org/10.1073/pnas.42.2.84.
http://dx.doi.org/10.1073/pnas.42.2.84...
), with progressive stratification over the year such as what frequently occurs in the tropics. These characteristics (i.e., effective retention of phosphorus, TS, turbidity and chlorophyll a) are relevant in a reservoir with catchment flow ranging from 2.1 to 6.5 m3/s, qualifying this environment as a reliable source of water for supply, and reducing costs of treatment, mainly if the water catchment is derived from the epilimnion. In addition, the water quality improvement promotes the multiple uses of the reservoirs (irrigation, fishing, leisure). However, the monitoring of water inflow and outflow are important to an efficient management plan, considering that changes in the HRT could influence in the mass balance in the ecosystem.

Monitoring water flux is essential to support strategies for water management mainly in the tropics (Fowe et al., 2015Fowe, T., Karambiri, H., Paturel, J.-E., Poussin, J.-C., & Cecchi, P., 2015. Water balance of small reservoirs in the Volta basin: a case study of Boura reservoir in Burkina Faso. Agric. Water Manage. 152, 99-109. http://dx.doi.org/10.1016/j.agwat.2015.01.006.
http://dx.doi.org/10.1016/j.agwat.2015.0...
). In this context, it is important to highlight that, by exclusively using input and output data, the calculations of the assimilation coefficients include characteristics of the elements circulation, regardless of the occurrence of thermal stratification of the aquatic environment. Thus, for example, in the case of phosphorus, occurrences of internal loading result in a decrease in the alpha value (i.e., lower retention capacity) or, in the case of predominance of this event, in negative values (i.e., the aquatic system is a source of phosphorus, resulting in higher outlet concentrations than inlet concentrations). On the other hand, if phosphorus is incorporated into the biomass (e.g., primary production), trapped in the hypolimnion or in the sedimentary surface of the environment with full circulation of the water column, the alpha value increases. A relevant difficulty factor in the application of this model is to define the period required for the longitudinal transport of the elements; in this case, the presence of thermal stratification may interfere with the calculation of the HRT. Specifically for the Itupararanga Reservoir, although the presence of the hypolimnion is frequent, its volume is very small in relation to the volumes of the epilimnion and metalimnion (Ryan & Harleman, 1971Ryan, P.J., & Harleman, D.R.F., 1971. Prediction of the annual cycle of temperature changes in stratified lake or reservoir: mathematical model and user’s manual. Cambridge, MA: MIT Department of Civil Engineering, Technical Report 137, 132 p.), not interfering substantially in the HRT calculation.

Changes in rainfall pattern would cause changes in HRT. In a land use change scenario, the input of pollutants will also increase due to the high runoff leading to the HRT decreasing, and consequently influencing the retention capability. The increase in rainfall events with high magnitude are a consequence of the human activities and may dramatically alter water quality. Changes in land use are one of the key factors affecting ecosystem services and human welfare due to its influence on the water quality (Costanza et al., 1997Costanza, R., d’Arge, R., de Groot, R., Farber, S., Grasso, M., Hannon, B., Limburg, K., Naeem, S., O’Neill, R.V., Paruelo, J., Raskin, R.G., Sutton, P., & van den Belt, M., 1997. The value of the world’s ecosystem services and natural capital. Nature 387(6630), 253-260. http://dx.doi.org/10.1038/387253a0.
http://dx.doi.org/10.1038/387253a0...
). In this context, the efficiency of retention capability of an aquatic ecosystem is closely connected with the current distribution of land use (Prokopová et al., 2019Prokopová, M., Salvati, L., Egidi, G., & Cudlín, O., 2019. Envisioning present and future land-use change under varying ecological regimes and their influence on landscape stability. Sustainability 11(17), 654. http://dx.doi.org/10.3390/su11174654.
http://dx.doi.org/10.3390/su11174654...
).

Reservoirs with low area:volume ratio and high HRT act as sink systems, producing a gradient of water quality improvement from headwater to the dam (Thornton, 1990Thornton, K.W., 1990. Sedimentary processes. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 43-69.; Straškraba, 1999Straškraba, M., 1999. Retention time as a key variable of reservoir limnology. In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 385-410.). Thus, due to the strong prevalence of chemical and biochemical reactions (e.g., phytoplankton and macrophyte growths) over the hydraulic flow in the Itupararanga reservoir, in the scenario of increased temperature and maintenance of current flows, the retention of elements may be even more accentuated because the reaction rates are sensitive to the temperature variation (Jørgensen & Fath, 2010Jørgensen, S.E., & Fath, B.D., 2010. Fundamentals of ecological modelling: application in environmental management and research. Amsterdam: Elsevier.). In this situation, the longitudinal gradient of the concentrations should be more evident than those usually observed (UFSCar, 2008Universidade Federal de São Carlos - UFSCar, 2008. Inventário limnológico para o monitoramento da qualidade da água e dos sedimentos do reservatório da UHE Itupararanga (Relatório Final). São Carlos: Fai-UFSCar.). However, in a scenario of increasing rainfall, the HRT can become more similar in the wet and dry seasons and the retentive capacity of an aquatic ecosystem may present some alterations (owing to the change in the thermal regime and the increase in internal loading), and the exportation of reactive elements would impact downstream environments, and their ecosystem services. Considering the changes in the water temperature and the shorter time of persistence of thermal stratification (i.e., increase in the period of complete circulation of the water column) owing to the increase in the flow rates, the toxic conditions prevalent in the hypolimnion may be exposed for a longer time. This could cause nutrient enrichment in the epilimnion, resulting in an increase in eutrophication, and leading to a decrease in the multiple uses of the reservoirs and changes in energy flows and trophic chains in the aquatic system. Finally, water intake from the bottom of the reservoir would damage the rivers downstream.

5. Conclusions

In the Itupararanga Reservoir, hydraulic characteristics determined a high capacity for the retention of elements, as demonstrated by the high reaction assimilation coefficient (α), due to high-rate constants (k), and low hydraulic flushing. Except for DO concentrations, all the selected values decreased from the inlet to the outlet. Seasonality was of minor importance to trapping elements, but retention capability is related to HRT, which is in turn linked to the hydrological, morphological, and biological characteristics of the ecosystem. However, the similarity of HRT during the wet and dry season indicated that hydraulic characteristics, such as low area:volume ratio were relevant for retention capability. In the tropics, small reservoirs with high HRT (> 120 days) can be used to retain elements considering the provisioning service (e.g., water supply). However, the sediment conditions must be recognized as an important source of pollutants if thermal stratification persists or if the oxidant conditions change, for example. Understanding the hydraulic characteristics of man-made reservoirs can help to provide informed management tools for these ecosystems that could improve the management of issues, such as phosphorus loading and eutrophication.

Acknowledgements

The authors are grateful to the Companhia Brasileira de Alumínio (CBA-Votorantim), currently Votorantim Energia, for providing data about the reservoir, for supporting part of the field sampling, and for providing part of the limnological data. We would also like to thank the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP process number: 08/55636-9), and Dr. Kevin Murphy (University of Glasgow) for his critical proof reading of the manuscript.

  • Cite as: Bottino, F. et al. Potential of the retention capability of a Neotropical reservoir (São Paulo State, Brazil). Acta Limnologica Brasiliensia, 2023, vol. 35, e4.

References

  • Abreu, M.C., & Tonello, K.C., 2017. Avaliação dos parâmetros hidrometeorológicos na bacia do rio Sorocaba/SP. Rev. Bras. Meteorol. 32(1), 99-109. http://dx.doi.org/10.1590/0102-778632120150164
    » http://dx.doi.org/10.1590/0102-778632120150164
  • Agência Nacional de Energia Elétrica - ANEEL, 2005. Atlas de energia elétrica (2. ed.). Brasília: ANEEL.
  • Akbarzadeh, Z., Maavara, T., Slowinski, S., & Van Cappellen, P., 2019. Effects of damming on river nitrogen fluxes: a global analysis. Global Biogeochem. Cycles 33(11), 1339-1357. http://dx.doi.org/10.1029/2019GB006222
    » http://dx.doi.org/10.1029/2019GB006222
  • American Public Health Association - APHA, 2005. Standard methods for examination of water and wastewater (25th ed.). Washington: APHA, AWWA, WEF.
  • Araújo, F.G., Azevedo, M.C.C., & Ferreira, M.N.L., 2011. Seasonal changes and spatial variation in the water quality of a eutrophic tropical reservoir determined by the inflowing river. Lake Reserv. Manage. 27(4), 343-354. http://dx.doi.org/10.1080/07438141.2011.627753
    » http://dx.doi.org/10.1080/07438141.2011.627753
  • Armengol, J., Garcia, J.C., Comerma, M., Romero, M., Dolz, J., Roura, M., Han, B.H., Vidal, A., & Šimek, K., 1999. Longitudinal processes in canyon type reservoirs: the case of Sau (N.E. Spain). In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 313-345.
  • Bartoszek, L., & Koszelnik, P., 2016. The qualitative and quantitative analysis of the coupled C, N, P and Si retention in complex of water reservoirs. Springerplus 5(1), 1157. PMid:27504255. http://dx.doi.org/10.1186/s40064-016-2836-7
    » http://dx.doi.org/10.1186/s40064-016-2836-7
  • Beghelli, F.G.S., Frascareli, D., Pompêo, M.L.M., & Moschini-Carlos, V., 2016. Trophic state evolution over 15 years in a tropical reservoir with low nitrogen concentrations and cyanobacteria predominance. Water Air Soil Pollut. 227(3), 95. http://dx.doi.org/10.1007/s11270-016-2795-1
    » http://dx.doi.org/10.1007/s11270-016-2795-1
  • Bianchini Junior, I., Fushita, Â.T., & Cunha-Santino, M.B., 2019. Evaluating the retention capacity of a new subtropical run-of-river reservoir. Environ. Monit. Assess. 191(3), 161. PMid:30771013. http://dx.doi.org/10.1007/s10661-019-7295-5
    » http://dx.doi.org/10.1007/s10661-019-7295-5
  • Bottino, F., 2011. Diversidade, biomassa e decomposição de macrófitas aquáticas no Reservatório Itupararanga - SP [Doctoral dissertation in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/T.18.2011.tde-08022012-104315
    » https://doi.org/10.11606/T.18.2011.tde-08022012-104315
  • Casali, S.P., 2014. A comunidade fitoplanctônica no reservatório de Itupararanga (Bacia do Rio Sorocaba, SP) [Doctoral dissertation in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/T.18.2014.tde-25092014-152955
    » https://doi.org/10.11606/T.18.2014.tde-25092014-152955
  • Chapra, S.C., 2008. Surface water-quality modeling. Long Grove: Waveland Press.
  • Cheng, F.Y., & Basu, N.B., 2017. Biogeochemical hotspots: role of small water bodies in landscape nutrient processing. Water Resour. Res. 53(6), 5038-5056. http://dx.doi.org/10.1002/2016WR020102
    » http://dx.doi.org/10.1002/2016WR020102
  • Cole, T.M., & Hannan, H.H., 1990. Dissolved oxygen dynamics. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 71-107.
  • Companhia Ambiental do Estado de São Paulo - CETESB, 2010. Relatório de Águas Interiores. Retrieved in 2020, October 20, from https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios
    » https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios
  • Companhia Ambiental do Estado de São Paulo - CETESB, 2011. Relatório de Águas Interiores. Retrieved in 2020, October 20, from https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios
    » https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios
  • Companhia Ambiental do Estado de São Paulo - CETESB, 2012. Relatório de Águas Interiores. Retrieved in 2020, October 20, from https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios
    » https://cetesb.sp.gov.br/aguas-interiores/publicacoes-e-relatorios
  • Companhia Brasileira de Alumínio - CBA, 1993. Relatório técnico administrativo das usinas hidroelétricas. São Paulo: Grupo Votarantim.
  • Cooke, G.D., Welch, E.B., Peterson, S.A., & Newroth, P.R., 1993. Restoration and management of lakes and reservoirs. Boca Raton: Lewis Publishers.
  • Costanza, R., d’Arge, R., de Groot, R., Farber, S., Grasso, M., Hannon, B., Limburg, K., Naeem, S., O’Neill, R.V., Paruelo, J., Raskin, R.G., Sutton, P., & van den Belt, M., 1997. The value of the world’s ecosystem services and natural capital. Nature 387(6630), 253-260. http://dx.doi.org/10.1038/387253a0
    » http://dx.doi.org/10.1038/387253a0
  • Cunha-Santino, M.B., Fushita, A.T., & Bianchini Junior, I., 2017. A modeling approach for a cascade of reservoirs in the Juquiá-Guaçu River (Atlantic Forest, Brazil). Ecol. Modell. 356, 48-58. http://dx.doi.org/10.1016/j.ecolmodel.2017.04.008
    » http://dx.doi.org/10.1016/j.ecolmodel.2017.04.008
  • Dai, A., Qian, T., Trenberth, K.E., & Milliman, J.D., 2009. Changes in continental freshwater discharge from 1948 to 2004. J. Clim. 22(10), 2773-2792. http://dx.doi.org/10.1175/2008JCLI2592.1
    » http://dx.doi.org/10.1175/2008JCLI2592.1
  • Davidson, E.A., & Janssens, I.A., 2006. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440(7081), 165-173. PMid:16525463. http://dx.doi.org/10.1038/nature04514
    » http://dx.doi.org/10.1038/nature04514
  • Ford, D.E., 1990. Reservoir transport processes. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 15-41.
  • Fowe, T., Karambiri, H., Paturel, J.-E., Poussin, J.-C., & Cecchi, P., 2015. Water balance of small reservoirs in the Volta basin: a case study of Boura reservoir in Burkina Faso. Agric. Water Manage. 152, 99-109. http://dx.doi.org/10.1016/j.agwat.2015.01.006
    » http://dx.doi.org/10.1016/j.agwat.2015.01.006
  • Fylstra, D., Lasdon, L., Watson, J., & Waren, A., 1998. Design and use of the microsoft excel solver. Interfaces 28(5), 29-55. http://dx.doi.org/10.1287/inte.28.5.29
    » http://dx.doi.org/10.1287/inte.28.5.29
  • Gonzaga, M., Cunha-Santino, M.B., & Bianchini Junior, I., 2007. Methodological test of efficiency of heterotrophic potential. Acta Sci. Biol. Sci. 29(2), 197-201. https://doi.org/10.4025/actascibiolsci.v29i2.526
    » https://doi.org/10.4025/actascibiolsci.v29i2.526
  • Hansen, E., Chan, K.-S., Jones, C.S., & Schilling, K., 2016. Assessing the relative importance of nitrogen-retention processes in a large reservoir using time-series modeling. J. Agric. Biol. Environ. Stat. 21(1), 152-169. http://dx.doi.org/10.1007/s13253-015-0218-1
    » http://dx.doi.org/10.1007/s13253-015-0218-1
  • Harper, D., 1982. Eutrophication of freshwaters. London: Chapman & Hall.
  • Hoeinghaus, D.J., Agostinho, A.A., Gomes, L.C., Pelicice, F.M., Okada, E.K., Latini, J.D., Kashiwaqui, E.A.L., & Winemiller, K.O., 2009. Effects of river impoundment on ecosystem services of large tropical rivers: embodied energy and market value of artisanal fisheries. Conserv. Biol. 23(5), 1222-1231. PMid:19459891. http://dx.doi.org/10.1111/j.1523-1739.2009.01248.x
    » http://dx.doi.org/10.1111/j.1523-1739.2009.01248.x
  • Horne, A.J., & Goldman, C.R., 1994. Limnology. New York: McGraw-Hill.
  • Hutchinson, G.E., & Löffler, H., 1956. The thermal classification of lakes. Proc. Natl. Acad. Sci. USA 42(2), 84-86. PMid:16589823. http://dx.doi.org/10.1073/pnas.42.2.84
    » http://dx.doi.org/10.1073/pnas.42.2.84
  • Instituto Nacional de Meteorologia - INMET, 2022. Normais climatológicas do Brasil período 1991-2020. Retrieved in 2022, May 20, from https://portal.inmet.gov.br/normais
    » https://portal.inmet.gov.br/normais
  • Jørgensen, S.E., & Fath, B.D., 2010. Fundamentals of ecological modelling: application in environmental management and research. Amsterdam: Elsevier.
  • Kawara, O., Yura, E., Fujii, S., & Matsumoto, T., 1998. A study on the role of hydraulic retention time in eutrophication of the Asahi River Dam reservoir. Water Sci. Technol. 37(2), 245-252. http://dx.doi.org/10.2166/wst.1998.0146
    » http://dx.doi.org/10.2166/wst.1998.0146
  • Kennedy, R.H., & Walker, W.W., 1990. Reservoir nutrient dynamics. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 109-131.
  • Kerimoglu, O., & Rinke, K., 2013. Stratification dynamics in a shallow reservoir under different hydro-meteorological scenarios and operational strategies. Water Resour. Res. 49(11), 7518-7527. http://dx.doi.org/10.1002/2013WR013520
    » http://dx.doi.org/10.1002/2013WR013520
  • Kõiv, T., Nõges, T., & Laas, A., 2011. Phosphorus retention as a function of external loading, hydraulic turnover time, area and relative depth in 54 lakes and reservoirs. Hydrobiologia 660(1), 105-115. http://dx.doi.org/10.1007/s10750-010-0411-8
    » http://dx.doi.org/10.1007/s10750-010-0411-8
  • Köppen, W., 1931. Grundriss der klimakunde. Berlin: De Gruyter. http://dx.doi.org/10.1515/9783111667751
    » http://dx.doi.org/10.1515/9783111667751
  • Kumwimba, M.N., Bao, L., Jie, Z., Li, X., Huang, J., Wang, W., Li, X., Su, J., Muyembe, D.K., Guide, A., & Dzakpasu, M., 2022. Nutrients retention of a series of small dam-impacted urban rivers in northern China. J. Environ. Chem. Eng. 10(3), 107967. http://dx.doi.org/10.1016/j.jece.2022.107967
    » http://dx.doi.org/10.1016/j.jece.2022.107967
  • Maavara, T., Parsons, C.T., Ridenour, C., Stojanovic, S., Dürr, H.H., Powley, H.R., & Van Cappellen, P., 2015. Global phosphorus retention by river damming. Proc. Natl. Acad. Sci. USA 112(51), 15603-15608. PMid:26644553. http://dx.doi.org/10.1073/pnas.1511797112
    » http://dx.doi.org/10.1073/pnas.1511797112
  • Némery, J., Gratiot, N., Doan, P.T.K., Duvert, C., Alvarado-Villanueva, R., & Duwig, C., 2016. Carbon, nitrogen, phosphorus, and sediment sources and retention in a small eutrophic tropical reservoir. Aquat. Sci. 78(1), 171-189. http://dx.doi.org/10.1007/s00027-015-0416-5
    » http://dx.doi.org/10.1007/s00027-015-0416-5
  • Pacheco, F.S., Soares, M.C.S., Assireu, A.T., Curtarelli, M.P., Roland, F., Abril, G., Stech, J.L., Alvalá, P.C., & Ometto, J.P., 2015. The effects of river inflow and retention time on the spatial heterogeneity of chlorophyll and water-air CO2 fluxes in a tropical hydropower reservoir. Biogeosciences 12(1), 147-162. http://dx.doi.org/10.5194/bg-12-147-2015
    » http://dx.doi.org/10.5194/bg-12-147-2015
  • Prokopová, M., Salvati, L., Egidi, G., & Cudlín, O., 2019. Envisioning present and future land-use change under varying ecological regimes and their influence on landscape stability. Sustainability 11(17), 654. http://dx.doi.org/10.3390/su11174654
    » http://dx.doi.org/10.3390/su11174654
  • Qin, L., Lei, P., Lei, Q., Liu, H., Li, X., Zhang, H., & Lindsey, S., 2020. Evaluating the effect of dam construction on the phosphorus fractions in sediments in a reservoir of drinking water source, China. Environ. Monit. Assess. 192(2), 99. PMid:31912244. http://dx.doi.org/10.1007/s10661-019-8053-4
    » http://dx.doi.org/10.1007/s10661-019-8053-4
  • Rosa, A.H., Silva, A.A.M.J., Melo, C.A., Moschini-Carlos, V., Guandique, E.G., Fraceto, L.F., & Lourenço, R.W., 2015. Diagnóstico ambiental e avaliação do uso e ocupação do solo visando a sustentabilidade da represa de Itupararanga, importante área da bacia do médio Tietê. In: Pompêo, M., Moschini-Carlos, V., Nishimura, P.Y., Silva, S.C., & Doval, J.C.L., eds. Ecologia de reservatórios e interfaces. São Paulo: USP, 212-231. Retrieved in 2022, May 20, from http://www.livrosabertos.sibi.usp.br/portaldelivrosUSP/catalog/book/35
    » http://www.livrosabertos.sibi.usp.br/portaldelivrosUSP/catalog/book/35
  • Ryan, P.J., & Harleman, D.R.F., 1971. Prediction of the annual cycle of temperature changes in stratified lake or reservoir: mathematical model and user’s manual. Cambridge, MA: MIT Department of Civil Engineering, Technical Report 137, 132 p.
  • Secchin, L.F., 2012. Caracterização ambiental e avaliação de cargas difusas da área de drenagem da represa Itupararranga, SP. [Master’s thesis in Hidráulica e Saneamento]. São Carlos: Universidade de São Paulo. https://doi.org/10.11606/D.18.2012.tde-30082012-150106
    » https://doi.org/10.11606/D.18.2012.tde-30082012-150106
  • Søndergaard, M., Jensen, J.P., & Jeppesen, E., 2003. Role of sediment and internal loading of phosphorus in shallow lakes. Hydrobiologia 506-509(1-3), 135-145. http://dx.doi.org/10.1023/B:HYDR.0000008611.12704.dd
    » http://dx.doi.org/10.1023/B:HYDR.0000008611.12704.dd
  • Straškraba, M., & Tundisi, J.G., 1999. Reservoir ecosystem functioning: theory and application. In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 565-583.
  • Straškraba, M., 1999. Retention time as a key variable of reservoir limnology. In: Tundisi, J.G., & Straškraba, M., eds. Theoretical reservoir ecology and its applications. Leiden: Backhuys, 385-410.
  • Teixeira, A.P., Bortolini, J.C., & Carneiro, F.M., 2022. Taxonomic and functional spatial distribution model of phytoplankton in tropical cascading reservoirs. Rev. Bras. Bot. Braz. J. Bot. 45(2), 791-805. http://dx.doi.org/10.1007/s40415-022-00810-7
    » http://dx.doi.org/10.1007/s40415-022-00810-7
  • Teodoru, C., & Wehrli, B., 2005. Retention of sediments and nutrients in the Iron Gate I Reservoir on the Danube River. Biogeochemistry 76(3), 539-565. http://dx.doi.org/10.1007/s10533-005-0230-6
    » http://dx.doi.org/10.1007/s10533-005-0230-6
  • Thornton, K.W., 1990. Sedimentary processes. In: Thornton, K.W., Kimmel, B.L., & Payne, F.E., eds. Reservoir limnology: ecological perspectives. New York: Wiley, 43-69.
  • Universidade Federal de São Carlos - UFSCar, 2008. Inventário limnológico para o monitoramento da qualidade da água e dos sedimentos do reservatório da UHE Itupararanga (Relatório Final). São Carlos: Fai-UFSCar.
  • Van Cappellen, P., & Maavara, T., 2016. Rivers in the anthropocene: global scale modifications of riverine nutrient fluxes by damming. Ecohydrol. Hydrobiol. 16(2), 106-111. http://dx.doi.org/10.1016/j.ecohyd.2016.04.001
    » http://dx.doi.org/10.1016/j.ecohyd.2016.04.001
  • Vo, N.X.Q., Doan, T.V., & Kang, H., 2014. Impoundments increase potential for phosphorus retention and remobilization in an urban stream. Environ. Eng. Res. 19(2), 175-184. http://dx.doi.org/10.4491/eer.2014.19.2.175
    » http://dx.doi.org/10.4491/eer.2014.19.2.175
  • Vollenweider, R.A., & Kerekes, J., 1982. Eutrophication of waters: monitoring, assessment and control. Paris: Organization for Economic Cooperation and Development.
  • Zhou, B., Xu, Y., Vogt, R.D., Lu, X., Li, X., Deng, X., Yue, A., & Zhu, L., 2016. Effects of land use change on phosphorus levels in surface waters - a case study of a watershed strongly influenced by agriculture. Water Air Soil Pollut. 227(5), 160. http://dx.doi.org/10.1007/s11270-016-2855-6
    » http://dx.doi.org/10.1007/s11270-016-2855-6

Edited by

Associate Editor: André Megali Amado.

Publication Dates

  • Publication in this collection
    20 Mar 2023
  • Date of issue
    2023

History

  • Received
    20 Oct 2022
  • Accepted
    16 Feb 2023
Associação Brasileira de Limnologia Av. 24 A, 1515, CEP: 13506-900 , Tel.:+55 (19) 3526-4225 - Rio Claro - SP - Brazil
E-mail: actalimno@gmail.com