Effects of soil composition and mineralogy on remote sensing of crop residue cover

The management of crop residues (non-photosynthetic vegetation) in agricultural fields influences soil erosion and soil carbon sequestration. Remote sensing methods can efficiently assess crop residue cover and related tillage intensity over many fields in a region. Although the reflectance spectra of soils and crop residues are often similar in the visible, near infrared, and the lower part of the shortwave infrared (400-1900 nm) wavelength region, specific diagnostic chemical absorption features are evident in the upper shortwave infrared (1900-2500 nm) region. Two reflectance band height indices used for estimating residue cover are the Cellulose Absorption Index (CAI) and the Lignin-Cellulose Absorption (LCA) index, both of which use reflectances in the upper shortwave infrared (SWIR). Soil mineralogy and composition will affect soil spectral properties and may limit the usefulness of these spectral indices in certain areas. Our objectives were to (1) identify minerals and soil components with absorption features in the 2000 nm to 2400 nm wavelength region that would affect CAI and LCA and (2) assess their potential impact on remote sensing estimates of crop residue cover. Most common soil minerals had CAI values ≤ 0.5, whereas crop residues were always > 0.5, allowing for good contrast between soils and residues. However, a number of common soil minerals had LCA values > 0.5, and, in some cases, the mineral LCA values were greater than those of the crop residues, which could limit the effectiveness of LCA for residue cover estimation. The LCA of some dry residues and live corn canopies were similar in value, unlike CAI. Thus, the Normalized Difference Vegetation Index (NDVI) or similar method should be used to separate out green vegetation pixels. Mineral groups, such as garnets and chlorites, often have wide ranges of CAI and LCA values, and thus, mineralogical analyses often do not identify individual mineral species required for precise CAI estimation. However, these methods are still useful for identifying mineral soils requiring additional scrutiny. Future advanced multi- and hyperspectral remote sensing platforms should include CAI bands to allow for crop residue cover estimation.     

Mitigating the effects of soil and residue water contents on remotely sensed estimates of crop residue cover

Crop residues on the soil surface decrease soil erosion and increase soil organic carbon and the management of crop residues is an integral part of many conservation tillage systems. Current methods of measuring residue cover are inadequate for characterizing the spatial variability of residue cover over large fields. The objectives of this research were to determine the effects of water content on the remotely sensed estimates of crop residue cover and to propose a method to mitigate the effects of water content on remotely sensed estimates of crop residue cover. Reflectance spectra of crop residues and soils were measured in the lab over the 400-2400 nm wavelength region. Reflectance of scenes with various residue cover fractions and water contents was simulated using a linear mixture model. Additional spectra of scenes with mixtures of crop residues and soil were also acquired in corn, soybean, and wheat fields with different tillage treatments and different water content conditions. Crop residue cover was linearly related to the cellulose absorption index (CAI), which was defined as the relative intensity of an absorption feature near 2100 nm. Water in the crop residue significantly attenuated CAI and changed the slope of the residue cover vs. CAI relationship. Without an appropriate correction, crop residue covers were underestimated as scene water content increased. Spectral vegetation water indices were poorly related to changes in the water contents of crop residues and soils. A new reflectance ratio water index that used the two bands located on the shoulders of the cellulose absorption feature to estimate scene water conditions was proposed and tested with data from corn, soybean, and wheat fields. The ratio water index was used to describe the changes in the slope of crop residue cover vs. CAI and improve the predictions of crop residue cover. These results indicate that spatial and temporal adjustments in the spectral estimates of crop residue cover are possible. Current mutispectral imaging systems will not provide reliable estimates of crop residue cover when scene water content varies. Hyperspectral data are not required, because the three narrow bands that are used for both CAI and the scene moisture correction could be incorporated in advanced multispectral sensors. Thus, regional surveys of soil conservation practices that affect soil carbon dynamics may be feasible using either advanced multispectral or hyperspectral imaging systems.     

Soil moisture estimation in a semiarid rangeland using ERS-2 and TM imagery

Soil moisture is important information in semiarid rangelands where vegetation growth is heavily dependent on the water availability. Although many studies have been conducted to estimate moisture in bare soil fields with Synthetic Aperture Radar (SAR) imagery, little success has been achieved in vegetated areas. The purpose of this study is to extract soil moisture in sparsely to moderately vegetated rangeland surfaces with ERS-2/TM synergy. We developed an approach to first reduce the surface roughness effect by using the temporal differential backscatter coefficient (Δσ_wet-dry⁰). Then an optical/microwave synergistic model was built to simulate the relationship among soil moisture, Normalized Difference Vegetation Index (NDVI) and Δσ_wet-dry⁰. With NDVI calculated from TM imagery in wet seasons and Δσ_wet-dry⁰ from ERS-2 imagery in wet and dry seasons, we derived the soil moisture maps over desert grass and shrub areas in wet seasons. The results showed that in the semiarid rangeland, radar backscatter was positively correlated to NDVI when soil was dry (m_v<10%), and negatively correlated to NDVI when soil moisture was higher (m_v>10%). The approach developed in this study is valid for sparse to moderate vegetated areas. When the vegetation density is higher (NDVI>0.45), the SAR backscatter is mainly from vegetation layer and therefore the soil moisture estimation is not possible in this study.     

Estimating fractional cover of photosynthetic vegetation, non-photosynthetic vegetation and bare soil in the Australian tropical savanna region upscaling the EO-1 Hyperion and MODIS sensors

Quantitative estimation of fractional cover of photosynthetic vegetation (f_PV), non-photosynthetic vegetation (f_NPV) and bare soil (f_BS) is critical for natural resource management and for modeling carbon dynamics. Accurate estimation of fractional cover is especially important for monitoring and modeling savanna systems, subject to highly seasonal rainfall and drought, grazing by domestic and native animals, and frequent burning. This paper describes a method for resolving f_PV, f_NPV and f_BS across the ~ 2 million km² Australian tropical savanna zone with hyperspectral and multispectral imagery. A spectral library compiled from field campaigns in 2005 and 2006, together with three EO-1 Hyperion scenes acquired during the 2005 growing season were used to explore the spectral response space for f_PV, f_NPV and f_BS. A linear unmixing approach was developed using the Normalized Difference Vegetation Index (NDVI) and the Cellulose Absorption Index (CAI). Translation of this approach to MODerate resolution Imaging Spectroradiometer (MODIS) scale was assessed by comparing multiple linear regression models of NDVI and CAI with a range of indices based on the seven MODIS bands in the visible and shortwave infrared region (SWIR) using synthesized MODIS surface reflectance data on the same dates as the Hyperion acquisitions. The best resulting model, which used NDVI and the simple ratio of MODIS bands 7 and 6 (SWIR3/SWIR2), was used to generate a time series of fractional cover from 16 day MODIS nadir bidirectional reflectance distribution function-adjusted reflectance (NBAR) data from 2000-2006. The results obtained with MODIS NBAR were validated against grass curing measurement at ten sites with good agreement at six sites, but some underestimation of f_NPV proportions at four other sites due to substantial sub-pixel heterogeneity. The model was also compared with remote sensing measurements of fire scars and showed a good matching in the spatio-temporal patterns of grass senescence and posterior burning. The fractional cover profiles for major grassland cover types showed significant differences in relative proportions of f_PV, f_NPV and f_BS, as well as large intra-annual seasonal variation in response to monsoonal rainfall gradients and soil type. The methodology proposed here can be applied to other mixed tree-grass ecosystems across the world.     

Multiple shadow fractions in spectral mixture analysis of a cotton canopy

Shadows are being used more frequently to estimate plant canopy biophysical characteristics. Typically, a zero value is assumed or a threshold value is derived from histogram analysis of imagery to determine the shadow endmember (EM). Here, two distinct shadow EMs were measured in situ for use in spectral mixture analysis of a cotton canopy on five dates in 2003. The four EMs used in the analysis were: sunlit green leaf, sunlit dry soil, self-shadowed leaf, shadowed dry soil. This 4-EM model was compared to a 3-EM model where a zero-value shade EM was used for unmixing with the two sunlit EMs. Multiple endmember spectral mixture analysis (MESMA) was used to allow EM composition to vary across each scene. The analysis and EMs were applied to fine-scale hyperspectral image data collected in the wavelength range, 440 to 810 nm. Ground data collected included percent cover, height, SPAD (a measure of leaf greenness), and chlorophyll a concentration. The normalized difference vegetation index (NDVI) was also compared to the unmixing results. Regression analysis showed that NDVI was equal to the 4-EM model for estimation of percent cover (r² = 0.95, RMSE = 6.6) although the NDVI y-intercept was closer to zero. The 4-EM model was best for estimating height (r² = 0.79, RMSE = 0.07 m) and chlorophyll a concentration (r² = 0.46, RMSE = 7.0 μg/cm²). The 3-EM model and NDVI performed poorly when estimating chlorophyll a concentration. Inclusion of two distinct shadow EMs in the model improved relationships to crop biophysical parameters and was better than assuming one, zero-value shade EM. Since MESMA operates at the pixel level and allows variable EM assignment to each pixel, mapping the spatial variability of shadows and other variables of interest is possible, providing a powerful input to canopy and ecosystem models as well as precision farming.     

A procedure for obtaining green plant cover: Relation to NDVI in a case study for barley

This paper describes a procedure for obtaining green plant cover from crop digital pictures. However, the appearance of shades caused by vegetation canopy characteristics, in addition to different coloration from several parts of the plant during the maturity process, makes evaluation difficult, especially when separating green vegetation from dry vegetation. In order to broach this problem, a methodology has been developed to obtain green plant cover. An experience in field radiometry is described where different average spectra have been captured for each experimental plot at the same time as pictures were taken, obtaining the Normalized Difference Vegetation Index (NDVI) for each one of the spectra. Finally we observe the relation between NDVI and green plant cover.     

Comparison of ERS-2 SAR and Landsat TM imagery for monitoring agricultural crop and soil conditions

Studies over the past 25 years have shown that measurements of surface reflectance and temperature (termed optical remote sensing) are useful for monitoring crop and soil conditions. Far less attention has been given to the use of radar imagery, even though synthetic aperture radar (SAR) systems have the advantages of cloud penetration, all-weather coverage, high spatial resolution, day/night acquisitions, and signal independence of the solar illumination angle. In this study, we obtained coincident optical and SAR images of an agricultural area to investigate the use of SAR imagery for farm management. The optical and SAR data were normalized to indices ranging from 0 to 1 based on the meteorological conditions and sun/sensor geometry for each date to allow temporal analysis. Using optical images to interpret the response of SAR backscatter (σ^o) to soil and plant conditions, we found that SAR σ^o was sensitive to variations in field tillage, surface soil moisture, vegetation density, and plant litter. In an investigation of the relation between SAR σ^o and soil surface roughness, the optical data were used for two purposes: (1) to filter the SAR images to eliminate fields with substantial vegetation cover and/or high surface soil moisture conditions, and (2) to evaluate the results of the investigation. For dry, bare soil fields, there was a significant correlation (r²=.67) between normalized SAR σ^o and near-infrared (NIR) reflectance, due to the sensitivity of both measurements to surface roughness. Recognizing the limitations of optical remote sensing data due to cloud interference and atmospheric attenuation, the findings of this study encourage further studies of SAR imagery for crop and soil assessment.     

Unmixing of soil types and estimation of soil exposure with simulated hyperspectral imagery

Hyperspectral determination of soil types has the potential to become an important addition to the methods used for classification and mapping of soils. In this study laboratory measured spectra of different soils, vegetation and crop residue were combined to simulate hyperspectral remote sensing imagery. The overall aim was to examine the spectral unmixing of these materials under laboratory conditions to better understand the limits to prediction of soil types and determination of cover fractions. Two different methods were utilized to mix spectra of the soil and vegetation and substantial differences were observed in the unmixing results from the different image types, particularly in mixed pixels. Results found pure soils were easily distinguished from each other when not mixed with vegetation, while some mixes of soil and vegetation were confused as pure soil spectra. The accuracy of abundance fractions retrieved in the unmixing process also varied substantially with soil type and vegetation cover.     

Dynamics of the relationship between NDVI and SWIR32 vegetation indices in southern Africa: implications for retrieval of fractional cover from MODIS data

Fractional cover of photosynthetic vegetation (F_PV), non-photosynthetic vegetation (F_NPV), and bare soil (F_BS) has been retrieved for Australian tropical savannah based on linear unmixing of the two-dimensional response envelope of the normalized difference vegetation index (NDVI) and short wave infrared ratio (SWIR)32 vegetation indices (VI) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance data. The approach assumes that cover fractions are made up of a simple mixture of green leaves, senescent leaves, and bare soil. In this study, we examine retrieval of fractional cover using this approach for a study area in southern Africa with a more complex vegetation structure. Region-specific end-members were defined using Hyperion images from different locations and times of the season. These end-members were applied to a 10-year time series of MODIS-derived NDVI and SWIR32 (from 2002 to 2011) to unmix F_PV, F_NPV, and F_BS. Results of validation with classified high-resolution imagery indicated major bias in estimation of F_NPV and F_BS, with regression coefficients for predicted versus observed data substantially less than 1.0 and relatively large intercept values. Examination with Hyperion images of the inverse relationship between the MODIS-equivalent SWIR32 index and the Hyperion-derived cellulose absorption index (CAI) to which it nominally approximates revealed: (1) non-compliant positive regression coefficients for certain vegetation types; and (2) shifts in slope and intercept of compliant regression curves related to day of year and geographical location. The results suggest that the NDVI–SWIR32 response cannot be used to approximate the NDVI–CAI response in complex savannah systems like southern Africa that cannot be described as simple mixtures of green leaves, dry herbaceous material high in cellulose, and bare soil. Methods that use a complete set of multispectral channels at higher spatial resolution may be needed for accurate retrieval of fractional cover in Africa.     

Characterizing vegetation indices derived from active and passive sensors

New ‘active’ sensors containing their own light source may provide consistent measures of plant and soil characteristics under varying illumination without calibration to reflectance. In 2006, an active sensor (Crop Circle) and various passive sensors were compared in a wheat (Triticum aestivum L., c.v. Chara) experiment in Horsham, VIC, Australia. The normalized difference vegetation index (NDVI) and soil-adjusted vegetation index (SAVI) were calculated from plot data with a range of canopy cover, leaf area and biomass. The active sensor NDVI and SAVI data were slightly less effective than corresponding passive sensor data at estimating green cover (r ²?=?0.80?0.90 vs. ～0.95). Passive sensor measurements showed strong non-linearity for estimating dry biomass and green leaf area index (GLAI), whereas SAVI calculated from the active sensor was linear (r ²?=?0.86 and 0.90). Scaling effects were not apparent when point, transect and plot areas were compared at the given level of spatial variation. Sensor height above the target confounded SAVI data probably due to differential irradiance from the light sources and the unbalanced effect of the ‘L’ factor within the algorithm. The active sensor was insensitive to light conditions (r ²?=?0.99 for cloudy vs. clear skies) and had no requirement for optical calibration.     

The impact of soil reflectance on the quantification of the green vegetation fraction from NDVI

The green vegetation fraction (Fg) is an important climate and hydrologic model parameter. A common method to calculate Fg is to create a simple linear mixing model between two NDVI endmembers: bare soil NDVI (NDVI_o) and full vegetation NDVI (NDVI_∞). Usually it is assumed that NDVI_o is close to zero (NDVI_o ∼ 0.05) and is generally chosen from the lowest observed NDVI values. However, the mean soil NDVI computed from 2906 samples is much larger (NDVI = 0.2) and is highly variable (standard deviation = 0.1). We show that the underestimation of NDVI_o yields overestimations of Fg. The largest errors occur in grassland and shrubland areas. Using parameters for NDVI_o and NDVI_∞ derived from global scenes yields overestimations of Fg that are larger than 0.2 for the majority of U.S. land cover types when pixel NDVI values are 0.2 < NDVI_pixel < 0.4. When using conterminous U.S. scenes to derive NDVI_o and NDVI_∞, the overestimation is less (0.10-0.17 for 0.2 < NDVI_pixel < 0.4). As a result, parts of the conterminous U.S. are affected at different times of the year depending on the local seasonal NDVI cycle. We propose using global databases of NDVI_o along with information on historical NDVI_pixel values to compute a statistically most-likely estimate of Fg. Using in situ measurements made at the Sevilleta LTER, we show that this approach yields better estimates of Fg than using global invariant NDVI_o values estimated from whole scenes. At the two studied sites, the Fg estimate was adjusted by 52% at the grassland and 86% at the shrubland. More significant advances will require information on spatial distribution of soil reflectance.     

Mapping of Sahelian vegetation parameters from ERS scatterometer data with an evolution strategies algorithm

The West African Sahel rainfall regime is known for its spatio-temporal variability at different scales which has a strong impact on vegetation development. This study presents results of the combined use of a simple water balance model, a radiative transfer model and ERS scatterometer data to produce map of vegetation biomass and thus vegetation cover at a spatial resolution of 25 km. The backscattering coefficient measured by spaceborne wind scatterometers over Sahel shows a marked seasonality linked to the drastic changes of both soil and vegetation dielectric properties associated to the alternating dry and wet seasons. For lack of a direct observation, METEOSAT rainfall estimates are used to calculate temporal series of soil moisture with the help of a water balance model. This a priori information is used as input of the radiative transfer model that simulates the interaction between the radar wave and the surface components (soil and vegetation). Then, an inversion algorithm is applied to retrieve vegetation aerial mass from the ERS scatterometer data. Because of the nonlinear feature of the inverse problem to be solved, the inversion is performed using a global stochastic nonlinear inversion method. A good agreement is obtained between the inverse solutions and independent field measurements with mean and standard deviation of −54 and 130 kg of dry matter by hectare (kg DM/ha), respectively. The algorithm is then applied to a 350,000 km² area including the Malian Gourma and Seno region and a Sahelian part of Burkina Faso during two contrasted seasons (1999 and 2000). At the considered resolution, the obtained herbaceous mass maps show a global qualitative consistency (r²=0.71) with NDVI images acquired by the VEGETATION instrument.     

Relationships between vegetation indices,fractional cover retrievals and the structure and composition of Brazilian Cerrado natural vegetation

This study explores the use of the relationship between the normalized difference vegetation index (NDVI) and the shortwave infrared ratio (SWIR32) vegetation indices (VI) to retrieve fractional cover over the structurally complex natural vegetation of the Cerrado of Brazil using a time series of imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS). Data from the EO-1 Hyperion sensor with 30 m pixel resolution is used to sample geographic and seasonal variation in NDVI, SWIR32, and the hyperspectral cellulose absorption index (CAI), and to derive end-member values for photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV), and bare soil (BS) from a suite of protected and/or natural vegetation sites across the Cerrado. The end-members derived from relatively pure 30 m pixels are then applied to a 500 m pixel resolution MODIS time series using linear spectral unmixing to retrieve PV, NPV, and BS fractional cover (F_PV, F_NPV, and F_BS). The two-way interaction response of MODIS-equivalent NDVI and SWIR32 was examined for regions of interest (ROI) collected within protected areas and nearby converted lands. The MODIS NDVI, SWIR32 and retrieved F_PV, F_NPV, and F_BS are then compared to detailed cover and structural composition data from field sites, and the influence of the structural and compositional variation on the VIs and cover fractions is explored. The hyperion ROI analysis indicated that the two-way NDVI–SWIR32 response behaved as an effective surrogate for the two-way NDVI–CAI response for the campo limpo/grazed pasture to cerrado sensu stricto woody gradient. The SWIR32 sensitivity to the NPV and BS variation increased as the dry season progressed, but Cerrado savannah exhibited limited dynamic range in the NDVI–CAI and NDVI–SWIR32 two-way responses compared to the entire landscape, which also comprises fallow croplands and forests. Validation analysis of MODIS retrievals with Quickbird-2 images produced an RMSE value of 0.13 for F_PV. However, the RMSE values of 0.16 and 0.18 for F_BS and F_NPV, respectively, were large relative to the seasonal and inter-annual variation. Analysis of site composition and structural data in relation to the MODIS-derived NDVI, SWIR32 and F_PV, F_NPV, and F_BS, indicated that the VI signal and derived cover fractions were influenced by a complex mix of structure and cover but included a strong year-to-year seasonal effect. Therefore, although the MODIS NDVI–SWIR32 response could be used to retrieve cover fractions across all Cerrado land covers including bare cropland, pastures and forests, sensitivity may be limited within the natural Cerrado due to sub-pixel heterogeneity and limited BS and NPV sensitivity.     

Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data

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

土壤背景对冠层NDVI的影响分析

归一化差值植被指数NDVI是植被遥感中应用最为广泛的指数之一, 但它受土壤背景等因素的干扰比较强烈。结合实测的土壤数据以及公式推导、PROSAIL 模型模拟等方法分析了这种影响。首先, 假定与土壤线性混合且叶片呈水平分布的植被冠层, 根据土壤与植被分别在红光、近红外波段处的反射率值、植被覆盖度等参数, 利用公式推导了土壤背景对不同覆盖度下冠层NDVI的影响。其次, 利用PROSAIL冠层光谱模拟模型, 模拟分析了土壤背景对不同LAI下冠层NDVI的影响。分析的结果表明:LAI 越小, 土壤背景的影响越大; 暗土壤背景下的冠层NDVI值大于亮土壤背景下冠层的NDVI值; 并且,暗土壤条件下,NDVI值对土壤亮度的变化更敏感,而亮土壤下,NDVI值则对LAI或覆盖度的变化更敏感。最后利用实测的不同土壤背景下的冬小麦冠层光谱数据, 验证了公式推导和模型模拟的结果。     

Nondestructive estimation of canopy chlorophyll content using Hyperion and Landsat/TM images

Remote sensing offers a nondestructive tool for the quick and precise estimation of canopy chlorophyll content that serves as an important indicator of the plant ecosystem. In this study, the canopy chlorophyll content of 26 samples in 2007 and 40 samples in 2008 of maize were nondestructively estimated by a set of vegetation indices (VIs; Normalized Difference Vegetation Index, NDVI; Green Chlorophyll Index, CI_green; modified soil adjust vegetation index, MSAVI; and Enhanced Vegetation Index, EVI) derived from the hyperspectral Hyperion and Thematic Mapper (TM) images. The PROSPECT model was used for sensitivity analysis among the indices and results indicated that CI_green had a large linear correlation with chlorophyll content ranging from 100–1000 mg m^?2. EVI showed a moderate ability in avoiding saturation and reached a saturation of chlorophyll content above 600 mg m^?2. Both of the other two indices, MSAVI and NDVI, showed a clear saturation at chlorophyll content of 400 mg m^?2, which demonstrated they may be inappropriate for chlorophyll interpretation at high values. A validation study was also conducted with satellite observations (Hyperion and TM) and in-situ measurements of chlorophyll content in maize. Results indicated that canopy chlorophyll content can be remotely evaluated by VIs with r ² ranging from the lowest of 0.73 for NDVI to the highest of 0.86 for CI_green. EVI had a greater precision (r ²=0.81) than MASVI (r ²=0.75) in canopy chlorophyll content estimation. The results agreed well with the sensitivity study and will be helpful in developing future models for canopy chlorophyll evaluation.     

Using lidar and effective LAI data to evaluate IKONOS and Landsat 7 ETM+ vegetation cover estimates in a ponderosa pine forest

Structural and functional analyses of ecosystems benefit when high accuracy vegetation coverages can be derived over large areas. In this study, we utilize IKONOS, Landsat 7 ETM+, and airborne scanning light detection and ranging (lidar) to quantify coniferous forest and understory grass coverages in a ponderosa pine (Pinus ponderosa) dominated ecosystem in the Black Hills of South Dakota. Linear spectral mixture analyses of IKONOS and ETM+ data were used to isolate spectral endmembers (bare soil, understory grass, and tree/shade) and calculate their subpixel fractional coverages. We then compared these endmember cover estimates to similar cover estimates derived from lidar data and field measures. The IKONOS-derived tree/shade fraction was significantly correlated with the field-measured canopy effective leaf area index (LAI_e) (r²=0.55, p<0.001) and with the lidar-derived estimate of tree occurrence (r²=0.79, p<0.001). The enhanced vegetation index (EVI) calculated from IKONOS imagery showed a negative correlation with the field measured tree canopy effective LAI and lidar tree cover response (r²=0.30, r=−0.55 and r²=0.41, r=−0.64, respectively; p<0.001) and further analyses indicate a strong linear relationship between EVI and the IKONOS-derived grass fraction (r²=0.99, p<0.001). We also found that using EVI resulted in better agreement with the subpixel vegetation fractions in this ecosystem than using normalized difference of vegetation index (NDVI). Coarsening the IKONOS data to 30 m resolution imagery revealed a stronger relationship with lidar tree measures (r²=0.77, p<0.001) than at 4 m resolution (r²=0.58, p<0.001). Unmixed tree/shade fractions derived from 30 m resolution ETM+ imagery also showed a significant correlation with the lidar data (r²=0.66, p<0.001). These results demonstrate the power of using high resolution lidar data to validate spectral unmixing results of satellite imagery, and indicate that IKONOS data and Landsat 7 ETM+ data both can serve to make the important distinction between tree/shade coverage and exposed understory grass coverage during peak summertime greenness in a ponderosa pine forest ecosystem.     

Soil moisture retrieval using thermal inertia, determined with visible and thermal spaceborne data, validated for European forests

Variations in soil moisture strongly affect surface energy balances, regional runoff, land erosion and vegetation productivity (potential crop yield). Hence, the detection of soil moisture content (SMC) is very valuable in the social, economic, humanitarian (food security) and environmental segments of society. A method to estimate SMC from optical and thermal spectral information of METEOSAT imagery based on thermal inertia (TI) is presented. Minimum and maximum TI values from time series are combined in the Soil Moisture Saturation Index (SMSI). To convert surface to soil profile values, a Markov type filter is used, based on a simple two layer water balance equation (the surface layer and the reservoir below) and an autocorrelation function. Ten-daily SMC values are compared with up-scaled (using AVHRR/NDVI) observations on 10 EUROFLUX sites in Europe for the 1997 growing season (March-October). Moreover, the thermal inertia approach is compared for 1997, with ERS Scatterometer data for eight EUROFLUX sites. METEOSAT pixels are up-scaled to accommodate the ERS Scatterometer spatial resolution. The regression coefficients (slope, intercept and R²) of the thermal inertia approach versus the up-scaled soil moisture observations from EUROFLUX sites vary between 0.811-1.148, − 0.0029-0.66 and 0.544-0.877, respectively, with a RRMSE range of 3.9% to 35.7%. The regression coefficients of the comparison of ERS Scatterometer derived Soil Water Index (SWI) versus the up-scaled Soil Moisture Saturation Index for the pooled case (binning the eight EUROFLUX sites) are 0.587, 0.105 and 0.441, respectively, with a RRMSE of 38%. A simple error propagation model applied for the thermal inertia approach reveals that the absolute and relative errors of the obtained soil moisture content is at least 0.010 m³ m^− 3 or 2.0% with a SMC of 0.203 m³ m^− 3. Recommendations are made to test and implement the TI methodology using NOAA/AVHRR imagery.     

Monitoring herbaceous fuel moisture content with SPOT VEGETATION time-series for fire risk prediction in savanna ecosystems

This paper evaluated the capacity of SPOT VEGETATION time-series to monitor herbaceous fuel moisture content (FMC) in order to improve fire risk assessment in the savanna ecosystem of Kruger National Park in South Africa. In situ herbaceous FMC data were used to assess the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Vegetation Dryness Index (VDI), Improved VDI (IVDI), and Accumulated Relative NDVI Decrement (ARND) during the dry season. The effect of increasing amounts of dead vegetation on the monitoring capacity of derived indices was studied by sampling mixed live and dead FMC. The IVDI was proposed as an improvement of the VDI to monitor herbaceous FMC during the dry season. The IVDI is derived by replacing NDVI with the integrated Relative Vegetation Index (iRVI), as an approximation of yearly herbaceous biomass, when analyzing the 2-dimensional space with NDWI. It was shown that the iRVI offered more information than the NDVI in combination with NDWI to monitor FMC. The VDI and IVDI exhibited a significant relation to FMC with R² of 0.25 and 0.73, respectively. The NDWI, however, correlated best with FMC (R² = 0.75), while the correlation of ARND and FMC was weaker (R² = 0.60) than that found for NDVI, NDWI, and IVDI. The use of in situ herbaceous FMC consequently indicated that NDWI is appropriate as spatio-temporal information source of herbaceous FMC variation which can be used to optimize fire risk and behavior assessment for fire management in savanna ecosystems.     

Neural network and crop residue index multiband models for estimating crop residue cover from Landsat TM and ETM+ images

Crop residues on the soil surface provide not only a barrier against water and wind erosion, but they also contribute to improving soil organic matter content, infiltration, evaporation, temperature, and soil structure, among others. In Argentina, soybean (Glycine max (L.) Merill) and corn (Zea mays L.) are the most important crops. The objective of this work was to develop and evaluate two different types of model for estimating soybean and corn residue cover: neural networks (NN) and crop residue index multiband (CRIM) index, from Landsat images. Data of crop residue were acquired throughout the summer growing season in the central plains of Córdoba (Argentina) and used for training and validating the models. The CRIM, a linear mixing model of composite soil and residue, and the NN design, included reflectance and digital numbers from a combination of different TM bands to estimate the fractional residue cover. The results show that both methodologies are appropriate for estimating the residue cover from Landsat data. The best developed NN model yielded R² = 0.95 when estimating soybean and corn residue cover fraction, whereas the best fit using CRIM yielded R² = 0.87; in addition, this index is dependent on the soil and residue lines considered.