SciELO - Scientific Electronic Library Online

 
vol.40 número1Las empresas forestales comunitarias en la región de la Mariposa Monarca, México: un enfoque empresarialEvaluación preliminar de la influencia de la disponibilidad de larvas en el comportamiento de forrajeo del carpintero negro índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google

Compartir


Bosque (Valdivia)

versión On-line ISSN 0717-9200

Bosque (Valdivia) vol.40 no.1 Valdivia  2019

http://dx.doi.org/10.4067/S0717-92002019000100071 

Artículo

Modelos de efectos mixtos de altura-diámetro para Drimys winteri en el sur (41-43° S) de Chile

Mixed-effects height-diameter models for Drimys winteri in the south (41-43° S) of Chile

Christian Salas-Eljatib a   b   * 

Patricio Corvalán c  

Nicolás Pino a  

Pablo J Donoso d  

Daniel P Soto e  

a Universidad Mayor, Facultad de Ciencias, Centro de Modelación y Monitoreo de Ecosistemas, José Toribio Medina 29, Santiago, Chile, tel.: +56(2) 25189719

b Universidad de La Frontera, Departamento de Ciencias Forestales, Laboratorio de Biometría, Temuco, Chile

c Universidad de Chile, Facultad de Ciencias Forestales y de la Conservación de la Naturaleza, Santiago, Chile, pcorvalanvera@gmail.com

d Universidad Austral de Chile, Facultad de Ciencias Forestales y Recursos Naturales, Instituto de Bosques y Sociedad, Valdivia, Chile, pdonoso@uach.cl

e Universidad de Aysén, Departamento de Ciencias Naturales y Tecnología, Coyhaique, Chile, daniel.soto@uaysen.cl

RESUMEN:

Drimys winteri es una especie muy abundante en los bosques naturales de Chile, conformando bosques secundarios de alto rendimiento entre las latitudes 41° y 43° S. La relación altura-diámetro es clave para realizar estimaciones de volumen, biomasa, y para entender características autoecológicas de las especies forestales. Sin embargo, esta relación alométrica ha sido poco investigada para esta especie. El objetivo del presente estudio fue modelar la relación altura-diámetro de D. winteri mediante modelos no-lineales de efectos mixtos. Se empleó una amplia base muestral de árboles entre las provincias de Llanquihue y Chiloé en el sur de Chile. Se analizaron diferentes modelos, y posteriormente se evaluó a qué parámetro(s) agregar efectos aleatorios. Finalmente, se modelaron dichos efectos en función de variables agregadas de rodal (e.g., densidad, área basal). Se propone un modelo que: posee un comportamiento biológicamente consistente; sus parámetros son estimados de forma eficiente al considerar la estructura jerárquica de los datos; e incluye el efecto de la competencia (representado por la densidad y diámetro del árbol de área basal media) y la calidad del sitio (representado por el incremento medio anual en altura dominante) de los bosques donde estos árboles crecen. Esto implica que la alometría altura-diámetro de D. winteri varía dependiendo de la competencia y productividad del sitio donde crece la especie, dos variables fundamentales de considerar en la toma de decisiones silvícolas.

Palabras clave: alometría; máxima verosimilitud; ecología forestal; dinámica de bosques; tipo forestal siempreverde

SUMMARY:

Drimys winteri is an abundant tree species in the natural forests of Chile, growing as secondary forest of high ecological and economical importance between 41? and 43? S. The height-diameter relationship is key for estimating volume, biomass and for understanding ecological features of tree species. However, there is a lack of research about allometric relationships for this species. We aim at modeling the height-diameter relationship of D. winteri by nonlinear mixed-effects models. We used an extensive sample of trees between Llanquihue and Chiloé Provinces in southern Chile. We analyzed several models, and afterwards, assessed to which parameters allocate the random effects. Finally, we modeled these random effects as a function of stand-level variables (e.g., density, basal area). The proposed model has a biologically consistent behavior; its parameters are efficiently estimated because we considered the hierarchical structure of the data and took into account the effect of competition (represented by density and diameter of the mean basal area tree) and site quality (represented by the mean annual increment in dominant height) of the forests where the trees are growing. Our findings showed that the height-diameter allometry of D. winteri varies depending on competition and site productivity, two key variables for silvicultural management.

Key words: allometry; maximum likelihood; forest ecology; forest dynamics; evergreen

INTRODUCCIÓN

Drimys winteri J.R et G. Forster, es una especie arbórea que posee una gran plasticidad ecológica que le permite crecer en la mayor parte de los bosques naturales de Chile, entre los 30° y los 50° S, tanto en la cordillera de la Costa como en la de Los Andes, preferentemente ocupando terrenos húmedos (Larjavaara y Muller-Landau 2013). Drimys winteri es una especie pionera que coloniza e invade terrenos desnudos después de cortas e incendios especialmente en la región sur-austral de Chile formando importantes bosques secundarios (Corvalán et al. 1987a). Específicamente en la isla de Chiloé y la provincia de Llanquihue en el sur de Chile, D. winteri es una de las especies nativas de mayor interés productivo, formando bosques secundarios de un alto potencial de manejo silvícola (Corvalán et al. 1987b, Navarro et al. 1999).

La relación altura-diámetro de los árboles es muy importante en ecología y manejo de bosques. Esta relación es el resultado de diversos factores que afectan el crecimiento de un árbol, como por ejemplo, el sitio donde crece y la competencia a la cual se ve expuesto. El error y el alto costo de medición de las alturas de los árboles es un tema clásico de la mensura forestal que se resuelve utilizando modelos de regresión empleando el diámetro (d) como variable predictora. Los errores de medición son especialmente comunes en rodales de alta densidad como en los renovales de D. winteri (St-Onge et al. 2004, Larjavaara y Muller-Landau 2013), por lo que es conveniente desarrollar modelos de altura-diámetro (h-d) a base de mediciones obtenidas en una submuestra de árboles (Salas y Real 2013, Mehtätalo et al. 2015).

La mayoría de los modelos de altura-diámetro para D. winteri han sido ajustados en forma local y a base de muestras restringidas. Los modelos de altura-diámetro se han clasificado en: “locales” a los que solamente ocupan como variable predictora a la variable diámetro del árbol; y “generalizados” a los que incluyen además alguna(s) variable(s) de rodal (Temesgen et al. 2007, Mehtätalo et al. 2015). Todos los modelos de h-d para D. winteri reportados en la literatura son locales y se han desarrollado aplicando modelos estadísticos tradicionales, i.e., regresión lineal mediante el método de los mínimos cuadrados (Corvalán et al. 1987b, Navarro et al. 1999, Reyes et al. 2009). Los modelos locales de h-d solo son apropiados para predicciones en condiciones de bosque y sitio restringidas. Por el contrario, los modelos generalizados de h-d, permiten su uso en diferentes condiciones de bosque y sitios, lo cual en el caso de D. winteri, que es una especie que crece naturalmente en rodales densos y con una gran amplitud de calidad de sitios, los hace particularmente de interés dado su plasticidad ecológica. Sin embargo, a la fecha no existen estudios al respecto para esta especie.

El efecto de variables de rodal y/o sitio en la modelación de la relación altura-diámetro ha sido deficientemente evaluado desde un punto de vista estadístico. Cuando un modelo generalizado de h-d se ajusta mediante el método de los mínimos cuadrados (e.g., en Staudhammer y LeMay 2000, Schröder y Álvarez 2001, Trincado y Leal 2006), no se considera la estructura jerárquica de los datos (i.e., árboles anidados a unidades de muestreo), induciendo por lo tanto una correlación que debe ser tomada en cuenta al momento de estimar los parámetros de los modelos y sus respectivos estimadores de varianza. Esto trae problemas de pseudoreplicación (Hurlbert 1984), que no han sido siempre bien reconocidos en la literatura. Una alternativa estadísticamente correcta para representar el efecto de las variables de rodal (e.g., densidad, área basal) y del sitio (e.g., índice de sitio, o nivel de crecimiento medio) en un modelo h-d, es ajustarlo en un contexto de efectos mixtos (Mehtätalo 2004, Robinson y Wykoff 2004, Mehtätalo et al. 2015). El objetivo del presente estudio es modelar el efecto de variables de rodal y sitio en la alometría altura-diámetro de D. winteri, mediante un enfoque de efectos mixtos. Los modelos que incluyen estas variables pueden ser de gran apoyo para la toma de decisiones silviculturales tales como la selección de sitio y el control de la competencia.

MÉTODOS

Datos

Se establecieron 122 unidades de muestreo (i.e., parcelas) en bosques no-intervenidos de D. winteri en las provincias de Chiloé y Llanquihue en el sur de Chile (figura 1). Dichas unidades fueron establecidas entre los 41° y 43° S por Corvalán et al. (1987a) y Donoso et al. (2007). La superficie de las parcelas varía entre 100 y 400 m2. En cada unidad de muestreo se midieron variables dendrométricas a todos los árboles con un diámetro a la altura del pecho (d) mayor a 5 cm, tales como: d, especie, posición sociológica, y otras variables. Además, a una submuestra aleatoria de árboles por parcela, se les midió la altura total (h) con hipsómetro. El tamaño de la submuestra de altura por parcela fue variable, oscilando entre 3 y 43 árboles. Además, a una submuestra de tres o cinco árboles dominantes (que pertenecen al dosel superior) de cada parcela se les determinó la edad, mediante análisis fustal (Corvalán et al. 1987a) o extracción de un tarugo de incremento (Donoso et al. 2007), respectivamente. La dispersión entre la altura y el diámetro muestra un amplio rango de variación de la especie en el área de estudio, demostrando la plasticidad ecológica de la especie (figura 2).

Figura 1 Área de estudio y distribución de las unidades de muestreo (i.e., parcelas) en bosques de Drimys winteri en las provincias de Llanquihue y Chiloé, sur de Chile. Nótese que debido a la escala del mapa, existen parcelas que están traslapadas 

Figura 2 Dispersión entre la altura y el diámetro para árboles de Drimys winteri en las provincias de Llanquihue y Chiloé, en el sur de Chile 

Con la información medida en cada unidad muestral, se calcularon variables a nivel agregado o de estado de rodal, como la densidad (N, número de árboles), área basal (G), diámetro del árbol de área basal media ( dg ), según lo propuesto por van Laar y Akça (2007). La altura dominante ( Hd ) fue estimada mediante el algoritmo de García y Batho (2005), el cual toma en consideración la superficie de las parcelas, lo cual afecta su estimación. A base de las mediciones de los tarugos de incremento de los árboles dominantes se determinó la edad dominante de la unidad de muestreo ( Ed ), como el promedio de las edades de los árboles dominantes. Finalmente, se obtuvo el incremento medio anual en altura dominante ( imaHd ), que es el cuociente entre Hd y Ed , como una aproximación a la productividad del sitio (van Laar y Akça 2007). El amplio rango de las variables a nivel individual y de parcela (cuadro 1), así como también su gran variabilidad demuestra lo exhaustivo del muestreo en terreno.

Cuadro 1 Estadística descriptiva para las principales variables a nivel de árbol y de rodal. Las variables de árbol son: diámetro (d) y altura total (h), mientras que las de rodal son densidad (N), área basal (G), diámetro del árbol de área basal media ( dg ), edad dominante ( Ed ), e incremento medio anual en altura dominante ( imaHd ). Descriptive statistics of tree- and stand-level variables. Tree-level variables are: diameter (d) and total height (h); meanwhile stand-level variables are: tree density (N), basal area (G), diameter of the mean basal area tree ( dg ), dominant age ( Ed ), and mean annual increment of dominant height ( imaHd ). 

Modelación

Primero se seleccionó una ecuación base de altura-diámetro, y en segunda instancia, se empleó dicha ecuación para ajustarlo mediante modelos de efectos mixtos. Para seleccionar la ecuación base, se ajustaron una serie de modelos no-lineales, dentro de los cuales destacan el de Richards (1959) y el de Michaelis-Menten (Bates y Watts 1980, Mehtätalo et al. 2015) mediante el método de mínimos cuadrados no-lineales. Luego de una evaluación predictiva y de las características biológicas de estos modelos, i.e., que al menos uno, o la mayoría, de sus parámetros tenga una interpretación directa dimensional (e.g., la asíntota, ver anexo), se seleccionó el modelo de Michaelis-Menten, que se expresa como

[1]

donde: h 𝑖 y 𝑑 𝑖 es la altura total y diámetro a los 1,3 m sobre el fuste para la i-ésima observación, respectivamente; α y β son parámetros; y 𝜀 𝑖 es el i-ésimo error aleatorio del modelo con distribución Gaussiana con media 0 y varianza 𝜎 𝜀 2 .

Este modelo ha sido ampliamente empleado para estudios de cinética enzimática en reacciones químicas y en estudios del crecimiento de plantas (Kobe 1999). Una de las particularidades de la ecuación [1] es que posee una interpretación biológica de sus parámetros, donde α representa a la asíntota (i.e., máxima altura) y β al diámetro donde se alcanza la mitad de la asíntota. En el anexo se detalla una demostración matemática al respecto.

En una segunda etapa, y empleando el modelo de Michaelis-Menten como base se ajustaron diferentes modelos de efectos mixtos, en donde se agregaron efectos aleatorios a cada uno de los parámetros-coeficientes de ecuación [1]. Los modelos de efectos mixtos permiten representar la estructura jerárquica de los datos (Pinheiro y Bates 2000, Trincado et al. 2007, Salas et al. 2007), ya que las alturas individuales son medidas dentro de unidades de muestreo, por lo tanto, existen datos anidados. Se evaluaron diferentes variantes de modelos de efectos mixtos, al agregar efectos aleatorios a uno y/o a todos o diferentes combinaciones de estos. El ajuste se realizó mediante máxima verosimilitud empleando el paquete nlme (Pinheiro et al. 2017) de R (R Core Team 2017). Luego de comparar los estadísticos de bondad de ajuste entre las diferentes variantes de modelos de efectos mixtos (cuadro 2), tales como AIC, BIC y test de razón de verosimilitud (Pinheiro y Bates 2000), se seleccionó el siguiente modelo de efectos mixtos como la mejor variante

[2]

donde: ℎ 𝑖𝑗 y 𝑑 𝑖𝑗 es altura y diámetro para el i-ésimo árbol en la j-ésima parcela, respectivamente; α y β son parámetros; 𝑎 𝑗 y 𝑏 𝑗 son efectos aleatorios a nivel de parcela y 𝜀 𝑖𝑗 es el error aleatorio del modelo.

El modelo [2] se ajustó además considerando una varianza de los residuales diferente para cada clase de la variable incremento medio anual (i.e., cuatro clases definidas por los cuartiles de la variable, ima Hd ) y fuente de datos, tal como lo aplicaron Salas et al. (2008).

Finalmente, se modelaron los efectos aleatorios en función de variables de estado de rodal que representan la densidad de los bosques (e.g., N, G, d g ) y la calidad del sitio (e.g., ima Hd ), debido a su potencial efecto en la alometría altura-diámetro. Para tales efectos, se obtuvieron parámetros del modelo de altura-diámetro a nivel de parcela como sigue,

[3]

[4]

donde: 𝛼 y 𝛽 son parámetros estimados; 𝑎 𝑗 y 𝑏 𝑗 son efectos aleatorios predichos por el modelo de efectos mixtos para la j-ésima parcela.

Se evaluaron diferentes modelos para predecir 𝛼 𝑗 y 𝛽 𝑗 de la forma general,

[5]

[6]

donde: 𝑓( ) es un modelo lineal o no-lineal; 𝜃 y ∅ son vectores de parámetros a ser estimados; 𝑋 𝛼 y 𝑋 𝛽 son matrices de variables predictoras para el modelo [5] y [6], respectivamente; y 𝜀 𝛼 𝑗 y 𝜀 𝛽 𝑗 son vectores de errores aleatorios para el modelo de 𝛼 𝑗 y 𝛽 𝑗 , respectivamente.

Estos modelos (ecuaciones [5] y [6]) de efectos aleatorios se denominaron sub-modelos, y se ajustaron mediante mínimos cuadrados, evaluando diferentes potenciales variables predictoras que puedan representar el efecto de la competencia y de la calidad del sitio de los bosques donde se desarrollan los árboles. La selección de los mejores sub-modelos, contempló evaluar el efecto de su predicción en el comportamiento del modelo de altura-diámetro. Para esto se desarrollaron análisis gráficos de comportamiento bajo diferentes niveles de las variables predictoras. Para todos los modelos y sub-modelos ajustados, se calcularon estadísticos de capacidades predictivas como la raíz cuadrada media de las desviaciones, la diferencia agregada, y diferencia absoluta (Salas et al. 2010).

RESULTADOS

La dispersión de la relación altura-diámetro muestra la alta variabilidad de D. winteri creciendo en condiciones naturales, así como demuestra la plasticidad ecológica de la especie o la variabilidad genética (figura 2). La mejor variante de modelos de efectos mixtos fue la que considera agregar efectos aleatorios a los dos parámetros de la ecuación de Michaelis-Menten (cuadro 2), es decir un modelo completo (ecuación [2]). Tal como lo sugieren Pinheiro y Bates (2000), se evaluó la sobre-parametrización del modelo de efectos mixtos al analizar la relación entre los efectos aleatorios, no encontrándose evidencia de esta.

Cuadro 2 Estadísticos de bondad de ajuste basados en máxima verosimilitud para tres variantes del modelo de efectos mixtos de la ecuación altura-diámetro de Michaelis-Menten (ecuación [1]). La variante α es cuando los efectos aleatorios son agregados a dicho parámetro del modelo base, y así sucesivamente para las otras dos variantes. Nótese que la variante α y β es representada en la ecuación [2]. AIC es el criterio de información de Akaike y BIC es el criterio de información Bayesiano. Maximum likelihood based statistics of the height-diameter Michaelis-Menten (equation [1]), for three mixed-effects model variants. The α variant add random effects into that parameter, and so on for the other variants. Notice that the α and β variant is represented by equation [2]. AIC is the Akaike’s information criterion and BIC is the Bayesian information criterion. 

La estructura del modelo de efectos mixtos propuesto, implica un efecto de las características de los bosques donde crecen los árboles en la asíntota del modelo (α en ecuación [1]) y en el parámetro que determina el diámetro cuando se alcanza la mitad de la asíntota (β en ecuación [1]). Así también el modelo propuesto toma en consideración la estructura jerárquica de los datos, por lo tanto, estimando eficientemente las varianzas de los parámetros estimados del modelo, lo cual es clave para la correcta inferencia estadística.

Luego de analizados una gran variedad de sub-modelos para 𝛼 𝑗 y 𝛽 𝑗 se seleccionó una misma relación funcional para ambos parámetros, la cual fue un modelo lineal del siguiente tipo

[7]

[8]

donde: 𝑖𝑚𝑎 𝐻𝑑 𝑗 , 𝑑 𝑔 𝑗 y 𝑁 𝑗 son el incremento medio anual en altura dominante, diámetro del árbol de área basal media, y densidad para la j-ésima parcela, respectivamente; α0, ..., α4 y β0, ..., β4 son parámetros, y 𝜀 𝛼𝑗 y 𝜀 𝛽𝑗 son errores aleatorios de los sub-modelos de 𝛼 𝑗 y 𝛽 𝑗 , respectivamente.

El ajuste de estos sub-modelos fueron estadísticamente significativos (P < 0,05), y con errores medios de 2,38 (m) y 1,56 (cm) para 𝛼 𝑗 y 𝛽 𝑗 , respectivamente (cuadro 3).

Cuadro 3 Parámetros estimados del modelo de altura-diámetro de efectos mixtos (ecuación [2]) y sub-modelos para predecir α j (ecuación [7]) y β j (ecuación [8]). Estimated parameters of the mixed-effects height-diameter model (equation [2]) and sub-models α j (equation [3]) and β j (equation [4]). 

Estos sub-modelos implican que la competencia y la calidad del sitio donde crecen los árboles afectan en forma simultánea a la relación alométrica altura-diámetro. La competencia es representada, en el modelo propuesto, por la densidad y el diámetro medio de los bosques, mientras que la calidad del sitio, por el incremento medio anual en altura dominante. Según los análisis de comportamiento, mayores densidades producen relaciones altura-diámetro de mayor magnitud, es decir, para un mismo diámetro de referencia, un árbol puede alcanzar una altura mayor (figura 3).

Figura 3 Comportamiento del modelo propuesto de altura-diámetro para Drimys winteri según diferentes niveles de competencia (representado por la densidad, N) y productividad de sitio (representado por el incremento medio anual en altura dominante, imaHd). Para las tres figuras se empleó un diámetro del árbol de área basal media (dg) igual a 15 cm. La curva “poblacional” se refiere al valor esperado del modelo sin considerar los efectos de la competencia y sitio 

Los estadísticos de validación del modelo propuesto, es decir: predicción de α y β de la ecuación [1] al emplear los valores respectivos para las variables de rodal y sitio de la j-ésima parcela, y luego calculó del valor predicho para el i-ésimo árbol dentro de la j-ésima parcela, muestran una disminución en relación al ajuste del modelo base mediante mínimos cuadrados. En particular la raíz cuadrada media de las diferencias (RMSD) disminuye un 6 % (cuadro 4).

Cuadro 4 Estadísticos de capacidades predictivas para modelo base (ecuación [1]) y modelo de efectos mixtos (ecuación [2]), con efectos aleatorios predichos según los sub-modelos αj (ecuación [7]) y βj (ecuación [8]). RMSD representa la raíz cuadrada media de las diferencias, DA a la diferencia agregada (o error medio), y DAA a la diferencia media absoluta. Statistics of prediction capabilities for the base model (equation [1]) and the mixed-effects model (equation [2]), with predicted random-effects from submodels αj (equation [3]) and βj (equation [4]). RMSD is the root mean square differences; DA is the aggregated difference (or average error); and DAA the aggregated of the absolute value differences.  

DISCUSIÓN

Las alturas y diámetros reportadas en esta investigación, se encuentran dentro del rango esperado para el área de estudio, con alturas máximas de 25 m y diámetros inferiores a 100 cm (Donoso et al. 2006). Si bien los valores máximos de estas variables son menores a los reportados para otras especies del bosque nativo chileno, por ejemplo, para especies del género Nothofagus (Salas y Real 2013), son valores típicos para D. winteri en bosques de segundo crecimiento (Corvalán et al. 1987b, Navarro et al. 1999, Soto y Donoso 2006). La gran variabilidad de las alturas que D. winteri puede alcanzar para un mismo diámetro demuestra la plasticidad de la especie para crecer en diferentes condiciones de sitio y de competencia.

En general, para representar el efecto de variables de rodal y/o de sitio en modelos de altura-diámetro, se han ajustado modelos mediante el método de mínimos cuadrados (Staudhammer y LeMay 2000, Trincado y Leal 2006), lo que tiende a sobreestimar el efecto o la importancia estadística de las variables de rodal debido a la pseudoreplicación (Hurlbert 1984). En cambio, en el presente estudio, se analizó el efecto de variables de rodal mediante modelos de efectos mixtos, lo cual permite asignar eficientemente el efecto de dichas variables en la jerarquía correspondiente. Hubo una mejora en el ajuste empleando efectos mixtos, donde los estadísticos de bondad y de validación disminuyeron en comparación al ajuste por mínimos cuadrados. Así también, al definir diferentes escenarios y condiciones para las variables predictoras (figura 2), se logró cubrir el rango de variabilidad de la alometría altura-diámetro.

La relación altura-diámetro de D. winteri se ve influenciada por dos grandes factores: la competencia y la calidad de sitio. El efecto de la competencia ha sido representado mayormente a nivel agregado por el número de árboles (Parresol 1992, Zeide y Vanderschaaf 2002, Soares y Tomé 2002) y el área basal (Parresol 1992, Staudhammer y LeMay 2000), y se basa en el hecho que los árboles que crecen en rodales más densos tienden a desarrollarse más en altura que en diámetro, para así competir de mejor forma por la luz. Por lo tanto, esto altera la alometría altura-diámetro, y el modelo propuesto (que depende del número de árboles) así lo representa. El efecto del tamaño medio de los árboles en un bosque, representado por dg, tiene que analizarse en relación al diámetro del árbol a ser evaluado. Esto es, si el árbol tiene un diámetro menor al dg del bosque, representaría a un árbol que está sometido a mayor competencia que un árbol de área basal media. Por lo tanto, debería considerarse su efecto en el presente modelo, como un tipo de índice de competencia a nivel individual (van Laar y Akça 2007). En los estudios de Mehtätalo (2004), Trincado y Leal (2006) y Mehtätalo et al. (2015) se han empleado con éxito también el dg como variable de rodal en modelos de h-d. Considerando todo lo anterior, el modelo de h-d propuesto considera el efecto de la competencia a nivel agregado y a nivel individual, con variables de fácil obtención.

El segundo factor, la calidad del sitio, en modelos h-d ha sido representada casi exclusivamente por la altura dominante (Hd) del bosque, variable de rodal que se debe calcular (o estimar) a base de mediciones (o estimaciones) de alturas individuales que normalmente provienen de un modelo h-d, por lo tanto, su uso es circular en el contexto del presente estudio. Es por esto que en el presente estudio deliberadamente se excluyó emplear la variable Hd. Como indicador indirecto de la productividad de un sitio, se empleó al imaHd, que incluso puede ser categorizado en sitios buenos, regulares y malos, para su uso más práctico. Debido a la mayor productividad del sitio, producto de mejores condiciones edáficas, topográficas y climáticas, existe una mayor disponibilidad de recursos para el crecimiento de un árbol, lo cual se traduce en un mayor desarrollo en altura (i.e., efecto positivo). Es importante destacar que el efecto combinado de ambos factores (i.e., competencia y sitio), implica que la competencia produce árboles de mayor altura en sitios de mayor productividad (Pretzsch 2009).

Finalmente, los resultados acá presentados proveen modelos de altura para diferentes sitios cubiertos por bosques secundarios dominados por D. winteri, y adicionalmente consideran el efecto de la densidad sobre el crecimiento en altura dentro de sitios de la misma calidad. Este es el primer estudio en Chile que evidencia que la altura de los árboles dominantes es influida por la densidad, con una relación negativa como es de esperar, ya que a menor densidad los árboles pueden asignar menos carbohidratos al crecimiento en altura en favor del desarrollo de las copas y los fustes (Kramer y Koslowski 1960, Lambers et al. 2008). Evidencias de este tipo han sido reportadas en otras especies, como por ejemplo en Pseudotsuga menziesii (Mirbel) Franco, en Estados Unidos (Salas et al. 2008). Los modelos acá provistos en consecuencia constituyen una herramienta de gran valor para mejorar las estimaciones de crecimiento y productividad de los bosques secundarios dominados por D. winteri en Chile, los cuales constituyen una reserva de bosques de alto potencial para el manejo y para el desarrollo forestal especialmente entre los 41° y 43° S (Navarro et al. 1999, Navarro y Cabello 2018).

CONCLUSIONES

La modelación de la alometría altura-diámetro mediante un enfoque de modelos no-lineales de efectos mixtos, permitió representar satisfactoriamente su comportamiento bajo diferentes factores que la afectan. Estos factores son la competencia y la calidad del sitio, los cuales fueron representados por variables predictoras, cuya significancia estadística fue correctamente estimada, dado que el modelo toma en consideración la estructura jerárquica de los datos. Además, el modelo propuesto posee una interpretación ecológica coherente con el conocimiento actual del crecimiento de árboles y, ofrece mejores capacidades predictivas que un modelo de altura-diámetro local.

AGRADECIMIENTOS

Al proyecto de investigación FIBN 027/2016 y al proyecto FIA 007/84.

REFERENCIAS

Bates DM, DG Watts. 1980. Relative curvature measures of nonlinearity. Journal of the Royal Statistical Society. Series B (Methodological) 42(1):1-25. DOI: 10.1080/03610919408813178 [ Links ]

Corvalán P, L Araya, S Blanco, F Cox. 1987a. El Canelo: una alternativa de desarrollo para la X región. Volumen III Metodología. Santiago, Chile. Fondo de Investigación Agraria y Universidad de Chile. 145 p. [ Links ]

Corvalán P, L Araya, R Calquin, V Loewe, S Niebuhr. 1987b. El Canelo: una alternativa de desarrollo para la X región. Volumen IV Resultados. Santiago, Chile. Fondo de Investigación Agraria y Universidad de Chile. 185 p. [ Links ]

Donoso C, B Escobar, P Donoso, F Utreras. 2006. Drimys winteri J.R et G. Forster. In Donoso C ed. Las especies arbóreas de los bosques templados de Chile y Argentina. Autoecología. Valdivia, Chile. Marisa Cuneo Ediciones. p. 220-232. [ Links ]

Donoso PJ, DP Soto, RA Bertín. 2007. Size-density relationships in Drimys winteri secondary forests of the Chiloe Island, Chile: Effects of physiography and species composition. Forest Ecology and Management 239(1-3):120-127. DOI: 10.1016/j.foreco.2006.11.015 [ Links ]

García O, Batho A. 2005. Top height estimation in lodgepole pine sample plots. Western Journal of Applied Forestry 20(1): 64-68. DOI: 10.1093/wjaf/20.1.64 [ Links ]

Hurlbert SH. 1984. Pseudoreplication and the design of ecological field experiments. Ecological Monographs 54(2):187-211. DOI: 10.2307/1942661 [ Links ]

Kobe RK. 1999. Light gradient partitioning among tropical tree species through differential seedling mortality and growth. Ecology 80(1):187-201. DOI: 10.2307/176989 [ Links ]

Kramer PJ, TT Kozlowski. 1960. Physiology of trees. New York, USA. McGraw-Hill. 642 p. [ Links ]

Lambers H, FS Chapin III, TL Pons. 2008. Plant Physiological Ecology. New York, USA. Springer. 605 p. [ Links ]

Larjavaara M, H Muller-Landau. 2013. Measuring tree height: a quantitative comparison of two common field methods in a moist tropical forest. Methods in Ecology and Evolution 4(9):793-801. DOI: 10.1111/2041-210X.12071 [ Links ]

Mehtätalo L. 2004. A longitudinal height-diameter model for Norway spruce in Finland. Canadian Journal of Forest Research 34:131-140. DOI: 10.1139/x03-207 [ Links ]

Mehtätalo L, S de Miguel, TG Gregoire. 2015. Modeling height-diameter curves for prediction. Canadian Journal of Forest Research 45:826-837. DOI: 10.1139/cjfr-2015-0054 [ Links ]

Navarro C, C Donoso, V Sandoval. 1999. Los renovales de canelo. In Donoso C, A Lara eds. Silvicultura de los bosques nativos de Chile. Santiago, Chile. Editorial Universitaria. p. 341-377. [ Links ]

Navarro C, J Cabello. 2018. Caracterización espacial de bosques de segundo crecimiento de Drimys winteri en Chile para la determinación exploratoria de áreas de gestión silvicultural. In Donoso P, A Promis, D Soto eds. Silvicultura en Bosques Nativos. Experiencias en silvicultura y restauración en Chile, Argentina y el oeste de los Estados Unidos. p. 135-156. [ Links ]

Parresol BR. 1992. Baldcypress height-diameter equations and their prediction confidence intervals. Canadian Journal of Forest Research 22(9):1429-1434. DOI: 10.1139/x92-191 [ Links ]

Pinheiro J, D Bates, S DebRoy, D Sarkar, the R Core Team. 2017. nlme: Linear and nonlinear mixed effects models. R package version 3.1-131. [ Links ]

Pinheiro JC, DM Bates. 2000. Mixed-Effects Models in S and Splus. New York, USA. Springer-Verlag. 528 p. [ Links ]

Pretzsch H. 2009. Forest Dynamics, Growth, and Yield. Berlin, Germany. Springer-Verlag. 664 p. [ Links ]

R Core Team. 2017. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org. [ Links ]

Reyes R, P Donoso, C Donoso, C Navarro. 2009. Crecimiento de renovales de Drimys winteri después de 16 años de aplicados distintos tratamientos de raleo en las cordilleras de Los Andes y de la Costa en Chile. Bosque 30(3):117-126. DOI: 10.4067/s0717-92002009000300001 [ Links ]

Richards FJ. 1959. A flexible growth function for empirical use. Journal of Experimental Botany 10(29):290-300. DOI: 10.1093/jxb/10.2.290 [ Links ]

Robinson AP, WR Wykoff. 2004. Imputing missing height measures using a mixed-effects modeling strategy. Canadian Journal of Forest Research 34(12):2492-2500. DOI: 10.1139/x04-137 [ Links ]

Salas C, L Ene, TG Gregoire, E Næsset, T Gobakken. 2010. Modelling tree diameter from airborne laser scanning derived variables: A comparison of spatial statistical models. Remote Sensing of Environment 114(6):1277-1285. DOI: 10.1016/j.rse.2010.01.020 [ Links ]

Salas C, P Real. 2013. Biometría de los bosques naturales de Chile: estado del arte. In Donoso P, A Promis eds. Silvicultura en los bosques nativos: avances en la investigación en Chile, Argentina y Nueva Zelanda. Valdivia, Chile. Marisa Cuneo Ediciones. p. 109-151. [ Links ]

Salas C, AR Stage, AP Robinson. 2008. Modeling effects of overstory density and competing vegetation on tree height growth. Forest Science 54(1):107-122. DOI: 10.1093/forestscience/54.1.107 [ Links ]

Soto DP, PJ Donoso. 2006. Patrones de regeneración en renovales de Drimys winteri en el centro-norte de la Isla de Chiloé: cambios de acuerdo al tamaño y la densidad relativa. Bosque 27(3): 241-249. DOI: 10.4067/S0717-92002006000300003 [ Links ]

Schröder J, JG Álvarez. 2001. Comparing the performance of generalized diameter-height equations for maritime pine in Northwestern Spain. Forstwissenschaftliches Centralblatt 120:18-23. DOI: 10.1007/BF02796077 [ Links ]

Soares P, M Tomé. 2002. Height-diameter equation for first rotation eucalypt plantations in Portugal. Forest Ecology and Management 166:99 -109. DOI: 10.1016/S0378-1127(01)00674-0 [ Links ]

St-Onge B, J Jumelet, M Cobello, C Véga. 2004. Measuring individual tree height using a combination of stereophotogrammetry and lidar. Canadian Journal of Forest Research 34(10):2122-2130. DOI: 10.1139/x04-093 [ Links ]

Staudhammer C, VM LeMay. 2000. Height prediction equations using diameter and stand density measures. The Forestry Chronicle 76(2):303-309. DOI: 10.5558/tfc76303-2 [ Links ]

Temesgen H, ME Goerndt, GP Johnson, DM Adams, RA Monserud. 2007. Forest measurement and biometrics in forest management: Status and future needs of the Pacific Northwest USA. Journal of Forestry 105(5):233-238. [ Links ]

Trincado G, C Leal. 2006. Ecuaciones locales y generalizadas de altura-diámetro para pino radiata (Pinus radiata). Bosque 27(1):23-34. DOI: 10.4067/S0717-92002006000100003 [ Links ]

Trincado G, CL VanderSchaaf, HE Burkhart. 2007. Regional mixed-effects height-diameter models for loblolly pine (Pinus taeda L.) plantations. European Journal of Forest Research 126:253-262. DOI: 10.1007/s10342-006-0141-7 [ Links ]

van Laar A, A Akça. 2007. Forest Mensuration. Dordrecht, Netherlands. Springer. 383 p. [ Links ]

Zeide B, C Vanderschaaf. 2002. The effect of density on the height-diameter relationship. In Outcalt K ed. Proceedings of the eleventh biennial southern silvicultural research conference. USDA For. Ser. Gen. Tech. Rep. SRS-48. p. 463-466. [ Links ]

ANEXO

.

Demostración matemática del significado de los parámetros del modelo de Michaelis-Menten

.

Bajo la hipótesis de que a mayor diámetro, mayor es la altura del árbol, se analiza el modelo mediante el límite como sigue

[9]

[10]

[11]

[12]

Por lo tanto, α es la asíntota del modelo, es decir, el máximo valor posible de h

.

Sí se re-escribe el modelo de Michaelis-Menten [9] como

[13]

y si en [13] 𝑑=𝛽, se tiene que

[14]

[15]

Por lo tanto, β es el diámetro donde la altura del modelo es la mitad de la asíntota

.

Recibido: 28 de Junio de 2018; Aprobado: 24 de Enero de 2019

*Autor de correspondencia: cseljatib@gmail.com

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons