Snow Cover Monitoring in Qinghai-Tibetan Plateau Based on Chinese Fengyun-3/VIRR Data

Snow cover monitoring in the Qinghai-Tibetan Plateau is very important to global climate change research. Because of the geographic distribution of ground meteorological stations in Qinghai-Tibetan Plateau is too sparse, satellite remote sensing became the only choice for snow cover monitoring in Qinghai-Tibetan Plateau. In this paper, multi-channel data from Visible and Infrared Radiometer (VIRR) on Chinese polar orbiting meteorological satellites Fengyun-3(FY-3) are utilized for snow cover monitoring, in this work, the distribution of snow cover is extracted from the normalized difference snow index(NDSI), and the multi-channel threshold from the brightness temperature difference in infrared channels. Then, the monitoring results of FY-3A and FY-3B are combined to generate the daily composited snow cover product. Finally, the snow cover products from MODIS and FY-3 are both verified by snow depth of meteorological station observations, result shows that the FY-3 products and MODIS products are basically consistent, the overall accuracy of FY-3 products is higher than MODIS products by nearly 1 %. And the cloud coverage rate of FY-3 products is less than MODIS by 2.64 %. This work indicates that FY-3/VIRR data can be reliable data sources for monitoring snow cover in the Qinghai-Tibetan Plateau.     

Snow permittivity retrieval inversion algorithm for estimating snow wetness

The environmental satellite (ENVISAT) advanced synthetic aperture radar (ASAR) offers the opportunity for monitoring snow parameters with dual polarization and multi-incidence angle. Snow wetness is an important index for indicating snow avalanche, snowmelt runoff modelling, water supply for irrigation and hydropower stations, weather forecasts and understanding climate change. We used a first-order scattering model that includes both volume and air/snow surface scattering based on a developed inversion model to estimate snow dielectric constant, which can be further related for estimating snow wetness. Comparison with field measurement showed that the correlation coefficient for snow permittivity estimated from ASAR data was observed to be 0.8 at 95% confidence interval and model bias was observed as 2.42% by volume at 95% confidence interval. The comparison of ASAR-derived snow permittivity with ground measurements shows the average absolute error 2.5%. The snow wetness range varies from 0 to 15% by volume.     

Synergy of in situ and space borne observation for snow depth mapping in the Swiss Alps

The Swiss Federal Institute for Snow and Avalanche Research in Davos (SLF) provides snow depth maps for Switzerland on a spatial resolution of 1 km × 1 km. These snow depth maps are derived from snow station measurements using a spatial interpolation method based on the dependency of snow depth and altitude. During a winter season the number of operating snow stations varies and the area-wide snow depth interpolation becomes increasingly difficult in spring. The objective of the study is to develop an operational and near-real time method to calculate snow depth maps using a combination of in situ snow depth measurements and the snow cover extent provided from space borne observations. The operational daily snow cover product obtained from the polar-orbiting NOAA-AVHRR satellite is used to gain an additional set of virtual snow stations to densify the in situ measurements for an improved spatial interpolation. The capacity of this method is demonstrated on selected days during winter 2005. Cross-validation tests are conducted to examine the quantitative accuracy of the synergetic interpolation method. The error estimators prove the decrease in error variance and increase of overall accuracy pointing out the high capacity of this new interpolation method that can be run in near real-time over a large horizontal domain at high horizontal resolution. A solid method for snow–no snow classification in the processing of the satellite data is essential to the quality of the snow depth maps.     

MTSAT-2静止气象卫星中国区域雪盖监测

利用MTSAT-2静止气象卫星数据开展了中国区域的雪盖监测研究,结合MODIS雪盖产品及站点雪深观测数据对判识结果进行对比分析和验证。首先,根据MTSAT-2静止气象卫星数据特点,进行角度效应校正及多时相数据合成,以减少云对图像的影响;其次,根据多个雪盖判识因子建立中国区域雪盖判识算法;最后,对比分析2011年1月份MTSAT-2和MODIS雪盖判识结果,并使用站点观测数据进行精度验证。研究表明:(1)MTSAT-2雪盖判识受云影响比例约30%,MODIS雪盖产品受云影响比例约60%,MTSAT-2去云效果明显。(2)无云情况下,MTSAT-2雪盖判识和MODIS雪盖产品判识精度均高于92%;有云覆盖时,MTSAT-2判识精度约65%,优于MODIS雪盖产品35%的判识精度。(3)MTSAT-2静止气象卫星在保持高积雪判识精度的前提下,可以更有效减少云对雪盖判识影响,实时获取更多地表真实信息。该研究对中国区域雪盖信息准确监测、气候变化研究以及防灾减灾等具有重要意义。     

Accuracy assessment of MODIS/Terra snow cover product for parts of Indian Himalayas

Snow cover mapping is important for snow and glacier-related research. The spatial and temporal distribution of snow cover area is a fundamental input to the atmospheric models, snowmelt runoff models and climate models, as well as other applications. Daily snow cover maps from Moderate Resolution Imaging Spectroradiometer (MODIS) Terra satellite were retrieved for the period between 2004 and 2007, and pixels in these images were classified as cloud, snow or snow-free. These images have then been compared with ground snow depth (SD) measurements from the four observatories located at different parts of Himalayas. Comparison of snow maps with in situ data showed good agreement with overall accuracies in between 78.15 and 95.60%. When snow cover was less, MODIS data were found to be less accurate in mapping snow cover region. As the SD increases, the accuracy of MODIS snow cover maps also increases.     

Potential of multispectral reflectance for assessment of snow geophysical parameters in Solang valley in the lower Indian Himalayas

ABSTRACT

Snow geophysical parameters such as wetness, density and permittivity are a significant input in hydrological models and water resource management. In this paper, we utilize the triangle method based on a feature space developed with the near-infrared (NIR) reflectance and the Normalized Differenced Snow Index (NDSI) for the estimation of surface snow wetness, permittivity and density. The triangular feature space based on NIR reflectance and NDSI is parameterized to yield a linear relationship between the snow wetness and the NIR reflectance. Snow density and permittivity are derived based on the least squares solution of empirical relations based on the observations of surface snow wetness. The proposed methodology was evaluated using Sentinel-2 data, and the modeled snow geophysical parameters were validated with respect to field measurements. Based on the results, it was inferred that the NIR reflectance varies linearly with the liquid water content in the snow. A good agreement was determined between the modeled and measured parameters for wet snow conditions as observed by the coefficient of determination of 0.968, 0.521 and 0.969 for the snow wetness, density and permittivity (real part), respectively. The proposed approach can be significantly utilized with unmanned aerial sensors for monitoring of physical properties of fresh or wet snow and is thus expected to contribute considerably in hydrological applications and avalanche studies.     

Accuracy assessment of four cloud-free snow cover products over the Qinghai-Tibetan Plateau

Four up-to-date daily cloud-free snow products – IMS (Interactive Multisensor Snow products), MOD-SSM/I (combination of the MODIS and SSM/I snow products), MOD-B (Blending method basing on the MODIS snow cover products) and TAI (Terra–Aqua–IMS) – with high-resolutions over the Qinghai-Tibetan Plateau (QTP) were comprehensively assessed. Comparisons of the IMS, MOD-SSM/I, MOD-B and TAI cloud-free snow products against meteorological stations observations over 10 snow seasons (2004–2013) over the QTP indicated overall accuracies of 76.0%, 89.3%, 92.0% and 92.0%, respectively. The K_hat values of the IMS, MOD-SSM/I, MOD-B and TAI products were 0.084, 0.463, 0.428 and 0.526, respectively. The TAI products appear to have the best cloud-removal ability among the four snow products over the QTP. Based on the assessment, an I-TAI (Improvement of Terra–Aqua–IMS) snow product was proposed, which can improve the accuracy to some extent. However, the algorithms of the MODIS series products show instability when identifying wet snow and snow under forest cover over the QTP. The snow misclassification is an important limitation of MODIS snow cover products and requires additional improvements.     

Integration of the MODIS Snow Cover Produced Into Snowmelt Runoff Modeling

Because of the difficulty of monitoring and measuring snow cover in mountainous watersheds, satellite images are used as an alternative to mapping snow cover to replace the ground operations in the watershed. Snow cover is one of the most important data in simulation snowmelt runoff. The daily snow cover maps are received from Moderate Resolution Imaging Spectroradiometer (MODIS), and are used in deriving the snow depletion curve, which is one of the input parameters of the snowmelt runoff model (SRM). Simulating Snowmelt runoff is presented using SRM model as one of the major applications of satellite images processing and extracting snow cover in the Ghara - Chay watershed. The first results of modeling process show that MODIS snow covered area product can be used for simulation and forecast of snowmelt runoff in Ghara - Chay watershed. The studies found that the SCA results were more reliable in the study area.     

An Introduction to MODISI and SCMOD Methods for Correction of the MODIS Snow Assessment Algorithm

Detection, monitoring and precise assessment of the snow covered regions is an important issue. Snow cover area and consequently the amount of runoff generated from snowmelt have a significant effect on water supply management. To precisely detect and monitor the snow covered area we need satellite images with suitable spatial and temporal resolutions where we usually lose one for the other. In this study, products of two sensors MODIS and ASTER both on board of TERRA platform having low and high spatial resolution respectively were used. The objective of the study was to modify the snow products of MODIS by using simultaneous images of ASTER. For this, MODIS snow index image with high temporal resolution were compared with that of ASTER, using regression and correlation analysis. To improve NDSI index two methods were developed. The first method generated from direct comparison of ASTER averaged NDSI with those of MODIS (MODISI). The second method generated by dividing MODIS NDSI index into 10 codes according to their percentage of surface cover and then compared the results with the difference between ASTER averaged and MODIS snow indices (SCMOD). Both methods were tested against some 16 MODIS pixels. It is found that the precision of the MODISI method was more than 96%. This for SCMOD was about 98%. The RMSE of both methods were as good as 0.02.     

Evaluating MODIS snow products for modelling snowmelt runoff: Case study of the Rio Grande headwaters

Snow-covered area (SCA) is a key variable in the Snowmelt-Runoff Model (SRM) and in other models for simulating discharge from snowmelt. Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM + ) or Operational Land Imager (OLI) provide remotely sensed data at an appropriate spatial resolution for mapping SCA in small headwater basins, but the temporal resolution of the data is low and may not always provide sufficient cloud-free dates. The coarser spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) offers better temporal resolution and in cloudy years, MODIS data offer the best alternative for mapping snow cover when finer spatial resolution data are unavailable. However, MODIS’ coarse spatial resolution (500 m) can obscure fine spatial patterning in snow cover and some MODIS products are not sensitive to end-of-season snow cover. In this study, we aimed to test MODIS snow products for use in simulating snowmelt runoff from smaller headwater basins by a) comparing maps of TM and MODIS-based SCA and b) determining how SRM streamflow simulations are changed by the different estimates of seasonal snow depletion. We compared gridded MODIS snow products (Collection 5 MOD10A1 fractional and binary SCA; SCA derived from Collection 6 MOD10A1 Normalised Difference Snow Index (NDSI) Snow Cover), and the MODIS Snow Covered-Area and Grain size retrieval (MODSCAG) canopy-corrected fractional SCA (SCA_MG), with reference SCA maps (SCA_REF) generated from binary classification of TM imagery. SCA_MG showed strong agreement with SCA_REF; excluding true negatives (where both methods agreed no snow was present) the median percent difference between SCA_REF and SCA_MG ranged between −2.4% and 4.7%. We simulated runoff for each of the four study years using SRM populated with and calibrated for snow depletion curves derived from SCA_REF. We then substituted in each of the MODIS-derived depletion curves. With efficiency coefficients ranging between 0.73 and 0.93, SRM simulation results from the SCA_MG runs yielded the best results of all the MODIS products and only slightly underestimated discharge volume (between 7 and 11% of measured annual discharge). SRM simulations that used SCA derived from Collection 6 NDSI Snow Cover also yielded promising results, with efficiency coefficients ranging between 0.73 and 0.91.In conclusion, we recommend that when simulating snowmelt runoff from small basins (<4000 km²) with SRM, we recommend that users select either canopy-corrected MODSCAG or create their own site-specific products from the Collection 6 MOD10A1 NDSI.     

MODIS大气校正精度评价及其对表层雪密度提取影响

以古尔班通古特沙漠为研究区,以中分辨率成像光谱仪（MODIS）为遥感数据源,结合ASD FieldSpec准同步实测积雪反射光谱数据对FLAASH大气校正能力进行了评价。研究表明： ①校正后的MODIS各波段积雪反射率与准同步实测积雪反射率波形相似, 在第1～7波段整体相关系数达0.82,表明FLAASH大气校正能极大地提高MODIS地物识别能力; ②校正后的MODIS 第6波段反射率和归一化差值积雪指数（NDSI）与实测雪密度呈线性相关,可用回归拟合构建MODIS雪密度遥感计算模式。     

Monitoring of Snow Cover Variability in Chenab Basin Using IRS AWiFS Sensor

Monitoring of seasonal snow cover is important for many applications such as melt runoff estimation, climate change studies and strategic requirements. Contribution of seasonal snow melt runoff of Chenab River is significant and important to meet hydrological requirement at foothills. Seasonal snow cover of Chandra, Bhaga, Miyar, Bhut, Warwan and Ravi, six major tributaries of Chenab River, becomes crucial to assess the water availability. In addition, altitudinal distribution of snow cover significantly influences the melt runoff which is highly sensitive to minor variations in atmospheric temperature. In this investigation, remote sensing based Normalized Difference Snow Index technique has been used to generate 10 daily snow cover product. Snow cover monitoring of all the sub-basins were carried out for 10 years from 2004–2005 to 2013–2014 during hydrological year (October to June) using Advanced Wide Field Sensor (AWiFS) of Indian remote sensing satellite (IRS). Accumulation and ablation patterns of snow cover have also been analyzed for the six sub-basins. Accumulation and ablation pattern of snow cover, from 2004 to 2014 which shows slightly increasing trend for all the sub-basins. Meteorological data of Kelong at Bhaga sub-basin was also analysed. Average monthly snow line altitude was estimated for all the sub-basins using hypsographic curve. Chandra and Bhaga sub-basins are at higher altitude and Ravi sub-basin is at lower altitude. It was also observed that areal extent of snow reaches to lower altitude during last 5 years, particularly in Ravi sub-basin.     

MODIS数据在积雪检测中的应用

积雪作为影响环境的一个因素,是非常重要的。自1999年Terra卫星升空以来,MODIS数据在环境监测的各个方面得到了广泛的应用。由于MODIS数据的高光谱、高空间分辨率、高时间分辨率等特征,越来越多地应用到积雪监测方面。本文就MODIS数据的积雪检测算法进行了探讨,对森林中雪的检测以及云和雪的区分进行了大量的研究。结果显示：MODIS数据对积雪检测非常有效。     

An Assessment of Snow Avalanche Paths and Forest Dynamics Using Ikonos Satellite Data

Abstract

Ikonos panchromatic and multispectral satellite data were acquired in October 2000 and August 2002 for a test area along US Highway 2, the southern border of Glacier National Park (GNP), Montana, USA. The research goals were to map snow avalanche paths and to characterize vegetation patterns in selected paths for longitudinal (i.e., source, track, and runout) and transverse (i.e., inner, flanking, outer) zones as part of a study of forest dynamics and nutrient flux from paths into terrestrial and aquatic systems. In some valleys, as much as 50 percent of the area may be covered by snow avalanche paths, and as such, serve as an important carbon source servicing terrestrial and aquatic ecosystems. Snow avalanches move woody debris down‐slope by snapping, tipping, trimming, and excavating branches, limbs, and trees, and by injuring and scaring trees that remain in‐place. Further, snow avalanches alter the vegetation structure on paths through secondary plant succession of disturbed areas. Contrast and edge enhancements, Normalized Difference Vegetation Index (NDVI), and the Tasseled Cap greenness and wetness transformations were used to examine vegetation patterns in selected paths that were affected by high magnitude snow avalanches during the winter of 2001-2002. Using image transects organized in longitudinal patterns in paths and in forests, and transects arranged in transverse patterns across the sampled paths, the Tasseled Cap transforms (and NDVI values) were plotted and assessed. Preliminary results suggest that NDVI patterns are different for paths and forests, and Tasseled Cap greenness and wetness patterns are different for longitudinal and transverse zones that describe the morphology of snow avalanche paths. The differentiation of paths from the background forest and the characterization of paths by morphometric zones through remote sensing has implications for mapping forest disturbances and dynamics over time and for large geographic areas and for modeling nutrient flux in terrestrial and aquatic systems.     

用被动微波AMSR数据反演地表温度及发射率的方法研究

 针对对地观测卫星多传感器的特点，提出了借助MODIS地表温度产品从被动微波数据中反演地表温度的方法。即利用MODIS地表温度产品和AMSR不同通道之间的亮度温度，建立地表温度的反演方程。该方法克服了以往需要测量同步数据的困难，为不同传感器之间的参数反演相互校正和综合利用多传感器的数据提供实际应用和理论依据。文中以MODIS地表温度产品作为评价标准，对方法进行检验，其平均误差为2~3℃。另外，微波的发射率是土壤水分反演的关键参数，在对微波地表温度反演的基础上，进一步对发射率进行了研究。     

Snow Depth Inversion Using the Localized HUT Model Based on FY-3B MWRI Data in the Farmland of Heilongjiang Province,China

Snow depth parameter inversion in the farmland using passive microwave remote sensing is of great significance to the agricultural production in Northeast China. Firstly, the Helsinki University of Technology (HUT) snow emission model was validated in the farmland based on microwave radiation imager (MWRI) onboard FengYun-3B satellite (FY-3B). The results showed that there was a big difference between the brightness temperature of HUT model simulation and MWRI for 18.7 GHz horizontal polarization (18.7 H) and 36.5 GHz horizontal polarization (36.5 H). To improve HUT model, the empirical parameter in the model was localized. Then the localized HUT (LHUT) model was built, where the extinction coefficient was calculated by the new extinction coefficient formula. Next, LHUT model was validated based on MWRI data and compared with HUT model. The results showed that LHUT underestimates slightly the brightness temperature with 0.91 and 4.19 K for 18.7 and 36.5 H respectively, and LHUT is superior to HUT model. Finally, the genetic algorithm (GA) was used to invert snow depth based on LHUT. The results showed that snow depth was underestimated with 6.79 cm based on LHUT. The inverted snow depth based on LHUT model is in better agreement with the measured snow depth.     

Qualitative and quantitative analysis of snow parameters using passive microwave remote sensing

Snow cover is an important variable for climatic and hydrologic models due to its effect on surface albedo, energy, and mass balance. Satellite observations successfully provide a global and comprehensive hemispheric-scale record of the short-term, as well as inter-seasonal variations in snow cover. Passive microwave sensors provide an excellent method to monitor temporal and spatial variations in large-scale snow cover parameters, overcoming problems of cloud cover. Using microwave remote sensing data, snow parameters (snow surface temperature, snow water equivalence, scattering index, emissivity, snow depth) have been retrieved to integrate with the snow cover simulation model developed by SASE for avalanche risk assessment on regional basis. Multispectral and multitemporal brightness temperature data obtained from the Special Sensor Microwave Imager (SSM/I), flown onboard the DMSP satellites, for the period November 2000 to April 2001 and from November 2001 to February 2002 have been analysed. A comparative data set on snow measurements and meteorological observations of a region covering large area of Pir-Panjal and the Greater Himalayan range, available on near real time basis from SASE field observatories were also used. Model calculations were carried out to study the effects of atmospheric transmission on the microwave radiation emitted from the snow covered and snow free ground and atmosphere. The sensitivity of combinations of the SSM/I channels at 19, 37 and 85 GHz, in both horizontal and vertical polarizations, in respect to snow depth, surface temperature of the snowpack have been carried out. Decision rule based algorithms are developed to identify snow cover and non-snow area.     

Cloud classification in NOAA AVHRR imageries using spectral and textural features

Clouds contribute significantly to the formation of many of the natural hazards. Hence cloud mapping and its classification becomes a major component of the various physical models which are used for forecasting natural hazards. The problem of cloud data classification from NOAA AVHRR (advance very high resolution radiometer) satellite imagery using image transformation techniques is considered in this paper. The singular value decomposition (SVD) scheme is used to extract the salient spectral and textural features attributed to satellite snow and cloud data in visible and IR channels. The goals of this paper are to discriminate between clear sky and clouds in an 8 × 8 pixel array of 1.1 km AVHRR data. If clouds are present then classify them as low, medium or high range. This scheme can effectively segregate clouds and non-cloud features in the visible and IR bands of the imagery. It can also classify clouds as low, medium or high range with a success rate of 70–90%. Computer-based snow and cloud discrimination and automatic cloud classification system will help the forecaster in various climatological applications, viz., energy balance estimation, precipitation forecasting, landslide forecasting, weather forecasting and avalanche forecasting etc.     

Snow cover area variability assessment in the upper part of the Satluj river basin in India

In high-altitude areas, snow cover plays a significant role in mountainous hydrology. Satluj, which is a snow-fed river, is a part of the Indus River system in the western Himalayas. Snow cover area (SCA) variability in this river basin affects the spatio-temporal flow availability and avalanche events. Keeping this in mind, the present study focuses on SCA variability and its relationship with various topographical features such as elevation, slope and aspect. The study has been carried out in the upper part of the Satluj River Basin on the basis of MODIS Terra (MOD10A2) data from 2001 to 2014. It has been noticed that the average annual SCA in this part of the Satluj River Basin varies from 44 to 56% with an average of about 48% of the total basin area of 16, 650 km². Further, snow accumulation and depletion curves have been suggested for assessing the SCA in the study area.     

Estimation of Snow Cover Area,Snow Physical Properties and Glacier Classification in Parts of Western Himalayas Using C-Band SAR Data

Snow physical properties, snow cover and glacier facies are important parameters which are used to quantify snowpack characteristics, glacier mass balance and seasonal snow and glacier melt. This study has been done using C-band synthetic aperture radar (SAR) data of Indian radar imaging satellite, radar imaging satellite-1 (RISAT)-1, to estimate the seasonal snow cover and retrieve snow physical properties (snow wetness and snow density), and glacier radar zones or facies classification in parts of North West Himalaya (NWH), India. Additional SAR data used are of Radarsat-2 (RS-2) satellite, which was used for glacier facies classification of Smudra Tapu glacier in Himachal Pradesh. RISAT-1 based snow cover area (SCA) mapping, snow wetness and snow density retrieval and glacier facies classification have been done for the first time in NWH region. SAR-based inversion models were used for finding out wet and dry snow dielectric constant, dry and wet SCA, snow wetness and snow density. RISAT-1 medium resolution scan-SAR mode (MRS) in HV polarization was used for first time in NWH for deriving time series of SCA maps in Beas and Bhagirathi river basins for years 2013–2014. The SAR-based inversion models were implemented separately for RISAT-1 quad pol. FRS2, for wet snow and dry snow permittivity retrieval. Masks for layover and shadow were considered in estimating final snow parameters. The overall accuracy in terms of R² value comes out to be 0.74 for snow wetness and 0.72 for snow density based on the limited ground truth data for subset area of Manali sub-basin of Beas River up to Manali for winter of 2014. Accuracy for SCA was estimated to be 95 % when compared with optical remote sensing based SCA maps with error of ±10 %. The time series data of RISAT-1 MRS and hybrid data in RH/RV mode based decompositions were also used for glacier radar zones classification for Gangotri and Samudra Tapu glaciers. The various glaciers radar zones or facies such as debris covered glacier ice, clean or bare glacier ice radar zone, percolation/refreeze radar zone and wet snow, ice wall etc., were identified. The accuracy of classified maps was estimated using ground truth data collected during 2013 and 2014 glacier field work to Samudra Tapu and Gangotri glaciers and overall accuracy was found to be in range of 82–90 %. This information of various glacier radar zones can be utilized in marking firn line of glaciers, which can be helpful for glacier mass balance studies.