Temperature Index Approach Research Articles

Estimation of component-wise runoff contribution using temperature index approach, in a snow- and glacier-fed transboundary Shyok River catchment of the Upper Indus Basin

Estimation of component-wise runoff contribution using temperature index approach, in a snow- and glacier-fed transboundary Shyok River catchment of the Upper Indus Basin

Estimating degree-day factors of snow based on energy flux components

Abstract. Meltwater from mountainous catchments dominated by snow and ice is a valuable source of fresh water in many regions. At mid-latitudes, seasonal snow cover and glaciers act like a natural reservoir by storing precipitation during winter and releasing it in spring and summer. Snowmelt is usually modelled either by energy balance or by temperature-index approaches. The energy balance approach is process-based and more sophisticated but requires extensive input data, while the temperature-index approach uses the degree-day factor (DDF) as a key parameter to estimate melt of snow and ice merely from air temperature. Despite its simplicity, the temperature-index approach has proved to be a powerful tool for simulating the melt process especially in large and data-scarce catchments. The present study attempts to quantify the effects of spatial, temporal, and climatic conditions on the DDF of snow in order to gain a better understanding of which influencing factors are decisive under which conditions. The analysis is based on the individual energy flux components; however, formulas for estimating the DDF are presented to account for situations where observed data are limited. A detailed comparison between field-derived and estimated DDF values yields a fair agreement with bias = 0.14 mm ∘C−1 d−1 and root mean square error (RMSE) = 1.12 mm ∘C−1 d−1. The analysis of the energy balance processes controlling snowmelt indicates that cloud cover and snow albedo under clear sky are the most decisive factors for estimating the DDF of snow. The results of this study further underline that the DDF changes as the melt season progresses and thus also with altitude, since melting conditions arrive later at higher elevations. A brief analysis of the DDF under the influence of climate change shows that the DDFs are expected to decrease when comparing periods of similar degree days, as melt will occur earlier in the year when solar radiation is lower, and albedo is then likely to be higher. Therefore, the DDF cannot be treated as a constant parameter especially when using temperature-index models for forecasting present or predicting future water availability.

Long-term (~40 years) mass balance appraisal and response of the Patsio glacier, in the Great Himalayan region towards climate change

Comprehensive assessment of the long-term mass balance (40 years) has been performed for Patsio glacier located in Bhaga river basin in the Great Himalayan region. The temperature-index approach has been employed for generating the time-series of annual specific mass balance of 25 years (1993–2018) using Snow and Avalanche Study Establishment (SASE) field observatory data in Patsio region. Results based on observational data showed nearly balanced condition until the year 1999 (+0.03 m.w.e/yr), followed by an accelerated rate of mass loss during the period 2000–2007 (−0.17 m.w.e./yr). This was evident due to the annual rise (+0.7°C) in the temperature and drop (−122 mm) in the amount of precipitation during 2000–2007. A slight drop in the mass balance (−0.14 m.w.e/yr) was observed during the past decade (2008–2018) which was well supported by in-situ climate observations and published geodetic results in the region. The altitudinal distribution of modelled mass balance was observed well in correlation with the glaciological measurements for the periods 2010–2011 and 2011–2012. Further, seasonal sensitivity characteristic (SSC) was developed for the glacier illustrating the dependence of specific mass balance on monthly anomalies in temperature and precipitation. The developed SSC was used with bias-corrected temperature and precipitation estimates from recently released ERA5 reanalysis to reconstruct annual specific mass balance during since 1979s.

Development and application of a spatially distributed snowmelt runoff model for limited data condition

Two main approaches are used to calculate snowmelt, namely, energy balance (physically based) and temperature index (conceptual). Energy balance requires comprehensive data to solve physical equations. In this study, a distributed conceptual model, “spatially distributed snowmelt runoff model”, based on modified temperature index approach was developed to estimate snow-related parameters like snow density and degree day factor in a distributed manner over space and time to meet the hydrological demands of data-scarce snow-dominated mountainous catchment. The model needs only daily temperature, precipitation, snow cover, and albedo as basic inputs. It gives pixel-wise output of snow density, snow depth, snow water equivalent, degree day, net radiation, snowmelt depth, rain-induced runoff, snowmelt runoff, actual evaporation, and infiltration, while routing the flow up to the catchment outlet to estimate discharge as time series along with contribution breakup between snowmelt and rain-induced runoff. The developed model was tested for a small-sized watershed Nuranang (58 km2) as well as for a medium-sized watershed Mago basin (844 km2) both located in the Tawang district of Arunachal Pradesh, India, and was successfully calibrated for 2004–2005 and validated for 2008–2009 for Nuranang watershed. In the case of Mago basin, a single year (2007) was used for calibration and validation was also performed for a single year (2009) as well due to the unavailability of data. Performance of the model was found satisfactory for both the basins with modelling efficiency (ME) greater than 0.6 and coefficient of residual mass (CRM) in between − 0.2 and 0.2 for both calibration and validation years in these data-scarce eastern Himalayan watersheds. Considering the satisfactory performance of the model in both the small- and medium-sized watersheds, the model can be very useful for any un-instrumented mountainous basins for pixel-wise snowmelt runoff computation and gridded output generation using degree day approach, where observed snow density data is not available.

Evaluation of SWAT Model performance on glaciated and non-glaciated subbasins of Nam Co Lake, Southern Tibetan Plateau, China

This paper presents an assessment of the Soil and Water Assessment Tool (SWAT) on a glaciated (Qugaqie) and a non-glaciated (Niyaqu) subbasin of the Nam Co Lake. The Nam Co Lake is located in the southern Tibetan Plateau, two subbasins having catchment areas of 59 km2 and 388 km2, respectively. The scores of examined evaluation indices (i.e., R2, NSE, and PBIAS) established that the performance of the SWAT model was better on the monthly scale compared to the daily scale. The respective monthly values of R2, NSE, and PBIAS were 0.94, 0.97, and 0.50 for the calibration period while 0.92, 0.88, and −8.80 for the validation period. Glacier melt contribution in the study domain was simulated by using the SWAT model in conjunction with the Degree Day Melt (DDM) approach. The conjunction of DDM with the SWAT Model ensued improved results during both calibration (R2=0.96, NSE=0.95, and PBIAS=−13.49) and validation (R2=0.97, NSE=0.96, and PBIAS=−2.87) periods on the monthly time scale. Average contribution (in percentage) of water balance components to the total streamflow of Niyaqu and Qugaqie subbasins was evaluated. We found that the major portion (99.45%) of the streamflow in the Niyaqu subbasin was generated by snowmelt or rainfall surface runoff (SURF_Q), followed by groundwater (GW_Q, 0.47%), and lateral (LAT_Q, 0.06%) flows. Conversely, in the Qugaqie subbasin, major contributor to the streamflow (79.63%) was glacier melt (GLC_Q), followed by SURF_Q (20.14%), GW_Q (0.13%), and LAT_Q (0.089%). The contribution of GLC_Q was the highest (86.79%) in July and lowest (69.95%) in September. This study concludes that the performance of the SWAT model in glaciated catchment is weak without considering glacier component in modeling; however, it performs reasonably well in non-glaciated catchment. Furthermore, the temperature index approach with elevation bands is viable in those catchments where streamflows are driven by snowmelt. Therefore, it is recommended to use the SWAT Model in conjunction with DDM or energy base model to simulate the glacier melt contribution to the total streamflow. This study might be helpful in quantification and better management of water resources in data scarce glaciated regions.

Assessing glacier melt contribution to streamflow at Universidad Glacier, central Andes of Chile

Abstract. Glacier melt is an important source of water for high Andean rivers in central Chile, especially in dry years, when it can be an important contributor to flows during late summer and autumn. However, few studies have quantified glacier melt contribution to streamflow in this region. To address this shortcoming, we present an analysis of meteorological conditions and ablation for Universidad Glacier, one of the largest valley glaciers in the central Andes of Chile at the head of the Tinguiririca River, for the 2009–2010 ablation season. We used meteorological measurements from two automatic weather stations installed on the glacier to drive a distributed temperature-index and runoff routing model. The temperature-index model was calibrated at the lower weather station site and showed good agreement with melt estimates from an ablation stake and sonic ranger, and with a physically based energy balance model. Total modelled glacier melt is compared with river flow measurements at three sites located between 0.5 and 50 km downstream. Universidad Glacier shows extremely high melt rates over the ablation season which may exceed 10 m water equivalent in the lower ablation area, representing between 10 and 13 % of the mean monthly streamflow at the outlet of the Tinguiririca River Basin between December 2009 and March 2010. This contribution rises to a monthly maximum of almost 20 % in March 2010, demonstrating the importance of glacier runoff to streamflow, particularly in dry years such as 2009–2010. The temperature-index approach benefits from the availability of on-glacier meteorological data, enabling the calculation of the local hourly variable lapse rate, and is suited to high melt regimes, but would not be easily applicable to glaciers further north in Chile where sublimation is more significant.

The Spatial and Temporal Contribution of Glacier Runoff to Watershed Discharge in the Yarkant River Basin, Northwest China

In this paper, a glacial module based on an enhanced temperature-index approach was successfully introduced into the Soil and Water Assessment Tool (SWAT) model to simulate the glacier runoff and water balance of a glacierized watershed, the mountainous region of the Yarkant River Basin (YRB) in Karakoram. Calibration and validation of the SWAT model were based on comparisons between the simulated and observed discharge with a monthly temporal resolution from 1961 to 2011 for the Kaqun hydrological station. The results reaffirmed the viability of the approach for simulating glacier runoff, as evidenced by a Nash–Sutcliff Efficiency (NSE) of 0.82–0.86 as well as a percentage bias (PBIAS) of −4.5% to 2.4%, for the calibration and validation periods, respectively. Over the last 50 years, the total discharge and glacier runoff both exhibited increasing trends with 0.031 × 109 m3·a−1 and 0.011 × 109 m3·a−1. The annual glacier runoff contribution to the streamflow was between 42.3% and 64.5%, with an average of 51.6%, although the glaciers accounted for only 12.6% of the watershed drainage area in the mountainous YRB. The monthly contribution of the glacier runoff ranged from 11.0% in April to 62.1% in August, and the glacier runoff from June to September accounted for about 86.3% of the annual glacier runoff. Runoff from the mountainous regions above 5000 m a.s.l. accounted for 70.5% of the total discharge, with glacier runoff contributions being approximately 46.4%.

Modification of the degree-day formula for diurnal meltwater generation and refreezing

The standard degree-day, temperature-index approach to calculating snowmelt generation and refreezing (the SDD method) is convenient and popularly used but seriously miscalculates the volumes of water that change phase on days when temperatures fluctuate either side of the freezing point. Additionally, the SDD method does not provide any estimate of the duration of daily melting and refreezing events. A modified version of the standard formula is introduced (the MDD method) that overcomes such problems by removing dependence on a single temperature index (the average daily temperature estimated over a 24-h period beginning at midnight) and instead transfers reliance onto daily air temperature extremes (maximum and minimum temperatures) at known times of occurrence. In this way, the modified formula retains the simplicity of the standard approach while targeting those segments of the diurnal air temperature curve that directly relate to periods of melting and freezing. Newly introduced temperature and time degree-day parameters allow the duration of melting and refreezing events to be estimated. The MDD method was evaluated for two sites in the snow-belt region of Canada where the availability of hourly records of daily temperature allowed the required MDD input parameters to be calculated reliably and thus used for comparative purposes. During testing, the MDD input parameters were obtained from daily temperature extremes and their times of occurrence, using two alternative approaches to synthetic air temperature curve generation, one linear, the other trigonometric. Very good agreement was obtained in both cases and confirms the value of the MDD approach. However, there is no significant benefit to be gained by using air temperature approximating functions more complicated than the linear method for supplementing the missing continuous air temperature measurements. Finally, the MDD approach is not seen as a replacement for the regular SDD method, so much as tool that can be applied when the SDD methodology is likely to become unreliable. This is best achieved by using a hybrid SDD-MDD algorithm that invokes the MDD approach only when the necessary conditions arise.

Assessing the benefit of snow data assimilation for runoff modeling in Alpine catchments

Abstract. In Alpine catchments, snowmelt is often a major contribution to runoff. Therefore, modeling snow processes is important when concerned with flood or drought forecasting, reservoir operation and inland waterway management. In this study, we address the question of how sensitive hydrological models are to the representation of snow cover dynamics and whether the performance of a hydrological model can be enhanced by integrating data from a dedicated external snow monitoring system. As a framework for our tests we have used the hydrological model HBV (Hydrologiska Byråns Vattenbalansavdelning) in the version HBV-light, which has been applied in many hydrological studies and is also in use for operational purposes. While HBV originally follows a temperature-index approach with time-invariant calibrated degree-day factors to represent snowmelt, in this study the HBV model was modified to use snowmelt time series from an external and spatially distributed snow model as model input. The external snow model integrates three-dimensional sequential assimilation of snow monitoring data with a snowmelt model, which is also based on the temperature-index approach but uses a time-variant degree-day factor. The following three variations of this external snow model were applied: (a) the full model with assimilation of observational snow data from a dense monitoring network, (b) the same snow model but with data assimilation switched off and (c) a downgraded version of the same snow model representing snowmelt with a time-invariant degree-day factor. Model runs were conducted for 20 catchments at different elevations within Switzerland for 15 years. Our results show that at low and mid-elevations the performance of the runoff simulations did not vary considerably with the snow model version chosen. At higher elevations, however, best performance in terms of simulated runoff was obtained when using the snowmelt time series from the snow model, which utilized data assimilation. This was especially true for snow-rich years. These findings suggest that with increasing elevation and the correspondingly increased contribution of snowmelt to runoff, the accurate estimation of snow water equivalent (SWE) and snowmelt rates has gained importance.

Backward snow depth reconstruction at high spatial resolution based on time‐lapse photography

AbstractWe report a methodology for reconstructing the daily snow depth distribution at high spatial resolution in a small Pyrenean catchment using time‐lapse photographs and snow depletion rates derived from an on‐site measuring meteorological station. The results were compared with the observed snow depth distribution, determined on a number of separate occasions using a terrestrial laser scanner (TLS). The time‐lapse photographs were projected onto a digital elevation model of the study site, and converted into snow presence/absence information. The melt‐out date (MOD; first occurrence of melt out after peak snow accumulation) was obtained from the projected photograph series. Commencing the backward reconstruction for each grid cell at the MOD, the method uses simulated snow depth depletion rates using a temperature index approach, which are extrapolated to the grid cells of the domain to arrive at the snow distribution of the previous day. Two variants of the reconstruction techniques were applied (1) using a spatially constant degree day factor (DDF) for calculating the daily expected snow depth depletion rate, and (2) allowing a spatially distributed DDF calculated from two consecutive TLS acquisitions compared to the snow depth depletion rate observed at the meteorological station. Validation revealed that both methods performed well (average R2 = 0.68; standard RMSE = 0.58), with better results obtained from the spatially distributed approach. Nevertheless, the spatially corrected DDF reconstruction, which requires TLS data, suggests that the constant DDF approach is an efficient, and for most applications sufficiently accurate and easily reproducible method. The results highlight the usefulness of time‐lapse photography for not only determining the snow covered area, but also for estimating the spatial distribution of snow depth. Copyright © 2016 John Wiley & Sons, Ltd.

A conceptual, distributed snow redistribution model

Abstract. When applying conceptual hydrological models using a temperature index approach for snowmelt to high alpine areas often accumulation of snow during several years can be observed. Some of the reasons why these "snow towers" do not exist in nature are vertical and lateral transport processes. While snow transport models have been developed using grid cell sizes of tens to hundreds of square metres and have been applied in several catchments, no model exists using coarser cell sizes of 1 km2, which is a common resolution for meso- and large-scale hydrologic modelling (hundreds to thousands of square kilometres). In this paper we present an approach that uses only gravity and snow density as a proxy for the age of the snow cover and land-use information to redistribute snow in alpine basins. The results are based on the hydrological modelling of the Austrian Inn Basin in Tyrol, Austria, more specifically the Ötztaler Ache catchment, but the findings hold for other tributaries of the river Inn. This transport model is implemented in the distributed rainfall–runoff model COSERO (Continuous Semi-distributed Runoff). The results of both model concepts with and without consideration of lateral snow redistribution are compared against observed discharge and snow-covered areas derived from MODIS satellite images. By means of the snow redistribution concept, snow accumulation over several years can be prevented and the snow depletion curve compared with MODIS (Moderate Resolution Imaging Spectroradiometer) data could be improved, too. In a 7-year period the standard model would lead to snow accumulation of approximately 2900 mm SWE (snow water equivalent) in high elevated regions whereas the updated version of the model does not show accumulation and does also predict discharge with more accuracy leading to a Kling–Gupta efficiency of 0.93 instead of 0.9. A further improvement can be shown in the comparison of MODIS snow cover data and the calculated depletion curve, where the redistribution model increased the efficiency (R2) from 0.70 to 0.78 (calibration) and from 0.66 to 0.74 (validation).

Bulk volumetric liquid water content in a seasonal snowpack: modeling its dynamics in different climatic conditions

We focus on the dynamics of volumetric liquid water content in seasonal snow covers. This is a key variable describing the fate of snowpacks during the melting season. However, its measurement and/or prediction by means of models at high spatial and temporal resolutions is still difficult due to both practical and theoretical reasons. To overcome these limitations in operational applications, we test the capability of a one-dimensional model to predict the dynamics of bulk volumetric liquid water content during a snow season. Multi-year data collected in three experimental sites in Japan are used as an evaluation. These sites are subjected to different climatic conditions. The model requires the calibration of one or two parameters, according to the degree of detail used. Either a simple temperature-index or a coupled melt–freeze temperature-index approach are considered to predict melting and/or melt–freeze dynamics of liquid water. Results show that, if melt–freeze dynamics are modeled, median absolute differences between data and predictions are consistently lower than 1 vol% at the sites where data of liquid water content are available. In addition, we find also that the model predicts correctly a dry condition in 80% of the observed cases at a site where calibration data are scarce. At the same site, observed isothermal conditions of the snow cover at 0 °C correspond to predictions of bulk volumetric liquid water content that are greater than 0.

Modified temperature index model for estimating the melt water discharge from debris-covered Lirung Glacier, Nepal

Abstract. In the Nepalese Himalayas, the complex topography, occurrence of debris covered glaciers, and limited data availability creates substantial difficulties for modelling glacier melt. The proper recognition of melt processes governs the accurate estimation of melt water from glacier dominated systems, even in the presence of debris-covered glaciers. This paper presents a glacier melt model developed for the Lirung sub-basin of Langtang valley, which has both a clean glacier area, 5.86 km2, and a debris-covered glacier area, 1.13 km2. We use a temperature index approach to estimate sub-daily melt water discharge for a two week period at the end of monsoon, and the melt factor is varied according to the aspect and distributed to each grid processed from the digital elevation model. The model uses easily available data and simple extrapolation techniques capable of generating melt with limited data. The result obtained from this method provides accurate estimate with an R2 value of 0.89, bias of 0.9% and Nash-Sutcliffe efficiency of 0.86, and suitable in Himalaya where data availability is major issue.

Effect of meteorological forcing and snow model complexity on hydrological simulations in the Sieber catchment (Harz Mountains, Germany)

Abstract. Detailed physically based snow models using energy balance approaches are spatially and temporally transferable and hence regarded as particularly suited for scenario applications including changing climate or land use. However, these snow models place high demands on meteorological input data at the model scale. Besides precipitation and temperature, time series of humidity, wind speed, and radiation have to be provided. In many catchments these time series are rarely available or provided by a few meteorological stations only. This study analyzes the effect of improved meteorological input on the results of four snow models with different complexity for the Sieber catchment (44.4 km2) in the Harz Mountains, Germany. The Weather Research and Forecast model (WRF) is applied to derive spatial and temporal fields of meteorological surface variables at hourly temporal resolution for a regular grid of 1.1 km × 1.1 km. All snow models are evaluated at the point and the catchment scale. For catchment-scale simulations, all snow models were integrated into the hydrological modeling system PANTA RHEI. The model results achieved with a simple temperature-index model using observed precipitation and temperature time series as input are compared to those achieved with WRF input. Due to a mismatch between modeled and observed precipitation, the observed melt runoff as provided by a snow lysimeter and the observed streamflow are better reproduced by application of observed meteorological input data. In total, precipitation is simulated statistically reasonably at the seasonal scale but some single precipitation events are not captured by the WRF data set. Regarding the model efficiencies achieved for all simulations using WRF data, energy balance approaches generally perform similarly compared to the temperature-index approach and partially outperform the latter.

Comparison of approaches for snowpack estimation in New York City watersheds

AbstractSnow is a substantial component of historical annual precipitation in New York City (NYC) water supply watersheds in the Catskill Mountains, and the pattern of snow accumulation and snowmelt has important implications for the management of the reservoirs and watersheds that are part of the NYC water supply. NYC currently estimates reservoir basin‐scale snowpack throughout the snow season by extrapolation from biweekly snow survey data. These estimates are complemented by the NOAA Snow Data Assimilation System (SNODAS) product. Snowpack models are used in short‐term projections to support reservoir operations and long‐term simulations to evaluate the potential effects of climate change, land use change, and watershed management on the water supply. We tested three snowpack estimation approaches compared with snow survey data: the lumped‐parameter temperature index approach from the Generalized Watershed Loading Function (GWLF) watershed model, a spatially distributed temperature index (SDTI) model, and the spatially distributed NOAA SNODAS product. Of the spatially distributed approaches, SNODAS estimated the spatial variability of snow water equivalent (SWE) among snow survey sites within a basin better than the SDTI model. All three snowpack estimation approaches, including the lumped‐parameter GWLF model, performed well in estimating basin‐wide SWE for most of the basins studied. Copyright © 2013 John Wiley & Sons, Ltd.

Performance of complex snow cover descriptions in a distributed hydrological model system: A case study for the high Alpine terrain of the Berchtesgaden Alps

[1] Runoff generation in Alpine regions is typically affected by snow processes. Snow accumulation, storage, redistribution, and ablation control the availability of water. In this study, several robust parameterizations describing snow processes in Alpine environments were implemented in a fully distributed, physically based hydrological model. Snow cover development is simulated using different methods from a simple temperature index approach, followed by an energy balance scheme, to additionally accounting for gravitational and wind-driven lateral snow redistribution. Test site for the study is the Berchtesgaden National Park (Bavarian Alps, Germany) which is characterized by extreme topography and climate conditions. The performance of the model system in reproducing snow cover dynamics and resulting discharge generation is analyzed and validated via measurements of snow water equivalent and snow depth, satellite-based remote sensing data, and runoff gauge data. Model efficiency (the Nash-Sutcliffe coefficient) for simulated runoff increases from 0.57 to 0.68 in a high Alpine headwater catchment and from 0.62 to 0.64 in total with increasing snow model complexity. In particular, the results show that the introduction of the energy balance scheme reproduces daily fluctuations in the snowmelt rates that trace down to the channel stream. These daily cycles measured in snowmelt and resulting runoff rates could not be reproduced by using the temperature index approach. In addition, accounting for lateral snow transport changes the seasonal distribution of modeled snowmelt amounts, which leads to a higher accuracy in modeling runoff characteristics.

Glacier mass balance of Norway 1961-2010 calculated by a temperature-index model

AbstractGlacier mass balance in Norway is only observed over a small portion (<15%) of the glacierized surface and only for short time periods (<10 years) for most sites. To provide a comprehensive overview of the temporal mass-balance evolution, we modeled surface mass balance for the glacierized area of mainland Norway from 1961 to 2010. The model is forced by operationally gridded daily temperature and precipitation fields which are available at 1 km horizontal resolution from 1957 until the present. The applied mass-balance model accounts for melting of snow and ice by using a distributed temperature-index approach. The precipitation input is corrected to obtain agreement between modeled and observed winter mass balance, and a melt factor and two radiation coefficients are optimized to the corresponding summer balance. The model results show positive trends of winter balance between 1961 and 2000 followed by a remarkable decrease in both summer and winter balances which resulted in an average annual balance of –0.86 ± 0.15 m w.e. a-1 between 2000 and 2010 after four decades of zero to slightly positive annual mass balances.

Snow cover variation and streamflow simulation in a snow-fed river basin of the Northwest Himalaya

Snowmelt is an important component of any snow-fed river system. The Jhelum River is one such transnational mountain river flowing through India and Pakistan. The basin is minimally glacierized and its discharge is largely governed by seasonal snow cover and snowmelt. Therefore, accurate estimation of seasonal snow cover dynamics and snowmelt-induced runoff is important for sustainable water resource management in the region. The present study looks into spatio-temporal variations of snow cover for past decade and stream flow simulation in the Jhelum River basin. Snow cover extent (SCE) was estimated using MODIS (Moderate Resolution Imaging Spectrometer) sensor imageries. Normalized Difference Snow Index (NDSI) algorithm was used to generate multi-temporal time series snow cover maps. The results indicate large variation in snow cover distribution pattern and decreasing trend in different sub-basins of the Jhelum River. The relationship between SCE-temperature, SCE-discharge and discharge-precipitation was analyzed for different seasons and shows strong correlation. For streamflow simulation of the entire Jhelum basin Snow melt Runoff Model (SRM) used. A good correlation was observed between simulated stream flow and in-situ discharge. The monthly discharge contribution from different sub-basins to the total discharge of the Jhelum River was estimated using a modified version of runoff model based on temperature-index approach developed for small watersheds. Stream power — an indicator of the erosive capability of streams was also calculated for different sub-basins.

Integrating glaciers raster-based modelling in large catchments hydrological balance: the Rhone case study

A raster-based glacier sub-model was successfully introduced in the distributed hydrological model FEST-WB to simulate the water balance and surface runoff of large Alpine catchments. The glacier model is based on temperature-index approach for melt, on linear reservoir for melt water propagation into the ice and on mass balance for accumulation; the initialization of the volume of ice on the basin was based on a formulation depending on surface topography. The model was first tested on a sub-basin of the Rhone basin (Switzerland), which is for 62% glaciated; the calibration and validation were based on comparison between simulated and observed discharge from 1999 to 2008. The model proved to be suitable to simulate the typical discharge seasonality of a heavily glaciated basin. The performance of the model was also tested by simulating discharge in the whole Swiss Rhone basin, in which glaciers contribution is not negligible, in fact, in summer, about the 40% of the discharge is due to glacier melt. The model allowed to take into account the volume of water coming from glaciers melt and its simple structure is suitable for analysis of the effects of climate change on hydrological regime of high mountain basins, with available meteorological forcing from current RCM. Copyright © 2012 John Wiley & Sons, Ltd.

Streamflow Modeling in a Highly Managed Mountainous Glacier Watershed Using SWAT: The Upper Rhone River Watershed Case in Switzerland

Streamflow simulation is often challenging in mountainous watersheds because of irregular topography and complex hydrological processes. Rates of change in precipitation and temperature with respect to elevation often limit the ability to reproduce stream runoff by hydrological models. Anthropogenic influence, such as water transfers in high altitude hydro- power reservoirs increases the difficulty in modeling since the natural flow regime is altered by long term storage of water in the reservoirs. The Soil and Water Assessment Tool (SWAT) was used for simulating streamflow in the upper Rhone watershed located in the south western part of Switzerland. The catchment area covers 5220 km 2 , where most of the land cover is dominated by forest and 14 % is glacier. Streamflow calibration was done at daily time steps for theperiodof2001-2005,andvalidated for2006-2010.Twodifferent approaches wereused for simulating snow and glacier melt process, namely the temperature index approach with and without elevation bands. The hydropower network was implemented based on theintake points that form part of the inter-reservoir network. Subbasins were grouped into two major categories