Subpixel mapping of snow cover in forests by optical remote sensing

Forest represents a challenging problem for snow-cover mapping by optical satellite remote sensing. To investigate reflectance variability and to improve the mapping of snow in forested areas, a method for subpixel mapping of snow cover in forests (SnowFrac) has been developed. The SnowFrac method is based on linear spectral mixing modelling of snow, trees and snow-free ground. The focus has been on developing a physically based reflectance model which uses a forest-cover map as prior information. The method was tested in flat terrain covered by spruce, pine and birch forests, close to the Jotunheimen region of South Norway. Experiments were carried out using a completely snow-covered Landsat Thematic Mapper (TM) scene, aerial photos and in situ reflectance measurements. A detailed forest model was photogrammetrically derived from the aerial photos. Modelled and observed TM reflectances were compared. In the given situation, the results demonstrate that snow and individual tree species, in addition to cast shadows on the snow surface from single trees, are the most influencing factors on visible and near-infrared reflectance. Modelling of diffuse radiation reduced by surrounding trees slightly improve the results, indicating that this effect is less important. The best results are obtained for pine forest and mixed pine and birch forest. Future work will focus on deriving a simplified reflectance model suitable for operational snow-cover mapping in forests.     

Snow-cover mapping in forests by constrained linear spectral unmixing of MODIS data

A snow-cover mapping method accounting for forests (SnowFrac) is presented. SnowFrac uses spectral unmixing and endmember constraints to estimate the snow-cover fraction of a pixel. The unmixing is based on a linear spectral mixture model, which includes endmembers for snow, conifer, branches of leafless deciduous trees and snow-free ground. Model input consists of a land-cover fraction map and endmember spectra. The land-cover fraction map is applied in the unmixing procedure to identify the number and types of endmembers for every pixel, but also to set constraints on the area fractions of the forest endmembers. SnowFrac was applied on two Terra Moderate Resolution Imaging Spectroradiometer (MODIS) images with different snow conditions covering a forested area in southern Norway. Six experiments were carried out, each with different endmember constraints. Estimated snow-cover fractions were compared with snow-cover fraction reference maps derived from two Landsat Enhanced Thematic Mapper Plus (ETM+) images acquired the same days as the MODIS images. Results are presented for non-forested areas, deciduous forests, coniferous forests and mixed deciduous/coniferous forests. The snow-cover fraction estimates are enhanced by increasing constraints introduced to the unmixing procedure. The classification accuracy shows that 96% of the pixels are classified with less than 20% error (absolute units) on 7 May 2001 when all forested and non-forested areas are included. The corresponding figure for 4 May 2000 is 88%.     

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.     

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.     

Annual changes in MODIS vegetation indices of Swedish coniferous forests in relation to snow dynamics and tree phenology

Remote sensing provides spatially and temporally continuous measures of forest reflectance, and vegetation indices calculated from satellite data can be useful for monitoring climate change impacts on forest tree phenology. Monitoring of evergreen coniferous forest is more difficult than monitoring of deciduous forest, as the new buds only account for a small proportion of the green biomass, and the shoot elongation process is relatively slow. In this study, we have analyzed data from 186 coniferous monitoring sites in Sweden covering boreal, southern-boreal, and boreo-nemoral conditions. Our objective was to examine the possibility to track seasonal changes in coniferous forests by time-series of MODIS eight-day vegetation indices, testing the coherence between satellite monitored vegetation indices (VI) and temperature dependent phenology. The relationships between two vegetation indices (NDVI and WDRVI) and four phenological indicators (length of snow season, modeled onset of vegetation period, tree cold hardiness level and timing of budburst) were analyzed.The annual curves produced by two curve fitting methods for smoothening of seasonal changes in NDVI and WDRVI were to a large extent characterized by the occurrence of snow, producing stable seasonal oscillations in the northern part and irregular curves with less pronounced annual amplitude in the southern part of the country. Measures based on threshold values of the VI-curves, commonly used for determining the timing of different phenological phases, were not applicable for Swedish coniferous forests. Evergreen vegetation does not have a sharp increase in greenness during spring, and the melting of snow can influence the vegetation indices at the timing of budburst in boreal forests. However, the interannual variation in VI-values for specific eight-day periods was correlated with the phenological indicators. This relation can be used for satellite monitoring of potential climate change impacts on northern coniferous spring phenology.     

Accuracy assessment of MODIS,NOAA and IRS data in snow cover mapping under Himalayan conditions

Snow cover information is an essential parameter for a wide variety of scientific studies and management applications, especially in snowmelt runoff modelling. Until now NOAA and IRS data were widely and effectively used for snow‐covered area (SCA) estimation in several Himalayan basins. The suit of snow cover products produced from MODIS data had not previously been used in SCA estimation and snowmelt runoff modelling in any Himalayan basin. The present study was conducted with the aim of assessing the accuracy of MODIS, NOAA and IRS data in snow cover mapping under Himalayan conditions. The total SCA was estimated using these three datasets for 15 dates spread over 4 years. The results were compared with ground‐based estimation of snow cover. A good agreement was observed between satellite‐based estimation and ground‐based estimation. The influence of aspect in SCA estimation was analysed for the three satellite datasets and it was observed that MODIS produced better results. Snow mapping accuracy with respect to elevation was tested and it was observed that at higher elevation MODIS sensed more snow and proved better at mapping snow under mountain shadow conditions. At lower elevation, IRS proved better in mapping patchy snow cover due to higher spatial resolution. The temporal resolution of MODIS and NOAA data is better than IRS data, which means that the chances of getting cloud‐free scenes is higher. In addition, MODIS has an automated snow‐mapping algorithm, which reduces the time and errors incorporated during processing satellite data manually. Considering all these factors, it was concluded that MODIS data could be effectively used for SCA estimation under Himalayan conditions, which is a vital parameter for snowmelt runoff estimation.     

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.     

Seasonal and inter-annual variation in view angle effects on MODIS vegetation indices at three forest sites

The bi-directional reflectance distribution function (BRDF) has been widely studied across different vegetation types. However, these studies generally report values for only one point in time. We were interested in the potential for seasonal and inter-annual variation in BRDF parameters. NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board the EOS satellites has now been collecting data for 10 years. Since BRDF parameters are reported for the individual spectral bands, these data can be used to examine intra-annual variation. However, MODIS BRDF parameters are not calculated for the various vegetation indices which are derived from the spectral bands. Our objective in this study was to use the 10 years of MODIS data now available to examine seasonal and inter-annual variation in the view angle sensitivity of three vegetation indices; the normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), and the photochemical reflectance index (PRI) at 3 flux tower sites (Harvard Forest, Howland Forest and Morgan Monroe State Forest). For these 3 sites, only EVI was significantly affected by view angle. There was also a substantial variation in the view angle sensitivity of EVI across seasons and this variation was different for backscatter vs. forward scatter data. It is possible that differences in the scattering of radiation between the spring and the fall are responsible for the seasonal difference in view angle sensitivity. There were also complimentary variations in forward and backscatter view angle sensitivity of EVI across years. The greater view angle sensitivity of EVI, as opposed to NDVI, suggests that greater care must be taken to correct for BRDF effects when using this vegetation index.     

View angle effects in the radiometric measurement of plant canopy temperatures

The thermal infrared sensor response from a wheat canopy was extremely non-Lambertian because of spatial variations in energy flow processes; the effective radiant temperature of the sensor varied as much as 13°C with changing view angle. This variation of sensor response was accurately quantified (root-mean-square of deviations between theoretical and measured responses reduced to 1.1°C) as a function of vegetation canopy geometry, vertical temperature distribution of canopy components, and sensor view angle. The results have important implications for optimizing sensor view angles for remote sensing missions.     

Retrieval of snow grain size over Greenland from MODIS

This paper presents a new automatic algorithm to derive optical snow grain size at 1 km resolution using Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. The retrieval is conceptually based on an analytical asymptotic radiative transfer model which predicts spectral bidirectional snow reflectance as a function of the grain size and ice absorption. The snow grains are modeled as fractal rather than spherical particles in order to account for their irregular shape. The analytical form of solution leads to an explicit and fast retrieval algorithm. The time series analysis of derived grain size shows a good sensitivity to snow melting and snow precipitation events. Pre-processing is performed by a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which includes gridding MODIS data to 1 km resolution, water vapor retrieval, cloud masking and an atmospheric correction. MAIAC cloud mask is a new algorithm based on a time series of gridded MODIS measurements and an image-based rather than pixel-based processing. Extensive processing of MODIS TERRA data over Greenland shows a robust discrimination of clouds over bright snow and ice. Because in-situ grain size measurements over Greenland were not available at the time of this work, the validation was performed using data of Aoki et al. (Aoki, T., Hori, M., Motoyoshi, H., Tanikawa, T., Hachikubo, A., Sugiura, K., et al. (2007). ADEOS-II/GLI snow/ice products — Part II: Validation results using GLI and MODIS data. Remote Sensing of Environment, 111, 274-290) collected at Barrow (Alaska, USA), and Saroma, Abashiri and Nakashibetsu (Japan) in 2001-2005. The retrievals correlate well with measurements in the range of radii ~ 0.1-1 mm, although retrieved optical diameter may be about a factor of 1.5 lower than the physical measured diameter. As part of validation analysis for Greenland, the derived grain size from MODIS over selected sites in 2004 was compared to the microwave brightness temperature measurements of SSM/I radiometer which is sensitive to the amount of liquid water in the snowpack. The comparison showed a good qualitative agreement, with both datasets detecting two main periods of snowmelt. Additionally, MODIS grain size was compared with predictions of the snow model CROCUS driven by measurements of the automatic weather stations of the Greenland Climate Network. We found that the MODIS value is on average a factor of two smaller than CROCUS grain size. This result agrees with the direct validation analysis indicating that the snow reflectance model may need a “calibration” factor of ~ 1.5 for the retrieved grain size to match the physical snow grain size. Overall, the agreement between CROCUS and MODIS results was satisfactory, in particular before and during the first melting period in mid-June. Following detailed time series analysis of snow grain size for four permanent sites, the paper presents maps of this important parameter over the Greenland ice sheet for the March-September period of 2004.     

View azimuth and zenith,and solar angle effects on wheat canopy reflectance

Since crop canopies are not lambertian reflectors, their reflectance varies with sun and view positions. It is not always possible or convenient to make reflectance measurements from the nadir position nor at the same time of day. Therefore, ways of estimating nadir reflectance from off-nadir views and for various solar zenith angles are needed. In this study, spectral measurements were made with a Mark II radiometer five times during the day on each of four dates from 15° interval zenith and 45° azimuth positions for wheat canopies during the development interval stem extension to watery ripeness of the grain. The ratio of off-nadir [R(Z_v,A_v)] to nadir [R(0)] radiance in NIR band (0.76–0.90 μm) was described by the regression equation: $\text{R(Z}{}_{\mathrm{v}}\text{,A}{}_{\mathrm{v}}\text{)}\text{R(0)}\text{= 1.0 + [β}{}_0\text{+ β}{}_1\text{sin}\text{(}\text{A}{}_{\mathrm{v}}\text{2}\text{) + β}{}_2\text{(1/}\text{cos}\text{Z}{}_{\mathrm{s}}\text{)]}\text{sin}\text{Z}{}_{\mathrm{v}}$ where A_v is view azimuth angle relative to the sum position, Z_s is solar zenith angle, and Z_v is view zenith angle. The coefficient of determination was 0.70 or higher. The equation describes the observations that 1) the ratio of off-nadir to nadir radiance increases or decreases as view zenith angle increases depending on view azimuth angle; backscattering is stronger than forwardscattering and the pattern is azimuthally symmetric about the principal plane of the sun; and 2) the rate of change in the radiance ratio increases with increasing solar zenith angle. The coefficients, β₀, β₁ and β₂, changed as the canopies grew. Although the equation needs to be more fully tested, it should help summarize and compare various angular observation data taken in crop fields.     

Spatiotemporal variability of snow cover on Tibet,China using MODIS remote-sensing data

ABSTRACT

Snow cover is an important component of the cryosphere, and the study on spatial and temporal variations of snow cover is essential for understanding the consequences and impacts of climate change and water resources management. In this study, the temporal variation of snow-covered area (SCA) and spatial variability of snow-cover frequency (SCF) on Tibet is analysed based on the Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra snow cover product (MOD10A2) from 2000 to 2015, and relationships with main climate variables are investigated. Results are as follows: (1) there is a very weak decreasing trend in annual mean SCA, and a slight increasing trend in autumn and winter and a slight decreasing trend in spring and more robust decreasing trend in summer for SCA are found. (2) The temporal variation of SCA is negatively correlated with temperature, whereas it is little correlated with corresponding precipitation. (3) The general trend of spatial SCF variation on Tibet, predominated by snow-cover variations in spring and autumn, tends to decrease in spring while it tends to increase in autumn. (4) The spatial variability of SCF is attributed to snow-cover variations in autumn and spring, which is more obvious in higher latitudes in autumn while it is more noticeable in lower-latitude southeastern plateau in spring. (5) The regions with higher variability of snow cover are main pastoral land and more prone to snow-related disaster in Tibet, becoming key zone of snow-cover monitoring and disaster prevention and mitigation.     

The effects of variation in snow properties on passive microwave snow mass estimation

Estimating snow mass at continental scales is difficult, but important for understanding land-atmosphere interactions, biogeochemical cycles and the hydrology of the Northern latitudes. Remote sensing provides the only consistent global observations, but with unknown errors. We test the theoretical performance of the Chang algorithm for estimating snow mass from passive microwave measurements using the Helsinki University of Technology (HUT) snow microwave emission model. The algorithm's dependence upon assumptions of fixed and uniform snow density and grainsize is determined, and measurements of these properties made at the Cold Land Processes Experiment (CLPX) Colorado field site in 2002–2003 used to quantify the retrieval errors caused by differences between the algorithm assumptions and measurements. Deviation from the Chang algorithm snow density and grainsize assumptions gives rise to an error of a factor of between two and three in calculating snow mass. The possibility that the algorithm performs more accurately over large areas than at points is tested by simulating emission from a 25 km diameter area of snow with a distribution of properties derived from the snow pit measurements, using the Chang algorithm to calculate mean snow-mass from the simulated emission. The snow mass estimation from a site exhibiting the heterogeneity of the CLPX Colorado site proves only marginally different than that from a similarly-simulated homogeneous site. The estimation accuracy predictions are tested using the CLPX field measurements of snow mass, and simultaneous SSM/I and AMSR-E measurements.     

Operational snow mapping using multitemporal Meteosat SEVIRI imagery

The Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation is the first geostationary satellite instrument with all visual and infrared channels that are important for snow mapping. In this paper, we present an algorithm for deriving snow cover maps from SEVIRI data that makes use of the unique combination of adequate spectral resolution and very high frequency. The short interval of 15 min between images makes it possible to extend traditional spectral classification with a detection of changes between images. This improves the detection of clouds and cloud shadows in instantaneous images, because these often display more variation in time than the surface. It therefore allows a more accurate mapping of surface snow cover, as is shown by a validation of the results with ground observations and other satellite data. The accurate classification of each single image allows the generation of temporal composite snow maps in near real-time, which is for example of interest for numerical weather prediction models. When compared to many in situ measurements from the winter of 2005/2006, the accuracy of the algorithm is 95%.     

A new water-resistant snow index for the detection and mapping of snow cover on a global scale

An up-to-date spatio-temporal change analysis of global snow cover is essential for better understanding of climate–hydrological interactions. The normalized difference snow index (NDSI) is a widely used algorithm for the detection and estimation of snow cover. However, NDSI cannot discriminate between snow cover and water bodies without use of an external water mask. A stand-alone methodology for robust detection and mapping of global snow cover is presented by avoiding external dependency on the water mask. A new spectral index called water-resistant snow index (WSI) with the capability of exhibiting significant contrast between snow cover and other cover types, including water bodies, was developed. WSI uses the normalized difference between the value and hue obtained by transforming red, green, and blue, (RGB) colour composite images comprising red, green, and near-infrared bands into a hue, saturation, and value (HSV) colour model. The superiority of WSI over NDSI is confirmed by case studies conducted in major snow regions globally. Snow cover was mapped by considering monthly variation in snow cover and availability of satellite data at the global scale. A snow cover map for the year 2013 was produced at the global scale by applying the random walker algorithm in the WSI image supported by the reference data collected from permanent snow-covered and non-snow-covered areas. The resultant snow-cover map was compared to snow cover estimated by existing maps: MODIS Land Cover Type Product (MCD12Q1 v5.1, 2012), Global Land Cover by National Mapping Organizations (GLCNMO v2.0, 2008), and European Space Agency’s GlobCover 2009. A significant variation in snow cover as estimated by different maps was noted, and was was attributed to methodological differences rather than annual variation in snow cover. The resultant map was also validated with reference data, with 89.46% overall accuracy obtained. The WSI proposed in the research is expected to be suitable for seasonal and annual change analysis of global snow cover.     

Evaluating MODIS data for mapping wildlife habitat distribution

Habitat distribution models have a long history in ecological research. With the development of geospatial information technology, including remote sensing, these models are now applied to an ever-increasing number of species, particularly those located in areas in which it is logistically difficult to collect habitat data in the field. Many habitat studies have used data acquired by multi-spectral sensor systems such as the Landsat Thematic Mapper (TM), due mostly to their availability and relatively high spatial resolution (30 m/pixel). The use of data collected by other sensor systems with lower spatial resolutions but high frequency of acquisitions has largely been neglected, due to the perception that such low spatial resolution data are too coarse for habitat mapping. In this study we compare two models using data from different satellite sensor systems for mapping the spatial distribution of giant panda habitat in Wolong Nature Reserve, China. The first one is a four-category scheme model based on combining forest cover (derived from a digital land cover classification of Landsat TM imagery acquired in June, 2001) with information on elevation and slope (derived from a digital elevation model obtained from topographic maps of the study area). The second model is based on the Ecological Niche Factor Analysis (ENFA) of a time series of weekly composites of WDRVI (Wide Dynamic Range Vegetation Index) images derived from MODIS (Moderate Resolution Imaging Spectroradiometer – 250 m/pixel) for 2001. A series of field plots was established in the reserve during the summer–autumn months of 2001–2003. The locations of the plots with panda feces were used to calibrate the ENFA model and to validate the results of both models. Results showed that the model using the seasonal variability of MODIS-WDRVI had a similar prediction success to that using Landsat TM and digital elevation model data, albeit having a coarser spatial resolution. This suggests that the phenological characterization of the land surface provides an appropriate environmental predictor for giant panda habitat mapping. Therefore, the information contained in remotely sensed data acquired with low spatial resolution but high frequency of acquisitions has considerable potential for mapping the habitat distribution of wildlife species.     

Accuracy of wet snow mapping using simulated Radarsat backscattering coefficients from observed snow cover characteristics

Wet snow cover mapping by means of airborne and spaceborne SAR is operational today and successfully applied in rugged high mountain terrain and in agricultural area. This paper proposes a numerical study to estimate the accuracy of wet snow mapping by using a radar backscattering model that simulates backscattering from a multi-layer snowpack for various snow cover conditions and for SAR parameters specific to Radarsat (C-HH). Field measurements carried out in numerous sites during the winters of 1994 to 1996 in several areas of Quebec (Canada) have allowed to choose some typical snow profiles and the corresponding snow/soil parameters. Results indicate that under the assumptions used in the model and the simulations, for the standard mode S1 of Radarsat (20 to 27.4) and in the case of wet snow cover with liquid water content of 1%, the optimum relative under-and over-estimation of wet snow pixels are of the order of 23.9% and 13.4%, respectively. For wet snow cover at 2%, the algorithm operates with a relative under-estimation of wet snow pixels around 8.5% and a relative over-estimation of the order of 1.7%. For wet snow with liquid water content of 4%, the relative under-and over-estimation of wet snow pixels is around 0.8% and 0.3%, respectively. They are negligible for wet snow with liquid water content higher than 4%. With the standard mode S7 of Radarsat (44.9 to 49.4), the wet snow mapping algorithm leads to a slightly lower performance than with the standard mode S1. The accuracy of the method for wet snow mapping demonstrates the high potential of SAR for snow monitoring. It is considered sufficient when the liquid water content of the snowpack is higher than 1% for actual snow conditions similar to those eight observed conditions used in this study.     

Blending MODIS and Landsat images for urban flood mapping

Satellite images provide important data sources for monitoring flood disasters. However, the trade-off between spatial and temporal resolutions of current satellite sensors limits their uses in urban flooding studies. This study applied and compared two data fusion models, the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM), in generating synthetic flooding images with improved temporal and spatial resolution for flood mapping. The synthetic images are produced in two scenarios: (1) for real-time prediction based on Landsat and MODIS images acquired before the investigated flooding; and (2) for post-disaster prediction based on images acquired after the flooding. The 2005 Hurricane Katrina in New Orleans was selected as a case study. The result shows that the Landsat-like images generated can be successfully applied in flood mapping. Particularly, ESTARFM surpasses STARFM in predicting surface reflectance in both real-time and post-flooding predictions. However, the flood mapping results from the Landsat-like images produced by both STARFM and ESTARFM are similar with overall accuracy around 0.9. Only for the flooding maps of real-time predictions does ESTARFM get a slightly higher overall accuracy than STARFM, indicating that the lower quality of the Landsat-like image generated by STARFM may not affect flood mapping accuracy, due to the marked contrast between land and water. This study suggests great potential of both STARFM and ESTARFM in urban flooding research. Blending multi-sources images could also support other disaster studies that require remotely sensed data with both high spatial and temporal resolution.     

Evaluation of MODIS land surface temperature with in-situ snow surface temperature from CREST-SAFE

This article presents the procedure and results of a temperature-based validation approach for the Moderate Resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature (LST) product provided by the National Aeronautics and Space Administration Terra and Aqua Earth Observing System satellites using in-situ LST observations recorded at the Cooperative Remote Sensing Science and Technology Center – Snow Analysis and Field Experiment (CREST-SAFE) during the years of 2013 (January–April) and 2014 (February–April). A total of 314 day-and-night clear-sky thermal images, acquired by the Terra and Aqua satellites, were processed and compared to ground-truth data from CREST-SAFE with a frequency of one measurement every 3 min. CREST-SAFE is a synoptic ground station, located in the cold county of Caribou in Maine, USA, with a distinct advantage over most meteorological stations because it provides automated and continuous LST observations via an Apogee Model SI-111 Infrared Radiometer. This article also attempts to answer the question of whether a single pixel (1 km²) or several spatially averaged pixels should be used for satellite LST validation by increasing the MODIS window size to 5 × 5, 9 × 9, and 25 × 25 windows.

Several trends in the MODIS LST data were observed, including the underestimation of daytime values and night-time values. Results indicate that although all the data sets (Terra and Aqua, diurnal and nocturnal) showed high correlation with ground measurements, day values yielded slightly higher accuracy (about 1°C), both suggesting that MODIS LST retrievals are reliable for similar land-cover classes and atmospheric conditions. Increasing the MODIS window size showed an overestimation of in-situ LST and some improvement in the daytime Terra and night-time Aqua biases, with the highest accuracy achieved with the 5 × 5 window. A comparison between MODIS emissivity from bands 31, 32, and in-situ emissivity showed that emissivity errors (relative error = ?0.30%) were insignificant.     

Using crowdsourced and weather station data to fill cloud gaps in MODIS snow cover datasets

We present the design, development, and testing of a new software package for generating snow cover maps. Using a custom inverse distance weighting method, we combine volunteer snow reports, cross-country ski track reports and station measurements to fill cloud gaps in the Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover product. The method is demonstrated by producing a continuous daily time step snow probability map dataset for the Czech Republic region. For validation, we checked the ability of our method to reconstruct MODIS snow cover under cloud by simulating cloud cover datasets and comparing estimated snow cover to actual MODIS snow cover. The percent correctly classified indicator showed accuracy between 80 and 90% using this method. The software is available as an R package. The output data sets are published on the HydroShare website for download and through a web map service for re-use in third-party applications.