Monitoring vegetation phenology using MODIS

Accurate measurements of regional to global scale vegetation dynamics (phenology) are required to improve models and understanding of inter-annual variability in terrestrial ecosystem carbon exchange and climate-biosphere interactions. Since the mid-1980s, satellite data have been used to study these processes. In this paper, a new methodology to monitor global vegetation phenology from time series of satellite data is presented. The method uses series of piecewise logistic functions, which are fit to remotely sensed vegetation index (VI) data, to represent intra-annual vegetation dynamics. Using this approach, transition dates for vegetation activity within annual time series of VI data can be determined from satellite data. The method allows vegetation dynamics to be monitored at large scales in a fashion that it is ecologically meaningful and does not require pre-smoothing of data or the use of user-defined thresholds. Preliminary results based on an annual time series of Moderate Resolution Imaging Spectroradiometer (MODIS) data for the northeastern United States demonstrate that the method is able to monitor vegetation phenology with good success.     

Satellite-based modeling of gross primary production in a seasonally moist tropical evergreen forest

A CO₂ eddy flux tower study has recently reported that an old-growth stand of seasonally moist tropical evergreen forest in Santarém, Brazil, maintained high gross primary production (GPP) during the dry seasons [Saleska, S. R., Miller, S. D., Matross, D. M., Goulden, M. L., Wofsy, S. C., da Rocha, H. R., de Camargo, P. B., Crill, P., Daube, B. C., de Freitas, H. C., Hutyra, L., Keller, M., Kirchhoff, V., Menton, M., Munger, J. W., Pyle, E. H., Rice, A. H., & Silva, H. (2003). Carbon in amazon forests: Unexpected seasonal fluxes and disturbance-induced losses. Science, 302, 1554-1557]. It was proposed that seasonally moist tropical evergreen forests have evolved two adaptive mechanisms in an environment with strong seasonal variations of light and water: deep roots system for access to water in deep soils and leaf phenology for access to light. Identifying tropical forests with these adaptive mechanisms could substantially improve our capacity of modeling the seasonal dynamics of carbon and water fluxes in the tropical zone. In this paper, we have analyzed multi-year satellite images from the VEGETATION (VGT) sensor onboard the SPOT-4 satellite (4/1998-12/2002) and the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra satellite (2000-2003). We reported temporal analyses of vegetation indices and simulations of the satellite-based vegetation photosynthesis model (VPM). The Enhanced Vegetation Index (EVI) identified subtle changes in the seasonal dynamics of leaf phenology (leaf emergence, leaf aging and leaf fall) in the forest, as suggested by the leaf litterfall data. The land surface water index (LSWI) indicated that the forest experienced no water stress in the dry seasons of 1998-2002. The VPM model, which uses EVI, LSWI and site-specific climate data (air temperature and photosynthetically active radiation, PAR) for 2001-2002, predicted high GPP in the late dry seasons, consistent with observed high evapotranspiration and estimated GPP from the CO₂ eddy flux tower.     

Land surface phenology from MODIS: Characterization of the Collection 5 global land cover dynamics product

Information related to land surface phenology is important for a variety of applications. For example, phenology is widely used as a diagnostic of ecosystem response to global change. In addition, phenology influences seasonal scale fluxes of water, energy, and carbon between the land surface and atmosphere. Increasingly, the importance of phenology for studies of habitat and biodiversity is also being recognized. While many data sets related to plant phenology have been collected at specific sites or in networks focused on individual plants or plant species, remote sensing provides the only way to observe and monitor phenology over large scales and at regular intervals. The MODIS Global Land Cover Dynamics Product was developed to support investigations that require regional to global scale information related to spatio-temporal dynamics in land surface phenology. Here we describe the Collection 5 version of this product, which represents a substantial refinement relative to the Collection 4 product. This new version provides information related to land surface phenology at higher spatial resolution than Collection 4 (500-m vs. 1-km), and is based on 8-day instead of 16-day input data. The paper presents a brief overview of the algorithm, followed by an assessment of the product. To this end, we present (1) a comparison of results from Collection 5 versus Collection 4 for selected MODIS tiles that span a range of climate and ecological conditions, (2) a characterization of interannual variation in Collections 4 and 5 data for North America from 2001 to 2006, and (3) a comparison of Collection 5 results against ground observations for two forest sites in the northeastern United States. Results show that the Collection 5 product is qualitatively similar to Collection 4. However, Collection 5 has fewer missing values outside of regions with persistent cloud cover and atmospheric aerosols. Interannual variability in Collection 5 is consistent with expected ranges of variance suggesting that the algorithm is reliable and robust, except in the tropics where some systematic differences are observed. Finally, comparisons with ground data suggest that the algorithm is performing well, but that end of season metrics associated with vegetation senescence and dormancy have higher uncertainties than start of season metrics.     

A Two-Step Filtering approach for detecting maize and soybean phenology with time-series MODIS data

The crop developmental stage represents essential information for irrigation scheduling/fertilizer management, understanding seasonal ecosystem carbon dioxide (CO₂) exchange, and evaluating crop productivity. In this study, we devised an approach called the Two-Step Filtering (TSF) for detecting the phenological stages of maize and soybean from time-series Wide Dynamic Range Vegetation Index (WDRVI) data derived from Moderate Resolution Imaging Spectroradiometer (MODIS) 250-m observations. The TSF method consists of a Two-Step Filtering scheme that includes: (i) smoothing the temporal WDRVI data with a wavelet-based filter and (ii) deriving the optimum scaling parameters from shape-model fitting procedure. The date of key crop development stages are then estimated by using the optimum scaling parameters and an initial value of the specific phenological date on the shape model, which are preliminary defined in reference to ground-based crop growth stage observations. The shape model is a crop-specific WDRVI curve with typical seasonal features, which were defined by averaging smoothed, multi-year WDRVI profiles from MODIS 250-m data collected over irrigated maize and soybean study sites.In this study, the TSF method was applied to MODIS-derived WDRVI data over a 6-year period (2003 to 2008) for two irrigated sites and one rainfed site planted to either maize or soybean as part of the Carbon Sequestration Program (CSP) at the University of Nebraska-Lincoln. A comparison of satellite-based retrievals with ground-based crop growth stage observations collected by the CSP over the six growing seasons for these three sites showed that the TSF method can accurately estimate the date of four key phenological stages of maize (V2.5: early vegetative stage, R1: silking stage, R5: dent stage and R6: maturity) and soybean (V1: early vegetative stage, R5: beginning seed, R6: full seed and R7: beginning maturity). The root mean square error (RMSE) of phenological-stage estimation for maize ranged from 2.9 [R1] to 7.0 [R5] days and from 3.2 [R6] to 6.9 [R7] days for soybean, respectively. In addition, the TSF method was also applied for two years (2001 and 2002) over eastern Nebraska to test its ability to characterize the spatio-temporal patterns of these key phenological stages over a larger geographic area. The MODIS-derived crop phenological stage dates agreed well with the statistical crop progress data reported by the United State Department of Agriculture (USDA) National Agricultural Statistics Service (NASS) for eastern Nebraska's three crop agricultural statistic districts (ASDs). At the ASD-level, the RMSE of phenological-stage estimation ranged from 1.6 [R1] to 5.6 [R5] days for maize and from 2.5 [R7] to 5.3 [R5] days for soybean.     

Mapping paddy rice agriculture in southern China using multi-temporal MODIS images

Information on the area and spatial distribution of paddy rice fields is needed for trace gas emission estimates, management of water resources, and food security. Paddy rice fields are characterized by an initial period of flooding and transplanting, during which period open canopy (a mixture of surface water and rice crops) exists. The Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the NASA EOS Terra satellite has visible, near infrared and shortwave infrared bands; and therefore, a number of vegetation indices can be calculated, including Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) that is sensitive to leaf water and soil moisture. In this study, we developed a paddy rice mapping algorithm that uses time series of three vegetation indices (LSWI, EVI, and NDVI) derived from MODIS images to identify that initial period of flooding and transplanting in paddy rice fields, based on the sensitivity of LSWI to the increased surface moisture during the period of flooding and rice transplanting. We ran the algorithm to map paddy rice fields in 13 provinces of southern China, using the 8-day composite MODIS Surface Reflectance products (500-m spatial resolution) in 2002. The resultant MODIS-derived paddy rice map was evaluated, using the National Land Cover Dataset (1:100,000 scale) derived from analysis of Landsat ETM+ images in 1999/2000. There were reasonable agreements in area estimates of paddy rice fields between the MODIS-derived map and the Landsat-based dataset at the provincial and county levels. The results of this study indicated that the MODIS-based paddy rice mapping algorithm could potentially be applied at large spatial scales to monitor paddy rice agriculture on a timely and frequent basis.     

A scalable approach to mapping annual land cover at 250 m using MODIS time series data: A case study in the Dry Chaco ecoregion of South America

Land use and land cover (LULC) maps from remote sensing are vital for monitoring, understanding and predicting the effects of complex human-nature interactions that span local, regional and global scales. We present a method to map annual LULC at a regional spatial scale with source data and processing techniques that permit scaling to broader spatial and temporal scales, while maintaining a consistent classification scheme and accuracy. Using the Dry Chaco ecoregion in Argentina, Bolivia and Paraguay as a test site, we derived a suite of predictor variables from 2001 to 2007 from the MODIS 250 m vegetation index product (MOD13Q1). These variables included: annual statistics of red, near infrared, and enhanced vegetation index (EVI), phenological metrics derived from EVI time series data, and slope and elevation. For reference data, we visually interpreted percent cover of eight classes at locations with high-resolution QuickBird imagery in Google Earth. An adjustable majority cover threshold was used to assign samples to a dominant class. When compared to field data, we found this imagery to have georeferencing error < 5% the length of a MODIS pixel, while most class interpretation error was related to confusion between agriculture and herbaceous vegetation. We used the Random Forests classifier to identify the best sets of predictor variables and percent cover thresholds for discriminating our LULC classes. The best variable set included all predictor variables and a cover threshold of 80%. This optimal Random Forests was used to map LULC for each year between 2001 and 2007, followed by a per-pixel, 3-year temporal filter to remove disallowed LULC transitions. Our sequence of maps had an overall accuracy of 79.3%, producer accuracy from 51.4% (plantation) to 95.8% (woody vegetation), and user accuracy from 58.9% (herbaceous vegetation) to 100.0% (water). We attributed map class confusion to limited spectral information, sub-pixel spectral mixing, georeferencing error and human error in interpreting reference samples. We used our maps to assess woody vegetation change in the Dry Chaco from 2002 to 2006, which was characterized by rapid deforestation related to soybean and planted pasture expansion. This method can be easily applied to other regions or continents to produce spatially and temporally consistent information on annual LULC.     

Systematic mapping of Leaf Area Index across Canada using 250-meter MODIS data

An approach to generate a 250-meter Canada wide Leaf Area Index (LAI) map using 250-meter MODIS data is described. The full processing chain is introduced. The approach is based on intercalibration of Landsat and MODIS vegetation indices (VI's) combined with LAI-VI's empirical relationships. Preliminary validation over two sites where field work was carried out shows promising results. Intercalibration of MODIS VI's before deriving LAI maps provides up to 65% improvement of the LAI overall accuracy.     

Mapping paddy rice agriculture in South and Southeast Asia using multi-temporal MODIS images

In this paper, we developed a new geospatial database of paddy rice agriculture for 13 countries in South and Southeast Asia. These countries have ∼ 30% of the world population and ∼ 2/3 of the total rice land area in the world. We used 8-day composite images (500-m spatial resolution) in 2002 from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the NASA EOS Terra satellite. Paddy rice fields are characterized by an initial period of flooding and transplanting, during which period a mixture of surface water and rice seedlings exists. We applied a paddy rice mapping algorithm that uses a time series of MODIS-derived vegetation indices to identify the initial period of flooding and transplanting in paddy rice fields, based on the increased surface moisture. The resultant MODIS-derived paddy rice map was compared to national agricultural statistical data at national and subnational levels. Area estimates of paddy rice were highly correlated at the national level and positively correlated at the subnational levels, although the agreement at the national level was much stronger. Discrepancies in rice area between the MODIS-derived and statistical datasets in some countries can be largely attributed to: (1) the statistical dataset is a sown area estimate (includes multiple cropping practices); (2) failure of the 500-m resolution MODIS-based algorithm in identifying small patches of paddy rice fields, primarily in areas where topography restricts field sizes; and (3) contamination by cloud. While further testing is needed, these results demonstrate the potential of the MODIS-based algorithm to generate updated datasets of paddy rice agriculture on a timely basis. The resultant geospatial database on the area and spatial distribution of paddy rice is useful for irrigation, food security, and trace gas emission estimates in those countries.     

Water content estimation from hyperspectral images and MODIS indexes in Southeastern Arizona

Hyperspectral water retrievals from AVIRIS data, equivalent water thickness (EWT), were compared to in situ leaf water content and LAI measurements at a semiarid site in southeastern Arizona. Retrievals of EWT showed good correlation with field canopy water content measurements. Statistical analysis suggested that EWT was significantly different among seven community types, from savanna to agriculture. Four band-ratio indexes (NDVI, EVI, NDWI, and NDII) were derived from MODIS showing strong spatial agreement between maps of AVIRIS EWT and MODIS indexes, and good statistical agreement for the range of habitats at the site. Temporal patterns of these four indexes in all vegetation communities except creosote bush and agriculture showed distinct seasonal patterns that responded to the timing and amount of precipitation. Moreover, these time series captured different ecological responses among the different vegetation communities.     

Effect of foliage spatial heterogeneity in the MODIS LAI and FPAR algorithm over broadleaf forests

This paper presents the analysis of radiative transfer assumptions underlying moderate resolution imaging spectroradiometer (MODIS) leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) algorithm for the case of spatially heterogeneous broadleaf forests. Data collected by a Boston University research group during the July 2000 field campaign at the Earth Observing System (EOS) core validation site, Harvard Forest, MA, were used for this purpose. The analysis covers three themes. First, the assumption of wavelength independence of spectral invariants of transport equation, central to the parameterization of the MODIS LAI and FPAR algorithm, is evaluated. The physical interpretation of those parameters is given and an approach to minimize the uncertainties in its retrievals is proposed. Second, the theoretical basis of the algorithm was refined by introducing stochastic concepts which account for the effect of foliage clumping and discontinuities on LAI retrievals. Third, the effect of spatial heterogeneity in FPAR was analyzed and compared to FPAR variation due to diurnal changes in solar zenith angle (SZA) to asses the validity of its static approximation.     

Optimal interpolation analysis of leaf area index using MODIS data

A simple data analysis technique for vegetation leaf area index (LAI) using Moderate Resolution Imaging Spectroradiometer (MODIS) data is presented. The objective is to generate LAI data that is appropriate for numerical weather prediction. A series of techniques and procedures which includes data quality control, time-series data smoothing, and simple data analysis is applied. The LAI analysis is an optimal combination of the MODIS observations and derived climatology, depending on their associated errors σ_o and σ_c. The “best estimate” LAI is derived from a simple three-point smoothing technique combined with a selection of maximum LAI (after data quality control) values to ensure a higher quality. The LAI climatology is a time smoothed mean value of the “best estimate” LAI during the years of 2002-2004. The observation error is obtained by comparing the MODIS observed LAI with the “best estimate” of the LAI, and the climatological error is obtained by comparing the “best estimate” of LAI with the climatological LAI value. The LAI analysis is the result of a weighting between these two errors. Demonstration of the method described in this paper is presented for the 15-km grid of Meteorological Service of Canada (MSC)'s regional version of the numerical weather prediction model. The final LAI analyses have a relatively smooth temporal evolution, which makes them more appropriate for environmental prediction than the original MODIS LAI observation data. They are also more realistic than the LAI data currently used operationally at the MSC which is based on land-cover databases.     

Evaluation of the onset of green-up in temperate deciduous broadleaf forests derived from Moderate Resolution Imaging Spectroradiometer (MODIS) data

Vegetation phenology is the chronology of periodic phases of development. It constitutes an efficient bio-indicator of impacts of climate changes and a key parameter for understanding and modelling vegetation-climate interactions and their implications on carbon cycling. Numerous studies were devoted to the remote sensing of vegetation phenology. Most of these were carried out using data acquired by AVHRR instrument onboard NOAA meteorological satellites. Since 1999, multispectral images were acquired over the whole earth surface every one to two days by MODIS instrument onboard Terra and Aqua platforms. In comparison with AVHRR, MODIS constitutes a significant technical improvement in terms of spatial resolution, spectral resolution, geolocation accuracy, atmospheric corrections scheme and cloud screening and sensor calibration. In this study, 250 m daily MODIS data were used to derive precise vegetation phenological dates over deciduous forest stands. Phenological markers derived from MODIS time-series and provided by MODIS Global Land Cover Dynamics product (MOD12Q2) were compared to field measurements carried out over the main deciduous forest stands across France and over five years. We show that the inflexion point of the asymmetric double-sigmoid function fitted to NDVI temporal profile is a good marker of the onset of green-up in deciduous stands. At plot level, the prediction uncertainty is 8.5 days and the bias is 3.5 days. MODIS Global Land Cover Dynamics MOD12Q2 provides estimates of onset of green-up dates which deviate substantially from in situ observations and do not perform better than the null model. RMSE values are 20.5 days (bias -17 days) using the onset of greenness increase and 36.5 days (bias 34.5 days) using the onset of greenness maximum. An improvement of prediction quality is obtained if we consider the average of MOD12Q2 onset of greenness increase and maximum as marker of onset of green-up date. RMSE decreases to 16.5 days and bias to 7.5 days.     

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.     

京津唐地区城市化对植被物候的影响研究

城市化过程使得土地利用、地表植被覆盖发生着显著的变化,如何定量化描述城市化对植被物候的影响越来越受到各方的关注。基于京津唐地区2001~2006年NDVI时间序列影像,得出了京津唐地区植被物候空间分布格局,计算出北京、天津、唐山3个核心城市城区的距离变量与平均植被物候,并分析了城市化对植被物候指标的影响趋势。结果表明:①2001~2006年,北京城市化使得城区及离城区较近的地方植被生长开始时间提前、结束时间推后、生长季周期变长、NDVI振幅减小;②天津和唐山的城市化使得城区及离城区较近的地方植被生长开始时间延后、结束时间提前、生长季周期变短\,NDVI振幅减小;③城市化对植被物候的影响与该地区城市扩张类型存在相关性关系。     

Disaggregation of MODIS surface temperature over an agricultural area using a time series of Formosat-2 images

The temporal frequency of the thermal data provided by current spaceborne high-resolution imagery systems is inadequate for agricultural applications. As an alternative to the lack of high-resolution observations, kilometric thermal data can be disaggregated using a green (photosynthetically active) vegetation index e.g. NDVI (Normalized Difference Vegetation Index) collected at high resolution. Nevertheless, this approach is only valid in the conditions where vegetation temperature is approximately uniform. To extend the validity domain of the classical approach, a new methodology is developed by representing the temperature difference between photosynthetically and non-photosynthetically active vegetation. In practice, both photosynthetically and non-photosynthetically active vegetation fractions are derived from a time series of Formosat-2 shortwave data, and then included in the disaggregation procedure. The approach is tested over a 16 km by 10 km irrigated cropping area in Mexico during a whole agricultural season. Kilometric MODIS (MODerate resolution Imaging Spectroradiometer) surface temperature is disaggregated at 100 m resolution, and disaggregated temperature is subsequently compared against concurrent ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data. Statistical results indicate that the new methodology is more robust than the classical one, and is always more accurate when fractional non-photosynthetically active vegetation cover is larger than 0.10. The mean correlation coefficient and slope between disaggregated and ASTER temperature is increased from 0.75 to 0.81 and from 0.60 to 0.77, respectively. The approach is also tested using the MODIS data re-sampled at 2 km resolution. Aggregation reduces errors in MODIS data and consequently increases the disaggregation accuracy.     

On intra-annual EVI variability in the dry season of tropical forest: A case study with MODIS and hyperspectral data

Considerable controversy is associated with dry season increases in the Enhanced Vegetation Index (EVI), observed using the Moderate Resolution Imaging Spectroradiometer (MODIS), compared with field-based estimates of decreasing plant productivity. Here, we investigate potential causes of intra-annual variability by comparing EVI from mature forest with field-measured Leaf Area Index (LAI) to validate space-based observations. EVI was calculated from 19 nadir and off-nadir Hyperion images in the 2005 dry season, and inspected for consistency with MODIS observations from 2004 to 2009. The objective was to evaluate the possible influence of the view-illumination geometry and of canopy foliage and leaf flush on the EVI. Spectral mixture models were used to evaluate the relationship between EVI and the shade fraction, a measure that varies with pixel brightness. MODIS LAI values were compared with LAI estimated using hemispherical photographs taken in two field campaigns in the dry season. To keep LAI and leaf flush conditions as constant variables and vary solar illumination, we used airborne Hyperspectral Mapper (Hymap) data acquired over mature forest from another region on the same day but with two distinct solar zenith angles (SZA) (29° and 53°). Results showed that intra-annual variability in MODIS and nadir Hyperion EVI in the dry season of tropical forest were driven by solar illumination effects rather than changes in LAI. The reflectance of the MODIS and Hyperion blue, red and near infrared (NIR) bands was higher at the end of the dry season because of the predominance of sunlit canopy components for the sensors due to decreasing SZA from June (44°) to September (26°). Because EVI was highly correlated with the reflectance of the NIR band used to generate it (r of + 0.98 for MODIS and + 0.88 for Hyperion), this vegetation index followed the general NIR pattern, increasing with smaller SZA towards the end of the dry season. Hyperion EVI was inversely correlated with the shade fraction (r = − 0.93). Changes in canopy foliage detected from MODIS LAI data were not consistent with LAI estimates from hemispherical photographs. Although further research is necessary to measure the impact of leaf flush on intra-annual EVI variability in the Querência region, analysis of Hymap data with fixed LAI and leaf flush conditions confirmed the influence of the illumination effects on the EVI.     

Tree density estimation in a tropical woodland ecosystem with multiangular MISR and MODIS data

In this paper we evaluate the potential of spectral, temporal and angular aspect of remotely sensed data for quantitative extraction of forest structure information in tropical woodlands. Moderate resolution imaging spectroradiometer (MODIS) multispectral data at 500-meter spatial resolution from different dates, multiangle imaging spectroradiometer (MISR) bidirectional reflectance factors (BRF) and normalized difference angular index (NDAI) derived from MISR data at 275-meter spatial resolution were used as input data. The number of trees per hectare bigger than 20cm in diameter at breast height was taken as variable of interest. Simple and multiple ordinary least square regressions and artificial neural networks (ANN) were tested to understand the relationships between the various sources of remotely sensed data and the output variable. An experimental design technique, followed by a classification of the input variables and a factor analysis were implemented in order to understand the structure, reduce the dimensionality of the data and avoid the overfitting of the neural network. The results show that there is a significant amount of independent information in the angular dimension, and this information is highly relevant to the estimation of tree densities in the study area. The MISR NDAI indexes improved the performance of the MISR BRF. The non-linear ANN outperformed the linear regressions. The best results were obtained with the ANN after selecting the input variables according to the results of the experimental design, the classification and the factor analysis, with a 0.71 correlation coefficient against the 0.58 of the best linear regression model.     

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%.     

Assessing the sensitivity of MODIS to monitor drought in high biomass ecosystems

Drought monitoring is important to analyse the influence of rainfall deficiency patterns on bushfire behaviour. Remote sensing provides tools for spatially explicit monitoring of drought across large areas. The objective of this study was to assess the performance of MODIS-based reflectance spectral indices to monitor drought across forest and woodland vegetation types in the fire prone Sydney Basin Bioregion, NSW, Australia. A time series of eight spectral indices were created from 2000 to 2009 to monitor inter-annual changes in drought and were compared to the Standardized Precipitation Index (SPI), a precipitation deficit/surplus indicator. A pixel-to-weather station paired correlation approach was used to assess the relationship between SPI and the MODIS-based spectral indices at different time scales. Results show that the Normalised Difference Infrared Index—band 6 (NDIIb6) provided the most suitable indicator of drought for the high biomass vegetation types considered. The NDIIb6 had the highest sensitivity to drought intensity and was highly correlated with SPI at all time scales analysed (i.e., 1, 3 and 6-month SPI) suggesting that variations in precipitation patterns have a stronger influence on vegetation water content than vegetation greenness properties. Spatial similarities were also found between patterns of NDIIb6-based drought maps and SPI values distribution. NDIIb6 outperformed the spectral index currently in use for operational drought monitoring systems in the region (Normalised Difference Vegetation Index, NDVI) and its implementation in existing drought-monitoring systems is recommended.     

利用MODIS卫星资料分析北京地区地表反照率时空分布及变化特征

利用2000～ 2005 年MOD IS 地表双向反照率(MOD43B3) , 植被指数(MOD13A 2) 和地表覆盖类型(MOD12Q 1) 资料, 分析了北京及周边地区地表反照率的时空分布, 计算了2000～ 2005 年平均地表反照率的时间变化。结果表明, 北京城区年平均地表反照率为0. 12, 山区森林(0. 11) 明显小于平原地区(0. 15) , 而永定河流域反照率较大(0. 18) , 这主要是因为永定河流域植被覆盖度较低(植被指数低) , 山区地表反照率季节性变化不明显。2000～ 2005 年6 年的统计回归显示北京平原地区反照率呈略有下降。