Agricultural Area Diversification and Crop Water Demand Analysis: A Remote Sensing and GIS Approach

Large scale adoption of input intensive rice–wheat cropping system in the centrally located Jalandhar district of Indian Punjab has led to over-exploitation of ground water resources, intensive use of chemical fertilizers and deterioration of soil health. To overcome these shortfalls, in the present study, agricultural area diversification plan has been generated from agricultural area and crop rotation maps derived from remote sensing data (IRS P6-AWiFS and RADARSAT ScanSAR) along with few agro-physical parameters in GIS environment. Cropping system indices (area diversity, multiple cropping and cultivated land utilization) were also worked out from remote sensing data .Analysis of remote sensing data (2004–05) revealed that rice and wheat individually remained the dominant crops, occupy 57.8% and 64.9% of total agricultural area (TAA), respectively. Therefore, in the diversified plan, it is suggested that at least 39% of the current 40% TAA under rice–wheat rotation should be replaced by other low water requiring, high value and soil enriching crops, particularly in coarse textured alluvial plain having good quality ground water zones with low annual rainfall(<700 mm). This will reduce water requirement to the tune of 15,660 cm depth while stabilizing the production and profitability by crop area diversification without further degradation of natural resources.     

Monitoring Horizontal and Vertical Cropping Pattern and Dynamics in Bihar over a Decade (2001–2012) Based on Time-Series Satellite Data

Increasing population and natural disasters like drought, flood, cyclone etc., has impacted global agriculture area and hence continuously modifying cropping pattern and associated statistics. The present study analysed agriculture dynamics over one of the densely populated and disaster prone state (Bihar) in India and derived vital statistics (single, double and triple cropping area, and monthly, seasonal, annual and long term status at the state and district level) for the years 2001–2012. The study used time-series MODIS vegetation index (EVI; MOD13A2, 1 km, 16 day, 2001–2012), MODIS annual Land Cover product (MCD12Q1, 500 m, 2001–2012) and Global Land Cover map (Scasia_V4, 1 km, 2000; Globcover_V2.2, 300 m, 2005/2006 and V2.3, 2009, 300 m), and extracted horizontal (i.e., area change) and vertical (i.e., cropping intensification) agriculture change pattern. The results were inter-compared, and validated using government reports as well as with high spatial resolution data (IRS-LISS III 23.5 m). From 2001–2006 to 2007–2012, the net horizontal and vertical change in agriculture area is +145.24 and +907.82 km², respectively, and net change in seasonal crop area (winter, summer and monsoon) is +959.21, +1009.84 and ?1061.64 km², respectively. The districts which are located along the eastern part of Ganges experienced maximum positive changes and the districts along Gandak river in the north-western part of the study area experienced maximum negative changes. Overall, the study has quantified and revealed interesting space–time agriculture change patterns over 12 years including impacts caused by droughts and floods in the study area.     

Agricultural drought assessment at disaggregated level using AWiFS/WiFS data of Indian Remote Sensing satellites

Spatial differences in drought proneness and intensity of drought caused by differences in cropping patterns and crop growing environments within a district indicate the need for agricultural drought assessment at disaggregated level. The objective of this study is to use moderate resolution satellite images for detailed assessment of the agricultural drought situation at different administrative units (blocks) within a district. Monthly time composite NDVI images derived from moderate resolution AWiFS (60 m) and WiFS (180 m) images from Indian Remote Sensing satellites were analysed along with ground data on rainfall and crop sown areas for the kharif seasons (June – November) of 2002 (drought year), 2004 (early season drought) and 2005 (good monsoon year). The impact of the 2002 meteorological drought on crop area in different blocks of the district was assessed. The amplitude of crop condition variability in a severe drought year (2002) and a good year (2005) was used to map the degree of vulnerability of different blocks in the district to agricultural drought. The impact of early season deficit rainfall in 2004 on the agricultural situation and subsequent recovery of the agricultural situation was clearly shown. Agricultural drought assessment at disaggregated level using moderate resolution images is useful for prioritizing the problem areas within a district to undertake, in season drought management plans, such as alternate cropping strategies, as well as for end of the season drought relief management actions. The availability of ground data on rainfall, cropping pattern, crop calendar, irrigation, soil type etc., is very crucial in order to interpret the seasonal NDVI patterns at disaggregated level for drought assessment. The SWIR band of AWiFS sensor is a potential data source for assessing surface drought at the beginning of the season.     

Impact of cropping pattern changes on the exploitation of water resources: A remote sensing and GIS approach

A remote sensing-based approach was applied to study the impact of changes in cropping system on the exploitation of water resources in two districts namely Ludhiana in central Punjab and Muktsar in South-Western Punjab. Rice-wheat remained dominant rotation in Ludhiana while cottonwheat rotation was replaced partially by rice-wheat in Muktsar within a span of over five years (1998–99 to 2003–04). The solo rice-wheat system in Ludhiana district has resulted in large-scale ground water exploitation as is evident from the faster decline in water table (up to 0.9m year⁻¹) and higher tube-wells density (440 per 1000 ha). As a result, nearly 60 per cent of the total area of Ludhiana district has the water table depth greater than 10m and in some blocks, it has reached to a depth of 22 m. In Muktsar district, corresponding rise in water table is 0.2m per year and tube well density is 114 per 1000 ha. Irrigation water associated with paddy cultivation in Ludhiana and adjoining areas moves laterally through the buried paleo-channels of Sutlaj in the deeper soil profile and gets accumulated in the basin lands of Muktsar and adjoining areas and causes an extra water flux and subsequent rise in water table, recorded at 3 to 7m depth. To minimize the hydrological imbalance of the state, it is suggested to diversify some of the area in the central districts from irrigation water intensive rice-wheat system to less water intensive cropping system.     

Results on potential use of simulated IRS-1C WiFS data for crop monitoring

Some of the basic requirements for cropping system analysis are updated information on crops grown, their phenological behaviour, method and duration of establishment and harvest, inter and intra crop variability, sequential cropping patterns. The next generation Indian Remote Sensing Satellite with high repeat cycle opens new possibility of crop surveys to derive such information. In this study, an attempt has been made to analyse cropping system at district level using simulated IRS-1C Wide Field Sensor (WiFS) data. Data acquired for nineteen dates during 1992–93 season for Bardhaman district, West Bengal has been used. It was feasible to derive accurate information on cropping pattern, crop rotation, crop duration, progress of harvest, crop growth profiles and annual crop acreage using multidate data. It was observed that even a seven to eight day interval of data acquisition during critical growth periods significantly affected classification and identification accuracy.     

Analysis of cropping pattern changes in Bathinda district,Punjab

A study was conducted in the Bathinda district of Punjab state for mapping the cropping pattern and crop rotation, monitoring long term changes in cropping pattern by using the satellite based remote sensing data along other spatial and non-spatial collateral data. Multi-date IRS LISS I and IRS WiFS sensor data have been used for this study. Cropping pattern maps and crop rotation maps were generated for the years 1988-89 and 1998-99. The present study has shown the increase of cropping intensity significantly, mainly due to increase in rice area. However, crop diversity has decreased mainly due to decline in the area under the minor crops like pearl millet, gram, rapeseed/ mustard. There is increase in area coverage of cotton-wheat and rice-wheat rotation, at the expense of the minor crops.     

Monitoring of deforestation and land use changes (1925–2012) in Idukki district,Kerala, India using remote sensing and GIS

The studies on forest cover change can reveal the status of forests and facilitate for its conservation planning. Idukki is the largest district in the state of Kerala having a total geographical area of 5019 km². The objectives of the present study are to map forest cover in Idukki district using multi-temporal remote sensing data (1975, 1990, 2001 and 2012) and topographical maps (1925), to analyze the trends in deforestation and land use changes. Overall statistics for the period of 1925 indicate that about 4675.7 km² (93.2 %) of the landscape was under forest. The forest cover in 2012 was estimated as 2613.4 km² (52.1 %). Recently, due to the implementation of policies and protection efforts, the rate of deforestation was greatly reduced. The commencement of hydroelectric projects during 1925–1990 responsible for an increase of area under water bodies by inundating other land uses. The long term analysis shows agricultural area been decreasing and commercial plantations been increasing in the district. There has been a significant increase in the area of plantations from 1236.2 km² (1975) to 1317.3 km² (2012).     

Mapping the Extent of Agroforestry Area in Different Altitudes of Tehri District,North Western Himalaya,India Through GIS and Remote Sensing Data

Assessment of area under agroforestry in Tehri district of North Western Himalaya, Uttarakhand, India has been done using GIS and remote sensing technology. The study district characterized by hilly terrain with varying elevations from 288 m to more than 2800 m and generally gentle slopes, valleys, flat land covers and agricultural terraces. High-resolution satellite imageries (spatial resolution 5.8 m) were used in this study for land uses and land covers classification. According to unsupervised classification, highest area was found under forest class (65.22%) followed by cropland (20.41%). Considerable area was also found under snow cover (9.45%) in the district. Area under agroforestry was estimated to be 5572.26 ha (1.53%) by this method, whereas it was estimated to be 7029.06 ha (1.93%) by supervised classification. Estimated cropland area comes out to be about 20.0%. An accuracy of 86.5% was found in this classification for agroforestry class. Highest area under agroforestry of 3707.36 ha was obtained in 1200–2000 m elevations followed by 2231.26 ha in 288–1200 m elevations. Negligible area was found on high elevation zones of more than 2800 m. The major agroforestry systems of dominated by Grewia oppositifolia (Bhimal), Celtis australis (Kharik) and Quercus leucotrichophora (Banj) were identified and mapped and remaining systems were grouped as others class. Estimated area under G. oppositifolia, C. australis and Q. leucotrichophora based systems come out to be 2330.82, 1456.80 and 1129.10 ha, respectively. These systems are multiple usufructs are food, fuelwood, fodder, fiber and small timber. It has been observed from the accuracy assessment that the estimates of area under agroforestry obtained under this study are reliable.     

Local variability in the timing and intensity of tropical dry forest deciduousness is explained by differences in forest stand age

Tropical Dry Forest deciduousness is a behavioral response to climate conditions that determines ecosystem-level carbon uptake, energy flux, and habitat conditions. It is regulated by factors related to stand age, and landscape scale variability in deciduous phenology may affect ecosystem functioning in forests throughout the tropics. This study determines whether observed phenological differences are explainable by forest age in the southern Yucatán Peninsula in Mexico, where forest clearing for shifting cultivation has created a mosaic of forest stands of varying age. Matched-pair statistical tests compare neighboring forest pixels of different age class (12–22 years versus 22+ years) and detect significant differences in Moderate Resolution Imaging Spectroradiometer (MODIS) Enhanced Vegetation Index (EVI)-derived metrics related to the timing and intensity of deciduousness during three dry seasons (2008–2011). In all seasons, young forests exhibit significantly more intense deciduousness, measured as total seasonal change of EVI normalized by annual maximum EVI (p < 0.001), and larger normalized EVI change during successive dry season months relative to start-of-dry-season EVI (p < 0.001), than neighboring older forests subject to similar environmental conditions.     

Use of high-resolution satellite data,GIS and NRCS-CN technique for the estimation of rainfall-induced run-off in small catchment of Jharkhand India

The present study demonstrates the use of NRCS-CN technique for rainfall-induced run-off estimation using high-resolution satellite data for small watershed of Palamu district, Jharkhand. The CN model was applied to the daily rainfall data of 15 years (1986–2000) along with use of large-scale thematic maps (1:10,000) pertaining to land use/land cover using IRS-P6 LISS-IV satellite data. The LU/LC map was spatially intersected with the hydrological soil group map to calculate the watershed area under different hydrological similar units for assigning CN values to compute discharge. The study showed that Daltonganj watershed exhibits an average run-off volume of 7,881,019 m³ from an average cumulative monsoon rainfall of 821 mm and the average actual direct run-off generated during the southwest monsoon season was 203 mm. The strong correlation between rainfall and run-off as well as between observed run-off and estimated run-off indicated high accuracy of run-off estimation by NRCS-CN technique.     

Analyzing ice dynamics using Sentinel-1 data at the Solheimajoküll Glacier,Iceland

ABSTRACT

The climate in southern Iceland has warmed over the last 70 years, resulting in accelerated glacier dynamics at the Solheimajoküll glacier. In this study, we compare glacier terminus locations from 1973 to 2018, to changes in climate across the study area, and we derive ice-surface velocities (2015–2018) from satellite remote-sensing imagery (Sentinel-1) using the offset-tracking method. There have been two regional temperature trends in the study period: cooling (1973–1979) and warming (1980–2018). Our results indicate a time lag of about 20 years between the onset of glacier retreat (?53 m/year since 2000) and the inception of the warming period. Seasonally, the velocity time series suggest acceleration during the summer melt season since 2016, whereas glacier velocities during accumulation months were constant. The highest velocities were observed at high elevations where the ice-surface slope is the steepest. We tested several scenarios to assess the hydrological time response to glacier accelerations, with the highest correlations being found between one and 30 days after the velocity estimates. Monthly correlation analyses indicated inter-annual and intra-annual variability in the glacier dynamics. Additionally, we investigate the linkage between glacier velocities and meltwater outflow parameters as they provide useful information about internal processes in the glacier. Velocity estimates positively correlate with water level and negatively correlate with water conductivity between April and August. There is also a disruption in the correlation trend between water conductivity and ice velocity in June, potentially due to a seasonal release of geothermal water.     

Application of remote sensing technology for Impact Assessment of Watershed Development Programme

The current study was taken up to investigate the utility of remote sensing and GIS tools for evaluation of Integrated Wasteland Development Programme (IWDP) implemented during 1997–2001 in Katangidda Nala watershed, Chincholi taluk, Gulbarga district, Karnataka. The study was carried out using IRS 1C, LISS III data of December 11, 1997 (pre-treatment) and November 15, 2002 (post-treatment) covering the watershed to assess the changes in land use / land cover and biomass that have changed over a period of five years (1997–2002). The images were classified into different land use/land cover categories using supervised classification by maximum likelihood algorithm. They were also classified into different biomass levels using Normalized Difference Vegetation Index (NDVI) approach. The results indicated that the area under agriculture crops and forest land were increased by 671 ha (5.7%) and 1,414 ha (11.94%) respectively. This is due to the fact that parts of wastelands and fallow lands were brought into cultivation. This increase in the area may be attributed to better utilization of surface and ground waters, adoption of soil and water conservation practices and changes in cropping pattern. The area under waste lands and fallow lands decreased by 1,667 ha (14.07%) and 467 ha (3.94%), respectively. The vegetation vigour of the area was classified into three classes using NDVI. Substantial increase in the area under high and low biomass levels was observed (502 ha and 19 ha respectively). The benefit-cost analysis indicates that the use of remote sensing and GIS was 2.2 times cheaper than the conventional methods. Thus, the repetitive coverage of the satellite data provides an excellent opportunity to monitor the land resources and evaluate the land cover changes through comparison of images for the watershed at different periods.     

Variations in the Ice Phenology and Water Level of Ayakekumu Lake,Tibetan Plateau,Derived from MODIS and Satellite Altimetry Data

The alpine lakes on the Tibetan Plateau (TP) are highly sensitive to variations in climate changes, and the lake ice phenology and water level are considered to be direct indicators of regional climate variability. In this study, we first used 14 years of moderate resolution imaging spectroradiometer snow cover products to analyse the freeze dates, ablation dates, and ice coverage durations. The lake level changes during 2002–2015 were estimated, derived from satellite altimetry and Hydroweb data. Unexpectedly, the freeze dates of lake ice greatly advanced, and the ablation dates were markedly delayed. The complete freezing duration lengthened by approximately 77 days. As a result of the warm-wet climate in the northern TP, the lake area expanded from 770 to 995 km² during 2002–2015, and the water levels rose by 4.2 m in total, at a rate of 0.3 m/year. The progressive expansion of Ayakekumu Lake profoundly affected the ice phenology. Larger water volume with larger thermal capacity likely led to the delaying of ablation dates, with the freezing point depression caused by decreasing salinity. Some new narrow and shallow bays located in southern and eastern Ayakekumu Lake were conducive to early freezing of ice. Additionally, the changes in air temperature, precipitation, potential evaporation, and sunshine duration may be related to the prolonged ice cover duration since 2002. In sum, accurate measurements of lake ice and water levels are critical for understanding the water resource balance and hydrologic cycle in arid or semi-arid regions of China.     

Assessing changes in the above ground biomass and carbon stocks of Lidder valley,Kashmir Himalaya,India

The changes in the land use and land cover (LULC), above ground biomass (AGB) and the associated above ground carbon (AGC) stocks were assessed in Lidder Valley, Kashmir Himalaya using satellite data (1980–2013), allometric equations and phytosociological data. Change detection analysis of LULC, comprising of eight vegetation and five non-vegetation types, indicated that 6% (74.5 km²) of the dense evergreen forest has degraded. Degraded forest and settlement increased by 20 and 52.8 km², respectively. Normalized difference vegetation index was assessed and correlated with the field-based biomass estimates to arrive at best-fit models for remotely sensed AGB estimates for 2005 and 2013. Total loss of 1.018 Megatons of AGB and 0.5 Megatons of AGC was estimated from the area during 33-year period which would have an adverse effect on the carbon sequestration potential of the area which is already facing the brunt of climate change.     

Change in regional land subsidence in Beijing after south-to-north water diversion project observed using satellite radar interferometry

ABSTRACT

Over-exploitation of groundwater has caused severe land subsidence in Beijing during the past two decades. Since the middle route of South-to-North Water Diversion Project (SNWDP), the biggest water diversion project in China, started to deliver water to Beijing in December 2014, the groundwater shortage has been greatly alleviated. This study aims to analyze the impact of SNWDP on the spatiotemporal evolution of land subsidence in Beijing. Change in surface displacement in Beijing after SNWDP was retrieved and the spatiotemporal patterns of the change were analyzed based on long time-series Envisat Advanced Synthetic Aperture Radar (ASAR) (2004–2010), Radarsat-2 (2011–2014), and Sentinel-1 (2015–2017) satellite datasets using Permanent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) techniques. Land subsidence unevenness index (LSUI) was proposed to represent the spatial unevenness of surface displacement. PS-Time approach was then adapted to examine the time series evolution of LSUI. The results showed that the InSAR measurements agree well with leveling measurements with R² over 0.96. Although the maximum annual displacement rate reached ?159.7 mm/year by 2017, over 57% of the area within 25 mm/year contour line showed decreasing or unchanged displacement rate after the south-north water delivered to Beijing. The settlement rate in Chaoyang-Dongbalizhuang (CD) subsidence center has decreased for 26 mm/year from 2011–2014 to 2015–2017. Only around 15% of the area experienced continued accelerating settlement rate through the three time periods, which was mainly located in the area with the compressible layer thickness over 190 m, while the magnitude of velocity increment considerably decreased after SNWDP. Land subsidence unevenness, represented by LSUI, developed more slowly after SNWDP than that during 2011–2014. However, LSUI at the edge of settlement funnel has kept developing and reached 1.7‰ in 2017. Decreasing groundwater level decline after SNWDP and the positive relationship (R² > 0.74) between land subsidence and groundwater level clearly showed impacts of SNWDP on the alleviating land subsidence. Other reasons include geological background, increasing precipitation, and strict water management policies implemented during these years.     

利用Landsat TM和IRS P6 LISS III卫星影像及GIS对耕地和耕种模式的时空动态分析(英文)

The study of the spatial patterns and temporal changes of cropland is important to understand the underlying factors and the functional effects of the agricultural landscape. On the other hand, crop dynamics mapping is essential to know the overall agro-spatial diversity of the area. Therefore, this paper addressed a spatio-temporal analysis of cropland and cropping pattern change in the Bogra district of Bangladesh over the last 16 years (between 1988/89 and 2004/05). In this paper, crop mapping from multi-temporal and multi-sensor satellite images was described. Landsat TM and IRS P6 LISS III satellite images were used with GIS for spatial dynamics of cropland and cropping pattern change analysis. First, seasonal cropland maps were derived from object-based classification of satellite images, then two-date classified image differencing with GIS overlay technique and decision rules were applied. Cropping pattern change was analyzed in a spatial and quantitative way for the 16 years and for this, Integrated Land and Water Information System (ILWIS) and Land Change Modular (LCM) of IDRISI Andes were used. The results showed that in the area, mono crop cultivation was found in summer, but in winter, areas under different crop cultivation had changed dramatically. Change analysis showed that the changes mainly occurred in the north northwest and southwest of the areas, and during the time the highest change area was found under the rice-potato pattern.      

Urban Built-up Area Assessment of Ranchi Township Using Cartosat-I Stereopairs Satellite Images

In the present study, the Cartosat-I digital elevation model (DEM) was utilized to deduce the vertical characteristics of Ranchi urban area and its relation to long term built-up expansion (1927–2010). The DEM represents moderate variation in terrain relief ranging from 595 m to 754 m with majority of area exhibiting upto 3° of slope and 3° to 6° indicating flat to undulating nature of terrain in Ranchi township. The DEM was used to generate location of sinks within urban area, which are generally delineated along the drainage channels, adjacent to high-rise built-up land and along the elevated road network. The pattern of urban sprawl over the eight decades (1927–2010) were examined with reference to terrain relief zones, which indicated that the built-up growth was mainly taken place over the elevation range of moderate (620–660 m) (67.0%) and high relief (660–680 m) (19.8%) zones. Although earlier preference for built-up development was more in high elevation zones (660–680 m), the low elevation zones (<600–620 m) are now preferred for multistoried built-up land development where better groundwater availability occur. The spatial pattern of vertical growth of built-up land was assessed using contour density obtained from Cartosat-I DEM. The results show that the high density contours predominately correspond to hilly area and high-rise buildings at majority of locations. The urban sprawl pattern and population trend exhibited rapid increase in vertical built-up growth after 1996 indicating beginning of urban densification in Ranchi township.     

The market value of cultural heritage in urban areas: an application of spatial hedonic pricing

The current literature often values intangible goods like cultural heritage by applying stated preference methods. In recent years, however, the increasing availability of large databases on real estate transactions and listed prices has opened up new research possibilities and has reduced various existing barriers to applications of conventional (spatial) hedonic analysis to the real estate market. The present paper provides one of the first applications using a spatial autoregressive model to investigate the impact of cultural heritage—in particular, listed buildings and historic–cultural sites (or historic landmarks)—on the value of real estate in cities. In addition, this paper suggests a novel way of specifying the spatial weight matrix—only prices of sold houses influence current price—in identifying the spatial dependency effects between sold properties. The empirical application in the present study concerns the Dutch urban area of Zaanstad, a historic area for which over a long period of more than 20 years detailed information on individual dwellings, and their market prices are available in a GIS context. In this paper, the effect of cultural heritage is analysed in three complementary ways. First, we measure the effect of a listed building on its market price in the relevant area concerned. Secondly, we investigate the value that listed heritage has on nearby property. And finally, we estimate the effect of historic–cultural sites on real estate prices. We find that, to purchase a listed building, buyers are willing to pay an additional 26.9 %, while surrounding houses are worth an extra 0.28 % for each additional listed building within a 50-m radius. Houses sold within a conservation area appear to gain a premium of 26.4 % which confirms the existence of a ‘historic ensemble’ effect.     

Reservoir capacity estimation using SARAL/AltiKa altimetry data coupled with Resourcesat P6-AWiFS and RISAT 1 microwave data

Reservoir water levels extracted from SARAL/AltiKa GDR data for the period 2013–2014 and water spread areas delineated from Resourcesat P6-AWiFS sensor and RISAT 1 microwave data corresponding to SARAL/AltiKa cycles were used for assessment of reservoir capacity in the Mayurakshi reservoir, Jharkhand state, India. It was found that the reservoir capacity based on the SARAL is around 474.62 Mm³ in comparison to in situ based estimate i.e. around 486.6 Mm³, indicating variation of <3%. Further, comparison of these estimates computed using SARAL and in situ with original reservoir capacity (547.59 Mm³) indicated loss of reservoir capacity is around 13.33 and 11.14%, respectively, within a span of 59 years. The hydrographic survey in the year 1999–2000 also proved that the storage capacity has reduced from 547.6 Mm³ in 1955 to 474.8 Mm³ indicating loss of nearly 13.3 % of total live capacity over period of 45 years.     

Frontal recession of Parkachik Glacier between 1971-2015, Zanskar Himalaya using remote sensing and field data

Parkachik Glacier is located in the Suru sub-basin of the Upper Indus River, Zanskar Himalaya. The Glacier has been analysed using Corona KH-4B (1971), Landsat-TM (1999), field survey (2015), Google Earth^TM (2015) and ASTER GDEM (2015) for frontal recession and area changes. Overall, from 1971 to 2015, the Glacier has retreated by 127 ± 0.09 m i.e. (0.75 ± 0.07%) at a rate of 2.9 ± 0.004 ma^?1 with a simultaneous decrease in area from 49.5 to 48.8 km² i.e. 740 ± 0.7 m² (1.5 ± 0.09%) at a rate of 74 ± 0.7 m²a^?1. However, during recent decade (1999–2015), the rate of glacier recession of 3.9 ± 0.004 ma^?1 with a corresponding area loss of 500 ± 0.74m² (1 ± 0.1%) was higher than the retreat rate of 2.3 ± 0.001 ma^?1 and an area loss of 240 ± 0.02m² (0.48 ± 0.08%) during 1971–1999. In the field, the evidences of glacier recession are present in the form of separated dead ice blocks from the main Glacier, recessional dumps/moraines, active ice calving activity and a small proglacial pond/lake at the terminus/snout of the Glacier. However, the recession over the studied period has been very slow and is controlled by its topographic configuration, particularly the large altitudinal range (6030–3620 m), almost northerly aspect and steep slope (average ~ 30°).