高光谱遥感积雪制图算法及验证

雪盖面积是高山地区和季节雪盖区水文和气象模型的重要输入因子。机载和星载遥感数据提取的雪盖面积是融雪径流模型的重要组成部分。对应不同传感器件的光谱特征，多种分类方法被相继提出。但是，缺乏相对独立的验证手段来评价各种分类方法，其主要原因是缺乏地面真实状态。针对该现状，本研究利用高光谱图像的细分光谱特征，建立高光谱影像及其对应“地面真相”的像对数据库来发展和验证积雪制图算法，并展示MODIS积雪制图算法验证和ASTER混合像元分解雪盖制图算法研究的应用实例。     

Atmospheric correction of AMSR-E brightness temperatures for dry snow cover mapping

Differences between the brightness temperatures (spectral gradient) collected by the Advanced Microwave Scanning Radiometer for EOS (AMSR-E) at 18.7 and 36.5 GHz are used to map the snow-covered area (SCA) over a region including the western U.S. The brightness temperatures are corrected to take into account for atmospheric effects by means of a simplified radiative transfer equation whose parameters are stratified using rawinsonde data collected from a few stations. The surface emissivity is estimated from the model, and the brightness temperatures at the surface are computed as the product of the surface temperature and the computed emissivity. The SCA derived from microwave data is compared with that obtained from the Moderate Resolution Imaging Spectroradiometer for both cases of corrected and noncorrected brightness temperatures. The improvement to the SCA retrievals based on the corrected brightness temperatures shows an average value around 7%.     

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.     

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.     

Snow Cover Area Mapping Using Synthetic Aperture Radar in Manali Watershed of Beas River in the Northwest Himalayas

The current study has used Synthetic Aperture Radar (SAR) satellite data to estimate the Snow Cover Area (SCA) in Manali watershed of Beas River in Northwest Himalayas of Himachal Pradesh, India. SAR data used in this study is of Radarsat-2 (RS2) and Environmental Satellite (ENVISAT), Advanced Synthetic Aperture Radar (ASAR). The SAR preprocessing was done with SAR image processing tools for converting raw SAR images into calibrated geo-coded backscatter images. Maps for forest, built area, layover and shadow were created and used for masking snow cover in these areas. The backscattering ratio of wet snow to reference image threshold method with value range from ?2 to ?3 db was used to estimate wet SCA for study area. In this technique, if the threshold is too high (≥-2 db) wet SCA is overestimated and if it is too low (≤-3db), this method underestimates the SCA. The wet SCA is under/over estimated (+6 % to?8 % on average) in late spring season due to the inherent terrain and SAR imaging effects of layover/foreshortening and shadow and also due to the masking of forest areas. Overall, the SCA derived from SAR data matches well when compared with total SCA derived from cloud free optical remote sensing data products, especially during wet season.     

HUT模型的改进及其雪深反演

基于风云-3B（FY-3B）卫星的微波成像仪（MWRI）数据对HUT模型（Helsinki university of technology snow emission model）进行验证,结果表明,无论是18.7 GHz还是36.5 GHz水平极化亮温,HUT模型模拟亮温都与MWRI亮温存在较大的偏差。因此,本文对消光系数进行了本地化改进,得到了改进的HUT模型（IMPHUT模型）。IMPHUT模型在18.7 GHz水平极化和36.5 GHz水平极化时的模拟亮温偏差分别为-0.91 K和-4.19 K,较原始的HUT模型模拟精度（偏差分别为14.03 K和-16.33 K）有很大提高。最后,利用遗传算法进行雪深反演,基于IMPHUT模型的雪深反演（偏差为-6.79 cm）优于HUT模型和Chang算法,反演与实测雪深具有较好的一致性。     

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

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

FY-3B微波成像仪海洋数据无线电干扰识别

中国风云3号B星(FY-3B)上的微波成像仪MWRI通过5个频率(10.65 GHz,18.7 GHz,23.8 GHz,36.5 GHz和89.0 GHz)的双极化通道对地球表面进行监测。研究表明,MWRI资料的低频波段数据中存在着无线电频率干扰(RFI)现象,这些污染信号对遥感数据和反演产品质量产生极大的影响。本文尝试使用多通道回归方法和双主成分分析(DPCA)方法识别MWRI的10.65 GHz水平通道亮温海洋区域中的RFI信号。结果表明,双主成分分析法可以有效地识别出海洋上的RFI信号。微波成像仪10.65 GHz水平通道亮温数据中的RFI信号主要分布在地中海等欧洲附近海域,也存在于美国、日本、澳大利亚等近岸地区。     

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.     

基于MEMLS模型的积雪深度反演方法

利用被动微波遥感反演积雪深度一直是积雪遥感领域中的研究热点。在现有的积雪深度反演算法中,NASA算法因其简洁、易于扩展的特点,成为应用最为广泛的算法。但NASA算法存在着一定不足:首先,由于NASA算法基于线性拟合得出,在应用到其他研究区域时需要对反演公式进行重新拟合,适用范围受到一定限制;其次,由于算法中引入的19GHz与37GHz的亮温差在雪深达到一定范围时会达到饱和,因此算法会低估积雪深度。本文针对现有反演算法的不足之处,结合蚁群智能算法的特点,发展了基于蚁群算法的积雪深度反演算法;此外,针对NASA算法中存在的雪深低估问题,引入了AMSR-E10.7GHz亮温数据,对算法进行了改进。利用MEMLS模型的模拟数据与AMSR-E辐射亮温数据对算法进行实验,并采用实测数据与AMSR-E雪水当量产品对算法的反演精度进行评价。结果表明,两种积雪深度反演算法均是可行的,反演精度与现有产品相比有较为明显的改进。     

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.     

用模型和车载微波辐射仪研究多频率多角度下玉米的散射和衰减特性

在使用被动微波技术反演土壤水分的过程中,为去除植被的影响,通常采用适用于低频的τ-ω模型。为准确评估较高频率下植被的散射和衰减特性,以玉米为例,采用基于光线跟踪原理的双矩阵(Matrix-Doubling)微波辐射模型,研究不同高度的作物在C(6.925GHz)、X(10.65GHz)和Ku(18.7GHz)波段下的单散射反照率和传输率。模型模拟的亮度温度跟车载微波辐射仪的野外实测数据接近。为验证模拟的玉米自身微波辐射,在玉米地上铺设了一层铝箔屏蔽地表的辐射。通过给验证后的模型输入不同的参数,建立一个亮度温度数据库,以模拟自然状态下不同高度玉米的亮度温度。然后把模型模拟的结果,跟相同环境下τ-ω模型得到的结果按最小二乘法进行匹配,从而获得不同高度的玉米在C、X和Ku波段上等效的单散射反照率和传输率。     

Comparison of SAR-Based Snow-Covered Area Estimation Methods for the Boreal Forest Zone

Spaceborne synthetic aperture radar data have been utilized for regional-scale snow-covered area (SCA) monitoring for several years. Different methods have been developed and demonstrated for different geographical regions. A method utilizing a single reference image for SCA estimation has been shown to function well on mountainous and nonforested regions. For the boreal forest zone, a method using two reference images and a forest compensation procedure has been previously utilized. The single-reference-image method is evaluated here for the boreal forest zone, and its performance is compared with the Helsinki University of Technology (TKK) SCA method that is specifically developed for boreal forest regions. The SCA evaluations are carried out using Radarsat-1 data for the snow-melt seasons of 2004–2007. The SCA estimation accuracies for the radar-based methods are determined using optical satellite-based SCA data as reference. The results show that SCA estimation using a single reference image is usable for the boreal forest zone, although the accuracy is significantly weaker than that of the TKK-developed boreal forest-specific SCA method. The best accuracy obtained shows a root-mean-square error (rmse) of 0.176 for the single-reference-image method and an rmse of 0.123 for the TKK SCA method.      

Microwave response of seasonal snow-cover measured by using a ground-based radiometer at 6.93 and 18.7 ghz frequencies and at dual polarization

Snow cover is an important variable for climatic and hydrologic models due to its effects on surface albedo, energy, and moisture budgets. Passive microwave sensors can be used to monitor temporal and spatial variations in large-scale snow cover parameters, avoiding problems of cloud cover and polar nights. In the present study, brightness temperature values were estimated (using calibration curves) for moist snow on natural and blackbody/metal surface. T_B response on snow depth, density, SWE and angular variation from nadir were measured and found that T_B decreases with increase of snow depth and with increase of angle from nadir. Empirical relations were used to estimate emissivity, dielectric constant and dielectric loss factor. It was observed that emissivity decreases with the increase of dielectric constant. The dielectric constant and dielectric loss factor both increases with the increase of density. Experiments were performed during winter of year 2005 at Dhundi and Solang (H.P.), India, using ground based passive microwave radiometer having 6.9 and 18.7 GHz antenna frequencies at dual polarization.     

Concurrent use of active and passive microwave remote sensing data for monitoring of rice crop

Estimation of crop area, growth and phenological information is very important for monitoring of agricultural crops. However, judicious combination of spatial and temporal data from different spectral regions is necessary to meet the requirement. This study highlights the use of active microwave QuikSCAT Ku-band scatterometer and Special Sensor Microwave/Imager (SSM/I) passive microwave radiometer data to derive information on important phenological phases of rice crop. The wetness index, a weekly composite product derived using brightness temperatures from 19, 37 and 85 GHz channels of SSM/I, was used to identify the puddling period. Ku-band scatterometer data provided the signal of transplanted rice seedlings since they acts as scatterers and increases the backscattering. Dual peak nature of temporal backscatter curve around the heading stage of rice crop was observed in Ku-band. The decrease of backscatter after first peak was associated with the threshold value of 60% crop canopy cover. The symmetric (Gaussian) and asymmetric (lognormal) curve fits were attempted to derive the date of initiation of the heading phase. The temporal signature from each of these sensors was found to complement each other in crop growth monitoring. Image showing pixel-wise timings of heading stage revealed the differences exists in various parts of the study area.     

SAR interferometric coherence analysis for snow cover mapping in the western Himalayan region

Abstract

Information of snow cover (SC) over Himalayan regions is very important for regional climatological and hydrological studies. Precise monitoring of SC in the Himalayan region is essential for water supply to hydropower stations, irrigation requirements, and flood forecasting. Microwave remote sensing has all weather, day and night earth observation capability unlike optical remote sensing. In this study, spaceborne synthetic aperture radar interferometric (InSAR) coherence analysis is used to monitor SC over Himalayan rugged terrain. The feasibility of monitoring SC using synthetic aperture radar (SAR) interferometry depends on the ability to maintain coherence over InSAR pair acquisition time interval. ERS-1/2 InSAR coherence and ENVISAT ASAR InSAR coherence images are analyzed for SC mapping. Data sets of winter and of snow free months of the Himalayan region are taken for interferogram generation. Coherence images of the available data sets show maximum decorrelation in most of the area which indicates massive snowfall in the region in the winter season and melting in the summer. Area showing coherence loss due to decorrelation is mapped as a snow-covered area. The result is validated with field observations of snow depth and it is found that standing snow is inversely related to coherence in the Himalayan region.     

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.     

青藏高原冰雪覆盖变化监测的卫星立体影像方法

针对青藏高原冰雪覆盖变化监测问题,以各拉丹冬冰川为例,提出了利用ALOS立体像对提取冰雪DEM,并以同源不同时相的冰雪DEM监测冰雪覆盖变化及体量变化的方法。依据DEM纠正裁切后的影像,采用ISODATA分类方法得到了冰雪覆盖面积变化量,由不同时相网格DEM高程点数据高差与面积求得体积变化,得出了2009年与2010年冬季冰雪量大小。结果显示2010年12月比2009年12月冰雪量减少了19.728 3 km³,雪覆盖面积减少了349.691 km²。文中还列出了遇到的技术困难及有待进一步研究的问题。     

Mapping floods due to Hurricane Sandy using NPP VIIRS and ATMS data and geotagged Flickr imagery

In this study, we present an approach to estimate the extent of large-scale coastal floods caused by Hurricane Sandy using passive optical and microwave remote sensing data. The approach estimates the water fraction from coarse-resolution VIIRS and ATMS data through mixed-pixel linear decomposition. Based on the water fraction difference, using the physical characteristics of water inundation in a basin, the flood map derived from the coarse-resolution VIIRS and ATMS measurements was extrapolated to a higher spatial resolution of 30 m using topographic information. It is found that flood map derived from VIIRS shows less inundated area than the Federal Emergency Management Agency (FEMA) flood map and the ground observations. The bias was mainly caused by the time difference in observations. This is because VIIRS can only detect flood under clear conditions, while we can only find some clear-sky data around the New York area on 4 November 2012, when most flooding water already receded. Meanwhile, microwave measurements can penetrate through clouds and sense surface water bodies under clear-or-cloudy conditions. We therefore developed a new method to derive flood maps from passive microwave ATMS observations. To evaluate the flood mapping method, the corresponding ground observations and the FEMA storm surge flooding (SSF) products are used. The results show there was good agreement between our ATMS and the FEMA SSF flood areas, with a correlation of 0.95. Furthermore, we compared our results to geotagged Flickr contributions reporting flooding, and found that 95% of these Flickr reports were distributed within the ATMS-derived flood area, supporting the argument that such crowd-generated content can be valuable for remote sensing operations. Overall, the methodology presented in this paper was able to produce high-quality and high-resolution flood maps over large-scale coastal areas.     

北极一年海冰表面积雪深度遥感反演与时序分析

北极海冰表面的积雪深度是重要的地球物理变量, 是研究物质与能量平衡、计算海冰厚度的重要参数。为减小不同被动微波传感器观测数据的系统误差, 对国防气象卫星计划(Defense Meteorological Satellite Program, DMSP)F17-SSMIS与F13-SSM/I重叠期亮度温度数据进行交叉定标, 建立4个频段48个月尺度定标模型, 并与传统年尺度定标模型进行比较和优选, 在此基础上估算并分析2003-2014年北极一年海冰表面积雪深度变化。结果表明:19H、19V、22V、37V频段1~5月的月尺度模型决定系数高于传统年尺度拟合模型; 2003-2014年, 北极一年海冰表面积雪深度总体呈现下降趋势, 同时积雪深度存在明显的周期性变化, 每年7月积雪深度最小, 9月最大; 东西伯利亚海、拉普捷夫海和巴伦支海海冰表面积雪深度呈现减少的趋势。