Acessibilidade / Reportar erro

Impact of large reservoirs on simulated discharges of Brazilian rivers

O impacto de grandes reservatórios nas vazões simuladas de rios brasileiros

ABSTRACT

Tens of thousands of dams were built around the world to reduce flood risks, produce energy, and maximize benefits of limited freshwater resources. In Brazil, the main and largest reservoirs are related to hydropower plants. Improving the understanding of reservoir dynamics is important not only to evaluate their impact in the flow regime of Brazilian rivers, but also to simulate the combined effect of constructing new dams and potential alterations under future climatic conditions. Here, we analyze how an ideal representation of reservoirs in terms of forced discharge would improve a previously calibrated hydrological model under the Brazilian domain. We forced the continental-scale version of the MGB model on observed reservoir outflows from 109 hydropower dams, which are part of the Brazilian National Interconnected System controlled by the National Electrical System Operator. Model simulated flows were replaced by the reservoir outflows in all dam locations and were compared to the original discharge in downstream gauges. The forced discharge simulation presented a mean improvement for Kling-Gupta Efficiency of 21%, when compared to the original model (naturalized flow). This analysis is a preliminary step towards an explicit representation of the reservoirs in the model, what will be conducted in a future study.

Keywords:
Hydrologic modeling; Reservoirs; Large-scale hydrology; MGB, Brazil

RESUMO

Milhares de barragens foram construídas no mundo para diminuir os riscos de inundação, gerar energia, e maximizar os benefícios dos recursos limitados de água doce. No Brasil, os principais e maiores reservatórios são do setor elétrico. Melhorar a compreensão da dinâmica dos reservatórios é importante não só para avaliar o seu impacto no regime de vazões dos rios brasileiros, mas também para simular o efeito combinado da construção de novas barragens e potenciais alterações sob condições climáticas futuras. Neste estudo, analisa-se como uma representação ideal de reservatórios em termos de vazões substituídas melhoraria um modelo hidrológico previamente calibrado sobre o domínio brasileiro. Forçou-se a versão em escala continental do modelo MGB com defluências observadas de 109 usinas hidrelétricas, que fazem parte do Sistema Interligado Nacional (SIN) - controlados pelo Operador Nacional do Sistema Elétrico (ONS). A simulação com vazões substituídas apresentou um Índice de Melhoria no Desempenho médio (IMD) para a Eficiência de Kling-Gupta (KGE) de 21%, em relação ao modelo original (i.e., vazões naturalizadas). Esta análise é um passo preliminar para uma representação explícita dos reservatórios no modelo, que será conduzida em um estudo futuro.

Palavras-chave:
Modelo hidrológico; Reservatórios; Hidrologia de grande escala; MGB; Brasil

INTRODUCTION

The hydrological cycle of the Earth is changing due to anthropogenic alterations in water resources management and land uses and cover. In one hundred years, the population has quadrupled to 7 billion people, and while half of them live in urban areas, rapid urbanization is occurring in many developed and developing regions of the world (Klein Goldewijk et al., 2011Klein Goldewijk, K., Beusen, A., Van Drecht, G., & De Vos, M. (2011). The HYDE 3.1 spatially explicit database of human‐induced global land‐use change over the past 12,000 years. Global Ecology and Biogeography, 20(1), 73-86. http://dx.doi.org/10.1111/j.1466-8238.2010.00587.x.
http://dx.doi.org/10.1111/j.1466-8238.20...
). For food, the irrigated agricultural area increased six times in the last century (Siebert et al., 2015Siebert, S., Kummu, M., Porkka, M., Döll, P., Ramankutty, N., & Scanlon, B. R. (2015). A global data set of the extent of irrigated land from 1900 to 2005. Hydrology and Earth System Sciences, 19(3), 1521-1545. http://dx.doi.org/10.5194/hess-19-1521-2015.
http://dx.doi.org/10.5194/hess-19-1521-2...
) and the water demand continues to rise, driving an ever-increasing agricultural demand (>70% of all fresh water). In Brazil, irrigation represents about 50% of the total estimated abstractions and 70% of the total estimated water consumption (Agência Nacional de Águas, 2018Agência Nacional de Águas – ANA. (2018). Conjuntura recursos hídricos Brasil 2018. Brasília: ANA.).

In order to maximize the benefits of limited freshwater resources and mitigate flood risks, tens of thousands of artificial dams have been built in the world's major river systems (Nilsson et al., 2005Nilsson, C., Reidy, C. A., Dynesius, M., & Revenga, C. (2005). Fragmentation and flow regulation of the world’s large river systems. Science, 308(5720), 405-408. PMid: http://dx.doi.org/10.1126/science.1107887.
http://dx.doi.org/10.1126/science.110788...
; Lehner et al., 2011Lehner, B., Liermann, C. R., Revenga, C., Vörösmarty, C., Fekete, B., Crouzet, P., Döll, P., Endejan, M., Frenken, K., Magome, J., Nilsson, C., Robertson, J. C., Rödel, R., Sindorf, N., & Wisser, D. (2011). High‐resolution mapping of the world’s reservoirs and dams for sustainable river‐flow management. Frontiers in Ecology and the Environment, 9(9), 494-502. http://dx.doi.org/10.1890/100125.
http://dx.doi.org/10.1890/100125...
). Hydropower generation is the main use for the largest reservoirs in Brazil, corresponding to 65% of the country’s installed capacity (Empresa de Pesquisa Energética, 2018Empresa de Pesquisa Energética – EPE. (2018). Balanço energético nacional 2018: ano base 2017. Rio de Janeiro: EPE.). Due to the severe droughts that hit Brazil’s main basins between 2012 and 2015, the average share of hydroelectric power generation declined to 76% due to nation scale depletion of reservoirs, with great use of thermal plants and large reservoir depletion. Just before, between 2000 and 2012, artificial dams yielded an average of 91% of the power generation (Zambon, 2015Zambon, R. C. (2015). A operação dos reservatórios e o planejamento da operação hidrotérmica do Sistema Interligado Nacional. Revista USP, 104, 133-144.).

Despite of the national relevance of large reservoirs for energy production and reduction of flood risk, it is fundamental to understand how these major structural interventions change natural hydrological regimes. In general, reservoirs store water and release it downstream according to the current demand. Regulating water flow results in attenuated hydrographs, higher evaporation due to the increase in water surface area, and also major changes in sediment and nutrient transport (Ligon et al., 1995Ligon, F. K., Dietrich, W. E., & Trush, W. J. (1995). Downstream ecological effects of dams. Bioscience, 45(3), 183-192. http://dx.doi.org/10.2307/1312557.
http://dx.doi.org/10.2307/1312557...
). The existence of dams also leads to fragmentation of rivers, reducing the natural connectivity within and between fluvial systems (Tischendorf & Fahrig, 2000Tischendorf, L., & Fahrig, L. (2000). On the usage and measurement of landscape connectivity. Oikos, 90(1), 7-19. http://dx.doi.org/10.1034/j.1600-0706.2000.900102.x.
http://dx.doi.org/10.1034/j.1600-0706.20...
; Moilanen & Hanski, 2001Moilanen, A., & Hanski, I. (2001). On the use of connectivity measures in spatial ecology. Oikos, 95(1), 147-151. http://dx.doi.org/10.1034/j.1600-0706.2001.950116.x.
http://dx.doi.org/10.1034/j.1600-0706.20...
), among other impacts.

An interesting way to assess the impact of dams on river dynamics is by making use of conceptual and physically-based hydrological models. In this context, reservoirs can be simulated explicitly by incorporating them in the model structure (e.g., solving the water balance for storage with stage-area-volume curves and operation rules) or not explicitly by using a simple substitution form, in which the simulated flows are replaced by reservoir operation data at the corresponding dam locations (i.e., the hydrological model is forced with observed reservoir outflows). In Brazil, reservoirs have been simulated in hydrological models either by explicit or substitution approaches (e.g. Fleischmann et al., 2019bFleischmann, A., Collischonn, W., & Paiva, R. C. D. (2019b). Estimating design hydrographs at the basin scale: from event-based to continuous hydrological simulation. Revista Brasileira de Recursos Hídricos, 24, 1-10., 2019cFleischmann, A., Collischonn, W., Paiva, R. C. D., & Tucci, C. (2019c). Modeling the role of reservoirs versus floodplains on large-scale river hydrodynamics. Natural Hazards, 99(2), 1075-1104. http://dx.doi.org/10.1007/s11069-019-03797-9.
http://dx.doi.org/10.1007/s11069-019-037...
), but the existing studies usually have a focus in specific river basins or particular dams. A comprehensive analysis regarding the impact of reservoirs at the national scale, on the other hand, has not been carried out so far.

Advances in modeling techniques, computer capabilities and big data has taken place in the last years (Siqueira et al., 2018Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842.), and assessing reservoir impacts over large spatial domains has become even more a reality. For instance, in Global Hydrology and Water Resources models (GHWRM), which aim to assess the anthropogenic impacts over the entire globe (see Bierkens, 2015Bierkens, M. F. P. (2015). Global hydrology 2015: state, trends, and directions. Water Resources Research, 51(7), 4923-4947. http://dx.doi.org/10.1002/2015WR017173.
http://dx.doi.org/10.1002/2015WR017173...
), numerous generic algorithms have been developed to explicitly represent the influence of lakes and reservoirs on river flows (Döll et al., 2003Döll, P., Kaspar, F., & Lehner, B. (2003). A global hydrological model for deriving water availability indicators: model tuning and validation. Journal of Hydrology (Amsterdam), 270(1-2), 105-134. http://dx.doi.org/10.1016/S0022-1694(02)00283-4.
http://dx.doi.org/10.1016/S0022-1694(02)...
; Haddeland et al., 2006Haddeland, I., Skaugen, T., & Lettenmaier, D. P. (2006). Anthropogenic impacts on continental surface water fluxes. Geophysical Research Letters, 33(8), L08406. http://dx.doi.org/10.1029/2006GL026047.
http://dx.doi.org/10.1029/2006GL026047...
; Hanasaki et al., 2006Hanasaki, N., Kanae, S., & Oki, T. (2006). A reservoir operation scheme for global river routing models. Journal of Hydrology (Amsterdam), 327(1-2), 22-41. http://dx.doi.org/10.1016/j.jhydrol.2005.11.011.
http://dx.doi.org/10.1016/j.jhydrol.2005...
; Wada et al., 2014Wada, Y., Wisser, D., & Bierkens, M. F. P. (2014). Global modeling of withdrawal, allocation and consumptive use of surface water and groundwater resources. Earth System Dynamics, 5(1), 15-40. http://dx.doi.org/10.5194/esd-5-15-2014.
http://dx.doi.org/10.5194/esd-5-15-2014...
; Zajac et al., 2017Zajac, Z., Revilla-Romero, B., Salamon, P., Burek, P., Hirpa, F. A., & Beck, H. (2017). The impact of lake and reservoir parameterization on global streamflow simulation. Journal of Hydrology (Amsterdam), 548, 552-568. PMid:28649141. http://dx.doi.org/10.1016/j.jhydrol.2017.03.022.
http://dx.doi.org/10.1016/j.jhydrol.2017...
; Shin et al. 2018Shin, S., Pokhrel, Y., & Miguez‐Macho, G. (2018). High Resolution Modeling of Reservoir Release and Storage Dynamics at the Continental Scale. Water Resources Research). However, generic operation rules used by GHWRMs are rather simplistic and often rely on global datasets (e.g., Lehner et al., 2011Lehner, B., Liermann, C. R., Revenga, C., Vörösmarty, C., Fekete, B., Crouzet, P., Döll, P., Endejan, M., Frenken, K., Magome, J., Nilsson, C., Robertson, J. C., Rödel, R., Sindorf, N., & Wisser, D. (2011). High‐resolution mapping of the world’s reservoirs and dams for sustainable river‐flow management. Frontiers in Ecology and the Environment, 9(9), 494-502. http://dx.doi.org/10.1890/100125.
http://dx.doi.org/10.1890/100125...
), which may hinder the representation of complex reservoir systems such as the Brazilian National Interconnected System (SIN). Exploring regionally available dam outflow data using a hydrological model with national coverage would be an important step to overcome such limitations, helping further development for robust operation rules and giving insights on the potential benefits of including reservoirs in the model structure. The main Brazilian dataset is the Reservoir Monitoring System (Sistema de Acompanhamento de Reservatórios) from the National Water Agency (ANA, in the Portuguese acronym), which provides daily time series of information as water level, percentage of active storage, inflows and outflows.

MGB (Collischonn, 2001Collischonn, W. (2001). Simulação hidrológica de grandes bacias. (Doctoral dissertation). Universidade Federal do Rio Grande do Sul, Porto Alegre.; Collischonn et al., 2007Collischonn, W., Allasia, D., Silva, B. C., & Tucci, C. M. (2007). The MGB-IPH model for large-scale rainfall: runoff modelling. Hydrological Sciences Journal, 52(5), 878-895. http://dx.doi.org/10.1623/hysj.52.5.878.
http://dx.doi.org/10.1623/hysj.52.5.878...
; Paiva et al., 2013Paiva, R. C. D., Buarque, D. C., Collischonn, W., Bonnet, M. P., Frappart, F., Calmant, S., & Bulhões Mendes, C. A. (2013). Large‐scale hydrologic and hydrodynamic modeling of the Amazon River basin. Water Resources Research, 49(3), 1226-1243. http://dx.doi.org/10.1002/wrcr.20067.
http://dx.doi.org/10.1002/wrcr.20067...
) is a large scale hydrological model which was already applied to the whole South America in a multi-basin manner with satisfactory performance (Siqueira et al., 2018Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842.). Although some studies have explored the explicit simulation of reservoirs in MGB (Larentis et al., 2008Larentis, D., Collischonn, W., & Tucci, C. E. M. (2008). Simulação da Qualidade de Água em Grandes Bacias: rio Taquari-Antas, RS. Revista Brasileira de Recursos Hídricos, 13(3), 5-22. http://dx.doi.org/10.21168/rbrh.v13n3.p5-22.
http://dx.doi.org/10.21168/rbrh.v13n3.p5...
; Collischonn et al., 2011Collischonn, B., Paiva, R. C. D., Collischonn, W., Meirelles, F. S. C., Schettini, E. B. C., & Fan, F. M. (2011). Modelagem hidrológica de uma bacia com uso intensivo de água: Caso do Rio Quaraí-RS. Revista Brasileira de Recursos Hídricos, 16(4), 119-134.; Fleischmann et al., 2015Fleischmann, A. S., Siqueira, V. A., Collischonn, W., & Fan, F. M. (2015, November 22-27). Desenvolvimento do módulo de reservatórios do modelo hidrológico MGB. In Associação Brasileira de Recursos Hídricos, Anais do 21º Simpósio Brasileiro de Recursos Hídricos. Brasília, Brasil,., 2019bFleischmann, A., Collischonn, W., & Paiva, R. C. D. (2019b). Estimating design hydrographs at the basin scale: from event-based to continuous hydrological simulation. Revista Brasileira de Recursos Hídricos, 24, 1-10., 2019cFleischmann, A., Collischonn, W., Paiva, R. C. D., & Tucci, C. (2019c). Modeling the role of reservoirs versus floodplains on large-scale river hydrodynamics. Natural Hazards, 99(2), 1075-1104. http://dx.doi.org/10.1007/s11069-019-03797-9.
http://dx.doi.org/10.1007/s11069-019-037...
), reservoir modeling is not yet consolidated in MGB. Then, the aim of this work is to evaluate how an ideal representation of SIN reservoirs in terms of forced discharge would improve a continental-scale, previously calibrated hydrological model.

The paper is divided as follows: firstly in the Material and Methods section, the model is described, together with reservoir data and performance evaluation criteria. Then results are shown for an exploratory analysis of the reservoirs behavior, followed by the model with forced discharges, and finally for the differences in discharges between natural and forced simulation scenarios.

MATERIAL AND METHODS

Model description

MGB is a conceptual, semi-distributed hydrological model designed to simulate land surface hydrological processes and flow propagation in tropical basins (Collischonn et al., 2007Collischonn, W., Allasia, D., Silva, B. C., & Tucci, C. M. (2007). The MGB-IPH model for large-scale rainfall: runoff modelling. Hydrological Sciences Journal, 52(5), 878-895. http://dx.doi.org/10.1623/hysj.52.5.878.
http://dx.doi.org/10.1623/hysj.52.5.878...
). It divides the basin into unit-catchments to capture the spatial variability of both meteorological and physical characteristics and further into hydrological response units (HRUs) to represent land use and soil heterogeneities. Water and energy budget are computed at the HRU level considering the soil as a bucket model. Evapotranspiration (from soil, plant and canopy) is calculated using the Penman-Monteith equation. Surface, subsurface and groundwater runoff generated within the unit-catchment are propagated to the stream network using linear reservoirs to represent the catchment damping effect. Flow routing in river channels can be computed using different approaches, ranging from simple methods such as the Muskingum-Cunge (Collischonn et al., 2007Collischonn, W., Allasia, D., Silva, B. C., & Tucci, C. M. (2007). The MGB-IPH model for large-scale rainfall: runoff modelling. Hydrological Sciences Journal, 52(5), 878-895. http://dx.doi.org/10.1623/hysj.52.5.878.
http://dx.doi.org/10.1623/hysj.52.5.878...
) to more complex ones like the full 1D St. Venant equations (Paiva et al., 2013Paiva, R. C. D., Buarque, D. C., Collischonn, W., Bonnet, M. P., Frappart, F., Calmant, S., & Bulhões Mendes, C. A. (2013). Large‐scale hydrologic and hydrodynamic modeling of the Amazon River basin. Water Resources Research, 49(3), 1226-1243. http://dx.doi.org/10.1002/wrcr.20067.
http://dx.doi.org/10.1002/wrcr.20067...
).

The continental MGB model version for South America (MGB-SA) (Siqueira et al., 2018)Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842. is used in this study. It uses the inertial equation of propagation in rivers described in Pontes et al., (2017)Pontes, P. R. M., Fan, F. M., Fleischmann, A. S., de Paiva, R. C. D., Buarque, D. C., Siqueira, V. A., Jardim, P. F., Sorribas, M. V., & Collischonn, W. (2017). MGB-IPH model for hydrological and hydraulic simulation of large floodplain river systems coupled with open source GIS. Environmental Modelling & Software, 94, 1-20. http://dx.doi.org/10.1016/j.envsoft.2017.03.029.
http://dx.doi.org/10.1016/j.envsoft.2017...
, with sufficient physical basis to represent river hydrodynamics and floodplains effects. This version has 33,749 unit-catchment with river reaches at every 15 km and adopts the Multi-Source Weighted Ensemble Precipitation (MSWEP) as the forcing data, a product that combines estimates of rainfall by satellite, reanalysis and observed precipitation (Beck et al., 2017Beck, H. E., van Dijk, A. I. J. M., Levizzani, V., Schellekens, J., Miralles, D. G., Martens, B., & de Roo, A. (2017). 3-hourly 0.25° global gridded precipitation (1979–2015) by merging gauge, satellite, and reanalysis data. Hydrology and Earth System Sciences, 21(1), 589-615. http://dx.doi.org/10.5194/hess-21-589-2017.
http://dx.doi.org/10.5194/hess-21-589-20...
). The model was manually calibrated for the period between 1990 and 2010 with >600 discharge gauge stations (including reservoir naturalized flows from ONS, i.e., without the effect of regulation), and was validated using multiple remote sensing databases. Performances for daily flow regarding the Kling-Gupta and Nash-Sutcliffe efficiencies were > 0.6, respectively, for 70% and 55% of the gauges (Siqueira et al., 2018Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842.). Here we used the same simulation period as the original model version, from 01/01/1990 to 31/12/2009.

Reservoir data and substitution of flows in the model

To evaluate the impact of reservoir regulation on water discharges, the model was forced with observed outflows from reservoirs controlled by the ONS. For that, simulated discharges at unit-catchments corresponding to reservoir locations were replaced by the observed outflows at each model time step (daily), and this simulation was called the F simulation. Reservoirs that meet one or more of the following criteria were not considered in this study: a) built near the end simulation date of the model; b) located in the same unit-catchment and upstream of another reservoir, since in the approach used here only one reservoir per unit-catchment is possible (the most downstream reservoir was simulated); c) outside the simulated spatial domain; and d) with insufficient flow data. Thus, it was included discharge data from 109 reservoirs (out of an original dataset of 150 reservoirs), obtained from the Reservoir Monitoring System (SAR, in the Portuguese acronym). Figure 1 shows the name and location of all reservoirs used in this study. Only the impact on the flows was analyzed.

Figure 1
Location of the 109 dams considered in this study. The labeled reservoirs are those whose results will be discussed in the text.

We carried out an exploratory analysis of the reservoir SAR data by a visual inspection, and also computed the correlation between the outflows and other variables: inflows, storage (water level and volume) and the day of the year (Julian day). For the latter case, the correlation was computed in terms of circular statistics to avoid the discontinuity of days 365 and 1.

Evaluation of model performance

The F simulation and naturalized scenarios (i.e., the naturalized flow simulation without reservoirs, called here N simulation) were compared with in situ observed flow data obtained from the National Water Agency (ANA) monitoring network.

The Kling-Gupta efficiency indexes (KGE, Gupta et al., 2009Gupta, H. V., Kling, H., Yilmaz, K. K., & Martinez, G. F. (2009). Decomposition of the mean squared error and NSE performance criteria: implications for improving hydrological modelling. Journal of Hydrology (Amsterdam), 377(1-2), 80-91. http://dx.doi.org/10.1016/j.jhydrol.2009.08.003.
http://dx.doi.org/10.1016/j.jhydrol.2009...
) and Nash-Sutcliffe (NSE, Nash & Sutcliffe, 1970Nash, J. E., & Sutcliffe, J. V. (1970). River flow forecasting through conceptual models part I: A discussion of principles. Journal of Hydrology (Amsterdam), 10(3), 282-290. http://dx.doi.org/10.1016/0022-1694(70)90255-6.
http://dx.doi.org/10.1016/0022-1694(70)9...
), correlation (r) and volume error (bias) were also used. Another metric used by Zajac et al. (2017)Zajac, Z., Revilla-Romero, B., Salamon, P., Burek, P., Hirpa, F. A., & Beck, H. (2017). The impact of lake and reservoir parameterization on global streamflow simulation. Journal of Hydrology (Amsterdam), 548, 552-568. PMid:28649141. http://dx.doi.org/10.1016/j.jhydrol.2017.03.022.
http://dx.doi.org/10.1016/j.jhydrol.2017...
is the skill score (IMD, in the Portuguese acronym). A positive IMD means that the performance of the model with reservoirs has improved, while negative values indicate the opposite trend. IMD can be applied to any performance index; an example is given in Equation 1 by computing IMD for the KGE.

I M D K G E = K G E F K G E N K G E o p t K G E N (1)

The skill score for the KGE (IMDKGE) is a function of the KGE for the reservoirs simulation (KGEF), the KGE for the naturalized flows simulation (KGEN), and the optimum KGE (KGEopt= 1).

Changes in flows were analyzed with the following metrics: normalized root mean square error (nRMSE – i.e., RMSE divided by average discharge), correlation (r), standard deviation, and absolute differences in Q10 and Q90, between the N and F simulation. The average absolute percentage difference between N and F simulations is computed annually for the main basins of South America (Amazon, Tocantins, São Francisco and Paraná).

RESULTS AND DISCUSSION

Exploratory analysis of the reservoirs’ behavior

The average behavior of some variables throughout the year was investigated. Three dams were chosen to cover different regions and volume capacity (Figure 1). Based on Figure 2b it is possible to infer that, for the Porto Estrela reservoir, there is a period with lower flow release between May and mid-October. For all analyzed reservoirs, the outflow can vary expressively for the same day, in different years. The dispersion is visually large, for both outflow and inflow (Figure 2a), due to the climate variability. The correlations were in average 0.2, indicating a certain relation between the day of the year and the outflow.

Figure 2
(a, b) daily inflow/outflow during the year, for three Brazilian hydroelectric reservoirs after their constructions. The average for all the SAR reservoir’s data is shown in red, all observations in blue; in the title, the reservoir’s name and the correlation between the inflow/outflow and the day of the year. Some points are not shown to better visualization of the pattern; (c) outflow versus storage (water elevation) of three Brazilian hydroelectric reservoirs, comparing the years seasons.

One can conclude the same from the outflow to the inflow. The largest difference is in the Sobradinho reservoir correlation. Its correlation for inflow is 0.79; and for outflow, 0.16, what shows its great storage capacity and the impact on downstream flow regime. In seasonal basins, a greater correlation between inflow and season is expected.

Regarding the relationship between daily water elevation and outflow (Figure 2c), it is hard to deduce any operating rule based on its behavior. That was expected because ONS uses a complex operation rule, not based only in the reservoir’s water elevation. Moreover, the analysis performed here considered only daily data; a future assessment based on monthly data could lead to different conclusions. In the figure, it is observed that the reservoirs of Porto Estrela and Sobradinho have their highest discharge values in summer and fall, while 14 de Julho, a run-of-the-river plant, in winter and spring.

The analysis of the relation between outflow and possible driving factors (e.g. storage, inflow, time of year) can be studied through the correlation coefficient (Figure 3). It was done for 109 reservoirs. The variable that most explains the outflow is the inflow. The mean correlation is 0.79. The second most relevant factor is storage, represented here by volume and water elevation, with mean correlations of 0.31 and 0.22, respectively. The circular correlation between outflow and day of the year presented a mean correlation of 0.19, being the least important factor.

Figure 3
Distribution of the correlation between outflow and four other variables

Forced Discharge Simulation

Figure 4 presents a spatial assessment of the IMD metric for KGE index. In most gauge stations (≅ 74%) the performance improved. A poorer performance was observed along the southeast and extreme west coast of Brazil (Paraíba do Sul and Alto Paraguai basins). This may be caused by uncertainties related to input data in all the domain. Furthermore, since the MGB-SA was calibrated for the situation without reservoirs and, although care has been taken to use natural observed discharge and gauge stations without reservoirs influence, the calibrated parameters (such as the soil storage capacity parameter or the relationship between storage and saturation in the soil parameter) may have been adjusted in the natural scenario with the effect of the reservoirs on the outflow. That is, the effect of the reservoirs in the simulation was possibly compensated during model parameterization. After forcing the model with a reservoir outflow, it would be necessary to recalibrate the model to remove this compensation. This is the main reason we may have lower results in the forced discharge scenario at some gauge stations.

Figure 4
Skill Score (IMD) for KGE in the forced discharge simulation (F). IMD ≈ 0 is when both simulations presented similar results. Black lines are the countries. Thicker drainage network with density > 30 000 km2. Less thick drainage network with density > 8000 km2.

Due to MGB-SA inability to perform well in the Brazilian Northeast (semi-arid climate, Kottek et al., 2006Kottek, M., Grieser, J., Beck, C., Rudolf, B., & Rubel, F. (2006). World map of the Köppen-Geiger climate classification updated. Meteorologische Zeitschrift (Berlin), 15(3), 259-263. http://dx.doi.org/10.1127/0941-2948/2006/0130.
http://dx.doi.org/10.1127/0941-2948/2006...
), there was a considerable improvement in the results downstream of São Francisco river, mainly due to the consideration of Sobradinho and Luiz Gonzaga reservoirs, close to the river mouth. The Três Marias reservoir discharge also has an important role there. Most of the discharge in São Francisco River is generated in these south headwaters, since the downstream reaches are located in a more arid region. The impact of Brazilian reservoirs on discharge can be seen even outside Brazil, affecting other countries such as Paraguay and Argentina, located in downstream parts of the Paraguay and Paraná rivers (Bonetto et al., 1987Bonetto, A. A., Castello, H. P., & Wais, I. R. (1987). Stream regulation in Argentina, including the superior Paraná and Paraguay rivers. Regulated Rivers: Research and Management, 1(2), 129-143. http://dx.doi.org/10.1002/rrr.3450010204.
http://dx.doi.org/10.1002/rrr.3450010204...
).

There was no clear relationship between IMD and reservoirs volumes (Figure 5). For volumes larger than 1000 hm3 there seems to be a trend towards IMD improvement with increasing volume. The F simulation showed improvement in 80% of the stations and, of these, 40% improved more than 20% (IMD> 0.2). Hydrographs downstream of three selected reservoirs are presented in Figure 6 to compare the F and N simulations (see Figure 1 for dams’ location). The reservoir of Irapé (Jequitinhonha River) presented significant improvement in the hydrograph, whereas, for the others, the model with naturalized flows was satisfactory. The 14 de Julho dam in Taquari-Antas River is a run-of-the-river plant, which almost does not alter the downstream flow regime. Great improvement is observed in the base flow downstream of Manso reservoir in the Cuiabá River, which was underestimated in the naturalized flows model.

Figure 5
IMD x Reservoir capacity. Only the first station downstream of each dam was represented here. In red, the reservoirs discussed in the text.
Figure 6
Hydrographs of the simulations with forced flows (F, in blue) and the original model (with naturalized flows, N, in red) of the MGB-SA. In black the ANA (Agência Nacional de Águas, 2018Agência Nacional de Águas – ANA. (2018). Conjuntura recursos hídricos Brasil 2018. Brasília: ANA.) observed flow gauge stations (UHE 14 de Julho Alça - 86470900; UHE Irapé Jusante - 1642044; Rosário Oeste - 1456008). The metrics in the titles correspond to the simulation with forced flows.

The calculated metrics for discharge gauge stations downstream of the 109 dams were summarized by the 10th percentile, median and 90th percentile (Table 1). F scenario had a performance superior than scenario N. The metric with the greatest improvement was the logarithm of the NSE (due to more uncertainties in the model base flow), followed by the NSE, KGE and r. The KGE showed an improvement of 21% for the F simulation.

Table 1
The 10th percentile (P10), 90th percentile (P90) and median of the performance indices for the two simulations (F and N); below, the skill score (IMD) for the F simulation (141 gauge stations). In grey the cell where the simulation improved.

Each metric depicts different flow signatures. The correlation reflects the seasonality, but is unable to evaluate the average flow. A NSE equal to zero shows that the simulations are as accurate as the mean of the observed data, and negative NSE means that the mean is a better estimate than the model. To further improve the assessment of hydrological models, the KGE was created (Gupta et al., 2009Gupta, H. V., Kling, H., Yilmaz, K. K., & Martinez, G. F. (2009). Decomposition of the mean squared error and NSE performance criteria: implications for improving hydrological modelling. Journal of Hydrology (Amsterdam), 377(1-2), 80-91. http://dx.doi.org/10.1016/j.jhydrol.2009.08.003.
http://dx.doi.org/10.1016/j.jhydrol.2009...
), which takes into account the correlation, bias and variability. This way, it is ensured that bias and variability are not correlated.

Differences in discharges between natural and forced simulation scenarios

Changes in daily discharges between simulations F and N were evaluated through the following metrics: nRMSE, and relative mean absolute differences in Q10 and Q90 (Figure 7a, b and c, respectively). The difference between both scenarios must be analyzed considering that MGB-SA estimates present relevant biases in certain basins, which is registered in the supplementary material of Siqueira et al. (2018)Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842. and is summarized in Figure 7d for the 109 dams. In this case, we computed the bias between scenario N and the observed outflow for each dam. For a future study, the impact of reservoirs on river discharges may be analyzed by firstly performing a bias correction on the model naturalized flows.

Figure 7
Maps showing the alteration in calculated flow regime in river reaches due to the presence of reservoirs: a) nRMSE; b) Q90; c) Q10. Results should be interpreted considering the MGB-SA bias (Figure d). Light gray lines represent the country’s borders.

The nRMSE calculated exceeded 80% in the lower part of the São Francisco River, while MGB-SA bias is lower than 50% for most gauges there (Figure 7d), indicating that although the model has a large bias in such semi-arid regions, important flow regulation is performed by the reservoirs. The Paraná River presents nRMSE higher than 60% in its upstream portion, nRMSE between 20 and 40% until the Brazilian border, and less than 20% in its downstream reaches. Compared to the Paraná River MGB-SA bias (smaller than 10% for most dams, Figure 7d), nRMSE indicates that considerable portions of the basin upstream from the Brazilian border are under flow regulation. The Amazon basin is little affected by reservoirs in terms of flow.

In general, the difference in Q10 was less than 20%, while the difference in Q90 was greater for the Paraná River and its tributaries (20-40%). However, we cannot say that this difference is due to the reservoirs’ presence, because the MGB-SA flow bias is, on average, greater than 20%.

CONCLUSIONS

The aim of this study was to evaluate how an ideal representation of reservoirs in terms of forced discharge would improve a continental-scale and previously calibrated hydrological model. For this purpose, we considered outflows of 109 SIN reservoirs in the MGB-SA model (Siqueira et al., 2018Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842.).

Through an exploratory analysis, it is identified how observed reservoirs outflows are related to inflows and water elevation. As most reservoirs are run-of-the-river (64 out of 109), there is a strong relation between inflow and outflow, however, no expressive correlation is observed between outflow and water storage. Thus, it is not trivial to build a general operation rule for all 109 SIN reservoirs.

The simulation with forced discharge (F simulation) presented better performance than the MGB-SA with naturalized flows (N simulation), showing a 21% improvement in KGE. Large rivers as São Francisco and Paraná were shown to be largely regulated, even considering the MGB-SA bias. In a next step, we intend to explicitly represent the reservoirs in the model and test operating rules. It is possible to develop a model of water uses, and then use the output schemes developed to meet the demand of water downstream. In addition, level-volume relationships provided by the ONS and the new ANA water use products can be used (Agência Nacional de Águas, 2018Agência Nacional de Águas – ANA. (2018). Conjuntura recursos hídricos Brasil 2018. Brasília: ANA.).

ACKNOWLEDGMENTS

The present work was carried out with the support of CNPq (National Council of Scientific and Technological Development /Brazil), under the grant numbers 132663/2017-1 and 141161/2017-5. The authors also acknowledges IPH and UFRGS for their support, and the anonymous reviewer’s suggestions that significantly improve the text.

REFERENCES

  • Agência Nacional de Águas – ANA. (2018). Conjuntura recursos hídricos Brasil 2018 Brasília: ANA.
  • Beck, H. E., van Dijk, A. I. J. M., Levizzani, V., Schellekens, J., Miralles, D. G., Martens, B., & de Roo, A. (2017). 3-hourly 0.25° global gridded precipitation (1979–2015) by merging gauge, satellite, and reanalysis data. Hydrology and Earth System Sciences, 21(1), 589-615. http://dx.doi.org/10.5194/hess-21-589-2017
    » http://dx.doi.org/10.5194/hess-21-589-2017
  • Bierkens, M. F. P. (2015). Global hydrology 2015: state, trends, and directions. Water Resources Research, 51(7), 4923-4947. http://dx.doi.org/10.1002/2015WR017173
    » http://dx.doi.org/10.1002/2015WR017173
  • Bonetto, A. A., Castello, H. P., & Wais, I. R. (1987). Stream regulation in Argentina, including the superior Paraná and Paraguay rivers. Regulated Rivers: Research and Management, 1(2), 129-143. http://dx.doi.org/10.1002/rrr.3450010204
    » http://dx.doi.org/10.1002/rrr.3450010204
  • Collischonn, B., Paiva, R. C. D., Collischonn, W., Meirelles, F. S. C., Schettini, E. B. C., & Fan, F. M. (2011). Modelagem hidrológica de uma bacia com uso intensivo de água: Caso do Rio Quaraí-RS. Revista Brasileira de Recursos Hídricos, 16(4), 119-134.
  • Collischonn, W. (2001). Simulação hidrológica de grandes bacias (Doctoral dissertation). Universidade Federal do Rio Grande do Sul, Porto Alegre.
  • Collischonn, W., Allasia, D., Silva, B. C., & Tucci, C. M. (2007). The MGB-IPH model for large-scale rainfall: runoff modelling. Hydrological Sciences Journal, 52(5), 878-895. http://dx.doi.org/10.1623/hysj.52.5.878
    » http://dx.doi.org/10.1623/hysj.52.5.878
  • Döll, P., Kaspar, F., & Lehner, B. (2003). A global hydrological model for deriving water availability indicators: model tuning and validation. Journal of Hydrology (Amsterdam), 270(1-2), 105-134. http://dx.doi.org/10.1016/S0022-1694(02)00283-4
    » http://dx.doi.org/10.1016/S0022-1694(02)00283-4
  • Empresa de Pesquisa Energética – EPE. (2018). Balanço energético nacional 2018: ano base 2017 Rio de Janeiro: EPE.
  • Fleischmann, A. S., Siqueira, V. A., Collischonn, W., & Fan, F. M. (2015, November 22-27). Desenvolvimento do módulo de reservatórios do modelo hidrológico MGB. In Associação Brasileira de Recursos Hídricos, Anais do 21º Simpósio Brasileiro de Recursos Hídricos Brasília, Brasil,.
  • Fleischmann, A., Collischonn, W., & Paiva, R. C. D. (2019b). Estimating design hydrographs at the basin scale: from event-based to continuous hydrological simulation. Revista Brasileira de Recursos Hídricos, 24, 1-10.
  • Fleischmann, A., Collischonn, W., Paiva, R. C. D., & Tucci, C. (2019c). Modeling the role of reservoirs versus floodplains on large-scale river hydrodynamics. Natural Hazards, 99(2), 1075-1104. http://dx.doi.org/10.1007/s11069-019-03797-9
    » http://dx.doi.org/10.1007/s11069-019-03797-9
  • Gupta, H. V., Kling, H., Yilmaz, K. K., & Martinez, G. F. (2009). Decomposition of the mean squared error and NSE performance criteria: implications for improving hydrological modelling. Journal of Hydrology (Amsterdam), 377(1-2), 80-91. http://dx.doi.org/10.1016/j.jhydrol.2009.08.003
    » http://dx.doi.org/10.1016/j.jhydrol.2009.08.003
  • Haddeland, I., Skaugen, T., & Lettenmaier, D. P. (2006). Anthropogenic impacts on continental surface water fluxes. Geophysical Research Letters, 33(8), L08406. http://dx.doi.org/10.1029/2006GL026047
    » http://dx.doi.org/10.1029/2006GL026047
  • Hanasaki, N., Kanae, S., & Oki, T. (2006). A reservoir operation scheme for global river routing models. Journal of Hydrology (Amsterdam), 327(1-2), 22-41. http://dx.doi.org/10.1016/j.jhydrol.2005.11.011
    » http://dx.doi.org/10.1016/j.jhydrol.2005.11.011
  • Klein Goldewijk, K., Beusen, A., Van Drecht, G., & De Vos, M. (2011). The HYDE 3.1 spatially explicit database of human‐induced global land‐use change over the past 12,000 years. Global Ecology and Biogeography, 20(1), 73-86. http://dx.doi.org/10.1111/j.1466-8238.2010.00587.x
    » http://dx.doi.org/10.1111/j.1466-8238.2010.00587.x
  • Kottek, M., Grieser, J., Beck, C., Rudolf, B., & Rubel, F. (2006). World map of the Köppen-Geiger climate classification updated. Meteorologische Zeitschrift (Berlin), 15(3), 259-263. http://dx.doi.org/10.1127/0941-2948/2006/0130
    » http://dx.doi.org/10.1127/0941-2948/2006/0130
  • Larentis, D., Collischonn, W., & Tucci, C. E. M. (2008). Simulação da Qualidade de Água em Grandes Bacias: rio Taquari-Antas, RS. Revista Brasileira de Recursos Hídricos, 13(3), 5-22. http://dx.doi.org/10.21168/rbrh.v13n3.p5-22
    » http://dx.doi.org/10.21168/rbrh.v13n3.p5-22
  • Lehner, B., Liermann, C. R., Revenga, C., Vörösmarty, C., Fekete, B., Crouzet, P., Döll, P., Endejan, M., Frenken, K., Magome, J., Nilsson, C., Robertson, J. C., Rödel, R., Sindorf, N., & Wisser, D. (2011). High‐resolution mapping of the world’s reservoirs and dams for sustainable river‐flow management. Frontiers in Ecology and the Environment, 9(9), 494-502. http://dx.doi.org/10.1890/100125
    » http://dx.doi.org/10.1890/100125
  • Ligon, F. K., Dietrich, W. E., & Trush, W. J. (1995). Downstream ecological effects of dams. Bioscience, 45(3), 183-192. http://dx.doi.org/10.2307/1312557
    » http://dx.doi.org/10.2307/1312557
  • Moilanen, A., & Hanski, I. (2001). On the use of connectivity measures in spatial ecology. Oikos, 95(1), 147-151. http://dx.doi.org/10.1034/j.1600-0706.2001.950116.x
    » http://dx.doi.org/10.1034/j.1600-0706.2001.950116.x
  • Nash, J. E., & Sutcliffe, J. V. (1970). River flow forecasting through conceptual models part I: A discussion of principles. Journal of Hydrology (Amsterdam), 10(3), 282-290. http://dx.doi.org/10.1016/0022-1694(70)90255-6
    » http://dx.doi.org/10.1016/0022-1694(70)90255-6
  • Nilsson, C., Reidy, C. A., Dynesius, M., & Revenga, C. (2005). Fragmentation and flow regulation of the world’s large river systems. Science, 308(5720), 405-408. PMid: http://dx.doi.org/10.1126/science.1107887
    » http://dx.doi.org/10.1126/science.1107887
  • Paiva, R. C. D., Buarque, D. C., Collischonn, W., Bonnet, M. P., Frappart, F., Calmant, S., & Bulhões Mendes, C. A. (2013). Large‐scale hydrologic and hydrodynamic modeling of the Amazon River basin. Water Resources Research, 49(3), 1226-1243. http://dx.doi.org/10.1002/wrcr.20067
    » http://dx.doi.org/10.1002/wrcr.20067
  • Pontes, P. R. M., Fan, F. M., Fleischmann, A. S., de Paiva, R. C. D., Buarque, D. C., Siqueira, V. A., Jardim, P. F., Sorribas, M. V., & Collischonn, W. (2017). MGB-IPH model for hydrological and hydraulic simulation of large floodplain river systems coupled with open source GIS. Environmental Modelling & Software, 94, 1-20. http://dx.doi.org/10.1016/j.envsoft.2017.03.029
    » http://dx.doi.org/10.1016/j.envsoft.2017.03.029
  • Shin, S., Pokhrel, Y., & Miguez‐Macho, G. (2018). High Resolution Modeling of Reservoir Release and Storage Dynamics at the Continental Scale. Water Resources Research
  • Siebert, S., Kummu, M., Porkka, M., Döll, P., Ramankutty, N., & Scanlon, B. R. (2015). A global data set of the extent of irrigated land from 1900 to 2005. Hydrology and Earth System Sciences, 19(3), 1521-1545. http://dx.doi.org/10.5194/hess-19-1521-2015
    » http://dx.doi.org/10.5194/hess-19-1521-2015
  • Siqueira, V. A., Paiva, R. C. D. D., Fleischmann, A. S., Fan, F. M., Ruhoff, A. L., Pontes, P. R. M., Collischonn, W. (2018). Toward continental hydrologic–hydrodynamic modeling in South America. Hydrology and Earth System Sciences, 22(9), 4815-4842.
  • Tischendorf, L., & Fahrig, L. (2000). On the usage and measurement of landscape connectivity. Oikos, 90(1), 7-19. http://dx.doi.org/10.1034/j.1600-0706.2000.900102.x
    » http://dx.doi.org/10.1034/j.1600-0706.2000.900102.x
  • Wada, Y., Wisser, D., & Bierkens, M. F. P. (2014). Global modeling of withdrawal, allocation and consumptive use of surface water and groundwater resources. Earth System Dynamics, 5(1), 15-40. http://dx.doi.org/10.5194/esd-5-15-2014
    » http://dx.doi.org/10.5194/esd-5-15-2014
  • Zajac, Z., Revilla-Romero, B., Salamon, P., Burek, P., Hirpa, F. A., & Beck, H. (2017). The impact of lake and reservoir parameterization on global streamflow simulation. Journal of Hydrology (Amsterdam), 548, 552-568. PMid:28649141. http://dx.doi.org/10.1016/j.jhydrol.2017.03.022
    » http://dx.doi.org/10.1016/j.jhydrol.2017.03.022
  • Zambon, R. C. (2015). A operação dos reservatórios e o planejamento da operação hidrotérmica do Sistema Interligado Nacional. Revista USP, 104, 133-144.

Publication Dates

  • Publication in this collection
    30 Mar 2020
  • Date of issue
    2020

History

  • Received
    20 June 2019
  • Reviewed
    18 Nov 2019
  • Accepted
    27 Dec 2019
Associação Brasileira de Recursos Hídricos Av. Bento Gonçalves, 9500, CEP: 91501-970, Tel: (51) 3493 2233, Fax: (51) 3308 6652 - Porto Alegre - RS - Brazil
E-mail: rbrh@abrh.org.br