*Article* **Hydrological Response to Drought Occurrences in a Brazilian Savanna Basin**

#### **Rubens Junqueira \*, Marcelo R. Viola, Jhones da S. Amorim and Carlos R. de Mello**

Programa de Pós-Graduação em Recursos Hídricos, Departamento de Recursos Hídricos, Universidade Federal de Lavras, Lavras 3037, Brazil; marcelo.viola@ufla.br (M.R.V.); jhonesamorim@gmail.com (J.d.S.A.); crmello@ufla.br (C.R.d.M.)

**\*** Correspondence: rubensjunqueira@live.com

Received: 29 August 2020; Accepted: 14 October 2020; Published: 16 October 2020

**Abstract:** The Brazilian savanna is one of the world's 25 biodiversity hotspots. However, droughts can decrease water availability in this biome. This study aimed to analyze meteorological and hydrological droughts and their influence on the hydrological behavior in a Brazilian savanna basin. For that, hydrological indicators were calculated to analyze the hydrological behavior in the Pandeiros river basin (PRB). The Standardized Precipitation Index (SPI) and Standardized Streamflow Index (SSI) were calculated for the hydrological year and rainy season from 1977 to 2018. The propagation of the meteorological to hydrological drought was studied by means of the Pearson coefficient of correlation between the SSI and SPI with 0, 3, 6, 9, and 12-month lags. A longer meteorological drought was observed from 2014/15 to 2017/18 which caused a reduction in the groundwater recharge, besides potentially reducing the ecological functions of the Brazilian savanna. This drought was intensified by an increase in the average annual temperature, resulting in the increasing of evapotranspiration. Regarding drought propagation, there is no significant difference among the coefficients of correlation from 0 to 6-month lags. For the lags of 9 and 12 months, the correlation decreases, indicating a greater influence of the current rainy season than the past ones.

**Keywords:** extreme events; hydrological indicators; SPI; water resources management

#### **1. Introduction**

The Brazilian savanna (known as "Cerrado") is the second-largest biome in South America, covering an area of 2 million km2. It is one of the world's 25 biodiversity hotspots due to the endemic species concentration and the high degree of threat [1,2]. Moreover, this biome plays an important role in providing water, maintaining its ecohydrological functionality for industry, agriculture, navigation, tourism, and hydroelectricity in several Brazilian and South American basins, including the São Francisco river basin (SFRB) [3]. Therefore, it is important to conduct hydrological studies for better assisting water resource management in the Brazilian savanna aiming to maintain its eco-hydrological services [4].

Hydrological indicators, such as the aquifer restitution rate (ARR) and the surface runoff rate (SRR), have been widely used to evaluate the hydrological behavior of Brazilian river basins [4–6]. Nevertheless, during extreme events, for example, floods and droughts, the basin may present different behavior. Extreme drought events can abnormally reduce the streamflow, agricultural production, lake and reservoir levels, and groundwater recharge [7].

Recently, several studies on drought have been developed in Brazil focused on meteorological droughts [8–10]. According to Van Loon [7], this drought is associated with a deficiency in precipitation, and possibly an increase in evapotranspiration. However, few studies have attempted to analyze the influence of meteorological drought on the streamflow (hydrological drought).

Although the origin of hydrological drought is commonly related to meteorological drought, other factors can influence it, highlighting lithology, vegetation, and human influence [11]. Junqueira et al. [12] reported a longer hydrological drought than meteorological ones from 2015 to 2017 due to the reduction in groundwater recharge in previous seasons in the Tocantins river basin, Brazil. As a result, the effects on irrigation, hydroelectricity, and urban supply were prolonged.

To analyze the occurrence, duration, and intensity of droughts, several indexes have been developed in recent years. The Standardized Precipitation Index (SPI) [13] has been widely used in Brazil and worldwide [12,14–16]. According to the World Meteorological Organization (WMO) [17], this index is considered standard due to its accuracy and simplicity. Several indexes have emerged from the SPI, such as the Standardized Streamflow Index (SSI) [11], a hydrological drought index with the same characteristics as the SPI. These indexes represent anomalies from a normal situation and allow for comparison in different regions [7].

To calculate a standardized index, a Probability Distribution Function (PDF) is required. McKee et al. [13] applied the two parameters of Gamma distribution for SPI. Nevertheless, to be more accurate in estimating droughts, Vicente-Serrano et al. [11] suggested it the most suitable PDF rather than adopting a single distribution for all situations. An inappropriate PDF may over- or underestimate the magnitude of the drought, as the extreme events are in the tail of the PDF [16]. Therefore, researchers have analyzed several PDFs to calculate standardized drought indexes and found different results worldwide [11,18,19].

In this context, this study aimed to analyze the occurrence, intensity, duration, and propagation of meteorological to hydrological droughts and their influence on the hydrological behavior in a Brazilian savanna basin.

#### **2. Materials and Methods**

#### *2.1. Study Area*

The Pandeiros river basin (PRB), located in the north Minas Gerais State, is inserted in the Brazilian savanna and has an area of 3220 km2. Due to its ecological relevance for the Brazilian savanna and for the SFRB, the Pandeiros River Environmental Protection Area was created through State Law n◦ 11,901 to protect the native fish species, which represent 70% of the reproduction and development fish from the middle São Francisco River [20,21]. Figure 1 shows the PRB location, the streamflow and rain gauge stations, and the Digital Elevation Model (DEM) ALOS (Advanced Land Observing Satellite) PALSAR (Phased Array L-band Synthetic Aperture Radar), with a spatial resolution of 12.5 m.

The elevation ranges from 496 to 847 m, with an average of 677 m. The climate, according to Köppen type-climate classification, is Aw (tropical with wet summers and dry winters) [22]. The average precipitation for the hydrological year (October to September) is 1085 mm, of which 92% occurs during the rainy period (October to March).

The predominant soil in the PRB is the Latosol, which covers 88.3% of its area [23], and the natural pasture, typical of the Brazilian savanna, is the land use predominant, covering 96.3% of the basin's area [24].

**Figure 1.** Locations of the PRB, streamflow and rain gauge stations, and DEM.

#### *2.2. Input Data*

Daily streamflow and precipitation datasets were obtained from the Brazilian National Water Agency (ANA) from October 1977 to September 2018 (41 hydrological years), using only those with up to 10% of gaps, aiming to meet the data quality standards for carrying out drought analysis [19]. The double mass curve was used to analyze the homogeneity and consistency of the precipitation series [25]. Then, gap filling was carried out on a monthly time scale using the Inverse Squared Distance Weighted method. The basin-scale precipitation was obtained using the area-weighted procedure with the eight rain gauge stations available (Figure 1) [26].

#### *2.3. Hydrological Behavior*

#### 2.3.1. Baseflow

To appraise the hydrological behavior of the PRB, the baseflow was taken using the recursive digital filters method, which divides the streamflow into direct surface runoff and baseflow (Equation (1)) [27].

$$\mathbf{y\_k} = \mathbf{f\_k} + \mathbf{b\_k} \tag{1}$$

where y is the total streamflow (m3 s<sup>−</sup>1); f is the direct surface runoff (m3 s<sup>−</sup>1); b is the baseflow (m3 s<sup>−</sup>1); and k is the time step (day).

For bk calculation, the methodology proposed by Eckhardt [27] for the general form of the recursive filter with one parameter is given by Equation (2).

$$\mathbf{b}\_{\mathbf{k}} = \frac{(1 - \text{BFI}\_{\text{max}}) \times \mathbf{a} \times \mathbf{b}\_{\mathbf{k}-1} + (1 - \mathbf{a}) \times \mathbf{BFI}\_{\text{max}} \times \mathbf{y}\_{\mathbf{k}}}{1 - \mathbf{a} \times \mathbf{BFI}\_{\text{max}}} \tag{2}$$

where BFImax is the maximum value of the baseflow index that can be assumed equal to 0.8 as recommended by Eckhardt [27] for perennial streams; and a corresponds to the groundwater recession constant (0.8925), dimensionless. This model assumes that the aquifer outflow is linearly proportional to its storage, which means an exponential recession of baseflow during the period without groundwater recharge [27].

#### 2.3.2. Hydrological Indicators

The hydrological indicators adopted to study the hydrological behavior in the PRB are presented in Table 1.



BFI is the long-term ratio between baseflow and total streamflow and the closer to 1 it is, the greater the contribution of baseflow to the streamflow. The ARR was obtained by the long-term ratio between baseflow and precipitation and indicates the precipitation rate that contributes to aquifer restitution. The evapotranspiration was calculated based on annual water balance (ET = precipitation − total streamflow, in mm), and then the ETR was calculated as ET divided by precipitation. The SRR consists of the long-term ratio between direct surface runoff (SR = total streamflow − baseflow) and precipitation. The α and Af are described in Equations (3) and (4). According to Silva et al. [28], α values close to zero indicate higher natural regularization capacity.

$$\mathbf{Q\_{l}} = \mathbf{Q\_{0}} \times \mathbf{e^{-\alpha t}}\tag{3}$$

$$\mathbf{A}\_{\mathbf{f}} = \frac{\mathbf{Q}\_{\mathbf{i}} \times 86.4}{\alpha \times \mathbf{A}\_{\mathbf{b}}} \tag{4}$$

where Q0 was the streamflow at the beginning of the recession, in m3 s−1; Qt is the streamflow after *t* days of the beginning of the recession, in m<sup>3</sup> s<sup>−</sup>1; Qi is the streamflow at the end of the hydrological year, in m<sup>3</sup> s<sup>−</sup>1; and Ab is the basin's area, in km2.

The Qmean was obtained from the average of the daily streamflow, while the Qmin and Qmax were obtained from the average of the minimum and maximum daily annual streamflow, respectively. To obtain Q7,10, 10 PDFs fitted using the L-moments were analyzed [29]: 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 PDF was then selected based on the best fitting according to the Anderson–Darling (AD) test. The Q90% and Q95% were obtained based on the flow-duration curve. The SY is the result of the calculated streamflows (Q90%, Q95%, etc.) divided by the basin's area, which allows comparing them among different regions or sub-basins.

#### *2.4. Meteorological and Hydrological Droughts*

To analyze the occurrence, intensity, and duration of the droughts in the PRB, the SPI [13] and SSI [11] were calculated by year, considering the hydrological year (October to September), and for half-year time scales, considering the rainy season (October to March).

For the SPI and SSI calculation, the basin-scale precipitation and streamflow series, respectively, were accumulated for the studied period and a PDF was fitted. Following that, they were transformed into a normal distribution, with mean and variance, respectively, equal to zero and one. The analysis of the most suitable PDF was performed in a similar way to the Q7,10. This procedure, according to Vicente-Serrano et al. [11], allows higher precision in obtaining drought indexes. Afterward, the indexes were classified following the WMO [17] classification (Table 2).


**Table 2.** Classification for the SPI and SSI values.

The Pearson correlation of the coefficient (r) with a statistical significance of 5% (α = 0.05) was used to analyze the correlation between the SPI and SSI on different time scales. In addition, to assess the propagation of meteorological to hydrological drought, the correlations between the 0, 3, 6, 9, and 12-month lags with both the SSI and SPI were performed on an annual scale.

#### **3. Results and Discussion**

#### *3.1. Hydrological Behavior*

Figure 2A presents the average monthly precipitation, streamflow, and baseflow for the PRB. Surface runoff occurs predominantly in the rainy season, whereas in the dry season, it is the baseflow. The BFI indicator was equal to 0.80, which is a characteristic of perennial rivers [27], showing the importance of baseflow over the total streamflow in the studied basin. The ARR for the PRB was 14.9%, a high value when considering that 81.4% of the precipitation returns to the atmosphere by evapotranspiration (ETR). Besides, only 3.7% of the total precipitation is converted to surface runoff (SRR). The high ARR and low SRR values are associated with the predominance of Latosol, which is a deep soil with a high infiltration capacity [30], occurring in a flatter topography (average slope of 5.9%). Moreover, the PRB is part of an environmental protection area and its vegetation has been preserved, which favors the soil-water infiltration and reduces direct surface runoff.

The SYmean obtained for the PRB was 6.46 L s−<sup>1</sup> km−2, while the SYmin and SYmax were 3.6 and 24.8 L s−<sup>1</sup> km−2, respectively. These values are low when compared to others obtained for Brazilian savanna basins such as by Rodrigues et al. [4] in the Manual Alves da Natividade river basin. However, the ETR obtained by Rodrigues et al. [4] was 69.1%, lower than that obtained in the PRB (81.4%). Therefore, the high evapotranspiration contributed to the reduction of the streamflow in the PRB.

Regarding Q7,10 (6.9 m3 s−1), for the streamflow that is used as a basis for granting water resources in the region, the PDF that presented the best fitting was Wakeby, with AD equal to 0.30. Amorim et al. [29] found a similar result in the Mortes River, Southeastern Brazil, demonstrating the satisfactory performance of this PDF for the Q7 series. The SY7,10 was equal to 2.1 L s−<sup>1</sup> km−2, similar to that obtained by Silva et al. [6] in the Minas Gerais State. According to the Euclydes et al. [31] classification, the PRB has a high natural regularization capacity once SY7,10 corresponds to 33% of

SYmean, which follows the α value close to zero (0.0073). The α is close to that obtained by Freitas and Bacellar [32] in sub-basins of the upper São Francisco River and by Junqueira et al. [12] in the Tocantins river basin. From the flow-duration curve (Figure 2B), Q90% and Q95% was equal to 9.7 and 7.7 m<sup>3</sup> s<sup>−</sup>1, respectively, which generated SY90 equal to 3.0 L s−<sup>1</sup> km−<sup>2</sup> and SY95 equal to 2.4 L s−<sup>1</sup> km−2, higher than SY7,10.

**Figure 2.** Hydrograph containing average monthly precipitation (P), total streamflow (Qt), and baseflow (Qb) (**A**); and permanent flow curve (**B**).

#### *3.2. Meteorological and Hydrological Droughts*

Among the PDFs analyzed, Weibull presented a better fit in two historical series, whereas Wakeby showed better performance in six. Although McKee et al. [13] have used Gamma for the SPI calculation, this PDF did not produce the best fit in either case. GPA presented unsatisfactory results for the basin-scale precipitation and streamflow series, not being approved in the AD test in both time scales. Table 3 presents the AD test results for the basin-scale precipitation and streamflow series for the hydrological years and rainy seasons obtained from each PDF.


**Table 3.** AD test results for the basin-scale precipitation and streamflow for the PRB.

Note: \* Best fit, Hy = Hydrological year, L3 = 3 months lag, L6 = 6 months lag, L9 = 9 months lag, L12 = 12 months lag, Rs = Rainy season.

The SPI and SSI results for the hydrological years and rainy seasons, calculated according to the most suitable PDF, can be observed in Figure 3. There is an agreement between the SPI values for the hydrological year and the rainy season (r = 0.96). A similar result is observed in the SSI, where r = 0.99. However, when comparing the SPI against the SSI at the two-time scales, the correlation decreases (r = 0.64). This is due to the complexity of the factors associated with the hydrological cycle in the basins, where the response to precipitation depends on factors such as vegetation, lithology, and topography [11].

**Figure 3.** SPI and SSI results for hydrological years and rainy seasons in the PRB.

The most significant meteorological droughts occurred during the hydrological years of 1986/87, 1994/95, 1995/96, 1997/98, 2007/08, and between 2014/15 and 2016/17 in the PRB. Similar behavior was observed for the rainy season, however, with some differences concerning intensity. Although the dry season accounted for only 8% of annual precipitation, the total precipitation in this period for the hydrological year of 1986/87 was 171.1 mm, approximately twice the average (89.5 mm). Therefore, the attenuation of drought on a hydrological year scale may have occurred due to increased precipitation in the dry season.

As a result of some years with below-average precipitation, the hydrological years between 2013/14 to 2017/18 presented a long and intense hydrological drought. In these years, there was a high reduction in the long-term mean streamflow, from 20.8 m3 s−<sup>1</sup> in the whole period to 10.3 m3 s−<sup>1</sup> for the drought period (2013/14 to 2017/18), which is lower than the minimum streamflow for the entire studied period (11.5 m3 s−1). Also, there was a significant reduction in ARR (14.9 to 9.2%), SSR (3.7 to 2.4%), SYmean (6.5 to 3.2 L s−<sup>1</sup> km<sup>−</sup>2), SYmin (3.6 to 1.5 L s−<sup>1</sup> km<sup>−</sup>2), and SYmax (24.8 to 17.4 L s−<sup>1</sup> km<sup>−</sup>2).

Besides the reduction in precipitation, the increase in temperature and, consequently, in evapotranspiration, contributed to the hydrological drought intensification, as highlighted by Van Loon [7] and Junqueira et al. [12]. According to data obtained from the Brazilian National Institute of Meteorology (INMET) using a weather station close to the basin (Januária station), the daily annual average temperature from 1977/78 to 2012/13 was 23.4 ◦C, however, from 2013/14 to 2017/18 there was an increase of 1.3 ◦C. Thus, the ETR in this period was equal to 88.4%, higher than the average observed in the studied period (81.4%).

Studies suggest an increase in the intensity and frequency of droughts by the end of the century worldwide due to atmospheric evaporative demand increases, reducing water availability, soil water

storage, and agricultural production [15,33,34]. Yet, Santos et al. [35] observed a reduction in the water table in the PRB wetlands, downstream of the streamflow gauge station of this study, due to a reduction in precipitation and an increase in temperature in recent years, affecting vegetation dynamics in the wetland.

During the period 2013/14 to 2017/18, similar events were found in other Brazilian regions such as the Doce river basin [14], Tocantins river basin [12], Northeastern Brazil [36], Amazon [36], Paraná river basin [37], the metropolitan region of São Paulo [38], and Ceará State [9], affecting the hydropower generation and urban supply. Besides, droughts can affect the ecological functions of the Brazilian savanna [34], as well as reduce the native fish species and increase exotic species [39].

According to Azevedo et al. [40], after a long period of drought, the Sobradinho hydropower plant reservoir, located in SFRB and responsible for almost 60% of the water resources in Northeastern Brazil, presented a reduction of up to 50% in surface water in 2015/16. In that year, hydropower generation was only 170 MW (its total installed power is 1050 MW).

Some researchers associate drought occurrence with anomalies in sea surface temperatures in the Pacific and the Atlantic Oceans [9,12,36]. Santos et al. [41] reported that anomalies in the El Niño-Southern Oscillation (ENSO) phenomenon are related to the cycle of meteorological droughts in the region in which the PRB is located, however, the Pacific Decadal Oscillation (PDO) influence has not been identified as having a clear influence on drought events in this region. El Niño (positive phase of the ENSO) occurred during all the years in which there was a meteorological drought in the PRB, except for 2007/08 and 2016/17, when the La Niña phenomenon (negative phase of the ENSO) took place for some months. Junqueira et al. [12] and Marengo et al. [36] also found the influence of the ENSO on the drought occurrences in the Tocantins river basin and Northeastern Brazil. According to Garreaud et al. [42], although the PDO-related precipitation and temperature anomalies have the same behavior as the ENSO in South America, their effects seem to be less intense in the PRB region.

#### *3.3. Drought Propagation*

Table 4 presents the correlation results between the SSI and SPI with a lag of 0, 3, 6, 9, and 12 months, both on an annual scale. There is no difference between the r values from a 0 to 6-month lag. However, a higher correlation was obtained for the 3-month lag, indicating a lag in the propagation of the meteorological drought to the hydrological one. A similar result was found by Rodrigues et al. [34] in sub-basins of the Tocantins River, where the correlation between the indices was higher in the period from 0 to 3 months.


**Table 4.** r results between the SSI and the SPI with a lag of 0, 3, 6, 9, and 12 months.

Furthermore, the correlation is lower for 9 and 12-month lags, that is, changes in streamflow are more influenced by the current than the past rainy season. According to Van Loon [7], rainfall deficit in the rainy season can influence following dry season conditions. In cases of longer meteorological drought, there may be a reduction in groundwater recharge, affecting the streamflow even after the meteorological conditions return to normal, as also reported by Jesus et al. [14] in the Doce river basin and by Junqueira et al. [12] in the Tocantins river basin. As an example, in the hydrological year 2017/18, the SPI was classified as near normal, however, the SSI was classified as moderately dry. This increase in precipitation from one year to the next attenuated the hydrological drought but was not enough to fill the recharge deficit from previous years.

According to Junqueira et al. [12], baseflow depends on the precipitation/recharge rate in previous months, as well as on aquifer characteristics. As the streamflow in the PRB is strongly influenced by baseflow (BFI = 0.80), the deficiency in groundwater recharge has a prolonged effect on the streamflow. This deficiency is confirmed by the Af values, where, from 2013/14 to 2017/18, the average was 21.1 mm,

less than half the historical average (48.9 mm). The recharging deficiency increased during the drought until it reached its lowest value at the end of 2016/17 (Af = 16.3 mm). After the 2017/18 hydrological year, there was a slight recovery of water storage in the aquifer (Af = 20.0 mm) due to the increase in precipitation, however, it remained below the historical average.

#### **4. Conclusions**

The SPI allowed us to identify the main droughts that hit the PRB (1986/87, 1994/95, 1995/96, 1997/98, 2000/01, 2007/08, 2014/15, 2015/16, and 2016/17). These events may be related to the occurrence of macro-scale climatic phenomena such as the ENSO, as observed in other Brazilian regions. Although the droughts observed in hydrological years presented a high correlation with the drought during the rainy season, the dry season can occasionally influence the drought intensity.

The hydrological drought from 2013/14 to 2017/18 occurred mainly due to a prolonged reduction in precipitation, reducing water availability, and affecting the fish reproduction in the PRB and SFRB and the ecological functions of the Brazilian savanna. In addition to the reduction in precipitation, the increase in temperature and, consequently, in evapotranspiration contributed to the intensification of the hydrological drought in the cited period.

Based on the hydrological indicators, a high influence of baseflow on the total streamflow (BFI = 0.8) and evapotranspiration (81.4% of the total precipitation) on the hydrological behavior of the PRB was observed. Moreover, the analysis of hydrological indicators highlighted the influence of droughts on groundwater recharge, evidenced by a reduction in Af values.

The propagation of meteorological to hydrological drought occurred in shorter lags (0 to 6 months), with greater correlation in 3-month lag. However, for the lags of 9 and 12 months, the correlation decreased, indicating that changes in streamflow are more influenced by the current than previous rainy seasons.

**Author Contributions:** Conceptualization, R.J., M.R.V. and J.d.S.A.; methodology, R.J., M.R.V. and J.d.S.A.; software, R.J. and J.d.S.A.; validation, R.J.; formal analysis, R.J. and M.R.V.; investigation, R.J.; resources, R.J., M.R.V. and J.d.S.A.; data curation, R.J. and J.d.S.A.; writing—original draft preparation, R.J.; writing—review and editing, R.J., M.R.V., J.d.S.A. and C.R.d.M.; visualization, R.J. and M.R.V.; supervision: M.R.V.; project administration, R.J., M.R.V. and J.d.S.A.; funding acquisition, R.J., M.R.V. and J.d.S.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was funded by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—CAPES (grant numbers 88882.446869/2019-01 and 88882.446854/2019-01) and the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) with the grant for authors, and Universidade Federal de Lavras (UFLA) through public notice PRPG/UFLA n◦ 046/2020 with the support for publication.

**Acknowledgments:** We greatly acknowledge the Brazilian National Water Agency (ANA) and the Brazilian National Institute of Meteorology (INMET) for providing the input data to develop this study.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
