Research Papers

Evaluating land use and aboveground biomass dynamics in an oil palm–dominated landscape in Borneo using optical remote sensing

[+] Author Affiliations
Minerva Singh

University of Oxford, School of Geography and the Environment, Oxford OX1 3QY, United Kingdom

University of Cambridge, Department of Plant Sciences, Forest Ecology and Conservation Group, Cambridge CB2 3EA, United Kingdom

Yadvinder Malhi

University of Oxford, Environmental Change Institute, Oxford University Centre for the Environment, Oxford OX1 3QY, United Kingdom

Shonil Bhagwat

University of Oxford, School of Geography and the Environment, Oxford OX1 3QY, United Kingdom

J. Appl. Remote Sens. 8(1), 083695 (Jan 10, 2014). doi:10.1117/1.JRS.8.083695
History: Received September 18, 2013; Revised November 22, 2013; Accepted December 10, 2013
Text Size: A A A

Open Access Open Access

Abstract.  The focus of this study is to assess the efficacy of using optical remote sensing (RS) in evaluating disparities in forest composition and aboveground biomass (AGB). The research was carried out in the East Sabah region, Malaysia, which constitutes a disturbance gradient ranging from pristine old growth forests to forests that have experienced varying levels of disturbances. Additionally, a significant proportion of the area consists of oil palm plantations. In accordance with local laws, riparian forest (RF) zones have been retained within oil palm plantations and other forest types. The RS imagery was used to assess forest stand structure and AGB. Band reflectance, vegetation indicators, and gray-level co-occurrence matrix (GLCM) consistency features were used as predictor variables in regression analysis. Results indicate that the spectral variables were limited in their effectiveness in differentiating between forest types and in calculating biomass. However, GLCM based variables illustrated strong correlations with the forest stand structures as well as with the biomass of the various forest types in the study area. The present study provides new insights into the efficacy of texture examination methods in differentiating between various land-use types (including small, isolated forest zones such as RFs) as well as their AGB stocks.

Figures in this Article

Over the past several decades, optical remote sensing (RS) data has been widely used for evaluating changes in land-use dynamics,1,2 biomass dynamics,3 and monitoring of degradation.4 At regional and national levels, vegetation indices have been commonly used to evaluate biophysical properties of forests. There are numerous examples of using optical RS for the study of land use, biomass, and degradation dynamics in tropical forests across the globe. The following sections will try to synthesize some of these.

Tangki and Chappell5 have evaluated the biomass variation in a selectively logged forest concession of Sabah by comparing field aboveground biomass (AGB) data with Normalized Difference Vegetation Index (NDVI) measures obtained from Landsat thematic mapper (TM) data. Results indicated that the NDVI explains up to 58% of the variation in AGB. However, the use of variables derived from optical RS data has had limited success in Borneo in predicting biomass.6 Ten vegetation indices were used in conjunction with field data to generate predictive relation for biomass estimations for a number of tropical sites, and these found weak correlations with field-based AGB measures.7

In addition to spectral-based methods using RS data, texture-based methods have provided valuable insights into the variation of forest structure, biophysical properties, and biomass dynamics in tropical forests of both Malaysia and Thailand.8 Research by Cutler et al.8 and Kuplich et al.9 using texture indices derived from high-resolution radar data has shown strong associations with AGB measures. Texture analysis of optical RS data has also been successfully used to quantify structural properties of forests.10 Wijaya et al.11 used a number of Landsat-derived measures, including vegetation indices and texture variables, to examine spatial variations in AGB and forest stand parameters in Kalimantan. Their research has shown that the biomass and forest stand values have stronger associations with texture measures derived using gray-level co-occurrence matrices (GLCMs).11 These findings have been corroborated by 12 using SPOT 5 data to derive spectral parameters and texture variables. Texture variables showed stronger associations with forest stand parameters, such as basal area, bole volume, and canopy height, than with spectral parameters.

These examples illustrate that the spectral characteristics of optical data (such as vegetation indices) have been widely used for studying biomass dynamics, but their performance has been limited to explain forest biophysical parameters in dense tropical forests. Texture analysis (such as those derived from GLCM) is a promising alternative to characterize the variation in AGB and forest stand parameters. However, these techniques have so far not been tested on landscapes dominated by oil palm plantations containing isolated forests such as riparian forest (RF) zones.

Logging, deforestation, and conversion of land to oil palm plantation have increased forest fragmentation in Borneo. Given the extent of forest loss and logging, it is essential to evaluate the ability of remnant forests to provide AGB storage and retention of tree biodiversity, especially forest fragments and riparian margins. A riparian forest zone (RF) is defined as the land adjacent to streams and rivers. Malaysia has 189 river systems, out of which 78 are located in Sabah. Riparian zones are legally protected and should be maintained for all permanent water courses. The width of these buffer zones varies according to state laws.13 In the riparian zones of the study area, a width of 30 m on each side of the waterway has been retained in all forest types.14 A significant body of literature indicates that the presence of RFs in a landscape can assist in biodiversity conservation.15,16 However, despite legal protection, RF zones are vulnerable to illegal logging.17 Currently, there are no frameworks in place to identify RF zones at the landscape scale and to monitor the variation in AGB stocks. Optical RS imageries, such as Landsat and SPOT, are widely available for free or at low cost. The motivation behind the present study is to identify techniques (such as vegetation and texture indices) which could be applied to these widely available dataset and would allow for monitoring the impacts of logging and the identification of isolated forest zones in mixed landscapes such as those found in Borneo. It is anticipated that this type of research could inform RS-based monitoring of oil palm tropical forest landscapes and smaller isolated forests located within these landscapes.

The objectives of the present study are: (a) to examine whether it is possible to distinguish between old growth primary forests, RFs of different logging intensities, and oil palm plantations using the spectral characteristics of Landsat TM and SPOT 5 data, (b) to examine whether texture-based measures can provide insights into the spatial variation of forest structure and biomass dynamics across the different forest types, and (c) to generate biomass estimation models of the study area using RS data. To accomplish these objectives, several image-processing techniques are employed to enable the detection and mapping of land-cover patterns and to generate biomass estimation models for the different land-use types.

Study Area

The research was carried out at the study site of the Stability of Altered Forest Ecosystems (SAFE) Project located in the Yayasan Sabah Concession area14,18 in Sabah, Malaysia (Fig. 1). The area consists of a mixed landscape that includes areas of a twice logged forest virgin jungle reserve (VJR), oil palm plantations (OP), and a 7200 ha heavily logged area known as the experimental area (EA), which has been ear marked for the conversion to oil palm plantation beginning in December 2011. The study area extends to the Maliau Conservation Basin (116.87 E, 4.82 N). It has been proposed to retain circular arrangements of forests (of varying forest cover). A number of RF sites will also be retained both in the proposed and in the existing oil palm plantations. In addition, RF zones will be retained in other land-use types.

Graphic Jump LocationF1 :

(a) Location of Sabah in Borneo. (b) Location of the study area within Sabah.18

Based on ground-truth data provided by SAFE, several different land-use types are present in the study area.14 These include old growth primary forests (OG), disturbed forests (which have been exposed to slight anthropogenic disturbance and clearance), VJR (which have been logged once at most), twice logged forests (LF), heavily logged forests (EA; which have undergone three rounds of logging and are in a state of severe degradation), RFs, and oil palm plantations (OP). Details of the structures of the various forest types in study area are presented in Table 1. These different forest types have significant differences in their stand and canopy structure (see Tables 2 and 3).

Table Grahic Jump Location
Table 1Structure of different forest types present in the study area (photographs taken at the SAFE project site, 2011).
Table Grahic Jump Location
Table 2Aboveground forest parameters across the riparian forests (RF).
Table Grahic Jump Location
Table 3Aboveground forest parameters across the non-RF (NRF) zones.

Data for forest mensuration parameters across different land-use types [diameter at breast height (DBH)] was obtained using 193 (25×25m2) vegetation plots spread across the different forest types in the study area.14 These vegetation plots have been set up following a fractal design.20 In the RFs of each of the different forest types, six plots (10×50m2) were set up in three riparian zones each. For each forest type, 18 plots were established from September to December 2011. Plot establishment followed a random stratified sampling strategy. The plots are located in all forest types present in the study area, but within each study area, they have been located randomly. In riparian plots, both DBH and tree height (H) were measured. Table 2 presents the forest mensuration variables measured from both riparian and non-riparian plots of different land-use types.

Tree heights for non-riparian plots were estimated using a DBH-height relationship provided by 21. The AGB of the trees in both riparian and non-riparian zones was calculated using the biomass equation recommended by 22Display Formula

AbovegroundBiomass(AGB)=0.0776×(ρ×DBH2×H)0.94,(1)
where ρ refers to the wood-specific gravity. Given the difference in the physiology of oil palm trees and regular forest trees, specific biomass equations are needed to calculate the AGB of oil palm plantations. The AGB of oil palm trees was calculated using the biomass equation recommended by 23Display Formula
(AGB)Trunk=100×π×(r×z)2×h×ρ,(2)
where r is the radius of the trunk (in cm) without frond bases, z is the ratio of the trunk diameter below the frond bases to the measured diameter above the frond bases (estimated to be 0.777 from the sampled trunks), and h is the height of the trunk (in m) to the base of the fronds. The trunk density ρ (in kgm3) is defined as follows: Display Formula
ρ=(0.0076x+0.083)/100,(3)
where x is the age of the oil palm plantation.

Satellite Data

Landsat TM and SPOT 5 from 2009 were used in this study. Landsat TM data has a spatial resolution of 30 m and has six spectral bands. The SPOT 5 data consisted of 10-m multispectral bands [green, red, near-infrared (NIR)] along with a 20-m short-wave infrared (SWIR) band, which was resampled to 10 m. Appropriate subsets covering the entire study area were clipped from the imagery data. These image subsets were co-registered and geo-referenced using appropriate ground-control points. A digital elevation model (DEM) data was obtained from the Shuttle Radar Topography Mission data.24 The DEM (originally at 90 m) was used to orthorectify the Landsat TM data.

Image Processing

A number of preprocessing and image-processing techniques were applied on both the Landsat TM and SPOT 5 datasets. Image preprocessing techniques included atmospheric and haze corrections. Unsupervised classification was carried out to isolate nonrelevant features (such as clouds and shadows) in both the satellite datasets. Supervised classification was carried out using the ground survey data collected while doing field survey. The maximum likelihood classification algorithm was implemented in ENVI(R) image processing software.

Vegetation Indices

The satellite data was converted to reflectance values, and these values were been used to calculate three vegetation indices: the Normalized Difference Vegetation Index (NDVI), the Soil-Adjusted Vegetation Index (SAVI), and the Normalized Difference Infrared Index (NDII). The NDVI is one of the most commonly used vegetation indices and is an indicator of the green vegetation in the study area. It is determined as follows: Display Formula

NDVI=(ρB4ρB3)(ρB4+ρB3)(4)
Display Formula
SAVI=1.5(ρB4ρB3)(ρB4+ρB3+0.5),(5)
where SAVI is a modification of NDVI and corrects for the influence of soil brightness when the vegetation cover is low.25Display Formula
NDII5=(ρB4ρB5)(ρB4+ρB5),(6)
where NDII5 is the vegetation index that makes use of the NIR and SWIR band information to distinguish between different land-use and forest type classes, as recommended by Souza et al.26

The NDII5 helps distinguish forest disturbances on the basis of difference in water content. Values derived using spectral characteristics, such as vegetation indices, and band reflectances of SPOT 5 data were used to distinguish between different land-use types. Individual pixels from different classes were randomly selected for each of the bands. Tukey tests were applied to the values obtained from the different land-use types in order to evaluate the class separability. Tukey tests consist of a multicomparison between the population means of the different land-use types.26

Texture Measures: GLCM

Texture analysis of an image can be carried out using statistical, spectral, and spatial techniques.27 In statistical-based texture analysis methods, information is obtained by measuring the spatial variation in an image’s tonal values.28 Image texture measures can be classified into two categories. The first category is occurrence, also known as first-order statistics. This relates to the frequency of tonal values in a specified neighborhood around each pixel.29 This does not take spatial relationships among pixels into account. The second category is co-occurrence, also known as second-order statistics. This measures the frequency of associations between brightness value pairs within a given area.28,30 The statistical texture indices derived using the GLCM algorithms are tabulated in Table 4.

Table Grahic Jump Location
Table 4Texture variables derived from GLCM (from 19).

The texture variables are calculated on the basis of the red band.31 According to GLCM analysis carried out by 19, the texture variables derived from the red band of optical RS data are the best explanatory variables for vegetation attributes. These texture-derived variables, especially mean, standard deviation, contrast, and entropy, have the ability to provide important information about the forest stand parameters and biomass dynamics and the variation and heterogeneity in these characteristics.

Land-Cover Map of the Study Area

The land-cover map of the study area (Fig. 2) was created using Landsat TM data and clearly delineates the areas of different logging intensities in the study area along with the areas of primary forest, oil palm plantations, and logged forests.

Graphic Jump LocationF2 :

Landsat TM-derived land-cover map of the study area. Gray lines represent the SAFE and Maliau basin conservation area (MBCA) boundaries.

While most of the field survey data were used for generating training pixels or regions of interest (ROIs), 30% of the field survey data was set aside for validation. A confusion matrix was used to describe the accuracy of the classification.32 The classified Landsat image had a classification accuracy of 77% and a kappa coefficient of 0.45.

Distinguishing Between Different Land-Use Classes Using Reflectance and Vegetation Indices
Efficacy of SPOT 5-based RS data in distinguishing between different land-use types

Band reflectance and spectral-based vegetation indices were used to distinguish between the different forest types in the study area (Fig. 3). All forest-use types show substantial differences in the shortwave band of the SPOT data. While other bands, such as infrared and red bands, can distinguish between the major forest classes, they cannot distinguish non-riparian zones from riparian zones and twice logged forests from heavily logged forests. All three bands can distinguish oil palm plantations from other forest types. However, the different logged forests show overlap in the red and NIR bands.

Graphic Jump LocationF3 :

SPOT-derived vegetation indices for different land-use types. OG: old growth forests; OP: oil palm plantations; RF: riparian forests; VJR: virgin jungle reserve; EA: experimental area; and LF: twice logged forests.

Vegetation indices derived from the higher resolution SPOT data were used to distinguish between different forest types present in the study area. The values for NDVI and NDII5 decrease sharply from primary forests to once logged forests and then increase from once logged to twice logged to heavily logged forests. The NDVI values based on SPOT 5 are able to distinguish between different forest types and forest transition systems in the study area. The NDVI values differ significantly between old growth primary forests, once logged forests, twice logged forests, and heavily logged forests. The oil palm plantations can also be distinguished from all other land-use types. However, NDVI was not able to distinguish between the RFs and the pristine, once logged, and twice logged forests. The NDVI values differ significantly between heavily logged forests and oil palm plantations. However, SAVI is not as effective in distinguishing between different forest types. The SAVI values varied significantly between the old growth pristine forests and once or slightly logged forests and old growth forests and heavily logged forests. However, SAVI values did not differ significantly between the once logged, twice logged, and heavily logged forests. The SAVI values (like NDVI) differ significantly between oil palm plantations and all other forest transition systems. The NDII5 is not able to distinguish oil palm plantations from other land-use types.

Efficacy of Landsat TM-based RS data in distinguishing between different land-use types in the area

Based on Landsat TM data, NDVI values differ significantly between oil palm plantations and other land-use types. However, NDVI values cannot distinguish between forests that have undergone different logging intensities. These results can be attributed to the saturation of the NDVI value at high biomass levels.33,34 This limits their ability to predict biomass in tropical forests such as those in Borneo.

Correlation Between VIs and Field Measures of AGB

Predictive models of biomass were generated using AGB data collected using field surveys in conjunction with vegetation indices and band reflectances derived from Landsat TM and SPOT 5 data. Reflectance values of the green (band 1), red (band 2), NIR (band 3), and SWIR (band 4) bands of the SPOT 5 data were correlated with field-based AGB data to derive biomass estimation models. The biomass estimation model is as follows: Display Formula

(Biomass)SPOT=106.37*(Band1)33.72*(Band2)+124.33*(Band3)+40.73*(Band4)130.71.(7)

A strong association was found between field-based AGB values and those derived using all of the SPOT bands. However, the SPOT biomass model in Eq. (7) tends to underestimate the AGB values. While most of the field-based AGB data were used to establish the coefficients in Eq. (7), a sample of the data was set aside for validation purposes. The results of validation are presented in Fig. 4. A biomass estimation model using the red, NIR, and SWIR bands also displays a high correlation (R2=0.8512, p<0.001) with field-based AGB values. Single-band biomass estimation models consisting of red (band 3) and NIR (band 4) bands have an R2 value of 0.80 and 0.833, respectively (p<0.001).

Graphic Jump LocationF4 :

Actual and expected regression models between AGB values derived using SPOT imagery and field observations. All four SPOT bands (green, red, NIR, and SWIR) were used for deriving SPOT-derived AGB.

Based on Landsat data, NDVI has a very weak correlation (R2=0.0562, p<0.01) with the field-based AGB values. Other vegetation indices, such as SAVI, also show a very weak correlation with the field-based AGB values. On the other hand, the reflectance values of band 3 show a strong association with the field-based AGB values (R2=0.85, p<0.01), whereas the correlation is weaker for band 4 (R2=0.59, p<0.01). The biomass estimation model based on band reflectance value can be expressed as Display Formula

Ba=2175.4(Band4)+2628.3(8)
Display Formula
Ba=2776(Band3)+2817.4,(9)
where Ba is the mean biomass of study area (in Mg/ha), which has been obtained from the reflectance values of bands 3 and 4.

Use of Texture-Based Variables in Evaluating Forest Stand Parameters and Biomass Dynamics

High-resolution satellite data is better suited for the calculation of texture indices than coarse- and medium-resolution data.28,35 For this reason, comparatively higher resolution SPOT data is used for obtaining the texture-based variables. Two vegetation attributes are employed: AGB and basal area. A statistical examination of the AGB values in the different land-use types indicates that they are influenced by the disturbance due to several rounds of logging and oil palm cultivation. Regression models were generated to describe the relationships between the vegetation attributes and the texture variables described in Table 4. The texture variables showed strong but varying associations with AGB and basal area of the different forest types (Table 5).

Table Grahic Jump Location
Table 5Ability of texture-based linear regression models to explain variation in vegetation attributes.

However, the results for oil palm plantations were not significant. The texture-based biomass estimation models (which were built using the variables mean, variance, and contrast) underestimated the AGB values for all the land-use types considered. The validation results for some of the land-use types are presented in Fig. 5.

Graphic Jump LocationF5 :

(a): Comparison between Field and Texture derived AGB for EA (b) Comparison between Field and Texture derived AGB for LF (c) Comparison between Field and Texture derived AGB for RF. In all three cases, it can be seen that while there is strong relationship between field and texture AGB, texture based AGB values underestimate field AGB values for all three land use types under consideration.

The present study compares a number of RS variables that can be used for the evaluation of land use and biomass dynamics. The research has accomplished two out of its three objectives. The first objective was to use vegetation indices and band reflectance values derived from Landsat and SPOT data to distinguish between different forest types in the study area. The results indicate that the vegetation indices derived from Landsat TM and SPOT 5 data have a limited potential to distinguish between the different land-use types. Vegetation indices derived from Landsat data (such as NDVI) can distinguish between major land-use types, such as old growth pristine forests and oil palm plantations, but not between forests with different logging intensities. The NDVI derived from higher resolution SPOT data can distinguish between pristine forests, forests of different logging intensities, and oil palm plantations. However, other SPOT-based vegetation indices also cannot distinguish between forests of different logging intensities. In addition, vegetation indices derived from both Landsat TM and SPOT 5 (such as NDVI) are not strongly correlated with the field-based AGB values. For the particular study area, the vegetation indices are insufficient for distinguishing between forests of different logging intensities and for generating biomass estimation models for individual land-use types. This confirms the growing body of literature which indicates that the use of spectral information is limited due to its saturation at high biomass values, thereby restricting its utility.7,11,12,26,31,36 As opposed to the vegetation indices, reflectance values of bands 3 and 4 of Landsat are strongly correlated with the field-based AGB values. Similarly, strong associations have been observed in a logged forest-pristine forest study site (similar to ours) located in the Danum Valley Conservation Area, which is close to the SAFE study area.5

The second and third objectives were to use texture-based variables to examine the structure and biomass parameters of different forest types in the study area and to evaluate the possibility of developing biomass estimation models for the different forest types. The use of the GLCM-based texture analysis methods provides insight into the biomass and structural dynamics across different land-use types. The texture analysis also helps identify that the different land-use types, including forests having undergone different logging rotations, vary in terms of their structure. Optical RS data can be used to identify (and delineate) small fragmented areas such as riparian zones and forests having different logging intensities. Statistically based texture variables have been widely cited in the literature for their strong association with variation in forest structure, forest stand parameters, and biomass dynamics.8,11,31 This accomplishes the second objective of the research. Table 5 shows that the regression models using texture variables and field-based AGB values vary for different forest types in terms of their strength. It can therefore be concluded that the use of texture analysis opens up the possibility of developing different biomass estimation models for different land-use types. This accomplished the third objective of the research.

In summary, the present study confirms the findings of previous research that vegetation indices derived from optical RS data have limited potential to distinguish between different forest types compared with texture-derived measures. The present study has determined that AGB is underestimated by both spectral and textural variables. A possible reason for this could be biomass saturation. This suggests that optical RS biomass estimation models need to be interpreted with care. Most importantly, however, the present study has shown that it is possible to identify isolated RF zones and to show that their structure and biomass dynamics differ from surrounding logged forests. To the best of our knowledge, only one field-based study has previously been carried out on the RFs of Malaysia.13 Arguably, texture-based analyses could be applied for studying the biomass and structure dynamics of isolated forest fragments.

The RS data offers the potential to study various forest properties including their structure, carbon dynamics, and assessment of degradation. The present study has shown a number of different techniques that may be used to overcome the shortcomings of optical RS data. Most importantly, the efficacy of these different techniques in examining the structural and biomass dynamics of tropical forests has been demonstrated, especially in the case of mixed land-use types. This research can be applied to the management of lowland Dipterocarp forests and to the evaluation of their carbon stocks. The present study has established that the band reflectance and texture measures derived from GLCM can be used to generate biomass estimation models that correlate strongly with the field-based AGB. The biomass estimation models derived from texture variables can be applied to similar mixed land-use types. Application of these biomass estimate models could allow for the assessment of timber stocks and for the assessment of the ecological status in terms of forest recovery and variation in carbon stocks. Most importantly, the use of biomass estimate models makes it possible to examine both the spatial variation in forest canopy structure and carbon stocks across a variety of different land-use types and disturbance gradients. Texture analysis of optical RS images is a very promising technique for a study area like Borneo, where traditional vegetation indices-based analysis may be limited due to data saturation for sites with high biomass density or sites having complex forest stand structure.37

Achard  F. et al., “Determination of deforestation rates of the world’s humid tropical forests,” Science. 297, (5583 ), 999 –1002 (2002), CrossRef. 0036-8075 
Hansen  M. C. et al., “A method for integrating MODIS and Landsat data for systematic monitoring of forest cover and change in the Congo Basin,” Remote Sens. Environ.. 112, (5 ), 2495 –2513 (2008), CrossRef. 0034-4257 
Asner  G. P. et al., “High-resolution forest carbon stocks and emissions in the Amazon,” Proc. Natl. Acad. Sci. U. S. A.. 107, (38 ), 16738 –16742 (2010), CrossRef. 0027-8424 
Margono  B. A. et al., “Mapping and monitoring deforestation and forest degradation in Sumatra (Indonesia) using Landsat time series data sets from 1990 to 2010,” Environ. Res. Lett.. 7, (3 ), 034010  (2012), CrossRef. 1748-9326 
Tangki  H., Chappell  N. A., “Biomass variation across selectively logged forest within a 225-km2 region of Borneo and its prediction by Landsat TM,” For. Ecol. Manage.. 256, (11 ), 1960 –1970 (2008), CrossRef. 0378-1127 
Foody  G. M. et al., “Mapping the biomass of Bornean tropical rain forest from remotely sensed data,” Global Ecol. Biogeogr.. 10, (4 ), 379 –387 (2001), CrossRef. 1466-8238 
Foody  G. M., Boyd  D. S., Cutler  M. E., “Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions,” Remote Sens. Environ.. 85, (4 ), 463 –474 (2003), CrossRef. 0034-4257 
Cutler  M. E. J. et al., “Estimating tropical forest biomass with a combination of SAR image texture and Landsat TM data: an assessment of predictions between regions,” ISPRS J. Photogramm. Remote Sens.. 70, , 66 –77 (2012), CrossRef. 0924-2716 
Kuplich  T. M., Curran  P. J., Atkinson  P. M., “Relating SAR image texture to the biomass of regenerating tropical forests,” Int. J. Remote Sens.. 26, (21 ), 4829 –4854 (2005), CrossRef. 0143-1161 
Kayitakire  F., Hamel  C., Defourny  P., “Retrieving forest structure variables based on image texture analysis and IKONOS-2 imagery,” Remote Sens. Environ.. 102, (3 ), 390 –401 (2006), CrossRef. 0034-4257 
Wijaya  A. et al., “Improved strategy for estimating stem volume and forest biomass using moderate resolution remote sensing data and GIS,” J. For. Res.. 21, (1 ), 1 –12 (2010), CrossRef. 1007-662X 
Castillo-Santiago  M. A., Ricker  M., de Jong  B. H. J., “Estimation of tropical forest structure from SPOT-5 satellite images,” Int. J. Remote Sens.. 31, (10 ), 2767 –2782 (2010), CrossRef. 0143-1161 
Azliza  M. et al., “Characterization of riparian plant community in lowland forest of peninsular Malaysia,” Int. J. Bot.. 8, (4 ), 181 –191 (2012), CrossRef. 1811-9700 
“Stability of Altered Forest Ecosystems SAFE,” http://www.safeproject.net/ (2011).
Lees  A. C., Peres  C. A., “Conservation value of remnant riparian forest corridors of varying quality for Amazonian birds and mammals,” Conserv. Biol.. 22, (2 ), 439 –449 (2008), CrossRef. 0888-8892 
MacDonald  M. A., “The role of corridors in biodiversity conservation in production forest landscapes: a literature review,” Tasforests. 14, , 41 –52 (2003).
Ibbotson  R., “Personal email communication about Sabah’s logging history” (2013).
Ewers  R. M. et al., “A large-scale forest fragmentation experiment: the Stability of Altered Forest Ecosystems Project,” Philos. Trans. R. Soc. B. 366, (1582 ), 3292 –3302 (2011), CrossRef. 0962-8436 
Gallardo-Cruz  J. A. et al., “Predicting tropical dry forest successional attributes from space: is the key hidden in image texture?” PLoS One. 7, (2 ), e30506  (2012), CrossRef. 1932-6203 
Marsh  C. J., Ewers  R. M., “A fractal-based sampling design for ecological surveys quantifying β-diversity,” Methods Ecol. Evol.. 4, (1 ), 63 –72 (2013), CrossRef.
Morel  A. C. et al., “Estimating aboveground biomass in forest and oil palm plantation in Sabah, Malaysian Borneo using ALOS PALSAR data,” For. Ecol. Manage.. 262, (9 ), 1786 –1798 (2011), CrossRef. 0378-1127 
Chave  J. et al., “Tree allometry and improved estimation of carbon stocks and balance in tropical forests,” Oecologia. 145, (1 ), 87 –99 (2005), CrossRef. 0029-8549 
Morel  A. C., Fisher  J. B., Malhi  Y., “Evaluating the potential to monitor aboveground biomass in forest and oil palm in Sabah, Malaysia, for 2000–2008 with Landsat ETM+ and ALOS-PALSAR,” Int. J. Remote Sens.. 33, (11 ), 3614 –3639 (2012), CrossRef. 0143-1161 
USGS, “US Geological Survey,” http://srtm.usgs.gov/ (2011).
Huete  A. R., “A soil-adjusted vegetation index (SAVI),” Remote Sens. Environ.. 25, (3 ), 295 –309 (1988), CrossRef. 0034-4257 
Souza  C. M.  Jr., Roberts  D. A., Monteiro  A., “Multitemporal analysis of degraded forests in the southern Brazilian Amazon,” Earth Interact.. 9, (19 ), 1 –25 (2005), CrossRef. 1087-3562 
Ploton  P., “Analyzing canopy heterogeneity of the tropical forests by texture analysis of very-high resolution images—a case study in the western ghats of 76 remote sensing of biomass—principles and applications. India,” Pondy Pap. Ecol.. 10, , 1 –71 (2010).
Jones  H. G., Vaughan  R. A., Remote Sensing of Vegetation: Principles, Techniques, and Applications. ,  Oxford University Press ,  New York  (2010).
St-Louis  V. et al., “High-resolution image texture as a predictor of bird species richness,” Remote Sens. Environ.. 105, (4 ), 299 –312 (2006), CrossRef. 0034-4257 
Estes  L. D. et al., “Habitat selection by a rare forest antelope: a multi-scale approach combining field data and imagery from three sensors,” Remote Sens. Environ.. 112, (5 ), 2033 –2050 (2008), CrossRef. 0034-4257 
Eckert  S., “Improved forest biomass and carbon estimations using texture measures from WorldView-2 satellite data,” Remote Sens.. 4, (4 ), 810 –829 (2012), CrossRef. 2072-4292 
Foody  G. M., “Status of land cover classification accuracy assessment,” Remote Sens. Environ.. 80, (1 ), 185 –201 (2002), CrossRef. 0034-4257 
Basuki  T. M. et al., “The potential of spectral mixture analysis to improve the estimation accuracy of tropical forest biomass,” Geocarto Int.. 27, (4 ), 329 –345 (2012), CrossRef. 1010-6049 
Sader  S. A. et al., “Tropical forest biomass and successional age class relationships to a vegetation index derived from Landsat TM data,” Remote Sens. Environ.. 28, , 143 –198 (1989), CrossRef. 0034-4257 
Beguet  B. et al., “Retrieving forest structure variables from very high resolution satellite images using an automatic method,” in  ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume I-7, 2012 XXII ISPRS Congress ,  Melbourne, Australia .
Lu  D. et al., “Application of spectral mixture analysis to Amazonian land-use and land-cover classification,” Int. J. Remote Sens.. 25, (23 ), 5345 –5358 (2004), CrossRef. 0143-1161 
Lu  D. et al., “Aboveground forest biomass estimation with Landsat and lidar data and uncertainty analysis of the estimates,” Int. J. For. Res.. 2012, , 436537  (2012), CrossRef. 0300-9831 

Minerva Singh is a graduate student at the Department of Plant Sciences, University of Cambridge, with Dr. David Coomes.

Yadvinder Malhi is a professor of ecosystem science at the School of Geography and the Environment. He is interested on interactions between forest ecosystems and the global atmosphere, with a particular focus on their role in global carbon, energy, and water cycles, and in understanding how the ecology of natural ecosystems may be shifted in response to global atmospheric change.

Shonil Bhagwat is the senior visiting research associate in the School of Geography and the Environment at University of Oxford. His research topics include conservation beyond protected-area boundaries, cultural, and spiritual values in conservation, and climate change and conservation.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Minerva Singh ; Yadvinder Malhi and Shonil Bhagwat
"Evaluating land use and aboveground biomass dynamics in an oil palm–dominated landscape in Borneo using optical remote sensing", J. Appl. Remote Sens. 8(1), 083695 (Jan 10, 2014). ; http://dx.doi.org/10.1117/1.JRS.8.083695


Figures

Graphic Jump LocationF1 :

(a) Location of Sabah in Borneo. (b) Location of the study area within Sabah.18

Graphic Jump LocationF2 :

Landsat TM-derived land-cover map of the study area. Gray lines represent the SAFE and Maliau basin conservation area (MBCA) boundaries.

Graphic Jump LocationF3 :

SPOT-derived vegetation indices for different land-use types. OG: old growth forests; OP: oil palm plantations; RF: riparian forests; VJR: virgin jungle reserve; EA: experimental area; and LF: twice logged forests.

Graphic Jump LocationF4 :

Actual and expected regression models between AGB values derived using SPOT imagery and field observations. All four SPOT bands (green, red, NIR, and SWIR) were used for deriving SPOT-derived AGB.

Graphic Jump LocationF5 :

(a): Comparison between Field and Texture derived AGB for EA (b) Comparison between Field and Texture derived AGB for LF (c) Comparison between Field and Texture derived AGB for RF. In all three cases, it can be seen that while there is strong relationship between field and texture AGB, texture based AGB values underestimate field AGB values for all three land use types under consideration.

Tables

Table Grahic Jump Location
Table 1Structure of different forest types present in the study area (photographs taken at the SAFE project site, 2011).
Table Grahic Jump Location
Table 2Aboveground forest parameters across the riparian forests (RF).
Table Grahic Jump Location
Table 3Aboveground forest parameters across the non-RF (NRF) zones.
Table Grahic Jump Location
Table 4Texture variables derived from GLCM (from 19).
Table Grahic Jump Location
Table 5Ability of texture-based linear regression models to explain variation in vegetation attributes.

References

Achard  F. et al., “Determination of deforestation rates of the world’s humid tropical forests,” Science. 297, (5583 ), 999 –1002 (2002), CrossRef. 0036-8075 
Hansen  M. C. et al., “A method for integrating MODIS and Landsat data for systematic monitoring of forest cover and change in the Congo Basin,” Remote Sens. Environ.. 112, (5 ), 2495 –2513 (2008), CrossRef. 0034-4257 
Asner  G. P. et al., “High-resolution forest carbon stocks and emissions in the Amazon,” Proc. Natl. Acad. Sci. U. S. A.. 107, (38 ), 16738 –16742 (2010), CrossRef. 0027-8424 
Margono  B. A. et al., “Mapping and monitoring deforestation and forest degradation in Sumatra (Indonesia) using Landsat time series data sets from 1990 to 2010,” Environ. Res. Lett.. 7, (3 ), 034010  (2012), CrossRef. 1748-9326 
Tangki  H., Chappell  N. A., “Biomass variation across selectively logged forest within a 225-km2 region of Borneo and its prediction by Landsat TM,” For. Ecol. Manage.. 256, (11 ), 1960 –1970 (2008), CrossRef. 0378-1127 
Foody  G. M. et al., “Mapping the biomass of Bornean tropical rain forest from remotely sensed data,” Global Ecol. Biogeogr.. 10, (4 ), 379 –387 (2001), CrossRef. 1466-8238 
Foody  G. M., Boyd  D. S., Cutler  M. E., “Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions,” Remote Sens. Environ.. 85, (4 ), 463 –474 (2003), CrossRef. 0034-4257 
Cutler  M. E. J. et al., “Estimating tropical forest biomass with a combination of SAR image texture and Landsat TM data: an assessment of predictions between regions,” ISPRS J. Photogramm. Remote Sens.. 70, , 66 –77 (2012), CrossRef. 0924-2716 
Kuplich  T. M., Curran  P. J., Atkinson  P. M., “Relating SAR image texture to the biomass of regenerating tropical forests,” Int. J. Remote Sens.. 26, (21 ), 4829 –4854 (2005), CrossRef. 0143-1161 
Kayitakire  F., Hamel  C., Defourny  P., “Retrieving forest structure variables based on image texture analysis and IKONOS-2 imagery,” Remote Sens. Environ.. 102, (3 ), 390 –401 (2006), CrossRef. 0034-4257 
Wijaya  A. et al., “Improved strategy for estimating stem volume and forest biomass using moderate resolution remote sensing data and GIS,” J. For. Res.. 21, (1 ), 1 –12 (2010), CrossRef. 1007-662X 
Castillo-Santiago  M. A., Ricker  M., de Jong  B. H. J., “Estimation of tropical forest structure from SPOT-5 satellite images,” Int. J. Remote Sens.. 31, (10 ), 2767 –2782 (2010), CrossRef. 0143-1161 
Azliza  M. et al., “Characterization of riparian plant community in lowland forest of peninsular Malaysia,” Int. J. Bot.. 8, (4 ), 181 –191 (2012), CrossRef. 1811-9700 
“Stability of Altered Forest Ecosystems SAFE,” http://www.safeproject.net/ (2011).
Lees  A. C., Peres  C. A., “Conservation value of remnant riparian forest corridors of varying quality for Amazonian birds and mammals,” Conserv. Biol.. 22, (2 ), 439 –449 (2008), CrossRef. 0888-8892 
MacDonald  M. A., “The role of corridors in biodiversity conservation in production forest landscapes: a literature review,” Tasforests. 14, , 41 –52 (2003).
Ibbotson  R., “Personal email communication about Sabah’s logging history” (2013).
Ewers  R. M. et al., “A large-scale forest fragmentation experiment: the Stability of Altered Forest Ecosystems Project,” Philos. Trans. R. Soc. B. 366, (1582 ), 3292 –3302 (2011), CrossRef. 0962-8436 
Gallardo-Cruz  J. A. et al., “Predicting tropical dry forest successional attributes from space: is the key hidden in image texture?” PLoS One. 7, (2 ), e30506  (2012), CrossRef. 1932-6203 
Marsh  C. J., Ewers  R. M., “A fractal-based sampling design for ecological surveys quantifying β-diversity,” Methods Ecol. Evol.. 4, (1 ), 63 –72 (2013), CrossRef.
Morel  A. C. et al., “Estimating aboveground biomass in forest and oil palm plantation in Sabah, Malaysian Borneo using ALOS PALSAR data,” For. Ecol. Manage.. 262, (9 ), 1786 –1798 (2011), CrossRef. 0378-1127 
Chave  J. et al., “Tree allometry and improved estimation of carbon stocks and balance in tropical forests,” Oecologia. 145, (1 ), 87 –99 (2005), CrossRef. 0029-8549 
Morel  A. C., Fisher  J. B., Malhi  Y., “Evaluating the potential to monitor aboveground biomass in forest and oil palm in Sabah, Malaysia, for 2000–2008 with Landsat ETM+ and ALOS-PALSAR,” Int. J. Remote Sens.. 33, (11 ), 3614 –3639 (2012), CrossRef. 0143-1161 
USGS, “US Geological Survey,” http://srtm.usgs.gov/ (2011).
Huete  A. R., “A soil-adjusted vegetation index (SAVI),” Remote Sens. Environ.. 25, (3 ), 295 –309 (1988), CrossRef. 0034-4257 
Souza  C. M.  Jr., Roberts  D. A., Monteiro  A., “Multitemporal analysis of degraded forests in the southern Brazilian Amazon,” Earth Interact.. 9, (19 ), 1 –25 (2005), CrossRef. 1087-3562 
Ploton  P., “Analyzing canopy heterogeneity of the tropical forests by texture analysis of very-high resolution images—a case study in the western ghats of 76 remote sensing of biomass—principles and applications. India,” Pondy Pap. Ecol.. 10, , 1 –71 (2010).
Jones  H. G., Vaughan  R. A., Remote Sensing of Vegetation: Principles, Techniques, and Applications. ,  Oxford University Press ,  New York  (2010).
St-Louis  V. et al., “High-resolution image texture as a predictor of bird species richness,” Remote Sens. Environ.. 105, (4 ), 299 –312 (2006), CrossRef. 0034-4257 
Estes  L. D. et al., “Habitat selection by a rare forest antelope: a multi-scale approach combining field data and imagery from three sensors,” Remote Sens. Environ.. 112, (5 ), 2033 –2050 (2008), CrossRef. 0034-4257 
Eckert  S., “Improved forest biomass and carbon estimations using texture measures from WorldView-2 satellite data,” Remote Sens.. 4, (4 ), 810 –829 (2012), CrossRef. 2072-4292 
Foody  G. M., “Status of land cover classification accuracy assessment,” Remote Sens. Environ.. 80, (1 ), 185 –201 (2002), CrossRef. 0034-4257 
Basuki  T. M. et al., “The potential of spectral mixture analysis to improve the estimation accuracy of tropical forest biomass,” Geocarto Int.. 27, (4 ), 329 –345 (2012), CrossRef. 1010-6049 
Sader  S. A. et al., “Tropical forest biomass and successional age class relationships to a vegetation index derived from Landsat TM data,” Remote Sens. Environ.. 28, , 143 –198 (1989), CrossRef. 0034-4257 
Beguet  B. et al., “Retrieving forest structure variables from very high resolution satellite images using an automatic method,” in  ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume I-7, 2012 XXII ISPRS Congress ,  Melbourne, Australia .
Lu  D. et al., “Application of spectral mixture analysis to Amazonian land-use and land-cover classification,” Int. J. Remote Sens.. 25, (23 ), 5345 –5358 (2004), CrossRef. 0143-1161 
Lu  D. et al., “Aboveground forest biomass estimation with Landsat and lidar data and uncertainty analysis of the estimates,” Int. J. For. Res.. 2012, , 436537  (2012), CrossRef. 0300-9831 

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.