# The Effect of Four Abiotic Factors on Macro-Anatomical Markers Development in Parkia biglobosa, Jack, R. Br., 1830 (Fabaceae) Crown

Beda Innocent Adji, Doffou Sélastique Akaffou, Kouadio Henri Kouassi, Yao Patrice Houphouët, Jerôme Duminil, Sylvie-Annabel Sabatier

► To cite this version:

Beda Innocent Adji, Doffou Sélastique Akaffou, Kouadio Henri Kouassi, Yao Patrice Houphouët, Jerôme Duminil, et al.. The Effect of Four Abiotic Factors on Macro-Anatomical Markers Development in Parkia biglobosa, Jack, R. Br., 1830 (Fabaceae) Crown. American Journal of Plant Sciences, 2021, 12 (04), pp.645-661. 10.4236/ajps.2021.124044. hal-03211834

HAL Id: hal-03211834

https://hal.inrae.fr/hal-03211834

Submitted on 29 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The Effect of Four Abiotic Factors on Macro-Anatomical Markers Development in Parkia biglobosa, Jack, R. Br., 1830 (Fabaceae) Crown

Beda Innocent Adji<fnref n="1" />,<fnref n="2" />,<fnref n="3" />*, Doffou Sélastique Akaffou<fnref n="1" />, Kouadio Henri Kouassi<fnref n="1" />, Yao Patrice Houphouët<fnref n="1" />, Jerôme Duminil<fnref n="4" />,<fnref n="5" />, Sylvie Sabatier<fnref n="2" />,<fnref n="3" />

**How to cite this paper:** Adji, B.I., Akaffou, D.S., Kouassi, K.H., Houphouët, Y.P., Duminil, J. and Sabatier, S. (2021) The Effect of Four Abiotic Factors on Macro-Anatomical Markers Development in Parkia biglobosa, Jack, R. Br., 1830 (Fabaceae) Crown. *American Journal of Plant Sciences*, **12**, 645-661.

https://doi.org/10.4236/ajps.2021.124044

Received: March 6, 2021

Accepted: April 26, 2021

Published: April 29, 2021

Copyright © 2021 by author(s) and
Scientific Research Publishing Inc.
This work is licensed under the Creative
Commons Attribution International
License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

<figure><img src="image_2.png" /></figure>

Open Access

## Abstract

Interpretation of primary growth markers (modules) is neglected in sustainable resource management processes, yet it opens up prospects for long time series on tree crown development, necessary for their characterization in the current context of climate change. This study aimed to assess the morphological variation of crown shoots in *Parkia biglobosa* in the face of a changing environment. Axis growth modules of 420 individuals of different ages were retrospectively analyzed in the presence or absence of shading during the wet and dry seasons in seven localities in Côte d'Ivoire. The results showed that the developmental environment of individuals did not influence the growth modules dimensions (P > 0.05). However, module size remained significantly different between locations (P < 0.05). The southern part of the gradient is still favourable and has priority for the establishment of permanent plots. The modules morphology differs from the youngest to the oldest individuals (P < 0.05). The rainy season remains the ideal period for the implementation of agroforestry reforestation programmes based on this species (P < 0.05). *Parkia biglobosa* is a monochasial sympod (Pseudo-monopod) with a relay axis that follows the Paul Champagnat architectural model in the young stage, and transits to the Wilhelm Troll model later in the adult and old stage. This information contributes to the understanding of the functioning of crown and the adaptation of this species to a varying environment. It could guide choice of suitable environment and ideal genotype for the implementation of a reforestation or agroforestry programme based on *Parkia biglobosa*.

<footnote n="1">¹Université Jean Lorougnon Guédé, Agroforestry UFR, Daloa, Côte d'Ivoire</footnote>

<footnote n="2">²CIRAD, UMR AMAP, Montpellier, France</footnote>

<footnote n="3">³AMAP, Université Montpellier, CIRAD, CNRS, INRAE, IRD, Montpellier, France</footnote>

<footnote n="4">⁴IRD, UMR-DIADE, Montpellier Cedex 5, France</footnote>

<footnote n="5">⁵DIADE, Université Montpellier, CIRAD, IRD, Montpellier, France</footnote>

<footnote>Email: *adjibedainnocent@gmail.com, *beda-innocent.adji@cirad.fr</footnote>
Keywords

*Parkia biglobosa*, Environment, Development, Growth Modules, Crown

**1. Introduction**

Architectural and retrospective analysis is very important for the monitoring and
sustainable management of trees [1] [2]. Architectural analysis allows us to in-
terpret the overall structure of the tree and to understand the morphological
mechanisms that gave birth to it; retrospective analysis allows us to trace the past
growth of trees via morphological markers left by growth stoppages due to the
dry seasons. In fact, the latter provides access to long-term growth series for un-
derstanding and interpreting the structural functioning of tree crowns, the choice
of resilient and adapted ideotypes, and the choice of environment conducive to
the establishment of permanent plots necessary for their sustainable management
[3] [4].

Several studies of sustainable management of plant genetic resources via arc-
hitectural and retrospective analysis have proven effective in Europe and Latin
America [5] [6] [7] [8]. These methods have not yet been applied to species in
Africa; indeed, studies on architectural and retrospective analysis for the sus-
tainable management of many overexploited species in Africa are almost absent.
We believe that these methods should be applied in Africa, specifically in West
Africa. This would allow a better understanding of the coordination over time of
growth processes (elongation, thickening) of shoots within the crown and the
associated morpho-anatomy essential for long-term retrospective analysis of sus-
tainable management, tree development and any analysis of the causes of ob-
served variations [1]. In view of the need to characterise plant responses to cur-
rent climate change, such retrospective analysis will be welcome as a comple-
ment to data from permanent plots.

Our study focused on a well-known tree from sub-Saharan Africa, *Parkia big-
lobosa* (Fabaceae), for several reasons: 1) it is an emerging tree still in the wild
and widely distributed in this area [9]; 2) it is a beloved and overexploited spe-
cies for its food, medicinal and socio-cultural function [10] [11]; 3) it is an agro-
forestry species with soil restorative properties (Leguminosae) [12]; 4) the genus
*Parkia* has anatomical features suitable for counting annual rings [13]; 5) this
specie present unusual morphological features (modules) suitable for determin-
ing growth modules or units, which makes it particularly interesting for a re-
trospective study of savanna tree shoot growth.

This study aimed to characterise the development of growth modules in the
crown of *Parkia biglobosa* in the face of a changing environment in Côte d'Ivoire
through retrospective analysis. This approach would allow us to understand: 1)
the variability of shoot size in relation to different habitats; 2) the adaptability of
the species to various environments; 3) to identify genotypes adapted and resis-
tant to climatic variations and 4) to guide the choice of priority and favourable environments for the establishment of permanent plots (species domestication) within the framework of possible “Néré” based agroforestry programmes in Côte d'Ivoire.

## 2. Material and Methods

### 2.1. Plant Material

The plant material is made up of axis of 420 individuals of different ages chosen randomly and arbitrarily in natural stands according to two types of habitats. The criteria for sorting the age categories, the dendrometric characteristics and the number of individuals per age category and per habitat type as well as the age types examined are recorded in Table 1.

### 2.2. Methods

#### 2.2.1. Study Sites

Study was carried out at seven locations in Côte d'Ivoire. The study sites are arranged along a South-North bioclimatic gradient. The localities were chosen according to the specific vegetation to the *Parkia biglobosa* species (Figure 1). The characteristics of the surveyed locations are recorded in Table 2.

#### 2.2.2. Observation Method

##### 1) Individuals selection and sorting by age category

A total of 420 individuals from natural regeneration were sampled and sorted by age category (young, adult and old). The age of the individuals was arbitrarily chosen on the basis of the size (height and diameter) of individuals raised in nurseries and of individuals present on old plots set up by the CNRA (*Centre National de Recherche Agronomique*) of Côte d'Ivoire in the 1980s in the Korhogo

<table><caption>Table 1. Age selection criteria, dendrometric characteristics and number of individuals used in relation to age and habitat in *Parkia biglobosa*.</caption><thead><tr><th rowspan="2">Ages</th><th rowspan="2">Habitats/environments</th><th colspan="3">H (m)</th><th colspan="3">D or DBH (cm)</th><th rowspan="2">Indiv-Nbr</th><th rowspan="2">Axis examined</th></tr><tr><th>Mini</th><th>Maxi</th><th>Mea</th><th>Mini</th><th>Maxi</th><th>Mea</th></tr></thead><tbody><tr><td>Young tree<br>(1 - 5 Yr)</td><td>Full sunlight</td><td>0.53</td><td>5.5</td><td>1.95 ± 0.11</td><td>0.56</td><td>18.96</td><td>4.35 ± 0.24</td><td>70</td><td>Main stem</td></tr><tr><td>H (m): 0.1 - 6;<br>D (cm): 1 - 19</td><td>Undergrowth/shading</td><td>0.42</td><td>5.8</td><td>2.56 ± 0.05</td><td>1.26</td><td>17.83</td><td>7.06 ± 1.22</td><td>70</td><td>Main stem</td></tr><tr><td>Adult tree<br>(6 - 20 Yr)</td><td>Full sunlight</td><td>7.5</td><td>11</td><td>9.68 ± 1.14</td><td>20.06</td><td>49.68</td><td>38.19 ± 3.02</td><td>70</td><td>Branch</td></tr><tr><td>H (m): 7 - 14;<br>D (cm): 20 - 50</td><td>Undergrowth/shading</td><td>8</td><td>13.7</td><td>12.02 ± 1.01</td><td>22.4</td><td>48.8</td><td>36.27 ± 1.53</td><td>70</td><td>Branch</td></tr><tr><td>Old tree<br>(21 Yr ≤)</td><td>Full sunlight</td><td>18</td><td>31.7</td><td>25.22 ± 3.21</td><td>61.3</td><td>178.44</td><td>85.51 ± 6.12</td><td>70</td><td>Branch</td></tr><tr><td>H (m): 15 ≤;<br>D (cm): 60≤</td><td>Undergrowth/shading</td><td>20.5</td><td>29.5</td><td>25.3 ± 3.32</td><td>60.1</td><td>115.61</td><td>76.98 ± 7.36</td><td>70</td><td>Branch</td></tr></tbody></table>

Yr = Year old, H = Height, D or DBH = Diameter or diameter at breast height, Mini = minimum, Maxi = Maximum, Moy = Mean, m = meter, cm = centimeter, Indiv-Nbr = number of individuals assessed.
<table><caption>Table 2. Characteristics of the surveyed localities [14]-[28].</caption><thead><tr><th>Sites/locations</th><th>GPS coordinates</th><th>Vegetation</th><th>Climate</th><th>Temperature (°C)</th><th>Rainfall (mm/year)</th><th>Soil type</th></tr></thead><tbody><tr><td>Ferké</td><td>5°23'43.39644"W; 9°36'1.87056"N</td><td>Grassy and little wooded savannah</td><td>Dry tropical</td><td>27 - 40</td><td>263 - 1200</td><td>Ferralitics soil (Ferrisols, Cambisols, Fluvisols, Luvisols), highly to moderately desaturated.</td></tr><tr><td>Korhogo</td><td>5°36'12.39612"W; 9°33'24.68988"N</td><td>Open forest (wooded savannah)</td><td>Dry tropical</td><td>26.6 - 35.7</td><td>817 - 1216</td><td>Ferruginous (90%) and Ferralitic (10%): superficial gravelly soil, deep gravel with a heavy texture, low in organic matter, highly desaturated.</td></tr><tr><td>Niakara</td><td>5°18'40.73544"W; 8°40'47.97912"N</td><td>Wooded and grassy savannah</td><td>Dry tropical</td><td>24.7 - 38</td><td>800 - 1230</td><td>Complex of slightly desaturated ferrallitic soils and eutrophic brown tropical soils derived from basic rocks.</td></tr><tr><td>Katiola</td><td>5°7'35.814"W; 8°13'53.94"N</td><td>Wooded and grassy savannah</td><td>Dry tropical</td><td>24 - 36</td><td>1100 - 1200</td><td>Moderately and highly desaturated ferrallitic soils</td></tr><tr><td>Bouaké</td><td>5°5'47.3289"W; 7°40'45.335"N</td><td>Clear forest (wooded savannah)</td><td>Wet tropical</td><td>23.6 - 34</td><td>1100 - 1200</td><td>Gravelly, moderately saturated, reworked, shallow ferralitic gravel from a granitic alteration material with a sandy-clay texture.</td></tr><tr><td>Toumodi</td><td>5°1'34.95576"W; 6°22'42.67848"N</td><td>Open forest (wooded savannah, grassland and gallery forests)</td><td>Wet tropical</td><td>26.6 - 30</td><td>1092 - 1200</td><td>Ferralitic soil on granitic bedrock (sandy-clayey soil), characterised by the weak differentiation and friable consistency of their horizons.</td></tr><tr><td>Daloa</td><td>6°26'9.19788"W; 6°54'32.058"N</td><td>Tropical rainforest</td><td>Wet tropical</td><td>21 - 34</td><td>1000 - 1900</td><td>Ferralitic, deep, acidic and desaturated in exchangeable bases, rich in organic matter.</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>Figure 1. Location of the study area.</figcaption></figure>

<footnote>*C = Celsius degree, mm = millimeter, W = west, N = North.</footnote>
locality. The characteristics of the age groups as well as the criteria for sorting individuals by age category are recorded in Table 1. 10 individuals were identified by chance and evaluated by habitat and age in each locality surveyed in the rainy season (optimal condition) and in the dry season (restrictive condition). We therefore have for each age category (Young, Adult and Old), 10 individuals identified and evaluated by habitat type (Undergrowth/Darkness and Full Sun/Sun) in each of the 7 localities above (10 × 3 × 2 × 7 = 420 individuals) (Table 1).

## 2) Choice of axis and habitat type

Observations were made on different types of axes depending on the accessibility of the crowns in two different habitats. For young trees, the axes assessed were the main trunks because of the easy access and the non-frequency of secondary branches on all individuals in this category. For mature and old trees, the axes assessed were only tertiary branches and short twigs due to accessibility. Assessments were carried out in situ for young trees; whereas for mature and old trees, branches were cut and transported to the laboratory for observations. Two types of environments were considered: undergrowth and full sun. The former refers to individuals living in a very shaded environment with a forest canopy or in an overcrowded environment or with superior shelter. The latter refers to individuals isolated in full sun or in an open environment in direct contact with sunlight.

## 3) Identification of primary markers (Modules) of growth by analysis

This analysis is made possible by the recognition of morphological markers that result from the functioning of meristems and that remain for several years. These morphological markers are the scars left by bud scales, cataphylls and deciduous leaves, naturally pruned twigs, a bend in the axis at the site of apical necrosis or trauma, the texture of the bark, etc. The growth unit or elongation unit or morphological unit or extension unit (GU) [6] [8] [29] [30], corresponding to the set of organs (phytomeres) initiated by the apical meristem during an uninterrupted phase of organogenesis and elongation is delimited at its base by a large number of scars of pairs of scales (7 to 16) left by the terminal bud. From this characteristic, it is possible to identify the successive growth units arranged along a leafy stem.

In general, two types of shoots are generated by the establishment of growth units by the terminal bud: monocyclic shoots (establishment of a single growth unit during a year or spring growth wave) and polycyclic shoots (establishment of several growth units during a year or spring and summer growth waves) [3] [6] [32]. If the first GU (spring growth wave) can be qualified as preformed (in the sense of [33]: presence of all the leaf organs of the future GU in the bud before its budburst), some organs in the GU are qualified as neoformed (prolongation of the setting of the leaf organs by the apical meristem in the same time after the preformed organs). In the tropics, the majority of shoots are polycyclic due to a rhythmicity adopted by the growth of species in this zone [33] [34]. This is the case for *Parkia biglobosa*.

In *Parkia biglobosa*, the growth units (GU) or macro-anatomical markers cor-
respond to modules (MOD), due to a defined and limited growth of the axes
generally caused by an apical trauma (death of apex), the axillary bud close to
the apex thus takes over, until it also dies, and so on. The structure of the axis
thus has a twisted appearance (Figure 2); we speak of a modular plant because of
a succession of structures or several superposed modules [6] [34]. Indeed, the
axis of these types of plants is a branched structure resulting from the function-
ing of several meristems and, depending on the number of branches formed, we
speak of a monochasial (one relay axis), dichasial (two relay axes) or polychasial
(several relay axes) sympod.

4) Parameters evaluated

The retrospective analysis focused on the first four growth modules (MODs) from the apical apex of the sampled axis. Morphological parameters assessed on each axis concerned the total number of modules carried by the observed axis, the length of each module, the diameter at the base of each module and the number of leaves or phytomers of the first four modules from the tip (apical start) of the sampled axes [1] [4] [35] [36]. The harvested stems and branches were split (longitudinal and transverse cuts) to identify growth stops (portions delimiting modules) imperceptible on the bark of the axes. Data were collected during the rainy season and this action was repeated on the same individuals during the dry season.

**2.2.3. Statistical Data Analysis**

Statistical analyses were first performed using one-dimensional descriptive
statistics and linkage analysis methods using XLSTAT 2020 version 7.5. The

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Images related to morphology and growth boundaries or growth stops (GS) of macro-anatomical markers or Growth modules (MOD) in *Parkia biglobosa*.</figcaption></figure>
difference between the evaluated parameters was carried out by means of one-
way analyses of variance (ANOVA) with SAS software version 9.4. The Student-
Newman-Keuls test at the 5% threshold was used for post hoc comparisons.

**3. Results**

**3.1. Architectural Development of the Individuals Observed**

According to the observations, the architectural development of the species is
sympodial monochasial with only delayed relaying (delayed growth and branch-
ing). The terminal bud dies after a certain period of elongation, the axillary bud
closest to the terminal bud takes over and so on. The plant appears as a pseu-
do-monopod with a twisted trunk structure. *Parkia biglobosa* can be classified in
the architectural model of Paul Champagnat (1921-2011) and also of Wilhelm
Troll (1897-1978) because of the way it builds its structure. Indeed, according to
our remarks, the development of the species is a repeated complex of unbranched
and discretely branched axes with short branches. The axes are vertical at their
base and their distal end is by large radius curvature, reduced to horizontal or
collapsed due to gravity and the flexibility of the wood. The leaves are spirally
inserted and the flowering is sometimes lateral, sometimes terminal. In the young
stage, the plant adopts the Champagnat pattern, while in adulthood and old age,
the Troll pattern is more expressive.

**3.2. Morphological Characteristics of Growth Modules**

Globally, the number of growth modules carried by the observed axes varied from 2 to 26 with an average of 10.08 ± 4.24; the average length of the Modules varied from 3.85 to 51.40 cm with an average of 18.57 ± 8.15 cm; the average diameter of the Modules varied from 0.43 to 8.81 cm with an average of 1.43 ± 0.82 cm. The average number of phytomers counted per Module varied from 5.5 to 33.5 with an average of 15.09 ± 5.35. The morphology and limits of the growth modules are presented in Figure 2.

**Figure 3** shows the average distribution of the number and length of phyto-
mers per growth module (MOD) as well as the growth stop or growth boundary
zones (GS) along a *Parkia biglobosa* axis.

**3.3. Influence of the Habitat on Growth Modules Morphology**

Analysis of the variance of modules morphological dimensions present on the
observed axes according to the two environments considered (Table 3), showed
that the morphology of the modules of the axes exposed to sunlight is statistical-
ly identical to the morphology of the modules of the axes benefiting from a su-
perior shelter (P > 0.05). The environment in which this species evolves does not
seem to have an effect on the dimensions of the modules. In *Parkia biglobosa*,
the type of habitat (bright or dark light) does not significantly influence the mor-
phology of the growth modules (P > 0.05).
<table><caption>**Table 3.** Comparison of modules morphological parameters according to the habitat type.</caption><thead><tr><th>Habitat</th><th>Tree-H (m)</th><th>Tree-Diam (cm)</th><th>Axis-Nbr-Mo</th><th>Mo-leng (cm)</th><th>Mo-Diam (cm)</th><th>Mo-Phyto-Nbr</th></tr></thead><tbody><tr><td><strong>Full sun</strong></td><td>10.41 ± 1.23 a</td><td>32.72 ± 3.46 a</td><td>9.8 ± 0.46 a</td><td>18.73 ± 0.88 a</td><td>1.43 ± 0.09 a</td><td>14.61 ± 0.53 a</td></tr><tr><td>Min-Max</td><td>0.53-31.7</td><td>0.56-178.44</td><td>2-26</td><td>3.85-51.4</td><td>0.43-8.81</td><td>5.5-30.5</td></tr><tr><td><strong>Shadowing</strong></td><td>10.43 ± 1.76 a</td><td>29.72 ± 3.45 a</td><td>10.64 ± 0.55 a</td><td>18.25 ± 1.08 a</td><td>1.42 ± 0.07 a</td><td>16.05 ± 0.79 a</td></tr><tr><td>Min-Max</td><td>0.42 - 28.3</td><td>1.26 - 113.44</td><td>5 - 19</td><td>6.25 - 41.5</td><td>0.5 - 2.77</td><td>6.5 - 33.5</td></tr><tr><td><strong>Pr > F</strong></td><td><strong>0.9771</strong></td><td><strong>0.5213</strong></td><td><strong>0.2685</strong></td><td><strong>0.7440</strong></td><td><strong>0.9280</strong></td><td><strong>0.1277</strong></td></tr></tbody></table>

Values with the same letters are not statistically different at the 5% level. Tree-H (m) = tree height in meters; Tree-Diam (cm) = diameter at breast or crown height of the tree in centimetres; Axis-Nbr-Mo = number of modules or growth units carried by the measured axis; Mo-leng (cm) = average length of modules or growth units in centimeters; Mo-Diam (cm) = average diameter of modules or growth units in centimeters; Mo-Phyto-Nbr = average number of phytomers carried by the modules or growth units; Min = minimum et Max = maximum.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Modelling of growth, rank and average distribution of the number of phytomeres by growth module along the observed axes in *Parkia biglobosa*. MOD/GU = Module or growth unit, GS = Growth stop or module limits, cm = Centimeter.</figcaption></figure>

### 3.4. Influence of the Surveyed Locations on Growth Modules Morphology

The comparison of the measured parameters (Table 4) according to the locations, showed a significant difference between the locations visited for each parameter of module morphology evaluated (P < 0.05). All module dimensions were neither increasing nor decreasing according to the arrangement of the localities along the defined South-North vegetation gradient. However, the modules were well developed in the localities of Bouaké and Katiola in the south of the defined gradient. In *Parkia biglobosa*, the morphology of the modules differed according to the characteristics of the study area (P < 0.05).

### 3.5. Individuals Age Influence on the Growth Modules Morphology

The results showed that the majority of the morphological parameters of the modules differed significantly between the ages (P < 0.05), except for the mean diameter of the modules, which remained statistically identical regardless of age (P > 0.05). The number of modules per axis was higher in mature trees, while the
<table><caption>**Table 4.** Variation of the morphological parameters of growth modules according to the localities surveyed.</caption><thead><tr><th>Localities</th><th>Tree-H (m)</th><th>Tree-Diam (cm)</th><th>Axis-Nbr-Mo</th><th>Mo-leng (cm)</th><th>Mo-Diam (cm)</th><th>Mo-Phyto-Nbr</th></tr></thead><tbody><tr><td><strong>Ferké</strong></td><td>15.51 ± 3.76 <strong>a</strong></td><td>55.28 ± 7.87 <strong>a</strong></td><td>9.39 ± 1.61 <strong>b</strong></td><td>16.32 ± 2.02 <strong>bc</strong></td><td>1.51 ± 0.08 <strong>ab</strong></td><td>16.75 ± 1.03 <strong>ab</strong></td></tr><tr><td>Min-Max</td><td>2.3 - 29.5</td><td>5.4 - 178.44</td><td>5 - 17</td><td>3.85 - 51.4</td><td>0.87 - 3.2</td><td>9.5 - 30.5</td></tr><tr><td><strong>Korhogo</strong></td><td>12.86 ± 2.09 <strong>ab</strong></td><td>38.41 ± 4.06 <strong>ab</strong></td><td>10.14 ± 0.55 <strong>abc</strong></td><td>13.81 ± 1.3 <strong>c</strong></td><td>1.13 ± 0.09 <strong>b</strong></td><td>12.36 ± 0.94 <strong>bc</strong></td></tr><tr><td>Min-Max</td><td>1.4 - 27.5</td><td>2.3 - 84.4</td><td>5 - 17</td><td>4.23 - 26</td><td>0.69 - 1.56</td><td>6.75 - 19.75</td></tr><tr><td><strong>Niakara</strong></td><td>15.98 ± 1.33 <strong>a</strong></td><td>46 ± 3.59 <strong>a</strong></td><td>11.32 ± 0.87 <strong>a</strong></td><td>19.69 ± 1.51 <strong>b</strong></td><td>1.49 ± 0.37 <strong>ab</strong></td><td>20.69 ± 0.62 <strong>a</strong></td></tr><tr><td>Min-Max</td><td>2.4 - 28.5</td><td>6.8 - 83</td><td>7 - 18</td><td>11.13 - 35.4</td><td>0.95 - 2.77</td><td>16 - 33.5</td></tr><tr><td><strong>Katiola</strong></td><td>7.37 ± 2.11 <strong>b</strong></td><td>21.37 ± 3.5 <strong>b</strong></td><td>10.63 ± 0.69 <strong>abc</strong></td><td>23.01 ± 1.35 <strong>ab</strong></td><td>1.67 ± 0.05 <strong>ab</strong></td><td>15.5 ± 0.69 <strong>b</strong></td></tr><tr><td>Min-Max</td><td>0.73 - 28.30</td><td>1.25 - 113.44</td><td>4 - 26</td><td>11.93 - 38.95</td><td>0.97 - 2.16</td><td>8 - 20.5</td></tr><tr><td><strong>Bouaké</strong></td><td>6.45 ± 2.6 <strong>b</strong></td><td>20.46 ± 7.86 <strong>b</strong></td><td>11.12 ± 1.21 <strong>a</strong></td><td>24.59 ± 2.13 <strong>a</strong></td><td>1.76 ± 0.14 <strong>a</strong></td><td>15.84 ± 1.46 <strong>b</strong></td></tr><tr><td>Min-Max</td><td>1.14 - 31.7</td><td>1.79 - 115.61</td><td>3 - 19</td><td>9.88 - 38.5</td><td>0.65 - 2.73</td><td>5.5 - 23.25</td></tr><tr><td><strong>Toumodi</strong></td><td>11.35 ± 0.02 <strong>ab</strong></td><td>29.74 ± 0.08 <strong>ab</strong></td><td>12.2 ± 0.52 <strong>a</strong></td><td>20.55 ± 1.15 <strong>b</strong></td><td>1.93 ± 0.03 <strong>a</strong></td><td>15.28 ± 0.37 <strong>b</strong></td></tr><tr><td>Min-Max</td><td>0.77 - 30</td><td>2.96 - 73</td><td>5 - 19</td><td>9.55 - 41.5</td><td>1.01 - 8.81</td><td>10 - 20.75</td></tr><tr><td><strong>Daloa</strong></td><td>0.64 ± 0.03 <strong>c</strong></td><td>1.33 ± 0.01 <strong>c</strong></td><td>6 ± 0.67 <strong>c</strong></td><td>14.51 ± 2 <strong>c</strong></td><td>0.62 ± 0.1 <strong>c</strong></td><td>8.9 ± 0.97 <strong>c</strong></td></tr><tr><td>Min-Max</td><td>0.42 - 0.87</td><td>0.56 - 2.04</td><td>2 - 13</td><td>9.3 - 31.6</td><td>0.43 - 0.97</td><td>6 - 12</td></tr><tr><td><strong>Pr > F</strong></td><td><strong>0.0001</strong></td><td><strong>0.0001</strong></td><td><strong>0.0009</strong></td><td><strong>0.0001</strong></td><td><strong>0.0089</strong></td><td><strong>0.0001</strong></td></tr></tbody></table>

Values with the same letters are not statistically different at the 5% level.

<table><caption>**Table 5.** Variation in modules morphological parameters according to the individuals sampled age.</caption><thead><tr><th>Ages</th><th>Tree-H (m)</th><th>Tree-Diam (cm)</th><th>Axis-Nbr-Mo</th><th>Mo-leng (cm)</th><th>Mo-Diam (cm)</th><th>Mo-Phyto-Nbr</th></tr></thead><tbody><tr><td><strong>Old tree</strong></td><td>25.23 ± 1.64 <strong>a</strong></td><td>82.36 ± 6.99 <strong>a</strong></td><td>10.61 ± 0.52 <strong>ab</strong></td><td>16.92 ± 0.97 <strong>b</strong></td><td>1.53 ± 0.06 <strong>a</strong></td><td>17.34 ± 0.90 <strong>a</strong></td></tr><tr><td>Min-Max</td><td>18 - 31.7</td><td>60.1 - 178.44</td><td>5 - 14</td><td>5.05 - 26.3</td><td>0.75 - 2.16</td><td>10.5 - 30.5</td></tr><tr><td><strong>Adult tree</strong></td><td>10.81 ± 1.26 <strong>b</strong></td><td>37.47 ± 2.32 <strong>b</strong></td><td>11.21 ± 0.48 <strong>a</strong></td><td>16.94 ± 1.18 <strong>b</strong></td><td>1.59 ± 0.16 <strong>a</strong></td><td>15.82 ± 0.59 <strong>a</strong></td></tr><tr><td>Min-Max</td><td>7.5 - 13.7</td><td>20 - 49.68</td><td>6 - 19</td><td>3.85 - 51.4</td><td>0.69 - 8.81</td><td>6.75 - 25.75</td></tr><tr><td><strong>Young tree</strong></td><td>2.15 ± 0.81 <strong>c</strong></td><td>5.24 ± 1.11 <strong>c</strong></td><td>9 ± 0.63 <strong>b</strong></td><td>20.52 ± 1.1 <strong>a</strong></td><td>1.26 ± 0.07 <strong>a</strong></td><td>13.56 ± 0.72 <strong>b</strong></td></tr><tr><td>Min-Max</td><td>0.42 - 5.8</td><td>0.56 - 18.96</td><td>2 - 26</td><td>9.2 - 41.5</td><td>0.43 - 2.77</td><td>5.5 - 33.5</td></tr><tr><td><strong>Pr > F</strong></td><td><strong>0.0001</strong></td><td><strong>0.0001</strong></td><td><strong>0.0174</strong></td><td><strong>0.0336</strong></td><td><strong>0.0755</strong></td><td><strong>0.0027</strong></td></tr></tbody></table>

Values with the same letters are not statistically different at the 5% level.

average module length was higher in younger trees (Table 5). Old and mature trees had higher numbers of internodes per module than younger trees (Table 5). In *Parkia biglobosa*, the age of the individuals influenced the morphological dimensions of the growth modules (P < 0.05).

## 3.6. Influence of Humidity on the Growth Modules Morphology

The average diameter of the growth modules did not vary significantly between seasons (Table 6); humidity did not have a major influence on the diameter of the growth modules (P > 0.05). In contrast, the number of modules per axis observed, the average length of growth modules and the average
<table><caption>Table 6. Comparison of modules morphological dimensions according to season collection.</caption><thead><tr><th>Seasons</th><th>Tree-H (m)</th><th>Tree-Diam (cm)</th><th>Axis-Nbr-Mo</th><th>Mo-leng (cm)</th><th>Mo-Diam (cm)</th><th>Mo-Phyto-Nbr</th></tr></thead><tbody><tr><td><strong>Rainy</strong></td><td>10.87 ± 1.04 <strong>a</strong></td><td>35.3 ± 2.79 <strong>a</strong></td><td>11.6 ± 0.38 <strong>a</strong></td><td>27.02 ± 1.07 <strong>a</strong></td><td>1.47 ± 0.13 <strong>a</strong></td><td>16.09 ± 0.72 <strong>a</strong></td></tr><tr><td>Min-Max</td><td>0.53 - 31.7</td><td>0.56 - 178.44</td><td>2 - 26</td><td>4.23 - 51.4</td><td>0.45 - 8.81</td><td>6.75 - 25.75</td></tr><tr><td><strong>Dryer</strong></td><td>9.94 ± 1.47 <strong>a</strong></td><td>27.91 ± 4.07 <strong>a</strong></td><td>8.47 ± 0.53 <strong>b</strong></td><td>18.1 ± 0.87 <strong>b</strong></td><td>1.39 ± 0.06 <strong>a</strong></td><td>14.04 ± 0.51 <strong>b</strong></td></tr><tr><td>Min-Max</td><td>0.42 - 30</td><td>0.9 - 85.5</td><td>4 - 18</td><td>3.85 - 41.5</td><td>0.43 - 3.2</td><td>5.5 - 33.5</td></tr><tr><td><strong>Pr > F</strong></td><td><strong>0.3201</strong></td><td><strong>0.143</strong></td><td><strong>0.0001</strong></td><td><strong>0.0043</strong></td><td><strong>0.5666</strong></td><td><strong>0.0201</strong></td></tr></tbody></table>

Values with the same letters are not statistically different at the 5% level.

number of phytomer carried by growth modules differed significantly between seasons (P < 0.05). The majority of the dimensions increased during the rainy season (Table 6). In *Parkia biglobosa*, humidity influenced the majority of the morphological dimensions of the growth modules (P < 0.05), except the mean diameter (P > 0.05).

## 4. Discussion

Growth dynamics in the crown of plants over the course of a season are the consequence of changes in architectural traits. How plants exploit and survive in a given environment, from an abiotic (climate) or biotic (interaction) point of view, is of key importance for their sustainable management. Sustainable management of plant genetic resources must take into account all possible aspects of scientific research, including architectural concepts (retrospective analysis) and the phenological cycle (diagnostic) of trees. This study addressed the effect of habitat, different localities along an environmental gradient, age of individuals and humidity on the morphology of growth modules carried by different axis types in *Parkia biglobosa* in Côte d'Ivoire.

### 4.1. Development of Modules and Architectural Characters of Individuals

The results showed that the boundaries of the Modules are represented by successively very short internodes (Figure 2 and Figure 3). These observations were found on Guayule (*Parthenium argentatum*) [37] [38] [39] and on *Parkia velutina* [1] etc. *Parkia biglobosa* has a sympodial architectural development with a single relaying axis (monochasial sympod), in contrast to Guayule (*Parthenium argentatum*) which is a sympod with several relay axes (polychasial sympod). Guayule belongs to the architectural model of Anthony Leeuwenberg like cassava. [34] indicated that *Parkia biglobosa*, of the Fabaceae family had an architectural development following the model of Wilhelm Troll, we think that the architectural model of this species also integrates that of Paul Champagnat according to our observations. This remains to be verified in further studies on the architectural analysis of *Parkia biglobosa*.
## 4.2. Habitat Influence on Growth Modules

In our study, the environment in which the individuals evolved did not significantly influence the morphological dimensions of the growth modules ($P > 0.05$). This is due to the fact that we considered individuals of all age types observed in this environment in a global manner. The finding might be different if the observations on morphological dimensions were made only and separately by age category. Alternatively, this fact is due to the nature of the species; indeed, some species respect their own genetic programme of architectural development regardless of environmental stress [34]. Moreover, this species is a slow-growing savannah [9], and therefore the immediate perception of the effect of light on shoot growth might be difficult. Normally in a given plant, individuals in the shade (undergrowth) should be less developed than those exposed to the sun (full sun) as in the studies of [3] on *Fagus sylvatica* and of [40]. However, plants growing in the shade are very often tall and etiolated because of the terminal bud's search for bright light for photosynthetic activity. Indeed, bright light is the key factor for growth, it guarantees a good growth in length and width of plants contrary to dark light [40] [41] [42].

## 4.3. Influence of Locality on Growth Modules

The variation in morphological parameters between locations is due to the specific and distinct environmental characteristics of each location surveyed (climate, soil type, rainfall, relative humidity etc.). We had thought that the morphological dimensions of the measured Modules would be ordered in an increasing or decreasing manner according to the layout of the surveyed localities. This was not the case; in fact, even if the surveyed localities are arranged according to a “pseudo-climatic” gradient of dense vegetation decreasing towards the north, each locality has its own distinctive characteristics. Moreover, in the current context of climate change, climatic and microclimatic parameters are very unstable, uncontrolled and vary greatly within a day. Several studies have shown that the depth and fertility of the soil, the environment and even the age of individuals influence the architectural development of species [30] [42] [43]. [44] and [45] have also indicated in their research that climate or an ecological gradient has an effect on plant morphology. Similarly, studies by [46] and [47] [48] indicated that the origin of differences in tree morphology is due to factors such as soil type and the genetic characteristics of the individuals sampled. On the other hand, the modules were well developed in the southern part of the refining gradient (Bouaké and Katiola). This area could be a favourable and priority area for establishing permanent plots and increasing the success rate of reforestation or agroforestry programmes using this species.

## 4.4. Influence of Individuals Age on the Growth Modules
Morphology

The morphological dimensions of the modules evaluated varied with the age of
the individuals (P < 0.05). The number of phytomeres per module is higher in old and mature trees, due to their age and the preformation of organs in the latent buds. In young trees, module length is higher due to the growth dynamics of young shoots (more intense meristematic activity). Indeed, some studies have shown that the age of individuals influences tree morphology [31] [43] [47] [48]. Age strongly modifies the morphological dimensions of plant organs regardless of location, due to very different functional capacities (number of photosynthetic leaves, evapotranspiration capacity, mineral storage, root suctions etc.) and structural dimensions (root and leaf sizes, space occupancy etc.). Young individuals with well-developed modules than the other two categories could be suitable ideotypes in the search for resilient genotypes adapted to global change for the establishment of plots based on this species by cuttings.

## 4.5. Influence of Season

In the rainy season the size of the modules increases. This period is favourable to the synchronised bud break of the terminal buds of the axes. In French Guiana, [1] observed the same phenomenon on *Parkia velutina*. The meristems of the buds, which have long remained dormant or in slow growth due to a long dry season, take advantage of this period to maximise their development. Indeed, water is vital for plants and its availability guarantees a good physiological and harmonious functioning of plants (metabolism, photosynthesis, organogenesis etc.). Whereas water stress slows down all physiological activity and causes their decline. [49] has indicated in his research the positive effects of water regime on plant growth and development. In the process of establishing a plot or implementing reforestation or agroforestry programmes based on *Parkia biglobosa*, the rainy season remains the best and most favourable period.

## 5. Conclusion and Perspectives

This study observed and identified the types of variation in the morphology of growth modules as a function of various abiotic parameters in *Parkia biglobosa*. It was found that the developmental environment of the plants did not significantly influence the dimensions of the growth modules (P > 0.05). While the morphological dimensions of the modules differed significantly between locations (P < 0.05). The age of the individuals influences the size of the growth modules (P < 0.05); with the exception of the average diameter of the modules, which remains identical regardless of age (P > 0.05). In the rainy season, the morphological dimensions of the modules increase significantly due to the favourable water regime for synchronised bud break of the terminal buds of the axes (P < 0.05). From this study, we conclude that *Parkia biglobosa* is a monochasial sympod with a relay axis that follows Paul Champagnat's architectural model at the young stage and transitions to Wilhelm Troll's model later at the adult and old stage. We will try to confirm this in a future study on the architectural development of this species. This information contributes to the under-
standing of the functioning of the crown and the adaptation of this species to
various environments. It could guide the choice of the ideal environment for the
establishment of a reforestation or agroforestry programme based on *Parkia big-
lobosa*.

Acknowledgements

This study was financed by the Ministry of Higher Education and Scientific Re-
search of Côte d'Ivoire, the French Development Agency and IRD (Institut de
Recherche pour le Developpement) in the framework of PRESeD-CI 2 (Renewed
Partnership for Research for Development in Côte d'Ivoire) and C2D (Debt Re-
duction Contract) of the AMRUGE-CI project (Support for the Modernization
and Reform of Universities and Grandes Ecoles of Côte d'Ivoire). The authors
are grateful to the Centre de Coopération International de Recherche Agrono-
mique pour le Développement (CIRAD) for providing the technical equipment
necessary to conduct the study.

Author's Declaration

Beda Innocent Adji wrote the article, took the experimental measurements and analyzed the data. Doffou Sélastique Akaffou and Sylvie Annabel Sabatier contributed to methodology development and corrected the article. Kouadio Henri Kouassi, Jerôme Duminil and Yao Patrice Houphouët coordinated the research project.

Data Availability

The datasets generated during and/or analyzed in the current study are available
from the corresponding author on reasonable request.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

[1] Nicolini, E., Beauchêne, J., Leudet., De la Vallée, B., Ruelle, J., Mangenet, T. and Heuret, P. (2012) Dating Branch Growth Units in a Tropical Tree Using Morphological and Anatomical Markers: The Case of *Parkia velutina* Benoist (Mimosoïdeae). *Annals of Forest Science*, **69**, 543-555. https://doi.org/10.1007/s13595-011-0172-1

[2] Sabatier, S., Caraglio, Y. and Drénou, C. (2014) L'architecture des arbres au service des forestiers. *Innovations Agronomiques*, **41**, 119-128.

[3] Nicolini, E. and Caraglio, Y. (1994) L'influence de divers caractères architecturaux sur l'apparition de la fourche chez le *Fagus sylvatica*, en fonction de l'absence ou de la présence d'un couvert. *Canadian Journal of Botany*, **72**, 1723-1734.
https://doi.org/10.1139/b94-213

[4] Nicolini, E. (2000) Nouvelles observations sur la morphologie des unités de croissance du hêtre (Fagus sylvatica L.). Symétrie des pousses, reflet de la vigueur des arbres. *Canadian Journal of Botany*, **78**, 77-87. https://doi.org/10.1139/cjb-78-1-77
[5] Reinhardt, D. and Kuhlemeier, C. (2002) Plant Architecture. *European Molecular Biology Organization (EMBO Reports)*, **3**, 846-851.
https://doi.org/10.1093/embo-reports/kvf177

[6] Barthélémy, D. and Caraglio, Y. (2007) Plant Architecture: A Dynamic, Multilevel and Comprehensive Approach to Plant Form, Structure and Ontogeny. *Annals of Botany*, **99**, 375-407. https://doi.org/10.1093/aob/mcl260

[7] Rosati, A., Paoletti, A., Caporali, S. and Perri, E. (2013) The Role of Tree Architecture in Super High Density Olive Orchards. *Scientia Horticulturae*, **161**, 24-29.
https://doi.org/10.1016/j.scienta.2013.06.044

[8] Millan, M. (2016) Analyse de la variabilité des traits architecturaux des formes de croissance dans les communautés végétales. Thèse de Botanique. Université de Montpellier, Montpellier, 178 p.

[9] Sina, S. (2006) Reproduction et Diversité Génétique chez Parkia biglobosa (Jacq.) G.Don. PhD Thesis, Wageningen University, Wageningen, 102 p.

[10] Ouédraogo, A.S. (1995) *Parkia biglobosa* (Leguminosae) en Afrique de l'Ouest: Bio-systématique et Amélioration. Thèse Doctorat Wageningen, University, Institute for Forestry and Nature Research, IBNDLO, 205p.

[11] Maisharou, A. and Larwanou, M. (2015) Market Potential of Non-Wood Forest Products in the Sahelian Countries. *International Journal of Forestry Review*, **17**, 125-135. https://doi.org/10.1505/146554815816007018

[12] Avana-Tientcheu, M.L.A., Keouna, S., Dongock, N.D. and Mouga, M.B. (2019) Structure des peuplements et potentiel de domestication de Parkia biglobosa dans la région de Tandjilé-Ouest (Tchad). *International Journal of Biological and Chemical Sciences*, **13**, 219-236. https://doi.org/10.4314/ijbcs.v13i1.19

[13] Détienne, P. (1995) Growth Ring Types and Occurrence in Some French Guianan Species. *Bois et Forêt des Tropiques*, **243**, 65-75.

[14] N'dri, K.E., Gore, B.B.N. and Akaffou, D.S. (2019) Agro-Morphological Characterization of *Myrianthus Arboreus* (P. Beauv) in Cote d'Ivoire. *International Journal of Advances in Scientific Research and Engineering*, **5**, 158-168.
https://doi.org/10.31695/IJASRE.2019.33450

[15] Louppe, D. and Ouattara, N. (1996) Les arboretums d'espèces locales en Nord Côte d'Ivoire. IDEFOR, 14 p. http://agritrop.cirad.fr/581418

[16] Hérault, B., Anatole, K.N., N'klo, O., Assandé, A., Fabrice, B., Brahima, C., Doua-Bi, Y., Koffi, Y., Koffi, K.J.C., Konaté, I., Tiéoulé, F., Wourro, F., Zo-Bi, I.C. and Louppe, D. (2020) The Long-Term Performance of 35 Tree Species of Sudanian West Africa in Pure and Mixed Plantings. *Forest Ecology and Management*, **468**, 118171.
https://doi.org/10.1016/j.foreco.2020.118171

[17] Djaha, J.B.A., N'da, A.A.A., Koffi, E.K., Ballo, C.K. and Coulibaly, M. (2012) Croissance et aptitude au greffage de deux génotypes d'anacardier (*Anacardium occidentale* L.) élites utilisées comme porte-greffe en Côte d'Ivoire. *International Journal of Biological and Chemical Sciences*, **6**, 1453-1466. http://ajol.info/index.php/ijbcs
https://doi.org/10.4314/ijbcs.v6i4.5

[18] Séka, K., Diallo, A.H., Kouassi, N.K. and Aké, S. (2009) Incidence du Yam mosaic virus (YMV) et du *Cucumber mosaic virus* (CMV) sur des variétés de *Dioscorea spp*. Cultivées dans les régions de Bouaké et de Toumodi en Côte d'Ivoire. *International Journal Biological Chemic. Science*, **3**, 694-703.
https://doi.org/10.4314/ijbcs.v3i4.47179

[19] N'cho, B.S. (1991) Modélisation de l'accès des racines de maïs (*Zea mays*) à l'azote. Expérimentation au champ en Centre Côte d'Ivoire. Mémoire de DIAT-ESAT, Mont-
pellier, 22 p.

[20] Ettien, D.J.B. (2004) Intensification de la production d'igname (*Dioscorea* spp) par la fertilisation minérale et l'identification de nouvelles variétés en zone forestière et savanicole de Côte d'Ivoire. Thèse de Doctorat unique en science de la terre, option Agro-pédologie, Université de Cocody, Abidjan, 187 p.

[21] Akassimadou, G. and Yao, K. (2014) La variabilité climatique en Côte d'Ivoire: entre perception sociales et réponse agricole. *Cahier Agriculture, 14*, 533-540.

[22] Traoré, k., Sorho, F., Dramane, D.D. and Sylla, M. (2013) Adventices hôtes alternatifs de virus en culture de Solanaceae en Côte d'Ivoire. *Agro Africa, 25*, 231-237.

[23] N'zi, J.C., Kouame, C. and N'guetta, A.S.P. (2019) Influence de quelques paramètres climatiques sur les effectifs de *Bemisia tabaci* sur la tomate (*Solanum lycopersicum* L.). *International Journal of Biological and Chemical Sciences, 13*, 338-352.
https://doi.org/10.4314/ijbcs.v13i1.27

[24] Coulibaly, L., Kouassi, K.H., SORO, G.E and Savane, I. (2016) Analyse du processus de savanisation du nord de la Côte d'Ivoire par télédétection: Cas du département de Ferkessédougou. *International Journal of Innovation and Applied Studies, 17*, 136-143. http://www.ijias.issr-journals.org/

[25] N'guessan, K.A., Diarassouba, N., Alui, K.A., Nangha, K.Y., Fofana, I.J. and Yao-Kouamé, A. (2015) Indicateurs de dégradation physique des sols dans le Nord de la Côte d'Ivoire: Cas de Boundiali et Ferkessédougou. *Afrique Science, 11*, 115-128.

[26] Coulibaly, T.J.H. (2009) Répartition spatiale, gestion et exploitation des eaux sou-terraines: cas du département de Katiola, région des savanes de Côte d'Ivoire. Thèse de Doctorat, Sciences de la Terre. Université Paris-Est, 142 p.
https://tel.archives-ouvertes.fr/tel-00638690

[27] Soro, T.D., Kouakou, B.D., Kouassi, E.A., Soro, G., Kouassi, A.M., Kouadio, K.E., Yéi, M.S.O. and Soro, N. (2013) Hydroclimatologie et dynamique de l'occupation du sol du bassin versant du Haut Bandama à Tortiya (Nord de la Côte d'Ivoire). *VertigO—La revue électronique en sciences de l'environnement, 13*, 23 p.
https://doi.org/10.4000/vertigo.14468

[28] Adjanohoun, E. (1964) Végétation des savanes et des rochers découverts en Côte d'ivoire Centrale. Rapport, Office de la recherche scientifique et technique outre-Mer, Centre ORSTOM d'Adiopodoumé-Abidjan, 238 p.

[29] Hallé, F and Martin, R. (1968) Etude de la croissance rythmique chez l'Hévéa (*He-vea brasiliensis* Mull. Arg.). *Adansonia, Série 2, 8*, 475-503.

[30] Sabatier, S. (1999) Variabilité morphologique et architecturale de deux espèces de Noyers: *Juglans regia* L., *Juglans nigra* L. et deux Noyers hybrides interspecifiques. Thèse de Doctorat, Biologie végétale, sciences et techniques du Languedoc, Université de Montpellier II, 143 p. https://tel.archives-ouvertes.fr/tel-00106305

[31] Sabatier, S. and Barthélémy, D. (1999) Growth Dynamics and Morphology of An-nual Shoots According to Their Architectural Position in Young *Cedrus atlantica* (Endl.) Manetti ex Carrière (Pinaceae). *Annals of Botany, 84*, 387-392.
https://doi.org/10.1006/anbo.1999.0939

[32] Thiébaut, B. (1982) Observations sur le développement de plantules de hêtre (*Fagus sylvatica*) cultivés en pépinières, orthotropie et plagiotropie. *Canadian Journal of Botany, 60*, 1292-1303. https://doi.org/10.1139/b82-165

[33] Hallé, F., Oldeman, R.A.A. and Tomlinson, P.B. (1978) Tropical Trees and Forests: An Architectural Analysis. Springer-Verlag, Berlin.
https://doi.org/10.1007/978-3-642-81190-6

[34] Hallé, F. and Keller, R. (2019) Mais d'où viennent les plantes? Actes sud., 179 p.
[35] Sabatier, S., Barthélémy, D., Ducousso, I. and Germain, E. (1998) Modalités d'allongement et morphologie des pousses annuelles chez le noyer commun, *Juglans regia* L. 'Lara' (Juglandaceae). *Canadian Journal of Botany*, **76**, 1253-1264.
https://doi.org/10.1139/b98-055

[36] Godin, C. and Caraglio, Y. (1998) A Multiscale Model of Plant Topological Structures. *Journal of Theoretical Biology*, **191**, 1-46.
https://doi.org/10.1006/jtbi.1997.0561

[37] Taurines, M., Brancheriau, L., Palu, S., Pioch, D., Tardan, E., Boutahar, N., Sartre, P. and Meunier, F. (2019) Determination of Natural Rubber and Resin Content of Guayule Fresh Biomass by Near Infrared Spectroscopy. *Industrial Crop & Products*, **134**, 177-184. https://doi.org/10.1016/j.indcrop.2019.03.073

[38] Jaeger, M., Sabatier, S.A., De Reffye, P., Alsater, A.A. and Brancheriau, L. (2020) Efficient Stochastic Functional Structural Sympodial Shrubs Modelling Base DonStructural Hierarchy. Application to Guayule. FSPM 2020: Towards Computable Plants, Leibniz, Germany, hal-03105340.

[39] Brancheriau, L., Sabatier, A.S., Jaeger, M., Hemery, N. and Mougani, C. (2020) Modelling the Guayule Plant Growth and Development with a Functional Structural Plant Model. FSPM towards Computable Plants. 9th International Conference on Functional-Structural Plant models, Hanovre, hal-03056119.

[40] Tousignant, M.E. and Delorme, M. (2006) Connaître le fonctionnement de la plante pour mieux gérer son environnement. Article paru dans Québec Vert (IQDHO). Adapté de Understanding Plant Growth: A Review of the Basics. Conférence présentée lors de l'Ohio International Short Course 2005 par Paul A. Thomas et Bodie V. Pennisi, de l'Université de Georgie, et bonifié par l'équipe de l'IQDHO, 3 p.

[41] Abidi, F. (2012) Effets de la qualité de la lumière sur l'élaboration de l'architecture du rosier buisson. Thèse de doctorat, l'Université de Tunis El Manar, Tunisie, 277 p.

[42] Calonnec, A. (2013) Façonner l'architecture végétale pour contrôler les maladies des plantes. *BIOFUTUR*, **343**, 37-42.

[43] Lauri, P.E., Kelner. J.J., Trottier, J.C and Costes, E. (2010) Insights into Secondary Growth in Perennial Plants: Its Unequal Spatial and Temporal Dynamics in the Apple (*Malus domestica*) Is Driven by Architectural Position and Fruit Load. *Annals of Botany*, **105**, 607-616. https://doi.org/10.1093/aob/mcq006

[44] Maranz, S. and Wiesman, Z. (2003) Evidence for Indigenous Selection and Distribution of the Shea Tree, *Vitellaria paradoxa*, and Its Potential Significance to Prevailing Parkland Savanna Tree Patterns in Sub-Saharan Africa North of the Equator. *Journal of Global Biogeography*, **30**, 1505-1516.
https://doi.org/10.1046/j.1365-2699.2003.00892.x

[45] Soloviev, P., Niang, T.D., Gaye, A. and Totte, A. (2004) Variabilité des caractères physicochimiques des fruits de trois espèces ligneuses de cueillette, récoltés au Sénégal, *Adansonia digitata*, *Balanites aegyptiaca* et *Tamarindus indica*. *Fruits*, **59**, 109-119. https://doi.org/10.1051/fruits:2004011

[46] Salazar, R. and Quesada, M. (1987) Provenance Variation in *Guzuma ulmifolia* L. in Costa Rica. *Commonwealth Forestry Review*, **66**, 317-324.

[47] Assogbadjo, A.E., Sinsin, B., Codjia, J.T.C. and Van Damme, P. (2005) Ecological Diversity and Pulp, Seed and Kernel Production of the Baobab (*Adansonia digitata*) in Benin. *Belgian Journal of Botany*, **138**, 47-56.

[48] Assogbadjo, A.E., Kyndt, T., Sinsin, B., Gheysen, G. and Van Damme, P. (2006) Patterns of Genetic and Morphometric Diversity in Baobab (*Adansonia digitata* L.) Populations across Different Climatic Zones of Benin (West Africa). *Annals of Botany*, **97**, 819-830. https://doi.org/10.1093/aob/mcl043
[49] Meridja, S. (2011) Approche biophysique des processus de développement et de croissance des couverts végétaux: Interaction avec le stress hydrique et optimisation des pratiques culturales en climat méditerranéen. Thèse en Sciences agricoles, AgroParisTech, Français, NNT: 2011AGPT0057, pastel-00776775, 180 p.