In the natural system of a watershed, a number of hydrological processes take place in different space-time frequencies and scales, influenced by physical characteristics related to soil, vegetation, geology, relief, and the drainage network, as well as by prevailing meteorological conditions and mostly the characteristics of the rainfall event. The spatio-temporal patterns that characterize the hydrological response of a watershed may be altered due to changes in climate. These alterations may vary in intensity according to climate characteristics, dry or wet lands, and the magnitude of the change. Global climate changes are expected to take place in the next decades as a result of population growth, land use change and climate change. The Intergovernmental Panel on Climate Change-IPCC Fourth (AR4) and Fifth (AR5) Assessment Reports (available from: www.ipcc.ch) highlight the regions on earth that are more vulnerable to climate change, changes such as the increase or decrease of the total annual precipitation and by dryness or wetness that on the long-term affect the basin-scale water budget. Magrin et al. (2014)  identified changes in stream flow and water availability that are projected to continue into the future in South America, affecting already vulnerable regions. Increasing runoffs in the La Plata River basin and decreasing ones in the Central Andes (Chile, Argentina) in the second half of the 20th century have been associated with changes in precipitation. The current strategy for estimating impacts of climate scenarios on river flow regimes is to use the Global Circulation Models (GCM) outputs as input for hydrological models. These studies may focus on hydrographic basin scale    or global scales     . Arnell and Gosling (2013)  have used 21 GCMs and the global hydrological model Mac-PDM0.09 at a spatial resolution of 0.5˚ × 0.5˚. A quarter of the globe had an increase in runoff in more than two-thirds of GCMs. Nakaegawa et al. (2013)  simulated 24 large river basins around the globe using the MRI-AGCM3.1 GCM (20 km spatial resolution) and TRIP river flow model (0.5˚ spatial resolution). The simulations projected annual runoff reduction in seven river basins and an increase in runoff in seventeen river basins.
In addition to runoff, other hydrological processes are influenced by climate change and the impacts have been studied under the various IPCC scenarios. For example, a seasonal impact resulted from IPCC scenarios was reported by  in South Korea,  in China and  in the River Rhine basin. Arnell and Gosling (2013)  have evaluated the changes in December, January, February (DJF) and June, July, August (JJA) runoff across the global domain. The results showed similar changes in the summer and winter seasons across Brazil. In the other parts of South America, the increase in runoff occurs in JJA. Nakaegawa et al. (2013)  verified that the Amazon River runoff would increase in the high-water season and decrease in low-water, characterizing an amplitude of the seasonal cycle. In the same study, the Parana River had different changes along the basin for the same season (upstream increases less than downstream).
Jung et al. (2013)  evaluated seasonal runoff impacts in South Korea using 13 GCMs with three emission scenarios and three hydrological models. The simulations showed the territory has different responses to scenarios, during the dry season, North is wetter and South is drier, whilst during the wet season, there is flow increase throughout the country. Similar results were obtained by  in the Wei River Basin, China. The hydrological model SWAT simulated scenarios A2 and B1 using precipitation from three GCMs and found a decrease in the low-flow (Q90 in the flow duration curve) and an increase in the high-flow (Q10 in the flow duration curve), suggesting more extreme events of droughts and floods. Demirel et al. (2013)  used three GCMs and three IPCC scenarios in the River Rhine basin upstream from the Dutch- German border. The results showed a shift of the low flows from winter to summer.
In Brazil, studies have identified discharge reduction for semiarid land in Northeast Brazil      . Northeast Brazil is included in those arid or semiarid lands particularly exposed to the impacts of climate change on freshwater  . In other regions of Brazil, other authors have found projected runoff decreasing in the Amazon region   , increasing and decreasing runoff (the signal depends on the GCM used) in the Southeast (Grande River)   , and increasing in hydrographic regions in the Southern Brazil region   . Risk of water supply shortages will increase owing to precipitation reduction and evapotranspiration increase in semiarid regions, thus affecting water supply for cities, hydropower generation and agriculture. All these works have assessed the impacts under IPCC AR4 scenarios.
In view of the observed runoff trends in South America and Brazil, and of the projected trends derived from global and regional models, the objective of this study is to assess the impacts of climate change on the water budget and seasonal stream flow in Brazilian hydrographic regions using the new regional climate change projections from the global HadGEM2-ES (Hadley Global Environment Model 2―Earth System) and MIROC5 (Model for Interdisciplinary Research On Climate) IPCC AR5 models simulated by the Eta regional model. The study is based on the downscaling of the climate simulations from two GCMs and two AR5 emission scenarios, RCP 4.5 and 8.5. The projections are assessed in three time slices: 2011-2040, 2041-2070 and 2071-2099. The study area comprises the entire Brazilian territory and parts of neighboring South American countries included in the delimitation of the Brazilian hydrographic regions. The impacts of the projected hydrological scenarios are assessed in terms of precipitation, soil moisture, evapotranspiration and runoff. Comparison of changes among the different basins may help to identify the most vulnerable river basins as well as the most resilient ones. The novelties of this work are the use of higher resolution downscaling simulations, at 20-km resolution, in the assessment of all major Brazilian river basins, and the use of AR5 IPCC scenarios. The higher spatial resolution improves the descriptions of topography and land cover in a river basin. Previous works have used coarser resolutions for their assessment, studied fewer river basins, and adopted AR4 scenarios.
Section 2 shows the methodology used, describing the study area, the climate and hydrological dataset, the hydrological model and the hydrological indicators selected for the analysis. Section 3 presents the verification of the discharges simulated with the models and the spatial patterns of changes in precipitation, runoff, evapotranspiration, soil moisture, water budget, seasonal runoff and extreme discharges. The last section shows the conclusion of the work taking into account the impacts on the hydrology for establishing adaptation measures.
2.1. Study Area
Brazil is the fifth largest country in the world with 8,511,965.0 km2. The Brazilian National Water Agency (ANA) established the National Hydrographic Division comprising 12 regions. The present study bases spatial discretization on the National Hydrographic Division, but with more details and considering the land beyond the Brazilian border, making the area studied a total of 11,535,645.3 km2. Figure 1 shows the seventeen hydrographic regions defined in the study with their respective drainage areas.
The study area has different climate regimes and land covers, such as tropical forests (Amazon), savannas (Cerrado), semi-arid (Northeast region), wetlands (Pantanal) and temperate climate in the South. Therefore, the major Brazilian river basins may present distinct hydrologic responses induced by climate changes. On the average along time, the entire study area has a contribution of 272,460.0 m3/s.
2.2. Climate and Hydrological Datasets
The downscaling of climate simulations and projections was provided by the Eta regional climate model (RCM). The Eta model    is an atmospheric regional model used operationally by the Brazilian Center for Weather Forecasts and Climate Studies for weather forecasts since 1996  and for seasonal forecasts since 2002  . The model was modified and validated to be suitable for climate changes studies 
Figure 1. Hydrographic regions in the study area.
  . The name of the model stands for the Greek letter η that defines the model vertical coordinate. This coordinate was developed by  to reduce the errors of the terrain-following coordinate (the σ coordinate) in the calculations of pressure gradient near sloping topography. As shown by  , the η coordinate is more suitable to simulate the large-scale circulation over the South America due to the presence of the Andes Cordillera. The current version of the Eta regional climate model has received upgrades in the dynamics and physics components  .
The downscaling of two global climate model simulations, HadGEM2-ES and MIROC5, by the Eta model has been evaluated  and future climate projections were assessed  . A 30-year baseline period was taken between 1961 and 1990. The future period was divided into three time-slices of 30 years: 2011-2040, 2041-2070, and 2071-2099, under RCP 4.5 and RCP 8.5 IPCC AR5 scenarios. The spatial resolution of the Eta model was 20 km in the horizontal with 38 layers in the vertical.
The observed rainfall from the Climatic Research Unit database CRU TS 2.1  was used for comparison between measured and simulated precipitation. The time series used from CRU TS 2.1 has a monthly time step and spatial resolution of 0.5 degrees. The CRU precipitation dataset was used to correct the bias error of RCM precipitation predictions prior to the input to the hydrological model, by using cumulative distribution functions (CDFs) according to  .
River discharges were used to verify the hydrological model simulation of the present-day climate. The 95 daily discharge time series were taken from stream gauges of the Brazilian National Water Agency’s (ANA) hydrometeorological network. Eleven of the 17 hydrographic regions correspond to a specific river basin: Solimões, Negro, Madeira, Tapajós, Xingu, Tocantins, Parnaíba, São Francisco, Uruguay, Paraná and Paraguay. For each of them, the discharges used for comparison in this study were obtained based on the most downstream stream gauge with available observed data, considering the ratio between drainage areas of the hydrographic region and the stream gauge. The other six hydrographic regions (Northern Amazon, Maranhão, Northeast Atlantic, East Atlantic, Southeast Atlantic and South Atlantic) comprise many river basins and observed data from a total of 84 stream gauges was used to obtain the average observed discharges for each of them, again taking into account the drainage area of each stream gauge. The discharge data from the stream gauges was used for obtaining mean annual flows, Q95, Q10 and a set of twelve values from January to December corresponding to mean monthly discharges.
2.3. Hydrological Model
A large scale distributed hydrological model called MGB-IPH (an acronym from the Portuguese for Large Basins Model and Institute of Hydraulic Research)  was used in the simulations. This model has been widely applied in distinct South-American watersheds, for different purposes, such as stream flow forecasting  , simulation of land-use changes  , evaluation of climate change impacts  , assessment of hydrological processes  . The MGB-IPH model may be distributed by cells or catchments and runs in daily time steps. The cell version of the model was used in this work, considering a regular square-based grid of 0.2˚ × 0.2˚. In this model, each cell is divided into hydrological response units (HRU), called patches, which are the combination of unique land-use, vegetation, and soil type. Each HRU has a uniform hydrological response to meteorological forcing.
MGB-IPH follows the Xinanjiang model formulation  to calculate the soil water budget. Three linear reservoirs are used to represent independent routing of surface, subsurface and groundwater flow through the cell. Potential evapotranspiration is calculated by the Penman-Monteith equation. The original version of the MGB-IPH simulates flow propagation in the rivers using the Muskingum-Cunge method; and further studies have included hydrodynamic modelling  . In the present study, however, river flow routing was eliminated and the MGB-IPH run as the global hydrological model Mac-PDM0.09, that is, the model does not route water from one cell to other  .
The hydrological model needs information about soil type and land cover to determine the HRU. The soil type was obtained from the Digital Soil Map of the World created by FAO on a scale of 1:5,000,000  and the land cover map was obtained from the Joint Research Centre  . In the present study, seventeen HRUs were considered in the MGB-IPH. The digital terrain model (DTM) was used for calculation of the time of concentration in each cell of the model. The elevation is represented by the SRTM DTM processed by the Hydrosheds  .
Some of the MGB-IPH model parameters are generally calibrated according to observed and calculated hydrograph comparisons, manually or using an optimization algorithm, as described in  . In the present study, the model parameters were not calibrated due to lack of available stream flow data for the entire area. The alternative was to use the values of the parameters obtained from earlier applications of the MGB-IPH in some Brazilian river basins    . As described in the prior section, the climate data simulated by Eta model is used as input for the MGB-IPH simulations, which comprise the baseline (1961-1990) and three time-slices (2011-2040, 2041-2070 and 2071-2099) for HadGEM2-ES (RCP 4.5 and 8.5) and MIROC5 (RCP 8.5). The results of the simulation using precipitation from climate models and baseline period were compared with the discharge estimated with time series from stream gauges as mentioned before.
2.4. Hydrological Indicators
The processes were analyzed in terms of anomalies of long-term annual average and seasonal average. Precipitation, average river discharge, evapotranspiration and soil moisture were calculated for each hydrographic region.
The seasonality was calculated using river discharges from December, January and February (DJF) for austral summer, and June, July and August (JJA) for austral winter. The goal is to verify if the seasonal changes are equivalent or different to annual average changes under the climate change scenarios. In addition, the seasonal amplitude was calculated using the difference between the maximum and minimum average monthly stream flow according to  .
The extreme values of river discharge are assessed using a duration curve for each cell model. The high-flow is represented by the Q10 (flow exceeded 10% of the time) and the low-flow is represented by the Q95 (flow exceeded 95% of the time).
3. Results and Discussion
Initially, simulations of the river discharge for the baseline period are evaluated. The projections are shown in terms of changes in runoff, water budget, seasonality and extreme discharges.
3.1. Discharge Verification
The hydrological model was simulated for the present-day climate (1961-1990) using variables calculated by climate models. River discharge ratio (RDR) and correlation coefficient (R2) were used to evaluate the hydrological model performance. RDR is the ratio between discharges simulated and observed. The correlation coefficient was determined using the mean monthly discharges observed and simulated by the MGB- IPH/Eta models. Table 1 shows the values obtained from the downscaling by the regional Eta model and MGB-IPH. From now on, the downscaling of the Eta model from the HadGEM2-ES will be referred to as Eta-HadGEM2-ES, and the downscaling from the MIROC5 will be referred to as Eta-MIROC5.
Analyzing the simulations of the MGB-IPH driven by Eta-HadGEM2-ES, fourteen of the hydrographic regions presented RDR within ±30% of ideal value RDR (100%). For the regions Xingu, Northern Amazon and Paraguay, however, the river discharges were overestimated and the corresponding RDR fell out of the ideal range. Similar performance was simulated by the MGB-IPH/Eta-MIROC5 model. The correlation coefficient calculated between mean monthly observed and simulated discharges roughly indicate the agreement between them in terms of reproducing the overall flow regime and seasonality in each hydrographic region. The results obtained were satisfactory, with 12 of the 17 regions with R2 greater than 0.7. Four regions presented R2 lower than 0.5, due to the difficulty in representing the seasonality of the discharge in the basins. The
Table 1. Statistics of observed and calculated river discharges.
a. annual average for baseline period.
extreme situation is verified in the Paraguay river basin where the flood and drought flows were completed inverted (negative R2) due to the strong water storage effect of the Pantanal wetland  . As mentioned before, the version of the MGB-IPH used in the present work does not represent the river flow routing, which prevent the simulation of the storage in the Pantanal.
The verification of the extreme stream flows Q95 and Q10 exhibited inferior quality in comparison to the verification of the mean annual stream flow. The RDR of the Q95 stream flow varied from 21.9% (Tapajós basin and MGB-IPH/Eta-HadGEM2-ES) to 203.7% (Northeast Atlantic and MGB-IPH/Eta-MIROC5). The overall result for the entire study area is an underestimate of the Q95 (46,350 m3/s and 41,870 m3/s, respectively, from MGB-IPH/Eta-HadGEM2-ES and MGB-IPH/Eta-MIROC5 against the observed value of 107,041 m3/s). The RDR of the Q10 stream flow varied from 64.0% (South Atlantic and MGB-IPH/Eta-MIROC5) to 203.0% (Paraguay basin and MGB-IPH/Eta- MIROC5). The overall result for the entire study area is an overestimate of the Q10 (541,091 m3/s and 550,439 m3/s, respectively, from MGB-IPH/Eta-HadGEM2-ES and MGB-IPH/Eta-MIROC5 against the observed value of 474,605 m3/s). The maximum RDR values of Q95 and Q10 probably are influenced, respectively, by the low observed discharge in Northeast Atlantic (Q95 = 42.2 m3/s) and the storage effect in Pantanal wetlands.
3.2. Changes of Precipitation and Runoff
The projections were produced after the definition of the hydrological model parameters. The simulations enabled evaluation of the sensitivity of the hydrographic regions in relation to the precipitation changes. The sensitivity of long-term stream flow to changes in long-term rainfall may be evaluated by using the concept of elasticity. The rainfall elasticity of stream flow is defined as the proportional change in mean annual stream flow divided by the proportional change in mean annual rainfall  .
El is the elasticity
Qsce is the MGB-IPH discharge of the scenario
Qbase is the MGB-IPH discharge of the baseline (1961-1990)
Psce is the Eta precipitation of the scenario
Pbase is the Eta precipitation of the baseline (1961-1990)
Mean values for the entire study area are shown in Table 2. The elasticity varied from 0.54 to 9.79, with a global weighted average of 1.84, characterizing median sensitivity. However, the elasticity is lowest in the Amazon and central-western, whilst Northeast Brazil exhibits the highest values (weighted average of 5.53 and variation from 1.51 to 9.79). This result is in agreement with the conclusions obtained by  in Australia, where the author observed that stream flow is more sensitive to rainfall changes in drier watersheds. Chiew (2006)  found values varying between 2.0 and
Table 2. Mean elasticity for all periods, scenarios and models.
3.5 in 70% of the 219 watersheds. IPCC (2014)  also notes that watershed sensitivity is highly influenced by the runoff coefficient (runoff-rainfall ratio). The lower the coefficient, the higher the sensitivity of the watershed to climate stimulus (precipitation, temperature and evaporation).
3.3. Changes of Evapotranspiration and Soil Moisture
The relation between soil moisture and evapotranspiration is indirectly due to the soil suction, which increases with decreased soil moisture content. Plant roots uptake less water when soil suction is greater.
The soil moisture in Figures 2-4 is calculated with the MGB-IPH and represents the proportion of water in relation to the soil capacity. Observing evapotranspiration maps (Figures 2-4), it is possible to verify that the anomaly of this variable behaves differently in the Amazon and in the Northeast. In the Amazon, the anomaly is positive and in the Northeast, it is negative. Since the precipitation reduces in both regions, evapotranspiration is expected to have negative anomaly in both regions. The positive anomaly of evapotranspiration in the Amazon may be related to the greater capacity of the soil to regulate the evapotranspiration. Despite the lower precipitation in the scenarios, the high soil capacity in the Amazon promotes increased evapotranspiration. In the Northeast, soil storage capacity is lower and, consequently, the moisture condition of the soil affects the plant transpiration. Figures 2-4 also show that soil moisture has a direct relation to precipitation in all regions.
Figure 2. Changes in water budget components with respect to baseline simulations of the MGB-IPH/Eta-HadGEM2-ES under RCP 4.5 scenario.
Figure 3. Changes in water budget components with respect to baseline simulations of the MGB-IPH/Eta-HadGEM2-ES under RCP 8.5 scenario.
Figure 4. Changes in water budget components with respect to baseline simulations of the MGB-IPH/Eta-MIROC5 under RCP 8.5 scenario.
3.4. Water Budget Effects
The climate change scenarios projected changes in the annual water budget, with a spatial variability as depicted in Figures 2-4. In general, there is a gradient in the patterns of changes projected for each scenario relative to the three time slices, i.e. the changes vary gradually along time. Mostly, a positive or negative change is intensified along time, but in some cases the change gradually dampens until the signal decreases. Because of this, our results are described focusing on the central time slice (2041-2070), with comments regarding the differences for the other two time slices.
For the Eta-HadGEM2-ES RCP 4.5 scenario and time slice 2041-2070, the Peruvian, Colombian and Bolivian Amazon show a strong decrease in precipitation (predominantly between 200 and 300 mm, but with peaks around 1000 mm in the highest altitude areas) followed by a decrease in soil moisture (mostly between 5% and 20%). Projected ET changes present a mixed signal in this region, varying from a decrease of up to 100 mm in the westernmost part to an increase of up to 200 mm in the central and eastern parts. These changes in precipitation, soil moisture and ET result in a decreased runoff (mostly between 200 to 400 mm, but reaching peaks greater than 1000 mm in the Andes). This pattern of impacts on the water budget is the same for the other two time slices, although weaker in the 2011-2040 time slice and stronger in the final time slice (2071-2099).
The results for the central part of Amazon and the central part of Brazil (west of the 50˚W) follow in general the previous patterns described for the Peruvian, Colombian and Bolivian regions of Amazon, but with weaker changes. There is a predominance of precipitation reduction, between 50 and 200 mm and soil moisture reduction, between 5% and 30%. The projected evapotranspiration tends to increase (around 50 and 200 mm) west of the 50˚W and negative to the east of this location. However, the general pattern for this entire region is of runoff reduction (mostly between 30 and 250 mm), but slightly stronger to the west of the 50˚W. Again, the results for the time slices 2011-2040 and 2071-2099 show weaker and stronger changes, respectively, relative to the 2041-2070 time slice.
Southern Brazil show increase of precipitation (between 50 and 400 mm), soil moisture (less than 10%), evapotranspiration (around 50 mm) and runoff (between 50 and 400 mm) for the 2041-2070 time-slice. This runoff increase occurs for all the time slices. Runoff increase diminishes, however, between the first and the last time slices. Therefore, according to the Eta-HadGEM2-ES RCP 4.5 simulations, there is initially an increase of runoff but this effect tends to reduce over time.
In the Northeast region of Brazil, results show a dual pattern. There is an area that follows the same pattern described for South Brazil, but the anomalies of the hydrological variables are weaker: increase of precipitation, soil moisture and evapotranspiration that also results in a small runoff increase (less than 50 mm). Again, change is dampened along the three time-slices. The larger area around this core part and also a ~100-km narrow strip along the coast follow the patterns of projected changes described for the central part of Brazil (east of the 50˚-meridian) for the following: precipitation, soil moisture, evapotranspiration and runoff reduction.
Considering the scenario results of the Eta-HadGEM2-ES RCP 8.5, in general it is clear the maintenance of the patterns described for the previous scenario, with an increase of the changes along time (from the first to the last time slice). Large portions of Brazil, including central and northern Amazon, central Brazil and part of southeast, present precipitation reduction between 300 and 800 mm and soil moisture reduction between 25% and 70% for the 2071-2099 time slice. The evapotranspiration increase (reduction) to the west (east) of the 50˚-meridian shown in Eta-HadGEM2-ES RCP 4.5 scenario is intensified in the Eta-HadGEM2-ES RCP 8.5 scenario.
The result of the projected water budget for this scenario is a runoff decrease ranging between 20 and 200 mm for most of the studied area in all of the three time slices. Some areas in the central Amazon present even higher reduction, with the runoff reaching up to 350 mm in the 2011-2040 time slice, and this pattern intensifies and expands spatially over time, in such a way that for the furthest time slice almost the entire Amazon basin presents a runoff reduction ranging between 300 and 800 mm, with peaks of up to 2000 mm in highest altitude areas.
In the South Brazil, the Eta-HadGEM2-ES RCP 8.5 scenario leads to precipitation increase greater than in the Eta-HadGEM2-ES RCP 4.5 scenario. However, the runoff increase in the RCP 8.5 scenario is slightly smaller than in the RCP 4.5 scenario, with an additional relevant difference: the runoff increase raises along the three time slices in the 8.5 scenario, while it reduces along time for the first emission scenario.
In general, the Eta-MIROC-ES RCP 8.5 scenario projected changes following the spatial patterns described for both Eta-HadGEM2 scenarios. However, the intensity of MIROC-ES RCP 8.5 may be considered as intermediate impact between HadGEM2-ES RCP 4.5 and RCP 8.5 scenarios.
3.5. Impacts on Seasonal Runoff and Extreme Discharges
One possible impact of the change of the climate is the shift and duration of the seasons. This may mean summer or winter starting earlier or later, being shorter or longer, wetter or drier. Figure 5 and Figure 6 show the anomalies of the mean runoff for austral summer (DJF) and winter (JJA), respectively. The seasonal anomalies follow a similar pattern to the runoff impacts with decreases in the Amazon, the Northeast and Central Brazil and an increase in the South. However, the anomalies are more remarkable in the summer. During the winter, the South Amazon is less impacted than during the summer. In the South region of Brazil, the anomaly in the winter is still positive, the intensity of the anomaly signal is lower and the impacted area is smaller, when compared with summer.
Table 3 exhibits the projected seasonal amplitude for the period 2041-2070 as a percentage of the baseline simulations. Negative values indicate reduction of the seasonal amplitude and positive values mean increase of the amplitude. The values indicate a relationship between mean annual runoff changes and the value of the seasonal amplitude change.
Figure 5. Changes of DJF (austral summer) mean runoff (mm/day) with respect to the baseline.
The extreme discharges may also change in the future climate scenarios. Figure 7 and Figure 8 verify that the pattern of impacts is similar to the seasonality impacts. South Brazil exhibits positive anomaly for high-flow and low-flow discharges in similar areas for the Eta-MIROC5 model. The Eta-HadGEM-ES model results also have positive anomaly in the South, but with restricted areas for low-flow discharge. The worst situation happens when the anomaly is positive for high-flow discharge and negative for low-flow discharge. In other words, there would be an intensification of extreme discharges, leading to a condition of more natural disasters related to floods and droughts in the same area. This characteristic may be observed in the Northwest Amazon and in a small area in the Northeast Atlantic.
Figure 6. Changes of JJA (austral winter) mean runoff (mm/day) with respect to the baseline.
The spatial pattern of the runoff anomalies is quite similar to the results of studies that use global hydrological modelling    , but disagree with the results of studies such as  , which found positive anomalies in the Amazon and parts of Northeast and  , who verified runoff decrease in South Brazil using climate and hydrological multimodel.
Analysis of hydropower-producing river basins shows the expressive impact arising from the stream flow regime change. The Parana River basin projections indicate discharge increase for Eta-MIROC5 and reduction for Eta-HadGEM2-ES. São Francisco and Tocantins River basins exhibit discharge reduction for both climate models and all
Table 3. Seasonal amplitude for the period 2041-2070 as percentage of the baseline simulations: (1) MGB-IPH/Eta-HadGEM2-ES RCP 4.5, (2) MGB-IPH/Eta-HadGEM2-ES RCP 8.5, (3) MGB- IPH/Eta-MIROC5 RCP 8.5.
scenarios. The most severe reduction is projected for the Tocantins River basin, reaching values between 50% and 60%. Those three river basins have about 73.0% of the installed power of all hydroelectric power plants of Brazil.
The indirect relation between the anomalies of precipitation and evapotranspiration has also been verified in Central Europe by  . The study of  in the Xingu River Basin used a numerical model to evaluate the influence of the climate on the water budget. The results suggest a decrease of about 2% in precipitation and +3% increase in evapotranspiration. Similar behavior was found with MGB-IPH in the future-climate scenarios simulations for the Amazon region: negative anomaly for precipitation and positive anomaly for evapotranspiration.
The seasonal amplitude of the discharge may be related to the relative variability in stream flow along the year. The results showed positive seasonal amplitude in hydrographic regions of South Brazil, meaning an increase of relative variability in stream flow. This is coincident with the increase of the average annual stream flow in the same region. Arnell and Gosling (2013)  have concluded the relative variability and annual stream flow change in parallel, which is in agreement with the results obtained in this study.
Figure 7. Changes of maximum discharge (Q10) (mm/day) with respect to the baseline.
Projections of changes in the hydrological regimes of the major Brazilian basins are assessed in this work. The projections agree partially with the conclusions on the runoff anomalies pointed out by earlier studies. For example, river discharge should increase in mid-high latitudes and that the sensitivity of the catchment influences the magnitude of the impact  . The hydrographic regions in South Brazil are located in mid-lati- tudes and exhibit positive anomaly. On the other hand, considering that the smaller the ratio runoff to precipitation, the greater the catchment sensitivity, Northeast Brazil should exhibit the most severe results as compared to the Amazon region.
Figure 8. Changes of minimum discharge (Q95) (mm/day) with respect to the baseline.
The assessment of the average annual runoff anomaly in the hydrographic regions shows projected reduction of water availability in almost the entire country. The exception is the South region of Brazil. The Northwest, Amazon and a small area along the Northeast Atlantic coast exhibited intensification of the extremes discharges, where the anomaly is positive for high-flow (Q10) and negative for low-flow (Q95). The increase of high-flow is worrying for regions such as Southeast Atlantic, South Atlantic, Uruguay and Paraná due to high population density, while the highest decreases of the low-flow were verified in the Northwest Amazon. Considering the extreme flood and drought events that affected those regions during the past 10 years, the projected changes may trigger the occurrence of natural disasters such as droughts and flash floods, which could increase vulnerability of people living in risk areas.
Considering the climate change impacts in the agriculture, possible adaptation measures are the expansion of irrigated land and enhancement of the irrigation methods  . In Brazil, regions with agricultural vocation such as Paraná basin and others with high dependency on irrigation such as São Francisco river basin must receive special attention to pursuit adaptation to climate change with respect to agriculture production.
Despite the uncertainties associated with climate and hydrological modelling, it is very likely that the hydrological characteristics in the river basin scale will change in the future  . Given this, some countries and regions in the world are developing procedures for adaptation and risk management practices that try to incorporate the likely changes in the hydrological regimes  . The results of our study point out a similar hydrological condition in Brazil compared to other countries and regions in the world in terms of hydrological regime change up to the end of the 21st century. Adaptation procedures are related to measures of water availability on one hand, and water consumption on the other. Water availability measures refer to the extension of the water infrastructure. However, these measures have physical, environmental and economic constraints. Demand-side measures refer to actions that emphasize the reduction of water loss in the water supply systems, enhancement of the irrigation systems and extension of the water reuse  .
This paper is a contribution of the Brazilian Network on Global Climate Change Research funded by CNPq Grant Number 550022/2014-7 and FINEP Grant Number 01.13.0353.00. This work was also partially funded by CNPq Grant Number 400792/2012-5. The authors acknowledge the Ministry of Science, Technology and Innovation of Brazil for supporting the work through Global Environmental Facility funding (UNDP BRA/10/G32) within the Third National Communication to the United Nations Frame- work Convention on Climate Change. The English text of this paper was revised by Sidney Pratt, Canadian, MAT (The Johns Hopkins University), RSAdip-TESL (Cambridge University).