A vegetation index based technique for spatial sharpening of thermal imagery

High spatial resolution (∼ 100 m) thermal infrared band imagery has utility in a variety of applications in environmental monitoring. However, currently such data have limited availability and only at low temporal resolution, while coarser resolution thermal data (∼ 1000 m) are routinely available, but not as useful for identifying environmental features for many landscapes. An algorithm for sharpening thermal imagery (TsHARP) to higher resolutions typically associated with the shorter wavebands (visible and near-infrared) used to compute vegetation indices is examined over an extensive corn/soybean production area in central Iowa during a period of rapid crop growth. This algorithm is based on the assumption that a unique relationship between radiometric surface temperature (T_R) relationship and vegetation index (VI) exists at multiple resolutions. Four different methods for defining a VI − T_R basis function for sharpening were examined, and an optimal form involving a transformation to fractional vegetation cover was identified. The accuracy of the high-resolution temperature retrieval was evaluated using aircraft and Landsat thermal imagery, aggregated to simulate native and target resolutions associated with Landsat, MODIS, and GOES short- and longwave datasets. Applying TsHARP to simulated MODIS thermal maps at 1-km resolution and sharpening down to ∼ 250 m (MODIS VI resolution) yielded root-mean-square errors (RMSE) of 0.67-1.35 °C compared to the ‘observed’ temperature fields, directly aggregated to 250 m. Sharpening simulated Landsat thermal maps (60 and 120 m) to Landsat VI resolution (30 m) yielded errors of 1.8-2.4 °C, while sharpening simulated GOES thermal maps from 5 km to 1 km and 250 m yielded RMSEs of 0.98 and 1.97, respectively. These results demonstrate the potential for improving the spatial resolution of thermal-band satellite imagery over this type of rainfed agricultural region. By combining GOES thermal data with shortwave VI data from polar orbiters, thermal imagery with 250-m spatial resolution and 15-min temporal resolution can be generated with reasonable accuracy. Further research is required to examine the performance of TsHARP over regions with different climatic and land-use characteristics at local and regional scales.     

Estimating fractional snow cover from MODIS using the normalized difference snow index

Snow-cover information is important for a wide variety of scientific studies, water supply and management applications. The NASA Earth Observing System (EOS) Moderate Resolution Imaging Spectroradiometer (MODIS) provides improved capabilities to observe snow cover from space and has been successfully using a normalized difference snow index (NDSI), along with threshold tests, to provide global, automated binary maps of snow cover. The NDSI is a spectral band ratio that takes advantage of the spectral differences of snow in short-wave infrared and visible MODIS spectral bands to identify snow versus other features in a scene. This study has evaluated whether there is a “signal” in the NDSI that could be used to estimate the fraction of snow within a 500 m MODIS pixel and thereby enhance the use of the NDSI approach in monitoring snow cover. Using Landsat 30-m observations as “ground truth,” the percentage of snow cover was calculated for 500-m cells. Then a regression relationship between 500-m NDSI observations and fractional snow cover was developed over three different snow-covered regions and tested over other areas. The overall results indicate that the relationship between fractional snow cover and NDSI is reasonably robust when applied locally and over large areas like North America. The relationship offers advantages relative to other published fractional snow cover algorithms developed for global-scale use with MODIS. This study indicates that the fraction of snow cover within a MODIS pixel using this approach can be provided with a mean absolute error less than 0.1 over the range from 0.0 to 1.0 in fractional snow cover.     

A comparison of multi-spectral, multi-angular, and multi-temporal remote sensing datasets for fractional shrub canopy mapping in Arctic Alaska

Shrub cover appears to be increasing across many areas of the Arctic tundra biome, and increasing shrub cover in the Arctic has the potential to significantly impact global carbon budgets and the global climate system. For most of the Arctic, however, there is no existing baseline inventory of shrub canopy cover, as existing maps of Arctic vegetation provide little information about the density of shrub cover at a moderate spatial resolution across the region. Remotely-sensed fractional shrub canopy maps can provide this necessary baseline inventory of shrub cover. In this study, we compare the accuracy of fractional shrub canopy (> 0.5 m tall) maps derived from multi-spectral, multi-angular, and multi-temporal datasets from Landsat imagery at 30 m spatial resolution, Moderate Resolution Imaging SpectroRadiometer (MODIS) imagery at 250 m and 500 m spatial resolution, and MultiAngle Imaging Spectroradiometer (MISR) imagery at 275 m spatial resolution for a 1067 km² study area in Arctic Alaska. The study area is centered at 69 °N, ranges in elevation from 130 to 770 m, is composed primarily of rolling topography with gentle slopes less than 10°, and is free of glaciers and perennial snow cover. Shrubs > 0.5 m in height cover 2.9% of the study area and are primarily confined to patches associated with specific landscape features. Reference fractional shrub canopy is determined from in situ shrub canopy measurements and a high spatial resolution IKONOS image swath. Regression tree models are constructed to estimate fractional canopy cover at 250 m using different combinations of input data from Landsat, MODIS, and MISR. Results indicate that multi-spectral data provide substantially more accurate estimates of fractional shrub canopy cover than multi-angular or multi-temporal data. Higher spatial resolution datasets also provide more accurate estimates of fractional shrub canopy cover (aggregated to moderate spatial resolutions) than lower spatial resolution datasets, an expected result for a study area where most shrub cover is concentrated in narrow patches associated with rivers, drainages, and slopes. Including the middle infrared bands available from Landsat and MODIS in the regression tree models (in addition to the four standard visible and near-infrared spectral bands) typically results in a slight boost in accuracy. Including the multi-angular red band data available from MISR in the regression tree models, however, typically boosts accuracy more substantially, resulting in moderate resolution fractional shrub canopy estimates approaching the accuracy of estimates derived from the much higher spatial resolution Landsat sensor. Given the poor availability of snow and cloud-free Landsat scenes in many areas of the Arctic and the promising results demonstrated here by the MISR sensor, MISR may be the best choice for large area fractional shrub canopy mapping in the Alaskan Arctic for the period 2000-2009.     

Toward improved daily snow cover mapping with advanced combination of MODIS and AMSR-E measurements

Taking three snow seasons from November 1 to March 31 of year 2002 to 2005 in northern Xinjiang, China as an example, this study develops a new daily snow cover product (500 m) through combining MODIS daily snow cover data and AMSR-E daily snow water equivalent (SWE) data. By taking advantage of both high spatial resolution of optical data and cloud transparency of passive microwave data, the new daily snow cover product greatly complements the deficiency of MODIS product when cloud cover is present especially for snow cover product on a daily basis and effectively improves daily snow detection accuracy. In our example, the daily snow agreement of the new product with the in situ measurements at 20 stations is 75.4%, which is much higher than the 33.7% of the MODIS daily product in all weather conditions, even a little higher than the 71% of the MODIS 8-day product (cloud cover of ~ 5%). Our results also indicate that i) AMSR-E daily SWE imagery generally agrees with MOD10A1 data in detecting snow cover, with overall agreement of 93.4% and snow agreement of 96.6% in the study area; ii) AMSR-E daily SWE imagery underestimates the snow covered area (SCA) due to its coarse spatial resolution; iii) The new snow cover product can better and effectively capture daily SCA dynamics during the snow seasons, which plays a significant role in reduction, mitigation, and prevention of snow-caused disasters in pastoral areas.     

Extending the MODIS 1 km ocean colour atmospheric correction to the MODIS 500 m bands and 500 m chlorophyll-a estimation towards coastal and estuarine monitoring

National and regional obligations to control and maintain water quality have led to an increase in coastal and estuarine monitoring. A potentially valuable tool is high temporal and spatial resolution satellite ocean colour data. NASA's MODIS-Terra and -Aqua can capture data at both 250 m and 500 m spatial resolutions and the existence of two sensors provides the possibility for multiple daily passes over a scene. However, no robust atmospheric correction method currently exists for these data, rendering them unusable for quantitative monitoring applications. Therefore, this paper presents an automatic and dynamic atmospheric correction approach allowing the determination of ocean colour. The algorithm is based around the standard MODIS 1 km atmospheric correction, includes cloud masking and is applicable to all of the visible 500 m bands. Comparison of the 500 m ocean colour data with the standard 1 km data shows good agreement and these results are further supported by in situ data comparisons. In addition, a novel method to produce 500 m chlorophyll-a estimates is presented. Comparisons of the 500 m estimates with the standard MODIS OC3M algorithm and to in situ data from a near-coast validation site are given.     

Monitoring snow cover variability in an agropastoral area in the Trans Himalayan region of Nepal using MODIS data with improved cloud removal methodology

Monitoring the extent and pattern of snow cover in the dry, high altitude, Trans Himalayan region (THR) is significant to understand the local and regional impact of ongoing climate change and variability. The freely available Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover images, with 500 m spatial and daily temporal resolution, can provide a basis for regional snow cover mapping, monitoring and hydrological modelling. However, high cloud obscuration remains the main limitation. In this study, we propose a five successive step approach — combining data from the Terra and Aqua satellites; adjacent temporal deduction; spatial filtering based on orthogonal neighbouring pixels; spatial filtering based on a zonal snowline approach; and temporal filtering based on zonal snow cycle — to remove cloud obscuration from MODIS daily snow products. This study also examines the spatial and temporal variability of snow cover in the THR of Nepal in the last decade. Since no ground stations measuring snow data are available in the region, the performance of the proposed methodology is evaluated by comparing the original MODIS snow cover data with least cloud cover against cloud-generated MODIS snow cover data, filled by clouds of another densely cloud-covered product. The analysis indicates that the proposed five-step method is efficient in cloud reduction (with average accuracy of > 91%). The results show very high interannual and intra-seasonal variability of average snow cover, maximum snow extent and snow cover duration over the last decade. The peak snow period has been delayed by about 6.7 days per year and the main agropastoral production areas of the region were found to experience a significant decline in snow cover duration during the last decade.     

面向亚像元级积雪信息提取的MODIS影像波段融合方法研究

MODIS影像因其共享性和时间序列的完整性而成为大区域积雪监测研究广泛使用的数据源,进行MODIS影像波段间融合,能够为积雪研究提供较高分辨率的影像数据源。为了充分利用MODIS影像250 m分辨率波段的空间和光谱信息,提取亚像元级的积雪面积,使用两种具有高光谱保真度的影像融合方法：基于SFIM变换和基于小波变换的融合方法,采取不同的波段组合策略,对MODIS影像bands 1~2和bands 3~7进行融合,并以Landsat TM影像的积雪分类图作为“真值”,对融合后影像进行混合像元分解得到的积雪丰度图的精度进行评价。结果表明：利用基于SFIM变换和小波变换方法融合后影像提取的积雪分类图精度较高,数量精度为75%,比未融合影像积雪分类图的精度提高了6%,表明MODIS影像波段融合是一种提取高精度积雪信息的有效方法。     

Mapping regional land cover with MODIS data for biological conservation: Examples from the Greater Yellowstone Ecosystem, USA and Pará State, Brazil

The paper investigated the application of MODIS data for mapping regional land cover at moderate resolutions (250 and 500 m), for regional conservation purposes. Land cover maps were generated for two major conservation areas (Greater Yellowstone Ecosystem—GYE, USA and the Pará State, Brazil) using MODIS data and decision tree classifications. The MODIS land cover products were evaluated using existing Landsat TM land cover maps as reference data. The Landsat TM land cover maps were processed to their fractional composition at the MODIS resolution (250 and 500 m). In GYE, the MODIS land cover was very successful at mapping extensive cover types (e.g. coniferous forest and grasslands) and far less successful at mapping smaller habitats (e.g. wetlands, deciduous tree cover) that typically occur in patches that are smaller than the MODIS pixels, but are reported to be very important to biodiversity conservation. The MODIS classification for Pará State was successful at producing a regional forest/non-forest product which is useful for monitoring the extreme human impacts such as deforestation. The ability of MODIS data to map secondary forest remains to be tested, since regrowth typically harbors reduced levels of biodiversity. The two case studies showed the value of using multi-date 250 m data with only two spectral bands, as well as single day 500 m data with seven spectral bands, thus illustrating the versatile use of MODIS data in two contrasting environments. MODIS data provide new options for regional land cover mapping that are less labor-intensive than Landsat and have higher resolution than previous 1 km AVHRR or the current 1 km global land cover product. The usefulness of the MODIS data in addressing biodiversity conservation questions will ultimately depend upon the patch sizes of important habitats and the land cover transformations that threaten them.     

Generation of dense time series synthetic Landsat data through data blending with MODIS using a spatial and temporal adaptive reflectance fusion model

Landsat imagery with a 30 m spatial resolution is well suited for characterizing landscape-level forest structure and dynamics. While Landsat images have advantageous spatial and spectral characteristics for describing vegetation properties, the Landsat sensor's revisit rate, or the temporal resolution of the data, is 16 days. When considering that cloud cover may impact any given acquisition, this lengthy revisit rate often results in a dearth of imagery for a desired time interval (e.g., month, growing season, or year) especially for areas at higher latitudes with shorter growing seasons. In contrast, MODIS (MODerate-resolution Imaging Spectroradiometer) has a high temporal resolution, covering the Earth up to multiple times per day, and depending on the spectral characteristics of interest, MODIS data have spatial resolutions of 250 m, 500 m, and 1000 m. By combining Landsat and MODIS data, we are able to capitalize on the spatial detail of Landsat and the temporal regularity of MODIS acquisitions. In this research, we apply and demonstrate a data fusion approach (Spatial and Temporal Adaptive Reflectance Fusion Model, STARFM) at a mainly coniferous study area in central British Columbia, Canada. Reflectance data for selected MODIS channels, all of which were resampled to 500 m, and Landsat (at 30 m) were combined to produce 18 synthetic Landsat images encompassing the 2001 growing season (May to October). We compared, on a channel-by-channel basis, the surface reflectance values (stratified by broad land cover types) of four real Landsat images with the corresponding closest date of synthetic Landsat imagery, and found no significant difference between real (observed) and synthetic (predicted) reflectance values (mean difference in reflectance: mixed forest x? = 0.086, σ = 0.088, broadleaf x? = 0.019, σ = 0.079, coniferous x? = 0.039, σ = 0.093). Similarly, a pixel based analysis shows that predicted and observed reflectance values for the four Landsat dates were closely related (mean r² = 0.76 for the NIR band; r² = 0.54 for the red band; p < 0.01). Investigating the trend in NDVI values in synthetic Landsat values over a growing season revealed that phenological patterns were well captured; however, when seasonal differences lead to a change in land cover (i.e., disturbance, snow cover), the algorithm used to generate the synthetic Landsat images was, as expected, less effective at predicting reflectance.     

Enhanced algorithm for estimating snow depth from geostationary satellites

Observations in the visible and infrared spectral bands from the Imager instrument onboard Geostationary Operational Environmental Satellite (GOES) have been used to derive snow depth. The technique makes use of correlation between depth of the snow pack and satellite-derived subpixel fractional snow cover. Previous efforts to infer snow depth from satellite data with this technique were focused on grasslands and croplands, where the snow depth/snow fraction relationship is most pronounced. In this paper we improve the retrieval algorithm to extend snow depth estimates to forested areas. The enhanced algorithm accounts for the tree cover fraction and for the type of forest, deciduous or coniferous.The developed technique was used to derive maps of snow depth over mid-latitude areas of North America during winter seasons of 2003-2004 and 2004-2005. Satellite-based snow depth maps were produced daily at 4 km spatial resolution. To validate the retrievals we compared them with surface observations of snow depth and with the snow depth analysis prepared at the NOAA National Operational Hydrological Remote Sensing Center (NOHRSC). The estimated retrieval error was about 30% for snow depths below 30 cm and increased to 50% for snow depths ranging from 30 to 50 cm. Snow depth retrievals were limited to scenes with less than 80% deciduous forest cover fraction and less than 50% needle leaf forest cover.     

Integrated assessment on multi-temporal and multi-sensor combinations for reducing cloud obscuration of MODIS snow cover products of the Pacific Northwest USA

MODIS (Moderate Resolution Imaging Spectroradiometer) snow cover products, of daily, freely available, worldwide spatial extent at medium spatial resolution, have been widely applied in regional snow cover and modeling studies, although high cloud obscuration remains a concern in some applications. In this study, various approaches including daily combination, adjacent temporal deduction, fixed-day combination, flexible multi-day combination, and multi-sensor combination are assessed to remove cloud obscuration while still maintain the temporal and spatial resolutions. The performance of the resultant snow cover maps are quantitatively evaluated against in situ observations at 244 SNOTEL stations over the Pacific Northwest USA during the period of 2006-2008 hydrological years. Results indicate that daily Terra and Aqua MODIS combination and adjacent temporal deduction can reduce cloud obscuration and classification errors although an annual mean of 37% cloud coverage remains. Classification errors in snow-covered months are actually small and tend to underestimate the snow cover. Primary errors of MODIS daily, fixed and flexible multi-day combination products occur during transient months. Flexible multi-day combination is an efficient approach to maintain the balance between temporal resolution and realistic estimation of snow cover extent since it uses two thresholds to control the combination processes. Multi-sensor combinations (daily or multi-day), taking advantage of MODIS high spatial resolution and AMSR-E cloud penetration ability, provide cloud-free products but bring larger image underestimation errors as compared with their MODIS counterparts.     

Cropland distributions from temporal unmixing of MODIS data

Knowledge of the distribution of crop types is important for land management and trade decisions, and is needed to constrain remotely sensed estimates of variables, such as crop stress and productivity. The Moderate Resolution Imaging Spectroradiometer (MODIS) offers a unique combination of spectral, temporal, and spatial resolution compared to previous global sensors, making it a good candidate for large-scale crop type mapping. However, because of subpixel heterogeneity, the application of traditional hard classification approaches to MODIS data may result in significant errors in crop area estimation. We developed and tested a linear unmixing approach with MODIS that estimates subpixel fractions of crop area based on the temporal signature of reflectance throughout the growing season. In this method, termed probabilistic temporal unmixing (PTU), endmember sets were constructed using Landsat data to identify pure pixels, and uncertainty resulting from endmember variability was quantified using Monte Carlo simulation. This approach was evaluated using Landsat classification maps in two intensive agricultural regions, the Yaqui Valley (YV) of Mexico and the Southern Great Plains (SGP). Performance of the mixture model varied depending on the scale of comparison, with R² ranging from roughly 50% for estimating crop area within individual pixels to greater than 80% for crop cover within areas over 10 km². The results of this study demonstrate the importance of subpixel heterogeneity in cropland systems, and the potential of temporal unmixing to provide accurate and rapid assessments of land cover distributions using coarse resolution sensors, such as MODIS.     

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.     

An application of MODIS data to snow cover monitoring in a pastoral area: A case study in Northern Xinjiang, China

Snow is an important land cover on the earth's surface. It is characterized by its changing nature. Monitoring snow cover extent plays a significant role in dynamic studies and prevention of snow-caused disasters in pastoral areas. Using NASA EOS Terra/MODIS snow cover products and in situ observation data during the four snow seasons from November 1 to March 31 of year 2001 to 2005 in northern Xinjiang area, the accuracy of MODIS snow cover mapping algorithm under varied snow depth and land cover types was analyzed. The overall accuracy of MODIS daily snow cover mapping algorithm in clear sky condition is high at 98.5%; snow agreement reaches 98.2%, and ranges from 77.8% to 100% over the 4-year period for individual sites. Snow depth (SD) is one of the major factors affecting the accuracy of MODIS snow cover maps. MODIS does not identify any snow for SD less than 0.5 cm. The overall accuracy increases with snow depth if SD is equal to or greater than 3 cm, and decreases for SD below 3 cm. Land cover has an important influence in the accuracy of MODIS snow cover maps. The use of MOD10A1 snow cover products is severely affected by cloud cover. The 8-day composite products of MOD10A2 can effectively minimize the effect of cloud cover in most cases. Cloud cover in excess of 10% occurs on 99% of the MOD10A1 products and 14.7% of the MOD10A2 products analyzed during the four snow seasons. User-defined multiple day composite images based on MOD10A1, with flexibilities of selecting composite period, starting and ending date and composite sequence of MOD10A1 products, have an advantage in effectively monitoring snow cover extent for regional snow-caused disasters in pastoral areas.     

Using MODIS snow cover maps in modeling snowmelt runoff process in the eastern part of Turkey

Water perhaps is the most valuable natural asset in the Middle East as it was a historical key for settlement and survival in Mesopotamia, “the land between two rivers”. At present, the Euphrates and Tigris are the two largest trans-boundary rivers in Western Asia where Turkey, Syria, Iran, Iraq and Saudi Arabia are the riparian countries. The Euphrates and Tigris basins are largely fed from snow precipitation whereby nearly two-thirds occur in winter and may remain in the form of snow for half of the year. The concentration of discharge mainly from snowmelt during spring and early summer months causes not only extensive flooding, inundating large areas, but also the loss of much needed water required for irrigation and power generation purposes during the summer season. Accordingly, modeling of snow-covered area in the mountainous regions of Eastern Turkey, as being one of the major headwaters of Euphrates-Tigris basin, has significant importance in order to forecast snowmelt discharge especially for energy production, flood control, irrigation and reservoir operation optimization.A pilot basin, located on the upper Euphrates River, is selected where five automated meteorological and snow stations are installed for real time operations. The daily snow cover maps obtained from Moderate Resolution Imaging Spectroradiometer MODIS at 500 m resolution are compared with ground information for the winter of 2002-2003 both during accumulation and ablation and at accumulation stage for the winter of 2003-2004. The snow presence on the ground is determined from the snow courses performed. Such measurements were made at 19 points in and around the upper Euphrates River in Turkey and at 20 points in the upper portion of the pilot basin for the winters of 2002-2003 and 2003-2004, respectively. Comparison of MODIS snow maps with in situ measurements over the snow season show good agreement with overall accuracies ranging between 62% and 82% considering the shift in the days of comparison. The main reasons to have disagreement between MODIS and in situ data are the high cloud cover frequency in the area and the current version of the MODIS cloud-mask that appears to frequently map edges of snow-covered areas and land surfaces. The effect of elevation and land cover types on validation of MODIS snow cover maps is also analyzed. In order to minimize the cloud cover and maximize the snow cover, MODIS-8 daily snow cover products are used in deriving the snow depletion curve, which is one of the input parameters of the snowmelt runoff model (SRM). The initial results of modeling process show that MODIS snow-covered area product can be used for simulation and also for forecasting of snowmelt runoff in basins of Turkey.     

Developing clear-sky,cloud and cloud shadow mask for producing clear-sky composites at 250-meter spatial resolution for the seven MODIS land bands over Canada and North America

A new technology was developed at the Canada Centre for Remote Sensing (CCRS) for generating Canada-wide and North America continental scale clear-sky composites at 250 m spatial resolution for all seven MODIS land spectral bands (B1–B7). The MODIS Level 1B (MOD02) swath level data are used as input to circumvent the problems with image distortion in the mid latitude and polar regions inherent to the global sinusoidal (SIN) projection utilized for the standard MODIS data products. The MODIS 500 m land bands B3 to B7 are first downscaled to 250 m resolution using an adaptive regression and normalization scheme for compatibility with the 250 m bands B1 and B2. A new method has been developed to produce the mask of clear-sky, cloud and cloud shadow at 250 m resolution. It shows substantial advantages in comparison with the MODIS 250 m standard cloud masks. The testing of new cloud mask showed that it is in reasonable agreement with the MODIS 1-km standard product once it is aggregated to 1-km scale, while the cloud shadow detection looks more reliable with the new methodology. Nevertheless, more quantitative analyses of the presented scene identification technique are required to understand its performance over the range of input scenes in various seasons. The new clear-sky compositing scheme employs a scene-dependent decision matrix. It is demonstrated that this new scheme provides better results than any others based on a single compositing criterion, such as maximum NDVI or minimum visible reflectance. To account for surface bi-directional properties, two clear-sky composites for the same time period are produced by separating backward scattering and forward scattering geometries, which separate pixels with the sun-satellite relative azimuth angles within 90°–270° and outside of this range. Comparison with Landsat imagery and with MODIS standard composite products demonstrated the advantage of the new technique for screening cloud and cloud shadow, and generating high spatial resolution MODIS clear-sky composites. The new data products are mapped in the Lambert Conformal Conic (LCC) projection for Canada and the Lambert Azimuthal Equal-Area (LAEA) projection for North America. Presently this activity is limited to MODIS/TERRA due to known problems with band-to-band registration and noisy SWIR channels on MODIS/AQUA.     

A prototype MODIS-SSM/I snow-mapping algorithm

Data in the wavelength range 0.545-1.652 w m from the Moderate Resolution Imaging Spectroradiometer (MODIS), launched aboard the Earth Observing System (EOS) Terra in December 1999, are used to map daily global snow cover at 500 m resolution. However, during darkness, or when the satellite's view of the surface is obscured by cloud, snow cover cannot be mapped using MODIS data. We show that during these conditions, it is possible to supplement the MODIS product by mapping the snow cover using passive microwave data from the Special Sensor Microwave Imager (SSM/I), albeit with much poorer resolution. For a 7-day time period in March 1999, a prototype MODIS snow-cover product was compared with a prototype MODIS-SSM/I product for the same area in the mid-western USA. The combined MODIS-SSM/I product mapped 9% more snow cover than the MODIS-only product.     

Validation of operational AVHRR subpixel snow retrievals over the European Alps based on ASTER data

Snow is of great economic and social importance for the European Alps. Accurate monitoring of the alpine snow cover is a key component in studying regional climate change as well as in daily weather forecasting and snowmelt run‐off modelling. These applications require snow cover information on a high temporal resolution in near‐real time. For the European Alps, operational snow cover fraction maps are generated on a daily basis using data from the Advanced Very High Resolution Radiometer (AVHRR) on board the National Oceanic and Atmospheric Administration (NOAA) platforms. Snow cover distribution is inherently discontinuous and heterogeneous in this mountainous region. We have therefore implemented a straightforward multiple endmember unmixing approach to estimate fractional snow cover. Subpixel proportions are difficult to validate because similar products are not available and appropriate ground‐based observations do not exist. In this study, we validate AVHRR subpixel snow retrievals using binary classified data sets from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) to establish absolute errors of our operational approach at three test sites. Our analysis indicates that the AVHRR subpixel maps compare well with the aggregated ASTER data, showing an overall correlation of 0.78 and providing subpixel estimates with a mean absolute error of 10.4% fractional snow cover. Discrepancies between AVHRR and ASTER snow fraction maps can be attributed to varying snow conditions, terrain effects and density in forest cover.     

东北地区MODIS亚像元积雪覆盖率反演及验证

以中巴资源卫星数据作为地面“真值”影像,根据东北地区地理环境与气候特点对Salomoson亚像元积雪覆盖率模型参数进行修正,反演东北地区MODIS像元积雪覆盖率,并用不同方案对模型的稳定性和精度进行分析。研究结果表明,经修正后的Salomoson亚像元积雪覆盖率反演模型对不同地貌--景观单元具有稳定性,其中较小的波动源于积雪物理性质差异、大气效应、积雪影像分类误差及影像配准误差。在东北平原区,NDSI值在0.52~0.65时,模型反演精度高,但反演雪盖率总体偏低,主要是由NDSI基于对波段反射率的非线性转换引起的;雪盖率高估的像元主要分布在城区外围以及农村居民点,而覆盖城区、乡、镇以及居民点之间道路的像元雪盖率误差小,其原因是人类活动频率影响像元内积雪组分与非积雪组分的光谱特性的差异程度。与MODIS雪产品进行对比分析,积雪覆盖率提供较传统雪盖制图更加丰富的信息,然而对林区冠层下积雪覆盖二者均未给出准确估计。