Tree size and fecundity influence ballistic seed dispersal of two dominant mast-fruiting species in a tropical rain forest

Seed production, seed dispersal, and seedling recruitment are integral to forest dynamics, especially in masting species. Often these are studied separately, yet scarcely ever for species with ballistic dispersal even though this mode of dispersal is common in legume trees of tropical African rain forests. Here, we studied two dominant main-canopy tree species, Microberlinia bisulcata and Tetraberlinia bifoliolata (Caesalpinioideae), in 25 ha of primary rain forest at Korup, Cameroon, during two successive masting events (2007/2010). In the vicinity of c. 100 and 130 trees of each species, 476/580 traps caught dispersed seeds and beneath their crowns c. 57,000 pod valves per species were inspected to estimate tree-level fecundity. Seed production of trees increased non-linearly and asymptotically with increasing stem diameters. It was unequal within the two species’ populations, and differed strongly between years to foster both spatial and temporal patchiness in seed rain. The M. bisulcata trees could begin seeding at 42–44 cm diameter: at a much larger size than could T. bifoliolata (25 cm). Nevertheless, per capita life-time reproductive capacity was c. five times greater in M. bisulcata than T. bifoliolata owing to former’s larger adult stature, lower mortality rate (despite a shorter life-time) and smaller seed mass. The two species displayed strong differences in their dispersal capabilities. Inverse modelling (IM) revealed that dispersal of M. bisulcata was best described by a lognormal kernel. Most seeds landed at 10–15 m from stems, with 1% of them going beyond 80 m (<100 m). The direct estimates of fecundity significantly improved the models fitted. The lognormal also described well the seedling recruitment distribution of this species in 121 ground plots. By contrast, the lower intensity of masting and more limited dispersal of the heavier-seeded T. bifoliolata prevented reliable IM. For this species, seed density as function of distance to traps suggested a maximum dispersal distance of 40–50 m, and a correspondingly more aggregated seedling recruitment pattern ensued than for M. bisulcata. From this integrated field study, we conclude that the reproductive traits of M. bisulcata give it a considerable advantage over T. bifoliolata by better dispersing more seeds per capita to reach more suitable establishment sites, and combined with other key traits they explain its local dominance in the forest. Understanding the linkages between size at onset of maturity, individual fecundity, and dispersal capability can better inform the life-history strategies, and hence management, of co-occurring tree species in tropical forests.

Improved methods for measuring forest landscape structure: LiDAR complements field-based habitat assessment

Conservation and monitoring of forest biodiversity requires reliable information about forest structure and composition at multiple spatial scales. However, detailed data about forest habitat characteristics across large areas are often incomplete due to difficulties associated with field sampling methods. To overcome this limitation we employed a nationally available light detection and ranging (LiDAR) remote sensing dataset to develop variables describing forest landscape structure across a large environmental gradient in Switzerland. Using a model species indicative of structurally rich mountain forests (hazel grouse Bonasa bonasia), we tested the potential of such variables to predict species occurrence and evaluated the additional benefit of LiDAR data when used in combination with traditional, sample plot-based field variables. We calibrated boosted regression trees (BRT) models for both variable sets separately and in combination, and compared the models’ accuracies. While both field-based and LiDAR models performed well, combining the two data sources improved the accuracy of the species’ habitat model. The variables retained from the two datasets held different types of information: field variables mostly quantified food resources and cover in the field and shrub layer, LiDAR variables characterized heterogeneity of vegetation structure which correlated with field variables describing the understory and ground vegetation. When combined with data on forest vegetation composition from field surveys, LiDAR provides valuable complementary information for encompassing species niches more comprehensively. Thus, LiDAR bridges the gap between precise, locally restricted field-data and coarse digital land cover information by reliably identifying habitat structure and quality across large areas.

Effects of elevation and land use on the biomass of trees, shrubs and herbs at Mount Kilimanjaro

The protection and sustainable management of forest carbon stocks, particularly in the tropics, is a key factor in the mitigation of global change effects. However, our knowledge of how land use and elevation affect carbon stocks in tropical ecosystems is very limited. We compared aboveground biomass of trees, shrubs and herbs for eleven natural and human-influenced habitat types occurring over a wide elevation gradient (866–4550 m) at the world's highest solitary mountain, Mount Kilimanjaro. Thanks to the enormous elevation gradient, we covered important natural habitat types, e.g., savanna woodlands, montane rainforest and afro-alpine vegetation, as well as important land-use types such as maize fields, grasslands, traditional home gardens, coffee plantations and selectively logged forest. To assess tree and shrub biomass with pantropical allometric equations, we measured tree height, diameter at breast height and wood density and to assess herbaceous biomass, we sampled destructively. Among natural habitats, tree biomass was highest at intermediate elevation in the montane zone (340 Mg ha−1), shrub biomass declined linearly from 7 Mg ha−1 at 900 m to zero above 4000 m, and, inverse to tree biomass, herbaceous biomass was lower at mid-elevations (1 Mg ha−1) than in savannas (900 m, 3 Mg ha−1) or alpine vegetation (above 4000 m, 6 Mg ha−1). While the various land-use types dramatically decreased woody biomass at all elevations, though to various degrees, herbaceous biomass was typically increased. Our study highlights tropical montane forest biomass as important aboveground carbon stock and quantifies the extent of the strong aboveground biomass reductions by the major land-use types, common to East Africa. Further, it shows that elevation and land use differently affect different vegetation strata, and thus the matrix for other organisms.

Using plant functional traits and phylogenies to understand patterns of plant community assembly in a seasonal tropical forest in Lao PDR

Plant functional traits reflect different evolutionary responses to environmental variation, and among extant species determine the outcomes of interactions between plants and their environment, including other plant species. Thus, combining phylogenetic and trait-based information can be a powerful approach for understanding community assembly processes across a range of spatial scales. We used this approach to investigate tree community composition at Phou Khao Khouay National Park (18°14’-18°32’N; 102°38’- 102°59’E), Laos, where several distinct forest types occur in close proximity. The aim of our study was to examine patterns of plant community assembly across the strong environmental gradients evident at our site. We hypothesized that differences in tree community composition were being driven by an underlying gradient in soil conditions. Thus, we predicted that environmental filtering would predominate at the site and that the filtering would be strongest on sandier soil with low pH, as these are the conditions least favorable to plant growth. We surveyed eleven 0.25 ha (50x50 m) plots for all trees above 10 cm dbh (1221 individual trees, including 47 families, 70 genera and 123 species) and sampled soils in each plot. For each species in the community, we measured 11 commonly studied plant functional traits covering both the leaf and wood economic spectrum traits and we reconstructed a phylogenetic tree for 115 of the species in the community using rbcL and matK sequences downloaded from Genebank (other species were not available). Finally we compared the distribution of trait values and species at two scales (among plots and 10x10m subplots) to examine trait and phylogenetic community structures. Although there was strong evidence that an underlying soil gradient was determining patterns of species composition at the site, our results did not support the hypothesis that the environmental filtering dominated community assembly processes. For the measured plant functional traits there was no consistent pattern of trait dispersion across the site, either when traits were considered individually or when combined in a multivariate analysis. However, there was a significant correlation between the degree of phylogenetic dispersion and the first principle component axis (PCA1) for the soil parameters.Moreover, the more phylogenetically clustered plots were on sandier soils with lower pH. Hence, we suggest that the community assembly processes across our sitemay reflect the influence ofmore conserved traits that we did not measure. Nevertheless, our results are equivocal and other interpretations are possible. Our study illustrates some difficulties in combining trait and phylogenetic approaches that may result from the complexities of integrating spatial and evolutionary processes that vary at different scales.

Effect sizes and standardization in neighbourhood models of forest stands: potential biases and misinterpretations

Effects of conspecific neighbours on survival and growth of trees have been found to be related to species abundance. Both positive and negative relationships may explain observed abundance patterns. Surprisingly, it is rarely tested whether such relationships could be biased or even spurious due to transforming neighbourhood variables or influences of spatial aggregation, distance decay of neighbour effects and standardization of effect sizes. To investigate potential biases, communities of 20 identical species were simulated with log-series abundances but without species-specific interactions. No relationship of conspecific neighbour effects on survival or growth with species abundance was expected. Survival and growth of individuals was simulated in random and aggregated spatial patterns using no, linear, or squared distance decay of neighbour effects. Regression coefficients of statistical neighbourhood models were unbiased and unrelated to species abundance. However, variation in the number of conspecific neighbours was positively or negatively related to species abundance depending on transformations of neighbourhood variables, spatial pattern and distance decay. Consequently, effect sizes and standardized regression coefficients, often used in model fitting across large numbers of species, were also positively or negatively related to species abundance depending on transformation of neighbourhood variables, spatial pattern and distance decay. Tests using randomized tree positions and identities provide the best benchmarks by which to critically evaluate relationships of effect sizes or standardized regression coefficients with tree species abundance. This will better guard against potential misinterpretations.

Canopy gaps promote selective stem-cutting by small mammals of two dominant tree species in an African lowland forest: the importance of seedling chemistry

Small mammals can impede tree regeneration by injuring seedlings and saplings in several ways. One fatal way is by severing their stems, but apparently this type of predation is not well-studied in tropical rain forest. Here, we report on the incidence of 'stem-cutting' to new, wild seedlings of two locally dominant, canopy tree species monitored in 40 paired forest understorey and gap-habitat areas in Korup, Cameroon following a 2007 masting event. In gap areas, which are required for the upward growth and sapling recruitment of both species, 137 seedlings of the long-lived, light-demanding, fast-growing large tropical tree (Microberlinia bisulcata) were highly susceptible to stem-cutting (83% of deaths) - it killed 39% of all seedlings over a c. 2-y period. In stark contrast, seedlings of the more shade-tolerant, slower-growing tree species (Tetraberlinia bifoliolata) were hardly attacked (4.3%). In the understorey, however, stem-cutting was virtually absent. Across the gap areas, the incidence of stem-cutting of M. bisulcata seedlings showed significant spatial variation that could not be explained significantly by either canopy openness or Janzen-Connell type effects (proximity and basal area of conspecific adult trees). To examine physical and chemical traits that might explain the species difference to being cut, bark and wood tissues were collected from a separate sample of seedlings in gaps (i.e. not monitored for stem-cutting). These analyses suggested that, compared with T. bifoliolata, the lower stem density, higher Mg and K and fatty acid concentrations in bark, and fewer phenolic and terpene compounds in M. bisulcata seedlings made them more palatable and attractive to small-mammal predators, likely rodents. We conclude that selective stem-cutting is a potent countervailing force to the current local canopy dominance of the grove-forming M. bisulcata by limiting the recruitment and abundance of its saplings. Given the ubiquity of gaps and ground-dwelling rodents in pantropical forests, it would be surprising if this form of lethal browsing was restricted to Korup.

Seven years of annual pollen influx at the forest limit in the Swiss Alps studied by pollen traps: relations to vegetation and climate

Annual pollen influx has been monitored in short transects across the altitudinal tree limit in four areas of the Swiss Alps with the use of modified Tauber traps placed at the ground surface. The study areas are Grindelwald (8 traps), Aletsch (8 traps), Simplon (5 traps), and Zermatt (5 traps). The vegetation around the traps is described. The results obtained are: (1) Peak years of pollen influx (one or two in seven years) follow years of high average air temperatures during June–November of the previous year for Larix and Picea, and less clearly for Pinus non-cembra, but not at all for Pinus cembra and Alnus viridis. (2) At the upper forest limit, the regional pollen influx of trees (trees absent within 100 m of the pollen trap) relates well to the average basal area of the same taxon within 10–15 km of the study areas for Pinus cembra, Larix, and Betula, but not for Picea, Pinus non-cembra, and Alnus viridis. (3) The example of Zermatt shows that pollen influx characterises the upper forest limit, if the latter is more or less intact. (4) Presence/absence of Picea, Pinus cembra, Larix, Pinus non-cembra, and Alnus viridis trees within 50–100 m of the traps is apparent in the pollen influx in peak years of pollen influx but not in other years, suggesting that forest-limit trees produce significant amounts of pollen only in some years. (5) Pollen influx averaged over the study period correlates well with the abundance of plants around the pollen traps for conifer trees (but not deciduous trees), Calluna, Gramineae, and Cyperaceae, and less clearly so Compositae Subfam. Cichorioideae and Potentilla-type. (6) Influx of extra-regional pollen derived from south of the Alps is highest in Simplon, which is open to southerly winds, slightly lower in Aletsch lying just north of Simplon, and lowest in Zermatt sheltered from the south by high mountains and Grindelwald lying north of the central Alps.

Stand scale variability of topsoil organic matter composition in a high-elevation Norway spruce forest ecosystem

Our knowledge about the effect of single-tree influence areas on the physicochemical properties of the underlying mineral soil in forest ecosystems is still limited. This restricts our ability to adequately estimate future changes in soil functioning due to forest management practices. We studied the stand scale spatial variation of different soil organic matter species investigated by 13C NMR spectroscopy, lignin phenol and neutral sugar analysis under an unmanaged mountainous high-elevation Norway spruce (Picea abies L.) forest in central Europe. Multivariate geostatistical approaches were applied to relate the spatial patterns of the different soil organic matter species to topographic parameters, bulk density, oxalate- and dithionite-extractable iron, pH, and the impact of tree distribution. Soil samples were taken from the mineral top soil. Generally, the stand scale distribution patterns of different soil organic matter compounds could be divided into two groups: Those compounds, which were significantly spatially correlated with topography/altitude and those with small scale spatial pattern (range ≤ 10 m) that was closely related to tree distribution. The concentration of plant-derived soil organic matter components, such as lignin, at a given sampling point was significantly spatially related to the distance of the nearest tree (p ≤ 0.05). In contrast, the spatial distribution of mainly microbial-derived compounds (e.g. galactose and mannose) could be attributed to the dominating impact of small-scale topography and the contribution of poorly crystalline iron oxides that were significantly larger in the central depression of the study site compared to crest and slope positions. Our results demonstrate that topographic parameters dominate the distribution of overall topsoil organic carbon (OC) stocks at temperate high-elevation forest ecosystems, particularly in sloped terrain. However, trees superimpose topography-controlled OC biogeochemistry beneath their crown by releasing litter and changing soil conditions in comparison to open areas. This may lead to distinct zones with different mechanisms of soil organic matter degradation and also stabilization in forest stands.

Primary lowland dipterocarp forest at Danum Valley, Sabah, Malaysia. Species composition and patterns in the understorey

The spatial distributions of species of tree ≥10 cm gbh were examined in two 4 ha plots and related to the local variation in topography and soil chemistry. The plots were similar in their species composition, particularly in terms of the densities of small trees, and they showed very similar edaphic characteristics. Size class distributions varied little within and between plots. Ordination of 0.25 ha subplots highlighted parallel gradients in the vegetation of both plots when the densities of trees ≥10 cm gbh were considered. Focusing on understorey trees in the 10-<50 cm gbh class at the 0.04 ha subplot scale showed a similar vegetation gradient in both plots closely associated with change from lower slope to ridge. No relationship with soil chemistry was found. On the ridges a special group of understorey species formed clumps and these species contributed importantly to the ordinations. Borneo has a regional history of occasionally severe droughts. It is suggested here that the observed patterns in the understorey are due to differential responses to low soil water supply, the ridges probably tending to dryness more than the lower slopes. Within the large and diverse family Euphorbiaceae, which dominates the understorey at Danum, there may be ecophysiological groupings of species. The long-term effects of disturbance interacting with local edaphic factors on forest structure and composition are discussed.

Primary forest dynamics in lowland dipterocarp forest at Danum Valley, Sabah, Malaysia, and the role of the understorey

Changes in species composition in two 4–ha plots of lowland dipterocarp rainforest at Danum, Sabah, were measured over ten years (1986 to 1996) for trees greater than or equal to 10 cm girth at breast height (gbh). Each included a lower–slope to ridge gradient. The period lay between two drought events of moderate intensity but the forest showed no large lasting responses, suggesting that its species were well adapted to this regime. Mortality and recruitment rates were not unusual in global or regional comparisons. The forest continued to aggrade from its relatively (for Sabah) low basal area in 1986 and, together with the very open upper canopy structure and an abundance of lianas, this suggests a forest in a late stage of recovery from a major disturbance, yet one continually affected by smaller recent setbacks. Mortality and recruitment rates were not related to population size in 1986, but across subplots recruitment was positively correlated with the density and basal area of small trees (10 to <50 cm gbh) forming the dense understorey. Neither rate was related to topography. While species with larger mean gbh had greater relative growth rates (rgr) than smaller ones, subplot mean recruitment rates were correlated with rgr among small trees. Separating understorey species (typically the Euphorbiaceae) from the overstorey (Dipterocarpaceae) showed marked differences in change in mortality with increasing gbh: in the former it increased, in the latter it decreased. Forest processes are centred on this understorey quasi–stratum. The two replicate plots showed a high correspondence in the mortality, recruitment, population changes and growth rates of small trees for the 49 most abundant species in common to both. Overstorey species had higher rgrs than understorey ones, but both showed considerable ranges in mortality and recruitment rates. The supposed trade–off in traits, viz slower rgr, shade tolerance and lower population turnover in the understorey group versus faster potential growth rate, high light responsiveness and high turnover in the overstorey group, was only partly met, as some understorey species were also very dynamic. The forest at Danum, under such a disturbance–recovery regime, can be viewed as having a dynamic equilibrium in functional and structural terms. A second trade–off in shade–tolerance versus drought–tolerance is suggested for among the understorey species. A two–storey (or vertical component) model is proposed where the understorey–overstorey species’ ratio of small stems (currently 2:1) is maintained by a major feedback process. The understorey appears to be an important part of this forest, giving resilience against drought and protecting the overstorey saplings in the long term. This view could be valuable for understanding forest responses to climate change where drought frequency in Borneo is predicted to intensify in the coming decades.

What drives European beech (Fagus sylvatica L.) mortality after forest fires of varying severity?

Predicting the timing and amount of tree mortality after a forest fire is of paramount importance for post-fire management decisions, such as salvage logging or reforestation. Such knowledge is particularly needed in mountainous regions where forest stands often serve as protection against natural hazards (e.g., snow avalanches, rockfalls, landslides). In this paper, we focus on the drivers and timing of mortality in fire-injured beech trees (Fagus sylvatica L.) in mountain regions. We studied beech forests in the southwestern European Alps, which burned between 1970 and 2012. The results show that beech trees, which lack fire-resistance traits, experience increased mortality within the first two decades post-fire with a timing and amount strongly related to the burn severity. Beech mortality is fast and ubiquitous in high severity sites, whereas small- (DBH <12 cm) and intermediate-diameter (DBH 12–36 cm) trees face a higher risk to die in moderate-severity sites. Large-diameter trees mostly survive, representing a crucial ecological legacy for beech regeneration. Mortality remains low and at a level similar to unburnt beech forests for low burn severity sites. Beech trees diameter, the presence of fungal infestation and elevation are the most significant drivers of mortality. The risk of beech to die increases toward higher elevation and is higher for small-diameter than for large-diameter trees. In case of secondary fungi infestation beech faces generally a higher risk to die. Interestingly, fungi that initiate post-fire tree mortality differ from fungi occurring after mechanical injury. From a management point of view, the insights about the controls of post-fire mortality provided by this study should help in planning post-fire silvicultural measures in montane beech forests.

Epiphytic orchid diversity and community composition in managed and natural moist evergreen Afromontane forest in SW Ethiopia

In SW Ethiopia, the moist evergreen Afromontane forest has become extremely fragmented and most of the remnants are intensively managed for coffee cultivation (Coffea arabica), with considerable impacts on biodiversity and ecosystem functioning. Because epiphytic orchids are potential indicators for forest quality and a proxy for overall forest biodiversity, we assessed the effect of forest management and forest fragmentation on epiphytic orchid diversity. We selected managed forest sites from both large and small forest remnants and compared their epiphytic orchid diversity with the diversity of natural unfragmented forest. We surveyed 339 canopy trees using rope climbing techniques. Orchid richness decreased and community composition changed, from the natural unfragmented forest, over the large managed forest fragments to the small managed forest fragments. This indicates that both forest management and fragmentation contribute to the loss of epiphytic orchids. Both the removal of large canopy trees typical for coffee management, and the occurrence of edge effects accompanying forest fragmentation are likely responsible for species loss and community composition changes. Even though some endangered orchid species persist even in the smallest fragments, large managed forest fragments are better options for the conservation of epiphytic orchids than small managed forests. Our results ultimately show that even though shade coffee cultivation is considered as a close-to-nature practice and is promoted as biodiversity conservation friendly, it cannot compete with the epiphytic orchid conservation benefit generated by unmanaged moist evergreen Afromontane forests.

Epiphytic orchid diversity found on trees and shrubs in coffee cultivated Afromontane rainforest, SW Ethiopia

The moist evergreen Afromontane forest of SW Ethiopia has become extremely fragmented and most remnants are intensively managed for cultivation of coffee (Coffea arabica). We investigated the distributions of epiphytic orchids in shade trees and their understory in forests with contrasting management intensity to determine biodiversity losses associated with coffee cultivation and to determine the capacity of coffee shrubs to act as refugia for orchid species. We studied epiphytic orchids in managed forests and natural forests and recorded orchid diversity and abundance in different tree zones of 339 trees and in the understory. Coffee management was associated with a downward shift of orchid species as orchid species were occurring in significantly lower tree zones in managed forest. The number of shrubs in the understory of managed forest was not higher than in natural forests, yet orchid abundance was higher in the understory of managed forests. Local extinctions of epiphytic orchids and species losses in the outer tree zones (a contraction of habitat) in managed forests are most likely driven by losses of large, complex-structured climax trees, and changes in microclimate, respectively. Coffee shrubs and their shade trees in managed forests are shown here to be a suitable habitat for only a limited set of orchid species. As farmers continue to convert natural forest into managed forest for coffee cultivation, further losses of habitat quality and collateral declines in regional epiphytic orchid diversity can be expected. Therefore, the conservation of epiphytic orchid diversity, as well as other components of diversity of the coffee forests, must primarily rely on avoiding coffee management intensification in the remaining natural forest. Convincing farmers to keep forest-climax trees in their coffee forest and to tolerate orchids on their coffee shrubs may also contribute to a more favorable conservation status of orchids in Ethiopian coffee agroecosystems.

Fine spatial pattern of an epiphytic lichen species is affected by habitat conditions in two forest types in the Iberian Mediterranean region

Persistence and abundance of species is determined by habitat availability and the ability to disperse and colonize habitats at contrasting spatial scales. Favourable habitat fragments are also heterogeneous in quality, providing differing opportunities for establishment and affecting the population dynamics of a species. Based on these principles, we suggest that the presence and abundance of epiphytes may reflect their dispersal ability, which is primarily determined by the spatial structure of host trees, but also by host quality. To our knowledge there has been no explicit test of the importance of host tree spatial pattern for epiphytes in Mediterranean forests. We hypothesized that performance and host occupancy in a favourable habitat depend on the spatial pattern of host trees, because this pattern affects the dispersal ability of each epiphyte and it also determines the availability of suitable sites for establishment. We tested this hypothesis using new point pattern analysis tools and generalized linear mixed models to investigate the spatial distribution and performance of the epiphytic lichen Lobaria pulmonaria, which inhabits two types of host trees (beeches and Iberian oaks). We tested the effects on L. pulmonaria distribution of tree size, spatial configuration, and host tree identity. We built a model including tree size, stand structure, and several neighbourhood predictors to understand the effect of host tree on L. pulmonaria. We also investigated the relative importance of spatial patterning on the presence and abundance of the species, independently of the host tree configuration. L. pulmonaria distribution was highly dependent on habitat quality for successful establishment, i.e., tree species identity, tree diameter, and several forest stand structure surrogates. For beech trees, tree diameter was the main factor influencing presence and cover of the lichen, although larger lichen-colonized trees were located close to focal trees, i.e., young trees. However, oak diameter was not an important factor, suggesting that bark roughness at all diameters favoured lichen establishment. Our results indicate that L. pulmonaria dispersal is not spatially restricted, but it is dependent on habitat quality. Furthermore, new spatial analysis tools suggested that L. pulmonaria cover exhibits a distinct pattern, although the spatial pattern of tree position and size was random.

Small-scale variation of vegetation in a mixed forest understorey is partly controlled by the effect of overstory composition on litter accumulation

We investigated how richness and composition of vascular plant species in the understory of a mixed hardwood forest stand varied with respect to the abundance and composition of the overstory. The stand is in central Spain and represents the southernmost range of distribution of several tree and herbaceous species in Europe. Understory species were identified in 46 quadrats (0.25 m2) where variables litter depth and light availability were measured. In addition, we estimated tree density, basal area, and percent basal area by tree species within 6-m-radius areas around each plot. Species richness and composition were studied using path analysis and scale-dependent geostatistical methods, respectively. We found that the relative abundance of certain trees species in the overstory was more important than total overstory abundance in explaining understory species richness. Richness decreased as soil litter depth increased, and soil litter increased as the relative proportion of Fagus sylvatica in the overstory increased, which accounted for a negative, indirect effect of Fagus sylvatica on richness. Regarding understory species composition, we found that some species distributed preferentially below certain tree species. For example, Melica uniflora was most frequent below Fagus sylvatica and Quercus petraea while the increasing proportion of Q. pyrenaica in the overstory favored the presence of Cruciata glabra, Arenaria montana, Prunus avium, Conopodium bourgaei, Holcus mollis, Stellaria media and Galium aparine in the understory. Overall, these results emphasize the importance of individual tree species in controlling the assemblage and richness of understory species in mixed stands. We conclude that soil litter accumulation is one way through which overstory composition shapes the understory community.