Matrices de transición, una herramienta para evaluar la sostenibilidad de sistemas agroforestales
Transition Matrix, a tool to assess the sustainability of agroforestry systems
Wilson Antonio Calero Borge1
Resumen
El trabajo se desarrolló en la Reserva Indígena de Talamanca, Costa Rica. Se midieron 64 parcelas de cacao (C) y banano (B), para determinar la sostenibilidad del aprovechamiento de madera que realizan los agricultores indígenas, midiendo el incremento de madera de laurel (Cordia alliodora) y cedro amargo (Cedrela odorata L.) de regeneración natural en estos sistemas agroforestales a través de matrices de transición como herramienta que permite calcular las tasas de paso o permanencia en diferentes clases diamétricas de individuos en crecimiento. Se realizó inventario de laurel y cedro amargo y se midieron árboles con d ≥ 4 cm. El volumen comercial (Vc) (d ≥ 45 cm) se proyectó con matrices de transición tipo Usher. Se encontró que el aprovechamiento de madera en los SAF evaluados es sostenible, puesto que hay incrementos de madera y de densidad de individuos en los cacaotales, esto indica el buen manejo que realizan los productores indígenas en los sistemas agroforestales. Se recomienda el uso de las matrices de transición como herramienta para evaluar la sostenibilidad de estos sistemas productivos.
Palabras clave: sistemas agroforestales; regeneración natural; matrices de transición; sostenibilidad.
Summary
This work was developed in the Indigenous Reserve of Talamanca, Costa Rica. 64 plots of cocoa (C) and banana (B) were measured to determine the sustainability of the use of wood by indigenous farmers, measuring the increase in wood of laurel (Cordia alliodora) and cedro amargo (Cedrela odorata L.) natural regeneration in these agroforestry systems through transition matrix as a tool that allows to calculate the rates of passage or permanence in different diametric classes of individuals growing. An inventory of laurel and cedro amargo was carried out and trees with d > 4 cm were measured. The commercial volume (Vc) (d > 45 cm) was projected with Usher-type transition matrix. It was found that the use of wood in the evaluated AFS is sustainable, since there are increases in wood and density of individuals in the cacao plantations, this indicates the good management carried out by indigenous producers in agroforestry systems. The use of transition matrix is recommended as a tool to evaluate the sustainability of these productive systems.
Key Words: agroforestry systems; natural regeneration; transition matrix, sustainability
I. Introducción
Los Sistemas Agroforestales (SAF) son formas de uso de la tierra que se caracterizan por contener en sus diversos componentes grandes cantidades de biomasa y abundantes acumulaciones de elementos como el carbono (Nair 1997). Así mismo, el componente arbóreo maderable presente en los SAF es un recurso valioso que brinda diferentes beneficios al agricultor, presta servicios ambientales como captura de carbono, protección del suelo, conservación de biodiversidad, entre otros (Barrance et al, 2003; Albertin y Nair, 2004).
El cacao y banano son cultivos que las comunidades indígenas, en este caso Talamanca en Costa Rica, han cultivado de antaño, bajo diferentes doseles de sombra y constituyen sus principales actividades económicas. El laurel (Cordia alliodora) y cedro amargo (Cedrela odorata) son las principales especies maderables en las fincas de indígenas de Talamanca (Suatunce, 2002). Estas especies son muy apreciadas por su madera, abundante regeneración natural, rápido crecimiento, sombra rala, y autopoda, en el caso del laurel. Los agricultores obtienen ingresos económicos al aprovecharlas en épocas de baja producción o bajos precios del cacao (Albertin y Nair, 2004; Russell y Francell, 2004).
Considerando la diversidad de especies en los SAF, el aprovechamiento de madera que realiza el productor indígena, el reclutamiento y la muerte de individuos, el presente trabajo propone el uso de matrices de transición como una herramienta para estimar la sostenibilidad de los SAF de indígenas de Talamanca.
II. Revisión de literatura
Regeneración y crecimiento del laurel en SAF
El éxito del establecimiento de la regeneración natural de laurel depende de la protección contra animales que depredan las semillas, en potreros, la regeneración natural está muy influenciada por los factores edáficos en sus primeras etapas y por manejo y uso anterior de la tierra en las sucesivas etapas de desarrollo de la planta (Camargo et al. 2000).
Modelos de crecimiento
Un modelo es una abstracción o una representación simplificada de algunos aspectos de la dinámica natural de un rodal y puede involucrar crecimiento, mortalidad y otros cambios en la estructura y composición. Un modelo puede predecir el crecimiento y rendimiento del rodal bajo diversas condiciones y se pueden agrupar asi; a nivel de rodal, de árboles individuales y por clase de tamaños (Vanclay, 1994). Existen tres variaciones de estas aproximaciones: las cadenas de Markov, las matrices de Usher y sus generalizaciones (Vanclay,1994; Moscovich, 2004).
Las cadenas de Markov
En las cadenas de Markov, un árbol en una clase diamétrica, durante el próximo período puede permanecer en la clase, crecer y pasar a la siguiente clase, ser cosechado o morir. Las probabilidades de paso son expresadas como una matriz y son usadas para predecir cambios en un intervalo de tiempo o sobre varios intervalos de tiempo ; donde
es la matriz que contiene las probabilidad de movimiento y
y
son los vectores que describen los estados inicial y final, respectivamente (Vanclay, 1994).
Las matrices de Usher
Las matrices de Usher intentan reducir el largo número de parámetros requeridos para formar una cadena de Markov, seleccionando el intervalo de tiempo y una amplitud de clase en la que un árbol durante ese intervalo no pueda crecer más que una clase, permite una sustancial reducción en el número de parámetros a ser estimados, entonces, la matriz puede ser reducida a cuatro vectores: crecimiento, fecundidad, cosecha y mortalidad (Vanclay, 1994).
Aplicaciones de los modelos matriciales de Usher
Diversas aplicaciones de estos modelos matriciales para árboles en sistemas agroforestales han sido desarrolladas. Somarriba (2001) modeló la dinámica poblacional de Acacia pennatula. La ecuación empleada fue , en donde
y
, denota los vectores de las distribuciones diamétricas en el periodo
y
, y
es el vector de reclutamiento. Los coeficientes de transición de la matriz
,
y
fueron calculados con base a la tasa anual de paso y tasa anual de permanencia, respectivamente, y ajustados por un factor de supervivencia por clase diamétrica.
En la zona de Talamanca, Costa Rica, en sistemas productivos de cacao y banano, para la especie Cordia alliodora en condiciones de regeneración natural, Suárez (2001) desarrolló un modelo de matrices de transición modificada (Somarriba, 2001) del tipo , donde los vectores
y
, denotan la distribución diamétrica entre los años sucesivos
y
y
es la matriz de los coeficientes de transición y de permanencia (Suárez, 2001).
III. Materiales y métodos
Población de estudio
La población de estudio consistió en los cacaotales y bananales de indígenas Bribri y Cabécar, en Talamanca, Costa Rica, para ello se midieron 64 parcelas que han sido medidas cada 5 años.
Modelación de la dinámica poblacional de laurel en sistemas agroforestales
Se proyectó la dinámica poblacional desarrollando un modelo de matrices de transición tipo Usher que proponen Somarriba (2001) y Suárez (2001). La población fue dividida en clases diamétricas de 5 cm. El diámetro mínimo de corta fue establecido en 45 cm.
La iteración se realizó por periodos de un año hasta una proyección de cinco años para determinar la capacidad de producción de madera y comportamiento de la población de laurel de los SAF en estudio. La ecuación matricial empleada fue la siguiente:
Donde los vectores y
denotan la distribución diamétrica (árboles ha-¹) entre los años sucesivos
y
,
es la matriz de coeficientes de transición y
es el reclutamiento de individuos en el periodo que se analiza. La ecuación
devuelve en cada iteración en la clase diamétrica superior la densidad de árboles que alcanzan el diámetro mínimo de corta.
Suponiendo una población de clases diamétricas (z₁, z₂... zl) en el año , en donde zl es la clase diamétrica de los árboles de cosecha, la proyección de la población en el año
con la ecuación
en forma matricial es:
+ R
=
*
Dónde:
- : Coeficientes de transición del periodo
al periodo
eij = pij · Si
- : Coeficientes de permanencia en el periodo
rii = qii · Si
El cálculo de estos coeficientes fue derivado de la estimación de una tasa anual de paso (pij ) de clase a clase
y una tasa anual de permanencia (qii ) en la misma clase
ajustados por un factor de supervivencia anual por clase (Si), entonces:
Donde denota el incremento corriente anual en diámetro en la clase
estimado con la primera derivada de la ecuación de crecimiento de Chapman Richards (Suárez, 2001) y
el intervalo de clase (5 cm).
Dónde:
Parámetros de la ecuación de regresión
: Base de los logaritmos
Edad para alcanzar la categoría diamétrica
- Entonces la tasa anual de permanencia se calculó de la siguiente manera:
qii = 1 - Pij
Dónde:
: Tasa anual de permanencia
: Tasa anual de paso de la clase
a la clase
Cálculo de la supervivencia
La supervivencia anual por clase () se estimó con base en la distribución diamétrica y las tasas de aprovechamiento anual por clase diamétrica para laurel. La ecuación que se empleó fue la siguiente:
Dónde:
: Es la supervivencia en la clase
: es el número de individuos en la clase diamétrica
: es la tasa de aprovechamiento anual para la clase diamétrica
.
Cálculo del reclutamiento anual
Se estimó comparando las curvas de distribución diamétrica de los inventarios sucesivos. La diferencia entre ambas curvas permitió obtener el número de individuos que se incorporaron a las clases diamétricas en el periodo, la ecuación fue la siguiente:
Dónde:
Tasa anual de reclutamiento (árboles ha-¹ año-¹)
Individuos en la clase diamétrica
por hectárea en el inventario actual.
Individuos en la clase diamétrica
por hectárea en el inventario anterior.
Cálculo de la tasa de aprovechamiento de laurel
El volumen de madera aprovechado anualmente en las parcelas se estimó mediante recorridos por la parcela, el propietario indicó los tocones aprovechados y el año de extracción. Se midió la altura de corte y el diámetro a todos los tocones aprovechados en el periodo. El diámetro se midió en dos direcciones tratando de encontrar los puntos intermedios del tocón debido a la presencia de aletones. Posteriormente, el volumen se estimó empleando la ecuación de Somarriba y Beer (1987):
Donde:
: Volumen del fuste con corteza (m³)
: Base de los logaritmos naturales
: Diámetro a la altura del pecho (cm)
Debido a que las alturas de corte son variables, se estandarizó el diámetro de los tocones a una altura de corte de 1,3 m (dap) con la ecuación de Pérez (1954) y Suárez (2001):
Donde:
dt: Diámetro a la altura del tocón (cm)
ht: Altura del tocón (cm)
d: Diámetro del tocón a 1,30 m de altura (cm)
Análisis de escenarios
El análisis de la dinámica poblacional de laurel se realizó a nivel de poblaciones y de parcelas, para evaluar el efecto como resultado de cambios en las prácticas de manejo, se simuló para cinco años el comportamiento con diferentes tasas de sobrevivencia y reclutamiento. Los escenarios fueron: 1) Proyección a cinco años con los datos observados (árboles ha-¹), 2) 100% de sobrevivencia en las clases diamétricas inferiores (dap < 45 cm) y 3) Incrementos del 50% en las tasas de reclutamiento.
Proyección del rendimiento de madera
Se realizó a partir de la distribución poblacional en clases diamétricas, con la siguiente ecuación se estimó la producción de madera (m³ ha-¹ año-¹) para cada sistema agroforestal.
Donde:
Incremento corriente anual del volumen comercial (m³ ha-¹ año-¹)
Vc i + 1: Volumen comercial proyectado en el año siguiente (m³ ha-¹ año-¹)
Vci: Volumen comercial proyectado en el año actual (m³ ha-¹ año-¹)
Evaluación de la sostenibilidad del sistema agroforestal
Se asumió lo siguiente: si el incremento corriente anual del volumen comercial () es mayor que cero el aprovechamiento es sostenible, si es igual a cero está en equilibrio con el crecimiento de la población, si el
es menor que cero el aprovechamiento es no sostenible. También se determinó el incremento periódico anual del volumen comercial
. La ecuación fue:
Dónde:
Incremento periódico anual del volumen comercial proyectado (m³ ha-¹ año-¹).
Vc f: Volumen comercial (m³ ha-¹) al final del periodo proyectado de cinco años.
Vc i: Volumen comercial al inicio del periodo proyectado (m³ ha-¹).
El análisis a nivel de parcelas se realizó calculando las tasas de sobrevivencia y reclutamiento anual, la distribución diamétrica fue la encontrada en cada parcela en este estudio.
Se realizaron iteraciones para cinco años y se encontraron los incrementos en densidad, si el incremento fue positivo se consideró que hay crecimiento de la población, si es negativo entonces hay disminución de la población, se encontraron las tasas anuales de incremento considerando el período de proyección, las ecuaciones fueron las siguientes:
Dónde:
: Cambio en la densidad de árboles en inventarios sucesivos (árboles ha-¹).
: Densidad de árboles a cinco años de proyección (árboles ha-¹)
: Inventario base en la parcela (árboles ha-¹)
Dónde:
: Incrementos periódicos anuales de la densidad (árboles ha-¹ año-¹)
: Cambios en la densidad de árboles en el periodo analizado (árboles ha-¹)
IV. Resultados y discusión
Aprovechamiento de madera de laurel en sistemas agroforestales de Talamanca
Las tasas anuales de corta de laurel (Cuadro 1) son mayores en los cacaotales; en los bananales se extrae el 38% de madera en relación a los cacaotales. En las laderas se extrae menos madera que en los valles (26% con relación al valle). El volumen comercial representó el 84% del volumen total anual extraído en los valles y el 80% en las laderas.
Se están aprovechando 0,92 m³ ha-¹ año-¹ lo que equivale a una extracción de 216 m³ en cinco años en 47 hectáreas, las tasas de extracción en los SAF de cacao y banano son de 1,16 m³ ha-¹ año-¹ y 0,67 m³ ha-¹ año-¹ respectivamente, son superiores a las reportadas por Suárez (2001) para SAF de cacao y banano siendo de 0,89 y 0,11 m³ ha-¹ año-¹ respectivamente.
Cuadro 1. Tasas anuales de corta de madera de laurel (m³ ha¯¹ año¯¹) en cacaotales y bananales en Talamanca, Costa Rica.
|
Clase diamétrica (cm) |
Cacao |
Banano |
||
|
Valle |
Ladera |
Valle |
Ladera |
|
|
m³ ha¯¹ año¯¹ |
m³ ha¯¹ año¯¹ |
m³ ha¯¹ año¯¹) |
m³ ha¯¹ año¯¹ |
|
|
4-10,9 |
0,00 |
0,00 |
0,00 |
0,00 |
|
11-15,9 |
0,00 |
0,00 |
0,02 |
0,00 |
|
16-20,9 |
0,01 |
0,04 |
0,02 |
0,00 |
|
21-25,9 |
0,00 |
0,01 |
0,00 |
0,00 |
|
26-30,9 |
0,02 |
0,01 |
0,02 |
0,00 |
|
31-35,9 |
0,00 |
0,00 |
0,00 |
0,00 |
|
36-40,9 |
0,18 |
0,02 |
0,03 |
0,12 |
|
41-45,9 |
0,07 |
0,07 |
0,05 |
0,00 |
|
46-50,9 |
0,20 |
0,01 |
0,09 |
0,00 |
|
51-55,9 |
0,33 |
0,00 |
0,05 |
0,25 |
|
56-60,9 |
0,18 |
0,00 |
0,00 |
0,61 |
|
61-65,9 |
0,34 |
0,00 |
0,09 |
0,00 |
|
66-70,9 |
0,32 |
0,00 |
0,00 |
0,00 |
|
71-75,9 |
0,14 |
0,00 |
0,00 |
0,00 |
|
76-80,9 |
0,21 |
0,00 |
0,00 |
0,00 |
|
81-85,9 |
0,00 |
0,00 |
0,00 |
0,00 |
|
>86 |
0,17 |
0,00 |
0,00 |
0,00 |
|
Volumen total |
2,17 |
0,16 |
0,37 |
0,98 |
|
V comercial (≥ 45 cm de dap) |
1,89 |
0,01 |
0,23 |
0,86 |
Hay mayor extracción de árboles en los cacaotales, las mayores frecuencias de extracción correspondieron a individuos con dap ≥ 35 cm en cacaotales y bananales. La extracción de individuos jóvenes (dap ≤ 15 cm) es alta y similar en ambos sistemas agroforestales (Figura 1).
Figura 1. Diámetros de árboles de laurel aprovechados anualmente en cacaotales y bananales de Talamanca, Costa Rica.
Proyección de la densidad, rendimiento e incremento del volumen comercial de madera de laurel en los sistemas agroforestales
Las probabilidades de paso son muy semejantes para valle o ladera, ya sea cacao o banano (Cuadro 2). Las tasas de aprovechamiento se consideraron como las únicas fuentes de disminución de la población de laurel y fueron de 0,52 y 0,26 árboles ha¯¹ año¯¹ y el reclutamiento fue de 3,26 y 1,80 árboles ha¯¹ año¯¹ para bananales y cacaotales respectivamente, en general los bananales tienen los mejores resultados en estas dos variables, presentan las menores tasas de aprovechamiento y las mayores tasas de reclutamiento.
Cuadro 2. Probabilidades de paso () y de permanencia (
) por categorías diamétricas para las poblaciones de laurel.
|
Clase diámetrica (cm) |
Cacao |
Banano |
|||||||
|
Valle |
Ladera |
Valle |
Ladera |
||||||
|
I |
j |
|
|
|
|
|
|
|
|
|
4 |
10 |
0,54 |
0,37 |
0,39 |
0,52 |
0,49 |
0,50 |
0,42 |
0,57 |
|
10 |
15 |
0,55 |
0,41 |
0,40 |
0,57 |
0,54 |
0,40 |
0,39 |
0,55 |
|
15 |
20 |
0,53 |
0,44 |
0,38 |
0,60 |
0,54 |
0,45 |
0,38 |
0,61 |
|
20 |
25 |
0,50 |
0,47 |
0,35 |
0,63 |
0,51 |
0,48 |
0,35 |
0,64 |
|
25 |
30 |
0,47 |
0,50 |
0,32 |
0,65 |
0,47 |
0,50 |
0,32 |
0,66 |
|
30 |
35 |
0,44 |
0,54 |
0,29 |
0,70 |
0,44 |
0,55 |
0,29 |
0,70 |
|
35 |
40 |
0,39 |
0,55 |
0,24 |
0,70 |
0,38 |
0,55 |
0,24 |
0,69 |
|
40 |
45 |
0,35 |
0,58 |
0,21 |
0,72 |
0,35 |
0,59 |
0,21 |
0,74 |
|
45 |
50 |
0,30 |
0,60 |
0,17 |
0,73 |
0,31 |
0,62 |
0,17 |
0,76 |
|
50 |
55 |
0,25 |
0,60 |
0,13 |
0,72 |
0,27 |
0,64 |
0,14 |
0,78 |
|
55 |
60 |
0,22 |
0,65 |
0,10 |
0,77 |
0,22 |
0,66 |
0,10 |
0,78 |
|
60 |
65 |
0,19 |
0,68 |
0,06 |
0,80 |
0,18 |
0,67 |
0,06 |
0,78 |
|
65 |
70 |
0,13 |
0,61 |
0,03 |
0,71 |
0,17 |
0,82 |
0,04 |
0,95 |
|
70 |
75 |
0,11 |
0,74 |
0,02 |
0,86 |
0,13 |
0,86 |
0,01 |
0,99 |
|
75 |
80 |
0,07 |
0,72 |
0,09 |
0,90 |
||||
|
80 |
85 |
0,05 |
0,94 |
0,05 |
0,94 |
||||
: coeficientes de paso;
: coeficientes de permanencia; i: clase diamétrica de permanencia; j: clase diamétrica de paso
En todas las poblaciones la curva se desplazará hacia la derecha haciendo moda en diámetros de 30 a 35 cm (Figura 2), lo que indica que no se producirán incrementos considerables, esto puede provocar una distribución poblacional donde las clases diamétricas inferiores tengan muy pocos individuos para mantener una adecuada estructura de individuos jóvenes y diámetros aprovechables.
Figura 2. Proyección de la población de laurel para cinco años en a) Cacao ladera, b) Cacao valle,
c) Banano ladera y d) Banano valle en Talamanca, Costa Rica.
La proyección del rendimiento de madera en volumen comercial a nivel de poblaciones fue creciente en todos los sistemas agroforestales, los mayores rendimientos se producirán en los sistemas agroforestales en valle, a cinco años de proyección los SAF en laderas estarán produciendo un 30% de madera menos en relación a los SAF en valle (Cuadro 3).
Cuadro 3. Proyección de las existencias (m³ ha-1) de madera comercial de laurel.
|
SAF |
Año base |
Año 1 |
Año 2 |
Año 3 |
Año 4 |
Año 5 |
|
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
|
|
CL |
30,40 |
33,80 |
35,09 |
37,96 |
41,02 |
44,31 |
|
CV |
33,11 |
34,89 |
40,15 |
46,07 |
52,64 |
59,85 |
|
BL |
28,70 |
31,34 |
34,39 |
37,59 |
41,14 |
44,86 |
|
BV |
32,80 |
35,52 |
41,50 |
48,15 |
55,50 |
63,65 |
CV: cacao en valle; CL: cacao en ladera; BV: banano valle; BL; banano ladera
El ICAVc fue creciente en todos los sistemas agroforestales a excepción del cacao en ladera que disminuyó en el tiempo, los BV fueron los que presentaron los mayores ICAVc. El incremento periódico anual del volumen comercial (IPAvc) será mayor en los sistemas agroforestales en valle (Cuadro 4). A nivel de población en cacaotales y bananales el aprovechamiento es sostenible puesto que hay incrementos en la producción de madera.
Cuadro 4. Proyección de los incrementos volumétricos de madera comercial (dap ≥ 45 cm) de laurel.
|
SAF |
ICAVc año1 |
ICAVc año 2 |
ICAVc año 3 |
ICAVc año 4 |
ICAVc año 5 |
IPAvc |
|
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
|
|
CL |
3,40 |
1,29 |
2,87 |
3,06 |
3,29 |
2,78 |
|
CV |
1,78 |
5,26 |
5,92 |
6,58 |
7,21 |
5,38 |
|
BL |
2,64 |
3,05 |
3,20 |
3,55 |
3,72 |
3,23 |
|
BV |
3,14 |
5,98 |
6,65 |
7,35 |
8,15 |
6,25 |
CV: cacao en valle; CL: cacao en ladera; BV: banano valle; BL: banano ladera
La densidad de laurel disminuirá en el 40 y 50% de las parcelas de banano y cacao respectivamente. Mientras que el 14 y el 17% de las parcelas de banano y cacao respectivamente tendrán altos incrementos de densidad; este resultado sugiere que el incremento de la densidad que se observa a nivel de poblaciones de SAF es producto del alto incremento de pocas parcelas en los sistemas agroforestales evaluados, mientras que en la mayoría de las parcelas la densidad está disminuyendo o tienen incrementos bajos (Figura 3).
Figura 3. Comportamiento de la densidad de laurel a nivel de parcelas con una proyección a cinco años.
Análisis de escenarios
Sin aprovechamiento de individuos ≤ 45 cm de dap
En este escenario, considerando la sobrevivencia del 100% de los individuos con dap ≤ 45 cm (Cuadro 5) se obtiene que en todos los sistemas agroforestales habrán incrementos en el ICAVc y el IPAvc, siendo los bananales los de mayores incrementos (13% más que los cacaotales), sin embargo, los incrementos son bajos con respecto a la proyección con los datos observados.
Cuadro 5. Proyección de los incrementos volumétricos IPAvc de madera comercial de laurel por 5 años
|
SAF |
Año 1 ICAVc |
Año 2 ICAVc |
Año 3 ICAVc |
Año 4 ICAVc |
Año 5 ICAVc |
IPAvc |
|
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
|
|
CL |
2,17 |
2,63 |
2,95 |
3,87 |
2,84 |
2,99 |
|
CV |
4,53 |
5,35 |
6,18 |
6.82 |
7,59 |
6,09 |
|
BL |
2,65 |
3,13 |
3,27 |
3,67 |
3,88 |
3,32 |
|
BV |
5,14 |
6,09 |
6,62 |
7,51 |
8,53 |
6,82 |
CV: cacao en valle; CL: cacao en ladera; BV: banano valle; BL; banano ladera
Incrementos del 50% de reclutamiento de individuos
En este escenario, se obtiene un incremento en el Vc de madera en todos los SAF evaluados siendo los bananales los que mayores incrementos presentan, no obstante, los sistemas en laderas no presentan incrementos cuando se compararan con las proyecciones de los datos observados (Cuadro 6).
Cuadro 6. Proyección de los incrementos volumétricos IPAvc de madera comercial de laurel por cinco años en Talamanca, Costa Rica.
|
SAF |
Año 1 ICAVc |
Año 2 ICAVc |
Año 3 ICAVc |
Año 4 ICAVc |
Año 5 ICAVc |
IPAvc |
|
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
(m³ ha-1) |
|
|
CL |
3,30 |
1,20 |
2,70 |
3,00 |
2,99 |
2,67 |
|
CV |
4,53 |
5,35 |
6,18 |
6,82 |
6,09 |
6,09 |
|
BL |
2,60 |
3,06 |
3,20 |
3,55 |
3,32 |
3,23 |
|
BV |
5,10 |
5,98 |
6,66 |
7,35 |
6,82 |
6,65 |
CV: cacao en valle; CL: cacao en ladera; BV: banano valle; BL; banano ladera
Proyección del rendimiento e incremento de madera y aplicación del modelo
Los incrementos de madera tanto en cacaotales como en bananales son inferiores (2,8; 5,4; 3,2 y 6,25 m³ ha-¹ año-1 para CL, CV, BL y BV respectivamente) a los reportados en otros estudios realizados con laurel. Rosero y Gewald (1979) reportan incrementos en volumen de 14,8 m³ ha-1 año-1; Somarriba (1990) reporta incrementos en volumen de 9 – 24 m³ ha-1 año-1. Sin embargo, los resultados obtenidos en este estudio son superiores a la proyección realizada por Suárez (2001) para estas mismas parcelas.
La distribución de la población de laurel con datos observados muestra que la clase diamétrica de 20-25 cm tiene la mayor densidad de individuos y luego desciende hacia los diámetros superiores, siendo esto una característica de las poblaciones de regeneración natural.
La poca presencia de individuos en las clases diamétricas inferiores a 20 cm (figura 2), puede deberse a factores ecológicos y de manejo, los factores ecológicos estarían asociados a: la biomasa en descomposición sobre el suelo la cual incide en la viabilidad de la semilla; la sombra del dosel de árboles y del cacao y el periodo de lluvias puesto que el laurel germina en época seca y estas arrastraran muchas semillas. Los factores de manejo están asociados a: las frecuencia de las chapias que realiza el agricultor en la cual son eliminadas muchas plantitas; al raleo para regular sombra para el cacao o banano y que implica eliminar árboles pequeños, a partir de dap > 12 cm este raleo disminuye. En la proyección de la población de laurel para cinco años se observa la misma tendencia (Figura 2) e incluso la curva hace moda en 25 cm.
En el análisis de los dos escenarios no se observan incrementos importantes en la producción de madera con respecto a los datos observados. Es probable que el periodo de análisis de cinco años es muy poco tiempo para que se exprese en el volumen de madera de laurel el efecto de mayores reclutamientos o el no aprovechamiento de individuos menores de 45 cm de dap. Con estas proyecciones, el sistema es sostenible porque hay incrementos de madera.
El modelo desarrollado puede considerarse como un reajuste al modelo de Suárez (2001), este consideró el reclutamiento como una nueva variable, calculó la sobrevivencia a partir de dos momentos de evaluación e incorporó datos de individuos desde las primeras clases diamétricas, por lo que el modelo gana en eficiencia.
V. Conclusiones
• El aprovechamiento de madera en estos sistemas agroforestales es sostenible, puesto que ocurren incrementos tanto en volumen de madera como en la densidad poblacional.
• Los sistemas agroforestales de cacao y banano orgánico de indígenas de Talamanca son fuentes importantes de madera.
• En el periodo proyectado con las matrices de transición, con las actuales tasas de extracción y reclutamiento, las poblaciones de laurel en los SAF aumentarán en (volumen y densidad).
• En los escenarios evaluados (disminución de la mortalidad e incrementos de las tasas de reclutamiento) los incrementos serán crecientes en todos los SAF.
• Los sistemas agroforestales evaluados son sostenibles con las actuales tasas de aprovechamiento y manejo.
• El modelo aplicado es una herramienta útil que permite evaluar la sostenibilidad de los sistemas productivos diversificados como son los SAF.
VI. Lista de referencias
Albertin, A; Nair, PK. (2004). Farmers Perspectives on the role of shade Trees in Coffee Productions Systems: An assessment from the Nicoya Peninsula, Costa Rica. Human Ecology 32(4): 443-463.
Barrance, AJ; Flores, L; Padilla, E; Gordon, JE; Schreckenberg, K. (2003). Trees and farming in the dry zone of southern Honduras I: campesino tree husbandry practices. Agroforestry Systems 59: 97-106.
Camargo, JC; Ibrahim, M; Somarriba, E; Finegan, B; Current, D. (2000). Factores ecológicos y socioeconómicos que influyen en la regeneración natural del laurel en sistemas silvopastoriles del trópico húmedo y subhúmedo de Costa Rica. Agroforestería en las Américas 7(26):46-49.
Moscovich, FA. (2004). Modelos de crecimiento y producción forestal. INTA. Centro Regional Misiones, Estación Experimental Montecarlo, Argentina. 42p.
Nair, PK. (1997). Agroforestería. Centro de Agroforestería para el desarrollo sostenible. Universidad Autónoma de Chapingo, México. Editor L. Krishnamurthy. 543 p.
Pérez, CF. (1954). Estudio forestal del Laurel Cordia alliodora (Ruiz y Pavón) Cham., en Costa Rica. Tesis Mag. Sc. Turrialba, Costa Rica, IICA. 195 p.
Rosero, P; Gewald, N. (1979). Growth of laurel (Cordia alliodora) in coffe and cacao plantations, and pastures, in the Atlantic region of Costa Rica. In De las Salas G ed. Proceeding of the Workshop Agro-forestry Systems in Latin America. Turrialba, Costa Rica CATIE. pp 205-208.
Russell, D; Franzel, S. (2004). Trees of prosperity: Agroforestry, markets and the African smallholder. Agroforestry Systems 61: 345-355.
Somarriba, E. (2001). Acacia pennatula en los potreros de la reserva natural Mesas de Moropotente, Estelí, Nicaragua. 28 p.
Somarriba, E. (1990). Sustainable timber production from uneven-aged shade stands of Cordia alliodora in small coffe farms. Agroforestry Systems 10: 253-263.
Somarriba, E; Beer, J. (1987). Dimensions, volumes and growth of Cordia alliodora in agroforestry systems. Forest Ecology and Management 18(2):113-126.
Suárez, A. (2001). Aprovechamiento sostenible de madera de Cordia alliodora y Cedrela odorata de regeneración natural en cacaotales y bananales de indígenas de Talamanca, Costa Rica. Tesis Mag. Sc. Turrialba, Costa Rica, CATIE. 77 p.
Suatunce, P. (2002). Diversidad de escarabajos estercoleros en bosques y en cacaotales de diferentes estructuras y composición florística, Talamanca, Costa Rica. Tesis Mag. Sc. Turrialba, Costa Rica, CATIE. 135 p.
Vanklay, J. (1994). Modelling forest growth and yield Applications to Mixed Tropical Forests. Cab International, Oxford, UK. 312 p.
1 Master en Agroforestería Tropical, Coordinador de Investigación y Postgrado de la Universidad de las Regiones Autónomas de la Costa Caribe Nicaragüense-Recinto Universitario Nueva Guinea. Correo: investigació[email protected] ORCID: https://orcid.org/0000-0003-4243-0490