PROPIEDADES FÍSICAS EN Guazuma crinita POR MEDIOS CONVENCIONALES Y ESPECTROSCOPIA INFRARROJA CERCANA

PHYSICAL PROPERTIES OF Guazuma crinita BY CONVENTIONAL METHODS AND NEAR INFRARED SPECTROSCOPY

La correcta utilización de la madera, implica el conocimiento de sus características anatómicas (fibras y ángulo microfibrilar), físicas (densidad y contracciones); mecánicas (rigidez), y químicas (contenido de celulosa y de lignina). Dicha información está directamente asociada a la calidad y el uso del producto final (Oliveira mencionado por Filgueira et al. 2015). Las técnicas tradicionales para determinar las propiedades de la madera son destructivas; sin embargo, en la actualidad se han desarrollado técnicas para aumentar la eficiencia y rapidez de resultados de estas propiedades, siendo la Espectroscopía de Infrarrojo Cercano (NIR) un ejemplo. Esta técnica es utilizada en diferentes países industrializados por su probada utilidad, sin embargo, en el Perú no se registra su uso con fines forestales.

La evaluación de las propiedades físicas de la madera se realiza utilizando métodos directos e indirectos. El método directo más común para determinar la densidad básica consiste en un procedimiento gravimétrico, método destructivo que requiere mucho tiempo y limita su aplicación en muestreos de campo (Rybníček et al. 2012). En los últimos años, la espectroscopia aparece como una técnica ideal para reemplazar a los métodos tradicionales de laboratorio, teniendo como ventajas realizar un gran número de mediciones en cortos intervalos de tiempo, análisis cuantitativos y cualitativos rápidos, poca preparación y no destrucción de las muestras (Tsuchikawa y Kobori 2015, Lang et al. 2015). La espectroscopía NIR es una técnica que forma parte del campo de estudio de la espectroscopía molecular, la cual analiza la interacción de la radiación electromagnética con la materia y cuyo rango espectral está comprendido entre los 750 - 2500 nm (Reis-Milagres et al. 2013).

Dentro de las propiedades de la madera, la densidad es sin duda la propiedad más empleada para la evaluación de la misma por ser de fácil determinación y presentar una estrecha relación con la mayor parte de sus demás características. Es un parámetro importante para la evaluación de la calidad de madera y para otras estimativas en el análisis de las plantaciones forestales (Suckow et al. 2009). Otro índice importante dentro de las propiedades físicas son las contracciones de la madera que permiten evaluar su estabilidad dimensional; ya que las maderas que se contraen menos son más estables (Gonçalves et al. 2009). La magnitud de la variación dimensional es normalmente mayor en maderas con mayor peso específico debido a la mayor cantidad de madera por unidad de volumen; además, maderas de mayor peso específico para un mismo contenido de humedad contienen más agua en la pared celular (Silva et al. 2006).

En el Perú es poca la información en términos de variación de las propiedades de las especies maderables utilizadas en plantaciones forestales siendo la especie Guazuma crinita una de ellas (Weber et al. 2011, Sotelo et al. 2008). En los últimos años se han realizado diversos trabajos con espectroscopia NIR en distintos materiales y con el fin de estimar distintas propiedades. En el campo forestal se ha utilizado la espectroscopía NIR en diversas líneas de investigación; así tenemos a Viana et al. (2010), Vimal et al. (2015), Cassia et al. (2016) y Lengowski et al. (2018) en la predicción de la composición química en madera y papel; en madera y carbón a Bolzon et al. (2016). En propiedades físicas y mecánicas de Mimosa tenuiflora a Diesel et al. (2014); en la influencia de la anisotropía en los resultados estadísticos, a Costa et al. (2018); y en la estimación de carbono en bambú a Dallagnol (2012). En la determinación de origen y trazabilidad de la madera a Silva et al. (2018) y Paredes et al. (2018); en discriminación de especies y clasificación de maderas de distinta procedencia con espectrómetro NIR portátil a Soares et al. (2017).

Uno de los géneros más estudiados por espectroscopia es el género Eucalyptus; así para densidad básica, Reis-Milagres et al. (2013) obtuvo coeficientes de determinación de 0,84 a 0,89; Viana et al. (2010), valores de correlación de 0,56 a 0,97 y Tyson et al. (2009) de 0,75; así mismo, diversos trabajos relacionados con las propiedades físicas utilizando NIRs demostraron buenas estimaciones de esta propiedad (Hein et al. 2009). Así mismo, Alves et al. (2012), en un estudio de la densidad básica de Pinus pinaster y Larix x eurolepis hallaron un RPD de 3,5 y 3,2 respectivamente, con coeficientes de determinación entre 0,82 y 0,92. Por otro lado, Santos et al. (2012), reportan para densidad básica de Acacia melanoxylon, RPDs de 2,8 a 4,1 y coeficientes de determinación de 0,84 a 0,94.

Los objetivos de este trabajo fueron determinar las propiedades físicas (densidad básica, contracción volumétrica, tangencial y radial) de G. crinita por el método destructivo para edades de 5 y 8 años; desarrollar modelos de calibración multivariada empleando ambas edades y evaluar la capacidad predictiva de cada modelo en la calibración y validación de los modelos.

La madera de árboles de Guazuma crinita Mart. proviene de plantaciones forestales de la Empresa Reforestadora Amazónica S.A.; sector Maquizapallo localizada en la ciudad de Puerto Inca, Huánuco - Perú; ubicada geográficamente a 9° 37’ 43’’, latitud S y 74° 58’ 28’’ longitud W; 145 m de altitud (Figura 1). La zona tiene una precipitación media anual de 2482 mm y temperatura mínima de 22 °C y máxima de 35 °C. Los suelos de la plantación son de tipo aluvial de textura franco y franco-limoso; pH 5,5, con contenido de materia orgánica 4%; CIC 10,8 meq/100 g y profundidad efectiva mayor a 90 cm. Las muestras de madera fueron colectadas en campo en noviembre del 2016. El espaciamiento de los árboles en la plantación fue de 3 m x 3 m.

Corte, colección y preparación de muestras de madera

Los arboles al momento del corte tenían cinco años (63 meses) y ocho años (103 meses). Se seleccionaron cinco árboles por edad; se midió la altura comercial, tumbados y divididos en tres partes (basal, media y apical) conforme a la metodología utilizada por Valente et al. (2013). Se tomó un metro de la parte inferior de cada troza de cada nivel y de cada una de ellas se obtuvo una rodaja de 5 cm de espesor, de donde se obtuvieron muestras diametrales de 2 mm de espesor y 2 cm de ancho para su evaluación por espectroscopia. De lo restante se obtuvo una porción central de 5 cm de espesor para la obtención de probetas para la determinación de las propiedades físicas; conforme a la Figura 2.

Figura 1:
Localización de las plantaciones de Guazuma crinita en Huánuco (Perú).

Las muestras para los ensayos destructivos fueron preparadas en el taller de trabajabilidad de la Universidad Nacional Agraria la Molina (UNALM), codificadas y guardadas en bolsas herméticas, buscando mantener la máxima condición de humedad, para luego ser medidas y analizadas en el Laboratorio de Ensayos Tecnológicos de la UNALM. Para la determinación de las propiedades físicas se utilizaron las normas: NTP 251.011 (2014), para determinación de densidad y NTP 251.012 (2015), para determinación de las contracciones: volumétrica, tangencial y radial.

Para los ensayos físicos se prepararon 45 probetas por árbol correspondiente a tres niveles axiales totalizando 225 probetas por edad; habiéndose determinado en ellas la densidad básica y contracciones. Las dimensiones de las muestras fueron de 3 cm x 3 cm x 10 cm (tangencial x radial x longitudinal); las que fueron pesadas y medidas con micrómetro de precisión 0,01 mm en estado húmedo y llevado a estufa (103 ± 2 °C) hasta estado anhidro, para determinar densidad básica y contracciones totales. Las mediciones en las muestras se realizaron en dos puntos: uno en dirección tangencial y otro en dirección radial, tomando como referencia el centro de las caras antes indicadas. Para la contracción volumétrica, el volumen fue determinado por el método de desplazamiento de agua (Arquímedes); y para determinar el volumen en estado anhidro las probetas fueron parafinadas para evitar la absorción de agua.

Las muestras consistieron en láminas de 2 mm de espesor obtenidas diametralmente a partir de las porciones centrales de las rodajas. Se obtuvo una lámina por cada uno de los tres niveles de altura de cada uno de los 10 árboles para las dos edades. Para la medición de los espectros se utilizó un espectrofotómetro NIR Bruker Tensor 37 que obtuvo espectros a cada 2 nm de longitud de onda, en el intervalo de 1000 a 2500 nm con resolución de 1 nm, totalizando 300 lecturas para ambas edades. Se tomaron 10 lecturas equidistantes (5 a cada lado de la medula) para cada muestra y se utilizó el promedio por nivel. Las muestras no recibieron ningún tratamiento químico y las condiciones de trabajo en laboratorio fueron de 20 ºC y 65% HR. La información de las mediciones y procesamiento de los espectros fue registrada por el software OPUS SPECTROSCOPY del equipo BRUKER y procesados por el software Unscrambler X Versión 10.4.

Figura 2:
Toma de muestra para determinación de propiedades físicas y espectroscopia.

Los valores de densidad básica (DB) y de contracción volumétrica, tangencial y radial; son mayores para la edad de 8 años (Tabla 1). Los resultados de densidad básica promedio (360 kg/m³ para 5 años y 400 kg/m³ para 8 años) clasifican la madera de G. crinita como madera de baja densidad, conforme a la clasificación de Souza et al. (2002). Dichos valores son similares comparados con 430 kg/m³ para la misma zona y edad de 8 años (Tuisima et al. 2017), 420 kg/m³ para 32 meses de edad (Weber y Sotelo 2008), y 410 kg/m³ para madera proveniente de bosque natural (Acevedo y Kikata 1994). A pesar que existen diferencias estadísticas significativas entre la densidad básica de 5 y 8 años, Zobel y Jett (1995) señalan que la edad no es determinante en esta propiedad, sino que depende de su carácter genético; mientras que en un estudio de once procedencias para la misma especie, Weber y Sotelo (2008) indican que las diferencias en esta propiedad dependen de la procedencia.

Tabla 1:
Valores medios de propiedades físicas de Guazuma crinita para 5 y 8 años.

Axialmente los valores de densidad básica y contracciones son mayores en la parte basal del árbol en ambas edades (Tabla 2) y disminuyen de la base a la parte apical del árbol, evidenciando diferencias estadísticas en los tres niveles evaluados. La variación longitudinal encontrada sigue uno de los patrones de variación presentado por Panshin y De Zeeuw (1980), así como patrones para G. crinita y Calycophyllum spruceanum reportados por Tuisima et al. (2017) y Weber y Sotelo (2008). La mayor densidad básica en la parte basal del árbol se debe a la mayor compactación por la superposición de células en este nivel del tronco; así como a una respuesta adaptativa al estrés de flexión producido por el viento (Ali et al. 2010, Weber y Sotelo 2008).

Tabla 2:
Valores medios de propiedades físicas de Guazuma crinita para 5 y 8 años en dirección axial.

Los valores de contracción para la especie están acorde con lo reportado por Tuisima et al. (2017) y considerados bajos según la clasificación de Durlo y Marchiori (1992). Las contracciones son mayores para la edad de 8 años y en ambas edades son mayores en la parte basal. Longitudinalmente se encontraron diferencias estadísticas significativas en contracción volumétrica para ambas edades en el nivel basal en relación a los otros dos niveles que fueron iguales. En G. crinita de 5 años, la contracción tangencial y radial presenta diferencias estadísticas significativas entre los 3 niveles cuyos valores disminuyen de la base al ápice. El mismo comportamiento se da para 8 años excepto para contracción tangencial entre el nivel intermedio y apical, los cuales no presentan diferencias significativas. Los coeficientes de variabilidad en el presente estudio son mayores a los de Tuisima et al. (2017).

Los parámetros estadísticos de los modelos NIR desarrollados para las propiedades físicas de G. crinita se presentan en la Tabla 3. Los resultados de calibración muestran que los modelos para densidad básica, contracción volumétrica y contracción radial alcanzaron valores elevados de coeficientes de determinación de calibración mientras que los valores de coeficientes de determinación para la predicción de las mismas variables, resultaron medios para la densidad básica y contracción volumétrica.

Tabla 3:
Parámetros estadísticos de los modelos NIR desarrollados para las propiedades físicas de Guazuma crinita de 5 y 8 años.

Los parámetros estadísticos de los modelos para densidad básica y contracción volumétrica obtenidos por NIR permiten una adecuada predicción. El para densidad básica (0,88) es mayor al reportado por Gouvea et al. (2011), quienes encontraron valores de de 0,84 y mayor a los encontrados por Schimleck et al. (2006) cuyo está en el rango de 0,64 a 0,76; así mismo, está dentro del rango obtenido por Reis-Milagres et al. (2013), quienes obtuvieron un de 0,80 a 0,98. El para contracción volumétrica (0,83) también resulta consistente para la predicción. El RPD_p para densidad básica de 1,55 es satisfactorio según Schimleck et al. (2005) quienes señalan que valores iguales o superiores a 1,5 son suficientes para la predicción. Lo mismo ocurre con el RPD_p de la contracción volumétrica (3,27); en tanto, Reis-Milagres et al. (2013) descartan que modelos con valores de RPD_p por debajo de 1 sean buenos en la predicción; agregando Fujimoto et al. (2008) que mientras los valores de RPD_p son más altos, la calibración es más confiable. Para contracción tangencial y radial los valores de son muy bajos por lo que la predicción no es útil; debido a que la data tomada para NIR, se realizó en 10 posiciones longitudinales sobre la cara transversal; en tanto, la medición de las contracciones se realizó en 3 posiciones entre la médula y la corteza: 0-33%, 33-66%, 66-100%; que sumado a la alta variabilidad en las contracciones radiales y tangenciales que superan el 20%, podría explicar la baja correlación en los resultados y corroborado por Costa et al. (2018) quien señala que los planos de corte radial y tangencial son los de mayor variabilidad.

Acorde con los resultados obtenidos por ensayos destructivos y NIR se concluye que:

La densidad básica promedio obtenida para las edades de 5 y 8 años fue de 360 y 400 kg/m³ respectivamente, lo que permite clasificarla como madera de baja densidad.

La densidad básica, contracción volumétrica, tangencial y radial son mayores para la edad de 8 años y axialmente disminuyen de la base al ápice del árbol para ambas edades con variaciones estadísticas significativas.

Buenas correlaciones por NIR fueron obtenidas para los modelos desarrollados para densidad básica y contracción volumétrica para ambas edades; en tanto, que la contracción tangencial y radial no proporciona buena predicción.