text in 
English (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
Introducción
A nivel mundial, se registra una extracción de 1 300 000 t de resina de pino, la cual se utiliza, principalmente, como resina, colofonia y trementina en el mercado (Cunningham, 2009). La producción resinífera se concentra en tres países: China, Brasil e Indonesia (92 %), mientras que 8 % del total procede de Argentina, Rusia, India, México, Portugal, España, Vietnam y Fiyi (Cunningham, 2009; Conafor, 2013). En México, el aprovechamiento de resina se realiza en los estados de Jalisco, Oaxaca, Estado de México y Michoacán, este último ocupa el primer lugar con 22 408 t anuales, que equivalen a 82.7 % de la producción nacional (Semarnat, 2016).
En el país, después de los productos maderables, la resina de pino es el ingreso más seguro y permanente para la población rural más pobre de las zonas con bosques de pino y pino-encino, en particular en las entidades antes citadas (Ayala, 2011).
La resinación de pino en México se ha practicado desde la época Prehispánica (Mas y Prado, 1981); sin embargo, su aprovechamiento con fines comerciales inició en 1920. En la actualidad, se realiza de acuerdo con la NOM-026-Semarnat-1996 y la normatividad vigente que comprende la Ley General de Desarrollo Forestal Sustentable y su Reglamento (Semarnat, 2016).
El sistema de resinación utilizado corresponde al Método Tradicional Francés modificado para México e implementado en 1937 por disposiciones oficiales; en el estado de Michoacán inició su aplicación en 1954 (Reyes-Ramos et al., 2019). Este sistema es aplicable para árboles mayores a 30 centímetros de diámetro (d), el cual determina el número de cortes (caras) que puede soportar el árbol, cuya apertura se inicia a partir de los 20 años de edad del arbolado (Romahn, 1992).
Históricamente, la producción de resina en Michoacán ha tenido una destacada participación con un valor entre 70 000 t anuales en la década de los años 30 (Cofom, 2007), y 22 408 t en 2016 (Semarnat, 2016); en el 2017, se registró un repunte en la producción de poco más de 90 %; de tal manera que, en la actualidad, el estado ocupa el primer lugar nacional con una producción de 24 141 t (Semarnat, 2017). Las especies de mayor producción resinera en la entidad son: Pinus oocarpa Schiede ex Schltdl., P. leiophylla Schiede ex Schltdl. & Cham., P. lawsonii Roezl ex Gordon, P. teocote Schltdl. & Cham., P. herrerae Martínez, Pinus tenuifolia L., P. montezumae Lamb., P. ponderosa Douglas ex C. Lawson, P. pringlei Shaw y P. pseudostrobus Lindl. (Instituto Nacional de Investigaciones Forestales, 1973).
Algunas industrias y organizaciones han retomado la tarea de incrementar la producción de resina de pino a través de reforestaciones productivas en bosques resineros y plantaciones forestales comerciales no maderables, según lo permite la normatividad vigente, como es el caso de la Comunidad Indígena de Nuevo San Juan Parangaricutiro (CINSJP), donde 70 % de la población se dedica a las actividades forestales, entre las que destacan aquellas relacionadas con la industria del aserrío, el transporte y la resinación. La superficie con valor comercial para madera y resina que posee la CINSJP es de 18 132 ha de bosques; de ella, entre el 40 y 50 % se destina a la producción de resina, cuya extracción se centra en tres especies: Pinus leiophylla (pino chino), P. pseudostrobus (pino canis o pino blanco) y P. oocarpa (pino trompillo o artigo) (Arias y Chávez, 2006). Además, la CINSJP tiene una planta destiladora de resina, con una capacidad de 750 t año-1 (Aguilar, 2008).
En los últimos años, se ha generado información cuantitativa para algunos recursos forestales no maderables (RFNM). Entre los estudios publicados, sobresalen algunos modelos para la estimación de la producción de hongos silvestres comestibles (Tricholoma magnivelare (Peck) Redhead y Cantharellus formosus Corner) en Tlaxcala (Zamora-Martínez et al., 2009; Velasco et al., 2010); rendimiento para agaves mezcaleros (Agave karwinskii Zucc. y A. potatorum Zucc.) en Oaxaca (Velasco et al., 2009b); palma camedor (Chamaedorea oblongata Mart.) en Escárcega, Campeche (Zamora-Martínez et al., 2011); agave lechuguilla (Agave lechuguilla Torr.) en San Luis Potosí (Velasco et al., 2009a); resina en especies de P. pseudostrobus y P. leiophylla en Michoacán (Zamora-Martínez et al., 2013) y P. oocarpa en Chiapas (Reynoso-Santos et al., 2018; Cadena-Iñiguez et al., 2019).
A la fecha, el aprovechamiento de la resina de pino se ha realizado de manera empírica, y no se cuenta con estudios sobre la validación de modelos predictivos que estimen el potencial de producción. Por ello, es necesario tener un método estandarizado con técnicas cuantitativas que permita al productor conocer el potencial de recolecta por unidad de superficie y a los técnicos forestales evaluar el proceso de extracción de una manera práctica y confiable. Si existiera una relación entre la dendrometría del árbol (altura, diámetro normal y porcentaje de copa viva) y el rendimiento de resina, estos parámetros podrían utilizarse para la generación de un modelo matemático, ya que requiere menos esfuerzo que las mediciones directas del flujo de resina.
La hipótesis del estudio plantea que ciertas características dasométricas del árbol están relacionadas con su potencial de rendimiento de resina, lo que permite obtener estimaciones más precisas. Debido a lo anterior, se planteó como objetivo ajustar un modelo de predicción para la producción de resina por árbol a partir de las variables como el diámetro normal (d, cm), diámetro de copa (dc, m) y altura total (At, m) para bosques de P. pseudostrobus de la CINJSP, Michoacán, México. Los resultados se discuten en términos de identificar el mejor modelo matemático que estime con mayor precisión el rendimiento de resina en los rodales de P. pseudostrobus.
Materiales y Métodos
El estudio se realizó en las áreas productoras de resina dentro de los bosques naturales de la CINSJP, Michoacán, la cual se ubica entre las coordenadas 19°25′00″de latitud norte y 102°17’46’’ de longitud oeste. El intervalo altitudinal donde se llevó a cabo la investigación fue de 2 226-2 785 m; el relieve es de tipo Sierra, de geología ígnea extrusiva y suelo de tipo Andosol; el clima es templado húmedo (Cw), con precipitación de 1 500 mm y temperatura de 15 °C (INEGI, 2010; Inegi, 2017).
A través de recorridos de campo, se eligieron árboles productores de resina en 12 parajes. Los individuos seleccionados fueron representativos de los cuarteles de resinación de cada localidad, a los cuales se les medió el d, dc y At. El rendimiento de resina se midió por cara (PredResina[g]) de 215 árboles (2 186 caras de resinación), con una periodicidad de 45 días, durante los periodos 2009-2010, 2010-2011, 2016-2017 y 2018-2019. Cabe mencionar que los individuos evaluados en los distintos años no fueron los mismos y presentaron diferente etapa de aprovechamiento.
Una vez capturada la información en una base de datos en Excel®, se auditaron para descartar los puntos aberrantes a la muestra y se agrupó el dc a cada metro de longitud (5, 6 y 18 m) para generar la variable de clasificación por dimensión o clase de copa (ClasCop), esto con la finalidad de reducir la variabilidad al incluirse en un ajuste bajo modelos de efectos mixtos. Posteriormente, se calculó la estadística básica de la muestra, en la cual los valores de Curtosis y asimetría debieron estar entre el intervalo de 3 y -3 para considerar la información con distribución normal (Martínez et al., 2006).
Construida y depurada la base de datos, se incluyó una variable combinada que se construyó al elevar el d al cuadrado y multiplicarlo por la At (d 2 At). Posterior a ello, se transformaron a logaritmo natural las variables de d 2 At y PredResina (ln[d 2 At] y ln[PredResina]). En un primer ajuste fueron incluidas a un modelo lineal (1) (Ayala et al., 1992) y dos de sus variantes, (2) y (3), mediante mínimos cuadrados ordinarios (MCO) en el programa de licencia libre R ® con el procedimiento lm (R Core Team, 2016).
Donde:
y = Producción de resina por cara (g)
d = Diámetro normal (cm)
At = Altura total (m)
ln = Logaritmo natural
a, b y c = Parámetros por estimar
e = Término de error
En un segundo enfoque de ajuste, y al no corregir el modelo por el efecto de la media como lo propone Ayala et al. (1992), se incluyeron las variables de clasificación o agrupación de ClasCop, altitud y mes de colecta dentro de un modelo de efectos mixtos (MEM:
Esta modelación bajo MEM con la inclusión de las variables de agrupación, asume varianza específica para cada grupo de clasificación y generaliza el ajuste para la información empleada (Santana y Mateos, 2014; R Core Team, 2016); además ajusta de forma simultánea los valores de los parámetros específicos de cada nivel especifico; por lo cual, el error, varianza y valor serán únicos para cada grupo (Pinheiro et al., 1994; Castedo et al., 2006). En la modelación se asumió una varianza cero y distribución normal, sin embargo, se incluyó una matriz varianza-covarianza para cada nivel como lo proponen Littell et al. (2006), aplicado mediante la técnica EBLUP´S por Corral et al. (2019). También, se verifico de forma gráfica la normalidad de los residuales mediante una prueba gráfica (SAS Institute Inc., 2014).
Como criterio paramétrico de selección del mejor modelo, se utilizaron los de Akaike (4) y Bayesiano (5) (AIC y BIC, respectivamente), además del valor de verosimilitud del modelo (logLik), los cuales evalúan la verosimilitud y parsimonia de cada expresión, donde el de menor valor representa el modelo más adecuado, además de la significancia de todos sus parámetros (p=0.05) (Bolker et al., 2009; Cayuela, 2014; Correa y Salazar, 2016; R Core Team, 2016).
Donde:
AIC y BIC = Criterios de información de Akaike y Bayesiano
logLik = Valor de log-verosimilitud del modelo
k = Valor de penalización por el exceso de parámetros
n = Número de parámetros del modelo
log = Logaritmo base diez
Resultados y Discusión
La estadística descriptiva de la muestra final obtuvo 1 761 registros de producción de resina por cara e indicó que el intervalo de los datos fue de 37 a 109 cm en el d, mientras que la At varió de 19 a 45 m y el dc promedio fue de 10 m. Además, las tres variables mostraron una distribución mesocúrtica y simétrica de la información. Sin embargo, el rendimiento de resina por cara presentó una distribución leptocúrtica (Curtosis=16.153) y una asimetría positiva (Coeficiente de Asimetría=3.339), situación que se corrigió con la transformación logarítmica de esta variable (Cuadro 1).
Cuadro 1 Estadística descriptiva de la muestra utilizada para ajustar los modelos predictivos de resina en Pinus pseudostrobus Lindl., en la Comunidad Indígena de Nuevo San Juan Parangaricutiro, Michoacán.
| Estadísticos/Variables | Diámetro normal (cm) | Altura total (m) | Diámetro de copa (m) | Rendimiento resina (g cara-1) | ln Rendimiento resina (g cara -1 ) |
|---|---|---|---|---|---|
| Mínimo | 37.00 | 19.00 | 4.65 | 250.00 | 5.52 |
| Máximo | 109.00 | 45.00 | 18.25 | 1 420.00 | 7.91 |
| Media | 66.43 | 35.38 | 10.38 | 518.81 | 6.15 |
| Desviación estándar | 13.849 | 5.481 | 2.837 | 293.245 | 0.415 |
| Varianza de la muestra | 191.803 | 30.042 | 8.047 | 85 992.900 | 0.172 |
| Curtosis | -0.876 | -0.054 | -0.569 | 16.153 | 1.689 |
| Coeficiente de Asimetría | 0.147 | -0.318 | 0.379 | 3.339 | 1.062 |
En el ajuste bajo MCO se observó que al incluir la variable combinada por si sola (sin la transformación logarítmica), el parámetro b resultaba no significativo (p>0.05), además solo la expresión (1) mostró resultados adecuados. En el enfoque de MEM, el valor del intercepto (a) resultó ser de cero y el parámetro relacionado a esta variable de igual forma no fue significativo en todas sus combinaciones de inclusión de las variables de clasificación en sus parámetros. De igual manera, el valor de AIC, BIC y logLik, fueron superiores al modelo ajustado por MCO, lo cual indicó que no hubo ninguna mejora de ajuste (Cuadro 2).
Cuadro 2 Modelos, parámetros estimados y estadísticos de bondad de ajuste para los modelos predictivos de resina en su forma logarítmica (ln(PredResina)) en Pinus pseudostrobus Lindl., en la Comunidad Indígena de Nuevo San Juan Parangaricutiro, Michoacán.
| Modelo | Variable de clasificación | Parámetro de efecto | Parámetro | Estimación | Eea | Valor t | Pr>|t| | AIC | BIC | logLik |
|---|---|---|---|---|---|---|---|---|---|---|
| (1) | a | 3.4370000 | 8.71E-01 | 3.948 | <0.0001 | 1 290.54 | 1 310.90 | -641.27 | ||
| b | 0.2502000 | 8.01E-02 | 3.124 | 0.00183 | ||||||
| c | -0.0000016 | 5.11E-07 | -3.050 | 0.00234 | ||||||
| (3) | a | 5.9686200 | 2.63E-01 | 22.733 | <0.0001 | 1 297.83 | 1 313.09 | -645.91 | ||
| b | 0.0153000 | 2.20E-02 | 0.694 | 0.48800 | ||||||
| (1) | Altitud | b | a | 5.1963720 | 8.13E-01 | 6.393 | <0.0001 | 1 095.79 | 1 131.40 | -540.89 |
| b | -0.0000010 | 5.00E-07 | -1.314 | 0.18890 | ||||||
| c | 0.1060970 | 7.28E-02 | 1.457 | 0.14540 | ||||||
| (3) | a | a | 5.7304590 | 2.81E-01 | 20.380 | <0.0001 | 1 068.47 | 1 098.99 | -528.23 | |
| b | 0.0527240 | 2.00E-02 | 2.634 | 0.00860 | ||||||
| (1) | Clase de copa | b | a | 3.4303940 | 8.98E-01 | 3.822 | 0.00010 | 1 325.09 | 1 360.70 | -655.55 |
| b | -0.0000020 | 5.00E-07 | -2.859 | 0.00430 | ||||||
| c | 0.2499670 | 8.24E-02 | 3.033 | 0.00250 | ||||||
| c | a | 3.3335510 | 8.75E-01 | 3.808 | 0.00010 | 1 327.47 | 1 363.08 | -656.74 | ||
| b | -0.0000020 | 5.00E-07 | -3.223 | 0.00130 | ||||||
| c | 0.2597660 | 8.05E-02 | 3.227 | 0.00130 | ||||||
| (3) | a | a | 6.0033800 | 3.42E-01 | 17.533 | <0.0001 | 1 307.76 | 1 338.29 | -647.88 | |
| b | 0.0119820 | 2.83E-02 | 0.423 | 0.67230 | ||||||
| (1) | Mes de colecta | b | a | 4.1428180 | 8.77E-01 | 4.722 | <0.0001 | 1 287.60 | 1 323.21 | -636.80 |
| b | -0.0000010 | 5.00E-07 | -2.149 | 0.03190 | ||||||
| c | 0.1846210 | 8.05E-02 | 2.294 | 0.02200 |
Eea = Error estándar aproximado; Pr>|t| = Probabilidad p=0.05; AIC y BIC = Criterios de información de Akaike y Bayesiano; logLik = Valor de log-verosimilitud del modelo; a, b y c = Parámetros estimados.
En el ajuste de MEM al emplear la variable ln(PredResina) e incluir los distintos efectos por variable de clasificación, se obtienen valores no significativos en algunos parámetros (i.e., [1] > variable de clasificación=altitud > Parámetro de efecto=b > parámetro no significativo=b y c; y [3] > variable de clasificación=Clase de copa > Parámetro de efecto=a > parámetro no significativo=b), y en el caso de la inclusión de los efectos por clase de copa, los valores de AIC, BIC y logLik no mejoraron (aumento del valor) por lo cual se descartaron para su elección (Cuadro 2). Al incluir las variables de clasificación en el parámetro b (relacionado a la variable independiente) de altitud y mes de colecta, evidenciaron mejoras en los valores de los criterios paramétricos de selección del mejor modelo en el que disminuyeron los valores de AIC, BIC y logLik, fue superior estadísticamente la expresión con efectos aleatorios por altitud de colecta de resina por cara en árboles de P. pseudostrobus (Cuadro 2).
El modelo seleccionado como el mejor ajustado por MEM=[3] > variable de clasificación=altitud > Parámetro de efecto=a (Cuadro 2) presentó una mejora estadística de 17.21 % (222.07), 16.16 % (211.90) y 17.63 % (113.04) en los valores de los criterios paramétricos de selección de AIC, BIC y logLik, respectivamente, y con respecto al ajuste tradicional por MCO. Además, al verificar de forma gráfica los supuestos de regresión de normalidad y homocedasticidad de los residuales, la gráfica de cuartiles tiende a la línea recta, lo que es deseable (Martínez et al., 2006); y los residuales no aumentan con respecto al incremento de la variable explicativa, por lo que no se observan problemas de incumplimiento de los supuestos de regresión (Martínez et al., 2006) (Figura 1 y 2). Por tal motivo, se afirma que el modelo es estadísticamente robusto y superior a los demás.
Figura 1 Cumplimiento de los supuestos de normalidad (a) y homocedasticidad de los residuales (b) en el modelo seleccionado para predecir la resina en Pinus pseudostrobus Lindl., en la Comunidad Indígena de Nuevo San Juan Parangaricutiro, Michoacán.
Figura 2 Proyección del rendimiento resinero (g cara-1 por árbol) por categoría diamétrica a una altura estándar de 25 m para Pinus pseudostrobus Lindl., en la Comunidad Indígena de Nuevo San Juan Parangaricutiro, Michoacán.
Una vez elegido el mejor modelo y mediante la función summary(coef()) del programa R ® , se extrajeron los valores específicos de la ordenada al origen (parámetro a) para cada nivel de agrupación en el que se obtuvieron valores de 5.440657, 5.945814 y 5.804905 en las altitudes de 2 300, 2 500 y 2 700 m, respectivamente; mientras que el parámetro relacionado con la variable d2At fue constante en todos los casos (b=0.05272409). Las expresiones propuestas y la forma de implementación fue la siguiente:
Modelo general:
Modelo para 2 300 m de altitud:
Modelo para 2 500 m de altitud:
Modelo para 2 700 m de altitud:
Donde:
y = Producción de resina por cara (g)
d = Diámetro normal (cm)
At = Altura total (m)
ln = Logaritmo natural
exp = Función exponencial
e = Término de error
Al sustituir los valores de los parámetros por altitud especifica en la expresión (3) y aplicarlos en árboles con una altura promedio de 25 m y categorías diamétricas de 30 a 110 cm, se observaron valores más altos en la mejor producción de resina (g cara-1): en los árboles de P. pseudostrobus fue a los 2 500 m de altitud, posterior 2 700 m y la más baja a 2 300 m. Además, en un árbol de 60 cm de d y una At de 25 m, al considerar por lo menos cinco colectas por año de producción, y que solo tenga una cara, el rendimiento varió de 2 103.90 g cara-1 por árbol a 2 300 m de altitud a 3 486.68 g cara-1 por árbol a 2 700 m (Figura 2).
A la luz de las evidencias observadas en la Figura 2 y con la finalidad de verificar si los resultados eran estadísticamente diferentes, se decidió realizar un análisis de varianza (ANOVA) entre las tres estimaciones de productividad de resina por categoría diamétrica a una altura estándar de 25 m, esto al considerarlas como poblaciones independientes. Para ello, se utilizó un ANOVA para comparar las medias mediante la prueba de t con un nivel de confiabilidad del 99 % (p=0.01) (Martínez et al., 2006) en el que se plantearon las hipótesis de igualdad de medias (hipótesis nula: Ho) y que por lo menos una es distinta a las demás (hipótesis alternativa: Ha).
Los resultados del ANOVA mostraron diferencias significativas entre las estimaciones (F=387.5, p=<0.001), por lo cual se rechaza Ho y se acepta Ha. Al realizar la separación de medias con la prueba de t y graficarla se ratificó que la producción de resina fue superior a una altitud de 2 500 m, seguida de los individuos en el intervalo que considera 2 700 m, y el área de menor productividad correspondió a la altitud de 2 300 m (Figura 3). En la Figura 3, se observa que el modelo general (producción media) sobrestima el rendimiento de 2 300 m y subestima la producción de resina de los intervalos altitudinales de 2 500 y 2 700 m.
Figura 3 Separación de medias (t) para las estimaciones realizadas entre las categorías diamétricas de 30 a 110 cm a una altura estándar de 25 m con los distintos modelos propuestos.
Con la validación del ANOVA, se observó que las pruebas de Shapiro-Wilk y de Bartlett indicaron falta de evidencia suficiente para rechazar la hipótesis nula (Moore, 2005), la normalidad de los residuos y que las varianzas son iguales, respectivamente, a un nivel de significancia de p=0.01 (99 %), debido a que W=0.9623 y p=0.3759, mientras que K 2 =3.9137 y p=0.2709. Lo anterior se ratifica con las gráficas q-q plot de normalidad (Figura 4a) y de homocedasticidad de residuales por grupo (Figura 4b).
Figura 4 Pruebas gráficas de validación del análisis de varianza (ANOVA) de normalidad (a) y homocedasticidad (b).
La distribución altitudinal muestreada para P. pseudostrobus en la CINSJP fue de 2 200-2 800 m, con intervalos altitudinales de 200 m, incluido el óptimo de desarrollo de la especie (2 400-2 600 m) que se consideró como el de la más alta producción de resina (g cara-1), la cual se obtuvo a 2 500 m de altitud, la siguiente a 2 700 m, y la más baja en altitudes de 2 300 m. En trabajos similares con P. oocarpa, los árboles seleccionados como buenos productores de resina se registraron en altitudes de 1 192-1 789 m (Reyes et al., 2019). Cabe mencionar que para esta especie, su óptimo desarrollo se registra en altitudes de 1 200-1 800 m (Dvorak et al., 2000a). Al respecto, existen pocos estudios que sustenten cual es la altitud ideal para tener una mejor producción de resina; no obstante, algunos autores (Nanos et al., 2001; Calama et al., 2010; Rodríguez-García et al., 2014) reconocen que la combinación de factores ambientales en espacio y tiempo, y su interacción con la dinámica del rodal, han sido identificadas como causas de variabilidad en el rendimiento de la resina. Estos resultados pueden explicarse por las diferentes condiciones de índice de sitio presentes a diferentes altitudes. El índice de sitio constituye el principal impulsor de los modelos empíricos de crecimiento y rendimiento, dedicados a la producción de madera, la cual está relacionada con variables dasométricas como el d y la At (Calama et al., 2010).
La producción de resina mensual promedio obtenida (0.250-1.420 kg, en promedio 0.518 kg) fue menor a lo citado por Plesníková (2014) en P. pringlei, con 0.3 a 4.3 kg de resina mensual, en una distribución altitudinal de 1 500 a 2 500 m (Dvorak et al., 2000b). Así como a lo documentado por Reyes-Ramos et al. (2019), quienes obtuvieron de 0.5 kg a 4.5 kg de resina mensual promedio para P. oocarpa. Sin embargo, la producción de resina mensual se comportó de manera similar a lo indicado por el Instituto Nacional de Investigaciones Forestales (1973), con una producción promedio anual para P. oocarpa de 5.5 kg (0.4 kg mes-1) y para Pinus pinaster Aiton, con valores estimados de 3.5-3.8 kg árbol-1 año-1 (Tadesee et al., 2001; Rodríguez-García et al., 2014).
Únicamente para las variables de clasificación de altitud y mes de colecta, los valores de los criterios paramétricos de selección del mejor modelo mostraron mejoras, y el mejor modelo fue el (3) al incluir la altitud y por ser estadísticamente robusto y óptimo a los demás; destacó que los árboles con una altura promedio de 25 m y categorías diamétricas de 30-110 cm presentaron mejor producción de resina (g cara-1), en comparación a otras altitudes (2 300 y 2 700 m).
En trabajos similares con otras especies de pino se ha señalado que la altura de fuste limpio y el número de caras vivas por individuo, son variables que determinan la producción de resina (Reyes et al., 2019). Así, Reynoso-Santos et al. (2018) consignan que los modelos desarrollados a partir de la variable combinada d y At, explican en una buena proporción la producción de resina de pino. Otras de las variables seleccionadas por el modelo fueron la clase de copa y la época de colecta, esto es comprensible, si se considera que al aumentar el tamaño de la copa y en función de la época o temporada de colecta puede aumentar o disminuir la producción de resina (Rodríguez-García et al., 2014).
De acuerdo con la Conafor (2013), la producción de resina se realiza en dos temporadas del año: la temporada alta de marzo a junio, y la baja de julio a febrero. La producción mensual en la estación seca es de 2.2 t mes-1 y de 1 t mes-1 en la temporada baja. Para Michoacán, la temporada seca (especialmente de marzo a junio) es la época más recomendable para la recolecta y, en promedio, un pino resinero produce un kilogramo de resina al mes; cantidad que puede variar en función del número de cortes verticales o caras que se realizan en el tronco (Zamora-Martínez et al., 2013).
Dvorak et al. (2002) registran para P. tecunumanii F. Schwerdtf. ex Eguiluz & J. P. Perry en áreas resinadas, una producción promedio de 0.599 kg árbol-1 mes-1; durante la época seca en especies de P. caribaea var. hondurensis (Sénécl.) W. H. Barrett & Golfari y P. oocarpa obtuvo una producción promedio mensual de 0.476 kg árbol-1 mes-1.
La producción de resina en las diferentes especies presentan rendimientos no constantes durante el transcurso del año, estos están condicionados a diversos factores como la influencia del clima (temperatura, precipitación, humedad, principalmente), factores naturales (procesos fisiológicos complejos tales como la asignación de carbohidratos para el crecimiento), o propios del árbol que se relacionan con las características genéticas de la especie, y entre ellas, la variabilidad dentro del árbol propiamente (Nanos et al., 2001; Tadesse et al., 2001; Calama et al., 2010), así mismo, la edad de los árboles, diámetro de copa, diámetro del fuste, distribución y densidad de los canales resiníferos, la fertilidad del suelo, entre otros (Tadesse et al., 2001; Reynoso-Santos et al., 2018; Reyes-Ramos et al., 2019).
Si se considera que a nivel nacional la producción media anual de resina es de 2.0 kg árbol-1, y para especies netamente resineras como P. oocarpa es de 3.5 kg (Reyes-Ramos et al., 2019), los árboles de P. pseudostrobus localizados en el estrato altitudinal de 2 500 m en la CINSJP son individuos que por su alta productividad en la zona, pueden considerarse como una variable de selección en la etapa inicial de un programa orientado al mejoramiento genético para la producción de resina.
Los modelos generados difieren a los elaborados para otros productos forestales no maderables tales como agaves mezcaleros (A. karwinskii y A. potatorum) (Velasco et al., 2009b); producción de hongos silvestres comestibles (Zamora-Martínez et al., 2009; Velasco et al., 2010); palma comedor en Escárcega, Campeche (Zamora-Martínez et al., 2011); agave lechuguilla (Agave lechuguilla) en San Luis Potosí (Velasco et al., 2009a); resina de P. pseudostrobus y P. leiophylla en Michoacán (Zamora-Martínez et al., 2013) y P. oocarpa en Chiapas (Reynoso-Santos et al., 2018; Cadena-Iñiguez et al., 2019), porque en este estudio se integró al análisis el intervalo altitudinal de distribución de la especie en Michoacán, ubicada en la Sierra Purépecha donde se localiza la CINSJP. Esta última se presenta en altitudes de 2 200 hasta 3 000 m (Madrigal, 1982).
Conclusiones
El modelo propuesto, al considerar variables dasométricas de fácil medición como el diámetro normal y la altura total del arbolado, así como los parámetros globales resultado de la regresión bajo el enfoque de modelos de efectos mixtos, es robusto y estadísticamente confiable para estimar la producción de resina por cara y posteriormente por individuo, en el cual al utilizar los valores aleatorios de sus parámetros por nivel de agrupación, puede emplearse de manera específica por gradiente altitudinal.
La altitud con producciones resineras aparentemente más altas se ubica entre 2 400 y 2 600 m, de tal manera que a medida que disminuye o incrementa el nivel sobre el mar de la topografía, el rendimiento es menor. Por lo tanto, esta información puede emplearse para la planeación de la gestión forestal no maderable en los bosques de la Comunidad Indígena de Nuevo San Juan Parangaricutiro, Michoacán, en donde crece y se desarrolla P. pseudostrobus. Además, constituye una herramienta para determinar el potencial de producción por unidad de superficie o en el establecimiento de plantaciones forestales comerciales no maderables para la producción de resina.
