The association between rainforest disturbance and recovery, tree community composition, and community traits in the Yangambi area in the Democratic Republic of the Congo

Despite their key role in biodiversity conservation, forests in the Congo Basin are increasingly threatened by human activities, but it remains challenging to assess the impact of forest degradation under a more or less intact canopy. Likewise, forest recovery following agricultural abandonment remains poorly understood in the Congo Basin. Here, we surveyed 125 vegetation quadrats across 25 forest inventory plots in the Yangambi area. We aimed to find associations between both selective logging and forest recovery, and a range of forest community and tree community trait characteristics, as compared to reference undisturbed old-growth forest. We found that plots in undisturbed old-growth forest harboured both more tree individuals and tree species with a higher wood density as compared to plots in disturbed old-growth forest. In addition, their tree community composition was significantly different, whereas species diversity recovered since relatively recent agricultural abandonment ( < 60 years), community composition and forest structure remained significantly different from the plots in undisturbed old-growth forest. Our study provides some insights into the rate of forest recovery in the Congo basin after agricultural abandonment and highlights the need of proper conservation of the remaining relatively undisturbed old-growth forests. Finally, we stress the need for more extensive vegetation surveys in the Congo Basin to further unravel the effects of anthropogenic disturbance.


Introduction
The African continent is home to the world's second largest area of tropical forests which are, together with earth's subtropical forests, by far the most floral-rich biomes on earth (Kier et al. 2005, Malhi et al. 2013).The African tropical forests appear to be less biodiverse than their Asian and American counterparts (Kier et al. 2005, Parmentier et al. 2007, Slik et al. 2015).Nevertheless, Sosef et al. (2017) recorded 22,577 plant species in tropical African forests, of which 30% are endemic, and concluded that much botanical exploration still remains to be done.The forests of the Congo Basin, which account for 89% of the African tropical forests, are also extremely important as carbon sinks, removing up to six times more carbon than the forests in the Amazon region (Harris et al. 2021, Malhi et al. 2013).Furthermore, the African carbon sink has been stable for the past three decades until 2015, while the Amazonian carbon sink has been declining since the 1990s (Hubau et al. 2020).
Despite their key role in biodiversity conservation and carbon sequestration, forests in tropical Africa are increasingly threatened by human activities (Tyukavina et al. 2018).Between 1990 and 2000, the annual net loss of forests in Africa corresponded to 3.3 million ha yr −1 , increasing to 3.4 million ha yr −1 between 2000 and 2010, and further to 3.9 million ha yr −1 between 2010 and 2020 (FAO & UNEP 2020).Africa is currently the only region where annual forest change rates are still increasing (FAO & UNEP 2020).The deforestation on the African continent is mainly caused by small-scale subsistence agriculture (Curtis et al. 2018, Jayathilake et al. 2021, Tegegne et al. 2016, Tyukavina et al. 2018), driven by very low agricultural yields (Carrasco et al. 2017, West et al. 2010), and a rapidly growing human population (Gerland et al. 2014).
Although deforestation is a very conspicuous and severe threat to biodiversity, degradation of tropical forests may even be more detrimental, based on the large spatial scales at which it occurs (Barlow et al. 2016).Forest degradation refers to within-forest disturbance and mainly includes selective logging, removal of understorey vegetation, and overexploitation of specific species (Sasaki & Putz 2009, Chaturvedi et al. 2012, Tyukavina et al. 2018).In tropical Africa, forest degradation is primarily driven by selective logging, mainly for the small-scale production of charcoal which is still the most commonly used energy source in many Central African countries (Hosonuma et al. 2012, Jayathilake et al. 2021, Tyukavina et al. 2018).The degradation rate of forests in the Congo Basin has been estimated to be 317,000 ha yr −1 between 2000 and 2005 (Ernst et al. 2013).Lambin et al. (2003) estimated that the annual rate of human-mediated forest degradation in Africa corresponds to approximately 50% of its annual rate of deforestation.This would imply that 1.95 million ha of African forest has been degraded between 2010 and 2020.However, Shapiro et al. (2021) recently reported that this assumption may drastically underestimate the actual rate.Using a forest condition index, they estimated that between 2000 and 2016, over 23 million ha of forest has been degraded in the Congo Basin alone.
Although remote-sensing techniques are constantly improving, it remains challenging to monitor the extent and impact of forest degradation due to the interference with the often untouched canopy in degraded forests (Matricardi et al. 2020, Peres et al. 2006).Evaluation of the extent and impact of forest degradation must therefore largely rely on vegetation surveys on the ground (Heym et al. 2021).Impacts of forest degradation on biodiversity have been frequently studied in the Neotropics (e.g.Carreño-Rocabado et al. 2012, Mirabel et al. 2020), Asia (e.g.Baboo et al. 2017, Millet et al. 2010), and other parts of Africa (e.g.Bobo et al. 2006, Oyugi et al. 2008), but similar studies are rare in the Congo Basin (but see Hall et al. 2003, Ifo et al. 2016), and almost lacking in the Democratic Republic of Congo (DR Congo) (but see Makana & Thomas 2006).Yet, as DR Congo harbours almost 60% of the African tropical forests (Malhi et al. 2013), it is essential to comprehend how plant community composition is affected by forest disturbance.
While conservation of existing tropical forests should still be the priority, the restoration of deforested land and degraded forests can aid in the recovery of tropical biodiversity and the recovery of ecosystem services (Brudvig 2011, FAO & UNEP 2020).Between 2000 and 2019, 0.7 million ha of forest has been actively or passively restored in Sub-Saharan Africa (FAO & UNEP 2020).By 2030, the African Forest Landscape Restoration Initiative (AFR100) aims to increase the land surface under restoration to 100 million ha.DR Congo has committed to restore 8 million ha (AFR100 2021).Because of the high variability in restoration outcomes, more research is crucial to clear the path for restoration (Chazdon 2008).Again, studies on the biodiversity consequences of tropical forest recovery are more commonly found in the Neotropics (e.g.Letcher & Chazdon 2009, Oberleitner et al. 2021) and tropical Asia (e.g.Hayward et al. 2021, Mahayani et al. 2020), but are rare in the Congo Basin (but see Bauters et al. 2019, Makelele et al. 2021).
In this study, we established a set of forest inventory plots in the Yangambi region in DR Congo and aimed to find associations between both anthropogenic disturbance through selective logging and forest recovery on the one side, and a range of forest and tree community characteristics on the other side, as compared to reference undisturbed old-growth forest.Based on 125 sampling quadrats across 25 forest inventory plots, we specifically quantified differences among (i) disturbed old-growth forest, (ii) regrowth forest, and (iii) undisturbed old-growth forest in terms of tree species diversity, community traits, community composition, forest structure, and canopy cover.We hypothesise that (i) selective logging results in the decreasing abundance of some specific tree species and in shifts in community composition, without having a profound effect on the overall tree species diversity, (ii) that selective logging results in a lower community-weighted mean wood density and in lower crown closure, and (iii) that forest recovery following agricultural abandonment will be fast for secondary tree species, with low wood density, but slower for primary forest tree species, with higher wood density, implying that the overall tree species diversity will have recovered, but that the community composition will be different from undisturbed old-growth forest.We established 25 forest inventory plots of 125 m × 125 m (1.56 ha), covering an area of ca.50 by 20 km, just North of the Congo River (Figure 1).Plots were assigned to three different forest categories by georeferencing (WGS84) them onto historical landuse maps, which were issued by Evrard (1954).Plots located on historical agricultural land from INERA were classified as regrowth forest (R).Although we could not exactly date the start of the forest recovery, our best estimate based on local sources is that the agricultural land was abandoned somewhere between 40 and 60 years ago.Plots located in historically old-growth forest areas but where we observed indications of considerable small-scale selective logging through the presence of tree stumps and other remnants of logging such as parts of trees or active logging were considered as disturbed old-growth forest (DO).These areas were typically located near regrowth forest and close to places with relatively higher population density.The remaining plots that were located in areas with historical old-growth forests, and where no signs of selective logging were present, were categorised as undisturbed old-growth forest (UO).The current UO and DO forests have previously been classified by Evrard (1954) as heterogenous forests with primitive character.These old vegetation maps have been shown to be very reliable (Vancutsem et al. 2009), and we may therefore assume there were no pre-existing differences between both forest categories.

Vegetation survey and plant trait quantification
Five quadrats (25 m × 25 m) were randomly established in each forest inventory plot, resulting in a total of 125 sampling quadrats.In each quadrat, the abundance and diameter at breast height (DBH) of all woody species with a DBH ≥ 5 cm were measured.The quadrats were geolocated (Garmin GPSMAP 64 series), and crown closure (%) was quantified from four readings in the cardinal directions with a model A convex spherical densiometer.
The specific leaf area (SLA) of all taxa identified at species level was assessed based on dry herbarium specimens from the Meise Botanic Garden herbarium following Pérez-Harguindeguy et al. (2013).Perez et al. (2020) provided evidence that the SLA can reliably be estimated from herbarium-based measurements.Leaf area was assessed using leaf disks (4.9 cm 2 ), to minimise damage to herbarium specimens, or otherwise by a whole leaf scan using ImageJ (Rueden et al. 2017).The obtained disks and whole leaves were weighted afterwards.Five replicates were used per species.In addition to the SLA, the wood density (WD) of 141 species was compiled from the literature (Table S1) (RMCA 2021, Carsan et al. 2012, Harja et al. 2021).

Data analyses
Statistical analyses were performed using R software (R Core Team 2020).The data of the five quadrats were summed into one value per plot.For each plot, the number of trees, species richness (Hill's N 0 ), species diversity (Hill's N 1 or the exponential of the Shannon entropy), and evenness (Hill's N 1 /log(Hill's N 0 )) were then calculated.These numbers were then standardised for sample completeness, as sample sizes differed between plots.Sample completeness represents the proportion of the total number of individuals in a unit that belongs to the species in the sample, contrary to sample size standardisation which allows less biased comparisons (Chao & Jost 2012).Extrapolation of the estimates was done to twice the reference sample size per plot, following Chao & Jost (2012).Calculations were based on the collected species abundance data with 9999 bootstraps using the iNEXT function of the iNEXT package (Hsieh et al. 2016).Subsequently, the different variables and the crown closure were tested for differences between the three forest categories using Kruskal-Wallis rank sum tests and pairwise Wilcoxon rank sum tests with correction for multiple testing, as the assumptions for parametric testing were not met.
Dissimilarity of the tree community composition across the forest categories was visualised using a non-metric multidimensional scaling (NMDS), based on the Bray-Curtis distance matrix, the most suitable distance for abundance data (Ricotta et al. 2017).A two-axis solution provided a stress level (0.13) well below the 0.20 threshold, indicating a good representation in reduced dimensions.The ordination was performed with the metaMDS function in the vegan package (Oksanen et al. 2020).Subsequently, a Mantel test with 9999 permutations was performed to test the effect of geographic distance on the dissimilarity of the tree community composition using the mantel function in the Vegan package (Oksanen et al. 2020).In addition to the NMDS, a permutational multivariate analysis of variance with 9999 permutations was performed to test for differences in community composition between the three different forest categories using the adonis function in the Vegan package (Oksanen et al. 2020).Pairwise comparisons of the different forest groups with 9999 permutations and Bonferroni correction were performed with the pairwise.adonis.function in the pairwiseAdonis package (Arbizu 2017).Tree species indicative for each of the three forest categories were identified based on their indicator value (Dufrene & Legendre 1997).This value was calculated using the multipatt function in the indicspecies package (De Cáceres & Legendre 2009), with 9999 permutations.
Community-weighted means (CWM) for each plot were calculated for the SLA and the WD, using the number of individuals as weight.The community-weighted means are known to be very responsive to environmental change and highly sensitive to disturbance (Carreño-Rocabado et al. 2012).Both variables met the assumptions for parametric testing, and to evaluate differences in community-weighted mean SLA and WD among forest categories, a one-way analysis of variance (ANOVA) was conducted followed by pairwise Tukey post hoc tests.
Forest structure was evaluated using the measured DBHs.Trees were categorised in ten size classes according to their DBH.An interval of 10 cm was used in the first nine classes, starting from 5 cm DBH.The tenth class had a wider interval, aggregating individuals with a DBH over 95 cm.Subsequently, the diameter frequency distribution was compared among the three forest categories with Pearson's chi-squared tests with simulated p-values based on 9999 replicates using the chisq.testfunction in the stats package (R Core Team 2020).The chi-square test was preferred over the frequently used Kolmogorov-Smirnov test, because the former is more flexible than the latter (Mitchell 1971).The diameter frequency distribution comparison was complemented by the calculation of the structural diversity, which is a single measure summarising the diameter distribution (McElhinny et al. 2005).Structural diversity was quantified by the first order Hill number which was standardised for sample completeness following the same methodology as outlined above.Calculation and extrapolation of structural diversity was done with the iNEXT function with 9999 bootstraps in the iNEXT package, as above.All figures were made using the ggplot2 package (Wickham 2016).

Community composition
The NMDS of the 25 forest inventory plots showed a spatial separation of the three forest categories (Figure 3).This was supported by the permutational multivariate analysis of variance which indicated differences between categories (R 2 = 0.253; P < 0.001).Specifically, community composition significantly differed between DO forests and UO forests (R 2 = 0.207; P < 0.001), and between R forests and UO forests (R 2 = 0.267; P < 0.001).The Mantel test between the geographic distance matrix and the Bray-Curtis distance matrix showed a significant but weak correlation (r = 0.200; P < 0.05).Across the 25 forest inventory plots, 27 out of the 167 at species level identified trees were indicative based on their indicator value (both P < 0.05 and indicator value > 0.20).Seventeen indicator species were indicative for UO forest (e.g.Cola griseiflora), four for DO forest (e.g.Leplaea cedrata), and six for R forest (e.g.Heisteria parvifolia) (Table S3).

Specific leaf area and wood density
The community-weighted means for SLA did not differ between forest categories (F = 1.38,P > 0.05); however, the communityweighted mean WD exhibited a significant difference (F = 9.03, P < 0.01) (Figure 4).In DO forests, community-weighted mean WD was 5.0 % lower (P < 0.05), whereas in R forests it was 6.2 % lower (P < 0.01) as compared to UO forests.

Discussion
Understanding how anthropogenic disturbance affects forest structure, tree species richness, and community composition is of major importance for biodiversity conservation in globally underexplored biodiversity hotspots such as the Congo Basin.The goal of our study was to find associations between both anthropogenic disturbance through selective logging and forest recovery on the one side, and a range of forest community and tree community trait characteristics on the other side, as compared to reference undisturbed old-growth forest.Although our three delineated forest categories certainly enveloped some within-variation in terms of time since abandonment (for the regrowth forest plots) and in terms of degree of disturbance (for the disturbed old-growth forest plots), we were able to identify differences in community composition, forest structure, and community traits among forest categories.

Selective logging
As hypothesised, our survey showed that there were clear differences in community composition and forest structure between plots in undisturbed old-growth forest and plots in disturbed old-growth forest.There were no differences in the number of woody species between any forest category; yet contrary to our expectation, species diversity was higher in plots in disturbed oldgrowth forest.This higher species diversity likely indicates that the species abundance distribution is less skewed in disturbed forests.Indeed, the evenness index showed that the species abundance is more evenly distributed in plots in disturbed forest than in plots in undisturbed old-growth forest.Combined with the observation that plots in disturbed old-growth forests contained significantly less tree individuals, these findings are in line with our hypothesis and might suggest that selective logging, in which exploitable species are selectively removed, is closely associated with these patterns.In the Yangambi area, timber is removed by local villagers using, for example, machetes, mainly for construction material and cooking purposes (Depecker pers. obs.).The stems that are removed are usually from species from which the wood is exploitable.We expect that selective logging at first only results in a shift in densities, but that it will ultimately result in the loss of certain species as selective loggers will specifically target species with stems that are exploitable.For example, Prioria oxyphylla and Scorodophloeus zenkeri were much less abundant in disturbed forests and are used as construction material.Moreover, in DR Congo the latter species is also exploited for its medicinal properties and is used to make charcoal (PROTA 2021).Other species that were highly reduced in numbers were Dialium pachyphyllum, Greenwayodendron suaveolens, and Strombosia pustulata, which all bear seeds that are dispersed by animals (PROTA 2021, Meunier et al. 2015).Animal-dispersed seeded tree species tend to have higher wood densities than abiotic-dispersed seeded tree species (Osuri et al. 2016), possibly indicating that loggers specifically target species with high WD, suggesting that species with animal-dispersed seeds are more vulnerable to selective logging.Careful monitoring of tree species is needed to unravel the specific effects of selective logging at the species level.Plots in disturbed old-growth forest on the other hand still contained many shade-tolerant tree species of which the wood is nonexploitable, such as species of the genus Cola which have a low DBH and highly irregular trunks.Cola acuminata also produces highly valued nuts with stimulating and medicinal properties  4. Comparisons between the community-weighted means of specific leaf area (SLA) and wood density (WD) across the three forest categories.Hinges represent the 25 th , 50 th , and 75 th percentiles, respectively.Whiskers extend to maximum 1.5 times the interquartile range.Letters code for significant differences between groups.( Meunier et al. 2015).Likewise, the abundance of most Diospyros spp. is indicative for old-growth forests, of which only Diospyros crassiflora is exploited for its wood (Meunier et al. 2015).Although a clear difference was found in community composition between forest categories in our study, the response in terms of species diversity and species richness was more variable.An earlier study on selective logging in the Congo Basin compared disturbed and undisturbed forests in terms of species diversity in two distinct locations (Makana & Thomas 2006).Although a significant influence of disturbance was detected in the first location, similar values for species richness and species diversity were detected in the second location.Furthermore, a clear effect of disturbance, primarily selective logging, on species diversity was detected in tropical forests in the Republic of the Congo (Ifo et al. 2016) and India (Baboo et al. 2017), but not in a tropical forest in Bolivia (Carreño-Rocabado et al. 2012).It is clear that solely assessing the possible impact of anthropogenic disturbance through species diversity can mask its impact and that community composition is a better indicator of the possible impacts of forest degradation (Edwards et al. 2014, Imai et al. 2014).Indeed, although several studies across different tropical regions have reported varying responses of species diversity to forest disturbance, they all reported significant shifts in the community composition following anthropogenic disturbance (Makana & Thomas 2006, Baboo et al. 2017, Carreño-Rocabado et al. 2012, Ifo et al. 2016) Because ecosystem functioning is mainly mediated by species' functional traits (Chaturvedi et al. 2021, de Bello et al. 2021), community-weighted mean functional traits can provide better ecologically information than only using taxonomic data (Katovai et al. 2012).We hypothesised that we would find lower community-weighted mean wood density and a lower canopy closure in plots in disturbed old-growth forest, as compared to plots in undisturbed old-growth forest.In terms of canopy closure, we did not find differences between any forest category.Consequently, the light availability in the plots remained unchanged, which presumably clarifies why we observed no significant differences in community-weighted mean specific leaf area between forest categories.It has indeed been established by Carreño-Rocabado et al. (2012) that increased light availability can induce changes in community-weighted mean functional traits in tropical forests.Furthermore, although disturbed and undisturbed old-growth forests were characterised by different indicator species assemblages, the majority can still be classified as shade-tolerant, which is commonly associated with low SLA (Poorter 2009).Nonetheless, the absence of differences between forest categories should be interpreted cautiously as SLA is highly responsive to environmental variation, such as sun exposure.In contrast to communityweighted mean SLA, community-weighted mean WD did significantly differ between plots in disturbed and plots in undisturbed old-growth forest.As anticipated, a decrease in communityweighted mean WD was observed in plots in disturbed old-growth forest, which indicates that trees with high WD and of high qualities as construction material have been selectively removed (Saranpää 2003).
The mean diameter distribution in the forests in the Yangambi region has previously been described by the reverse J-shaped model (Fayolle et al. 2014).In our study, the reverse J-shape was more pronounced in plots in undisturbed old-growth forests than in plots in disturbed old-growth forests which contained significantly less trees in the first size class, suggesting less regeneration (Gebeyehu et al. 2019).Similar results of reduced regeneration in forests were obtained in other tropical forests and were assigned to disturbance (Makana & Thomas 2006, Baboo et al. 2017, Hall et al. 2003).

Forest recovery following agricultural abandonment
Based on the historical land-use maps, our results indicated that after 40 to 60 years there were no significant differences between plots in regrowth forests and plots in undisturbed old-growth forest in terms of species diversity, which was hypothesised.They only differed in terms of the number of tree individuals and community-weighted mean WD.These findings are in line with other studies in the Congo Basin (Bauters et al. 2019, Makelele et al. 2021) and in other tropical regions (Letcher & Chazdon 2009, Mahayani et al. 2020).In contrast to species diversity, we found that the community composition in plots in regrowth forest was significantly different from those in undisturbed old-growth forest.These findings are in line with what was expected and might suggest that community composition needs more time to converge in similarity to plots in undisturbed old-growth forest.Similar to our results, the community composition of young regrowth forest (25-30 years) differed from old-growth forest in the Maringa-Lopori-Wamba forest landscape in DR Congo, whereas old regrowth forests (150-300 years) were very close to old-growth forests (Bauters et al. 2019).Likewise, a recent study in the Yoko forest reserve, at approximately 30 km Southeast of Kisangani, showed that the community composition of 60-year-old regrowth forests was still significantly different from local old-growth forests (Makelele et al. 2021).Our results are furthermore in line with studies in the Neotropics (Oberleitner et al. 2021, Poorter et al. 2021) and tropical Asia (Hayward et al. 2021).
Despite having a similar community composition as plots in disturbed old-growth forests, plots in regrowth forests were characterised by a different indicator species assemblage, including Myrianthus arboreus, a late pioneer species (PROTA 2021).Additionally, Combretum lokele and Dacryodes edulis were also characteristic for regrowth forests.Both species are commonly found in villages and are indicative for the important anthropogenic influence in these forests (Meunier et al. 2015).
The forest structure in the plots in regrowth forests was comparable to plots in disturbed old-growth forest, but significantly different from plots in undisturbed old-growth forest.These differences were not caused by the lower abundance of large-sized trees in regrowth forests, as would be expected because these are slower to recover (Bauters et al. 2019, Makelele et al. 2021).On the contrary, the observed differences can be explained by the higher frequency of small-sized trees in plots in undisturbed old-growth forests.This possibly implies that the tree regeneration could be hampered.Similar to other studies (Bauters et al. 2019, Makelele et al. 2021), we suggest that forests in the Yangambi landscape initially recover rapidly to a close canopy forest with shadetolerant understorey species following agricultural abandonment, but that the recruitment of primary forest species is slow, and that further forest succession may be hampered by ongoing anthropogenic disturbance in these regrowth forests.

Limitations
We acknowledge that our study is not flawless and that there are some limitations that have to be considered.Our plots are to a certain extent spatially clustered, which was not possible to avoid due to the location of the disturbance zone in the Yangambi area.Furthermore, due to reasons of feasibility, it was not possible to incorporate additional control plots to account for potential spatial variability (Davies & Gray 2015).The Mantel test was used to test the effect of geographic distance on the dissimilarity of the tree community composition.It should be noted that the Mantel test showed a significant correlation between geographic distance and dissimilarity of community composition.This distance decay of similarity is commonly reported in tropical forests (Ganivet et al. 2020, Liebsch et al. 2008, Wittmann et al. 2006).Pre-existing spatial variation seems to confound with the main treatment (logging), but the correlation is rather weak.Although we cannot conclude that the observed differences between sites are a direct effect of selective logging per se, we can interpret the results in light of our ecological understanding of the area and its activities (Oksanen 2004, Davies & Gray 2015).

Conclusion
We found that the tree community composition and forest structure were significantly different between forest categories and that the response of species diversity was more variable.Based on our knowledge of the activities in the Yangambi region, our results suggest that these differences in community composition found between UO forest and DO forest can be attributed to selective logging.Our data furthermore suggest that 40 to 60 years after agricultural abandonment, species diversity is able to recover, but community composition and forest structure are not.Given the limitations of our sampling design, we suggest vegetation surveys across larger geographic areas to deal with the spatial clustering of our survey plots within one disturbance category.Based on the known correlation between forest productivity and SLA, and between carbon storage and WD, our results also indicate that disturbed old-growth forest and regrowth forest are as productive as undisturbed old-growth forest but are not able to store as much carbon (Chave et al. 2009, Madani et al. 2018).Conservation of undisturbed old-growth forests is vital, although also regrowth of forests can play a crucial role in maintaining biodiversity and providing ecosystem services (Crouzeilles et al. 2016, Van de Perre et al. 2018).
Supplementary material.For supplementary material accompanying this paper visit https://doi.org/10.1017/S0266467422000347Data Availability Statement.The data that support the findings of this study will be made openly available in Zenodo

Figure 1 .
Figure 1.Location of the 25 forest inventory plots across the three forest categories in the Yangambi region (DR Congo).

Figure 2 .Figure 3 .
Figure2.Comparisons between tree abundance and tree species diversity metrics across old-growth undisturbed forest (OU), old-growth disturbed forest (OD), and regrowth forest (R) in the Yangambi region (DR Congo).Hinges represent the 25th, 50 th , and 75th percentiles, respectively.Whiskers extend to maximum 1.5 times the interquartile range.Letters code for significant differences between groups.(a) Tree abundance, (b) species richness, (c) species diversity, (d) evenness.

Figure
Figure4.Comparisons between the community-weighted means of specific leaf area (SLA) and wood density (WD) across the three forest categories.Hinges represent the 25 th , 50 th , and 75 th percentiles, respectively.Whiskers extend to maximum 1.5 times the interquartile range.Letters code for significant differences between groups.

Figure 5 .
Figure 5. Observed frequencies among the ten DBH size categories compared among the three forest categories.Size classes one to nine have a 10 cm interval, while the tenth has a wider interval, grouping individuals with a DBH of over 95 cm.