Spatial structure,sampling design and scale in remotely-sensed imagery of a California savanna woodland

Abstract

This article describes research related to sampling techniques for establishing linear relations between land surface parameters and remotely-sensed data. Predictive relations are estimated between percentage tree cover in a savanna environment and a normalized difference vegetation index (NDVI) derived from the Thematic Mapper sensor. Spatial autocorrelation in original measurements and regression residuals is examined using semi-variogram analysis at several spatial resolutions. Sampling schemes are then tested to examine the effects of autocorrelation on predictive linear models in cases of small sample sizes. Regression models between image and ground data are affected by the spatial resolution of analysis. Reducing the influence of spatial autocorrelation by enforcing minimum distances between samples may also improve empirical models which relate ground parameters to satellite data.     

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.     

Relationships between percent vegetation cover and vegetation indices

In this paper, percent vegetation cover is estimated from vegetation indices using simulated Advanced Very High Resolution Radiometer (AVHRR) data derived from in situ spectral reflectance data. Spectral reflectance measurements were conducted on grasslands in Mongolia and Japan. Vegetation indices such as the normalized difference, soil-adjusted, modified soil-adjusted and transformed soil-adjusted vegetation indices (NDVI, SAVI, MSAVI and TSAVI) were calculated from the spectral reflectance of various vegetation covers. Percent vegetation cover was estimated using pixel values of red, green and blue bands of digitized colour photographs. Relationships between various vegetation indices and percent vegetation cover were compared using a second-order polynomial regression. TSAVI and NDVI gave the best estimates of vegetation cover for a wide range of grass densities.     

Land cover assessment with MODIS imagery in southern African Miombo ecosystems

Global land use and land cover products in highly dynamic tropical ecosystems lack the detail needed for natural resource management and monitoring at the national and provincial level. The MODIS sensor provides improved opportunities to combine multispectral and multitemporal data for land use and land cover mapping. In this paper we compare the MODIS Global Land Cover Classification Product with recent land use and land cover maps at the national level over a characteristic location of Miombo woodlands in the province of Zambezia, Mozambique. The performances of three land cover-mapping approaches were assessed: single-date supervised classification, principal component analysis of band-pair difference images, and multitemporal NDVI analysis. Extensive recent field data were used for the definition of the test sites and accuracy assessment. Encouraging results were achieved with the three approaches. The classification results were refined with the help of a digital elevation model. The most consistent results were achieved using principal component analysis of band-pair difference images. This method provided the most accurate classifications for agriculture, wetlands, grasslands, thicket and open forest. The overall classification accuracy reached 90%. The multitemporal NDVI provided a more accurate classification for the dense forest cover class. The selection of the right image dates proved to be critical for all the cases evaluated. The flexibility of these alternatives makes them promising options for rapid and inexpensive land cover mapping in regions of high environmental variability such as tropical developing countries.     

Discriminating between cool season and warm season grassland cover types in northeastern Kansas

This study assesses the ability of multitemporal Landsat Thematic Mapper (TM) data and the normalized difference vegetation index (NDVI) to spectrally separate grazed cool season and warm season grassland cover types in Douglas County, Kansas. Biophysical data collected during the summer of 1997 suggest that differences in the per cent of total living vegetation cover, per cent of senescent vegetation, and proportion of forb cover between the two grassland cover types could make cool season and warm season grassland cover types spectrally distinct. The results show that the two grassland cover types were spectrally different in several spring (May) and mid-summer (July) bands, but not in any fall (September) bands. Furthermore, the two grassland cover types could be discriminated with a high level of accuracy. Accuracy assessments of the three single dates showed that the mid-summer (July) image and NDVI discriminated between the grassland cover types most accurately (81.8%). The multitemporal TM and NDVI data did not improve the spectral discrimination of the two grassland cover types over the mid-summer image or NDVI and had classification accuracy levels of 63.6% and 68.2%, respectively.     

Smoke effects on NDVI determination of savannah vegetation types

Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data collected on 25 and 27 August and 1 September 1995, in central Brazil, were used to describe smoke effects on Normalized Difference Vegetation Index (NDVI) determination of savannah vegetation types. Image-derived empirical relationships were used to predict the response of the red band (sensitive to smoke) from the shortwave infrared (SWIR) bands (insensitive to smoke), so as to minimize smoke effects on NDVI. Results showed that smoke effects were target dependent and produced NDVI differences up to 0.35 between the three dates for dense vegetation covers such as savannah woodland and riparian forest. Although this approach did not consider the smoke influence on the near-infrared (NIR) interval, the small SWIR-derived NDVI differences observed between the dates for a given land cover, under distinct smoke conditions, indicated its applicability in the Brazilian cerrado.     

Boosted decision tree classifications of land cover over Turkey integrating MODIS,climate and topographic data

This study investigates the impact of using different combinations of Moderate Resolution Imaging Spectroradiometer (MODIS) and ancillary datasets on overall and per-class classification accuracies for nine land cover types modified from the classification system of the International Geosphere Biosphere Programme (IGBP). Twelve land cover maps were generated for Turkey using boosted decision trees (BDTs) based on the stepwise addition of 14 explanatory variables derived from a time series of 16-day MODIS composites between 2000 and 2006 (Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and four spectral bands) and ancillary climate and topographic data (minimum and maximum air temperature, precipitation, potential evapotranspiration, aspect, elevation, distance to sea and slope) at 500-m resolution. Evaluation of the 12 BDTs indicated that the BDT built as a function of all the MODIS and climate variables, aspect and elevation produced the highest degree of overall classification accuracy (79.8%) and kappa statistic (0.76) followed by the BDTs that additionally included distance to sea (DtS), and both DtS and slope. Based on an independent validation dataset derived from a pre-existing national forest map and Landsat images of Turkey, the highest overall accuracy (64.7%) and kappa coefficient (0.58) among the 12 land cover maps was achieved by using MODIS-derived NDVI time series only, followed by NDVI and EVI time series combined; NDVI, EVI and four MODIS spectral bands; and the combination of all MODIS and climate data, aspect, elevation and distance to sea, respectively. The largest improvements in producer's accuracies were observed for grasslands (+50%), barrenlands (+46%) and mixed forests (+39%) and in user's accuracies for grasslands (+53%), shrublands (+30%) and mixed forests (+28%), in relation to the lowest producer's accuracy. The results of this study indicate that BDTs can increase the accuracy of land cover classifications at the national scale.     

Satellite remote sensing of rangelands in Botswana II. NOAA AVHRR and herbaceous vegetation

NOAA-7 Advanced Very High Resolution Radiometer (AVHRR) global-area coverage (GAC) data for the visible and near-infrared bands were used to investigate the relationship between the normalized difference vegetation index (NDVI) and the herbaceous vegetation in three representative rangeland types in eastern Botswana. Regressions between Landsat MSS band-7/band-5 ratios and field measurements of the cover of the live parts of herbaceous plants, above-ground biomass of live herbaceous plants and bare ground were used in conjunction with MSS data in order to interpolate the field data to 144 km² areas for comparison with blocks of nine AVHRR GAC pixels. NOAA NDVI data were formed into 10-day composites in order to remove cloud cover and extreme off-nadir viewing angles. Both individual NDVI composite data and multitemporal integrations throughout the period May 1983-April 1984 were compared with the field data.

In multiple linear regressions, the cover and biomass of live herbaceous plants and bare ground measurements accounted for 42, 56 and 19 per cent respectively of the variation in NDVI. When factors were included in I he regression models to specify the site and date of acquisition of the data, between 93 and 99 per cent of the variation in NDVI was accounted for. The total herbaceous biomass at the end of the season was positively related to integrated NDVI, up lo the maximum biomass observed in a 12km × 12km area (590kgha^?1)- These results give a different regression of herbaceous biomass values on integrated AVHRR NDVI to that reported by Tucker et at. (1985 b) for Senegalese grasslands. The effect of the higher cover of the tree canopy in Botswana on this relationship and on the detection of forage available to livestock is discussed.     

Multi-year monitoring of rice crop phenology through time series analysis of MODIS images

Precise phenological calendars, for each cultivated species and variety, are necessary both to highlight anomalous agronomic situations and to feed crop models. This study, conducted in the Italian rice area, focuses on the evaluation of the contribution of remote sensing satellite data to providing phenological information on rice cropping systems. A time series of 5 years (2001–2005) of the Moderate Resolution Imaging Spectrometer (MODIS) Normalized Difference Vegetation Index (NDVI) 16-day composites was analysed with the TIMESAT program in order to automatically retrieve key phenological information such as the start of season (emergence), peak (heading) and end of season (maturity). The procedure involved two steps: (1) interpolation and smoothing of MODIS NDVI temporal profile and (2) the analysis of a temporal signal for the extraction of the phenological metrics. The remote sensing estimates were evaluated using information regarding cultivated variety, sowing dates, management and production directly acquired from rice farmers. A good correlation (R ²?=?0.92, n?=?24) has been observed between estimates derived from satellites and estimates produced with the traditional Growing Degree Days (GDD) method based on thermal unit accumulation. Improved estimates of the maturity stage were obtained using a procedure that integrates satellite and GDD methods; however its application requires spatially distributed information on the cultivated varieties. Satellite derived maps of the retrieved phenological parameters showed an intra-seasonal pattern related to different cultivated varieties. Inter-seasonal analysis allowed the anomalous behaviour of the year 2003 to be highlighted, characterized by rapid growth at the beginning of the spring and an early senescence. The results confirm the potential of remotely sensed data for the monitoring of crop status and for the forcing of crop models in a spatially distributed way.     

Analysis of NDVI and scaled difference vegetation index retrievals of vegetation fraction

The normalized difference vegetation index (NDVI) is the most widely used vegetation index for retrieval of vegetation canopy biophysical properties. Several studies have investigated the spatial scale dependencies of NDVI and the relationship between NDVI and fractional vegetation cover, but without any consensus on the two issues. The objectives of this paper are to analyze the spatial scale dependencies of NDVI and to analyze the relationship between NDVI and fractional vegetation cover at different resolutions based on linear spectral mixing models. Our results show strong spatial scale dependencies of NDVI over heterogeneous surfaces, indicating that NDVI values at different resolutions may not be comparable. The nonlinearity of NDVI over partially vegetated surfaces becomes prominent with darker soil backgrounds and with presence of shadow. Thus, the NDVI may not be suitable to infer vegetation fraction because of its nonlinearity and scale effects. We found that the scaled difference vegetation index (SDVI), a scale-invariant index based on linear spectral mixing of red and near-infrared reflectances, is a more suitable and robust approach for retrieval of vegetation fraction with remote sensing data, particularly over heterogeneous surfaces. The proposed method was validated with experimental field data, but further validation at the satellite level would be needed.     

Relations between AVHRR NDVI and ecoclimatic parameters in China

Based on monthly AVHRR NDVI data and weather records collected from 160 stations throughout China for 10 years, the relationships between NDVI and two ecoclimatic parameters (growing degree-days (GDD) and rainfall) were analysed. The results indicate that a significant correlation exists between NDVI and the two ecoclimatic parameters; the NDVI/GDD correlation was stronger than the NDVI/rainfall correlation. The NDVI/rainfall correlation coefficient exhibits a clear structure in terms of spatial distribution. Further, the results indicate that for natural vegetation, the NDVI/rainfall correlation coefficient increases in order from evergreen forest, to deciduous forest, to shrubs and desert vegetation, to steppe and savanna. The correlation coefficients associated with cultural vegetation type depend on a number of factors including annual rainfall, seasonal variation in precipitation, type and intensity of irrigation practice and other environmental factors.     

An assessment of AVHRR/NDVI-ecoclimatological relations in Nebraska,U.S.A.

Research designed to better define relations between 1-km multitemporal AVHRR-derived NDVI data and selected climatological parameters, soil hydrological properties and land cover characteristics is summarized. Bi-weekly maximum value composite NDVI data and concurrently measured meteorological data acquired in 1990 and 1991 for Nebraska were utilized to study relations between NDVI and accumulated growing degree days,soil temperature, precipitation and potential evapotranspiration. Temporal change in NDVI was found to be closely linked with the temperature regime. NDVI-precipitation and NDVI-potential evapotranspiration relations exhibited time lags, although the length of lag varied with land cover type, precipitation, and soil hydrologic properties. NDVI response to precipitation was stronger in natural grasslands and grassland/wet meadows than in areas of irrigated cropland and mixed crop/ grass. NDVI-climate relations were strongest where vegetation was developed on soils with low root zone available water capacity and high permeability. Relations derived by using NDVI values over 3 pixel by 3 pixel windows showed little difference from those using single 1 km pixel. This may reflect both the relatively homogeneous land cover characteristics of the study area and the effect of off-nadir viewing geometry on AVHRR data acquisition.     

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.     

A new approach for using time-series remote-sensing images to detect changes in vegetation cover and composition in drylands: a case study of eastern Kenya

Vegetation cover and composition are important aspects of the dryland environment because they provide livelihood to humans and also protect soil resources against erosion. Currently, scientists are advancing various techniques for detecting vegetation degradation in the drylands and the possibilities for its control. This study contributed through the testing of time-series mixed-effects modelling of the normalized difference vegetation index (NDVI) and rainfall relationship to trace the footprints of vegetation dynamics in the drylands. The approach aimed at providing guidelines for quick diagnosis of the changes in vegetation cover and composition to trigger necessary action. The mixed-effects technique used in this study is a novel regression approach for simultaneous modelling of the NDVI–rainfall relationship in different dominant vegetation types. Its time-series application with Advanced Very High Resolution Radiometer (AVHRR) and Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI images between 1982 and 2008 was tested in eastern Kenya. The results show how the original dominant vegetation types had been converted to cereal croplands, open grasslands, or reduced to bare ground in a span of 27 years. In some places, it shows how the changes in vegetation composition resulted in the overall loss of vegetation cover. Field validation positively confirmed these observations; thus, indicating that the method was a promising tool for tracing vegetation dynamics in the drylands. In spite of its success, the method was found to be only useful in detecting changes in large areas with dominant vegetation types. The technique can therefore be recommended for regional analysis, and can be used as a first approximation to guide more detailed subsequent analysis.     

NDVI from low altitude aircraft and composited NOAA AVHRR data for scaling Arctic ecosystem fluxes

The relationship between normalized difference vegetation index (NDVI) patterns obtained from high spatial resolution aircraft and low spatial resolution satellite data (Advanced Very High Resolution Radiometer (AVHRR)) was investigated with the intent of using multilevel data to scale carbon flux models in Arctic tundra ecosystems. Despite variable illumination conditions during the aircraft missions and maximum value compositing of the AVHRR data, the difference between 3?km average aircraft and AVHRR NDVI values was generally constant along each flight transect. However, the magnitude of the offset differed between flight dates and small lakes had a greater effect on area averaged aircraft NDVI values than on the satellite values. A cloud index was calculated using incident solar radiation measured by the aircraft and this index was used to identify periods when the aircraft NDVI values may have been biased by cloud cover. Removal of NDVI values based on a cloud index threshold did not appear to be justified given the marginal improvement in the relationship between the two NDVI datasets. If the systematic difference between AVHRR and aircraft NDVI values can be determined, then the scaling of carbon flux models based on the NDVI should be a viable approach in Arctic ecosystems.     

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

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

Classification management for grassland using MODIS data: a case study in the Gannan region,China

Classification of grasslands is a convenient method to measure and manage the sustainability of Chinese grasslands. In this study, a timely and reliable procedure was examined using remote-sensing (RS) techniques. Linear regression analysis between field survey data and Moderate-Resolution Imaging Spectroradiometer (MODIS) data showed that among 17 vegetation indices (VIs) evaluated, the enhanced vegetation index (EVI) was the best VI to simulate forage dry biomass and cover in the Gannan region. The results of precision estimation of the models showed that power and logarithm regression satisfactorily simulated grassland dry biomass and grassland cover, respectively. The index of classification management of grasslands (ICGs) was used to subdivide grasslands into conservation grasslands and moderately productive grasslands in the Gannan region, where no grasslands fell into intensively productive grasslands. Conservation grasslands accounted for 2.04% of the available grasslands, whereas moderately productive grasslands were 97.96% of the available grasslands, and this is related to the history of the grasslands’ use and the per capita income in the Gannan region. This study proposes that the area of conservation grasslands and that of moderately productive grasslands are determined by increases in per capita income and changes in the human use of grasslands.     

Land surface phenology, climatic variation, and institutional change: Analyzing agricultural land cover change in Kazakhstan

Kazakhstan is the second largest country to emerge from the collapse of the Soviet Union. Consequent to the abrupt institutional changes surrounding the disintegration of the Soviet Union in the early 1990s, Kazakhstan has reportedly undergone extensive land cover/land use change. Were the institutional changes sufficiently great to affect land surface phenology at spatial resolutions and extents relevant to mesoscale meteorological models? To explore this question, we used the NDVI time series (1985-1988 and 1995-1999) from the Pathfinder Advanced Very High Resolution Radiometer (AVHRR) Land (PAL) dataset, which consists of 10 days maximum NDVI composites at a spatial resolution of 8 km. Daily minimum and maximum temperatures were extracted from the NCEP Reanalysis Project and 10 days composites of accumulated growing degree-days (AGDD) were produced. We selected for intensive study seven agricultural areas ranging from regions with rain-fed spring wheat cultivation in the north to regions of irrigated cotton and rice in the south. We applied three distinct but complementary statistical analyses: (1) nonparametric testing of sample distributions; (2) simple time series analysis to evaluate trends and seasonality; and (3) simple regression models describing NDVI as a quadratic function of AGDD.The irrigated areas displayed different temporal developments of NDVI between 1985-1988 and 1995-1999. As the temperature regime between the two periods was not significantly different, we conclude that observed differences in the temporal development of NDVI resulted from changes in agricultural practices.In the north, the temperature regime was also comparable for both periods. Based on extant socioeconomic studies and our model analyses, we conclude that the changes in the observed land surface phenology in the northern regions are caused by large increases in fallow land dominated by weedy species and by grasslands under reduced grazing pressure. Using multiple lines of evidence allowed us to build a case of whether differences in land surface phenology were mostly the result of anthropogenic influences or interannual climatic fluctuations.     

数据融合技术在提高NPP估算精度中的应用

针对现有遥感数据不能同时满足在时间和空间上精确监测植被动态变化的问题,提出利用时空适应性反射率融合模型（STARFM）的方法对MODIS-NDVI和TM-NDVI影像数据进行融合处理获得30 m较高时空分辨率的融合NDVI影像,进而将多种尺度的MODIS-NDVI和融合NDVI数据分别输入到CASA模型,对锡林浩特地区进行植被净初级生产力（NPP）的多尺度估算。将不同尺度的NPP估算结果与地上生物量地面实测值进行验证比较,结果表明：随着输入NDVI空间分辨率的提高,NPP估算值与实测地上生物量之间的相关性也逐渐增大,[r]最大值达到了0.915。此外以融合NDVI影像作为输入数据之一的NPP估算值与实测地上生物量的相关性均比未融合NDVI的相关性高,说明融合NDVI估算NPP的效果较未融合NDVI好,并且以融合NDVI影像作为模型输入数据可提高NPP估算精度。