Pakistan Water and Power Development Authority (WAPDA) Climate Change Impact on Water Resources of Pakistan Glacier Monitoring & Research Centre Muhammad Arshad Pervez Project Director (GMRC)
Outline of Presentation Introduction Significance of Upper Indus Basin GMRC and its Present Role in UIB Monitoring Activities of Glacier Monitoring & Research Centre Overall Progress and Future Plans 2
Introduction Pakistan is an agriculture based country and is highly reliant on the Indus irrigation system fed by melt water originating from large scale blue gold glaciers reserves. 3
Significance of Upper Indus Basin Glaciers melting under climate change is emerging as one of the biggest challenge for planners and managers of water resources and agriculture around the world The Indus basin in Pakistan contains the largest continuous irrigation system in the world, and millions of people are dependent on the food produced in this area The Karakoram mountain range plays a crucial role in sustaining the Indus basin s water resources through glacial and snow melt hydrological processes Upper Indus Basin (UIB) 4
Flows (MAF) Contribution of Snow and ice 70 60 50 40 30 85% Flows Ice & Snow Contributions Contribution of Snow & Ice in river inflows 70% for entire Indus Basin and 85% for Indus Main Stem Any changes in glaciers behavior will have direct impacts on agriculture and food security in Pakistan 20 10 0 75% 65% 50% Therefore it has become more important to have better knowledge about the status of cryosphere Indus at Tarbela Kabul at Nowshera Jhelum at Mangla Chenab at Marala
Role of GMRC-WAPDA in UIB Monitoring High-Altitude Hydro-met Network Hydrologic Modelling Field Investigations 20 AWS UBC Watershed Model Glacier field monitoring - Stake surveys - Snout surveys 15 + 15 AWS 13 ARL GR4JSG & SRM Isotope analysis Permafrost monitoring GLOF Studies 6
Temperature Trend in Upper Indus Basin (Man Kendall) (1995-2015) Month Shangla Rattu Rama Shogran Ushkore Hushey Kelash Zani Pass Yasin Khod Pass Shendure Ziarat Burzil Khunjerab January 0.62 0.15-0.34 1.22 0.95 0.53 0.46 0.93 0.70 0.30 0.86 0.81 1.09 2.17 February 0.62 0.00 0.76 0.00 0.84 1.19-0.63 0.98 0.25-0.69 0.37 2.12 1.75 1.89 March 2.20-0.04 0.04 2.09 1.33 1.82 1.37 2.44 0.82 0.99 1.40 1.49 1.12 2.80 April 0.89-0.08-0.15 1.22 0.70-0.53 0.35 0.06 0.08 0.30 0.31 0.18-0.53 0.70 May 0.34 0.45 0.72 0.24 0.74 0.77 0.07 0.07 1.11-0.20 0.73 1.31 0.04 0.35 June 0.14 0.50 0.59-0.92 0.38 0.08 0.91 0.06 0.86-0.45 0.79 0.77 0.19 0.53 July -0.96-0.83-0.38-2.21 0.00-0.46-1.65 0.71-0.05-0.40-1.47-1.85-0.07 0.14 August -0.42-0.23-0.33-0.62-0.74 0.21-2.00-2.06-1.13 0.00-0.86-2.57 1.09-0.28 September 0.14-2.08-2.16-2.01-2.10-1.79-2.10-1.44-1.90-0.69-1.65-2.40-1.86-0.98 October -0.21-0.30-0.41-0.24 0.77-0.04 0.00 0.79 0.12 0.00-0.43 0.45 1.05 1.40 November -0.55 0.61 0.37 0.61 0.39 0.00 0.00 0.11-0.04 0.82 0.18 0.99 0.28 1.54 December -0.62-0.19 0.12 0.31 0.07 0.95-1.23 0.67-0.25 0.93 0.43 0.99 0.00 1.26 ANNUAL 0.48-1.21-0.98-0.31 0.14 0.00-1.30-0.40 0.04 0.89 0.43 1.40 0.63 1.54
Temperature Change in UIB ( o C ) (1995-2015) Month Shangla Rattu Rama Shogran Ushkore Hushey Kelash Zani Pass Yasin Khod Pass Shendure Ziarat Burzil Khunjerab January 0.69 0.45-0.36 1.33 0.95 0.90 0.63 0.70 0.67 0.75 0.90 0.65 1.29 2.59 February 1.80 0.15 1.03 0.10 1.02 1.07-0.95 1.74 0.43-1.35 0.29 1.65 1.83 3.05 March 2.18 0.00 0.00 2.03 1.38 1.36 2.19 3.30 0.85 1.43 2.17 1.15 0.84 2.45 April 1.41 0.00-0.30 2.97 0.83-0.53 0.43 0.13 0.14 0.17 0.51 0.13-0.40 1.13 May 2.18 1.29 1.54 0.59 1.45 1.24 0.47 0.49 1.42-0.75 1.08 1.24 0.05 0.71 June 0.30 0.67 1.05-0.54 0.80 0.13 0.82 0.22 1.64-0.50 1.91 1.43 0.21 0.44 July -0.49-0.68-0.74-0.78 0.00-0.48-0.95 0.47-0.13-0.32-1.56-1.10-0.10 0.08 August -0.03-0.09-0.19-0.35-0.57 0.24-1.37-1.12-1.60 0.00-0.64-1.83 0.47-0.38 September 0.08-2.30-2.22-1.44-2.22-2.34-1.97-1.34-2.08-1.88-1.95-1.71-1.27-1.30 October -0.41-0.45-0.62-0.56 0.75 0.00 0.00 1.43 0.36 0.27-0.41 0.23 1.27 1.94 November -1.15 1.00 0.27 0.53 0.46-0.06 0.00 0.93-0.05 2.98 0.33 1.10 0.38 1.96 December -0.80-0.24 0.09 0.45 0.32 1.40-2.38 1.23-0.10 2.75 0.72 1.06-0.09 1.12 ANNUAL 0.47-0.44-0.56-0.16 0.14 0.01-0.53-0.46 0.05 0.53 0.20 0.68 0.19 0.92
Conclusion of Temperature Analysis There is no significant change in the mean annual temperatures of the stations located in the Upper Indus Basin Most of the significant falling trends were seen in the months of July and August which are glacier melt periods, Similarly September also showed significant falling trends for most of the stations. For most of the stations, winter months showed rising trends in which some were significant for January and February.
Flow Trends of Indus River at Tarbela 1962-2015 Month First year Last Year N Man-Kendall Statistic Total Change (MAF) January 1962 2015 54 1.35 0.04 February 1962 2015 54 2.91 0.11 March 1962 2015 54 2.22-0.09 April 1962 2015 54 0.51-0.45 May 1962 2015 54 3.04-0.70 June 1962 2015 54-1.63-1.01 July 1962 2015 54-1.42-1.91 August 1962 2015 54-1.49 0.14 September 1962 2015 54-0.12-0.31 October 1962 2015 54 2.10 0.18 November 1962 2015 54 4.13 0.22 December 1962 2015 54 3.00 0.13 Annual 1962 2015 54-0.67-0.13
Monthly Mean Annual Flow Volume Flow Trends of Indus River at Tarbela 1962-2015 8.00 7.00 6.00 5.00 4.00 3.00 2.00 1.00 0.00 1960 1965 1970 1975 1980 1985 Year Data Sen's estimate
Conclusion of Flow Analysis The flow data of 1995-2005 showed an insignificant falling trend in the mean annual analysis. The winter months showed a rising trend with November as significant. Summer months showed an insignificant falling trend. For the months of February, March and May, there was a significant rising trend for the flows
UIB Glaciers Selected for Mass Balance Studies 13
Ablation Gradient at Passu Glacier # Description August 1990 August 2011 August 2012 August 2013 Changes in Snout at Passu Glacier from 2011 to 2016 1 Net ablation at PS profile (NAS) cm/day 2 Net ablation at PG Profile (NAG) cm/day 10.3 9.14 6.0 7.5 6.7 6.91 4.9 6.1 3 ABL = (NAS-NAG) cm/day 3.6 2.23 1.1 1.4 4 Elevation PS (ES) m.a.s.l 2850 2677 2635 2634 5 Elevation PG (EG) m.a.s.l 3500 3436 3428 3426 6 ELV (EG-ES) m 650 759 793 729 7 Ablation Gradient cm/day/100- m between the PS & PG profiles 0.554 0.294 0.14 0.18 Findings Glacier has receded in 2012 relative to its position in 2011 A part of glacier at its terminus detached from the main body of the glacier and melted during 2012 The glacier snout position was unchanged in 2013 The comparison of the snout survey from 2012-13 to 2015-16 showed that terminus of the Passu glacier has advanced about 20 m 14
Ablation Rate (cm/day) Map of Barpu Glacier Melt Rates for Barpu Glacier (Jul-Sep, 2014) 8.0 7.0 6.0 5.0 R² = 0.7541 Barpu Glacier Snout Barpu Glacier Data shows that the ablation rate drastically reduces with an increase in the thickness of the debris cover, especially over 2 cm 4.0 3.0 2.0 1.0 0.0 R² = 0.9018 0 2 4 6 8 10 12 14 Debris Thickness (cm) Relationship between debris thickness and ablation rates 15
Melt Rates for Biafo Glacier (Jul-Sep, 2015) Net Ablation Sr. No. Profiles Month Rates 1 July (11-31) 6.8 Mango Profile 2 August 5.2 3 July (11-31) 6.5 Baintha Profile 4 August 5.0 Findings Ablation rate at Baintha profile is same i.e., 6 cm/day, as compared with the ablation rate of the study in 1985 The spatial and temporal ablation rates didn t show any significant difference at both profiles instead of very higher difference in temperatures (on-ice) At Mango, for example, the average temperature (on-ice) was 25% higher but the ablation rate was only 4.5% more than the Baintha during July Similarly, the glacier velocity at Baintha also was same as measured during 1985 These latest results (i.e. from the data collected during 2015) reveals that no significant change in ablation rate and glacier velocity has occurred as compared with the results of 1985 16
Hydrological modelling studies on Source Platform with collaboration of CSIRO,AUSTRALIA Correlation 0.83 Volume Bias -7.42 Results of Flows calibration at Dainyor Bridge, Hunza River Basin 17
Elevation (m.a.s.l) 122.4 m 139.3 m 114 m Elevation (m.a.s.l) 81.8 m 105.5 m 86.9 m Raikot Glacier Depth Profiles (2016) 3400 3380 Profile 1 PR1-1 RC-3 PR1-2 3360 3340 3320 3300 3280 3260 3240 Glacier Surface Glacier Bed 0 100 200 300 400 500 600 700 Distance from Glacier's Left Edge (m) 3550 Profile 2 PR2-1 RC-6 PR2-2 3500 3450 3400 3350 Glacier Surface Series2 0 200 400 600 800 Distance from Glacier's Left Edge (m) 18
Glacier Snout Surveys Snout Surveys planned: 50 Conducted for 37 glaciers so far Monitoring and mapping of fluctuations in their terminus They will serve as a strong control points on the ground for the verification of remotely sensed and interpreted images 19
Isotope Study for Indus River above Tarbela Samples have been collected and will be sent to China for lab analysis 20
Permafrost Monitoring A permafrost monitoring equipment has been installed in Deosai DCP Station. 21
GLOF Study Glacier Monitoring Research Centre (GMRC) was continuously monitoring Khurdopin Glacier and the Lake area through field visits and RS analysis from 2015 and two reports were generated in this regard GMRC was the only entity to reach the lake physically. 22
1 st Visit Observations The Khurdopin Glacier has advanced and raised the river bed that has resulted the formation of glacial lake in the upstream side by blocking the flow of Virjerab Glacier Ice dammed lake was formed which was nearly 40-50m long, 8-10m wide 4-6m deep With the help of RS/GIS it was observed that due to further surging of glacier the length of lake was increased about 400m till mid of June The water continuously draining out through the snout of the glacier and the lake disappeared in July-August 2017 Recommendations Another assessment of Khurdopin glacier was suggested in the month of Sep 2017 so that the probable impact of newly developed lake may be predicted Ice blocks were falling from Khurdopin Glacier into the lake blocking location but the flow of the river overtopped these ice blocks and washed away them within hour or two. Awareness sessions can be helpful to reduce the mental stress of local people and calamity from the GLOF event If the Khurdopin Glacier continues to advance then there is a chance of the large glacial lake development that will be hazardous to downstream inhabitants, vegetated areas and public infrastructure 23
2 nd Visit Observations GMRC Team revisited the lake for the assessment in Sep-Oct 2017 The lake, developed due to surging of glacier as observed in May-June 2017, has been drained The Khurdopin Glacier was further advanced towards the opposite mountain as compared with the previous observations made in May-June 2017 and blocked the river path completely at some locations The water continuously draining out through the snout of the glacier and the lake disappeared in July-August 2017 Recommendations Drainage of the lake by controlled breaching for the formation of spill channel to increase the drainage rate thereby reducing the risk of high flow Dam Breach Study to assess the situation D/s in case of a GLOF event Safe evacuation of the inhabitants to avoid any loss of lives The advancement of Khurdopin Glacier continued in 2018 as well and due to complete blockage of the stream, the lake was again started developing in November 2017. It was apprehended that blockage of crevasses due to falling debris may block the flow subsequently may result in development of big lake which may cause high flood event in coming summer season As a result of glacier surge the process of deformation has been accelerated and large numbers of crevasses have been developed in the glacier Continuous monitoring of lake formation using RS/GIS and physical verification at lake site The breakage of glacier and water flow through crevasses may result in loss of mass at snout due to the detachment of lowermost portion of glacier snout Early Warning Arrangement 24
RS/GIS Analysis Development of Lake from 13.05.17 to 29.03.18 using Landsat 8 Satellite Images downloaded from http://earthexplorer.usgs.gov/ an approximate length of lake 2.5 3.5 Km and ~ 175m in width, which is very close to what the local of the area informed telephonically. 25
NDMA Initiative NDMA took initiative and form a ToE. Two experts Mr. Mujahid Tanveer Ahmed (Director Glaciology) and Mr. Syed Hammad Ali (Glaciologist) from GMRC WAPDA As a result of that meeting, a Heli visit was arranged by NDMA. observed that water level of the lake has been lowered due partial drainage from the lake but still there is every chance of increase of lake level as the upstream discharge increases. It will be useful to make arrangements on urgent bases to clear the path of water through controlled breaching at blocked positions, if Khurdopin glacier continue surge in the same direction it may result in formation of big Glacial Lake. In this regard April to June 2018 may be critical months and situation may face a GLOF event most probably in July-August, 2018 26
Overall Progress and Future Plans High-Altitude Weather Stations It is planned to further enhance the AWS network at high altitude. ARL For Glacier Fed Streams 13 ARL are planned to install in the glacier fed streams for the estimation of future water availability Glacier Field Monitoring Three (3) of the major Glaciers have been investigated for stake ablation out of five (5) selected for the Mass-Balance studies and Thirty Seven (37) out of fifty (50) of them for the monitoring and mapping of fluctuations in their terminus. Hydrological Modeling Flows for Hunza River Basin has been calibrated. The flows for other sub-catchments will be calibrated and future climate scenarios will be used to estimate the availability of flows from the sub-catchments of UIB. (C-4 funding) 27
Existing and Proposed GMRC AWS and ARL Stations 28
THANKS 29