Study Area

This study was conducted at the La Mancha Coastal Research Centre (CICOLMA, from its name in Spanish: Centro de Investigaciones Costeras La Mancha) located in the state of Veracruz, Mexico (19°35′50″ N; 96°22′45″ W; Figure 1). Climate is AW₂ with a mean annual temperature of 25.6℃ and a mean annual precipitation of 1,286 mm/year. There are two seasons based on the precipitation regime: a rainy season from June to September when 78% of the annual precipitation falls and a dry season from October to May, within which, from November to February there is a relatively humid and cold period characterized by strong northeasterly winds, locally known as nortes (Castillo-Campos & Travieso-Bello, 2006).

Figure 1.

Sampling plot location in Centro de Investigaciones Costeras La Mancha, municipality of Actopan, Veracruz, Mexico. The figure shows the polygons for each site sampled: secondary forest, dune scrub, and protected forest.

For this study, we selected three locations (Figure 1) that correspond to distinct types of habitats or landscape units (hereafter called “sites”) that differ in their physiognomy and floristic composition. All three are located in CICOLMA, which is a protected area (70 ha): (a) a protected tropical semideciduous forest (PF; protected forest); (b) a SF, and (c) coastal DS. The semideciduous forest (selva mediana subcaducifolia, sensu Miranda & Hernandez, 1963) is a type of seasonally dry tropical forest (sensu Pennington, Lavin, & Oliveira-Filho, 2009), in this case an old-growth forest that is more than 200 years old and has been protected since 1977 (Castillo-Campos & Travieso-Bello, 2006). The SF used to be a man-made pasture 5.9 ha in area that was abandoned at the end of 1995, when it was incorporated into CICOLMA. This pasture was used as a single management unit (i.e., never subdivided), where the African grass Panicum maximum Jacq. was planted to raise cattle for at least 20 years before it was abandoned. After 17 years of succession, this SF has a closed canopy that is 10 to 15 m tall (Guerrero-Reyna, 2016). Within CICOLMA, an area of 24.5 ha is covered by a sand dune system with coastal scrub vegetation. Cover is predominantly grasses with small, scattered patches of shrubs and trees growing in the dune slacks (Moreno-Casasola & Travieso-Bello, 2006).

The scrub site is part of a large and relatively recent parabolic semistabilized sand dune, with an estimated age of no more than a few centuries (Kellman & Roulet, 1990). It is 10 to 15 m and up to 20 m higher in elevation than the surrounding areas, where the other two sampling sites (SF and PF) are located. The SF and PF are located on an ancient dune system that is much older as it is late glacial in origin, and they have more fertile soils with a higher organic matter content. The soil is 80% sand, while the younger and the poorer soils of the DS are >90% sand. However, in the deepest dune slacks of the DS where there are dense woody patches, the soil properties are more similar to those of the PF and SF than to those of the surrounding dune slopes and crests (Suárez, Campos, & Cruz, 2015). Several authors (see review by Moreno-Casasola & Travieso-Bello, 2006) have proposed that the DS is undergoing a long-term primary succession process; based on their review, the current patches of scrub vegetation would be expected to expand by nucleation and ultimately become a semideciduous tropical forest.

Vegetation Sampling

A 6 ha plot (200 × 300 m) was delimited in each of the three sites, and within each of these plots, thirty 5 × 5 m quadrats were placed randomly. All woody plants rooted within the quadrats with a stem diameter < 5.0 cm and a height ≥ 20 cm were recorded. Don David Díaz, a local eco-guide, has worked with plant ecologists and botanists at CICOLMA for many years and assisted with the species identification of seedlings and saplings in the field. Collected specimens were also shown to botanists at the XAL herbarium, including Dr. Gonzalo Castillo-Campos, head of the Flora of Veracruz project (Instituto de Ecología, A.C. [INECOL], 2016), who is an expert on the local woody flora. Taxonomic nomenclature follows The Plant List (2013).

The percentage of tree canopy cover and leaf area index (LAI) above each quadrat was estimated using the Digital Plant Canopy Imager CI-110 ( http://www.cid-inc.com/ci-110.php), which takes hemispherical photos and processes each digital image with specialized software. For each 25 m² quadrat, the slope of the terrain was measured with an inclinometer and percentage ground cover was visually estimated for grasses, dicotyledonous herbs, leaf litter, and bare ground. Additionally, a 10 × 10 m quadrat with the 5 × 5 m quadrat at its center was set up to record, identify, and measure the diameter at breast height (DBH) (cm) of all rooted woody plants with a DBH ≥ 5 cm, which we regarded as adults (see Norden et al., 2009b). These were classified as local variables in further analyses and are considered proxies for the microenvironmental conditions in each quadrat that could affect the AR.

Data Analysis

Total species richness and abundance per site and their respective means per quadrat were estimated. We compared richness and abundance per quadrat among sites using generalized linear models, with a negative binomial error type because of the overly dispersed nature of our data. Then, when significant differences were found, a post hoc contrast test (Tukey) was used to determine which sites were different. This analysis was run in the statistical programming language R (version 2.14.10; R Development Core Team, 2015). The diversity of both the AR and adult communities in each site was analyzed by way of diversity profiles, estimating Hill numbers for observed richness (q = 0), for typical diversity (q = 1), and for the diversity of the most abundant species (q = 2), using the R software package iNEXT (Hsieh et al., 2016).

The floristic composition of the AR and its variation within and among sites was analyzed with a principal components analysis (PCA) ordination run in PC-ORD, version 6 (McCune & Grace, 2002). For each species, abundance data per quadrat was also log-transformed (log + 1). Quadrats with no woody plants inside (12 from the DS) were excluded from the ordination. The option “variance/covariance” in PC-ORD was used to emphasize variation in species composition. The species ordination (PCA) was complemented with a correlation analysis of the PCA scores on each axis and the local variables of each quadrat mentioned earlier. A cutoff value of r² > .2 was used to decide which local variables to include in the ordination graph or biplot. The PCA ordination of samples and important (i.e., r² > .2) local variables are shown on the same figure, drawing local variables as vectors, following McCune and Grace (2002) and Peck (2010). Also, a variance partitioning analysis of the AR species composition was run using the local variables grouped into four categories as explanatory matrices: topography (slope), adults (richness, abundance, and biomass), canopy (cover and LAI), and ground cover (percentages of grass, herbs, litter, and bare ground). This analysis was run in R with the Vegan library (varpart function; Oksanen et al., 2017) for the three sites and for each one independently. Additionally, the Jaccard similarity index was used to compare both the adult and AR communities among the three sampling sites.

For every species recorded in this study, its dispersal syndrome was determined based on fruit and seed morphology and an extensive review of the published literature (mainly INECOL, 2016, Castillo-Campos & Travieso-Bello, 2006). There were only three categories: endozoochory (species with fleshy fruit dispersed by vertebrate frugivores, mainly birds and bats), anemochory (species dispersed by wind), and other types of dispersal (including: barochory, dispersal by ungulate herbivores—mainly domestic animals, and unknown dispersal methods). Growth form was also assessed based mainly on the maximum height attained by each species and fell into four categories: (a) tall trees (adults taller than 15 m); (b) medium-sized trees (taller than 8 m but shorter than 15 m in height); (c) small trees, treelets, and shrubs (adults less than 8 m in height, regardless of their branching pattern); and (d) lianas. For simplicity, palm species were regarded as trees. We compared the variation in the percentage of growth forms and dispersal syndromes among sites for both the AR and adult communities in terms of richness and abundance.

Richness, Density, and Composition of the Advanced Regeneration

The site with the richest AR community was the SF with 67 morphospecies, followed by the PF with 55, and then the DS with only 27 species. Mean species richness (χ² test, p < .001) and mean abundance per quadrat (χ² test, p < .001) were significantly different among sites. The post hoc test (Tukey test, p < .05) confirmed that the three sites were significantly different in both variables (Table 1); richness was highest in SF followed by PF and then DS, while abundance was highest in PF, followed by SF and then DS.

Table 1.

Species Richness and Abundance per Quadrat (25 m²) in the Advance Regeneration Stage for Each Site.

The PCA ordination explained 65% of the total variance in floristic composition (36% along Axis 1 and 29% along Axis 2). Axis 1 separates the quadrats of the PF (highest values on the right) from those of the SF and the DS (Figure 2). The species whose presence and abundance per quadrat were most highly correlated (Pearson’s correlation) with the ordination scores for Axis 1 were B. alicastrum (r = .94), N. salicifolia (.71), Leucaena leucocephala (−.32), and C. corymbosa (−.33). The first two species had higher abundances and were recorded in more quadrats within the PF, while the other two were more abundant in the DS and the SF. Along Axis 2, the SF quadrats had the highest variation in ordination scores. Most of the PF quadrats are in the middle of this axis, while the DS quadrats are located at the bottom (Figure 2). The species with the highest correlations along Axis 2 were B. macrocarpa (.70), C. corymbosa (.69), C. uragoga (.71), and P. erythrocarpa (.70). Of all our local variables (i.e., environmental proxies and adult vegetation attributes) recorded per quadrat, only three were significantly correlated with the PCA scores from the species ordination: grass cover (%), litter ground cover per quadrat (%), and woody plant canopy cover (%) above each quadrat. Grass cover was negatively correlated with PCA scores along both ordination axes. Most of the DS quadrats and some SF quadrats had the highest values of grass cover (lower left corner of the plot). The woody canopy cover above each quadrat was positively correlated with PCA scores on both axes (upper right in the ordination plot; Figure 2).

Figure 2.

Principal components analysis (PCA) ordination of sample units from each site (PF = protected forest; DS = dune scrub; SF = secondary forest) and local variables shown as vectors that had the highest correlation with PCA scores (r² > .2 with at least one of the axes, following Peck, 2010); G = grass cover; L = leaf litter cover; C = canopy cover above quadrat tree (see Methods section).

Regarding the variance partition analysis of AR composition among and within sites, we found that the variation explained by all of the local variables together was highest for the DS (.821), followed by PF (.343; see online Appendix B). For the SF, the variables combined as adult attributes and those combined as ground cover type accounted for most of the variation, although with relatively low explanatory power (.210). For the PF, the combination of topography, adults, and canopy were the most explanatory (.322); and in the DS, it was the combination of adults, canopy, and cover type (.813).

Comparison of Adult Vegetation and AR

Adult vegetation was richest in the PF (53 spp.) followed by the SF (40 spp.) and the DS, which had less than one third of the species (13 spp.) recorded in the first two sites. In the SF, adult vegetation had the highest density (890.0 ind./ha) and basal area (23.3 m²/ha), with similar but lower values in the PF (670.0 ind./ha and 20.9 m²/ha) and much lower values in the DS (173.3 ind./ha and 4.4 m²/ha), where nearly 60% of its area (6 ha) was covered by grasses and forbs with no woody plant cover.

Plant diversity was much lower in the DS for both the AR and the adult vegetation in comparison with the PF and the SF (Figure 3). Species accumulation curves show that the AR of the 17-year-old SF had a slightly higher richness than the PF did (Figure 3(a)), while the opposite was true when comparing the community of adults (i.e., plants with DBH ≥ 5 cm) between the two sites (Figure 3(c)). However, the estimated 95% confidence intervals (CI) of these curves as well as the diversity profiles (Hill numbers) suggest that there was no statistical difference in the observed richness (q0 ± 95% CI) between these sites for either the AR or adult vegetation (Figure 3(b) and (d)). Nevertheless, the numbers of typical (q1) and very abundant (q2) species in the AR were higher in the SF (Figure 3(b); q1 = 43.2 spp. ± 4.1, CI 95%; q2 = 31.7 ± 3.8) than in the PF (q1 = 29.7 ± 4.4; q2 = 17.7 ± 2.5), where just two species, B. alicastrum and N. salicifolia, were notoriously abundant (70.5% and 11.4% of AR abundance, respectively). In the adult vegetation, diversity was very similar between the SF and the PF, with around 30 typical species (q1) and around 21 very abundant species (q2) in both sites (Figure 3(d)).

Figure 3.

Species accumulation curves for the advance regeneration (a) and the adult community (c) in each of the sites sampled: dune scrub (DS), protected forest (PF), and secondary forest (SF) in Centro de Investigaciones Costeras La Mancha, Veracruz, Mexico, and their diversity profiles (b and d, respectively) showing Hill’s numbers (q0 = observed richness, q1 = number of typical species, q2 = dominant species) and their 95% confidence intervals.

Overall, pairwise comparisons between sites for both of the stages sampled (AR and adult) had low similarity values, with the highest recorded for the AR between the SF and the PF (.341), which had the highest number of shared species (31) between them (Table 2), representing almost half of the total number of species recorded in the SF and more than half of those recorded in the PF. The latter was followed closely by the similarity between the AR of the SF and that of the DS (.306). The similarity of the AR to the adults of the SF (.277) and the adults of the PF (.270) was also comparatively high. The least similar samples were the AR of the PF and the adults of the DS, which only shared five species. It is noteworthy that similarities among sites were generally higher for AR than for adults and that the AR of a given site was less similar to the adults growing in the same site than with the AR of the other two sites.

Table 2.

Similarity Between Advance Regeneration and Adult Communities Within and Among the Three Habitats Sampled.

Regarding dispersal syndrome, the most prevalent type of dispersal in the three sites and for both the AR and the adult vegetation was endozoochorous, followed by anemochorous dispersal (Figure 4). A relatively high proportion of adult plants was anemochorous in the SF and the DS; while in the AR of both of these sites, there was a much higher abundance of endozoochorous seedlings and saplings. Regarding the growth form of the woody plants, the absence of lianas in the DS is notable, but they were recorded in the SF and PF as adults, saplings, and seedlings, although their richness and abundance were low in both sites. Small trees (treelets) and shrubs, which in this study were grouped into a single growth form category (see Methods section), were the most prevalent type of species in the three sites for both AR and adults and accounted for a greater proportion of species in AR than in the adults. In the adult vegetation of the PF and the SF, the most abundant growth form was medium-sized trees, followed closely by tall trees in the PF and by shrubs and treelets in the SF. In the AR of the SF and the DS, shrubs and treelets were the most abundant growth form (Figure 4(d)); while in the PF, tall trees were by far the most abundant, mainly due to the high abundance of B. alicastrum seedlings in this site.

Figure 4.

Proportional richness (% species) and abundance (% individuals) of the advance regeneration (Adv. Regen.) and adult vegetation (adults), grouped by dispersal syndrome (left panels) and growth form (right panels), for the dune scrub (DS), secondary forest (SF), and the protected forest (PF). Growth forms: tall trees > 15 m in height; medium trees 8 to 15 m; treelets < 8 m and shrubs; and lianas. Dispersal syndrome: endozoochorous species (Zoo); anemochorous (Ane); and other type of dispersal (see Methods section).

Successional Trajectories Deduced From the AR Stage

The analysis of species accumulation curves suggests that differences in the observed richness of adult woody plants (DBH ≥5 cm in this study) were not statistically significant between the secondary and the preserved forests. Nevertheless, diversity profiles show that the AR in the SF was more diverse than it was in the preserved forest. This is a remarkable result, since the SF site was used continuously for raising cattle 20 years prior to its abandonment, yet 17 years after the exclusion of cattle, not only is its richness in adults comparable to that of the preserved forest but its AR is also more diverse than that of the preserved forest owing to the mixture of early, mid- and late successional species in the site. The main difference between the SF and the preserved forest was their floristic composition, most notably in the dominant species of adult plants.

It has been pointed out that forest recovery within seasonally dry tropical regions should be much slower than it is in more humid regions (Janzen, 1988; and see reviews in Quesada et al., 2009, and Poorter et al., 2016). Janzen (1988) also mentioned that abandoned pastures in the dry tropics could be colonized and dominated over several years by a few species of fast growing, wind dispersed trees, which would impede or slow down the establishment of tree species dispersed by forest animals that usually avoid open areas. Consequently, secondary succession could be deflected toward a poorer composition that is different from that of the original forest. The fast recovery of species richness does not imply the complete recovery of composition (see Barlow et al., 2007, for a case of quick recovery in richness but not in composition in the wet tropics).

If we look at the large adult plants of the 17-year-old SF, we can see that the closed secondary canopy is dominated by just a few wind dispersed species (Cedrela, Tecoma, Diphysa, and Gliricidia) and some that are dispersed by livestock (Enterolobium, Guazuma, and Leucaena). Bursera simaruba, the only endozoochorous tree (dispersed by birds) that could be regarded as dominant in the secondary canopy, is widely planted as a living fence in the area, so there were plenty of sources of this species’ seeds nearby when the pasture was abandoned. These observations would lead us to believe that secondary succession in the dry tropics is very slow and also that a first wave of a few colonizer species that do not depend on frugivorous fauna for dispersal could be stopping or at least slowing down the establishment of late successional species dispersed by birds and bats. However, when we look at the smaller woody plants (DBH < 5 cm) that are part of the AR stage in this SF, two things stand out: there are very few or no immature individuals of the first colonizers of the site in its understory, and yet there are plenty of seedlings, saplings, and juveniles of several bird- and bat-dispersed species (e.g., B. alicastrum, N. salicifolia, Ehretia tinifolia, and Diospyros acapulcensis, among others) growing under the light shade of the current secondary tree canopy. This explains in part why several quadrats of the SF had a high degree of similarity with those of the primary forest (Figure 2), and also why overall similarity among sites (Table 2) was highest in the AR between these two sites. Shade tolerant species that are already present in the AR but have slow growth rates and long life cycles will be recruited as adults in the future, increasing the similarity of its floristic composition to that of the nearby preserved forest. SF in the dry tropics may appear to be dominated by a small subset of fast growing, easily dispersed species of early secondary succession phases, but a careful examination of the understory reveals that mid- to late successional species have already established there and their growth is facilitated by the light shade of the first successional wave of early colonizers (see also Lebrija-Trejos, Perez-García, Meave, Bongers, & Poorter, 2010b; Norden et al., 2009b).

These trends were evident on just over 5 ha of the large 6 ha plot of SF that we studied; however, within the same plot, there was a compact patch of almost 1 ha that had a very different structure and composition. This patch was almost exclusively dominated by L. leucocephala, whose short, slender trees (3–4 m tall, DBH < 8 cm) formed a thin, deciduous tree canopy, under which the African grass P. maximum was still growing profusely, covering 100% of the ground. The only other tree species in this patch were Diphysa americana, Tecoma stans, and B. simaruba, whose stunted and widely scattered individuals were indicative of the stressful conditions for woody plants within this patch. If our vegetation sampling had been restricted to this patch, the most likely conclusion would have been that both the cultivated grass and the ruderal Leucaena tree had prevented the establishment of other woody plants, thus inhibiting succession. However, the generous size of the plot (6 ha) together with its topographic heterogeneity allowed us to discern that the main factor explaining the composition of this patch was not biotic but rather abiotic and directly linked to hydric stress. It is important to note that the importance of topography in the spatial variation of vegetation was not detected the way we measured it in the field as the slope of the 5 × 5 m quadrat (see DS topography discussion below). It was only by the direct observation of the entire plot and its internal topographic variation that we were able to determine that the area still covered in grass had the highest elevation within this plot. Because it was 2 to 3 m above the level of the rest of the plot, the water table is deeper down than in the remaining area where conditions were more mesic, and the rich, dense SF was developed. We would like to make it clear that the entire plot was used during its cattle ranching stage as a single undivided grazing unit for 20 years, and all of it was abandoned at the same time. Dispersal failure does not explain the lack of woody species in this patch given that other areas of the same plot were similar distances from seed sources and were just as attractive or unattractive to dispersal agents as this patch was. The C4 African grass P. maximum and the deep-rooted, ruderal tree L. leucocephala are both well known for their drought tolerance (Centro Internacional de Agricultura Tropical, 1992), and this explains their dominance in the patch.

The sharp contrast in the spatial variability of succession within this large plot serves as a warning regarding drawing conclusions based on the data from the small plots or sampling areas that are the norm in tropical secondary succession studies. It would also be prudent to exercise caution when interpreting the results of studies based on chronosequences since they assume that the structural (density, biomass, richness) and compositional characteristics of SF are mainly explained by or correlated with the time since abandonment (or elapsed successional time). By averaging the values between different plots of similar age, these studies ignore the spatial heterogeneity and local factors that can exert a strong influence on secondary succession (see also Chazdon, 2014; Power et al., 2009). Other studies have also found a high degree of spatial heterogeneity in the attributes of the AR of SF (Almazán-Núñez, Arizmendi, Eguiarte, & Corcuera, 2012; Capers et al., 2005), but our study is the first to report such striking spatial heterogeneity in a tropical seasonally dry region and within a single large area (6 ha) that was abandoned all at once and that has similar distances to seed sources.

For the DS, previous research reported frequent visits by several frugivorous birds to the woody patches of this habitat, mainly from the adjacent-preserved forest, dispersing the seeds of zoochorous woody plants between the two habitats (Ortíz-Pulido, Laborde, & Guevara, 2000; Vázquez-Balbuena, 2014). However, our results strongly suggest that few of the dispersed seeds establish successfully or are recruited as adults or juveniles in the DS. In spite of the proximity of the adjacent preserved forest and the fact that both habitats had been protected since 1977, AR in the DS was very poor and sparse in woody species, even when compared with the 17-year-old SF. Our data indicate that there is strong niche limitation in the DS and that the limited seedling establishment of forest woody plants is not a result of dispersal failure (seed limitation sensu Levine & Murrel, 2003). The heterogeneous topography of the DS has important consequences in the spatial distribution of woody plants; however, this was not as important in our results (online Appendix B) as expected. This is due, in large part, to the way we measured it in the field by noting the slope of the sampled quadrat, thus we need a better estimator of variation in topography that reflects local conditions, and based on our observations and vegetation results, we suggest that elevation and the depth of the water table should be taken into account more quantitatively. The structural attributes of established adult vegetation in the dunes had more explanatory power than the slope of the sampled quadrat not only in the DS but also in the PF. This suggests that a better descriptor of canopy variation, such as for instance the degree of deciduousness of the trees on top of the site, should be used in future studies.

The largest trees in the DS were Bursera and Cedrela, but they were restricted to the deepest and therefore wettest dune slacks and they were considerably smaller than the largest conspecifics growing in the adjacent forest, highlighting the stressful hydric conditions of this habitat. Soils in this site are poorer and comprised almost of pure sand (Suárez et al., 2015), but the site is also more than 10 m higher in elevation than surrounding areas and therefore highly prone to desiccation. The DS plot is dominated by clumps of highly drought tolerant woody plants that, additionally, are able to grow in very poor sandy soils. These species include Diphysa, Randia, Tecoma, and Casearia, which are in turn surrounded by areas covered by native grasses (Andropogon spp., Trachypogon plumosus, and Schizachyrium scoparium, among others) growing on the dune crests and flanks, where even cacti (Opuntia stricta) are found (Moreno-Casasola & Travieso-Bello, 2006). Our results show that even though some woody plants are able to establish in the dune slacks, only a few of them are recruited as adults and most of them are shrubs or small trees (Figure 4). The absence or extremely low abundance in the DS of medium to tall tree species that are common in the old-growth forest of the region explains in great part the DS’s much lower values of diversity and similarity when compared to the other two sites. Therefore, the successional pathway to an old-growth semideciduous tropical forest, as suggested in the literature (Moreno Casasola & Travieso-Bello, 2006), is still quite far in the future for this system. It is possible that several centuries of soil formation in the site with the accumulation of leaf litter from drought resistant woody plants are needed before the current scrub patches can expand through a nucleation process and ultimately coalesce with other patches to form a continuous arboreal cover in the site.

Seedling Carpets of B. alicastrum and the Absence of a Keystone Animal

A noteworthy result of this study is the detection of dense seedling carpets of B. alicastrum beneath the crowns of conspecific adults in the preserved forest, since these carpets are usually found in more humid tropical forest (Álvarez-Loayza & Terborgh, 2011). Almost three decades ago, Delfosse (1990) working in the same forest patch noted that seedling carpets of this canopy tree species were only evident during mast seeding years and that they were ephemeral, since most if not all of them would die off during the dry season. Another peculiarity of this forest at the time was the very low density of tree seedlings and saplings in the understory relative to that of other tropical forests (Blain & Kellman, 1991). The organism mainly responsible for these phenomena was the red land crab, Gecarcinus lateralis Frem., which was very abundant in our study site prior to 2000 (with 50–100 thousand crabs/ha in the 1980s; Delfosse, 1990). The crabs were observed to feed heavily on the medium-sized to large seeds of woody plants (Capistrán-Barradas, Moreno-Casasola, & Defeo, 2006; Delfosse, 1990).

These crabs are now absent from the site and the seedling density we recorded beneath the crowns of large B. alicastrum adults was very high and similar to that recorded in 1987 by Delfosse (1990) inside crab exclosures constructed beneath the crowns of large B. alicastrum adults. The absence of this crab has changed the understory dramatically, as the forest floor is now covered year round with a continuous layer of decaying leaf litter and has higher densities of seedlings and saplings (see also Sherman, 2002). To the best of our knowledge, this is the first report of the absence of this crab from this study site and is also the first record of the defaunation (sensu Dirzo et al., 2014) of an invertebrate in a seasonally dry tropical forest (SDTF) having remarkable effects on its understory. That being said, the numerical dominance of B. alicastrum in the AR stage of the preserved forest does not mean that the future canopy of this site will be dominated by this tree; seedlings that form part of dense carpets are rarely or never recruited as juveniles (Álvarez-Loayza & Terborgh, 2011). Future studies, and specifically the resampling of our 25 m² permanent quadrats (with all individuals recorded in the AR), will tell us which woody species are thriving and which are not in this recently defaunated (no land crabs) forest.