Aquino Ramirez M MC Forestales 2014
Aquino Ramirez M MC Forestales 2014
Aquino Ramirez M MC Forestales 2014
CAMPUS MONTECILLO
ESTIMACIN DE VOLUMEN,
BIOMASA Y CARBONO AREO EN
TRES ESPECIES ARBREAS DE
PLUMA HIDALGO, OAXACA
T E S I S
PRESENTADA COMO REQUISITO PARCIAL
PARA OBTENER EL GRADO DE:
MAESTRO EN CIENCIAS
2014
La presente tesis titulada: ESTIMACIN DE VOLUMEN, BIOMASA Y CARBONO
AREO EN TRES ESPECIES ARBREAS DE PLUMA HIDALGO, OAXACA,
realizada por el alumno: MARTN AQUINO RAMREZ, bajo la direccin del Consejo
Particular indicado, ha sido aprobada por el mismo y aceptada como requisito parcial para
obtener el grado de:
MAESTRO EN CIENCIAS
EN CIENCIAS FORESTALES
CONSEJO PARTICULAR
CONSEJERO:
ASESOR:
ASESOR:
Al Consejo Nacional de Ciencia y Tecnologa (CONACYT) por la beca otorgada para los
estudios de maestra.
Al Colegio de Postgraduados (COLPOS) por la oportunidad de seguir estudiando en esta
excelente institucin acadmica.
Al Dr. Alejandro Velzquez Martnez por su amistad brindada y toda la confianza
depositada en m, asi como el gran apoyo a la elaboracin y direccin de esta investigacin,
adems de su infinita paciencia mostrada conmigo.
Al Dr. Juan Francisco Castellanos Bolaos por su amistad, excelentes consejos, estimulo
acadmico y brindarme la oportunidad de adentrarme a la investigacin forestal en el
estado de Oaxaca, asi como el apoyo para la culminacin exitosa de esta investigacin.
Al Dr. Hctor de los Santos Posadas por su apoyo, enseanzas y sugerencias para mejorar
sustancialmente este trabajo de investigacin.
Al Dr. Jorge D. Etchevers Barra por sus valiosas aportaciones a este trabajo y apoyo en la
fase de laboratorio.
A la M.C. Haidie Hervert Zamora por sus comentarios, revisiones y sugerencias para
enriquecer el escrito de esta investigacin.
Al Ing. Alberto Perez Mariscal por la autorizacin y la oportunidad de desarrollar la
presente investigacin en la finca La Cabaa en Pluma Hidalgo, Oaxaca.
A la brigada en la fase de campo, Sr. Rogelio, Sr. Cano, Sr. Agustn, por las enseanzas
en el bosque tropical.
A mis amigos de Oax por su aliento a continuar estudiando: Elizabeth, Maribel, Gaby,
Xchitl, Jess, Rene, sin olvidar a los amigos y compaeros del CP: Lissette Hervert,
Wenceslao Santiago, Rolando Tlaxcala, Esteban Blas, Luis ngel Martnez, Magnolia
Mndez, Berenice List y German Bautista.
A los profesores del CP que influyeron en mi formacin acadmica en esta aventura,
Doctores, Patricia Hernndez, Armando Gmez, Antonio Trinidad, Juan Ignacio
Valdez, Miguel A. Lpez y Sergio Prez.
iii
DEDICATORIA
Dedico esta tesis a mi familia, especialmente a mi madre Margarita Ramrez Lpez por
sus consejos, amor y confianza, y a mi padre Martn Aquino Santiago por ensearme
sentido de compromiso conmigo mismo, a los dos por su incondicional apoyo.
A mis hermanos: Reyna alhel y ngel por su cario, confianza y el apoyo brindado
durante mi formacin profesional.
iv
CONTENIDO
v
2.4.1 Biomasa area y densidad especfica .................................................................... 31
2.4.2 Grupos de especies................................................................................................ 33
2.4.3 Modelos de biomasa ............................................................................................. 35
2.4.4 Particin de biomasa ............................................................................................. 38
2.5 CONCLUSIONES ...................................................................................................... 42
2.6 AGRADECIMIENTOS .............................................................................................. 42
2.7 LITERATURA CITADA............................................................................................ 42
CAPTULO III ..................................................................................................................... 47
VOLUMEN COMERCIAL Y VARIACIN EN LA CONCENTRACIN DE CARBONO
EN TRES ESPECIES ARBREAS TROPICALES DE OAXACA, MXICO ................. 47
3.1 RESUMEN .................................................................................................................. 47
3.2 INTRODUCCIN ...................................................................................................... 48
3.3 MATERIALES Y MTODOS ................................................................................... 50
3.3.1 rea de estudio ..................................................................................................... 50
3.3.2 Modelos de regresin para estimar el volumen comercial.................................... 50
3.3.3 Determinacin de la concentracin de carbono .................................................... 52
3.3.4 Anlisis estadstico de la concentracin de carbono ............................................. 53
3.3.5 Modelo para estimar la cantidad de carbono ........................................................ 55
3.3.6 Factores de expansin de biomasa y carbono ....................................................... 55
3.4 RESULTADOS Y DISCUSIN................................................................................. 56
3.4.1 Variaciones en la concentracin de carbono ......................................................... 56
3.4.2 Estimacin del volumen comercial ....................................................................... 61
3.4.3 Estimacin de carbono areo total ........................................................................ 64
3.4.4 Estimacin de biomasa y carbono a partir del volumen ....................................... 65
3.5 CONCLUSIONES ...................................................................................................... 70
3.6 LITERATURA CITADA............................................................................................ 70
CAPITULO IV ..................................................................................................................... 73
4.1 CONCLUSIONES GENERALES .............................................................................. 73
4.2 RECOMENDACIONES GENERALES ..................................................................... 73
vi
NDICE DE CUADROS
Cuadro 1.1.Promedio de la biomasa viva y total de los principales ecosistemas terrestres del
mundo. (Fuente: Houghton et al., 2009). ............................................................ 7
Cuadro 1.2. Principales ecuaciones alomtricas desarrolladas para la estimacin de biomasa
en los diferentes tipos de vegetacin en Mxico. .............................................. 15
Cuadro 2.1 Modelos ajustados para estimacin de biomasa total (Bt). ................................ 30
Cuadro 2.2 Principales parmetros dendromtricos, pesos secos y densidad especfica de la
madera de los rboles derribados....................................................................... 32
Cuadro 2.3 Resultados de la prueba de adicionalidad para la agrupacin de especies con el
modelo: Bt 0 1Z1 2 Z 2 0 1Z1 2 Z 2 Dn2 H . ................... 34
Cuadro 2.4 Comparativo y estadsticos de bondad de ajuste de los modelos para estimar la
Bt. ....................................................................................................................... 35
Cuadro 2.5 Parmetros y estadsticos de los modelos ajustados para estimar la Bt. ............ 36
Cuadro 2.6 Resultados de ajuste para estimar la particin de biomasa del rbol en las tres
especies arbreas. .............................................................................................. 39
vii
NDICE DE FIGURAS
Figura 2.1 Distribucin de la Bt por componentes estructurales (hojas, ramas, corteza, fuste,
tocn) en las diferentes categoras diamtricas de los rboles muestreados. .... 32
Figura 2.2 Densidad especfica de la madera (g cm-3) en las tres especies arbreas. Valores
con diferente letra son significativamente diferentes (Tukey's HSD; p
0.0001). .............................................................................................................. 33
Figura 2.3 Biomasa total (Bt) en funcin al Dn2H en las tres especies arbreas. ................ 35
Figura 2.4 Ajuste de las ecuaciones para la estimacin de Bt en los dos grupos de especies
tropicales. ........................................................................................................... 37
Figura 2.5 Error relativo en la estimacin de la Bt en los diferentes dimetros del rbol. ... 37
Figura 2.6 Modelos ajustados para estimar biomasa en los componentes estructurales del
rbol: a) tocn, b) fuste, c) corteza, d) ramas y e) hojas. ................................... 40
Figura 2.7 Distribucin de biomasa entre componentes estructurales del rbol. ................. 41
viii
Figura 3.7. Modelo ajustado para estimar el carbono areo total (Ct) en los dos grupos de
especies. ............................................................................................................. 65
Figura 3.8. Estimacin de la biomasa area total (Bt) a partir del volumen comercial (vol)
en las tres especies arbreas. ............................................................................. 67
Figura 3.9. Estimacin del carbono areo total (Ct) a partir del volumen comercial (vol) en
las tres especies arbreas. .................................................................................. 67
Figura 3.10. FEB y FEC (kgm-3) variables para los diferentes tamaos del rbol en a)
Grupo 1 y b) Grupo 2. ....................................................................................... 69
ix
ESTIMACIN DE VOLUMEN, BIOMASA Y CARBONO AREO EN TRES
ESPECIES ARBREAS DE PLUMA HIDALGO, OAXACA
RESUMEN GENERAL
x
ABOVEGROUND VOLUME, BIOMASS AND CARBON ESTIMATION IN THREE
TREE SPECIES IN PLUMA HIDALGO, OAXACA
SUMMARY
Pools and storage potential of biomass and carbon are unknown in most tropical forests of
Mexico, so it is important to make accurate estimates throughout the development of
allometric equations in order to know both spatial and temporal variations in the forest
ecosystems. This study was conducted in Pluma Hidalgo, Oaxaca, Mexico, in the
physiographic region of the Sierra Madre del Sur, in order to develop equations to estimate
commercial wood volume, above-ground biomass and total above-ground carbon in three
tree species: Cupania dentata DC., Alchornea latifolia Sw., and Inga punctata Willd.
Samples of the structural aerial components were taken throughout destructive sampling on
48 trees for measurement of biomass and carbon concentrations. Two different groups were
obtained throughout an additive test, where tree species were grouped by using a regression
model. In the adjusted models for volume, above-ground biomass and above-ground total
carbon, the independent variables were the DBH and total tree height; in addition, indicator
variables were included to define the group effect, these models obtained high values of
2
determination coefficient adjustment ( R Adj > 0.97). Differences in the amount of above-
ground total biomass among groups were attributed to wood specific density and tree
morphology. Biomass of the tree structural components (stump, stem, bark, branches, and
leaves) were estimated as an equation system throughout partitioning factors using SUR
(Seemingly Unrelated Regressions) regression to estimate regression parameters; getting a
logical consistency among the functions of the components of biomass, because additive
property is obtained. Inter and intra specific average carbon concentration (CC), in the three
species showed significant differences (p0.0001). In addition, equations were developed
to estimate both biomass expansion factors (BEF) and carbon expansion factors (FEC)
variables, throughout the ratio of the models of total aerial biomass and commercial volume
with bark, in order to use in forest inventory. Greatest biomass accumulation of the three
species is concentrated mainly in descendent order in the branches, stem, bark, stump, and
leaves.
xi
CAPTULO I
INTRODUCCIN GENERAL
1
carbn y gas natural (utilizados para generar energa para mover a vehculos, electricidad
para las industrias, hogares, etc.), de incendios, deforestacin y actividades agrcolas; el
metano y el xido nitroso, resultado de mltiples actividades humanas (IPCC, 2001; Karl y
Trenberth, 2003).
En el 2011, las concentraciones de estos GEI fueron de 391 ppm1 de dixido de carbono,
1803 ppb de metano y 324 ppb de xido nitroso excediendo los niveles pre-industriales
alrededor del 40%, 150% y 20 % respectivamente (IPCC, 2013). Aproximadamente la
mitad de las emisiones actuales estn siendo absorbidas por el ocano y los ecosistemas
terrestres pero esta absorcin es sensible al clima y a las concentraciones de CO2 en la
atmsfera (Betts et. al., 1997; Cao and Woodward, 1998b; Cox et al., 2000), el incremento
atmosfrico del CO2 puede afectar la eficiencia de absorcin por parte de estos y por
consiguiente, al clima (Friedlingstein et al., 2001).
Motivado por este rpido incremento del CO2 atmosfrico, a partir de la revolucin
industrial, se han venido realizando diversas investigaciones por cientficos internacionales,
analizado el papel de los componentes individuales del sistema de la tierra en el ciclo
global del carbono (Falkowski et al., 2000). Las emisiones anuales netas de CO2
procedentes de la combustin de combustibles fsiles y la produccin de cemento fueron de
8.3 (7.6-9.0) GtC2 ao-1promediadas entre 2002-2001y de 9.5 (8.7-10.3) GtC ao-1 en 2011,
un 54% por encima del nivel correspondiente a 1990. Las emisiones anuales netas de CO2
antropognicas por el cambio del uso del suelo fueron de 0.9 (0.1-1.7) GtC ao-1 en
promedio durante 2001-2011 y de 1750 al 2011, las emisiones de CO2de la combustin de
combustibles fsiles y produccin de cemento se han disparado a 375 (345-405) GtC a la
atmsfera, mientras la deforestacin y el cambio del uso del suelo son estimados en 180
(100-260) GtC, resultando emisiones acumulativas antropognicas de 555 (470-640) GtC.
De esta acumulacin antropognica de emisiones de CO2, 240 (230-250) GtC estn
acumuladas en la atmsfera, 155 (125-185) GtC han sido absorbida por el ocano y 160
(70-250) GtC han sido acumulados en ecosistemas terrestres naturales (IPCC, 2013).
1
ppm: partes por milln o ppb: partes por billn, es la proporcin del nmero de molculas de gas con respecto al nmero de
molculas de aire seco. Por ejemplo, 300 ppm corresponde en promedio 300 molculas de gas por un milln de molculas de
aire seco.
2
1 GtC: una gigatonelada de carbono, es igual a 1015 gramos de carbono. Esto corresponde a 3.667 GtCO2.
2
Ante tal situacin, se han realizado una serie de acuerdos internacionales entre los pases
para reducir sus emisiones de CO2 o buscar sumideros para reducir las tasas de
concentracin en la atmsfera de este gas de GEI entre los principales se encuentra en
primer lugar el protocolo de Kyoto en 1997 en el Marco de la Convencin de las Naciones
Unidas para el Cambio Climtico. En l se da reconocimiento a la importancia de los
ecosistemas terrestres en el ciclo global del carbono, donde se recomienda la proteccin,
mejora y cuantificacin de los sumideros de la biosfera terrestre por las emisiones de CO2
antropognicas (Cramer et al., 2001).
Los ecosistemas terrestres desempean un papel importante en el ciclo global del carbono,
son un gran reservorio y regulador de las concentraciones de bixido de carbono (CO2) y
metano (CH4). Estos ecosistemas terrestres ejercen influencias sobre el clima mediante el
intercambio de carbono, energa e intercambio de agua (IPCC, 2001), desempeando un
papel fundamental en la mitigacin del cambio climtico contribuyendo a la reduccin de
las emisiones de CO2.
En los ltimos 200 aos las actividades humanas han alterado significativamente el ciclo de
carbono (Falkowski et al., 2000). La vegetacin y el suelo presente en los ecosistemas
terrestres desempean un papel importante en el ciclo del carbono y en las variaciones de
CO2 atmosfrico (Cao y Woodward, 1998a; Cramer et al., 2001), requiriendo una mejor
comprensin y definicin de la dinmica del balance de carbono a diferentes escalas, tanto
para la interpretacin de las variaciones observadas en el intercambio carbono entre la
atmsfera-biosfera, como para la evaluacin de las polticas para mitigar las emisiones
antropognicas de CO2 (Fung et al., 1997). Por lo anterior, es necesario realizar
investigacin respecto al potencial que se tienen de los sumideros (o fuentes) de los
ecosistemas terrestres para CO2 atmosfrico requiriendo la cuantificacin del intercambio
neto de carbono entre la vegetacin, suelo y la atmsfera (Cao y Woodward, 1998a).
Debido a que los ecosistemas presentan heterogeneidad, se torna complejo, el obtener una
alta precisin en la cuantificacin del ciclo de carbono en los ecosistemas terrestres y la
prediccin de los impactos por el cambio climtico global.
3
La fotosntesis de las plantas, la respiracin y la descomposicin del carbono orgnico del
suelo, estn involucrados en el ciclo de carbono en los ecosistemas terrestres. Por lo tanto,
la vegetacin y el suelo poseen un papel importante en el ciclo global del carbono y las
variaciones de CO2 atmosfrico. En consecuencia, los cambios en el tamao de la
vegetacin y en los reservorios de carbono en el suelo, pueden afectar significativamente el
CO2 atmosfrico (Cao y Woodward, 1998a).
4
materia orgnica muerta es menor que la productividad primaria neta (PNE >0) (Brown,
1997a).
1.1.3 Biomasa
5
Los ecosistemas terrestres cuentan casi exclusivamente con la energa del sol (Kimmins
,1997) para realizar sus funciones fisiolgicas como el crecimiento y el metabolismo de sus
organismos residentes(no as muchas plantaciones comerciales que son fertilizadas; estos
son energa que se agrega al sistema).Las plantas son relativamente poderosos factores de
biomasa gracias a la luz solar que reciben, abasteciendo de tal forma a organismos
superiores de cadenas de alimentacin con energa y bloques estructurales de construccin
de vida (Gough, 2012).Las plantas son los productores primarios terrestres, constituyndose
en los organismos que manufacturan mediante la fotosntesis, molculas nuevas orgnicas
como los carbohidratos y lpidos para formar materiales orgnicos naturales. Estos nuevos
bloques de compuestos orgnicos proveen energa que generalmente es accesible para los
organismos hetertrofos, mismos que consumen las tasas de molculas orgnicas
producidas. De esta forma, los productores primarios son un vehculo esencial para trasferir
la energa del sol a los consumidores incluyendo a los humanos, con combustible para
conducir su metabolismo mientras proveen compuestos esenciales que contienen carbono
que forman los bloques de clulas vivas (Gough, 2012). La biomasa de plantas y la materia
orgnica del suelo contiene cerca de 2200 Gt de carbono a una profundidad de 2 m, dos
veces ms que la cantidad que se encuentra en la atmsfera, los flujos de carbono entre el
sistema de vegetacin-suelo y la atmsfera son 10 veces ms grande que las emisiones de
CO2 de combustibles fsiles (Schimel et al., 1995).
La biomasa como peso seco representa alrededor del 50% de carbono en el ciclo de
carbono, la biomasa de plantas leosas se centra principalmente en los bosques ya que estos
sostienen alrededor del 70-90% de la biomasa area y subterrnea terrestre (Cairns et al.,
1997). La cantidad de biomasa en los ecosistemas terrestres del mundo varan desde 385 a
650 Pg3C (Houghton et al., 2009) (Cuadro 1.1).
3
1 Pg: un pentagramo, es igual a 1015g o 1 Gt.
6
Cuadro 1.1.Promedio de la biomasa viva y total de los principales ecosistemas
terrestres del mundo. (Fuente: Houghton et al., 2009).
Promedio de la
Total de
6 densidad de
Tipo de ecosistema rea (10 ha) biomasa
biomasa
(Pg)
(Mgha-1)
Bosques tropicales 1750, 1850 680, 350 390, 190
Bosques templados 1040 280 270
Bosques boreales 1370 110 83
Bosques templados+bosques 2410, 2450 390, 185 160, 75
boreales
Tundra rtica 560 4 7
Terrenos con arbustos mediterrneos 280 34 120
Terrenos de cultivo 1350 8 6
Sabanas y pastizales tropicales 2760 160 57
Pastizales tropicales 1500 12 8
Desierto 2770 20 7
Hielo 1550 0 0
Total 14930, 15070 1300, 733 87,51
La materia orgnica del suelo, de manera global, sostiene de dos a tres veces ms carbono
que la biomasa, sin embargo, mucho del carbono en el suelo esta fsicamente y
qumicamente protegido y no es fcilmente oxidado. En contraste, la biomasa
particularmente la area, es vulnerable al fuego, tala, cambio del uso del suelo, tormentas,
plagas etc. y por tal razn este carbono es fcilmente desprendido a la atmsfera (Davidson
y Janssens, 2006).
7
comparar atributos estructurales y funcionales de un ecosistema forestal, a travs de un
amplio rango de condiciones ambientales (Brown, 1997b).
8
resultado de los disturbios y su recuperacin. Por otro lado, a escala del paisaje (1-100 km)
conformado por diferentes rodales, la densidad de la biomasa vara a travs del espacio
consecuencia de las diferencias en las edades (Houghton et al., 2009).Por lo anterior, se han
desarrollado mtodos indirectos para estimar la densidad de biomasa. En donde, la tcnica
emprica ms comn se basa en ecuaciones alomtricas, las cuales se desarrollan a partir de
la toma de muestras destructivas, que permite la estimacin de la densidad de biomasa de
manera ms fcil, con el uso de variables dasomtricas, tales como el dimetro normal y la
altura total del rbol (Brown, 1997b).
( )
9
encuentra arriba del fuste comercial, corteza, ramas gruesas y delgadas, as como el follaje
(Somogyi et al., 2007).
( ) ( ) ( ) ( )
( ) ( )
( ) ( )
( ) ( ) ( ( )) ( ( )) ( ).
10
donde es la biomasa total o de algn componente del rbol, una variable del rbol,
parmetro del modelo y trmino del error. Las principales variables utilizadas
comnmente del rbol como se mencion anteriormente, son el Dn, Dn2, H, Dn2H,
adems, de la edad y longitud de la copa viva (LCL). El modelo 1) corresponde a un
modelo de regresin lineal mltiple que puede ser ajustado por el procedimiento de
estimacin de mnimos cuadrados ordinarios; el modelo 2) ecuaciones de regresin no
lineal necesita de procedimientos interactivos para la estimacin de los parmetros,
mientras que para el modelo 3) las ecuaciones de regresin no lineal normalmente son
transformadas a ecuaciones de regresin lineal (error aditivo) mediante el uso de logaritmos
en ambos lados de la ecuacin. De esta forma los parmetros de la ecuacin pueden
fcilmente ser estimadas por el procedimiento de mnimos cuadrados ordinarios (Parresol,
1999).
De acuerdo con Bi et al. (2004) cuando se desarrollan ecuaciones para los componentes de
biomasa del rbol frecuentemente estos son estimados por separado sin tomar en cuenta dos
situaciones: 1) la correlacin inherente entre los componentes de la biomasa medidos en la
misma muestra de rboles y 2) la restriccin lgica entre la suma de las predicciones de los
componentes del rbol y la prediccin total del rbol.
11
El primer procedimiento implica el mejor ajuste de las funciones de regresin para cada
componente de manera independiente, as la funcin de regresin para la biomasa total es
igual a la suma de cada una de las funciones de regresin individuales de los k
componentes, el sistema el sistema de ecuaciones es el siguiente:
( )
( )
.
.
.
( )
( ) () ( )
Donde la covarianza ( ) ( ) ( )
12
Bajo el supuesto de independencia de los componentes, la varianza de es la suma de
las varianzas de las .
( ) ()
El tercer procedimiento descrito por Parresol (1999) asegura la aditividad de las ecuaciones
de biomasa forzando los coeficientes de regresin de cada componente para estimar la
biomasa total; esta aditividad se logra mediante el ajuste de contrastes de los coeficientes de
regresin. La regresin lineal generalizada (SUR por sus siglas en ingls Seemingly
unrelated regression) puede ser utilizada si se asumen correlaciones contemporneas entre
ecuaciones de componentes para la estimacin de los parmetros y con ello se logra obtener
una varianza baja. La estructura de un sistema de ecuaciones de modelos de biomasa
ajustado por SUR est definido por:
( )
( )
.
.
.
( )
( )
donde para cada ecuacin se tiene un vector de componentes de biomasa , una matriz
que contiene las variables independientes , un vector de parmetros y un vector de error
aleatorio .
13
ajuste en la mayora de las ecuaciones desarrolladas recure al uso de logaritmos (natural o
base 10) para transformar el modelo potencial a un lineal (Cuadro 2), logrando con ello
reducir la heterocedasticidad, as mismo, en otros trabajos utilizan algoritmos de
optimizacin lineal para resolver el problema de mnimos cuadrados para modelos no
lineales. Para mejorar la precisin en las estimaciones de la cantidad de biomasa es
indispensable el desarrollo de nuevas ecuaciones a nivel de especie o grupos de especies
donde se tome en cuenta la variabilidad de las condiciones ambientales y estructurales por
el sitio, adems de la localizacin geogrfica. De la misma manera, es importante tomar en
cuenta las variaciones inter-especficas e intra-especficas en las concentraciones de
carbono en las especies arbreas para desarrollar modelos para estimar la cantidad de
carbono en el arbolado.
14
Cuadro 1.2. Principales ecuaciones alomtricas desarrolladas para la estimacin de biomasa en los diferentes tipos de
vegetacin en Mxico.
Tipo de Dn
Autor Especies N Ecuacin
vegetacin (cm)
Pinus devoniana Lindl., P. Plantacin mixta en B=Exp(-1.705+1.936 ln(Dn))
Mndez-Gonzlez et al.,2011 13-27 40
pseudostrobus Lindl. Guanajuato B=Exp(1.250+0.157Dn)
Avendao-Hernndez et al.,2009 Abies religiosa(H.B.K.) Schl. et Cham. Bosque de Abies 6.5-79 26 B=0.0713 Dn2.5104
15
Pinus durangensis Martinez, Pinus
arizonica Engelmann.
Plantacin de hule en
Rojo-Martnez et al.,2005 Hevea brasiliensis Mll. Arg. 10-40 28 ln(B)=3.1426+2.69273 ln(Dn)
Papaloapan, Oaxaca
Bosque templado de
Figueroa-Navarro et al.,2010 Pinus patula Schltdl. et Cham. 2.6-56.9 18 B=5.338+0.018635(Dn2 H)
Pinus patula
Caesalpinia eriostachys Benth.,
Thouinia paucidentata Radlk.,
Apoplanesia paniculata Presl.,
Lonchocarpus constrictus Pitt.,
Caesalpinia coriaria (Jacq.) Willd.,
Ipomea Wolcottiana Rose., log(B)=-0.5352+0.9996 log(BA)
Lonchocarpus lanceolatus Benth., log(B)=-0.4492+0.9972 log(BA)+0.7127
Martinez-Yrizar et al.,1993 Opuntia excelsa Sanchez-Mejorada., Selva baja caducifolia >3 191 log(WSG)
Lysiloma microphylla Benth., log(B)=-0.7590+0.9011 log(BA)+0.5715
Heliocarpus pallidus Rose., Bursera log(WSG)+0.5654 log(H)
excelsa (HBK) Engl., Caesalpinia
sclerocarpa Standl., Cordia
elaeagnoides DC., Jatropha
malacophylla Standl., Guettarda
elliptica Sw.
Bosque templado de
Daz-Franco et al.,2007 Pinus patula Schl. et Cham 6.7-64.1 25 B=0.0357Dn2.6916
Pinus patula
Bosque templado de >7.5- B=22.3476-4.9470
Pimienta De la Torre et al.,2007 Pinus cooperi Blanco 633
pinos <67.5 (Dn)+0.4911(Dn2)+0.0039(Dn2H)
ln (B)=-1.62+2.12 ln(Dn)
ln (B)=-3.51+2.67 ln(Dbase)
ln (B)=-4.04+4.00 ln(Dn)-0.29 (ln(Dn))2
ln (B)=-4.60+2.93(Dbase)
Guazuma ulmifolia, Trichos permum
+0.45(ln(Dbase))2-0.12(ln(Dbase ))3
Douterlungne et al.,2013 mexicanum,Inga vera, Ochroma Bosque neotropical 0-39.8 103
ln (B)=-2.45+2.30ln(Dn)
pyramidale
ln (B)=-4.56+2.82ln(Dbase)+0.29(ln(Dbase))2
-0.09(ln(Dbase))3
ln (B)=-2.82+2.42 ln(Dn)
ln (B)=-3.75+2.65ln(Dbase)-0.01(ln(Dbase))3
N: nmero de observaciones; B:biomasa total; Dn: dimetro normal; WSG: densidad bsica de la madera; Dbase: dimetro de la base;
ln=logaritmo natural; log: logaritmo base 10.
16
1.1.5 Caractersticas de las especie bajo estudio
Se distribuye en la vertiente del Golfo de Mxico en Tamaulipas, Puebla, San Luis Potos,
Veracruz, Tabasco y Quintana Roo. En la vertiente del Pacfico, se le encuentra en Jalisco,
Oaxaca y Chiapas. Altitud: 0 a 1,600 m pero preferentemente alrededor de los 400 a 800 m.
A menudo abundante a lo largo de pendientes abruptas de barrancas y caadas. Prospera en
acahuales, relictos de selva media, claros de bosque, cafetales y a lo largo de caminos.
Crece en suelos derivados de material calizo o gneo y metamrfico, pero siempre en suelos
17
con buen drenaje y buena capacidad de retencin de humedad. Suelos: negro, rendzina con
caliza, caf-oscuro pedregoso, oscuro arcilloso. Es rbol perennifolio, de 8 a 25 m (hasta 30
m) de altura, con un dimetro normal de hasta 60 cm. Copa redonda o irregular. Hojas
simples, alternas, grandes, elpticas; lminas de 7 a 18 cm de largo por 9 a 28 cm de ancho,
ovadas a orbiculares, margen crenado; verde oscuras en el haz y verde plidas en el envs
con nerviacin amarillenta prominente. Especie dioica. Los frutos son cpsulas biloculares
de hasta 1.3 cm de ancho, perianto, estilo y estigmas persistentes. Semillas 2 de 6 mm de
dmetro, ligeramente comprimidas, maduran de marzo a mayo. Su madera relativamente
suave se utiliza en construccin rurales ligeras, se emplea mucho como sombra de
plantaciones de caf en la zona de Pluma Hidalgo, Oaxaca (Pennington y Sarukhn, 2005).
18
1.2 OBJETIVOS
1.2.1 General
1.2.2 Especficos
1.3 HIPTESIS
El volumen comercial, la biomasa area total y el carbono areo total en las tres
especies tropicales pueden ser modelados cada uno de ellos, bajo un mismo modelo
mediante el uso de variables indicadoras y con ello obtener modelos ms prcticos
y precisos.
Las diferencias en la cantidad de biomasa y carbono de las tres especies son
influidas por la densidad de madera presente.
La concentracin de carbono vara en funcin del componente estructural del rbol
(fuste, ramas, corteza) y a la especie arbrea.
19
1.4 LITERATURA CITADA
Betts, R. A., P. M. Cox, S. E. Lee and F. I. Woodward. 1997. Contrasting physiological and
structural vegetation feedbacks in climate change simulations. Nature 387:796-799.
Bi, H., J. Turner and M. J. Lambert. 2004. Additive biomass equations for native eucalypt
forest trees of temperate Australia. Trees18:467-479.
Brown, S. 1997a. Los bosques y el cambio climtico: el papel de los terrenos forestales
como sumideros de carbono. In Actas del XI Congreso Mundial Forestal: Recursos
Forestales y Arboles pp. 13-22.
Brown, S. 1997b. Estimating biomass and biomass change of tropical forests: a primer.
FAO Forestry Paper 134. Food and Agriculture Organization of the United Nations
(FAO), Rome, Italy, pp. 55.
Brown, S., A. J. R. Gillespie and A. E. Lugo. 1989. Forest with applications to forest
inventory data. Forest Science, 35(4):881-902.
Cao, M. And F. I. Woodward. 1998 a. Net primary and ecosystem production and carbon
stocks of terrestrial ecosystems and their responses to climate change. Global
Change Biology, 4:185-198.
20
Dixon, R. K., S. Brown, R. A. Houghton, A. M. Solomon, M. C. Trexler and J.
Wisniewski. 1994. Carbon pools and flux of global forest ecosystems. Science
(Washington), 263(5144):185-1991.
FAO. 2010. Evaluacin de los recursos forestales mundiales. FAO, Roma, pp.44.
FAO. 2011. Situacin de los bosques en el mundo. 2011. Estudio FAO Montes 163. Roma
(Italia).Disponible en: http://www.fao.org/docrep/013/i1757s/i1757s.pdf.
Helgerson, O. T., K. Cromack, S. Stafford, R. E. Miller and R. Slagle. 1988. Equations for
estimating aboveground components of young Douglas-fir and red alder in a coastal
Oregon plantation. Can. J. For. Res. 18:1082-1085.
Houghton, J. T., L.G. Meira Filho, B.A. Callander, N. Harris, A. Kattenberg, and K.
Maskell (eds.). 1996. Climate Change 1995: The Science of Climate Change:
contribution of working group I to the second assessment report of the
Intergovernmental Panel on Climate Change. Vol. 2. Cambridge University Press,
Houghton, R.A. 2005. Aboveground forest biomass and the global carbon balance. Global
Change Biology 11: 945-958.
Houghton, R.A., F. Hall., and S. J. Goetz. 2009. Importance of biomass in the global
carbon cycle. Journal of Geophysical Research: Biogeosciences, 114(G2):1-13.
IPCC (Intergovernmental Panel on Climate Change). 2001. Climate Change 2001: The
Scientific Basis, Cambridge University Press, Cambridge, 892 pp.
21
IPCC (Intergovernmental Panel on Climate Change). 2013: Summary for Policymakers. In:
Climate Change 2013: The Physical Science Basis. Contribution of Working Group
I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change
[Stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, A.
Nauels, Y. Xia, V. Bex and P.M. Midgley (eds.)]. Cambridge University Press,
Cambridge, United Kingdom and New York, NY, USA.
IPCC, 2012: Glossary of terms. In: Managing the Risks of Extreme Events and Disasters to
Advance Climate Change Adaptation [Field, C.B., V. Barros, T.F. Stocker, D. Qin,
D.J. Dokken, K.L. Ebi, M.D. Mastrandrea, K.J. Mach, G.-K. Plattner, S.K. Allen,
M. Tignor, and P.M. Midgley (eds.)]. A Special Report of Working Groups I and II
of the Intergovernmental Panel on Climate Change (IPCC). Cambridge University
Press, Cambridge, UK, and New York, NY, USA, pp. 555-564
IPCC. 2007. Fourth assessment report, climate change 2007. Synthesis report: summary for
policymakers. Ginebra (Suiza). IPCC (disponible tambin en www.ipcc.
ch/pdf/assessment-report/ar4/syr/ar4_syr_spm.pdf).
Kaitaniemi, P., 2004. Testing allometric scaling laws. J. Theor. Biol. 228:149-153.
Kimmins, J.P.1997.Balancing Act. Environmental issues in forestry. 2nd ed., UBC Press,
Vancouver, Canada.Pp.33-97.
Navar, N. 2009. Allometric equations for tree species and carbon stocks for forests of
northwestern Mexico. For. Ecol. Manage. 257:427-434.
Parresol, B. R. 1999. Assessing tree and stand biomass: a review with examples and critical
comparison. For Sci. 45:573-593.
22
Pennington, T. D., y J. Sarukhn.2005. rboles de Mxico. Manual para la identificacin
de las principales especies. Universidad Nacional Autnoma de Mxico. Fondo de
Cultura Econmica. Mxico.523 pp.
Pilli, R., T. Anfodillo and M. Carrer. 2006. Towards a functional and simplified allometry
for estimating forest biomass. For. Ecol. Manage.237:583-593.
Zianis, D.2008. Predicting mean aboveground forest biomass and its associated variance.
For. Ecol. Manage. 256:1400-1407.
23
CAPTULO II
24
2.2 INTRODUCCIN
Los trabajos realizados para generar ecuaciones de las especies nativas del trpico
mexicano son escasos, aun cuando, muchas de estas especies, son utilizadas para programas
de plantaciones forestales comerciales y de restauracin, as como para establecerse en los
sistemas agroforestales. Algunos estudios relevantes para generar ecuaciones alomtricas
en la estimacin de biomasa, son los realizados en el bosque tropical seco en la costa de
Jalisco (Martinez-Yrizar et al, 1992), la pennsula de Yucatn (Cairns et al, 2003) y en
Sinaloa (Navar, 2009); mientras que en el bosque tropical hmedo son los desarrollados en
Oaxaca (Acosta-Mireles et al., 2002) y en Chiapas (Douterlungne et al., 2013); as como en
plantaciones comerciales de hule (Rojo-Martnez et al.,2005; Monroy y Navar, 2004),
donde ajustaron distintos modelos lineales y no lineales, incluyendo como variables
25
predictoras al dimetro normal, altura total, dimetro de la base, rea basal y gravedad
especifica de la madera.
Los objetivos de este estudio fueron: 1) desarrollar modelos alomtricos para estimar la
cantidad de biomasa area total en Cupania dentata DC., Alchornea latifolia Sw. e Inga
punctata Willd., especies arbreas nativas de la Sierra Sur de Oaxaca, de manera conjunta
mediante el uso de variables indicadoras y, 2) desarrollar un sistema de identificacin de la
distribucin de la biomasa en los distintos componentes del rbol.
26
2.3.2 Muestreo destructivo de rboles
El fuste fue seccionado en trozas de 1 m de longitud para obtener el peso verde; as mismo,
cada componente del rbol se pes in situ por separado para obtener su peso verde con una
bscula digital marca TECNOCOR modelo IPEN de 100 kg, con una precisin de 0.5 kg.
El peso seco de cada componente fue obtenido mediante la relacin peso seco: peso fresco,
para lo cual se obtuvieron submuestras de madera consistentes en secciones transversales
(rodajas) de fuste de aproximadamente 4-6 cm de grosor.
Las primeras dos rodajas fueron obtenidas a una altura de 0.30 y 1.30 m sobre el nivel suelo
y despus, en la parte superior de cada troza, separando la corteza de cada una de ellas; as
mismo, se obtuvieron 4 muestras para las ramas y 4 muestras de 0.5 kg para las hojas.
Todas las submuestras se pesaron inmediatamente para evitar prdidas de humedad, se
colocaron en bolsas de papel con sus respectivas etiquetas, se transportaron y se pusieron
96 h en una estufa de secado de madera (Nyle modelo L500) a 70 C o hasta lograr un peso
constante. El peso verde y el peso seco de las submuestras pequeas se obtuvieron en una
bscula electrnica (TECNOCOR modelo COR-30) con una precisin de 5 g. La relacin
peso seco: peso fresco promedio se utiliz para obtener la biomasa (peso seco) en el tocn,
las ramas y las hojas, mientras que para el fuste y la corteza dicha relacin se obtuvo de
cada una de las submuestras de las trozas.
27
La biomasa total area (Bt) de cada rbol (kg) corresponde a la suma de la biomasa de cada
uno de sus componentes, expresado as:
La densidad especfica de la madera ( ) fue obtenida por la relacin del peso seco con el
volumen verde de la muestra de madera (g cm-3), mientras que el volumen verde se
determin por el mtodo de desplazamiento de agua, previamente saturadas las secciones
para evitar su contraccin. El peso anhidro se obtuvo colocando las muestras en una estufa
de secado por 48 h o hasta a alcanzar un peso constante.
La biomasa total de las tres especies se analiz de manera conjunta para conocer la
posibilidad de agruparlas como una sola especie o grupos de especies, y poder definir en
una sola ecuacin la contribucin relativa de las mismas, para ello se desarroll una prueba
de adicionalidad simple por medio de un anlisis de regresin lineal incorporando variables
indicadoras (variables dummy) (Draper y Smith, 1981; Montgomery y Peck, 1982).
Bt 0 1 ( Dn2 H )
1, si pertenece a I. punctata
Z2
0, de otra forma
donde Bt es la biomasa rea total del rbol (kg), Dn es el dimetro normal (cm), H es la
altura total (m), Z1 y Z2 corresponde a las variables indicadoras,0, 0 son los parmetros de
regresin, 1,2, 1, 2 los parmetros de regresin debidos a la adicionalidad y el error
28
aleatorio. Se consider a Cupania dentata como la especie base en el modelo por presentar
el mayor nmero de observaciones. Para agrupar a las especies como un solo grupo o
grupos de especies, se consider la homogeneidad de interceptos y de pendientes
(H0=1=2=1=2=0).
Una vez agrupadas las especies se ajustaron los modelos alomtricos mediante
procedimientos de regresin. Tres modelos utilizados en la estimacin de la biomasa area
total (Brown et al., 1989; Parresol, 1999, 2001) fueron seleccionados, los cuales expresan la
Bt como una funcin del dimetro normal y la altura del rbol:
Bt f ( Dn, H )
En el Cuadro 2.1 se presentan los tres modelos los cuales incluyen a las variables
indicadoras que definen el efecto del grupo de especies. Con la finalidad de comparar estos
2
modelos y seleccionar el mejor, se consider el coeficiente de determinacin ajustado( R Adj )
n
bi bi n p
2
2
RAdj 1 i n1
b bi 2 n 1
i 1
i
b b
n
2
RMSE i i (n p)
i 1
29
b b n p
n
2
i i
CV % i 1
100
b
n
bi bi
E %
100
n
i 1 bi
b b
n
i i
S i 1
n
Modelo Ecuacin
Bt: Biomasa rea total del rbol (kg); Dn: dimetro normal (cm); H altura total (m); 0 , 0
3 parmetros de regresin de C. dentata; 1 , 1 parmetros de adicionalidad debido a A.
latifolia; 2 , 2 parmetros de adicionalidad debido a I. punctata; Z1, Z2 variables
indicadoras; Exp base de los logaritmos neperianos y el error aleatorio.
La distribucin de biomasa en los componentes estructurales del rbol, se ajust como una
funcin directa de la Bt y un factor de particin (F) expresado como:
Bcomponente f componenteBt, F
El cual permiti derivar el siguiente sistema de ecuaciones de particin, con las siguientes
funciones:
Btocn g Bt,
30
B fuste hBt,
Bcorteza iBt,
Bramas j Bt,
Bhojas k Bt,
Los parmetros de regresin en los modelos de particin de biomasa fueron estimados con
el mtodo SUR (Seemingly Unrelated Regressions), ya que se supone que existe una
correlacin entre las ecuaciones individuales de los componentes, para as lograr una
varianza baja en el sistema y por lo tanto un mejor ajuste que aumente la sensibilidad de las
pruebas de hiptesis. La aditividad en el sistema se garantiz con la definicin de
restricciones al modelo mediante el desarrollo de contrastes de los coeficientes de regresin
(Parresol, 1999). El procedimiento MODEL de SAS/ETS (SAS Institute Inc., 2004) se
us para el ajuste de los modelos.
31
tabla FITECMA (Sotomayor-Castellanos, 2008) como densidad baja (0.201-0.40 g cm-3); y
las ltimas dos especies, como densidad media (0.401-0.60 g cm-3).
Promedio, Mnimo-Mximo
Especie
H (m) Dn (cm) Peso seco del rbol (kg) (g cm-3)
H: altura total del rbol (m); Dn: dimetro normal (cm); densidad especfica de la
madera.
2000
Hojas
1600 Ramas
Biomasa (kg)
Corteza
1200 Fuste
Tocn
800
400
0
< 7.5
10
15
20
25
30
< 7.5
10
15
20
30
35
45
50
< 7.5
10
15
20
25
30
32
(Tamarit-Urias y Lpez-Torres, 2007). Sin embargo, el valor promedio para I. punctata es
inferior a 0.63 0.10 g cm-3, valor reportado por Segura et al. (2006), diferencia que
probablemente se deba a las prcticas de manejo implementadas en el sistema agroforestal
de caf con esta especie en esta zona, similares a las implementadas en San Ramn,
Nicaragua (Segura et al., 2006).
Figura 2.2 Densidad especfica de la madera (g cm-3) en las tres especies arbreas.
Valores con diferente letra son significativamente diferentes (Tukey's
HSD; p 0.0001).
2.4.2 Grupos de especies
Los resultados obtenidos en la prueba de adicionalidad (Cuadro 2.3), indican que las
especies C. dentata e I. punctata pueden agruparse como una sola especie (Grupo 1), es
decir, no existen diferencias significativas entre interceptos ( 0 = 1 , p = 0.6578), ni entre
33
parmetros de la pendiente es estadsticamente significativa ( 0 , p .0001), aunque no
de las especies, indican que hay mayor biomasa area por rbol en las especies de C
.dentata e I. punctata, que para la especie de A. latifolia, cuando presentan el mismo
dimetro y altura (Figura 2.3), estas diferencias son estadsticamente significativas. Por lo
tanto, la cantidad de biomasa asociada a la densidad especfica de la madera difiere entre
grupos de especies, aun cuando los rboles se encuentran bajo las mismas condiciones (luz,
temperatura, tipo de suelo, agua). Las especies del Grupo 1, forman parte del estrato
inferior (12-20 m) de la estructura del bosque tropical de Pluma Hidalgo, mientras la
especie del Grupo 2 se encuentra en el estrato de rboles altos (20-35 m) (Acosta, 1997).
Parmetros estimados
Especie
Smbolo Valor Valor t Pr > |t|
0 -14.7815 -1.15 0.2561
C. dentata
0 0.0395 21.76 <0.0001
1 17.3487 1.05 0.3020
A. latifolia
1 -0.0190 -10.15 <0.0001
2 8.0586 0.45 0.6578
I. punctata
2 0.0009 0.36 0.7239
34
2000
1200
800
A. latifolia
400 C. dentata
I. punctata
0
0 20000 40000 60000 80000
Dn H (cm m)
2 2
Figura 2.3 Biomasa total (Bt) en funcin al Dn2H en las tres especies arbreas.
El uso de variables indicadoras en un modelo lineal, tiene la ventaja de que se pueden
realizar de forma directa pruebas de hiptesis con el mtodo de la suma extra de cuadrados
(Montgomery, 2005), y con ello, comparar los parmetros de los modelos de estimacin de
biomasa debido al efecto de la especie.
Los modelos ajustados para la estimacin de la Bt en los dos grupos (Cuadro 2.4),
presentaron altos valores en el coeficiente de determinacin ajustado ( R Adj
2
>0.98). Por esta
razn, la variacin explicada por el modelo, cuando se incluye adems del Dn, la H o una
combinacin (Dn2H) es alta; es decir, la precisin en la estimacin de la biomasa mejora
significativamente. Estos mismos resultados han sido obtenidos en otros trabajos para la
estimacin de biomasa en rboles tropicales (Brown et al. 1989; Brown, 1997; Arajo et al.
1999; Nogueira et al., 2008; Feldpausch et al, 2012).
RMSE(kg)
2
Modelo CV (%) S (%) E (kg) R Adj
35
RMSE: raz de los cuadrados medios del error; CV: coeficiente de variacin; S sesgo
promedio; E : error relativo estndar promedio; R Adj
2
coeficiente de determinacin
ajustado.
Los dos grupos de especies presentan diferencias estadsticas, por lo que los modelos
generales (Cuadro 2.1), tienden a reducirse y, el efecto de la variable indicadora Z, es
significativo (Cuadro 2.5). De acuerdo a los valores bajos de RMSE, CV, E y mejor
comportamiento, se seleccion el Modelo 3 de la forma Schumacher y Hall (Figura 2.4),
aun cuando el S fue ligeramente mayor al del Modelo 1; as mismo se pueden considerar
satisfactorios los estadsticos de ajuste de los Modelos 1 y 2. Todos los coeficientes de
regresin en los modelos fueron altamente significativos, excepto 0 del Modelo 1 (Cuadro
2.5).
Cuadro 2.5 Parmetros y estadsticos de los modelos ajustados para estimar la Bt.
1 0 32.36 <0.0001
0.0392 0.0012
1 -0.0185 0.0011 -15.49 <0.0001
36
2400
2000 Grupo 1
Biomasa area total (kg) Grupo 2
1600
1200
( )
800
400
( )
0
0 10 20 30 40 50 60
Dn (cm)
Figura 2.4 Ajuste de las ecuaciones para la estimacin de Bt en los dos grupos de
especies tropicales.
En los tres modelos, los valores de error relativo ms altos corresponden a los rboles
menores a 20 cm de dimetro (Figura 2.5), parte de esta variacin en la biomasa puede
deberse a que en estas categoras diamtricas los rboles crecen bajo una fuerte
competencia por la luz con los rboles vecinos (Nvar, 2009).
2
Modelo 1
Modelo 2
1
Modelo 3
Error relativo
-1
-2
0 10 20 30 40 50 60
Dn (cm)
37
El modelo seleccionado con variables indicadoras es equivalente a tener dos modelos
especficos para cada grupo (Figura 2.4), los cuales son estadsticamente diferentes. Es
probable que las diferencias en los pesos secos sean atribuibles al efecto de la densidad de
la madera (Kuyah, et al., 2012; Fayolle et al., 2013) y a la arquitectura del rbol (Basuki et
al., 2009). As mismo, el ajuste de un modelo con variables indicadoras provoca mayor
robustez en el modelo en comparacin con el ajuste individual. Adems, el desarrollo de
ecuaciones especficas para especies locales es importante para mejorar la precisin en la
estimacin de biomasa, porque las condiciones ambientales o estructurales influyen en la
alometra del rbol, dimetro: altura, de la misma manera que la variacin originada por la
localizacin geogrfica (Feldpausch et al., 2011; Banin et al., 2012; Fayolle et al., 2013).
En todas las ecuaciones de particin de biomasa (Cuadro 2.6), se incorpor el efecto de las
variables indicadoras, por lo tanto, es necesario sumar o restar el efecto de dichas variables
sobre los parmetros de regresin para conocer la prediccin de biomasa en algn
componente estructural especfico del rbol (Figura 2.6). Los valores del coeficiente de
determinacin ajustado en cada modelo fueron altos ( R Adj
2
>0.90), excepto para la corteza
38
Cuadro 2.6 Resultados de ajuste para estimar la particin de biomasa del rbol en las
tres especies arbreas.
Todos los parmetros ajustados fueron significativos (p 0.0001); EE: error estndar de
los parmetros; RMSE: raz de los cuadrados medios del error; R Adj 2
: coeficiente de
determinacin ajustado; 11 11 : parmetros de regresin de C. dentata; 12 12 :
parmetros de adicionalidad debido a A. latifolia; 13 13 : parmetros de adicionalidad
debido a I. punctata;, Z1,Z2 variables indicadoras.
39
Figura 2.6 Modelos ajustados para estimar biomasa en los componentes estructurales
del rbol: a) tocn, b) fuste, c) corteza, d) ramas y e) hojas.
40
La distribucin de biomasa entre los componentes fue muy variable debido al tamao del
rbol. C. dentata destina en promedio la mayor cantidad de biomasa a las ramas 51.64 % y
al fuste 40.69 %, seguidos de la corteza 5.59 %, el tocn 4.12% y 2.68% en las hojas. En A.
latifolia las ramas y el fuste tuvieron la mayor distribucin de biomasa con 44.07 % y
38.53% respectivamente, seguidos de la corteza 8.17 %, las hojas 6.51% y 2.70 en el
tocn; y en I. punctata tambin las ramas y el fuste presentaron la mayor cantidad de
biomasa con 55.73 % y 33.65%, seguidos del tocn, las hojas y la corteza con 3.77 %, 3.53
% y 3.29 %, respectivamente.
41
2.5 CONCLUSIONES
Se logr separar a las especies en dos grupos por medio de la adicionalidad en un modelo
de regresin con variables indicadoras. La diferencia estadstica entre grupos de especies es
atribuible a la densidad especfica de la madera y a la morfologa del rbol.
En todos los casos, las ramas son el componente estructural del rbol que acumula mayor
biomasa seguidas del fuste, la corteza, el tocn y las hojas, y esto es consistente en C.
dentata, A. latifolia e I. punctata.
2.6 AGRADECIMIENTOS
Este estudio fue financiando por la Fundacin Produce Oaxaca A.C. a travs del proyecto
Almacenamiento de carbono y conservacin de la biodiversidad vegetal en el rea de
Pluma Hidalgo-Loxicha, Oaxaca. Agradecimientos especiales a todas las personas de la
finca La Cabaa por los esfuerzos dedicados, as como a los revisores que hicieron posible
la elaboracin de este trabajo.
Acosta, C.S. 1997. Afinidades fitogeogrficas del bosque mesfilo de montaa de la zona
de Pluma Hidalgo, Oaxaca, Mxico. Polibotnica.6:25-39.
42
tropical forest architecture? Testing environmental, structural and floristic drivers.
Global Ecol. Biogeogr. 21:1179-1190.
Basuki, T. M., P. E. van Laake, A.K. Skidmore, and Y. A. Hussin. 2009. Allometric
equations for estimating the above-ground biomass in tropical lowland Dipterocarp
forests. For. Ecol. Manage. 257:1684-1694.
Bi, H., J. Turner, and M. J. Lambert. 2004. Additive biomass equations for native eucalypt
forest trees of temperate Australia.Trees.18:467-479.
Brown, S. 1997. Estimating biomass and biomass change of tropical forests: A primer.
FAO. For. Pap.134, Rome, pp 55.
Cairns, M.A., I. Olmsted, J. Granados, and J. Argaez. 2003. Composition and aboveground
tree biomass of a dry semi-evergreen forest on Mexicos Yucatan Peninsula. For.
Ecol. Manage.186: 125-132.
Castellanos, F. J., A. Velzquez M., J.J. Vargas H., C. Rodrguez F., y A. M. Fierros G.
1996. Produccin de biomasa en un rodal de Pinus patula Schl. et Cham.
Agrociencia 30: 123-128.
Clutter, J.L., J.C. Forston, L.V. Pienaar, G.H. Brister, and R.L. Bailey. 1983. Timber
Management: a Quantitative Approach. John Wiley and Sons. New York. 333p.
43
Draper, N. R., and H. Smith. 1981. Applied Regression Analisys. 2da. ed. John Wiley &
Sons. New York. 709 p.
Fayolle, A., J.L. Doucet, J.F. Gillet, N. Bourland, and P. Lejeune. 2013. Tree allometry in
Central Africa: Testing the validity of pantropical multi-species allometric equations
for estimating biomass and carbon stocks. For. Ecol. Manage.305:29-37.
Feldpausch, T.R., L. Banin, O.L. Phillips, T.R. Baker, S.L. Lewis, C.A. Quesada, K.
Affum-Baffoe, E. Arets, N.J. Berry, and M. Bird. 2011. Height-diameter allometry
of tropical forest trees. Biogeosciences 8:1081-1106.
INEGI. 2005. Conjunto de datos vectoriales del uso del suelo y vegetacin, Serie III
(continuo nacional), escala 1:250 000. Instituto Nacional de Estadstica, Geografa e
Informtica, Aguascalientes, Mxico.
44
INEGI. 2009. Prontuario de informacin geogrfica municipal de los Estados Unidos
Mexicanos. http://www3.inegi.org.mx/sistemas/mexicocifras/datos
geograficos/20/20071.pdf. Consultado el 18 de septiembre 2013.
Nvar, J. 2009. Allometric equations for tree species and carbon stocks for forests of
northwestern Mexico. For. Ecol. Manage.257: 427434.
Nvar, J., E. Mndez, A. Njera, J. Graciano, V. Dale, and B. Parresol. 2004. Biomass
equations for shrub species of Tamaulipan thornscrub of Northeastern Mexico.
Journal of Arid Environments 59(4): 657-674.
Nvar-Chidez, J.J. 2010. Biomass allometry for tree species of Northwestern Mexico.
Tropical and Subtropical Agroecosystems. 12: 507-519.
Parresol, B. 1999. Assessing tree and stand biomass: a review with examples and critical
comparisons. For. Sci. 45:573-593.
45
Rutishauser, E., F. Nooran, Y. Laumonier, J. Halperin, Rufiie, K. Hergoualch, and
L.Verchot. 2013. Generic allometric models including height best estimate forest
biomass and carbon stocks in Indonesia. For. Ecol. Manage. 307:219-225.
SAS Institute Inc. 2004. SAS/ETS 9.1 Users Guide. Cary, NC. 2416 p.
Segura, M., M. Kanninen, and D. Surez. 2006. Allometric models for estimating
aboveground biomass of shade trees and coffee bushes grown together. Agrofor.
Syst. 68(2):143-150.
46
CAPTULO III
3.1 RESUMEN
El objetivo de este estudio fue medir la concentracin de carbono y desarrollar modelos que
permitan estimar el volumen comercial, la biomasa y carbono areo total en tres especies
arbreas tropicales de la Sierra Sur de Oaxaca. La concentracin de carbono (CC) se
determin en una muestra de 27 rboles de las especies Cupania dentata DC., Alchornea
latifolia Sw. e Inga punctata Willd., en el fuste, las ramas y las hojas. La CC promedio
inter-especfica e intra-especfica en las tres especies tropicales mostraron diferencias
significas. La concentracin de carbono promedio ponderado (CCPP) en C. dentata fue de
46.290.29%, mientras que en A. latifolia de 45.920.33%20.33% y en I. punctata de
45.430.3%. Los modelos de volumen, biomasa y carbono fueron ajustados mediante
regresin lineal y no lineal con valores del coeficiente de determinacin ajustado superiores
a 0.95. Las ecuaciones para estimar los factores de expansin de biomasa (FEB) y los
factores de expansin de carbono (FEC) variables fueron desarrollados a partir de la razn
del modelo de biomasa area total y el volumen comercial con corteza del rbol. Es
importante considerar las variaciones inter-especficas e intra-especficas de la CC en los
rboles tropicales para mejorar la precisin en las estimaciones del contenido de carbono.
Los factores de expansin variables de biomasa y carbono variables propuestos permiten la
estimacin rpida del contenido de biomasa y carbono presente en los rboles, cuando solo
se dispone del volumen comercial.
47
3.2 INTRODUCCIN
A pesar de que el clima en nuestro planeta es dinmico y cambiante por su ciclo natural,
actualmente, parte de esa variacin, es consecuencia de las actividades humanas; tales
como, el aumento de las concentraciones de los Gases de Efecto Invernadero (GEI)
principalmente en el dixido de carbono (CO2), metano (CH4) y xido nitroso (N2O), los
cuales han provocado cambios importantes en las escalas climticas de la tierra (IPCC,
2007). La emisin de CO2 a la atmsfera, gas de efecto invernadero, se ha intensificado
por la quema de combustibles fsiles y el cambio del uso del suelo (Dixon et al., 1994).
Este continuo incremento en la concentracin atmosfrica de CO2, prev que conducir a
cambios significativos en el clima (Houghton et al., 1996; Friedlingstein et al., 2001).
48
Chidez, 2010). Los FEB dependen principalmente de la especie, caractersticas
ambientales, edad, estructura del rodal y calidad de la estacin (Aalde et al., 2006; Silva-
Arredondo y Nvar-Chidez, 2009); as mismo, los FEC pueden ser afectados por el
crecimiento en la distribucin de la biomasa y la concentracin de carbono (Bert y Danjon,
2006; Fukatsu, et al, 2008).
Considerando lo anterior, los objetivos planteados para este estudio fueron: a) medir y
comparar la concentracin de carbono inter-especfica e intra-especfica en las especies
arbreas nativas de la Sierra Sur de Oaxaca, Cupania dentata DC., Alchornea latifolia Sw.
e Inga punctata Willd., y b) desarrollar ecuaciones para estimar los factores de expansin
de biomasa y de carbono variables a nivel de rbol para su aplicabilidad a inventarios
forestales.
49
3.3 MATERIALES Y MTODOS
El fuste comercial fue definido hasta que el tronco se empez a bifurcar en la copa de rbol,
el cual una vez derribado, fue seccionado en trozas de 1 m de longitud. En cada troza se
midieron los dimetros superior e inferior, as como la longitud con cinta mtrica. El
volumen comercial fue obtenido por el mtodo de trozas traslapadas (Bailey, 1995). Antes
del derribo, se procedi a la medicin del dimetro normal (Dn) a una altura de 1.30 m
sobre el nivel del suelo con cinta diamtrica; posterior al derribo, se midi la altura total del
rbol (H) con cinta mtrica.
50
El volumen de las tres especies se analiz de manera conjunta para conocer la posibilidad
de agruparlas como una sola especie o grupos de especies, y poder definir en una sola
ecuacin el comportamiento de las mismas, para ello se desarroll una prueba de
adicionalidad en el modelo de la variable combinada incorporando variables indicadoras
(variables dummy). El modelo fue el siguiente, considerando a I. punctata como especie
base:
Vol 0 1Z1 2 Z 2 0 1Z1 2 Z 2 Dn2 Hk
1, si pertenece a C. dentata
Z1
0, de otra forma
1, si pertenece a A. latifol ia
Z2
0, de otra forma
Los modelos fueron ajustados por la tcnica de regresin lineal y no lineal de acuerdo a su
estructura, incorporando variables indicadoras:
Variable combinada.. Bt 0 1Z1 2 Z 2 0 1Z1 2 Z 2 Dn2 H
Variable combinada logartmica Bt Exp0 1Z12Z2 Dn2 H 0 1Z1 2Z2
las variables indicadoras; k=/40000, Exp base de los logaritmos neperianos ( =2.7182) y
es el error aleatorio. Para comparar estos modelos y seleccionar el mejor, se consider el
2
coeficiente de determinacin ajustado ( R Adj ) que expresa la variabilidad explicada por el
51
modelo tomando en cuenta el nmero de parmetros, la raz de los cuadrados medios del
error (RMSE) que define la precisin de las estimaciones y el coeficiente de variacin en %
(CV), mismos que se calcularon de la siguiente forma:
n
yi y i n p
2
2
R Adj 1 i n1
yi y i n 1
2
i 1
n
RMSE y
i 1
i yi (n p)
2
y i y i
2
n p
CV % i 1
100
y
Donde yi , y i , y i es el volumen comercial del fuste del rbol observado, predicho y medio;
De los 48 rboles derribados para obtener la biomasa area total en el captulo II, solo se
seleccionaron 27 para realizar el anlisis de la concentracin de carbono. Para ellos se
conformaron tres categoras diamtricas para cada una de las especies arbreas: la primera
corresponde a menores a 10 cm, la segunda, mayores a 10 cm y menores a 20 cm y la
ltima correspondiente a mayores a 20 cm de dimetro. En cada categora diamtrica se
seleccionaron tres rboles de manera aleatoria, de donde se obtuvieron tres muestras
(probetas) de las ramas, tres muestras de 0.5 kg en hojas y una muestra en el tronco a una
altura de 1.30 m sobre el nivel de suelo (dimetro normal), los cuales permitieron preparar
una muestra compuesta para cada uno de los componentes estructurales del rbol. De todas
las muestras obtenidas en ramas y tronco, se extrajo una submuestra de 20 g que incluy la
albura, duramen y corteza, mientras que las hojas, incluyeron flores y frutos, los cuales
fueron finamente picados con herramientas de carpintera. Posteriormente, para lograr la
homogeneidad en el material de las submuestras, se llevaron a un molino de material
52
vegetal Thomas-Wliley Modelo 4 hasta que pasaran una malla de 20 mm y por ltimo, para
conseguir partculas ms finas del material vegetal se colocaron en un molino General
Electric Mod. 5KH39QN5525 hasta pasar por una malla de 1 mm. Todas las submuestras
fueron llevadas a una estufa de secado a 70 C por 24 horas, o hasta alcanzar un peso
constante, para eliminar la humedad absorbida del ambiente. En todas las submuestras se
obtuvo una porcin de 30 mg, las cuales fueron pesadas en una balanza analtica con una
precisin de 0.0001g, posteriormente, se colocaron en un analizador automtico de carbono
Shimadzu TOC 5050-A para determinar la concentracin de carbono por la tcnica de
combustin seca en el Laboratorio de Fertilidad de Suelo y Qumica Ambiental del Colegio
de Postgraduados. Las curvas de calibracin de la concentracin de carbono fueron
construidas con sacarosa, despus de cada 10 muestras determinadas aproximadamente 20
mg de sacarosa fue medido para compensar el error de medicin.
n
CCPP Ci Pi
i 1
53
Cuadro 3.1 Factores y niveles para el anlisis de varianza en la concentracin de
carbono (CC).
Factor Niveles
E1 C. dentata
Especie (E) E2 A. latifolia
E3 I. punctata
CET Tronco
Componente estructural (CE) CER Ramas
CEH Hojas
CD1 Menor a 10 cm
Categora diamtrica (CD) CD2 Entre 10-20 cm
CD3 Mayor a 20 cm
54
3.3.5 Modelo para estimar la cantidad de carbono
carbono (%). La suma del contenido carbono en todos los componentes del rbol
corresponde al carbono total almacenado en el rbol. Se seleccion el mismo modelo
utilizado para estimar la biomasa area total (Bt) en el captulo II, para ajustar el contenido
de carbono en cada una de las especies arbreas tropicales, expresado de la siguiente forma:
Ct Exp 0 Dn1 H 2
Para estimar la biomasa area total del rbol (Bt) o el carbono areo total (kg) a partir del
volumen comercial (m3) se ajustaron dos modelos lineales que expresan el vol como una
funcin directa de la Bt y el Ct, bajo este supuesto; el primer modelo multiplica el vol por
un factor de expansin de biomasa (FEB) variable, expresado de la siguiente forma:
Bt 0 1 I1 2 I 2 0 1 I1 2 I 2 vol
Ct 0 1 I1 2 I 2 0 1 I1 2 I 2 vol
55
1, si pertenece a Alchornea latifolia
I1
0, de otra forma
1, si pertenece a I. punctata
I2
0, de otra forma
Donde Bt, Ct, vol es la biomasa area total del rbol (kg), carbono areo total (kg) y
volumen de fuste comercial (m3) respectivamente, 0 , 1 , 2 , 0 , 0 , 1 , 2 , 0 , 1 , 2
56
(p>0.05), con valores en el fuste (45.56 %), ramas (46.39 %) y hojas (46.07 %) muy
similares.
50
Fuste
Ramas a
a) 49 Hojas ab
48 abc abc
CC (%) bc
47 c c
c
c
46
45
44
C1 C2 C3
Categoras diamtricas
Fuste
49 Ramas a
Hojas
b) 48
ab ab
47 ab
ab
CC (%)
46 b b
b b
45
44
43
42
C1 C2 C3
Categoras diamtricas
49
Fuste
a
48 Ramas
c)
Hojas
47 a a
a a
CC (%)
46 a
a a
a
45
44
43
42
C1 C2 C3
Categoras diamtricas
57
As mismo, en I. punctata la CC promedio en hojas (46.4%) y fuste (45.48%) no son
estadsticamente diferentes (p>0.05), ni tampoco hay diferencias significativas entre la CC
promedio de las ramas (45.19%) y el fuste. La CC fue significativamente diferente entre las
especies bajo estudio (F2,78=25.19, p0.001) al excluir el factor componente estructural y
categora diamtrica. La CC promedio de A. latifolia fue de 46.010.21 % y en I. punctata
de 45.690.20% los cuales no fueron estadsticamente diferentes (p>0.05) de acuerdo a la
prueba de comparacin de medias de Tukey's HSD; mientras que la CC promedio en C.
dentata fue mayor a las otras dos especies con 46.710.19% y por lo tanto estadsticamente
diferente (p<0.05) (Figura 3.3). As mismo, el contenido de carbono promedio ponderado
(CCPP) para C. dentata fue de 46.290.29 %, en A. latifolia de 45.920.33% y en I.
punctata de 45.430.3 (Cuadro 3.2). Los valores promedios de las concentraciones de
carbono obtenidas en los componentes del rbol y los promedios ponderados por especie en
este estudio fueron significativamente diferentes (p0.0001) al valor de 50%
frecuentemente utilizado en la literatura para estimar la cantidad de carbono en la biomasa
seca.
49
Fuste
Ramas a
48 Hojas
47 a
b a a
CC (%)
b
46 a ab
b
45
44
43
C. dentata A. latifolia I. puctata
Especies
58
Cuadro 3.2. Concentracin de carbono (CC, %) entre componentes estructurales del
rbol y concentracin de carbono promedio ponderado (CCPP, %) en las
tres especies tropicales.
Componente estructural
Especie CD CCPP
Tronco Ramas Hojas
CD1 45.89 45.62 46.90 45.87
C. dentata CD2 46.93 46.60 47.92 46.84
CD3 46.05 46.12 48.32 46.15
46.290.32 46.110.28 47.710.42 46.290.29
CD: categora diamtrica; CD1: rboles con un dimetro menor a 10 cm, CD2 arboles con
dimetro de 10-20 cm, CD3 rboles con dimetros mayores a 20 cm; CCPP concentracin
de carbono promedio ponderado; EE error estndar.
Los valores obtenidos de la concentracin de carbono en las tres especies son similares a
los reportados en otros estudios para latifoliadas, por ejemplo, Elias y Potvin (2003),
registraron que la concentracin de C tuvo una variacin entre 44.4% y 49.4% en 32
especies tropicales de Panam; Figueroa et al. (2005), determinaron la concentracin de C
en 6 especies tropicales de la Sierra Norte de Oaxaca con valores de 47.3 a 51.3%; por otro
lado, Gayoso y Guerra (2005) determinaron valores de 34.86% a 48.31% para 16 especies
de bosques nativos de Chile. La concentracin de C en latifoliadas es inferior comparado
con las conferas. Se ha encontrado que los valores en conferas varan de 47.21 %-55.20%
mientras que en latifoliadas de 46.27%- 49.97% (Lamlom y Savidge, 2003). Por otro lado,
la concentracin de C puede estar en un intervalo de 47 a 59 % como una funcin de los
compartimentos del rbol o de las especies (Elias y Potvin, 2003; Lamlom y Savidge,
2003).
59
Comparar las concentraciones de carbono con otros estudios puede ser complicado debido a
las diferencias en el tipo de vegetacin, tcnicas de manejo, mtodo de muestreo y la
tcnica empleada para la determinacin de carbono; as como edad del rbol y condiciones
climticas (Bert y Danjon, 2006). Las variaciones en las concentraciones de C en los
diferentes componentes estructurales del rbol y entre especies, se debe a que la biomasa de
las plantas est conformada de molculas orgnicas con proporciones variables de carbono;
es decir, la biomasa de la madera est conformada principalmente de la lignina (1835%) y
la holocelulosa (6575%), esta ltima compuesta de -celulosa y hemicelulosa; adems, en
la madera se encuentran otros materiales como extractos orgnicos y minerales orgnicos
(cenizas) que representan aproximadamente de 4-10% de la biomasa; los extractos se
componen de grasas, ceras, alcaloides, protenas, compuestos fenlicos simples y
compuestos, azcares simples, pectinas, mucilagos, gomas, resinas, terpenos, almidones,
glucsidos, saponinas y aceites esenciales (Bert y Danjon, 2006).
Por otro lado, se realiz una prueba de t-student para muestras pareadas para comparar la
cantidad de carbono de los rboles derribados en el captulo II considerando la variabilidad
de la concentracin de carbono intra- especfica e inter- especfica, y el valor genrico de
50% frecuentemente utilizado en la estimacin del contenido de carbono. El resultado
obtenido de la prueba, muestra diferencias significativas entre los valores de la cantidad de
60
carbono considerando la variabilidad intra- especifica e inter-especfica de la concentracin
de C y la cantidad de carbono comparado con el valor de 50% sin importar la especie ni
componente estructural del rbol (t=-4.7936, gl=47, p0.0001). El error relativo promedio
de la cantidad de C en C. dentata es de 7.880.24 %, en A. latifolia de 8.680.26% y del
10.490.22 en I. punctata cuando se realiza estimaciones de la cantidad de C ignorando la
variacin intra-especfica en la concentracin de carbono (Figura 3.4), en otras palabras,
utilizar el valor genrico de 50% en las estimaciones de C genera una sobre-estimacin en
la cantidad de carbono en las especies bajo estudio.
11
Error realtivo (%)
10
7
C. dentata A. latifolia I. punctata
Especies
61
rechazar estos resultados sera necesario obtener una muestra ms amplia de rboles de las
tres especies, cubriendo un amplio intervalo de dimetros y realizar las mismas pruebas de
hiptesis.
Promedio, Mnimo-Mximo
Especie N
Dn (cm) H (m) Vol (m3)
C. dentata 17 20.64(5.92-54.75) 14.36(8.12-27.20) 0.1427(0.0029-0.4528)
A. latifolia 17 15.55(3.34-29.60) 13.99(6.90-20.15) 0.2812(0.0099-1.4143)
I. punctata 14 14.28(4.71-27.82) 14.52 (9.70-20.20) 0.1217(0.0089-0.4791)
N: nmero de observaciones; Dn: dimetro normal (cm); H: altura total del rbol (m);
Parmetros estimados
Especie
Smbolo Valor EE Pr > |t|
0 0.005186 0.0264 0.8457
I. punctata
0 0.3952 0.0562 <0.0001
1 0.00466 0.0376 0.9021
C. dentata
1 -0.05931 0.0740 0.4279
2 0.0467 0.0343 0.1816
A. latifolia
2 -0.15151 0.0575 0.0117
Los resultados de ajuste de los modelos seleccionados para estimar el volumen de fuste
comercial en los dos grupos de especies, presentan en general buenos estadsticos, es decir,
un coeficiente de determinacin ajustado superior a 0.91 (Cuadro 3.5); sin embargo, en los
estadsticos de RMSE y CV, el modelo Schumacher-Hall presenta los valores ms bajos,
por lo cual fue seleccionado.
62
Cuadro 3.5. Resultados de ajuste de los modelos con variables indicadoras para
estimar volumen comercial.
0 0.026474ns 0.0143
Variable
0 0.000026*** 2.568E-6 0.0775 41.73 0.9193
combinada
1 -6.53E-6*** 2.534E-6
Variable
0 8.5438*** 0.4077
combinada 0 0.795344*** 0.0444 0.0641 34.52 0.9448
logartmica 1 -0.00151 ns 0.00958
0 8.29357*** 0.3194
no significativo p>0.05.
Como se puede observar en estas dos ecuaciones, el volumen del Grupo 1 es mayor al del
Grupo 2 cuando presentan la misma altura y el mismo dimetro (Figura 3.5), estos
resultados son congruentes con los coeficientes de forma obtenidos por especie, es decir,
los rboles del Grupo1 son ms cilndricos que los del Grupo 2.
63
2
A. latifolia
1.6
Volumen (m3)
C. dentata
1.2 I. punctata
0.8
0.4
0
0 10 20 30 40 50 60
Dn (cm)
Figura 3.5. Modelo ajustado para estimar volumen de fuste comercial (m3) para los
dos grupos de especies.
3.4.3 Estimacin de carbono areo total
Figura 3.6. Resultados de ajuste de los modelos para estimar el contenido de carbono
areo (Ct).
Parmetros estimados 2
RMSE (kg) R Adj
Smbolo Valor EE
0 4.177307*** 0.1310
0 2.4872*** 0.0462
7.56 0.9969
1 -0.2113*** 0.0079
2 0.506852*** 0.0366
EE: error estndar de los parmetros; RMSE: raz de los cuadrados medios del error; R Adj
2
64
Los parmetros estimados de adicionalidad del modelo son equivalentes a tener dos
ecuaciones, uno para cada grupo de especies de la siguiente forma:
Los resultados obtenidos del ajuste del modelo confirman las conclusiones obtenidas en el
captulo II, el Grupo 1 presenta mayor cantidad de carbono (kg) que el Grupo 2 cuando
presentan el mismo dimetro y altura (Figura 3.7).
1000
Carbono areo total (kg)
A. latifolia
800 C. dentata
I. punctata
600
400
200
0
0
30 10
40 5020 60
Dn(cm)
Figura 3.7. Modelo ajustado para estimar el carbono areo total (Ct) en los dos grupos
de especies.
3.4.4 Estimacin de biomasa y carbono a partir del volumen
El modelo ajustado para estimar la biomasa area total y carbono areo total a partir del
volumen comercial presentaron buenos estadsticos de ajuste (Cuadro 3.7 y 3.8). Para
estimar la Bt (kg) por rbol se proponen los siguientes modelos para los grupos de especies:
65
Mientras que para estimar el Ct (kg) por rbol son los siguientes:
Como se puede observar en los valores de los parmetros de los modelos (Cuadro 3.7 y
3.8), estos nos proporcionan un factor de expansin de biomasa y de carbono variable
directamente. Para el Grupo1 el FEB es de 1358.5 kg de peso seco de la biomasa area
total por rbol y un FEC es de 631.25 de peso seco de carbono areo total por rbol por
cada metro cubico del volumen del fuste comercial; mientras que los valores de FEB y FEC
en el Grupo 2 son de 998.56 y 462.36 kgm-3 respectivamente. Las diferencias en los valores
del FEB y FEC se observan en las Figuras 3.8 y 3.9.
EE: error estndar de los parmetros; RMSE: raz de los cuadrados medios del error; R Adj
2
66
Cuadro 3.8. Resultados de ajuste de la estimacin del carbono a partir del volumen
comercial agrupando las especies.
Parmetros estimados RMS 2
Ecuacin R Adj
(kg)
Smbolo Valor EE Pr > |t|
0 -9.91 5.763 0.0924
Bt 0 ( 0 1 I1 )vol 0 631.24 34.52 <0.0001 30.36 0.9506
1 -168.87 34.51 <0.0001
2
EE: error estndar de los parmetros; RMSE: raz de los cuadrados medios del error; R Adj :
coeficiente de determinacin ajustado; 0 , 0 ,parmetros de regresin; 1 parmetro de
adicionalidad debido al grupo 2.
2000
Biomasa area total (kg)
A. latifolia
1600 C. dentata
I. punctata
1200
800
400
0
0 0.5 1 1.5 2 2.5
Vol (m3)
Figura 3.8. Estimacin de la biomasa area total (Bt) a partir del volumen comercial
(vol) en las tres especies arbreas.
1600
Carbon areo total (kg)
A. latifolia
C. dentata
1200 I. punctata
800
400
0
0 0.5 1 1.5 2 2.5
3
Vol (m )
Figura 3.9. Estimacin del carbono areo total (Ct) a partir del volumen comercial
(vol) en las tres especies arbreas.
67
Adems, a partir de la relacin del modelo de biomasa rea total (Bt) y el volumen de
comercial (vol) tambin es posible obtener una ecuacin de factor de expansin variable
(FEB), expresado de la siguiente forma:
Con las ecuaciones anteriores se pueden construir tablas de FEB y FEC para un
determinado dimetro y altura total del rbol (Anexo 3.1), que al multiplicar el valor
correspondiente de FEB y FEC por el volumen del rbol obtenemos la biomasa area total
y carbono areo total respectivamente. Los FEB y FEC variables pueden resultar ms
precisos en las estimaciones de biomasa que los FEB y FEC constantes, porque toman en
cuenta la variacin en la alometra de los rboles (Jalkanen et al., 2005). En los dos grupos
de especies, los FEB y FEC incrementan con el tamao del rbol (Figura 3.10), debido a
que en las especies de rboles tropicales que crecen en espacios abiertos acumulan mayor
biomasa en las ramas, es decir, la expansin de la copa es ms rpida que la del fuste
cuando el dimetro aumenta, por esta razn los valores de FEB y FEC se incrementan en
los dimetros mayores.
68
1200
a)
1000
600
FEB
400 FEB_estimados
FEC
200 FEC_estimados
0
0 10 20 30 40 50 60
Dn (cm)
2000
b)
1600
FEB-FEC (kgm-3)
1200
800 FEB
FEB_estimados
400 FEC
FEC_estimados
0
0 5 10 15 20 25 30 35
Dn (cm)
Figura 3.10. FEB y FEC (kgm-3) variables para los diferentes tamaos del rbol en a)
Grupo 1 y b) Grupo 2.
69
3.5 CONCLUSIONES
La mayor concentracin de carbono se encontr en las hojas y las ramas seguido del fuste
en C. dentata e I. punctata, por lo contrario en A. latifolia los valores de las
concentraciones fueron similares. La concentracin promedio ponderada entre especies
mostr diferencias significativas.
Bailey, R. L. 1995. Upper-stem volume from stem-analysis data and overlapping bolt
method. Canadian J. Forest Res. 25:170-173.
Basuki, T. M., P. E. van Laake, A.K. Skidmore, and Y. A. Hussin. 2009. Allometric
equations for estimating the above-ground biomass in tropical lowland Dipterocarp
forests. For. Ecol. Manage. 257:1684-1694.
Bert, D., and F. Danjon. 2006. Carbon concentration variations in the roots, stem and crown
of mature Pinus pinaster (Ait.). For. Ecol. Manage. 222:279-295.
Brown, S. 1997. Estimating biomass and biomass change of tropical forests: a primer. FAO
Forestry Paper 134. Food and Agriculture Organization of the United Nations
(FAO), Rome, Italy, pp. 55.
Cao, M., and F. I. Woodward. 1998 a. Net primary and ecosystem production and carbon
stocks of terrestrial ecosystems and their responses to climate change. Global
Change Biology 4: 185-198.
70
Clutter, J.L., J.C. Forston, L.V. Pienaar, G.H. Brister, and R.L. Bailey. 1983. Timber
Management: a Quantitative Approach. John Wiley and Sons. New York. 333p.
Houghton, J. T., L.G. Meira Filho, B.A. Callander, N. Harris, A. Kattenberg, and K.
Maskell (eds.). 1996. Climate Change 1995: The Science of Climate
Change: contribution of working group I to the second assessment report of the
Intergovernmental Panel on Climate Change. Vol. 2. Cambridge University
Press.
71
Jalkanen A, R. Mkip, G. Stahl, A. Lehtonen, and H. Petersson. 2005. Estimation of the
biomass stock of trees in Sweden: comparison of biomass equations and age-
dependent biomass expansion factors. Ann. For. Sci. 62:845-851.
Nvar-Chidez, J. J. 2009. Alometric equations and expansion factors for tropical dry trees
of eastern Sinaloa, Mexico. Tropical and Subtropical Agroecosystems 10:45-
52.
72
CAPITULO IV
Los resultados obtenidos en esta investigacin permiten proponer una metodologa para
desarrollar modelos para estimar biomasa area total y un sistema de particin de biomasa
en tres especies tropicales del municipio de Pluma Hidalgo, Oaxaca. As mismo para
estimar el volumen y el carbono areo se desarrollaron los modelos correspondientes. Es
importante destacar, que las ecuaciones propuestas son una herramienta prctica para su
uso en campo, debido a que sus variables independientes son dimetro normal y la altura
total del rbol. Para mejorar las estimaciones del contenido de carbono en las especies bajo
estudio es necesario considerar las variaciones inter-especficas e intra-especficas en la
concentracin de carbono.
73
ANEXOS
2.1 ndice de Valor de Importancia (IVI) para las especies presentes en la zona de estudio.
Homalium trichostemon S.F. Blake Palo piedra 1.196 3.359 1.826 6.381
Palo mazo 1.914 2.037 1.826 5.777
Coralillo 0.957 1.932 1.826 4.716
Saurauia serrata DC. Mamellito 1.675 0.264 2.740 4.679
Cecropia obtusifolia Bertol. Guarumbo 1.435 1.573 1.370 4.378
Juan diego 0.957 1.695 1.370 4.022
Diphysa robinioides Benth. Cuachupil 1.914 0.412 1.370 3.695
Lonchocarpus guatemalensis Benth. Palo gusano 1.675 0.164 1.826 3.665
Abutilon purpusii Standl. Yaco 0.957 1.290 1.370 3.617
Cedrela odorata L.
Cedro rojo 1.914 0.284 1.370 3.568
74
Garbanso 0.718 0.106 1.370 2.193
Cuil cuadrado 0.957 0.111 0.913 1.981
Bursera simaruba Mulato 0.478 0.558 0.913 1.950
Macuil 0.239 1.248 0.457 1.944
Caulote 0.718 0.264 0.913 1.895
Llaco de piedra 0.957 0.467 0.457 1.881
Naranjal 0.718 0.177 0.913 1.808
Cola pava 0.718 0.104 0.913 1.735
Guanabo 0.718 0.089 0.913 1.719
Chachalaco 0.718 0.073 0.913 1.704
Palo de ropa 0.478 0.072 0.913 1.464
Huaje verde 0.478 0.053 0.913 1.445
Mango 0.239 0.586 0.457 1.282
Palo naba 0.239 0.580 0.457 1.276
Cucharillo 0.239 0.436 0.457 1.132
Hierba blanca 0.478 0.110 0.457 1.045
Palo gallina 0.478 0.106 0.457 1.041
Cuasia 0.478 0.063 0.457 0.998
Zorrillo 0.478 0.055 0.457 0.990
Llagozan 0.239 0.284 0.457 0.980
Cuilocho 0.239 0.145 0.457 0.841
Cuil de cuero 0.239 0.126 0.457 0.822
Floripundio 0.239 0.109 0.457 0.805
Ciruelo 0.239 0.093 0.457 0.789
Llaco de cuero 0.239 0.041 0.457 0.737
San jose 0.239 0.041 0.457 0.737
Huaje 0.239 0.023 0.457 0.719
Cacao 0.239 0.019 0.457 0.715
Palo cano 0.239 0.014 0.457 0.710
64 100 100 100 300
75
3.1 Factores de expansin de biomasa (FEB) y Factores de expansin de carbono
variables
817.36 1069.46 1251.59 1399.32 1525.81 1637.61 1738.50 1830.91 1916.49 1996.42
30 872.43 1141.51 1335.90 1493.58 1628.60 1747.93 1855.62 1954.25 2045.59 2130.91
35 921.86 1206.19 1411.60 1578.22 1720.89 1846.98 1960.77 2064.99 2161.51 2251.67
40 966.95 1265.18 1480.64 1655.41 1805.05 1937.31 2056.66 2165.98 2267.22 2361.78
45 1008.54 1319.61 1544.33 1726.61 1882.70 2020.64 2145.13 2259.15 2364.75 2463.38
50 1047.26 1370.27 1603.62 1792.90 1954.98 2098.22 2227.49 2345.89 2455.54 2557.95
55 1083.57 1417.77 1659.21 1855.06 2022.75 2170.96 2304.71 2427.22 2540.66 2646.63
60 1117.81 1462.58 1711.65 1913.68 2086.68 2239.57 2377.55 2503.92 2620.95 2730.27
65 1150.27 1505.04 1761.34 1969.25 2147.26 2304.60 2446.58 2576.62 2697.05 2809.54
70 1181.16 1545.46 1808.64 2022.13 2204.92 2366.48 2512.28 2645.81 2769.48 2884.99
25 374.53 492.91 578.82 648.70 708.67 761.77 809.75 853.74 894.53 932.67
30 400.22 526.72 618.52 693.20 757.29 814.02 865.29 912.31 955.89 996.65
35 423.31 557.11 654.21 733.20 800.98 860.99 915.22 964.94 1011.05 1054.15
40 444.39 584.85 686.78 769.70 840.86 903.86 960.78 1012.99 1061.39 1106.64
45 463.85 610.46 716.86 803.41 877.68 943.44 1002.86 1057.35 1107.87 1155.10
50 481.98 634.32 744.87 834.81 911.99 980.31 1042.06 1098.68 1151.17 1200.24
55 498.99 656.70 771.16 864.27 944.17 1014.91 1078.83 1137.45 1191.79 1242.60
60 515.04 677.83 795.97 892.08 974.55 1047.56 1113.54 1174.04 1230.13 1282.58
65 530.26 697.86 819.49 918.44 1003.35 1078.52 1146.45 1208.74 1266.49 1320.48
70 544.76 716.94 841.89 943.55 1030.77 1108.00 1177.78 1241.78 1301.11 1356.58
76
FEB variables para el Grupo 2 (A. latifolia)
Altura total (m)
5 10 15 20 25 30 35 40 45 50
5 356.34 466.25 545.65 610.06 665.21 713.95 757.93 798.22 835.53 870.38
10 409.14 535.33 626.49 700.44 763.76 819.72 870.23 916.48 959.32 999.33
15 443.58 580.39 679.23 759.40 828.05 888.72 943.48 993.63 1040.07 1083.45
20 469.76 614.65 719.32 804.22 876.92 941.18 999.16 1052.27 1101.45 1147.39
25
Dimetro normal (cm)
491.12 642.60 752.03 840.80 916.81 983.98 1044.60 1100.13 1151.55 1199.58
30 509.30 666.38 779.86 871.92 950.74 1020.40 1083.26 1140.84 1194.17 1243.97
35 525.19 687.18 804.20 899.12 980.40 1052.24 1117.06 1176.44 1231.43 1282.79
40 539.36 705.71 825.89 923.38 1006.85 1080.62 1147.20 1208.17 1264.64 1317.39
45 552.17 722.48 845.51 945.31 1030.76 1106.29 1174.45 1236.87 1294.68 1348.68
50 563.89 737.81 863.45 965.37 1052.64 1129.77 1199.37 1263.12 1322.16 1377.31
55 574.70 751.96 880.01 983.89 1072.83 1151.44 1222.38 1287.35 1347.52 1403.72
60 584.76 765.12 895.41 1001.10 1091.60 1171.58 1243.76 1309.87 1371.09 1428.28
65 594.16 777.42 909.81 1017.20 1109.15 1190.42 1263.76 1330.94 1393.14 1451.25
70 603.00 788.99 923.35 1032.34 1125.66 1208.14 1282.57 1350.74 1413.88 1472.85
25 223.72 294.43 345.75 387.49 423.32 455.03 483.69 509.97 534.34 557.12
30 232.19 305.58 358.84 402.16 439.34 472.26 502.00 529.28 554.56 578.21
35 239.60 315.33 370.29 415.00 453.36 487.33 518.02 546.17 572.26 596.66
40 246.21 324.03 380.50 426.45 465.87 500.77 532.31 561.24 588.05 613.12
45 252.19 331.90 389.75 436.81 477.19 512.94 545.25 574.87 602.34 628.02
50 257.67 339.11 398.21 446.29 487.55 524.07 557.08 587.35 615.41 641.65
55 262.72 345.76 406.02 455.04 497.11 534.35 568.01 598.87 627.49 654.24
60 267.42 351.94 413.29 463.19 506.01 543.91 578.17 609.59 638.71 665.94
65 271.82 357.73 420.08 470.80 514.33 552.86 587.68 619.61 649.22 676.90
70 275.96 363.18 426.48 477.97 522.16 561.27 596.63 629.04 659.10 687.20
77