South Asian river basins in a 1.5 °C warmer world

In 2015, with the signing of the “Paris Agreement”, 195 countries committed to limiting the increase in global temperature to less than 2 °C with respect to pre-industrial levels and to aim at limiting the increase to 1.5 °C by 2100. The regional ramifications of those thresholds remain however largely unknown and variability in the magnitude of change and the associated impacts are yet to be quantified. We provide a regional quantitative assessment of the impacts of a 1.5 versus a 2 °C global warming for a major global climate change hotspot: the Indus, Ganges, and Brahmaputra river basins (IGB) in South Asia, by analyzing changes in climate change indicators based on 1.5 and 2 °C global warming scenarios. In the analyzed ensemble of general circulation models, a global temperature increase of 1.5 °C implies a temperature increase of 1.4–2.6 (μ = 2.1) °C for the IGB. For the 2.0 °C scenario, the increase would be 2.0–3.4 (μ = 2.7) °C. We show that climate change impacts are more adverse under 2 °C versus 1.5 °C warming and that changes in the indicators’ values are in general linearly correlated to average temperature increase. We also show that for climate projections following Representative Concentration Pathways 4.5 and 8.5, which may be more realistic, the regional temperature increases and changes in climate change indicators are much stronger than for the 1.5 and 2 °C scenarios.


Introduction
Climate change is identified as one of the major global issues, threatening the planet at various levels (IPCC 2014). Multiple attempts to streamline global policy on climate change mitigation have been made over the past decades, and the BParis Agreement^which was signed at the 21st Conference of the Parties in 2015 is considered a major breakthrough in formulating adequate measures to tackle climate change (Clemencon 2016;Savaresi 2016). During that conference, 195 countries agreed on Ba long-term goal of keeping the increase in global average temperature to well below 2°C above pre-industrial levels,^and Bto aim to limit the increase to 1.5°C, since this would significantly reduce risks and the impacts of climate change.^In response to this development, the Intergovernmental Panel on Climate Change (IPCC) has published a Special Report on global warming of 1.5°C above pre-industrial levels (IPCC 2018). However, scientific evidence of the impacts of a 1.5°C global warming, and the differences in impacts between 1.5°C and a 2°C global warming, is still sparse. Therefore, the scientific community has been mobilized to provide this scientific evidence (NCC 2016). A global study of differences in impacts at 1.5 and 2.0°C forms a starting point for such assessments and highlights the importance of regional differentiation (Schleussner et al. 2016b). The aim of our study is to provide a regional quantitative assessment of the impacts of a 1.5°C versus a 2°C global warming for a major global climate change hotspot: the Indus, Ganges, and Brahmaputra river basins in South Asia (De Souza et al. 2015).
Developing countries and their inhabitants are considered to be more vulnerable to climate change than western countries (Mertz et al. 2009). The glacier-fed Indus, Ganges, and Brahmaputra (IGB) river basins in South Asia are home to a population around 660 million people (Shrestha et al. 2015), which is projected to increase to almost 1 billion by 2050 (UN 2015). With the increasing population, demands for water, food, and energy will increase strongly too. The region is facing enormous challenges in adapting to climate change, which is impacting many different sectors, in addition to conventional development challenges that lead to adaptation deficits (Kilroy 2015;Vinke et al. 2017).
The regional climate is dominated by the Asian summer monsoon, protruding from the south-east to the north-west of the region, and bringing large amounts of precipitation during June to September. The largest amounts of precipitation fall close to the Bay of Bengal, and at the southern slopes of the Hindu-Kush-Himalayas where the steep rise in topography causes orographic precipitation (Galewsky 2009) (Fig. 1). Under climate change, changes in monsoon duration and intensity are expected (Turner and Annamalai 2012;Sperber et al. 2013), and over a historical period increases in dry and wet spells are observed during the monsoon season (Singh et al. 2014). Increases in frequency and magnitude of heavy precipitation events will increase the occurrence of floods (Mirza 2011;Hirabayashi et al. 2013), river bank erosion, and other natural hazards like landslides, debris flows, and glacier-lake outburst floods (GLOFs) (Tariq et al. 2014). Besides, increases in precipitation intensity lead to increasing sedimentation, affecting hydropower production and the lifetime of hydropower infrastructure.
Droughts are expected to increase globally under climate change in combination with an increase in human water use (Trenberth et al. 2013;Wanders and Wada 2014). This affects agricultural production and therefore regional food security (Wheeler and von Braun 2013;Biemans et al. 2013), and may lead to insufficient flows to sustain ecosystem services, or the drying up of springs (Tambe et al. 2012). In the Indus and Ganges basins, large portions of irrigation water come from nonrenewable groundwater already in the present situation (Gleeson et al. 2012;Gleeson and Wada 2013). Heat waves also affect agricultural production, can trigger GLOFs, and affect human health (Pal and Eltahir 2015), and are projected to be severe in South Asia (Im et al. 2017). Changes in temperature and humidity affect the transmissivity of vector-borne diseases like malaria (Dhiman et al. 2011).
The Indus, Ganges, and Brahmaputra rivers depend significantly on glacier and snow melt water (Immerzeel et al. 2010;Lutz et al. 2014), and climate change affects cryospheric water reserves stored in the Hindu-Kush, Karakoram, and Himalayan mountain ranges. Glaciers have been retreating in recent decades (Bolch et al. 2012;Brun et al. 2017), and ice masses will reduce further under climate change (Kraaijenbrink et al. 2017), with temperature increases being stronger at high altitude compared to lower land areas (Pepin et al. 2015;Palazzi et al. 2016). The changes in the cryosphere will lead to changes in downstream water availability, and timing and magnitude of flow (Lutz et al. 2016a).
The Ganges-Brahmaputra delta is additionally threatened by increases in flooding events related to the subsidence of land (Higgins et al. 2014), rising sea levels (Nicholls and Cazenave 2010;Kay et al. 2015), and increases in tropical cyclone activity over the Bay of Bengal during the postmonsoon season (Balaguru et al. 2014).
Many of these impacts can be directly linked to changes in air temperature and precipitation, or climate change indicators derived from those variables. The aim of this study is to provide a quantitative overview of the impacts of climate change for the 1.5°C global warming scenario, the 2.0°C global warming scenario, and an ensemble of downscaled region-representative climate projections for t h e e n d o f t h e t w e n t y -f i r s t c e n t u r y, u n d e r t h e Representative Concentration Pathways (van Vuuren et al. 2011) RCP4.5 and RCP8.5. To this end, we analyze a set of ten climate change indicators, which can be derived from changes in air temperature and precipitation, and compare these for the four scenarios (i.e., 1.5°C, 2.0°C, end of century RCP4.5 and end of century RCP8.5). The ten indicators are chosen in such a way that they can be directly or indirectly related to the most important sectors in the region.

Downscaled climate change scenarios
We use an ensemble of 2 × 4 downscaled general circulation models (GCMs), which is representative for the complete set of climate change projections under RCP4.5 and RCP8.5 in the Coupled Model Intercomparison Project Phase 5 (CMIP5) multi-model ensemble (Taylor et al. 2012), over the IGB region (Table 1). The GCMs have been selected in such a way that they reflect the full ranges of projected climate change in terms of changes in climatic means and climatic extremes. Besides, only climate models with sufficient skill in simulating historical climate in the region have been included in the ensemble (Lutz et al. 2016b). The eight selected GCMs' outputs for precipitation, mean air temperature, maximum air temperature, and minimum air temperature have been downscaled and bias corrected by applying the empiricalstatistical quantile mapping method (e.g., Piani et al. 2010;Themeßl et al. 2011a), to the future GCM projections using a historical reference climate dataset for the IGB basins. This reference dataset covers 1981-2010 at a daily time step with a 10 × 10-km spatial resolution (Lutz and Immerzeel 2015). The downscaled GCMs cover 2011-2100 at the same time step and spatial resolution. The quantile mapping methodology is robust and well established, and has been proven to perform well over mountainous regions (Themeßl et al. 2011b;Immerzeel et al. 2013). It has been successfully applied for the Indus, Ganges, and Brahmaputra river basins before (Wijngaard et al. 2018). We follow the methodology described by Themeßl et al. (2011a). However, instead of calculating empirical cumulative distribution functions (ecdfs) for each day of the year, we construct ecdfs for each month of the year (moy, 12 months). We do this for each 10 × 10km grid cell (i), from the daily values spanning 1981-2010 of the reference climate dataset (obs), which is based on observations.
For the historical GCM runs for 1981-2010, the same procedure is followed to generate ecdfs for each month of the year, after the GCM grids have been smoothed from their original spatial resolution to 10 × 10 km using a bilinear interpolation. Both sets of ecdfs (i.e., for the reference climate dataset (obs) and the historical GCM runs (GCM)) are used to construct a correction function (CF): Fig. 1 a Study area. b Average annual precipitation for reference period (1981. c Average air temperature for reference period (1981 where P t,i is the probability of a future GCM value, which is derived as: The correction function CF is used to downscale and bias-correct the future GCM runs spanning 2011-2100 at a daily time step. For each day and each grid cell, this is done as follows: Similar as described by Themeßl et al. (2011a), we include frequency adaptation, to account for cases where the dry-day frequency in the GCM is greater than in the observations, and the generation of new extremes by linear extrapolation of the correction value (CF) at the highest and lowest quantiles of the calibration range (i.e., the 1981-2010 reference period). See Themeßl et al. (2011a) for further details. In this way, transient gridded meteorological time series from 2011 to 2100 at 10 × 10-km spatial resolution and daily time step are constructed for each of the selected GCM runs.
The quantile mapping is validated by a decadal Bleave one out^cross-validation approach (e.g., Themeßl et al. 2011b, a). This is done by constructing the ecdfs for the reference climate dataset and historical GCM runs for 20 out of 30 years of the reference period (e.g., 1981-2000) and subsequently using those ecdfs and the derived correction function to downscale the historical GCM runs for the remaining 10 years of the reference period (e.g., 2001-2010). This procedure is repeated to downscale the two remaining 10-year periods of the reference period: ecdfs for 1991-2010 are used to downscale 1981-1990 and ecdfs for 1981-1990 and 2001-2010 combined are used to downscale 1991-2000. The downscaled GCM data of precipitation for 1981-2010 is validated to the reference dataset, by comparing the mean precipitation, number of dry days, and 99th percentile value of precipitation. For air temperature, we validate the mean and 1st and 99th percentile values of mean, maximum, and minimum air temperature. This comparison is done for the entire distributions of daily values spanning 30 years covering 1981-2010. Because the reference dataset covers 27,164 grid cells, it is computationally not feasible to do this comparison for each individual grid cell. To overcome this, we do the comparison for 100 randomly sampled grid cells in each of the three river basins.
Extraction and analysis of GCM time slices of 1.5 and 2.0°C global temperature increase with respect to the preindustrial era Because of the large computational effort required for the bias-correction and downscaling described in the BDownscaled climate change scenarios^section, the downscaling has been limited to the IGB region, and the period 2011-2100 which is adjacent to the reference period . Since the 1.5 and 2°C scenarios concern global averages, analyses comparing global and regional future climate to the preindustrial area were done using the raw (i.e., non-downscaled, at their original spatial resolutions) global daily grids of the eight selected GCM runs as available from the CMIP5 archive. We analyze the global temperature change by taking Table 1 Analyzed GCM runs listing time slices for which global temperature changes (ΔT) of 1.5 and 2.0°C compared to pre-industrial (PI) levels are reached, the corresponding regional changes in air temperature (ΔT IGB) and precipitation (ΔP IGB), and the projected regional changes in air temperature (ΔT IGB EOC) and precipitation (ΔP IGB EOC), at the end of the 21 st century (EOC 2071-2100), compared to pre-industrial levels GCM the average of 30-year time slices and comparing them to the average temperature during the pre-industrial period (1861-1890, PI). We start by comparing 1891-1920 to PI and iteratively move the 30-year time slice 1 year forward in time to identify for each GCM the earliest time slices where the global temperature increase thresholds of 1.5 and 2.0°C compared to PI levels are reached (Table 1). For those time slices, we analyze for each GCM run how the regional temperature and precipitation have changed for the IGB basins with respect to the PI period, to identify how the regional temperature changes correspond to global temperature change (Table 1). Furthermore, we analyze for each GCM run how regional temperature levels for the IGB basins will change between the end of the twenty-first century (EOC, 2071-2100) and the PI era. Besides changes in mean air temperature, we analyze regional precipitation changes for the mentioned time slices, for each GCM run ( Table 1).
Analysis of future changes in climate change indicators for 1.5 and 2.0°C global temperature increase scenarios We analyze changes in a set of ten climate change indicators (Table 2) for the 1.5°C time slice, 2.0°C time slice, and end of century projections for RCP4.5 and RCP8.5, compared to the reference period (REF, 1981(REF, -2010. These analyses are done by comparing the downscaled bias-corrected future GCM data to the reference climate dataset. In contrast to the analyses described in the BExtraction and analysis of GCM time slices of 1.5 and 2.0°C global temperature increase with respect to the preindustrial era^section, here, we compare the future to the reference period, and not to the pre-industrial period, to focus on changes that are yet to come, and not on changes that have already happened between the pre-industrial period and the reference period. The eight GCMs' time slices for the 1.5°C scenario are averaged to get an ensemble mean and the same is done for the time slices for the 2.0°C. For the end of century projections for RCP4.5 and RCP8.5, respectively, the four members of each RCP are averaged for 2071-2100 to get their ensemble means. The indicators have been selected in such a way that they can be directly or indirectly linked to the most important climate change impacts in the region, as identified in this study and by others (Kilroy 2015;Vinke et al. 2017) (Table 2), and can be derived from daily precipitation and daily mean, and maximum and minimum air temperature.
In this study, monsoon precipitation is defined as the precipitation falling during the months of June to September. The variability of precipitation is calculated as the coefficient of variation (CV), being the ratio of the standard deviation to the mean. To calculate the inter-annual variability, the mean and standard deviation are calculated from annual precipitation sums within a 30-year time slice. To calculate the intraannual variability, the CV is calculated for each year within a 30-year time slice using the mean and standard deviation of daily precipitation sums within that year. Subsequently, the 30 CV values are averaged. To calculate the maximum 5-day precipitation sum, the maximum sum of precipitation falling during five consecutive days within a climatological period of 30 years (i.e., a full 30-year time slice) is calculated. The reference evapotranspiration is calculated using the Modified Hargreaves equation, which requires daily maximum, minimum, and mean air temperature as input (Droogers and Allen 2002). The obtained values are averaged over the years within a 30-year period. Consecutive dry days periods are calculated for a climatological period of 30 years as the number of consecutive days with precipitation < 1 mm day −1 . Hot nights are defined as days with minimum air temperature > 21°C. The number of hot nights within each individual year is calculated, and the annual values are averaged for a 30-year period.  Figure 2 shows that for both precipitation and air temperature the biases of the uncorrected GCM data are drastically reduced by the quantile mapping approach and the downscaled data matches well with the observed dataset. For daily mean air temperature values, the remaining biases are very close to zero for means and extremes. The same holds for daily maximum and minimum air temperatures ( Supplementary Fig.  S1). At the high ends of the distributions of precipitation values, in particular for the 99th percentile values of precipitation, some bias remains. This is most likely due to the generation of new extremes by extrapolating the quantile mapping error function linearly outside the range of the empirical cumulative distribution functions. However, the remaining biases are in the order of several percents and the largest bias is still below 10% (RCP 4.5 CSIRO−Mk3−6−0 r4i1p1 model in the Ganges basin). We therefore conclude that the downscaled GCM data is suitable for our purpose to estimate future changes in climate change indicators.

Validation of GCM downscaling and bias-correction
Future climate change in the Indus, Ganges, and Brahmaputra river basins The first striking result is the difference in warming of the IGB regional average compared to the global average. A global average warming of 1.5°C with respect to PI levels would mean a warming of 1.4-2.6 (μ = 2.1)°C for the IGB (Table 1). At 2.0°C global average warming, the IGB would have warmed up by 2.0-3.4 (μ = 2.7)°C. This difference can be largely attributed to elevationdependent warming in the upstream IGB basins, i.e., the stronger warming of areas at high altitude compared to low-lying areas (Pepin et al. 2015;Palazzi et al. 2016    Our ensemble of representative scenarios for the IGB for RCP4.5 and RCP8.5, which may include more likely scenarios, indicates a regional temperature increase of 2.5-4.0 (μ = 3.4)°C and 4.4-6.8 (μ = 5.8)°C respectively by the end of the century, compared to PI levels.
For the 1.5°C scenario, temperature in the IGB is increasing ≈ 1.15°C with respect to the reference period (REF) (Fig. 3), meaning that temperature has already risen by ≈ 0.95°C between REF and PI. For the 2.0°C scenario, temperature increase in the IGB is ≈ 1.5°C with respect to REF. Temperature increases differ spatially, and are strongest in the lower Indus basin and on the Tibetan Plateau. Precipitation most likely increases for the entire IGB, with strongest relative increases in the most southern part of the Indus basin and southwestern part of the Ganges basin (Fig. 3). The southern part of the Indus basin is very arid (Fig. 1b), and thus even large relative increases do not result in large amounts of precipitation in the future. Furthermore, the range in the precipitation projections is large, especially for the areas where the mean of the projections indicates large precipitation increases.
For monsoon precipitation, the spatial patterns are similar to the changes in average annual totals of precipitation (Fig. 4a-c). Striking is the small decrease in monsoon precipitation in the eastern part of the Brahmaputra basin for the 1.5°C scenario. The intra-annual variability in precipitation decreases in the mountains, which may indicate increasing winter precipitation from westerly flows (Fig. 4d-f). For most of the lower basins, an increase in intra-annual variability is projected, which can be linked to the projected increases in monsoon precipitation. The inter-annual variability of precipitation likely decreases in areas with low inter-annual variability and increases in areas with already high inter-annual variability ( Fig. 4g-i).
Extreme precipitation (Fig. 4j-o) shows a very clear trend with increases for most of the IGB domain. Only the west of the Indus basin shows decreases in P 99 , which are stronger for the 1.5°C scenario than the 2.0°C scenario. The largest increases are projected for the Ganges basin (Table 3) and the very upstream parts of the Indus basin. Largest increases in 5-day sums of precipitation are projected for the Himalayan mountain range in the upper Indus basin. For the remaining areas, 5-day sums increase in general by 50-100%, with a few Fig. 3 Projected changes in mean annual air temperature (a, b) and mean annual precipitation (c, d) for the 1.5°C and 2.0°C scenarios with respect to the reference period (REF, 1981(REF, -2010. Grid values indicate the ensemble mean value and contour lines the ensemble ranges exceptions of decreases in 5-day sums. Reference evapotranspiration increases for most of the IGB domain as a result of increasing temperatures (Fig. 4p-r). The increases are mostly around + 0.1 mm day −1 , which does not make a large difference for the plains where reference evapotranspiration is in the order of 3-4 mm day −1 , but is substantial for the mountainous upstream basins where reference evapotranspiration is in the order of 1 mm day −1 . Stronger absolute increases are projected for the Ganges-Brahmaputra delta and western Indus basin, whereas a decrease in reference evapotranspiration is projected for the Indus delta.
Interesting is that the length of dry periods decreases for the southern Indus basin, which is currently very prone to droughts (Fig. 4s-u). The length of dry periods also decreases in the Karakoram mountain range, which may be related to increases in winter precipitation. In the Brahmaputra basin, a decrease in dry periods' length is projected, but for the Ganges basin, the length of dry periods increases (Table 3). These contrasting trends may be related to the contrasting trends in intra-annual variability of precipitation in those two basins (Fig. 4d-f).
The number of hot nights increases for all areas, except for the areas at high altitude, where the threshold for hot nights is not reached. The most drastic increases are in areas in the foothills of the mountains, where the threshold was not reached in the reference period, but will be reached in the future (Fig. 4v-x). We analyzed the changes in the number of hot nights specifically for five major cities in the region, with large, growing, populations, located in different parts of the region (Fig. 5). Remarkable is the observation that the difference between the 1.5 and 2.0°C scenarios is rather small. Increases are particularly drastic for the RCP8.5 EOC scenario. For relatively cool cities like Kathmandu, the occurrence of hot nights in the future will be a new phenomenon. For cities already facing severe heat stress, like Karachi, the number of hot nights per year can increase to up to 300 nights per year (Fig. 5d). Urban heat island effects have not been taken into account in our analysis, meaning that the found numbers may be underestimates. Figure 6 shows the difference between the precipitation and ET REF for both the kharif (monsoon) and rabi (winter) seasons. This is shown for two irrigated command areas located in the Indus and Ganges basins (Fig. 1a). The relation between P and ET REF is an indicator of water availability for the crops in the irrigated command area. The change in precipitation for the command area in the Indus basin is small and the climate effect on ET REF is small, P-ET REF is not expected to change much. The deficit likely increases slightly for the rabi season and decreases slightly for the kharif season, with roughly the same magnitude for each of the four scenarios. For the command area in the Ganges basin, water availability is likely to increase in the kharif season and this is mostly a result of a strong climate signal in the precipitation. The projected precipitation increases for the RCP4.5 and RCP8.5 EOC projections are significantly stronger than for the 1.5 and 2.0°C scenarios (Table 3), resulting in stronger increases in P-ET REF .
For the rabi season, the changes in the Ganges basin are similar to changes in the Indus basin. For both command areas, the uncertainty is much higher during the kharif season compared to the rabi season. This is because precipitation amounts are much higher during the kharif season compared to the rabi season, and therefore also the uncertainty in the precipitation amounts.

Impacts of changes
From the results, it is clear that significant changes in the region are imminent for each of the scenarios. The likely increases in monsoon precipitation and precipitation extremes could possibly lead to more frequent and more severe flooding in the region in the future. The changes in precipitation sums and extremes both show a reasonable linear correlation with temperature increase (Fig. 7b-g). Strongest increases are projected for the Ganges basin, indicating that future flooding could be a major issue in this basin, even for the 1.5°C scenario, while flooding is a major issue already in the present situation (e.g., the large floods during the 2017 monsoon season). A study comparing the impacts of 1.5 and 2.0°C global temperature increase for water availability and hydrological extremes in the Brahmaputra river concludes indeed that flooding events will be more frequent and more severe at 2.0°C global warming compared to 1.5°C (Mohammed et al. 2017). The increases in frequency and magnitude of precipitation extremes, and possible associated increases in floods, sedimentation, landslides, and GLOFs negatively affect human lives and livelihoods and ecosystems. Also, adverse effects on infrastructure can possibly be expected, including hydropower infrastructure, which is rapidly developing in the region (Siddiqi et al. 2012;Molden et al. 2014).
The overall likely increase in annual precipitation totals could imply overall increases in water availability. However, because interannual and intra-annual variability in precipitation in general increase for the Indus and Ganges basins, the increase in water availability could possibly be not directly beneficial. In the Ganges basin, there could possibly be more periods of water excess and more periods of water shortage, as indicated by the increase in extreme precipitation and increase in dry periods (Fig. 7f-h). Although the dry periods show an overall decrease in the Indus and Brahmaputra basin (Fig. 7h), Fig. 4 Values of climate change indicators for the reference period (REF, 1981(REF, -2010, and projected changes for the 1.5°C scenario (middle column) and 2.0°C scenario (right column) with respect to REF. The maps for the scenarios show the ensemble mean value this is partly misleading because in the Indus the decrease in dry periods is mainly observed in the upstream mountainous parts of the basin (Fig. 4t, u), where drought is less of an issue compared to the downstream areas, where most of the agricultural areas and people are located. Nevertheless, the length of dry periods in the downstream areas shows only minor increases for the 1.5 and 2.0°C scenarios, in line with results from a global study on changes in drought risk at 1.5 and 2.0°C scenarios (Lehner et al. 2017). Despite minor changes in the number of dry periods, the likely increases in reference evapotranspiration indicate increasing crop water demands, which could put further stress on water resources.
As the Indus, Ganges, and Brahmaputra rely to varying extent on water supplied by upstream melt water, future changes in snow and ice have to be considered here as well. With likely increases in dry periods and more erratic precipitation, water supplied by snow and glacier melt may become more important in the future. Kraaijenbrink et al. (2017) show that the ice volume in Asia will decrease by ≈ 36% under the 1.5°C scenario and ≈ 51-64% under RCP4.5 and RCP8.5 respectively. In particular, in the Indus basin, with high dependency on melt water (Lutz et al. 2014), this will affect future water availability. Further increases in the number of hot nights could affect human health further in a region where heat stress is already a major problem. Densely populated areas which have not encountered severe heat until now will likely also face this issue in the future, under each of the scenarios.
Concerning differences in impacts at 1.5 and 2.0°C global temperature increase, we can draw the general conclusion that changes in climate change indicators other than air temperature correlate linearly with temperature increase (Fig. 7). For some indicators, we even observe a somewhat flattening trend with temperature increase (Fig. 7d-j). This is especially pronounced for the change in interannual variability of precipitation, which shows little difference between the 1.5 and 2.0°C scenarios (Fig. 7d). In this case, for the Brahmaputra basin, the change becomes even negative for the 2.0°C scenario. As the change in interannual variability in this study depends on change in the Fig. 4 (continued) mean (and therefore total) and standard deviation of precipitation, this may be related to the strong increases in monsoon precipitation (Table 3,   of dry periods between the Indus and Brahmaputra on one hand and the Ganges on the other (Fig. 7h). On the whole, we can conclude that the regional impacts of climate change will be more severe for the 2.0°C than the 1.5°C scenario.
Much more important however are the large differences in changes between the 1.5°C and 2.0°C scenarios on the one hand, and the more likely RCP4.5 and RCP8.5 scenarios on the other hand. As regional temperature increases under the RCP4.5 and RCP8.5 end of century scenarios are significantly higher than for the 1.5°C and 2.0°C scenarios (Table 1), the changes in other climate change indicators, and their impacts in different sectors are also much stronger (Table 3). Current global warming is already close to 1°C with respect to pre-industrial temperature levels. Even considering the effects of emission mitigation policies, the chances that global temperature increase will be kept below 2.0°C are small. We therefore stress that the 1.5°C (or even the 2.0°C) scenario is not suitable for climate change adaptation planning, and it is therefore undesirable that the 1.5°C scenario will become a new paradigm. We emphasize that the formulation of robust climate change adaptation policies should be based on the full range of climate change projections.

Uncertainties and limitations
By using an ensemble of GCMs, we include the uncertainty in future climate in our assessment. This uncertainty is rather large, as illustrated in Figs. 3 and 5, and Tables 1 and 3. The changes in climate change indicators associated with a 1.5 or 2.0°C global temperature increase have wide ranges. This implies that care must be taken in interpreting the results. In our  Fig. 1 analysis, we mostly discuss changes based on the ensemble mean, and it has to be taken into account that there are large uncertainty ranges associated with the ensemble mean values. The uncertainty in future climate can be largely attributed to the general poor skill of GCMs in simulating the important regional climatological features of this region (Sperber et al. 2013;Palazzi et al. 2014;Ramesh and Goswami 2014;Sperber and Annamalai 2014;Mishra 2015). We reduced the uncertainty by only selecting GCMs with sufficient skill in simulating historical climate for the IGB region, but still large discrepancies in GCMs' ability to simulate the important features of the regional climate are evident (Lutz et al. 2016b).
Our ensemble of GCMs has been selected based on the range in the projections over the entire IGB. This scale issue may mean that the full range of projections is not captured for each separate area within the IGB. Because of the regional scale of the study, we analyze climate data at 10 × 10-km spatial resolution. Climatic processes at smaller scales, which in many cases dictate the actual impact, may be insufficiently represented.
Our study forms a first-order quantitative assessment of the difference in regional impacts of climate change at 1.5 and 2.0°C global temperature increase scenarios. We link quantitative changes in a set of ten climate change indicators to expected impacts for different sectors. Although we show that the changes in the analyzed climate change indicators increase in general linear with air temperature, on the impact side many processes may respond in a non-linear fashion to changes in climate change indicators. To quantify and verify the suggested impacts and differences in impacts of these scenarios for the different sectors requires additional sector-specific modeling studies in the fields of, i.e., hydrology, food and energy production, disease transmission, and natural hazards. It also requires inclusion of the projected future socioeconomic changes in the region. With an expected strong continued population growth, pressure on natural resources in the region will increase, as well as the number of people who are prone to climate change induced natural disasters. Our analysis of ten climate change indicators based on air temperature and precipitation does not include future changes in sea level rise, although this is an important climate hazard in the Ganges-Brahmaputra delta. A recent study investigating the differences in impacts from sea level rise resulting from 1.5, 2.0, and 3.0°C global temperature rise concludes for this delta that in unprotected areas the flooding depth is projected to double at 3.0°C compared to 1.5°C. In protected areas, this increase is approximately 50%. The average area inundated would increase 2.5 times (Brown et al. 2018).

Conclusions
South Asian river basins are vulnerable to climate change. A global temperature increase of 1.5°C with respect to pre-industrial levels would imply a ≈ 2.1°C temperature increase for the Indus, Ganges, and Brahmaputra river basins. Under a 2.0°C global temperature increase scenario, these river basins would warm up by ≈ 2.7°C. Our analysis of changes in climate change indicators under these scenarios shows that the region could face strong negative impacts of climate change for the occurrence of floods, agricultural production, hydropower production, and human health. The changes in climate change indicators are in general linearly increasing with temperature increase.
However, more likely scenarios for RCP4.5 and RCP8.5 imply regional temperature increases of ≈ 3.4 and ≈ 5.8°C at the end of the twenty-first century, and therefore possibly much stronger regional impacts. We stress that the 1.5°C (or even the 2.0°C) scenario is not suitable for robust climate change adaptation planning, and emphasize that the formulation of robust climate change adaptation policies should be based on the full range of climate change projections.