Introducción
En el manejo forestal se requiere el uso y diseño de herramientas cuantitativas para la toma de decisiones, los modelos de crecimiento y rendimiento forestal son los instrumentos más utilizados para lograr tal fin. Estos modelos tienen como función determinar la evolución en el tiempo de una o varias variables dendrométricas (del árbol individual o de rodal) que se definen en el sistema a estudiar. Los modelos a nivel de rodal son los más adecuados para la planeación del manejo de rodales coetáneos. Sin embargo, en algunos casos se requiere información más detallada para la toma de decisiones. Los modelos de distribución diamétrica permiten realizar una desagregación del volumen total para conocer la distribución de productos maderables, en tanto que, los modelos de árbol individual permiten estimar variables de crecimiento a ese nivel, así como conocer el número de árboles de cada clase diamétrica (Quevedo, Moret y Jerez, 2003; Gómez-García, Crecente-Campo, Pérez-Rodríguez y Diéguez-Aranda, 2013; Hirigoyen y Rachid, 2014)
La estimación de las distribuciones diamétricas se realiza a través de funciones de densidad de probabilidades (fdp) (Hafley y Schreuder, 1977), siendo las distribuciones Beta, Johnson SB, Weibull, lognormal, gamma y normal las comúnmente usadas para describir distribuciones diamétricas de rodales en bosques coetáneos e incoetáneos (Bailey y Dell, 1973; Hafley y Schreuder, 1977; Maltamo, Puumalainen y Päivinen, 1995).
La función Weibull derivada por Waloddi Weibull (1939) en estudios de resistencia de materiales, ha sido la más utilizada en la ciencia forestal por su relativa facilidad de aplicación, y fue introducida por Bailey y Dell (1973) para modelar distribuciones diamétricas de Pinus taeda, P. echinata, P. banksiana y Pseudotsuga menziesii. Entre las ventajas de la fdp Weibull están su manejo matemático sencillo y la flexibilidad para adoptar diferentes formas, desde la “J” invertida hasta distribuciones en forma de campana, con diferentes grados de sesgo; además, se puede integrar analíticamente para generar una función de distribución acumulativa de forma cerrada (Bailey y Dell, 1973; Maldonado y Navár, 2002) para determinar el número de árboles entre una clase diamétrica inferior y una superior. Esta función puede caracterizarse de manera dinámica a partir de la predicción de los parámetros, esto es, después de estimar los parámetros de localización, escala y forma (a, b y c, respectivamente) con alguno de los métodos diseñados para ello (Torres-Rojo, Magaña-Torres y Acosta-Mireles, 2000), se pueden relacionar con variables del rodal usando modelos lineales o no lineales. Con este método, se han obtenido resultados de ajuste estadísticamente adecuados que además permiten obtener directamente los parámetros de las distribuciones futuras proyectadas y facilita la planeación de los regímenes de cortas intermedias (Reynolds, Burk y Huang, 1988; García, Cañadas y Montero, 2002). Otra técnica para estimar los parámetros de la función Weibull es la mezcla entre métodos de estimación de parámetros y percentiles; como el método de recuperación de parámetros que consiste en recobrar los parámetros de la distribución a partir de sus momentos no-centrales y momentos centrales, al tener como base un conjunto de percentiles los cuales son fáciles de modelar de manera dinámica (Hyink y Moser, 1983).
En México se han generado modelos basados en la fdp Weibull para estimar distribuciones diamétricas en rodales coetáneos e incoetáneos y plantaciones de coníferas, entre las que destacan Pinus caribaea var. hondurensis Barr y Golf., P. durangensis Martínez, P. cooperi Blanco, P. engelmannii Carr., P. arizonica Engl. y P. patula (Montero y Fierros, 2000; Maldonado y Návar, 2002; Magaña, Torres, Rodríguez, Aguirre y Fierros, 2008; Santiago-García et al., 2014; Quiñonez, De los Santos, Cruz, Velázquez y Ramírez, 2015).
Objetivos
Predecir el rendimiento maderable por clase diamétrica para Pinus patula Schiede ex Schlechtendal & Chamisso con el uso de la función de densidad de probabilidades Weibull y las técnicas de predicción de parámetros y predicción de percentiles.
Materiales y métodos
Área de estudio
El estudio se llevó a cabo en rodales coetáneos de Pinus patula en bosques del predio de Ixtlán de Juárez, Oaxaca, México (Fig. 1), entre las coordenadas 17°18’16” y 17°30’00” latitud norte y 96°31’38” y 96°22’00” longitud oeste, con una altitud media de 2780 m snm. La región se localiza en la provincia fisiográfica denominada Sistema Montañoso del Norte de Oaxaca. Dicha provincia abarca la mitad septentrional de Oaxaca y áreas adyacentes de Puebla y Veracruz. Los tipos de climas predominantes en la zona son templado subhúmedo y templado subhúmedo con lluvias en verano, con temperatura media anual que va de 14 ºC a 18 °C. Los grupos de suelo más comunes en el área corresponden a acrisol, luvisol y cambisol. El principal tipo de vegetación presente corresponde a bosques de pino-encino (Servicios técnicos forestales de Ixtlán de Juárez [STF], 2015).
Datos dasométricos
Para el análisis se utilizó una base de datos de 132 distribuciones diamétricas, provenientes de dos mediciones realizadas en 66 sitios permanentes de investigación silvícola de forma cuadrada (400 m2) en rodales coetáneos de Pinus patula, los sitios se establecieron en 2015 cubriendo diferentes intervalos de densidad, calidad de sitio (20 m - 38 m) y edad (de 5 años a 80 años), vueltos a medir en 2016. Dentro de cada sitio se midió el diámetro normal (Dn) de todos los árboles, la altura (At) de al menos ocho árboles por sitio, de los cuales, cuatro se identificaron como dominantes de acuerdo con la definición de altura dominante que corresponde a los 100 árboles de mayor altura y diámetro por hectárea (Assmann, 1970; Alder, 1980).
Las variables de estado del rodal fueron: altura total promedio (At, m), altura dominante (Hd, m), área basal (Ab, m2 ha-1), número de árboles (Na ha-1), diámetro promedio (Dm, cm), diámetro medio cuadrático (Dq, cm), diámetro mínimo (Dmin, cm) y volumen (V, m3 ha-1). El diámetro medio cuadrático corresponde al diámetro del árbol de área basal media:
La altura total para los individuos no medidos fue obtenida con el modelo altura-diámetro ajustado por López-Villegas et al. (2017):
El volumen del fuste total con corteza para cada árbol (v, m3) se calculó con la ecuación de doble entrada ajustada por Rodríguez-Justino (2017):
En la tabla 1 se presentan las estadísticas descriptivas de los datos empleados en el ajuste de los sistemas de predicción implícita.
Variable | Observaciones | Media | Desviación estándar | Mínimo | Máximo |
---|---|---|---|---|---|
At | 9.57 | 4.87 | 1.88 | 36.74 | |
Dn | 8015 | 11.88 | 6.96 | 2.10 | 58.30 |
E † | 14.81 | 11.54 | 5.00 | 80.00 | |
Hd | 15.39 | 6.27 | 4.55 | 35.41 | |
Ab | 22.59 | 13.74 | 1.35 | 68.26 | |
Na | 1517.99 | 1011.85 | 300.00 | 6050.00 | |
Dm | 132 | 13.93 | 6.49 | 5.80 | 39.84 |
Dq | 14.93 | 7.16 | 5.87 | 43.52 | |
Dmin | 5.82 | 2.77 | 2.10 | 17.00 | |
V | 180.66 | 165.86 | 3.93 | 681.11 | |
E § | 18.72 | 16.87 | 5.00 | 75.25 |
At: altura total del árbol individual (m); Dn: diámetro normal con corteza (cm); E: edad (años); Hd: altura dominante (m); Ab: área basal (m2 ha-1); Na: número de árboles por hectárea (Na ha-1); Dm: diámetro promedio (cm); Dq: diámetro medio cuadrático (cm); Dmin: diámetro mínimo (cm); y V: volumen (m3 ha-1).
† Edad del árbol individual.
§ Edad promedio de la distribución diamétrica
La calidad del sitio forestal (índice de sitio, IS) y su correspondiente altura dominante se estimaron mediante el modelo polimórfico de Levakovic II, a una edad base de 40 años (Pérez-López et al., 2017; Santiago-García et al., 2017):
Función de distribución Weibull
La distribución Weibull triparamétrica está definida por la función de densidad probabilística siguiente (Bailey y Dell, 1973; Clutter, Forston, Pienaar, Brister y Bailey, 1983; Cao, 2004):
donde f(x) es la probabilidad asociada con cada posible valor de la variable aleatoria x (diámetro normal); a, b y c son los parámetros de localización, escala y forma, respectivamente.
La distribución por categoría diamétrica se determinó con la distribución acumulada de forma cerrada (Bailey y Dell, 1973; Clutter, Forston, Pienaar, Brister y Bailey, 1983):
La porción de la población comprendida dentro de un intervalo determinado está dada por la ecuación (Clutter, Forston, Pienaar, Brister y Bailey, 1983; Rennolls, Geary y Rollinson, 1985):
donde P es la proporción de árboles en la categoría diamétrica correspondiente, L y U son los límites inferior y superior de la categoría diamétrica, respectivamente, X es la categoría diamétrica y exp indica la función exponencial; el resto se definió previamente.
El producto de la proporción de árboles (P) en cada categoría diamétrica (X) y la densidad del rodal determinan el número de árboles por clase diamétrica, es decir, cuántos árboles de cierto diámetro se tendrán en una edad determinada (Quiñonez, De los Santos, Cruz, Velázquez y Ramírez, 2015).
Predicción de parámetros de la fdp Weibull
Los parámetros de la fdp Weibull se estimaron de acuerdo con la distribución diamétrica observada en cada sitio de muestreo; estos valores se obtuvieron mediante el procedimiento CAPABILITY de SAS/ETS® 9.3, el cual utiliza métodos basados en máxima verosimilitud para la estimación de parámetros (SAS Institute Inc., 2011).
Una vez estimados los parámetros fueron relacionados con los atributos dasométricos del rodal: área basal (Ab), diámetro mínimo (Dmin), diámetro promedio (Dm), diámetro medio cuadrático (Dq), número de árboles (Na), edad (E), altura dominante (Hd) e índice de sitio (IS) por medio de regresión lineal, con el procedimiento STEPWISE (selección por pasos) de SAS (SAS Institute Inc., 2011) (Maldonado y Návar, 2002). En este caso se consideraron dos ecuaciones para cada parámetro (Tabla 2).
Parámetro | Ecuación | Clave |
---|---|---|
Localización (a) |
|
a.E1 |
|
a.E2 | |
Escala (b) |
|
b.E1 |
|
b.E2 | |
Forma (c) |
|
c.E1 |
|
c.E2 |
Dmin: diámetro mínimo (cm); Dm: diámetro promedio (cm); Dq: diámetro medio cuadrático; IS: índice de sitio; Ca: coeficiente de asimetría de la distribución diamétrica; ln: logaritmo natural y βi: parámetros a estimar.
El coeficiente de asimetría (Ca) se calculó con la expresión:
donde:
n: número de observaciones
s es la desviación estándar
Predicción de percentiles
Para este método se requieren al menos dos percentiles de la distribución diamétrica del rodal; los cuales se ajustaron como funciones del diámetro medio cuadrático del rodal (Dq), en este estudio los percentiles 55 y 85 resultaron adecuados. Así, las ecuaciones de predicción de percentiles tomaron la forma siguiente:
donde:
p55: percentil 55 de la distribución diamétrica del rodal
p85: percentil 85 de la distribución diamétrica del rodal
Dq: diámetro medio cuadrático del rodal
α y β: parámetros por estimar
Los parámetros de la distribución Weibull fueron recuperados por el método de momentos, considerando el criterio de Pienaar y Rheney (1993). Las ecuaciones utilizadas se muestran a continuación:
parámetro de localización (a):
si Dmin < 5cm, entonces
A partir del parámetro de localización fue posible estimar el parámetro de forma (c) y de escala (b) respectivamente:
donde:
Dmin: diámetro mínimo del rodal (cm)
Estimación del número de árboles
Para la predicción del rendimiento maderable con distribuciones diamétricas fue necesario realizar la estimación del número de árboles sobrevivientes por hectárea a través del tiempo, esta variable fue estimada con la ecuación de mortalidad siguiente (Santiago-García et al., 2017):
donde Na 2 es el número de árboles vivos en la edad de proyección E 2 , Na 1 corresponde al número de árboles en la edad inicial E 1 .
Modelo de proyección del diámetro mínimo y diámetro medio cuadrático
El diámetro mínimo del rodal (percentil cero “p0”) es una variable de estado indispensable para realizar estimaciones de distribuciones diamétricas con el sistema basado en la fdp Weibull, debido a la estrecha relación que presenta con el parámetro de localización (a) (Santiago-García et al., 2014; Quiñonez, De los Santos, Cruz, Velázquez y Ramírez, 2015). Para la proyección del diámetro mínimo se utilizó la expresión anamórfica basada en el modelo de Chapman Richards (Amaro, Reed, Tomé y Temido, 1998):
donde:
Dmin 2: diámetro mínimo del rodal en E 2
Dmin 1: diámetro mínimo del rodal en E 1
E 1 : edad inicial del rodal (años)
E 2 : edad de proyección del rodal (años)
β i : parámetros por estimar
La proyección del diámetro medio cuadrático (Dq) del rodal se realizó mediante el ajuste del modelo compatible propuesto por Santiago-García et al. (2014) para predecir y proyectar el crecimiento de esta variable en el tiempo:
donde:
Dq 1 : diámetro medio cuadrático del rodal en la edad inicial E 1
Dq 2 : diámetro medio cuadrático del rodal a la edad de proyección E 2
Na 1 : número de árboles por hectárea en E 1
Na 2 : número de árboles por hectárea en E 2
Hd 1 : altura dominante en E 1
Hd 2 : altura dominante en E 2
β i : parámetros por estimar
Ajuste de los sistemas e indicadores estadísticos
A través del método de mínimos cuadrados ordinarios, se ajustaron los modelos de predicción de parámetros, en tanto que, el modelo de proyección del diámetro mínimo y las ecuaciones de predicción de los percentiles 55 y 85 de la distribución diamétrica, se ajustaron a través de métodos iterativos adecuados para modelos de regresión no lineal, en este caso, se empleó el método de Gauss-Newton. Por su parte, el modelo compatible de crecimiento para diámetro medio cuadrático fue ajustado a través de regresión aparentemente no relacionada (SUR, por sus siglas en inglés), con el paquete SAS/ETS®, mediante el procedimiento MODEL (SAS Institute Inc., 2011).
La técnica de estimación de parámetros SUR produce estimadores consistentes y eficientes en presencia de correlaciones contemporáneas en un sistema de ecuaciones, además, permite compatibilidad total, de forma que los parámetros comunes de las ecuaciones toman los mismos valores mientras se cumpla con el criterio de minimización de cuadrados de los residuos (Borders, Souter, Bailey y Ware, 1987; Borders y Patterson, 1990; Galán, De los Santos y Valdez, 2008).
La bondad de ajuste de los modelos estudiados se midió a través del análisis
numérico con la obtención de la suma de cuadrados del error
(SCE), la raíz del error cuadrático medio
(RECM), el coeficiente de determinación (R
2
), el coeficiente de determinación ajustado por el número de parámetros
(R
2
adj
), el sesgo promedio absoluto (
Bondad de ajuste de la fdp Weibull
Para verificar la hipótesis de que las distribuciones diamétricas de las parcelas siguen una distribución Weibull, se utilizó el estadístico de Kolmogorov-Smirnov (KS), este se basa en las diferencias absolutas entre las distribuciones estimadas y observadas, para esto existen valores tabulados que permiten decidir si la diferencia máxima entre las distribuciones es significativa (Sokal y Rohlf, 1979). El cálculo de la diferencia máxima se realizó a través de la expresión siguiente:
donde:
D max : diferencia máxima entre las distribuciones acumuladas
F(x): distribución acumulada teórica
S(x): distribución acumulada observada
Esto implica que dicha prueba solo es válida si los parámetros de la fdp Weibull son conocidos para cada distribución diamétrica (Torres-Rojo, Acosta y Magaña, 1992). Para esta prueba se utilizaron distintos niveles de significancia (α = 0.05, α = 0.10 y α = 0.20) para contrastar el número de parcelas no ajustadas a la fdp Weibull de acuerdo con cada criterio.
Resultados
Predicción de parámetros de la fdp Weibull
En la tabla 3 se presentan los resultados obtenidos de la prueba KS con diferentes niveles de significancia. De esta manera, se logró conocer la bondad de ajuste de la fdp Weibull a las estructuras diamétricas de los rodales de Pinus patula estudiados.
α | Número de distribuciones rechazadas |
Proporción* (%) |
---|---|---|
0.05 | 8 | 6.06 |
0.10 | 13 | 9.85 |
0.20 | 21 | 15.91 |
* Proporción calculada a partir de 132 distribuciones diamétricas
Pece, G. de Benítez y J. de Galíndez (2000) indican que el nivel de α = 0.20 es el más exigente, porque reduce las desviaciones mínimas permitidas para el no rechazo de la concordancia. Con este nivel de significancia se trabajaron las 111 distribuciones diamétricas restantes y se procedió al ajuste de las ecuaciones de predicción de los parámetros de la fdp Weibull en relación con las variables de estado del rodal.
El ajuste de las ecuaciones probadas fue estadísticamente adecuado (Tabla 4 y Tabla 5), al generar, en la mayoría de los casos, coeficientes de
determinación (R
2
y R
2
adj
) altos, suma de cuadrados del error (SCE), raíz del
error cuadrático medio (RECM), sesgo promedio absoluto (
Ecuación | SCE | RECM | R 2 | R 2 adj | AIC | |
---|---|---|---|---|---|---|
a.E1 | 112.279 | 1.010 | 0.859 | 0.859 | 0.041 | 3.271 |
a.E2 | 36.248 | 0.582 | 0.955 | 0.953 | 3.7E-15 | -116.226 |
b.E1 | 217.900 | 1.427 | 0.916 | 0.914 | 3.5E-15 | 82.870 |
b.E2 | 184.300 | 1.313 | 0.929 | 0.927 | 2.8E-15 | 64.281 |
c.E1 | 41.568 | 0.620 | 0.235 | 0.220 | 0.023 | -93.228 |
c.E2 | 36.358 | 0.583 | 0.330 | 0.312 | 0.024 | -115.889 |
SCE: suma de cuadrados del error; RECM: raíz del error cuadrático
medio; R2: coeficiente de determinación;
R2-adj: coeficiente de determinación ajustado por
el número de parámetros;
Ec | Estimación |
Error
estándar |
P > |t| | Ec | Estimación |
Error
estándar |
P > |t| | |||
---|---|---|---|---|---|---|---|---|---|---|
β 0 | 1.0755 | 0.137 | <0.0001 | |||||||
a.E1 | β 0 | 0.8817 | 0.015 | <0.0001 | a.E2 | β 1 | -0.7646 | 0.022 | <0.0001 | |
β 2 | 1.7209 | 0.065 | <0.0001 | |||||||
β 3 | -0.8566 | 0.058 | <0.0001 | |||||||
β 0 | -4.1956 | 0.514 | <0.0001 | β 0 | -6.0631 | 0.486 | <0.0001 | |||
b.E1 | β 1 | 2.9085 | 0.205 | <0.0001 | b.E2 | β 1 | 2.6496 | 0.183 | <0.0001 | |
β 2 | -1.4985 | 0.136 | <0.0001 | β 2 | -0.0360 | 0.003 | <0.0001 | |||
β 3 | 0.8536 | 0.028 | <0.0001 | β 3 | 0.8074 | 0.023 | <0.0001 | |||
β 0 | -0.3400 | 0.062 | <0.0001 | β 0 | -0.2908 | 0.060 | <0.0001 | |||
c.E1 | β 1 | 0.3886 | 0.069 | <0.0001 | c.E2 | β 1 | 0.3272 | 0.067 | <0.0001 | |
β 2 | 0.0454 | 0.004 | <0.0001 | β 2 | 0.0559 | 0.004 | <0.0001 | |||
β 3 | -0.3277 | 0.084 | 0.0002 |
Desde el punto de vista estadístico, las ecuaciones a.E2, b.E2 y c.E2 predicen satisfactoriamente los parámetros de localización, escala y forma de la fdp Weibull. Las ecuaciones ajustadas en conjunto con las funciones de altura y volumen hicieron posible predecir el rendimiento maderable en términos implícitos, que consiste en determinar las frecuencias de árboles por categoría diamétrica (Fig. 2 y Fig. 3), es importante resaltar que, para este método, el diámetro medio cuadrático se derivó de la proyección del área basal y número de árboles por hectárea del sistema explícito para P. patula (Santiago-García et al., 2017).
Predicción de percentiles
Las ecuaciones ajustadas para proyectar el diámetro mínimo, diámetro medio
cuadrático y las ecuaciones de predicción de los percentiles 55 y 85 de la
distribución diamétrica, generaron indicadores de bondad de ajuste
satisfactorios (Tabla 6 y Tabla 7), es decir, coeficientes de
determinación (R
2
y R
2
adj
) altos, suma de cuadrados del error (SCE), raíz del error cuadrático
medio (RECM), sesgo promedio absoluto (
Ecuación | SCE | RECM | R 2 | R 2 adj |
|
AIC |
---|---|---|---|---|---|---|
p55 | 690.100 | 2.304 | 0.908 | 0.908 | 0.090 | 222.333 |
p85 | 552.800 | 2.062 | 0.959 | 0.959 | 0.034 | 193.050 |
Dmin 2 | 7.141 | 0.334 | 0.986 | 0.986 | 0.003 | -142.767 |
Dq 1 | 258.400 | 2.009 | 0.925 | 0.924 | 0.043 | 96.080 |
Dq 2 | 29.382 | 0.672 | 0.991 | 0.991 | 0.059 | -47.412 |
SCE: suma de cuadrados del error; RECM: raíz del error cuadrático
medio; R2: coeficiente de determinación;
R2
adj: coeficiente de determinación ajustado por el
número de parámetros;
Ecuación | Parámetro | Estimación | Error estándar | P > |t| |
---|---|---|---|---|
p55 | α | 0.7532 | 0.065 | <0.0001 |
β | 1.0813 | 0.028 | <0.0001 | |
p85 | α | 1.1051 | 0.063 | <0.0001 |
β | 1.0563 | 0.018 | <0.0001 | |
Dmin 2 | β 1 | 0.1911 | 0.050 | 0.0003 |
β 2 | 1.9149 | 0.453 | <0.0001 | |
β 0 | 2.2347 | 0.073 | <0.0001 | |
Dq 1 y Dq 2 | β 1 | -0.00201 | 0.000 | <0.0001 |
β 2 | 0.0419 | 0.003 | <0.0001 |
La predicción de percentiles y el método de momentos, así como las ecuaciones de proyección del crecimiento en diámetro mínimo y diámetro medio cuadrático del rodal, permitieron recuperar los parámetros de la fdp Weibull en el tiempo, esto fue posible debido a que dichos parámetros dependen de las variables de estado, por lo que, al proyectarlas en el tiempo, se recupera o conoce el valor de los parámetros a una edad de interés (Prodan, Peters, Cox y Real, 1997; Poudel y Cao, 2013; Santiago-García et al., 2014), es importante mencionar que por la estrecha relación entre momentos (o percentiles en su caso) y los atributos del rodal, esta técnica ha brindado modelos con ajustes satisfactorios (Borders y Patterson, 1990).
Con las ecuaciones ajustadas fue posible determinar la proporción de árboles para cada categoría diamétrica en el rodal (Fig. 4). Para estimar el rendimiento maderable con el método de recuperación de parámetros mediante predicción de percentiles se requieren los mismos modelos adicionales (ecuaciones de altura-diámetro y volumen individual) usados en el sistema de predicción de parámetros, al interrelacionar estas ecuaciones fue posible la predicción del rendimiento maderable por categoría diamétrica (Fig. 3 y Fig. 4).
Comparación de los métodos
Al comparar gráficamente el rendimiento maderable derivado de un modelo de totalidad del rodal (predicción explícita) diseñado para Pinus patula (Santiago-García et al., 2017) y los rendimientos predichos implícitamente (Fig. 5), es evidente la diferencia entre las estimaciones, porque el rendimiento maderable total predicho por el método de predicción de percentiles es mayor al volumen explícito, mientras que el método de predicción de parámetros predice el rendimiento maderable de forma similar al sistema explícito.
En una predicción implícita del rendimiento existen más errores acumulados que en una predicción explícita (porque se usan más ecuaciones para realizar la predicción implícita), esta última se usa para corregir las estimaciones de la primera. El procedimiento de corrección consiste en distribuir las diferencias de ambas proyecciones en área basal o volumen de acuerdo con la participación de cada categoría diamétrica en estas variables. El factor de corrección (FC) se calculó como:
donde:
VE : volumen total explícito (m3 ha-1)
VI : volumen total implícito (m3 ha-1) (Magaña, Torres, Rodríguez, Aguirre y Fierros, 2008)
Una vez obtenido este valor, se multiplica por el volumen de cada categoría diamétrica.
El FC (Tabla 8) se estimó para ambas predicciones implícitas y se obtuvieron las curvas de rendimiento volumétrico corregidas por dicho factor (Fig. 6); de esta manera, la estimación es más realista y ambos sistemas tienden a ser similares.
Edad | VTE | VTI | FC | Edad | VTE | VTI | FC | |||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
S1 | S2 | S1 | S2 | S1 | S2 | S1 | S2 | |||||
5 | 0.23 | 12.79 | 21.41 | 0.02 | 0.01 | 45 | 651.38 | 666.12 | 983.48 | 0.98 | 0.66 | |
10 | 25.64 | 65.18 | 77.62 | 0.39 | 0.33 | 50 | 685.82 | 716.61 | 1091.49 | 0.96 | 0.63 | |
15 | 123.51 | 143.73 | 168.84 | 0.86 | 0.73 | 55 | 712.71 | 755.29 | 1182.55 | 0.94 | 0.60 | |
20 | 253.98 | 242.23 | 287.85 | 1.05 | 0.88 | 60 | 733.94 | 783.39 | 1253.83 | 0.94 | 0.59 | |
25 | 374.71 | 344.17 | 428.00 | 1.09 | 0.88 | 65 | 750.87 | 801.91 | 1314.62 | 0.94 | 0.57 | |
30 | 472.72 | 440.66 | 575.75 | 1.07 | 0.82 | 70 | 764.53 | 811.99 | 1355.94 | 0.94 | 0.56 | |
35 | 548.67 | 527.68 | 722.81 | 1.04 | 0.76 | 75 | 775.65 | 814.84 | 1382.69 | 0.95 | 0.56 | |
40 | 606.81 | 603.16 | 859.50 | 1.01 | 0.71 | 80 | 784.78 | 811.46 | 1396.14 | 0.97 | 0.56 |
VTE: volumen total explícito; VTI: volumen total implícito; FC: factor de corrección: S1: predicción de parámetros, S2: predicción de percentiles.
Con el factor de corrección, las estimaciones son equivalentes, aunque con el sistema basado en predicción de percentiles, el volumen estimado por categoría diamétrica es ligeramente mayor que el estimado por el sistema basado en predicción de parámetros.
Discusión
Predicción de parámetros de la fdp Weibull
La ecuación para predecir el parámetro de localización (a) seleccionada por su bondad de ajuste fue:
Esta ecuación indica la relación existente entre el parámetro de localización (a), diámetro promedio (Dm) y el diámetro medio cuadrático (Dq), este resultado coincide con el registrado por De la Fuente et al. (1998), para rodales coetáneos de Pinus rudis Endl. en Pueblos Mancomunados (Oaxaca), y por Torres-Rojo, Magaña-Torres y Acosta-Mireles (2000), para P. montezumae Lamb. en San Juan Tetla (Puebla), ya que estos autores utilizaron tanto el Dm como el Dq para estimar el parámetro a, además, estas variables son fundamentales para caracterizar la distribución diamétrica del rodal.
Con respecto al parámetro de escala (b), la ecuación de predicción seleccionada fue:
esta es parcialmente similar a la propuesta por Quiñonez, De los Santos, Cruz, Velázquez y Ramírez (2015), quienes para estimar las distribuciones diamétricas en masas mezcladas de especies del género Pinus, sugieren como variables predictoras de este parámetro el diámetro mínimo, el diámetro medio cuadrático y el diámetro promedio del rodal.
El parámetro de forma (c) fue predicho de mejor manera con la ecuación:
Haan (1986) indicó que la técnica de momentos relaciona el coeficiente de asimetría (Ca) con el parámetro c, de acuerdo con esta condición Návar-Cháidez (2009) propuso predecir este parámetro en función del Ca; Magaña, Torres, Rodríguez, Aguirre y Fierros (2008) obtuvieron que para P. rudis Endl. en Aloapan (Oaxaca), la predicción del parámetro c es satisfactoria al considerar como variables independientes el diámetro medio cuadrático, la calidad de sitio (IS), el número de árboles y la altura promedio. En tanto que Quiñonez, De los Santos, Cruz, Velázquez y Ramírez (2015) sugieren como variables independientes el diámetro promedio, el diámetro medio cuadrático y el diámetro mínimo.
Predicción de percentiles
En la predicción de la distribución diamétrica, la elección de los percentiles es crucial y depende de la muestra, por lo que una pareja de percentiles al estar más centrada ayuda a reducir el sesgo en la predicción (Santiago-García et al., 2014). En este estudio, los percentiles 55 y 85 permitieron calcular el parámetro c de la fdp Weibull y obtener estimaciones adecuadas de la distribución diamétrica de los rodales de Pinus patula. Pienaar y Rheney (1993) y Montero y Fierros (2000) usaron los percentiles 24 y 93 para calcular el parámetro c. García, Cañadas y Montero (2002) señalan que con la pareja de percentiles 40 y 82 pueden conseguirse ajustes eficaces. En tanto que, Santiago-García et al. (2014) obtuvieron resultados satisfactorios con los percentiles 50 y 90 de la distribución diamétrica.
Comparación de los métodos
En la estimación del rendimiento maderable, el sistema de predicción de percentiles generó valores más altos que el sistema de predicción de parámetros, esto, debido a que las distribuciones diamétricas del rodal obtenidas con predicción de percentiles abarcaron más valores con probabilidad de ocurrencia (Fig. 4). En el caso del método de predicción de parámetros, el intervalo de las clases diamétricas estimadas para P. patula abarcó de 5 cm a 60 cm en un intervalo de edad de 5 años a 80 años, en tanto que, con el método de predicción de percentiles se obtuvieron valores de diámetro menores a 2.5 cm e iguales a 90 cm en el mismo intervalo de edad. Embrechts, Klüppelberg y Mikosch (1997) argumentan que las distribuciones subexponenciales como la Weibull son de colas pesadas, porque tienden a 0 lentamente. Por tanto, tienen más probabilidad de encontrar valores extremos repartidas en ellas. En este caso, la probabilidad mayor ocurrió para las categorías diamétricas superiores en el sistema de predicción de percentiles, en consecuencia, se obtuvieron valores más altos de rendimiento maderable.
Las estimaciones implícitas de rendimiento maderable se compararon contra un modelo de volumen explícito de P. patula existente en el área de estudio (Santiago-García et al., 2017). Se considera que este tipo de modelos tienen menor error en la estimación del rendimiento de rodales coetáneos, porque el nivel de resolución es directamente en m3ha-1 y se reconoce que representan un buen compromiso entre precisión y generalidad (García, 1988; Vanclay, 1994; Diéguez-Aranda, Castedo, Álvarez y Rojo, 2006). Por esta razón, las estimaciones realizadas con el método de predicción de parámetros de la fdp Weibull en función de las variables del rodal son más realistas que el sistema basado en predicción de percentiles, ya que las estimaciones fueron similares a las obtenidas con el modelo explícito.
En el presente estudio se estimó, con el sistema implícito de predicción de parámetros, que el volumen del fuste total con corteza (Tabla 8) a los 15 años fue de 143.73 m3 ha-1; en tanto que, con el sistema de predicción de percentiles, se estimó 168.84 m3 ha-1; para el área basal se estimó 16.99 m2 ha-1 y 19.21 m2 ha-1 para los sistemas de predicción de parámetros y predicción de percentiles, respectivamente, al considerar una densidad de 1326 árboles por hectárea. Al respecto, Arteaga (2003) obtuvo que el rendimiento en volumen de una plantación de P. patula de 14.6 años en Perote, Veracruz fue de 32.39 m3 ha-1 de fuste total sin corteza con una densidad de 976 árboles ha-1, área basal de 9.21 m2 ha-1 y altura media de 6.53 m. Por su parte, Santiago-García et al. (2014) estimaron, para los 15 años, cerca de 150 m3 ha-1 de fuste total con corteza, de acuerdo con una densidad inicial de 1400 árboles por hectárea a los cinco años de edad e índice de sitio de 29 m a la edad base de 40 años para rodales coetáneos de P. patula en el ejido La Mojonera, Municipio de Zacualtipán, Hidalgo.
Los sistemas de crecimiento implícitos presentados en este artículo son instrumentos que apoyan la planificación forestal, al permitir definir la distribución de productos maderables, los regímenes silvícolas que generen estructuras residuales deseables, o los tipos, intensidad y frecuencia de aclareos. Estos sistemas pueden complementarse con una ecuación de ahusamiento-volumen, que permita realizar la distribución de productos a nivel árbol individual.
Conclusiones
Los sistemas estudiados permiten obtener predicciones del rendimiento maderable por categoría diamétrica a nivel de unidad de superficie, son sencillos de implementar y permiten caracterizar los productos que se tendrán en el bosque en un determinado momento. Con los modelos de proyección de las variables dasométricas del rodal se recuperaron los parámetros de la función de distribución de probabilidades Weibull a diferentes edades de interés. Para este estudio, el sistema basado en predicción de parámetros arrojó resultados satisfactorios. Estos modelos permiten representar la estructura diamétrica de un rodal en términos numéricos a través de una tabla de rendimiento, la cual muestra el cambio de estructura de un rodal a medida que cambia la edad, por lo que constituyen una herramienta valiosa para planificar el manejo de los rodales de Pinus patula en Ixtlán de Juárez, Oaxaca, México. La precisión e intervalo de validez de los sistemas propuestos puede mejorarse al aumentar el número de mediciones y reajustar el sistema con cada remedición.