Rajesh Kumar

Four years of mass balance on Chhota Shigri Glacier, Himachal Pradesh, India, a new benchmark glacier in the western Himalaya

Little is known about the Himalayan glaciers, although they are of particular interest in terms of future water supply, regional climate change and sea-level rise. In 2002, a long-term monitoring programme was started on Chhota Shigri Glacier (32.28 N, 77.58 E; 15.7 km 2 , 6263-4050 m a.s.l., 9 km long) located in Lahaul and Spiti Valley, Himachal Pradesh, India. This glacier lies in the monsoon-arid transition zone (western Himalaya) which is alternately influenced by Asian monsoon in summer and the mid-latitude westerlies in winter. Here we present the results of a 4 year study of mass balance and surface velocity. Overall specific mass balances are mostly negative during the study period and vary from a minimum value of -1.4 m w.e. in 2002/03 and 2005/06 (equilibrium-line altitude (ELA) � 5180 m a.s.l.) to a maximum value of +0.1 m w.e. in 2004/05 (ELA 4855 m a.s.l.). Chhota Shigri Glacier seems similar to mid-latitude glaciers, with an ablation season limited to the summer months and a mean vertical gradient of mass balance in the ablation zone (debris-free part) of 0.7 m w.e. (100 m) -1 , similar to those reported in the Alps. Mass balance is strongly dependent on debris cover, exposure and the shading effect of surrounding steep slopes.

Applicability of Landsat 8 data for characterizing glacier facies and supraglacial debris

Abstract The present work evaluates the applicability of operational land imager (OLI) and thermal infrared sensor (TIRS) on-board Landsat 8 satellite. We demonstrate an algorithm for automated mapping of glacier facies and supraglacial debris using data collected in blue, near infrared (NIR), short wave infrared (SWIR) and thermal infrared (TIR) bands. The reflectance properties in visible and NIR regions of OLI for various glacier facies are in contrast with those in SWIR region. Based on the premise that different surface types (snow, ice and debris) of a glacier should show distinct thermal regimes, the ‘at-satellite brightness temperature’ obtained using TIRS was used as a base layer for developing the algorithm. This base layer was enhanced and modified using contrasting reflectance properties of OLI bands. In addition to facies and debris cover characterization, another interesting outcome of this algorithm was extraction of crevasses on the glacier surface which were distinctly visible in output and classified images. The validity of this algorithm was checked using field data along a transect of the glacier acquired during the satellite pass over the study area. With slight scene-dependent threshold adjustments, this work can be replicated for mapping glacier facies and supraglacial debris in any alpine valley glacier.

A lake detection algorithm (LDA) using Landsat 8 data: A comparative approach in glacial environment

Abstract Glacial lakes show a wide range of turbidity. Owing to this, the normalized difference water indices (NDWIs) as proposed by many researchers, do not give appropriate results in case of glacial lakes. In addition, the sub-pixel proportion of water and use of different optical band combinations are also reported to produce varying results. In the wake of the changing climate and increasing GLOFs (glacial lake outburst floods), there is a need to utilize wide optical and thermal capabilities of Landsat 8 data for the automated detection of glacial lakes. In the present study, the optical and thermal bandwidths of Landsat 8 data were explored along with the terrain slope parameter derived from Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model Version2 (ASTER GDEM V2), for detecting and mapping glacial lakes. The validation of the algorithm was performed using manually digitized and subsequently field corrected lake boundaries. The pre-existing NDWIs were also evaluated to determine the supremacy and the stability of the proposed algorithm for glacial lake detection. Two new parameters, LDI (lake detection index) and LF (lake fraction) were proposed to comment on the performances of the indices. The lake detection algorithm (LDA) performed best in case of both, mixed lake pixels and pure lake pixels with no false detections (LDIxa0=xa00.98) and very less areal underestimation (LFxa0=xa00.73). The coefficient of determination (R2) between areal extents of lake pixels, extracted using the LDA and the actual lake area, was very high (0.99). With understanding of the terrain conditions and slight threshold adjustments, this work can be replicated for any mountainous region of the world.

Remote sensing flow velocity of debris-covered glaciers using Landsat 8 data

Changes in ice velocity of a glacier regulate its mass balance and dynamics. The estimation of glacier flow velocity is therefore an important aspect of temporal glacier monitoring. The utilisation of conventional ground-based techniques for detecting glacier surface flow velocity in the rugged and alpine Himalayan terrain is extremely difficult. Remote sensing-based techniques can provide such observations on a regular basis for a large geographical area. Obtaining freely available high quality remote sensing data for the Himalayan regions is challenging. In the present work, we adopted a differential band composite approach, for the first time, in order to estimate glacier surface velocity for non-debris and supraglacial debris covered areas of a glacier, separately. We employed various bandwidths of the Landsat 8 data for velocity estimation using the COSI-Corr (co-registration of optically sensed images and correlation) tool. We performed the accuracy assessment with respect to field measurements for two glaciers in the Indian Himalaya. The panchromatic band worked best for non-debris parts of the glaciers while band 6 (SWIR – short wave infrared) performed best in case of debris cover. We correlated six temporal Landsat 8 scenes in order to ensure the performance of the proposed algorithm on monthly as well as yearly timescales. We identified sources of error and generated a final velocity map along with the flow lines. Over- and underestimates of the yearly glacier velocity were found to be more in the case of slow moving areas with annual displacements less than 5 m. Landsat 8 has great capabilities for such velocity estimation work for a large geographic extent because of its global coverage, improved spectral and radiometric resolutions, free availability and considerable revisit time.

Automated detection and temporal monitoring of crevasses using remote sensing and their implications for glacier dynamics

Abstract Detailed studies on temporal changes of crevasses and their linkage with glacier dynamics are scarce in the Himalayan context. Observations of temporally changing surficial crevasse patterns and their orientations are suggestive of the processes that determine seasonal glacier flow characteristics. In the present study, on a Himalayan valley glacier, changing crevasse patterns and orientations were detected and mapped on Landsat 8 images in an automated procedure using the ratio of Thermal Infrared Sensor (TIRS) band 10 to Optical Land Imager (OLI) shortwave infrared (SWIR) band 6. The ratio was capable of mapping even crevasses falling under mountain shadows. Differential GPS observations suggested an average error of 3.65% and root-mean-square error of 6.32m in crevasse lengths. A year-round observation of these crevasses, coupled with field-based surface velocity measurements, provided some interesting interpretations of seasonal glacier dynamics.

Development of a Glacio-hydrological Model for Discharge and Mass Balance Reconstruction

The reconstruction of glacio-hydrological records for the data deficient Himalayan catchments is needed in order to study the past and future water availability. The study provides outcomes of a glacio-hydrological model based on the degree-day approach. The model simulates the discharge and mass balance for glacierised Shaune Garang catchment. The degree-day factors for different land covers, used in the model, were estimated using daily stake measurements on Shaune Garang glacier and they were found to be varying between 2.6u2009±u20090.4 and 9.3u2009±u20090.3xa0mm °C−1day−1. The model is validated using observed discharge during ablation season of 2014 with coefficient of determination (R2) 0.90 and root mean square error (RMSE) 1.05 m3 sec−1. The model is used to simulate discharge from 1985 to 2008 and mass balance from 2001 to 2008. The model results show significant contribution of seasonal snow and ice melt in total discharge of the catchment, especially during summer. We observe the maximum discharge in July having maximum contribution from snow and ice melt. The annual melt season discharge shows following a decreasing trend in the simulation period. The reconstructed mass balance shows mass loss of 0.89xa0m we per year between 2001 and 2008 with slight mass gain during 2000/01 and 2004/05 hydrological years.

Dynamics of suspended sediment load with respect to summer discharge and temperatures in Shaune Garang glacierized catchment, Western Himalaya

The observed and predicted rise in temperature will have deleterious impact on melting of snow and ice and form of precipitation which is already evident in Indian Himalayan Region. The temperature-dependent entities like discharge and sediment load will also vary with the observed and predicted rise posing environmental, social and economic threat in the region. There is little known about sediment load transport in relation to temperature and discharge in glacierized catchments in Himalaya mainly due to the scarcity of ground-based observation. The present study is an attempt to understand the suspended sediment load and transportation in relation to variation in discharge and temperature in the Shaune Garang catchment. The result shows strong dependence of sediment concentration primarily on discharge (R2u2009=u20090.84) and then on temperature (R2u2009=u20090.79). The catchments with similar geological and climate setting were observed to have comparatively close weathering rate. The sediment load was found to be higher in the catchments in eastern and central part of Indian Himalayan Region in comparison with western part due to dominance of Indian Summer Monsoon leading to high discharge. The annual physical weathering rate in Shaune Garang catchment was found to be 411xa0txa0km−2xa0year−1 which has increased from 327xa0txa0km−2xa0year−1 in around three decades due to rise in temperature causing increase in discharge and proportion of debris-covered glacierized area.

Ambient black carbon, PM 2.5 and PM 10 at Patna: Influence of anthropogenic emissions and brick kilns

Particulate Matters like Black Carbon, PM2.5 and PM10 present in the atmosphere not only poses a threat to human health but also contributes to near-term regional and global atmospheric warming. There has been a large concern of this kind of pollutants in Indo-Gangetic Plains (IGP). Hence, an attempt has been made to see the impact in one of the highly developing city of IGP as Patna. This city has been ranked among the top 100 air polluted cities in the world. To establish the role of BC, PM2.5, and PM10 on local air quality, continuous measurements were conducted at seven locations of Patna from January to December 2015. The seasonal mass concentration of BC were 13.92±3.48μgm-3 in the winter, 9.65±3.0μgm-3 in the pre-monsoon, 5.83±1.90μgm-3 in the monsoon and 7.86±3.66μgm-3 in the post-monsoon. Similarly, the seasonal concentrations of PM2.5 (PM10) were 68.86±18.83μgm-3 (108.13±21.49μgm-3) in the winter; 64.62±18.76μgm-3 (93.45±18.42μgm-3) in the pre-monsoon; 37.83±11.27μgm-3 (62.82±14.81μgm-3) in the monsoon and 40.14±16.66μgm-3 (64.72±22.40μgm-3) in the post-monsoon. About 76.67% of PM2.5 and 87.78% of PM10 concentrations were greater than NAAQ Standards in the winter on a daily basis and 46.74% and 36.96% in the pre-monsoon season. The backward trajectory analysis was also carried out through HYSPLIT model which suggests that the additional source of these pollutants during the winter and pre-monsoon season from the northwest and northern region of Patna. The ratios of PM10/PM2.5 observed at brick kilns cluster monitoring locations during the brick manufacturing period were significantly higher (0.87-4.48μgμgm-3) than other monitoring sites and increase the level of these pollutants over the city.

Dynamical Characteristics of Atmospheric Aerosols over IG Region

The dynamical characteristics of atmospheric aerosols over the Indo-Gangetic (IG) region are primarily dependent on the geographical settings and meteorological conditions. Detailed analysis of multi satellite data and ground observations have been carried out over three different cities i.e. Kanpur, Greater Noida and Amritsar during 2010-2013. Level-3 Moderate Resolution Imaging Spectroradiometer (MODIS) terra daily global grid product with spatial resolution of 1° × 1° shows the mean AOD at 500 nm wavelength value of 0.73, 0.70 and 0.67 with the standard deviation of 0.43, 0.39 and 0.36 respectively over Amritsar, Greater Noida and Kanpur. Our detailed analysis shows characteristic behavior of aerosols from west to east in the IG region depending upon the proximity of desert regions of Arabia. We have observed large influx of dusts from the Thar desert and Arabia peninsula during pre-monsoon season (April–June), highly affecting Amritsar which is close to the desert region.

Heterogeneity in topographic control on velocities of Western Himalayan glaciers

Studies of the seasonal and annual patterns of glacier velocities improve our understanding of the ice volume, topography, responses to climate change, and surge events of glaciers. Such studies are especially relevant and equally rare for the Himalayan glaciers, which supply many rivers that sustain some of the most heavily populated mountainous regions in the world. In particular, the control of the hypsometric distribution of geomorphometric parameters, such as slope, aspect, and curvature, on the dynamics of Himalayan glaciers have never been studied so far, at the river basin scale. Here, we present the degree to which topographic and hypsometric parameters affect the seasonal and annual average flow velocities of 112 glaciers in the Baspa River basin in the Western Indian Himalaya by analysing Global Land Ice Velocity Extraction from Landsat 8 (GoLIVE) datasets for the years 2013–2017. We observe, (i) significant heterogeneity in topographic controls on the velocities of these glaciers, (ii) elevation and the seasons play important roles in regulating the degree to which morphometric parameters (slope, aspect, and curvature) affect these velocities, (iii) a possible polythermal regime promoting both sliding and deformational forms of motion in a majority of these glaciers, and (iv) a detailed analysis of complex topographic controls within various elevation zones using a novel hypso-morphometric approach. These findings can help us to better model the dynamics of Himalayan glaciers and their responses to the future climatic scenarios. The inferences also suggest the need to incorporate dynamic topography in glacio-hydrological models in the wake of constant glacial evolutions.