INTRODUCCIÓN

Vochysia guatemalensis (Donn. Smith) es una de las principales especies forestales nativas utilizadas en la reforestación comercial en Costa Rica. Para el 2004, el área plantada en la región Huetar Norte, alcanzó el 25% del total reforesta- do en todo el país con esta especie (^{Solís y Moya 2004}), lo cual correspondía a 947 hectáreas; 10 años más tarde en el 2014 alcanzó 2116 hectáreas, que corresponde a un 64% del área plantada con esta especie en todo el país (INEC 2014). Debido a su facultad de acumular aluminio en el follaje, así como en otros tejidos del árbol (Chenery y Spor- ne 1976, Foy et al. 1978, citados por ^{Pérez 1993}), se le considera una especie hyperacumuladora o tolerante al aluminio (^{Camacho et al. 2017}). Por esta cualidad y por su capacidad para ocupar terrenos con pastos degradados (^{Montagnini y Sancho 1994}), se le reconoce un gran potencial para crecer mejor que otras especies forestales en los suelos ácidos y poco fértiles de la región tropical. ^{Camacho et al. (2016}) afirman que los árboles de esta especie tienen una ventaja compe- titiva, por encima de Tectona grandis y Gmelina arborea, debido a su mejor adaptabilidad a suelos bajos en pH y altos en aluminio. Por otro lado, V. guatemalensis presenta mayores tasas de creci- miento y productividad cuando se le compara con otras especies forestales nativas del país (^{Delgado 2002}, ^{Arias et al. 2011}). Según lo anterior, ha dado lugar para que sea propuesta como una especie complementaria a las especies exóticas que dominan el mercado en los programas de reforestación locales.

Crecimiento en altura y clasificación de sitios

Los estudios para evaluar el desarrollo en altura de las especies forestales nativas han alcanzado en los últimos años, mayor importan- cia práctica, debido a su relación directa con el cálculo del volumen total de madera (^{Andrade 2002}). En el caso de V. guatemalensis los pocos estudios que existen por lo general, han repor- tado alturas hasta edades muy tempranas, como por ejemplo, el estudio presentado por ^{ACEN (1992}), en donde se hace referencia a los datos de crecimiento en altura del proyecto de Ensayos Forestales (TRIAL) desarrollado por la OET. Esta deficiencia en la información limita las pro- yecciones del volumen para esta especie a edades más avanzadas (^{Arias et al. 2011}).

En Costa Rica, los trabajos sobre calidad de sitio y clasificación de sitios para V. guate- malensis son muy escasos. De hecho, solamente se conocen unos pocos estudios publicados, uno de ellos realizado por ^{Barraza y Díaz (1999}), quienes de manera preliminar clasificaron 23 localidades en las cuales se encontraba cre- ciendo la especie. Para construir el sistema de clasificación de sitios los autores definieron una edad base de 11 años y establecieron 3 clases de calidad de sitio: Baja (IS entre 24 y 27), Media (IS entre 27 y 33) y Alta (IS entre 33 y 36). Por su parte, COSEFORMA (2001) publicó una guía sobre V. Guatemalensis, en la cual se presenta un esquema para la clasificación de las plantaciones a través de la estratificación de las tasas de cre- cimiento diamétrico, lo cual permitió construir 3 escenarios de crecimiento en diámetro para esta especie. No obstante, este esquema de clasifica- ción, basado en el crecimiento diamétrico, tiene su mayor aplicación en la definición de propues- tas para el manejo silvícola de esta especie, más que para una clasificación de sitios.

Más recientemente, ^{Mora et al. (2015}) utilizaron una clasificación de sitios preliminar para evaluar los factores edáficos que influyen en el crecimiento de esta especie en la región Huetar Norte de Costa Rica. En la presente investigación, se utilizaron variables físicas y químicas de 2 tipos de suelos (Inceptisoles y Ultisoles), con el fin de explicar el comportamiento del índice de sitio (IS) y su dependencia de aquellas variables edáficas que mostraron algún tipo de relación con la calidad de los terrenos (fuera esta positiva o negativa) para V. guatemalensis.

Se propuso evaluar los patrones de creci- miento de la altura dominante y desarrollar un sistema para clasificar la calidad de los terrenos plantados con V. guatemalensis, al utilizar el crecimiento de la altura dominante como un indi- cador de la calidad de sitio. Esta última constituye un concepto integrador, que resume las variadas condiciones que brindan los distintos lugares para el crecimiento de los árboles (^{Prodan et al. 1997}). Bajo este punto de vista, la medida de la calidad de sitio refleja el potencial de productividad o capacidad productiva de un sitio, por lo que se convierte en una expresión de la capacidad de un área de terreno para producir madera o cualquier otro tipo de vegetación (^{Davis y Johnson 1987}).

El objetivo del presente trabajo fue selec- cionar un modelo adecuado, de acuerdo con los criterios señalados por ^{Curtis (1964}) y ^{Bailey y Clutter (1974}), ^{Mora et al. (2014}), que represente acertadamente el crecimiento de la altura domi- nante para Vochysia guatemalensis y que permita construir curvas actualizadas de índice de sitio (IS), con el fin de facilitar la toma de decisiones relacionadas con el manejo silvicultural de esta especie, en la regiones Huetar Norte y Huetar Caribe de Costa Rica.

MATERIALES Y MÉTODOS

Área de estudio

Para llevar a cabo este trabajo fue nece- sario compilar la información existente acerca del crecimiento de V. guatemalensis, proveniente de diferentes organizaciones públicas y privadas, como FUNDECOR, EARTH, COSEFORMA y TEC. El área de estudio abarca las regiones Huetar Norte y Huetar Caribe de Costa Rica, donde las Parcelas Permanentes de Muestreo (PPM) de COSEFORMA-TEC se distribuyen principalmente en los cantones de Los Chiles y San Carlos; las de FUNDECOR se concentran en los cantones de Sarapiquí y Pocosí; y las de la EARTH se ubican en los cantones de Guácimo y Siquirres (Figura 1).

Fig. 1 Distribución de las parcelas permanentes y el área de cobertura del estudio para cada organización involucrada, regiones Huetar Norte y Huetar Caribe, Costa Rica. 

Recopilación de los datos

El proceso de compilación dio como resultado una base de datos compuesta por los

promedios de 43 parcelas permanentes de mues- treo (PPM), las cuales produjeron 397 pares de datos altura dominante-edad. El área total del estudio abarca 4233 km2 (Cuadro 1).

La mayoría de las parcelas fueron medidas desde 1 año hasta 16,5 años (esta edad se consi- dera adecuada, ya que los turnos de corta para esta especie no superan los 12 años). En todas las parcelas se registró la edad para cada medición, se midió el diámetro (Dap) en centímetros y la altura total (Htot) en metros; además, se calculó la altura dominante como el promedio en altura total de los árboles más altos en un número equi- valente al 10% de los árboles de la parcela.

Por tratarse de una base de datos cons- truida a partir del trabajo de 4 instituciones, fue necesario llevar a cabo una revisión y depuración de la misma en busca de inconsistencias. Se trató de corregirlas en todos los casos, cuando esto fue posible, y se removieron de la base de datos cuando no se pudo aplicar la corrección deseada.

Identificación de los patrones de crecimiento mostrados por la altura dominante

Con el propósito de encontrar evidencia empírica de la existencia de clases de calidad de sitio en las localidades, donde se encuentra cre- ciendo V. guatemalensis, fue necesario seguir los siguientes pasos:

Se graficaron los datos de altura dominan- te en metros (definida como el promedio en altura total del 10% de los árboles más altos de la parcela) versus la edad en años, con el fin de estudiar su comportamiento y detectar las diferencias existentes entre las parcelas a la edad de 8 años. Esto permitió identificar, de manera preliminar, 3 clases de calidad de sitio basados en el creci- miento de la altura dominante a la edad mencionada. Esta clasificación constituye un paso preliminar hacia la estratificación definitiva de la muestra.

Posteriormente, se graficaron los datos de las 3 clases por separado, para realizar el análisis del crecimiento en altura y estu- diar el comportamiento de varios modelos lineales y no lineales con cada uno de los grupos. El proceso permitió detectar la existencia de 2 patrones de crecimiento (2 maneras de crecer) para V. guatemalensis.

Se procedió a conformar un primer grupo con las clases I y II (llamada Muestra 1), ya que mostraron un comportamiento muy similar. El segundo grupo quedó confor- mado por los datos de la clase III (Muestra 2), los cuales difieren en su comportamien- to con respecto a los datos de la Muestra 1.

Durante esta fase, se optó por cambiar la edad de referencia de 8 a 9 años, debido a que esta última estaba bien representada en los datos, mostraba una desviación están- dar menor y se encontraba más próxima al turno de rotación.

Para determinar si ambos grupos o mues- tras provienen de una misma población, o si representan poblaciones diferentes, se aplicó la prueba inferencial de ^{Levene (1960}).

Posteriormente, para dilucidar si ambas muestras son o no estadísticamente diferentes, se aplicó la prueba t para la igualdad de medias.

Con el fin de poner en evidencia las discre- pancias encontradas para ambas muestras, se procedió a crear un gráfico altura domi- nante-edad que muestra las diferencias en el comportamiento de ambos grupos de datos.

Para apoyar la teoría estadística, se aportó evidencia empírica al graficar el incremen- to corriente anual (ICA) y el incremento medio anual (IMA) de la altura dominante contra la edad.

Finalmente, ante la imposibilidad de encontrar un modelo lo suficientemente flexible para representar ambos patrones de crecimiento, se resolvió trabajar con cada grupo por separado y utilizar un modelo diferente para ajustar las curvas de índice de sitio, utilizando un enfoque polimórfico en ambos casos. Los 2 juegos de curvas se representan en un mismo gráfico, con el fin de contrastar la forma de las curvas y las tendencias que ambas muestras representan.

Selección del modelo de mejor ajuste para la construcción de las curvas de índice de sitio

Para cada grupo de datos se probaron los mismos modelos de regresión, lo cual permitió seleccionar el modelo que mejor ajustó los datos, tanto para la Muestra 1 como para la Muestra 2. Para evaluar cada modelo se tomaron en cuenta los siguientes estadísticos: r = coeficiente de correlación; R2 = coeficiente de determinación; R2adj = coeficiente de determinación ajustado; MSE = cuadrado medio del error; SEE = error estándar de las estimaciones; y MAE = error medio absoluto. Además, se graficaron los valo- res estimados versus los valores observados, con el fin de evaluar gráficamente el ajuste de los datos.

En una primera fase, se trabajó con los siguientes modelos lineales: Logarítmico, Semi- logarítmico, Cuadrático, Doble Recíproco, Raíz cuadrada de Y-Logaritmo de X y la ecuación de ^{Schumacher (1939}) en su forma lineal. También se probaron 2 modelos no lineales, pero linea- lizables: Prodan (1951), citado por ^{Thompson et al. 1984}, y ^{Bailey y Clutter (1974}), los cuales, a través de una serie de transformaciones, es posi- ble ajustarlos mediante técnicas de cuadrados mínimos lineales.

En una segunda fase se trabajó con los siguientes modelos no lineales: Mistcherlich o Monomolecular, Weibull, Gompertz y Chapman- Richards. Todos estos modelos se pueden obtener a partir de una misma función, la ecuación de von Bertanlanffy también fue ensayada con los datos de la muestra. A esta lista se agregaron 6 modelos más, los cuales constituyen variantes de algunos de los modelos anteriores. A partir del modelo de Gompertz se obtiene el modelo Gom- pertz modificado; de la ecuación ^{de Richards (1959)} se obtienen el Modelo Natural de Cre- cimiento cuando c = 1 y la Función Logística Generalizada cuando c = 1; de esta última se puede obtener la Función Logística o Autocata- lítica cuando d = 1. También se probaron 2 ver- siones de ambas funciones logísticas, reescritas y recompuestas por ^{Sit y Poulin-Costello (1994}). Finalmente, se incluyó el modelo de ^{Schumacher (1939}) en su forma no lineal.

Desarrollo teórico-simbólico de los modelos utilizados para la construcción de las curvas de índice de sitio

Debido a las discrepancias irreconcilia- bles observadas en los patrones de crecimiento que presenta V. guatemalensis, en las regiones Huetar Norte y Huetar Caribe de Costa Rica, se procedió a ajustar modelos separados para cada tendencia. El modelo de Bailey y Clutter se ajustó para el grupo de parcelas que compo- nen la Muestra 1 (clases I y II) y el modelo de Prodan se utilizó para los datos de la Muestra 2 (clase III). Con el fin de elevar el valor didácti- co del actual trabajo, las funciones utilizadas y su desarrollo teórico-simbólico, se presentan a continuación.

Modelo de Bailey y Clutter

Para los datos altura-edad correspondien- tes a las plantaciones de V. Guatemalensis repre- sentadas en la Muestra 1 (clases I y II), el modelo utilizado fue la versión generalizada del modelo de ^{Schumacher (1939}), modificado por ^{Bailey y Clutter (1974}). Su desarrollo simbólico se presen- ta a continuación:

Ln(H) = a-b(1/E)^c (1)

Donde, a, b y c son los coeficientes a ser ajustados H es la altura dominante (m), E representa la edad (en años) y Ln es el logarit- mo natural.

Este modelo, al igual que el de Schu- macher, asume que: “el logaritmo de la altura dominante es una función lineal del recíproco de la edad”, sin embargo, se diferencia en que este incorpora un tercer coeficiente c, que actúa como exponente de la variable independiente. Esta con- dición le confiere varios atributos importantes: por un lado, se vuelve más flexible, al ser capaz de asumir formas diferentes, lo cual le permite ajustar datos de muy diversos orígenes; por otro lado, la forma de las curvas se mantiene aunque cambie la edad base (^{Bailey y Clutter 1974}). En todo caso, según ^{Alder (1980}), para la mayoría de las especies, los valores apropiados de c varían entre 0,2 y 2,0.

Es necesario aclarar que este modelo no lineal puede ser ajustado linealmente por medio de un proceso de 2 etapas (^{Bailey y Clutter 1974}): Etapa 1. Se ajusta la primera derivada de la ecuación de regresión para estimar ĉ. La derivada de esta ecuación es de la forma:

Ln(I/H) = Ln (-cb)+(c+1)Ln(1/E) (2)

Donde, Ln(-cb) corresponde al punto de intercepción y (c+1) equivale a la pendiente del modelo; I es el incremento corriente anual en altura dominante; H y E fueron previamente definidos. Como esta igualdad es lineal en sus parámetros, los procedimientos normales de regresión van a producir una estimación de ĉ+1, del cual ĉ puede ser derivado al restar el valor de 1 a la pendiente de esta ecuación (^{Bailey y Clut- ter 1974}). Para aplicar este método se requiere información producida a partir de mediciones sucesivas realizadas en parcelas permanentes, o bien, por análisis fustal, pues es necesario contar con datos del incremento corriente anual en altura para cada una de las edades (^{Mora y Meza 2003}).

Etapa 2. Se utiliza ĉ en la ecuación básica (1) y se ajustan los datos por medio de la regre- sión lineal simple,

Ln(H) = a-b(1/E)^ĉ

Ahora bien, antes de construir las curvas de índice de sitio, es necesario recordar la defini- ción del índice de sitio, la cual dice que: ”cuando la edad (E) asume el valor de la edad base (Eb), la altura dominante (H) es igual al índice de sitio (IS)”. Al tener en cuenta este enunciado, se procedió a sustituir estas 2 variables en la ecua- ción anterior, por lo que quedó expresada de la siguiente forma,

Ln(IS) = a-b(1/Eb)^ĉ (3)

Para construir curvas polimórficas de índice de sitio es necesario poner a variar el parámetro b, al cual llamaremos “coeficiente de sitio” (bi), ya que este varía con el índice de sitio. De este modo, despejando bi en la ecuación (3), se tiene,

-bi = (a-Ln(IS))(Eb)^ĉ (4) Sustituyendo en (1), se puede tener,

Ln(H) = a+{(a-Ln(IS))(Eb^ĉ)}(1/Ei)^ĉ

= a+(a-Ln(IS)){Eb/Ei}^ĉ (5)

Se eliminan los logaritmos naturales (base

e) y se obtiene,

H = e^a (IS/e^a)^{Eb/Ei}^ĉ (6)

Donde e es la base de los logaritmos naturales (e = 2,71828); IS es el índice de sitio a escoger; a y ĉ fueron definidos en la ecuación básica; Ei asume todos los valores de la edad; y Eb es la edad base (o edad clave), que en este caso es igual a 9 años.

Modelo de Prodan

En el caso de los datos que componen la Muestra 2 (clase III), se utilizó el modelo de Prodan. El desarrollo simbólico del mismo se presenta a continuación:

H = (E2) / (a+bE+cE2) (7)

Donde a, b y c son los coeficientes a ser ajustados; H es la altura dominante (m) y E repre- senta la edad (en años). Para ajustar este modelo, mediante mínimos cuadrados, es necesario remo- ver la variable E2 del numerador de la fracción, de manera que al realizar el despeje esta variable pasa al lado contrario del igual,

E2/H = a+bE+cE2 (8)

Con el fin de facilitar la manipulación de esta función, se hace necesario crear una nueva variable, la cual llamaremos Z; de forma tal que, Z = (E2/H). Entonces,

Z = a+bE+cE2 (9)

Tomando en cuenta la definición del índice de sitio, señalada en la página anterior, se sabe que cuando E = Eb, entonces H = IS. De esta manera, sustituyendo en (8), se tiene,

Eb2/IS = a+bEb+cEb2 (10)

Para construir curvas polimórficas con este modelo b pasa a depender del índice de sitio, por eso se denomina “coeficiente de sitio” (bi) y el mismo puede ser despejado de la ecuación anterior de la siguiente manera,

bi = (Eb/IS)-(a/Eb)-cEb (11)

Al sustituir bi en la fórmula (7) se obtiene la siguiente función,

H = Ei2/(a+(Eb/IS)-(a/Eb)-(cEb) Ei+(cEi2)) (12)

Donde a y c son los mismos coeficientes de la ecuación original; Ei es cualquier edad; y Eb es la edad clave (o edad base), que para fines del presente estudio es igual a 9 años.

RESULTADOS Y DISCUSIÓN

Definición de las categorías de crecimiento en altura dominante

El análisis de los datos permitió identificar 2 patrones de crecimiento claramente diferencia- dos; esto por cuanto los datos definen 2s zonas en el gráfico (Figura 2). En la franja superior se encuentran los valores que representan los mejo- res sitios y corresponden a la Muestra 1 (clases I y II), y en la banda inferior se localizan las parcelas de peor crecimiento, correspondientes a la Muestra 2 (clase III). Mientras que los datos de la Muestra 1 presentan una mayor pendiente y dan muestras de un crecimiento sostenido, los datos de la Muestra 2 presentan una pendiente más leve y una tendencia a buscar la horizon- talidad a una edad muy temprana; esto sugiere que pronto podrían desacelerar aun más su lento crecimiento

Fig. 2 Dispersión de los datos que muestran las discrepancias observadas en el patrón de crecimiento de la altura dominante para Vochysia guatemalensis, en las regiones Huetar Norte y Huetar Caribe de Costa Rica. 

En las primeras edades no se distingue con claridad un patrón diferenciado en la forma de crecer de V. guatemalensis, pero después de la edad de 8 años empieza a visualizarse una ten- dencia más clara en los promedios de la Muestra 2 (clase III), los cuales comienzan a separarse de la masa de datos que conforman la Muestra 1 (clases I y II). A partir de esa edad se observa que las parcelas de la Muestra 2 comienzan a mostrar una reducción en la pendiente y la evidencia sugiere que muy pronto podrían dejar de crecer y hasta podrían alcanzar su asíntota a una edad muy temprana. Los sitios donde se encuentran creciendo las parcelas de la Muestra 1, se pueden considerar de muy buena calidad para esta espe- cie, mientras que los sitios que corresponden a la Muestra 2 se pueden considerar como malos o muy malos.

Bajo este esquema, a la edad base de 9 años la Muestra 1 está delimitada entre las alturas 19,0 m y 28,13 m; la Muestra 2, que representa a los sitios de peor calidad, queda definida entre las alturas 12,0 m y 18,0 m.

Como se puede ver en la Figura 2, a la edad de referencia, la separación entre los 2 gru- pos de datos es muy leve, ya que la brecha entre ambas muestras se reduce a solamente un metro de distancia. Sin embargo, a la edad de 18 años el panorama cambia. La curva que representa el límite inferior de la Muestra 2 cierra con una altura de 20,0 m y la curva del límite superior de termina en 27,85 m a la edad señalada. Por otro lado, la curva que representa el límite inferior de la Muestra 1 cierra con una altura de 33,70 m a la edad de 18 años, mientras que la curva que actúa como límite superior alcanza una altura de 43,11 m a dicha edad. De esta manera, entre el límite inferior de la Muestra 1 y el límite superior de la Muestra 2, existe una brecha de casi 6,0 m de diferencia (5,85 m) a la edad de 18 años, lo que confirma el crecimiento diferencia- do para ambas muestras.

Teoría estadística que respalda las diferencias encontradas en la muestra

Prueba de Levene para la igualdad de varianzas

Al aplicar la prueba de Levene se pudo determinar que esta no rechazó la hipótesis de igualdad de varianzas, lo cual indica con una alta probabilidad que las muestras proceden de una misma población. Como se puede ver en el Cuadro 2, el p-valor de 0,633 es mayor que el nivel de significancia de 0,05, escogido para trabajar en este caso; lo que indica que no se puede rechazar la hipótesis nula.

Prueba t para la igualdad de medias

La prueba t rechazó la hipótesis nula de igualdad de medias, ya que el p-valor (Sig. bila- teral) de 0,000 es menor al nivel de significancia de 0,05, elegido para trabajar con estos datos (Cuadro 3). Además, el intervalo de confianza no incluyó el valor de cero, lo cual es una razón más para desechar la hipótesis nula. La evidencia indica que realmente existen diferencias signifi- cativas en el promedio de las alturas dominantes para ambas muestras.

Representación de los patrones de crecimiento de la especie

Para sustentar los resultados estadísticos obtenidos a través de la prueba de Levene y de la prueba t, se analizó el comportamiento de la variable altura dominante, en su desarrollo a tra- vés del tiempo, para las 2 muestras. Para esto se procedió a graficar el incremento corriente anual (ICA) y el incremento medio anual (IMA) de la altura dominante versus la edad.

La Figura 3 presenta la relación entre el Incremento Corriente Anual (ICA) y la Edad para V. guatemalensis al crecer en 2 calidades de sitio. En estas curvas se manifestaron diferencias atribuibles a la capacidad productiva de los terre- nos. En el mejor sitio, la curva de la Muestra 1 inició desde una posición superior cercana a los 3,0 m, por lo que alcanzó un ICAmáx de 2,77 m a la edad de 3 años. Por su parte, la curva de la Muestra 2 arranca lentamente, con un ICA de 1,16 m y solamente logra alcanzar un ICAmáx de 1,98 m a la edad de 6 años.

Fig. 3.  Relación ICA(Hdom)-Edad para las muestras 1 y 2 de Vochysia guatemalensis, región Huetar Norte y Huetar Caribe, Costa Rica. 

Las curvas representadas en la Figura 4 muestran la relación que existe entre el Incre- mento Medio Anual (IMA) y la Edad, para V. guatemalensis en 2 condiciones de calidad de sitio. Se puede ver claramente la posición superior que mantiene la primera curva, per- teneciente a la Muestra 1, a lo largo de todo el recorrido sobre la segunda curva perteneciente a la Muestra 2. La primera presenta una forma más aplanada, alcanzando el IMA máximo a la edad de seis años; por su parte, la segunda curva no solo se aparta de la forma mostrada por la primera, sino que alcanza su IMA máximo a la edad de 8 años.

Fig. 4.  Relación IMA (Hdom)-Edad para las muestras 1 y 2 de Vochysia guatemalensis, región Huetar Norte y Huetar Caribe, Costa Rica. 

Con estas curvas se pone en evidencia las diferencias atribuibles a la calidad de sitio. En el mejor de los sitios, la curva de la Muestra 1 arranca en una posición superior con 2,44 m, alcanzando un IMAmáx de 2,60 m a la edad de 6 años. La curva de la Muestra 2 inició y lamen- tablemente, con un IMA de 1,32 m, y solamente logró llegar a un IMAmáx de 1,77 m a la edad de 8 años.

El comportamiento de las curvas, de ICA (Figura 3) e IMA (Figura 4) de la Muestra 2 es típico de las condiciones de una calidad de sitio mala: a) arrancan lentamente, b) les toma más tiempo alcanzar su pico máximo, c) nunca alcanzan los niveles de las otra y d) presentan una declinación abrupta en ambos incrementos, al alcanzar valores muy bajos a edades tempranas. Todo esto viene a ratificar las discrepancias que se encontraron al graficar las alturas dominantes versus la edad, para ambas muestras, representa- das en la Figura 2.

Construcción de las curvas de índice de sitio para la Muestra 1: modelo de Bailey y Clutter

Al correr el modelo de Bailey y Clutter con los datos de la Muestra 1, se obtuvieron los siguientes resultados: el valor de b es igual a -5,18103 (la pendiente de este modelo siempre va a ser negativa); por su parte, a es igual a 5,47261 y se encuentra dentro del rango de 2 a 7 indicado por ^{Alder (1980}). El valor del parámetro ĉ corres- ponde a 0,34326 y se encuentra entre 0,2 y 2 que son los valores propuestos por (Alder 1980) para la mayoría de las especies.

La Ecuación (1) se utilizó para ajustar los datos de la Muestra 1; a continuación se presentan los parámetros calculados durante la regresión:

Ln(Hdom) = 5.47261 -5.18103*(1/Edad)^0.34326

(13)

Al incorporar dichos parámetros a la ecua- ción, con la cual se construyeron las curvas de índice de sitio (Ecuación (6)), queda expresada de la siguiente manera: se obtuvo la siguiente función, con la cual se construyó un sistema de curvas para la Muestra 2.

H = e^5.47261  (IS/e^5.47261)^{9/Ei}^0.34326

(14)

H = Ei2/(0,465523 + (9/IS) -(0,465523/9)

-(0,01836639)Ei + (0,0183663Ei2))

(16)

Construcción de las curvas de índice de sitio para la Muestra 2: Modelo de Prodan

Al correr el modelo de Prodan con los datos de la Muestra 2, se obtuvieron los siguien- tes coeficientes, a = 0,465523; b = 0,538814; y c = 0,0183663. Ahora bien, al sustituir estos valo- res en la Ecuación (9), se obtiene la siguiente expresión,

Z = 0,465523 + 0,538814E + 0,0183663E2

(15)

Finalmente, al utilizar estos coeficientes y sustituyendo el valor de bi en la Ecuación (12),

Con las Ecuaciones (14) y (16) se logró generar 2 conjuntos de curvas polimórficas de índice de sitio, las cuales cubren todo el intervalo de edades y de alturas dominantes incorporadas en los datos de ambas muestras.

En el Cuadro 4 se presentan los esta- dísticos obtenidos como resultado de ajustar ambos modelos a los datos que conforman la Muestra 1 y la Muestra 2. En términos gene- rales, los resultados de los análisis estadísticos son bastante buenos, por lo que se indica que los modelos escogidos resultaron ser apropia- dos en ambos casos.

Curvas de índice de sitio diferenciadas según los patrones de crecimiento observados

Dadas las diferencias irreconciliables, encontradas en los patrones de crecimiento de la altura dominante representada en la Figura 2, se desarrollaron 2 juegos de curvas independientes con el fin de construir un sistema de clasificación de terrenos basado en los índices de sitio para V. guatemalensis. Para lograr esto, se construyeron 3 curvas de índice de sitio para las clases de creci- miento superior Muestra 1 y 3 curvas de índice de sitio para la clase de crecimiento inferior Muestra 2. Los índices de sitio, o índices de clasificación de sitios, fueron definidos utilizando la edad de 9 años como edad de referencia.

Se observó (Figura 5), que las localida- des de mejor crecimiento están representadas por los índices de sitio 19, 24 y 30; mientras que el grupo de parcelas que se encontraban creciendo en los peores sitios, fueron repre- sentadas por los índices de sitio 10, 14 y 18. Se puede ver en la figura anterior que a la edad de 9 años la curva de IS = 18, que correspon- de a la Muestra 2, se encuentra a muy poca distancia de la curva de IS = 19 que pertenece a la Muestra 1. Luego, estas curvas se irán separando conforme aumenta la edad, hasta el punto de que ambos patrones de crecimiento se vuelven incompatibles.

Fig. 5.  Curvas de índice de sitio diferenciadas para las 2 condiciones observadas en las plantaciones de V. guatemalensis, en las regiones Huetar Norte y Huetar Caribe de Costa Rica. 

Estas diferencias en el crecimiento pudie- ron deberse a factores edáficos que afectan el desarrollo de la especie en los sitios de mala cali- dad (por ejemplo un pH alto y una alta concen- tración de bases), o bien, en el caso de los sitios buenos, debido a factores edáficos que favorecen su crecimiento (por ejemplo una acidez elevada y una alta concentración de hierro y aluminio), tal y como lo establecieron (^{Mora et al. 2015}, Cama- cho et al. 2016, ^{Camacho et al. 2017}).

CONCLUSIONES

Al analizar el desempeño de la altura dominante en relación con la edad, para V. guatemalensis en las regiones Huetar Norte y Huetar Caribe de Costa Rica, se pudo detectar la existencia de 2 patrones de crecimiento cla- ramente diferenciados, lo cual permitió definir zonas de crecimiento en el gráfico altura dominante-edad.

Los sitios correspondientes a la Muestra 1 se pueden considerar de muy buena calidad. El crecimiento de las plantaciones en estos sitios es de muy bueno a excelente. Por otro lado, los sitios correspondientes a la Muestra 2 se pueden considerar de calidad de mala a muy mala para dicha especie.

La prueba de Levene indicó, con un 95% de confianza, que todos los datos provienen de una misma población, mientras que la prueba t confirmó que existen diferencias significativas en los promedios de las alturas dominantes de ambas muestras.

Los gráficos de ICAHdom e IMAHdom ver- sus Edad corroboran las diferencias encontradas al aplicar las pruebas basadas en estadística infe- rencial. Ambos gráficos demostraron la superio- ridad de las plantaciones de la Muestra 1 sobre los datos de la Muestra 2, debido a que las plantacio- nes de esta última crecen de muy mala manera.

Las curvas de índice de sitio se constru- yeron al respetar las discrepancias encontradas en la muestra y para eso se utilizó una edad base de 9 años. Las localidades de mejor crecimiento están representadas por los índices de sitio 19, 24 y 30; mientras que las de peor crecimiento son representadas por los índices de sitio 10, 14 y 18. Las 2 ecuaciones presentaron un ajuste bastante aceptable, la ecuación para las Clases I y II (Bailey y Clutter) presentó los siguientes estadís- ticos: r = -0,94; R2ajust =0,89; MSE = 0,0565; SEE

= 0,2378; y la ecuación de la Clase III (Prodan): r = 0,93; R2ajust = 0,87; MSE = 0,3886; SEE = 0,6234.

RECOMENDACIONES

Se recomienda establecer la especie sola- mente en los índices de sitio 19 a 30, ya que bajo estas condiciones los árboles muestran su mejor desarrollo.

Es conveniente extender el muestreo a otras áreas y otras regiones donde se encuentra creciendo la especie, de manera que se pueda diversificar la información y aumentar el tamaño de la muestra.

Con el sistema de clasificación de sitios propuesto en este trabajo, se recomienda hacer acopio de toda la información relacionada con los diferentes tipos de suelo sobre los cuales se encuentra creciendo la especie, para establecer correlaciones entre las características físicas y químicas de los mismos con el crecimiento en altura dominante de V. guatemalensis.