Introducción
El suelo brasileño y las condiciones climáticas, en combinación con incentivos fiscales otorgados en la década de 1960, así como la inversión en investigación y tecnología han permitido a la silvicultura nacional brasileña de los géneros Eucalyptus y Pinus alcanzar altos niveles de productividad, convirtiéndose en un modelo de éxito mundial. Recientemente, el interés en la información sobre propagación y desempeño silvícola de especies nativas ha aumentado debido a la necesidad de restaurar o promover la regeneración natural de 18 millones de hectáreas para 2030, meta asumida por Brasil en el el Glasgow en la 26a Conferencia de las Naciones Unidas sobre el Cambio Climático (Ministério do Meio Ambiente, 2021). Los costos de la restauración forestal son altos y para aumentar el éxito de las plantaciones se requiere una base de datos sobre el desempeño de las especies nativas en condiciones de campo bajo diferentes condiciones edafoclimáticas, seleccionando las especies apropiadas (Morais Júnior et al., 2019). Esta base de datos debe incluir los requerimientos para el buen desarrollo de cada especie. En el campo, las plántulas recién plantadas pueden beneficiarse de la fertilización al promover el crecimiento temprano y aumentar las tasas de supervivencia donde la baja fertilidad del suelo y la alta competencia herbácea pueden inhibir el crecimiento (Earnshaw, Baribault y Jacobs, 2016). La dosificación de fertilizantes y la formulación química influyen en su efectividad. El crecimiento inicial de los árboles y la producción de biomasa están más comúnmente limitados por el fósforo (P) (Netzer, Pozzi, Dubbert y Herschbach, 2019), especialmente en suelos degradados en regiones tropicales.
Las especies forestales nativas mostraron una gran variación en el comportamiento de la fertilización con fosfato. Se ha demostrado que las especies pioneras responden a la fertilización y, con el avance del grupo sucesional, el estímulo al crecimiento proporcionado por la fertilización es menos pronunciado y a veces inexistente (Santos, Resende, Neto y Corte, 2008). Se encontraron resultados positivos de la adición de P para Senna macranthera (Fonseca, Paiva y Neves, 2011), Anadenanthera macrocarpa (Gonçalves, Paiva, Neves y Gomes, 2012), Mimosa caesalpiniaefolia (Gonçalves, Paiva, Lima-Neves y Gomes, 2013), Dalbergia nigra (Gonçalves, Paiva, Neves, Klippel y Caldeira, 2014), Schinopsis brasiliensis (Oliveira, Teixeira-Rios, Melo y Yano-Melo, 2015) y Cassia grandis (Freitas, Paiva, Leite y Oliveira Neto, 2017). Souza et al. (2012) observaron que las plántulas de Peltophorum dubium presentaban altos requerimientos nutricionales para nitrógeno (N), utilizando P de manera más eficiente. Sin embargo, Santos, Resende, Neto y Corte (2008) encontraron que el P tenía poca influencia en el crecimiento y la nutrición de las plántulas de Calophyllum brasiliensis y Myroxylon peruiferum.
La mayoría de los estudios sobre el efecto de la adición de P en especies forestales nativas se han realizado en condiciones de vivero, descuidando el objetivo principal de la respuesta de campo. Martinotto, Martinotto, Coelho, Azevedo y Figueiredo (2012), evaluando la supervivencia y el crecimiento inicial de seis especies del bosque de Cerrado en consorcio con yuca (Manihot esculenta), encontraron que la fertilización con fosfato favoreció solo a la especie Sclerolobium paniculatum. Por lo tanto, para repetir con especies nativas el éxito alcanzado con especies exóticas, es necesario dominar las etapas de crecimiento y los requisitos de las especies, incluyendo fertilización, encalado, espaciado, manejo y otros.
Vinhático, vinhático-do-campo, candeia e pau-amarelo son algunos nombres populares de Plathymenia reticulata Benth (Fabaceae-Mimosoideae), especie arbórea nativa que crece en una amplia área de distribución natural en el Cerrado, el Bosque Atlántico, Caatinga, Amazon y otros biomas de América del Sur (Carvalho, 2009; Morim, 2015). Descrita como pionera (Morais Junior et al., 2019) o secundaria temprana, esta especie puede presentar individuos adultos de 5 m a 30 m de altura con diámetro que varía de 30 cm a 150 cm, dependiendo del bioma (Carvalho, 2009).
El crecimiento de P. reticulata se describe como relativamente rápido, ocurre en suelos de alta y baja fertilidad, tiene potencial para su uso en la recuperación de áreas degradadas, forestación urbana y sector maderero (Lorenzi, 2008; Carvalho, 2009). P. reticulata se usa ampliamente en la medicina popular para el tratamiento de infecciones, hemorragias e intoxicaciones por serpientes (Moura et al., 2016). La versatilidad de P. reticulata la hace interesante para su uso en programas de restauración forestal y plantaciones comerciales, sin embargo, los trabajos que evalúan su crecimiento en el campo aún son incipientes.
Objetivos
El objetivo de este estudio fue evaluar el efecto de diferentes dosis de P en el rendimiento silvícola de P. reticulata en un área del Bosque Atlántico en Fazenda São Marcos en Sooretama, norte del estado de Espírito Santo, Brasil. Para esto, se probaron las siguientes hipótesis: (i) P. reticulata tiene potencial para su uso en programas de reforestación en el Bosque Atlántico, presentando un rápido crecimiento inicial; (ii) una mayor disponibilidad de P en el sustrato mejora el crecimiento inicial y a mediano plazo de un rodal monoespecífico de P. reticulata.
Materiales y métodos
Área de estudio
El experimento se estableció en un área de 1.31 hectáreas del proyecto Biomas, ubicado en Fazenda São Marcos en Sooretama, norte del estado de Espírito Santo, Brasil, desde abril de 2013 hasta abril de 2016, en las coordenadas 19 ° 12'55 "S y 40 ° 03 ' 09” W y 43 m de altitud. El sitio del experimento se inserta en el bioma del Bosque Atlántico con vegetación primaria clasificada como bosque de ombrófilos densos de tierras bajas. El clima de la región es cálido y semihúmedo Aw (Köppen), con una temperatura media anual de 23.5 °C y una precipitación anual de 1276 mm (Alvares et al., 2013) (Fig. 1).
Efecto de P sobre el rendimiento silvícola de Plathymenia reticulata
El experimento se realizó en un diseño de bloques al azar con cinco dosis de P y tres repeticiones/bloques. Los tratamientos consistieron en dosis de 0 g, 27 g, 47 g, 67 g y 87 g de P2O5 por cada árbol. Cada bloque se subdividió en 5 parcelas de 12 m × 33 m (396 m² por parcela). Cada tratamiento consistió en un conjunto de árboles dispuestos en cuatro filas que contenían 11 plantas cada una, sin embargo, para las mediciones se consideró un borde simple, siendo la parcela útil las dos líneas centrales con nueve árboles, con un total de 18 árboles útiles por parcela.
El suelo del área se caracteriza por Argisol amarillo distrófico (Sistema Brasileiro de Classificação de Solos [SiBCS] 2013). Antes de plantar, el área estaba cubierta de pasto con Pennisetum clandestinum. El control de las plantas espontáneas se llevó a cabo aplicando herbicida de glifosato en toda el área del experimento para facilitar las actividades de plantación y realización del experimento. Previo al experimento, se hizo una caracterización de atributos químicos del suelo (Tabla 1), de acuerdo con Donagema, Campos, Calderano, Teixeira y Viana (2011).
Prof. | pH | P* | K | Ca + Mg | Al 3+ | H + Al | SB | t | T | m | V | M.O |
(cm) | CaCl2 | mg dm-3 | mg dm-3 | -------------------- | cmolc dm-3 | ----------------------- | ----- % ---- | g dm-3 | ||||
0-12 | 4.6 | 1.5 | 0.1 | 1.0 | 0.2 | 2.9 | 1.7 | 1.9 | 4.7 | 9.4 | 37.3 | 19.9 |
12-30 | 4.6 | 0.4 | 0.1 | 1.8 | 0.2 | 3.1 | 1.9 | 2.1 | 4.9 | 11.3 | 37.3 | 11.3 |
30-60 | 4.6 | 0.3 | 0.1 | 1.5 | 0.3 | 3.1 | 1.9 | 1.9 | 4.9 | 17.8 | 34.2 | 17.8 |
* Extraído por el método Melhich-1.
La corrección del suelo para la siembra consistió en la aplicación manual de 1.0 Mg ha-1 de piedra caliza dolomítica. Posteriormente, el suelo sebarbechó con un subsolador-descompactador a una profundidad de 60 cm para facilitar el crecimiento del sistema de raíces de la planta. El espaciamiento entre plántulas adoptado fue de 3 m × 3 m y para su implantación se excavaron hoyos o pocillos de 30 cm × 30 cm × 30 cm.
Se usó superfosfato simple como fuente de P. Se agregaron al suelo 5 g de sulfato de zinc por pocillo, junto con la fertilización con P. Las plántulas se compraron en el vivero Vale Natural Reserve en Sooretama y tenían una altura promedio 30 cm. Las plántulas se plantaron tres meses después de la aplicación de dosis de P en los pozos preexcavados, agregando polímero de hidrocloruro hidratado. Debido a las condiciones climáticas del sitio durante el período de siembra, fue necesario colocar mantillo alrededor de las plántulas para retener la humedad y suavizar la temperatura del suelo y, en consecuencia, la evapotranspiración.
Se aplicaron nitrógeno y potasio en dosis de 20 g por plántula, de la siguiente manera: 35% de la dosis total tres meses después de la siembra; 35% seis meses después de plantar y 30% un año después de plantar. Además de la primera fertilización de cobertura, también se aplicaron 10 g de bórax por plántula, por medio de hoyuelos que la rodeaban.
A los 3 meses, 9 meses y 15 meses de la siembra, se midieron las siguientes variables: altura total (H, m) desde el suelo hasta el final del brote apical y diámetro a la altura del suelo (DAS, mm); esta última medición se hizo a una altura de 7 cm. La primera medición de H, a los 3 meses, se realizó con cinta métrica y las dos mediciones posteriores se hicieron con regla telescópica; el DAS se midió con un calibrador digital. A los 36 meses, se midieron la altura total (H, cm), el diámetro a la altura del pecho (DAP, mm; medido a 1.3 m sobre el nivel del suelo) con la ayuda de una cinta métrica; el área basal (G, m2 ha-1), mediante la suma de áreas seccionales a 1.3 m, extrapoladas a hectáreas; y supervivencia (%), contando el número de plántulas vivas en los tiempos de evaluación en función del número de plántulas plantadas.
Para evaluar el estado nutricional de los árboles de P. reticulata en diferentes tratamientos, se recolectaron hojas a los 36 meses. La colección se zigzague en la gráfica útil y las hojas se recolectaron en cada punto cardinal en la parte intermedia de la corona. Las muestras de hojas por parcela formaron una muestra compuesta. Las muestras se tomaron para su secado (en un horno de circulación forzada de aire a 60 °C hasta masa constante), molienda (en el molino Willey) y almacenamiento (en bolsas de papel) y luego se enviaron al laboratorio para análisis químico, de acuerdo con la metodología descrita por Silva (2009).
Se recolectaron muestras simples de hojarasca acumulada, totalizando 10 muestras simples por parcela, en los mismos tiempos y puntos de recolección de hojas. La operación se realizó con la ayuda de una plantilla metálica de 0.25 m × 0.25 m. Después de secar (en un horno de circulación forzada de aire a 60 °C hasta una masa constante) y medir la masa seca, las 10 muestras de cada parcela se unieron para formar una muestra compuesta. Luego, la muestra fue triturada en un molino tipo Willey, tamizada en una malla de 1.0 mm de abertura, almacenada en bolsas de papel y enviada para análisis químico (Tabla 2), según Silva (2009). El suelo se caracterizó por atributos químicos al final del experimento a los 36 meses (Tabla 3), de acuerdo con Donagema et al. (2011).
P 2 O 5 | N | P | K | Ca | Mg | S |
Gramos por árbol | --------------------------------- | g kg-1---------------------------------- | ||||
0 | 14.96 | 0.39 | 0.98 | 9.13 | 1.50 | 1.53 |
27 | 15.96 | 0.43 | 1.04 | 9.69 | 1.77 | 1.40 |
47 | 16.61 | 0.55 | 1.04 | 9.55 | 1.71 | 1.64 |
67 | 15.56 | 0.53 | 1.59 | 7.11 | 1.54 | 2.09 |
87 | 12.60 | 0.54 | 1.69 | 6.70 | 1.54 | 1.43 |
Prof. | pH | P* | K | Ca + Mg | Al 3+ | H + Al | SB | t | T | m | V | M.O |
(cm) | H2O | mg dm-3 | mg kg-3 | -------------------- | cmolc dm-3 | ----------------------- | ----- % ---- | g dm-3 | ||||
0-5 | 6.19 | 2.6 | 74.2 | 2.38 | 0.00 | 1.61 | 2.71 | 2.71 | 4.32 | 0.00 | 62.57 | 2.17 |
5-10 | 5.98 | 2.4 | 63.8 | 1.88 | 0.00 | 1.95 | 2.17 | 2.18 | 4.12 | 0.00 | 52.8 | 2.12 |
10-20 | 5.68 | 2.13 | 52.53 | 1.38 | 0.08 | 2.43 | 1.62 | 1.70 | 4.05 | 5.01 | 39.83 | 2.00 |
* Extraído por el método Melhich-1
Los datos se sometieron a análisis de regresión para verificar la dosis óptima de P para cada variable, utilizando la primera derivada de los estimadores β0 y β1. Para elegir la ecuación, se consideró la importancia de los parámetros, la prueba t, la importancia biológica y el coeficiente de determinación (R2). El modelo adoptado fue el lineal y = β0 + β 1x . Los análisis se realizaron con el programa R (R Core Team, 2018) versión 3.5.0.
Modelos analizados y estructuras de varianza y autocorrelación
Se ajustó el modelo logístico comúnmente utilizado para modelar el crecimiento en altura de árboles individuales (Ec. 1).
Donde:
h = altura (m)
id = edad actual (meses)
β i= coeficientes del modelo
Como el crecimiento en altura varía de un árbol a otro, el modelo logístico se ajustó para efectos fijos y mixtos (efectos fijos y aleatorios), para ver si la inclusión del efecto aleatorio (árbol) daría lugar a ganancias en términos de precisión en la estimación de la altura futura de los árboles de P. reticulata. El modelo logístico se ajustó en la forma mixta de la siguiente manera:
Donde:
h ij = altura (m) del i-ésimo árbol en la j-ésima medida
id = edad del j-ésimo árbol en la i-ésima medida
β0, β1, β2 = efectos fijos de los parámetros
e ij = error aleatorio, que se asume que se distribuye como eij ~ N (0, Ri)
ui = efectos aleatorios de los parámetros con una matriz de varianza-covarianza (D), es decir, ui ~ N (0, D)
La matriz de varianza-covarianza (Ri) se puede determinar como:
Donde:
G i = matriz diagonal para describir la heterocedasticidad del grupo de efectos aleatorios
En los modelos de crecimiento de árboles individuales, la variación de las alturas generalmente aumenta a medida que los árboles crecen, este efecto se conoce como heterocedasticidad (Pinheiro y Bastes, 2000; Campos y Leite, 2013). En este caso, para corregir dicho efecto, se analizaron las funciones varPower (potencia de una covariable) y varIdent (varianza constante por grupo). Además, para corregir una eventual autocorrelación entre los datos observados de medidas repetidas a lo largo del tiempo en el mismo individuo, se probaron dos estructuras de autocorrelación: autorregresiva de primer orden (AR1) y autorregresiva contínua de primer orden (CAR1).
El modelo logístico en la forma fija se ajustó mediante el método Gauss Newton utilizando el software R Core Team (2018), utilizando la función nls. Para ajustar los modelos no lineales mixtos, se utilizó el algoritmo de máxima verosimilitud restringido, disponible en el paquete nlme (Pinheiro y Bates, 2000).
Estimación de la precisión de las estimaciones
La evaluación de las estimaciones generadas por los modelos mixtos y el modelo fijo se realizó con base en las siguientes estadísticas: sesgo del coeficiente de correlación (r), parcialidad [P (%)], raíz de la desviación cuadrática media [RMSE (%)], criterio de información de Akaike (AIC) y criterio bayesiano (BIC) (Tabla 4). Para complementar las estadísticas, se prepararon gráficos que relacionaban los valores observados y los estimados por los modelos.
Resultados
A los 3 meses, los resultados mostraron un efecto negativo de las dosis de P sobre el crecimiento, tanto en H como en DAS de plántulas de P. reticulata. El crecimiento máximo observado, 76.4 cm en H y 9.9 mm en DAS se obtuvo para la dosis de 0 g de P2O5 (Fig. 2a, b) y, con cada aumento de dosis de P, se observó una reducción promedio de 6.6 cm en H y de 0.9 mm en DAS.
*y ** significativo al nivel de 5% y 1% de probabilidad por la prueba F respectivamente
ns no significativo a 5% de probabilidad, por la prueba t.
Se observó la misma tendencia a los 9 meses de crecimiento, el aumento de las dosis de P redujo linealmente el incremento de H y DAS de las plántulas de P. reticulata (Fig.2c, d). Los valores máximos de H y DAS fueron 2.40 m y 3.47 mm, respetivamente, observados en la dosis de 0 g de P2O5, con una reducción promedio de 31.6 cm y 3.8 mm en cada aumento de dosis de P.
A los 15 meses de crecimiento, el soporte de P. reticulata todavía tuvo un efecto negativo de la fertilización con fosfato, la mayor altura observada fue 4.16 m, se observó para la dosis de 0 g de P2O5, con una tasa de crecimiento de reducción en el orden de 20.4 cm por cada aumento de dosis de P (Fig.2e). También se encontró que a los 15 meses de crecimiento las diferencias presentadas para DAS en función de la adición de P tendían a desaparecer, presentando un promedio general de 5.42 cm (Fig.2f).
Con respecto a los 36 meses de crecimiento, no hubo efecto de la adición de P sobre el crecimiento del árbol de P. reticulata, el rodal presentó un promedio de H de 5.9 m y DAP de 9.51 cm (Fig. 3a, b). Sin embargo, se encontró que la fertilización con fosfato promovió una reducción lineal significativa en la supervivencia de los árboles, lo que refleja la reducción del área basal. Se observaron los mejores resultados para la supervivencia (92.92%) y el área basal (7.34 m2 ha-1) para la dosis de 0 g de P2O5 (Fig. 3c, d).
*y ** significativo al nivel de 5% y 1% de probabilidad por la prueba F respectivamente
ns no significativo a 5% de probabilidad, por la prueba t.
No se observó ningún efecto de la adición de P en la plantación sobre la nutrición de los árboles de P. reticulata a los 36 meses de crecimiento, excepto por el contenido de P de la hoja, que mostró una tendencia de crecimiento lineal creciente con el mejor resultado siendo 0.82 g kg-1, a 87 g de P2O5 (Fig. 4b). En general, los árboles de P. reticulata presentaron el siguiente orden de absorción N > Ca > K > Mg > S > P.
*y ** significativo al nivel de 5% y 1% de probabilidad por la prueba F respectivamente
ns no significativo a 5% de probabilidad, por la prueba t.
Los parámetros estimados del modelo logístico de efectos fijos para la estimación de altura futura, así como los errores estándar respectivos y las pruebas de significación se pueden observar en la tabla 5. Se puede observar que todos los parámetros del modelo fueron significativos (p <0.01) por la prueba t.
Estimar | Error estándar | tc | Pr>|tc| |
6.48689 | 0.86446 | 7.504 | <0.001 |
21.5194 | 2.11995 | 10.151 | <0.001 |
0.23403 | 0.02098 | 11.154 | <0.001 |
La tabla 6 muestra los indicadores estadísticos utilizados para evaluar la calidad de ajuste de los modelos de efectos mixtos y su forma original. Los ajustes de modelos mixtos con diferentes estructuras de varianza-covarianza y autocorrelación dieron como resultado estimaciones RMSE (%) que van desde 20.58 a 36. Se observó sesgo y los coeficientes de correlación fueron mayores a 0.88.
Modelo | F. variación | F. autocorrelación | RMSE (%) | P (%) | r | AIC | BIC |
1 | - | - | 34.86 | 0.380 | 0.88 | 1012.08 | 1028.82 |
2 | varIdent | AR1 | 20.58 | -0.26 | 0.96 | 868.47 | 893.54 |
3 | varIdent | CAR1 | 20.07 | -0.24 | 0.96 | 874.61 | 912.2118 |
4 | varPower | AR1 | 36.02 | -6.05 | 0.88 | 1074.61 | 1012.2118 |
5 | varPower | CAR1 | 36.02 | -6.05 | 0.88 | 1074.61 | 1012.2118 |
En general, los modelos mixtos (modelos 2 y 3) obtuvieron una mayor precisión en el ajuste y la homogeneidad de los valores estadísticos en comparación con el modelo logístico en su forma fija (Modelo 1). Sin embargo, los modelos 4 y 5 dieron como resultado valores más discrepantes en comparación con otros modelos, considerando la RMSE (%). Con respecto a las estadísticas de las estimaciones que utilizan los modelos mixtos, los ajustes que utilizan la función varIdent fueron superiores en comparación con la función varPower.
Al observar los gráficos de la relación entre la altura observada y la estimada (Fig. 4), es evidente que el modelo logístico ajustado solo con efecto fijo (Fig. 5b) presentó una mayor dispersión en relación con la línea ideal (1:1), que representa menor precisión. En el modelo mixto ajustado (Fig. 5 a), los puntos fueron más agregados y uniformes, lo que corrobora los resultados de los indicadores estadísticos que se muestran en la tabla 6.
Discusión
En general, se confirma la primera hipótesis planteada, P. reticulata mostró una alta supervivencia y un rápido crecimiento inicial, lo que indica que la especie tiene un alto potencial para su uso en programas de reforestación en el Bosque Atlántico. También hubo una tendencia de mejora en los atributos químicos del suelo después de 36 meses de asentamiento, con énfasis en los aumentos de P, Ca + Mg, suma base, CTC efectivo, saturación de bases y materia orgánica (Tabla 3), indicando el comienzo de los procesos de ciclo de nutrientes.
Hay pocas plantaciones de la especie a nivel experimental, comercial y de recomposición, lo que dificulta una comparación más segura. Morais Junior et al. (2019), con el objetivo de identificar el potencial de 25 especies de árboles que compensan el carbono en un área degradada del Bosque Atlántico, en el estado de Minas Gerais, descubrieron que Plathymenia foliolosa Benth (sinónimo de P. reticulata) tenía un alto potencial para tales plantaciones debido a su alta tasa de supervivencia (75%), su crecimiento inicial en H (5.34 m), DAS (13.43 cm) y su incremento promedio de carbono (29.1 kg año-1) a los 29 meses. Los resultados de crecimiento existentes de P. reticulata son similares a los obtenidos en este estudio.
Sin embargo, el rendimiento silvícola de P. reticulata fue influenciado negativamente por la adición de P, lo que contradice la segunda hipótesis, ya que las especies de leguminosas de la sucesión ecológica temprana generalmente presentan una alta demanda nutricional. Aunque las diferencias en H y DAS fueron más prominentes en los primeros meses de implantación y desaparecieron a los 36 meses, el aumento en la fertilización con fosfato promovió una reducción significativa en el número de individuos sobrevivientes, lo que refleja la disminución en el área basal baja.
Millones de años de evolución han ayudado a las plantas a desarrollar una serie de mecanismos genéticos, morfológicos, fisiológicos y bioquímicos de adaptación a los hábitats pobres de P. Estos mecanismos dan como resultado una mayor adquisición de P y un uso interno más eficiente (Rao, Miles, Beebe y Horst, 2016) e incluyen la modificación de la arquitectura de la raíz, el desarrollo de un sistema radical extenso, raíces más largas y raíces más delgadas, acidificación de rizosfera, exudación de ácidos orgánicos, expresión de transportadores de P de alta afinidad, eficiencia de P de la homeostasis citoplasmática y removilización interna de P, secreción de fosfatasa y simbiosis micorrízica (Makoi, Chimphango y Dakora, 2010; Balemi y Negisho, 2012; Folli-Pereira, Meira-Haddad, Bazzolli y Kasuya, 2012). Por lo tanto, se sugiere que P. reticulata tiene una alta adaptabilidad a las condiciones restrictivas de P, ya que originalmente el suelo tenía en promedio 0.73 mg dm-3 de P en las capas analizadas (Tabla 1), un valor que se considera muy bajo. Los niveles de P de las hojas en los árboles a menudo se usan como un indicador del suministro de P (Fäth, Kohlpaintner, Göttlein y Mellert, 2019), sin embargo, la sutil tendencia al alza en las hojas de P. reticulata encontradas a los 36 meses en función del aumento de la fertilización con fosfato no fue acompañada por ningún efecto positivo del crecimiento del asentamiento. Esto es aún más evidente cuando se observan los niveles de P que se encuentran en la camada (Tabla 2), que presentó el más bajo entre los macronutrientes, lo que indica que la especie tiene una baja demanda o alta eficiencia de utilización de este nutriente.
Silva et al. (2006), evaluando la eficiencia en el uso de nutrientes por tres especies de árboles pioneros que se encuentran con frecuencia en pasturas degradadas donde el P siempre es limitante, después de encalado y la aplicación de P, verificaron diferencias entre las especies estudiadas. Vismia japurensis utilizó la mayoría de los nutrientes analizados de manera más eficiente, especialmente N, P y Mg, con altos valores de tasa fotosintética, bajos contenidos foliares de P, y la limitación de estos nutrientes no fue limitante para el buen desempeño de esta especie, como es el caso de P. reticulata en este estudio. Por lo tanto, a pesar de los notables avances en la silvicultura nativa, se necesitan estudios para establecer condiciones nutricionales óptimas para cada especie cultivada en suelos con baja disponibilidad nutricional.
Con respecto al modelo de crecimiento en altura para la especie estudiada en este trabajo, es evidente la ventaja del modelo logístico ajustado en forma mixta en comparación con el ajuste tradicional. Esto se debe al modelo mixto propuesto que presenta el árbol como un componente aleatorio. En este caso, considerar el árbol en el análisis significa asignar curvas de crecimiento a cada uno de los árboles estudiados, reduciendo así el error cometido. Es razonable decir que este es el factor principal para reducir la RMSE de 34.86% (modelo 1) a 20.07% (modelo 3).
Una de las ventajas de los modelos que incluyen efectos aleatorios en su estructura en comparación con los modelos de efectos fijos es la reducción del error estándar residual. Obras como Mendonça, Pádua Chaves y Calegario (2015) obtuvieron una ganancia de precisión de 53.89% en modelos hipsométricos; Carvalho et al. (2011), tuvieron una reducción en el error estándar residual (Syx) de 26.31% a 4.01% en predicción volumétrica. Gouveia, Silva, Ferreira, Gadelha y Lima-Filho (2015), también trabajando con modelos volumétricos, obtuvieron una reducción en el error estándar residual de 26 veces.
Los modelos de regresión tradicionales requieren supuestos de cumplimiento como la independencia entre las observaciones y la homogeneidad de la varianza, que no siempre se satisfacen (Gouveia et al., 2015). Una alternativa a estos modelos es la sustitución por modelos mixtos, que permiten generalizar las estructuras de correlación espacio-temporal y la varianza no constante, extrayendo mejor las relaciones de variables independientes con datos agrupados en comparación con los modelos de regresión tradicionales (Fausto, Carneiro, Antunes, Pinto y Colosimo, 2008; Gouveia et al., 2015). Esto demuestra una vez más la ventaja de este tipo de modelos aplicados a los datos con medidas repetidas en el tiempo para evaluar el crecimiento de las especies forestales, como se estudia en esta investigación.
Conclusiones
La especie Plathymenia reticulata tiene potencial para su uso en programas de reforestación bajo condiciones edafoclimáticas probadas y condiciones similares.
La especie no es exigente en P, la adición de este nutriente promueve efectos negativos en su crecimiento a corto y mediano plazos.
Con respecto al modelado de crecimiento en altura, el modelo logístico ajustado en la forma mixta fue superior al tradicional.