INTRODUCCIÓN
La industria forestal es una herramienta necesaria para hacer que la producción forestal sostenible sea un objetivo viable. También constituye un agente para el desarrollo sostenible por derecho propio, ya que aumenta los beneficios socioeconómicos del sector.
Pupo et al. (2018), plantean que actualmente es un objetivo de la política forestal crecer cada año en la producción de madera, fundamentalmente de plantaciones, aliviando la presión existente sobre los bosques naturales ya que en su gran mayoría las especies de alto valor comercial se encuentran distribuidas en estos bosques.
La demanda de materia prima continúa en aumento y la oferta de madera aún más limitada, jugando un papel fundamental la diversificación de especies de rápido crecimiento que, además de adaptarse a las condiciones edafoclimáticas, son de alta productividad y tienen madera de buena calidad, donde se destaca el género Eucalyptus(França et al., 2019).
Sin embargo, de acuerdo con Beltrame et al. (2016), la mayor limitación técnica que presenta este género es la magnitud de las tensiones de crecimientos que, probablemente, resulte el indicador de mayor significación, puesto que altos niveles de tensión (particularmente las tensiones longitudinales) ocasionan los defectos que más perjudican el rendimiento y la obtención de madera aserrada de calidad.
Las tensiones de crecimiento comienzan a liberarse desde el mismo momento de la tala y se manifiestan como rajaduras en los extremos de las trozas. En las piezas aserradas, las tensiones remanentes también provocan rajaduras y alabeos (França et al. 2019 y Silva et al., 2019). Por lo que el objetivo general del presente es incrementar los rendimientos de madera aserrada de Eucalyptus pellita y Eucalyptus saligna a partir de la utilización de la combinación de tratamientos para reducir el efecto de las tensiones de crecimiento y los esquemas de aserrado.
MATERIALES Y MÉTODOS
Metodología utilizada
La investigación se llevó a cabo en el aserradero de banda Rigo Fuentes, perteneciente a la Empresa Agroforestal La Palma, específicamente en la Unidad Empresarial de Base Agroindustrial "Los Palacios" ubicada en el Consejo Popular de San Diego de los Baños, Municipio Los Palacios, Provincia de Pinar del Río, Cuba.
El municipio Los Palacios, se localiza en la parte sur-este de la provincia de Pinar del Río. Limita al Norte con los municipios de La Palma, al Sur con las aguas del Golfo de Batabanó, al Oeste con el municipio Consolación del Sur y al Este con San Cristóbal, el territorio se extiende desde las elevaciones de la Cordillera de Guaniguanico hasta la costa sur bañada por el mar Caribe y posee una extensión territorial de 764,51 km2.
La Unidad Empresarial de Base Agroindustrial ¨Los Palacios¨ radica en la zona montañosa del municipio y cuenta con un patrimonio forestal de 12249,60 ha, con 199 trabajadores, distribuidos en seis brigadas para atender todo lo referido a la Silvicultura, al aprovechamiento y la industria forestal.
Durante la materialización de la investigación se utilizan dos especies de rápido crecimiento del género Eucalyptus: Eucalyptus pellita y Eucalyptus saligna en plantaciones que abarcan 373,6 ha de las 1451,70 ha que constituyen plantaciones establecidas del Patrimonio Forestal del Municipio Los Palacios.
Determinación del número de muestras utilizadas en la investigación
Se empleó el muestreo aleatorio simple a partir de la utilización de 100 árboles, teniendo como viable de interés el diámetro (D1.30); según la metodología propuesta por Aldana et al. (2010) y Álvarez et al. (2020); que utilizan las siguientes ecuaciones Ecuación 1 y Ecuación 2:
Donde: n = número de muestras; t 2 = Estadígrafo t de student al cuadrado; LE= límite del error de muestreo admitido.
La muestra queda constituida por 200 trozas de 4,0 m de longitud de cada especie, con un diámetro entre 21 y 30 cm, este número fue dividido en grupos de 50 para aplicar los diferentes tratamientos para reducir los índices de rajaduras de las trozas e incrementar los rendimientos de madera aserrada
Determinación de la eficiencia del proceso de aserrado de Eucalyptus pellita y Eucalyptus saligna en dependencia de los esquemas de aserrado y tratamientos para reducir las tensiones de crecimiento
A partir de la combinación de tres tratamientos para reducir el efecto de las tensiones de crecimiento sobre la calidad de las trozas; combinado con la aplicación de diferentes métodos de aserrado, propuestos por Pedro et al. (2014) citado por Juizo et al. (2018) y Carvalho et al. (2019) se pretende incrementar la eficiencia del proceso de aserrado de E. pellita y E. saligna. Los tratamientos son los siguientes:
Tratamiento de trozas almacenadas con riego (T1)
Este tratamiento consiste en mantener las trozas sometidas a riego durante 4 veces al día, con el objetivo de mantener un elevado gradiente de humedad en la madera para disminuir el efecto de las tensiones internas de crecimientos y a su vez hacer más fácil el aserrado de las trozas, según (Pupo et al., 2018).
Tratamiento de trozas almacenadas con anillo en los extremos (T2)
Consiste en hacer un anillo de 30.0cm en cada extremo de las trozas, profundizando 1/3 del diámetro, teniendo en cuenta lo planteado por Pupo et al. (2018).
Anillado de árboles en pie (T4)
Consiste en hacer un anillo de 30 cm de ancho por encima de la altura de talado de los árboles, profundizando 1/3 del diámetro con el objetivo de cortar la circulación de la savia eliminando el cambium del árbol, hasta llegar a la profundidad deseada, para lo cual fue necesario esperar tres meses. (Neto, 2017 y Pupo et al., 2020). A continuación, se obtiene una troza de cuatro metros de longitud ubicada por encima de la altura de talado y que contenga el anillado realizado.
Las trozas fueron trasladadas del área de tala para la Unidad Empresarial de Base Agroindustrial "Los Palacios" donde fueron montados todos los tratamientos. Como índice de determinación del efecto de los métodos aplicados para la reducción de las tensiones de crecimiento sobre la madera de Eucalyptus sp se toma como referencia el Índice de rajadura.
Determinación del índice de rajadura
En el desarrollo de los tratamientos se siguió de manera sistemática la cuantificación de las fendas o grietas producidas en las trozas debido a la liberación de tensiones de crecimiento en cuanto a longitud (LG), profundidad (PG), ancho (AG) y número (NG). Estos datos fueron procesados con las facilidades del sistema Microsoft Excel donde se determinaron las medias de cada uno de los parámetros cuantificados (Figura 1), así como el índice de rajadura de las trozas (IRT), según el tratamiento y la especie, utilizando la ecuación utilizada por Beltrame et al., (2015). (Ecuación 3).
Dónde: IR - índice de rajadura de las trozas, %; ai - abertura máxima de la rajadura (i = 1...n), cm; Ci - tamaño de la rajadura (medula-corteza), cm; D - diámetro medio de corte, cm.
Para la construcción de los modelos que permitieron la predicción del Índice de Rajaduras se tomó la muestra patrón (T3), teniendo como variables dependiente el IR y las variables independientes el diámetro en el centro de la troza y la conicidad de las mismas. Basándose la selección de ambos variables, principalmente en fundamentos conceptuales o teóricos (Hair et al., 1999). Para el ajuste de los modelos de estimación del Índice de rajadura se utilizó la regresión lineal múltiple, mediante el método por pasos.
Después de concluido el tiempo establecido en los tratamientos aplicados, las muestras fueron trasladadas al aserrío Martín González donde se llevó a cabo la segunda fase de la investigación.
Debido a la particularidad del Eucalyptus sp en cuanto a las elevadas tensiones internas de crecimiento y sus posteriores defectos en la madera aserrada se aplicaron a las muestras dos esquemas de aserrado, denominados (A)- aserrado tangencial con cortes sucesivos y el (B)- aserrado tangencial con cortes alternativos, tal y como son representados en la (Figura 1); coincidiendo con Carvalho et al. (2019).
Se procesaron 25 trozas por tratamiento para definir los rendimientos y calidad de la madera aserrada para las especies de Eucalyptus investigadas. Los rendimientos de madera aserrada se obtienen a partir de la siguiente expresión matemática, coincidiendo con Leyva et al., (2017) Ecuación 4.
Donde: Rv - Rendimiento volumétrico total, (%); Vma- Volumen de madera aserrada, (m3); VT- Volumen de madera en trozas (m3)
El volumen de madera aserrada total en cada troza en los aserraderos, se determinó sobre la base de las mediciones lineales obtenidas de madera aserrada de acuerdo con las expresiones que se exponen a continuación Ecuación 5:
Donde: - volumen de madera aserrada de una troza, m3;,, - ancho, grueso y longitud de la pieza i obtenida de una troza o grupo de troza, m; n- número de piezas aserradas de una troza
El cómputo del volumen en bruto de cada troza se efectuó a través de la fórmula de Smalian, teniendo en consideración estudios desarrollados por (Casagrande et al., 2019), como se muestra en la siguiente ecuación Ecuación 6:
Donde: Vmb: volumen de la madera en bolo (m3); g1: área basal en el extremo fino (m2); g2: área basal en el extremo grueso (m2); L: longitud de la troza (m).
Análisis estadísticos
Las variables longitud (LG), profundidad (PG), ancho (AG) y número de grietas (NG) fueron analizadas para conocer si se cumplían los supuestos teóricos de distribución normal y homogeneidad de varianzas. Se realizó un ANOVA de dos factores, teniendo en cuenta el comportamiento de cada variable dependiente en cada métodos o tratamientos aplicados en las especies investigadas, con excepción del caso específico del número de grietas que no cumple con los supuestos teóricos ya que no sigue una distribución normal y se le aplicaron transformaciones logarítmicas.
A través del test de Duncan a un nivel de 5 % de probabilidad de error se determinó la existencia de diferencias significativas entre los métodos o tratamientos aplicados en cada una de las variables cuantificadas.
Se realizó además un análisis de varianza para definir diferencias entre las especies investigadas y los métodos o tratamientos aplicados para reducir las magnitudes de las rajaduras (número, ancho, longitud y profundidad de las grietas) y el Índice de Rajadura de las trozas (IRT) a través de una prueba HSD Tukey considerando un valor de p= 0.05
Se utilizó un análisis t de muestras pareadas en el proceso de validación de los modelos de regresión obtenidos para la `predicción del Índice de rajadura en las trozas investigadas.
RESULTADOS Y DISCUSIÓN
Análisis de los tratamientos para reducir el efecto de las tensiones internas de crecimiento en la madera de Eucalyptus pellita y Eucalyptus saligna
En la Figura 2, se muestra el comportamiento del número de grietas por tratamiento en dependencia de la especie. Observando que la muestra patrón es la que presenta un mayor número de grietas, considerando que a estas trozas no se le ha aplicado ningún tratamiento para reducir el efecto de las tensiones de crecimiento sobre el desarrollo de grietas en los extremos; lo cual afecta directamente la efectividad del proceso de transformación mecánica de estas trozas en volumen y calidad de la madera procesada; coincidiendo con Braz et al. (2017) (Figura 2).
En las Figuras 3 y 4, se aprecia que los mejores resultados se obtuvieron en los tratamientos de trozas almacenadas con riego (T1) y en árboles anillados en pie (T4) para el E. saligna y en árboles anillados en pie (T8) para E. pellita (Figura 3 y Figura 4).
Los resultados expuestos en las Figuras 3 y 4, posibilitaron exponer una reducción de las dimensiones de las magnitudes de las grietas, lo cual es consecuencia de una mayor liberación de las tensiones de crecimiento; coincidiendo con Müller et al. (2017) y Silva et al. (2017).
Estos resultados fueron demostrados a través de un ANOVA de dos factores teniendo en cuenta el comportamiento de cada variable dependiente en cada prueba o tratamiento realizado para cada especie.
Al analizar el comportamiento de cada variable de forma independiente se resume que en la variable número de grietas (NG) se aprecia la interacción entre los factores (tratamientos y especies) y los mejores tratamientos sin diferencia significativa entre ellos, son los tratamientos de trozas almacenadas con riego (T1) y árboles anillados en pie (T4) para E. saligna, donde se obtuvo valores de 3,25 cm y 2,0 cm respectivamente y en árboles anillados en pie (T4) para E. pellita con un valor promedio de (1,66cm) (Tabla 1).
En el caso de la variable longitud de grietas (LG) no hubo diferencias significativas ente las especies ni entre los tratamientos. El tratamiento de trozas almacenadas con riego (T1) con 7,04cm en E. saligna y 7,29 cm en E. pellita no difiere estadísticamente del tratamiento trozas almacenadas con anillo en los extremos (T2) con 7,26 cm en E. saligna y 7,35 cm en E. pellita, indicando esto que la liberación de las tensiones de crecimientos fue más evidente en la aparición de las grietas en los extremos de las trozas.
Sin embargo, se precia la diferencia significativa entre el tratamiento de árboles anillados en pie (T4) con los demás tratamientos y el patrón (T3) considerando este, como el más efectivo para reducir las grietas en los extremos de las trozas para ambas especies con valores mínimos obtenidos de 4,65 cm en E. saligna y de 4,67 cm en E. pellita (Tabla 1).
Especies | Tratamientos | NG | LGcm | AGcm | PGcm |
Eucalyptus saligna | 1 | 3,25 d | 7,04 b | 0,018 a | 0,045 a |
- | 2 | 6,25 c | 7,26 b | 0,017 a | 0,043 a |
- | 3 | 9,00 ab | 10,28 a | 0,019 a | 0,048 a |
- | 4 | 2,00 d | 4,65 c | 0,010 b | 0,010 b |
Eucalyptus pellita | 1 | 8,00 bc | 7,29 b | 0,020 a | 0,049 a |
- | 2 | 7,50 bc | 7,35 b | 0,018 a | 0,043 a |
- | 3 | 10,25 a | 10,48 a | 0,020 a | 0,044 a |
- | 4 | 1,66 d | 4,67 c | 0,011 b | 0,010 b |
Nota: Valores con la misma letra demuestran que no existen diferencias significativas entre los tratamientos para una P < 0,05.
Lo cual es consecuencia de que al realizar el anillado se provoca la muerte de los árboles, perdiendo humedad, reduciendo el efecto de las tensiones de secado además otro motivo es la liberación parcial de las tensiones de crecimientos provocando una disminución en la magnitud de las grietas.
Silva et al. (2017), obtienen valores medios de agrietamientos aplicando el tratamiento de árboles anillados en pie en E. europhylla de (26,73 cm) mucho más altos que los obtenidos en la presente investigación. El ancho (AG) y la profundidad (PG) de las grietas muestran similar comportamiento que en la variable anterior. Con valores mínimos obtenidos de (0,011cm) en E. saligna y de (0,010cm) en E. pellita en el tratamiento de árboles anillados en pie, reafirmando la efectividad del mismo.
Análisis de la determinación del Índice de Rajadura
Mediante el análisis ANOVA de dos factores, se comprobó que no existen diferencias significativas entre las especies sobre el Índice de rajadura; considerando además la existencia de un efecto significativo del factor tratamiento sobre el Índice de rajadura (Tabla 2).
Análisis De Varianza | |||||
Fuente | Suma de Cuadrados | Gl | Cuadrado Medio | F | P |
A: Prueba | 0,415184 | 3 | 0,138395 | 27,40 | 0,0000 |
B: Especie | 0,00300312 | 1 | 0,0030012 | 0,59 | 0,4482 |
A x B | 0,403344 | 3 | 0,0134448 | 2,66 | 0,0709 |
Residuo | 0,121225 | 24 | 0,00505104 | - | - |
Total | 0,579747 | 31 | - | - | - |
Pupo et al. (2018) y Pupo et al. (2020), establecen que la presencia de rajaduras es consecuencia de la presencia de las tensiones de crecimiento y es por ello que a partir de la aplicación de diferentes tratamientos se pretende reducir el efecto de las mismas, las que se desarrollan con marcada significación en el género Eucalyptus, coincidiendo en este sentido con Silva et al. (2019).
En la Tabla 3, se comprueba la existencia de diferencias significativas entre los tratamientos aplicados, donde el T4 (anillado de árboles en pie) es el de mejor comportamiento con el menor índice de rajadura (0,15) en E. saligna y (0,26) E. pellita, aunque en los demás tratamientos aplicados se aprecia una tendencia al decrecimiento del mismo con relación al patrón (Tabla 3) y coincidiendo con Beltrame et al. (2015).
TRATAMIENTOS | IRT (%) E. saligna | IRT (%) E. pellita |
---|---|---|
1 | 0.32 b | 0,37 b |
2 | 0.29 c | 0.29 c |
3 | 0.56 a | 0.48 a |
4 | 0.15 d | 0.26 da |
Valores con la misma letra demuestran que no existen diferencias significativas entre los tratamientos para una P > 0,05.
En estudios realizados en diferentes clones de Eucalyptus sp por Beltrame et al. (2015), el valor promedio obtenido es de 0.46 % aunque en algunos de ellos los resultados fueron entre 0.14 y 0.34 %, muy similares a los obtenidos en la presente investigación.
Modelación matemática de la influencia de las variables diámetro y conicidad sobre la predicción del Índice de Rajadura de las trozas
Es muy importante la determinación del efecto de las características del árbol sobre los valores de madera, con la finalidad de proveer información para la selección de los modelos, por lo que se utilizan las variables de las trozas diámetro y conicidad (CON) y para determinar el grado de esta relación entre las variables empleadas en la estimación de los índices de rajaduras en las trozas se construyó la matriz de correlación que se muestra en la (Tabla 4 y Tabla 5).
Como se aprecia en las Tablas 4 y 5 los coeficientes de correlación de las variables empleadas en la estimación del índice de rajadura en las dos especies estudiadas son significativos, lo cual demuestra el alto grado de relación existente entre estas (Tabla 4 y Tabla 5).
- | Diámetro | IR | CON |
---|---|---|---|
DIÁMETRO Correlación de Pearson | 1 | -.857** | -.0.45 |
Sig. (bilateral) | - | .000 | .757 |
N | 50 | 50 | 60 |
IR Correlación de Pearson | -.857** | - | .201 |
Sig. (bilateral) | .00 | - | - |
N | 60 | 50 | 50 |
CON Correlación de Pearson | -.045 | .184 | 1 |
Sig. (bilateral) | .757 | .201 | - |
N | 50 | 50 | 50 |
** La correlación es significativa en el nivel 0.01 (dos colas)
- | Diámetro | IR | Conicidad |
---|---|---|---|
DIÁMETRO Correlación de Pearson | 1 | -.845** | -.757** |
Sig. (bilateral) | - | .000 | .000 |
N | 34 | 34 | 60 |
IR Correlación de Pearson | -.845** | 1 | .579** |
Sig. (bilateral) | .000 | - | .000 |
N | 34 | 34 | 34 |
CON Correlación de Pearson | -.757** | .579** | 1 |
Sig. (bilateral) | .000 | .000 | - |
N | 34 | 34 | 34 |
Nota: ** La correlación es significativa en el nivel 0.01 (dos colas)
Utilizando las variables independientes (Diámetro y Conicidad de las trozas) de mayor relación con la variable dependiente (IR) se obtienen ecuaciones de predicción del índice de rajadura de las trozas para las especies E. saligna y E. pellita y los modelos comparados se observan en la Tabla 6, con sus respectivos coeficientes de correlación r, determinación (r2) y error estándar de estimación (s), obteniendo un conjunto de resultados que denotan que existe una correlación lineal entre el índice de rajadura de las trozas, y las variables en estudio, reportando que existe una r cuadrado de 0,715 para E. saligna y de 0,806 para E. pellita; lo cual es considerado como fuerte. A partir de la aplicación del método de eliminación en el análisis de regresión lineal, la variable conicidad queda excluida para la construcción del modelo de predicción (Tabla 6).
Especie | Modelos de regresión | R | R2 | S |
---|---|---|---|---|
E. saligna | IRT= 1.474-0.046(Diámetro) | 0.845 | 0.715 | 0.09155 |
E. pellita | IRT= -0.606+0.047(Diámetro) | 0.898 | 0.806 | 0.07087 |
Las dos ecuaciones construidas para determinar el índice de rajadura en las trozas tienen elevados coeficientes de correlación, determinación y bajos errores estándar de estimación. Sin embargo, el mejor comportamiento lo posee el modelo que corresponde al E. pellita, que tiene mejor ajuste (las variables independientes explican el 89 % de la variable dependiente) y el menor valor del error estándar de estimación.
Validación de los modelos de regresión obtenidos
En las Tablas 7 y 8, se representa la comparación entre la predicción de la variable IR en un 50 % de nuevas muestras para validar con los valores reales de esa variable obteniéndose una exactitud en el comportamiento del modelo propuesto al no encontrarse diferencia significativa entre el valor real y estimado de los modelos para E. saligna y E. pellita (Tabla 7) y (Tabla 8).
- | Diferencias emparejadas | - | ||||||
- | Media | Desviación estándar | Media de error estándar | 95 % de intervalo de confianza de la diferencia | t | gl | Si. (bilateral) | |
Inferior | Superior | |||||||
Par 1 IR- Unstandarized Predicted Value | ,00173255 | ,079642299 | ,01592860 | .03114246 | ,03460756 | ,109 | 24 | ,914 |
- | Diferencias emparejadas | - | ||||||
- | Media | Desviación estándar | Media de error estándar | 95% de intervalo de confianza de la diferencia | t | gl | Si. (bilateral) | |
Inferior | Superior | |||||||
Par 1 IR- ajust | ,06958 | ,22310 | ,04554 | ‘022562 | ,16279 | ,1,506 | 23 | ,146 |
La validación del modelo garantiza su utilización con precisión en la práctica social para la predicción y control de la incidencia de la formación de grietas y su implicación sobre la calidad de la madera en bolo y su influencia sobre los rendimientos y calidad de la madera aserrada.
Análisis de los rendimientos de madera aserrada en función de los métodos de aserrado propuestos
Teniendo en consideración los métodos de aserrado experimentados se pudo establecer que en los resultados expuestos en la Tabla 9, los mejores rendimientos se obtuvieron para ambas especies con la utilización del método de aserrado con cortes tangencial alternativo y en el tratamiento de árboles anillados en pie donde se obtienen los valores menores de índice de rajadura (Tabla 9).
Tratamientos | IR- E. saligna % | IR E. pellita % | Rend- E. saligna CTCS % | Rend- E. saligna CTA % | Rend-E. pellita CTCS% | Rend- E. pellita CTA% |
Almacenamiento con riego | 0.32 | 0.37 | 46,34 | 49,98 | 45,75 | 46,95 |
Anillamiento de trozas en los extremos | 0.29 | 0.29 | 47,25 | 50,65 | 47,88 | 51,55 |
Patrón | 0.56 | 0.48 | 42,19 | 44,14 | 43,21 | 42,76 |
Anillamiento de árboles en pie |
0.15 | 0.26 | 48,19 | 51,87 | 49,10 | 52,55 |
Nota: IR- Índice de rajadura, CTCS- Corte tangencial con cortes sucesivos, CTA-Corte tangencial alternativo.
Anjos y Fonte (2017) en E. grandis y E. dunnii obtuvo valores de rendimiento entre 41.07 % y 43.85 %, 36.71 % y 41.02 % respectivamente, también inferiores a los obtenidos, aunque en E. saligna obtuvo valores entre 43.87 % y 47.46 %, este último muy similar a los resultados obtenidos en ambas especies con la utilización del esquema de aserrado tangencial con cortes sucesivos, aplicando el tratamiento de anillado en los extremos de las trozas.
También Carvalho et al., (2019), en Eucalyptus sp obtienen valores de rendimientos de 47.08 % similar al resultado anterior y de 53.43 % superior a los obtenidos. Melo et al, (2016), en estudios realizados con especies nativas obtienen resultados de rendimiento de 49.0 %, similar a los obtenidos en esta investigación.
Demostrando la necesidad de tener en consideración estos planteamientos para obtener los mejores resultados en relación a la eficiencia del proceso de trasformación primaria en los aserraderos para las trozas de E. saligna y E. pellita.
El modelo de simulación se sustenta en la solución del problema de corte de dos dimensiones en un espacio irregular que podría extenderse a una tercera dimensión. En su implementación computacional el algoritmo utiliza la información comercial de los pedidos de madera requeridos para optimizar el patrón de aserrado para una troza en específico.