Vegetation water content during SMEX04 from ground data and Landsat 5 Thematic Mapper imagery

Vegetation water content is an important parameter for retrieval of soil moisture from microwave data and for other remote sensing applications. Because liquid water absorbs in the shortwave infrared, the normalized difference infrared index (NDII), calculated from Landsat 5 Thematic Mapper band 4 (0.76-0.90 μm wavelength) and band 5 (1.55-1.65 μm wavelength), can be used to determine canopy equivalent water thickness (EWT), which is defined as the water volume per leaf area times the leaf area index (LAI). Alternatively, average canopy EWT can be determined using a landcover classification, because different vegetation types have different average LAI at the peak of the growing season. The primary contribution of this study for the Soil Moisture Experiment 2004 was to sample vegetation for the Arizona and Sonora study areas. Vegetation was sampled to achieve a range of canopy EWT; LAI was measured using a plant canopy analyzer and digital hemispherical (fisheye) photographs. NDII was linearly related to measured canopy EWT with an R² of 0.601. Landcover of the Arizona, USA, and Sonora, Mexico, study areas were classified with an overall accuracy of 70% using a rule-based decision tree using three dates of Landsat 5 Thematic Mapper imagery and digital elevation data. There was a large range of NDII per landcover class at the peak of the growing season, indicating that canopy EWT should be estimated directly using NDII or other shortwave-infrared vegetation indices. However, landcover classifications will still be necessary to obtain total vegetation water content from canopy EWT and other data, because considerable liquid water is contained in the non-foliar components of vegetation.     

A comparison of empirical and neural network approaches for estimating corn and soybean leaf area index from Landsat ETM+ imagery

Plant foliage density expressed as leaf area index (LAI) is used in many ecological, meteorological, and agronomic models, and as a means of quantifying crop spatial variability for precision farming. LAI retrieval using spectral vegetation indices (SVI) from optical remotely sensed data usually requires site-specific calibration values from the surface or the use of within-scene image information without surface calibrations to invert radiative transfer models. An evaluation of LAI retrieval methods was conducted using (1) empirical methods employing the normalized difference vegetation index (NDVI) and a new SVI that uses green wavelength reflectance, (2) a scaled NDVI approach that uses no calibration measurements, and (3) a hybrid approach that uses a neural network (NN) and a radiative transfer model without site-specific calibration measurements. While research has shown that under a variety of conditions NDVI is not optimal for LAI retrieval, its continued use for remote sensing applications and in analysis seeking to develop improved parameter retrieval algorithms based on NDVI suggests its value as a “benchmark” or standard against which other methods can be compared. Landsat-7 ETM+ data for July 1 and July 8 from the Soil Moisture EXperiment 2002 (SMEX02) field campaign in the Walnut Creek watershed south of Ames, IA, were used for the analysis. Sun photometer data collected from a site within the watershed were used to atmospherically correct the imagery to surface reflectance. LAI validation measurements of corn and soybeans were collected close to the dates of the Landsat-7 overpasses. Comparable results were obtained with the empirical SVI methods and the scaled SVI method within each date. The hybrid method, although promising, did not account for as much of the variability as the SVI methods. Higher atmospheric optical depths for July 8 leading to surface reflectance errors are believed to have resulted in overall poorer performance for this date. Use of SVIs employing green wavelengths, improved method for the definition of image minimum and maximum clusters used by the scaled NDVI method, and further development of a soil reflectance index used by the hybrid NN approach are warranted. More importantly, the results demonstrate that reasonable LAI estimates are possible using optical remote sensing methods without in situ, site-specific calibration measurements.     

Evaluating image-based estimates of leaf area index in boreal conifer stands over a range of scales using high-resolution CASI imagery

Leaf area index (LAI) is an important surface biophysical parameter as a measure of vegetation cover, vegetation productivity, and as an input to ecosystem process models. Recently, a number of coarse-scale (1-km) LAI maps have been generated over large regions including the Canadian boreal forest. This study focuses on the production of fine-scale (≤30-m) LAI maps using the forest light interaction model-clustering (FLIM-CLUS) algorithm over selected boreal conifer stands and the subsequent comparison of the fine-scale maps to coarse-scale LAI maps synthesized from Landsat TM imagery. The fine-scale estimates are validated using surface LAI measurements to give relative root mean square errors of under 7% for jack pine sites and under 14% for black spruce sites. In contrast, finer scale site mean LAI ranges between 49% and 86% of the mean of surface estimates covering only part of the sites and 54% to 110% of coarse-scale site mean LAI. Correlations between fine-scale and coarse-scale estimates range from near 0.5 for 30-m coarse-scale images to under 0.3 to 1-km coarse-scale images but increase to near 0.90 after imposing fine-scale zero LAI areas in coarse-scale estimates. The increase suggests that coarse-scale image-based LAI estimates require consideration of sub-pixel open areas. Both FLIM-CLUS and coarse-scale site mean LAI are substantially lower than surface estimates over northern sites. The assumption of spatially random residuals in regression-based estimates of LAI may not be valid and may therefore add to local bias errors in estimating LAI remotely. Differences between fine-scale airborne LAI maps and 30-m-scale Landsat TM LAI maps suggests that, for sparse boreal conifer stands, LAI maps produced from Landsat TM alone may not always be sufficient for validation of coarser scale LAI maps. In addition, previous studies may have used biased LAI estimates over the study site. Fine-scale spatial LAI maps offer one means of assessing and correcting for effects of sub-pixel open area patches and for characterising the spatial pattern of residuals in coarse-scale LAI estimates in comparison to the true distribution of LAI on the surface.     

Retrieval of the canopy leaf area index in the BOREAS flux tower sites using linear spectral mixture analysis

This paper reports on the use of linear spectral mixture analysis for the retrieval of canopy leaf area index (LAI) in three flux tower sites in the Boreal Ecosystem-Atmosphere Study (BOREAS) southern study area: Old Black Spruce, Old Jack Pine, and Young Jack Pine (SOBS, SOJP, and SYJP). The data used were obtained by the Compact Airborne Spectrographic Imager (CASI) with a spatial resolution of 2 m in the winter of 1994. The convex geometry method was used to select the endmembers: sunlit crown, sunlit snow, and shadow. Along transects for these flux tower sites, the fraction of sunlit snow was found to have a higher correlation with the field-measured canopy LAI than the fraction of sunlit crown or the fraction of shadow. An empirical equation was obtained to describe the relation between canopy LAI and the fraction of sunlit snow. There is a strong correlation between the estimated LAI and the field-measured LAI along transects (with R² values of 0.54, 0.71, and 0.60 obtained for the SOBS, SYJP, and SOJP sites, respectively). The estimated LAI for the whole tower site is consistent with that obtained by the inversion of a canopy model in our previous study where values of 0.94, 0.92, and 0.63 were obtained for R² for the SOBS, SYJP and SOJP sites, respectively.The CASI 2-m summer data over the SOBS site was also employed to investigate the possibility of deriving canopy LAI from the summer data using linear mixture analysis. At a spatial resolution of 10 m, the correlation between the field-measured LAI and the estimated LAI along transects is small at R² less than 0.3, while R² increases to 0.6 at a spatial resolution of 30 m. The difficulty in canopy LAI retrieval from the summer data at a spatial resolution of 10 m is likely due to the variation of the understory reflectance across the scene, although spatial misregistration of the CASI data used may also be a possible contributing factor.     

SMEX02: Field scale variability, time stability and similarity of soil moisture

Evaluation of air- or space-borne remote sensors measuring soil moisture requires strategic ground-based sampling. As part of the Soil Moisture Experiment 2002 (SMEX02), daily surface soil moisture sampling at 90-140 locations were conducted in four fields in Walnut Creek watershed, Iowa. Various combinations of soils, vegetation, and topography characterize the fields. Depending on the field's characteristics and soil moisture content, 3-32 independent measurements were necessary to capture the field mean volumetric soil moisture with a ±2% bias and 95% confidence interval. Validation of the retrieved soil moisture products from the aircraft microwave instruments using the average of 14 samples per field is more appropriate for dry (<10% volumetric soil moisture) or wet (>25% volumetric soil moisture) range than for intermediate soil moisture range. Time stability analysis showed that an appropriately selected single sampling point could provide similar accuracy across a range of soil moisture conditions. Analyses based on landscape position (depression, hilltop, steep slope, and mild slope) showed that locations with mild slopes consistently exhibit time stable features. Hilltop and steep slope locations consistently underestimated mean field soil moisture. Soils parameters could not be used to identify time stable features as sampling locations with relatively high sand content consistently underestimated the field mean while those locations with relatively high clay content consistently overestimated the field mean. However, the slope position characterization of time stable features was enhanced using soils properties. The mild slope locations having the best time-stable features are those with moderate to moderately high clay content as compare to the field average (28-30% clay).     

Retrieval of soil moisture from passive and active L/S band sensor (PALS) observations during the Soil Moisture Experiment in 2002 (SMEX02)

The Soil Moisture Experiments in 2002 (SMEX02) were conducted in Iowa between June 25th and July 12th, 2002. A major aim of the experiments was examination of existing algorithms for soil moisture retrieval from active and passive microwave remote sensors under high vegetation water content conditions. The data obtained from the passive and active L and S band sensor (PALS) along with physical variables measured by in situ sampling have been used in this study to demonstrate the sensitivity of the instrument to soil moisture and perform soil moisture retrieval using statistical regression and physical modeling techniques. The land cover conditions in the region studied were predominantly soybean and corn crops with average vegetation water contents ranging from 0 to ∼5 kg/m². The PALS microwave sensitivity to soil moisture under these vegetation conditions was investigated for both passive and active measurements. The performance of the PALS instrument and retrieval algorithms has been analyzed, indicating soil moisture retrieval errors of approximately 0.04 g/g gravimetric soil moisture. Statistical regression techniques have been shown to perform satisfactorily with soil moisture retrieval error of around 0.05 g/g gravimetric soil moisture. The retrieval errors were higher for the corn than for the soybean fields due to the higher vegetation water content of the corn crops. However, the algorithms performed satisfactorily over the full range of vegetation conditions.     

Retrieving leaf area index using a genetic algorithm with a canopy radiative transfer model

Leaf area index (LAI) is an important structural property of vegetation canopy and is also one of the basic quantities driving the algorithms used in regional and global biogeochemical, ecological and meteorological applications. LAI can be estimated from remotely sensed data through the vegetation indices (VI) and the inversion of a canopy radiative transfer (RT) model. In recent years, applications of the genetic algorithms (GA) to a variety of optimization problems in remote sensing have been successfully demonstrated. In this study, we estimated LAI by integrating a canopy RT model and the GA optimization technique. This method was used to retrieve LAI from field measured reflectance as well as from atmospherically corrected Landsat ETM+ data. Four different ETM+ band combinations were tested to evaluate their effectiveness. The impacts of using the number of the genes were also examined. The results were very promising compared with field measured LAI data, and the best results were obtained with three genes in which the R² is 0.776 and the root-mean-square error (RMSE) 1.064.     

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.     

Airborne discrete-return LIDAR data in the estimation of vertical canopy cover, angular canopy closure and leaf area index

Remote sensing of forest canopy cover has been widely studied recently, but little attention has been paid to the quality of field validation data. Ecological literature has two different coverage metrics. Vertical canopy cover (VCC) is the vertical projection of tree crowns ignoring within-crown gaps. Angular canopy closure (ACC) is the proportion of covered sky at some angular range around the zenith, and can be measured with a field-of-view instrument, such as a camera. We compared field-measured VCC and ACC at 15° and 75° from the zenith to different LiDAR (Light Detection and Ranging) metrics, using several LiDAR data sets and comprehensive field data. The VCC was estimated to a high precision using a simple proportion of canopy points in first-return data. Confining to a maximum 15° scan zenith angle, the absolute root mean squared error (RMSE) was 3.7-7.0%, with an overestimation of 3.1-4.6%. We showed that grid-based methods are capable of reducing the inherent overestimation of VCC. The low scan angles and low power settings that are typically applied in topographic LiDARs are not suitable for ACC estimation as they measure in wrong geometry and cannot easily detect small within-crown gaps. However, ACC at 0-15° zenith angles could be estimated from LiDAR data with sufficient precision, using also the last returns (RMSE 8.1-11.3%, bias -6.1-+4.6%). The dependency of LiDAR metrics and ACC at 0-75° zenith angles was nonlinear and was modeled from laser pulse proportions with nonlinear regression with a best-case standard error of 4.1%. We also estimated leaf area index from the LiDAR metrics with linear regression with a standard error of 0.38. The results show that correlations between airborne laser metrics and different canopy field characteristics are very high if the field measurements are done with equivalent accuracy.     

Mapping leaf chlorophyll and leaf area index using inverse and forward canopy reflectance modeling and SPOT reflectance data

Reflectance data in the green, red and near-infrared wavelength region were acquired by the SPOT high resolution visible and geometric imaging instruments for an agricultural area in Denmark (56°N, 9°E) for the purpose of estimating leaf chlorophyll content (C_ab) and green leaf area index (LAI). SPOT reflectance observations were atmospherically corrected using aerosol data from MODIS and profiles of air temperature, humidity and ozone from the Atmospheric Infrared Sounder (AIRS), and used as input for the inversion of a canopy reflectance model. Computationally efficient inversion schemes were developed for the retrieval of soil and land cover-specific parameters which were used to build multiple species and site dependent formulations relating the two biophysical properties of interest to vegetation indices or single spectral band reflectances. Subsequently, the family of model generated relationships, each a function of soil background and canopy characteristics, was employed for a fast pixel-wise mapping of C_ab and LAI.The biophysical parameter retrieval scheme is completely automated and image-based and solves for the soil background reflectance signal, leaf mesophyll structure, specific dry matter content, Markov clumping characteristics, C_ab and LAI without utilizing calibration measurements.Despite the high vulnerability of near-infrared reflectances (ρ_nir) to variations in background properties, an efficient correction for background influences and a strong sensitivity of ρ_nir to LAI, caused LAI-ρ_nir relationships to be very useful and preferable over LAI-NDVI relationships for LAI prediction when LAI > 2. Reflectances in the green waveband (ρ_green) were chosen for producing maps of C_ab.The application of LAI-NDVI, LAI-ρ_nir and C_ab-ρ_green relationships provided reliable quantitative estimates of C_ab and LAI for agricultural crops characterized by contrasting architectures and leaf biochemical constituents with overall root mean square deviations between estimates and in-situ measurements of 0.74 for LAI and 5.0 μg cm^− 2 for C_ab.The results of this study illustrate the non-uniqueness of spectral reflectance relationships and the potential of physically-based inverse and forward canopy reflectance modeling techniques for a reasonably fast and accurate retrieval of key biophysical parameters at regional scales.     

A consistency analysis of surface reflectance and leaf area index retrieval from overlapping clear-sky Landsat ETM+ imagery

In this study, the consistency of systematic retrievals of surface reflectance and leaf area index was assessed using overlap regions in adjacent Landsat Enhanced Thematic Mapper-Plus (ETM+) scenes. Adjacent scenes were acquired within 7-25 days apart to minimize variations in the land surface reflectance between acquisition dates. Each Landsat ETM+ scene was independently geo-referenced and atmospherically corrected using a variety of standard approaches. Leaf area index (LAI) models were then applied to the surface reflectance data and the difference in LAI between overlapping scenes was evaluated. The results from this analysis show that systematic LAI retrieval from Landsat ETM+ imagery using a baseline atmospheric correction approach that assumes a constant aerosol optical depth equal to 0.06 is consistent to within ±0.61 LAI units. The average absolute difference in LAI retrieval over all 10 image pairs was 26% for a mean LAI of 2.05 and the maximum absolute difference over any one pair was 61% for a mean LAI of 1.13. When no atmospheric correction was performed on the data, the consistency in LAI retrieval was improved by 1%. When a scene-based dense, dark vegetation atmospheric correction algorithm was used, the LAI retrieval differences increased to 28% for a mean LAI of 2.32. This implies that a scene-based atmospheric correction procedure may improve the absolute accuracy of LAI retrieval without having a major impact on retrieval consistency. Such consistency trials provide insight into the current limits concerning surface reflectance and LAI retrieval from fine spatial resolution remote sensing imagery with respect to the variability in clear-sky atmospheric conditions.     

Retrieval of crop chlorophyll content and leaf area index from decompressed hyperspectral data: the effects of data compression

The objective of this study is to evaluate whether the retrieval of the leaf chlorophyll content and leaf area index (LAI) for precision agriculture application from hyperspectral data is significantly affected by data compression. This analysis was carried out using the hyperspectral data sets acquired by Compact Airborne Spectrographic Imager (CASI) over corn fields at L'Acadie experimental farm (Agriculture and Agri-Food Canada) during the summer of 2000 and over corn, soybean and wheat fields at the former Greenbelt farm (Agriculture and Agri-Food Canada) in three intensive field campaigns during the summer of 2001. Leaf chlorophyll content and LAI were retrieved from the original data and the reconstructed data compressed/decompressed by the compression algorithm called Successive approximation multi-stage vector quantization (SAMVQ) at compression ratios of 20:1, 30:1, and 50:1. The retrieved products were evaluated against the ground-truth.In the retrieval of leaf chlorophyll content (the first data set), the spatial patterns were examined in all of the images created from the original and reconstructed data and were proven to be visually unchanged, as expected. The data measures R², absolute RMSE, and relative RMSE between the leaf chlorophyll content derived from the original and reconstructed data cubes, and the laboratory-measured values were calculated as well. The results show the retrieval accuracy of crop chlorophyll content is not significantly affected by SAMVQ at the compression ratios of 20:1, 30:1, and 50:1, relative to the observed uncertainties in ground truth values. In the retrieval of LAI (the second data set), qualitative and quantitative analyses were performed. The results show that the spatial and temporal patterns of the LAI images are not significantly affected by SAMVQ and the retrieval accuracies measured by the R², absolute RMSE, and relative RMSE between the ground-measured LAI and the estimated LAI are not significantly affected by the data compression either.     

Vegetation water content mapping using Landsat data derived normalized difference water index for corn and soybeans

Information about vegetation water content (VWC) has widespread utility in agriculture, forestry, and hydrology. It is also useful in retrieving soil moisture from microwave remote sensing observations. Providing a VWC estimate allows us to control a degree of freedom in the soil moisture retrieval process. However, these must be available in a timely fashion in order to be of value to routine applications, especially soil moisture retrieval. As part of the Soil Moisture Experiments 2002 (SMEX02), the potential of using satellite spectral reflectance measurements to map and monitor VWC for corn and soybean canopies was evaluated. Landsat Thematic Mapper and Enhanced Thematic Mapper Plus data and ground-based VWC measurements were used to establish relationships based on remotely sensed indices. The two indices studied were the Normalized Difference Vegetation Index (NDVI) and the Normalized Difference Water Index (NDWI). The NDVI saturated during the study period while the NDWI continued to reflect changes in VWC. NDWI was found to be superior based upon a quantitative analysis of bias and standard error. The method developed was used to map daily VWC for the watershed over the 1-month experiment period. It was also extended to a larger regional domain. In order to develop more robust and operational methods, we need to look at how we can utilize the MODIS instruments on the Terra and Aqua platforms, which can provide daily temporal coverage.     

Optimal interpolation analysis of leaf area index using MODIS data

A simple data analysis technique for vegetation leaf area index (LAI) using Moderate Resolution Imaging Spectroradiometer (MODIS) data is presented. The objective is to generate LAI data that is appropriate for numerical weather prediction. A series of techniques and procedures which includes data quality control, time-series data smoothing, and simple data analysis is applied. The LAI analysis is an optimal combination of the MODIS observations and derived climatology, depending on their associated errors σ_o and σ_c. The “best estimate” LAI is derived from a simple three-point smoothing technique combined with a selection of maximum LAI (after data quality control) values to ensure a higher quality. The LAI climatology is a time smoothed mean value of the “best estimate” LAI during the years of 2002-2004. The observation error is obtained by comparing the MODIS observed LAI with the “best estimate” of the LAI, and the climatological error is obtained by comparing the “best estimate” of LAI with the climatological LAI value. The LAI analysis is the result of a weighting between these two errors. Demonstration of the method described in this paper is presented for the 15-km grid of Meteorological Service of Canada (MSC)'s regional version of the numerical weather prediction model. The final LAI analyses have a relatively smooth temporal evolution, which makes them more appropriate for environmental prediction than the original MODIS LAI observation data. They are also more realistic than the LAI data currently used operationally at the MSC which is based on land-cover databases.     

Modeling and assimilation of root zone soil moisture using remote sensing observations in Walnut Gulch Watershed during SMEX04

Soil moisture status in the root zone is an important component of the water cycle at all spatial scales (e.g., point, field, catchment, watershed, and region). In this study, the spatio-temporal evolution of root zone soil moisture of the Walnut Gulch Experimental Watershed (WGEW) in Arizona was investigated during the Soil Moisture Experiment 2004 (SMEX04). Root zone soil moisture was estimated via assimilation of aircraft-based remotely sensed surface soil moisture into a distributed Soil-Water-Atmosphere-Plant (SWAP) model. An ensemble square root filter (EnSRF) based on a Kalman filtering scheme was used for assimilating the aircraft-based soil moisture observations at a spatial resolution of 800 m × 800 m. The SWAP model inputs were derived from the SSURGO soil database, LAI (Leaf Area Index) data from SMEX04 database, and data from meteorological stations/rain gauges at the WGEW. Model predictions are presented in terms of temporal evolution of soil moisture probability density function at various depths across the WGEW. The assimilation of the remotely sensed surface soil moisture observations had limited influence on the profile soil moisture. More specifically, root zone soil moisture depended mostly on the soil type. Modeled soil moisture profile estimates were compared to field measurements made periodically during the experiment at the ground based soil moisture stations in the watershed. Comparisons showed that the ground-based soil moisture observations at various depths were within ± 1 standard deviation of the modeled profile soil moisture. Density plots of root zone soil moisture at various depths in the WGEW exhibited multi-modal variations due to the uneven distribution of precipitation and the heterogeneity of soil types and soil layers across the watershed.     

Analysis of leaf area index products from combination of MODIS Terra and Aqua data

A prototype product suite, containing the Terra 8-day, Aqua 8-day, Terra-Aqua combined 8- and 4-day products, was generated as part of testing for the next version (Collection 5) of the MODerate resolution Imaging Spectroradiometer (MODIS) leaf area index (LAI) products. These products were analyzed for consistency between Terra and Aqua retrievals over the following data subsets in North America: single 8-day composite over the whole continent and annual time series over three selected MODIS tiles (1200 × 1200 km). The potential for combining retrievals from the two sensors to derive improved products by reducing the impact of environmental conditions and temporal compositing period was also explored. The results suggest no significant discrepancies between large area (from continent to MODIS tile) averages of the Terra and Aqua 8-day LAI and surface reflectances products. The differences over smaller regions, however, can be large due to the random nature of residual atmospheric effects. High quality retrievals from the radiative transfer based algorithm can be expected in 90-95% of the pixels with mostly herbaceous cover and about 50-75% of the pixels with woody vegetation during the growing season. The quality of retrievals during the growing season is mostly restricted by aerosol contamination of the MODIS data. The Terra-Aqua combined 8-day product helps to minimize this effect and increases the number of high quality retrievals by 10-20% over woody vegetation. The combined 8-day product does not improve the number of high quality retrievals during the winter period because the extent of snow contamination of Terra and Aqua observations is similar. Likewise, cloud contamination in the single-sensor and combined products is also similar. The LAI magnitudes, seasonal profiles and retrieval quality in the combined 4-day product are comparable to those in the single-sensor 8-day products. Thus, the combined 4-day product doubles the temporal resolution of the seasonal cycle, which facilitates phenology monitoring in application studies during vegetation transition periods. Both Terra and Aqua LAI products show anomalous seasonality in boreal needle leaf forests, due to limitations of the radiative transfer algorithm to model seasonal variations of MODIS surface reflectance data with respect to solar zenith angle. Finally, this study suggests that further improvement of the MODIS LAI products is mainly restricted by the accuracy of the MODIS observations.     

An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest

Canopy leaf area index (LAI), defined as the single-sided leaf area per unit ground area, is a quantitative measure of canopy foliar area. LAI is a controlling biophysical property of vegetation function, and quantifying LAI is thus vital for understanding energy, carbon and water fluxes between the land surface and the atmosphere. LAI is routinely available from Earth Observation (EO) instruments such as MODIS. However EO-derived estimates of LAI require validation before they are utilised by the ecosystem modelling community. Previous validation work on the MODIS collection 4 (c4) product suggested considerable error especially in forested biomes, and as a result significant modification of the MODIS LAI algorithm has been made for the most recent collection 5 (c5). As a result of these changes the current MODIS LAI product has not been widely validated. We present a validation of the MODIS c5 LAI product over a 121 km² area of mixed coniferous forest in Oregon, USA, based on detailed ground measurements which we have upscaled using high resolution EO data. Our analysis suggests that c5 shows a much more realistic temporal LAI dynamic over c4 values for the site we examined. We find improved spatial consistency between the MODIS c5 LAI product and upscaled in situ measurements. However results also suggest that the c5 LAI product underestimates the upper range of upscaled in situ LAI measurements.     

Prototyping of Land-SAF leaf area index algorithm with VEGETATION and MODIS data over Europe

The Satellite Application Facility on Land Surface Analysis (Land-SAF) aims to provide land surface variables for the meteorological and environmental science communities from EUMETSAT satellites. This study assesses the performance of a simplified (i.e. random distribution of vegetation is assumed) version of the Land-SAF algorithm for the estimation of Leaf Area Index (LAI) when prototyped with VEGETATION (processed in CYCLOPES program) and MODIS reflectances. The prototype estimates of LAI are evaluated both by comparison with validated CYCLOPES and MODIS LAI products derived from the same sensors and directly through comparison with ground-based estimates. Emphasis is given on evaluating the impact of the algorithm and input data on LAI retrieval discrepancies. Analysis is achieved over Europe for the 2000-2003 period. The results demonstrate the capacity of the Land-SAF algorithm to retrieve consistent LAI estimates from multiple optical sensors even when their reflectances present systematic differences. High spatial and temporal consistencies between Land-SAF prototype estimates and existing LAI products are found. The differences between Land-SAF and CYCLOPES LAI are lower than their uncertainties (RMSE (relative RMSE) within 0.4 (30%)). Land-SAF prototype estimates and MODIS LAI show larger discrepancies mainly due to differences in the vegetation structure representation and algorithm assumptions (RMSE ranging from 0.2 (30%) up to 0.8 (40%)). Land-SAF prototype provides higher LAI values than MODIS for herbaceous canopies (i.e. shrubs, grasses and crops) and lower values for woody biomes (i.e. savannas and forests). Direct validation indicates that LAI estimates from prototyping of the Land-SAF algorithm with CYCLOPES and MODIS reflectances achieve similar performances (differences with ground measurements are lower than 0.5 LAI units in 60% and 50% of the cases, respectively) as CYCLOPES and MODIS LAI products. Results from this prototyping exercise appear useful for improved retrieval of LAI and constitute a step forward for refinement, validation and consolidation of the Land-SAF algorithm.     

Reflectance seasonality and its relation to the canopy leaf area index in an eastern Siberian larch forest: Multi-satellite data and radiative transfer analyses

Reliable monitoring of seasonality in the forest canopy leaf area index (LAI) in Siberian forests is required to advance the understanding of climate-forest interactions under global environmental change and to develop a forest phenology model within ecosystem modeling. Here, we compare multi-satellite (AVHRR, MODIS, and SPOT/VEGETATION) reflectance, normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and LAI with aircraft-based spectral reflectance data and field-measured forest data acquired from April to June in 2000 in a larch forest near Yakutsk, Russia. Field data in a 30 × 30-m study site and aircraft data observed around the field site were used. Larch is a dominant forest type in eastern Siberia, but comparison studies that consider multi-satellite data, aircraft-based reflectance, and field-based measurement data are rarely conducted. Three-dimensional canopy radiative transfer calculations, which are based on Antyufeev and Marshak's [Antyufeev, V.S., & Marshak, A.L. (1990). Monte Carlo method and transport equation in plant canopies, Remote Sensing of Environment, 31, 183-191] Monte Carlo photon transport method combined with North's [North, P.R. (1996). Three-dimensional forest light interaction model using a Monte Carlo method, IEEE Transactions on Geoscience and Remote Sensing, 34(4), 946-956] geometric-optical hybrid forest canopy scene, helped elucidate the relationship between canopy reflectance and forest structural parameters, including several forest floor conditions. Aircraft-based spectral measurements and the spectral response functions of all satellite sensors confirmed that biases in reflectance seasonality caused by differences in spectral response functions among sensors were small. However, some reflectance biases occur among the near infrared (NIR) reflectance data from satellite products; these biases were potentially caused by absolute calibration errors or cloud/cloud shadow contamination. In addition, reflectance seasonality in AVHRR-based NIR data was very small compared to other datasets, which was partially due to the spring-to-summer increase in the amount of atmospheric water vapor. Radiative transfer simulations suggest that bi-directional reflectance effects were small for the study site and observation period; however, changes in tree density and forest floor conditions affect the absolute value of NIR reflectance, even if the canopy leaf area condition does not change. Reliable monitoring of canopy LAI is achieved by minimizing these effects through the use of NIR reflectance difference, i.e., the difference in reflectance on the observation day from the reflectance on a snow-free/pre-foliation day. This may yield useful and robust parameters for multi-satellite monitoring of the larch canopy LAI with less error from intersensor biases and forest structure/floor differences. Further validation with field data and combined use of other index (e.g. normalized difference water index, NDWI) data will enable an extension of these findings to all Siberian deciduous forests.     

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.