Next Article in Journal
The Urban Environment Can Modify Drought Stress of Small-Leaved Lime (Tilia cordata Mill.) and Black Locust (Robinia pseudoacacia L.)
Previous Article in Journal
Effect of Organic Layer Thickness on Black Spruce Aging Mistakes in Canadian Boreal Forests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geospatial Estimation of above Ground Forest Biomass in the Sierra Madre Occidental in the State of Durango, Mexico

by
Pablito M. López-Serrano
1,
Carlos A. López Sánchez
2,*,†,
Raúl Solís-Moreno
3,† and
José J. Corral-Rivas
2,†
1
Doctorado Institucional en Ciencias Agropecuarias y Forestales, Universidad Juárez del Estado de Durango, Boulevard Guadiana 501, Fraccionamiento Ciudad Universitaria, 34120 Durango, Mexico
2
Instituto de Silvicultura e Industria de la Madera, Universidad Juárez del Estado de Durango, Boulevard Guadiana 501, Fraccionamiento Ciudad Universitaria, 34120 Durango, Mexico
3
Facultad de Ciencias Forestales, Universidad Juárez del Estado de Durango, Río Papaloapan, Valle del Sur, 34120 Durango, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Forests 2016, 7(3), 70; https://doi.org/10.3390/f7030070
Submission received: 19 November 2015 / Revised: 8 March 2016 / Accepted: 10 March 2016 / Published: 15 March 2016

Abstract

:
Combined use of new geospatial techniques and non-parametric multivariate statistical methods enables monitoring and quantification of the biomass of large areas of forest ecosystems with acceptable reliability. The main objective of the present study was to estimate the aboveground forest biomass (AGB) in the Sierra Madre Occidental (SMO) in the state of Durango, Mexico, using the M5 model tree (M5P) technique and the analysis of medium-resolution satellite-based multi-spectral data, and field data collected from a network of 201 permanent forest growth and soil research sites (SPIFyS). Research plots were installed by systematic sampling throughout the study area in 2011. The digital levels of the images were converted to apparent reflectance (ToA) and surface reflectance (SR). The M5P technique that constructs tree-based piecewise linear models was used. The fitted model with SR and tree abundance by species group as predictive variables (ASG) explained 73% of the observed AGB variance (the root mean squared error (RMSE) = 39.40 Mg·ha−1). The variables that best discriminated the AGB, in order of decreasing importance, were the normalized difference vegetation index (NDVI), tree abundance of other broadleaves species (OB), Band 4 of Landsat 5 TM (Thematic Mapper) satellite and tree abundance of pines (Pinus). The results demonstrate the potential usefulness of the M5P method for estimating AGB based in the surface reflectance values (SR).

1. Introduction

The Sierra Madre Occidental (SMO) mountain range is of great ecological interest because of its environmental heterogeneity, which is attributed to the broad physiographical and climatic diversity in the area [1]. Moreover, the SMO is home to pine and oak species that are economically important in ecosystems in Mexico and other parts of the world [2]. The SMO crosses several states in western Mexico, including the state of Durango (the SMO occupies 71.5% of the surface area of the state). The state of Durango generates between 25% and 30% of the national timber production, producing a total of 1.5 million·m3 of roundwood per year, and boasts forest reserves that are important sources of environmental services [3]. Studies that attempt to estimate forest biomass in this type of ecosystem are expensive due to its large coverage and difficult access for direct estimation of biomass. Thus, the emergence of geospatial techniques is becoming increasingly relevant for estimating and monitoring forest biomass in short periods of time because of its low cost and acceptable accuracy [4,5,6,7]. Because of the macrospatial scale and high heterogeneity of these ecosystems, the quantitative data obtained often do not comply with the underlying assumptions of simple statistical analysis (homogeneity and normality of distribution), so other techniques such as logistic regression and non-parametric classification methods are often applied [8,9,10]. The M5 model tree (M5P) technique is a reconstruction of M5 algorithm for inducing trees of regression models [11]. M5P is used for numeric prediction and at each leaf it stores a linear regression model that predicts the class value of instances that reach the leaf. To determine which attribute is the best to split the portion of the training data that reaches a particular node, the splitting criterion is used. The standard deviation of the training class is treated as a measure of the error at that node and each attribute at that node is tested by calculating the expected reduction in error. The attribute that is chosen for splitting maximizes the expected error reduction at that node. The main objective of the present study was to estimate the aboveground forest biomass (AGB) in the SMO in the state of Durango, Mexico, using the M5P technique and the analysis of medium-resolution satellite-based multi-spectral data, and field data collected from a network of 201 permanent forest growth and soil research sites (SPIFyS).

2. Material and Methods

2.1. Study Area

The study area is a mountainous zone in the state of Durango (Mexico) that forms part of the Sierra Madre Occidental (Figure 1).
The area occupied by the state of Durango represents 6.3% of the land in Mexico. The total area covered by the state is 12.3 million ha, of which 9.1 million ha (74.35% of the land in the state) is forestland managed by 11 Regional Forest Management Units (UMAFORES). A large part of the forestland (4.9 million ha) is occupied by temperate forest and is subjected to precipitation levels of between 800 and 1200 mm per year, with frost occurring in winter as a result of the combination of low temperatures and humid winds from the Pacific Ocean; a smaller area of the land (0.5 million), affected by warmer climate, is occupied by forest classified as rainforest [12]. The mean elevation in this zone is 2650 m above sea level. These forests have rich biodiversity and include at least 27 coniferous tree species (of which 20 are Pinus species) and 43 species of Quercus; the predominant forest stands comprise pines and oaks, often mixed with Arbutus and Juniperus, among other tree species [13]. These unique forests are irregular and have been subject to selective harvesting for almost a century to provide a mix of services to local communities. This irregularity refers to the spatial arrangement of trees (vertical and horizontal irregularity) and the variation in the age structure of trees and stands. This structure is the result of the management history, which has depended on land ownership, as well as the economic and social changes that have taken place in the state, and also natural conditions [14].

2.2. Field Data

The dasometric data were obtained from 201 permanent forest growth and soil monitoring plots (SPIFyS) in the SMO in the state of Durango. The plots were installed during the winter of 2011 using the protocol developed by [15]. The data of these permanent sample plots are used to monitor the growth and yield of Durango’s forests. The plots cover the main forest types and the current diameter distributions of commercial forests in Durango. The plots are 50 × 50 m in size (distance was corrected by the slope) and are distributed by systematic sampling (with some exceptions), with a variable grid ranging from 3 to 5 kilometers, depending on the size of the “Ejidos”. Ejidos are communal groups that live in rural areas and whose lands are managed with some level of governmental control. The sampling plots are intended to be re-measured at five-year intervals. Among other variables, tag number, species code, breast height diameter (measured in cm at 1.3 m above ground level), total tree height (m), height to the live crown (m), azimuth (°) and radius (m) from the center of the plot of all trees equal or larger than 7.5 centimeters (cm) in diameter were recorded. The database used here includes measurement data from 31,979 trees.
The aboveground biomass in each of the SPIFyS plots was estimated using specific allometric equations developed by [16] for the same study area. Depending on the species, the goodness of fit statistics ranged between 0.82 and 0.97 of the coefficient of determination (R2) and the root mean square error (RMSE) between 22.68 and 133.68 kg.
The main descriptive statistics for the total aboveground biomass per hectare in the study sites are summarized in Table 1.

2.3. Tree Abundance by Species Group (ASG)

The tree abundance (number of trees per group of species per plot) was estimated for posterior analysis in this study. A total of seventy-two different tree species were grouped in four groups of species as they present similar growth patterns: (P) Pinus species (16); (OC) other conifers species (12); (Q) oaks species (26); and (OB) other broadleaves species (18).

2.4. Source of Spectral Data

The data used for the study were obtained from six Landsat 5 TM (Thematic Mapper) satellite images captured between March and May 2011 and covering all of the SMO within the state of Durango (path/row: 30/44, 31/42, 31/43, 31/44, 32/42 and 32/43) [17]. This satellite platform, of medium spatial resolution, operates in seven bands of the electromagnetic spectrum: blue (bandwidth 0.45–0.52 µm), green (bandwidth 0.52–0.60 µm), red (0.63–0.69 µm), near infrared (0.78–0.89 µm), mid infrared (1.55–1.75 µm) and far infrared (2.08–2.35 µm). These bandwidths correspond, respectively, to bands 1, 2, 3, 4, 5 and 7 of the Landsat TM5 satellite [18]. Band 6, designed for the thermal mapping and soil moisture, was not considered because of its lower (120 m) spatial resolution.
The satellite images were digitally pre-processed by radiometric correction techniques, according to the procedures suggested by [19,20]. The images are produced by USGS with a rectification using a cubic convolution geometric correction for discrete data (level L1T), with a root mean square error (RMSE) of less than 1 pixel, thus making them suitable for digital image processing [21]. The digital levels (DLs) were converted to radiance values to generate images that were calibrated with the minimal radiance (Lmin) and maximal radiance (Lmax) values for each band of the sensor [22]. The radiance was subsequently converted to apparent reflectance (Top of Atmosphere (ToA)) with the aim of converting the original values of each image into standard physical variables that are comparable over time for the same sensor [19]. This process was carried out with IDRISI® Selva software [23] and the ATMOS algorithm, which fits the radiometric effect on considering the solar elevation, yielding an image with reflectance values (0–1).
Furthermore, the same images were downloaded from the National Landsat Archive Processing System (NLAPS), corresponding to the product Landsat 4–5 Thematic Mapper level 1 of reflectance on surfaces (SR), radiometrically and atmospherically corrected, and processed through the Standard Landsat Product Generation System (LPGS) using the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) algorithm [17].
Bands 1, 2, 3, 4, 5 and 7 (Band 1 to Band 7) of Landsat TM5 were used; Band 6 was not used, because of its thermal characteristics [17]. The Normalized Difference Vegetation Index (NDVI) was also calculated, with the aim of compensating the factors that influence the images in relation to biomass estimation, such as the illumination conditions, the slope and the orientation of the surface. The use of NDVI calculated for ToA and SR, as a predictor variable to model AGB has been successfully reported in previous studies [24,25,26].
N D V I = ( N I R R ) ( N I R + R )
where NIR is the spectral band in the near infrared region (Band 4) and R the band in the red region (Band 3).

2.5. Integration of Data Files

Once the images were obtained by the previously mentioned processes (ToA and SR spectral bands), a mosaic was constructed with six of the scenes covering the SMO. Posterior geolocation of the SPIFyS in the mosaic enabled extraction of the information at the pixel level by bilinear interpolation. ArcGIS 10® software [27] was used for this extraction. ToA and SR values were integrated in a database together with the extracted total aboveground biomass (Mg·ha−1) as inputs for the model.

2.6. Fitted Model

We used a machine learning technique to estimate the AGB at stand level. M5P technique combines a conventional decision tree with the possibility of linear regression functions at the nodes. First, a decision-tree induction algorithm is used to build a tree, but instead of maximizing the information gain at each inner node, a splitting criterion is used that minimizes the intra-subset variation in the class values down each branch. The splitting procedure in M5P stops if the class values of all instances that reach a node vary very slightly, or only a few instances remain. Second, the tree is pruned back from each leaf. When pruning, an inner node is turned into a leaf with a regression plane. Third, to avoid sharp discontinuities between the subtrees, a smoothing procedure is applied that combines the leaf model prediction with each node along the path back to the root, smoothing it at each of these nodes by combining it with the value predicted by the linear model for that node. Techniques devised by [11] for their classification and regression trees system are adapted in order to deal with enumerated attributes and missing values. All enumerated attributes are turned into binary variables so that all splits in M5P are binary. As to missing values, M5P uses a technique called “surrogate splitting” that finds another attribute to split on in place of the original one and uses it instead [11,28,29].
In this study, in a first stage, the six spectral bands of the Landsat 5 TM sensor (1, 2, 3, 4, 5 and 7) and the NDVI were analyzed with the algorithms ToA and SR to estimate AGB. In a second stage, its spectral variables (SR) were evaluated with variables that incorporate aspects of forest structure (ASG). All analyses were performed with M5P technique implemented into the WEKA open source software [30].
To compare the performance of the models, the coefficient of determination (R2), the root mean squared error (RMSE) and the root relative squared error (RRSE) were used as goodness-of-fit criteria for evaluating model performance and were expressed as follows:
R 2 = 1 i = 1 n ( y i y ^ i ) 2 i = 1 n ( y i y ¯ i ) 2
R M S E = i = 1 n ( y i y ^ i ) 2 n p
R R S E = i = 1 n ( y i y ^ i ) 2 i = 1 n ( y ^ i y ¯ i ) 2
where, yi, y ^ i and y ¯ i are the observed, estimated and mean values of AGB, respectively; n is the total number of observations used to fit the model; and p is the number of model parameters.
The selected model was applied for mapping AGB in the SMO area using ArcGIS 10® software [27].

3. Results

The decision tree generated by the M5P technique for ToA, SR and SR with ASG variables were implemented in WEKA software, using the pixel level values extracted from the images of the 201 SPIFyS plots is shown in Figure 2.
In accordance with the hierarchical structure of the decision trees, the following variables that best discriminated or predicted the AGB, in order of decreasing importance were, for ToA: Band 7, Band 3, Band 1, NDVI and Band 5; for SR: NDVI, Band 1 and Band 7; and for SR with ASG: NDVI, OB, Band 4 and tree abundance of pines. Categorization of the trees continued following the path determined by the responses to the questions at the internal nodes, until reaching a terminal node, where the predetermined label will be that assigned to the classification pattern—in this case, the pixel values for AGB estimation. The Table 2 show the goodness-of-fit statistics derived from the M5P technique with ToA values explained 54% (R2) of the observed variability in the AGB of the 201 research plots, with a RMSE of 50.47 Mg·ha−1. SR model explained 69% (RMSE = 42.17 Mg·ha−1), and when including the ASG variables the explanation of the variance increases to 73% (RMSE = 39.40 Mg·ha−1).
Graphical analysis of the residual values and the observed values plotted against the predicted values of AGB did not reveal any important problems in relation to heterogeneity of the variance or lack of normal distribution of the residuals, with the exception of a slight trend of underestimation for high AGB (Figure 3).
The spatial distribution of the estimated AGB (Mg·ha−1) in the SMO area obtained by the application of the classification rules included in the regression tree model (M5P) for SR variables is shown in Figure 4. The lighter color pixels represent the lowest amounts of AGB, below 75 Mg·ha−1, whereas the dark green pixels represent the largest amounts of AGB, which consistently correspond to the most dense areas of temperate forest. Calculated mean amount of AGB for the study area was around 106 Mg·ha−1.
The total AGB content estimations from the MP5 technique for SR variables for the analyzed forest management units are shown in Table 3. The highest mean value of AGB was observed in the UMAFOR 1006 (Municipally of San Dimas) with 148.98 Mg·ha−1 and a total estimation of 64,033,008.59 Mg. This zone encompasses the largest area of forestland and therefore the largest amount of AGB. On the other hand, the lowest amount of AGB was observed in the UMAFOR 1001 with a mean estimate of 78.66 Mg·ha−1, making it the forest region with the lowest density out of the eleven forest management units considered in this study.

4. Discussion

The results of the present study demonstrate that the data acquired by a medium spatial resolution (Landsat) sensor are potentially useful for estimating AGB in structurally complex forests, such as those in the SMO in the state of Durango (Mexico), with satisfactory results and low cost. The deterministic predictors were the bands belonging to the blue, green and most red, near and mid infrared spectral regions. This finding was similar to that reported by [31], who demonstrated that reflectance in the red and near infrared regions yielded good predictions for AGB estimation in forest zones of the Yellowstone National Park, USA. In the present study, the model tended to underestimate AGB values above approximately 250 Mg·ha−1. This might possibly be due to the saturation of NDVI, which is the most influential variable in predicting biomass for high values. In this sense, several studies have similarly found that the NDVI loses its sensitivity to dense vegetation because of the saturation in red and near infrared wavelength in measuring and monitoring plant growth, vegetation cover and biomass production from satellite data [32,33,34]. Models fitted ToA, SR, and SR with ASG, respectively, showed an increasing capacity to overcome these NDVI saturation problems.
Furthermore, [35] concluded that the vegetation indices or individual bands, which include one or more bands in the infrared spectrum, provide satisfactory descriptions of zones occupied by conifer or broadleaf species. Moreover, in the case of SR, which was better than ToA, the contribution of Band 1 in two terminal nodes of the M5P model is associated with the structural variability of the canopy [36]. Günlü et al. [37] found that the reflectance from Landsat TM satellite Band 1 was the best predictor of AGB (R2 = 0.465, RMSE = 91,836 t·ha−1), given the structural conditions of the canopy and understory, as the reflectance from this band increased as the AGB increased. In the present study, the result of the M5P analysis with SR spectral bands (R2 = 0.69, RMSE = 42.17 Mg·ha−1) was higher to that reported by Houghton et al. [38], who analyzed data from the MODIS sensor (resolution, 500 m) and forest inventory data using the non-parametric Random Forest (R2 = 0.61) method to map forest biomass in Russia.
In a recent study, Tian et al. [39] used the non-parametric k-nearest neighbours (k-NN) technique to produce an optimized model (R2 = 0.59, RMSE = 24.92 ton·ha−1) from Landsat-TM images of a sample of 133 plots, with topographic correction based on sun-canopy-sensor (SCS + C). Likewise, Hall, [30] used Landsat 5 (TM) images rectified by SCS+C radiometric correction and compared the performance of the k-NN method and support vector machine (SVM) method for estimating AGB. They found that k-NN performed better (R2 = 0.54; RMSE = 26.62 ton·ha−1) than SVM (R2 = 0.51; RMSE = 27.45 ton·ha−1).
In general, most previous studies report significant relationships between AGB and the reflectance values yielded by each sensor. The reliability was within the range reported in diverse research studies that estimate AGB from medium resolution spectral data, from Landsat and SPOT, which often yield R2 values between 0.50 and 0.70 with absolute errors of the estimates of between 30 and 60 Mg·ha−1 [5,40,41,42,43,44,45,46,47]. The present study also shows that incorporation of spectral data and tree abundance estimated by species group in mixed and uneven-aged forests (SR with ASG), such as the SMO, can increase the level of estimation of the AGB (R2 = 0.73; RMSE = 39.40 Mg·ha−1). In this sense, previous studies have reported significant variations in forest biomass estimation between different ecological zones, tree species, ages, density and management types [48,49,50].
In other studies [51,52,53], several authors have concluded that the spectral data derived after atmospheric and topographic correction may improve the accuracy of the biomass estimation, irrespective of the statistical method used. As the areas being monitored are mountainous zones, the quality of the data is negatively affected by the reflectance between sunny and shaded slopes. Interactive parameter fitting in the topographical correction methods may improve the quality of the spectral data and of the AGB estimates [52,53].

5. Conclusions

In the present study, we estimated the AGB in the SMO in the state of Durango, Mexico, using the M5P technique and the analysis of medium-resolution satellite-based multi-spectral data, and field data collected from a network of 201 SPIFyS.
The findings show that the M5P method is potentially useful for estimating forest biomass. Data from the infrared channel of the Landsat TM5 sensor proved best for discriminating or predicting AGB.
The surface reflectance values (SR) in comparison with atmospheric correction from the sensor (ToA), was best for the estimation of AGB.
The results of this study indicate that performing atmospheric corrections and considering variables related to forest structure (SR with ASG variables) can help to solve problems of saturation of NDVI for high values of biomass.

Acknowledgments

This research was supported by SEP-PROMEP (project: Seguimiento y Evaluación de Sitios Permanentes de Investigación Forestal y el Impacto Socio-económico del Manejo Forestal en el Norte de México).

Author Contributions

Pablito M. López-Serrano and Carlos A. López-Sánchez conceived, designed and performed the experiments and wrote the manuscript. Raúl Solís-Moreno and José J. Corral-Rivas analyzed the data and revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. González Elizondo, M.S.; González Elizondo, M.; Tena Flores, J.A.; Ruacho González, L.; López-Enríquez, I.L. Vegetación de la Sierra Madre Occidental, México: una síntesis. Acta Botánica Mex. 2012, 100, 351–403. [Google Scholar]
  2. Sánchez, O.; Vega, E.; Peters, E.Y.; Monrroy, V.O. Conservación de Ecosistemas Templados de Montaña en México; Instituto Nacional de Ecología (INE–SEMARNAT): México D.F., México, 2003. [Google Scholar]
  3. Secretaría de Medio Ambiente y Recursos Naturales. Anuario Estadístico de la Producción Forestal 2011; México, 2011; pp. 11–24. [Google Scholar]
  4. Foody, G.M.; Boyd, D.S.; Cutler, M.E.J. Predictive relations of tropical forest biomass from Landsat TM data and their transferability between regions. Remote Sens. Environ. 2003, 85, 463–474. [Google Scholar] [CrossRef]
  5. Hall, R.J.; Skakun, R.S.; Arsenault, E.J.; Case, B.S. Modeling forest stand structure attributes using Landsat ETM+ data: Application to mapping of aboveground biomass and stand volume. For. Ecol. Manag. 2006, 225, 378–390. [Google Scholar] [CrossRef]
  6. Fuchs, H.; Magdon, P.; Kleinn, K.; Flessa, H. Estimating aboveground carbon in a catchment of the Siberian forest tundra: Combining satellite imagery and field inventory. Remote Sens. Environ. 2009, 113, 518–531. [Google Scholar] [CrossRef]
  7. Verbesselt, J.; Hyndman, R.; Newnham, G.; Culvenor, D. Detecting trend and seasonal changes in satellite image time series. Remote Sens. Environ. 2010, 114, 106–115. [Google Scholar] [CrossRef]
  8. Barajas, F.H. Comparación Entre Análisis Discriminante no-Métrico y Regresión Logística Multinomial. Tesis de Maestría, Facultad de Ciencias, Universidad Nacional de Colombia, Medellín, Colombia, 2007. [Google Scholar]
  9. Gibbons, J.D.; Chakraborti, S. Nonparametric Statistical Interference; Marcel Denker, Inc.: New York, NY, USA, 2003; p. 645. [Google Scholar]
  10. Karjalainen, M.; Kankare, V.; Vastaranta, M.; Holopainen, M.; Hyyppa, J. Prediction of plot-level forest variables using TerraSAR-X stereo SAR data. Remote Sens. Environ. 2012, 117, 338–347. [Google Scholar] [CrossRef]
  11. Quinlan, J.R. Learning with Continuous Classes. In 5th Australian Joint Conference on Artificial Intelligence; Word Scientific: Singapore, 1992; pp. 343–348. [Google Scholar]
  12. SRNyMA-CONAFOR. Plan Estrátegico Forestal 2030; Gobierno del estado de Durango: Durango, México, 2007. [Google Scholar]
  13. Zhao, X.; Corral-Rivas, J.J.; Zhang, C.; Temesgen, H.; Gadow, K. Forest observational studies-an essential infrastructure for sustainable use of natural resources. Forest Ecosyst. 2014. [Google Scholar] [CrossRef]
  14. Wehenkel, C.; Corral-Rivas, J.J.; Hernández-Díaz, J.C.; Gadow, K. Estimating Balanced Structure Areas in multi-species forests on the Sierra Madre Occidental, Mexico. Ann. For. Sci. 2011, 68, 385–394. [Google Scholar] [CrossRef]
  15. Corral Rivas, J.; Vargas, B.; Wehenkel, C.; Aguirre, O.; Álvarez, J.; Rojo, A. Guía para el Establecimiento de Sitios de Inventario Periódico Forestal y de Suelos del Estado de Durango; Facultad de Ciencias Forestales. Universidad Juárez del Estado de Durango: Durango, México, 2009. [Google Scholar]
  16. Vargas-Larreta, B.; González-Herrera, L.; López-Sánchez, C.A.; Corral-Rivas, J.J.; López-Martínez, J.O.; Aguirre-Calderón, C.G.; Álvarez-González, J.J. Biomass equations and Carbon content of the temperate forests of northwestern México. Biomass Bioenerg. 2015, in press. [Google Scholar]
  17. United States Geological Survey. Available online: http://glovis.usgs.gov (accessed on 11 January 2011).
  18. NASA. Landsat 7 Science Data Users Handbook, 2011. Available online: http://landsathandbook.gsfc.nasa.gov/pdfs/Landsat7_Handbook.pdf (accessed on 11 January 2011).
  19. Greenle, P.T. Remote Sensing science applications in arid environments. Remote Sens. Environ. 1993, 23, 143–154. [Google Scholar]
  20. Chuvieco, E. Teledetección Ambiental: La observación de la tierra desde el espacio, 3rd ed.; Ariel.: Barcelona, España, 2010; pp. 262–306. [Google Scholar]
  21. Keys, R.G. Cubic convolution interpolation for digital image processing. IEEE Trans. Acoust. Speech Signal Process. 1981, 29, 1153–1160. [Google Scholar] [CrossRef]
  22. Eastman, J.R. IDRISI, Selva., Guía para SIG y Procesamiento de Imágenes; Clark Labs Clark University: Worcester, MA, USA, 2012. [Google Scholar]
  23. Eastman, J.R. IDRISI version Selva, software for GIS and image processing; Clark University: Worcester, MA, USA, 2012.
  24. Gasparri, N.I.; Parmuchi, M.G.; Bono, J.; Karszenbaum, H.; Montenegro, C.L. Assessing multi-temporal Landsat 7 ETM+ images for estimating above-ground biomass in subtropical dry forests of Argentina. J. Arid Environ. 2010, 74, 1262–1270. [Google Scholar] [CrossRef]
  25. Liang, S.; Li, X.; Wang, J. Advanced Remote Sensing: Terrestrial Information Extraction and Applications; Academic Press: Oxford, UK, 2012. [Google Scholar]
  26. Zhu, X.; Liu, D. Improving forest aboveground biomass estimation using seasonal Landsat NDVI time-series. ISPRS J. Photogramm. Remote Sens. 2015, 102, 222–231. [Google Scholar] [CrossRef]
  27. ArcGIS for Desktop, version 10; software for analyzing geographic information systems; ESRI: Redlands, CA, USA, 1999–2012.
  28. Breiman, L.; Friedman, J.; Olshen, R.; Sone, C. Classification and Regression Trees; Wadsworth International Group: Belmont, CA, USA, 1984. [Google Scholar]
  29. Wang, Y.; Witten, I.H. Induction of model trees for predicting continuous classes. In Proceedings of the Poster Papers of the 9th European Conference on Machine Learning, Hamilton, New Zealand, 23 October 1997.
  30. Hall, M. Correlation-based Feature Selection for Machine Learning. Ph.D. Thesis, University of Waikato, Hamilton, New Zealand, 1999. [Google Scholar]
  31. Jakubauskas, M.E. Thematic Mapper characterization of lodge pole pine serals in Yellowstone National Park, USA. Remote Sens. Environ. 1996, 56, 118–132. [Google Scholar] [CrossRef]
  32. Huete, A.R.; Liu, H.; van Leeuwen, W.J. The Use of Vegetation Indices in Forested Regions: Issues of Linearity and Saturation. In Proceedings of 1997 IEEE International Geoscience and Remote Sensing Symposium, Singapore, 3–8 August 1997; European Space Agency Publications Division: Noordwijk, The Netherlands, 1997; pp. 1966–1968. [Google Scholar]
  33. Lu, D.; Mausel, P.; Brondızio, E.; Moran, E. Relationships between forest stand parameters and Landsat TM spectral responses in the Brazilian Amazon Basin. For. Ecol. Manag. 2004, 198, 149–167. [Google Scholar] [CrossRef]
  34. Mutanga, O.; Skidmore, A.K. Narrow band vegetation indices overcome the saturation problem in biomass estimation. Int. J. Remote Sens. 2004, 25, 3999–4014. [Google Scholar] [CrossRef]
  35. Fassnacht, K.S.; Gower, S.T.; MacKenzie, M.D.; Nordheim, E.V.; Lillesand, T.M. Estimating the leaf area index of North Central Wisconsin forests using the Landsat Thematic Mapper. Remote Sens. Environ. 1997, 61, 229–245. [Google Scholar] [CrossRef]
  36. Palestina, R.A.; Equihua, M.; Pérez-Maqueo, O.M. Influencia de la complejidad estructural del dosel en la reflectancia de datos Landsat TM. Madera y Bosques 2015, 21, 63–75. [Google Scholar]
  37. Günlü, A.; Ercanli, İ.; Başkent, E.Z.; Çakır, G. Estimating aboveground biomass using Landsat TM imagery: A case study of Anatolian Crimean pine forests in Turkey. Ann. For. Res. 2014, 57, 289–298. [Google Scholar]
  38. Houghton, R.A.; Butman, D.; Bunn, A.G.; Krankina, O.N.; Schlesinger, P.; Stone, T.A. Mapping Russian forest biomass with data from satellites and forest inventories. Environ. Res. Lett. 2007, 4, 1–7. [Google Scholar] [CrossRef]
  39. Tian, X.; Zengyuan, L.; Zhongbo, S.; Erxue, C.; Christiaan, T.; Xin, L.; Guo, Y.; Longhui, L.; Feilong, L. Estimating montane forest above-ground biomass in the upper reaches of the Heihe River Basin using Landsat-TM data. Int. J. Remote Sens. 2014, 35, 7339–7362. [Google Scholar] [CrossRef]
  40. Guoa, Y.; Tianb, X.; Lib, Z.; Linga, F.; Chenb, E.; Yanb, M.; Lib, C. Comparison of estimating forest above-ground biomass over montane area by two non-parametric methods. ISPRS J. Photogramm. Remote Sens. 2014, 92, 137–146. [Google Scholar]
  41. Foody, G.M.; Cutler, M.E.; Mcmorrow, J.; Pelz, D.; Tangki, H.; Boyd, D.S.; Douglas, I. Mapping the biomass of Bornean tropical rain forest from remotely sensed data Global. Ecol. Biogeogr. 2001, 10, 379–387. [Google Scholar] [CrossRef]
  42. Tomppo, E.; Nilsson, M.; Rosengren, M.; Aalto, P.; Kennedy, P. Simultaneous use of Landsat-TM and IRS-1C WiFS data in estimating large area tree stem volume and aboveground biomass. Remote Sens. Environ. 2002, 82, 156–171. [Google Scholar] [CrossRef]
  43. Zheng, D.; Rademacher, J.; Chen, J.; Crow, T.; Bresee, M.; le Moine, J.; Ryu, S.R. Estimating aboveground biomass using Landsat 7 ETM+ data across a managed landscape in northern Wisconsin, USA. Remote Sens. Environ. 2004, 93, 402–411. [Google Scholar] [CrossRef]
  44. Muukkonen, P.; Heiskanen, J. Estimating biomass for boreal forests using ASTER satellite data Combined with standwise forest inventory data. Remote Sens. Environ. 2005, 99, 434–447. [Google Scholar] [CrossRef]
  45. Chen, W.; Blain, D.; Li, J.; Keohler, K.; Fraser, R.; Zhang, K.; Leblanc, S.; Olthof, I.; Wang, J.; McGovern, M. Biomass measurements and relationships with Landsat-7/ETM+ and JERS-1/SAR data over Canada's western sub-arctic and low arctic. Int. J. Remote Sens. 2009, 30, 2355–2376. [Google Scholar] [CrossRef]
  46. Castillo Santiago, M.A.; Ricker, M.; Jong, B.H.J. Estimation of tropical forest structure from SPOT-5 satellite images. Int. J. Remote Sens. 2010, 31, 2767–2782. [Google Scholar] [CrossRef]
  47. Tian, X.; Zhongbo, S.; Erxue, C.; Zengyuan, L.; Christiaan, T.; Jianping, G.; Qisheng, H. Estimation of forest above-ground biomass using multi-parameter remote sensing data over a cold and arid area. Int. J. Appl. Earth Obs. Geoinform. 2012, 14, 160–168. [Google Scholar] [CrossRef]
  48. Henry, M.; Picard, N.; Trotta, C.; Manlay, R.J.; Valentini, R.; Bernoux, M.; Saint-André, L. Estimating tree biomass of sub-Saharan African forests: A review of available allometric equations. Silva. Fenn. 2011, 45, 477–569. [Google Scholar] [CrossRef]
  49. De-Miguel, S.; Pukkala, T.; Assaf, N.; Shater, Z. Intra-specific differences in allometric equations for aboveground biomass of eastern Mediterranean Pinus. Brutia. Ann. For. Sci. 2014, 71, 101–112. [Google Scholar] [CrossRef]
  50. Zou, W.-T.; Zeng, W.-S.; Zhang, L.-J.; Zeng, M. Modeling Crown Biomass for Four Pine Species in China. Forests 2015, 6, 433–449. [Google Scholar] [CrossRef]
  51. Richter, R.; Kellenberger, T.; Kaufmann, H. Comparison of topographic correction methods. Remote. Sens. 2009, 3, 184–196. [Google Scholar] [CrossRef] [Green Version]
  52. Hantson, S.; Chuvieco, E. Evaluation of different topographic correction methods for Landsat imagery. Int. J. Appl. Earth Obs. Geoinform. 2011, 5, 691–700. [Google Scholar] [CrossRef]
  53. Balthazar, V.; Vanacker, V.; Lambin, E.F. Evaluation and parameterization of ATCOR3 topographic correction method for forest cover mapping in mountain areas. Int. J. Appl. Earth Obs. Geoinform. 2012, 18, 436–450. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Forests 07 00070 g001
Figure 2. Decision tree obtained using the M5P technique with ToA (upper), SR (middle) and SR with ASG variables (bottom). (Band 1 to 7): of Landsat 5 TM (Thematic Mapper) satellite, (NDVI): normalized difference vegetation index, (OB) tree abundance of other broadleaves species, and (Pinus) tree abundance of pines.
Figure 2. Decision tree obtained using the M5P technique with ToA (upper), SR (middle) and SR with ASG variables (bottom). (Band 1 to 7): of Landsat 5 TM (Thematic Mapper) satellite, (NDVI): normalized difference vegetation index, (OB) tree abundance of other broadleaves species, and (Pinus) tree abundance of pines.
Forests 07 00070 g002
Figure 3. Graphs showing the distribution of the residuals and of the observed AGB values with ToA (upper), SR (middle) and SR with ASG variables (bottom).
Figure 3. Graphs showing the distribution of the residuals and of the observed AGB values with ToA (upper), SR (middle) and SR with ASG variables (bottom).
Forests 07 00070 g003
Figure 4. Spatial distribution of the total AGB in the SMO, state of Durango, Mexico.
Figure 4. Spatial distribution of the total AGB in the SMO, state of Durango, Mexico.
Forests 07 00070 g004
Table 1. Descriptive statistics of the total aboveground biomass per hectare in the 201 permanent forest growth and soil monitoring plots (SPIFyS).
Table 1. Descriptive statistics of the total aboveground biomass per hectare in the 201 permanent forest growth and soil monitoring plots (SPIFyS).
VariableMeanStandard DeviationMinimum ValueMaximum Value
Number of stems per ha645271.842242264
Stand basal area (m2·ha−1)23.448.068.2154.83
Dominant height (m)17.475.086.8630.60
Stand biomass (Mg·ha−1)141.6475.0127.73469.42
Table 2. Summary of the goodness-of-fit statistics for estimation of the AGB.
Table 2. Summary of the goodness-of-fit statistics for estimation of the AGB.
StatisticsToASRSR with ASG
R20.540.690.73
RMSE50.4742.1739.40
RRSE67.4556.3652.66
Table 3. Estimation of AGB for the regional forest management units in the SMO, state of Durango, Mexico.
Table 3. Estimation of AGB for the regional forest management units in the SMO, state of Durango, Mexico.
UMAFORMean AGBSurface AreaTotal AGB
(Mg·ha−1)(ha)(Mg)
100178.66423,990.0033,350,319.02
100285.58351,498.0030,079,977.18
1003107.10126,054.0013,500,870.12
100499.54318,104.0031,663,478.70
1005125.38424,753.0053,256,210.80
1006148.98429,806.0064,033,008.59
100788.29253,619.0022,393,159.18
1008111.12373,308.0041,482,686.80
1009120.90162,075.0019,594,636.15
1010111.49358,944.0040,017,365.50
101184.97262,488.0022,303,475.59
TOTAL105.643,484,639.00283,143,798.25

Share and Cite

MDPI and ACS Style

López-Serrano, P.M.; López Sánchez, C.A.; Solís-Moreno, R.; Corral-Rivas, J.J. Geospatial Estimation of above Ground Forest Biomass in the Sierra Madre Occidental in the State of Durango, Mexico. Forests 2016, 7, 70. https://doi.org/10.3390/f7030070

AMA Style

López-Serrano PM, López Sánchez CA, Solís-Moreno R, Corral-Rivas JJ. Geospatial Estimation of above Ground Forest Biomass in the Sierra Madre Occidental in the State of Durango, Mexico. Forests. 2016; 7(3):70. https://doi.org/10.3390/f7030070

Chicago/Turabian Style

López-Serrano, Pablito M., Carlos A. López Sánchez, Raúl Solís-Moreno, and José J. Corral-Rivas. 2016. "Geospatial Estimation of above Ground Forest Biomass in the Sierra Madre Occidental in the State of Durango, Mexico" Forests 7, no. 3: 70. https://doi.org/10.3390/f7030070

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop