Ubicación y caracterización del área de estudio

La investigación se realizó entre enero de 2010 y diciembre de 2017, en la Reserva Nacional Tambopata (RNT) del Centro de Investigaciones de Tambopata (Tambopata Research Center/TRC), ubicada a 250 msnm, en las coordenadas 433162 E - 8548037 N, cerca al límite del Parque Nacional Bahuaja Sonene, región Madre de Dios, Perú (Figura 1). El área de estudio abarca una superficie de 1 613 ha, con una precipitación media anual de 2 925 mm y una temperatura media anual de 25 °C.

Figura 1 Ubicación del área de estudio y distribución de parcelas en la RNT, Madre de Dios, Perú. 

^{Ribeiro et al. (1999)} clasificó el área en cuatro tipos de bosques: Bosque de Aguajal (Ba), con suelos inundados durante casi todo el año y predominando Mauritia flexuosa L.f. (aguaje); Bosque de Bajío (Bb), con suelos ocasionalmente inundados y con pocas especies emergentes como Dipteryx micrantha Harms (shihuahuaco) y Ceiba pentandra (L.) Gaertn. (ceiba); Bosque Sucesional (Bs), situado en la llanura del río y dominado por taxones de rápido crecimiento como Erythrina poeppigiana (Walp.) O.F.Cook y E. ulei Harms, Triplaris americana L. (tangarana) y parches de Guadua sp. (paca); y Bosque de Tierra Firme (Bt), con suelos arcillosos, con especies emergentes como Hymenaea courbaril L. (azúcar huayo) y Eschweilera coriacea (DC.) S.A. Mori.

Colecta de datos y procesamiento de muestras

Las parcelas fueron instaladas siguiendo la orientación del sistema de trochas del Centro de Investigación. Se evaluaron 20 parcelas de 50 × 30 m, cinco para cada tipo de bosque, que representan 3.5 veces el área de 20 × 20 m recomendada para una subparcela (^{Phillips et al., 2016}). En cada parcela se registraron todos los individuos arbóreos con DAP (diámetro a 1.30 m del suelo) ≥ 10 cm. Los individuos registrados (n= 1 958) se colectaron, herborizaron e identificaron a partir de visitas virtuales a colecciones botánicas: Trópicos (https://www.tropicos.org/), The Plant List (http://www.theplantlist.org/) y GBIF (https://www.gbif.org/); Así como de consultas a especialistas y visitas a herbarios físicos próximos al área de estudio: Herbario Alwyn Gentry (HAG), herbario Vargas (CUZ) y herbario Selva Central de Oxapampa (HOXA).

Posteriormente, se ordenaron por familias de acuerdo al sistema propuesto por el grupo para la filogenia de Angiospermas (^{APG IV, 2016}); los ejemplares se depositaron en el Herbario de la RNT. Enseguida, se procedió a identificar taxones indicadores por tipo de bosque, a partir de la fidelidad y afinidad de la especie por el tipo de hábitat (^{Dufrêne y Legendre, 1997}), cuyos individuos se monitorearon mensualmente con binoculares para evaluar las siguientes fenofases: botón floral, flor, fruto inmaduro y fruto maduro. En la toma de datos fenológicos, se utilizó una modificación de la escala de ^{Fournier (1974)} para ausencia/presencia de fenofases, en la que 0 indica ausencia y 1 presencia. Los datos de precipitación y temperatura para el periodo de estudio se obtuvieron de la Estación Meteorológica de la RNT y la Estación Meteorológica de Puerto Maldonado (Madre de Dios).

Análisis de datos

Se construyeron tres matrices: una florística, una fenológica y una ambiental. A partir de la matriz florística, se llevó a cabo un análisis de ordenamiento NMDS (Non-metric Multidimensional Distance Scaling) (^{Kruskal, 1964}) que utiliza la matriz de disimilitud entre parcelas, aplicando la distancia de Bray-Curtis. Además, se aplicó un análisis PERMANOVA (Permutational Multivariate Analysis of Variance) para verificar diferencias significativas entre los grupos formados y se usó el método de rarefacción para comparar la riqueza de especies entre tipos de bosques (hábitats); se estimó el máximo número de especies a partir del índice de Chao2, como medida de eficiencia del muestreo (^{Colwell, 2013}). Para verificar las especies indicadoras de tipos de hábitats, se hizo un análisis ISA (Indicator Species Analysis) (^{Dufrêne y Legendre, 1997}), que determina la significancia de las especies en análisis a partir de una prueba de Monte Carlo para probar la hipótesis nula de que las especies en estudio no tenían un valor indicador (IndVal).

Posteriormente, para verificar la influencia de las variables ambientales sobre las fenológicas, se aplicó un análisis CCA (Canonical Correlation Analysis) (^{McCune y Grace, 2002}). Se obtuvo la relación de los promedios de precipitación y temperatura con los promedios de las fenofases para las especies indicadoras por tipo de bosque, para la media de cada mes entre los años de 2010-2017.Por último, se verificó la correlación de esos meses con las fenofases en estudio para las especies indicadoras. Se utilizó la prueba de Monte Carlo con 1 000 permutaciones para determinar la significancia de la ordenación. Todos los análisis se efectuaron con el ambiente estadístico R 3.6.3 (^{R Core Team, 2021}) y el programa Fitopac 2.1 (^{Shepherd, 2010}).

Composición florística

En los cuatro tipos de bosque se registraron 1 958 individuos, agrupados en 57 familias, 173 géneros y 300 especies. Las familias con mayor número de géneros fueron Fabaceae (19), Euphorbiaceae (12), Malvaceae (12), Moraceae (12) y Annonaceae (11). Las familias con más especies fueron Fabaceae (36), Moraceae (27) y Annonaceae (18); y las que agruparon más individuos fueron Arecaceae (607), Euphorbiaceae (144), Fabaceae (142) y Myristicaceae (122) (Cuadro 1). Además se registraron 68 especies para el Bosque de Aguajal, 142 para el Bosque de Bajío, 82 para el Bosque Sucesional y 158 para el Bosque de Tierra Firme; cabe señalar que una especie se observó en más de un tipo de bosque.

Los bosques de la RNT con mayor similitud fueron el Bosque de Tierra Firme y el Bosque Sucesional. Por otro lado, los más disímiles fueron el Bosque de Aguajal y el Bosque Sucesional. Tales diferencias se reflejan en el gráfico de la Figura 2 del NMDS con stress de Bray-Curtis = 0.08 y varianza significativa (p<0.001).

Ba= Bosque de Aguajal; Bb= Bosque de Bajío; Bt= Bosque de Tierra Firme; Bs= Bosque Sucesional. Cada punto representa una parcela de 50 × 30 m en la RNT, Madre de Dios.

Figura 2 NMDS de especies de los cuatro bosques en estudio. 

Asimismo, el método de rarefacción-extrapolación, basado en las muestras mostró diferencia entre la riqueza específica y la heterogeneidad ambiental de las comunidades, el cual permitió constatar mayor riqueza en el Bosque de Tierra Firme y de Bajío (Figura 3).

Figura 3 Curva de rarefacción de la riqueza estimada usando el índice de Chao2 para los bosques de la RNT, Madre de Dios. 

Análisis de especies indicadoras

El ISA identificó ocho especies como indicadoras para los respectivos bosques; las más destacadas: Mauritia flexuosa y Euterpe precatoria Mart., para el Bosque de Aguajal; Iriartea deltoidea Ruiz & Pav., para el Bosque de Bajío; Erytrhina ulei, para el Bosque sucesional y Quararibea malacocalyx A. Robyns & S. Nilsson, para el Bosque de Tierra Firme (Cuadro 3).

El análisis de CCA, basado en los registros de las fenofases de las especies indicadoras y en las variables ambientales de temperatura y precipitación demostró ser significativo (MC = 0.001); con la variable precipitación fuertemente correlacionada al eje 1 (-0.99); en cuanto a la variable temperatura al eje 2 (-0.84). En el gráfico que tiene 22 % de variación acumulada (Figura 4), se observa que ambas variables ambientales estuvieron muy correlacionadas con los meses de noviembre a marzo; precisamente, los meses con las mayores temperaturas y precipitaciones que corresponden al verano neotropical. Asimismo, las fenofases botón (2.65) y flor (2.55) presentaron los valores canónicos más altos para el eje 1; el fruto maduro (-1.61) y fruto inmaduro (1.52) registraron mayores valores para el eje 2 (Figura 5).

Círculos mayores = Meses con mayor influencia de la precipitación; Escala de colores =Influencia de la variable temperatura, donde variaciones de azul y verde, meses con menores temperaturas; amarillo y naranja, meses con mayores temperaturas.

Figura 4 CCA de las variables ambientales en los bosques evaluados en la RNT, Madre de Dios. 

Círculos mayores = Meses con mayor influencia de la precipitación; Escala de colores = Influencia de la variable temperatura, donde variaciones de azul y verde, meses con menores temperaturas; amarillo y naranja, meses con mayores temperaturas.

Figura 5 CCA de las fenofases en la RNT, Madre de Dios. 

El análisis de las fenofases de las especies indicadoras (Cuadro 4) evidenció que en el Bosque de Aguajal la fenofase botón floral se inició en octubre y alcanzó su máximo valor en noviembre (21 %), y los mayores valores de fruto maduro se alcanzaron en diciembre (16 %) y enero (15 %). En el Bosque de Bajío el valor más alto de la fenofase botón floral (17 %) se obtuvo en septiembre y de fruto maduro en diciembre (29 %). En el Bosque Sucesional el máximo valor de la fenofase botón floral (14 %) se registró en septiembre y de fruto maduro (23 %) en diciembre; en el Bosque de Tierra Firme el valor más alto de la fenofase botón floral se observó en septiembre (35 %), mientras que la de fruto maduro en enero (24 %).

Composición florística

Es importante conocer las dimensiones de la heterogeneidad florística para describir y comprender los patrones dinámicos de las especies y su interacción con el ambiente. El presente estudio evidencia que las especies indicadoras analizadas con frecuencia se citan en estudios florísticos locales (^{Baez y Garate, 2017}; ^{Dueñas y Garate, 2018}; ^{Alvarez-Montalván et al., 2021}). Las familias, como Fabaceae, Moraceae y Arecaceae se consignan como representativas de los bosques amazónicos del Perú (^{Vásquez et al., 2010}; ^{Ureta, 2015}, ^{Alvarez-Montalván et al., 2021}).

De la misma manera, el agrupamiento específico de individuos por tipo de bosque permite establecer similitudes cercanas entre hábitats. Es así que el análisis de agrupamiento NMDS evidenció que hábitats con mayores números de individuos arbóreos, como el Bosque de Bajío y Bosque de Tierra Firme, presentaban más similitud que los bosques con menor cantidad de individuos. Por otro lado, se ha indicado que el análisis tiene la desventaja de generar más de una respuesta, por lo que resulta necesario llevar a cabo procedimientos para determinar el menor valor de tensión (^{Rocha-Loredo et al., 2010}). En ese sentido, la confiabilidad del análisis realizado válido; puesto que, el valor de stress estuvo por debajo de 0.1.

Adicionalmente, la comparación directa de la riqueza observada por tipo de hábitats mostró cómo habitas no correlacionados por el análisis de agrupamiento presentaron más especies, como el caso de los Bosque de Bajío y el Bosque Sucesional. Sin embargo, ninguno de los tipos de bosques alcanzó la estabilidad en el número de taxones indicada por el coeficiente indicador de Chao2. En complemento, ^{Samaniego et al. (2015)} con el mismo análisis de rarefacción, no obtuvieron diferencias en la riqueza. En contraste, ^{Silva et al. (2016)} indicaron que la estabilidad de la curva de rarefacción en el Bosque Nacional Tapajós se observó, especialmente, para los bosques secundarios.

Análisis de las respuestas fenológicas de las especies indicadoras

Las relaciones entre los diferentes tipos de hábitats y los parámetros meteorológicos evidenciaron mayores variaciones de las respuestas fenológicas en el Bosque de Tierra Firme, lo que indicaría más sensibilidad a los cambios de temperatura y precipitación. En cambio, al correlacionar los parámetros ambientales en estudio para los Bosques de Aguajal, de Bajío y Sucesional estos serían influenciados más por la temperatura que la precipitación. La ausencia de precipitaciones o sequías extremas y la temperatura, dependiendo de los periodos de duración (^{Siegmund et al., 2016}), afectarían a los individuos arbóreos del Bosque de Tierra Firme que se encuentran en la Amazonía (^{Zhao et al., 2017}). Uno de los indicadores inmediatos es la pérdida de verdor de los bosques y la senescencia de hojas; lo que conlleva a la acumulación de material orgánico, descomposición y aumento de CO₂.Al respecto, en el estudio de ^{Pires et al. (2018)}, realizado en una reserva de Mata Atlántica de Río de Janeiro, se determinó una relación débil entre las condiciones climáticas y la dinámica de formación y maduración de frutos, la que se relacionó con la fenodinámica de la floración. Por el contrario, en el presente estudio, se observó que la presencia de frutos inmaduros y maduros disminuyó, de manera ostensible, conforme se aproximaban los meses de junio a septiembre, época en que la precipitación alcanzó sus menores valores.

El comportamiento fenológico de las plantas, como indicador del cambio climático, está siendo investigado y monitoreado en diversos lugares del mundo (^{Workie y Debella, 2018}). En estos estudios el comportamiento fenológico resultará determinante en el análisis e interpretación de los componentes climáticos en el tiempo y el espacio (^{Pires et al., 2018}), tal como se observa en el estudio que aquí se documenta. Es por ello que, al margen de otros factores ambientales no considerados, el efecto de la temperatura y precipitación resultaron significativos en el establecimiento de las fenofases de las especies evaluadas.

Por otra parte, los bosques tropicales se caracterizan por presentar una fenología con alta diversidad debido a que carecen de una estación fría bien definida que afecte y restringa su crecimiento. Por esto, es posible observar especies con flores y frutos todo el año, otras que florecen y fructifican una vez al año, o una sola vez en varios años (^{Pires et al., 2018}). De igual forma, la mayoría de las plantas leñosas tropicales producen nuevas hojas y flores de manera repentina y no continuamente, por lo que sus cambios fenológicos representan adaptaciones a factores bióticos y abióticos (^{Van Schaik et al., 2003}).

Es importante destacar que, aunque la contribución de los resultados del presente estudio es limitada, constituye una base de apoyo al conocimiento científico local, y refleja la necesidad de integrar mayores números de especies, con el fin de no discriminar posibles taxones que podrían ser fundamentales para el entendimiento de la fenodinámica del área de estudio. Por lo tanto, es imperativo la conservación de la reserva y sus hábitats.