GIS based mapping of land cover changes utilizing multi-temporal remotely sensed image data in Lake Hawassa Watershed,Ethiopia

Classifying multi-temporal image data to produce thematic maps and quantify land cover changes is one of the most common applications of remote sensing. Mapping land cover changes at the regional level is essential for a wide range of applications including land use planning, decision making, land cover database generation, and as a source of information for sustainable management of natural resources. Land cover changes in Lake Hawassa Watershed, Southern Ethiopia, were investigated using Landsat MSS image data of 1973, and Landsat TM images of 1985, 1995, and 2011, covering a period of nearly four decades. Each image was partitioned in a GIS environment, and classified using an unsupervised algorithm followed by a supervised classification method. A hybrid approach was employed in order to reduce spectral confusion due to high variability of land cover. Classification of satellite image data was performed integrating field data, aerial photographs, topographical maps, medium resolution satellite image (SPOT 20 m), and visual image interpretation. The image data were classified into nine land cover types: water, built-up, cropland, woody vegetation, forest, grassland, swamp, bare land, and scrub. The overall accuracy of the LULC maps ranged from 82.5 to 85.0 %. The achieved accuracies were reasonable, and the observed classification errors were attributable to coarse spatial resolution and pixels containing a mixture of cover types. Land cover change statistics were extracted and tabulated using the ERDAS Imagine software. The results indicated an increase in built-up area, cropland, and bare land areas, and a reduction in the six other land cover classes. Predominant land cover is cropland changing from 43.6 % in 1973 to 56.4 % in 2011. A significant portion of land cover was converted into cropland. Woody vegetation and forest cover which occupied 21.0 and 10.3 % in 1973, respectively, diminished to 13.6 and 5.6 % in 2011. The change in water body was very peculiar in that the area of Lake Hawassa increased from 91.9 km² in 1973 to 95.2 km² in 2011, while that of Lake Cheleleka whose area was 11.3 km² in 1973 totally vanished in 2011 and transformed into mud-flat and grass dominated swamp. The “change and no change” analysis revealed that more than one third (548.0 km²) of the total area was exposed to change between 1973 and 2011. This study was useful in identifying the major land cover changes, and the analysis pursued provided a valuable insight into the ongoing changes in the area under investigation.     

Mapping the spatio-temporal distribution of key vegetation cover properties in lowland river reaches,using digital photography

The presence of vegetation in stream ecosystems is highly dynamic in both space and time. A digital photography technique is developed to map aquatic vegetation cover at species level, which has a very high spatial and a flexible temporal resolution. A digital single-lens reflex (DSLR) camera mounted on a handheld telescopic pole is used. The low-altitude (5 m) orthogonal aerial images have a low spectral resolution (red-green-blue), high spatial resolution (～1.9 pixels cm^?2, ～1.3 cm length) and flexible temporal resolution (monthly). The method is successfully applied in two lowland rivers to quantify four key properties of vegetated rivers: vegetation cover, patch size distribution, biomass and hydraulic resistance. The main advantages are that the method is (i) suitable for continuous and discontinuous vegetation covers, (ii) of very high spatial and flexible temporal resolution, (iii) relatively fast compared to conventional ground survey methods, (iv) non-destructive and (v) relatively cheap and easy to use, and (vi) the software is widely available and similar open source alternatives exist. The study area should be less than 10 m wide, and the prevailing light conditions and water turbidity levels should be sufficient to look into the water. Further improvements of the image processing are expected in the automatic delineation and classification of the vegetation patches.     

Aerial and terrestrial-based monitoring of channel erosion,headcutting, and sinuosity

Headcuts are points of accelerated channel erosion that frequently have ecological consequences. A particularly large and dynamic headcut in southwest Wyoming has affected natural and anthropogenic resources for decades. To better understand and address this issue, we undertook a review of the headcut’s upstream retreat, followed by photogrammetric monitoring of the present condition for erosion monitoring. Aerial photography shows the Bitter Creek headcut retreated >?200 m upstream in 68 years (1948–2016) at ~?1.4 m year^?1. Following installation of a concrete slab structure in the mid-1970s, headcut retreat slowed to ~?0.5 m year^?1. Channel sinuosity downstream of the headcut is greater than upstream, which we attribute to the presence of the headcut, given that there are no major changes in valley geometry, geology, or soils through this reach. Both aerial and terrestrial-based image platforms were used to collect stereo imagery and create 3D photogrammetric models of the headcut in 2016. From these two models, we measured soil loss downstream of the headcut at ~?126 m³ m^?1 valley length. Since 1954, soil loss within the channel has been ~?98 m³ year^?1 or ~871 t ha^?1 year^?1since then. Models created from aerial- and terrestrial-based images differed in volumetric estimates by 2%, indicating that either method could be used for this type of monitoring. The ground-based imagery model showed more detail, especially on vertical and overhanging surfaces, while the aerial imagery model produced a more realistic orthomosaic and efficiently covered a larger area. Ground-based image acquisition took longer and was more costly per unit area, but is an efficient method for small project areas, or areas where aerial imagery cannot be safely or practically acquired. Historical imagery and photogrammetric modeling proved very useful in elucidating stream dynamics associated with this large, dynamic headcut.     

Sampling for landscape elements—a case study from Lower Saxony,Germany

The estimation of coverage, i.e., the proportion of the total area in a study region covered by a given target class, is essential to many aspects of environmental monitoring. We analyze and compare the efficiency of different sample-based approaches for the estimation of coverage of different land cover classes from aerial imagery in a case study in Lower Saxony, Germany on the basis of the estimated standard errors. A complete delineation of vegetation classes in n?=?279 aerial photo plots of 400 × 400 m thrown onto the study region of 1,117.7 km² in accordance with a systematic grid is compared to different configurations of line intercept sampling and clusters of points. The observation designs under study are characterized by different complexity and total size of the observation units and therefore also to the efforts related to yield a single observation. Especially for those classes that cover a relatively large proportion of the sampling frame, our results show that difference in performance between the different designs are negligible. A cluster of four transects of 200 m each allows estimating the area of land cover classes with high coverage with nearly similar precision as a complete mapping of fixed area plots of 16 ha each. Clusters of points show unexpected high precision for the estimated coverage of land cover classes with relatively high coverage.     

Relationship between land use classification and grass shrimp Palaemonetes spp. population metrics in coastal watersheds

Estuaries in the southeastern USA have experienced increased loading of contaminants from nonpoint source runoff as well as changes in habitat (e.g., loss of wetlands) due to urbanization. These changes may pose significant risks to estuarine fauna, including crustaceans. Several studies have shown relationships between land use classification and levels of stress in estuarine populations. The grass shrimp of the genus Palaemonetes is one of the dominant species found in estuarine tidal creeks, accounting for more than 50 % of all macropelagic fauna. Grass shrimp populations were sampled monthly for 3 years at six estuarine creeks on Kiawah Island, SC. Creek watersheds were estimated using National Aerial Photograph Program color infrared and low-altitude true color aerial photography combined with in situ differentially corrected global positioning system mapping of engineered features. Land classifications delineated included water, marsh, buildings, roads, and lawns. Pairwise comparisons for grass shrimp densities among sites showed significant differences on an annual and seasonal basis. Significant relationships (p?<?0.05) between land class variables and grass shrimp density were identified both annually and seasonally. These findings suggest an influence of land use on Palaemonetes spp. populations.     

A new multiscale approach for monitoring vegetation using remote sensing-based indicators in laboratory, field, and landscape

Remote sensing is an important tool for studying patterns in surface processes on different spatiotemporal scales. However, differences in the spatiospectral and temporal resolution of remote sensing data as well as sensor-specific surveying characteristics very often hinder comparative analyses and effective up- and downscaling analyses. This paper presents a new methodical framework for combining hyperspectral remote sensing data on different spatial and temporal scales. We demonstrate the potential of using the “One Sensor at Different Scales” (OSADIS) approach for the laboratory (plot), field (local), and landscape (regional) scales. By implementing the OSADIS approach, we are able (1) to develop suitable stress-controlled vegetation indices for selected variables such as the Leaf Area Index (LAI), chlorophyll, photosynthesis, water content, nutrient content, etc. over a whole vegetation period. Focused laboratory monitoring can help to document additive and counteractive factors and processes of the vegetation and to correctly interpret their spectral response; (2) to transfer the models obtained to the landscape level; (3) to record imaging hyperspectral information on different spatial scales, achieving a true comparison of the structure and process results; (4) to minimize existing errors from geometrical, spectral, and temporal effects due to sensor- and time-specific differences; and (5) to carry out a realistic top- and downscaling by determining scale-dependent correction factors and transfer functions. The first results of OSADIS experiments are provided by controlled whole vegetation experiments on barley under water stress on the plot scale to model LAI using the vegetation indices Normalized Difference Vegetation Index (NDVI) and green NDVI (GNDVI). The regression model ascertained from imaging hyperspectral AISA-EAGLE/HAWK (DUAL) data was used to model LAI. This was done by using the vegetation index GNDVI with an R ² of 0.83, which was transferred to airborne hyperspectral data on the local and regional scales. For this purpose, hyperspectral imagery was collected at three altitudes over a land cover gradient of 25 km within a timeframe of a few minutes, yielding a spatial resolution from 1 to 3 m. For all recorded spatial scales, both the LAI and the NDVI were determined. The spatial properties of LAI and NDVI of all recorded hyperspectral images were compared using semivariance metrics derived from the variogram. The first results show spatial differences in the heterogeneity of LAI and NDVI from 1 to 3 m with the recorded hyperspectral data. That means that differently recorded data on different scales might not sufficiently maintain the spatial properties of high spatial resolution hyperspectral images.     

Growing stock and woody biomass assessment in Asola-Bhatti Wildlife Sanctuary,Delhi, India

Biomass is an important entity to understand the capacity of an ecosystem to sequester and accumulate carbon over time. The present study, done in collaboration with the Delhi Forest Department, focused on the estimation of growing stock and the woody biomass in the so-called lungs of Delhi—the Asola-Bhatti Wildlife Sanctuary in northern Aravalli hills. The satellite-derived vegetation strata were field-inventoried using stratified random sampling procedure. Growing stock was calculated for the individual sample plots using field data and species-specific volume equations. Biomass was estimated from the growing stock and the specific gravity of the wood. Among the four vegetation types, viz. Prosopis juliflora, Anogeissus pendula, forest plantation and the scrub, the P. juliflora was found to be the dominant vegetation in the area, covering 23.43 km² of the total area. The study revealed that P. juliflora forest with moderate density had the highest (10.7 m³/ha) while A. pendula forest with moderate density had the lowest (3.6 m³/ha) mean volume. The mean woody biomass was also found to be maximum in P. juliflora forest with moderate density (10.3 t/ha) and lowest in A. pendula forest with moderate density (3.48 t/ha). The total growing stock was estimated to be 20,772.95 m³ while total biomass worked out to be 19,366.83 t. A strong correlation was noticed between the normalized difference vegetation index (NDVI) and the growing stock (R ²?=?0.84)/biomass (R ²?=?0.88). The study demonstrated that growing stock and the biomass of the woody vegetation in Asola-Bhatti Wildlife Sanctuary could be estimated with high accuracy using optical remote sensing data.     

Environmental change monitoring in the arid and semi-arid regions: a case study Al-Basrah Province, Iraq

In recent years, land use/cover dynamic change has become a key subject urgently to be dealt with in the study of global environmental change. This research utilizes the integrated remote sensing and geographic information systems (GIS) in the southern part of Iraq (Basrah Province was taken as a case) to monitor, map, and quantify the environmental change using a 1:250,000 mapping scale. Remote sensing and GIS software were used to classify Landsat TM in 1990 and Landsat ETM+ in 2003 imagery into five land use and land cover (LULC) classes: vegetation land, sand land, urban area, unused land, and water bodies. Supervised classification and normalized difference buildup index, normalized difference vegetation index, normalized difference bare land index, the normalized differential water index, crust index (CI) algorithms, and change detection techniques were adopted in this research and used, respectively, to retrieve its class boundary. An accuracy assessment was performed on the 2003 LULC map to determine the reliability of the map. Finally, GIS software was used to quantify and illustrate the various LULC conversions that took place over the 13-year span of time. The results showed that the urban area, sand lands, and bare lands had increased by the rate of 1.2%, 0.8%, and 0.4% per year, with area expansion from 3,299.1, 4,119.1 km², and 3,201.9 km² in 1990 to 3,794.9, 4,557.7, and 3,351.7 km² in 2003, respectively. While the vegetation cover and water body classes were about 43.5% in 1990, the percentage decreased to about 39.6% in 2003. This study demonstrates the effectiveness of the remote sensing and GIS technologies in detecting, assessing, mapping, and monitoring the environmental changes.     

Impacts of ecological water conveyance on groundwater dynamics and vegetation recovery in the lower reaches of the Tarim River in northwest China

The ecological water conveyance project (EWCP) in the lower reaches of the Tarim River provided a valuable opportunity to study hydro-ecological processes of desert riparian vegetation. Ecological effects of the EWCP were assessed at large spatial and temporal scales based on 13 years of monitoring data. This study analyzed the trends in hydrological processes and the ecological effects of the EWCP. The EWCP resulted in increased groundwater storage—expressed as a general rise in the groundwater table—and improved soil moisture conditions. The change of water conditions also directly affected vegetative cover and the phenology of herbs, trees, and shrubs. Vegetative cover of herbs was most closely correlated to groundwater depth at the last year-end (R?=?0.81), and trees and shrubs were most closely correlated to annual average groundwater depth (R?=?0.79 and 0.66, respectively). The Normalized Difference Vegetation Index (NDVI) responded to groundwater depth on a 1-year time lag. Although the EWCP improved the NDVI, the study area is still sparsely vegetated. The main limitation of the EWCP is that it can only preserve the survival of existing vegetation, but it does not effectively promote the reproduction and regeneration of natural vegetation.     

Preliminary investigation of submerged aquatic vegetation mapping using hyperspectral remote sensing

The use of airborne hyperspectral remote sensing imagery for automated mapping of submerged aquatic vegetation (SAV) in the tidal Potomac River was investigated for near to real-time resource assessment and monitoring. Airborne hyperspectral imagery and field spectrometer measurements were obtained in October of 2000. A spectral library database containing selected ground-based and airborne sensor spectra was developed for use in image processing. The spectral library is used to automate the processing of hyperspectral imagery for potential real-time material identification and mapping. Field based spectra were compared to the airborne imagery using the database to identify and map two species of SAV (Myriophyllum spicatum and Vallisneria americana). Overall accuracy of the vegetation maps derived from hyperspectral imagery was determined by comparison to a product that combined aerial photography and field based sampling at the end of the SAV growing season. The algorithms and databases developed in this study will be useful with the current and forthcoming space-based hyperspectral remote sensing systems.     

Assessing and monitoring semi-arid shrublands using object-based image analysis and multiple endmember spectral mixture analysis

Arid and semi-arid shrublands have significant biological and economical values and have been experiencing dramatic changes due to human activities. In California, California sage scrub (CSS) is one of the most endangered plant communities in the US and requires close monitoring in order to conserve this important biological resource. We investigate the utility of remote-sensing approaches—object-based image analysis applied to pansharpened QuickBird imagery (QBPS/OBIA) and multiple endmember spectral mixture analysis (MESMA) applied to SPOT imagery (SPOT/MESMA)—for estimating fractional cover of true shrub, subshrub, herb, and bare ground within CSS communities of southern California. We also explore the effectiveness of life-form cover maps for assessing CSS conditions. Overall and combined shrub cover (i.e., true shrub and subshrub) were estimated more accurately using QBPS/OBIA (mean absolute error or MAE, 8.9 %) than SPOT/MESMA (MAE, 11.4 %). Life-form cover from QBPS/OBIA at a 25?×?25 m grid cell size seems most desirable for assessing CSS because of its higher accuracy and spatial detail in cover estimates and amenability to extracting other vegetation information (e.g., size, shape, and density of shrub patches). Maps derived from SPOT/MESMA at a 50?×?50 m scale are effective for retrospective analysis of life-form cover change because their comparable accuracies to QBPS/OBIA and availability of SPOT archives data dating back to the mid-1980s. The framework in this study can be applied to other physiognomically comparable shrubland communities.     

Mapping broom snakeweed through image analysis of color-infrared photography and digital imagery

A study was conducted on a south Texas rangeland area to evaluate aerial color-infrared (CIR) photography and CIR digital imagery combined with unsupervised image analysis techniques to map broom snakeweed [Gutierrezia sarothrae (Pursh.) Britt. and Rusby]. Accuracy assessments performed on computer-classified maps of photographic images from two sites had mean producer's and user's accuracies for broom snakeweed of 98.3 and 88.3%, respectively; whereas, accuracy assessments performed on classified maps from digital images of the same two sites had mean producer's and user's accuracies for broom snakeweed of 98.3 and 92.8%, respectively. These results indicate that CIR photography and CIR digital imagery combined with image analysis techniques can be used successfully to map broom snakeweed infestations on south Texas rangelands.     

Post-classification comparison of land cover using multitemporal Landsat and ASTER imagery: the case of Kahramanmaraş, Turkey

This study assessed land cover (LC) changes in Kahramanmara? (K.Mara?) and its environs by using multitemporal Landsat and ASTER imagery, respectively belong to 1989, 2000 and 2004. A priori defined nine land cover classes in the classification scheme were urban and built-up, forest, sparsely vegetated areas, grassland, vegetated stream beds, unvegetated stream beds, bare areas, crop fields, and water bodies. Individual classifications were employed using the combination of both unsupervised and supervised classification methods. Iterative Self Organizing Data Analysis (ISODATA) was used to reduce spectral variation in the scenes arising from complex pattern of crop fields. Maximum Likelihood classifier was used in the LC classification of the individual images. Image pairs of consecutive dates were compared by overlaying the thematic LC maps and cross-tabulating the LC statistics. Urbanization and expansion of agriculture were the major reasons for the dramatic LC conversions. The amount of conversion from crop fields to water occurred as large as 927.67 ha, accounting for 73% of the total land-to-water conversion. Conversions to agriculture have mainly been occurred from grasslands and sparsely vegetated areas as large as 1,314.95 and 1,325.84 ha, respectively. Urban coverage doubled in this period as a result of 1,443.45 ha of increase. Urban area increased in the second period from 2,920 to 3,526 ha. Conversions to agriculture occurred at high amounts. A total of 1,075.79 ha area changed from sparsely vegetated areas to crop fields. A landscape-level environmental monitoring scheme based on satellite remote sensing was proposed for effective environmental resource management.     

Weed mapping in cotton using ground-based sensors and GIS

Site-specific weed management presupposes the careful monitoring and mapping of weed infestation areas. Cut-edge sensor technologies coupled with geographical information systems (GIS) provide the means for reliable decision-making concerning weed management even in sub-field level. In present research, two different spectral sensing systems were engaged in order to digitally map weed patches as grown in four different cotton fields in Central Greece. The systems used were a set of two Crop Circle multispectral sensors ACS-430 and a digital camera Nikon D300S. The spaces between cotton rows were scanned and photographed with the two systems accordingly. Raw recorded data were stored and analyzed in GIS environment producing spatially interpolated maps of red-edge normalized difference vegetation index (NDVI) and weed cover percentage values. Both mapping approaches were satisfactorily related to weed distribution as occurred in the fields; however, the photographic method tended to underestimate weed populations. Correlation of red-edge NDVI and weed cover values, at the points where photographs were taken, as revealed by Pearson’s correlation coefficient was high (r?>?0.83) and statistically significant at the 0.01 level. A first-degree linear equation adequately modeled (R²?>?0.7) the between value pair relations, strengthening the validity of the two methodologies in spatially monitoring weed patches. The methodologies and the technologies used in the study can be used for yearly mapping weed flora in cotton cultivation and potentially constitute a means of rationalizing herbicide application in terms of doses and spatio-temporal decision-making.     

Vegetation NDVI Linked to Temperature and Precipitation in the Upper Catchments of Yellow River

Vegetation in the upper catchment of Yellow River is critical for the ecological stability of the whole watershed. The dominant vegetation cover types in this region are grassland and forest, which can strongly influence the eco-environmental status of the whole watershed. The normalized difference vegetation index (NDVI) for grassland and forest has been calculated and its daily correlation models were deduced by Moderate Resolution Imaging Spectroradiometer products on 12 dates in 2000, 2003, and 2006. The responses of the NDVI values with the inter-annual grassland and forest to three climatic indices (i.e., yearly precipitation and highest and lowest temperature) were analyzed showing that, except for the lowest temperature, the yearly precipitation and highest temperature had close correlations with the NDVI values of the two vegetation communities. The value of correlation coefficients ranged from 0.815 to 0.951 (p?<?0.01). Furthermore, the interactions of NDVI values of vegetation with the climatic indicators at monthly interval were analyzed. The NDVI of vegetation and three climatic indices had strong positive correlations (larger than 0.733, p?<?0.01). The monthly correlations also provided the threshold values for the three climatic indictors, to be used for simulating vegetation growth grassland under different climate features, which is essential for the assessment of the vegetation growth and for regional environmental management.     

Applying Satellite Imagery to Triage Assessment of Ecosystem Health

Considerable evidence documents that certain changes in vegetation and soils result in irreversibly degraded rangeland ecosystems. We used Advanced Very High Resolution Radiometer (AVHRR) imagery to develop calibration patterns of change in the Normalized Difference Vegetation Index (NDVI) over the growing season for selected sites for which we had ground data and historical data characterizing these sites as irreversibly degraded. We used the NDVI curves for these training sites to classify and map the irreversibly degraded rangelands in southern New Mexico. We composited images into four year blocks: 1988–1991, 1989–1992, and 1990–1993. The overlap in pixels classified as irreversibly degraded ranged from 42.6% to 84.3% in year block comparisons. Quantitative data on vegetation composition and cover were collected at 13 sites within a small portion of the study area. Wide coverage reconnaissance of boundaries between vegetation types was also conducted for comparisons with year block maps. The year block 1988–1991 provided the most accurate delineation of degraded areas. The rangelands of southern New Mexico experienced above average precipitation from 1990–1993. The above average precipitation resulted in spatially variable productivity of ephemeral weedy plants on the training sites and degraded rangelands which resulted in much smaller areas classified as irreversibly degraded. We selected imagery for a single year, 1989, which was characterized by the absence of spring annual plant production in order to eliminate the confounding effect of reflectance from annual weeds. That image analysis classified more than 20% of the rangelands as irreversibly degraded because areas with shrub-grass mosaic were included in the degraded classification. The single year image included more than double the area classified as irreversibly degraded by the year blocks. AVHRR imagery can be used to make triage assessments of irreversibly degraded rangeland but such assessment requires understanding productivity patterns and variability across the landscapes of the region and careful selection of the years from which imagery is chosen.     

Assessment of some remote sensing techniques used to detect land use/land cover changes in South-East Transilvania,Romania

This paper assesses the image differencing technique for the Normalized Difference Vegetation Index (NDVI), the second principal component (PC2), and the TM 4 band (TM 4), as well as the post-classification comparison (PCC) in order to analyze the land use/land cover changes in the South-East Transilvania, Romania. The analysis was performed using two frames from Landsat 5 TM satellite images acquired on August 5, 1993 and July 24, 2009. After applying the NDVI, PC2, and TM 4 image differencing techniques, the images obtained were transformed into change/no change maps. The thresholds identified to highlight the changes were set at 0.6 s for NDVI and 0.7 s for PC2 and TM 4. Before applying the PCC technique, the satellite images were classified through the supervised classification method. The overall accuracy obtained was 85.91 % and the kappa statistics 0.8249 for 1993, 88.18 % and 0.8497 for 2009, respectively. The assessment of the changes detection methods in the studied area shows that the first place is occupied by NDVI image differencing with an overall accuracy of 83.80 %, followed by PCC method with 83.20 %, PC2 difference with an overall accuracy of 81.60 %, and TM 4 difference with an overall accuracy of 79.40 %.     

Evaluating the difference between the normalized difference vegetation index and net primary productivity as the indicators of vegetation vigor assessment at landscape scale

Both the net primary productivity (NPP) and the normalized difference vegetation index (NDVI) are commonly used as indicators to characterize vegetation vigor, and NDVI has been used as a surrogate estimator of NPP in some cases. To evaluate the reliability of such surrogation, here we examined the quantitative difference between NPP and NDVI in their outcomes of vegetation vigor assessment at a landscape scale. Using Landsat ETM+ data and a process model, the Boreal Ecosystem Productivity Simulator, NPP distribution was mapped at a resolution of 90 m, and total NDVI during the growing season was calculated in Heihe River Basin, Northwest China in 2002. The results from a comparison between the NPP and NDVI classification maps show that there existed a substantial difference in terms of both area and spatial distribution between the assessment outcomes of these two indicators, despite that they are strongly correlated. The degree of difference can be influenced by assessment schemes, as well as the type of vegetation and ecozone. Overall, NDVI is not a good surrogate of NPP as the indicators of vegetation vigor assessment in the study area. Nonetheless, NDVI could serve as a fairish surrogate indicator under the condition that the target region has low vegetation cover and the assessment has relatively coarse classification schemes (i.e., the class number is small). It is suggested that the use of NPP and NDVI should be carefully selected in landscape assessment. Their differences need to be further evaluated across geographic areas and biomes.