Crispin Ilunga-Mulala, Hulda Hatakiwe, Hans Beeckman, Olivier J. Hardy, Gauthier Ligot, Dieu-Merci Assumani, Léopold Ndjele & Nils Bourland

Influence of spacing and seed trees on the growth of Pericopsis elata saplings during the first twenty months of a planting trial

Résumé

Influence de l’écartement et du semencier sur la croissance des plants de Pericopsis elata durant les vingt premiers mois d’un essai de plantation

Description du sujet. Pericopsis elata (Fabaceae) est un arbre héliophile des forêts tropicales africaines qui produit du bois de grande valeur économique. Ses populations naturelles souffrent de surexploitation et présentent un déficit de régénération naturelle. Les plantations pourraient augmenter sa production et limiter la pression sur les forêts naturelles. Cependant, les connaissances sur l'influence de l'écartement et du semencier sur la croissance des arbres font encore défaut pour les protocoles de plantation.

Objectifs. Cette étude évalue l'influence de la densité de plantation et du semencier sur les croissances en diamètre et en hauteur des plants de P. elata, 20 mois après la plantation.

Méthode. Les graines ont été collectées sur 19 semenciers dans un peuplement naturel de 400 ha près de Kisangani (RDC). En octobre 2017, les plantules ont été transplantées dans une plantation suivant le dispositif de Nelder en trois répétitions, pour comparer la croissance de 540 plants en plein découvert suivant 10 densités différentes. Les accroissements en diamètre et en hauteur ont été analysés en fonction de la compétition locale et du semencier, à l'aide des modèles additifs généralisés.

Résultats. Vingt mois après la plantation, les croissances étaient maximales, en hauteur (165,6 ± 39,2 cm·an^-1) à une densité intermédiaire de 47 000 tiges·ha^-1, en diamètre (24,6 ± 8,3 mm·an^-1) à une faible densité de 4 200 tiges·ha^-1. L’origine génétique des plants a affecté leur croissance, les coefficients de variation entre semenciers ont atteint 7,6 % et 6,4 % pour les accroissements en diamètre et en hauteur, respectivement. Le diamètre initial a influencé positivement la croissance diamétrique mais la hauteur initiale n'a pas affecté la croissance en hauteur.

Conclusions. Le dispositif de Nelder s'est avéré utile pour confirmer la dépendance de la croissance des plants vis-à-vis de l’écartement et mettre en évidence l’influence de facteurs génétiques sur la croissance.

Abstract

Description of the subject. Pericopsis elata (Fabaceae) is a long-lived light-demanding tree from African rain forests that produces timber of high economic value. Natural populations suffer from overexploitation and a deficit of natural regeneration. Plantations could increase its production and limit the pressure on natural forests. However, we lack knowledge on the influence of spacing and seed tree origin on tree growth to optimize plantation protocols.

Objectives. This study evaluated the impact of sapling density and seed tree origin on the growth in diameter and height of P. elata plants, 20 months after plantation.

Method. Seeds were collected on 19 mother trees in a 400 ha natural stand near Kisangani (DRC). In October 2017, seedlings were transplanted in a Nelder design plantation with three replicates, to compare sapling growth at 10 contrasting plant densities under full sun. Diameter and height increments of 540 plants were analyzed according to local competition and maternal origin using generalized additive models.

Results. Twenty months after planting, height growth peaked at an intermediate density of 47,000 stems·ha^-1 (165.6 ± 39.2 cm·year^-1), while diameter growth peaked at a lower density of 4,200 stems·ha^-1 (24.6 ± 8.3 mm·year^-1). The maternal origin of saplings affected their growth, the coefficients of variation among mother trees reaching 7.6% and 6.4% for diameter and height increments, respectively. The initial diameter affected positively diameter growth while the initial height did not influence height growth.

Conclusions. The Nelder design proved useful to assess how plant growth depends both on spacing and genetic factors.

Table of content

1. Introduction

1The African tropical forest plays a major role in the world, through its multiple ecosystem services (Makana et al., 2011; Fétéké et al., 2015; Tieguhong et al., 2017). However, it is under increasing pressure from shifting agriculture and timber extraction (Bourbier et al., 2013), especially for industrial purposes. Timber is an important source of income for the Congo Basin countries (Toledo et al., 2011). However, its exploitation can have negative impacts on the natural regeneration of some species (Petrokofsky et al., 2015). There is a need to undertake in-depth studies on these commercially important species, to better understand their ecology (Gond et al., 2013), in order to develop sustainable management practices. Efforts to restore degraded forest landscapes are also emerging, such as the AFR100 initiative, seeking to restore one million km² in Africa (https://afr100.org).

2Effective management of production forests requires prior knowledge of species dynamics (Obiang et al., 2010; Vlam et al., 2014), in particular for the most sought-after species with high market value. For several decades, regulatory measures have been taken by the Congo Basin producer countries (Ouédraogo et al., 2014; Ligot et al., 2019) to ensure the sustainability of the resource.

3Some authors argue that the establishment of plantations in these countries would guarantee the continuity of timber production after logging, in particular for long-lived light-demanding species (Ngueguim et al., 2010; Petrokofsky et al., 2015). Pericopsis elata (Harms) Meeuwen, a legume tree traded under the name Afrormosia or African teak, is a good example: its demand is strong on the international markets (Pérez et al., 2005; Doucet et al., 2016) and its natural regeneration is problematic under dense canopy (Bourland et al., 2012a). Having been completely overexploited in West Africa at the end of the 20^th century, P. elata has been listed on CITES Appendix II Annotation #17 (logs, sawn wood, veneer sheets, plywood and transformed wood) and is considered "endangered A1cd" by the IUCN (https://www.iucnredlist.org/species/33191/9759606). Nowadays, the main producing countries are DRC, Cameroon and the Republic of the Congo (Bourland et al., 2012b).

4Natural populations of P. elata generally display Gaussian distribution of stem diameter, a characteristic of species with low natural regeneration (Bourland et al., 2012a; Ouédraogo et al., 2014; Vlam et al., 2014). This population structure does not guarantee the long-term availability of P. elata stems of exploitable diameter (Halpin & Lorimer, 2017).

5The introduction of the species in large-scale plantations may be a solution to compensate its overexploitation in Central Africa, at least if the wood conserves its good mechanical properties in plantations. Among various aspects, it requires to identify the optimal plantation density, hence to better understand how between tree interferences affect tree growth in both height and diameter. We expect that competition must impact growth. However, tree species usually react to light competition by investing more energy into height growth than diameter growth, which is interesting to obtain large stems (Brunner & Nigh, 2000), so that maximal height growth may occur at an intermediate density. Genetic factors can also affect tree growth, so that the choice of seed trees can also matter when planning plantations (Sotelo et al., 2008). However, this aspect has rarely been investigated in hardwood tropical tree species, especially in Africa (Koskela et al., 2014).

6The present study aims to determine the factors influencing the early growth rate of P. elata saplings under full light conditions and identify the density that optimizes the growth. More specifically, it addresses the following questions:

7– at which density and level of interference with neighbouring plants do P. elata saplings maximize their growth in stem diameter and height?

8– how much of the variability in growth is explained by mother tree, indicating a genetic determinism?

2. Materials and methods

2.1. Study location and environmental conditions

9The study took place in the Tshopo province, Democratic Republic of the Congo (DRC), near the city of Kisangani. Seeds were harvested in a 400 ha natural stand of P. elata located in Biaro forest (25°20' E, 0°11' N), 30 km southeast of Kisangani. The experimental design (plantation) is located on the outskirts of the city of Kisangani (25°15' 58'' E, 0°30' 42'' N, elevation 495 m), in the concession of a logging company (Compagnie Forestière et de Transformation, CFT).

10The climate of Kisangani is Af type according to the Köppen classification with an average annual temperature of 25 °C, a maximum temperature in February (32 °C) and an average annual rainfall of 1,672 mm (De Ridder et al., 2014). The driest month has a precipitation of more than 60 mm and the thermal amplitude is less than 5 °C (Dudu, 1991). The region is characterized by a long rainy season from August to December, interrupted by two dry seasons in January-February and June-July and the mean relative atmospheric humidity is between 90% and 95% (De Ridder et al., 2014).

11The Kisangani region is formed by recent superficial horizons and old ferralitic layers of soils rich in sand and clay (Mosango, 1983). The geological formation of the Biaro forest consists of sandstone, red clay, marl and limestone, with macro- and microporosity in the upper part of the soils (De Ridder et al., 2014). Fertility is low (cation exchange capacity – CEC – varies between 2 and 8 meq per 100 g) and the soil is acidic (pH varies between 3.5 and 5.5) (De Ridder et al., 2014).

2.2. Seed sampling and nursery

12Seeds were harvested between 17 September and 12 December 2016 in the 400 ha plot of Biaro forest, where all 196 P. elata trees (diameter at breast height ≥ 10 cm) had been inventoried. Seed trees were retained only if they produced mature seeds and were sufficiently isolated to avoid any confusion of seed provenance (19 seed trees eventually contributed to the seed pool used). One week after harvesting, seeds were sowed in polyethylene bags containing soil with natural manure, and an ID label indicating the seed tree. The nursery was located next to the botanical garden of the Science Faculty campus of the University of Kisangani, under mild shading. After germination, we randomized the positions of the seedlings in the nursery and followed their growth for one year before transplantation.

2.3. Transplantation following the Nelder design

13We installed three plantation trials in October 2017 according to the design Ia proposed by Nelder (1962), whereby plants are located along a set of concentric circles ensuring a gradient of sapling density (Figure 1). Each unit (Nelder 1, Nelder 2 and Nelder 3) is referred to as a replicate in the following sections. The experiment was carried out over an area of 135 x 45 = 6,075 m² (including areas around each Nelder design) and each replicate covers an area of 908 m². This area was cleaned from former vegetation (fallow land), allowing to prevent the destruction of the designs by uncontrolled bush fires.

14The replicates were placed next to each other to maximize the similarity of environmental conditions (edaphic, topographic and sun exposure). Each Nelder design consisted of 12 concentric circles, each one with 18 equally spaced saplings (Figure 1). The saplings located on the first and last circles formed the buffer zones, placed respectively at a distance of 0.6 and 15.6 m from the centre of the circles, so that a total of 180 saplings out of 216 were monitored for each replicate. The local spacing between neighbouring saplings ranged from 0.2 x 0.2 m to 3.8 x 3.8 m (Table 1).

15The theoretical area S available for each sapling is related to the density of each circle:

16with l the radius of the circle where the target sapling is located, m the radius of the largest neighbouring circle, k the radius of the smallest neighbouring circle and 18 the number of saplings per circle (Figure 2A).

17Saplings from the nursery were transplanted on October 2017, when they had reached heights ranging from 40 cm to 160 cm (median around 65 cm). As available sample sizes varied widely among maternal families, we could not conceive a well-balanced design. In Nelder 1 and Nelder 2, where 19 and 10 families were represented, we ensured that saplings from the same mother tree were never contiguous. This was not feasible in Nelder 3, represented by seven families but with two of them constituting 88% of the saplings. Within each Nelder design, we ensured that represented families were well spread among the different circles to avoid any confounding effect between maternal tree and local density. Initial sapling height did not differ significantly among replicates. Dead saplings were replaced immediately when found.

2.4. Data collection and analyses

18The total height and diameter at 10 cm from the collar were measured after two months and then at three months intervals on all saplings between December 2017 and June 2019, using a ruler and a calliper graduated in cm and mm, respectively. The total height was measured from the collar to the terminal bud of the longest branch. Weeding was carried out before each measurement period.

19To avoid a bias resulting from the stress experienced during transplantation, which could slow down initial growth, the growth rates in diameter and height were measured by the increments that occurred between the 5^th and 20^th month after transplantation, multiplied by 12/15 to obtain annual growth rates. Saplings replaced (2.8%) as a result of mortality caused by vandalism or planting stress were not included in the statistical growth analyses.

20The degree of interference between neighbouring saplings is expected to decay from circle C1 to C12 as the spacing between plants increases (Table 1) but it should also depend on the respective sizes of adjacent plants. To take this into account, we computed a competition index (CI_i) for each sapling i as a function of the distance and diameter of neighbouring saplings situated within a circle of radius K (Figure 2B), according to the relationship (Hegyi, 1974; Vanclay et al., 2013):

21where D_i and D_j are the respective diameters of the sapling subject i and competitor j and L_ij is the distance between i and j (L_ij ≤ K). The competition radius K, which is the maximal distance at which we consider that saplings can interfere, was calculated as the mean height of all the saplings of a given replicate, reaching 2.7 m, 2.4 m and 1.9 m, for Nelder 1, 2 and 3, respectively (Uhl et al., 2015; Figure 2B).

22The medians of the diameter and height annual growth rates were compared between replicates using Kruskal-Wallis tests. The coefficient of variation (CV) was used to calculate the variability between circles. We modelled growth rates with generalized additive models (GAM) of the Gaussian type (Zuur et al., 2009). These models combine linear and non-linear explanatory variables influencing the response variable (Uhl et al., 2015).

23The models used are of the form:

24where ΔH_qpi and ΔD_qpi are the response variables, respectively the annual increments in height and diameter in a Nelder design q for sapling i from a mother tree p, a_q is the effect of the design (fixed effect), b the slope describing the impact of the initial height (InH_i) or initial diameter (InD_i) of the sapling (at the 5^th month after transplantation), f(CI_i) is a nonparametric smoothing function describing the influence of the competition index on sapling i, α_p is the mother tree effect (random effect following a centred normal distribution of standard deviation d) and ε_i is the residual error of the model (centred normal distribution of standard deviation σ). We estimated the CV expressing the relative effect of the mother trees by the ratio of d over the mean growth rate, where d is estimated by the GAM. To assess if the individual-based competition index (CI_i) better explains growth rates than a more simple model accounting only from the local sapling density, we also ran models similar to (3) and (4) but replacing f(CI_i) by f(CN_i), where CN_i denotes the circle number (1 to 12) where sapling i is located.

25Mother trees represented by < 10 individuals (seven out of 19 mother trees) were not included in the models, removing only 4% of the saplings. To find the optimal model, an iterative approach consisting of adding or removing each explanatory variable in turn was used. The best models were selected on the basis of the Akaike Information Criterion (AIC). The 'mgcv' package version 1.8-28 (Wood, 2011) for the statistical software R (R Core Team, 2018) was used for the GAM. The distances between saplings were computed with the GeoGebra Classique software.

3. Results

3.1. Height and diameter growth trajectories

26Sapling height and diameter varied with time and densities and showed weak growth rates during the first five months after transplantation (Figure 3). Growth rates seemed fairly constant between the 8^th and 17^th month, but tended to be lower between March and June, both in 2018 and 2019 (Figure 3). At 20 months, sapling height ranged from 24 to 451 cm, with a mean and standard deviation of 152.7 ± 82.8 cm (CV = 54.2%), while the diameter ranged from 3 to 72 mm, with a mean and standard deviation of 19.2 ± 11.1 mm (CV = 57.6%).

3.2. Competition Index (CI)

27The competition index averaged 37.6 m^-1 with a standard deviation of 49.1 m^-1. Differences between replicates were not significant (Kruskal-Wallis; X² = 2.14; p = 0.34). The CI was highest on circle C1 (128.7 ± 63.9 m^-1) and decreased accordingly to circle density, reaching zero for the last circle C10 (Figure 4A). The competition index varied also significantly within each circle (except for C10) with a CV ranging between 21% and 50%.

3.3. Annual height increment

28The overall annual height increment was 136.5 ± 48.7 cm⋅year^-1. It was largest in Nelder 1 (152.9 ± 48.1 cm⋅year^-1) followed by Nelder 2 (132.4 ± 49.0 cm⋅year^-1) and then Nelder 3 (124.2 ± 44.4 cm⋅year^-1). These differences were significant (Kruskal-Wallis; X² = 35.35; p = 2.11·10^-8). It increased from circle C1 to C3 or C4 and then decreased towards C10 (Figure 4B). Highest increments occurred on circles C3 (165.6 ± 39.2 cm⋅year^-1), C4 (164.2 ± 42.0 cm⋅year^-1) and C5 (156.6 ± 39.8 cm⋅year^-1). Variability was high within all Nelder circles (CV between 19 and 47%).

29The GAM explained 40% of the deviance in sapling height increment, and leads to residuals showing satisfactory distribution (Appendix Figure A1). Annual height increment significantly depended on competition (p < 2⋅10^-16), showing a non-linear and non-monotonic relationship (Table 2). It peaked when the competition index was close to 50 m^-1 (Figure 5A). Height increment was also significantly influenced by sapling mother trees (p = 7.56⋅10^-16), with a standard deviation of annual increment of 10.4 cm⋅year^-1 (Table 2). Hence, reported to the mean annual increment, the variation due to mother trees reached a coefficient of variation of CV = 10.4 / 136.5 = 7.6%. The GAM confirmed that height increment significantly differed among Nelder replicates (Table 2), while the initial height had no significant impact and was removed from the final model. When we used the circle number instead of the local competition index to run the model, the explained variance dropped from R² = 40.2% to 33.0% and the AIC increased by 53 units (5,276 to 5,329).

3.4. Annual diameter increment

30The overall annual diameter increment was 18.4 ± 7.7 mm⋅year^-1. The largest mean was observed in Nelder 1 (21.4 ± 9.1 mm⋅year^-1) followed by Nelder 2 (17.8 ± 7.6 mm⋅year^-1) and then Nelder 3 (16.0 ± 4.8 mm⋅year^-1). These differences were significant (Kruskal-Wallis; X² = 43.21; p = 4.15⋅10^-10). The increments tended to increase from circle C1 to C7 and then to decrease until C10 (Figure 4C). High averages were located on C5 (20.2 ± 5.8 mm⋅year^-1), C6 (20.5 ± 6.7 mm⋅year^-1), C7 (24.6 ± 8.3 mm⋅year^-1) and C8 (20.5 ± 6.7 mm⋅year^-1). Variability was high within all circles (CV between 22 and 50%).

31The GAM explained 56% of the deviance in annual diameter increment (residuals shown in Appendix Figure A2). Diameter increment decreased non-linearly but monotically with the competition index (p < 2⋅10^-16; Table 2; Figure 5B). It also depended on the initial diameter (p = 1.38⋅10^-12) and varied significantly among Nelder designs (Table 2). Finally, it was also significantly influenced by the mother trees of the saplings (p = 6.73⋅10^-5; estimated standard deviation among maternal families reaching 1.2 mm⋅year^-1, CV = 6.4%; Table 2). Using the circle number instead of the local competition index, the explained variance dropped from R² = 56.1% to 50.5% and the AIC increased by 64 units (3,197 to 3,261).

4. Discussion

32Our results highlight that the density at planting (competition) and the genetic origin (i.e. mother tree) of the saplings influence their growth. For both height and diameter growth, our results show that the competition index (CI), which accounts of the relative heights of the neighbouring plants, better explains the growth rate than the circle number, which reflects only the local density of saplings. The impact of density or competition on growth rates is far from linear and shows a mode in the case of height growth (Figures 4 and 5), justifying the use of a smoothing function (GAM) to describe it. Hereafter, we discuss the different effects observed.

4.1. Determinants of saplings height growth

33Density seems at this stage of the evolution of the plantation to stimulate height growth of P. elata saplings, with an optimum corresponding to a spacing of c. 0.5 – 0.6 m between saplings (Figure 4; Table 1). Similar observations were reported by Fayolle et al. (2015) when comparing the performance of non-competitive and competitive stems.

34The mean height increment (135.3 ± 50.1 cm·year^-1) was greater than that found by Ouédraogo et al. (2014) in Cameroon (80.32 ± 36.05 cm·year^-1) in logging gaps, probably due to higher sun exposure. At the densities with the highest increments, there was a high variability in the competition index.

35The competition for access to light would have favoured the scrolling of saplings (Brunner & Nigh, 2000; Ngueguim et al., 2012; Parrott et al., 2012; Kuehne et al., 2013; Fayolle et al., 2015), explaining why size is negatively correlated with the competition index until some density threshold (Vanclay et al., 2013; Uhl et al., 2015). Above 46,000 stems·ha^-1, competition no longer favours height growth but impedes it. This is in line with the observations of Gourlet-Fleury (1998) showing that dominated trees do not develop their crowns sufficiently, with the result that their growth is limited. It is likely that the optimal density maximizing height growth rate will move toward lower densities as trees age and increasingly compete for resources. Therefore, the optimal density at 20 months should not be considered as an optimal operational density for P. elata plantations.

36The growth rate was larger in Nelder 1 than Nelder 2 or 3. Whether subtle differences of soil fertility between the three Nelder plots might explain these differences remains to be investigated. By contrast, the initial height did not influence the growth of P. elata saplings.

4.2. Determinants of saplings diameter growth

37The mean diameter increment (18.4 ± 7.7 mm·year^-1) found in this study is larger than the values reported in Cameroonian logging gaps under silvicultural treatments (clearings) at 60 months (9.1 ± 4.6 mm·year^-1 in Fayolle et al., 2015; 8.0 ± 4.0 mm·year^-1 in Ouédraogo et al., 2014). The latter studies occurred in the concession of Pallisco Company, certified FSC (Forest Stewardship Council), which practices controlled felling. Under these conditions, the average size of the gaps is relatively small, which severely limits the amount of sunshine for the saplings (including gaps that have been completely cleared). Hence, difference of light exposure might explain these differences in diameter growth.

38As for height growth, diameter growth varied among Nelder plots, was significantly influenced by the mother tree, and strongly responded to sapling density (competition). However, unlike height growth, diameter growth was positively affected by the initial diameter of the saplings after transplantation, and it increased monotonically as competition decreased. Thus, saplings grow faster in diameter at lower densities (Lemoine & Sartolou, 1980; Kuehne et al., 2013; Uhl et al., 2015). In addition, local environmental conditions influenced by the transport and use of water and nutrients by saplings, are likely to influence their development (Beckage & Clark, 2003; Rio et al., 2014; Colombaroli et al., 2016; Kearsley et al., 2017). In future analyses, it may be appropriate to integrate the environmental changes that occur as the seedlings grow, for example due to the faster accumulation of litter in the centre of repetitions, compared to their extremities.

4.3. Influence of mother tree on growth

39Our results corroborate the few studies conducted on tropical hardwood species showing that a significant part of the variability in tree growth can be explained by genetic factors (e.g., León et al., 2017; Luechanimitchit et al., 2017; Sawitri et al., 2020), a prerequisite for tree breeding. Kuehne et al. (2013) found that the influence of the mother tree on growth is most noticeable during the juvenile stages. This is the well-known "early selection" principle applied by forest geneticists right from the nursery (Nanson, 2004). We used open-pollinated seeds but as P. elata is able to self-fertilize (D.-M. Assumani & O.J. Hardy, unpublished), we could not estimate the heritability of growth rates by simply assuming that we had half-sib families. Saplings will need to be genotyped using genetic markers to assess their actual relatedness and to estimate the heritability of measured traits. Among the genetic causes influencing growth, inbreeding depression may also be an important factor in species endogamous like P. elata. It results from the expression of deleterious recessive genes in homozygotes, diminishing the fitness of the more consanguineous individuals (Hedrick et al., 2015). The vast majority of tropical trees have developed mechanisms preventing or strongly limiting inbreeding (e.g., self-incompatibility system) but P. elata is an exception (Micheneau et al., 2011). Hence, by identifying self-fertilized saplings using genetic markers, our plantation trials will also allow us to assess whether inbreeding depression is a feature to take into account.

4.4. Study improvement options

40Knowledge of the growth dynamics of saplings can lead to better management of stands rich in P. elata stems, in plantations or in natural forests where P. elata is very gregarious. Our results are based on seedlings aged about 30 months (7 to 10 months in the nursery plus 20 months of planting) and we need to wait a few more years before being able to make robust recommendations regarding optimal tree spacing for plantations. Moreover, studying the influence of environmental conditions (climatic, topographic, etc.) for the development of P. elata could be relevant (Addo-Danso et al., 2012). Nevertheless, our preliminary results show the interest of trials based on Nelder design, which are also appropriate to assess the performances of different seed families.

5. Conclusions

41Our study of the growth dynamics of P. elata saplings allows us to assess the optimal plantation density (i.e. best increments in both height and diameter). Density had a non-linear relationship with diameter and height growth. As expected, the sapling density that maximizes the increment in height is larger than the one maximizing the increment in diameter. At 20 months, the best compromise is found at the C5 arc corresponding to some 14,079 stems·ha^-1, leading to mean increments of 157 ± 40 cm⋅year^-1 in height and 20 ± 6 mm⋅year^-1 in diameter. Pericopsis elata could therefore be planted at a spacing of 0.8 x 0.8 m to promote sapling growth in the first 20 months. However, this optimal density will need to be re-evaluated on longer term. The genetic origin of the saplings influenced the increments in height and diameter, indicating that selection of plus seed trees should be efficient. On longer term, the design will make it possible to assess the impact of inbreeding on the growth of saplings. Our plantation trials based on Nelder design prove interesting to assess the effects of both the density and the seed tree on tree growth. Such design, requiring a limited surface, could be replicated in other areas and for other commercial species that may be of interest for restoring forest landscapes.

Acknowledgements

42We are grateful to the “FOrmation, Recherche, Environnement dans la TShopo” (FORETS) project for financially supporting the first author during his studies. This project, funded with resources from the XI European Development Fund, is led by CIFOR in partnership with the company Resources and Synergies Development (R&SD) and the University of Kisangani. Our thanks also go to the Belgian Directorate-general Development Cooperation and Humanitarian Aid funded project N3-GeForCo led by the Royal Museum for Central Africa, at the origin of the experiment, as well as the ARES “Projets de synergie-2020/Régénération assistée de Pericopsis elata (Harms) Meeuwen : étude par une approche de traits fonctionnels foliaires” grant. We thank the company R&SD DRC, through its Manager G. Assimbo, for logistical assistance and the company “Compagnie Forestière et de Transformation” (CFT) through its Director K. Ammacha, for hosting the tree nursery and the Nelder wheels, providing unlimited access to the site and other logistical assistance. We are grateful to C. Lubwilu Lolo, P. Mayele Maloko, E. Ngombo Luzolo (CFT) and E. Tudico and his team (R&SD Construction sia) for technical assistance on the field. O.H. is supported by a Belgian Fund for Scientific Research, F.R.S-FNRS.

43N.B. and H.B. conceived and implemented the Nelder designs. D-M.A. collected the seeds in the Biaro forest, set up and monitored the saplings in the nursery under N.B. and O.H.’s supervision. H.H. and N.B. led plantation works. H.H. and C.I. collected data on the field and at the laboratory. C.I., O.H. and G.L. performed the statistical analyses. C.I. wrote the manuscript with the support of N.B., O.H., G.L., L.N. and H.B. The final version of the manuscript was read and corrected by all the authors.

Bibliographie

To cite this article

About: Crispin Ilunga-Mulala

About: Hulda Hatakiwe

About: Hans Beeckman

About: Olivier J. Hardy

About: Gauthier Ligot

About: Dieu-Merci Assumani

About: Léopold Ndjele

About: Nils Bourland