An operational hydrological forecast model was set-up based on the Snowmelt-Runoff Model (SRM) in order to forecast Kharif flows from Upper Jhelum catchment. Zone-wise degree-day factor functions were derived by diagnostic calibration and are applied according to a defined temperature rule when melting starts. While predicting the depletion of snow-covered area by SRM's modified depletion curve approach, scenario runs with temperature and precipitation of past years are carried out which are evaluated statistically to forecast the seasonal flow volume.
Irrigated agriculture provides 90 % of Pakistan's food requirements,
22 % of its Gross Domestic Product and provides employment to 60 % of
the population of the country (Euroconsult, 2011). Pakistan ranks 5th in
rice export and 8th in wheat production in the world (FAO, 2011).
Consequently, agriculture's share in water usage is about 97 %, which is
well above the global average of about 70 % (Akram, 2009). The Indus
Basin Irrigation System, comprising a canal command area of more than
14 million ha, is one of the largest contiguous irrigation systems of the
world. On average, annual canal water diversion from River Indus and its
major tributaries amounts to about 130 km
Flows in the Indus Basin mainly originate from the high mountain ranges of the Western Himalaya – Karakoram – Hindu Kush region. Though the flow regime also comprises glacier melt and runoff from rainfall during monsoon season, the predominant contribution comes from the melting of seasonal snow accumulated during the preceding winter and spring (Archer et al., 2010). Depending on the altitude of the catchment, flows originating from snowmelt start to rise in March to April. Average monthly flow volumes reach their peak in May to July.
Expected flows during the forthcoming Kharif cropping season (April–September) are forecasted by end of March. Based on these forecasts, the Indus River System Authority (IRSA) decides the provincial shares and the provincial irrigation departments subsequently determine the seasonal water allocation to the different canal command areas. A reliable seasonal forecast of the water resources to be expected from snowmelt is therefore of utmost importance for the agricultural production.
One of only two major reservoirs in Pakistan, the Mangla reservoir, is
located at the foothills of the Western Himalayas on Jhelum river (Fig. 1).
After a dam raising project completed in 2009, the storage capacity is 9.1 km
Landsat TM scene of study area.
The catchment of Jhelum River upstream of Mangla dam has an area of about
33 500 km
Snowmelt is a complex physical process that is primarily driven by the energy flux entering the snowpack. Two categories of snowmelt models exists: energy balance models attempt to quantify the components of the heat balance equation while temperature index models use air temperature as a principal predictor of melt rates. Temperature index models have been widely used as they offer generally a good model performance while data requirements are comparatively low which is of particular importance in remote and rugged mountainous catchments like in the Himalaya – Karakoram – Hindu Kush region. Among the class of temperature index models SRM (Martinec, 1975) has become very popular in conjunction with the increasing availability of remote sensing snow cover data, as it circumvents the error prone snow accumulation approach used by other models.
An Excel® version (Bogacki and Hashmi, 2013) of WinSRM (Martinec et al., 2011) was developed in order to allow for changes in the underlying code like the handling of heavy rainfalls that had to be adapted to the catchment characteristics. Furthermore, Excel® is well known to most engineers who can easily adapt the layout or data structure to their specific needs.
Distribution of monthly inflow to Mangla reservoir 2000–2014.
SRM uses air temperature, snow covered area, and precipitation as daily input variables. The watershed is divided into elevation zones using their mean hypsometric elevation as reference for zone-wise degree-day calculation.
There are three high elevation climate stations in the Pakistani part of the
Upper Jhelum catchment operated by WAPDA's Pakistan Water and
Power Development Authority Global Summary Of the Day. Download at:
The MODIS (Terra) Snow Cover Daily product Moderate-resolution
imaging spectro-radiometer
As spatial interpolation of daily precipitation station data in mountainous
regions is particularly difficult, the remote sensing based RFE 2.0 South
Asia RainFall Estimates version 2.0 created by the NOAA Climate
Prediction Center's FEWS-NET group (2011–2015) sponsored by USAID. Download at:
The most important model parameter controlling daily snowmelt is the
degree-day factor [cm
The degree-day factor generally increases during melt season as the snowpack becomes “ripe” due to accumulation of energy by the time (e.g. Martinec et al., 2011; Hock, 2003). This process happens later at higher elevation zones as temperatures are lower but the increase is quicker than in the lower zones, as energy input by solar radiation is longer in effect until actual melting starts. In order to develop a generalised rule as needed in the forecasting procedure, zone-wise degree-day factor functions (Ismail et al., 2015) where developed. The relation of degree-day factors versus time, i.e. the 10-days periods, was plotted for each elevation zone for all calibrated years, using the first period when the increase commences as starting point. Figure 3 shows the evolution of degree-day factors by time for elevation zones 6 (2500–3000 m a.s.l.) and 9 (4000–4500 m a.s.l.) respectively. Finally a degree-day factor function was developed by linear regression for each elevation zone of Upper Jhelum catchment (Table 1).
Increase of degree-day factors with time (10-days periods) after melting start for elevation zones 6 and 9. Degree-day factors are obtained by diagnostic calibration.
The start of melting in each zone differs from year to year, depending on
the actual temperature development during spring. Based on average 10-days
period temperatures, a forecast rule was developed controlling when to start
the respective degree-day factor function in a particular zone. These
threshold temperature values decrease from 6
Zone-wise degree-day factors depending on 10-days periods after melting start.
Comparison of annual and Kharif (April–September) mean absolute volumetric error [%] of 10-daily flow forecasts using (a) scenario approach, (b) NCEP, and (c) ECMWF 15-days weather forecasts.
The other model parameters required by SRM like temperature lapse-rate, recession coefficient, runoff coefficient for snow, lag-time, etc., were applied basin-wide and kept constant for all years. The values of these parameters were determined according to the methods described by Martinec et al. (2011) and slightly adjusted to achieve a good fit over the whole calibration period. Only the procedure of adjusting the recession coefficient for heavy rainfalls, which is hard coded in WinSRM had to be adapted to the catchment characteristics, as otherwise peak runoffs were heavily overestimated.
The main objective of the Upper Jhelum snowmelt runoff model is the
operational forecast of water availability in the Kharif cropping season
(April–September) but the model is also utilised for 10-daily forecasts
in order to fine-tune irrigation water distribution during the season.
Because medium-term weather forecasts were deemed more promising than
long-term seasonal, in a first analysis 10-daily flow hindcasts were carried
out using 15-days temperature and precipitation forecasts from
NCEP National Centers for Environmental Prediction, USA
( European
Centre for Medium-Range Weather Forecasts ( TIGGE (
When forecasting the three model variables snow-covered area, temperature, and precipitation have to be predicted. In order to estimate the future depletion of the snow-covered area SRM uses so called “modified depletion curves” which are derived from the conventional depletion curves of each elevation zone by replacing the time scale with the cumulative daily snow-melt depth (Martinec et al., 2011). The decline of the modified depletion curves depends on the initial accumulation of snow and represents the actual snow-water equivalent. When initial snow depth is low the modified depletion curve declines faster than in years when a lot of snow has accumulated.
Results of validation of final flow forecast model (dashed line)
compared to observed inflows to Mangla reservoir (solid line) for the year
2011. The graph also shows at the bottom rainfall and snowmelt depth [cm d
End of March, when the seasonal forecast is carried out, an elevation zone showing already some decline in snow-covered area, and hence having also some cumulated degree-days, is chosen as “key zone”. Comparing the relation of decline in snow-covered area versus cumulated degree-days with a statistical analysis of the modified depletion curves of previous years the actual amount of snow is estimated and the future depletion anticipated accordingly, while assuming similar snow conditions for all elevation zones.
While the snow-covered area and its depletion is calculated only once for the actual season to be forecasted, scenario runs are carried out with historic temperature and precipitation data-sets of various years. This results in an ensemble of predicted seasonal flows representing historic weather conditions from which by statistical analysis a forecast of “most likely” (median) as well as expected flows under “dry” or “wet” conditions can be derived.
The final flow forecast model, i.e. using fixed model parameters as well as
the developed degree-day factor functions and the start rule as described in
Sect. 3.4, was validated for the series of years 2002–2011. The relative
volume difference
The years 2002 and 2011 that had not been used for parameter calibration and
determination of the degree-day factor functions resulted in an
The forecast capability of the developed model was evaluated by hindcasts
for the years 2000–2011 using the same procedures as developed for the
operational seasonal forecasts. Meanwhile, real forecasts exist for the
Kharif seasons of the years 2012–2015. Results are compared with IRSA's
forecasts that are based on a statistical approach and with forecasts from
the UBC University of British Columbia Watershed Model
Statistics of relative error of Kharif flow volume forecasts by IRSA's statistical approach, UBC watershed model, and SRM. SRM 2000–2011 are hindcasts, all other forecasts. The series include the approx. 200 years flood in 2010.
The quite good accuracy of seasonal flow prediction in the Upper Jhelum catchment might be due to a combination of favourable factors. One is definitely SRM's good approximation of the actual snow situation at the start of the melting season by utilising the MODIS snow-cover remote-sensing data in conjunction with the modified depletion curve approach. Once the amount of available snow is known, it is most likely, that all snow is melted during the Kharif season.
In addition, the variability of the second important flow component, i.e. rainfall from westerly disturbances in spring and monsoon events in summer, is sufficiently considered by the scenario approach, which transforms the inter-annual variation of precipitation into a series of seasonal flows that can be statistically evaluated. The use of the gridded precipitation product RFE, although it sometimes over- or underestimates a single rainfall event considerably, might also be favourable in that context as the spatial distribution seems to be more fitting as any interpolation of station data in such mountainous regions.
In order to further improve the seasonal forecast capability of the scenario
approach, research is under way on the early identification of cold/warm
respectively dry/wet years e.g. correlated to the ENSO El Niño
Southern Oscillation
Considering that the error in predicting the Kharif flow volume is less than 10 % in 14 out of 16 years with a mean absolute error in this period of 6.4 % it can be stated that the combination of SRM with remote sensing data and the scenario approach has proven to be a reliable procedure for operational seasonal flow forecasting in the Upper Jhelum catchment. Furthermore this approach could be easily improved if the meteorological characteristics of the forthcoming Kharif season can be teleconnected to global climatic conditions like the ENSO status during the winter.
The study was carried out as part of the project “Upgrading of Tools, Water Resources Database, Management Systems and Models” under sub component “B1” of World Bank financed “PK Water Sector Capacity Building & Advisory Services Project” (WCAP): P110099. The authors thank National Engineering Services Pakistan (Pvt.) Ltd. (NESPAK), Lahore and AHT Group AG, Essen, Germany that they could be part of the project team. They are highly grateful to civil engineering department of University of Engineering and Technology Lahore and civil engineering department of Hochschule Koblenz, University of Applied Sciences, for supporting subsequent research and wish also to express their high gratitude to Indus River System Authority (IRSA) and WAPDA's Snow and Ice Hydrology Project for sharing their forecast results.