Acessibilidade / Reportar erro

Streamflow regionalization for the Mortes River Basin upstream from the Funil Hydropower Plant, MG

Regionalização de vazões para a bacia hidrográfica do Rio das Mortes a montante da Usina Hidrelétrica do Funil, MG

Abstract

Maximum and minimum streamflow are fundamental for water resource management, especially for water rights. However, lack of monitoring and scarce streamflow data limit such studies. Streamflow regionalization is a useful tool to overcome these limitations. The study developed models for regionalization of the maximum and minimum reference streamflows for the Mortes River Basin (MRB) (Water Resources Planning and Management Unit - GD2), Southern Minas Gerais State. The study used long-term streamflow historical series provided by the Brazilian National Water Agency (ANA). Previous exploratory analysis was performed, and it was observed that the streamflow series are stationary according to the Mann-Kendall test. The estimation of the streamflow for different return periods (RP) was performed by fitting Probability Density Functions (PDFs) that were tested by the Anderson-Darling (AD) test. The Generalized Extreme Values (GEV) and Wakeby were the most appropriate PDFs for maximum and minimum streamflows, respectively. The streamflow models were fitted using a power regression procedure, considering the drainage area of the watersheds as inputs. The fittings reached the coefficient of determination (R2) greater than 0.90. Thus, the streamflow regionalization models demonstrated good performance and are a potential tool to be used for water resource management in the studied basin.

Keywords:
probability density functions; statistical hydrology; water resource management

Resumo

Vazões máximas e mínimas são fundamentais para a gestão dos recursos hídricos, especialmente para outorga de uso da água. Entretanto, a falta de monitoramento e a escassez de dados de vazão são fatores limitantes para tais estudos. Para contornar tal dificuldade, uma das ferramentas mais utilizadas é a regionalização de vazões. O objetivo neste estudo foi desenvolver modelos de regionalização para vazões máximas e mínimas de referênciapara a bacia hidrográfica do Rio das Mortes (MRB) (Unidade de Planejamento e Gestão dos Recursos Hídricos - GD2), sul do estado de Minas Gerais. As séries históricas de vazão de longo termo utilizadas foram disponibilizada pela Agência Nacional das Águas (ANA). Uma análise exploratório preliminar foi realizada, bem como a análise da estacionariedade das séries por meio do teste de tendências de Mann-Kendall. As Funções Densidade de Probabilidade (PDFs) com diferentes números de parâmetros foram utilizadas para estimar as vazões mínimas e máximas associadas a diferentes tempos de retorno e a performance foi avaliada por meio do teste de Anderson-Darling. A Generalized Extreme Values e Wakeby foram as distribuições mais apropriadas para vazões máximas e mínimas, respectivamente. O modelo de regressão para vazões foi ajustado utilizando como variável independente a área de drenagem das estações fluviométricas, sendo que os coeficientes de determinação obtidos foram superiores a 0,90. Desta forma, os modelos de regionalização demonstraram uma boa performance e são uma ferramenta potencial para ser utilizada na gestão dos recursos hídricos.

Palavras-chave:
funções densidade de probabilidade; gestão dos recursos hídricos; hidrologia estatística

1. INTRODUCTION

To perform studies regarding floods and droughts, analyses of different scenarios of water uses, and the streamflow prediction, among other possible applications, it is important to understand the behavior of streamflow in basins as accurately as possible (Fan et al., 2014FAN, F. M.; PONTES, P.; PAIVA, R.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. https://dx.doi.org/10.21168/rbrh.v19n4.p137-147
https://dx.doi.org/10.21168/rbrh.v19n4.p...
). This understanding may be performed by analyzing the minimum reference streamflows and the frequency (or return period) of the historical floods, using historical series.

Despite the importance of these studies for water resource management, streamflow data have often been a limiting factor for most of the region of interest (Zaman et al., 2012ZAMAN, M. A.; RAHMAN, A.; HADDAD, K. Regional flood frequency analysis in arid regions: A case study for Australia. Journal of Hydrology, v. 475, p. 74-83, 2012. https://dx.doi.org/10.1016/j.jhydrol.2012.08.054
https://dx.doi.org/10.1016/j.jhydrol.201...
). Concerning the number of fluviometric gauge stations, Brazil has only a few stations, and some of them present short or even insufficient streamflow data, especially in small-medium basins (Beskow et al., 2014BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SIMÕES, M. C.; CALDEIRA, T. L.; NUNES, G. S. Índices de sazonalidade para regionalização hidrológica de vazões de estiagem no Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 18, n. 7, p. 748-754, 2014. https://dx.doi.org/10.1590/S1415-43662014000700012
https://dx.doi.org/10.1590/S1415-4366201...
).

A common practice adopted to overcome this problem is streamflow regionalization (Agarwal et al., 2016AGARWAL, A.; MAHESWARAN, R.; SEHGAL, V.; KHOSA, R.; SIVAKUMAR, B.; BERNHOFER, C. Hydrologic regionalization using wavelet-based multiscale entropy method. Journal of Hydrology, v. 538, p. 22-32, 2016. https://doi.org/10.1016/j.jhydrol.2016.03.023
https://doi.org/10.1016/j.jhydrol.2016.0...
), which is still considered a challenging issue for hydrological science (Samuel et al., 2011SAMUEL, J.; COULIBALY, P.; METCALFE, R. A. Estimation of continuous streamflow in Ontario ungauged basins: comparison of regionalization methods. Journal of Hydrologic Engineering, v. 16, n. 5, p. 447-459, 2011. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000338
https://doi.org/10.1061/(ASCE)HE.1943-55...
; Agarwal et al., 2016AGARWAL, A.; MAHESWARAN, R.; SEHGAL, V.; KHOSA, R.; SIVAKUMAR, B.; BERNHOFER, C. Hydrologic regionalization using wavelet-based multiscale entropy method. Journal of Hydrology, v. 538, p. 22-32, 2016. https://doi.org/10.1016/j.jhydrol.2016.03.023
https://doi.org/10.1016/j.jhydrol.2016.0...
). In general, regionalization is a technique that transfers hydrologic information from a monitored drainage basin to sub-basins with scarce data (Mwale et al., 2011MWALE, D.; GAN, T. Y.; DEVITO, K. J.; SILINS, U.; MENDOZA, C.; PETRONE, R. Regionalization of runoff variability of Alberta, Canada, by wavelet, independent component, empirical orthogonal function, and geographical information system analyses. Journal of Hydrologic Engineering, v. 16, n. 2, p. 93-107, 2011. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000284
https://doi.org/10.1061/(ASCE)HE.1943-55...
). The quality of the regionalization depends on the hydrological homogeneity between the sites of interest and the monitored site (Razavi and Coulibali, 2013RAZAVI, T.; COULIBALY, P. Streamflow prediction in ungauged basins: review of regionalization methods. Journal of Hydrologic Engineering, v. 18, n. 8, p. 958-975, 2013. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000690
https://doi.org/10.1061/(ASCE)HE.1943-55...
). Beskow et al. (2014)BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SIMÕES, M. C.; CALDEIRA, T. L.; NUNES, G. S. Índices de sazonalidade para regionalização hidrológica de vazões de estiagem no Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 18, n. 7, p. 748-754, 2014. https://dx.doi.org/10.1590/S1415-43662014000700012
https://dx.doi.org/10.1590/S1415-4366201...
emphasized that clustering into homogeneous regions increases the success potential of the regionalization study.

In the quantile regionalization method, streamflows are calculated for different return periods (RP) based on the frequency distribution for J stations (QTJ) within the region. From these results, the parameters of the model are estimated by multiple regression procedures taking QTJ values and exploratory variables, such as drainage area, length of the main watercourse, slope of the basin, and rainfall of the respective stations (Naghettini and Pinto, 2007NAGHETTINI, M.; PINTO, É. J. de A. Hidrologia estatística. Brasília: CPRM, 2007.). Multiple regression is one of the most used methods to obtain a model for the regionalization of the streamflows (Pruski et al., 2012PRUSKI, F. F.; NUNES, A. A.; REGO, F. S.; SOUZA, M. F. Extrapolação de equações de regionalização de vazões mínimas: Alternativas para atenuar os riscos. Water Resources and Irrigation Management, v. 1, n. 1, p. 51-59, 2012.).

The regionalization method based on the probability density functions (PDF) parameters requires the same frequency distribution for all fluviometric stations. Further, the parameters are related to the physiographic and/or climatic characteristics of the respective basin. The PDF parameters can be derived based on statistical criteria for their fitting, and then, they can be obtained with a single function for the entire region (Cassalho et al., 2018CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
).

Other approaches to regionalization have been used. McIntyre et al. (2005)MCINTYRE, N.; LEE, H.; WHEATER, H.; YOUNG, A.; WAGENER, T. Ensemble predictions of runoff in ungauged catchments. Water Resources Research, v. 41, n. 12, 2005. https://doi.org/10.1029/2005WR004289
https://doi.org/10.1029/2005WR004289...
presented a method based on ensemble modeling in which the relationship between the model parameters and the basin characteristics (e.g., area, perimeter) should be dealt with as a response surface of likelihood. This surface must be integrated into a representative range of the basin characteristics to generate a continuous joint distribution function. Therefore, the model parameters are kept as an ensemble and its interdependencies are not neglected. Despite some advantages, McIntyre et al. (2005)MCINTYRE, N.; LEE, H.; WHEATER, H.; YOUNG, A.; WAGENER, T. Ensemble predictions of runoff in ungauged catchments. Water Resources Research, v. 41, n. 12, 2005. https://doi.org/10.1029/2005WR004289
https://doi.org/10.1029/2005WR004289...
pointed to an underestimation of peak flows in a study carried out in the United Kingdom.

A “similarity” approach is when the parameters from a gauged basin are calibrated and transferred to the ungauged basin. In this procedure, an appropriate calibrated basin is selected through physical similarities between it and the ungauged basin. According to Choubin et al. (2019)CHOUBIN, B.; SOLAIMANI, K.; REZANEZHAD, F.; ROSHAN, M. H.; MALEKIAN, A.; SHAMSHIRBAND, S. Streamflow regionalization using a similarity approach in ungauged basins: Application of the geo-environmental signatures in the Karkheh River Basin, Iran. Catena, v. 182, p. 104128, 2019. https://doi.org/10.1016/j.catena.2019.104128
https://doi.org/10.1016/j.catena.2019.10...
, accessing some geo-environmental information may be the main limitation of the similarity approach, mainly if rainfall-runoff models are used to estimate the streamflow in ungauged basins. Gibbs et al. (2012)GIBBS, M. S.; MAIER, H. R.; DANDY, G. C. A generic framework for regression regionalization in ungauged catchments. Environmental Modelling & Software, v. 27, p. 1-14, 2012. https://doi.org/10.1016/j.envsoft.2011.10.006
https://doi.org/10.1016/j.envsoft.2011.1...
point out that for regions where the number of appropriated calibrated basins is low, the similarity approach is less attractive.

Given that the regionalization of quantiles obtained from regression using as independent variables the physiographic characteristics has presented satisfactory results in Brazil (Maciel et al., 2019MACIEL, A. L.; VIEIRA, E. M.; MONTE MOR, R. C. Regionalização e espacialização de vazões de permanência: Estudo aplicado na bacia rio Piracicaba - MG. Revista brasileira de climatologia, v. 24, p. 114-133, 2019. http://dx.doi.org/10.5380/abclima.v24i0.58420
http://dx.doi.org/10.5380/abclima.v24i0....
; Silva et al., 2006SILVA, A. M. D.; OLIVEIRA, P. M. D.; MELLO, C. R. D.; PIERANGELI, C. Vazões mínimas e de referência para outorga na região do Alto Rio Grande, Minas Gerais. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 10, n. 2, p. 374-380, 2006. http://dx.doi.org/10.1590/S1415-43662006000200019
http://dx.doi.org/10.1590/S1415-43662006...
), this method has been applied frequently in studies in several places in this country. In addition, in cases that basin characteristics are available in GIS databases, the regression approach may be more appropriate (Gibbs et al., 2012GIBBS, M. S.; MAIER, H. R.; DANDY, G. C. A generic framework for regression regionalization in ungauged catchments. Environmental Modelling & Software, v. 27, p. 1-14, 2012. https://doi.org/10.1016/j.envsoft.2011.10.006
https://doi.org/10.1016/j.envsoft.2011.1...
).

Grande River Basin (GRB) is located in southeastern Brazil. This basin corresponds to approximately 12% of the hydropower produced in Brazil (Nobrega et al., 2011NÓBREGA, M. T.; COLLISCHONN, W.; TUCCI, C. E. M.; PAZ, A. R. Uncertainty in climate change impacts on water resources in the Rio Grande Basin, Brazil. Hydrology and Earth System Sciences, v. 15, n. 2, p. 585, 2011. https://doi.org/10.5194/hess-15-585-2011
https://doi.org/10.5194/hess-15-585-2011...
). Mortes River Basin (MRB) is one of the sub-basins of the GRB, and plays an important hydrological function for the water resources management, especially feeding the Funil Hydropower Plant (Oliveira, 2013OLIVEIRA, V. A. Regionalização de vazões nas regiões das Unidades de Planejamento E Gestão De Recursos Hídricos GD1 e GD2, Minas Gerais. 2013. 99 f. Dissertação (Mestrado em Recursos Hídricos em Sistemas Agrícolas) - Universidade Federal de Lavras, Lavras, 2013.). In addition, MRB develops several mining activities, focused on metallurgy and limestone, and agriculture, mainly coffee production (Minas Gerais, 2013). Despite the importance of MRB and a lot of water rights being discussed in the basin, there are few studies in this sense in the region.

Silva et al. (2006)SILVA, A. M. D.; OLIVEIRA, P. M. D.; MELLO, C. R. D.; PIERANGELI, C. Vazões mínimas e de referência para outorga na região do Alto Rio Grande, Minas Gerais. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 10, n. 2, p. 374-380, 2006. http://dx.doi.org/10.1590/S1415-43662006000200019
http://dx.doi.org/10.1590/S1415-43662006...
performed hydrological regionalization in the GRB. However, PDFs more appropriate to simulate the frequency distribution, with 4 and 5 parameters, were not tested, which limited their study. A more appropriate study is then necessary. Moreover, PDFs with a greater number of parameters can be considered as a fundamental flood distribution and have been used successfully in hydrological modeling (Busaba Bodin et al., 2016BUSABA BODHIN, P.; AM SEO, Y.; PARK, J. S.; KUMPHON, B. O. LH-moment estimation of Wakeby distribution with hydrological applications. Stochastic Environmental Research and Risk Assessment, v. 30, n. 6, p. 1757-1767, 2016. https://doi.org/10.1007/s00477-015-1168-4
https://doi.org/10.1007/s00477-015-1168-...
).

Therefore, a robust statistical study was proposed for the streamflows and, thus, for the regionalization process. The objectives of the study were: i) to compare the performance of PDFs with 2 to 5 parameters using the L-moments method (LMM); and, ii) to develop the streamflow regionalization of the maximum and minimum streamflows (Q7,10 and Q90) for the Mortes River Basin, southeast Brazil.

2. MATERIALS AND METHODS

2.1. Study Area

The study was carried out in the Mortes River Basin (MRB), affluent of the Grande River Basin, in southeastern Brazil. Regarding the Minas Gerais Water Resources Planning, this basin is known as UPGRH-GD2 (Figure 1). MRB is located in southern Minas Gerais State, with an area of approximately 10,540 km². This basin includes the Jacaré, Cervo, and Mortes Rivers, which are the main tributaries to the Funil hydropower plant reservoir (Eduardo et al., 2016EDUARDO, E. N.; MELLO, C. R.; VIOLA, M. R.; OWENS, P. R.; CURI, N. Hydrological simulation as subside for management of surface water resources at the Mortes River Basin. Ciência e Agrotecnologia, v. 40, n. 4, p. 390-404, 2016. http://dx.doi.org/10.1590/1413-70542016404009516
http://dx.doi.org/10.1590/1413-705420164...
), whose installed power is 180 MW.

The predominant climates are the Cwa and Cwb, according to the Köppen classification, with a mean annual precipitation of 1500 mm and an average annual temperature of 18°C (Eduardo et al., 2016EDUARDO, E. N.; MELLO, C. R.; VIOLA, M. R.; OWENS, P. R.; CURI, N. Hydrological simulation as subside for management of surface water resources at the Mortes River Basin. Ciência e Agrotecnologia, v. 40, n. 4, p. 390-404, 2016. http://dx.doi.org/10.1590/1413-70542016404009516
http://dx.doi.org/10.1590/1413-705420164...
). For both climate types, summer is characterized by moderate temperatures and approximately 80% of the total precipitation falls in this season, while the winter is cold and dry.

Figure 1.
Study area location and fluviometric gauges stations.

The terrain is gently undulating to undulating relief, with a predominance of Cambisols in the headwaters and Latosols and Argisols in regions where the terrain is smooth. The altitude ranges from 770 to 1498 m. The predominant vegetation comprises high-altitude fields in the upstream regions with the occurrence of the seasonal semideciduous forest (IMGA, 2013INSTITUTO MINEIRO DE GESTÃO DAS ÁGUAS. Resumo Executivo do Plano Diretor de Recursos Hídricos da Bacia do Rio das Mortes. Belo Horizonte, 2013., Scolforo and Carvalho, 2006SCOLFORO, J. R. S.; CARVALHO, L. M. T. Mapeamento da flora nativa e dos reflorestamentos de Minas Gerais. Lavras: UFLA, 2006. v. 1. 288 p.).

2.2. Streamflow Database

The fluviometric gauge stations available in the Brazilian National Water Agency (ANA) Hidroweb platform in MRB are shown in Figure 1. Because of the availability of monitored data, the streamflow stations were selected based on a minimum of 10 years of daily observations (Beskow et al. 2014BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SIMÕES, M. C.; CALDEIRA, T. L.; NUNES, G. S. Índices de sazonalidade para regionalização hidrológica de vazões de estiagem no Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 18, n. 7, p. 748-754, 2014. https://dx.doi.org/10.1590/S1415-43662014000700012
https://dx.doi.org/10.1590/S1415-4366201...
; Caldeira et al. 2015CALDEIRA, T. L.; BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SOUZA, M. R.; GUEDES, H. A. S. Modelagem probabilística de eventos de precipitação extrema no estado do Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 19, n. 3, p. 197-203, 2015. http://dx.doi.org/10.1590/1807-1929/agriambi.v19n3p197-203
http://dx.doi.org/10.1590/1807-1929/agri...
; Cassalho et al., 2018CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
). The maximum and minimum annual streamflow historical series were defined taking the maximum and minimum streamflow rates observed in each hydrological year (October to September) of the fluviometric gauge stations.

The historical series must be independent, randomness, homogeneity, and stationarity to fit the PDFs (Cassalho et al., 2018CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
). Trend analysis was performed using the Mann-Kendall test (Mann, 1945MANN, H. B. Non-parametric tests against trend. Econometrica, v. 13, p. 245-259, 1945. ; Kendall, 1975KENDALL, M. Rank Correlation Methods. London: Charles Griffin, 1975. ). The primary advantages of this test are that it does not require a Gauss distribution and is less influenced by abrupt changes (Salviano et al., 2016SALVIANO, M. F.; GROPPO, J. D.; PELLEGRINO, G. Q. Análise de tendências em dados de precipitação e temperatura no Brasil. Revista Brasileira de Meteorologia, v. 31, n. 1, p. 64-73, 2016. http://dx.doi.org/10.1590/0102-778620150003
http://dx.doi.org/10.1590/0102-778620150...
).

The test is based on accepting (or not accepting) the null hypothesis H 0 with a probability significant level α. H0 hypothesis assumes that the data (x1, ..., xn) of the series do not exhibit any trends or significant temporal variations.

The Mann-Kendall’s ZMK index follows a normal distribution and can be obtained from Equations 1, 2 and 3. A decreasing trend is observed when ZMK is negative and positive values result in an increasing trend (Salviano et al., 2016SALVIANO, M. F.; GROPPO, J. D.; PELLEGRINO, G. Q. Análise de tendências em dados de precipitação e temperatura no Brasil. Revista Brasileira de Meteorologia, v. 31, n. 1, p. 64-73, 2016. http://dx.doi.org/10.1590/0102-778620150003
http://dx.doi.org/10.1590/0102-778620150...
). In the cited equations, X is the streamflow value in the series, and the variable S can be compared to a normal distribution when n ≥ 8. The indices i and J represent the position of x in the series.

z M K = S - 1 V a r S ; p a r a S > 0 0 ; p a r a S = 0 S + 1 V a r ( S ) ; p a r a S < 0 (1)

S = i = 1 n - 1 j = i + 1 n s g n ( x j - x i ) (2)

s g n ( x j - x i ) = + 1 ; s e x j > x i 0 ; s e x j = x i - 1 ; s e x j < x i (3)

A significance level of 5% was adopted. Thus, the series shows no trend if the absolute value of ZMK is less than or equal to 1.96. This level of significance was chosen because it is the most commonly used to avoid Type I and Type II errors.

2.3. Probability Distribution Functions and Regionalization

Before fitting the PDFs, an exploratory analysis of the data from each station was performed to identify possible outliers, which were identified using a box plot. The parameters of the PDFs were estimated using the Moments-L (LMM). According to Cassalho et al. (2018)CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
, this method is less sensitive to the presence of outliers and less susceptible to bias.

The PDF has an impact on the estimates of quantiles, which may lead to an overestimation or underestimation (Rahman et al., 2015RAHMAN, A.; ZAMAN, M. A.; HADDAD, K.; EL ADLOUNI, S.; ZHANG, C. Applicability of Wakeby distribution in flood frequency analysis: a case study for eastern Australia. Hydrological processes, v. 29, n. 4, p. 602-614, 2015. https://doi.org/10.1002/hyp.10182
https://doi.org/10.1002/hyp.10182...
). The PDFs used in this study were the ones available in the studies by Cassalho et al. (2018)CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
, Naghettini and Pinto (2007)NAGHETTINI, M.; PINTO, É. J. de A. Hidrologia estatística. Brasília: CPRM, 2007. and Mello and Silva (2013)MELLO, C. R. de; SILVA, A. M. da. Hidrologia: princípios e aplicações em sistemas agrícolas. Lavras: UFLA, 2013., i.e., the 3-parameter log-normal (LN3), 3-parameter Pearson (PE3), the Gumbel extreme value, the Generalized Extreme Values (GEV), Gamma, Weibull, the 4-parameter Kappa, the 5-parameter Wakeby, the generalized logistic (GLO) and the generalized Pareto (GPA).

The goodness-of-fit values of the PDFs were evaluated using the Anderson-Darling (AD) test. This test performs well when used to evaluate asymmetric distributions and is one of the most commonly used in hydrology (Baldassarre et al., 2009BALDASSARRE, G. D.; LAIO, F.; MONTANARI, A. Design flood estimation using model selection criteria. Physics and Chemistry of the Earth, v. 34, n. 10-12, p. 606-611, 2009. http://dx.doi.org/10.1016/j.pce.2008.10.066
http://dx.doi.org/10.1016/j.pce.2008.10....
; Rahman et al., 2015RAHMAN, A.; ZAMAN, M. A.; HADDAD, K.; EL ADLOUNI, S.; ZHANG, C. Applicability of Wakeby distribution in flood frequency analysis: a case study for eastern Australia. Hydrological processes, v. 29, n. 4, p. 602-614, 2015. https://doi.org/10.1002/hyp.10182
https://doi.org/10.1002/hyp.10182...
). Another advantage of this test is the ability to compare two or more PDFs (Naghettini and Pinto, 2007NAGHETTINI, M.; PINTO, É. J. de A. Hidrologia estatística. Brasília: CPRM, 2007.).

The homogeneous regions were identified based on statistical criteria, comparing the observed frequencies made dimensionless by dividing the long-term mean streamflows of the respective series (Naghettini and Pinto, 2007NAGHETTINI, M.; PINTO, É. J. de A. Hidrologia estatística. Brasília: CPRM, 2007.). This criterion is necessary to identify stations with the same hydrological behavior so that the data can be extrapolated to other regions without monitoring.

Morphometric data used as inputs for the regression models was the drainage area of each station as it can be easily obtained from topographic maps or remote sensing. The return periods of 2, 5, 10, 20, 50, 100, 500, and 1000 years for maximum streamflow and Q7,10 and Q90 for the minimum streamflow are used to perform the regionalization models.

3. RESULTS AND DISCUSSION

Based on the Mann-Kendall test, at the significance level of 5%, 14 series of maximum streamflow and 17 series of Q7 could be considered as stationary and adequate for fitting the PDFs and, subsequently, for regionalization (Table 1). As discussed by Ishak et al. (2013)ISHAK, E. H.; RAHMAN, A.; WESTRA, S.; SHARMA, A.; KUCZERA, G. Evaluating the non-stationarity of Australian annual maximum flood. Journal of Hydrology, v. 494, p. 134-145, 2013. https://doi.org/10.1016/j.jhydrol.2013.04.021
https://doi.org/10.1016/j.jhydrol.2013.0...
, this step is essential to avoid mistakes that result from the changes caused by non-stationarity of the series in the estimation of streamflow for different return periods, which may compromise the application of the regionalization technique.

Table 1.
Mann-Kendall test results, best PDF fitted and Q90 for fluviometric gauge stations in MRB.

The data demonstrate that GLO, GEV, Gamma, and LN3 PDFs had good performances for most of the maximum and minimum series (Table 2). Furthermore, Weibull distribution also showed a good performance for several Q7 series that it is in accordance with the studies of Lopes et al. (2017)LOPES, T R.; ZOLIN, C. A.; PRADO, G.; PAULINO, J.; ALMEIDA, F. T. Regionalization of maximum and minimum flow in the Teles Pires basin, Brazil. Engenharia Agrícola, v. 37, n. 1, p. 54-63, 2017. http://dx.doi.org/10.1590/1809-4430-eng.agric.v37n1p54-63/2017
http://dx.doi.org/10.1590/1809-4430-eng....
and Barbosa et al. (2005)BARBOSA, S. E. S.; BARBOSA JÚNIOR, A. R.; SILVA, G. Q.; CAMPOS, E. N. B.; RODRIGUES, V. C. Geração de modelos de regionalização de vazões máximas, médias de longo período e mínimas de sete dias para a Bacia do Rio do Carmo, Minas Gerais. Engenharia Sanitária e Ambiental, v.10, n. 1, p. 64-71, 2005. https://doi.org/10.1590/S1413-41522005000100008
https://doi.org/10.1590/S1413-4152200500...
. This demonstrates that such PDFs are the most recommended for MRB, being able to estimate the studied streamflow in this basin.

Table 2.
Numbers of historical series fitted and the best fit for each PDFs according to the AD results.

GEV PDF had a good fit for almost all the series, except for two minimum series (Table 2). Cassalho et al. (2018)CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
, in a regionalization study for the São Gonçalo Mirim Basin, also identified GEV as the best PDF to estimate maximum streamflow, demonstrating a good performance of this distribution.

For maximum streamflow estimates, some studies indicate that the best fits were obtained with the 3-parameter Pearson PDF. However, the LN3P (AD = 0.406) showed slightly inferior performance than GEV (AD = 0.405) (Table 1), and both showed better performance than Pearson PDF. A similar result was found by Lopes et al. (2017)LOPES, T R.; ZOLIN, C. A.; PRADO, G.; PAULINO, J.; ALMEIDA, F. T. Regionalization of maximum and minimum flow in the Teles Pires basin, Brazil. Engenharia Agrícola, v. 37, n. 1, p. 54-63, 2017. http://dx.doi.org/10.1590/1809-4430-eng.agric.v37n1p54-63/2017
http://dx.doi.org/10.1590/1809-4430-eng....
in a regionalization study conducted in the Teles Pires River Basin, in which LN3P PDF showed better adjustment for maximum flows.

Regarding PDFs with more parameters, the 4-parameter Kappa presented adequate results for most of the historical series (57% for the maximum and 76% for minimum streamflows). In addition, Kappa PDF produced the best fit for three series of both maximum and minimum streamflows. This higher performance was also observed by Beskow et al. (2014)BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SIMÕES, M. C.; CALDEIRA, T. L.; NUNES, G. S. Índices de sazonalidade para regionalização hidrológica de vazões de estiagem no Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 18, n. 7, p. 748-754, 2014. https://dx.doi.org/10.1590/S1415-43662014000700012
https://dx.doi.org/10.1590/S1415-4366201...
in a study of rainfall series in the state of Rio Grande do Sul, Brazil, applying Kappa PDF. Despite the good performance, Kappa is still little-used in Brazil. The Cassalho et al. (2018)CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
study was one of the first studies that applied Kappa PDF in a maximum streamflow frequency analysis in Brazil.

The best performing PDF was the 5-parameter Wakeby, which had a good fit for 64% of the maximum series and 88% of the Q7 series. This PDF had the least overall values for the AD test, indicating better estimations of the streamflow in MRB. The results obtained by Yerdelen et al. (2010)YERDELEN, C.; KARIMI, Y.; KAHYA, E. Frequency analysis of mean monthly streamflow in Çoruh Basin, Turkey. Fresenius Environmental Bulletin, v. 19, n. 7, p. 1300-1311, 2010. on streamflow regionalization in the Çoruh Basin, Turkey, showed that Wakeby PDF appeared to be a robust distribution in two of the five regions where the Çoruh Basin was subdivided, which were more heterogeneous. This flexibility of the Wakeby PDF application was also highlighted by Rahman et al. (2015)RAHMAN, A.; ZAMAN, M. A.; HADDAD, K.; EL ADLOUNI, S.; ZHANG, C. Applicability of Wakeby distribution in flood frequency analysis: a case study for eastern Australia. Hydrological processes, v. 29, n. 4, p. 602-614, 2015. https://doi.org/10.1002/hyp.10182
https://doi.org/10.1002/hyp.10182...
.

The historical series were grouped based on an analysis of the empirical frequencies made dimensionless dividing them by their respective mean streamflow. Using these frequencies as criteria and observing their behavior, homogeneous regions were obtained for the series of maximum streamflow and Q7 (Figure 2). The other historical series were not used since they did not present similar hydrological behavior. It has been demonstrated that the method based on the empirical frequencies was able to identify the homogeneous region within MRB.

In a regionalization study conducted in the Doce River Basin, Ribeiro et al. (2005)RIBEIRO, C. B. de M.; MARQUES, F. de A.; SILVA, D. D. da. Estimativa e regionalização de vazões mínimas de referência para a bacia do rio Doce. Engenharia na agricultura, v. 13, n. 2, p. 103-117, 2005. established the minimum threshold of six stations per region. This criterion has also been adopted in other studies in Brazil (Cassalho et al., 2018CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
https://doi.org/10.1007/s11269-017-1810-...
; Lopes et al., 2017LOPES, T R.; ZOLIN, C. A.; PRADO, G.; PAULINO, J.; ALMEIDA, F. T. Regionalization of maximum and minimum flow in the Teles Pires basin, Brazil. Engenharia Agrícola, v. 37, n. 1, p. 54-63, 2017. http://dx.doi.org/10.1590/1809-4430-eng.agric.v37n1p54-63/2017
http://dx.doi.org/10.1590/1809-4430-eng....
). Therefore, the number of stations selected for the homogeneous region (eight series for maximum streamflow and seven series for Q7) is following the abovementioned criteria.

Figure 2.
A) homogeneous regions for the series of Q7 and B) homogeneous regions for the series of maximum streamflow.

Since the Weibull PDF was the only one that had a good fit for all Q7 series within the homogeneous region, this PDF was used to calculate Q7,10, which consists of calculating the minimum streamflow of 7-day consecutive duration (Q7) series with a return period of 10 years. Q7,10 is extremely important for water resources rights in MRB, which is entirely within Minas Gerais State. The minimum streamflow in 90% of the time (Q90%) was obtained from the historical series positioned in descending order, calculating the excess frequency of 90%.

The models fitted for the maximum streamflow associated with the return periods analyzed in this study and the minimum reference flows in the homogeneous region are shown in Table 3. The contribution areas for each station ranged from slightly larger than 20 km² to larger than 14,000 km² (Table 1). This range is relevant for the application of the regionalization models because, according to Pruski et al. (2012)PRUSKI, F. F.; NUNES, A. A.; REGO, F. S.; SOUZA, M. F. Extrapolação de equações de regionalização de vazões mínimas: Alternativas para atenuar os riscos. Water Resources and Irrigation Management, v. 1, n. 1, p. 51-59, 2012., the estimates of streamflow should not be performed for areas outside the intervals used to fit the regression model.

For the maximum streamflow, the PDF parameters could not be regionalized with the best PDF (GEV) due to the coefficient of determination (R2) for the shape parameter (ҡ) was not satisfactory. Therefore, the quantiles were calculated based on the distribution that best fitted to the series in each station (Table 1). In this case, the regionalization was conducted based on the return periods along with the drainage area of the sub-basins.

The only physiographic variable used as input in the models was the drainage area (Table 3). The R2 values indicate the suitability of the models to estimate the maximum streamflow, considering the return periods, and minimum reference streamflow using only the drainage area.

Table 3.
Regionalization models and their respective coefficient of determination.

The result presented in Table 3 is similar to what was observed in the study of Lopes et al. (2017)LOPES, T R.; ZOLIN, C. A.; PRADO, G.; PAULINO, J.; ALMEIDA, F. T. Regionalization of maximum and minimum flow in the Teles Pires basin, Brazil. Engenharia Agrícola, v. 37, n. 1, p. 54-63, 2017. http://dx.doi.org/10.1590/1809-4430-eng.agric.v37n1p54-63/2017
http://dx.doi.org/10.1590/1809-4430-eng....
and Pruski et al. (2012)PRUSKI, F. F.; NUNES, A. A.; REGO, F. S.; SOUZA, M. F. Extrapolação de equações de regionalização de vazões mínimas: Alternativas para atenuar os riscos. Water Resources and Irrigation Management, v. 1, n. 1, p. 51-59, 2012. in Brazilian basins, where the regressions that considered the drainage area provided better performance than another criterion. However, Silva et al. (2006)SILVA, A. M. D.; OLIVEIRA, P. M. D.; MELLO, C. R. D.; PIERANGELI, C. Vazões mínimas e de referência para outorga na região do Alto Rio Grande, Minas Gerais. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 10, n. 2, p. 374-380, 2006. http://dx.doi.org/10.1590/S1415-43662006000200019
http://dx.doi.org/10.1590/S1415-43662006...
observed that, although the drainage area of the sub-basins has been the only significant morphometric variable for prediction of Q7,10, the inclusion of the annual precipitation improved the models for streamflow estimates in the upper Grande River Basin.

It was observed that as the return period increases, there is a reduction in the goodness-of-fit of the models, as indicated by the R² values (Table 3). This feature is associated with the length of the series, which were mostly short; It is therefore associated with uncertainty regarding the PDF fitting. Nevertheless, the use of equations may assist in the design of hydraulic structures that use the return periods considered here as safety margins.

The models fitted for Q90 and Q7,10 had very high R² values, which indicates reliability for their use for minimum reference streamflow estimates in MRB, mainly for studies involving water rights use concessions in regions with no fluviometric gauge stations.

Despite the good results obtained in the streamflow estimates through the regionalization process, the monitoring of this variable in the field and the quality of the measurements remains a challenge for the hydrological studies in this region. In addition, the regionalization process presents some uncertainties, as does any other technique used for data simulation.

4. CONCLUSIONS

The analysis of the empirical dimensionless frequencies was able to identify homogeneous regions for maximum and minimum streamflows in MRB. Thus, a hydrological behavior was observed mainly in the main channel of the Mortes River.

Although the distributions with fewer parameters were fitted in most of the maximum streamflows series, the PDF with four and five parameters showed greater plausibility in estimating the quantiles according to AD, except GPA PDF. Regarding minimum streamflows, Wakeby PDF was the one that adjusted to the highest number of series and best performance among the analyzed PDFs.

The area demonstrates a satisfactory input for streamflow regionalization, with coefficients of determination R2 greater than 0.9 for all the studied streamflow. Hydrological homogeneity may be the factor that favored this result.

5. ACKNOWLEDGEMENTS

This work was supported by CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior) and UFLA (Universidade Federal de Lavras).

6. REFERENCES

  • AGARWAL, A.; MAHESWARAN, R.; SEHGAL, V.; KHOSA, R.; SIVAKUMAR, B.; BERNHOFER, C. Hydrologic regionalization using wavelet-based multiscale entropy method. Journal of Hydrology, v. 538, p. 22-32, 2016. https://doi.org/10.1016/j.jhydrol.2016.03.023
    » https://doi.org/10.1016/j.jhydrol.2016.03.023
  • BALDASSARRE, G. D.; LAIO, F.; MONTANARI, A. Design flood estimation using model selection criteria. Physics and Chemistry of the Earth, v. 34, n. 10-12, p. 606-611, 2009. http://dx.doi.org/10.1016/j.pce.2008.10.066
    » http://dx.doi.org/10.1016/j.pce.2008.10.066
  • BARBOSA, S. E. S.; BARBOSA JÚNIOR, A. R.; SILVA, G. Q.; CAMPOS, E. N. B.; RODRIGUES, V. C. Geração de modelos de regionalização de vazões máximas, médias de longo período e mínimas de sete dias para a Bacia do Rio do Carmo, Minas Gerais. Engenharia Sanitária e Ambiental, v.10, n. 1, p. 64-71, 2005. https://doi.org/10.1590/S1413-41522005000100008
    » https://doi.org/10.1590/S1413-41522005000100008
  • BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SIMÕES, M. C.; CALDEIRA, T. L.; NUNES, G. S. Índices de sazonalidade para regionalização hidrológica de vazões de estiagem no Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 18, n. 7, p. 748-754, 2014. https://dx.doi.org/10.1590/S1415-43662014000700012
    » https://dx.doi.org/10.1590/S1415-43662014000700012
  • BUSABA BODHIN, P.; AM SEO, Y.; PARK, J. S.; KUMPHON, B. O. LH-moment estimation of Wakeby distribution with hydrological applications. Stochastic Environmental Research and Risk Assessment, v. 30, n. 6, p. 1757-1767, 2016. https://doi.org/10.1007/s00477-015-1168-4
    » https://doi.org/10.1007/s00477-015-1168-4
  • CALDEIRA, T. L.; BESKOW, S.; MELLO, C. R.; FARIA, L. C.; SOUZA, M. R.; GUEDES, H. A. S. Modelagem probabilística de eventos de precipitação extrema no estado do Rio Grande do Sul. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 19, n. 3, p. 197-203, 2015. http://dx.doi.org/10.1590/1807-1929/agriambi.v19n3p197-203
    » http://dx.doi.org/10.1590/1807-1929/agriambi.v19n3p197-203
  • CASSALHO, F.; BESKOW, S.; MELLO, C. R.; MOURA, M. M.; KERSTNER, L.; ÁVILA, L. F. At-site flood frequency analysis coupled with multiparameter probability distributions. Water resources management, v. 32, n. 1, p. 285-300, 2018. https://doi.org/10.1007/s11269-017-1810-7
    » https://doi.org/10.1007/s11269-017-1810-7
  • CHOUBIN, B.; SOLAIMANI, K.; REZANEZHAD, F.; ROSHAN, M. H.; MALEKIAN, A.; SHAMSHIRBAND, S. Streamflow regionalization using a similarity approach in ungauged basins: Application of the geo-environmental signatures in the Karkheh River Basin, Iran. Catena, v. 182, p. 104128, 2019. https://doi.org/10.1016/j.catena.2019.104128
    » https://doi.org/10.1016/j.catena.2019.104128
  • EDUARDO, E. N.; MELLO, C. R.; VIOLA, M. R.; OWENS, P. R.; CURI, N. Hydrological simulation as subside for management of surface water resources at the Mortes River Basin. Ciência e Agrotecnologia, v. 40, n. 4, p. 390-404, 2016. http://dx.doi.org/10.1590/1413-70542016404009516
    » http://dx.doi.org/10.1590/1413-70542016404009516
  • FAN, F. M.; PONTES, P.; PAIVA, R.; COLLISCHONN, W. Avaliação de um método de propagação de cheias em rios com aproximação inercial das equações de Saint-Venant. Revista Brasileira de Recursos Hídricos, v. 19, n. 4, p. 137-147, 2014. https://dx.doi.org/10.21168/rbrh.v19n4.p137-147
    » https://dx.doi.org/10.21168/rbrh.v19n4.p137-147
  • GIBBS, M. S.; MAIER, H. R.; DANDY, G. C. A generic framework for regression regionalization in ungauged catchments. Environmental Modelling & Software, v. 27, p. 1-14, 2012. https://doi.org/10.1016/j.envsoft.2011.10.006
    » https://doi.org/10.1016/j.envsoft.2011.10.006
  • INSTITUTO MINEIRO DE GESTÃO DAS ÁGUAS. Resumo Executivo do Plano Diretor de Recursos Hídricos da Bacia do Rio das Mortes. Belo Horizonte, 2013.
  • ISHAK, E. H.; RAHMAN, A.; WESTRA, S.; SHARMA, A.; KUCZERA, G. Evaluating the non-stationarity of Australian annual maximum flood. Journal of Hydrology, v. 494, p. 134-145, 2013. https://doi.org/10.1016/j.jhydrol.2013.04.021
    » https://doi.org/10.1016/j.jhydrol.2013.04.021
  • KENDALL, M. Rank Correlation Methods. London: Charles Griffin, 1975.
  • LOPES, T R.; ZOLIN, C. A.; PRADO, G.; PAULINO, J.; ALMEIDA, F. T. Regionalization of maximum and minimum flow in the Teles Pires basin, Brazil. Engenharia Agrícola, v. 37, n. 1, p. 54-63, 2017. http://dx.doi.org/10.1590/1809-4430-eng.agric.v37n1p54-63/2017
    » http://dx.doi.org/10.1590/1809-4430-eng.agric.v37n1p54-63/2017
  • MACIEL, A. L.; VIEIRA, E. M.; MONTE MOR, R. C. Regionalização e espacialização de vazões de permanência: Estudo aplicado na bacia rio Piracicaba - MG. Revista brasileira de climatologia, v. 24, p. 114-133, 2019. http://dx.doi.org/10.5380/abclima.v24i0.58420
    » http://dx.doi.org/10.5380/abclima.v24i0.58420
  • MANN, H. B. Non-parametric tests against trend. Econometrica, v. 13, p. 245-259, 1945.
  • MCINTYRE, N.; LEE, H.; WHEATER, H.; YOUNG, A.; WAGENER, T. Ensemble predictions of runoff in ungauged catchments. Water Resources Research, v. 41, n. 12, 2005. https://doi.org/10.1029/2005WR004289
    » https://doi.org/10.1029/2005WR004289
  • MELLO, C. R. de; SILVA, A. M. da. Hidrologia: princípios e aplicações em sistemas agrícolas. Lavras: UFLA, 2013.
  • MWALE, D.; GAN, T. Y.; DEVITO, K. J.; SILINS, U.; MENDOZA, C.; PETRONE, R. Regionalization of runoff variability of Alberta, Canada, by wavelet, independent component, empirical orthogonal function, and geographical information system analyses. Journal of Hydrologic Engineering, v. 16, n. 2, p. 93-107, 2011. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000284
    » https://doi.org/10.1061/(ASCE)HE.1943-5584.0000284
  • NAGHETTINI, M.; PINTO, É. J. de A. Hidrologia estatística. Brasília: CPRM, 2007.
  • NÓBREGA, M. T.; COLLISCHONN, W.; TUCCI, C. E. M.; PAZ, A. R. Uncertainty in climate change impacts on water resources in the Rio Grande Basin, Brazil. Hydrology and Earth System Sciences, v. 15, n. 2, p. 585, 2011. https://doi.org/10.5194/hess-15-585-2011
    » https://doi.org/10.5194/hess-15-585-2011
  • OLIVEIRA, V. A. Regionalização de vazões nas regiões das Unidades de Planejamento E Gestão De Recursos Hídricos GD1 e GD2, Minas Gerais. 2013. 99 f. Dissertação (Mestrado em Recursos Hídricos em Sistemas Agrícolas) - Universidade Federal de Lavras, Lavras, 2013.
  • PRUSKI, F. F.; NUNES, A. A.; REGO, F. S.; SOUZA, M. F. Extrapolação de equações de regionalização de vazões mínimas: Alternativas para atenuar os riscos. Water Resources and Irrigation Management, v. 1, n. 1, p. 51-59, 2012.
  • RAHMAN, A.; ZAMAN, M. A.; HADDAD, K.; EL ADLOUNI, S.; ZHANG, C. Applicability of Wakeby distribution in flood frequency analysis: a case study for eastern Australia. Hydrological processes, v. 29, n. 4, p. 602-614, 2015. https://doi.org/10.1002/hyp.10182
    » https://doi.org/10.1002/hyp.10182
  • RAZAVI, T.; COULIBALY, P. Streamflow prediction in ungauged basins: review of regionalization methods. Journal of Hydrologic Engineering, v. 18, n. 8, p. 958-975, 2013. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000690
    » https://doi.org/10.1061/(ASCE)HE.1943-5584.0000690
  • RIBEIRO, C. B. de M.; MARQUES, F. de A.; SILVA, D. D. da. Estimativa e regionalização de vazões mínimas de referência para a bacia do rio Doce. Engenharia na agricultura, v. 13, n. 2, p. 103-117, 2005.
  • SALVIANO, M. F.; GROPPO, J. D.; PELLEGRINO, G. Q. Análise de tendências em dados de precipitação e temperatura no Brasil. Revista Brasileira de Meteorologia, v. 31, n. 1, p. 64-73, 2016. http://dx.doi.org/10.1590/0102-778620150003
    » http://dx.doi.org/10.1590/0102-778620150003
  • SAMUEL, J.; COULIBALY, P.; METCALFE, R. A. Estimation of continuous streamflow in Ontario ungauged basins: comparison of regionalization methods. Journal of Hydrologic Engineering, v. 16, n. 5, p. 447-459, 2011. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000338
    » https://doi.org/10.1061/(ASCE)HE.1943-5584.0000338
  • SCOLFORO, J. R. S.; CARVALHO, L. M. T. Mapeamento da flora nativa e dos reflorestamentos de Minas Gerais. Lavras: UFLA, 2006. v. 1. 288 p.
  • SILVA, A. M. D.; OLIVEIRA, P. M. D.; MELLO, C. R. D.; PIERANGELI, C. Vazões mínimas e de referência para outorga na região do Alto Rio Grande, Minas Gerais. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 10, n. 2, p. 374-380, 2006. http://dx.doi.org/10.1590/S1415-43662006000200019
    » http://dx.doi.org/10.1590/S1415-43662006000200019
  • YERDELEN, C.; KARIMI, Y.; KAHYA, E. Frequency analysis of mean monthly streamflow in Çoruh Basin, Turkey. Fresenius Environmental Bulletin, v. 19, n. 7, p. 1300-1311, 2010.
  • ZAMAN, M. A.; RAHMAN, A.; HADDAD, K. Regional flood frequency analysis in arid regions: A case study for Australia. Journal of Hydrology, v. 475, p. 74-83, 2012. https://dx.doi.org/10.1016/j.jhydrol.2012.08.054
    » https://dx.doi.org/10.1016/j.jhydrol.2012.08.054

Publication Dates

  • Publication in this collection
    08 June 2020
  • Date of issue
    2020

History

  • Received
    11 Nov 2019
  • Accepted
    12 Mar 2020
Instituto de Pesquisas Ambientais em Bacias Hidrográficas Instituto de Pesquisas Ambientais em Bacias Hidrográficas (IPABHi), Estrada Mun. Dr. José Luis Cembranelli, 5000, Taubaté, SP, Brasil, CEP 12081-010 - Taubaté - SP - Brazil
E-mail: ambi.agua@gmail.com