nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
La estimación de biomasa de los árboles es determinante en los estudios de productividad, ciclo de nutrientes y cuantificación de las contribuciones al secuestro de carbono por parte de los ecosistemas forestales (Tesfaye, Bravo-Oviedo, Bravo y Ruiz-Peinado, 2015; Meng, et al., 2017). Así mismo, el conocimiento de las cantidades de biomasa sirve como un indicador ecológico para la sostenibilidad y para caracterizar la estructura y función de los bosques (Chave et al., 2005). En los últimos años, se ha incrementado la necesidad de disponer de información más eficiente sobre la cantidad de almacenamiento de biomasa y de los flujos de carbono forestal en los diferentes ecosistemas terrestres, debido a su influencia para mitigar el cambio climático global.
Para estimar la biomasa aérea en árboles se han generado ecuaciones alométricas, las cuales se derivan de la aplicación de un muestreo destructivo y del uso de técnicas de regresión (Acosta-Mireles, Vargas-Hernández, Velázquez-Martínez y Etchevers-Barra, 2002; Vargas-Larreta et al., 2017; Quiñonez-Barraza et al., 2019). Las ecuaciones de biomasa se fundamentan en las relaciones alométricas empíricas de los atributos del árbol que son fácilmente medibles como el diámetro normal, la altura total y la densidad específica de la madera (Temesgen, Affleck, Poudel, Gray y Sessions, 2015; Eker, Poudel y Özçelik, 2017). A nivel local o del sitio, el uso de ecuaciones resulta la opción más práctica en tiempo y costo para evaluar la cantidad de biomasa aérea a partir de datos de inventario forestal (Riofrío, Herrero, Grijalva y Bravo, 2015; Daryaei y Sohrabi, 2016). Diferentes estudios señalan que las ecuaciones específicas por especie mejoran las estimaciones de biomasa que las ecuaciones generales, porque se toman en cuenta las condiciones ambientales, la densidad de la madera y la variabilidad genética en donde se desarrollaron los árboles; factores que influyen significativamente en las relaciones alométricas (Skovsgaard y Vanclay, 2013; Lin et al., 2017; Mensah, Veldtman y Seifert, 2017).
Al estimar la biomasa de más de un componente del árbol es importante considerar la propiedad de aditividad puesto que se toman en cuenta las correlaciones inherentes entre los componentes medidos y la biomasa total, mejorando la precisión en las estimaciones (Carvalho y Parresol, 2003; Dong, Zhang y Li, 2014). Con la aditividad se asegura que la suma de los valores estimados de las ecuaciones de los componentes estructurales sea igual al valor estimado obtenido de la ecuación de la biomasa aérea total (Bi, Turner y Lambert, 2004).
En México, se han realizado considerables esfuerzos para generar ecuaciones para estimar la biomasa arbórea en las distintas especies de coníferas y latifoliadas (Douterlungne, Herrera-Gorocica, Ferguson, Siddique y Soto-Pinto, 2013; Aquino-Ramírez, Velázquez-Martínez, Castellanos-Bolaños, De los Santos-Posadas y Etchevers-Barra, 2015; Rojas-García, De Jong, Martínez-Zurimendí y Paz-Pellat, 2015; Quiñonez-Barraza et al., 2019). La mayoría de las ecuaciones del género Pinus son restrictivas a especies con fines maderables (Chávez-Aguilar et al., 2016; Vargas-Larreta et al., 2017). Sin embargo, no existen ecuaciones para las especies de pinos piñoneros a pesar de su amplia distribución y plasticidad ambiental. Por lo cual, la información que se tiene de la capacidad de almacenamiento de biomasa y carbono de los bosques piñoneros es escasa.
Los bosques de Pinus cembroides Zucc. son característicos de las zonas áridas y de distribución latitudinal muy amplia sobre el altiplano del norte y centro de México (Luna-Cavazos, Romero-Manzanares y García-Moya, 2008). P. cembroides es considerada como la principal especie de pinos piñoneros del país, crece sobre laderas de montaña en suelos delgados y pedregosos con climas de tipo templado seco a templado subhúmedo, donde la época de sequía puede durar de siete a ocho meses (Constante-García, Villanueva-Díaz, Cenaro-Paredes, Cornejo-Oviedo y Valencia-Manzo, 2009). A pesar de su limitada capacidad productividad puede constituir una fuente complementaria de captura y almacenamiento de carbono (Pompa-García y Yerena-Yamallel, 2014), además suministra bienes a los pobladores locales como leña y madera de cortas dimensiones, resina, árboles de navidad, alimento e ingresos económicos por la venta de su semilla (piñón). Por su amplia adaptabilidad y resistencia a condiciones extremas ha sido utilizada para proyectos de reforestación y restauración ecológica (Ríos, De Hoogh y Návar, 2008). No obstante, las poblaciones de esta especie se ven amenazadas o en peligro de desaparecer por las actividades antropogénicas como el cambio de uso del suelo, el pastoreo, los incendios, la extracción ilegal de la madera y principalmente por la recolección excesiva de la semilla (Rosas-Chavoya, Granados-Sánchez, Granados-Victorino y Esparza-Govea, 2016). Por lo anterior, son necesarias acciones que promuevan la conservación y el manejo sustentable de P. cembroides. Una alternativa viable es su incorporación a los programas de pago de servicios ambientales como la captura de carbono, puesto que estos bosques pueden considerarse como una posible fuente con potencial de mitigación ante el cambio climático (Pompa-García y Yerena-Yamallel, 2014), aunque los montos económicos por la venta de bonos de carbono serían mínimos, debido a su baja productividad. De manera que se requiere información detallada de la cantidad de los diferentes componentes estructurales del árbol y total para estimar con precisión los almacenes de carbono en los bosques de P. cembroides.
Objetivos
Desarrollar un sistema de ecuaciones de aditividad para estimar la biomasa de los componentes estructurales por separado y la biomasa aérea total a nivel de árbol individual para Pinus cembroides Zucc. en el bosque de la comunidad agraria de Sombrerete, Cadereyta de Montes, Querétaro, México.
Materiales y métodos
Área de estudio
El estudio se realizó en el mes de octubre de 2017, en una superficie de 8 ha de las 1333 ha de los bosques naturales de la Comunidad Agraria de Sombrerete, Municipio de Cadereyta de Montes, Querétaro, México, ubicada entre los 20º 49' 80" y 20º 50' 0.05" N; y los 99º 39' 88.0" y 99º 40' 13.0" O, a 2600 m snm. El clima es de tipo semiseco (BS1kw) con lluvias en verano. La temperatura media anual se encuentra entre 12 ºC y 18 ºC, mientras que la temperatura del mes más frío llega a ser de -3 ºC y la del mes más cálido de 22 ºC, con una precipitación media anual de 600 mm (García, 2004). Los análisis de laboratorio realizados indicaron que el tipo de suelo predominante es el Phaeozem con una textura franco-arcillo-arenosa, pH de 7.3 y un contenido de materia orgánica de 5.04%. El tipo de vegetación corresponde a un bosque de pino dominado por P. cembroides y asociados a este, Juniperus flaccida Schltdl., Yucca queretaroensis Piña, Agave lechuguilla Torr. y Opuntia spp.
Obtención de muestras y medición de la biomasa
El método utilizado para la medición de biomasa aérea fue similar al propuesto por Picard, Saint-Andre y Henry (2012), que consiste en cortar a cada árbol a ras de suelo y separarlo en componentes estructurales. Un total de 21 individuos fueron derribados de distintas categorías diamétricas (CD) (5 cm a 40 cm). La distribución de los árboles fue de dos individuos para las CD de 5 cm, 10 cm y 30 cm y tres individuos para la CD de 15 cm; mientras que en las CD de 20 cm y 25 cm fue de cinco individuos y un individuo para las CD de 35 cm y 40 cm. Los árboles seleccionados para su derribo estaban libres de plagas y enfermedades y no mostraban algún daño físico o mecánico. Previo al derribo, el diámetro normal (Dn) fue medido con cinta diamétrica a una altura de 1.30 m desde el nivel del suelo, mientras que la altura total (H) fue medida con cinta métrica posterior al derribo como lo sugieren Acosta-Mireles et al. (2002) y Carrillo-Anzures, Acosta-Mireles, Jiménez-Cruz, González-Molina y Etchevers-Barra, (2016). Los componentes estructurales en que se dividió cada árbol fueron: fuste, ramas y ramillas con hojas (diámetro < 1 cm).
El fuste con corteza fue seccionado en trozas de 1.0 m de longitud desde la base hasta la punta o hasta la primera bifurcación del árbol, se obtuvieron muestras de ramas y muestras de ramillas con hojas, para ello la copa se dividió en tres secciones transversales: baja, media y alta (Zhang, Wang, Wang y Quan, 2009). De acuerdo con el procedimiento sugerido por Acosta-Mireles et al. (2002), Aquino-Ramírez et al. (2015) y Carrillo-Anzures et al. (2016) todos los componentes estructurales se pesaron por separado directamente en campo para obtener los pesos frescos totales (kg) con una báscula digital de 200 kg (TECNOCOR modelo IPEN) y una precisión de ± 0.5 kg. Después, una rodaja (3 cm - 5 cm de grosor) fue obtenida en cada troza del fuste y tres rodajas al azar para las ramas en cada sección de la copa. De igual manera, para las ramillas con hojas se obtuvieron tres muestras al azar de 0.5 kg por sección. Las rodajas y muestras se pesaron inmediatamente con una báscula electrónica de 10 kg (TECNOCOR modelo PEN-20) con una precisión de ± 5 g y puestas en bolsas de papel debidamente identificadas; luego fueron transportadas al laboratorio y colocadas por 96 h en una estufa de secado de madera (Nyle modelo L500) a 75 °C o hasta lograr un peso constante. Finalmente, el peso seco de cada componente del árbol fue calculado extrapolando la relación peso seco a peso fresco de las muestras. La biomasa aérea total (Bt) de cada árbol correspondió a la suma total de los pesos secos (kg) de cada componente estructural (fuste, ramas y ramillas con hojas).
Ecuaciones de biomasa
Los siguientes dos modelos alométricos ampliamente utilizados, fueron seleccionados para estimar la biomasa de los componentes de P. cembroides (Ec. 1 y Ec. 2). Estos modelos relacionan a la biomasa como una función del Dn y la H o la combinación de ambas variables (Dn 2 H) (Bi, Turner y Lambert, 2004; Zhao, Kane, Markewitz, Teskey y Clutter, 2015):
donde:
y = |
biomasa de cualquiera de los componentes estructurales del árbol (kg) |
Dn = |
diámetro normal (cm) |
H = |
altura total (m) |
β 0 , β 1 y β 2 = |
parámetros de regresión |
ε = |
error aleatorio |
El primer sistema de ecuaciones de aditividad (S1) fue basado en la Ecuación 1 y formulado de la siguiente forma:
Como consecuencia de la aditividad del sistema, la estimación de la Bt del S1 será el resultado de la sumatoria de la estimación de la biomasa de los tres componentes, expresada de la siguiente manera:
Mientras que el segundo sistema de ecuaciones de aditividad (S2) se fundamentó en la Ecuación 2 y se especificó de la siguiente forma:
Y la aditividad de la Bt fue expresado como:
donde:
β f , = |
biomasa en fuste (kg) |
β r , = |
biomasa en ramas (kg) |
β rh = |
biomasa en ramillas con hojas (kg) |
β t , = |
biomasa total (kg) en kg respectivamente; |
β ij = |
coeficientes de regresión |
ε 1 - ε 4 = |
términos del error |
En los dos sistemas de ecuaciones los errores son aditivos y están inherentemente correlacionados, los cuales son expresados como un vector
con un valor esperado
y una matriz de varianza y covarianza
(Bi et al., 2004; Wang, Bi, Ximenes, Ramos y Li, 2017). La propiedad de aditividad se obtuvo en los dos sistemas de ecuaciones al expresar a la biomasa aérea total (Bt) como la suma de la biomasa de cada uno de los componentes estructurales del árbol (B f , B r , B rh ) (Parresol, 2001).
Ajuste y comparación
Los dos sistemas propuestos fueron ajustados simultáneamente usando el método generalizado de momentos (GMM, por sus siglas en inglés) con el procedimiento PROC MODEL del programa SAS/ETS® (SAS Institute Inc., 2004). Con este método se estiman parámetros en ecuaciones no lineales con términos de error aditivo más eficientes en condiciones heterocedásticas sin ninguna especificación de su naturaleza (Greene, 1999), especialmente en muestras pequeñas, donde el principal problema de estimación es la obtención de funciones de ponderación de los componentes de biomasa que describan con precisión el diseño de la matriz de varianza del error (Bi et al., 2004; Bi et al., 2015). Para evaluar y comparar las ecuaciones desarrolladas en cada sistema se obtuvieron cuatro estadísticos: el coeficiente de determinación ajustado (R
2
adj
), la raíz de los cuadrados medios del error (RMSE), el sesgo promedio (
Resultados
Variación en la distribución de biomasa aérea
Los árboles derribados de P. cembroides en el presente estudio registraron un intervalo de diámetros de 3 cm a 40 cm con una media de 20.2 cm ± 8.6 cm y alturas de 1.2 m a 9.4 m con un valor medio de 6.9 m ± 1.9 m; mientras que la biomasa total fue de 2.8 kg por árbol a 700.1 kg por árbol con un valor promedio de 180.8 kg ± 163.8 kg. Por otra parte, la distribución de biomasa entre los componentes estructurales (fuste, ramas y ramillas con hojas) mostraron importante variación debido a los tamaños del árbol (Fig. 1). El fuste y las ramas fueron los principales componentes del árbol donde se acumuló la mayor cantidad de biomasa y la de menor cantidad fueron las ramillas con hojas. La proporción en cuanto a la cantidad de biomasa del fuste manifestó una amplitud de valores entre 26.1% y 68.7% con un promedio de 43.2%. En cambio, la proporción de la cantidad de biomasa en ramas fue mayor que en el fuste, siendo de 18.8% en las categorías diamétricas de 10 cm y de 64.6% en las categorías diamétricas de 40 cm, con una media de 41.2%; algo similar sucedió con la proporción de la biomasa en las ramillas con hojas y se presentó una variación de 6.3% a 32.1%, en las categorías diamétricas de 5 cm a 10 cm con un valor medio de 15.5%.
Evaluación de las ecuaciones de biomasa
El valor estimado, error estándar y la significancia de cada parámetro para los dos sistemas de ecuaciones de aditividad (S1 y S2) son presentados en la Tabla 1. Todos los parámetros fueron significativamente diferentes de cero (p < 0.001). Los coeficientes de R
2
adj
obtenidos en las ecuaciones de distribución de biomasa se encontraron en un intervalo de 0.85% - 0.97% y para las ecuaciones de biomasa aérea total fueron superiores a 0.97% (Tabla 1). La ecuación para predecir la biomasa de ramillas con hojas, que incluyó al Dn
2
H como variable predictora fue la que mostró el valor más bajo de R
2
adj
(0.85) y valores altos de
Tabla 1 Parámetros estimados y estadísticos de bondad de ajuste de las ecuaciones de aditividad para estimar la biomasa de los componentes estructurales y la biomasa aérea total de Pinus cembroides.
| Sistema de ecuaciones | Componente del árbol | Parámetro | Estimador | Error estándar | Valor de t | Pr > |t| | RMSE (kg) |
|
|
|
| S1 | Fuste |
|
2.766848 | 0.2151 | 12.86 | <0.0001 | 9.4955 | 0.5266 | 1.72 | 0.9641 |
|
|
0.848754 | 0.0233 | 36.43 | <0.0001 | ||||||
| Ramas |
|
4.182507 | 0.2549 | 16.41 | <0.0001 | 23.8367 | -0.0715 | 30.27 | 0.9408 | |
|
|
1.0513449 | 0.0276 | 38.05 | <0.0001 | ||||||
| Ramillas con hojas |
|
2.838132 | 0.3376 | 8.41 | <0.0001 | 7.1995 | -0.4596 | 21.41 | 0.8597 | |
|
|
0.738552 | 0.0357 | 20.71 | <0.0001 | ||||||
| Total | 25.3352 | -0.0035 | 5.00 | 0.9761 | ||||||
| S2 | Fuste |
|
-2.007449 | 0.3337 | 5.68 | <0.0001 | 8.6137 | -0.0001 | 4.75 | 0.9704 |
|
|
1.776093 | 0.1028 | 17.28 | <0.0001 | ||||||
|
|
0.364976 | 0.2675 | 1.36 | <0.0001 | ||||||
| Ramas |
|
-5.193296 | 0.4592 | 11.31 | <0.0001 | 22.4340 | 0.0103 | 16.58 | 0.9476 | |
|
|
2.36942 | 0.1743 | 13.59 | <0.0001 | ||||||
|
|
1.11539 | 0.4266 | 2.61 | 0.0176 | ||||||
| Ramillas con hojas |
|
2.07281 | 0.4878 | 4.25 | 0.0005 | 5.1074 | 0.0005 | 19.12 | 0.9294 | |
|
|
2.463058 | 0.1935 | 12.73 | <0.0001 | ||||||
|
|
-1.17492 | 0.4276 | -2.75 | 0.0132 | ||||||
| Total | 23.0285 | 0.0116 | 5.76 | 0.9802 |
S1=sistema de ecuaciones de aditividad basado en el modelo
Los resultados de bondad de ajuste indicaron los mejores estadísticos en las ecuaciones del S2 que las obtenidas del S1. Las ecuaciones del S2 explicaron la mayor variabilidad de los datos observados con coeficientes de R
2
adj
superiores a 0.92, así mismo, valores bajos de RMSE (< 25 kg),
Figura 2 Comparación entre los valores observados y las curvas de predicción desarrolladas a partir de las ecuaciones alométricas de aditividad para a) la biomasa de fuste (Bf), b) biomasa de ramas (Br), c) biomasa de ramillas con hojas (Brh) y d) biomasa total (Bt) de Pinus cembroides.
Discusión
Los dos conjuntos de ecuaciones de aditividad obtenidos en este estudio proporcionan una herramienta operacional para estimar con precisión la biomasa aérea de los árboles individuales de P. cembroides. La variación explicada por las ecuaciones de aditividad fue alta, esto se debe a que en su estructura se incluyeron como variables independientes al Dn y a la H. El uso de estas variables en la predicción de la biomasa ha sido discutido ampliamente por diferentes autores, debido a que pueden mejorar significativamente los resultados de ajuste y la capacidad de predicción (Bi et al., 2004; Li y Zhao, 2013), en lugar de emplear únicamente al Dn como variable predictora. Además, al incluir a la H en las ecuaciones se tiene información relacionada sobre la competencia y/o fertilidad del sitio ocasionando estimaciones menos sesgadas (Tesfaye et al., 2015).
El uso del Dn como única variable predictora para la estimación de biomasa a nivel del árbol, ha mostrado buenos estadísticos de ajuste para diferentes especies de coníferas (Návar, González, Graciano, Dale y Parresol, 2004; Acosta-Mireles et al; 2002; Díaz-Franco et al., 2007). Sin embargo, la inclusión de la H como segunda variable predictora tiene la ventaja de ser un factor predictivo que debe tomarse en cuenta en la cuantificación de la biomasa, debido a que se espera que en árboles con el mismo Dn, pero distinto ahusamiento, tengan alometrías diferentes (Nogueira, Fearnside, Nelson, Barbosa y Keizer, 2008; Temesgen et al., 2015). Además, la altura es importante para explicar las diferencias específicas en la arquitectura o la plasticidad fenotípica de los árboles dependiendo de las condiciones de crecimiento (Archibald y Bond, 2003). Una de las principales razones para omitir la altura del árbol es su dificultad para medirla en condiciones de dosel cerrado (Díaz-Franco et al., 2007; Hunter, Keller, Victoria y Morton, 2013; Larjavaara y Muller-Landau, 2013). No obstante, existen diferentes métodos que garantizan su correcta medición (Sullivan et al., 2018).
El sistema S2 basado en el modelo
fue el que obtuvo mejores estadísticos de bondad de ajuste, expresando la relación entre la biomasa de los componentes con las variables dasométricas del árbol, por lo tanto, fue seleccionado para estimar la biomasa de los componentes específicos y del árbol individual de P. cembroides. Este modelo ha sido considerado como uno de los más flexibles ya que mejora la precisión de las estimaciones de la Bt y de los componentes de la copa, especialmente la biomasa de ramas y de las hojas (Zhou, Brandle, Schoeneberger y Awada, 2007; Dong, Zhang y Li, 2015).
Las ecuaciones de aditividad del S2 explican de 92% a 97% de la variabilidad de la biomasa de los componentes y 98% de la variabilidad de la biomasa aérea total, valores que son diferentes a los señalados por Návar (2011) en modelos de potencia para estimar la partición de biomasa en P. cembroides, Pinus pinceana Gordon, Pinus pseudostrobus Brongn, Pinus brutia Ten y Cupressus arizonica Greene, incluyendo únicamente al Dn como variable predictora en especies de bosques templados secos del noreste de México, con coeficientes de determinación de 0.87, 0.60 y 0.81 para la biomasa del fuste, biomasa de hojas y ramas, y en la biomasa total, respectivamente. Estas diferencias pueden atribuirse a la alta variabilidad que se registra al integrar en un solo modelo de distribución de biomasa a cinco especies de coníferas, representadas por un número reducido de árboles por especie. Por el contrario, los resultados en este estudio son semejantes a los hallados por Chávez-Pascual, Rodríguez-Ortiz, Enríquez-del Valle, Velasco-Velasco y Gómez-Cárdenas (2017) en Pinus oaxacana Mirov bajo manejo en Ixtlán Oaxaca, México, quienes ajustaron modelos de la variable combinada generalizada para estimar la biomasa del fuste con corteza, biomasa de ramas y biomasa de hojas con coeficientes de determinación de 0.97, 0.93 y 0.98, respectivamente; y un valor de 0.97 para el R 2 adj en la ecuación de la biomasa total. Así mismo, a los obtenidos por Vargas-Larreta et al. (2017) quienes utilizaron ecuaciones de aditividad con la misma estructura a los de este estudio, con coeficientes de determinación de 0.85 a 0.95 para la biomasa área total y coeficientes de 0.63 a 0.97 en la distribución de biomasa en 11 especies de coníferas del noroeste de México.
La aditividad o la compatibilidad es una característica deseable al desarrollar un sistema de ecuaciones de biomasa para eliminar la inconsistencia entre la suma de los valores predichos para los diferentes componentes (fuste, ramas y follaje) y la predicción total del árbol (Bi et al., 2004; Dong et al., 2015), propiedad que frecuentemente es ignorada. En México son limitados los estudios que toman en cuenta estas relaciones de aditividad, por ejemplo, Návar et al. (2004), desarrollaron ecuaciones aditivas ponderadas para el cálculo de la biomasa en plantaciones de Pinus durangensis Martínez, Pinus cooperi Blanco y Pinus engelmannii Carr. de la Sierra Madre Occidental en Durango, México, de igual manera, Monroy-Rivera y Návar-Cháidez (2004) generaron ecuaciones de aditividad para componentes de biomasa en plantaciones comerciales de Hevea brasiliensis en Veracruz, México y Aquino-Ramírez et al. (2015) propusieron un sistema de ecuaciones de partición de biomasa con aditividad para tres especies de árboles tropicales en Pluma Hidalgo Oaxaca, México. Recientemente, Vargas-Larreta et al. (2017) ajustaron simultáneamente sistemas de ecuaciones para la biomasa de los componentes y la biomasa total para asegurar la aditividad, en 17 especies forestales en los bosques templados del noroeste de México, así mismo, Flores-Medina et al. (2018) desarrollaron ecuaciones alométricas considerando la aditividad para estimar biomasa por fracción de grosor y total en cuatro especies en etapa de regeneración del suroeste de Durango, México.
La proporción de biomasa en las ramas se incrementó conforme aumentaron los tamaños del árbol. En las categorías diamétricas de 25 cm a 40 cm la asignación de biomasa superó más de 40% a la biomasa del fuste (Fig. 1). Estos resultados son consistentes con el estudio realizado por Grier, Elliott y McCullough, (1992) en un bosque mixto de Piñón-Juniperus en el norte de Arizona, encontraron que la asignación de biomasa en las ramas de Pinus edulis Engelm fue de 43.7% en rodales jóvenes con Dn promedio de 16 cm y se incrementó a 51.1% en rodales maduros con Dn promedio de 25 cm. Los patrones de asignación de biomasa en los pinos piñoneros son diferentes a la mayoría de las especies del género Pinus donde la mayor cantidad de biomasa se concentra principalmente en el fuste (Chávez-Aguilar et al., 2016; Vargas-Larreta et al., 2017). Lo anterior se atribuye a que los árboles que crecen en espacios abiertos tienden a presentar copas extendidas (Návar, 2009), porque las ramas no tienen una buena poda natural lo que favorece un crecimiento libre y, en consecuencia, el árbol destina una mayor proporción de biomasa a los tejidos relacionados con el soporte del follaje y a las estructuras reproductivas (Loveall y Harrington, 2008).
Las ecuaciones de biomasa con aditividad presentadas en este estudio aportarán información más eficiente para mejorar la cuantificación de la biomasa y carbono, con el propósito de ayudar en el manejo sustentable de los bosques piñoneros de P. cembroides. Se recomienda la aplicación de las ecuaciones dentro del intervalo de diámetros en el que fueron construidas porque las relaciones alométricas no son constantes y pueden cambiar de diferentes maneras según la especie, el tamaño, calidad del sitio y la edad (Riofrío et al., 2015), lo que ocasionaría resultados con menor precisión.
Conclusiones
El sistema de ecuaciones alométricas (S2) propuesto que integró al diámetro normal y la altura total del árbol como variables predictoras permite estimaciones robustas de la biomasa de los componentes estructurales y la biomasa aérea total de Pinus cembroides, asegurando la aditividad de los componentes.
En los árboles de P. cembroides la proporción de biomasa por componentes fue de 43.2%, 41.2% y 15.5% para fuste, ramas y ramillas con hojas, sin embargo, en las categorías diamétricas de 25 cm a 40 cm la acumulación de biomasa fue mayor en ramas que la del fuste con más de 40%.
Las ecuaciones desarrolladas en este estudio pueden aplicarse de manera sencilla y precisa a datos de árboles individuales e inventarios forestales para estimaciones de biomasa y carbono, lo que contribuirá al desarrollo de prácticas de manejo sustentable y estrategias de conservación de la especie, donde la producción de madera no es la función principal.
