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.     

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.     

Analysis of climate change impacts on lake ice phenology in Canada using the historical satellite data record

Variability and trends in lake ice dynamics (i.e. lake ice phenology) are related to climate conditions. Climate influences the timing of lake ice melt and freeze onset, ice duration, and lake thermal dynamics that feedback to the climate system initiating further change. Phenology records acquired in a consistent manner and over long time periods are required to better understand variability and change in climate conditions and how changes impact lake processes. In this study, we present a new technique for extracting lake ice phenology events from historical satellite records acquired by the series of Advanced Very High Resolution Radiometer (AVHRR) sensors. The technique was used to extend existing in-situ measurements for 36 Canadian lakes and to develop records for 6 lakes in Canada's far north. Comparison of phenology events obtained from the AVHRR record and in-situ measurements show strong agreement (20 lakes, 180 cases) suggesting, with high confidence especially in the case of break-up dates, the use of these data as a complement to ground observations. Trend analysis performed using the combined in-situ and AVHRR record ∼ 1950-2004 shows earlier break-up (average — 0.18 days/year) and later freeze-up (average 0.12 days/year) for the majority of lakes analyzed. Less confidence is given to freeze-up date results due to lower sun elevation during this period making extraction more difficult. Trends for the 20 year record in the far north showed earlier break-up (average 0.99 days/year) and later freeze-up (average 0.76 days/year). The established lake ice phenology database from the historical AVHRR image archive for the period from 1985 to 2004 will to a certain degree fill data gaps in the Canadian in-situ observation network. Furthermore, the presented extraction procedure is not sensor specific and will enable continual data update using all available satellite data provided from sensors such as NOAA/AVHRR, MetOp/AVHRR, MODIS, MERIS and SPOT/VGT.     

A model for the estimation of the daily global sunshine duration from meteorological geostationary satellite data

In this research a simple model was proposed for the estimation of daily global sunshine duration and constructing spatially continous sunshine duration map from meteorological geostationary satellite sensor data. First cloud cover index and then daily mean cloud cover index value for each pixel were computed using a time‐series of Meteosat C3D visible type images. The statistical relationship between daily mean cloud cover index and measured bright sunshine duration values was tested and found to be linear. Using regression parameters daily sunshine duration of each pixel was calculated and then daily sunshine duration and monthly mean daily sunshine duration maps were constructed. It was shown that, by using the suggested model it is possible to calculate daily sunshine duration values over a large area where the sunshine duration data are not available and construct spatially continious long‐term sunshine duration maps.     

气象卫星资料虚拟存储系统设计与实现

利用服务器连接（SAS）技术和存储局域网（SAN）技术构建气象卫星资料存储系统。实现群集环境中计算机系统间的存储资源共享，设计虚拟存储管理软件对不同的物理存储设备进行逻辑管理，可以简化存储设备的使用，可以提高系统运行效率，国家卫星气象中心的气象卫星资料存储系统的运行结果表明，虚拟存储是解决群体环境中不同物理存储设备有效使用的一条切实可行的技术路线。     

Calculation of the sensible heat flux of the global ocean using satellite data

The sensible heat flux of the global ocean is derived using satellite data of Special Sensor Microwave/Imager (SSM/I) precipitable water, Advanced Very High Resolution Radiometer (AVHRR) sea-surface temperature (SST) and scatterometer wind speed. Prior to heat flux derivation, the air temperature over the sea surface in the global ocean is obtained with an iterative solving technique applied to a simplified equation that specifies the relationship among boundary-layer parameters. It is found that a bias exists between the calculated air temperature and the climatology data, which is corrected by a linear model based on the climatology of air temperature. Using the corrected air temperature and the bulk formula, we calculate the sensible heat flux from January 1992 to October 1998. The heat flux calculation is consistent with previous results. An error analysis suggests that based on the bulk formula, the uncertainty caused by the air temperature is comparable with the error components from the SST and wind. Empirical orthogonal function (EOF) analysis is used to extract the temporal and spatial characteristics of the calculated sensible heat flux. The first EOF is characterized by a winter and summer oscillation with an annual cycle, the second shows another annual cycle in spring and autumn oscillation, the third EOF is related to the El Niño and La Niña cycle, reflecting a certain relationship between El Niño Southern Oscillation (ENSO) events and the global ocean sensible heat flux, and the fourth EOF indicates an increasing trend with a quasi-biennial oscillation and atmospheric influences.     

Daily global solar radiation mapping of Turkey using Meteosat satellite data

Many studies have indicated that the estimation of solar irradiation at ground level using meteorological satellite data has been an alternative and easy method compared to classical methods. In the present work, the incident of solar radiation over Turkey has been estimated at ground level between July 1997 and December 1998. Statistical regressions between ground data and digital satellite data, measured in the visible band (0.4–1.1?µm) by Meteosat radiometer, have been determined and these regression parameters have been used to estimate solar radiation at ground level. This is the so-called statistical method, which uses a simple model because satellites measure only a few parameters among the many that govern radiative transfers.

The visible image (C3D) data used in the present work was Meteosat Wefax type. While pursuing our studies the mean daily sum of global solar radiation over Turkey has been determined to be 18.44?MJ?m^?2?d^?1 with a correlation coefficient of 0.96. The rms error for the mean daily sum has been evaluated as 1.92?MJ?m^?2?d^?1. The monthly mean daily sum of solar radiation has been determined with an rms error of 1.82?MJ?m^?2?d^?1 in two years. During this period the maximum value of the daily sum has been found to occur in June 1998 as 28.47?MJ?m^?2?d^?1, whereas the minimum has been found to occur in December 1998 as 7.35?MJ?m^?2?d^?1. The evaluation procedure, results and possible sources of error are suggested and possible ways of improving the method are described and discussed.     

Hurricane Georges and vegetation change in Puerto Rico using AVHRR satellite data

This study assesses natural disturbances at Puerto Rico resulting from hurricane Georges in September 1998. The study was done using data from the National Oceanic and Atmospheric Administration's Advanced Very High Resolution Radiometer (AVHRR) 14 satellite. Specifically, the Luquillo Experimental Forest, the Guanica Dry Forest and five cities were examined using AVHRR data. It was found that although there is probably no significant relationship between Georges and the before/after temperature data, there is a statistically significant relationship between distance to the location affected by the hurricane and the extent of changes in NDVI, a fact that suggests that it is possible to measure hurricane impacts on vegetation by using AVHRR data.     

极轨气象卫星有效载荷航空校飞数据地理定位

航空校飞数据的地理定位方法用来确定数据中每一个像素的地理经纬度,它是整个数据处理过程的基础。介绍了一种用于航空校飞试验数据的地理定位方法。首先给出了问题的描述;其次建立了航迹和姿态模型;然后给出了地理定位方法的算法;最后给出计算结果及误差分析,验证了方法的可行性以及高精度。地理定位方法不仅仅可以用于极轨气象卫星的校飞试验数据处理,同时也适用于其他遥感卫星的校飞试验数据处理,具有非常重要的意义。     

Monitoring post-burn recovery of chaparral vegetation in southern California using multi-temporal satellite data

Monitoring the post-burn recovery condition of chaparral vegetation in southern California is important for managers to determine the appropriate time to conduct controlled burns. Due to the difficulty of monitoring post-fire recovery over large areas and the absence of detailed fire records in many areas, we examined the possibility of using satellite observations to establish the postfire recovery stage of chamise chaparral stands in this region. SPOT XS data collected on three dates between 1986 and 1992 were analysed to determine if temporal changes in a spectral vegetation index tracked the expected post-fire recovery trajectory of the above-ground biomass of chamise chaparral stands of varying post-fire ages. Results of the study indicated that neither the normalized difference vegetation index nor the soil adjusted vegetation index followed the expected post-fire recovery patterns in these stands. These findings are explained by interannual variations in precipitation having a larger than expected effect on the growth of this drought-resistant evergreen community, with changes in green leaf area dominating the temporal variations in the spectral vegetation indices.     

Retrieval of leaf area index in different vegetation types using high resolution satellite data

With the successful launch of the IKONOS satellite, very high geometric resolution imagery is within reach of civilian users. In the 1-m spatial resolution images acquired by the IKONOS satellite, details of buildings, individual trees, and vegetation structural variations are detectable. The visibility of such details opens up many new applications, which require the use of geometrical information contained in the images. This paper presents an application in which spectral and textural information is used for mapping the leaf area index (LAI) of different vegetation types. This study includes the estimation of LAI by different spectral vegetation indices (SVIs) combined with image textural information and geostatistical parameters derived from high resolution satellite data. It is shown that the relationships between spectral vegetation indices and biophysical parameters should be developed separately for each vegetation type, and that the combination of the texture indices and vegetation indices results in an improved fit of the regression equation for most vegetation types when compared with one derived from SVIs alone. High within-field spatial variability was found in LAI, suggesting that high resolution mapping of LAI may be relevant to the introduction of precision farming techniques in the agricultural management strategies of the investigated area.     

A simple single layer model to estimate transpiration from vegetation using multi-spectral and meteorological data

A methodology is developed here to model evapotranspiration (λE_c ) from the canopy layer over large areas by combining satellite and ground measurements of biophysical and meteorological variables. The model developed here follows the energy balance approach, where λE_c is estimated as a residual when the net radiation (Rn), sensible heat flux (H) and ground flux (G) are known. Multi-spectral measurements from the NOAA Advanced Very High Resolution Radiometer (AVHRR) were used along with routine meteorological measurements made on the ground to estimate components of the energy balance. The upwelling long wave radiation, and H from the canopy layer were modelled using the canopy temperature, obtained from a linear relation between the Normalized Difference Vegetation Index (NDVI) and surface temperature. This method separates flux measurements from the canopy and bare soil without the need for a complex two layer model. From theoretical analysis of canopy reflectance, leaf area, and canopy resistance, a model is developed to scale the transpiration estimates from the full canopy to give an area averaged estimate from the mean NDVI of the study area. The model was tested using data collected from the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE), and the results show that the modelled values of total surface evapotranspiration from the soil and canopy layers vary from the ground measurements by less than 9 per cent.     

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.     

Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data

A vegetation index (VI) model for predicting evapotranspiration (ET) from data from the Moderate Resolution Imaging Spectrometer (MODIS) on the EOS-1 Terra satellite and ground meteorological data was developed for riparian vegetation along the Middle Rio Grande River in New Mexico. Ground ET measurements obtained from eddy covariance towers at four riparian sites were correlated with MODIS VIs, MODIS land surface temperatures (LSTs), and ground micrometeorological data over four years. Sites included two saltcedar (Tamarix ramosissima) and two Rio Grande cottonwood (Populus deltoides ssp. Wislizennii) dominated stands. The Enhanced Vegetation Index (EVI) was more closely correlated (r=0.76) with ET than the Normalized Difference Vegetation Index (NDVI; r=0.68) for ET data combined over sites and species. Air temperature (T_a) measured over the canopy from towers was the meteorological variable that was most closely correlated with ET (r=0.82). MODIS LST data at 1- and 5-km resolutions were too coarse to accurately measure the radiant surface temperature within the narrow riparian corridor; hence, energy balance methods for estimating ET using MODIS LSTs were not successful. On the other hand, a multivariate regression equation for predicting ET from EVI and T_a had an r²=0.82 across sites, species, and years. The equation was similar to VI-ET models developed for crop species. The finding that ET predictions did not require species-specific equations is significant, inasmuch as these are mixed vegetation zones that cannot be easily mapped at the species level.     

Objective assessment of the NOAA global vegetation index data product

Abstract

For the last 10 years the U.S. National Oceanic and Atmospheric Administration has produced an experimental Global Vegetation Index (GVI) data set for terrestrial vegetation research. These data, sampled from advanced very high resolution radiometer (AVHRR) observations, have served as a primary stimulus for global-scale vegetation research but have, so far, not been adequately evaluated. This study reviews the GVI production procedures and compares the resultant observations with a more comprehensive compilation of the AVHRR data being produced at the NASA Goddard Space Flight Center. There are many aspects of the GVI production procedures which could be improved to achieve the desired objectives. In particular, the mapping and sampling procedures employed provide measurements which only approximate the observed GAC measurements. The GVI NDVI record varies more than ±NDVI units (～ 7 per cent of signal) from the GAC record and, in general, seriously underestimates the GAC NDVI measurements. The NDVI portion of the GVI record is compromised through use of digital numbers rather than calibrated reflectance. NDVI measurements from the calibrated channels of the GVI data set produces values that compare favourably with the GAC measurements, but with considerable residual variance. Calculation of a 3 by 3 pixel average of the GVI NDVI measurements reduces residual variance between the data sets to ±0.018 NDVI units (～3 per cent of signal). Decay of sensor calibration and orbital overpass time, experienced by all the AVHRR sensors, as well as differences in these parameters between the sensors are not addressed but the results suggest the importance of accounting for these factors. These results indicate that GVI data sets, following adequate reprocessing, provide reasonable estimates of major regional contrasts in vegetation activity but should not be employed to evaluate local or minor trends.     

North American vegetation dynamics observed with multi-resolution satellite data

We investigated normalized difference vegetation index data from the NOAA series of Advanced Very High Resolution Radiometers and found regions in North America that experienced marked increases in annual photosynthetic capacity at various times from 1982 to 2005. Inspection of these anomalous areas with multi-resolution data from Landsat, Ikonos, aerial photography, and ancillary data revealed a range of causes for the NDVI increases: climatic influences; severe drought and subsequent recovery; irrigated agriculture expansion; insect outbreaks followed by logging and subsequent regeneration; and forest fires with subsequent regeneration. Vegetation in areas in the high Northern Latitudes appear to be solely impacted by climatic influences. In other areas examined, the impact of anthropogenic effects is more direct. The pattern of NDVI anomalies over longer time periods appear to be driven by long-term climate change but most appear to be associated with climate variability on decadal and shorter time scales along with direct anthropogenic land cover conversions. The local variability of drivers of change demonstrates the difficulty in interpreting changes in NDVI and indicates the complex nature of changes in the carbon cycle within North America. Coarse scale analysis of changes could well fail to identify the important local scale drivers controlling the carbon cycle and to identify the relative roles of disturbance and climate change. Our results document regional land cover land use change and climatic influences that have altered continental scale vegetation dynamics in North America.     

A global poverty map derived from satellite data

A global poverty map has been produced at 30 arcsec resolution using a poverty index calculated by dividing population count (LandScan 2004) by the brightness of satellite observed lighting (DMSP nighttime lights). Inputs to the LandScan product include satellite-derived land cover and topography, plus human settlement outlines derived from high-resolution imagery. The poverty estimates have been calibrated using national level poverty data from the World Development Indicators (WDI) 2006 edition. The total estimate of the numbers of individuals living in poverty is 2.2 billion, slightly under the WDI estimate of 2.6 billion. We have demonstrated a new class of poverty map that should improve over time through the inclusion of new reference data for calibration of poverty estimates and as improvements are made in the satellite observation of human activities related to economic activity and technology access.     

A blended global snow product using visible,passive microwave and scatterometer satellite data

A joint US Air Force/National Aeronautics and Space Administration (NASA) blended global snow product that uses Earth Observation System Moderate Resolution Imaging Spectroradiometer (MODIS), Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) and Quick Scatterometer (QuikSCAT or QSCAT) data has been developed. Existing snow products derived from these sensors have been blended into a single, global, daily, user-friendly product by using a newly developed Air Force Weather Agency (AFWA)/NASA Snow Algorithm (ANSA). This initial blended snow product uses minimal modelling to expeditiously yield improved snow products, which include, or will include, snow-cover extent, fractional snow cover, snow water equivalent (SWE), onset of snowmelt and identification of actively melting snow cover. The blended snow products are currently 25-km resolution. These products are validated with data from the lower Great Lakes region of the USA, from Colorado obtained during the Cold Land Processes Experiment (CLPX), and from Finland. The AMSR-E product is especially useful in detecting snow through clouds; however, passive microwave data miss snow in those regions where the snow cover is thin, along the margins of the continental snowline, and on the lee side of the Rocky Mountains, for instance. In these regions, the MODIS product can map shallow snow cover under cloud-free conditions. The confidence for mapping snow-cover extent is greater with the MODIS product than with the microwave product when cloud-free MODIS observations are available. Therefore, the MODIS product is used as the default for detecting snow cover. The passive microwave product is used as the default only in those areas where MODIS data are not applicable due to the presence of clouds and darkness. The AMSR-E snow product is used in association with the difference between ascending and descending satellite passes or diurnal-amplitude variations (DAV) to detect the onset of melt, and a QSCAT product will be used to map areas of snow that are actively melting.     

A global NDVI and EVI reference data set for land-surface phenology using 13 years of daily SPOT-VEGETATION observations

Time series of vegetation indices (VIs) obtained by remote sensing are widely used to study phenology on regional and global scales. The aim of the study is to design a method and to produce a reference data set describing the seasonal and inter-annual variability of the land-surface phenology on a global scale. Specific constraints are inherent in the design of such a global reference data set: (1) the high diversity of vegetation types and the heterogeneous conditions of observation, (2) a near-daily resolution is needed to follow the rapid changes in phenology, (3) the time series used to depict the baseline vegetation cycle must be long enough to be representative of the current vegetation dynamic and encompass anomalies, and (4) a spatial resolution consistent with a land-cover-specific analysis should be privileged. This study focuses on the SPOT (Satellite Pour l’Observation de la Terre)-VEGETATION sensor and its 13-year time series of reflectance values. Five steps addressing the noise and the missing data in the reflectance time series were selected to process the daily multispectral reflectance observations. The final product provides, for every pixel, three profiles for 52 × 7-day periods: a mean, a median, and a standard deviation profile. The mean and median profiles represent the reference seasonal pattern for variation of the vegetation at a specific location whereas the standard deviation profile expresses the inter-annual variability of VIs. A quality flag at the pixel level demonstrated that the reference data set can be considered as a reliable representation of the vegetation phenology in most parts of the Earth.     

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.