Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions

The full realization of the potential of remote sensing as a source of environmental information requires an ability to generalize in space and time. Here, the ability to generalize in space was investigated through an analysis of the transferability of predictive relations for the estimation of tropical forest biomass from Landsat TM data between sites in Brazil, Malaysia and Thailand. The data sets for each test site were acquired and processed in a similar fashion to facilitate the analyses. Three types of predictive relation, based on vegetation indices, multiple regression and feedforward neural networks, were developed for biomass estimation at each site. For each site, the strongest relationships between the biomass predicted and that measured from field survey was obtained with a neural network developed specifically for the site (r>0.71, significant at the 99% level of confidence). However, with each type of approach problems in transferring a relation to another site were observed. In particular, it was apparent that the accuracy of prediction, as indicated by the correlation coefficient between predicted and measured biomass, declined when a relation was transferred to a site other than that upon which it was developed. Part of this problem lies with the observed variation in the relative contribution of the different spectral wavebands to predictive relations for biomass estimation between sites. It was, for example, apparent that the spectral composition of the vegetation indices most strongly related to biomass differed greatly between the sites. Consequently, the relationship between predicted and measured biomass derived from vegetation indices differed markedly in both strength and direction between sites. Although the incorporation of test site location information into an analysis resulted in an increase in the strength of the relationship between predicted and actual biomass, considerable further research is required on the problems associated with transferring predictive relations.     

Improving tropical forest mapping using multi-date Landsat TM data and pre-classification image smoothing

The present study explores the possibility of using Landsat imagery for mapping tropical forest types with relevance to forest ecosystem services. The central part in the classification process is the use of multi-date image data and pre-classification image smoothing. The study argues that multi-date imagery contains information on phenological and canopy structural properties, and shows how the use of multi-date imagery has a significant impact on classification accuracy. Furthermore, the study shows the value of applying small kernel smoothing filters to reduce in-class spectral variability and enhance between-class spectral separability. Making use of these approaches and a maximum likelihood algorithm, six tropical forest types were classified with an overall accuracy of 90.94%, and with individual forest classes mapped with accuracies above 75.19% (user's accuracy) and above 74.17% (producer's accuracy).     

Image segmentation for humid tropical forest classification in Landsat TM data

Humid tropical forest types have low spectral separability in Landsat TM data due to highly textured reflectance patterns at the 30m spatial resolution. Two methods of reducing local spectral variation, low-pass spatial filtering and image segmentation, were examined for supervised classification of 10 forest types in TM data of Peruvian Amazonia. The number of forest classes identified at over 90% accuracy increased from one in raw imagery to three in filtered imagery, and six in segmented imagery. The ability to derive less generalised tropical forest classes may allow greater use of classified imagery in ecology and conservation planning.     

Mapping forest vegetation using Landsat TM imagery and a canopy reflectance model

Estimates of mean tree size and cover for each forest stand from an invertible forest canopy reflectance model are part of a new forest vegetation mapping system. Image segmentation defines stands which are sorted into general growth forms using per-pixel image classifications. Ecological models based on terrain relations predict species associations for the conifer, hardwood, and brush growth forms. The combination of the model-based estimates of tree size and cover with species associations yields general-purpose vegetation maps useful for a variety of land management needs. Results of timber inventories in the Tahoe and Stanislaus National Forests indicate the vegetation maps form a useful basis for stratification. Patterns in timber volumes for the strata reveal that the cover estimates are more reliable than the tree size estimates. A map accuracy assessment of the Stanislaus National Forest shows high overall map accuracy and also illustrates the problems in estimating tree size.     

Detection of forest harvest type using multiple dates of Landsat TM imagery

A simple and relatively accurate technique for classifying time-series Landsat Thematic Mapper (TM) imagery to detect levels of forest harvest is the topic of this research. The accuracy of multidate classification of the normalized difference vegetation index (NDVI) and the normalized difference moisture index (NDMI) were compared and the effect of the number of years (1–3, 3–4, 5–6 years) between image acquisition on forest change accuracy was evaluated. When Landsat image acquisitions were only 1–3 years apart, forest clearcuts were detected with producer's accuracy ranging from 79% to 96% using the RGB-NDMI classification method. Partial harvests were detected with lower producer's accuracy (55–80%) accuracy. The accuracy of both clearcut and partial harvests decreased as time between image acquisition increased. In all classification trials, the RGB-NDMI method produced significantly higher accuracies, compared to the RGB-NDVI. These results are interesting because the less common NDMI (using the reflected middle infrared band) outperformed the more popular NDVI. In northern Maine, industrial forest landowners have shifted from clearcutting to partial harvest systems in recent years. The RGB-NDMI change detection classification applied to Landsat TM imagery collected every 2–3 years appears to be a promising technique for monitoring forest harvesting and other disturbances that do not remove the entire overstory canopy.     

Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon

The complicated forest stand structure and associated abundant tree species in the Amazon often induce difficulty in estimating aboveground biomass (AGB) using remotely sensed data. This paper explores AGB estimation using Landsat Thematic Mapper (TM) data in the eastern and western Brazilian Amazon, and discusses the impacts of forest stand structure on AGB estimation. Estimating AGB is still a challenging task, especially for the sites with complicated biophysical environments. The TM spectral responses are more suitable for AGB estimation in the sites with relatively simple forest stand structure than for the sites with complicated forest stand structure. Conversely, textures appear more important than spectral responses in AGB estimation in the sites with complicated forest stand structure. A combination of spectral responses and textures improves AGB estimation performance. Different study areas having various biophysical conditions affect AGB estimation performance.     

Estimation of timber volume in a coniferous plantation forest using Landsat TM

Optimisation of economic return from forests requires that comprehensive forest inventory data be available to support the design of harvesting strategies. Such inventory data can potentially be obtained by remote sensing. This study investigates the accuracy with which wood volume (m3 ha -1) in a plantation forest can be calculated from Landsat TM data at the pixel and foreststand spatial scales. Wood volumes were estimated from regression analysis, nonparametric line-fitting, and an N-dimensional K-nearest-neighbour classification scheme. At the pixel scale, relations between Landsat data and measured wood volume were found to be significant but weak, with r2 values of 0.3, and with correspondingly poor estimates of wood volume (root-mean-square errors rmses of 100 m3 ha -1). By averaging the pixel-scale estimates, wood volume estimates of acceptable accuracy were obtained for forest-stand areas of about 40 ha (rmses of 46 m3 ha -1). Parametric regression performed slightly better overall than non-parametric line fitting techniques for estimating wood volume. Estimates of similar accuracy to those obtained by regression were also given by NK-classification at the pixel-scale, provided K was large ( 15), although the classifier produced biased results at the forest-stand scale. It is concluded that Landsat TM only provides an acceptable data source for estimating wood volumes in plantation forests for areas of about 40 ha and larger. The very low dynamic range in the Landsat data is probably a significant factor limiting its use for inventory at more detailed scales.     

Assessment of tropical forest degradation by selective logging and fire using Landsat imagery

Many studies have assessed the process of forest degradation in the Brazilian Amazon using remote sensing approaches to estimate the extent and impact by selective logging and forest fires on tropical rain forest. However, only a few have estimated the combined impacts of those anthropogenic activities. We conducted a detailed analysis of selective logging and forest fire impacts on natural forests in the southern Brazilian Amazon state of Mato Grosso, one of the key logging centers in the country. To achieve this goal a 13-year series of annual Landsat images (1992-2004) was used to test different remote sensing techniques for measuring the extent of selective logging and forest fires, and to estimate their impact and interaction with other land use types occurring in the study region. Forest canopy regeneration following these disturbances was also assessed. Field measurements and visual observations were conducted to validate remote sensing techniques. Our results indicated that the Modified Soil Adjusted Vegetation Index aerosol free (MSAVI_af) is a reliable estimator of fractional coverage under both clear sky and under smoky conditions in this study region. During the period of analysis, selective logging was responsible for disturbing the largest proportion (31%) of natural forest in the study area, immediately followed by deforestation (29%). Altogether, forest disturbances by selective logging and forest fires affected approximately 40% of the study site area. Once disturbed by selective logging activities, forests became more susceptible to fire in the study site. However, our results showed that fires may also occur in undisturbed forests. This indicates that there are further factors that may increase forest fire susceptibility in the study area. Those factors need to be better understood. Although selective logging affected the largest amount of natural forest in the study period, 35% and 28% of the observed losses of forest canopy cover were due to forest fire and selective logging combined and to forest fire only, respectively. Moreover, forest areas degraded by selective logging and forest fire is an addition to outright deforestation estimates and has yet to be accounted for by land use and land cover change assessments in tropical regions. Assuming that this observed trend of land use and land cover conversion continues, we predict that there will be no undisturbed forests remaining by 2011 in this study site. Finally, we estimated that 70% of the total forest area disturbed by logging and fire had sufficiently recovered to become undetectable using satellite data in 2004.     

Comparison of two types of forest disturbance using multitemporal Landsat TM/ETM+ imagery and field vegetation data

Forest disturbances influence many landscape processes, including changes in microclimate, hydrology, and soil erosion. We analyzed the spectral response and temporal progress of two types of disturbances of spruce forest (bark beetle outbreak and clear-cuts) in the central part of Šumava Mountains at the border between the Czech Republic and Germany, Central Europe. The bark beetle (Ips typographus [L.]) outbreak in this region in the last 20 years resulted in regional-scale spruce forest decay. Clear-cutting was done here to prevent further bark-beetle propagation in the buffer zones.The aim of the study is to identify the differences in spectral response between the two types of forest disturbances and their temporal dynamics. General trends were analyzed throughout the study area, with sampled disturbance areas selected to assess the relationship between field vegetation data and their spectral response. Thirteen Landsat TM/ETM+ scenes from 1985 to 2007 were used for the assessment. The following spectral indices were estimated: NDMI, Tasseled Cap (Brightness, Greenness, Wetness), DI, and DI′. The DI′, Wetness, and Brightness indices show the highest sensitivity to forest disturbance for both disturbance types (clear-cuts and bark beetle outbreak). The multitemporal analysis distinguished three different stages of development. The highest spectral differences between the clear-cuts and the bark beetle disturbances were found in the period between 1996 and 2004 with increased levels of forest disturbance (repeated measures ANOVA, Scheffé post hoc test; p ≤ 0.05). Clear-cut disturbance resulted in significantly higher spectral differences from the original forest and occurred as a more discrete event in comparison to bark beetle outbreak.     

Predicting forest successional stages using multitemporal Landsat imagery with forest inventory and analysis data

Forest succession is an important ecological process that has profound biophysical, biological and biogeochemical implications in terrestrial ecosystems. Therefore, information on forest successional stages over an extensive forested landscape is crucial for us to understand ecosystem processes, such as carbon assimilation and energy interception. This study explored the potential of using Forest Inventory and Analysis (FIA) plot data to extract forest successional stage information from remotely sensed imagery with three widely used predictive models, linear regression (LR), decision trees (DTs) and neural networks (NNs). The predictive results in this study agree with previous findings that multitemporal Landsat Thematic Mapper (TM) imagery can improve the accuracy of forest successional stage prediction compared to models using a single image. Because of the overlap of spectral signatures of forests in different successional stages, it is difficult to accurately separate forest successional stages into more than three broad age classes (young, mature and old) with reasonable accuracy based on the age information of FIA plots and the spectral data of the plots from Landsat TM imagery. Given the mixed spectral response of forest age classes, new approaches need to be explored to improve the prediction of forest successional stages using FIA data.     

Land-cover classification in a moist tropical region of Brazil with Landsat TM imagery

This research aims to improve land-cover classification accuracy in a moist tropical region in Brazil by examining the use of different remote sensing-derived variables and classification algorithms. Different scenarios based on Landsat Thematic Mapper (TM) spectral data and derived vegetation indices and textural images, and different classification algorithms - maximum likelihood classification (MLC), artificial neural network (ANN), classification tree analysis (CTA), and object-based classification (OBC), were explored. The results indicated that a combination of vegetation indices as extra bands into Landsat TM multispectral bands did not improve the overall classification performance, but the combination of textural images was valuable for improving vegetation classification accuracy. In particular, the combination of both vegetation indices and textural images into TM multispectral bands improved overall classification accuracy by 5.6% and kappa coefficient by 6.25%. Comparison of the different classification algorithms indicated that CTA and ANN have poor classification performance in this research, but OBC improved primary forest and pasture classification accuracies. This research indicates that use of textural images or use of OBC are especially valuable for improving the vegetation classes such as upland and liana forest classes having complex stand structures and having relatively large patch sizes.     

Characterizing landscape changes in central Rondonia using Landsat TM imagery

An analysis of landscape changes in a region of pioneer settlements in central Rondonia, western Brazilian Amazon, was derived from Landsat TM data. Total deforested area increased from 206 x 103 ha in 1977, to 565 x 103 ha in 1985 and to 1210 x 103 ha, or 35.5% of the region, in 1995. Eighty-one per cent of the total 1995 deforestation had occurred in regions within 12.5km from areas of pioneer colonization deforested by 1977. Deforested area exceeded 79% in regions within 12.5km from the region's first road.     

Lineament mapping in a tropical environment using Landsat imagery

Remote sensing has proved to be a useful tool in lineament identification and mapping. This study demonstrates the use of multispectral Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM?+) satellite data obtained over two acquisition dates in 1990 and 2002 for lineament interpretation in a Malaysian tropical environment. A digital elevation model (DEM) was generated to improve the interpretation. We found that most of the major orientations in the field station could be successfully detected from the remotely sensed imagery. The results from the study show that the remote sensing technique is capable of extracting lineament trends in an inaccessible tropical forest.     

Using Landsat TM data to estimate carbon release from burned biomass in an Alaskan spruce forest complex

Fire disturbance in boreal forests can release carbon to the atmosphere stored in both the aboveground vegetation and the organic soil layer. Estimating pyrogenic emissions of carbon released during biomass burning in these forests is useful for understanding and estimating global carbon budgets. In this work, we have developed a method to estimate carbon efflux for the burned black spruce in an Alaskan forest by combining information derived from Landsat Thematic Mapper (TM) data and field measurements. We have used the spatial and spectral information of TM data to identify and measure two important factors: pre-burn black spruce density and burn severity. Field measurements provided estimates of aboveground and ground layer carbon per unit area for the pre-burn Landsat spectral classes, and percentage of carbon consumed for the post-burn Landsat spectral classes. Carbon release estimates for the burned black spruce were computed using field data and the co-occurrence of the pre-burn and post-burn spectral classes. The estimated carbon released was 39.9tha-1, which is 57% greater than an estimate computed using AVHRR data and estimates of pre-burn biomass and carbon fractions consumed that were not site specific or spatially varying. We conclude that the spectral bands and spatial resolution of Landsat TM data provide the potential for improved estimates of pyrogenic carbon efflux relative to the coarser spectral and spatial resolution of other multispectral sensors.     

Estimation of aboveground grassland phytomass with a growth model using Landsat TM and climate data

Estimation of aboveground phytomass in meadow grasslands was carried out using multitemporal satellite data of fine resolution but of low frequency from Landsat TM observations in 1984-1990. We developed two growth models for the estimation of the first-cut yield, based on an exponential plant growth of which initial phytomass was determined using NDVI or TM2/TM3 on the Landsat observation date together with the effective cumulative temperature during the growth period after the observation date. Validations of these models with different data sets of Landsat TM in 1990-1994 indicated that the measured and estimated yields agreed well, suggesting the great potential of applying fine resolution satellite data coupled with a growth model to phytomass study, in spite of the low frequency of observation.     

Mapping montane tropical forest successional stage and land use with multi-date Landsat imagery

When mapping land cover with satellite imagery in montane tropical regions, varying illumination angles and ecological zones can obscure the differences between spectral responses of old-growth forest, secondary forest and agricultural lands. We used multi-date, Landsat Thematic Mapper (TM) imagery to map secondary forests, agricultural lands and old-growth forests in the Talamanca Mountain Range in southern Costa Rica. With stratification by illumination and ecological zone, the overall accuracy for this classification was 87% with a Kappa coefficient of 0.83. We also examined spectral responses to forest successional stage, ecological zone and aspect illumination for the TM Tasselled Cap indices, TM (2 x 6)/7, TM 4/5 and TM difference bands, and whether using digital data from multiple decades improved classification accuracy. Digital maps of ecological zones should be useful for large-scale mapping of land use and forest successional stage in complex montane regions such as those in Central America.     

Comparison of data gap-filling methods for Landsat ETM+ SLC-off imagery for monitoring forest degradation in a semi-deciduous tropical forest in Mexico

A number of methods to overcome the 2003 failure of the Landsat 7 Enhanced Thematic Mapper (ETM+) scan-line corrector (SLC) are compared in this article in a forest-monitoring application in the Yucatan Peninsula, Mexico. The objective of this comparison is to determine the best approach to accomplish SLC-off image gap-filling for the particular landscape in this region, and thereby provide continuity in the Landsat data sensor archive for forest-monitoring purposes. Four methods were tested: (1) local linear histogram matching (LLHM); (2) neighbourhood similar pixel interpolator (NSPI); (3) geostatistical neighbourhood similar pixel interpolator (GNSPI); and (4) weighted linear regression (WLR). All methods generated reasonable SLC-off gap-filling data that were visually consistent and could be employed in subsequent digital image analysis. Overall accuracy, kappa coefficients (κ), and quantity and allocation disagreement indices were used to evaluate unsupervised Iterative Self-Organizing Data Analysis (ISODATA) land-cover classification maps. In addition, Pearson correlation coefficients (r) and root mean squares of the error (RMSEs) were employed for estimates agreement with fractional land cover. The best results visually (overall accuracy > 85%, κ < 9%, quantity disagreement index < 5.5%, and allocation disagreement index < 12.5%) and statistically (r > 0.84 and RMSE < 7%) were obtained from the GNSPI method. These results suggest that the GNSPI method is suitable for routine use in reconstructing the imagery stack of Landsat ETM+ SLC-off gap-filled data for use in forest-monitoring applications in this type of heterogeneous landscape.     

Comparison of automatic thresholding methods for snow-cover mapping using Landsat TM imagery

Knowledge of snow cover is essential to understanding the global water and energy cycle. Thresholding the normalized difference snow index (NDSI) image is a method frequently used to map snow cover from remotely sensed data. However, the threshold is dependent on the scenario and needs to be determined accordingly. In this study, nine automatic thresholding methods were tested on the NDSI. Comparisons of the automatic thresholding methods, optimal threshold, and support vector machine (SVM) classification show that Otsu's and Nie's methods appear to be the most robust among the nine automatic thresholding methods, achieving comparable accuracies with the latter two approaches. In addition, NDSI from the digital number (DN) can be an efficient substitution for NDSI obtained from atmospherically or topographically corrected data, with similar accuracy.     

Accuracy of forest mapping based on Landsat TM data and a kNN-based method

A multi-source forest inventory (MSFI) method has been developed for use in the Norwegian National Forest Inventory (NFI). The method is based on a k-nearest neighbour rule and uses field plots from the NFI, land cover maps, and satellite image data from Landsat Thematic Mapper. The inventory method is used to produce maps of selected forest variables and to estimate the selected forest variables for large areas such as municipalities. In this study, focus has been on the qualitative variables ‘dominating species group’ and ‘development class’ because these variables are of central interest to forest managers. A mid-summer Landsat 5 TM scene was used as image data, and all NFI plots inside the scene were used as a reference dataset. The relationship between the spectral bands and the forest variables was analysed, and it was found that the levels of association were low. A leave-one-out method based on the reference dataset was used to estimate the pixel-level accuracies. They were found to be relatively low with 63% agreement for species groups. An independent control survey was available for a municipality and estimates from the MSFI were compared to it. The levels of error were quite high. It was concluded that the large area estimates were biased by the reference dataset.     

On the estimation of biomass of submerged vegetation using Landsat thematic mapper (TM) imagery: A case study of the Honghu Lake,PR China

The visible and near infrared bands of the Landsat thematic mapper (TM) were used in an empirical assessment of submerged vegetation biomass in Honghu Lake in the middle reaches of the Yangtse River, in the People's Republic of China. The method used here was based on eigenvector rotation of the four bands to enhance submerged vegetation biomass variations. Field measurement of spectral reflectance of submerged vegetation was taken for various biomass and vegetation types to determine the possibility of estimating submerged vegetation biomass using remote sensing. The locations of sample points were determined by global positioning system (GPS) and field biomass was obtained at the same time as the TM image. Regression analyses were performed between the principal components and biomass, and a marked linear relationship between submerged vegetation biomass and first two principal components (PC) was revealed. This was used to determine the total biomass of submerged vegetation.