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.     

Cross-scalar satellite phenology from ground, Landsat, and MODIS data

Phenological records constructed from global mapping satellite platforms (e.g. AVHRR and MODIS) hold the potential to be valuable tools for monitoring vegetation response to global climate change. However, most satellite phenology products are not validated, and field checking coarse scale (≥ 500 m) data with confidence is a difficult endeavor. In this research, we compare phenology from Landsat (field scale, 30 m) to MODIS (500 m), and compare datasets derived from each instrument. Landsat and MODIS yield similar estimates of the start of greenness (r² = 0.60), although we find that a high degree of spatial phenological variability within coarser-scale MODIS pixels may be the cause of the remaining uncertainty. In addition, spatial variability is smoothed in MODIS, a potential source of error when comparing in situ or climate data to satellite phenology. We show that our method for deriving phenology from satellite data generates spatially coherent interannual phenology departures in MODIS data. We test these estimates from 2000 to 2005 against long-term records from Harvard Forest (Massachusetts) and Hubbard Brook (New Hampshire) Experimental Forests. MODIS successfully predicts 86% of the variance at Harvard forest and 70% of the variance at Hubbard Brook; the more extreme topography of the later is inferred to be a significant source of error. In both analyses, the satellite estimate is significantly dampened from the ground-based observations, suggesting systematic error (slopes of 0.56 and 0.63, respectively). The satellite data effectively estimates interannual phenology at two relatively simple deciduous forest sites and is internally consistent, even with changing spatial scale. We propose that continued analyses of interannual phenology will be an effective tool for monitoring native forest responses to global-scale climate variability.     

Evaluation of optical satellite remote sensing for rice paddy phenology in monsoon Asia using a continuous in situ dataset

In monsoon Asia, optical satellite remote sensing for rice paddy phenology suffers from atmospheric contaminations mainly due to frequent cloud cover. We evaluated the quality of satellite remote sensing of paddy phenology: (1) through continuous in situ observations of a paddy field in Japan for 1.5 years, we investigated phenological signals in the reflectance spectrum of the paddy field; (2) we tested daily satellite data taken by Terra/Aqua MODIS (MOD09 and L1B products) with regard to the agreement with the in situ data and the influence of cloud contamination. As a result, the in situ spectral characteristics evidently indicated some phenological changes in the rice paddy field, such as irrigation start, padding, heading, harvest and ploughing. The Enhanced Vegetation Index (EVI) was the best vegetation index in terms of agreement with the in situ data. More than 65% of MODIS observations were contaminated with clouds in this region. However, the combined use of Terra and Aqua decreased the rate of cloud contamination of the daily data to 43%. In conclusion, the most robust dataset for monitoring rice paddy phenology in monsoon Asia would be daily EVI derived from a combination of Terra/MODIS and Aqua/MODIS.     

Sensitivity of vegetation phenology detection to the temporal resolution of satellite data

Vegetation phenology derived from satellite data has increasingly received attention for applications in environmental monitoring and modelling. The accuracy of phenological estimates, however, is unknown at the regional and global level because field validation data are insufficient. To assess the accuracy of satellite‐derived phenology, this study investigates the sensitivity of phenology detection to both the temporal resolution of sampling and the number of consecutive missing values (usually representing cloud cover) in the time series of satellite data. To do this, time series of daily vegetation index data for various ecosystems are modelled and simulated using data from Moderate‐Resolution Imaging Spectroradiometer (MODIS) data. The annual temporal data are then fitted using piecewise logistic functions, which are employed to calculate curvature change rate for detecting phenological transition dates. The results show that vegetation phenology can be estimated with a high precision from time series with temporal resolutions of 6–16 days even if daily data contains some uncertainties. If the temporal resolution is no coarser than 16 days for time series sampled using an average composite, the absolute errors are less than 3 days. On the other hand, the phase shift of temporal sampling is shown to have limited impacts on phenology detection. However, the accuracy of phenology detection may be reduced greatly if missing values in the time series of 16‐day MODIS data occur around the onsets of phenological transition dates. Even so, the probability that the absolute error in phenological estimates is greater than 5 days is less than 4% when only one period is missing in the time series of 16‐day data during vegetation growing seasons; this probability increases to 20% if there are two consecutive missing values.     

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.     

Assessing viewing and illumination geometry effects on the MODIS vegetation index (MOD13Q1) time series: implications for monitoring phenology and disturbances in forest communities in Queensland,Australia

Time series analysis of satellite data can be used to monitor temporal dynamics of forested environments, thus providing spatial data for a range of forest science, monitoring and management issues. The moderate resolution imaging spectroradiometer (MODIS) vegetation index (MOD13Q1) product has potential for monitoring forest dynamics and disturbances. However, the suitability of the product to accurately measure temporal changes due to phenology and disturbances is questionable as the effects of variable solar and viewing geometry have not been removed from these data. This study aimed to examine the impact that viewing and illumination geometry differences had on MOD13Q1 vegetation index values, and their subsequent ability to map changes arising from phenology and disturbances in a number of forest communities in Queensland, Australia. MOD13Q1 normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) were compared to normalized NDVI and EVI (NDVI_normalized and EVI_normalized), which were derived from the reflectance modelled from a bidirectional reflectance distribution function (BRDF)/albedo parameters product (MCD43A1) using fixed viewing and illumination geometry. Time series plots of the vegetation index values from a number of pixels representing different forest types and known disturbances showed that the NDVI_normalized time series was more effective at capturing the changes in vegetation than the NDVI. MOD13Q1 NDVI showed higher seasonal amplitude, but was less accurate at capturing phenology and disturbances compared to the NDVI_normalized. The EVI was less affected by variable viewing and illumination geometry in terms of amplitude, but was affected in terms of time shift in periodicities providing erroneous information on phenology. More studies in different environments with appropriate vegetation phenology reference data will be needed to confirm these observations.     

Phenology of vegetation in Southern England from Envisat MERIS terrestrial chlorophyll index (MTCI) data

Given the close association between climate change and vegetation response, there is a pressing requirement to monitor the phenology of vegetation and understand further how its metrics vary over space and time. This article explores the use of the Envisat MERIS terrestrial chlorophyll index (MTCI) data set for monitoring vegetation phenology, via its estimates of chlorophyll content. The MTCI was used to construct the phenological profile of and extract key phenological event dates from woodland and grass/heath land in Southern England as these represented a range of chlorophyll contents and different phenological cycles. The period 2003–2008 was selected as this was known to be a period with temperature and phenological anomalies. Comparisons of the MTCI-derived phenology data were made with ground indicators and climatic proxy of phenology and with other vegetation indices: MERIS global vegetation index (MGVI), MODIS normalized difference vegetation index (NDVI) and MODIS enhanced vegetation index (EVI). Close correspondence between MTCI and canopy phenology as indicated by ground observations and climatic proxy was evident. Also observed was a difference between MTCI-derived phenological profile curves and key event dates (e.g. green-up, season length) and those derived from MERIS MGVI, MODIS NDVI and MODIS EVI. The research presented in this article supports the use of the Envisat MTCI for monitoring vegetation phenology, principally due to its sensitivity to canopy chlorophyll content, a vegetation property that is a useful proxy for the canopy physical and chemical alterations associated with phenological change.     

基于AVHRR和MODIS数据的全球植被物候比较分析

AVHRR和MODIS卫星数据在全球和区域尺度植物物候对气候变化响应研究中起着重要的作用,然而两种传感器在全球尺度物候监测的一致性有待验证。首先利用时间序列谐波分析法(HANTS)对2005年全球GIMMS AVHRR NDVI和MODIS 13A2 数据进行滤波处理;然后基于改进的动态阈值方法,提取全球植被的返青期(SOS)、枯黄期(EOS)和生长季长度(DOS);最后分区域比较和评估两种传感器提取物候参数的潜力。研究结果表明：2005年全球大部分地区植被在第100～140 d开始生长,到第260～300 d逐渐停止生长,生长季长度集中在130～180 d,并且和区域研究结果具有一致性;两种传感器提取的植被关键物候期的空间变化趋势是一致的,随着纬度升高,返青期呈现推迟趋势,枯黄期呈现提早趋势,生长季长度呈现缩短趋势;AVHRR和MODIS提取落叶林和草地的SOS、EOS和DOS在欧亚大陆和北美洲区域的相关系数大部分达到0.9以上。     

Comparison of multiple models for remote sensing of carbon exchange using MODIS data in conifer-dominated forests

Boreal forests in the northern hemisphere provide important sinks for storing carbon dioxide (CO₂). However, the size and distribution of these sinks remain uncertain. In particular, many remote-sensing models show a strong bias in the simulation of carbon fluxes for evergreen needleleaf forest. The objective of this study is to improve these predictive models for accurately quantifying temporal changes in the net ecosystem exchange (NEE) of conifer-dominated forest solely based on satellite remote sensing, including the Moderate Resolution Imaging Spectroradiometer (MODIS) Terra daytime land-surface temperature (LST), night-time LST′, enhanced vegetation index (EVI), land–surface water index (LSWI), fraction of absorbed photosynthetically active radiation (FPAR), and leaf area index (LAI). Considering that the component fluxes, gross primary production (GPP), and ecosystem respiration (R_e), are strongly influenced by vegetation phenology, seasonality information was extracted from time-series MODIS EVI data based on non-linear least-squares fits of asymmetric Gaussian model functions with a software package for analysing the time-series of satellite sensor data (TIMESAT). The results indicated that models directly incorporating phenological information failed to improve their performance for temperate deciduous forest. Instead, three methods to retrieve the component fluxes – GPP and R_e – including direct estimates, models incorporating the phenological information, and models developed based on the threshold value (LST 273 K), were explored respectively. All methods improved NEE estimates markedly and models developed based on the threshold value performed best, and provided a future framework for accurate remote sensing of NEE in evergreen forest.     

Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data

A semi-physical fusion approach that uses the MODIS BRDF/Albedo land surface characterization product and Landsat ETM+ data to predict ETM+ reflectance on the same, an antecedent, or subsequent date is presented. The method may be used for ETM+ cloud/cloud shadow and SLC-off gap filling and for relative radiometric normalization. It is demonstrated over three study sites, one in Africa and two in the U.S. (Oregon and Idaho) that were selected to encompass a range of land cover land use types and temporal variations in solar illumination, land cover, land use, and phenology. Specifically, the 30 m ETM+ spectral reflectance is predicted for a desired date as the product of observed ETM+ reflectance and the ratio of the 500 m surface reflectance modeled using the MODIS BRDF spectral model parameters and the sun-sensor geometry on the predicted and observed Landsat dates. The difference between the predicted and observed ETM+ reflectance (prediction residual) is compared with the difference between the ETM+ reflectance observed on the two dates (temporal residual) and with respect to the MODIS BRDF model parameter quality. For all three scenes, and all but the shortest wavelength band, the mean prediction residual is smaller than the mean temporal residual, by up to a factor of three. The accuracy is typically higher at ETM+ pixel locations where the MODIS BRDF model parameters are derived using the best quality inversions. The method is most accurate for the ETM+ near-infrared (NIR) band; mean NIR prediction residuals are 9%, 12% and 14% of the mean NIR scene reflectance of the African, Oregon and Idaho sites respectively. The developed fusion approach may be applied to any high spatial resolution satellite data, does not require any tuning parameters and so may be automated, is applied on a per-pixel basis and is unaffected by the presence of missing or contaminated neighboring Landsat pixels, accommodates for temporal variations due to surface changes (e.g., phenological, land cover/land use variations) observable at the 500 m MODIS BRDF/Albedo product resolution, and allows for future improvements through BRDF model refinement and error assessment.     

Green leaf phenology at Landsat resolution: Scaling from the field to the satellite

Despite the large number of in situ, plot-level phenological measurements and satellite-derived phenological studies, there has been little success to date in merging these records temporally or spatially. In this research, we bridge this scale gap through higher resolution satellite records (Landsat) and quantify the accuracy of satellite-derived metrics with direct field measurements. We compiled fifty-seven Landsat scenes from southern New England (P12 R51) from 1984 to 2002. Green vegetation areal abundance for each scene was derived from spectral mixture analysis and a single set of endmembers. The leaf area signal was fit with a logistic-growth simulating sigmoid curve to derive phenological markers (half-maximum leaf-onset and offset). Spring leaf-onset dates in homogenous stands of deciduous forests displayed significant and persistent local variability. The local variability was validated with multiple springtime ground observations (r² = 0.91). The highest degree of verified small-scale variation occurred where contiguous forests displayed leaf-onset gradients of 10-14 days over short distances (< 500 m). These dramatic gradients occur in of low-relief (< 40 m) upland regions. The patterns suggest that microclimates resulting from springtime cold-air drainage may be influential in governing the start of leaf growth; every 4.16 m loss in elevation delayed spring leaf onset by 1 day. These microclimates may be of crucial importance in interpreting in situ records and interpolating phenology from satellite data. Regional patterns from the Landsat analyses suggest topographic, coastal, and land-use controls on phenology. Our results indicate that deciduous forests in the Providence, RI metropolitan area leaf out 5-7 days earlier than comparable rural areas. The platform-independent curve-fit methodology may be extended across platforms and field data. The methodologically consistent approach, in tandem with Landsat data, allows an effective scaling from plot to satellite phenological observations.     

Interannual vegetation phenology estimates from global AVHRR measurements: Comparison with in situ data and applications

In a previous paper, we described a procedure to correct the directional effects in AVHRR reflectance time series. The corrected measurements show much less high frequency variability than their original counterparts, which makes them suitable to study vegetation dynamics. The time series are used here to estimate the start and ending dates of the growing season for 18 years from 1982 to 1999. We focus on the interannual variations of these phenological parameters.A database of in situ phenology observations is used to quantify the accuracy of the satellite-based estimates. Although based on a limited sampling of the Northern mid and high latitudes, the comparison indicates that i) the satellite phenological product contains meaningful information on interannual onset anomalies; ii) there is a higher degree of consistency over regions covered by Broadleaf Forests, Grasses and cereal Crops than over those covered by Needleleaf Forests or Savannas; and iii) the satellite phenological product is of lower quality in regions with mountainous terrain. In favorable conditions, interannual variations of the onset are captured with an accuracy of a few days.As this satellite-derived product captures the interannual onset variability at ground-truth sites, we confidently use it to larger scales studies. Mapped at a continental scale, the onset anomalies show coherent patterns at the regional (≈ 1000 km) scale for the mid and high latitudes of the Northern hemisphere, which is consistent with a meteorological forcing. In the tropics, there is larger spatial heterogeneity, which suggests more complex controls of the phenology. The relation between vegetation phenology and climate is further investigated over Europe by comparing the variability of the satellite-derived vegetation onset and that of the winter North Atlantic Oscillation index, at a fine spatial scale. The strong correlations observed confirm that this climate forcing parameter explains most of the onset variability over a large fraction of Northern Europe (earlier onsets for positive winter NAO), with lower impact towards the south and opposite effects around the Mediterranean basin. The NAO has a predictive character as the January-February NAO index is strongly correlated with the vegetation onset that occurs around April in Northern Europe.     

Monitoring boreal forest leaf area index across a Siberian burn chronosequence: a MODIS validation study

Landscapes containing differing amounts of ecological disturbance provide an excellent opportunity to validate and better understand the emerging Moderate Resolution Imaging Spectrometer (MODIS) vegetation products. Four sites, including 1‐year post‐fire coniferous, 13‐year post‐fire deciduous, 24‐year post‐fire deciduous, and >100 year old post‐fire coniferous forests, were selected to serve as a post‐fire chronosequence in the central Siberian region of Krasnoyarsk (57.3°N, 91.6°E) with which to study the MODIS leaf area index (LAI) and vegetation index (VI) products. The collection 4 MODIS LAI product correctly represented the summer site phenologies, but significantly underestimated the LAI value of the >100 year old coniferous forest during the November to April time period. Landsat 7‐derived enhanced vegetation index (EVI) performed better than normalized difference vegetation index (NDVI) to separate the deciduous and conifer forests, and both indices contained significant correlation with field‐derived LAI values at coniferous forest sites (r ² = 0.61 and r ² = 0.69, respectively). The reduced simple ratio (RSR) markedly improved LAI prediction from satellite measurements (r ² = 0.89) relative to NDVI and EVI. LAI estimates derived from ETM+ images were scaled up to evaluate the 1 km resolution MODIS LAI product; from this analysis MODIS LAI overestimated values in the low LAI deciduous forests (where LAI<5) and underestimated values in the high LAI conifer forests (where LAI>6). Our results indicate that further research on the MODIS LAI product is warranted to better understand and improve remote LAI quantification in disturbed forest landscapes over the course of the year.     

Comparison of multiple models for estimating gross primary production using MODIS and eddy covariance data in Harvard Forest

Gross primary production (GPP) defined as the overall rate of fixation of carbon through the process of vegetation photosynthesis is important for carbon cycle and climate change research. Three models, the Vegetation Photosynthesis Model (VPM), the Temperature and Greenness (TG) model and the Vegetation Index (VI) model have been compared for the estimation of GPP in Harvard Forest from 2003 to 2006 using climate variables acquired by eddy covariance (EC) measurements and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite images. All these models provide more reliable estimates of GPP than that of MODIS GPP product. High Pearsons correlation coefficients r equal to 0.94, 0.92 and 0.90 are observed for the VPM, the TG and the VI model, respectively. Relationships between GPP and land surface temperature (LST, R² = 0.72), and vapor pressure deficit (VPD, R² = 0.45) indicate that climate variables are important for GPP estimation. Due to proper characterization of temperature, water stress and leaf age by three scalars, VPM best follows the seasonal variations of GPP. By incorporation of the MODIS surface reflectance and LST product, the TG model is the most suitable choice for areas without prior knowledge as it is based entirely on remote sensing observations. Results from the VI model demonstrate the possibility of using a single vegetation index for light use efficiency (LUE) estimation in deciduous forest that is of high spatial heterogeneity. The validation and comparison of models will be helpful in development of future GPP models using combinations of climate variables and/or remote sensing observations.     

Fire-induced changes in green-up and leaf maturity of the Canadian boreal forest

Recent studies of vegetation phenology of northern forests using satellite data suggest that the observed earlier spring increase and peak amplitude of the normalized difference vegetation index (NDVI) are a result of climate warming. In addition to undergoing an increase in temperature, the northern forests of Canada have also seen a dramatic increase in area burned by wildfire over the same time period. Using the Canadian Large Fire Database, we analyzed the impact fire had on the phenological dates derived from fitting a logistical model to yearly data from 2004 for several different subsets of both AVHRR-NDVI and MODIS LAI in wildfire dominated terrestrial ecozones. Fire had a significant but complex effect on estimated phenology dates. The most recently burned areas (1994–2003) had later green-up dates in two ecozones for AVHRR data and all ecozones for MODIS. However, older forested (not burned during 1980–2003) had estimated green-up dates 1 to 9 days earlier than the entire forested area in the MODIS LAI data. These data corroborate studies in Canada and demonstrate that fire history is influencing boreal forest phenology and growing season LAI.     

Feasibility of using Geostationary Ocean Colour Imager (GOCI) data for land applications after atmospheric correction and bidirectional reflectance distribution function modelling

Korea's Geostationary Ocean Colour Imager (GOCI) has very high temporal resolution as well as wide spatial coverage. There is thus great interest in testing its applicability for monitoring land areas in addition to ocean areas. GOCI has eight spectral bands, from blue to near-infrared. These bands can be sensitive to vegetation change, but their wavelength ranges are slightly different from those of the extensively studied Moderate Resolution Imaging Spectroradiometer (MODIS). This study examines whether GOCI data can be applied for land monitoring and how GOCI data should be processed so as to reflect the spectral characteristics of land surfaces as detected by polar-orbit satellite sensors. Several image processing steps were performed for the GOCI data, including atmospheric correction and semi-empirical bidirectional reflectance distribution function modelling, before the results were compared with the MODIS land-surface product. Among the four GOCI normalized difference vegetation index (NDVI) products tested in this study, the GOCI NDVI with viewing-angle-adjusted reflectance showed the best agreement with MODIS NDVI calculated from normalized reflectance, with the lowest root mean square error of 0.126. Additionally, its temporal trends over forest and mixed vegetation areas were similar to those of MODIS NDVI during the study period from September to December.     

Development of an approach for generation of temporally complete daily nadir MODIS reflectance time series

Consistent, spatially and temporally complete reflectance time series are required for reliable terrestrial monitoring. The Moderate Resolution Imaging Spectroradiometer (MODIS), like other polar-orbiting wide field of view satellite sensors, can provide global observations on a nearly daily basis, but the sparseness of valid observations due to cloud, residual atmospheric effects, and sensor anomalies, may result in gaps in the derived reflectance time series. This paper presents an approach for the generation of temporally complete daily MODIS 500 m nadir view BRDF-adjusted reflectance (NBAR) time series. The research is illustrated and assessed quantitatively using two years of cloud and snow screened, daily MODIS Terra and Aqua reflectance data at four sites in Africa, and demonstrated for phenology monitoring using NBAR derived NDVI time series. The components of the approach include: 1) an outlier detection algorithm to remove residual anomalous daily observations undetected in the upstream processing, 2) the dynamic generation of NBAR time series on a daily basis when seven or more observations are available for a day under consideration over a 16-day period, and 3) the means to gap fill the NBAR time series where less than 7 observations are available. The MODIS Ross-Thick/Li-Sparse-Reciprocal BRDF model is used with a rolling approach whereby a 16-day BRDF inversion window is moved on a daily overlapping basis to provide more reliable outlier detection and daily NBAR. NBAR gap filling in periods of missing observations is investigated using static land cover specific archetype BRDF parameters and using BRDF parameters defined adaptively from the temporally closest 16-day periods with 7 or more observations. Scaling factor estimators using ordinary least squares (OLS) and median-based robust least squares regression are investigated, and the robust method is demonstrated to provide on average temporally more coherent gap filled NBAR values. For regions with persistent clouds, the utility of the adaptive NBAR gap filling method is demonstrated to be severely limited due to the decreased likelihood that the surface BRDF at each gap can be described reliably. The reliability of the NBAR gap filling methodology is evaluated statistically using a cross-validation approach. For the small number of study site considered, the adaptive method is shown to provide more accurate results than the archetype method when there are more than an average of ~ 4-5 observations per 16-day window, or when a gap day is on average less than about 30 days from a 16-day period with 7 or more observations. The resulting gap free daily NBAR time series and derived daily NBAR NDVI generated by the approach is shown to capture phenological variations in a coherent temporally consistent manner, suggesting that it is a fruitful avenue for future research and validation.     

A comparison of Landsat TM and MODIS vegetation indices for estimating forage phenology in desert bighorn sheep (Ovis canadensis nelsoni) habitat in the Sonoran Desert,USA

Sonoran Desert bighorn sheep (Ovis canadensis mexicana) occupy rugged upland areas that experience irregular periods of vegetation growth associated with precipitation events. These episodic and often spatially limited events provide important forage and preformed water resources that may be important drivers of animal movement and habitat use. Habitat-use models that incorporate forage phenology would broaden our understanding of desert bighorn ecology and have considerable potential to inform conservation efforts for the species. Field-based methods are of limited utility to characterize vegetation phenology across large areas. Vegetation indices (VI) derived from satellite imagery are a viable alternative, but may be confounded by areas of high relief and shadow effects that can degrade VI values. The varying spatial and temporal resolutions of readily available satellite sensors, such as the Landsat thematic mapper (TM) and moderate-resolution imaging spectrometer (MODIS), present additional challenges. In this study, we sought to minimize degrading effects of terrain on TM- and MODIS-based estimates of vegetation phenology. We compared effects of high topographic relief on time series MODIS- and TM-based VI such as the normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) using VI departures from average (DA) in shaded and unshaded areas. Sun elevation angle negatively impacted TM-derived NDVI and EVI values in areas of steep terrain. In contrast, MODIS-derived NDVI values were insensitive to sun elevation and terrain effects, whereas MODIS-derived EVI was degraded in areas of steep terrain. Time series MODIS NDVI and EVI DA values differed significantly during months of low sun elevation angle. Average MODIS EVI departure values were ≥20% lower than NDVI under these conditions, confounding time series estimates of plant phenology. Our best results were obtained from MODIS 16-day composited NDVI. These remote-sensing-based VI estimates of seasonal plant phenology and productivity can be used to inform models of habitat use and movements of desert bighorn over large areas.     

Characterization of seasonal variation of forest canopy in a temperate deciduous broadleaf forest, using daily MODIS data

In this paper, we present an improved procedure for collecting no or little atmosphere- and snow-contaminated observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. The resultant time series of daily MODIS data of a temperate deciduous broadleaf forest (the Bartlett Experimental Forest) in 2004 show strong seasonal dynamics of surface reflectance of green, near infrared and shortwave infrared bands, and clearly delineate leaf phenology and length of plant growing season. We also estimate the fractions of photosynthetically active radiation (PAR) absorbed by vegetation canopy (FAPAR_canopy), leaf (FAPAR_leaf), and chlorophyll (FAPAR_chl), respectively, using a coupled leaf-canopy radiative transfer model (PROSAIL-2) and daily MODIS data. The Markov Chain Monte Carlo (MCMC) method (the Metropolis algorithm) is used for model inversion, which provides probability distributions of the retrieved variables. A two-step procedure is used to estimate the fractions of absorbed PAR: (1) to retrieve biophysical and biochemical variables from MODIS images using the PROSAIL-2 model; and (2) to calculate the fractions with the estimated model variables from the first step. Inversion and forward simulations of the PROSAIL-2 model are carried out for the temperate deciduous broadleaf forest during day of year (DOY) 184 to 201 in 2005. The reproduced reflectance values from the PROSAIL-2 model agree well with the observed MODIS reflectance for the five spectral bands (green, red, NIR₁, NIR₂, and SWIR₁). The estimated leaf area index, leaf dry matter, leaf chlorophyll content and FAPAR_canopy values are close to field measurements at the site. The results also showed significant differences between FAPAR_canopy and FAPAR_chl at the site. Our results show that MODIS imagery provides important information on biophysical and biochemical variables at both leaf and canopy levels.     

Reflectance seasonality and its relation to the canopy leaf area index in an eastern Siberian larch forest: Multi-satellite data and radiative transfer analyses

Reliable monitoring of seasonality in the forest canopy leaf area index (LAI) in Siberian forests is required to advance the understanding of climate-forest interactions under global environmental change and to develop a forest phenology model within ecosystem modeling. Here, we compare multi-satellite (AVHRR, MODIS, and SPOT/VEGETATION) reflectance, normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and LAI with aircraft-based spectral reflectance data and field-measured forest data acquired from April to June in 2000 in a larch forest near Yakutsk, Russia. Field data in a 30 × 30-m study site and aircraft data observed around the field site were used. Larch is a dominant forest type in eastern Siberia, but comparison studies that consider multi-satellite data, aircraft-based reflectance, and field-based measurement data are rarely conducted. Three-dimensional canopy radiative transfer calculations, which are based on Antyufeev and Marshak's [Antyufeev, V.S., & Marshak, A.L. (1990). Monte Carlo method and transport equation in plant canopies, Remote Sensing of Environment, 31, 183-191] Monte Carlo photon transport method combined with North's [North, P.R. (1996). Three-dimensional forest light interaction model using a Monte Carlo method, IEEE Transactions on Geoscience and Remote Sensing, 34(4), 946-956] geometric-optical hybrid forest canopy scene, helped elucidate the relationship between canopy reflectance and forest structural parameters, including several forest floor conditions. Aircraft-based spectral measurements and the spectral response functions of all satellite sensors confirmed that biases in reflectance seasonality caused by differences in spectral response functions among sensors were small. However, some reflectance biases occur among the near infrared (NIR) reflectance data from satellite products; these biases were potentially caused by absolute calibration errors or cloud/cloud shadow contamination. In addition, reflectance seasonality in AVHRR-based NIR data was very small compared to other datasets, which was partially due to the spring-to-summer increase in the amount of atmospheric water vapor. Radiative transfer simulations suggest that bi-directional reflectance effects were small for the study site and observation period; however, changes in tree density and forest floor conditions affect the absolute value of NIR reflectance, even if the canopy leaf area condition does not change. Reliable monitoring of canopy LAI is achieved by minimizing these effects through the use of NIR reflectance difference, i.e., the difference in reflectance on the observation day from the reflectance on a snow-free/pre-foliation day. This may yield useful and robust parameters for multi-satellite monitoring of the larch canopy LAI with less error from intersensor biases and forest structure/floor differences. Further validation with field data and combined use of other index (e.g. normalized difference water index, NDWI) data will enable an extension of these findings to all Siberian deciduous forests.