Frédéric Mortier

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−1 greater than or equal to 80 cm DBH) and thinning (11–41 trees ha−1 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−1 and 79.7 ± 45.9 m3 ha−1, respectively. We found that (i) natural control forest was increasing in AGB (2.58 ± 1.73 Mg dry mass ha−1 yr−1) and decreasing in timber stock (−0.33 ± 1.57 m3 ha−1 yr−1); (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−1 yr−1) and tripled in the logged + thinned plots (8.03 ± 1.41 Mg ha−1 yr−1); (iii) the timber stock recovered slowly (mean recovery after 24 years: 41%), at a rate of 0.75 ± 0.51 m3 ha−1 yr−1 in the logged plots, and 0.81 ± 0.74 m3 ha−1 yr−1 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 km2 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.

Legacy of logging roads in the Congo Basin: How persistent are the scars in forest cover?

Logging roads in the Congo Basin are often associated with forest degradation through fragmentation and access for other land uses. However, in concessions managed for timber production, secondary roads are usually closed after exploitation and are expected to disappear subsequently. Little is known about the effectiveness of this prescription and the factors affecting vegetation recovery rate on abandoned logging roads. In a novel approach we assessed logging roads as temporary elements in the forest landscape that vary in persistence depending on environmental conditions. We analyzed road persistence during the period 1986�2013 in adjacent parts of Cameroon, Central African Republic and Republic of Congo. Three successive phases of road recovery were identified on LANDSAT images: open roads with bare soil, roads in the process of revegetation after abandonment and disappeared roads no longer distinguishable from the surrounding forest. Field based inventories confirmed significant differences between all three categories in density and richness of woody species and cover of dominant herbs. We used dead-end road segments, built for timber exploitation, as sampling units. Only 6% of them were identified as being re-opened. Survival analyses showed median persistence of four years for open roads before changing to the revegetating state and 20 years for revegetating roads before disappearance. Persistence of revegetating roads was 25% longer on geologically poor substrates which might result from slower forest recovery in areas with lower levels of soil nutrient content. We highlight the contrast amongst forests growing on different types of substrate in their potential for ecosystem recovery over time after roads have been abandoned. Forest management plans need to take these constraints into account. Logging activities should be concentrated on the existing road network and sites of low soil resource levels should be spared from business-as-usual exploitation.

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.

Relationships between demography and gene flow and their importance for the conservation of tree populations in tropical forests under selective felling regimes

Determining how tropical tree populations subject to selective felling (logging) pressure may be conserved is a crucial issue for forest management and studying this issue requires a comprehensive understanding of the relationships between population demography and gene flow. We used a simulation model, SELVA, to study (1) the relative impact of demographic factors (juvenile mortality, felling regime) and genetic factors (selfing, number and location of fathers, mating success) on long-term genetic diversity; and (2) the impact of different felling regimes on population size versus genetic diversity. Impact was measured by means of model sensitivity analyses. Juvenile mortality had the highest impact on the number of alleles and genotypes, and on the genetic distance between the original and final populations. Selfing had the greatest impact on observed heterozygote frequency and fixation index. Other factors and interactions had only minor effects. Overall, felling had a greater impact on population size than on genetic diversity. Interestingly, populations under relatively low felling pressure even had a somewhat lower fixation index than undisturbed populations (no felling). We conclude that demographic processes such as juvenile mortality should be modelled thoroughly to obtain reliable long-term predictions of genetic diversity. Mortality in selfed and outcrossed progenies should be modelled explicitly by taking inbreeding depression into account. The modelling of selfing based on population rate appeared to be oversimplifying and should account for inter-tree variation. Forest management should pay particular attention to the regeneration capacities of felled species.

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.

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 models 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.

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.

Supervised component generalized linear regression using a PLS-extension of the Fisher scoring algorithm

In the current estimation of a GLM model, the correlation structure of regressors is not used as the basis on which to lean strong predictive dimensions. Looking for linear combinations of regressors that merely maximize the likelihood of the GLM has two major consequences: (1) collinearity of regressors is a factor of estimation instability, and (2) as predictive dimensions may lean on noise, both predictive and explanatory powers of the model are jeopardized. For a single dependent variable, attempts have been made to adapt PLS regression, which solves this problem in the classical Linear Model, to GLM estimation. In this paper, we first discuss the methods thus developed, and then propose a technique, Supervised Component Generalized Linear Regression (SCGLR), that combines PLS regression with GLM estimation in the multivariate context. SCGLR is tested on both simulated and real data.

Hierarchical Bayesian Model for Gaussian, Poisson and Ordinal Random Fields

As most georeferenced data sets are multivariate and concern variables of different kinds, spatial mapping methods must be able to deal with such data. The main difficulties are the prediction of non Gaussian variables and the modelling of the dependence between processes. The aim of this paper is to propose a new approach that permits simultaneous modelling of Gaussian, count and ordinal spatial processes. We consider a hierarchical model implemented within a Bayesian framework. The method used for Gaussian and count variables is based on the generalized linear mixed models. Ordinal variable is taken into account through a generalization of the ordinal probit model. We use a moving average approach to model the spatial dependence between the processes. The proposed model is applied to pedological data collected in French Guiana.