The objective of this work is to assess the impacts of IPCC AR5 climate change scenarios on water resources and hydrological processes across the entire Brazilian territory. Hydrological simulations are carried out in total drainage area of about 11,535,645 km 2 and average stream flow of about 272,460 m 3/s. The study area consists of different climates and land covers such as the Amazon Forest, Northeast Semiarid, Brazilian Savannah, Pantanal wetlands and temperate climate in the South. The atmospheric forcing to drive the large-scale hydrological model MGB-IPH is derived from the downscaling of two global climate models, HadGEM2-ES and MIROC5, by the Eta Regional Climate Model, at 20 km resolution. The Eta model provided the downscaling of the baseline (1961-1990) and three time-slices (2011-2040, 2041-2070 and 2071-2099). These projections adopted two emission scenarios, the RCP 4.5 and RCP 8.5. The change in the average and extremes of precipitation, evapotranspiration, rates of river discharge and soil moisture were assessed. The simulations showed the response of the hydrographic regions due to change of precipitation and potential evapotranspiration in the scenarios. Water availability decreases in almost the entire study area (exception for the South) and the major basins for hydroelectric power generation are affected. 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 ( Q 10) and negative for low-flow ( Q 95). The results highlight the most climatic sensitive regions in Brazil in terms of hydrological variables and water resources.
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) [
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 [
Jung et al. (2013) [
In Brazil, studies have identified discharge reduction for semiarid land in Northeast Brazil [
An increase in precipitation and a reduction in evapotranspiration from land use changes have been associated with the trend in stream flow of the Parana-Plata river [
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.
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.
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.
The downscaling of climate simulations and projections was provided by the Eta regional climate model (RCM). The Eta model [
[
The downscaling of two global climate model simulations, HadGEM2-ES and MIROC5, by the Eta model has been evaluated [
Precipitation, air temperature, net radiation, relative humidity, wind speed, and surface pressure are the atmospheric variables that drive the hydrological model. Precipitation input is on a daily basis, while the other variables are considered as mean monthly values for calculation of the evapotranspiration using the Penman-Monteith method. The mean monthly values are taken from the average of each time-slice of thirty years. For more information on the characteristics of the Eta regional model and the HadGEM2-ES and MIROC5 global models, please refer to [
The observed rainfall from the Climatic Research Unit database CRU TS 2.1 [
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.
A large scale distributed hydrological model called MGB-IPH (an acronym from the Portuguese for Large Basins Model and Institute of Hydraulic Research) [
MGB-IPH follows the Xinanjiang model formulation [
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 [
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 [
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).
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.
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.
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
Hydrographic Region | Observed (m3/s)a | MGB-IPH/Eta-HadGEM2-ES | MGB-IPH/Eta-MIROC5 | ||||
---|---|---|---|---|---|---|---|
Baseline (m3/s) | RDR (%) | R2 | Baseline (m3/s) | RDR (%) | R2 | ||
Solimões | 101,916.2 | 83,292.6 | 81.7 | 0.41 | 78,409.7 | 76.9 | 0.46 |
Negro | 43,349.3 | 31,981.8 | 73.8 | 0.98 | 31,618.6 | 72.9 | 0.99 |
Madeira | 28,452.0 | 28,368.5 | 99.7 | 0.78 | 29,589.4 | 104.0 | 0.80 |
Xingu | 9130.8 | 15293.7 | 167.5 | 0.98 | 16,346.3 | 179.0 | 0.98 |
Tapajós | 13167.4 | 13,246.8 | 100.6 | 0.96 | 14,708.7 | 111.7 | 0.96 |
Northern Amazon | 21,092.6 | 29,212.0 | 138.5 | 0.92 | 30,650.5 | 145.3 | 0.94 |
Tocantins | 13,040.5 | 10,471.0 | 80.3 | 0.97 | 13,636.6 | 104.6 | 0.96 |
Maranhão | 3561.5 | 3627.3 | 101.8 | 0.49 | 4400.3 | 123.5 | 0.52 |
Parnaíba | 894.2 | 707.1 | 79.1 | 0.65 | 848.7 | 94.9 | 0.61 |
Northeast Atlantic | 671.8 | 826.7 | 123.0 | 0.49 | 842.5 | 125.4 | 0.45 |
São Francisco | 2892.9 | 2726.1 | 94.2 | 0.99 | 2971.9 | 102.7 | 0.99 |
East Atlantic | 1343.9 | 1571.1 | 116.9 | 0.84 | 1666.1 | 123.9 | 0.78 |
Southeast Atlantic | 2747.4 | 2449.1 | 89.1 | 0.99 | 2728.0 | 99.3 | 0.98 |
South Atlantic | 6194.4 | 4885.4 | 78.9 | 0.91 | 4416.6 | 71.3 | 0.92 |
Uruguay | 8749.6 | 7790.8 | 89.0 | 0.88 | 7481.2 | 85.5 | 0.76 |
Parana | 12,013.2 | 11,913.5 | 99.2 | 0.92 | 13,442.6 | 111.9 | 0.96 |
Paraguay | 3242.3 | 4527.7 | 139.6 | −0.65 | 4614.6 | 142.3 | −0.55 |
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 [
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.
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 [
where
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
Hydrographic region | Elasticity |
---|---|
Solimões | 0.84 |
Negro | 0.54 |
Madeira | 0.83 |
Xingu | 0.91 |
Tapajós | 1.04 |
Northern Amazon | 0.73 |
Tocantins | 1.50 |
Maranhão | 1.51 |
Parnaíba | 9.79 |
Northeast Atlantic | 4.15 |
São Francisco | 4.52 |
East Atlantic | 8.25 |
Southeast Atlantic | 1.10 |
South Atlantic | 1.48 |
Uruguay | 2.31 |
Parana | 1.39 |
Paraguay | 1.58 |
3.5 in 70% of the 219 watersheds. IPCC (2014) [
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.
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.
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.
The extreme discharges may also change in the future climate scenarios.
The spatial pattern of the runoff anomalies is quite similar to the results of studies that use global hydrological modelling [
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
Hydrographic region | (1) | (2) | (3) |
---|---|---|---|
Solimões | −16.4 | 13.1 | 26.5 |
Negro | 10.9 | 4.6 | −10.2 |
Madeira | −22.4 | −31.1 | −15.9 |
Xingu | −18.9 | −34.4 | −17.9 |
Tapajós | −17.3 | −37.6 | −21.7 |
Northern Amazon | −22.6 | −15.0 | −24.1 |
Tocantins | −28.5 | −48.2 | −39.4 |
Maranhão | −66.0 | −86.6 | −59.4 |
Parnaíba | −12.1 | −65.6 | −59.1 |
Northeast Atlantic | −6.4 | 1.7 | −48.6 |
São Francisco | −29.3 | −62.7 | −41.2 |
East Atlantic | −10.3 | −52.3 | −2.2 |
Southeast Atlantic | 29.7 | 29.4 | 35.7 |
South Atlantic | 89.4 | 90.7 | 89.9 |
Uruguay | 70.2 | 27.1 | 33.6 |
Parana | −28.9 | −44.8 | 7.8 |
Paraguay | −25.9 | −37.8 | −6.2 |
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 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) [
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 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 [
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 [
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).
Neto, A.R., da Paz, A.R., Marengo, J.A. and Chou, S.C. (2016) Hydrological Processes and Climate Change in Hydrographic Regions of Brazil. Journal of Water Resource and Protection, 8, 1103- 1127. http://dx.doi.org/10.4236/jwarp.2016.812087