An in situ study of the effects of surface anisotropy on the remote sensing of burned savannah

This Letter presents field‐based evidence of the perturbing effects of surface anisotropy on the remote sensing of burned savannah. The analysis is based on bidirectional spectral reflectance data collected at different solar illumination angles and convolved to Moderate‐resolution Imaging Spectroradiometer (MODIS) reflective bands. Results from a grass savannah site show that burning reduces the anisotropy of the surface compared to its pre‐burn state. In contrast, at a shrub savannah site, burning reduces or increases surface anisotropy. Spectral indices defined from 1.240 µm and 2.130 µm reflectance, and 1.640 µm and 2.130 µm reflectance, provided stronger diurnal separation between burned and unburned areas than individual reflectance bands but do not eliminate anisotropic effects. The Normalized Difference Vegetation Index (NDVI) provides weak diurnal separation relative to these near‐ and mid‐infrared based indices. Implications of the findings are discussed for burned area mapping.     

Assessing the feasibility of sub-pixel burned area mapping in miombo woodlands of northern Mozambique using MODIS imagery

The goal of this study was to evaluate the feasibility of sub-pixel burned area detection in the miombo woodlands of northern Mozambique, using imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS). Multitemporal Landsat-7 ETM+ data were acquired to produce a high spatial resolution map of areas burned between mid-August and late September 2000, and a field campaign was conducted in early November 2000 to gather ground truth data. Mapping of burned areas was performed with an ensemble of classification trees and yielded a kappa value of 0.896. This map was subsequently degraded to a spatial resolution of 500 m, to produce an estimate of burned area fraction, at the MODIS pixel size. Correlation analysis between the sub-pixel burned area fraction map and the MODIS reflective channels 1-7 yielded low but statistically significant correlations for all channels. The better correlations were obtained for MODIS channels 2 (0.86 µm), 5 (1.24 µm) and 6 (1.64 µm). A regression tree was constructed to predict sub-pixel burned area fraction as a function of those MODIS channels. The resulting tree has nine terminal nodes and an overall root mean square error of 0.252. The regression tree analysis confirmed that MODIS channels 2, 5, and 6 are the best predictors of burned area fraction. It may be possible to improve these results considering, as an alternative to individual channels, some appropriate spectral indices used to enhance the burnt scar signal, and by including MODIS thermal data in the analysis. It may also be possible to improve the accuracy of sub-pixel burned area fraction using MODIS imagery by allowing the regression tree to automatically create linear combinations of individual channels, and by using ensembles of trees.     

Temporal burn scar evolution in tallgrass prairie based on field spectroscopy

In situ field spectroscopy samples were used to simulate several Moderate Resolution Imaging Spectroradiometer (MODIS) bands and indices commonly used for burned area detection. Each band or index was tested for its ability to differentiate between burned and unburned tallgrass prairie during several time periods from spring (when burning took place) to late summer (peak biomass) with three analysis of variance tests. The normalized difference vegetation index (NDVI), global environmental monitoring index (GEMI), global environmental monitoring index – burn scar (GEMI-B), and normalized burn ratio (NBR) indices, as well as MODIS band 7 (longwave mid-infrared; LWMIR), showed virtually no promise for differentiating burned from unburned areas for more than several days after the burn. Others, including the burned area index (BAI), Mid-infrared burn index (MIRBI), and MODIS bands 3 (red), 4 (near-infrared; NIR), 5 (longwave near-infrared; LWNIR), and 6 (shortwave mid-infrared; SWMIR) were able to differentiate between burned and unburned areas well into the growing season – in some cases, even through its entire length. The performance of particular bands and indices often depended on grazing, vegetation phenology, ash/char/soil reflectance, and factors that influenced pre-burn biomass.     

Evaluation of linear spectral unmixing and ΔNBR for predicting post‐fire recovery in a North American ponderosa pine forest

The Differenced Normalized Burn Ratio (ΔNBR) is widely used to map post‐fire effects in North America from multispectral satellite imagery, but has not been rigorously validated across the great diversity in vegetation types. The importance of these maps to fire rehabilitation crews highlights the need for continued assessment of alternative remote sensing approaches. To meet this need, this study presents a first preliminary comparison of immediate post‐fire char (black ash) fraction, as measured by linear spectral unmixing, and ΔNBR, with two quantitative one‐year post‐fire field measures indicative of canopy and sub‐canopy conditions: % live tree and dry organic litter weight (gm^?2). Image analysis was applied to Landsat 7 Enhanced Thematic Mapper (ETM+) imagery acquired both before and immediately following the 2000 Jasper Fire, South Dakota. Post‐fire field analysis was conducted one‐year post‐fire. Although the immediate post‐fire char fraction (r ² = 0.56, SE = 28.03) and ΔNBR (r ² = 0.55, SE = 29.69) measures produced similarly good predictions of the % live tree, the standard error in the prediction of litter weight with the char fraction method (r ² = 0.55, SE = 4.78) was considerably lower than with ΔNBR (r ² = 0.52, SE = 8.01). Although further research is clearly warranted to evaluate more field measures, in more fires, and across more fire regimes, the char fraction may be a viable approach to predict longer‐term indicators of ecosystem recovery and may potentially act as a surrogate retrospective measure of the fire intensity.     

Comparison of burned area estimates derived from SPOT-VEGETATION and Landsat ETM+ data in Africa: Influence of spatial pattern and vegetation type

An algorithm for burned area mapping in Africa based on classification trees was developed using SPOT-VEGETATION (VGT) imagery. The derived 1 km spatial resolution burned area maps were compared with 30 m spatial resolution maps obtained with 13 Landsat ETM+ scenes, through linear regression analysis. The procedure quantifies the bias in burned area estimation present in the low spatial resolution burned area map. Good correspondence was observed for seven sites, with values of the coefficient of determination (R²) ranging from 0.787 to 0.983. Poorer agreement was observed in four sites (R² values between 0.257 and 0.417), and intermediate values of R² (0.670 and 0.613) were obtained for two sites. The observed variation in the level of agreement between the Landsat and VGT estimates of area burned results from differences in the spatial pattern and size distribution of burns in the different fire regimes encompassed by our analysis. Small and fragmented burned areas result in large underestimation at 1 km spatial resolution. When large and compact burned areas dominate the landscape, VGT estimates of burned area are accurate, although in certain situations there is some overestimation. Accuracy of VGT burned area estimates also depends on vegetation type. Results showed that in forest ecosystems VGT maps underestimate substantially the amount of burned area. The most accurate estimates were obtained for woodlands and grasslands. An overall linear regression fitted with the data from the 13 comparison sites revealed that there is a strong relationship between VGT and Landsat estimates of burned area, with a value of R² of 0.754 and a slope of 0.803. Our findings indicate that burned area mapping based on 1 km spatial resolution VGT data provides adequate regional information.     

Monitoring boreal forest leaf area index across a Siberian burn chronosequence: a MODIS validation study

Landscapes containing differing amounts of ecological disturbance provide an excellent opportunity to validate and better understand the emerging Moderate Resolution Imaging Spectrometer (MODIS) vegetation products. Four sites, including 1‐year post‐fire coniferous, 13‐year post‐fire deciduous, 24‐year post‐fire deciduous, and >100 year old post‐fire coniferous forests, were selected to serve as a post‐fire chronosequence in the central Siberian region of Krasnoyarsk (57.3°N, 91.6°E) with which to study the MODIS leaf area index (LAI) and vegetation index (VI) products. The collection 4 MODIS LAI product correctly represented the summer site phenologies, but significantly underestimated the LAI value of the >100 year old coniferous forest during the November to April time period. Landsat 7‐derived enhanced vegetation index (EVI) performed better than normalized difference vegetation index (NDVI) to separate the deciduous and conifer forests, and both indices contained significant correlation with field‐derived LAI values at coniferous forest sites (r ² = 0.61 and r ² = 0.69, respectively). The reduced simple ratio (RSR) markedly improved LAI prediction from satellite measurements (r ² = 0.89) relative to NDVI and EVI. LAI estimates derived from ETM+ images were scaled up to evaluate the 1 km resolution MODIS LAI product; from this analysis MODIS LAI overestimated values in the low LAI deciduous forests (where LAI<5) and underestimated values in the high LAI conifer forests (where LAI>6). Our results indicate that further research on the MODIS LAI product is warranted to better understand and improve remote LAI quantification in disturbed forest landscapes over the course of the year.     

Analysis of the relationship between spatial pattern and spectral detectability of areas burned in southern Africa using satellite data

Fires in Africa affect atmospheric emissions and carbon sequestration, landscape patterns, and regional and global climatic conditions. Studies of these effects require accurate estimation of the extent of measurable fire events. The goal of this study was to assess the influence of burned area spatial patterns on the spectral detectability of burned areas. Six Landsat‐7 ETM+ images from the southern Africa were used for burned area mapping and spatial pattern analysis, while contemporaneous MODIS 500 m spatial resolution images were used to measure the spectral detectability of burned areas. Using a 15 by 15 km sample quadrats analysis, we showed that above a burned area proportion threshold of approximately 0.5 the spectral detectability of burned areas increase due to the decrease in the number of mixed pixels. This was spatially related to the coalescence of burned patches and the decrease in the total burned area perimeter. Simple burned area shapes were found at the Botswana site, where the absence of tree cover and the presence of bright surfaces (soil and dry grass) enhanced the spectral contrast of the burned surfaces, thus enabling better estimates of burned area extent. At the Zambia and Congo sites, landscape fragmentation due to human activity and the presence of a tree vegetation layer, respectively, contribute to the presence of small burned area patches, which may remain undetectable using moderate spatial resolution satellite imagery, leading to less accurate burned area extent estimates.     

Comparison of Landsat 8 OLI and Landsat 7 ETM+ for estimating grassland LAI using model inversion and spectral indices: case study of Mpumalanga,South Africa

The leaf area index (LAI) is the key biophysical indicator used to assess the condition of rangeland. In this study, we investigated the implications of narrow spectral response, high radiometric resolution (12 bits), and higher signal-to-noise ratio of the Landsat 8 Operational Land Imager (OLI) sensor for the estimation of LAI. The Landsat 8 LAI estimates were compared to that of its predecessors, namely Landsat 7 Enhanced Thematic Mapper Plus (ETM+) (8 bits). Furthermore, we compared the radiative transfer model (RTM) and spectral indices approaches for estimating LAI on rangeland systems in South Africa. The RTM was inverted using artificial neural network (ANN) and lookup table (LUT) algorithms. The accuracy of the models was higher for Landsat 8 OLI, where ANN (root mean squared error, RMSE = 0. 13; R² = 0. 89), LUT (RMSE = 0. 25; R² = 0. 50), compared to Landsat 7 ETM+, where ANN (RMSE = 0. 35; R² = 0. 60), LUT (RMSE = 0. 38; R² = 0. 50). Compared to an empirical approach, the RTM provided higher accuracy. In conclusion, Landsat 8 OLI provides an improvement for the estimation of LAI over Landsat 7 ETM+. This is useful for rangeland monitoring.     

An active-fire based burned area mapping algorithm for the MODIS sensor

We present an automated method for mapping burned areas using 500-m Moderate Resolution Imaging Spectroradiometer (MODIS) imagery coupled with 1-km MODIS active fire observations. The algorithm applies dynamic thresholds to composite imagery generated from a burn-sensitive vegetation index and a measure of temporal texture. Cumulative active fire maps are used to guide the selection of burned and unburned training samples. An accuracy assessment for three geographically diverse regions (central Siberia, the western United States, and southern Africa) was performed using high resolution burned area maps derived from Landsat imagery. Mapped burned areas were accurate to within approximately 10% in all regions except the high-tree-cover sub-region of southern Africa, where the MODIS burn maps underestimated the area burned by 41%. We estimate the minimum detectable burn size for reliable detection by our algorithm to be on the order of 120 ha.     

Regionally adaptable dNBR-based algorithm for burned area mapping from MODIS data

Recent advances in instrument design have led to considerable improvements in wildfire mapping at regional and global scales. Global and regional active fire and burned area products are currently available from various satellite sensors. While only global products can provide consistent assessments of fire activity at the global, hemispherical or continental scales, the efficiency of their performance differs in various ecosystems. The available regional products are hard-coded to the specifics of a given ecosystem (e.g. boreal forest) and their mapping accuracy drops dramatically outside the intended area. We present a regionally adaptable semi-automated approach to mapping burned area using Moderate Resolution Imaging Spectroradiometer (MODIS) data. This is a flexible remote sensing/GIS-based algorithm which allows for easy modification of algorithm parameterization to adapt it to the regional specifics of fire occurrence in the biome or region of interest. The algorithm is based on Normalized Burned Ratio differencing (dNBR) and therefore retains the variability of spectral response of the area affected by fire and has the potential to be used beyond binary burned/unburned mapping for the first-order characterization of fire impacts from remotely sensed data. The algorithm inputs the MODIS Surface Reflectance 8-Day Composite product (MOD09A1) and the MODIS Active Fire product (MOD14) and outputs yearly maps of burned area with dNBR values and beginning and ending dates of mapping as the attributive information. Comparison of this product with high resolution burn scar information from Landsat ETM+ imagery and fire perimeter data shows high levels of accuracy in reporting burned area across different ecosystems. We evaluated algorithm performance in boreal forests of Central Siberia, Mediterranean-type ecosystems of California, and sagebrush steppe of the Great Basin region of the US. In each ecosystem the MODIS burned area estimates were within 15% of the estimates produced by the high resolution base with the R² between 0.87 and 0.99. In addition, the spatial accuracy of large burn scars in the boreal forests of Central Siberia was also high with Kappa values ranging between 0.76 and 0.79.     

Mapping fire scars in a southern African savannah using Landsat imagery

The spectral, spatial and temporal characteristics of the Landsat data record make it appropriate for mapping fire scars. Twenty-two annual fire scar maps from 1972–2002 were produced from historical Landsat imagery for a semi-arid savannah landscape on the South Africa–Botswana border, centred over Madikwe Game Reserve (MGR) in South Africa. A principal components transformation (PCT) helped differentiate the spectral signal of fire scars in each image. A simple, nonparametric, supervised classification (parallelepiped) of the PCT data differentiated burned and unburned areas. During most years, fire occurrences and the percentage of area burned annually were lowest in Botswana, highest in MGR, and intermediate in South Africa outside MGR. These fire scar maps are aiding MGR managers, who are endeavouring to restore a more active fire regime following decades of fire exclusion.     

Comparisons of land cover and LAI estimates derived from ETM+ and MODIS for four sites in North America: a quality assessment of 2000/2001 provisional MODIS products

The MODIS land science team produces a number of standard products, including land cover and leaf area index (LAI). Critical to the success of MODIS and other sensor products is an independent evaluation of product quality. In that context, we describe a study using field data and Landsat ETM+ to map land cover and LAI at four 49-km² sites in North America containing agricultural cropland (AGRO), prairie grassland (KONZ), boreal needleleaf forest, and temperate mixed forest. The purpose was to: (1) develop accurate maps of land cover, based on the MODIS IGBP (International Geosphere-Biosphere Programme) land cover classification scheme; (2) derive continuous surfaces of LAI that capture the mean and variability of the LAI field measurements; and (3) conduct initial MODIS validation exercises to assess the quality of early (i.e., provisional) MODIS products. ETM+ land cover maps varied in overall accuracy from 81% to 95%. The boreal forest was the most spatially complex, had the greatest number of classes, and the lowest accuracy. The intensive agricultural cropland had the simplest spatial structure, the least number of classes, and the highest overall accuracy. At each site, mapped LAI patterns generally followed patterns of land cover across the site. Predicted versus observed LAI indicated a high degree of correspondence between field-based measures and ETM+ predictions of LAI. Direct comparisons of ETM+ land cover maps with Collection 3 MODIS cover maps revealed several important distinctions and similarities. One obvious difference was associated with image/map resolution. ETM+ captured much of the spatial complexity of land cover at the sites. In contrast, the relatively coarse resolution of MODIS did not allow for that level of spatial detail. Over the extent of all sites, the greatest difference was an overprediction by MODIS of evergreen needleleaf forest cover at the boreal forest site, which consisted largely of open shrubland, woody savanna, and savanna. At the agricultural, temperate mixed forest, and prairie grassland sites, ETM+ and MODIS cover estimates were similar. Collection 3 MODIS-based LAI estimates were considerably higher (up to 4 m² m⁻²) than those based on ETM+ LAI at each site. There are numerous probable reasons for this, the most important being the algorithms' sensitivity to MODIS reflectance calibration, its use of a prelaunch AVHRR-based land cover map, and its apparent reliance on mainly red and near-IR reflectance. Samples of Collection 4 LAI products were examined and found to consist of significantly improved LAI predictions for KONZ, and to some extent for AGRO, but not for the other two sites. In this study, we demonstrate that MODIS reflectance data are highly correlated with LAI across three study sites, with relationships increasing in strength from 500 to 1000 m spatial resolution, when shortwave-infrared bands are included.     

Comparison of MODIS and Landsat TM5 images for mapping tempo–spatial dynamics of Secchi disk depths in Poyang Lake National Nature Reserve,China

Landsat has successfully been applied to map Secchi disk depth of inland water bodies. Operational use for monitoring a dynamic variable like Secchi disk depth is however limited by the 16‐day overpass cycle of the Landsat system and cloud cover. Low spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) image captured twice a day could potentially overcome these problems. However, its potential for mapping Secchi disk depth of inland water bodies has so far rarely been explored. This study compared two image sources, MODIS and Landsat Thematic Mapper (TM), for mapping the tempo–spatial dynamics of Secchi disk depth in Poyang Lake National Nature Reserve, China. Secchi disk depths recorded at weekly intervals from April to October in 2004 and 2005 were related to 5 Landsat TM and 22 MODIS images respectively. Two multiple regression models including the blue and red bands of Landsat TM and MODIS respectively explained 83% and 88% of the variance of the natural logarithm of Secchi disk depth. The standard errors of the predictions were 0.20 and 0.37 m for Landsat TM and MODIS‐based models. A high correlation (r = 0.94) between the predicted Secchi disk depth derived from the two models was observed. A discussion of advantages and disadvantages of both sensors leads to the conclusion that MODIS offers the possibility to monitor water transparency more regularly and cheaply in relatively big and frequently cloud covered lakes as is with Poyang Lake.     

Model development and comparison to predict softwood and hardwood per cent cover using high and medium spatial resolution imagery

The recent availability of high spatial resolution multispectral scanners provides an opportunity to adapt existing methods and test models to derive spatially explicit forest type and per cent cover information at the Landsat pixel level. A regression modelling methodology was applied for scaling‐up high resolution (IKONOS) to medium spatial resolution satellite imagery (Landsat) to predict softwood and hardwood forest type and density (per cent cover) in a northern Maine study area. Regression relationships (63 different models) were developed and compared. The model variables included vegetation indices and several date (season) combinations of Landsat Enhanced Thematic Mapper Plus (ETM+) imagery (August, September, October and May).

A model incorporating all variables from four dates of Landsat ETM+ imagery produced the highest coefficient of variation in predicting both softwood (0.655) and hardwood cover (0.66). The addition of vegetation indices with the six ETM+ reflected bands did not significantly improve or detract from the regression relationships for any of the multi‐date or single date models examined. A two‐date combination of October and May variables provided an acceptable (and arguably more cost‐effective) model as the adjusted R ² value was 0.645 for softwood and 0.649 for hardwood. A significant result was that all single‐date models produced inferior results with a sharp drop in adjusted R ², compared with the multi‐date seasonal models. This research has demonstrated that the regression models including multi‐date variables produce good results and can provide spatially explicit forest type and stand structure data that has been difficult or infeasible to obtain from medium spatial resolution imagery using traditional classification methods.     

Application of multitemporal Landsat data to monitor land cover changes in the Eastern Nile Delta region,Egypt

Due to the progressive increase in development of desert land in Egypt, the demand for efficient and accurate land cover change information is increasing. In this study, we apply the methodology of post‐classification change detection to map and monitor land cover change patterns related to agricultural development and urban expansion in the desert fringes of the Eastern Nile Delta region. Using a hybrid classification approach, we employ multitemporal Landsat Thematic Mapper (TM)/Enhanced Thematic Mapper Plus (ETM+) images from 1984, 1990 and 2003 to produce three land cover/land‐use maps. Post‐classification comparison of these maps was used to obtain ‘from–to’ statistics and change detection maps. The change detection results show that agricultural development increased by 14% through the study period. The average annual rate of land reclamation during 1990–2003 (4511 ha a^?1) was comparable to that during 1984–1990 (4644 ha a^?1), reflecting a systematic national plan for desert reclamation that went into effect. We find that the increase in urbanization (by ca 21 300 ha) during 1990–2003 was predominantly due to encroachment into traditionally cultivated land at the fringes of urban centres. Our results accurately quantify the land cover changes and delineate their spatial patterns, demonstrating the utility of Landsat data in analysing landscape dynamics over time. Such information is critical for making efficient and sustainable policies for resource management.     

Spatio-temporal distribution of burned areas by ecoregions in Mexico and Central America

Fire activity in Mexico and Central America, and its associated emissions, has impacts across multiple scales. On the local-to-regional scale, fire activity impacts land use, productivity, and biodiversity. On the regional-to-global scale, fire activity impacts hydrological, biogeochemical, and atmospheric processes. A consistent, reliable, large-scale characterization of the spatial and temporal distribution of fire burned area is required to assess its environmental impacts and to support natural resources’ management. The spatial and temporal distributions of fire burned areas in ecoregions of Mexico and Central America are evaluated in this study for the period 2001–2014, using the satellite Moderate Resolution Imaging Spectroradiometer (MODIS) MCD45 Burned Area data set. The methodology combines the 500 m burned area product with a MODIS land cover product and a map of North American land cover to calculate the spatiotemporal variability of fire activity as a function of land-use type.

The total burned area over Mexico and Central America over the period 2001–2014 was found to be 614,243.5 km², but with significant interannual variability over the 14 years included in the study. Indeed, the minimum burned area over the period was 9892.25 km² in 2014 and the maximum was 37,669.50 km² in 2011, a fourfold increase. Burned areas were found to be concentrated in northern Mexico and on the Pacific coast, mainly from October to June. Agricultural burned areas accounted for 37% and 43% of total detected burns in Mexico and Central America, respectively. The largest extent of burned surface occurs in May for most land-cover types. The maximum density of burned areas occurred in the tropical dry forests ecoregion during the dry season. Both in Mexico and Central America, burned area anomalies have significant anti-correlation with precipitation anomalies.     

Evaluation of novel thermally enhanced spectral indices for mapping fire perimeters and comparisons with fire atlas data

We evaluated the potential of two novel thermally enhanced Landsat Thematic Mapper (TM)‐derived spectral indices for discriminating burned areas and for producing fire perimeter data (as a potential surrogate to digital fire atlas data) within two wildland fires (1985 and 1993) in ponderosa pine (Pinus ponderosa) forests of the Gila Wilderness, New Mexico, USA. Image‐derived perimeters (manually produced and classified from an index image) were compared to fire perimeters recorded within a digitized fire atlas. For each fire, the highest spectral separability was achieved using the newly proposed Normalized Burn Ratio‐Thermal (NBRT₁) index (M = 1.18, 1.76, for the two fires respectively). Correspondence between fire atlas and manually digitized fire perimeters was high. Landsat imagery may be a useful supplement to existing historical fire perimeters mapping methods, but the timing of the post‐fire image will strongly influence the separability of burned and unburned areas.     

Comparing MODIS and ETM+ data for regional and global land classification

Nearly simultaneous reflectance data sets from the Landsat 7 Enhanced Thematic Mapper Plus (ETM+), at 30-m resolution, and the Terra satellite instrument MODIS, at 500-m resolution, are compared for their ability to map fractional coverage of surface types over large areas. Lower spatial resolution MODIS classification results are generally comparable those of ETM+, with discrepancies for some regions with mixed surface types. Analysis of laboratory and field spectra suggests an ambiguity, the “brightness ambiguity”, which can prevent accurate area estimation of pixels having two or more surface types. This ambiguity, plus general mathematical inversion issues, can account for the discrepancy. Thus, occasional high-resolution measurements, as from Landsat 7, are necessary to refine estimations of large area surface types from MODIS and similar lower spatial resolution instruments.     

Multitemporal assessment of the geomorphologic evolution of the Restinga of Marambaia,Rio de Janeiro,Brazil

The Restinga of Marambaia is an emerged sand bar located between the Sepetiba Bay and the South Atlantic Ocean, on the south‐east coast of Brazil. The objective of this study was to observe the geomorphologic evolution of the coastal zone of the Restinga of Marambaia using multitemporal satellite images acquired by multisensors from 1975 to 2004. The images were digitally segmented by a region growth algorithm and submitted to an unsupervised classification procedure (ISOSEG) followed by a raster edit based on visual interpretation. The image time‐series showed a general trend of decrease in the total sand bar area with values varying from 80.61 km² in 1975 to 78.15 km² in 2004. The total area calculation based on the 1975 and 1978 Landsat MSS data was shown to be super‐estimated in relation to the Landsat TM, Landsat ETM+, and CBERS‐2 CCD data. These differences can also be associated to the relatively poorer spatial resolution of the MSS data, nominally 79 m, against the 20 m of the CCD data and 30 m of the TM and ETM+ data. For the estimates of the width in the central portion of the sand bar the variation was from 158 m (1975) to 100 m (2004). The formation of a spit in the northern region of the study area was visually observed. The area of the spit was estimated, with values varying from 0.82 km² (1975) to 0.55 km² (2004).     

Leaf area index mapping in northern Canada

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

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