基于ART模型的MODIS积雪反照率反演研究

积雪反照率是研究局地或全球的能量收支平衡和气候变化中的重要参数,遥感反演为积雪反照率的获取提供了便利的手段。积雪反照率大小主要取决于积雪的自身物理属性(雪粒径、形状和污染物等因子)以及天气状况,遥感反演反照率大多基于双向反射模型(BRDF),积雪BRDF模型常使用积雪辐射传输模型获得。采用考虑了雪粒径、粒子形状以及污染物影响的渐进辐射传输理论(ART)模型,建立了MODIS积雪反照率反演算法,得到了MODIS 8d合成积雪反照率产品。将此算法应用于具有均一积雪地表的格陵兰岛地区,并使用GC-Net实测数据进行了验证,反演的总均方根误差(RMSE)为0.018,相关系数(r)为0.83,结果表明考虑了积雪特性的ART模型能够较好地反演积雪反照率,而且反演需要的参数较少。     

Spatial relationships between snow contaminant content, grain size, and surface temperature from multispectral images of Mt. Rainier, Washington (USA)

We use multispectral MODIS/ASTER Airborne Simulator (MASTER) data collected at Mt. Rainier, Washington (USA) to map spatial covariance between snowpack properties and to evaluate techniques for quantitative estimation of reflectance, grain size, and temperature. The late-August MASTER images reveal a distinct pattern of snow contaminant content, grain size, and temperature related to a recent snowfall and late-summer melting. Spatial correlation between grain size and temperature patterns suggests that rapid destructive metamorphism of the fresh snow occurred when temperatures were near 0 °C. We use 10 specific locations to evaluate hemispherical-directional reflectance factor (HDRF), grain size, and temperature retrievals. We map relative snow contaminant content using visible (0.4-0.8 μm) HDRF spectra. Atmospheric correction and topographic modeling limit the accuracy of HDRF estimates. We use MASTER-derived spectra near 1.8 and 2.2 μm to estimate optical grain size (by comparison to modeled layers of ice spheres) and physical grain size (by comparison to measured spectra with known physical grain size and by correlation to ground measurements). Estimated physical grain sizes were less than estimated optical grain sizes. Differing definitions of optical and physical grain sizes could contribute to this discrepancy. Limitations at 1.8 and 2.2 μm, including reduced discrimination between larger grain radii (>∼500 μm physical, >∼200 μm optical) and low signal-to-noise ration with atmospheric effects and decreasing solar irradiance, suggest that grain size retrieval may be improved at other wavelengths (e.g., 1.1 μm). Accounting for uncertainty in emissivity, atmospheric correction, and detector noise, we estimate systematic errors in our radiant temperatures at <1.8 °C. This study shows both strengths and limitations for coregistered visible, short-wave infrared, and thermal infrared images to estimate snowpack properties and reveal their spatial coherence.     

MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland

In situ measurements of snow albedo at five stations along a north-south transect in the dry-snow facies of the interior of Greenland follow the theoretically expected dependence of snow albedo with solar zenith angle (SZA). Greenland Climate Network (GC-Net) measurements from 1997 through 2007 exhibit the trend of modest surface brightening with increasing SZA on both diurnal and seasonal timescales. SZA explains up to 50% of seasonal albedo variability. The two other environmental factors considered, temperature and cloudiness, play much less significant roles in seasonal albedo variability at the five stations studied. Compared to the 10-year record of these GC-Net measurements, the five-year record of MODIS satellite-retrieved snow albedo shows a systematic negative bias for SZA larger than about 55°. Larger bias of MODIS snow albedo exists at more northerly stations. MODIS albedos successfully capture the snow albedo dependence on SZA and have a relatively good agreement with GC-Net measurements for SZA < 55°. The discrepancy of MODIS albedo with in situ albedo and with theory is determined mainly by two related factors, SZA and retrieval quality. While the spatiotemporal structure, especially zonal features, of the MODIS-retrieved albedo may be correct for large SZA, the accuracy deteriorates for SZA > 55° and often becomes physically unrealistic for SZA > 65°. This unphysical behavior biases parameterizations of surface albedo and restricts the range of usefulness of the MODIS albedo products. Seasonal-to-interannual trends in surface brightness in Greenland, and in polar (i.e., large SZA) regions in general, and model simulations of these trends, should be evaluated in light of these limitations.     

Accuracy assessment of the MODIS 16-day albedo product for snow: comparisons with Greenland in situ measurements

The accuracy of the Moderate Resolution Imaging Spectroradiometer (MODIS) 16-day albedo product (MOD43) is assessed using ground-based albedo observations from automatic weather stations (AWS) over spatially homogeneous snow and semihomogeneous ice-covered surfaces on the Greenland ice sheet. Data from 16 AWS locations, spanning the years 2000-2003, were used for this assessment. In situ reflected shortwave data were corrected for a systematic positive spectral sensitivity bias of between 0.01 and 0.09 on a site-by-site basis using precise optical black radiometer data. Results indicate that the MOD43 albedo product retrieves snow albedo with an average root mean square error (RMSE) of ±0.07 as compared to the station measurements, which have ±0.035 RMSE uncertainty. If we eliminate all satellite retrievals that rely on the backup algorithm and consider only the highest quality results from the primary bidirectional reflectance distribution function (BRDF) algorithm, the MODIS albedo RMSE is ±0.04, slightly larger than the in situ measurement uncertainty. There is general agreement between MODIS and in situ observations for albedo <0.7, while near the upper limit, a −0.05 MODIS albedo bias is evident from the scatter of the 16-site composite.     

Commentary on Wang and Zender—MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland

The recent paper by Wang and Zender [Wang, X., & Zender, C. S. (2010). MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland. Remote Sensing of Environment.] draws erroneous conclusions about solar zenith angle biases at high latitudes by not making appropriate use of the extensive quality flags available with the MODIS BRDF/Albedo. Coarse resolution MODIS white-sky albedo data are compared with actual blue-sky field albedometer measurements from the Greenland GC-Net. By utilizing large area averages of the MODIS data product that combine both high quality and poor quality data indiscriminately, the authors erroneously conclude that the accuracy deteriorates for solar zenith angle (SZA) > 55° and often becomes physically unrealistic for SZA > 65°. Once the quality flags are considered, however, the comparisons demonstrate that the MODIS product performs quite well out to the recommended limit for product use of 70° SZA. This verifies the conclusions of an earlier more rigorous evaluation performed by Stroeve et al. [Stroeve, J., Box, J. E., Gao, F., Liang, S., Nolin, A., & Schaaf, C. B. (2005). Accuracy assessment of the MODIS 16-day albedo product for snow: comparisons with Greenland in situ measurements. Remote Sensing of Environment.]. With over a decade of observations and products now available from the MODIS instrument, these data are increasingly being used to evaluate and tune climate and biogeochemical models. However, such use should take into account the documented quality and limitations of the satellite-derived product.     

Snow optical properties for different particle shapes with application to snow grain size retrieval and MODIS/CERES radiance comparison over Antarctica

We investigated the single scattering optical properties of snow for different ice particle shapes and degrees of microscopic scale roughness. These optical properties were implemented and tested in a coupled atmosphere-snow radiative transfer model. The modeled surface spectral albedo and radiance distribution were compared with surface measurements. The results show that the reflected radiance and irradiance over snow are sensitive to the snow grain size and its vertical profile. When inhomogeneity of the particle size distribution in the vertical is taken into account, the measured spectral albedo can be matched, regardless of the particle shapes. But this is not true for the modeled radiance distribution, which depends a lot on the particle shape. The usual “equivalent spheres” assumption significantly overestimates forward reflected radiances, and underestimates backscattering radiances, around the principal plane. On average, the aggregate shape assumption has the best agreement with the measured radiances to a mean bias within 2%.The snow optical properties with the aggregate assumption were applied to the retrieval of snow grain size over the Antarctic plateau. The retrieved grain sizes of the top layer showed similar and large seasonal variation in all years, but only small year to year variation. Using the retrieved snow grain sizes, the modeled spectral and broadband radiances showed good agreements with MODIS and CERES measurements over the Antarctic plateau. Except for the MODIS 2.13 μm channel, the mean relative model-observation differences are within few percent. The modeled MODIS radiances using measured surface reflectance at Dome C also showed good agreement in visible channels, where radiation is not sensitive to snow grain size and the measured surface bidirectional reflectance is applicable over the Antarctic plateau. But modeled radiances using local, surface-measured reflectance in the near infrared yielded large errors because of the high sensitivity to the snow grain size, which varies spatially and temporally. The CERES broadband shortwave radiance is moderately sensitive to the snow grain size, comparable to the MODIS 0.86 μm channel. The variation of broadband snow reflectance due to the seasonal variation in snow grain size is about 5% in a year over the Antarctic plateau. CERES broadband radiances simulated with grain sizes retrieved using MODIS are about 2% larger than those observed.     

Validation of the Climate-SAF surface broadband albedo product: Comparisons with in situ observations over Greenland and the ice-covered Arctic Ocean

This paper describes a validation study performed by comparing the Climate-SAF Surface Albedo Product (SAL) to ground truth observations over Greenland and the ice-covered Arctic Ocean. We compare Advanced Very High Resolution Radiometer (AVHRR)-based albedo retrievals to data from the Greenland Climate Network (GCN) weather stations and the floating ice station Tara for polar summer 2007. The AVHRR dataset consists of 2755 overpasses. The overpasses are matched to in situ observations spatially and temporally. The SAL algorithm presented here derives the surface broadband albedo from AVHRR channels 1 and 2 using an atmospheric correction, temporal sampling of an empirical Bidirectional Reflectance Distribution Function (BRDF), and a narrow-to-broadband conversion algorithm. The satellite product contains algorithms for snow, sea ice, vegetation, bare soil, and water albedo. At the Summit and DYE-2 stations on the Greenland ice sheet, instantaneous SAL RMSE is 0.073. The heterogeneous surface conditions at satellite pixel scale over the stations near the Greenland west coast increase RMSE to > 0.12. Over Tara, the instantaneous SAL RMSE is 0.069. The BRDF sampling approach reduces RMSE over the ice sheet to 0.053, and to 0.045 over Tara. Taking into account various sources of uncertainty for both satellite retrievals and in situ observations, we conclude that SAL agrees with in situ observations within their limits of accuracy and spatial representativeness.     

Retrieval of subpixel snow covered area, grain size, and albedo from MODIS

We describe and validate a model that retrieves fractional snow-covered area and the grain size and albedo of that snow from surface reflectance data (product MOD09GA) acquired by NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). The model analyzes the MODIS visible, near infrared, and shortwave infrared bands with multiple endmember spectral mixtures from a library of snow, vegetation, rock, and soil. We derive snow spectral endmembers of varying grain size from a radiative transfer model specific to a scene's illumination geometry; spectra for vegetation, rock, and soil were collected in the field and laboratory. We validate the model with fractional snow cover estimates from Landsat Thematic Mapper data, at 30 m resolution, for the Sierra Nevada, Rocky Mountains, high plains of Colorado, and Himalaya. Grain size measurements are validated with field measurements during the Cold Land Processes Experiment, and albedo retrievals are validated with in situ measurements in the San Juan Mountains of Colorado. The pixel-weighted average RMS error for snow-covered area across 31 scenes is 5%, ranging from 1% to 13%. The mean absolute error for grain size was 51 µm and the mean absolute error for albedo was 4.2%. Fractional snow cover errors are relatively insensitive to solar zenith angle. Because MODSCAG is a physically based algorithm that accounts for the spatial and temporal variation in surface reflectances of snow and other surfaces, it is capable of global snow cover mapping in its more computationally efficient, operational mode.     

In-situ measured spectral directional emissivity of snow and ice in the 8–14 μm atmospheric window

The directional emissivity of snow and ice surfaces in the 8–14 μm thermal infrared (TIR) atmospheric window was determined from spectral radiances obtained by field measurements using a portable Fourier transform infrared spectrometer in conjunction with snow pit work. The dependence of the directional emissivity on the surface snow type (grain size and shape) was examined. We obtained emissivity spectra for five different surface types, i.e., fine dendrite snow, medium granular snow, coarse grain snow, welded sun crust snow, and smooth bare ice. The derived emissivities show a distinct spectral contrast at wavelengths λ = 10.5–12.5 μm which is enhanced with increasing the snow grain size. For example, emissivities at both 10.5 μm and 12.5 μm for the nadir angle were 0.997 and 0.984 for the fine dendrite snow, 0.996 and 0.974 for the medium granular snow, 0.995 and 0.971 for the coarse grain snow, 0.992 and 0.968 for the sun crust, and 0.993 and 0.949 for the bare ice, respectively. In addition, the spectral contrast exhibits a strong angular dependence, particularly for the coarser snow and bare ice, e.g., the emissivity at λ = 12.5 μm for the off-nadir angle of 75° reaches down to 0.927, 0.896, and 0.709 for the coarse grain snow, sun crust, and bare ice cases, respectively. The angular dependent emissivity spectra of the bare ice were quite consistent with the spectra predicted by the Fresnel reflectance theory. The observed results firmly demonstrate that the directional emissivity of snow in the TIR can vary depending upon the surface snow type. The high variability of the spectral emissivity of snow also suggests the possibility to discriminate between snow and ice types from space using the brightness temperature difference in the atmospheric window.     

Algorithm for retrieval of the effective snow grain size and pollution amount from satellite measurements

We present an algorithm for retrieval of the effective Snow Grain Size and Pollution amount (SGSP) from satellite measurements. As well as our previous version (Zege et al., 2008, 1998), the new algorithm is based on the analytical solution for snow reflectance within the asymptotic radiative transfer theory. The SGSP algorithm does not use any assumptions on snow grain shape and allows for the snow pack bidirectional reflectance distribution function (BRDF). The algorithm includes a new atmospheric correction procedure that allows for snow BRDF. This SGSP algorithm has been thoroughly validated with computer simulations. Its sensitivity to the atmosphere model has been investigated. It is shown that the inaccuracy of the snow characteristic retrieval due to the uncertainty in the aerosol and molecular atmosphere model is negligible, as compared to that due to the measurement errors at least for aerosol loads typical for polar regions. The significant advantage of the SGSP over conventional algorithms, which use a priori assumptions about particle shape and (or) not allow for the BRDF of the individual snow pack, is that the developed retrieval still works at low sun elevations, which are typical for polar regions.     

积雪属性非均匀性对被动微波遥感积雪的影响

积雪属性的非均匀性在水平方向上表现为像元内积雪未完全覆盖和雪深分布的不均匀,在垂直方向上表现为积雪剖面上粒径和密度的不一致导致的积雪分层现象。这些积雪属性的非均匀性对被动微波遥感反演雪深或雪水当量带来很大的不确定性,并且给反演结果的验证带来不确定性。通过野外积雪的微波辐射特性观测、遥感积雪产品对比分析、积雪辐射传输模型模拟对这些问题进行阐述和探讨,为今后积雪微波遥感反演算法发展和结果验证提供参考。     

Semi-empirical approach for estimating broadband albedo of snow

Snow is a medium that exhibits highly anisotropic reflectance throughout the solar spectrum. The anisotropic nature of snow shows more variability in snow metamorphic processes for wavelengths beyond 1.0 μm than in the visible spectrum. This behavior poses challenges for the development of a model that can retrieve broadband albedo from reflectance measurements throughout the snow season. In this paper, a semi-empirical model is presented to estimate near infrared (0.8-2.5 μm) albedo of snow. This model estimates spectral albedo at a wavelength of 1.240 μm using only three variables: solar zenith angle, scattering angle and measured reflectance, which is used to retrieve near infrared albedo. To form a base for such a model, quantification of reflectance patterns and variability in varying snow condition, i.e. snow grain size, and sun-sensor geometry are prerequisite. In this study the DIScrete Ordinate Radiative Transfer (DISORT) model is used to simulate bi-directional reflectance. The performance of the developed model is evaluated by using DISORT simulated spectral albedo for various snow grain sizes and solar zenith angles, as well as the Moderate Resolution Imaging Spectroradiometer (MODIS) and in-situ measurements. The developed model is shown to be capable of estimating spectral albedo at 1.240 μm with acceptable accuracy. The mean error (ME), mean absolute error (MAE), and root mean squared error (RMSE) in the estimates are found to be 0.053, 0.055 and 0.064, respectively, for a wide range of sun-sensor geometries and snow grain sizes. The model shows better accuracy for spectral albedo estimates than for those computed using the Lambertian reflectance assumption for snow, reducing the error in the range and standard deviation by 75% and 65%, respectively. Applying the model to MODIS, the retrieved albedo is found to be in good quantitative agreement (r = 0.82) with in-situ measurements. These improvements in albedo estimation should allow more accurate use of remote sensing measurements in climate and hydrological models.     

Sizing snow grains using backscattered solar light

In this article, we describe a technique to determine dry snow grain size from optical observations. The method is based on analysis of the snow reflectance in the near-infrared region, in particular, the Medium Resolution Imaging Spectrometer (MERIS) band at 865 nm, which is common to many spaceborne optical sensors, is used. In addition, the algorithm is applied to the Moderate Resolution Imaging Spectroradiometer (MODIS) 1240 nm band. It is found that bands located at 1020 and 1240 nm are the most suitable for snow grain size remote-sensing applications. The developed method is validated using MODIS observations over flat snow deposited on a lake ice in Hokkaido, Japan.     

Estimating northern hemisphere snow water equivalent for climate research through assimilation of space-borne radiometer data and ground-based measurements

The key variable describing global seasonal snow cover is snow water equivalent (SWE). However, reliable information on the hemispheric scale variability of SWE is lacking because traditional methods such as interpolation of ground-based measurements and stand-alone algorithms applied to space-borne observations are highly uncertain with respect to the spatial distribution of snow mass and its evolution. In this paper, an algorithm assimilating synoptic weather station data on snow depth with satellite passive microwave radiometer data is applied to produce a 30-year-long time-series of seasonal SWE for the northern hemisphere. This data set is validated using independent SWE reference data from Russia, the former Soviet Union, Finland and Canada. The validation of SWE time-series indicates overall strong retrieval performance with root mean square errors below 40 mm for cases when SWE < 150 mm. Retrieval uncertainty increases when SWE is above this threshold. The SWE estimates are also compared with results obtained by a typical stand-alone satellite passive microwave algorithm. This comparison demonstrates the benefits of the newly developed assimilation approach. Additionally, the trends and inter-annual variability of northern hemisphere snow mass during the era of satellite passive microwave measurements are shown.     

Evaluation of the HUT modified snow emission model over lake ice using airborne passive microwave measurements

The algorithms designed to estimate snow water equivalent (SWE) using passive microwave measurements falter in lake-rich high-latitude environments due to the emission properties of ice covered lakes on low frequency measurements. Microwave emission models have been used to simulate brightness temperatures (Tbs) for snowpack characteristics in terrestrial environments but cannot be applied to snow on lakes because of the differing subsurface emissivities and scattering matrices present in ice. This paper examines the performance of a modified version of the Helsinki University of Technology (HUT) snow emission model that incorporates microwave emission from lake ice and sub-ice water. Inputs to the HUT model include measurements collected over brackish and freshwater lakes north of Inuvik, Northwest Territories, Canada in April 2008, consisting of snowpack (depth, density, and snow water equivalent) and lake ice (thickness and ice type). Coincident airborne radiometer measurements at a resolution of 80 × 100 m were used as ground-truth to evaluate the simulations.The results indicate that subsurface media are simulated best when utilizing a modeled effective grain size and a 1 mm RMS surface roughness at the ice/water interface compared to using measured grain size and a flat Fresnel reflective surface as input. Simulations at 37 GHz (vertical polarization) produce the best results compared to airborne Tbs, with a Root Mean Square Error (RMSE) of 6.2 K and 7.9 K, as well as Mean Bias Errors (MBEs) of −8.4 K and −8.8 K for brackish and freshwater sites respectively. Freshwater simulations at 6.9 and 19 GHz H exhibited low RMSE (10.53 and 6.15 K respectively) and MBE (−5.37 and 8.36 K respectively) but did not accurately simulate Tb variability (R = −0.15 and 0.01 respectively). Over brackish water, 6.9 GHz simulations had poor agreement with airborne Tbs, while 19 GHz V exhibited a low RMSE (6.15 K), MBE (−4.52 K) and improved relative agreement to airborne measurements (R = 0.47). Salinity considerations reduced 6.9 GHz errors substantially, with a drop in RMSE from 51.48 K and 57.18 K for H and V polarizations respectively, to 26.2 K and 31.6 K, although Tb variability was not well simulated. With best results at 37 GHz, HUT simulations exhibit the potential to track Tb evolution, and therefore SWE through the winter season.     

Fractional snow cover mapping through artificial neural network analysis of MODIS surface reflectance

Accurate areal measurements of snow cover extent are important for hydrological and climate modeling. The traditional method of mapping snow cover is binary where a pixel is considered either snow-covered or snow-free. Fractional snow cover (FSC) mapping can achieve a more precise estimate of areal snow cover extent by estimating the fraction of a pixel that is snow-covered. The most common snow fraction methods applied to Moderate Resolution Imaging Spectroradiometer (MODIS) images have been spectral unmixing and an empirical Normalized Difference Snow Index (NDSI). Machine learning is an alternative for estimating FSC as artificial neural networks (ANNs) have been successfully used for estimating the subpixel abundances of other surfaces. The advantages of ANNs are that they can easily incorporate auxiliary information such as land cover type and are capable of learning nonlinear relationships between surface reflectance and snow fraction. ANNs are especially applicable to mapping snow cover extent in forested areas where spatial mixing of surface components is nonlinear. This study developed a multilayer feed-forward ANN trained through backpropagation to estimate FSC using MODIS surface reflectance, NDSI, Normalized Difference Vegetation Index (NDVI) and land cover as inputs. The ANN was trained and validated with higher spatial-resolution FSC maps derived from Landsat Enhanced Thematic Mapper Plus (ETM+) binary snow cover maps. Testing of the network was accomplished over training and independent test areas. The developed network performed adequately with RMSE of 12% over training areas and slightly less accurately over the independent test scenes with RMSE of 14%. The developed ANN also compared favorably to the standard MODIS FSC product. The study also presents a comprehensive validation of the standard MODIS snow fraction product whose performance was found to be similar to that of the ANN.     

Comparison of satellite-derived and in-situ observations of ice and snow surface temperatures over Greenland

The most practical way to get spatially broad and continuous measurements of the surface temperature in the data-sparse cryosphere is by satellite remote sensing. The uncertainties in satellite-derived LSTs must be understood to develop internally-consistent decade-scale land surface temperature (LST) records needed for climate studies. In this work we assess satellite-derived “clear-sky” LST products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), and LSTs derived from the Enhanced Thematic Mapper Plus (ETM+) over snow and ice on Greenland. When possible, we compare satellite-derived LSTs with in-situ air temperature observations from Greenland Climate Network (GC-Net) automatic weather stations (AWS). We find that MODIS, ASTER and ETM+ provide reliable and consistent LSTs under clear-sky conditions and relatively-flat terrain over snow and ice targets over a range of temperatures from ? 40 to 0 °C. The satellite-derived LSTs agree within a relative RMS uncertainty of ~ 0.5 °C. The good agreement among the LSTs derived from the various satellite instruments is especially notable since different spectral channels and different retrieval algorithms are used to calculate LST from the raw satellite data. The AWS record in-situ data at a “point” while the satellite instruments record data over an area varying in size from: 57 × 57 m (ETM+), 90 × 90 m (ASTER), or to 1 × 1 km (MODIS). Surface topography and other factors contribute to variability of LST within a pixel, thus the AWS measurements may not be representative of the LST of the pixel. Without more information on the local spatial patterns of LST, the AWS LST cannot be considered valid ground truth for the satellite measurements, with RMS uncertainty ~ 2 °C. Despite the relatively large AWS-derived uncertainty, we find LST data are characterized by high accuracy but have uncertain absolute precision.     

Development of a winter snow water equivalent algorithm using in situ passive microwave radiometry over snow-covered first-year sea ice

A snow water equivalent (SWE) algorithm has been developed for thin and thick snow using both in situ microwave measurements and snow thermophysical properties, collected over landfast snow covered first-year sea ice during the Canadian Arctic Shelf Exchange Study (CASES) overwintering mission from December 2003 to May 2004. Results showed that the behavior of brightness temperatures (T_bs) in thin snow covers was very different from those in a thick snowpack. Microwave SWE retrievals using the combination of T_b 19 GHz and air temperature (multiple regression) over thick snow are quite accurate, and showed very good agreement with the physical data (R² = 0.94) especially during the cooling period (i.e., from freeze up to the minimum air temperature recorded) where the snow is dry and cold. Thin snow SWE predictions also showed fairly good agreement with field data (R² = 0.70) during the cold season. The differences between retrieved and in situ SWE for both thin and thick snow cover are mainly attributable to the variations in air temperature, snow wetness and spatial heterogeneity in snow thickness.     

被动微波遥感在青藏高原积雪业务监测中的初步应用

积雪范围、积雪深度和雪水当量等参数的遥感监测与反演对气候模式的建立以及积雪灾害的评估具有重要意义。被动微波遥感在这些参数的反演方面具有明显优势,但目前尚未应用到青藏高原地区的积雪遥感业务监测上来。2001年10月至2002年4月,利用SSM/I数据对青藏高原地区的积雪范围和积雪深度进行了实时监测,为西藏、青海遥感应用部门提供逐日的雪深分布图。对这次监测的总效果进行了分析和评价,并对发生在青海省内一次较大的降雪过程进行了遥感分析,结果表明:SSM/I反演的积雪范围变化趋势与MODIS结果总体上较为一致;SSM/I的雪深监测结果为当地遥感部门对大于10 cm的雪深做出正确判断提供了重要信息,是对雪灾定位的重要信息源。     

A satellite snow depth multi-year average derived from SSM/I for the high latitude regions

The hydrological cycle for high latitude regions is inherently linked with the seasonal snowpack. Thus, accurately monitoring the snow depth and the associated aerial coverage are critical issues for monitoring the global climate system. Passive microwave satellite measurements provide an optimal means to monitor the snowpack over the arctic region. While the temporal evolution of snow extent can be observed globally from microwave radiometers, the determination of the corresponding snow depth is more difficult. A dynamic algorithm that accounts for the dependence of the microwave scattering on the snow grain size has been developed to estimate snow depth from Special Sensor Microwave/Imager (SSM/I) brightness temperatures and was validated over the U.S. Great Plains and Western Siberia.

The purpose of this study is to assess the dynamic algorithm performance over the entire high latitude (land) region by computing a snow depth multi-year field for the time period 1987–1995. This multi-year average is compared to the Global Soil Wetness Project-Phase2 (GSWP2) snow depth computed from several state-of-the-art land surface schemes and averaged over the same time period. The multi-year average obtained by the dynamic algorithm is in good agreement with the GSWP2 snow depth field (the correlation coefficient for January is 0.55). The static algorithm, which assumes a constant snow grain size in space and time does not correlate with the GSWP2 snow depth field (the correlation coefficient with GSWP2 data for January is − 0.03), but exhibits a very high anti-correlation with the NCEP average January air temperature field (correlation coefficient − 0.77), the deepest satellite snow pack being located in the coldest regions, where the snow grain size may be significantly larger than the average value used in the static algorithm. The dynamic algorithm performs better over Eurasia (with a correlation coefficient with GSWP2 snow depth equal to 0.65) than over North America (where the correlation coefficient decreases to 0.29).