Estimating the characteristics of runoff inflow into Lake Gojal in ungauged, highly glacierized upper Hunza River Basin, Pakistan

Abstract

Motivated by the potential flood outburst of Lake Gojal in the ungauged highly glacierized (27%) upper Hunza River Basin (HRB) in Pakistan that was dammed by a massive landslide on 4 January 2010, we attempt to analyze the characteristics of water inflow to the lake employing remote sensing data, two hydrological models, and sparsely observed data. One of the models (Model I) is a monthly degree-day model, while another (Model II) is the variable infiltration capacity (VIC) model. The mixture of glacier runoff output from Model I and runoff over unglacierized areas calculated by Model II has a similar seasonal variation pattern as that estimated from data recorded at a downstream station. This suggests that glacier runoff is the main source (87%) of runoff inflow into the lake. A sensitivity analysis suggests that the water inflow to the lake is highly sensitive to an increase in air temperature. Runoff in May is predicted to sharply increase by 15% to more than two-fold if the air temperature increases by 1 to 7, but it is predicted to increase only from 9% to 34% if the precipitation increases by 10% to 40%. The results suggested that the water inflow into Lake Gojal will not sharply rise even if there is heavy rain, and it needs to be in caution if the air temperature sharply increases. Analysis on long-term air temperature record indicates that the water inflow into the lake in May 2010 was probably less than average owing to the relatively low air temperature. Consequently, the flood outburst did not occur before the completion of the spillway on 29 May 2010.

Keywords

ungauged / glacierized / runoff modeling / sensitivity analysis

References

[1] Archer D R. Contrasting Hydrological Regimes in the Upper Indus Basin. Journal of Hydrology, 2003, 274: 198-210.

[2] Cherkauer K A, Lettenmaier D P. Simulation of Spatial Variability in Snow and Frozen Soil. Journal of Geophysical Research, 2003, 108 D22 8858

[3] Ding Y J. Precipitation Conditions for the Development of the Present Glaciers. Geo Journal, 1991, 25(2/3): 243-248.

[4] Gao X, Ye B S, Zhang S Q, . Glacier Runoff Variation and Its Influence on River Runoff during 1961–2006 in the Tarim River Basin, China. Science China Series D: Earth Sciences, 2010, 53: 880-891.

[5] Hewitt K. Tributary Glacier Surges: An Exceptional Concentration at Panmah Glacier, Karakoram Himalaya. Journal of Glaciology, 2007, 53(181): 181-188.

[6] Hock R. Glacier Melt: A Review on Processes and Their Modelling. Progress in Physical Geography, 2005, 29(3): 362-391.

[7] Kargel J S, Fink W, Furfaro R, . Hunza Landslide and Monsoon Flooding in Pakistan Call for International Attention to Transboundary Natural Hazards. American Geophysical Union Fall Meeting 2010. American Geophysical Union, San Francisco, California, USA, Abstract #NH52A-0, 2010

[8] Kloke E J, Oerlemans J. Model Study of the Spatial Distribution of the Energy and Mass Balance of Morteratschgletscher, Switzerland. Journal of Glaciology, 2002, 48(163): 506-518.

[9] Liang X, Lettenmaier D P, Wood E F. A Simple Hydrologically Based Model of Land Surface Water and Energy Fluxes for General Circulation Models. Journal of Geophysical Research, 1994, 99(D7): 14415-14428.

[10] Liang X, Wood E F, Lettenmaier D P. Surface Soil Moisture Parameterization of the VIC-2L Model: Evaluation and Modification. Global and Planetary Change, 1996, 13: 195-206.

[11] Liu G. The Climate of the Batura Glacier and Its Adjacent Areas. Lanzhou Institute of Glaciology and Geocryology, Chinese Academy of Sciences, Professional Papers on the Batura Glacier, Karakorum Mountains, 1980 Beijing: Science Press, 99-110.

[12] Lohmann D, Nolte-Holube R, Raschke E. A Large Scale Horizontal Routing Model to be Coupled to Land Surface Parameterization Schemes. Tellus, 1996, 48A: 708-721.

[13] Lohmann D, Raschke E, Nijssen B, . Regional Scale Hydrology: I. Formulation of the VIC-2L Model Coupled to a Routing Model. Hydrological Sciences, 1998, 43: 131-141.

[14] Shi Y F, Zhang X S. Miller K J. Some Studies of the Batura Glacier in the Karakoram Mountains. International Karakoram Project, 1984 Cambridge: Cambridge University Press, 51-63.

[15] Slater A G, Bohn T J, McCreight J L, . A Multimodel Simulation of Pan-Arctic Hydrology. Journal of Geophysical Research, 2007, 112 G04S45

[16] Su F, Adam J C, Bowling L C, . Streamflow Simulations of the Terrestrial Arctic Domain. Journal of Geophysical Research, 2005, 110 D08112

[17] Su F, Adam J C, Trenberth K E, . Evaluation of Surface Water Fluxes of the Pan-Arctic Land Region with a Land Surface Model and ERA-40 Reanalysis. Journal of Geophysical Research, 2006, 111 D05110

[18] Wake C P. Glaciochemical Investigations as a Tool to Determine the Spatial Variation of Snow Accumulation in the Central Karakoram, Northern Pakistan. Annals of Glaciology, 1989, 13: 279-284.

[19] Wang W, Li K, Gao J. Monitoring Glacial Shrinkage Using Remote Sensing and Site-Observation Method on Southern Slope of Kalik Mountain, Eastern Tian Shan, China. Journal of Earth Science, 2011, 22(4): 503-514.

[20] Wu L H, Li H L, Wang L. Application of a Degree-Day Model for Determination of Mass Balance of Urumqi Glacier No. 1, Eastern Tianshan, China. Journal of Earth Science, 2011, 22(4): 470-481.

[21] Yatagai A, Arakawa O, Kamiguchi K, . A 44-Year Daily Gridded Precipitation Dataset for Asia Based on a Dense Network of Rain Gauges. SOLA, 2009, 5: 137-140.

[22] Zhang S Q, Ye B S, Liu S Y, . A Modified Monthly Degree-Day Model for Evaluating Glacier Runoff Changes in China. Part I: Model Development. Hydrological Processes, 2012, 26(11): 1686-1696.

[23] Zhang S Q, Gao X, Ye B S, . A Modified Monthly Degree-Day Model for Evaluating Glacier Runoff Changes in China. Part II: Application. Hydrological Processes, 2012, 26(11): 1697-1706.

[24] Zhao Q D, Ye B S, Ding Y J, . Hydrological Process of a Typical Catchment in Cold Region: Simulation and Analysis. Journal of Gaciology and Geocryology, 2011, 33(3): 595-605.