Anthropogenic alterations of natural ecosystems creates fragmented landscapes with isolated remnants of natural vegetation separated by a matrix of contrasting environmental conditions (Saunders et al. 1991), where the edges of fragments are transition zones with biotic and abiotic conditions differing from those of the original ecosystem (Murcia 1995, Suzán et al. 1999, Broadbent et al. 2008). Environmental conditions in edges of fragments often affect several parameters of the life cycle in the native biota, which is a phenomenon commonly known as “border or edge effects” (Saunders et al. 1991, Murcia 1995, Honnay et al. 2005).
Changes in the environmental conditions in the edges of vegetation fragments are associated to factors such as the increase in the incidence of solar radiation, higher temperatures and dryness modifying the composition of the soil microbiota and the nutrients cycles, as well as increasing environmental stress conditions for the macroorganisms (Didham et al. 1996, Gehlhausen et al. 2000). Modifications in the edge environmental conditions will depend on the resulting surrounding matrix (i.e. grasslands, agricultural fields, urban settlements) (Jules & Shahani 2003). Additionally, colonization of vegetation fragments by exogenous species and the restrictive factors associated to the degree of fragmentation (i.e. distances among vegetation patches) creates adverse conditions for the persistence of mature community species (Turner 1996, Laurance et al. 1998, Benitez-Malvido 1998, Honnay et al. 1999, Zhu et al. 2004).
The reduction in the regeneration ability of species from mature communities is an important effect of habitat fragmentation (Benitez-Malvido 1998). For example, in a dominant palm species in Mexican tropical rain forests (Astrocaryum mexicanum), reduction in fragment sizes increased the competition among pollinators and non-pollinators floral visitors, because pollinator density is positively correlated to fragment size (Aguirre et al. 2011). Other examples such as Samanea saman, a common species in the tropical dry forests of Costa Rica, decreases in the germination percentages and seedling vigor, results from isolation among conspecific individuals increasing geitonogamic crosses and incompatible pollen deposition (Cascante et al. 2002). In Chilean and Argentinean forest, habitat fragmentation affects the survival of seedlings and decreases the fruit size of the common creeper Lapageria rosea (Henríquez 2004).
Alterations caused by the vegetation fragmentation on the life cycles of mature communities species, can be evaluated comparing the age or size structure of populations, in areas with different disturbance levels such as cores and edges of the vegetation gaps fragments, since the structure of ages and/or sizes reflects predominant patterns of growth, reproduction and survival (Silvertown 1987, Silvertown & Lovett-Doust 1993). Frequently, species from mature vegetation exhibit a decrease in the recruitment rates, and an increase in the mortality in the border of fragments (Benitez-Malvido 1998). In contrast, the recruitment of species from early successional stages, often increases in the fragment borders (Laurance et al. 1998).
The Tropical Deciduous Forest (TDF) is one of the must threatened ecosystems in the world, as consequence of the deforestation associated to extensive grazing, agriculture expansion and wood extraction (Janzen 1988, Maass et al. 2002). In Mexico, around 60 % of the tropical communities correspond to TDF (Trejo 1996, Trejo & Dirzo 2000), which is an ecosystem characterized by high biological diversity and an important number of endemic species (Ceballos & García 1995, Lott & Atkinson 2002, Rzedowski 2006). Nevertheless, the areas with TDF in the country face serious conservation problems, remaining only 26 % of its original distribution in México for 2002 (Challenger et al. 2009, Koleff et al. 2012), and only a 10 % of their surface remains in conservation zones (Sánchez-Azofeifa et al. 2005, Koleff et al. 2012). Therefore, the implementation of suitable practices of handling and conservation of areas with TDC becomes fundamental. Additionally, it is important to consider that several species of trees that characterize tropical communities are particularly vulnerable to the effects of the fragmentation (Cascante et al. 2002), because many species have low densities, are pollinated by animals, and they have out-crossing pollination and complex self-incompatibility systems (Bawa 1974, Zapata & Arroyo 1978, Hamrick & Murawsky 1990).
In this study we evaluated the population structure and the reproductive success of Bursera fagaroides (Kunth) Engl. and B. palmeri S. Watson in the borders and core areas of TDC fragments in Querétaro, México. Bursera fagaroides and B. palmeri are two important species in mature vegetation of TDF in central México. Both species are particularly vulnerable to fragmentation because are dioecious (Rzedowski & Guevara-Féfer, 1992; Rzedowski et al., 2005), and in consequence, its reproductive success may be related to changes of density by the fragmentation and the consequent alteration of environmental conditions in the edges of TDF fragments.
Material and methods
Study area. Study area encompasses nine fragments of TDF distributed in a polygon with coordinates: (20° 26’ N and 20° 47’ N, 100° 26’ W and 100° 37’ W; Figure 1), located in the municipalities of Querétaro, Qro., and Apaseo el Grande, Guanajuato, México. Climates in the study area are semidry semiwarm and semidry temperate, and dominant soils are vertisol and phaeozem (INEGI, 2010). Dominant vegetation type is TDF surrounded by xeric shrub-lands (Zamudio et al. 1992). Average altitudes of fragments varied between 1,847 and 2,002 m a.s.l.
Species studied. Bursera fagaroides is a dioecious small tree or shrub occasionally hermaphrodite that reach 10 m in height and a trunk wit 30 cm in diameter. It has an exfoliation external crust of yellowish color and produces an off-white latex (Rzedowski & Guevara-Féfer 1992). It’s distributed between southeast USA to Oaxaca, in Mexico), and is considered an important element of the TDF, secondary vegetation from TDF and xerophytic shrub-lands, with a recent severe abundant reduction in the TDF (Rzedowski et al. 2004). The phenophase starts in April and May, with fructification between November and December and foliage present from June to November (Rzedowski & Guevara-Féfer 1992).
Bursera palmeri is a dioecious small tree or shrub that reach 8 m in height, which produce an aromatic resin, with a trunk reaching up to 30 cm in diameter with a green to reddish-green exfoliating cortex (Rzedowski & Guevara-Féfer 1992). In Mexico it’s distributed in TDF from the following states: Aguascalientes, Durango, Guanajuato, Jalisco, Michoacán, Querétaro and Zacatecas. This species inhabit at altitudes higher than 1,400 m a.s.l., and its abundance shows a significant decrease due to extraction for fuel (Rzedowski et al. 2005). The flowering period range between May and June, and the species preserve the foliage between May and November (Rzedowski & Guevara-Féfer 1992).
Fragments selected for the study. In order to evaluate de edge effect in fragmented vegetation we selected nine TDF fragments using orthophotos scale 1:10,000 (INEGI 2008), and field survey between May and September 2012. The fragments were selected considering the presence of the two Bursera species and a qualitative conservation status that considered the presence of edges clearly distinguishable of the surrounded matrix, appreciable abundance of primary tree species with DBH ≥ 15 cm (e.g. Cedrela dugesii, Ceiba aesculifolia, Erythrina coralloides, Eysenhardtia polystachia and Ipomoea murucoides, Lysiloma microphylla, Senna polyantha), and an average canopy height of 8 m. Each fragment was mapped, and the plots for environmental contrasting conditions (edge and interior of the fragment) were randomly selected by using a network of cells (50 × 20 m) with the program ArcMAp v. 10.1, and defining their central geographic coordinates to locate Whittaker modified plots of 25 × 50 m (Stohlgren et al. 1995) in the field. When selected coordinates of two or more plots were too close, within the same condition (interior or edge) in a fragment, we maintained a minimum distance of 200 m between centers of plots. The total area and surrounding matrix for each fragment was calculated with the software ESRI® ArcMapTM 10.0. The surrounding matrix considered a buffer outside area of 500 m wide according to Charbonneau and Fahrig (2004) and Houlahan et al. (2006).
For each condition (edge and interior of the fragments) we measured respectively 16 and 41 Whitakker plots. The number of plots per fragment differs according to the fragment size (Table 1). The plots in the edges were located in a buffer zone (0-50 m) from the edge to inside of fragments, while plots in the interior were located close to the center of each fragment (Figure 1), excluding the possible border effects (Williams-Linera 1990a, Chen et al. 1992).
Fragment | Area (has) | Condition | PAI (m2 / m2) |
Texture | Apparent density (g/cm3) |
Soil Compaction |
---|---|---|---|---|---|---|
1 | 7.68 | Edge (1) | 2.2 ± 0.19 | SaCL | 0.52 ± 0.16 | Very low |
Interior (2) | 2.83 ± 0.12 | CL | 0.5 ± 0.16 | Very low | ||
2 | 7.72 | Edge (2) | 1.75 ± 0.14 | SCL | 0.61 ± 0.12 | Very low |
Interior (1) | 2 ± 0.14 | SCL | 0.6 ± 0.14 | Very low | ||
3 | 9.39 | Edge (1) | 1.75 ± 0.22 | SCL | 0.61 ± 0.14 | Very low |
Interior (1) | 1.73 ± 0.3 | SCL | 0.93 ± 0.2 | Very low | ||
4 | 11.43 | Edge (1) | 0.86 ± 0.44 | SCL, FCL | 0.55 ± 0.09 | Very low |
Interior (4) | 1.32 ± 0.15 | CL | 0.62 ± 0.16 | Very low | ||
5 | 26.57 | Edge (1) | ND | FCL | 0.69 ± 0.09 | Very low |
Interior (2) | ND | FCL | 0.72 ± 0.19 | Very low | ||
6 | 26.93 | Edge (1) | 2.36 ± 0.17 | FCL | 0.89 ± 0.1 | Very low |
Interior (3) | 2.82 ± 0.31 | FCL | 0.82 ± 0.04 | Very low | ||
7 | 44.19 | Edge (1) | 2.1 ± 0.13 | CL, SC | 0.16 ± 0.61 | Very low |
Interior (5) | 1.73 ± 0.26 | SCL | 0.43 ± 0.05 | Very low | ||
8 | 104.08 | Edge (2) | 1.5 ± 0.11 | SCL, SC, CoSC, FCL, CL | 0.79 ± 0.04 | Very low |
Interior (12) | 1.82 ± 0.14 | SaCL, SCL | 0.69 ± 0.09 | Very low | ||
9 | 249.59 | Edge (6) | 1.75 ± 0.15 | SC, SaCL, CL, SCL | 0.36 ± 0.33 | Very low |
Interior (11) | 1.89 ± 0.14 | SaCL, SCL, CL | 0.66 ± 0.06 | Very low |
Environmental conditions in the edges and interior of the fragments. Plant area index (PAI) and soil compaction were obtained for each Whitakker plot in the fragments. PAI was selected as an estimator of the canopy cover and solar radiation within the plots, and soil compaction as an indirect estimator of the effect of anthropogenic activities in the area.
Plant Area Index (PAI) of the fragments.- A couple of Licor LAI 2000 (PCA, LI-COR®, Lincoln, NE, USA) instruments were used to calculate the PAI of the plots according to Asner et al. (2003) and de Wasseige et al. (2003), including leafs and branches for each measurement. All the measurements were conducted between 5.00 and 7.00 h with clear skies, in the rain season (October-November) 2012. One instrument sensor installed at 1 m height, beneath the vegetation canopy with a 90° cap, and the second instrument in an open area measuring automatically every 30 seconds as an outside-canopy control. In the fragment, three PAI measurements were conducted with a systematic sampling across the central longitudinal axis of each Whittaker plot, and with 10 m in distance between measurements. Each PAI measurement consists on the average of three beneath canopy measurements, and one outside measurement from the control sensor. A Wilcoxon T test was used to compare plots in edges and interior of TDF fragments (Sokal & Rohlf 1995).
Soil compaction.- Soil compaction was estimated by means of its apparent density and texture (Siebe et al. 1996), and classified according to texture ranges suggested by the Mexican system of integrated differential diagnosis of soils (Soto et al. 2003, Córdoba-Athanasiadis 2010). A total of 342 samples, 246 from the edges and 96 from the fragment interiors plots, were obtained with metal cylinder cores (124.35 cm3). In each Whittaker plot six equidistant samples were collected. In order to measure apparent density, the soil samples were weighted in fresh and after being dried in an oven (105 °C, 24 h), and the cylinder and aluminum foil weights were subtracting from the total weight of each sample. Soils densities for each condition (edge and interior of fragments) were compared with a Wilcoxon T test (Sokal & Rohlf 1995).
Population structure. The population structure of Bursera fagaroides and B. palmeri in the edge and interior of fragments was estimated within Whittaker modified plots (Stohlgren et al. 1995), which include a nested design with four sub-plot levels (A, B, C and D). In the subplots A (0.5 × 2 m), we register frequencies of individuals with height ≤ 50 cm (juvenile type 1); for the subplots B (2 × 5 m) we register individuals with height > 50 cm, and diameter at breast height (DBH) between 1 and 5 cm (juvenile type 2); in the subplots C (5 × 20 m) individuals between 5 and 10 cm in DBH (juvenile 2), and in subplots D (20 × 50 m) individuals > 10 cm in DBH (adults) Seedling recruitment was not estimated because the majority of seedling die in few weeks.
Comparisons of densities of individuals between the edge and interior of the fragments were analyzed by a bi-factorial analysis of variance (ANOVA) with a randomized block design, with condition (edge and interior) and size categories (juveniles type 1 and 2 and adults) as factors and fragments as blocks. Density data were square root transformed (Sokal & Rohlf 1995). Within each fragment, population structures in edge and interior were compared by a G heterogeneity test (Sokal & Rohlf 1995).
Additionally, in order to explore the contribution of surrounding matrix in the population structure of Bursera species, we conducted a Spearman correlation analysis with values for fragments of the proportion of surrounding matrix corresponding to the area covered by principal components of matrix (i.e. roads, paddocks, rainfed agriculture, rural settlements and areas with secondary vegetation), and values of density and proportion of individuals of two Bursera species in the three categories (juvenile 1, juvenile 2 and adults).
Reproductive individual size and reproductive potential in the fragments. Sex, height and two perpendicular canopy diameters of reproductive trees in the “D” subplots within each plot were recorded. The volume of reproductive individuals was calculated as a unique estimator of individual size with the inverted truncated cone formula:
Where h is the cone height, R is the mean radius of the circumference corresponding to the canopy and r is the radius of the circumference corresponding to the DBH.
The volume of reproductive individuals was compared between the edge and interior of the fragments with a Student “t” test, with volume data previously log-transformed (Sokal & Rohlf 1995).
Reproductive potential (PR) in the edge and interior of each fragment was estimated as:
Where d and d represent densities for female and male reproductive individuals in 1,000 m2; dAd adult densities/1,000m2; and ns the number of sexes in the Plot (equal to two if both sexes are present in the plot). PR values ranges between cero (no reproductive individuals) and one (all adult individuals reproductive, and both sexes in the plot). PR values were compared with a Wilcoxon T test (Sokal & Rohlf 1995).
For each species we compared adult densities between edge and interior of fragments, and between sexes with a bifactorial ANOVA in a randomized block design (fragments as blocks). Density data were transformed by the square root function (Sokal & Rohlf 1995).
Seeds viability and germinative ability. Seeds were collected between April and May 2013 for Bursera fagaroides, and between January and February 2013 for Bursera palmeri. In each fragment, for each species and each condition, seeds from at least five adult trees with similar DBH (15 cm < DBH < 20 cm) were collected. Percentages of viable seeds for each species and condition were calculated by a flotation test from a random sample of 100 seeds. Additionally, in B. fagaroides viability was also estimated with a tetrazolium chloride (TC) test, from seeds selected after the flotation test. In each condition for each fragment 30 seeds were selected and then scarified and exposed to a 0.2 % solution of TC for 24 h (Henríquez 2004). TC evaluation was not conducted in B. palmeri because the insufficient amount of seeds. The proportion of viable seeds in each condition was compared with a paired t test with proportions transformed with the arc-sine function (Sokal & Rohlf 1995).
Mature seeds of B. fagaroides selected after a flotation test were used to compare germination percentages between the edge and interior of fragments (100 seeds each) in laboratory conditions. Seeds were disinfected with sodium hypochlorite (30 %, 5 min) and then washed in distilled water for 5 min (Godínez-Álvarez et al. 2008). Seeds were revised daily and germination was considered after root emergence. Seed germination percentages between conditions were compared with a pair t test after the transformation of data with the arc-sine function (Sokal & Rohlf 1995).
All the statistical analysis was conducted with JMP 8.0 software.
Results
Fragment description. The arboreal stratum in most of the fragments was dominated by species such as Bursera fagaroides, Ipomoea murucoides, Lysiloma microphyllum, Myrtillocactus geometrizans y Senna polyantha, other common species were, Bursera palmeri, Ceiba aesculifolia, Celtis caudata, Erythrina coralloides, Prosopis laevigata. Evidence of wood extraction for the arboreal species was very scarce or null in some fragments. The shrub stratum was highly disturbed by cattle and goat grazing, and the dominant species were Acacia farnesiana, Celtis pallida, Forestiera phillyreoides and Karwinskia humboldtiana. The more abundant surrounding matrices were rainfed agriculture (28.8 %) and paddocks (27.2 %), while roads and rural areas were the least commons.
The Plant Area index (PAI) values for the fragment edges range between 0.86 and 2.36 m2/m2 with a mean and standard error (SE) = 1.78 ± 0.17 m2/m2, while in the interior range between 1.32 and 2.83 m2/m2, with mean and SE = 2.02 ± 0.19 m2/m2 (Table 1), but no significant differences were detected for PAI values between the edge and interior of the fragments (W = 12.0, d.f. = 7, P = 0.10).
Apparent densities of soil according to the texture ranks suggested by the integrated differential diagnosis had small values, indicating non-compaction in the fragments studied (Table 1). Average and SE soil densities for all edge fragments were 0.64 ± 0.05 g/cm3, and for the interiors 0.71 ± 0.04 g/cm3. No significant differences between the edge and interior of fragments were detected (W = 8.5, d.f. = 7, P = 0.27).
Population structure. Population structure of Bursera fagaroides and B. palmeri showed higher densities of individuals for all size classes in the edge of fragments (F = 43.92, d.f. = 1, P < 0.0001 and F = 31.39, d.f. = 1, P < 0.0001, respectively). In B. fagaroides differences in densities among size classes were not significant (P = 0.92), and the interaction between condition (edge, interior) and size classes was also no significant (P > 0.53). In contrast, for B. palmeri densities of juvenile 1 size class was smaller compared to other size classes (F = 10.11, d.f. = 2, P = 0.0003), and the interaction between condition and size classes was significant (F = 5.12, d.f.= 2, P = 0.0105), due to the similarities on juvenile 1 for both conditions (Table 2). On the other hand, B. fagariodes had higher densities than B. palmeri (F = 18.83, d.f. = 1, P = 0.0002). However, even that total densities of B. fagaroides individuals were higher in the edges (F = 38.8, d.f. = 1, P < 0.0001), densities of this species in the interior were intermediate between densities of B. palmeri in the edges and interior of the fragments, leading to a significant interaction between the factors species and condition ( F = 7.63, d.f. = 1, P = 0.0108).
Bursera fagaroides | Bursera palmeri | |||
---|---|---|---|---|
Class | Edge | Interior | Edge | Interior |
Juvenile 1 | 33.1 (± 12.8) | 4.0 (± 1.3) | 1.9 (± 0.8) | 1.1 (± 0.7) |
Juvenile 2 | 28.2 (± 9.2) | 4.5 (± 1.5) | 8.7 (± 2.8) | 1.9 (± 0.7) |
Adults | 40.2 (± 12.7) | 2.3 (± 0.6) | 11.6 (± 3.2) | 1.9 (± 0.9) |
TOTAL | 101.5 (± 24.7) | 10.9 (± 2.6) | 22.1 (± 5.9) | 4.9 (± 1.9) |
REPRODUCTIVE | ||||
Adults | 2.8 (± 0.6) | 0.9 (± 0.3) | 0.8 (± 0.3) | 0.3 (± 0.3) |
Adults | 4.7 (± 1.1) | 1.3 (± 0.3) | 0.9 (± 0.3) | 1.0 (± 0.5) |
For both species, population structure was different between the edges and interior for most of the fragments (12.74 ≤ G ≤ 152.13, d.f. = 2, P < 0.05) (Figure 2). In the interior of the fragments, the higher proportion of individuals was concentrated in juvenile 1 and juvenile 2 classes, while in the edges the biggest proportion on individuals were concentrated on juvenile 2 and/or adult classes (Figure 2, Table 1). The comparisons of population structure between the edge and interior of different fragments in both species, suggest that juvenile trees grow slowly in the interior of fragments, leading to a concentration of individuals on juvenile 1 and juvenile 2 classes, while in the edge of fragments, juvenile trees seem grow faster, changing toward upper size classes in smaller times, and leading to a concentration in juvenile 2 and adult categories.
Regarding the relationship between characteristics of surrounding matrix of fragments and the population structure, no significant associations were detected between the area covered by predominant components of the surrounding matrix, and the density and proportion of individuals of both Bursera species (P > 0.05).
Reproductive individual size and reproductive potential in the fragments. The size of the reproductive trees (average volume ± SE) for both species tend to be higher in the interior than the edges of fragments, for Bursera fagaroides 35.0 ± 2.8 m3 vs 29.8 ± 1.9 m3 respectively, and for B. palmeri 36.0 ± 3.1 m3 vs 27.7 ± 2.6 m3 respectively. However, differences among conditions were non-significant: B. fagaroides: t = 1.67, d.f. = 203, P = 0.10; and B. palmeri: t = 1.45, d.f. = 104, P = 0.14). The reproductive potential of B. fagaroides was higher in the interior of the fragments than in the edges (W = 20.5, d.f. = 1, P = 0.012; Table 3), and for B. palmeri a similar but no significant trend was detected (W = 14.0, d.f. = 1, P = 0.10; Table 3).
Busera fagaroides | Bursera palmeri | |||||
---|---|---|---|---|---|---|
Fragment | Condition | PR | Vb / VbCT (%) | Gm (%) | PR | Vb (%) |
1 | Edge | 0.41 | 52 / 37 | 82 | 0.00 | 43 |
Interior | 0.89 | 28 / 78 | 76 | 1.00 | 25 | |
2 | Edge | 0.09 | 26 / 66 | 89 | 0.04 | 39 |
Interior | 1.00 | 18 / 53 | 87 | 0.33 | 18 | |
3 | Edge | 0.09 | 32 / 56 | 88 | 0.05 | 38 |
Interior | 0.50 | 56 / 72 | 85 | 0.57 | 41 | |
4 | Edge | 0.14 | 73 / ND | 77 | 0.50 | 74 |
Interior | 0.50 | 57 / ND | 70 | 0.00 | 43 | |
5 | Edge | 0.25 | 69 / ND | 98 | 0.25 | 65 |
Interior | 0.50 | 27 / ND | 93 | 0.00 | 33 | |
6 | Edge | 1.00 | 54 / 44 | 89 | 0.21 | 62 |
Interior | 0.88 | 59 / 78 | 80 | 0.00 | 63 | |
7 | Edge | 0.12 | 35 / 36 | 79 | 0.00 | 42 |
Interior | 0.57 | 20 / 51 | 75 | 0.50 | 44 | |
8 | Edge | 0.14 | 67 / 49 | 87 | 0.19 | 60 |
Interior | 0.17 | 55 / 64 | 84 | 0.89 | 57 | |
9 | Edge | 0.52 | 43 / 38 | 93 | 0.56 | 60 |
Interior | 1.00 | 60 / 49 | 93 | 0.90 | 57 | |
Mean (EE) | Edge | 0.31 | 50.11 / 46.54 | 86.89 | 0.20 | 53.67 |
(0.10) | (5.74) / (4.23) | (2.22) | (0.07) | (4.41) | ||
Interior | 0.67 | 42.22 / 63.69 | 82.56 | 0.47 | 42.33 | |
(0.10) | (6.11) / (4.82) | (2.66) | (0.14) | (5.06) |
Density of reproductive individuals were higher in the edge than in the interior of fragments for Bursera fagaroides (F = 20.15, d.f. = 1, P = 0.0002), but not for B. palmeri (F = 1.29, d.f. = 1, P = 0.27). In two species density of male individuals tends to be higher than females, both in edge and interior of fragments (Table 2), but differences between sexes were not significant (B. fagaroides: F = 3.16, d.f. = 1, P = 0.09; and B. palmeri: F = 3.07, d.f. = 1, P = 0.09), like the interaction between conditions and sex factors (B. fagaroides: F = 0.43, d.f. = 1, P = 0.52; and B. palmeri: F = 1.20, d.f. = 1, P = 0.29).
Seeds viability and germinative ability. Seed viability estimated by flotation in Bursera fagaroides not showed differences among the edge and interior of fragments (t = 1.22, d.f. = 8, P = 0.26; Table 3), while in B. palmeri viability was higher in the edge (t = 2.40, d.f. = 8, P = 0.04; Table 3). In contrast, viability estimated with tetrazolium chloride for B. fagaroides was significantly higher for seeds collected in the interior of fragments (t = 2.58, d.f. = 6, P = 0.04; Table 3). Finally, germination percentages in B. palmeri were consistently higher in the edge than interior of fragments (Table 3), leading to significant differences between conditions (t = 4.49, d.f. = 8, P = 0.002).
Discussion
In this study we expected that the population structure and the reproductive success of Bursera fagaroides and B. palmeri, could reflect changes associated to the modification of the environmental conditions in the edges of fragments of the Tropical Deciduous forest (TDF), since both species are dioecious, have low densities and they are used frequently for fuel (Rzedowski et al. 2005), leading to variations in density and size structure between the edge and the interior of fragments, as a result of smaller rates of survival, growth and reproduction in the edges.
In contrast to our expectations, Bursera. fagaroides and B. palmeri seems to be favored by the environmental conditions that predominate in the edges of fragments, where the density of both species was greater. In this sense, probably the variation between edge and interior of fragments in the solar radiation that penetrates through canopy has influenced in the density. Similar results reported by Ortiz-Pulido and Rico-Gray (2006) from “El Morro de la Mancha” (Veracruz) in B. fagaroides shows that the presence of the species, their germination and growth, were favored by the TDF moderate shaded environment, contrasting with open exposed sites, or in heavy shaded environments within the Tropical semideciduous forests. However, our estimations of the PAI in the edge and inside fragments had similar values, suggesting two possibilities: 1) that the solar incidence of radiation is not an important factor in the variation of the density of individuals, or 2) that the differences in the incidence of the solar radiation between both conditions are currently minimum because the rapid recovery of the canopy. With respect to our estimation of the PAI, it is important to mention that the real value of LAI can be overestimated due to the contribution of trunks and branches of the canopy. Nevertheless, this bias for our instruments could be considered as non-significant (Cutini et al. 1998). On the other hand, is also important to consider that use of PAI to estimate the differences in the solar radiation that penetrates the canopy between edge and the interior of fragments, could be a coarse method to detect a subtle variation in the quantity and quality of light affecting growth of plants under the canopy, compared with other methods such as measuring of the photosynthetically active radiation (PAR). However, measures of LAI and PAR generally show a high correlation (e.g., Buckley et al. 1999, Bellow & Nair 2003, Hyer & Goetz 2004), so we considered that PAI estimations derived from LAI measures have the sufficient resolution to detect differences in radiation passing through the canopy.
In addition to the potential effect of differences in the solar radiation between edge and interior of fragments on species densities, other factors must be considered. For example, a lower density of individuals in the interior of fragments may be the result of higher rates of seed predation than in edges, whether seed predators concentrate their activity away of edges to avoid their own predation (e.g., Guzmán-Guzman & Williams-Linera 2006, Penido et al. 2015). On the other hand, higher density of plants in the edge, in contrast with the interior of fragments, may be the result of higher rates of removal and deposition of fruits that are more conspicuous for frugivore birds in the edges (Galetti et al. 2003). However, this seems unlikely for Bursera fagaroides, because germination rate of its seeds are reduced when pass through the digestive tract of birds (Ortiz-Pulido & Rico-Gray 2006).
Soil compaction is a non- relevant factor for understanding differences in densities and population structure in the edges and interior of the fragments. Soil compaction values were similar for both conditions and exhibited very low values according to the integrated differential diagnosis system for mexican soils (Soto et al. 2003, Córdoba-Athanasiadis 2010). Even that paddocks were a common surrounding matrix in the regions, and cattle was frequently observed for small periods of time, this activity seems not to affect significantly the soil compaction. It is important to consider that we used an indirect estimator of soil compaction, which probably did not detected differences between edge and interior of fragments. Nevertheless, the measures of apparent density of soil are frequently used as reliable estimators of compaction (Carter 1990, Håkansson & Lipiec 2000, Blanco-Sepúlveda 2009). Even more, compaction generally has detrimental effects on the growth of plants (Kozlowski 1999), which apparently did not happening in the studied fragments as is suggested by high densities of both Bursera species in the edges.
Previously, we mentioned that the similarity in values of the PAI between the interior and the edge of fragments could be the result of a relatively fast recovery of the canopy under the initial conditions of greater incidence of radiation in the edges. A recovery of the relatively accelerated canopy can be possible if the incidence of greater solar radiation leads to an increase of the recruitment rates, survival and growth of some common arboreal species in the mature and/or secondary communities of the TDF, such as Bursera fagaroides and B. palmeri (Rzedowski & Guevara-Féfer 1992). Generally, the arboreal species that comprise the mature vegetal communities are affected by modifications from environmental conditions after ecosystem fragmentation, showing higher mortality rates and smaller recruitment in the edges (e.g., Williams-Linera 1990b, Benitez-Malvido 1998, Nascimento & Laurance 2004). Nevertheless, has been documented that recruitment and growth of some species from mature communities increased in conditions of moderate disturbance such as in the edges fragments, becoming, dominant elements of the secondary communities (Chen et al. 1992, Oliveira-Filho et al. 1997, McDonald & Urban 2004, Laurance et al. 2006).
Differences in the population structure of Bursera fagaroides and B. palmeri between the edge and the interior of TDF fragments, is congruent with the idea that environmental conditions in the edges of fragments are or were favorable for increasing growth rates. Nevertheless, since the size or age structure of a population depends on the patterns of recruitment and survival (Silvertown 1987, Caswell 2001), a higher growth rate in the edges is not sufficient to explain the differences observed in the population structure of both species between edges and interior of the fragments. In this sense, one more reasonable explanation is than a higher individual growth rate in the edges has been accompanied with high recruitment and survival rates, as is suggested by the considerably high density of individuals of all size categories in the edges of fragments. Nevertheless, it is necessary to demonstrate by long term demographic studies that the recruitment, survival and growth rates of B. fagaroides and B. palmeri increased under conditions of greater light intensity. Additionally, it is necessary to evaluate the importance of other factors that contribute in the differences of density and variations of the population structure like the seeds and seedlings predation (Benitez-Malvido 1998, Cascante et al. 2002, Guzmán-Guzman & Williams-Linera 2006, Penido et al. 2015), changes in the density of common species in the mature community, which are less tolerant to the light, high temperatures and drying (Benitez-Malvido 1998, Laurance et al. 2006), and the modification of the abundance of generalist and specialist dispersers (Restrepo et al. 1999, Albrecht et al. 2013).
The increase in the recruitment rate of Bursera fagaroides and B. palmeri in the edge of fragments, indicates differences in juvenile densities between the edge and the interior of the fragments, and does not seem to be a consequence of the increasing densities of reproductive adult individuals, since in the edge of fragments the density of reproductive individuals was slightly higher than in the interior, which is reflected in the reproductive potential, showing smaller values for the edges compared to the interior of the fragments. Also, the similarity in the proportion of masculine and feminine reproductive individuals in the edge and inside fragments, as well as the similarity in the percentage of viability and germination, does not provide indications that these factors have important effects in the difference of density for juvenile individuals between both conditions. Even when germination percentages were significantly greater in the edges than in the interiors, the values obtained were high in both conditions (87 % and 83 %, respectively). A possible explanation for this small but significant difference could be a greater temperature fluctuation in the edge condition that increases the seeds germination, which has been demonstrated for several species of Bursera under experimental conditions (Bonfil-Sánders et al. 2008).
Finally, it is important to consider that our evaluations for densities of reproductive individuals, sex rates and the viability and germination for both species, represent only one snapshot in the process that have determined the variation in the density and the population structure between the edge and interior of fragments. In this way, it is necessary to consider the temporary variation to improve our estimation of the edge effects on the reproduction of studied species. For example, although we could not currently detect a bias between sexes, which, by the way, may depend on several factors as differences between male and female plants in the age or size to the sexual maturity, or differences in the flowering frequency between plants of different sexes (Meagher 1980), to long term, bias in sexual rates between edge and interior of fragments may be expected, because in contrast with male plants, survival rates of female plants can be decreased in stressed environments as consequence of a less resources allocation to the survival and a high allocation to the reproduction (e.g., Agren 1988, Gehring & Linhart 1993, Yu & Lu 2011). Additionally, it is important to consider other factors that can influence juvenile densities in the edge and the interior of the fragments. For example, the high production of seeds in the edge of fragments as consequence of a greater environmental stress on the female individuals, which has been observed in Bursera fagaroides (Ortiz-Pulido & Pavón 2010); and the ability of seeds to remain viable in seeds banks, and to germinate massively under conditions of greater luminosity. With respect to this second point, the time that seeds can remain viable is an aspect still little explored in species of the genus Bursera. Nevertheless, the available information suggests that the period in which the seeds remain viable in natural conditions, can be relatively short (Bonfil-Sánders et al. 2008, Sánchez-Martínez et al. 2011, Morgan & Jose 2013).