Introducción

En ecología, una de las preguntas más importantes es: ¿cuáles son los factores que explican la alta diversidad de especies de plantas en los bosques tropicales? Este tema representa una de las principales preguntas que la ecología tropical debe responder en los próximos años (Pennisi, 2005); por ello, se han planteado varias hipótesis para resolver esta interrogante. Una de ellas predice una relación positiva entre la diversidad de especies de plantas y la precipitación (^{Gentry, 1982a}, 1988; ^{Wright, 1992}; ^{Specht & Specht, 1993}; ^{Clinebell, Phillips, Gentry, Stark, & Zuuring, 1995}; ^{Aplet, Hughes, & Vitousek, 1998}). Esto se ha podido confirmar tanto a escalas continentales como locales (^{Bjorholm, Svenning, Skov, & Baslev, 2005}; ^{Slik et al., 2009}; Eiserhardt, Svenning, Kissling, & ^{Balslev, 2011}). Otra hipótesis postula que la diversidad florística está relacionada con la fertilidad del suelo; sin embargo, los estudios muestran evidencias contradictorias, así que la diversidad puede tanto aumentar con la fertilidad del suelo (^{Gentry, 1988}; Clinebell et al., 1995; ^{Duivenvoorden, 1994}; ^{Tuomisto, Zuquim, & Cárdenas, 2014}) o disminuir (^{Huston, 1980}; ^{Tuomisto et al., 2002}; Slik et al., 2009; ^{Peñas-Claros et al., 2012}). De igual manera, puede presentarse un pico intermedio donde existe una correlación positiva en varias partes del gradiente edáfico, negativa en otras (^{Ashton, 1992}), o simplemente no existe una relación entre la fertilidad de los suelos y la riqueza de las comunidades de plantas (^{Quinto-Mosquera & Moreno-Hurtado, 2014}). Los que mantienen que hay una relación proponen que las diferencias se deben a la heterogeneidad de los suelos y la metodología de los muestreos. En el neotrópico, los trabajos que relacionan las variables edáficas con la composición florística de las comunidades de palmas se han centrado en la Amazonia y en el noroccidente de Panamá (Vormisto, Svenning, ^{Hall, & Balslev, 2004}; ^{Andersen, Turner, & Dalling, 2010}; ^{Kristiansen et al., 2012}; ^{Cámara-Leret, Tuomisto, Ruokolainen, Balslev, & Kristiansen, 2017}; ^{Olivares, Svenning, Va Bodegom, Valencia, & Balslev, 2017}; ^{Stevenson, Aldana, Cardena, & Negret, 2018}). Además, se ha determinado que, en los bosques tropicales, la presencia y ausencia de algunas especies de plantas está directamente relacionada con la disponibilidad o concentración de ciertos nutrientes en el suelo (^{John et al., 2007}).

Por ejemplo, se ha reportado que existen asociaciones pronunciadas de muchas especies de árboles (^{Condit, Engelbrecht, Pino, Pérez, & Turner, 2013}) y palmas (^{Cámara-Leret et al., 2017}) a diferentes concentraciones de fósforo, lo que sugiere un papel predominante de este nutriente en la distribución y abundancia de las especies de plantas tropicales (^{Vitousek, 1984}; ^{Swaine, 1996}; ^{Potts, Ashton, Kaufman, & Plotkin, 2002}). El Chocó biogeográfico es uno de los lugares con mayor diversidad de plantas del mundo (^{Gentry, 1982b}, 1986, ¹⁹⁹³; ^{Forero & Gentry, 1989}; ^{Faber-Langendoen & Gentry, 1991}; ^{Galeano, 2001}). Su gran diversidad se ha explicado por la alta precipitación en la región (que pueden alcanzar los 11 000 mm anuales en la parte central) y por las bajas concentraciones de nutrientes (como Ca, Mg, K, P, Na) (^{Gentry, 1986}; Faber-Langendoen & Gentry, 1991). En general, los suelos se caracterizan por ser muy pobres en nutrientes, son de extrema acidez, alta saturación de aluminio y de gran concentración de materia orgánica (Faber-Langendoen & Gentry, 1991; ^{Quinto-Mosquera & Moreno-Hurtado, 2016}). Sus alto niveles de acidez se deben a la alta precipitación y lixiviación de nutrientes como Ca, Mg, K, y Na (Quinto-Mosquera & Moreno-Hurtado, 2016).

Las palmas son un ícono de los bosques tropicales donde, a menudo, son muy diversas localmente hiperdominates (^{Henderson, 2002}; ^{Dransfield et al., 2008}; ⁵⁷; ^{Balslev, Bernal, & Fay, 2016}). Presentan muchas formas de crecimiento (Balslev et al., 2011) y juegan un rol determinante en el funcionamiento de los ecosistemas (^{Kahn & de Granville, 1992}; Henderson, 2002). En los bosques de tierras bajas del Chocó, las comunidades de palmas están representadas por un conjunto de especies de géneros comunes como Geonoma y Bactris que son los más diversos (^{Henderson, Galeano, & Bernal, 1995}; Galeano & Bernal, 2010; ^{Balslev, Pedersen, Navarrete, & Pintaud, 2015}). Estas comunidades de palmas son reconocidas por su alta diversidad y densidad de especies (Tabla 1).

Adicionalmente, tienen multiplicidad de usos en esta zona (^{Ledezma-Renteria & Galeano, 2014}; Cámara-Leret, Copete, ^{Balslev, Soto-Gómez, & Macia, 2016}; ^{Schneider, Cámara-Leret, Barfod, & Weckerle, 2017}) y ha sido catalogada como la segunda familia más usada después de Lauraceae (^{Galeano, 2000}). Por estas razones, las palmas son un excelente grupo para estudiar cómo influyen los factores edáficos sobre la composición florística en la localidad antes descrita.

En las comunidades de palmas, las especies de sotobosque suelen ser más diversas y abundantes que las del dosel, patrón que algunas comunidades de árboles también siguen (^{Gentry & Dodson, 1987}; ^{Duivenvoorden, 1994}; ^{Balslev et al., 2011}). Sin embargo, son escasos los estudios ecológicos que han analizado las preferencias edáficas de ambas formas de crecimiento (^{Duque, Sánchez, Cavelier, & Duivenvoorden, 2002}; ^{Costa, Jean-Louis, Albertina, & Ocirio, 2009}; ^{Kristiansen et al., 2012}; ^{Cámara-Leret et al., 2017}). En el Chocó, pocas investigaciones han evaluado la influencia de los nutrientes sobre la composición florística. Tampoco existen estudios a nivel de comunidad de palmas o de especies individuales. Entre los trabajos antes realizados, se destaca el de Quinto-Mosquero & Moreno-Hurtado (2014), quienes realizaron un inventario florístico en dos parcelas permanentes de una hectárea en Salero (Chocó central) y encontraron que la diversidad de flores está positivamente correlacionada con la concentración de aluminio (Al) y negativamente con la topografía.

En este estudio se exploró la relación entre la composición florística de comunidades de palmas y la concentración de nutrientes en el suelo. Para ello, utilizamos datos florísticos y edáficos recolectados en 20 transectos distribuidos en el Chocó norte (Colombia) y Chocó sur (Ecuador). Específicamente, se analizó (1) cuánto varía la concentración de los nutrientes del suelo entre las dos localidades; (2) cuál es la relación entre los nutrientes del suelo y la composición florística de palmas; (3) cuáles nutrientes explican mejor la variación en la composición florística de las comunidades de palmas; y (4) cómo varía la abundancia de las especies más dominantes de sotobosque y de dosel a lo largo del gradiente de concentración de los nutrientes más importantes (bases intercambiables y fósforo).

MATERIALES Y MÉTODOS

Área de estudio: El Chocó biogeográfico se extiende desde el extremo sur de Panamá hasta el norte de Ecuador y cubre aproximadamente 1 500 km (^{Conservation International, 2014}). Presenta una alta precipitación de ~ 11 000 mm en el centro del Chocó biogeográfico, la cual va disminuyendo hasta 2 000 - 3 000 mm en los extremos norte y sur. La temperatura media anual es de 27 °C y la humedad relativa es 90 % (^{Eslava, 1992}; ^{Poveda, Rojas, Rudas, & Rangel, 2004}). El trabajo se realizó en Colombia en el departamento de Chocó, Municipio Acandí, corregimiento de Capurganá (2 983 mm) y en Ecuador en la provincia de Santo Domingo de los Tsáchilas (3 340 mm).

Muestreo de campo y análisis de suelos: En el 2011 y 2017, recolectamos datos sobre la distribución y abundancia de especies de palmas en dos localidades del Chocó: Capurganá y Santo Domingo (Tabla 2). Utilizamos imágenes satelitales para encontrar sitios apropiados para establecer 20 transectos en bosques de tierras bajas no inundables sin aparente influencia humana (10 en Capurganá y 10 en Santo Domingo). Cada transecto tuvo un tamaño de 5 × 500 m y fue dividido en 100 subunidades de 5 × 5 m. Estos se ubicaron estratificadamente a una distancia de 500 m cada uno siguiendo el protocolo establecido por ^{Balslev et al. (2010}). Se registró la localización geográfica en el punto medio de cada transecto con un GPS utilizando el sistema de coordenadas Universal Transversal Mercator (UTM). En cada subunidad, identificamos y contamos todos los individuos de palmas presentes de todos los estados de crecimiento (plántulas, juveniles, subadultos y adultos) siguiendo las definiciones de categorías de crecimiento propuestas por Balslev et al. (2010). Las muestras de las especies recolectadas fueron depositadas en los herbarios de la Universidad Nacional de Colombia sede Medellín (MEDEL) y la Pontificia Universidad Católica del Ecuador (QCA).

Las muestras de suelos fueron recolectadas a lo largo de cada transecto a 0, 250 y 500 m, correspondiendo a las subunidades 1, 50 y 100. En total, se analizaron 57 muestras de suelo, tres por cada transecto (Capurganá: 28, Santo Domingo: 29). Para tres de ellos, solo se analizaron dos muestras de suelo. Cada una fue obtenida al mezclar cinco submuestras recolectadas en cada esquina y en el centro de la subunidad a una profundidad de 0 - 10 cm. Primero, se retiró la capa de materia orgánica, las piedras, raíces y material diferente al suelo. Las muestras fueron secadas al aire libre y almacenadas en bolsas ziploc en campo. Seguidamente, se analizaron en el laboratorio de suelos de la Universidad de Aarhus de Dinamarca. La extracción de nutrientes (Al, Ca, Mg, K, Na y P) se realizó con la técnica ^{Mehlich-III (Mehlich, 1984}) y las concentraciones se analizaron con espectrometría de emisión con plasma de acoplamiento inductivo (ICP-OES). Para los análisis, se utilizó el promedio de la concentración de nutrientes de cada transecto.

Análisis de los datos: Para calcular la similitud florística entre transectos, se utilizó el índice de Sørensen, para datos de presencia - ausencia, y el índice de Steinhaus, para datos de abundancia. Para calcular este último, se utilizó la fórmula 2 W / (A + B), donde W es la abundancia de especies compartidas entre dos transectos, A es la suma de las abundancias de todas las especies en el primer transecto y B es la suma de las abundancias de todas las especies en el segundo transecto (^{Legendre & Legendre, 2012}). Para explorar las diferencias en la composición florística entre los sitios de muestreo, se realizó una ordenación usando la prueba de escalamiento multidimensional no métrico (NMDS), la cual mapea los valores de disimilitud florística entre sitios en un espacio de ordenación con un número predefinido de dimensiones (^{Minchin, 1987}). También se produjo una ordenación bidimensional utilizando la función metaMDS en el paquete R _VEGAN (^{Oksanen et al., 2016}) y se proyectaron las variables de suelo en la ordenación utilizando la función envfit. Asimismo, se calcularon los coeficientes de correlaciones al cuadrado (R²) para cada predictor en relaciones con los ejes NMDS y se determinó su importancia o probabilidad estadística con una prueba de permutación de Monte Carlo utilizando 999 permutaciones aleatorias.

Distancias geográficas: Para medir la distancia geográfica entre sitios y transectos, se utilizaron las coordenadas geográficas (UTM). Por su parte, la matriz de distancia geográfica se calculó utilizando la distancia euclidiana con datos transformados a logaritmo en base 10 (^{Andersen et al., 2010}). Adicionalmente, se aplicó una transformación logarítmica con el fin de hacer más simétrica la distribución de los datos y con el fin de que se acercara más a la distribución normal.

Test de Mantel y Mantel parcial: Se realizó esta prueba para calcular si la disimilitud en composición florística está correlacionada con la disimilitud en la concentración de nutrientes de suelo. Este se aplicó tanto con datos de presencia - ausencia como con datos de abundancias. Es importante aclarar que transformamos todas las variables edáficas al logaritmo en base 10 antes de los análisis. Cada una de ellas se estandarizó (promedio = 0, varianza = 1) antes de combinarlas en una matriz de distancia euclidiana. Asimismo, se realizó un test parcial de Mantel para controlar la autocorrelación espacial y así evitar inflar la importancia de las variables ambientales. Se utlizó el estadístico de Mantel (r_M) y establecimos la significancia estadística de las correlaciones en P < 0.001 con un test de Monte Carlo con 999 permutaciones aleatorias. Por último, para observar la relación entre la abundancia de las palmas del sotobosque, dosel y las bases intercambiables, se realizó un análisis de correlación y regresión lineal simple. Todos las pruebas con Mantel y Mantel parcial fueron realizadas con la función mantel en R, del paquete _VEGAN (^{Oksanen et al., 2016}) y, en general, todos los análisis fueron realizados con el programa estadístico R versión 3.3.2. (^{R Development Core Team, 2017}).

RESULTADOS

Patrones generales florísticos y edáficos: Un total de 9 158 individuos y 33 especies de palmas fueron registrados en 20 transectos que cubren una área de 5 ha (Tabla 3). La abundancia de palmas en los trasenctos varió en un rango de 96 - 1 670 individuos (media ± DS: 868 ± 319) y 7 - 15 especies (media ± DS: 6 ± 2). Las especies más frecuentes fueron Synechanthus warscewiczianus (presente en 18 de 20 transectos), Socratea exorrhiza (17 / 20) y Geonoma cuneata (16 / 20). Las más abundantes fueron Welfia regia (2 060 individuos), Phytelephas macrocarpa (1 726) y Synechanthus warscewiczianus (1 088). Cuatro especies (Attalea allenii, Desmoncus cirrhifer, Geonoma calyptrogynoidea y Geonoma undata) fueron encontradas en solo un transecto y siete (Astrocaryum standleyanum, Bactris barronis, Bactris gasipaes var. chichagui, Chamaedorea pinnatifrons, Geonoma interrupta, Iriartea deltoidea, Phytelephas macrocarpa) se identificaron en al menos 10 transectos (Tabla 3). De las 33 especies registradas, 15 fueron exclusivas de Santo Domingo y 11 de Capurganá. Únicamente siete se comparten entre los dos sitios. Las especies de dosel más abundantes fueron Welfia regia y Socratea exorrhiza y las del sotobosque fueron Phytelephas macrocarpa y Synechanthus warscewiczianus (Tabla 3). Los valores de disimilitud florística para el índice de Sørensen fueron 0.04 - 0.90 (media = 0.48), y para el índice de Steinhaus fueron 0.19 - 0.99 (media = 0.77).

Por su parte, las concentraciones de nutrientes no presentaron mayores diferencias entre regiones (t = 1.34, P = 0.21). Sin embargo, los suelos de Capurganá presentaron concentraciones más altas que los de Santo Domingo (Tabla 4). La concentración de bases intercambiables (Ca, K, Mg) fue cuatro veces mayor en Capurganá que en Santo Domingo y presentó una alta variación entre transectos en ambas regiones. Los nutrientes que más contribuyeron a la concentración de bases intercambiables fueron Ca y Mg (Tabla 4). Cabe agregar que la concentración de los nutrientes como Ca, Mg y Na fueron mayores en Capurganá, mientras que el Al, pH, K, P fueron mayores en Santo Domingo (Tabla 4).

Los nutrientes que más aportaron a la variación de la composición florística en las regiones fueron Ca, Mg, Na, K y P (Tabla 5) según el análisis de NMDS. Los valores de correlaciones del NMDS muestran que la composición florística y abundancia de palmas está fuertemente correlacionada con todas las variables edáficas, con excepción del Al y el pH (Tabla 5, Fig. 1). Además, el pH, P y K muestran un patrón similar y son más importantes en Santo Domingo (Fig. 1). El pH y el Al fueron las variables que menos se correlacionaron con el resto de las variables, mientras que el Na y Mg presentaron las mayores correlaciones (r = 0.90), seguidos por Ca y Mg (r = 0.88) (Tabla 6).

Fig. 1 Escalamiento multidimensional no métrico (NMDS) que muestra las proyecciones de las variables edáficas en 20 transectos en el Chocó biogeográfico, basado en la disimilitud en la composición florística medida con del índice de Steinhaus (Stress = 0.071). Fig. 1. Non-metric multidimensional scaling (NMDS) with projection of soil variables in 20 transects in the Chocó region, based on the dissimilarity of floristic composition as measured with the Steinhaus index (Stress = 0.071). 

Relación entre la composición florística y edáfica: La correlación entre las matrices de disimilitud en composición florística y edáfica, según los datos de presencia - ausencia, fue 0.69 y, según los datos de abundancia, fue 0.63 (P = 0.001 para ambos casos). Por otro lado, la correlación entre la disimilitud florística y geográfica fue estadísticamente significativa, tanto para los datos de presencia - ausencia (r_M = 0.75, P = 0.001) como de abundancia (r_M = 0.71, P = 0.001) (Tabla 7A). Todas estas relaciones entre la distancia geográfica y las variables edáficas fueron estadísticamente significativas excepto para el pH y el Al (Tabla 7B). El test de Mantel parcial entre las matrices de disimilitud florística y edáfica, que controlaba el efecto de la distancia geográfica, también arrojó resultados estadísticamente significativos para datos de presencia - ausencia (r_M = 0.31, P = 0.001) y de abundancia (r_M = 0.24, P = 0.002) (Tabla 7A).

Relación entre bases intercambiables y abundancia: Se registraron 24 especies de palmas del sotobosque (Fig. 2A, Fig. 2B, Fig. 2C, Fig. 2D, Fig. 2E y Fig. 2F) y nueve del dosel (Fig. 2G, Fig. 2H, Fig. 2I, Fig. 2J, Fig. 2K y Fig. 2L). Al considerar ambas localidades conjuntamente, se encontró una relación significativa entre las bases intercambiables (Ca + K+ Mg) y la abundancia de las palmas del sotobosque (r² = 0.37, P = 0.004). Sin embargo, al analizar por separado ambas zonas, se encontró un patrón diferente. En Capurgná, por un lado, las bases intercambiables no presentaron una relación estadísticamente significativa con la abundancia de palmas del sotobosque (Fig. 3A) o dosel (Fig. 3B). Por otro lado, en Santo Domingo las bases intercambiables sí presentaron una relación estadísticamente significativa con la abundancia de palmas de sotobosque (Fig. 3C), pero no con las de dosel (Fig. 3D).

Fig. 2 Ejemplo de especies de palmas de sotobosque (A - F) y dosel (G - L) que crecen en los bosques del Chocó biogeográfico. (A) Phytelephas macrocarpa, (B) Chamaedorea allenii, (C) Bactris barronis, (D) Geonoma cuneata, (E) Chamaedorea deneversiana, (F) Aiphanes tricuspidata, (G) Attalea colenda, (H) Astrocaryum standleyanum, (I) Bactris gasipaes var. chichagui, (J) Iriartea deltoidea, (K) Socratea exorrhiza, (L) Wettinia quinaria. Fig. 2. Examples of palm species of the understory (A - F) and canopy (G - L) that grow in the Chocó region’s forests. (A) Phytelephas macrocarpa, (B) Chamaedorea allenii, (C) Bactris barronis, (D) Geonoma cuneata, (E) Chamaedorea deneversiana, (F) Aiphanes tricuspidata, (G) Attalea colenda, (H) Astrocaryum standleyanum, (I) Bactris gasipaes var. chichagui, (J) Iriartea deltoidea, (K) Socratea exorrhiza, (L) Wettinia quinaria. 

Fig. 3 Relación entre el número de individuos de palmas de (A) sotobosque de Capúrganá, (B) dosel de Capurganá, (C) sotobosque de Santo Domingo, (D) dosel de Santo Domingo con el promedio de concentración de bases intercambiables (Ca + Mg + K) en el Chocó biogeográfico. Fig. 3. Relationship between the number of palm individuals in the (A) understory of Capurganá, (B) canopy of Capurganá, (C) understory of Santo Domingo, (D) canopy of Santo Domingo with the concentration of mean exchangeable bases (Ca + Mg + K) in the Chocó region. 

Respuesta de las palmas más abundantes a las bases intercambiables y P: Las especies de sotobosque Geonoma cuneata, G. interrupta y Phytelephas macrocarpa fueron las más abundantes a bajas concentraciones de P, mientras que Synechanthus warscewiczianus fue la más abundante en suelos con altas concentraciones de P (Fig. 4A). A su vez, S. warscewiczianus mostró preferencia por suelos con bajas concentraciones de bases intercambiables, mientras que G. cuneata, G. interrupta y P. macrocarpa abundaron en los transectos con altas concentraciones de bases intercambiables (Fig 4B). Entre las palmas del dosel, se encontró que Socratea exorrhiza crecía más en todo el gradiente de P y bases intercambiables (Fig. 4C, Fig. 4D, mientras que Iriartea deltoidea lo hacía en suelos con concentración intermedias de P y bases intercambiables (Fig. 4C, Fig. 4D).

Fig. 4 Abundancia de las palmas del sotobosque y del dosel a lo largo del gradiente de concentración de fósforo (P) y bases intercambiables (Ca + K + Mg) en 20 transectos en el Chocó biogeografico: (A - B) sotobosque. (C - D) dosel. Fig. 4. Abundance of understory and canopy palms along the concentration gradient of phosphorus (P) and exchangeable bases (Ca + K + Mg) in 20 transects in the Chocó region: (A - B) understory. (C - D) canopy. 

DISCUSIÓN

La concentración de nutrientes de los suelos de Capurganá y Santo Domingo presentó una variación que está dentro de los rangos resgistrados anteriormente en el Chocó biogeográfico (^{Faber-Langendoen & Gentry, 1991}; ^{Galeano, 2001}; ^{Quinto-Mosquera & Moreno-Hurtado, 2016}) y en la Amazonia que también posee los suelos de tierras firmes y un promedio similar de precipitación (^{Poulsen, Tuomisto, & Balslev, 2006}; ^{Kristiansen et al., 2012}). Una explicación a la alta variación en la concentración de nutrientes de los suelos en esta investigación es que los suelos en el Chocó suelen ser más ricos en nutrientes en lugares donde la precipitación es menor (Quinto-Mosquera & Moreno-Hurtado, 2016) y, en este estudio, se muestrearon transectos con un rango de precipitación de 2 983 - 3 340 mm, donde presentaban mayores promedios de concentración las bases intercambibles en Capurganá que en Santo Domingo. Mientras tanto, en el caso del P, por ejemplo, la concentración promedio en Santo Domingo fue más del doble que en Capurganá (Tabla 4).

Nuestros resultados apoyan la hipótesis general de que la composición florística de los bosques tropicales está relacionada con los nutrientes del suelo. Si bien esto se ha demostrado en bosques amazónicos y en Panamá (^{Poulsen et al., 2006}; ^{Kristiansen et al., 2012}; ^{Garibaldi, Nieto-Ariza, Macía, & Cayuela, 2014}; ^{Tuomisto et al., 2014}; ^{Cámara-Leret et al., 2017}; ^{Stevenson et al., 2018}; ^{Muscarella et al., 2018}), nuestro estudio es el primero en resaltar que, en las tierras bajas del Chocó, la composición de las comunidades de palmas está relacionada con los nutrientes: los suelos con concentraciones bajas en bases intercambiables son más ricos en especies.

Asimismo, en esta investigación, las variables edáficas mayormente correlacionadas con la composición florística, según el test de Mantel, fueron Mg, Ca, Na, y P. Esto sugiere que, al igual que en Amazonia (^{Peñas-Claros et al., 2012}, ^{Cámara-Leret et al., 2017}), en las tierras bajas del Chocó biogeográfico las bases intercambiables y el P son los factores importantes que determinan la composición de las comunidades de palmas. El rol del P como determinante de la distribución de árboles también se ha comprobado en las Dipterocarpaceae en Borneo (^{Paoli, Curran, & Zak, 2006}) y Entandrophragma en África (^{Hall, Mckenna, Ashton, & Gregoire, 2004}). Las bases intercambiables son más importantes para explicar la composición de especies de palmas en Capurgana, mientras que, en Santo Domingo, el P explica la riqueza de especies (Fig. 1). Nuestros resultados son consistentes con los de ^{Garibaldi et al. (2014}) quienes, en la Península Azuero de Panamá, encontraron que la fertilidad y el P presentan una alta correlación con la composición florística.

Por otra parte, se encontró una pequeña relación positiva entre las bases intercambiables y la abundancia de palmas de sotobosque en las tieras bajas del Chocó biogeográfico. Estos resultados son el reflejo de las altas abundancias de algunas especies de palmas del sotobosque como Phytelephas macrocarpa que es una especie indicadora de suelos ricos en nutrientes (^{Kahn & Mejia, 1990}). Aunque no se detectó una diferencia significativa entre las bases intercambiables y las palmas del dosel, esto puede deberse a que esas plantas tienden a ser más generalistas. Por ejemplo, Socratea exorrhiza (una palma de dosel) estuvo presente en casi todos los transectos y con abundancias muy semejantes entre cada uno. Lo mismo se ha reportado en la Amazonia (^{Cámara-Leret et al., 2017}), donde Socratea exorrhiza crece tanto en suelos pobres como ricos en nutrientes. Lastimosamente, nuestros datos no permiten explorar la respuesta de especies de dosel con baja abundancia, como Oenocarpus bataua, la cual crece en suelos pobres en nutrientes en la Amazonia occidental (Cámara-Leret et al., 2017).

Las bajas correlaciones entre la fertilidad y palmas, tanto de sotobosque como de dosel, pueden deberse al que las muestras se limitaron a una parte del gradiente climático del Chocó (Capurganá al norte y Santo Domigo al sur). También se observó que algunas especies fueron más abundantes en sitios con determinada concentración de nutrientes, por ejemplo, la palma de dosel Iriartea deltoidea fue más abundante en suelos medianamente fértiles y con concentración medias de P. Esto es similar a lo reportado en la Amazonia occidental (^{Cámara-Leret el al., 2017}). Otras palmas como Geonoma cuneata son más abundantes en suelos pobres en bases intercambiables y con baja concentración de P. Igualmente, Wettia quinaria, Welfia regia y Oenocarpus bataua eran más frecuentes en suelos pobres en nutrientes y en lugares con altas precipitaciones. Por el contrario, especies como Phytelephas macrocarpa y Astrocayrum standleyanum fueron muy abundantes en sitios ricos en nutrientes.