The effect of the temporal resolution of NDVI data on season onset dates and trends across Canadian broadleaf forests

Satellite remote sensing has the potential to contribute to plant phenology monitoring at spatial and temporal scales relevant for regional and global scale studies. Historically, temporal composites of satellite data, ranging from 8 days to 16 days, have been used as a starting point for satellite-derived phenology data sets. In this study we assess how the temporal resolution of such composites affects the estimation of the start of season (SOS) by: 1) calibrating a relationship between satellite derived SOS with in situ leaf unfolding (LU) of trembling aspen (Populus tremuloides) across Canada and 2) quantifying the sensitivity of calibrated satellite SOS estimates and trends, over Canadian broadleaf forests, to the temporal resolution of NDVI data. SOS estimates and trends derived from daily NDVI data were compared to SOS estimates and trends derived from multiday NDVI composites that retain the exact date of the maximum NDVI value or that assume the midpoint of the multiday interval as the observation date. In situ observations of LU dates were acquired from the PlantWatch Canada network. A new Canadian database of cloud and snow screened daily 1-km resolution National Oceanic and Atmospheric Administration advanced very high resolution radiometer surface reflectance images was used as input satellite data. The mean absolute errors of SOS dates with respect to in situ LU dates ranged between 13 and 40 days. SOS estimates from NDVI composites that retain the exact date of the maximum NDVI value had smaller errors (~ 13 to 20 days). The sensitivity analysis reinforced these findings: SOS estimates from NDVI composites that use the exact date had smaller absolute deviations from the LU date (0 to − 5 days) than the SOS estimates from NDVI composites that use the midpoint (− 2 to − 27 days). The SOS trends between 1985 and 2007 were not sensitive to the temporal resolution or compositing methods. However, SOS trends at individual ecozones showed significant differences with the SOS trends from daily NDVI data (Taiga plains and the Pacific maritime zones). Overall, our results suggest that satellite based estimates of vegetation green-up dates should preferably use sub-sampled NDVI composites that include the exact observation date of the maximum NDVI to minimize errors in both, SOS estimates and SOS trend analyses. For trend analyses alone, any of the compositing methods could be used, preferably with composite intervals of less than 28 days. This is an important finding, as it suggests that existing long-term 10-day or 15-day NDVI composites could be used for SOS trend analyses over broadleaf forests in Canada or similar areas. Future studies will take advantage of the growing in situ phenology networks to improve the validation of satellite derived green-up dates.     

Satellite passive microwave remote sensing for monitoring global land surface phenology

Vegetation phenology characterizes seasonal life-cycle events that influence the carbon cycle and land-atmosphere water and energy exchange. We analyzed global phenology cycles over a six year record (2003-2008) using satellite passive microwave remote sensing based Vegetation Optical Depth (VOD) retrievals derived from daily time series brightness temperature (T_b) measurements from the Advanced Microwave Scanning Radiometer on EOS (AMSR-E) and other ancillary data inputs. The VOD parameter derives vegetation canopy attenuation at a given microwave frequency (18.7 GHz) and varies with canopy height, density, structure and water content. An error sensitivity analysis indicates that the retrieval algorithm can resolve the VOD seasonal cycle over a majority of global vegetated land areas. The VOD results corresponded favorably (p < 0.01) with vegetation indices (VIs) and leaf area index (LAI) information from satellite optical-infrared (MODIS) remote sensing, and phenology cycles determined from a simple bioclimatic growing season index (GSI) for over 82% of the global domain. Lower biomass land cover classes (e.g. savannas) show the highest correlations (R = 0.66), with reduced correspondence at higher biomass levels (0.03 < R < 0.51) and higher correlations for homogeneous land cover areas (0.41 < R < 0.83). The VOD results display a unique end-of-season signal relative to VI and LAI series, and may reflect microwave sensitivity to the timing of vegetation biomass depletion (e.g. leaf abscission) and associated changes in canopy water content (e.g. dormancy preparation). The VOD parameter is independent of and synergistic with optical-infrared remote sensing based vegetation metrics, and contributes to a more comprehensive view of land surface phenology.     

Evaluation of compositing period and AVHRR and MERIS combination for improvement of spring phenology detection in deciduous forests

Accurate and precise detection of phenology events is needed to assess trends in seasonal vegetation development indicative of climate or other environmental change processes. In this research, detection accuracy of start of season (SOS) phenology for deciduous forest across Eastern Canada was assessed using satellite time series and in situ PlantWatch observations. Several aspects were evaluated regarding performance of phenology information extraction: 1) effect of compositing period, 2) individual performance of the Advanced Very High Resolution Radiometer (AVHRR) and the Medium Resolution Imaging Spectrometer (MERIS) sensors, and 3) performance for these sensors combined. The AVHRR and MERIS sensors were used as they are overlapping operational missions with planned future continuity. Three approaches to utilizing the multi-sensor data were tested: 1) inter-calibrating NDVI data between sensors and using the multi-sensor data stream to detect SOS, 2) combining independently derived SOS estimates from AVHRR and MERIS based on a weighted average, and 3) combining approaches 1 and 2. Comparison with in situ observations of leaf out and first bloom showed that combining independent SOS estimates from AVHRR and MERIS was better than using the inter-calibrated multi-sensor data. Combining SOS estimates from both sensors reduced error by 1-2 days compared to the individual sensor results. Composite periods from 7 to 11 days produced the best results for leaf out with a mean absolute error (MAE) of 5 days. Results for first bloom were not as good as those for leaf out, producing a MAE of 6.5 days. For first bloom, compositing periods greater than 11 days did not increase error at the same rate as seen for leaf out. However, the larger MAE observed for first bloom may have masked this effect.     

基于Landsat8热红外遥感数据的山地地表温度地形效应研究

地表温度是影响地表能量收支平衡的重要参量,能够综合反演地表的水热交换过程。虽然当前在基于地表温度开展全球或者区域尺度的地表能量平衡研究方面取得一系列的进展,但是面向山地区域尺度的类似研究仍然面临较大的挑战。为分析山地复杂地形对山地地表温度时空分布的影响规律,基于具有较高空间分辨率的Landsat 8热红外数据,以我国西南典型山地为研究对象,定量反演该区域的地表温度空间分布状况,结合SRTM90DEM数据,选择从海拔、坡度和坡向3个关键地形因子角度分析山地地表温度的地形效应特征。结果发现:山地地表温度随地形因子均呈现出十分显著的变化特征。总体而言,地表温度均随着海拔和坡度的升高而降低,而在坡向方面,南坡的温度相比北坡的温度要高。在地形效应分析的基础上,通过开展1km空间尺度地形和地表温度的空间统计分析发现,山地1km尺度下地表温度存在较大的空间异质性,且其影响不可忽略。研究结果表明:开展山地地表水热过程遥感动态监测需高空间分辨率地表温度作为数据支持,以准确描述山地地形因素对地表能量交换过程的影响。     

Monitoring elevation variations in leaf phenology of deciduous broadleaf forests from SPOT/VEGETATION time-series

In mountain forest ecosystems where elevation gradients are prominent, temperature gradient-based phenological variability can be high. However, there are few studies that assess the capability of remote sensing observations to monitor ecosystem phenology along elevation gradients, despite their relevance under climate change. We investigated the potential of medium resolution remotely sensed data to monitor the elevation variations in the seasonal dynamics of a temperate deciduous broadleaf forested ecosystem. Further, we explored the impact of elevation on the onset of spring leafing. This study was based on the analysis of multi-annual time-series of VEGETATION data acquired over the French Pyrenees Mountain Region (FPMR), in conjunction with simultaneous ground-based observations of leaf phenology made for two dominant tree species in the region (oak and beech). The seasonal variations in the perpendicular vegetation index (PVI) were analyzed during a five-year period (2002 to 2006). The five years of data were averaged into a one sole year in order to fill the numerous large spatio-temporal gaps due to cloud and snow presence - frequent in mountains - without altering the temporal resolution. Since a VEGETATION pixel (1 km²) includes several types of land cover, the broadleaf forest-specific seasonal dynamics of PVI was reconstructed pixel-by-pixel using a temporal unmixing method based on a non-parametric statistical approach. The spatial pattern of the seasonal response of PVI was clearly consistent with the relief. Nevertheless the elevational or geographic range of tree species, which differ in their phenology sensitivity to temperature, also has a significant impact on this pattern. The reduction in the growing season length with elevation was clearly observable from the delay in the increase of PVI in spring and from the advance of its decrease in the fall. The elevation variations in leaf flushing timing were estimated from the temporal change in PVI in spring over the study area. They were found to be consistent with those measured in situ (R2 > 0.95). It was deduced that, over FPMR, the mean delay of leaf flushing timing for every 100 m increase in elevation was estimated be approximately 2.3 days. The expected estimation error of satellite-based leaf unfolding date for a given elevation was approximately 2 days. This accuracy can be considered as satisfactory since it would allow us to detect changes in leafing timing of deciduous broadleaf forests with a magnitude equivalent to that due to an elevation variation of 100 m (2.3 days on average), or in other words, to that caused by a variation in the mean annual air temperature of 0.5 °C. Although averaging the VEGETATION data over five years led to a loss of interannual information, it was found to be a robust approach to characterise the elevation variations in spring leafing and its long-term trends.     

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.     

The effects of aggregated land cover data on estimating NPP in northern Wisconsin

Ecosystem models are routinely used to estimate net primary production (NPP) from the stand to global scales. Complex ecosystem models, implemented at small scales (< 10 km²), are impractical at global scales and, therefore, require simplifying logic based on key ecological first principles and model drivers derived from remotely sensed data. There is a need for an improved understanding of the factors that influence the variability of NPP model estimates at different scales so we can improve the accuracy of NPP estimates at the global scale. The objective of this study was to examine the effects of using leaf area index (LAI) and three different aggregated land cover classification products-two factors derived from remotely sensed data and strongly affect NPP estimates-in a light use efficiency (LUE) model to estimate NPP in a heterogeneous temperate forest landscape in northern Wisconsin, USA. Three separate land cover classifications were derived from three different remote sensors with spatial resolutions of 15, 30, and 1000 m. Average modeled net primary production (NPP) ranged from 402 gC m^− 2 year^− 1 (15 m data) to 431 gC m^− 2 year^− 1 (1000 m data), for a maximum difference of 7%. Almost 50% of the difference was attributed each to LAI estimates and land cover classifications between the fine and coarse scale NPP estimate. Results from this study suggest that ecosystem models that use biome-level land cover classifications with associated LUE coefficients may be used to model NPP in heterogeneous land cover areas dominated by cover types with similar NPP. However, more research is needed to examine scaling errors in other heterogeneous areas and NPP errors associated with deriving LAI estimates.     

Temporal and spatial dynamics of phenology along the North–South Transect of Northeast Asia

Vegetation phenology is sensitive to climate change and, as such, is often regarded as an indicator of climate change. It is a common practice to extract vegetation phenological indicators based on satellite remote sensing data. In this study, we used the Advanced Very High Resolution Radiometer (AVHRR) Global Inventory Modeling and Mapping Study (GIMMS) Third-Generation normalized difference vegetation index (NDVI_3G) to investigate temporal and spatial changes in phenology in Northeast Asia. Based on the maximum rate of change in the NDVI and dynamic threshold, we used the Asymmetric Gaussian model, Double Logistic method, and Savitzky-Golay filter to extract the start of the growing season (SOS), the end of the growing season (EOS), and the length of the growing season (LOS), respectively, along the North–South Transect of Northeast Asia (NSTNEA) from 1982 to 2014. We then compared the differences in SOS, EOS, and LOS and considered their spatio-temporal dynamics and relationship with temperature. The results show that the Asymmetric Gaussian model has the highest stability among the three methods. Dynamic thresholds corresponding to the maximum change rate of NDVI were mainly between 0.5 and 0.6. From 1982 to 2014, the SOS in the NSTNEA region occurred approximately 0.19 days earlier each year; the trends in EOS and LOS were not significant. In general, temperature and latitude have a strong linear relationship, both of which significantly impact vegetation phenology in the NSTNEA region. In addition, elevation also significantly impacts on vegetation phenology in the NSTNEA region.     

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.     

Extracting grassland vegetation phenology in North China based on cumulative SPOT-VEGETATION NDVI data

Plant phenology is one of the main indicators of climate or other environmental processes. This paper assesses the detection accuracy of start of season (SOS) and end of season (EOS) for grassland vegetation in north China from 2001 to 2010 using SPOT-VEGETATION normalized difference vegetation index (NDVI) data sets and in situ observations. The cumulative NDVI is calculated and fitted using a logistic model to identify phenological transition dates. The curvature of the fitted logistic models predicts phenological transition dates that correspond to the times at which the curvature in the yearly integrated NDVI exhibits local minimums or maximums. Validating with in situ observations, phenological dates are extracted from satellite time series data and are accurate to within 10 days. The spatial trends of SOS and EOS are very similar for 2001–2010. SOS mainly occurs from the day of year (DOY) 110 to DOY 170, and EOS occurs from DOY 240 to DOY 300. SOS displays a marked delay from south to north, while EOS gradually advances, indicating regional differences in climate and terrain. However, the effect of latitude and longitude on the average EOS of alpine grasslands is not significantly different, while SOS at low latitude and high longitude is 10 days earlier than at high-latitude and high-longitude regions. We detected an overall advance in SOS of 3.1 days over 10 years, and a 1.3-day delay in EOS. However, the amplitude is low (about 5 days) and the changes in most regions are not significant (close to zero). The results in this paper are concordant with many reported studies that explored the phenology of grasslands in North China, which is an important component of global grasslands science.     

Detection of foliage conditions and disturbance from multi-angular high spectral resolution remote sensing

Disturbance of forest ecosystems, an important component of the terrestrial carbon cycle, has become a focus of research over recent years, as global warming is about to increase the frequency and severity of natural disturbance events. Remote sensing offers unique opportunities for detection of forest disturbance at multiple scales; however, spatially and temporally continuous mapping of non-stand replacing disturbance remains challenging. First, most high spatial resolution satellite sensors have relatively broad spectral ranges with bandwidths unsuitable for detection of subtle, stress induced, features in canopy reflectance. Second, directional and background reflectance effects, induced by the interactions between the sun-sensor geometry and the observed canopy surface, make up-scaling of empirically derived relationships between changes in spectral reflectance and vegetation conditions difficult. Using an automated tower based spectroradiometer, we analyse the interactions between canopy level reflectance and different stages of disturbance occurring in a mountain pine beetle infested lodgepole pine stand in northern interior British Columbia, Canada, during the 2007 growing season. Directional reflectance effects were modelled using a bidirectional reflectance distribution function (BRDF) acquired from high frequency multi-angular spectral observations. Key wavebands for observing changes in directionally corrected canopy spectra were identified using discriminant analysis and highly significant correlations between canopy reflectance and field measured disturbance levels were found for several broad and narrow waveband vegetation indices (for instance, r²_NDVI = 0.90; r²_CHL3 = 0.85; p < 0.05). Results indicate that multi-angular observations are useful for extraction of disturbance related changes in canopy reflectance, in particular the temporally and spectrally dense data detected changes in chlorophyll content well. This study will help guide and inform future efforts to map forest health conditions at landscape and over increasingly coarse scales.     

NDVI-derived land cover classifications at a global scale

Phenological differences among vegetation types, reflected in temporal variations in the Normalized Difference Vegetation Index (NDVI) derived from satellite data, have been used to classify land cover at continental scales. Extending this technique to global scales raises several issues: identifying land cover types that are spectrally distinct and applicable at the global scale; accounting for phasing of seasons in different parts of the world; validating results in the absence of reliable information on global land cover; and acquiring high quality global data sets of satellite sensor data for input to land cover classifications. For this study, a coarse spatial resolution (one by one degree) data set of monthly NDVI values for 1987 was used to explore these methodological issues. A result of a supervised, maximum likelihood classification of eleven cover types is presented to illustrate the feasibility of using satellite sensor data to increase the accuracy of global land cover information, although the result has not been validated systematically. Satellite sensor data at finer spatial resolutions that include other bands in addition to NDVI, as well as methodologies to better identify and describe gradients between cover types, could increase the accuracy of results of global land cover data sets derived from satellite sensor data.     

Validating satellite phenology through intensive ground observation and landscape scaling in a mixed seasonal forest

Phenology is a key component of monitoring terrestrial ecosystem variations in response to global climate change. Satellite-measured land surface phenology (LSP) has been widely used to assess large scale phenological patterns and processes. However, the accuracy of LSP is rarely validated with spatially compatible field data due to the significant spatiotemporal scale mismatch. In this study, we employ intensive field observations specifically designed to address this deficiency. High density/frequency spring phenological observations were collected in a mixed seasonal forest during 2008 and 2009. A landscape up-scaling approach was used to derive landscape phenology (LP) indices from plot-level observations in order to validate Moderate-resolution Imaging Spectroradiometer (MODIS) based LSP. Results show that the MODIS Enhanced Vegetation Index (EVI) derived start of spring season (SOS) measure was able to predict LP full bud burst date with absolute errors less than two days. In addition, LSP derived SOS captured inter-annual variations and spatial differences that agreed with ground observations. Comparison of complete time series of LP and LSP revealed that fundamental differences exist between the two observation means, e.g., LP development had increased influence on LSP during the course of spring onset. Therefore, inferring continuous LP processes directly from LSP measures could be problematic. However, using LSP derived from techniques such as logistic curve modeling for extracting seasonal markers appears more promising. This study contributes to a more explicit understanding of the linkages between remotely sensed phenology and traditionally observed (ground-based) phenology.     

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.     

Land surface phenology dynamics and climate variations in the North East China Transect (NECT), 1982–2000

Land surface phenology dynamics reflect the response of the Earth's biosphere to inter‐ and intra‐annual dynamics of the Earth's climate and hydrologic regimes. Investigations of land surface phenology dynamics and its relation to long‐term climate variation could help us to detect the response of regional vegetation to climate variation. The present study developed a new algorithm for detecting regional land surface phenology dynamics (ARLSPD) and demonstrated it in detecting the vegetation response to inter‐annual climate variability in the North East China Transect (NECT), a mid‐latitude semi‐arid terrestrial transect with strong gradients in environmental conditions and vegetation formations. The spatial–temporal patterns of greenup‐onset date, maturity date, and senescence date during the period of 1982–2000 are presented. The resultant spatial–temporal patterns of land surface phenology were quite consistent with the land‐cover characteristics, moisture, and temperature gradients. The relations between inter‐annual variations in phenology and seasonal climate were investigated. It was found that besides human disturbance, land surface phenology depended primarily on the combined effects of preseason temperature and precipitation. The relative influence of preseason temperature and precipitation on land surface phenology was changing, which led to the different responses of land surface dynamics to climate variation along the moisture gradient in the NECT. In the arid and semi‐arid region of NECT, the dates of onset for phonological events in temperate typical grassland were most significantly related to the precipitation during the preceding 2–4 months. Temperature‐induced drought stress during the preceding months could delay greenup onset in cropland/grassland mosaic, and advance senescence in temporal typical grassland, and in cropland/grassland mosaic. The regional phenology algorithm, theoretically also applicable for complex ecosystems characterized by annual multiple growth cycles, is expected to couple with large‐scale biogeochemical models to regulate dynamically land surface phenology.     

基于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以上。     

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.     

Global interannual variability in terrestrial ecosystems: sources and spatial distribution using MODIS-derived vegetation indices,social and biophysical factors

Variations in global vegetation activity were measured at a global scale, from 2000 to 2006, based on the Enhanced Vegetation Index (EVI) extracted from the 1km resolution Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. Interannual variations in phenology and/or in annual integrated vegetation index are mapped using change metrics. The relationships between interannual variability and climate, ecosystem disturbances and land use are also examined. Around 14% of the study area experienced high interannual variability in land surface attributes over the six years. The ecosystems most subjected to large fluctuations in surface conditions were the boreal ecosystems, temperate ecozones, and subtropical and tropical steppes. These changes were largely related to rainfall variability and were associated with mean annual rainfall, agriculture, fire regimes and population density. Large population concentrations were mostly found in more stable ecozones. Rainfall and natural fire regimes explained more than half of the land surface variability in Australia. An additional global analysis on trends in vegetation activity also shows that around 4.5% of the vegetated surface of the Earth, excluding deserts and frequently cloudy regions, experienced a continuous decrease in vegetation activity over the six years. This represents more than twice the area experiencing a greening trend over this time period and concerns mostly tropical, subtropical and temperate forest ecozones. The total change in vegetation activity at a global scale and per year amounted to –2% of the annual integrated EVI aggregated across all ecosystems of the study area, on average for the years 2001–2006.     

AERONET-based surface reflectance validation network (ASRVN) data evaluation: Case study for railroad valley calibration site

The AERONET-based Surface Reflectance Validation Network (ASRVN) is an operational processing system developed for validation of satellite derived surface reflectance products at regional and global scales. The ASRVN receives 50 × 50 km² subsets of MODIS data centered at AERONET sites along with AERONET aerosol and water vapor data, and performs an atmospheric correction. The ASRVN produces surface bidirectional reflectance factor (BRF), albedo, parameters of the Ross-Thick Li-Sparse (RTLS) BRF model, as well as Hemispherical-Directional Reflectance Factor (HDRF), which is required for comparison with the ground-based measurements. This paper presents a comparison of ASRVN HDRF with the ground-based HDRF measurements collected during 2001-2008 over a bright calibration Railroad Valley, Nevada site as part of the MODIS land validation program. The ground measurements were conducted by the Remote Sensing Group (RSG) at the University of Arizona using an ASD spectrometer. The study reveals a good agreement between ASRVN and RSG HDRF for both MODIS Terra and Aqua with rmse ~ 0.01-0.025 in the 500 m MODIS land bands B1-B7. Obtained rmse is below uncertainties due to the spatial and seasonal variability of the bright calibration 1 km² area. While two MODIS instruments have a similar rmse in the visible bands, MODIS Aqua has a better agreement (lower rmse) with the ground data than MODIS Terra at wavelengths 0.87-2.1 μm. An independent overall good agreement of two MODIS instruments with the ground data indicates that the relative calibration of MODIS Terra and Aqua at medium-to-bright reflectance levels for the stated time period is significantly better than uncertainties of the ASRVN and ground data.