10 January 2014 Evaluating land use and aboveground biomass dynamics in an oil palm–dominated landscape in Borneo using optical remote sensing
Author Affiliations +
J. of Applied Remote Sensing, 8(1), 083695 (2014). doi:10.1117/1.JRS.8.083695
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.
Singh, Malhi, and Bhagwat: Evaluating land use and aboveground biomass dynamics in an oil palm–dominated landscape in Borneo using optical remote sensing

1.

Introduction

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 Ref. 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.

2.

Rationale and Objectives

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.

3.

Methodology

3.1.

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.

Fig. 1

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

JARS_8_1_083695_f001.png

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 1

Structure of different forest types present in the study area (photographs taken at the SAFE project site, 2011).

Land-use categoriesIllustrative photograph
Old Growth Forests (OG) These are pristine lowland forests located in the MBCA. The slope of these forest locations does not exceed 20 deg. These forests are characterized by a closed canopy and the presence of very large Dipterocarp trees. Some of the trees have heights of >50  m and diameter at breast height (DBH) of >80  cm. Tree species survey was carried out and it revealed that these forests also host IUCN red listed Endangered Dipterocarp species such as Shorea johorensis. These forests also have a very thick understory
Oil Palm Plantations (OP) The SAFE project site contains oil palm plantations with ages ranging from 5- to 8-years old. Oil palm plantations are characterized by a homogenous canopy structure and negligible tree species diversity. The picture shows the aerial view of the oil palm plantations in Sabah19
Heavily logged forests (EA) These are highly degraded forests, which have undergone several rounds of logging. These are characterized by the presence of large open areas and broken canopy. Virtually no large trees are observed, and the area is dominated by small successional trees and ginger shrubs
Riparian forests (RF) These are forest zones adjacent to rivers and streams. These are dominated by thin and tall vegetations and trees of the Maccaranga genera. In some of the riparian zones, the evidence of illegal logging has been observed based on to the presence of logging tracts. These forests are relatively intact compared with the surrounding heavily logged and twice logged forests

Table 2

Aboveground forest parameters across the riparian forests (RF).

(OG)RF(VJR)RF(LF)RF(EA)RF(OP)RF
Basal area (m2/ha)56.28±9.2755±8.153749.75±9.1929.146±5.5234.183±3.205
Basal area of trees with DBH>10  cm (m2/ha)55.8±8.87854.486±7.9848±9.7425.289±5.1832.44±3.18
Tree height (m)19.7±0.8318.9±1.3122.7±1.449.2±0.339.7±0.25
Stem density (/ha)6677146298401056
Stem density of trees with DBH>10  cm (/ha)488481456440601

OG: Old growth forests; VJR: virgin jungle reserve; LF: logged forest; EA: experimental area; OP: oil palm plantations.

Table 3

Aboveground forest parameters across the non-RF (NRF) zones.

(OG)NRF(VJR)NRF(LF)NRF(EA)NRF
Basal area of trees with DBH >10  cm (m2/ha)65.39±3.1NA32.13±13.4317.14±2.17
Stem density of trees with DBH >10  cm (/ha)820NA592417

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 Ref. 21. The AGB of the trees in both riparian and non-riparian zones was calculated using the biomass equation recommended by Ref. 22

(1)

AbovegroundBiomass(AGB)=0.0776×(ρ×DBH2×H)0.94,
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 Ref. 23

(2)

(AGB)Trunk=100×π×(r×z)2×h×ρ,
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:

(3)

ρ=(0.0076x+0.083)/100,
where x is the age of the oil palm plantation.

3.2.

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.

3.3.

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.

3.4.

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:

(4)

NDVI=(ρB4ρB3)(ρB4+ρB3)

(5)

SAVI=1.5(ρB4ρB3)(ρB4+ρB3+0.5),
where SAVI is a modification of NDVI and corrects for the influence of soil brightness when the vegetation cover is low.25

(6)

NDII5=(ρB4ρB5)(ρB4+ρB5),
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

3.5.

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 4

Texture variables derived from GLCM (from Ref. 19).

Texture variableFormulaDescription
MeanMEAN=∑ij=0N−1iPijMean of the probability values from the GLCM. It is directly related to the image spectral heterogeneity
VarianceVAR=∑ij=0N−1Pij(i−MEAN)2Measure of the global variation in the image. Large values denote high levels of spectral heterogeneity
CorrelationCOR=∑ij=0N−1Pij[ij−MEANij−MEAN)VAR]Measure of the linear dependency between neighbouring pixels
ContrastCONT=∑ij=0N−1Pij(i−j)2Quadratic measure of the local variation in the image. High values indicate large differences between neighbouring pixels
DissimilarityDISS=∑ij=0N−1Pij(i−j)Linear measure of the local variation in the image
HomogeneityHOM=∑ij=0N−1PPij1+(i−j)2Measure of the uniformity of tones in the image. A concentration of high values along the GLCM diagonal denotes to a high homogeneity
Angular second momentASM=∑ij=0N−1Pij2Measure of the order in the image. It is related to the energy required for arranging the elements in the system
EntropyENT=∑ij=0N−1PijInPijMeasure of the disorder in the image. It is inversely related to ASM

The texture variables are calculated on the basis of the red band.31 According to GLCM analysis carried out by Ref. 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.

4.

Results

4.1.

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.

Fig. 2

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

JARS_8_1_083695_f002.png

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.

4.2.

Distinguishing Between Different Land-Use Classes Using Reflectance and Vegetation Indices

4.2.1.

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.

Fig. 3

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.

JARS_8_1_083695_f003.png

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.

4.2.2.

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.

4.3.

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:

(7)

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

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).

Fig. 4

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.

JARS_8_1_083695_f004.png

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

(8)

Ba=2175.4(Band4)+2628.3

(9)

Ba=2776(Band3)+2817.4,
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.

4.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 5

Ability of texture-based linear regression models to explain variation in vegetation attributes.

Forest typeAGBBasal area
EAR2=0.86 (p<0.01)R2=0.92 (p<0.01)
LFR2=0.88 (p<0.01)R2=0.72 (p<0.01)
VJRR2=0.90 (p<0.01)R2=0.98 (p<0.01)
RFR2=0.91 (p<0.01)R2=0.84 (p<0.01)
OPNANA

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.

Fig. 5

(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.

JARS_8_1_083695_f005.png

5.

Discussion

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.

6.

Conclusions and Future Directions

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

References

1. 

F. Achardet al., “Determination of deforestation rates of the world’s humid tropical forests,” Science 297(5583), 999–1002 (2002), http://dx.doi.org/10.1126/science.1070656.SCIEAS0036-8075Google Scholar

2. 

M. C. Hansenet 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), http://dx.doi.org/10.1016/j.rse.2007.11.012.RSEEA70034-4257Google Scholar

3. 

G. P. Asneret al., “High-resolution forest carbon stocks and emissions in the Amazon,” Proc. Natl. Acad. Sci. U. S. A. 107(38), 16738–16742 (2010), http://dx.doi.org/10.1073/pnas.1004875107.PNASA60027-8424Google Scholar

4. 

B. A. Margonoet 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), http://dx.doi.org/10.1088/1748-9326/7/3/034010.ERLNAL1748-9326Google Scholar

5. 

H. TangkiN. A. Chappell, “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), http://dx.doi.org/10.1016/j.foreco.2008.07.018.FECMDW0378-1127Google Scholar

6. 

G. M. Foodyet al., “Mapping the biomass of Bornean tropical rain forest from remotely sensed data,” Global Ecol. Biogeogr. 10(4), 379–387 (2001), http://dx.doi.org/10.1046/j.1466-822X.2001.00248.x.GEBIFS1466-8238Google Scholar

7. 

G. M. FoodyD. S. BoydM. E. Cutler, “Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions,” Remote Sens. Environ. 85(4), 463–474 (2003), http://dx.doi.org/10.1016/S0034-4257(03)00039-7.RSEEA70034-4257Google Scholar

8. 

M. E. J. Cutleret 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), http://dx.doi.org/10.1016/j.isprsjprs.2012.03.011.IRSEE90924-2716Google Scholar

9. 

T. M. KuplichP. J. CurranP. M. Atkinson, “Relating SAR image texture to the biomass of regenerating tropical forests,” Int. J. Remote Sens. 26(21), 4829–4854 (2005), http://dx.doi.org/10.1080/01431160500239107.IJSEDK0143-1161Google Scholar

10. 

F. KayitakireC. HamelP. Defourny, “Retrieving forest structure variables based on image texture analysis and IKONOS-2 imagery,” Remote Sens. Environ. 102(3), 390–401 (2006), http://dx.doi.org/10.1016/j.rse.2006.02.022.RSEEA70034-4257Google Scholar

11. 

A. Wijayaet 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), http://dx.doi.org/10.1007/s11676-010-0001-7.1007-662XGoogle Scholar

12. 

M. A. Castillo-SantiagoM. RickerB. H. J. de Jong, “Estimation of tropical forest structure from SPOT-5 satellite images,” Int. J. Remote Sens. 31(10), 2767–2782 (2010), http://dx.doi.org/10.1080/01431160903095460.IJSEDK0143-1161Google Scholar

13. 

M. Azlizaet al., “Characterization of riparian plant community in lowland forest of peninsular Malaysia,” Int. J. Bot. 8(4), 181–191 (2012), http://dx.doi.org/10.3923/ijb.2012.181.191.IJBNCT1811-9700Google Scholar

14. 

“Stability of Altered Forest Ecosystems SAFE,” http://www.safeproject.net/ (2011).Google Scholar

15. 

A. C. LeesC. A. Peres, “Conservation value of remnant riparian forest corridors of varying quality for Amazonian birds and mammals,” Conserv. Biol. 22(2), 439–449 (2008), http://dx.doi.org/10.1111/j.1523-1739.2007.00870.x.CBIOEF0888-8892Google Scholar

16. 

M. A. MacDonald, “The role of corridors in biodiversity conservation in production forest landscapes: a literature review,” Tasforests 14, 41–52 (2003).Google Scholar

17. 

R. Ibbotson, “Personal email communication about Sabah’s logging history” (2013).Google Scholar

18. 

R. M. Ewerset al., “A large-scale forest fragmentation experiment: the Stability of Altered Forest Ecosystems Project,” Philos. Trans. R. Soc. B 366(1582), 3292–3302 (2011), http://dx.doi.org/10.1098/rstb.2011.0049.PTRBAE0962-8436Google Scholar

19. 

J. A. Gallardo-Cruzet al., “Predicting tropical dry forest successional attributes from space: is the key hidden in image texture?” PLoS One 7(2), e30506 (2012), http://dx.doi.org/10.1371/journal.pone.0030506.1932-6203Google Scholar

20. 

C. J. MarshR. M. Ewers, “A fractal-based sampling design for ecological surveys quantifying β-diversity,” Methods Ecol. Evol. 4(1), 63–72 (2013), http://dx.doi.org/10.1111/j.2041-210x.2012.00256.x.Google Scholar

21. 

A. C. Morelet 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), http://dx.doi.org/10.1016/j.foreco.2011.07.008.FECMDW0378-1127Google Scholar

22. 

J. Chaveet al., “Tree allometry and improved estimation of carbon stocks and balance in tropical forests,” Oecologia 145(1), 87–99 (2005), http://dx.doi.org/10.1007/s00442-005-0100-x.OECOBX0029-8549Google Scholar

23. 

A. C. MorelJ. B. FisherY. Malhi, “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), http://dx.doi.org/10.1080/01431161.2011.631949.IJSEDK0143-1161Google Scholar

24. 

USGS, “US Geological Survey,” http://srtm.usgs.gov/ (2011).Google Scholar

25. 

A. R. Huete, “A soil-adjusted vegetation index (SAVI),” Remote Sens. Environ. 25(3), 295–309 (1988), http://dx.doi.org/10.1016/0034-4257(88)90106-X.RSEEA70034-4257Google Scholar

26. 

C. M. Souza Jr.D. A. RobertsA. Monteiro, “Multitemporal analysis of degraded forests in the southern Brazilian Amazon,” Earth Interact. 9(19), 1–25 (2005), http://dx.doi.org/10.1175/EI132.1.EINTFP1087-3562Google Scholar

27. 

P. Ploton, “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).Google Scholar

28. 

H. G. JonesR. A. Vaughan, Remote Sensing of Vegetation: Principles, Techniques, and Applications, Oxford University Press, New York (2010).Google Scholar

29. 

V. St-Louiset al., “High-resolution image texture as a predictor of bird species richness,” Remote Sens. Environ. 105(4), 299–312 (2006), http://dx.doi.org/10.1016/j.rse.2006.07.003.RSEEA70034-4257Google Scholar

30. 

L. D. Esteset 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), http://dx.doi.org/10.1016/j.rse.2008.01.004.RSEEA70034-4257Google Scholar

31. 

S. Eckert, “Improved forest biomass and carbon estimations using texture measures from WorldView-2 satellite data,” Remote Sens. 4(4), 810–829 (2012), http://dx.doi.org/10.3390/rs4040810.17LUAF2072-4292Google Scholar

32. 

G. M. Foody, “Status of land cover classification accuracy assessment,” Remote Sens. Environ. 80(1), 185–201 (2002), http://dx.doi.org/10.1016/S0034-4257(01)00295-4.RSEEA70034-4257Google Scholar

33. 

T. M. Basukiet al., “The potential of spectral mixture analysis to improve the estimation accuracy of tropical forest biomass,” Geocarto Int. 27(4), 329–345 (2012), http://dx.doi.org/10.1080/10106049.2011.634928.1010-6049Google Scholar

34. 

S. A. Saderet 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), http://dx.doi.org/10.1016/0034-4257(89)90112-0.RSEEA70034-4257Google Scholar

35. 

B. Beguetet 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.Google Scholar

36. 

D. Luet al., “Application of spectral mixture analysis to Amazonian land-use and land-cover classification,” Int. J. Remote Sens. 25(23), 5345–5358 (2004), http://dx.doi.org/10.1080/01431160412331269733.IJSEDK0143-1161Google Scholar

37. 

D. Luet al., “Aboveground forest biomass estimation with Landsat and lidar data and uncertainty analysis of the estimates,” Int. J. For. Res. 2012, 436537 (2012), http://dx.doi.org/10.1155/2012/436537.IJVNAP0300-9831Google Scholar

Biography

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.

Minerva Singh, Yadvinder Malhi, Shonil Bhagwat, "Evaluating land use and aboveground biomass dynamics in an oil palm–dominated landscape in Borneo using optical remote sensing," Journal of Applied Remote Sensing 8(1), 083695 (10 January 2014). http://dx.doi.org/10.1117/1.JRS.8.083695
JOURNAL ARTICLE
14 PAGES


SHARE
KEYWORDS
Vegetation

Remote sensing

Biological research

Earth observing sensors

Landsat

Data modeling

Reflectivity

Back to Top