Abstract. We estimated the black carbon (BC) concentration over the Hindu Kush Himalayan region (HKH), its impact on snow albedo reduction, and sensitivity on annual glacier runoff over the identified glaciers. These estimates were based on free-running aerosol simulations (freesimu) and constrained aerosol simulations (constrsimu) from an atmospheric general circulation model, combined with numerical simulations of a glacial mass balance model. BC concentration estimated from freesimu performed better over higher altitude (HA) HKH stations than that over lower altitude (LA) stations. The estimates from constrsimu mirrored the measurements well when implemented for LA stations. Estimates of the spatial distribution of BC concentration in the snowpack (BCc) over the HKH region led to identifying a hot-spot zone located around Manora Peak. Among glaciers over this zone, BCc (>60 µg kg−1) and BC-induced snow albedo reduction (≈5 %) were estimated explicitly being high during the pre-monsoon for Pindari, Poting, Chorabari, and Gangotri glaciers (which are major sources of fresh water for the Indian subcontinent). The rate of increase of BCc in recent years (i.e., over the period 1961–2010) was, however, estimated to be the highest for the Zemu Glacier. Sensitivity analysis with a glacial mass balance model indicated the increase in annual runoff from debris-free glacier areas due to BC-induced snow albedo reduction (SAR) corresponding to the BCc estimated for the HKH glaciers was 4 %–18 %, with the highest being for the Milam and Pindari glaciers. The rate of increase in annual glacier runoff per unit BC-induced percentage SAR was specifically high for Milam, Pindari, and Sankalpa glaciers. The source-specific contribution to atmospheric BC aerosols by emission sources led to identifying the potential emission source being primarily from the biofuel combustion in the Indo-Gangetic Plain south of 30∘ N, but also from open burning in a more remote region north of 30∘ N.

Accurate long-term inflow forecasts are essential for optimal planning of hydropower production. In snow-rich regions, where spring snowmelt is often the largest reservoir of water, inflow forecasts may be improved by assimilating snow observations to achieve more accurate initial states for the hydrological models prior to the prognosis. In this study, we test whether an ensemble Kalman based approach is useful for this purpose for a mountainous catchment in Norway. For 15 years, annual snow observations near peak accumulation at three locations were assimilated into a distributed hydrological model. After the update, the model was run for a 4-month forecasting period with inflows compared to a base case scenario that omitted the snow observations. The assimilation framework improved the forecasts in several years, and in two of the years, the improvement was very large compared to the base case simulation. At the same time, the filter did not degrade the forecasts largely, indicating that though the updating might slightly degrade performance in some years, it maintains the potential for large improvements in others. Thus, the framework proposed here is a viable method for improving snow-related deficiencies in the initial states, which translates to better forecasts.

Abstract. Light absorbing impurities in snow and ice (LAISI) originating from atmospheric deposition enhance snowmelt by increasing the absorption of shortwave radiation. The consequences are a shortening of the snow duration due to increased snowmelt and, at the catchment scale, a temporal shift in the discharge generation during the spring melt season.In this study, we present a newly developed snow algorithm for application in hydrological models that allows for an additional class of input variable: the deposition mass flux of various species of light absorbing aerosols. To show the sensitivity of different model parameters, we first use the model as a 1-D point model forced with representative synthetic data and investigate the impact of parameters and variables specific to the algorithm determining the effect of LAISI. We then demonstrate the significance of the radiative forcing by simulating the effect of black carbon (BC) deposited on snow of a remote southern Norwegian catchment over a 6-year period, from September 2006 to August 2012. Our simulations suggest a significant impact of BC in snow on the hydrological cycle. Results show an average increase in discharge of 2.5, 9.9, and 21.4 %, depending on the applied model scenario, over a 2-month period during the spring melt season compared to simulations where radiative forcing from LAISI is not considered. The increase in discharge is followed by a decrease in discharge due to a faster decrease in the catchment's snow-covered fraction and a trend towards earlier melt in the scenarios where radiative forcing from LAISI is applied. Using a reasonable estimate of critical model parameters, the model simulates realistic BC mixing ratios in surface snow with a strong annual cycle, showing increasing surface BC mixing ratios during spring melt as a consequence of melt amplification. However, we further identify large uncertainties in the representation of the surface BC mixing ratio during snowmelt and the subsequent consequences for the snowpack evolution.

Discharge over the Narayani river catchment of Nepal was simulated using Statkraft's Hydrologic Forecasting Toolbox (Shyft) forced with observations and three global forcing datasets: (i) ERA-Interim (ERA-I), (ii) Water and Global Change (WATCH) Forcing Data ERA-I (WFDEI), and (iii) Coordinated Regional Climate Downscaling Experiment with the contributing institute Rossy Centre, Swedish Meteorological and Hydrological Institute (CORDEX-SMHI). Not only does this provide an opportunity to evaluate discharge variability and uncertainty resulting from different forcing data but also it demonstrates the capability and potential of using these global datasets in data-sparse regions. The fidelity of discharge simulation is the greatest when using observations combined with the WFDEI forcing dataset (hybrid datasets). These results demonstrate the successful application of global forcing datasets for regional catchment-scale modeling in remote regions. The results were also promising to provide insight of the interannual variability in discharge. This study showed that while large biases in precipitation can be reduced by applying a precipitation correction factor (p_corr_factor), the best result is obtained using bias-corrected forcing data as input, i.e. the WFDEI outperformed other forcing datasets. Accordingly, the WFDEI forcing dataset holds great potential for improving our understanding of the hydrology of data-sparse Himalayan regions and providing the potential for prediction. The use of CORDEX-SMHI- and ERA-I-derived data requires further validation and bias correction, particularly over the high mountain regions.