Water content estimation from hyperspectral images and MODIS indexes in Southeastern Arizona

Hyperspectral water retrievals from AVIRIS data, equivalent water thickness (EWT), were compared to in situ leaf water content and LAI measurements at a semiarid site in southeastern Arizona. Retrievals of EWT showed good correlation with field canopy water content measurements. Statistical analysis suggested that EWT was significantly different among seven community types, from savanna to agriculture. Four band-ratio indexes (NDVI, EVI, NDWI, and NDII) were derived from MODIS showing strong spatial agreement between maps of AVIRIS EWT and MODIS indexes, and good statistical agreement for the range of habitats at the site. Temporal patterns of these four indexes in all vegetation communities except creosote bush and agriculture showed distinct seasonal patterns that responded to the timing and amount of precipitation. Moreover, these time series captured different ecological responses among the different vegetation communities.     

Remote sensing of vegetation water content from equivalent water thickness using satellite imagery

Vegetation water content (VWC) is one of the most important parameters for the successful retrieval of soil moisture content from microwave data. Normalized Difference Infrared Index (NDII) is a widely-used index to remotely sense Equivalent Water Thickness (EWT) of leaves and canopies; however, the amount of water in the foliage is a small part of total VWC. Sites of corn (Zea mays), soybean (Glycine max), and deciduous hardwood woodlands were sampled to estimate EWT and VWC during the Soil Moisture Experiment 2005 (SMEX05) near Ames, Iowa, USA. Using a time series of Landsat 5 Thematic Mapper, Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Advanced Wide Field Sensor (AWiFS) imagery, NDII was related to EWT with R² of 0.85; there were no significant differences among land-cover types. Furthermore, EWT was linearly related to VWC with R² of 0.87 for corn and 0.48 for soybeans, with a significantly larger slope for corn. The 2005 land-cover classification product from the USDA National Agricultural Statistics Service had an overall accuracy of 92% and was used to spatially distribute VWC over the landscape. SMEX05 VWC versus NDII regressions were compared with the regressions from the Soil Moisture Experiment 2002 (SMEX02), which was conducted in the same study area. No significant difference was found between years for corn (P = 0.13), whereas there was a significant difference for soybean (P = 0.04). Allometric relationships relate the size of one part of a plant to the sizes of other parts, and may be the result from the requirements of structural support or material transport. Relationships between NDII and VWC are indirect, NDII is related to canopy EWT, which in turn is allometrically related to VWC.     

Comparison of principal component inversion with VI-empirical approach for LAI estimation using simulated reflectance data

A simulated canopy reflectance dataset for a total of six channels in visible, near-infrared (NIR) and shortwave-infrared (SWIR) region, corresponding to Landsat Thematic Mapper (TM) was generated using the PROSAIL (PROSPECT+SAIL) model and a range of Leaf Area Index (LAI), soil backgrounds, leaf chlorophyll, leaf inclination and viewing geometry inputs. This dataset was used to develop and evaluate approaches for LAI estimation, namely, standard two-band nonlinear empirical vegetation index (VI)–LAI formulation (using Normalized Difference Vegetation Index/simple ratio (NDVI/SR)) and a multi-band principal component inversion (PCI) approach. The analysis indicated that the multi-band PCI approach had a smaller rms error (RMSE=0.380) than the NDVI and SR approaches (RMSE=2.28, 0.88), for an independently generated test dataset.     

Analysis of vegetation recovery surrounding a restored wetland using the normalized difference infrared index (NDII) and normalized difference vegetation index (NDVI)

Watershed restoration efforts seek to rejuvenate vegetation, biological diversity, and land productivity at Cienega San Bernardino, an important wetland in southeastern Arizona and northern Sonora, Mexico. Rock detention and earthen berm structures were built on the Cienega San Bernardino over the course of four decades, beginning in 1984 and continuing to the present. Previous research findings show that restoration supports and even increases vegetation health despite ongoing drought conditions in this arid watershed. However, the extent of restoration impacts is still unknown despite qualitative observations of improvement in surrounding vegetation amount and vigor. We analyzed spatial and temporal trends in vegetation greenness and soil moisture by applying the normalized difference vegetation index (NDVI) and normalized difference infrared index (NDII) to one dry summer season Landsat path/row from 1984 to 2016. The study area was divided into zones and spectral data for each zone was analyzed and compared with precipitation record using statistical measures including linear regression, Mann–Kendall test, and linear correlation. NDVI and NDII performed differently due to the presence of continued grazing and the effects of grazing on canopy cover; NDVI was better able to track changes in vegetation in areas without grazing while NDII was better at tracking changes in areas with continued grazing. Restoration impacts display higher greenness and vegetation water content levels, greater increases in greenness and water content through time, and a decoupling of vegetation greenness and water content from spring precipitation when compared to control sites in nearby tributary and upland areas. Our results confirm the potential of erosion control structures to affect areas up to 5 km downstream of restoration sites over time and to affect 1 km upstream of the sites.     

Characterization of pasture biophysical properties and the impact of grazing intensity using remotely sensed data

Remote sensing has the potential of improving our ability to map and monitor pasture degradation. Pasture degradation is one of the most important problems in the Amazon, yet the manner in which grazing intensity, edaphic conditions and land‐use age impact pasture biophysical properties, and our ability to monitor them using remote sensing is poorly known. We evaluate the connection between field grass biophysical measures and remote sensing, and investigate the impact of grazing intensity on pasture biophysical measures in Rondônia, in the Brazilian Amazon. Above ground biomass, canopy water content and height were measured in different pasture sites during the dry season. Using Landsat Thematic Mapper (TM) data, four spectral vegetation indices and fractions derived from spectral mixture analysis, i.e., Non‐Photosynthetic Vegetation (NPV), Green Vegetation (GV), Soil, Shade, and NPV + Soil, were calculated and compared to field grass measures. For grazed pastures under dry conditions, the Normalized Difference Infrared Index (NDII5 and NDII7), had higher correlations with the biophysical measures than the Normalized Difference Vegetation Index (NDVI) and the Soil‐Adjusted Vegetation Index (SAVI). NPV had the highest correlations with all field measures, suggesting this fraction is a good indicator of pasture characteristics in Rondônia. Pasture height was correlated to the Shade fraction. A conceptual model was built for pasture biophysical change using three fractions, i.e., NPV, Shade and GV to characterize possible pasture degradation processes in Rondônia. Based upon field measures, grazing intensity had the most significant impact on pasture biophysical properties compared to soil order and land‐use age. The impact of grazing on pastures in the dry season could be potentially measured by using remotely sensed measures such as NPV.     

Upscaling ground observations of vegetation water content, canopy height, and leaf area index during SMEX02 using aircraft and Landsat imagery

Microwave-based remote sensing algorithms for mapping soil moisture are sensitive to water contained in surface vegetation at moderate levels of canopy cover. Correction schemes require spatially distributed estimates of vegetation water content at scales comparable to that of the microwave sensor footprint (10¹ to 10⁴ m). This study compares the relative utility of high-resolution (1.5 m) aircraft and coarser-resolution (30 m) Landsat imagery in upscaling an extensive set of ground-based measurements of canopy biophysical properties collected during the Soil Moisture Experiment of 2002 (SMEX02) within the Walnut Creek Watershed. The upscaling was accomplished using expolinear relationships developed between spectral vegetation indices and measurements of leaf area index, canopy height, and vegetation water content. Of the various indices examined, a Normalized Difference Water Index (NDWI), derived from near- and shortwave-infrared reflectances, was found to be least susceptible to saturation at high levels of leaf area index. With the aircraft data set, which did not include a short-wave infrared water absorption band, the Optimized Soil Adjusted Vegetation Index (OSAVI) yielded best correlations with observations and highest saturation levels. At the observation scale (10 m), LAI was retrieved from both NDWI and OSAVI imagery with an accuracy of 0.6, vegetation water content at 0.7 kg m⁻², and canopy height to within 0.2 m. Both indices were used to estimate field-scale mean canopy properties and variability for each of the intensive soil-moisture-sampling sites within the watershed study area. Results regarding scale invariance over the SMEX02 study area in transformations from band reflectance and vegetation indices to canopy biophysical properties are also presented.     

Leaf area index mapping in northern Canada

Leaf area index (LAI) is an important structural vegetation parameter that is commonly derived from remotely sensed data. It has been used as a reliable indicator for vegetation's cover, status, health and productivity. In the past two decades, various Canada-wide LAI maps have been generated by the Canada Centre for Remote Sensing (CCRS). These products have been produced using a variety of very coarse satellite data such as those from SPOT VGT and NOAA AVHRR satellite data. However, in these LAI products, the mapping of the Canadian northern vegetation has not been performed with field LAI measurements due in large part to scarce in situ measurements over northern biomes. The coarse resolution maps have been extensively used in Canada, but finer resolution LAI maps are needed over the northern Canadian ecozones, in particular for studying caribou habitats and feeding grounds.

In this study, a new LAI algorithm was developed with particular emphasis over northern Canada using a much finer resolution of remotely sensed data and in situ measurements collected over a wide range of northern arctic vegetation. A statistical relationship was developed between the in situ LAI measurements collected over vegetation plots in northern Canada and their corresponding pixel spectral information from Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data. Furthermore, all Landsat TM and ETM+ data have been pre-normalized to NOAA AVHRR and SPOT VGT data from the growing season of 2005 to reduce any seasonal or temporal variations. Various spectral vegetation indices developed from the Landsat TM and ETM?+?data were analysed in this study. The reduced simple ratio index (RSR) was found to be the most robust and an accurate estimator of LAI for northern arctic vegetation. An exponential relationship developed using the Theil–Sen regression technique showed an R ² of 0.51 between field LAI measurement and the RSR. The developed statistical relationship was applied to a pre-existing Landsat TM 250 m resolution mosaic for northern Canada to produce the final LAI map for northern Canada ecological zones. Furthermore, the 250 m resolution LAI estimates, per ecological zone, were almost generally lower than those of the CCRS Canada-wide VGT LAI maps for the same ecozones. Validation of the map with LAI field data from the 2008 season, not used in the derivation of the algorithm, shows strong agreement between the in situ LAI measurement values and the map-estimated LAI values.     

Quantitative accuracy assessment of contextually filtered classified images

Relative radiance recorded by the Landsat Thematic Mapper (TM) for l–1O-year-old plantations of oil palm (Elaeis guineensis Jacq.) in Sabah, Malaysia, was negatively correlated with stand age. Remotely sensed response is determined by biophysical variables related to age since field planting, notably leaf area, canopy architecture and progressive masking of ground cover vegetation. The relation was asymptotic. Age and age-related variables will be most accurately inferred from TM data for stands under 5-years old, especially using short wave infrared and thermal bands. Narrower age classes will be required to represent younger stands in image classification.     

Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery

This paper compares the normalized difference vegetation index (NDVI) and percent impervious surface as indicators of surface urban heat island effects in Landsat imagery by investigating the relationships between the land surface temperature (LST), percent impervious surface area (%ISA), and the NDVI. Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data were used to estimate the LST from four different seasons for the Twin Cities, Minnesota, metropolitan area. A map of percent impervious surface with a standard error of 7.95% was generated using a normalized spectral mixture analysis of July 2002 Landsat TM imagery. Our analysis indicates there is a strong linear relationship between LST and percent impervious surface for all seasons, whereas the relationship between LST and NDVI is much less strong and varies by season. This result suggests percent impervious surface provides a complementary metric to the traditionally applied NDVI for analyzing LST quantitatively over the seasons for surface urban heat island studies using thermal infrared remote sensing in an urbanized environment.     

On intra-annual EVI variability in the dry season of tropical forest: A case study with MODIS and hyperspectral data

Considerable controversy is associated with dry season increases in the Enhanced Vegetation Index (EVI), observed using the Moderate Resolution Imaging Spectroradiometer (MODIS), compared with field-based estimates of decreasing plant productivity. Here, we investigate potential causes of intra-annual variability by comparing EVI from mature forest with field-measured Leaf Area Index (LAI) to validate space-based observations. EVI was calculated from 19 nadir and off-nadir Hyperion images in the 2005 dry season, and inspected for consistency with MODIS observations from 2004 to 2009. The objective was to evaluate the possible influence of the view-illumination geometry and of canopy foliage and leaf flush on the EVI. Spectral mixture models were used to evaluate the relationship between EVI and the shade fraction, a measure that varies with pixel brightness. MODIS LAI values were compared with LAI estimated using hemispherical photographs taken in two field campaigns in the dry season. To keep LAI and leaf flush conditions as constant variables and vary solar illumination, we used airborne Hyperspectral Mapper (Hymap) data acquired over mature forest from another region on the same day but with two distinct solar zenith angles (SZA) (29° and 53°). Results showed that intra-annual variability in MODIS and nadir Hyperion EVI in the dry season of tropical forest were driven by solar illumination effects rather than changes in LAI. The reflectance of the MODIS and Hyperion blue, red and near infrared (NIR) bands was higher at the end of the dry season because of the predominance of sunlit canopy components for the sensors due to decreasing SZA from June (44°) to September (26°). Because EVI was highly correlated with the reflectance of the NIR band used to generate it (r of + 0.98 for MODIS and + 0.88 for Hyperion), this vegetation index followed the general NIR pattern, increasing with smaller SZA towards the end of the dry season. Hyperion EVI was inversely correlated with the shade fraction (r = − 0.93). Changes in canopy foliage detected from MODIS LAI data were not consistent with LAI estimates from hemispherical photographs. Although further research is necessary to measure the impact of leaf flush on intra-annual EVI variability in the Querência region, analysis of Hymap data with fixed LAI and leaf flush conditions confirmed the influence of the illumination effects on the EVI.     

Cover The wake of Hurricane Felix

Analytical canopy reflectance (CR) models have reached the level of adequacy that makes it possible to estimate vegetation parameters by inversion of such models. The growing efficiency of algorithms and the increasing power of computers urge the development of procedures for the estimation of vegetation phytometrical parameters on large areas using satellite data and inversion of theoretical CR models. In this article, clusterization of a Landsat Thematic Mapper (TM) quarter scene is performed in the space of spectral signatures, and the CR model is inverted for these clusters. Optical parameters of the atmosphere which are needed for the atmospheric correction are estimated on the same image. The estimated Leaf Area Index (LAI) pattern is in good accordance to the land use map. Estimated LAI and chlorophyll content of forests are systematically biased.     

Estimating live fuel moisture content from remotely sensed reflectance

Fuel moisture content (FMC) is used in forest fire danger models to characterise the moisture status of the foliage. FMC expresses the amount of water in a leaf relative to the amount of dry matter and differs from measures of leaf water content which express the amount of water in a leaf relative to its area. FMC is related to both leaf water content and leaf dry matter content, and the relationships between FMC and remotely sensed reflectance will therefore be affected by variation in both leaf biophysical properties. This paper uses spectral reflectance data from the Leaf Optical Properties EXperiment (LOPEX) and modelled data from the Prospect leaf reflectance model to examine the relationships between FMC, leaf equivalent water thickness (EWT) and a range of spectral vegetation indices (VI) designed to estimate leaf and canopy water content. Significant correlations were found between FMC and all of the selected vegetation indices for both modelled and measured data, but statistically stronger relationships were found with leaf EWT; overall, the water index (WI) was found to be most strongly correlated with FMC. The accuracy of FMC estimation was very low when the global range of FMC was examined, but for a restricted range of 0-100%, FMC was estimated with a root-mean-square error (RMSE) of 15% in the model simulations and 51% with the measured data. The paper shows that the estimation of live FMC from remotely sensed vegetation indices is likely to be problematic when there is variability in both leaf water content and leaf dry matter content in the target leaves. Estimating FMC from remotely sensed data at the canopy level is likely to be further complicated by spatial and temporal variations in leaf area index (LAI). Further research is required to assess the potential of canopy reflectance model inversion to estimate live fuel moisture content where a priori information on vegetation properties may be used to constrain the inversion process.     

双层植被结构冠层光谱特性的理论模拟

利用双层冠层反射率模型ACRM，模拟不同叶面积指数LAI、含水量Cw和结构参数N下，波长是820nm和1600nm冠层反射率的角度分布。结果说明，该模型能准确模拟出“热点”效应 |冠层反射率角度分布对LAI的敏感性小于LAI和Cw以及LAI和N的共同作用。其中LAI和Cw共同作用对波长1 600 nm反射率角度分布的影响非常显著，而LAI和N的共同作用在820 nm略微大于1 600 nm。另外，提取冠层含水量的土壤可调节水分指数SAWI受冠层结构的影响也较大。今后在模型选取中应该更好的考虑冠层结构影响。     

Dynamic analysis of Robinia pseudoacacia forest health levels from 1995 to 2013 in the Yellow River Delta,China using multitemporal Landsat imagery

The Robinia pseudoacacia forest in the Yellow River Delta (YRD), China, was planted in the 1970s and has continuously suffered dieback and mortality since the 1990s. Timely and accurate information on forest growth and forest condition and its dynamic change as well is essential for assessing and developing effective management strategies. In this study, multitemporal Landsat imagery was used to analyze and monitor changes of the R. pseudoacacia forest in the YRD from 1995 to 2013. To do so, Landsat image band reflectance, three fraction images calculated by using a multiple endmember spectral mixture analysis (MESMA) method, and four vegetation indices (VIs) were used to discriminate three health levels of R. pseudoacacia forest in years 1995, 2007, and 2013 with a random forest (RF) classifier. The four VIs include a difference infrared index (DII) developed in this study, normalized difference vegetation index, soil-adjusted vegetation index, and normalized difference infrared index (NDII), all of which were computed from Landsat Thematic Mapper and Operational Land Imager multispectral (MS) bands. The dynamic changes of the forest health levels during the periods of 1995–2007 and 2007–2013 were analysed. The analysis results demonstrate that three fraction images created by MESMA method and four VIs were powerful in separating the three forest health levels. In addition to the Landsat MS bands, the additional three fraction images increased the classification accuracy by 14?20%; if coupled with the four VIs, the overall accuracy was further increased by 5?6%. According to the importance values calculated by RF classifier for all input features, the DII vegetation index was the second effective feature, outperforming NDII. From 1995 to 2013, a total of 2615 ha of forest in the study area suffered from mortality or loss.     

Land cover mapping in support of LAI and FPAR retrievals from EOS-MODIS and MISR: Classification methods and sensitivities to errors

Land cover maps are used widely to parameterize the biophysical properties of plant canopies in models that describe terrestrial biogeochemical processes. In this paper, we describe the use of supervised classification algorithms to generate land cover maps that characterize the vegetation types required for Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FPAR) retrievals from MODIS and MISR. As part of this analysis, we examine the sensitivity of remote sensing-based retrievals of LAI and FPAR to land cover information used to parameterize vegetation canopy radiative transfer models. Specifically, a decision tree classification algorithm is used to generate a land cover map of North America from Advanced Very High Resolution Radiometer (AVHRR) data with 1 km spatial resolution using a six-biome classification scheme. To do this, a time series of normalized difference vegetation index data from the AVHRR is used in association with extensive site-based training data compiled using Landsat Thematic Mapper (TM) and ancillary map sources. Accuracy assessment of the map produced via decision tree classification yields a cross-validated map accuracy of 73%. Results comparing LAI and FPAR retrievals using maps from different sources show that disagreement in land cover labels generally do not translate into strong disagreement in LAI and FPAR maps. Further, the main source of disagreement in LAI and FPAR maps can be attributed to specific biome classes that are characterized by a continuum of fractional cover and canopy structure.     

Measurement and comparison of Leaf Area Index estimators derived from satellite remote sensing techniques

Leaf Area Index (LAI) is an important biophysical characteristic of vegetation that is directly related to rates of atmospheric gas exchange, biomass partitioning, and productivity. Mapping and monitoring LAI over scales from landscapes to regions is essential for understanding medium-scale biophysical properties and how these properties affect biogeochemical cycling, biomass accumulation, and primary productivity. This study developed and verified several models to estimate LAI using in situ field measurements, Landsat Thematic Mapper imagery, vegetation indices, simple and multiple regression, and artificial neural networks (ANNs). It was shown that while multiple band regression and regression with individual vegetation indices can estimate LAI, the most accurate way to estimate regional scale LAI is to train an ANN using in situ LAI data and remote sensing brightness values.     

Leaf area index for biomes of the Eastern Arc Mountains: Landsat and SPOT observations along precipitation and altitude gradients

Understanding of mechanisms underlying carbon flux dynamics in the Eastern Arc Mountains and their catchment areas is lacking, due to data shortage (e.g. biome specific canopy structure) and spatial heterogeneity of tropical ecosystems. This study focuses on documenting leaf area index (LAI) for the main biomes in the Eastern Arc Mountains and their surroundings. In situ optical instruments, i.e. hemispherical photography and a SunScan device, were used to acquire ground LAI measurements. Spectral vegetation indices (VIs) extracted from Landsat Enhanced Thematic Mapper (ETM +) and Système Probatoire d'Observation de la Terre (SPOT) reflectance data were used, along with mean annual precipitation (MAP), as explanatory variables of LAI variation. The results indicate that LAI significantly increases with increasing MAP for woody biomes. Implementing long-term MAP as a second predictor variable into the VI–LAI models significantly improved LAI predictions by up to 10% using the normalised difference vegetation index (NDVI), modified soil adjusted vegetation index (MSAVI 2) and 2-band enhanced vegetation index (EVI 2). Varying forest disturbances and agricultural management practises may have contributed to observed discrepancies of LAI with MAP across biomes. The importance of altitudinal gradients is yet to be explained fully with more study required. However, LAI appears to be higher in low-altitude forests compared to forests at higher altitudes. Our results indicate that SPOT and Landsat-derived VIs, in combination with long-term MAP, may be a suitable tool to develop landscape maps of LAI in Eastern Africa. This study also presents the in situ LAI measurements for further validation of global products for areas that are currently under-represented in Earth Observation (EO) global validation networks.     

A non-parametric,supervised classification of vegetation types on the Kaibab National Forest using decision trees

Traditional land classification techniques for large areas that use Landsat Thematic Mapper (TM) imagery are typically limited to the fixed spatial resolution of the sensors (30 m). However, the study of some ecological processes requires land cover classifications at finer spatial resolutions. We model forest vegetation types on the Kaibab National Forest (KNF) in northern Arizona to a 10-m spatial resolution with field data, using topographical information and Landsat TM imagery as auxiliary variables. Vegetation types were identified by clustering the field variables total basal area and proportion of basal area by species, and then using a decision tree based on auxiliary variables to predict vegetation types. Vegetation types modelled included pinyon-juniper, ponderosa pine, mixed conifer, spruce- and deciduous-dominated mixes, and openings. To independently assess the accuracy of the final vegetation maps using reference data from different sources, we used a post-stratified, multivariate composite estimator. Overall accuracy was 74.5% (Kappa statistic = 49.9%). Sources of error included differentiating between mixed conifer and spruce-dominated types and between openings in the forest and deciduous-dominated mixes. Overall, our non-parametric classification method successfully identified dominant vegetation types on the study area at a finer spatial resolution than can typically be achieved using traditional classification techniques.