A new data fusion model for high spatial- and temporal-resolution mapping of forest disturbance based on Landsat and MODIS

Investigating the temporal and spatial pattern of landscape disturbances is an important requirement for modeling ecosystem characteristics, including understanding changes in the terrestrial carbon cycle or mapping the quality and abundance of wildlife habitats. Data from the Landsat series of satellites have been successfully applied to map a range of biophysical vegetation parameters at a 30 m spatial resolution; the Landsat 16 day revisit cycle, however, which is often extended due to cloud cover, can be a major obstacle for monitoring short term disturbances and changes in vegetation characteristics through time.The development of data fusion techniques has helped to improve the temporal resolution of fine spatial resolution data by blending observations from sensors with differing spatial and temporal characteristics. This study introduces a new data fusion model for producing synthetic imagery and the detection of changes termed Spatial Temporal Adaptive Algorithm for mapping Reflectance Change (STAARCH). The algorithm is designed to detect changes in reflectance, denoting disturbance, using Tasseled Cap transformations of both Landsat TM/ETM and MODIS reflectance data. The algorithm has been tested over a 185 × 185 km study area in west-central Alberta, Canada. Results show that STAARCH was able to identify spatial and temporal changes in the landscape with a high level of detail. The spatial accuracy of the disturbed area was 93% when compared to the validation data set, while temporal changes in the landscape were correctly estimated for 87% to 89% of instances for the total disturbed area. The change sequence derived from STAARCH was also used to produce synthetic Landsat images for the study period for each available date of MODIS imagery. Comparison to existing Landsat observations showed that the change sequence derived from STAARCH helped to improve the prediction results when compared to previously published data fusion techniques.     

Possibility investigation of tree diversity mapping using Landsat ETM+ data in the Hyrcanian forests of Iran

Lack of data often limits understanding and management of biodiversity in forested areas. Remote sensing imagery has considerable potential to aid in the monitoring and prediction of biodiversity across many spatial and temporal scales. In this paper, we explored the possibility of defining relationships between species diversity indices and Landsat ETM+ reflectance values for Hyrcanian forests in Golestan province of Iran. We used the COST model for atmospheric correction of the imagery. Linear regression models were implemented to predict measures of biodiversity (species richness and reciprocal of Simpson indices) using various combinations of Landsat spectral data. Species richness was modeled using the band set ETM5, ETM7, DVI, wetness and variances of ETM1, ETM2 and ETM5 (adjusted R² = 0.59, RMSE = 1.51). Reciprocal of Simpson index was modeled using the band set NDVI, brightness, greenness, variances of ETM2, ETM5 and ETM7 (adjusted R² = 0.459 RMSE = 1.15). The results demonstrated that spectral reflectance from Landsat can be used to effectively model tree species diversity. Predictive map derived from the presented methodology can help evaluate spatial aspects and monitor tree species diversity of the studied forest. The methodology also facilitates the evaluation of forest management and conservation strategies in northern Iran.     

A vegetation index based technique for spatial sharpening of thermal imagery

High spatial resolution (∼ 100 m) thermal infrared band imagery has utility in a variety of applications in environmental monitoring. However, currently such data have limited availability and only at low temporal resolution, while coarser resolution thermal data (∼ 1000 m) are routinely available, but not as useful for identifying environmental features for many landscapes. An algorithm for sharpening thermal imagery (TsHARP) to higher resolutions typically associated with the shorter wavebands (visible and near-infrared) used to compute vegetation indices is examined over an extensive corn/soybean production area in central Iowa during a period of rapid crop growth. This algorithm is based on the assumption that a unique relationship between radiometric surface temperature (T_R) relationship and vegetation index (VI) exists at multiple resolutions. Four different methods for defining a VI − T_R basis function for sharpening were examined, and an optimal form involving a transformation to fractional vegetation cover was identified. The accuracy of the high-resolution temperature retrieval was evaluated using aircraft and Landsat thermal imagery, aggregated to simulate native and target resolutions associated with Landsat, MODIS, and GOES short- and longwave datasets. Applying TsHARP to simulated MODIS thermal maps at 1-km resolution and sharpening down to ∼ 250 m (MODIS VI resolution) yielded root-mean-square errors (RMSE) of 0.67-1.35 °C compared to the ‘observed’ temperature fields, directly aggregated to 250 m. Sharpening simulated Landsat thermal maps (60 and 120 m) to Landsat VI resolution (30 m) yielded errors of 1.8-2.4 °C, while sharpening simulated GOES thermal maps from 5 km to 1 km and 250 m yielded RMSEs of 0.98 and 1.97, respectively. These results demonstrate the potential for improving the spatial resolution of thermal-band satellite imagery over this type of rainfed agricultural region. By combining GOES thermal data with shortwave VI data from polar orbiters, thermal imagery with 250-m spatial resolution and 15-min temporal resolution can be generated with reasonable accuracy. Further research is required to examine the performance of TsHARP over regions with different climatic and land-use characteristics at local and regional scales.     

A comparison of multi-spectral, multi-angular, and multi-temporal remote sensing datasets for fractional shrub canopy mapping in Arctic Alaska

Shrub cover appears to be increasing across many areas of the Arctic tundra biome, and increasing shrub cover in the Arctic has the potential to significantly impact global carbon budgets and the global climate system. For most of the Arctic, however, there is no existing baseline inventory of shrub canopy cover, as existing maps of Arctic vegetation provide little information about the density of shrub cover at a moderate spatial resolution across the region. Remotely-sensed fractional shrub canopy maps can provide this necessary baseline inventory of shrub cover. In this study, we compare the accuracy of fractional shrub canopy (> 0.5 m tall) maps derived from multi-spectral, multi-angular, and multi-temporal datasets from Landsat imagery at 30 m spatial resolution, Moderate Resolution Imaging SpectroRadiometer (MODIS) imagery at 250 m and 500 m spatial resolution, and MultiAngle Imaging Spectroradiometer (MISR) imagery at 275 m spatial resolution for a 1067 km² study area in Arctic Alaska. The study area is centered at 69 °N, ranges in elevation from 130 to 770 m, is composed primarily of rolling topography with gentle slopes less than 10°, and is free of glaciers and perennial snow cover. Shrubs > 0.5 m in height cover 2.9% of the study area and are primarily confined to patches associated with specific landscape features. Reference fractional shrub canopy is determined from in situ shrub canopy measurements and a high spatial resolution IKONOS image swath. Regression tree models are constructed to estimate fractional canopy cover at 250 m using different combinations of input data from Landsat, MODIS, and MISR. Results indicate that multi-spectral data provide substantially more accurate estimates of fractional shrub canopy cover than multi-angular or multi-temporal data. Higher spatial resolution datasets also provide more accurate estimates of fractional shrub canopy cover (aggregated to moderate spatial resolutions) than lower spatial resolution datasets, an expected result for a study area where most shrub cover is concentrated in narrow patches associated with rivers, drainages, and slopes. Including the middle infrared bands available from Landsat and MODIS in the regression tree models (in addition to the four standard visible and near-infrared spectral bands) typically results in a slight boost in accuracy. Including the multi-angular red band data available from MISR in the regression tree models, however, typically boosts accuracy more substantially, resulting in moderate resolution fractional shrub canopy estimates approaching the accuracy of estimates derived from the much higher spatial resolution Landsat sensor. Given the poor availability of snow and cloud-free Landsat scenes in many areas of the Arctic and the promising results demonstrated here by the MISR sensor, MISR may be the best choice for large area fractional shrub canopy mapping in the Alaskan Arctic for the period 2000-2009.     

Cross-scalar satellite phenology from ground, Landsat, and MODIS data

Phenological records constructed from global mapping satellite platforms (e.g. AVHRR and MODIS) hold the potential to be valuable tools for monitoring vegetation response to global climate change. However, most satellite phenology products are not validated, and field checking coarse scale (≥ 500 m) data with confidence is a difficult endeavor. In this research, we compare phenology from Landsat (field scale, 30 m) to MODIS (500 m), and compare datasets derived from each instrument. Landsat and MODIS yield similar estimates of the start of greenness (r² = 0.60), although we find that a high degree of spatial phenological variability within coarser-scale MODIS pixels may be the cause of the remaining uncertainty. In addition, spatial variability is smoothed in MODIS, a potential source of error when comparing in situ or climate data to satellite phenology. We show that our method for deriving phenology from satellite data generates spatially coherent interannual phenology departures in MODIS data. We test these estimates from 2000 to 2005 against long-term records from Harvard Forest (Massachusetts) and Hubbard Brook (New Hampshire) Experimental Forests. MODIS successfully predicts 86% of the variance at Harvard forest and 70% of the variance at Hubbard Brook; the more extreme topography of the later is inferred to be a significant source of error. In both analyses, the satellite estimate is significantly dampened from the ground-based observations, suggesting systematic error (slopes of 0.56 and 0.63, respectively). The satellite data effectively estimates interannual phenology at two relatively simple deciduous forest sites and is internally consistent, even with changing spatial scale. We propose that continued analyses of interannual phenology will be an effective tool for monitoring native forest responses to global-scale climate variability.     

Effects of spatial and spectral resolution in estimating ecosystem α-diversity by satellite imagery

Remote sensing represents a powerful tool to derive quantitative and qualitative information about ecosystem biodiversity. In particular, since plant species richness is a fundamental indicator of biodiversity at the community and regional scales, attempts were made to predict species richness (spatial heterogeneity) by means of spectral heterogeneity. The possibility of using spectral variance of satellite images for predicting species richness is known as Spectral Variation Hypothesis. However, when using remotely sensed data, researchers are limited to specific scales of investigation. This paper aims to investigate the effects of scale (both as spatial and spectral resolution) when searching for a relation between spectral and spatial (related to plant species richness) heterogeneity, by using satellite data with different spatial and spectral resolution. Species composition was sampled within square plots of 100 m² nested in macroplots of 10,000 m². Spectral heterogeneity of each macroplot was calculated using satellite images with different spatial and spectral resolution: a Quickbird multispectral image (4 bands, spatial resolution of 3 m), an Aster multispectral image (first 9 bands used, spatial resolution of 15 m for bands 1 to 3 and 30 m for bands 4 to 9), an ortho-Landsat ETM+ multispectral image (bands 1 to 5 and band 7 used; spatial resolution, 30 m), a resampled 60 m Landsat ETM+ image.Quickbird image heterogeneity showed a statistically highly significant correlation with species richness (r = 0.69) while coarse resolution images showed contrasting results (r = 0.43, r = 0.67, and r = 0.69 considering the Aster, Landsat ETM+, and the resampled 60 m Landsat ETM+ images respectively). It should be stressed that spectral variability is scene and sensor dependent. Considering coarser spatial resolution images, in such a case even using SWIR Aster bands (i.e. the additional spectral information with respect to Quickbird image) such an image showed a very low power in catching spectral and thus spatial variability with respect to Landsat ETM+ imagery. Obviously coarser resolution data tend to have mixed pixel problems and hence less sensitive to spatial complexity. Thus, one might argue that using a finer pixel dimension should inevitably result in a higher level of detail. On the other hand, the spectral response from different land-cover features (and thus different species) in images with higher spectral resolution should exhibit higher complexity.Spectral Variation Hypothesis could be a basis for improving sampling designs and strategies for species inventory fieldwork. However, researchers must be aware on scale effects when measuring spectral (and thus spatial) heterogeneity and relating it to field data, hence considering the concept of scale not only related to a spatial framework but even to a spectral one.     

MODIS tree cover validation for the circumpolar taiga-tundra transition zone

A validation of the 2005 500 m MODIS vegetation continuous fields (VCF) tree cover product in the circumpolar taiga-tundra ecotone was performed using high resolution Quickbird imagery. Assessing the VCF's performance near the northern limits of the boreal forest can help quantify the accuracy of the product within this vegetation transition area. The circumpolar region was divided into 7 longitudinal zones and validation sites were selected in areas of varying tree cover where Quickbird imagery is available in Google Earth. Each site was linked to the corresponding VCF pixel and overlaid with a regular dot grid within the VCF pixel's boundary to estimate percent tree crown cover in the area. Percent tree crown cover was estimated using Quickbird imagery for 396 sites throughout the circumpolar region and related to the VCF's estimates of canopy cover for 2000-2005. Regression results of VCF inter-annual comparisons (2000-2005) and VCF-Quickbird image-interpreted estimates indicate that: (1) Pixel-level, inter-annual comparisons of VCF estimates of percent canopy cover were linearly related (mean R² = 0.77) and exhibited an average root mean square error (RMSE) of 10.1% and an average root mean square difference (RMSD) of 7.3%. (2) A comparison of image-interpreted percent tree crown cover estimates based on dot counts on Quickbird color images by two different interpreters were more variable (R² = 0.73, RMSE = 14.8%, RMSD = 18.7%) than VCF inter-annual comparisons. (3) Across the circumpolar boreal region, 2005 VCF-Quickbird comparisons were linearly related, with an R² = 0.57, a RMSE = 13.4% and a RMSD = 21.3%, with a tendency to over-estimate areas of low percent tree cover and anomalous VCF results in Scandinavia. The relationship of the VCF estimates and ground reference indicate to potential users that the VCF's tree cover values for individual pixels, particularly those below 20% tree cover, may not be precise enough to monitor 500 m pixel-level tree cover in the taiga-tundra transition zone.     

View angle effects on the discrimination of soybean varieties and on the relationships between vegetation indices and yield using off-nadir Hyperion data

Imaging spectrometry has the potential to provide improved discrimination of crop types and better estimates of crop yield. Here we investigate the potential of Hyperion to discriminate three Brazilian soybean varieties and to evaluate the relationship between grain yield and 17 narrow-band vegetation indices. Hyperion analysis focused on two datasets acquired from opposite off-nadir viewing directions but similar solar geometry: one acquired on 08 February 2005 (forward scattering) and the other on 14 January 2006 (back scattering). In 2005, the soybean canopies were observed by Hyperion at later reproductive stages than in 2006. Additional Hyperion datasets were not available due to cloud cover. To further examine the impact of viewing geometry within the same season, Hyperion data were complemented by 250 m Moderate Resolution Imaging Spectroradiometer (MODIS) images (bands 1 and 2) acquired in consecutive days (05-06 February 2005) with opposite viewing geometries (− 42° and + 44°, respectively). MODIS data analysis was used to keep reproductive stage as a constant factor while isolating the impact of viewing geometry. For discrimination purposes, multiple discriminant analysis (MDA) was applied over each dataset using surface reflectance values as input variables and a stepwise procedure for band selection. All possible Hyperion band ratios and the 17 narrow-band vegetation indices with soybean grain yield were evaluated across years through Pearson's correlation coefficients and linear regression. MODIS-derived Normalized Difference Vegetation Index (NDVI) and Simple Ratio (SR) were evaluated within the same growing season. Results showed that: (1) the three soybean varieties were discriminated with highest accuracy in the back scattering direction, as deduced from MDA classification results from Hyperion and MODIS data; (2) the highest correlation between Hyperion vegetation indices and soybean yield was observed for the Normalized Difference Water Index (NDWI) (r = + 0.74) in the back scattering direction and this result was consistent with band ratio analysis; (3) higher Hyperion correlation results were observed in the back scattering direction when compared to the forward scattering image. For the same reproductive stage, stronger shadowing effects were observed over the MODIS red band in the forward scattering direction producing lower and lesser variable reflectance for the sensor. As a result, the relationship between MODIS-derived NDVI and soybean yield improved from the forward (r of + 0.21) to the back scattering view (r of + 0.60). The same trend was observed for SR that increased from + 0.22 to + 0.58.     

Assessment of QuickBird high spatial resolution imagery to detect red attack damage due to mountain pine beetle infestation

High spatial resolution remotely sensed data has the potential to complement existing forest health programs for both strategic planning over large areas, as well as for detailed and precise identification of tree crowns subject to stress and infestation. The area impacted by the current mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in British Columbia, Canada, has increased 40-fold over the previous 5 years, with approximately 8.5 million ha of forest infested in 2005. As a result of the spatial extent and intensity of the outbreak, new technologies are being assessed to help detect, map, and monitor the damage caused by the beetle, and to inform mitigation of future beetle outbreaks. In this paper, we evaluate the capacity of high spatial resolution QuickBird multi-spectral imagery to detect mountain pine beetle red attack damage. ANOVA testing of individual spectral bands, as well as the Normalized Difference Vegetation Index (NDVI) and a ratio of red to green reflectance (Red-Green Index or RGI), indicated that the RGI was the most successful (p < 0.001) at separating non-attack crowns from red attack crowns. Based on this result, the RGI was subsequently used to develop a binary classification of red attack and non-attack pixels. The total number of QuickBird pixels classified as having red attack damage within a 50 m buffer of a known forest health survey point were compared to the number of red attack trees recorded at the time of the forest health survey. The relationship between the number of red attack pixels and observed red attack crowns was assessed using independent validation data and was found to be significant (r² = 0.48, p < 0.001, standard error = 2.8 crowns). A comparison of the number of QuickBird pixels classified as red attack, and a broader scale index of mountain pine beetle red attack damage (Enhanced Wetness Difference Index, calculated from a time series of Landsat imagery), was significant (r² = 0.61, p < 0.001, standard error = 1.3 crowns). These results suggest that high spatial resolution imagery, in particular QuickBird satellite imagery, has a valuable role to play in identifying tree crowns with red attack damage. This information could subsequently be used to augment existing detailed forest health surveys, calibrate synoptic estimates of red attack damage generated from overview surveys and/or coarse scale remotely sensed data, and facilitate the generation of value-added information products, such as estimates of timber volume impacts at the forest stand level.     

Remote sensing of woody shrub cover in desert grasslands using MISR with a geometric-optical canopy reflectance model

A new method is described for the retrieval of fractional cover of large woody plants (shrubs) at the landscape scale using moderate resolution multi-angle remote sensing data from the Multiangle Imaging SpectroRadiometer (MISR) and a hybrid geometric-optical (GO) canopy reflectance model. Remote sensing from space is the only feasible method for regularly mapping woody shrub cover over large areas, an important application because extensive woody shrub encroachment into former grasslands has been seen in arid and semi-arid grasslands around the world during the last 150 years. The major difficulty in applying GO models in desert grasslands is the spatially dynamic nature of the combined soil and understory background reflectance: the background is important and cannot be modeled as either a Lambertian scatterer or by using a fixed bidirectional reflectance distribution function (BRDF). Candidate predictors of the background BRDF at the Sun-target-MISR angular sampling configurations included the volume scattering kernel weight from a Li-Ross BRDF model; diffuse brightness (ρ0) from the Modified Rahman-Pinty-Verstraete (MRPV) BRDF model; other Li-Ross kernel weights (isotropic, geometric); and MISR near-nadir bidirectional reflectance factors (BRFs) in the blue, green, and near infra-red bands. The best method was multiple regression on the weights of a kernel-driven model and MISR nadir camera blue, green, and near infra-red bidirectional reflectance factors. The results of forward modeling BRFs for a 5.25 km² area in the USDA, ARS Jornada Experimental Range using the Simple Geometric Model (SGM) with this background showed good agreement with the MISR data in both shape and magnitude, with only minor spatial discrepancies. The simulations were shown to be accurate in terms of both absolute value and reflectance anisotropy over all 9 MISR views and for a wide range of canopy configurations (r² = 0.78, RMSE = 0.013, N = 3969). Inversion of the SGM allowed estimation of fractional shrub cover with a root mean square error (RMSE) of 0.03 but a relatively weak correlation (r² = 0.19) with the reference data (shrub cover estimated from high resolution IKONOS panchromatic imagery). The map of retrieved fractional shrub cover was an approximate spatial match to the reference map. Deviations reflect the first-order approximation of the understory BRDF in the MISR viewing plane; errors in the shrub statistics; and the 12 month lag between the two data sets.     

AERONET-based surface reflectance validation network (ASRVN) data evaluation: Case study for railroad valley calibration site

The AERONET-based Surface Reflectance Validation Network (ASRVN) is an operational processing system developed for validation of satellite derived surface reflectance products at regional and global scales. The ASRVN receives 50 × 50 km² subsets of MODIS data centered at AERONET sites along with AERONET aerosol and water vapor data, and performs an atmospheric correction. The ASRVN produces surface bidirectional reflectance factor (BRF), albedo, parameters of the Ross-Thick Li-Sparse (RTLS) BRF model, as well as Hemispherical-Directional Reflectance Factor (HDRF), which is required for comparison with the ground-based measurements. This paper presents a comparison of ASRVN HDRF with the ground-based HDRF measurements collected during 2001-2008 over a bright calibration Railroad Valley, Nevada site as part of the MODIS land validation program. The ground measurements were conducted by the Remote Sensing Group (RSG) at the University of Arizona using an ASD spectrometer. The study reveals a good agreement between ASRVN and RSG HDRF for both MODIS Terra and Aqua with rmse ~ 0.01-0.025 in the 500 m MODIS land bands B1-B7. Obtained rmse is below uncertainties due to the spatial and seasonal variability of the bright calibration 1 km² area. While two MODIS instruments have a similar rmse in the visible bands, MODIS Aqua has a better agreement (lower rmse) with the ground data than MODIS Terra at wavelengths 0.87-2.1 μm. An independent overall good agreement of two MODIS instruments with the ground data indicates that the relative calibration of MODIS Terra and Aqua at medium-to-bright reflectance levels for the stated time period is significantly better than uncertainties of the ASRVN and ground data.     

Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data

A semi-physical fusion approach that uses the MODIS BRDF/Albedo land surface characterization product and Landsat ETM+ data to predict ETM+ reflectance on the same, an antecedent, or subsequent date is presented. The method may be used for ETM+ cloud/cloud shadow and SLC-off gap filling and for relative radiometric normalization. It is demonstrated over three study sites, one in Africa and two in the U.S. (Oregon and Idaho) that were selected to encompass a range of land cover land use types and temporal variations in solar illumination, land cover, land use, and phenology. Specifically, the 30 m ETM+ spectral reflectance is predicted for a desired date as the product of observed ETM+ reflectance and the ratio of the 500 m surface reflectance modeled using the MODIS BRDF spectral model parameters and the sun-sensor geometry on the predicted and observed Landsat dates. The difference between the predicted and observed ETM+ reflectance (prediction residual) is compared with the difference between the ETM+ reflectance observed on the two dates (temporal residual) and with respect to the MODIS BRDF model parameter quality. For all three scenes, and all but the shortest wavelength band, the mean prediction residual is smaller than the mean temporal residual, by up to a factor of three. The accuracy is typically higher at ETM+ pixel locations where the MODIS BRDF model parameters are derived using the best quality inversions. The method is most accurate for the ETM+ near-infrared (NIR) band; mean NIR prediction residuals are 9%, 12% and 14% of the mean NIR scene reflectance of the African, Oregon and Idaho sites respectively. The developed fusion approach may be applied to any high spatial resolution satellite data, does not require any tuning parameters and so may be automated, is applied on a per-pixel basis and is unaffected by the presence of missing or contaminated neighboring Landsat pixels, accommodates for temporal variations due to surface changes (e.g., phenological, land cover/land use variations) observable at the 500 m MODIS BRDF/Albedo product resolution, and allows for future improvements through BRDF model refinement and error assessment.     

Leaf chlorophyll content retrieval from airborne hyperspectral remote sensing imagery

Hyperspectral remote sensing has great potential for accurate retrieval of forest biochemical parameters. In this paper, a hyperspectral remote sensing algorithm is developed to retrieve total leaf chlorophyll content for both open spruce and closed forests, and tested for open forest canopies. Ten black spruce (Picea mariana (Mill.)) stands near Sudbury, Ontario, Canada, were selected as study sites, where extensive field and laboratory measurements were carried out to collect forest structural parameters, needle and forest background optical properties, and needle biophysical parameters and biochemical contents chlorophyll a and b. Airborne hyperspectral remote sensing imagery was acquired, within one week of ground measurements, by the Compact Airborne Spectrographic Imager (CASI) in a hyperspectral mode, with 72 bands and half bandwidth 4.25-4.36 nm in the visible and near-infrared region and a 2 m spatial resolution. The geometrical-optical model 4-Scale and the modified leaf optical model PROSPECT were combined to estimate leaf chlorophyll content from the CASI imagery. Forest canopy reflectance was first estimated with the measured leaf reflectance and transmittance spectra, forest background reflectance, CASI acquisition parameters, and a set of stand parameters as inputs to 4-Scale. The estimated canopy reflectance agrees well with the CASI measured reflectance in the chlorophyll absorption sensitive regions, with discrepancies of 0.06%-1.07% and 0.36%-1.63%, respectively, in the average reflectances of the red and red-edge region. A look-up-table approach was developed to provide the probabilities of viewing the sunlit foliage and background, and to determine a spectral multiple scattering factor as functions of leaf area index, view zenith angle, and solar zenith angle. With the look-up tables, the 4-Scale model was inverted to estimate leaf reflectance spectra from hyperspectral remote sensing imagery. Good agreements were obtained between the inverted and measured leaf reflectance spectra across the visible and near-infrared region, with R² = 0.89 to R² = 0.97 and discrepancies of 0.02%-3.63% and 0.24%-7.88% in the average red and red-edge reflectances, respectively. Leaf chlorophyll content was estimated from the retrieved leaf reflectance spectra using the modified PROSPECT inversion model, with R² = 0.47, RMSE = 4.34 μg/cm², and jackknifed RMSE of 5.69 μg/cm² for needle chlorophyll content ranging from 24.9 μg/cm² to 37.6 μg/cm². The estimates were also assessed at leaf and canopy scales using chlorophyll spectral indices TCARI/OSAVI and MTCI. An empirical relationship of simple ratio derived from the CASI imagery to the ground-measured leaf area index was developed (R² = 0.88) to map leaf area index. Canopy chlorophyll content per unit ground surface area was then estimated, based on the spatial distributions of leaf chlorophyll content per unit leaf area and the leaf area index.     

Quantification of two decades of shallow-water coral reef habitat decline in the Florida Keys National Marine Sanctuary using Landsat data (1984–2002)

The loss of coral reef habitats has been witnessed at a global scale including in the Florida Keys and the Caribbean. In addition to field surveys that can be spatially limited, remote sensing can provide a synoptic view of the changes occurring on coral reef habitats. Here, we utilize an 18-year time series of Landsat 5/TM and 7/ETM+ images to assess changes in eight coral reef sites in the Florida Keys National Marine Sanctuary, namely Carysfort Reef, Grecian Rocks, Molasses Reef, Conch Reef, Sombrero Reef, Looe Key Reef, Western Sambo and Sand Key Reef. Twenty-eight Landsat images (1984–2002) were used, with imagery gathered every 2 years during spring, and every 6 years during fall. The image dataset was georectified, calibrated to remote sensing reflectance and corrected for atmospheric and water-column effects. A Mahalanobis distance classification was trained for four habitat classes (‘coral’, ‘sand’, ‘bare hardbottom’ and ‘covered hardbottom’) using in situ ground-truthing data collected in 2003–2004 and using the spectral statistics from a 2002 image. The red band was considered useful only for benthic habitats in depths less than 6 m. Overall mean coral habitat loss for all sites classified by Landsat was 61% (3.4%/year), from a percentage habitat cover of 19% (1984) down to 7.6% (2002). The classification results for the eight different sites were critically reviewed. A detailed pixel by pixel examination of the spatial patterns across time suggests that the results range from ecologically plausible to unreliable due to spatial inconsistencies and/or improbable ecological successions. In situ monitoring data acquired by the Coral Reef Evaluation and Monitoring Project (CREMP) for the eight reef sites between 1996 and 2002 showed a loss in coral cover of 52% (8.7%/year), whereas the Landsat-derived coral habitat areas decreased by 37% (6.2%/year). A direct trend comparison between the entire CREMP percent coral cover data set (1996–2004) and the entire Landsat-derived coral habitat areas showed no significant difference between the two time series (ANCOVA; F-test, p = 0.303, n = 32), despite the different scales of measurements.     

Remote detection of Trichodesmium blooms in optically complex coastal waters: Examples with MODIS full-spectral data

Remote detection of the Trichodesmium spp. cyanobacteria blooms on the west Florida shelf (WFS) has been problematic due to optical complexity caused by sediment resuspension, coastal runoff, and bottom interference. By combining MODIS data measured by the ocean bands and land bands, an approach was developed to identify surface mats of Trichodesmium on the WFS. The approach first identifies possible bloom patches in MODIS FAI (floating algae index) 250 m resolution imagery derived from the Rayleigh-corrected reflectance at 667, 859, and 1240 nm. Then, spectral analysis examines the unique reflectance characteristics of Trichodesmium at 469, 488, 531, 551, and 555 nm due to specific optical properties (absorption, backscattering, and fluorescence) of the unusual pigments in Trichodesmium. These spectral characteristics (i.e., high-low-high-low-high reflectance at 469-488-531-551-555 nm, respectively) differentiate Trichodesmium mats unambiguously from other features observed in the FAI imagery, such as Sargassum spp. Tests in other coastal locations show that the approach is robust and applicable to other optically complex waters. Results shown here can help study Trichodesmium bloom dynamics (e.g., initiation and bloom formation) and may also help design future sensors to better detect and quantify Trichodesmium, an important N₂ fixer in the global oceans.     

An assessment of photosynthetic light use efficiency from space: Modeling the atmospheric and directional impacts on PRI reflectance

Estimation of photosynthetic light use efficiency (ε) from satellite observations is an important component of climate change research. The photochemical reflectance index, a narrow waveband index based on the reflectance at 531 and 570 nm, allows sampling of the photosynthetic activity of leaves; upscaling of these measurements to landscape and global scales, however, remains challenging. Only a few studies have used spaceborne observations of PRI so far, and research has largely focused on the MODIS sensor. Its daily global coverage and the capacity to detect a narrow reflectance band at 531 nm make it the best available choice for sensing ε from space. Previous results however, have identified a number of key issues with MODIS-based observations of PRI. First, the differences between the footprint of eddy covariance (EC) measurements and the MODIS footprint, which is determined by the sensor's observation geometry make a direct comparison between both data sources challenging and second, the PRI reflectance bands are affected by atmospheric scattering effects confounding the existing physiological signal. In this study we introduce a new approach for upscaling EC based ε measurements to MODIS. First, EC-measured ε values were “translated” into a tower-level optical PRI signal using AMSPEC, an automated multi-angular, tower-based spectroradiometer instrument. AMSPEC enabled us to adjust tower-measured PRI values to the individual viewing geometry of each MODIS overpass. Second, MODIS data were atmospherically corrected using a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which uses a time series approach and an image-based rather than pixel-based processing for simultaneous retrievals of atmospheric aerosol and surface bidirectional reflectance (BRDF). Using this approach, we found a strong relationship between tower-based and spaceborne reflectance measurements (r² = 0.74, p < 0.01) throughout the vegetation period of 2006. Swath (non-gridded) observations yielded stronger correlations than gridded data (r² = 0.58, p < 0.01) both of which included forward and backscatter observations. Spaceborne PRI values were strongly related to canopy shadow fractions and varied with different levels of ε. We conclude that MAIAC-corrected MODIS observations were able to track the site-level physiological changes from space throughout the observation period.     

Generation of high-resolution East Antarctic landfast sea-ice maps from cloud-free MODIS satellite composite imagery

A method to generate high spatio-temporal resolution maps of landfast sea ice from cloud-free MODIS composite imagery is presented. Visible (summertime) and thermal infrared (wintertime) cloud-free 20-day MODIS composite images are used as the basis for these maps, augmented by AMSR-E ASI sea-ice concentration composite images (when MODIS composite image quality is insufficient). The success of this technique is dependent upon efficient cloud removal during the compositing process. Example wintertime maximum (~ 374,000 km²) and summertime minimum (~ 112,000 km²) fast-ice maps for the entire East Antarctic coast are presented. The summertime minimum map provides the first high-resolution indication of multi-year fast-ice extent, which may be used to help assess changes in Antarctic sea-ice volume. The 2σ errors in fast-ice extent are estimated to be ± 2.98% when ≥ 90% of the fast-ice pixels in a 20-day period are classified using the MODIS composite, or ± 8.76 otherwise (when augmenting AMSR-E or the previous/next MODIS composite image is used to classify > 10% of the fast ice). Imperfect composite image quality, caused by persistent cloud, inaccurate cloud masking or a highly dynamic fast-ice edge, was the biggest impediment to automating the fast-ice detection procedure.     

Regional scale impacts of Tamarix leaf beetles (Diorhabda carinulata) on the water availability of western U.S. rivers as determined by multi-scale remote sensing methods

Tamarix leaf beetles (Diorhabda carinulata) have been widely released on western U.S. rivers to control introduced shrubs in the genus Tamarix. Part of the motivation to control Tamarix is to salvage water for human use. Information is needed on the impact of beetles on Tamarix seasonal leaf production and subsequent water use over wide areas and multiple cycles of annual defoliation. Here we combine ground data with high resolution phenocam imagery and moderate resolution (Landsat) and coarser resolution (MODIS) satellite imagery to test the effects of beetles on Tamarix evapotranspiration (ET) and leaf phenology at sites on six western rivers. Satellite imagery covered the period 2000 to 2010 which encompassed years before and after beetle release at each study site. Phenocam images showed that beetles reduced green leaf cover of individual canopies by about 30% during a 6–8 week period in summer, but plants produced new leaves after beetles became dormant in August, and over three years no net reduction in peak summer leaf production was noted. ET was estimated by vegetation index methods, and both Landsat and MODIS analyses showed that beetles reduced ET markedly in the first year of defoliation, but ET recovered in subsequent years. Over all six sites, ET decreased by 14% to 15% by Landsat and MODIS estimates, respectively. However, results were variable among sites, ranging from no apparent effect on ET to substantial reduction in ET. Baseline ET rates before defoliation were low, 394 mm yr^? 1 by Landsat and 314 mm yr^? 1 by MODIS estimates (20–25% of potential ET), further constraining the amount of water that could be salvaged. Beetle–Tamarix interactions are in their early stage of development on this continent and it is too soon to predict the eventual extent to which Tamarix populations will be reduced. The utility of remote sensing methods for monitoring defoliation was constrained by the small area covered by each phenocam image, the low temporal resolution of Landsat, and the low spatial resolution of MODIS imagery. Even combined image sets did not adequately reveal the details of the defoliation process, and remote sensing data should be combined with ground observations to develop operational monitoring protocols.     

An improved high spatial and temporal data fusion approach for combining Landsat and MODIS data to generate daily synthetic Landsat imagery

Because of low temporal resolution and cloud influence, many remote-sensing applications lack high spatial resolution remote-sensing data. To address this problem, this study introduced an improved spatial and temporal data fusion approach (ISTDFA) to generate daily synthetic Landsat imagery. This algorithm was designed to avoid the weaknesses of the spatial and temporal data fusion approach (STDFA) method, including the sensor difference and spatial variability. A weighted linear mixed model was used to adjust the spatial variability of surface reflectance. A linear-regression method was used to remove the influence of differences in sensor systems. This method was tested and validated in three study areas located in Xinjiang and Anhui province, China. The other two methods, the STDFA and the Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM), were also applied and compared in those three study areas. The results showed that the ISTDFA algorithm can generate daily synthetic Landsat imagery accurately, with correlation coefficient r equal to 0.9857 and root mean square error (RMSE) equal to 0.0195, which is superior to the STDFA method. The ISTDFA method had higher accuracy than ESTARFM in areas greater than 200 × 200 MODIS pixels while the ESTARFM method had higher accuracy than the ISTDFA method in small areas. The correlation coefficient r had a negative power relation with ratio of land-cover change pixels. A land-cover change of 20.25% pixels can lead to a reduced correlation coefficient r of 0.295 in the blue band. The accuracy of the ISTDFA method indicated a logarithmic relationship with the size of the applied area, so it is recommended for use in large-scale areas.