SciELO - Scientific Electronic Library Online

 
vol.101 issue3Structural attributes and habitat of Juniperus jaliscana in Talpa de Allende, Jalisco, Mexico author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Botanical Sciences

On-line version ISSN 2007-4476Print version ISSN 2007-4298

Bot. sci vol.101 n.3 México Jul./Sep. 2023  Epub July 31, 2023

https://doi.org/10.17129/botsci.3204 

Ecology

Landscape heterogeneity drives spatial distribution of palm community in a Neotropical rainforest reserve affected by defaunation

La heterogeneidad del paisaje regula la distribución espacial de la comunidad de palmas en un bosque húmedo neotropical afectado por defaunación

Katia V. Miranda-Gallegos1  , Writing - original draft, Investigation, Formal analysis
http://orcid.org/0000-0003-3440-7952

Armando Navarrete-Segueda1  2  *  , Investigation, Formal analysis, Writing - original draft, Writing – review & editing, Project administration
http://orcid.org/0000-0001-6614-8751

Jorge Cortés-Flores3  , Writing – review & editing, Formal analysis
http://orcid.org/0000-0001-9578-7827

M. Lourdes González-Arqueros4  , Writing – review & editing, Formal analysis
http://orcid.org/0000-0001-8775-2886

Estefanía E. Acosta-Pérez5  , Investigation
http://orcid.org/0000-0002-6975-7007

Guillermo Ibarra-Manríquez1  , Writing – review & editing, Project administration
http://orcid.org/0000-0002-3739-8660

1Instituto de Investigaciones en Ecosistemas y Sustentabilidad, Universidad Nacional Autónoma de México, Morelia, Michoacán, Mexico.

2Escuela Nacional de Estudios Superiores, Unidad Morelia, Universidad Nacional Autónoma de México, Morelia, Michoacán, Mexico.

3Jardín Botánico, Instituto de Biología, Sede Tlaxcala, Universidad Nacional Autónoma de México, Tlaxcala, México.

4 CONAHCYT -Instituto de Investigaciones en Ciencias de la Tierra, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, Mexico.

5Herbario IBUG, Instituto de Botánica, Departamento de Botánica y Zoología, Universidad de Guadalajara, Zapopan, Jalisco, Mexico.


Abstract

Background:

Our study presents the first assessment of the relationships between environmental heterogeneity and the palm community at the landscape scale at the northern boundary of tropical rainforest distribution in America.

Question:

How does landscape-level heterogeneity (edaphic, topographic as well as the density, diversity, and tree biomass) shapes the palm community in a protected reserve affected by defaunation?

Study site and date:

The study was carried out in Los Tuxtlas Tropical Biology Station, Veracruz, Mexico (2020).

Methods:

We delineated landscape units (LUs) based on soil and geomorphology to characterize the biophysical mosaic. We built rank-abundance curves to compare the palm richness and their abundance among LUs. Permutational multivariate analysis of variance (PERMANOVA) and nonmetric multidimensional scaling (NMDS) were applied to detect shifts in the palm community across LUs.

Results:

Integrating abiotic and biotic components by delimiting LUs allowed us to evaluate the multifactorial effect of environmental heterogeneity on the palm community. The density and composition of the palms changed significantly among the LUs, influenced by soil coarse fragments, elevation, slope inclination, and annual precipitation. Tree richness and density had a significant effect on the palm community only when analyzed in combination with the soil-topography variables.

Conclusions:

The results reveal that there is high edaphic and topographic heterogeneity that influences palm species distribution at the landscape scale. The environmental heterogeneity approach allows capturing novel abiotic variation to explain the distribution patterns of palm species and their coexistence in highly diverse tropical communities.

Keywords: Landscape units; Los Tuxtlas Biosphere Reserve; Soil rooting depth; Species diversity

Resumen

Antecedentes:

Se evalúa la influencia de la heterogeneidad ambiental en la comunidad de palmas a escala del paisaje en el límite norte de la distribución de la selva tropical húmeda en América.

Preguntas:

¿Cómo la heterogeneidad a nivel de paisaje (edáfica, topográfica, así como la densidad, diversidad y biomasa arbórea) influye en la distribución de las especies de palma en una reserva forestal afectada por la defaunación?

Sitio y años de estudio:

El estudio se realizó en la Estación de Biológica Tropical Los Tuxtlas, Veracruz, México, en 2020.

Métodos:

Se delinearon unidades de paisaje (UP) basadas en suelo y geomorfología, se obtuvieron datos de la riqueza y abundancia de los árboles para caracterizar la heterogeneidad ambiental. Se analizó como influye la heterogeneidad ambiental en la composición, riqueza y abundancia de las especies de palma.

Resultados:

La abundancia y composición de las especies de palmas difirió entre las UP, en ambos casos fueron influenciadas por la presencia de fragmentos gruesos del suelo, altitud, inclinación de la pendiente y precipitación anual. La riqueza y densidad de árboles afectaron significativamente a la comunidad de palmas, pero sólo cuando se analizaron en combinación con variables edáficas y topográficas.

Conclusiones:

Los resultados revelan que existe una marcada heterogeneidad edáfica y topográfica que influye en la distribución de las especies de palmas a escala de paisaje, por lo que el enfoque de la heterogeneidad ambiental permite asociar la variación abiótica con los patrones de distribución de las especies y su coexistencia en comunidades tropicales altamente diversas.

Palabras clave: Diversidad de especies; Profundidad efectiva del suelo; Unidades de paisaje; Reserva de la Biosfera de Los Tuxtlas

Palms are one of the most diverse and abundant plant groups in tropical rainforest (TRF) and their diversity depends on multiple environmental factors which operate at different temporal and spatial scales (Eiserhardt et al. 2011). Edaphic conditions (particularly soil water availability, effective soil depth and texture) as well as topographic factors (e.g., relief position, elevation, and slope), have been shown a significant influence on the structure and composition of palm communities (Clark 1995, Svenning 2001a, b, Emilio et al. 2014, Muscarella et al. 2019). In turn, the relationship with biotic components is not clear since has been found contrasting results about effect of variables as forest canopy openness and tree density on palm species richness (Cintra et al. 2005, Rodrigues et al. 2014). In the past decades, the studies of Svenning (2001b) and Rodrigues et al. (2014) proved that the spatial covariation of biophysical components (soil, topography, and tree structure) influences the palm community. This knowledge represents the basis for understanding the influences of forest and landscape structure components on the palm species distribution and their coexistence. However, the effect of these biophysical drivers is still poorly understood in most Neotropical forests.

At the landscape scale, the spatial patterns of the palm community have been explored mostly by isolating abiotic and biotic components (Eiserhardt et al. 2011). However, soil, topography, microclimate, and tree structure are inextricably related, and their variation shapes the habitat mosaics that influence the distribution and density of plant species (Tylianakis et al. 2008, Baldeck et al. 2013, Li et al. 2018, Rodrigues et al. 2019, Schmitt et al. 2021). Consequently, the understanding of spatial patterns of palm species as well as identifying their predictors, requires the multivariate and spatially explicit assessment.

Neotropics harbor great plant diversity on a relatively small surface area (Raven et al. 2020). However, agricultural expansion in Latin America's tropical forests reduced forest cover at a rate of 40,000 km2 per year in the decade 2000-2010 alone (FAO & UNEP 2020). This loss of forest cover is one of the main drivers of biodiversity decline in tropical ecosystems (Barlow et al. 2018, Winkler et al. 2021). In response, throughout the Neotropics, geographical areas oriented to biodiversity conservation, habitat protection, and maintenance of ecosystem processes and services have been delineated (Cazalis et al. 2020). These protected areas, commonly considered intact, have been suffering an intense loss of fauna (Dirzo et al. 2014, Benítez-López et al. 2019).

Los Tuxtlas Tropical Biology Station (LTBS) reserve, which is located at the Sierra de Los Tuxtlas (Veracruz state, Mexico), represent the northern boundary of the TRF distribution in America and supports an old-growth TRF. Nevertheless, among the 40 eco-regions recognized in Mexico, the Sierra de Los Tuxtlas had the greatest forest cover decline during the period 2001-2014, and the remaining forest cover represents less than 10 % of its original cover (Bonilla-Moheno & Aide 2020, Von Thaden et al. 2020). Although there are still remnants, such as the LTBS, where high plant diversity continues to exist (Ibarra-Manríquez & Sinaca-Colín 1995), defaunation caused by human beings has eliminated or diminished populations of several medium and large mammal species (Dirzo & Miranda 1991, Martínez-Ramos et al. 2016). Such forest fauna decline has the potential to produce a variety of effects on plant populations, gradually modifying the plant communities’ composition and structure (Gardner et al. 2019, Villar et al. 2020, Valverde et al. 2021).

Regarding the effect of defaunation on plant ecology of LTBS, there is an unresolved debate about its influence on the trajectory of regeneration and plant diversity. Martínez-Ramos et al. (2016) proposed that loss of fauna favored demographic explosion of a hyper-dominant palm (Astrocaryum mexicanum Liebm. ex Mart.), which, in turn, reduced tree diversity and modified the forest’s composition. However, this conclusion has been challenged by Arroyo-Rodríguez & Melo (2016), firstly because it was based on a few small plots that are close together, and secondly, because it did not consider the role of environmental factors in shaping mesoscale plant spatial patterns. One approach to integrate these last factors is through the stratification of the landscape into discrete units based on biophysical components.

Landscape stratification based on soil and geomorphology has proven useful in capturing landscape heterogeneity and explaining variation in plant communities throughout Neotropical landscapes (e.g.,Clark et al. 1995, Denslow et al. 2019, Navarrete-Segueda et al. 2017). The effectiveness of this approach is based on the fundamental role of topography as a determinant of habitat differentiation and ecological processes at the landscape scale (Zonneveld 1989, Zinck et al. 2016). As a result, discrete landscape units (LUs) based on landforms can be detected, in which environmental conditions are relatively homogeneous within each type of landscape unit, while adjacent LUs differ in habitat characteristics (Dehn et al. 2001, MacMillan & Shary 2009). Consequently, we explored factors (edaphic, topographic as well as the tree density, diversity, and biomass) driving variation in the composition and distribution of the palm species in the LTBS, using the LUs delineated by Navarrete-Segueda et al. (2021).

Studies on the growth pattern of palms in LTBS have shown that continuous production of leaves from a single terminal stem allows age to be determined based on palm height. In A. mexicanum and Chamaedorea spp., stems >1 m indicate that individuals are mature of > 30 years (Piñero et al. 1984, Martínez-Ramos et al. 1988, Oyama et al. 1992, Martínez-Ramos et al. 2016). Therefore, by focusing on palms of this height, we can record the diversity and the density of individuals that were present in the decades before the defaunation occurred in the reserve.

The objectives of the study were to: i) document the environmental heterogeneity of major LUs of the LTBS, ii) identify major gradients of topography and soil that influence the composition and density of palm species, iii) evaluate the variation of both palm community attributes among types of LUs, and iv) assess the effects of abiotic (topography and soil) and biotic (density, biomass, and tree diversity) factors on the spatial variation of palm species composition and density from a multivariate perspective. We expected that LUs captured the variation in relief, soil parameters, and forest structure of the study area. Thus, we also expected significant differences among LUs in the palm community composition and density. From an abiotic perspective, the slope, due to its effect on habitat differentiation and its synergistic relationship with fertility and soil water storage (Zinck et al. 2016), was expected to regulate the spatial variation of density and diversity of palms. With respect to biotic variables (individuals’ density as well as the biomass and diversity of trees), we expected a negative effect on palm species (presence and density) across LUs (Rodrigues et al. 2014).

Material and methods

Study area. We carried out the study at LTBS (18° 35' N, 95° 05' W), which is located in a volcanic complex of the Upper Tertiary and Middle Quaternary (< 7 Ma), that is composed of basaltic andesite and basalt (Verma et al. 1993) and bordered by lava flows from the Holocene (Nelson & González-Caver 1992). The landscape in LTBS is topographically heterogeneous, with slopes ranging from 2 to 45°, covering an area of 640 ha, where the elevation ranges between 150 and 700 m asl (Ibarra-Manríquez et al. 1997). The annual average temperature is 24.5 °C and the mean annual precipitation (MAP) is 3,928 mm, with a dry season from March to May, where it rains less than 10 % of the total annual precipitation (SMN 2010). The most abundant tree species are Damburneya ambigens (S.F. Blake) Trofimov, Guarea glabra Vahl, Omphalea oleifera Hemsl., and Pseudolmedia glabrata (Liebm.) C.C. Berg (Ibarra-Manríquez et al. 1997). Additionally, 13 species of palms have been recorded; the most abundant include A. mexicanum, Chamaedorea alternans H. Wendl., C. pinnatifrons (Jacq.) Oerst., and C. tepejilote Liebm. Although there are few palm species, their high density makes them an important element of this TRF (Ibarra-Manríquez et al. 1997, Martínez-Ramos et al. 2016). Despite being surrounded by cattle pasture, the protected status of the LTBS has protected the forest from fires, logging, and plant extractions for more than 50 years, such that it is representative of conserved TRF in the region (Figure 1).

Figure 1 (A) Study area and location of plots on the land units (High Altitude Slopes Strongly Inclined, HSSI; Low Altitude Slopes Strongly Inclined, LSSI; Lava Flows, LF; Piedmont, PD; Steeply Dissected Slopes, SDS). (B) Satellite image of the study area on Google Earth. Data: Google ©2021, Maxar Technologies/TerraMetrics. 

Environmental heterogeneity characterization. We stratified the study area based on the five LUs delimited by Navarrete-Segueda et al. (2021). These discrete units are based on the systematic integration of environmental components such as geology, relief, and soil (Zinck et al. 2016). This approach has been used to stratify the environmental mosaic and the potential habitats in the study area (Navarrete-Segueda et al. 2021) as well as in other Neotropical rainforests (Clark et al. 1995, Brown et al. 2013, Denslow et al. 2019). Four of the five LUs were associated with cinder cones (late Tertiary and early Quaternary), that differ in inclination, elevation, microclimate, and their soils, which were formed from volcanic ash. The fifth landscape unit corresponds to a landscape associated with Holocene volcanic events and is characterized by shallow stony soils (Table 1). The selected LUs are representative nearly of 75 % of the landforms of the LTBS (640 ha). In each of the LUs, three plots of 20 × 50 m (0.1 ha) were established, covering a total sample area of 1.5 ha. Plots on steep slopes were oriented transverse to the slope to minimize the effect of topographic variation. Estimated inter‐site distances ranged between 400 and 4,550 m.

Table 1 Changes in environmental and structural variables of tree community among LUs [mean and standard error (S.E.)]. Available water holding capacity (AWHC). High altitude slopes strongly inclined (HSSI), Low altitude slopes strongly inclined (LSSI), Lava flows (LF), Piedmont (PD), Steeply dissected slopes (SDS). Means with the same superscript are not statistically different among geopedologic land units. 

  PD LSSI HSSI SDS LF
Abiotic factors
Slope (°) 10.3(0.9)ab 19.7(2.9)b 13.0(1.7)ab 29.7(3.9)c 7.7(1.5)a
Elevation (m asl) 199.3(34.6)c 273.7(29.45)ac 481.7(22.7)b 419.3(21.1)b 348.3(31.2)a
Aspect 183.6(46)a 135.4(58.8)a 82.3(13.5)a 219.5(76.1)a 90.2(13.9)a
Annual precipitation (mm) 4878.0(10.6)b 4860.3(29.9)b 4715.7(15.3)a 4845.7(38.0)b 4762.3(8.4)a
Clay (%) 23.9(3.9)ab 27.5(1.7)b 22.0 (1.2)ab 21.7(1.0)ab 11.7(6.0)a
AWHC (L m-2) 88.7(14.7)b 94.14(20.7)b 77.7(13.5)b 123.9(30.2)b 4.0(0.6)a
Coarse fragments (%) 4.9(2.2)b 3.1(2.1)b 0.3(0.3)b 5.5(3.9)b 87.5(2.5)a
Rooting depth (cm) 51.7(1.7)c 47.3(4.3)c 74.3(14.8)b 57.0(12.3)bc 18.7(2.4)a
Biotic factors
Tree richness (0D) 24.7(1.8)a 23.7(0.9)a 27.3(3.8)a 23.3(2.0)a 20.0(3.0)a
Tree biomass (t 0.1ha-1) 33.7(2.7)a 45.9(17.6)a 45.8(9.9)a 28.4(29.1)a 27.5(6.3)a
Tree common species (1D) 19.8(1.5)a 20.2(1.3)a 22.0(3.8)a 18.3(0.3)a 14.5(2.5)a
Tree density (ind. 0.1ha-1) 41.0(1.5)a 41.7(2.8)a 44.7(5.7)a 55.3(12.6)a 38.3(4.7)a

In each plot, we described soil profiles, based on Schoeneberger et al. (2012), within the rooting depth (Siebe et al. 1996). The available water holding capacity (AWHC) of each horizon was estimated and interpreted according to Eckelmann et al. (2005) and Siebe et al. (1996). This variable integrates the soil texture, organic matter content, bulk density, percentage of coarse fragments, and the thickness of each horizon estimated in the field. In addition, texture, rooting depth, and percentage of coarse fragments were selected as indicators of nutrient storage and soil fertility (Bünemann et al. 2018). We also integrated elevation and slope as relief variables since they have been shown to influence the spatial variation in vegetation structure and diversity (Baldeck et al. 2013). To obtain the abiotic components related to relief, we parameterized the topographic data of slope and elevation, based on a digital elevation model (DEM) constructed using contour lines to 10 m resolution of INEGI (2010), processed in ArcMap 10.1 (ESRI®). We extracted the average value and standard deviation of these variables for each plot.

We generated climatic data because in the Sierra de Los Tuxtlas there is a low density of meteorological stations (Gutiérrez-García & Ricker 2011). To accurately interpolate and strengthen the climate surface, we included weather data from the surrounding area. We employed eight meteorological stations, with elevations ranging from 4 to 1,200 m. Yearly average precipitation (mm yr-1) and temperature (°C) were generated with Inverse Distance Weighting (IDW) interpolation through ArcMap (ESRI®). IDW interpolation is a mathematical (deterministic) function that assumes that closer values are more related than farther values (Chen & Liu 2012).

In each 0.1 ha plot (20 × 50 m), we recorded tree individuals with a diameter at breast height (dbh) ≥ 10 cm and identified them to the species level. Using dbh data, we calculated the biomass per tree species (t dry mass 0.1 ha-1) based on the allometric equation [1] applied to tropical regions with a MAP between 2,000- and 4,000-mm year-1 (Brown 1997, Rügnitz et al. 2008):

AGBt Kg dry mass=exp-2289+2649×Indbh-0.021×Indbh2 

To estimate AGB per ha we used the equation [2]:

AGB Mg dry mass 0.1 ha-1= at/1000

Where ∑at is the sum of the dry mass of all trees in the plot (AGBt) of 0.1 ha and 1000 is the factor to convert kg into Mg.

The use of specific allometric equations based on forest type has proven to be an efficient method that captures 95 % of biomass variation (Brown 2002, Berenguer et al. 2015). In each plot of 20 × 50 m (0.1 ha). We recorded and identified all the palms with a height greater than 1 m, except the adult individuals of Reinhardtia gracilis (H. Wendl.) Drude ex Dammer, which are smaller than 1 m. However, it has been shown that this species responds to soil conditions (Mendoza & Franco 1998) and is an important element of the palm community in the study area (Ibarra-Manriquez & Sinaca-Colín 1995). Therefore, individuals of R. gracilis were sampled only when flowers or fruits were observed.

Data analysis. We carried out generalized linear models (GLM) and multiple comparisons of means with post hoc Tukey tests to assess differences in abiotic and biotic components among LUs using ‘multcomp’ package in R (R Core Team 2022). We apply error distribution of Poisson to count data and Gaussian to continuous data (after normality test), and log and identity link functions, respectively (Crawley 2007). We applied a principal component analysis (PCA) to describe major gradients in the soil and topography to reduce the number and collinearity of variables (Laurance et al. 1999). Therefore, we expected that PCA would allow us to detect and classify meaningful variables related to the palm community parameters. The analysis was performed with the "scale" argument to standardize the data, results were centered by mean and, variables were presented by their correlations. The PCA was performed using the dudi.pca function of package ˈade4ˈ (Dray & Dufour 2007).

Data of individuals and richness were used to obtain mean (± S.E.) values of density (number of palms in 0.1 ha) and effective number of species (0D and 1D) for each landscape unit, following Jost (2006). We carried out GLM and multiple comparisons of means to assess differences in 0D, 1D, and density of palms (0.1 ha) among LUs.

We constructed rank-abundance curves to assess changes in the dominant species of palms among the LUs, lumping the data of the three plots per unit. The curves were constructed following Magurran (2004). Linear models were fit for rank- abundance curves and differences between units were tested using an analysis of covariance (ANCOVA) considering as an explanatory variable the rank (rank; as numeric) and landscape unit (Landscape unit; as categorical). To compare the slope of rank-abundance curves according to LUs, abundance of species was transformed to Log10. Differences among slopes were evaluated between LUs (Izsák 2006).

The effects of environmental heterogeneity on palm community were analyzed by a permutational multivariate analysis of variance (PERMANOVA) (Anderson & Braak 2003). This analysis was performed with the adonis2 function of ˈveganˈ package of R. The environmental data were normalized with the scale function. PERMANOVA is especially useful for partitioning multivariate data in response to complex designs (e.g., LUs) (Anderson 2017). Values of P were obtained using 9,999 permutations. Non-metric multidimensional scaling (NMDS) techniques were used to visualize the PERMANOVA results, which allowed us to examine the multivariate relationship of palm species (presence and abundance) with biophysical factors among LUs. The NMDS analysis were executed using the function metaMDS (ˈveganˈ package) in R (Oksanen et al. 2019). We implemented NMDS with Bray-Curtis dissimilarities to density and with the Jaccard method to the presence-absence data. Since NMDS uses rank order information, this analysis represents a highly flexible method to explore the proximities and resemblances of the species composition and the structure of palm community in relation whit environmental heterogeneity (McCune & Grace 2002).

Results

Environmental heterogeneity among the LUs. The GLM results showed that LUs are different in several of the analyzed environmental factors (Table 1). The slope presented a high variation. The units with low inclination were Piedmont (PD) and Lava Flows (LF) (< 11° each). These units differed significantly from Low Altitude Slopes Strongly Inclined (LSSI) and Steeply Dissected Slopes (SDS) (> 19° slope) (F4, 10 = 56.86, P < 0.000). The PD unit was found at lower elevations than LF (Table 1; Figure 1). High Altitude Slopes Strongly Inclined (HSSI) and SDS have a higher elevation, corresponding to the highest parts of the volcanoes. The LF unit, whose soils were formed by basaltic rock flows, had high surface stoniness (> 85 %), leading to soils with shallow rooting depth and low water storage (Table 1). In the other units, because the soils were formed from volcanic ash, the rooting depth and the soil water storage were greater, while the presence of coarse fragments was low. There were no significant differences in relief aspect and the density, biomass, and tree diversity among the LUs (F4,10= 36844, P = 0.301; F4,10 = 8.115, P = 0.0875; F4,10 = 3330134039, P = 0.78; F4,10 = 99.68, P = 0.133, respectively) (Table 1).

The first axis of the PCA (42 % of the variance) showed a gradient associated with the AWHC for plants and the coarse fragments in the soil (Table 2; Figure S1, Supplementary material). Coarse fragments in the soil are related to the LF unit, while, at the opposite extreme are the LUs recognized in soils formed from volcanic ash, with less coarse fragments content and larger stores of water in the soil, mainly SDS (Table 1; Table 2). On the second axis (26 % variance), the LUs are differentiated by a gradient determined by the elevation and the annual precipitation (Table 2). At one end of this axis are the PD unit and one plot corresponding to LSSI, located at low elevation and with high precipitation, separated from HSSI that had the opposite trend in these two environmental variables (Figure S1, Supplementary material).

Table 2 Loadings of the soil properties and topography that determine the first two axes of the Principal Component Analysis (PCA). Significance levels are based on a Pearson's correlation between soil properties and topography and PCA axes: * P < 0.05, ** P < 0.01, *** P < 0.001. 

Variables Axis 1 Axis 2
Annual precipitation (mm) -0.37ns -0.87**
Aspect (°) -0.34ns -0.14ns
Available water holding capacity (L m2) -0.94*** 0.14ns
Clay (%) -0.60ns -0.35ns
Coarse fragments (%) 0.92*** 0.02ns
Elevation (m a.s.l) -0.06ns 0.94***
Rooting depth (cm) -0.76* 0.54ns
Slope (°) -0.63* 0.05ns
Cumulative percentage of explained variance (%) 42 68

Palm composition and density in the context of environmental heterogeneity. A total of 10 palm species and 2,761 individuals were recorded. No significant difference was found in species richness (0D) (F4, 10= 0.81, P = 0.92) and the number of common species (1D) (F4, 10 = 7.63, P = 0.153) of palms among the LUs (Table S1, Supplementary material). We found that palm density differed among LUs, except between LSSI and SDS (F4, 10 = 61.98, P < 0.000). Palm density was highest in the LF unit and lowest in HSSI (Table S1, Supplementary material), which was clearly due to a low density of A. mexicanum and C. alternans in the latter unit (Figure 2; Table S2, Supplementary material).

Figure 2 Species-rank curves based on species abundance for palm community present in the land units of Los Tuxtlas Tropical Biology Station reserve, SE Mexico. In the X-axis, species are ranked in decreasing order of their abundance, while in the Y-axis, the number of palms is indicated. HSSI (High Altitude Slopes Strongly Inclined), LSSI (Low Altitude Slopes Strongly Inclined), LF (Lava Flows), PD (Piedmont), SDS (Steeply Dissected Slopes). 

Rank-abundance curves showed that the most abundant species in all the LUs was A. mexicanum (Figure 2). The second most abundant species differed among LUs; C. alternans in LF and PD, C. ernesti-augusti H. Wendl SDS and LSSI, R. gracilis in HSSI. We found significant differences among rank-abundance curves of LUs with contrasting soil properties (Table S3, Supplementary material). The comparison showed that HSSI curve differs from LF, and PD curves (F1= 0.455, P < 0.000 and F1= 0.3143, P = 0.013, respectively). Simultaneously, we found difference between SDS and LSSI curves (F1= 0.207, P = 0.036). We did not find significative differences among other LUs curves (Table S3, Supplementary material).

There were distinct distribution patterns in palm assemblages among the LUs. Results from the PERMANOVA showed significant effects of slope of relief (P = 0.009), elevation (P = 0.002), MAP (P = 0.04), soil coarse fragments (P = 0.002), and tree D1 (P = 0.008; Table 3) on palm species presence. These results clustered the plots with low inclination (PD) and separated them from the LUs with abrupt relief and higher elevation, as visualized by NMDS (stress value 0.15) in Figure 3A. Species that responded to this gradient are Bactris mexicana Mart., which was registered mainly in low-inclination plots, while C. ernesti-augusti and C. elatior Mart. were found mostly in plots with abrupt relief (Table S2, Supplementary material). We found that 33 % of the variation in palm species presence was explained by topography since relief elevation and slope explain 21 and 12 %, respectively (Table 3). Coarse fragments in soil and tree diversity (1D) accounted for 25 and 13 % of the variation, respectively.

Table 3 Results of Permutational Multivariate Analysis of Variance (PERMANOVA) for testing the effect of environmental variables on presence and density (number of individuals) recorded in 15 vegetation sampling plots (0.1 ha each one) located in Los Tuxtlas Tropical Biology Station reserve. 

Source of variation DF Sums of squares F R2 P
PRESENCE
Slope 1 0.06 5.23 0.13 0.009**
Elevation 1 0.10 8.69 0.21 0.002**
Mean annual precipitation 1 0.04 3.52 0.08 0.040*
Coarse fragments 1 0.12 10.46 0.25 0.002**
Tree density 1 0.00 0.35 0.01 0.761
Tree 1D 1 0.06 5.62 0.13 0.008**
Residual 8 0.09 0.19
Total 14 0.49 1.00
DENSITY
Coarse fragments 1 0.32 4.76 0.19 0.003**
Elevation 1 0.21 3.16 0.12 0.030*
Rooting depth 1 0.06 0.83 0.03 0.517
Slope 1 0.21 3.19 0.12 0.034*
Tree 0D 1 0.23 3.45 0.13 0.028*
Tree density 1 0.16 2.34 0.09 0.084
Residual 8 0.53 0.31
Total 14 1.71 1

Figure 3 Non-metric multidimensional scaling (NMDS) of (A) presence and (B) abundance of the palm community in the study landscape. Species presence and abundance in 20 × 50 m plots were used. Both species scores and plot scores for the first two ordination axes are plotted. Topography, soil, and tree community attributes in the plots were fitted to ordination to test the relationship between species composition and abundance, and environmental heterogeneity. MAP (Mean annual precipitation), HSSI (High Altitude Slopes Strongly Inclined), LSSI (Low Altitude Slopes Strongly Inclined), LF (Lava Flows), PD (Piedmont), SDS (Steeply Dissected Slopes). Initial letter of the genera: A (Astrocaryum), B (Bactris), C (Chamaedorea), and R (Reinhardtia). 

Results from the PERMANOVA showed significant effects of relief slope (P = 0.034), elevation (P = 0.03), soil coarse fragments (P = 0.003), and tree D0 (P = 0.028; Table 3) on palm species abundance. The NMDS in Figure 3B (stress value < 0.2) shows that, except for LSSI, the plots were grouped consistently with the LUs. Coarse fragments in the soil explained 19 % of the palm density variation (Table 3). C. tepejilote had the highest density in plots with high coarse fragments content (LF), while R. gracilis was most abundant in plots with deeper soils (and HSSI) (Table S2, Supplementary material). Similar to the palm species presence patterns, the density of B. mexicana, C. ernesti-augusti, and C. elatior responded to elevation and slope, which accounted for 24 % of variation in palm species density among the plots. Tree richness (D0) accounted for 13 % of density variation, with mainly, R. gracilis and C. alternans at the opposite extremes concerning tree richness and soil rooting depth interaction (Figure 3B).

Discussion

Environmental heterogeneity and the LUs. The delimited LUs based on geomorphology and soil captured the spatial variation of the abiotic component, detecting two main gradients: i) the coarse fragments and the rooting depth in the soil and ii) elevation and precipitation. The first gradient separates deep soils derived from volcanic ash (SDS, PD, LSSI, and HSSI) from shallow soils formed from basalt (LF) (Figure S1, Supplementary material). This gradient is an indicator of the water store and nutrients in the soil (Tetegan et al. 2015, Drobnik et al. 2018). Variation in elevation and soil properties related to topographic position have been reported as predictors of palm density and diversity in other tropical forests (Clark et al. 1995, Cámara-Leret et al. 2017). We found that elevation and its relationship with precipitation influenced the spatial distribution of the palm species. In this regard, Perrigo et al. (2020) and Svenning (2001a) indicate that in landscapes with significant elevation variation (e.g., > 100 m difference), as occurs in our study area, precipitation may vary even over small distances. The variation in this climatic variable was particularly relevant in LF, with lower soil water available for plants, since palm individuals could experience water stress during three months in which evapotranspiration exceeds rainfall. This finding may indicate that precipitation variations due to climate change can have the potential to affect the distribution of palm species at the landscape scale, especially those with lower soil AWHC.

Palm composition and density in the context of environmental heterogeneity. The ordination analysis showed that LUs are a good stratification approach to identify the effect of environmental heterogeneity on palm community composition in Los Tuxtlas reserve. Site distributions in the ordination reflected palm species distribution between units. This result is in line with other studies that proved the efficiency of sampling based on soil and topography to capture and explain spatial variation of tree, understory plants, and palm communities that are structured by species responses to environmental factors in volcanic Neotropical landscapes (Clark et al. 1995, Denslow et al. 2019, Navarrete Segueda et al. 2021).

The altitude and the slope had a strong effect on the presence of B. mexicana and C. ernesti-augusti and was evident from the decrease in the density of A. mexicanum and C. alternans in the higher elevation units (SDS and HSSI; Figure 3B). This result coincides with those of Cepeda-Cornejo & Dirzo (2010), for C. alternans and C. elatior that showed specific affinity with the LUs having lower precipitation and well-drained soils, either due to a low percentage of clay or high stoniness (LSSI and LF; Figure 3A). According to studies of Clark et al. (1995), Svenning (2001a), and Muscarella et al. (2019), this result may be because soil water conditions are regulator of the composition of TRF palm species. Variables that differentiated PD from LF-elevation and coarse fragments in soil (Figure S1, Supplementary material) influence the distribution of C. ernesti-augusti and C. tepejilote, which are practically absent in PD and relatively abundant in LF. Bacon & Bailey (2006) and Miceli-Méndez et al. (2013) also reported higher density of C. tepejilote and C. ernesti-augusti in sites with shallow and stony soils, respectively. Renninger et al. (2013) suggest that palms are more efficient in resource allocation and distribution of aerial and root biomass in the face of constraints such as shallow soils with low AWHC. This strategy has been shown by C. tepejilote (Oyama & Dirzo 1988), which may have an advantage over tree species, which have deeper roots. The allocation of biomass to the stem and roots has been reported as an important indicator of competitive interactions in the face of restrictive factors (Qi et al. 2019).

Four species were recorded in all the LUs (A. mexicanum, C. alternans, C. pinnatifrons and R. gracilis) (Figure 2). The wide mesoscale distribution of C. pinnatifrons has been reported by Svenning (2001b), who found that the density of this species is associated with deep, well-drained soils. This result contrasts with our findings, since soil depth was not a factor that influenced the density of this species. In the case of A. mexicanum, its presence of has been widely reported in this TRF (Popma et al. 1988, Ibarra-Manríquez et al. 1997, Hernández-Ruedas et al. 2018) and is usually related to its high seed production capacity in the face of disturbances (Martínez-Ramos et al. 2016). However, the hyperdominance of A. mexicanum decreased significantly in the highest elevation and rooting depth units (SDS and HSSI) (Figure 2), contrasting with the proposal of Martínez-Ramos et al. (2016). This divergence in results is likely attributable to the fact that the plots of Martínez-Ramos et al. (2016) were all located on the PD unit, which excluded the effect of the elevation on this species. Our results support the argument of Arroyo-Rodríguez & Melo (2016) that it is important evaluate palm populations throughout the landscape in the Sierra de Los Tuxtlas, since the soil and topographic position have a clear relationship with A. mexicanum distribution throughout the analyzed forest.

The low density of A. mexicanum in SDS could also be a consequence of the increase in light in the lower strata of the forest. According to Clark et al. (1996) and Alves et al. (2010), the lateral incidence of light at sites with abrupt relief can be equivalent to light levels in clearings or forest edges. This effect limits the density of palm species with low tolerance to prolonged light exposure, which includes A. mexicanum (Martínez-Ramos et al. 2016, Hernández-Ruedas et al. 2018). Reduced exposure to light would also explain the increase of individuals of this species in the LUs with less sloping relief (LF and PD).

As expected, slight variations in the tree community, together with the abiotic factors, influenced the palm community (Figure 3). This effect is significant only when soil-topography variables are included in the analysis (Table 3). Interspecific competition with the tree community could explain this relationship. The interaction between soil and tree community had a positive effect on palm density, especially in landscape unit whit high stoniness, where palm species were favored due to their more efficient use of sources, as soil water, compared to tree species. This result lends support to the interpretations of Martínez-Ramos et al. (2016) by indicating that a higher density of C. alternans, C. tepejilote, and A. mexicanum (Figure 3B) alters the composition and structure of tree community. However, our results indicate that the composition of the palm community is driven by environmental heterogeneity.

The increase in the populations of C. elatior and C. ernesti-augusti in SDS could be related to the decrease in the density of A. mexicanum and C. pinnatifrons. Particularly, it has been reported that A. mexicanum in the LTBS is a highly competitive species that affects density of another palm and tree species (Noriega‐Piña et al. 2021). Cepeda-Cornejo & Dirzo (2010) suggested that there is a microhabitat separation, with C. pinnatifrons distributed at low elevation sites, while C. ernesti-augusti typically occurs at slightly higher elevations. This relationship suggests that the spatial variation of the most abundant palm species, in the context of environmental heterogeneity, affects the density of other palms. This result is highly relevant because palm community is an important structural component in the understory layer (plants up to 10 m) (Ibarra-Manríquez 1988). However, further studies are required to evaluate the density-dependent effects on palm populations, considering the density of trees-relief-light interaction.

Our results show that there is high edaphic and topographic variation that influences the distribution of species at the landscape scale. However, geopedological heterogeneity has not been considered to explain the distribution patterns and species diversity in the TRF studied. To integrate the complex variation of both abiotic and biotic ecosystem components, we proved that landscape stratification into spatially explicit LUs is highly effective. We found that palm species density differs among LUs, and our study supports the idea that abiotic and biotic components simultaneously regulate the palm composition and abundance. However, the response to each factor differs among species. Therefore, further studies are needed to fill the information gap on these species’ functional attributes and ecological niches in the context of environmental heterogeneity. While the composition of adult palms was related mainly to abiotic factors at the landscape level, there are other factors such as biotic interactions (e.g., herbivory, pollination, seed dispersal, seed predation) that influence palm composition, and operate at different scales and stages of the life cycle of palms. The systematic integration of landscape elements allows us to address the multifactorial effect of the distribution and density of palm species. Therefore, under this approach, the analysis of other plant groups together (e.g., trees or lianas) can provide useful information to understand species coexistence in plant communities and critical factors for their conservation, as well as to understand the importance of environmental filters in the face of environmental changes and their effect on the distribution of species.

Supplementary material

Supplemental data for this article can be accessed here: https://doi.org/10.17129/botsci.3204

Supplemental material

Acknowledgements

A. Navarrete-Segueda thanks the support received from the Consejo Nacional de Humanidades Ciencias y Tecnologías through the 2022 (1) postdoctoral fellowship grant of the “Estancias Posdoctorales por México - Académica”. Santiago Sinaca-Colín, Santiago Xolo and Armando Xolo helped with the field work. The authors want to thank to Estación de Biología Tropical Los Tuxtlas, mainly Rosamond Coates, for all the facilities to realize field surveys. The authors declare that they have no conflict of interest.

Literature cited

Alves LF, Vieira SA, Scaranello MA, Camargo PB, Santos FAM, Carlos AJ, Martinelli LA. 2010. Forest structure and live aboveground biomass variation along an elevational gradient of tropical Atlantic moist forest (Brazil). Forest Ecology and Management 260: 679-691. DOI: https://doi.org/10.1016/j.foreco.2010.05.023 [ Links ]

Anderson M, Braak CT. 2003. Permutation tests for multi-factorial analysis of variance. Journal of Statistical Computation and Simulation 73: 85-113. DOI: https://doi.org/10.1080/00949650215733 [ Links ]

Anderson MJ. 2017. Permutational Multivariate Analysis of Variance (PERMANOVA). Wiley StatsRef: Statistics Reference Online 1-15. DOI: https://doi.org/10.1002/9781118445112.stat07841 [ Links ]

Arroyo-Rodríguez V, Melo FPL. 2016. Commentary: Anthropogenic disturbances jeopardize biodiversity conservation within tropical rainforest reserves. Frontiers in Ecology and Evolution 4: 73. DOI: https://doi.org/10.3389/fevo.2016.00073 [ Links ]

Bacon CD, Bailey CD. 2006. Taxonomy and conservation: A case study from Chamaedorea alternans. Annals of Botany 98: 755-763. DOI: https://doi.org/10.1093/aob/mcl158 [ Links ]

Baldeck CA, Harms KE, Yavitt JB, John R, Turner BL, Valencia R, Navarrete H, Davies SJ, Chuyong GB, Kenfack D, Thomas DW, Madawala S, Gunatilleke N, Gunatilleke S, Bunyavejchewin S, Kiratiprayoon S, Yaacob A, Supardi MNN, Dalling JW. 2013. Soil resources and topography shape local tree community structure in tropical forests. Proceedings of the Royal Society B: Biological Sciences 280: 20130548. DOI: https://doi.org/10.1098/rspb.2012.2532 [ Links ]

Barlow J, França F, Gardner TA, Hicks CC, Lennox GD, Berenguer E, Castello L, Economo EP, Ferreira J, Guénard B, Gontijo C, Isaac V, Lees AC, Parr CL, Wilson SK, Young PJ, Graham NAJ. 2018. The future of hyperdiverse tropical ecosystems. Nature 559: 517-526. DOI: https://doi.org/10.1038/s41586-018-0301-1 [ Links ]

Benítez-López A, Santini L, Schipper AM, Busana M, Huijbregts MAJ. 2019. Intact but empty forests? Patterns of hunting-induced mammal defaunation in the tropics. PLOS Biology 17: e3000247. DOI: https://doi.org/10.1371/journal.pbio.3000247 [ Links ]

Berenguer E, Gardner TA, Ferreira J, Aragão LEOC, Camargo PB, Cerri CE, Durigan M, Oliveira-Junior RC, Vieira ICG, Barlow J. 2015. Developing cost-effective field assessments of carbon stocks in human-modified tropical forests. Plos One 10: e0133139. DOI: https://doi.org/10.1371/journal.pone.0133139 [ Links ]

Bonilla-Moheno M, Aide TM. 2020. Beyond deforestation: Land cover transitions in Mexico. Agricultural Systems 178: 102734. DOI: https://doi.org/10.1016/j.agsy.2019.102734 [ Links ]

Brown S. 1997. Estimating biomass and biomass change of tropical forests: A primer. Rome: Food and Agriculture Organization of the United Nations. ISBN: 92-5-103955-0 [ Links ]

Brown S. 2002. Measuring carbon in forests: Current status and future challenges. Environmental Pollution 116: 363-372. DOI: https://doi.org/10.1016/S0269-7491(01)00212-3 [ Links ]

Brown C, Burslem DFRP, Illian JB, Bao L, Brockelman W, Cao M, Chang LW, Dattaraja HS, Davies S, Gunatilleke CVS, Gunatilleke IAUN, Huang J, Kassim AR, LaFrankie JV, Lian J, Lin L, Ma K, Mi X, Nathalang A, Noor S, Ong P, Sukumar R, Su SH, Sun IF, Suresh HS, Tan S, Thompson J, Uriarte M, Valencia R, Yap SL, Ye W, Law R. 2013. Multispecies coexistence of trees in tropical forests: Spatial signals of topographic niche differentiation increase with environmental heterogeneity. Proceedings of the Royal Society B: Biological 280: 20130502. DOI: https://doi.org/10.1098/rspb.2013.0502 [ Links ]

Bünemann EK, Bongiorno G, Bai Z, Creamer RE, De Deyn G, de Goede R, Fleskens L, Geissen V, Kuyper TW, Mäder P, Pulleman M, Sukkel W, van Groenigen JW, Brussaard L. 2018. Soil quality - A critical review. Soil Biology and Biochemistry 120: 105-125. DOI: https://doi.org/10.1016/j.soilbio.2018.01.030 [ Links ]

Cámara-Leret R, Tuomisto H, Ruokolainen K, Balslev H, Kristiansen M. 2017. Modelling responses of western Amazonian palms to soil nutrients. Journal of Ecology 105: 367-381. DOI: https://doi.org/10.1111/1365-2745.12708 [ Links ]

Cazalis V, Princé K, Mihoub J-B, Kelly J, Butchart SHM, Rodrigues ASL. 2020. Effectiveness of protected areas in conserving tropical forest birds. Nature Communications 11: 4461. DOI: https://doi.org/10.1038/s41467-020-18230-0 [ Links ]

Cepeda-Cornejo V, Dirzo R. 2010. Sex-Related differences in reproductive allocation, growth, defense and herbivory in three dioecious Neotropical palms. PLOS ONE 5: e9824. DOI: https://doi.org/10.1371/journal.pone.0009824 [ Links ]

Chen FW, Liu CW. 2012. Estimation of the spatial rainfall distribution using inverse distance weighting (IDW) in the middle of Taiwan. Paddy Water Environment 10: 209-222. DOI: https://doi.org/10.1007/s10333-012-0319-1 [ Links ]

Cintra R, Ximenes A de C, Gondim FR, Kropf MS. 2005. Forest spatial heterogeneity and palm richness, abundance and community composition in Terra Firme forest, Central Amazon. Revista Brasileira Botanica 28: 75-84. DOI: https://doi.org/10.1590/S0100-84042005000100007 [ Links ]

Clark DA, Clark DB, Sandoval RM, Castro MVC. 1995. Edaphic and human effects on landscape-scale distributions of tropical rain forest palms. Ecology 76: 2581-2594. DOI: https://doi.org/10.2307/2265829 [ Links ]

Clark DB, Clark DA, Rich PM, Weiss S, Oberbauer SF. 1996. Landscape-scale evaluation of understory light and canopy structures: Methods and application in a Neotropical lowland rain forest. Canadian Journal of Forest Research 26: 747-757. DOI: https://doi.org/10.1139/x26-084 [ Links ]

Crawley MJ. 2007. The R Book. Chichester: John Wiley & Sons Ltd. ISBN: 978-0-470-51024-7 [ Links ]

Dehn M, Gärtner H, Dikau R. 2001. Principles of semantic modeling of landform structures. Computational Geosciences 27: 1005-1010. DOI: https://doi.org/10.1016/S0098-3004(00)00138-2 [ Links ]

Denslow JS, Chaverri GL, Vargas O. 2019. Patterns in a species‐rich tropical understory plant community. Biotropica 51: 664-673. DOI: https://doi.org/10.1111/btp.12704 [ Links ]

Dirzo R, Miranda A. 1991. Altered patterns of herbivory and diversity in the forest understory: A case study of the possible consequences of contemporary defaunation. In: Price PW, Lewinsohn TM, Fernandes GW, Benson WW, eds. Plant-animal interactions: Evolutionary ecology in tropical and temperate regions. New York: Wiley and Sons, pp. 273-287. DOI: https://doi.org/10.1086/417402 [ Links ]

Dirzo R, Young HS, Galetti M, Ceballos G, Isaac NJB, Collen B. 2014. Defaunation in the Anthropocene. Science 345: 401-406. DOI: https://doi.org/10.1126/science.1251817 [ Links ]

Dray S, Dufour A-B. 2007. The ade4 package: Implementing the duality diagram for ecologists. Journal of Statistical Software 22: 1-20. DOI: https://doi.org/10.18637/jss.v022.i04 [ Links ]

Drobnik T, Greiner L, Keller A, Grêt-Regamey A. 2018. Soil quality indicators-From soil functions to ecosystem services. Ecological Indicators 94: 151-169. DOI: https://doi.org/10.1016/j.ecolind.2018.06.052 [ Links ]

Eiserhardt WL, Svenning J-C, Kissling WD, Balslev H. 2011. Geographical ecology of the palms (Arecaceae): Determinants of diversity and distributions across spatial scales. Annals of Botany 108: 1391-1416. DOI: https://doi.org/10.1093/aob/mcr146 [ Links ]

Eckelmann W, Sponagel H, Grottenthaler W. 2005. Bodenkundliche Kartieranleitung. KA5. Verbesserte und erweiterte-Auflage. 5th ed. Hannover: Bundesanstalt für geowissenschaften und rohstoffe in zusammenarbeit mit den staatlichen geologischen diensten. ISBN 978-3-510-95920-4 [ Links ]

Emilio T, Quesada CA, Costa FRC, Magnusson WE, Schietti J, Feldpausch TR, Brienen RJW, Baker TR, Chave J, Álvarez E, Araújo A, Bánki O, Castilho CV, Honorio EN, Killeen TJ, Malhi Y, Oblitas-Mendoza EM, Monteagudo A, Neill D, Alexander-Parada G, Peña-Cruz A, Ramirez-Angulo H, Schwarz M, Silveira M, ter Steege H, Terborgh JW, Thomas R, Torres-Lezama A, Vilanova E, Phillips, OL. 2014. Soil physical conditions limit palm and tree basal area in Amazonian forests. Plant Ecology & Diversity 7: 215-229. DOI: https://doi.org/10.1080/17550874.2013.772257 [ Links ]

FAO [Food and Agriculture Organization of the United Nations], UNEP [United Nations Environmental Programme]. 2020. The state of the world’s forests 2020. Forests, biodiversity and people. Rome: Food and Agriculture Organization of the United Nations. ISBN 978-92-5-132419-6. DOI: https://doi.org/10.4060/ca8642en [ Links ]

Gardner CJ, Bicknell JE, Baldwin-Cantello W, Struebig MJ, Davies ZG. 2019. Quantifying the impacts of defaunation on natural forest regeneration in a global meta-analysis. Nature Communications 10: 4590. DOI: https://doi.org/10.1038/s41467-019-12539-1 [ Links ]

Gutiérrez-García G, Ricker M. 2011. Climate and climate change in the region of Los Tuxtlas (Veracruz, Mexico): A statistical analysis. Atmosfera 24: 347-373. [ Links ]

Hernández-Ruedas MA, Arroyo-Rodríguez V, Morante-Filho JC, Meave JA, Martínez-Ramos M. 2018. Fragmentation and matrix contrast favor understory plants through negative cascading effects on a strong competitor palm. Ecological Applications 28: 1546-1553. DOI: https://doi.org/10.1002/eap.1740 [ Links ]

Ibarra-Manríquez G. 1988. The palms of a tropical rain forest in Veracruz, México. Principes 32: 147-155. [ Links ]

Ibarra-Manríquez G, Sinaca-Colín S. 1995. Lista florística comentada de la Estación de Biología Tropical "Los Tuxtlas", Veracruz, México. Revista de Biología Tropical 43: 75-115. [ Links ]

Ibarra-Manríquez G, Martínez Ramos M, Dirzo R, Núñez-Farfán J. 1997. La vegetación. In: González-Soriano E, Dirzo R, Vogt R, eds. Historia Natural de Los Tuxtlas. Ciudad de México: Universidad Nacional Autónoma de México, pp. 62-85. ISBN: 968-36-5646-3 [ Links ]

INEGI [Instituto Nacional de Estadística y Geografía]. 2010. Conjuntos de datos vectoriales de la carta topográfica escala 1:20,000. [ Links ]

Izsák J. 2006. Some practical aspects of fitting and testing the Zipf-Mandelbrot model. Scientometrics 67: 107-120. DOI: https://doi.org/10.1007/s11192-006-0052-x [ Links ]

Jost L. 2006. Entropy and diversity. Oikos 113: 363-375. DOI: https://doi.org/10.1111/j.2006.0030-1299.14714.x [ Links ]

Laurance WF, Fearnside PM, Laurance SG, Delamonica P, Lovejoy TE, Rankin-de Merona JM, Chambers JQ, Gascon C. 1999. Relationship between soils and Amazon forest biomass: A landscape-scale study. Forest Ecology and Management 118: 127-138. DOI: https://doi.org/10.1016/S0378-1127(98)00494-0 [ Links ]

Li Y, Shipley B, Price JN, Dantas V de L, Tamme R, Westoby M, Siefert A, Schamp BS, Spasojevic MJ, Jung V, Laughlin DC, Richardson SJ, Bagousse-Pinguet YLe, Schöb C, Gazol A, Prentice HC, Gross N, Overton J, Cianciaruso MV, Louault F, Kamiyama C, Nakashizuka T, Hikosaka K, Sasaki T, Katabuchi M, Frenette-Dussault C, Gaucherand S, Chen N, Vandewalle M, Batalha MA. 2018. Habitat filtering determines the functional niche occupancy of plant communities worldwide. Journal of Ecology 106: 1001-1009. DOI: https://doi.org/10.1111/1365-2745.12802 [ Links ]

MacMillan RA, Shary PA. 2009. Landforms and landform elements in geomorphometry. In: Hengl T, Reuter HI, eds. Geomorphometry: Concepts, Software, Applications. The Netherlands: Elsevier, pp. 227-254. DOI: https://doi.org/10.1016/S0166-2481(08)00009-3 [ Links ]

McCune B, Grace JB. 2002. Analysis of ecological communities. 2nd ed. Gleneden Beach: MjM Software Design. ISBN: 0-9721290-0-6 [ Links ]

Magurran AE. 2004. Measuring biological diversity. Oxford: Blackwell Publishing. ISBN: 0-632-05633-9 [ Links ]

Martínez-Ramos M, Álvarez-Buylla E, Sarukhán J, Piñero D. 1988. Treefall age determination and gap dynamics in a tropical forest. Journal of Ecology 76: 700-716. DOI: https://doi.org/10.2307/2260568 [ Links ]

Martínez-Ramos M, Ortiz-Rodríguez IA, Piñero D, Dirzo R, Sarukhán J. 2016. Anthropogenic disturbances jeopardize biodiversity conservation within tropical rainforest reserves. Proceedings of the National Academy of Sciences 113: 5323-5328. DOI: https://doi.org/10.1073/pnas.1602893113 [ Links ]

Mendoza A, Franco M. 1998. Sexual reproduction and clonal growth in Reinhardtia gracilis (Palmae), an understory tropical palm. American Journal of Botany 85: 521-527. DOI: https://doi.org/10.2307/2446436 [ Links ]

Miceli Méndez CL, Sánchez Molina DF, López Mendoza S, Reyes Escutia F de J. 2013. Palma Cola de Pescado (Chamaedorea ernesti-augusti). 1st ed. Tuxtla Gutiérrez: Universidad de Ciencias y Artes de Chiapas. ISBN: 978-607-8240-11-1 [ Links ]

Muscarella R, Bacon CD, Faurby S, Antonelli A, Kristiansen SM, Svenning JC, Balslev H. 2019. Soil fertility and flood regime are correlated with phylogenetic structure of Amazonian palm communities. Annals of Botany 123: 641-655. DOI: https://doi.org/10.1093/aob/mcy196 [ Links ]

Navarrete-Segueda A, Martínez-Ramos M, Ibarra-Manríquez G, Cortés-Flores J, Vázquez-Selem L, Siebe C. 2017. Availability and species diversity of forest products in a Neotropical rainforest landscape. Forest Ecology and Management 406: 242-250. DOI: https://doi.org/10.1016/j.foreco.2017.08.037 [ Links ]

Navarrete-Segueda A, Cortés-Flores J, Cornejo-Tenorio G, González-Arqueros ML, Torres-García M, Ibarra-Manríquez G. 2021. Timber and non-timber forest products in the northernmost Neotropical rainforest: Ecological factors unravel their landscape distribution. Journal of Environmental Management 279: 111819. DOI: https://doi.org/10.1016/j.jenvman.2020.111819 [ Links ]

Nelson S, González-Caver E. 1992. Geology and K-Ar dating of the Tuxtla Volcanic Field, Veracruz, Mexico. Volcanology 55: 85-96. DOI: https://doi.org/10.1007/BF00301122 [ Links ]

Noriega‐Piña K, Piñero D, Valverde T, Martínez‐Ramos M. 2021. Competitive effects of a dominant palm on sapling performance in a Neotropical rainforest. Biotropica 53: 1558-1568. DOI: https://doi.org/10.1111/btp.13002 [ Links ]

Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Szoecs E, Wagner H. 2019. Vegan: Community Ecology Package. R Package. version 2.5-6 296. [ Links ]

Oyama K, Dirzo R. 1988. Biomass allocation in the dioecious tropical palm Chamaedorea tepejilote and its life history consequences. Plant Species Biology 3: 27-33. DOI: https://doi.org/10.1111/j.1442-1984.1988.tb00168.x [ Links ]

Oyama K, Dirzo R, Ibarra-Manríquez G. 1992. Population structure of the dominant palms species in the understory of a mexican lowland rain forest. Tropics 2: 23-28. DOI: https://doi.org/10.3759/tropics.2.23 [ Links ]

Perrigo A, Hoorn C, Antonelli A. 2020. Why mountains matter for biodiversity. Journal of Biogeography 47: 315-325. DOI: https://doi.org/10.1111/jbi.13731 [ Links ]

Piñero D, Martínez-Ramos M, Sarukhán J. 1984. A population model of Astrocaryum mexicanum and a sensitivity analysis of its finite rate of increase. Journal of Ecology 72: 977-991. DOI: https://doi.org/10.2307/2259545 [ Links ]

Popma J, Bongers F, Martínez-Ramos M, Veneklaas E. 1988. Pioneer species distribution in treefall gaps in Neotropical rain forest: A gap definition and its consequences. Journal of Tropical Ecology 4: 77-88. DOI: https://doi.org/10.1017/S0266467400002534 [ Links ]

Qi Y, Wei W, Chen C, Chen L. 2019. Plant root-shoot biomass allocation over diverse biomes: A global synthesis. Global Ecology and Conservation 18: e00606. DOI: https://doi.org/10.1016/j.gecco.2019.e00606 [ Links ]

R Core Team. 2022. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.orgLinks ]

Raven PH, Gereau RE, Phillipson PB, Chatelain C, Jenkins CN, Ulloa Ulloa C. 2020. The distribution of biodiversity richness in the tropics. Science Advances 6: eabc6228. DOI: https://doi.org/10.1126/sciadv.abc6228 [ Links ]

Renninger HJ, McCulloh KA, Phillips N. 2013. A comparison of the hydraulic efficiency of a palm species (Iriartea deltoidea) with other wood types. Tree Physiology 33: 152-160. DOI: https://doi.org/10.1093/treephys/tps123 [ Links ]

Rodrigues AC, Villa PM, Neri AV. 2019. Fine-scale topography shape richness, community composition, stem, and biomass hyperdominant species in Brazilian Atlantic Forest. Ecological Indicators 102: 208-217. DOI: https://doi.org/10.1016/j.ecolind.2019.02.033 [ Links ]

Rodrigues LF, Cintra R, Castilho CV, de Sousa Pereira O, Pimentel TP. 2014. Influences of forest structure and landscape features on spatial variation in species composition in a palm community in central Amazonia. Journal of Tropical Ecology 30: 565-578. DOI: https://doi.org/10.1017/S0266467414000431 [ Links ]

Rügnitz MT, Chacón ML, Porro R. 2008. Guía para la determinación de carbono en pequeñas propiedades rurales. Lima: Centro Mundial Agroflorestal (ICRAF) y Consórcio Iniciativa Amazônica (IA). ISBN: 978-92-9059-254-9 [ Links ]

Schmitt S, Tysklind N, Derroire G, Heuertz M, Hérault B. 2021. Topography shapes the local coexistence of tree species within species complexes of Neotropical forests. Oecologia 196: 389-398. DOI: https://doi.org/10.1007/s00442-021-04939-2 [ Links ]

Schoeneberger PJ, Wysocki DA, Benham EC, Soil Survey Staff. 2012. Fieldbook for describing and sampling soils. Version 3.0. Lincoln: Natural Resources Conservation Service, National Soil Survey Center. IBSN: 9780160915420 [ Links ]

SMN [Servicio Metereológico Nacional] . 2010. Información climatológica por estado. Serv. Metereológico Nac. https://shre.ink/QD6k (accessed March 25, marzo 2019). [ Links ]

Siebe C, Jahn R, Stahr K. 1996. Manual para la descripción y evaluación ecológica de suelos en el campo. Chapingo: Sociedad Mexicana de la Ciencia del Suelo. ISBN: 9686201189 [ Links ]

Svenning JC. 2001a. On the role of microenvironmental heterogeneity in the ecology and diversification of Neotropical rain-forest palms (Arecaceae). The Botanical Review 67: 1-53. DOI: https://doi.org/10.1007/BF02857848 [ Links ]

Svenning JC. 2001b. Environmental heterogeneity, recruitment limitation and the mesoscale distribution of palms in a tropical montane rain forest (Maquipucuna, Ecuador). Journal of Tropical Ecology 17: 97-113. DOI: https://doi.org/10.1017/S0266467401001067 [ Links ]

Tetegan M, Richer-de Forges AC, Verbeque B, Nicoullaud B, Desbourdes C, Bouthier A, Arrouays D, Cousin I. 2015. The effect of soil stoniness on the estimation of water retention properties of soils: A case study from central France. Catena 129: 95-102. DOI: https://doi.org/10.1016/j.catena.2015.03.008 [ Links ]

Tylianakis JM, Didham RK, Bascompte J, Wardle DA. 2008. Global change and species interactions in terrestrial ecosystems. Ecology Letter 11: 1351-1363. DOI: https://doi.org/10.1111/j.1461-0248.2008.01250.x [ Links ]

Valverde J, da Silva-Carvalho C, Jordano P, Galetti M. 2021. Large herbivores regulate the spatial recruitment of a hyperdominant Neotropical palm. Biotropica 53: 286-295. DOI: https://doi.org/10.1111/btp.12873 [ Links ]

Verma SP, Salazar VA, Negendank JFW, Milán M, Navarro LI, Besch T. 1993. Características petrográficas y geoquímicas de elementos mayores del Campo Volcánico de Los Tuxtlas. Geofísica Internacional 32: 237-248. DOI: https://doi.org/10.22201/igeof.00167169p.1993.32.2.558 [ Links ]

Villar N, Siqueira T, Zipparro V, Farah F, Schmaedecke G, Hortenci L, Brocardo CR, Jordano P, Galetti M. 2020. The cryptic regulation of diversity by functionally complementary large tropical forest herbivores. Journal of Ecology 108: 279-290. DOI: https://doi.org/10.1111/1365-2745.13257 [ Links ]

Von Thaden JJ, Laborde J, Guevara S, Mokondoko-Delgadillo P. 2020. Dinámica de los cambios en el uso del suelo y cobertura vegetal en la Reserva de la Biosfera Los Tuxtlas (2006-2016). Revista Mexicana de Biodiversidad 91: 913190. DOI: https://doi.org/10.22201/ib.20078706e.2020.91.3190 [ Links ]

Winkler K, Fuchs R, Rounsevell M, Herold M. 2021. Global land use changes are four times greater than previously estimated. Nature Communications 12: 2501. DOI: https://doi.org/10.1038/s41467-021-22702-2 [ Links ]

Zinck JA, Metternicht G, Bocco G, Del Valle HF. 2016. Geopedology. Cham: Springer. ISBN: 978-3-319-19158-4. DOI: https://doi.org/10.1007/978-3-319-19159-1 [ Links ]

Zonneveld IS. 1989. The land unit -A fundamental concept in landscape ecology, and its applications. Landscape Ecology 3: 67-86. DOI: https://doi.org/10.1007/BF00131171 [ Links ]

Received: June 21, 2022; Accepted: October 19, 2022; Published: June 05, 2023

*Author for correspondence: anavarretesegueda@gmail.com

Associate editor: Alejandro Hurtado Zavala

Author contributions: K. V. Miranda-Gallegos, writing-original draft, field work, formal analysis; A. Navarrete-Segueda, writing-original draft, writing-review and editing, field work, formal analysis, supervision and project administration; J. Cortés-Flores, writing-review and editing, formal analysis; M. L. González-Arqueros, writing-review and editing, formal analysis; E. E. Acosta-Pérez, field work; G. Ibarra-Manríquez, writing-review and editing, supervision and project administration.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License