40 years of forest dynamics and tree demography in an intact tropical forest at M'Baïki in central Africa.

A vast silvicultural experiment was set up in 1982 nearby the town of M'Baïki in the Central African Republic to monitor the recovery of tropical forests after disturbance. The M'Baïki experiment consists of ten 4-ha Permanent Sample Plots (PSPs) that were assigned to three silvicultural treatments in 1986 according to a random block design. In each plot, all trees with a girth at breast height greater than 30 cm were spatially located, numbered, measured, and determined botanically. Girth, mortality and newly recruited trees, were monitored almost annually over the 1982-2022 period with inventory campaigns for 35 years. The data were earlier used to fit growth and population models, to study the species composition dynamics, and the effect of silvicultural treatments on tree diversity and aboveground biomass. Here, we present new information on the forest stand structure dynamics and tree demography. The data released from this paper cover the three control plots and constitute a major contribution for further studies about the biodiversity of intact tropical forests.

Pólya-splitting distributions as stationary solutions of multivariate birth-death processes under extended neutral theory.

Multivariate count distributions are crucial for the inference of ecological processes underpinning biodiversity. In particular, neutral theory provides useful null distributions allowing the evaluation of adaptation or natural selection. In this paper, we build a broader family of multivariate distributions: the Polya-splitting distributions. We show that they emerge naturally as stationary distributions of a multivariate birth-death process. This family of distributions is a consistent extension of non-zero sum neutral models based on a master equation approach. It allows considering both total abundance of the community and relative abundances of species. We emphasize that this family is large enough to encompass various dependence structures among species. We also introduce the strong closure under addition property that can be useful to generate nested multi-level dependence structures. Although all Pólya splitting distributions do not share this property, we provide numerous example verifying it. They include the previously known example with independent species, and also new ones with alternative dependence structures. Overall, we advocate that Polya-splitting distribution should become a part of the classic toolbox for the analysis of multivariate count data in ecology, providing alternative approaches to joint species distribution framework. Comparatively, our approach allows to model dependencies between species at the observation level, while the classical JSDM's model dependencies at the latent process strata.

Unveiling African rainforest composition and vulnerability to global change.

Africa is forecasted to experience large and rapid climate change1 and population growth2 during the twenty-first century, which threatens the world's second largest rainforest. Protecting and sustainably managing these African forests requires an increased understanding of their compositional heterogeneity, the environmental drivers of forest composition and their vulnerability to ongoing changes. Here, using a very large dataset of 6 million trees in more than 180,000 field plots, we jointly model the distribution in abundance of the most dominant tree taxa in central Africa, and produce continuous maps of the floristic and functional composition of central African forests. Our results show that the uncertainty in taxon-specific distributions averages out at the community level, and reveal highly deterministic assemblages. We uncover contrasting floristic and functional compositions across climates, soil types and anthropogenic gradients, with functional convergence among types of forest that are floristically dissimilar. Combining these spatial predictions with scenarios of climatic and anthropogenic global change suggests a high vulnerability of the northern and southern forest margins, the Atlantic forests and most forests in the Democratic Republic of the Congo, where both climate and anthropogenic threats are expected to increase sharply by 2085. These results constitute key quantitative benchmarks for scientists and policymakers to shape transnational conservation and management strategies that aim to provide a sustainable future for central African forests.

Spatial validation reveals poor predictive performance of large-scale ecological mapping models.

Mapping aboveground forest biomass is central for assessing the global carbon balance. However, current large-scale maps show strong disparities, despite good validation statistics of their underlying models. Here, we attribute this contradiction to a flaw in the validation methods, which ignore spatial autocorrelation (SAC) in data, leading to overoptimistic assessment of model predictive power. To illustrate this issue, we reproduce the approach of large-scale mapping studies using a massive forest inventory dataset of 11.8 million trees in central Africa to train and validate a random forest model based on multispectral and environmental variables. A standard nonspatial validation method suggests that the model predicts more than half of the forest biomass variation, while spatial validation methods accounting for SAC reveal quasi-null predictive power. This study underscores how a common practice in big data mapping studies shows an apparent high predictive power, even when predictors have poor relationships with the ecological variable of interest, thus possibly leading to erroneous maps and interpretations.

A map of African humid tropical forest aboveground biomass derived from management inventories.

Forest biomass is key in Earth carbon cycle and climate system, and thus under intense scrutiny in the context of international climate change mitigation initiatives (e.g. REDD+). In tropical forests, the spatial distribution of aboveground biomass (AGB) remains, however, highly uncertain. There is increasing recognition that progress is strongly limited by the lack of field observations over large and remote areas. Here, we introduce the Congo basin Forests AGB (CoFor-AGB) dataset that contains AGB estimations and associated uncertainty for 59,857 1-km pixels aggregated from nearly 100,000 ha of in situ forest management inventories for the 2000 - early 2010s period in five central African countries. A comprehensive error propagation scheme suggests that the uncertainty on AGB estimations derived from c. 0.5-ha inventory plots (8.6-15.0%) is only moderately higher than the error obtained from scientific sampling plots (8.3%). CoFor-AGB provides the first large scale view of forest AGB spatial variation from field data in central Africa, the second largest continuous tropical forest domain of the world.

Effects of logging on roadless space in intact forest landscapes of the Congo Basin.

Forest degradation in the tropics is often associated with roads built for selective logging. The protection of intact forest landscapes (IFL) that are not accessible by roads is high on the biodiversity conservation agenda and a challenge for logging concessions certified by the Forest Stewardship Council (FSC). A frequently advocated conservation objective is to maximize the retention of roadless space, a concept that is based on distance to the nearest road from any point. We developed a novel use of the empty-space function - a general statistical tool based on stochastic geometry and random sets theory - to calculate roadless space in a part of the Congo Basin where road networks have been expanding rapidly. We compared the temporal development of roadless space in certified and uncertified logging concessions inside and outside areas declared IFL in 2000. Inside IFLs, road-network expansion led to a decrease in roadless space by more than half from 1999 to 2007. After 2007, loss leveled out in most areas to close to 0 due to an equilibrium between newly built roads and abandoned roads that became revegetated. However, concessions in IFL certified by FSC since around 2007 continuously lost roadless space and reached a level comparable to all other concessions. Only national parks remained mostly roadless. We recommend that forest-management policies make the preservation of large connected forest areas a top priority by effectively monitoring - and limiting - the occupation of space by roads that are permanently accessible.

Wood Specific Gravity Variations and Biomass of Central African Tree Species: The Simple Choice of the Outer Wood.

CONTEXT: Wood specific gravity is a key element in tropical forest ecology. It integrates many aspects of tree mechanical properties and functioning and is an important predictor of tree biomass. Wood specific gravity varies widely among and within species and also within individual trees. Notably, contrasted patterns of radial variation of wood specific gravity have been demonstrated and related to regeneration guilds (light demanding vs. shade-bearing). However, although being repeatedly invoked as a potential source of error when estimating the biomass of trees, both intraspecific and radial variations remain little studied. In this study we characterized detailed pith-to-bark wood specific gravity profiles among contrasted species prominently contributing to the biomass of the forest, i.e., the dominant species, and we quantified the consequences of such variations on the biomass. METHODS: Radial profiles of wood density at 8% moisture content were compiled for 14 dominant species in the Democratic Republic of Congo, adapting a unique 3D X-ray scanning technique at very high spatial resolution on core samples. Mean wood density estimates were validated by water displacement measurements. Wood density profiles were converted to wood specific gravity and linear mixed models were used to decompose the radial variance. Potential errors in biomass estimation were assessed by comparing the biomass estimated from the wood specific gravity measured from pith-to-bark profiles, from global repositories, and from partial information (outer wood or inner wood). RESULTS: Wood specific gravity profiles from pith-to-bark presented positive, neutral and negative trends. Positive trends mainly characterized light-demanding species, increasing up to 1.8 g.cm-3 per meter for Piptadeniastrum africanum, and negative trends characterized shade-bearing species, decreasing up to 1 g.cm-3 per meter for Strombosia pustulata. The linear mixed model showed the greater part of wood specific gravity variance was explained by species only (45%) followed by a redundant part between species and regeneration guilds (36%). Despite substantial variation in wood specific gravity profiles among species and regeneration guilds, we found that values from the outer wood were strongly correlated to values from the whole profile, without any significant bias. In addition, we found that wood specific gravity from the DRYAD global repository may strongly differ depending on the species (up to 40% for Dialium pachyphyllum). MAIN CONCLUSION: Therefore, when estimating forest biomass in specific sites, we recommend the systematic collection of outer wood samples on dominant species. This should prevent the main errors in biomass estimations resulting from wood specific gravity and allow for the collection of new information to explore the intraspecific variation of mechanical properties of trees.

Statistical ecology comes of age.

The desire to predict the consequences of global environmental change has been the driver towards more realistic models embracing the variability and uncertainties inherent in ecology. Statistical ecology has gelled over the past decade as a discipline that moves away from describing patterns towards modelling the ecological processes that generate these patterns. Following the fourth International Statistical Ecology Conference (1-4 July 2014) in Montpellier, France, we analyse current trends in statistical ecology. Important advances in the analysis of individual movement, and in the modelling of population dynamics and species distributions, are made possible by the increasing use of hierarchical and hidden process models. Exciting research perspectives include the development of methods to interpret citizen science data and of efficient, flexible computational algorithms for model fitting. Statistical ecology has come of age: it now provides a general and mathematically rigorous framework linking ecological theory and empirical data.

Tropical forest recovery from logging: a 24 year silvicultural experiment from Central Africa.

Large areas of African moist forests are being logged in the context of supposedly sustainable management plans. It remains however controversial whether harvesting a few trees per hectare can be maintained in the long term while preserving other forest services as well. We used a unique 24 year silvicultural experiment, encompassing 10 4 ha plots established in the Central African Republic, to assess the effect of disturbance linked to logging (two to nine trees ha⁻¹ greater than or equal to 80 cm DBH) and thinning (11-41 trees ha⁻¹ greater than or equal to 50 cm DBH) on the structure and dynamics of the forest. Before silvicultural treatments, above-ground biomass (AGB) and timber stock (i.e. the volume of commercial trees greater than or equal to 80 cm DBH) in the plots amounted 374.5 ± 58.2 Mg ha⁻¹ and 79.7 ± 45.9 m³ ha⁻¹, respectively. We found that (i) natural control forest was increasing in AGB (2.58 ± 1.73 Mg dry mass ha⁻¹ yr⁻¹) and decreasing in timber stock (-0.33 ± 1.57 m³ ha⁻¹ yr⁻¹); (ii) the AGB recovered very quickly after logging and thinning, at a rate proportional to the disturbance intensity (mean recovery after 24 years: 144%). Compared with controls, the gain almost doubled in the logged plots (4.82 ± 1.22 Mg ha⁻¹ yr⁻¹) and tripled in the logged + thinned plots (8.03 ± 1.41 Mg ha⁻¹ yr⁻¹); (iii) the timber stock recovered slowly (mean recovery after 24 years: 41%), at a rate of 0.75 ± 0.51 m³ ha⁻¹ yr⁻¹ in the logged plots, and 0.81 ± 0.74 m³ ha⁻¹ yr⁻¹ in the logged + thinned plots. Although thinning significantly increased the gain in biomass, it had no effect on the gain in timber stock. However, thinning did foster the growth and survival of small- and medium-sized timber trees and should have a positive effect over the next felling cycle.

Geological substrates shape tree species and trait distributions in African moist forests.

BACKGROUND: Understanding the factors that shape the distribution of tropical tree species at large scales is a central issue in ecology, conservation and forest management. The aims of this study were to (i) assess the importance of environmental factors relative to historical factors for tree species distributions in the semi-evergreen forests of the northern Congo basin; and to (ii) identify potential mechanisms explaining distribution patterns through a trait-based approach. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed the distribution patterns of 31 common tree species in an area of more than 700,000 km(2) spanning the borders of Cameroon, the Central African Republic, and the Republic of Congo using forest inventory data from 56,445 0.5-ha plots. Spatial variation of environmental (climate, topography and geology) and historical factors (human disturbance) were quantified from maps and satellite records. Four key functional traits (leaf phenology, shade tolerance, wood density, and maximum growth rate) were extracted from the literature. The geological substrate was of major importance for the distribution of the focal species, while climate and past human disturbances had a significant but lesser impact. Species distribution patterns were significantly related to functional traits. Species associated with sandy soils typical of sandstone and alluvium were characterized by slow growth rates, shade tolerance, evergreen leaves, and high wood density, traits allowing persistence on resource-poor soils. In contrast, fast-growing pioneer species rarely occurred on sandy soils, except for Lophira alata. CONCLUSIONS/SIGNIFICANCE: The results indicate strong environmental filtering due to differential soil resource availability across geological substrates. Additionally, long-term human disturbances in resource-rich areas may have accentuated the observed patterns of species and trait distributions. Trait differences across geological substrates imply pronounced differences in population and ecosystem processes, and call for different conservation and management strategies.

A hierarchical bayesian model for spatial prediction of multivariate non-gaussian random fields.

As most georeferenced data sets are multivariate and concern variables of different types, spatial mapping methods must be able to deal with such data. The main difficulties are the prediction of non-Gaussian variables and the modeling of the dependence between processes. The aim of this article is to present a new hierarchical Bayesian approach that permits simultaneous modeling of dependent Gaussian, count, and ordinal spatial fields. This approach is based on spatial generalized linear mixed models. We use a moving average approach to model the spatial dependence between the processes. The method is first validated through a simulation study. We show that the multivariate model has better predictive abilities than the univariate one. Then the multivariate spatial hierarchical model is applied to a real data set collected in French Guiana to predict topsoil patterns.

Finding confidence limits on population growth rates: bootstrap and analytic methods.

When predicting population dynamics, the value of the prediction is not enough and should be accompanied by a confidence interval that integrates the whole chain of errors, from observations to predictions via the estimates of the parameters of the model. Matrix models are often used to predict the dynamics of age- or size-structured populations. Their parameters are vital rates. This study aims (1) at assessing the impact of the variability of observations on vital rates, and then on model's predictions, and (2) at comparing three methods for computing confidence intervals for values predicted from the models. The first method is the bootstrap. The second method is analytic and approximates the standard error of predictions by their asymptotic variance as the sample size tends to infinity. The third method combines use of the bootstrap to estimate the standard errors of vital rates with the analytical method to then estimate the errors of predictions from the model. Computations are done for an Usher matrix models that predicts the asymptotic (as time goes to infinity) stock recovery rate for three timber species in French Guiana. Little difference is found between the hybrid and the analytic method. Their estimates of bias and standard error converge towards the bootstrap estimates when the error on vital rates becomes small enough, which corresponds in the present case to a number of observations greater than 5000 trees.

A spatial statistical model for landscape genetics.

Landscape genetics is a new discipline that aims to provide information on how landscape and environmental features influence population genetic structure. The first key step of landscape genetics is the spatial detection and location of genetic discontinuities between populations. However, efficient methods for achieving this task are lacking. In this article, we first clarify what is conceptually involved in the spatial modeling of genetic data. Then we describe a Bayesian model implemented in a Markov chain Monte Carlo scheme that allows inference of the location of such genetic discontinuities from individual geo-referenced multilocus genotypes, without a priori knowledge on populational units and limits. In this method, the global set of sampled individuals is modeled as a spatial mixture of panmictic populations, and the spatial organization of populations is modeled through the colored Voronoi tessellation. In addition to spatially locating genetic discontinuities, the method quantifies the amount of spatial dependence in the data set, estimates the number of populations in the studied area, assigns individuals to their population of origin, and detects individual migrants between populations, while taking into account uncertainty on the location of sampled individuals. The performance of the method is evaluated through the analysis of simulated data sets. Results show good performances for standard data sets (e.g., 100 individuals genotyped at 10 loci with 10 alleles per locus), with high but also low levels of population differentiation (e.g., FST<0.05). The method is then applied to a set of 88 individuals of wolverines (Gulo gulo) sampled in the northwestern United States and genotyped at 10 microsatellites.