Recibido, Febrero 2018; Aceptado, Julio 2018 Articulo original

Estimación del potencial de captura de carbono en dos especies del género

Polylepis, Recuay - Ancash.
Estimation of carbon capture potential in two species of the genus Polylepis, Recuay - Ancash.
1 a 1 b
Eleuterio Ramírez Apolinario Peláez Peláez Freddy

RESUMEN
El objetivo de esta investigación fue estimar el potencial de captura de carbono en especies de
Polylepis incana y Polylepis serícea, que crecen en bosques nativos del caserío de Poccrac, distrito de
Ticapampa, Ancash. Las ecuaciones alométricas para ambas especies, relacionando la biomasa (M)
como variable dependiente con el diámetro de altura de pecho como variable independiente (DAP). La
medida del DAP de los árboles de Polylepis talados, están dentro del rango que va desde 5 cm a 20cm.
En el análisis de regresión realizado con el programa estadístico SPSS V 22.0 ,encontramos para la
especie Polylepis incana un coeficiente de determinación de r2=0.929 y para la especie Polylepis
2
sericea un r =0.94, lo cual indica un alta relación entre las variables biomasa y DAP, mostrando buena
bondad de ajuste para los parámetros que constituyen la ecuación alométrica de ambas especies. Los
resultados obtenidos constituyen un importante aporte metodológico para las estimaciones de captura
de carbono en zonas destinadas a la conservación y una herramienta para valorar la función de captura
de CO2 en bosques andinos.

Palabras clave: biomasa, ecosistemas andinos, ecuaciones alométricas.

ABSTRACT
The objective of this research was to estimate the potential for carbon capture in Polylepis incana and
Polylepis sericea species, which grow in native forests of the village of Poccrac, district of Ticapampa,
Ancash. He allometric equations for both species, relating the biomass (M) as a dependent variable with
the breast height diameter as an independent variable (DAP). The DAP measurement of the cut
Polylepis trees is within the range from 5 cm to 20 cm. In the regression analysis performed with the
statistical program SPSS V 22.0, we find for the Polylepis incana species a coefficient of determination
of r2 = 0.929 and for the Polylepis sericea species a r2 = 0.94, which indicates a high relation between the
variables biomass and DAP, showing good goodness of fit for the parameters that constitute the
allometric equation of both species. The results obtained constitute an important methodological
contribution for the estimates of carbon capture in areas destined for conservation and a tool to assess
the CO2 capture function in Andean forests.

Keywords: biomass, Andean ecosystems, allometric equations.

1
Universidad Nacional de Trujillo – Escuela de Postgrado. Trujillo - Perú.
a
Físico
b
Biólogo

Rev. Tayacaja 1 (1): 17-27, enero - julio 2018 17

Estimación del potencial de captura de carbono Ramírez E.

INTRODUCCIÓN troncos, flores y frutos (Rodríguez et al., 2006)..

Siendo los arboles los que retienen en mayor
La actividad antrópica en la zona rural, el mal cantidad este elemento por unidad de área en
manejo y uso de los recursos naturales, las altas comparación con otros tipos de vegetación
emisiones del dióxido de carbono, vienen (IPCC, 2007).
ocasionando cambios en el planeta que son En el Perú, los bosques andinos son remanentes
evidentes; el incremento de la temperatura del de grandes extensiones de bosques de “quenual”
planeta, las concentraciones de los gases de (Polylepis spp.) y “colle” (Buddleia spp.), así
efecto invernadero, las pérdidas de masas de también existen 19 especies del género Polylepis
nieve y hielo, entre otros, constituyen efectos de las 28 que se tienen registradas y reportadas a
palpables . A lo largo de la extensión de nuestro nivel mundial, se habla de un centro de
territorio nacional existen cadenas de montañas diversidad del género ya que concentra gran
que sobrepasan altitudes 3000 msnm, que cantidad de especies. Un Decreto Supremo
permiten el desarrollo de bosques en dictado en el año 2006, cuyo objetivo es
ecosistemas diversos y paisajísticamente identificar y priorizar elementos de la Diversidad
extraordinarios los cuales son social y Biológica que se encuentran en amenaza por lo
económicamente importantes . La mayoría de que deben ser conservados y requieren atención
los ecosistemas podrían verse afectados por prioritaria, brinda la Lista Roja de Especies de
estas alteraciones que experimenta nuestro Flora y Fauna Amenazadas elaborada por la
clima.( Bernier y Schoene, 2007) Unión Mundial para la Conservación, del cual se
La vegetación de los Andes Centrales está extrae que en la categoría de Peligro Crítico se
dominada por zonas agrícolas, pastizales y encuentran las especies P. incana y P. racemosa,
zonas arbustivas. Los bosques de Polylepis en Peligro comprenden P. microphylla, P.
representan la vegetación natural de una gran multijuga, P. subsericans y P. tomentella, en la
parte de los Andes, a altitudes entre 3 500 m. a 5 categoría Vulnerable están las especies P.
200 m. las especies del género ocupan una gran besseri, P. pepei , P.rugulosa, P. serícea, P.
variedad de hábitats, desde el límite superior de tarapacana y P. weberbaueri, por último en la
los bosques de neblina hasta los volcanes áridos categoría Casi Amenazado se encuentra P. pauta.
del Altiplano. Sin embargo, durante milenios las Es decir 12 de las 19 especies presentes en el país
actividades y las necesidades humanas en los requieren planes de conservación debido a las
Andes han destruido estos bosques, amenazas que atraviesan, sin embargo, el vacío
restringiéndolos a hábitats especiales y de información existente sobre estas especies
modificando su composición florística y dificulta su conservación (Llerena et al. 2014;
faunística. Mendoza y Cano, 2011).
El CO2 presente en la atmosfera se captura en los En las Región Ancash, especialmente en la
procesos metabólicos de las plantas a través de la Cordillera Blanca, la vegetación mayormente se
fotosíntesis; en este proceso, el carbono presenta en forma de pajonales y pequeños
atmosférico secuestrado se expresa en términos arbustos, siendo el género Polylepis, el elemento
de biomasa constituido por follaje, ramas, raíces, arbóreo dominante.

18 Rev. Tayacaja 1 (1): 17-27, enero - julio 2018

Estimación del potencial de captura de carbono Ramírez E.

Su tala, (a pesar de ser una especie protegida) Avendaño et al., (2009), sostiene que: “la
está generando pérdida de la diversidad estimación de la biomasa en zonas que están
biológica, alteración del ciclo hidrológico y dedicadas a la conservación requieren de la
cambio micro meteorológico, afectando a la implementación de un método de medición con
población local y a la vez acelerando la datos que se ajusten a la estructura particular de
degradación del ecosistema. Su conservación cada bosque sin la necesidad de generar tala de
representa una prioridad por su extraordinaria árboles”. Ziani y Mencuccini, (2004); proponen
riqueza y endemismo, y porque varias de sus métodos indirectos a partir de mediciones
especies constituyentes están siendo dasométricas precisas como por ejemplo el
severamente amenazados (Aubad et al., 2010) diámetro de altura de pecho( DAP), como una
Frente a esta problemática de debe buscar alternativa adecuada para la estimación de la
mecanismos de solución diferentes a los biomasa en áreas donde existen restricciones
utilizados en la actualidad, que pasa solo por para la tumba de árboles y se requieren resultados
prohibir la tala de esta especie, actividad que no con bajo margen de error.
puede ser controlado en su totalidad debido a la Entonces para esta investigación se ha utilizado
extensión del área geográfica por donde se modelos alométricos como metodología que
extienden los pequeños bosques de Polylepis y permite obtener una estimación confiable y
por la necesidad específica de cada población directa de la biomasa y por ende del carbono en
que se ubica en las cercanías de estos bosques. los sistemas vegetales, porque tal estimación es
Una alternativa viable según Ruignitz et al, individual y mediante la sumatoria se puede
(2009), sería realizar intervenciones productivas obtener una estimación global para cada sistema.
de servicios ambientales, que tengan como meta (Etchevers et.al. 2002).
la estimación de captura de carbono en bosques
de productores familiares o comunidades MATERIALES Y MÉTODOS
rurales, como una forma de generar beneficios
financieros mediante el acceso al mercado de Objeto de estudio
bonos de carbono. Restaurando las áreas Biomasa de plantas de Polylepis spp de un
degradadas por medio de plantaciones y/o bosque interandino ubicado en la zona de
regeneración natural, y por la extracción de Poccrac, Disrito de Ticapampa en Ancash.
madera. En ambos casos se pretende almacenar Especies: Polylepis incana y Poyilepis sericea.
el carbono a través del crecimiento de árboles y, Selección de la muestra
al extraer la madera, convertirla en productos La selección de las plantas a talar se realizó
durables. El carbono acumulado se mantendrá empleando un muestreo no probabilístico de tipo
durante la vida útil del producto. Al extraer la intencional, es decir se seleccionó el tipo de
madera, la regeneración actuará almacenando planta de acuerdo a sus características
carbono por el crecimiento. Los sistemas dasometricas como: diámetro de altura de pecho
forestales y agroforestales pueden capturar en (DAP), altura de copa, fuste recto y de fácil
sus diferentes almacenes de 80 a 350 toneladas acceso. Las plantas taladas estuvieron en un
de carbono por hectárea. (Camilo y Díaz 2010). rango de: 5cm ≤ DAP ≤ 20cm.

Rev. Tayacaja 1 (1): 17-27, enero - julio 2018 19

Estimación del potencial de captura de carbono Ramírez E.

Para ambas especies se realizó 11 talas. Éstas vs DAP). Mediante análisis de regresión se forma
fueron hechas en el periodo que va de Mayo a una ecuación que se denomina alométrica, que
Diciembre que es el periodo de carencia de puede ser extrapolada a toda la población en
lluvias en la zona sierra de Ancash. estudio.
Área de estudio Ecuación alométrica
La recolección de datos se realizó en remanentes La relación M vs DAP, tiene comportamiento
de un bosque natural, ubicado en las quebradas polinomial de la forma:
A
del Caserío de Poccrac, Distrito de Ticapampa- M= B (DAP) ............................................. 1
R e c u a y. D e p a r t a m e n t o d e A n c a s h .
Geográficamente situado entre las coordenadas: Dónde: M= Biomasa aérea; DAP= Diámetro de
Latitud 09°45l28ll Sur.Longitud:77°26l 40ll altura de pecho (1.30 m de altitud de la planta
Oeste. Y una altitud de 3465msnm (Código respecto al suelo); A y B son parámetros para
Ubigeo: 021710-RENIEC). estimar
La expresión (1) se linealiza aplicando
Métodos y técnicas logaritmo natural a ambos lados de la ecuación:
En este trabajo se utilizó una metodología basada LNM = LNB + A LN (DAP)............................2
en la correlación entre el DAP y la biomasa, ya
que según Solano et al. (2014), ésta constituye la Es decir la expresión (2), coincide con una
variable que frecuentemente más se relaciona ecuación lineal de la forma:
con la biomasa y es la más fácil de ser medida en Y= a + b X
campo, ya que la altura y otras variables son Donde: Y = LNM ; X = LN (DAP); b = LNB
menos precisas de identificar y medir. Además, (intercepto de línea recta); a=A (pendiente de
según el mismo autor, los modelos alometricos línea recta); Las contantes: a, y b. (Son
con variables logarítmicas y con exponentes parámetros para estimar)
cuadrados, son los que mejor grado de
significancia tienen; es decir son los modelos Estimación de los parámetros de la ecuación
más confiables y que mejor se ajustan para alométrica
calcular la biomasa aérea de las especies
forestales. Método de mínimos cuadrados
Para estimar los valores de los parámetros a y b de
Estimación de biomasa la recta que mejor se ajusta a los datos
Para estimar la biomasa se utilizó el método experimentales, se utilizó el método de mínimos
indirecto, también llamado no destructivo, que cuadrados. El valor de éstos parámetros se
consiste en cortar y pesar una cantidad definida determina de la siguiente manera:
de árboles, el peso (M) de estas plantas se
relaciona con una variable dasométrica de fácil
nå xy - å xå y å y - aå x
medición, que en nuestro caso es el diámetro de a= b =
altura de pecho (DAP), luego se observa el nå x 2 - ( å x ) 2 n
comportamiento de la relación en un gráfico (M

20 Rev. Tayacaja 1 (1): 17- 27, enero - julio 2018

Estimación del potencial de captura de carbono Ramírez E.

Bondad de ajuste Técnica de recolección y pesaje de muestras

El criterio más sencillo para determinar la En la toma de datos se procedió a medir el DAP de
bondad de ajuste para una recta de regresión lo da la planta antes de talar, luego de la tala se
el coeficiente de correlación de Pearson: separaron las partes: fuste, ramas y hojas y se
pesaron separados. Luego del talado y el pesado
r=
å ( x - x)( y - y) en campo, se procedió a trasladar las plantas a un
2 2 almacén para su secado respectivo.
å ( x - x) å ( y - y )
El secado se realizó de manera natural y en
sombra, para ello se monitoreo cada semana el
Coeficiente de determinación: avance de secado hasta que alcanzó un peso
El coeficiente de determinación expresa la constante (que aproximadamente sucede a las 8
cantidad de variación de la variable dependiente semanas), a partir de este valor se considera el
que es explicada por la variable independiente. peso en seco de la planta. Datos consignados en
Cuanto más próximo a uno es el valor de r2 la tabla 1, para la especie Polylepis Incana y la
mejor fue realizado el ajuste. Machado et tabla 2 para la especie Polylepis Sericea.
al.(2002).
2
é ù
r =ê
2 å ( x - x )( y - y ) ú
ê 2 2 ú
ëê ( x - x) å ( y - y ) ûú

RESULTADOS
Tabla1. Datos del diámetro de altura de Pecho (DAP) y de la biomasa aérea de Polylepis Incana

N° 1 2 3 4 5 6 7 8 9 10 11

DAP (cm) 11,2 10,5 16,0 17,0 18,0 9,1 6,5 7,9 12,0 13,5 6,0

M (Kg) 16,5 12,9 25,1 28,5 35,5 14,0 5,62 6,4 18,9 28,5 3,2

Tabla 2. Datos del diámetro de altura de pecho (DAP) y de la biomasa aérea de Polylepis Sericea

N° 1 2 3 4 5 6 7 8 9 10 11
DAP (cm) 17,5 12,0 13,0 10,5 18,0 9,3 6,3 8,1 12,6 14,0 5,8

M (Kg) 29,0 17,0 24,0 13,1 36,0 14,2 5,5 6,6 20,4 28,0 3,0

Rev. Tayacaja 1 (1): 17-27, enero - julio 2018 22

Estimación del potencial de captura de carbono Ramírez E.

Para elaborar los modelos alometricos para ambas especies, se procedió a organizar un grafico con la
información de biomasa (M) y el diámetro de altura de pecho (DAP). En ella se observa un
comportamiento de tipo polinomial de la forma M= B (DAP) A , tal y como se observa en la Figura 1
para Polylepis incana y figura 2, para Polylepis serícea.

Figura 1. Comportamiento polinomial de la relación biomasa aérea y DAP, para Polylepis incana

Figura 2. Comportamiento polinomial de la relación biomasa aérea y DAP, para Polylepis Sericea

La ecuación polinomial fue transformado a una forma lineal aplicando logaritmo natural a ambos lados
de la ecuación, luego los parámetros A y B, fueron estimados con el método de mínimos cuadrados,
resultando valores ubicados en la tabla 3 y en la tabla 4.

23 Rev. Tayacaja 1 (1): 17-27, enero - julio 2018

Estimación del potencial de captura de carbono Ramírez E.

Tabla 3. Datos obtenidos al linealizar la función polinomial, aplicando logaritmo natural a las variables,
diámetro de altura de pecho (DAP) y biomasa aérea de la especie de Polylepis Incana.

N° 1 2 3 4 5 6 7 8 9 10 11

LNDAP (cm) 2,416 2,351 2,773 2,833 2,890 2,208 1,872 2,067 2,485 2,603 1,792

LNM (Kg) 2,803 2,557 3,223 3,350 3,570 2,630 1,726 1,856 2,939 3,350 1,163

Tabla 4. Datos obtenidos al linealizar la función polinomial, aplicando logaritmo natural a las variables,
diámetro de altura de pecho (DAP) y biomasa aérea de la especie de Polylepis Sericea

N° 1 2 3 4 5 6 7 8 9 10 11

LNDAP (cm) 2,862 2,485 2,565 2,351 2,890 2,230 1,841 2,092 2,534 2,640 1,758

LNM (Kg) 3,367 2,833 3,178 2,573 3,584 2,653 1,716 1,887 3,015 3,332 1,099

Al relacionar el Ln de la biomasa y el Ln del DAP , en ambas especies, resulta una relación lineal
positiva entre la biomasa y la variable DAP, se puede observar que el modelo funciona bien para las
variables consideradas.

Figura 3. Comportamiento lineal de la relación biomasa aérea y DAP, para Polylepis Incana .

Rev. Tayacaja 1 (1): 17-27, enero - julio 2018 24

Estimación del potencial de captura de carbono Ramírez E.

Figura 4. Comportamiento lineal de la relación biomasa aérea y DAP, para Polylepis Sericea .
Tabla 5. Resultado de análisis de regresión lineal con el programa estadístico SPSS V 22.0
Coeficientes

Nombre Constante LNDAP N r r2 r2 corregido Error tipo residual

Polylepis incana -2,137 2,0078 11 0,968 0,936 0,929 0,043

Polylepis sericea -2.1757 2,0250 11 0,970 0,941 0,934 1 0,04

DISCUSIÓN El valor alto del coeficiente de determinación

(r2) indica una gran variabilidad en los datos y
Para el análisis de la regresión lineal de las
que el modelo para ambas especies es muy
ecuaciones alométricas se utilizó el programa
adecuado y significativo. El estudio está
estadístico SPSS V 22.0; la cual arrojó como
restringido para diámetros de fuste en un
resultado valores de r2 = 0.929 para Polylepis
intervalo de, 5cm ≤ DAP≤ 20cm, por lo tanto la
incana y r2=0.934 para Polylepis sericea; esto
validez de las ecuaciones se dan para este rango
indica que para ambas especies existe una fuerte
de medición, para la misma especie y que
correlación entre las variables estimadas. El
presenten características similares de
93.6% de la variación de la variable biomasa de
crecimiento en la zona.
Polylepis incana se debe a la presencia de la
variable DAP en el modelo . El 6.4% se debe a El DAP como variable independiente es un buen
otros factores intervinientes. Para el caso de predictor de la biomasa aérea de las especies
Polylepis sericea , el 94.1% de la variación de la estudiadas, ya que las ecuaciones alométricas
variable biomasa se debe a la presencia de la obtenidas dieron r2 cercanos a la unidad. Aunque
variable DAP en el modelo. El 5.9% se debe a varios estudios dicen encontrar una buena
otros factores intervinientes. relación de la biomasa cuando se usan
simultáneamente las variables diámetro y altura,
25 Rev. Tayacaja 1 (1): 17-27, enero - julio 2018
Estimación del potencial de captura de carbono Ramírez E.

en la mayoría de los casos con solo emplear el relaciones alométricas en seis especies
DAP también se logra una buena estimación, tal arbóreas en oaxaca, mexico. agrociencia
y como lo afirman Brown (2001) y Acosta et al. 36. 752-736
(2002). Las constantes de la ecuación alométrica
Aubad, J., Aragón, P. y Rodríguez, 2010. Human
para Polylepis Incana reportado por la siguiente
access and landscape structure effects on
investigación, es comparable a lo encontrado
Andean forest bird richness, Acta
para la misma especie por Schlegel et al.
Oecologica 36: 396-402
(2001), quien establece que, los modelos que
mejor estiman biomasa son los de tipo alomerico, Brown, S. (2001). Estimating biomass and
de la forma Y=a Xb con transformación biomass change of tropical forest: a
logarítmica que incluye variable de DAP. La primer. Roma. FAO forestry paper, 134.
validez de las ecuaciones alometricas generadas
Bernier, P. y Schoene, D. 2007. La adaptación de
en este estudio debe restringirse al rango de
los bosques y su ordenación al cambio
variación de las dimensiones de DAP para un
climático.
intervalo de 5cm ≤ DAP≤20cm
Camilo, P.; y Diaz, T. 2010. Estimación del
CONCLUSIONES
carbono contenido en la biomasa forestal
La ecuación alométrica calculada para la aérea de dos bosques andinos en los
estimación de la biomasa aérea es altamente Departamentos de Santander y
confiable, ya que el valor del r2 para las dos Cundinamarca. Colombia. UDFJC.
ecuaciones es superior a 0.9, lo cual muestra alta
Etchevers.J; Vargas J; M. Acosta; A. Velásquez
dependencia entre las variables analizadas en los
.2002. Estimación de la biomasa, aérea
gráficos de correlación, lo que quiere decir que
mediante el uso de relaciones alometricas
puede ser usada para estimaciones de este
en seis especies arbóreas en Oaxaca-
componente restringidos solo para Polylepis
México Agrociencia. Vol.36.Nº006
Incana y Polylepis Sericea , debido a que la
ecuación fue desarrollada bajo las características IPCC (Intergovernmental Panel on Climate
propias de estas especies. La información Change). 2007. Climate Change (2014):
obtenida en el presente trabajo, podrá ser The Physical Science Basis. Summary for
utilizada como información de entrada para Policymakers. Paris, Francia. s.e. p 2.
modelización de escenarios de cambio climático
Llerena, P.;S.Yalle ;E. Silvestre. 2014. Los
y como linea base en futuros estudios de
bosques y el cambio climático en el Pero:
conservación de la especie .
Situación y perspectivas. Organización de
REFERENCIAS BIBLIOGRÁFICAS la Naciones Unidas para la alimentación y
la agricultura. Lima –Perú.
Avendaño, R., Galindo, A. And Angulo, A. 2009.
Ecología y educación ambiental . Machado, A., Conceicao, B., Figueiredo,J 2002.
Universidad autónoma de Sinaloa – Modelagem do volume individual para
Culiacan. Mexico DF. diferentes idades e regimes de desbaste
em plantacoes de pinus oocarpa . Revista
Acosta M. , Vargas H., Velasquez M. y
Ciencias Exactas e Naturais, Vol. n°4
Etchevers B. 2002. estimación de la
Jul/dez. : 185-197.
biomasa aérea mediante el uso de
Rev. Tayacaja 1 (1): 17-27, enero - julio 2018 26
Estimación del potencial de captura de carbono Ramírez E.

Mendoza, Wilfredo y Cano, Asunción. 2011. 311-332.

Diversidad del género Polylepis (
CORRESPONDENCIA:
Rosaceae , Sanguisorbeae ) en los Andes
peruanos. Revista Peruana de Biología. Mag. Eleuterio Ramirez Apolinario
Vol. 18, no. 2, p. 197–200
e-mail: [email protected]
Pearson, T. Walker, y Brown, S. 2005.
Sourcebook for land use, land-use
change and forestry projects. Winrock
International and the BioCarbon Fund
of the World Bank. 64 p.

Rodriguez, R. Jimenez, J. Aguirre, O. yTreviño,

E. 2006. Estimación del carbono
almacenado en un bosque de niebla en
Tamaulipas, México. CIENCIA UANL
9: 179 p.

Rugnitz,M. T.; Chacón, M,L; Porro R. 2009.

Guía para la determinación de carbono
en pequeñas propiedades rurales 1ª
edición, Lima Perù: Centro Mundial
Agroforestal (ICRAF)/ Consorcio
Iniciativa Amazónica (IA)

Schlegel, B., Gayoso, J. y Guerra, J. 2001 .

Medición de la capacidad de captura de
carbono en bosques de Chile y
promoción en el mercado mundial.
Manual de procedimientos para
inventarios .Universidad Austral de
Chile. Valdivia, Chile.

Solano,D.; C.vega, V. H. Heras, K.Cueva.

2014.Generación de modelos
alometricos para determinar biomasa
aérea a nivel de especies, mediante el
método destructivo de baja intensidad
para el estrato de bosqque seco
p l u v i e s t a c i o n a l d e l
Ecuador.CEDAMAZ. Ecuador

Zianis D., Mencuccini M., 2003. On simplifying

allometric analyses of forest biomass.
Forest Ecology and Management, 187:

27 Rev. Tayacaja 1 (1): 17-27, enero - julio 2018