Projection of Future Precipitation and Temperature Change over the Transboundary Koshi River Basin Using Regional Climate Model PRECIS ()
1. Introduction
The Himalayan region is highly influenced by monsoon rainfall. About 80% of the annual total rainfall occurs during the monsoon season (June-September) in the eastern parts of the Himalayas [1] , but the role of monsoon decreases as it proceeds west and the influence of westerly winter disturbances becomes more prominent. Over western Himalaya, the contribution of winter precipitation is as much as 50% of the total [2] . Large populations across this region rely heavily on monsoon rainfall for their agricultural production, hydro-electricity generation and other livelihood activities [3] . In the global level, the increased emissions of greenhouse gases and short-lived climate pollutants (SLCP) are modifying short-term and seasonal climate variability [4] . Such changes can have profound impact on sectors like agriculture, water resources and human health in the Himalayan region [5] [6] . Variability of future climate is a matter of great concern for resources management and planning [7] [3] . Understanding the nature of climate change has thus become essential. Many climate modelling groups have performed future climate simulations under IPCC’s different emission scenarios [8] . The major sources of water in the Hindu Kush Himalaya (HKH)―snow and ice reserves are considered vulnerable to climate change [9] [10] . These models indicate that future rainfall and temperature is likely to change but the magnitude of these changes will be different in various ecological regimes. In the case of the Himalayas, too, the impact is likely to be different across the region given the sharp altitudinal difference across the mountains extending from east to west. As a result of the differences in precipitation climatology, the hydrological regime of the region is also quite different from east to west. Eastern Himalaya has significantly longer snow melt period compared to western Himalaya and global warming may increase already observed glacier retreat [11] . Eastern river basin―Brahmaputra―is heavily influenced by the summer monsoon whereas western river basin―Indus―relies on melting of snow and glaciers [12] [5] .
There is scant available literature on the Himalayan region that uses climate change scenarios based on high-resolution climate models [13] some previous studies on climate change scenarios include [14] [15] [16] [17] were developed over the sub-continental domain. For impact assessment, particularly in the water sector, scenarios developed over river basins have been found to be useful [18] . This study attempts to examine future changes in climate over the Koshi river basin using three simulation outputs from the regional climate model PRECIS based on the SRES A1B scenario [19] .
The Koshi River is the largest tributary of the Ganges and is one of the key transboundary river basins in the Hindu Kush Himalayan (HKH) region. Tributaries of Koshi river originates from the high-altitude Himalayas including trans-Himalaya (Tibetan plateau) and flow through Nepal’s mid-mountains, hills and plains (lowland region of Nepal and India). Koshi basin covers approximately 87,970 square kilometres; 32% of this area lies in China, 45% in Nepal and 22% in India [20] . Physiographic distribution based on altitude and major river system is shown in Figure 1. The basin is particularly prone to natural hazards such as glacial lake outburst floods (GLOF), landslide and debris flow, droughts and floods [21] . Neupane et al. [22] have suggested that climate change is likely to lead to an increase in water scarcity in the basin, which will ultimately influence the agro-based livelihoods of the basin. The temperature analysis by Agarwal et al. [23] suggested that seasonal as well as mean annual minimum and maximum temperature will increase in the future. The maximum temperature was reported to be increasing at the rate of 0.058˚C/year in the Nepal part of the Koshi basin [23] . The precipitation trend is spatially less homogenous compared to temperature [24] . To better understand the future
Figure 1. Koshi basin location map with three physiographic zones (trans-Himalaya, middle mountains and southern plain) and major river system.
precipitation and temperature patterns, this study attempts to examine future climate change scenarios in the basin using PRECIS outputs.
In the context of climate change, the hazards may increase in magnitude and frequency [25] [3] . Additionally, the changes in physical processes and management practices in upstream areas might affect water availability in downstream areas [26] .
Due to the complex topography of the Himalaya, rainfall distribution over the region is highly heterogeneous [27] [28] [29] [30] . The sudden rise of the Himalayas from the Terai plains within a span of a few kilometres results in various types of climate [31] [32] [33] , and can influence the atmospheric flow pattern through orographic effects. Understanding spatial and temporal distribution of climate variables for the present and future is essential for understanding future impacts on different sectors, including water resources, agriculture and other sectors. The Koshi basin covers three different physiographic zones: trans-Himalaya in the northern Tibetan part of China; middle mountains and high Himalayas (hereafter referred to as middle mountains) in the centre (Nepal); and the southern plains of Nepal and India, as shown in Figure 1.
2. Data
Complex climate regimes and geographical features in the region poses a considerable difficulty in simulating observed climate scenarios. Panday et al [13] suggests that high resolution dynamical downscaling requires techniques that improves better representation of topography to simulate monsoon and other hydrological process in the region.
As Kripalani et al. [34] have noted, only a few GCMs (out of 22) effectively captured the spatial variability of the Indian summer monsoon. Among these models, HadCM3 (along with six other models) has simulated the inter-annual summer monsoon mean and variability well. The three QUMP simulations namely, Q0, Q1 and Q14 were able to simulate gross features of the Indian summer monsoon. These three output have been validated for their skill in simulating the present climate over India [35] . Further, outputs have been used to project the future climate over India [15] , Himalayan region [17] and Indus basin [18] . In this paper, the outputs of these three simulations were used to project the changes in future climate from the base period 1961-1990 over the Koshi basin. These simulations, run at IITM (Q0, Q1 and Q14), have been configured to 1.5˚N to 38˚N and 56˚E to 103˚E of whose output were provided to International Centre for Integrated Mountain Development (ICIMOD).
For the current study, APHRODITE’s (Asian Precipitation - Highly - Resolved Observational Daily Integration Towards Evaluation http://www.chikyu.ac.jp/precip/) gridded data has been taken as the base data for evaluation of the seasonal rainfall simulation. These gridded datasets are publicly available at a resolution of 0.25˚ for the monsoon Asian domain (60˚E - 150˚E and 15˚N - 55˚N). APHRODITE’s rainfall data is based on rain-gauge observations and depict the areal distribution and variability of the rainfall over the Himalayas [36] . APHRODITE rainfall product has been widely used for scientific assessment in the Himalayan region [37] . It is found to be better than other gridded products such as TRMM and gives better precipitation estimates [38] [39] . For quantitative and qualitative comparison, we have used version APHRO_V1101, which includes more rain-gauge data than the previous version [40] .
To evaluate the temperature outputs of the PRECIS model, we have taken APHRODITE’s mean temperature publicly available for monsoon Asia domain from 1961-2007 [41] and maximum and minimum temperature from Princeton University Hydroclimatology Group Bias Corrected Meteorological Forcing Datasets [42] available globally from 1947-2007 at 1.0˚ × 1.0˚ resolution. Further, the Princeton data is of coarse resolution so for the present study, we have used daily temperature range from this dataset and applied it to APHRODITE's mean daily temperature to build maximum and minimum temperatures of 0.25˚ resolution.
3. Results and Discussion
3.1. Precipitation
To evaluate the model in representing the regional climatological features, summer monsoon (June-September) rainfall characteristics have been studied. Table 1 provides seasonal rainfall statistics for the three simulations for the baseline period (1961-1990) over the whole basin, the trans-Himalayan part, the middle mountain area and the southern plains. The rainfall statistics are computed from area averaged daily rainfall and temperature over transboundary Koshi basin and its three sub-divisions. Rainfall simulations for the whole basin compared with observed (APHRODITE) resulted in a substantial wet bias over
the basin. Of the three simulations, annual total for Q1 simulation was found to be 1632 mm with 201 mm standard deviation closest to the observed value of 1024 mm. Q0 and Q14 simulations estimated annual total of 1663 mm and 1676 mm with standard deviations of 186 mm and 162 mm, respectively. Breaking down the basin into three zones shows that Q1 simulation with 1151 mm is closest to the observed value 1085 mm over the southern plains. A substantial wet bias is simulated by PRECIS over the trans-Himalayan area (rain shadow area). The estimation here is as much as three times the observed (Table 1(d)) value.
Ensemble average of the three simulations shows an overestimation (bias) of 1.75 mm per day. Several previous studies have reported a wet bias over the South Asia domain [14] [16] [18] . A PRECIS simulation output study over India showed Q0 and Q14 simulations to have a wet bias and Q1 a dry bias [15] , whereas other similar studies over India showed a wet bias [14] .
Monthly precipitation distribution for the whole basin and its breakdown into the three physiographic zones is provided in Figure 2. The model captures the monthly accumulation pattern. For July and August, the model results and the observed data were quite close, but for the rest of year, the model overestimated the observed total. Over the middle mountains and the southern plains, the observed data exceeded the model total for July and August, and for the remaining months, the model totals were higher.
Rainfall distribution for the monsoon season (June-September) is provided in Figure 3. Central parts of the basin (middle mountains) where the rainfall maximum is observed and the north-south gradient is well captured by the
Figure 2. Annual cycle of mean rainfall. (a) Whole Koshi basin (b) southern plains (c) middle mountains (d) trans-Himalaya.
Figure 3. Seasonal distribution of rainfall (June-September) in mm. (a) APHRODITE (b) Q0 PRECIS simulation, (c) Q1 simulation (d) Q14 simulation.
PRECIS simulation. However there is a wet bias extending slightly to the north part of the basin. Winter season spatial rainfall distribution (figure not shown) also shows over estimation by the simulation while the north-south gradient in rainfall is captured.
3.2. Rainfall Frequency and Intensity
A day with rainfall greater than or equal to 1.0 mm is considered a rainy day. Rainfall frequency is higher over the middle mountain areas and progressively decreases from north to south (Figure 4). The PRECIS simulations overestimated the number of rainy days but the distribution pattern is similar.
Figure 4. Number of rainy days as simulated by PRECIS compared with the observed. (a) APHRODITE (b) Q0 simulation (c) Q1 simulation (d) Q14 simulation.
APHRODITE estimated 169 annual rainy days for the baseline period 1961-1990, whereas PRECIS estimated 200, 189, and 203 days for the Q0, Q1, and Q14 simulation respectively (Table 1(a)). Rainfall intensity pattern on a rainy day (Figure 5) is higher over the central parts of the basin and lower in the northern parts. However, in the southern plains, rainfall intensity is between 6 - 10 mm/day in PRECIS simulations compared to the 8 - 12 mm/day for APHRODITE.
Figure 5. Simple daily rainfall intensity on a rainy day (mm/day). (a) APHRODITE (b) Q0 (c) Q1 (d) Q14 PRECIS simulation.
3.3. Temperature
Seasonal and annual temperature (˚C) and its standard deviation is provided in Table 2 for the whole basin and the different physiographic zones viz., southern plains, middle mountains and trans-Himalaya. The annual cycle of mean monthly temperature is shown in Figure 6. The annual march of temperature is
Table 2. Seasonal and annual average temperature (˚C) and standard deviation for the whole basin (top left), southern plains (top right), middle mountains (bottom left) and trans-Himalayan area (bottom right).
Figure 6. Annual cycle of mean monthly temperature (a) Koshi basin (b) southern plains (c) middle mountains (d) trans-Himalaya.
well captured by the model. The Q0 and Q1 simulation values were quite close to the observed, while the Q14 simulation shows a warm bias during April and May. In the middle mountain areas, the PRECIS underestimates the observed data throughout the year, although the gap is little for April and May. In the southern plains, PRECIS overestimates the pre-monsoon season, whereas the observed data overestimate the rest of the time period. In the trans-Himalayan zone, the simulated and observed data match well, although PRECIS slightly underestimates the observed data, except for the monsoon season. The temperature is underestimated from July for the three simulations. These biases may be due to a coarse network of monitoring stations over the high-altitude mountainous areas, as the stations are generally located in the valley bottoms.
Summer and winter distributions of air temperature as simulated by PRECIS are given in Figure 7 and Figure 8 respectively. Temperature distribution over the Koshi basin closely follows the topography and is in general agreement with the observed, but with a cold bias over the trans-Himalaya and a slight warm bias towards the southern plains.
3.4. Highest Maximum and Lowest Minimum Temperatures
Baseline simulations of the extreme maximum and minimum temperatures are shown in Figure 9 and Figure 10 respectively. They closely match the observed pattern. Both the highest maximum and lowest minimum show a warm bias over all parts of the basin. The model simulates the highest temperature exceeding 50˚C over the southern plains (mostly over India) whereas in the observed, it is seen over the southernmost part only.
4. Projection of the Future Climate
The continuous output from 1961 to 2098 from PRECIS based on IPCC SRES A1B scenario are divided into baseline period (1961-1990), and future projections from 2011 to 2098 are divided into three time slices 2011-2040 (near future), 2041-2070 (middle of the first century) and 2071-2098 (end of the century).
Figure 7. Spatial distribution of summer season mean temperature (˚C) (a) APHRODITE compared with three PRECIS simulations (b) Q0 (c) Q1 and (d) Q14.
Figure 8. Winter season mean surface air temperature in ˚C. (a) APHRODITE (b) PRECIS Q0 simulation (c) Q1 simulation (d) Q14 simulation.
Their results are discussed below.
4.1. Projected Change in Normal Rainfall
Projected changes in summer rainfall for the whole basin and the three physiographic zones for the three time slices in the future viz., 2011-2040, 2041-2070, and 2071-2098 are given in Figure 11. Over the whole Koshi basin, all three simulations show decrease in rainfall in the near future and then gradual increase till the end of the century. However over the southern plains, except for Q14 in the 20’s time slice, there is progressive increase in summer rainfall. Over the middle mountain zone, rainfall is projected to decrease in the near future and then progressively increase from the present baseline towards the end of the century. The rainfall is projected to decrease from the baseline period over the trans-Himalayan zone.
Figure 9. Highest maximum temperature (˚C). (a) APHRODITE (b) Q0 simulation (c) Q1 simulation (d) Q14 simulation.
Figure 10. Lowest minimum temperature (˚C). (a) APHRODITE (b) Q0 simulation (c) Q1 simulation (d) Q14 simulation.
Spatial distribution of projected changes in summer and winter season rainfall for three simulations are provided in Figure 12 and Figure 13 respectively. During summer season, Q0 and Q1 simulation projected increase in rainfall, mostly up to 10% over the plains area and southern part of the middle mountains, and decrease over the northern part of the middle mountains and the trans-Himalayan area. Q14 simulation projected decrease in rainfall over almost
Figure 11. Projected changes in summer rainfall for three time slices in the future (a) whole basin (b) southern plain (c) middle mountain, and (d) trans-Himalaya.
all parts of the basin, mostly up to 10% in the south and more in the north. Spatially, there is general agreement between all three simulations on rainfall projection for the 50’s time slice. During the 80’s, Q1 simulation projected a slightly larger area of decrease in rainfall over the northern part of the middle mountains and the trans-Himalaya. For the 20’s, all three simulations show decrease in rainfall during winter in the entire basin, with more than 50% decrease over the southern part of the basin by Q1 scenario. However during the 50’s, the result is mixed: Q0 projects increase in the entire basin and Q1 projects decrease, whereas Q14 projects increase over the middle mountains and trans-Himalaya and decrease in the south. Towards the end of the 21st century, there is projected increase in rainfall over the north and decrease in the south.
4.2. Rainfall Frequency and Intensity
Projections of future changes in rainfall frequency or the number of rainy days and simple daily rainfall intensity (mm/day) are shown in Figure 14 and Figure 15 respectively. Except for a small area over the eastern part of the middle mountain area, all three simulations show decrease in rainfall frequency by up to five days. Over the southern plains, further decrease in the number of rainy days is projected for all three simulations. During the 80’s, there is decrease by 10 - 25 days, 5 - 20 days and up to 15 days for Q0, Q1 and Q14 simulations respectively.
Figure 12. Projected change in mean summer monsoon rainfall (%) with respect to baseline 1961-1990 in the 2020s (left column), 2050s (middle column), and 2080s (right column) for the three simulations Q0 (top row), Q1 (middle row), and Q14 (bottom row).
Figure 13. Projected change in mean winter season rainfall (%) with respect to baseline 1961-1990 in the 2020s (left column), 2050s (middle column), and 2080s (right column) for the three simulations Q0 (top row), Q1 (middle row), and Q14 (bottom row).
Figure 14. Projected change in the number of rainy days in a year with respect to baseline 1961-1990 in the 2020s (left column), 2050s (middle column), and 2080s (right column) for the three simulations Q0 (top row), Q1 (middle row), and Q14 (bottom row).
Figure 15. Projected change in simple daily rainfall intensity (mm/day) with respect to baseline 1961-1990 in the 2020s (left column), 2050s (middle column), and 2080s (right column) for the three simulations Q0 (top row), Q1 (middle row), and Q14 (bottom row).
On the other hand, a slight increase (up to 2 mm/day) in rainfall intensity is projected over the southern plains, while a decrease over northern area (middle mountains and trans-Himalayan area) is projected for the near future, as shown in Figure 15. But towards the end of the century, with the exception of a small area over the middle mountains and trans-Himalayan border, rainfall intensity is projected to increase by about 2 mm/day compared to the baseline, indicating some increase in high intensity rainfall in the future.
4.3. Projected Changes in Temperature
Simulations of surface air temperature show progressive rise in both the maximum and minimum temperatures through the given time slices. Table 3 provides averages of the three simulations for summer maximum and winter minimum temperature anomalies against the baseline over the whole basin and its three physiographic zones for the three time slices. These temperature anomalies are the thirty years averages difference between baseline and the projection period. The change in temperature from the baseline (1961-1990) for the whole Koshi basin and its different physiographic zones are shown in Figure 16. Again progressive warming is projected for both summer and winter seasons. Average warming (Q0, Q1 & Q14) from the base period for the three time slices over whole the Koshi basin are 1.2˚C, 2.8˚C and 3.7˚C for maximum temperature and 2.0˚C, 3.7˚C and 5.3˚C for minimum temperature during winter season. In all the simulations and over all the physiographic zones, more warming is projected for winter than for summer. Towards the end of the century, there is more warming over the trans-Himalaya during summer and over the southern plains during winter. Spatial distribution of annual mean temperature is given in Figure 17. Although there is no definitive pattern in warming, progressive warming from the near future (1˚C - 2˚C), middle of the century (2˚C - 4˚C) and in the end of the century (3˚C - 6˚C) is simulated. Rajbhandari [20] has reported similar projections over the Koshi basin using CMIP5 GCMs.
4.4. Extreme Temperature
Three simulation results indicate that daily extreme temperature may intensify in the future. The spatial pattern of the changes in the highest maximum temperature suggests warming of 1˚C - 2˚C (Figure 18) in the 2011-2040 period,
Table 3. Projected seasonal maximum and minimum temperature anomalies (˚C) against the baseline over the Koshi basin and its physiographic zones for three time slices.
Figure 16. Changes in future projection of the temperature (˚C) from the baseline over the whole Koshi basin (top row), southern plains (second row), middle mountains (third row) and trans-Himalaya (fourth row). Left panel―summer season maximum temperature and right panel―winter season minimum temperature.
Figure 17. Projected change in annual mean surface temperature (˚C) with respect to baseline 1961-1990 in the 2020s (left column), 2050s (middle column), and 2080s (right column) for the three simulations Q0 (top row), Q1 (middle row), and Q14 (bottom row).
Figure 18. Projected change in highest maximum temperature (˚C) with respect to baseline 1961-1990 in the 2020s (left column), 2050s (middle column), and 2080s (right column) for the three simulations Q0 (top row), Q1 (middle row), and Q14 (bottom row).
which may exceed up to 6˚C over the higher Himalayan parts of the middle mountains in the distant future.
5. Conclusions
Presented below is the assessment of expected future changes in rainfall and temperature over the Koshi basin under the IPCC’s SRES A1B global warming scenario based on the PRECIS model.
・ The PRECIS model showed ability to simulate climate scenario for both seasonal rainfall and temperature reasonably well but there were quantitative biases. There was an overestimation of rainfall for both summer and winter seasons and a cold bias for temperature.
・ Rainfall is expected to decrease in the near future (2011-2040) and then progressively increase towards the end of the 21st century.
・ Over the southern plains, summer monsoon is expected to increase from the base period towards the end of the 21st century.
・ Over the trans-Himalayan area, summer monsoon rainfall is projected to decrease.
・ The projected change in the number of rainy days differs among the ensemble; however, in future, the number of rainy days may decrease and there may be increase in rainfall intensity over the southern plains.
・ Progressive increase in maximum and minimum temperatures is projected for the future. All three simulations indicated a significant rise in temperature (between 4˚C - 6˚C) throughout the basin.
・ Analysis of extreme temperature indicates that daily maximum temperature may be more intense in the future.
The above assessment is indicative of the range of expected changes in rainfall and temperature from the PRECIS baseline. It is important to note that projections contain significant uncertainties. The study was limited as it was based on a few simulations from a single high-resolution regional climate model.
Acknowledgements
This study was conducted under the Koshi Basin Programme at ICIMOD, funded by the Department of Foreign Affairs and Trade (DFAT) of Australia. This study was partially supported by the core funds of ICIMOD (contributed by the governments of Afghanistan, Australia, Austria, Bangladesh, Bhutan, China, India, Myanmar, Nepal, Norway, Pakistan, Switzerland and the United Kingdom). The views and interpretations in this publication are those of the authors and not necessarily attributable to ICIMOD.