# Effet de la gestion des feuilles d'élagage du palmier à huile sur le stock de carbone et les propriétés physico-chimiques du sol dans les palmerai es villageoises du Bénin

Hervé Nonwègnon Sayimi Aholoukpè (1), Guillaume Lucien Amadji (2), Didier Blavet (3), Jean-Luc Chotte (3), Philippe Deleporte (4), Bernard Dubos (5), Albert Flori (5), Christophe Jourdan (4)

Reçu le 29 avril 2015, accepté le 12 novembre 2015.

**Description du sujet.** Au Bénin, la culture du palmier à huile est largement répandue en milieu villageois. L'étude porte sur l'impact du recyclage des feuilles d'élagage sur les propriétés du sol des palmerai es.

**Objectifs.** Évaluer l'effet de deux modes de gestion des feuilles élaguées dans la palmerai e sur le stock de carbone et les propriétés physico-chimiques du sol.

**Méthode.** Six palmerai es villageoises de deux classes d'âges (7-12 ans et 13-24 ans) ont été retenues dans le sud-est du Bénin afin de comparer les propriétés du sol de deux zones : les zones de restitution totale des feuilles d'élagage (RT) et celles de restitution nulle des feuilles (RN) représentant le témoin. Dans chaque palmerai e, le sol sous chaque modalité de recyclage des feuilles a été prélevé, sur fosses, jusqu'à 50 cm de profondeur à proximité de six palmiers. Les propriétés du sol ont été déterminées par des méthodes d'analyse conventionnelles et par la spectroscopie proche infrarouge.

**Résultats.** Les stocks de carbone et d'azote sous RT à l'âge adulte sont de 58 Mg·ha⁻¹ et 4,3 Mg·ha⁻¹, respectivement et sont significativement supérieurs aux stocks dans tous les autres traitements (34,6 Mg C·ha⁻¹ et 2,9 Mg N·ha⁻¹). Le recyclage des feuilles a amélioré la fertilité du sol jusqu'à 20 cm de profondeur dans les palmerai es adultes. La porosité et la densité apparente du sol ont été aussi améliorées.

**Conclusions.** La gestion des feuilles élaguées dans les palmerai es modifie les propriétés du sol qui sont nettement améliorées en surface après 10 ans de recyclage des feuilles.

**Mots-clés.** Elaeis guineensis, pratique culturale, élimination des feuilles, propriétés physico-chimiques du sol, carbone, azote, fertilité du sol, Bénin.

## Effect of management of pruned fronds of oil palm on soil carbon stock and soil physico-chemical properties in Beninese smallholder plantations

**Description of the subject.** In Benin, oil palm plantations are widely managed by smallholder farmers. The study addresses the impact on soil properties of the recycling of pruned fronds from these trees.

**Objectives.** The study aims to assess the effect of two modes of management of pruned fronds on the carbon stock and physico-chemical properties of soil.

**Method.** Six oil palm plantations were selected to represent two age classes (7-12 years and 13-24 years) in southeastern Benin. They were used to compare the soil properties of two areas: the areas of total recycling of pruned fronds (TR) and the areas of no recycling of the fronds (NR), which represented the control. In each plantation, soil samples were collected under the two treatments, from holes of up to 50 cm depth, near six oil palm trees. The soil properties were determined using conventional laboratory methods and near-infrared spectroscopy method.

**Results.** Carbon and nitrogen stocks under TR of adult plantations were 58 Mg·ha⁻¹ and 4.3 Mg·ha⁻¹, respectively, and were significantly higher than the stocks in all other treatments (34.6 Mg C·ha⁻¹ and 2.9 Mg N·ha⁻¹). Recycling of pruned fronds

<footnote>(1) Centre de Recherches agricoles Plantes Pérennes. INRAB. BP 01 Pobè (Bénin). E-mail : aholoukpeherve@yahoo.fr</footnote>

<footnote>(2) Université d'Abomey-Calavi. Faculté des Sciences agronomiques. BP 526. Cotonou (Bénin).</footnote>

<footnote>(3) IRD. UMR Éco&Sols. 2, Place Viala. F-34060 Montpellier Cedex 2 (France).</footnote>

<footnote>(4) CIRAD. UMR Éco&Sols. 2, Place Viala. F-34060 Montpellier Cedex 2 (France).</footnote>

<footnote>(5) CIRAD. UPR Systèmes de Pérennes. Avenue Agropolis. F-34398 Montpellier Cedex 5 (France).</footnote>
improved soil fertility at a depth of up to 20 cm in the adult plantations. The porosity and bulk density of the soil were also improved.

**Conclusions.** The management of the pruned fronds in the smallholders’ oil palm plantations modified the soil properties, which were significantly improved in the top soil after 10 years of frond recycling.

**Keywords.** *Elaeis guineensis*, cultural practices, defoliation, physico-chemical properties of soil, carbon, nitrogen, soil fertility, Benin.

# 1. INTRODUCTION

Le palmier à huile (*Elaeis guineensis*) cultivé sous les latitudes tropicales et subtropicales humides, est l'oléagineux ayant le rendement en corps gras le plus élevé (Barcelos, 1998 ; Corley et al., 2003). Originaire d'Afrique, le palmier à huile a connu une forte expansion dans le monde, en particulier en Asie du Sud-Est, où la Malaisie et l'Indonésie détiennent aujourd'hui plus de 80 % des superficies mondiales plantées (Ahoyo, 2008). Ce succès asiatique s'explique par des conditions climatiques globalement plus favorables que celles de son aire d'origine (Nodichao, 2008). En Afrique de l'Ouest, où il est cultivé depuis la fin du 19e siècle, les conditions climatiques sont parfois marginales et limitent fortement les rendements, comme c'est le cas au Bénin. Dans ce pays où la culture n'est maintenant pratiquée que par des villageois, la station de recherche sur le palmier à huile située à Pobè a développé un programme spécifique de mise au point de matériel végétal adapté à la sècheresse (Nouy et al., 1999). Parallèlement, les programmes de recherche ont étudié l'amélioration de la productivité par une fertilisation minérale adaptée, la mise au point de techniques culturales et la gestion de la matière organique du sol dans le but de réduire l'impact des déficits hydriques (Djègui et al., 1992 ; Adjè, 2004).

Dans une palmeraie, les principales sources de matière organique qui peuvent être recyclées tout au long du cycle de culture sont constituées des rafles après égrappage des régimes, des sous-produits d'usinage (boues et effluents) et des feuilles d'élagage. La biomasse des troncs des palmiers ne peut être recyclée qu'au moment de la replantation. Tous ces produits contribuent à augmenter le stock de matière organique du sol et à améliorer les propriétés physiques et chimiques du sol, ce qui a un effet positif sur le rendement du palmier (Khalid et al., 2000a ; Khalid et al., 2000b ; Corley et al., 2003). Parmi les éléments nutritifs restitués par le recyclage des feuilles, l'azote et le potassium sont deux éléments majeurs en raison des quantités libérées dans le sol et en raison de leur importance pour la nutrition minérale du palmier à huile (Khalid et al., 1999 ; Dubos et al., 2011).

Dans les palmeraies villageoises du Bénin en phase productive, les biomasses qui sont recyclées directement au champ sont les résidus des cultures vivrières associées aux jeunes palmiers jusqu'à 5 ans,

puis les feuilles d'élagage des palmiers (Djègui et al., 1992 ; Aholoukpè et al., 2013). Ces feuilles représentent en moyenne 8 t et 4 t de biomasse sèche par hectare par an respectivement dans une palmeraie adulte et pré-adulte (Aholoukpè, 2013). Elles peuvent être :

- soit disposées en tas sur le sol où elles se décomposeront ; elles forment des andains discontinus situés à mi-distance sur la ligne séparant deux palmiers consécutifs (Figure 1) ;

- soit exportées hors de la parcelle pour un usage domestique (combustible de cuisine, palissade, etc.).

C'est ce dernier usage qui est le plus fréquent (Aholoukpè et al., 2013). Si l'exportation de la biomasse du palmier à huile a un effet négatif sur le maintien de la matière organique et sur la fertilité du sol (Khalid et al., 2000a), très peu d'études ont toutefois évalué le potentiel de la pratique de recyclage des résidus de récolte du palmier à huile dans l'amélioration des stocks organiques et la restitution des nutriments du sol (Djègui, 1992 ; Law et al., 2009 ; Nelson et al., 2011).

Dans une première étude sur la station de recherche de Pobè, Djègui (1995) a travaillé sur la variabilité spatiale de la matière organique du sol liée d'une part aux pratiques culturales propres aux cultures vivrières conduites sur billons dans les jeunes plantations et, d'autre part, au mode de plantation en ligne avec restitution en andains des feuilles d'élagage. Ces essais conduits en station n'ont cependant pas pris en compte l'effet sur le long terme de la restitution des feuilles d'élagage dans les conditions des palmeraies villageoises où la pratique de fertilisation disparait au moment de la mise en production. L'objectif général de l'étude présentée dans cet article est d'évaluer l'effet du recyclage ou non de la biomasse des feuilles sur les stocks de matière organique et la fertilité du sol en milieu réel, en étudiant des palmeraies villageoises soumises pendant des durées plus ou moins longues aux différentes pratiques. De façon spécifique, deux questions de recherche ont été étudiées :

- effet du mode de gestion des feuilles élaguées et de l'âge des palmeraies sur les stocks de carbone et d'azote du sol ;

- effet du mode de gestion des feuilles élaguées sur les propriétés physico-chimiques des sols des palmeraies adultes.
<figure><img src="image_1.png" /></figure>

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

**Figure 1.** Espace de Voronoï autour du palmier étudié (a) et triangle de Voronoï (b) présentant les dimensions de la zone Z3s de prélèvement du sol sous andain et dans l'interligne — Complete Voronoï space around the studied palm tree (a) and Voronoï triangle (b) with the dimensions of Z3s zone for the soil sampling under the pruned fronds pile and inside the interrow.

a. L'espace de Voronoï est l'hexagone défini par six palmiers qui entourent le palmier central étudié. Il est constitué de 12 triangles rectangles égaux (Nodichao, 2008). Pour les six palmiers choisis par palmeraie, il a été installé six espaces de Voronoï dont le palmier central (Pc) est chaque fois utilisé pour les observations — *The Voronoï space is the hexagon delimited by six palm trees which surround the studied palm tree. This space is split into 12 equivalent right triangles (Nodichao, 2008). For the six palm trees chosen by palm plantation, six Voronoï spaces have been set up and the central tree (Pc) was each time used for the observations* ; b. Le triangle de Voronoï est composé de 5 zones dont les zones Z1, Z2 et Z3s délimitent l'espace simplifié (Dassou, 2011). Les observations ont été faites dans la zone 3 simplifiée (Z3s) — *The Voronoï triangle includes 5 zones of which Z1, Z2 and Z3s delimit the simplified space (Dassou, 2011). The observations were done inside the simplified zone 3 (Z3s).*

# 2. MATÉRIEL ET MÉTHODES

## 2.1. Zone d'étude

L'étude a été conduite dans le département du Plateau situé au sud-est du Bénin (2°30' - 2°45' E et 6°35' -

7°45' N). Le climat est de type subéquatorial, caractérisé par une alternance de deux saisons pluvieuses (avril à juillet et mi-septembre à mi-novembre) et deux saisons sèches (aout à mi-septembre et mi-novembre à mars). Dans la zone d'étude, on enregistre une hauteur moyenne de 1300 mm de pluie par an et un
déficit hydrique<fnref n="1" /> du sol estimé entre 600 et 800 mm, ce qui limite la productivité du palmier à huile (Nouy et al., 1999). Les sols étudiés sont largement représentés dans le département et sont les plus cultivés avec le palmier à huile (Aholoukpè et al., 2013). Ils sont communément appelés « terre de barre » et ont été décrits par Azontondé (1991) et Amadji et al. (2008). Selon la nomenclature française, ce sont les sols ferrallitiques faiblement désaturés appauvris en argile (Volkoff et al., 1976). Dans la classification de la FAO, il s'agit de ferralsols rhodiques, molliques sous forêt (teneur en C > 0,6 %) et ochriques (teneur en C < 0,6 %) sous culture (Djègui et al., 1992).

Dans la zone d'étude, la végétation est dominée par des savanes arborées et arbustives, des palmeraies plantées avec du matériel végétal sélectionné (« palmeraies sélectionnées ») et des palmeraies issues de régénération naturelle (« palmeraies naturelles »), des jachères arbustives, des îlots de forêts semi-décidues et des galeries forestières (Djègui, 1992).

## 2.2. Sélection des palmeraies

Les deux communes produisant le plus de palmiers à huile dans le département (Adja-Ouèrè et Sakété) (Adjadi, 2012; Akpo et al., 2012) ont été retenues pour le choix de six palmeraies villageoises distantes entre elles d'un minimum de 10 km et dont quelques caractéristiques sont présentées dans le **tableau 1**. S'agissant d'une étude en milieu villageois, il n'a pas été possible d'identifier des palmeraies d'âges identiques à l'échelle du département. De ce fait, les palmeraies étudiées ont été regroupées en classes d'âges définies dans une étude typologique préalable (Aholoukpè et al., 2013). Les classes des palmeraies pré-adultes (7-12 ans) et adultes (13-24 ans) ont été retenues pour l'étude. Elles ont été plantées selon le même dispositif en triangles équilatéraux de 9 m de coté (143 palmiers par hectare) et ont été conduites selon des itinéraires techniques similaires : précédents culturaux constitués d'anciennes cultures annuelles, association de cultures vivrières au jeune âge, très rares apports d'engrais minéraux à partir de 5 ans, élagage à partir de 7 ans et dépôt en andains des feuilles d'élagage. Dans ces palmeraies, on a distingué les zones où sont déposées toutes les feuilles d'élagage entre deux palmiers sur la ligne de plantation (modalité RT, restitution totale des feuilles) et les zones non concernées par le recyclage dans les interlignes de plantation (modalité RN, restitution nulle). Les andains de feuilles couvrent une surface assimilée à un rectangle centré sur la ligne séparant deux palmiers consécutifs et dont la superficie moyenne est

d'environ 7 m². Ils représentent 10 % de la superficie d'une palmeraie (Aholoukpè, 2013). La modalité RN est assimilable à la pratique d'exportation totale des feuilles d'élagage rencontrée dans les palmeraies villageoises de la zone d'étude (Aholoukpè et al., 2013). Les classes d'âges correspondent respectivement à une moyenne de 3 et 10 ans de recyclage des feuilles élaguées.

L'étude comprend deux phases :

- dans un premier temps, les six palmeraies ont servi à étudier les stocks de carbone et d'azote du sol sur 0-50 cm de profondeur. Le dispositif a permis d'étudier l'effet des deux facteurs « mode de gestion des feuilles d'élagage » (RT ou RN) et « classe d'âge » (pré-adulte : 7-12 ans ou adulte : 13-24 ans) et leur interaction ;

- dans un second temps, les conclusions de la première phase de l'étude nous ont amenés à compléter les observations des trois palmeraies de la classe d'âge adulte par des analyses physico-chimiques (pH, cations échangeables, capacité d'échange cationique, densité apparente et porosité du sol) en distinguant les modalités RT et RN.

## 2.3. Échantillonnage de sol et mesures de terrain

Dans chaque palmeraie étudiée, les observations ont été faites pour six palmiers tirés de façon aléatoire sur une superficie de 1 ha. La surface plantée a été structurée selon un dispositif de Voronoï (Figure 1) qui définit des secteurs triangulaires équivalents autour du palmier (Nodichao, 2008). Pour comparer les modalités RT et RN, nous avons restreint les zones d'études aux rectangles Z3s du dispositif de Voronoï simplifié proposé par Dassou (2011). Des fosses de dimensions 1,80 m x 0,70 m x 1 m (longueur, largeur et profondeur) ont été creusées dans chacune des zones (andains et interligne) des six palmiers de la parcelle. Les analyses physico-chimiques et la détermination de la densité apparente du sol ont été faites sur les couches de profondeurs 0-5, 5-10, 10-20, 20-30 et 30-50 cm. Aucun profil de sol n'a donc été pris en compte dans cet échantillonnage. Pour chaque fosse, les quatre faces ont été utilisées pour prélever les échantillons de sol qui ont servi à constituer des échantillons composites pour les analyses chimiques de chaque fosse. Au total, 360 échantillons composites de sol ont été réalisés sur la base des cinq profondeurs de prélèvement, des deux fosses par palmier, des six palmiers par palmeraie et des six palmeraies retenues pour l'étude. Chaque face de la fosse a aussi été utilisée pour déterminer la densité apparente moyenne du sol au moyen d'un cylindre de 3 cm de diamètre et de 5 cm de hauteur.

## 2.4. Étude de laboratoire

Les sols ont été séchés à l'air et tamisés à 2 mm pour les analyses de laboratoire. L'étude au laboratoire a conduit

<footnote n="1">¹ Une estimation du déficit hydrique est obtenue à partir du bilan simplifié basé sur la pluviosité mensuelle, une réserve hydrique du sol fixée à 150 mm et une ETP de 100 à 150 mm par mois selon la saison (Nouy et al., 1999).</footnote>
<table><caption>**Tableau 1.** Caractéristiques des palmerai es étudiées — *Characteristics of the studied palm plantations.*</caption><thead><tr><th>N° palmeraie</th><th>Coordonnées géographiques du centre de la parcelle</th><th>Superficie (ha)</th><th>Âge de la plantation (an)</th><th>Classe d'âge</th><th>Durée de gestion des feuilles (an)</th></tr></thead><tbody><tr><td>1</td><td>6°49'26,1" N 2°38'25,2" E</td><td>2</td><td>11</td><td>Pré-adulte</td><td>4</td></tr><tr><td>2</td><td>6°47'59,4" N 2°37'01,6" E</td><td>2</td><td>9</td><td>Pré-adulte</td><td>2</td></tr><tr><td>3</td><td>6°49'39,1" N 2°38'51,9" E</td><td>3</td><td>10</td><td>Pré-adulte</td><td>3</td></tr><tr><td>4</td><td>6°49'45,1" N 2°38'29,3" E</td><td>4</td><td>17</td><td>Adulte</td><td>10</td></tr><tr><td>5</td><td>6°45'35,0" N 2°37'38,2" E</td><td>7</td><td>17</td><td>Adulte</td><td>10</td></tr><tr><td>6</td><td>6°59'06,3" N 2°40'34,7" E</td><td>2,5</td><td>18</td><td>Adulte</td><td>11</td></tr></tbody></table>

à l'obtention des teneurs des propriétés du sol sur les 360 échantillons. La méthode analytique adoptée a consisté en l'utilisation de la spectrométrie proche infrarouge pour la prédiction des teneurs avec des modèles de régression établis sur la base de résultats d'analyses physico-chimiques classiques faites sur un nombre plus réduit d'échantillons. Le principe a été basé sur les étapes suivantes :

- choix des échantillons de calibration et de validation des modèles de prediction. Le passage de tous les échantillons au spectromètre a permis d'acquérir leurs données spectrales qui ont été analysées séparément avec le logiciel WinISI III-v1.61e (Foss NIRSystems/ Tecator International, LLC, Silver Spring, MD, USA). Une analyse en composantes principales (ACP) a été réalisée sur les absorbances. Une centaine d'échantillons (exactement 101 échantillons dans cette étude) les plus représentatifs ont été sélectionnés par le logiciel pour constituer le lot d'échantillons de calibration. La sélection se base sur le concept du voisinage (distance entre deux échantillons voisins) dans l'espace des composantes principales. Par la suite, 30 échantillons ont été sélectionnés sur les 259 échantillons restants pour constituer le lot d'échantillons de validation ;

- analyses physico-chimiques des échantillons de sol. Les échantillons sélectionnés pour la calibration et la validation des modèles de prédiction (131 échantillons) ont été analysés au laboratoire selon les méthodes conventionnelles décrites par Pansu et al. (2006). Les teneurs en carbone et azote organiques ont été déterminées par la méthode du CHN par combustion sèche totale du sol sous oxygène avec comme gaz vecteur l'hélium pur à 870 °C. Le pH du sol a été déterminé au moyen d'un pHmètre après une mise en solution aqueuse du sol pendant 2 h, dans un rapport sol/eau de 1/2,5. Une extraction à l'acétate d'ammonium (1N, pH = 7) dans un rapport sol/solution de 1/20, puis un lavage du sol par le chlorure de potassium (1N) pour déplacer en solution les ions $NH_4^+$, ont permis de déterminer les cations échangeables, calcium (Ca), magnésium (Mg) et potassium (K) et la capacité d'échange

cationique (CEC). Les cations échangeables sont dosés dans l'extrait d'acétate d'ammonium au moyen d'un spectromètre d'émission plasma (ICP). Les ions $NH_4^+$ déplacés en solution représentent la CEC et sont dosés par colorimétrie ;

- prédiction des propriétés du sol. Les résultats des 131 analyses conventionnelles ont servi à établir des modèles de régression pls (partial least square) pour prédire les propriétés du sol sur les autres échantillons (Barthès et al., 2006; Brunet et al., 2007). La précision des modèles de prédiction par la SPIR a été appréciée en calibrant le modèle à partir des analyses des 101 échantillons de calibration et en le validant sur les données d'analyses des 30 échantillons restants. Les paramètres suivants ont été examinés : $R^2_{Cal}$ (coefficient de détermination de la calibration des modèles), $R^2_{val}$ (coefficient de validation des modèles) et $RPD_{cv}$ (ratio de déviation de la performance de validation croisée). La prédiction est d'autant meilleure que $R^2_{Cal}$ et $R^2_{val}$ sont proches de 1 et selon Chang et al. (2001), la validation croisée des modèles est bonne si $RPD_{cv} > 2$, elle est moyenne si $1,4 < RPD_{cv} < 2$ et mauvaise si $RPD_{cv} < 1,4$. Le **tableau 2** présente les valeurs obtenues pour les modèles retenus pour la prédiction des propriétés de nos échantillons.

## 2.5. Traitement des données

**Calcul de la densité apparente, de la porosité et des stocks de carbone et d'azote du sol.** La densité apparente du sol désigne la masse volumique du sol. Elle a été déterminée par le rapport de la masse sèche du sol par le volume du cylindre de prélèvement. Les valeurs de densité apparente du sol ont servi à calculer la porosité du sol selon la formule :

$$ \text{Porosité (\%)} = (1-\text{Da}/\text{Dr})*100 $$

où Da est la densité apparente du sol et Dr la densité réelle du sol. Pour Dr, nous avons retenu la valeur de 2,65 dont Aho et al. (1997) ont montré qu'elle représentait bien celle des sols ferrallitiques.
<table><caption>**Tableau 2.** Paramètres de calibration, de validation croisée et de validation externe des modèles de prédiction des propriétés du sol par la méthode de la spectroscopie proche infrarouge basée sur des régressions pls (partial least square) — Parameters of calibration, cross validation and external validation of the models of soil properties prediction established by partial least square regressions (pls) on near infrared spectroscopy data.</caption><thead><tr><th rowspan="2">Propriété du sol</th><th rowspan="2">Meilleur pré-<br>traitement<br>spectral</th><th colspan="3">Calibration<br>pls</th><th colspan="3">Validation croisée<br>pls</th><th colspan="4">Validation externe pls</th><th rowspan="2">Précision de<br>la prédiction<br>selon Chang<br>et al. (2001)</th></tr><tr><th>N1</th><th>SE<sub>CV</sub></th><th>R<sup>2</sup><sub>CV</sub></th><th>RPD<sub>CV</sub></th><th>N2</th><th>SEP</th><th>R<sup>2</sup><sub>val</sub></th><th>RPD<sub>val</sub></th></tr></thead><tbody><tr><td>C<sub>org</sub> (g·kg<sup>-1</sup>)</td><td>Std MSC 0011</td><td>95</td><td>1,50</td><td>0,97***</td><td>4,9</td><td>30</td><td>2,74</td><td>0,91***</td><td>3,0</td><td>Bonne</td></tr><tr><td>N<sub>tot</sub> (g·kg<sup>-1</sup>)</td><td>Std MSC 0011</td><td>94</td><td>0,11</td><td>0,95***</td><td>3,8</td><td>30</td><td>0,20</td><td>0,89***</td><td>2,6</td><td>Bonne</td></tr><tr><td>pH<sub>eau</sub></td><td>Detrend 1441</td><td>99</td><td>0,28</td><td>0,81***</td><td>1,9</td><td>30</td><td>0,24</td><td>0,72***</td><td>1,9</td><td>Moyenne</td></tr><tr><td>Ca (meq·100 g<sup>-1</sup>)</td><td>Std MSC 1441</td><td>94</td><td>0,44</td><td>0,94***</td><td>3,0</td><td>30</td><td>0,68</td><td>0,87***</td><td>2,5</td><td>Bonne</td></tr><tr><td>Mg (meq·100 g<sup>-1</sup>)</td><td>Std MSC 1441</td><td>95</td><td>0,29</td><td>0,89***</td><td>2,5</td><td>30</td><td>0,39</td><td>0,78***</td><td>2,1</td><td>Bonne</td></tr><tr><td>K (meq·100 g<sup>-1</sup>)</td><td>SNV 1441</td><td>93</td><td>0,02</td><td>0,88***</td><td>2,1</td><td>30</td><td>0,01</td><td>0,80***</td><td>2,1</td><td>Bonne</td></tr><tr><td>S (meq·100 g<sup>-1</sup>)</td><td>Std MSC 0011</td><td>95</td><td>0,76</td><td>0,93***</td><td>3,0</td><td>30</td><td>0,88</td><td>0,90***</td><td>2,9</td><td>Bonne</td></tr><tr><td>CEC (meq·100 g<sup>-1</sup>)</td><td>SNV 0011</td><td>94</td><td>0,87</td><td>0,91***</td><td>2,8</td><td>30</td><td>1,19</td><td>0,77***</td><td>2,1</td><td>Bonne</td></tr></tbody></table>

Pré-traitement spectral — spectra pre-treatments : Std MSC 0011 = standard multiplicative scatter correction ; 1er chiffre : ordre de la dérivée, 2e chiffre : nombre de longueurs d'ondes sur lesquelles la dérivation a été faite, 3e et 4e chiffres : niveau de lissage du modèle — 1er digit: derivative order, 2nd digit: number of wavelengths over which the derivative is applied, 3rd and 4th digits: model smoothing levels ; SNV : standard normal variate ; N1 : nombre d'échantillons utilisés pour la calibration des modèles après élimination des échantillons déviants. Le nombre d'échantillons retenus au départ pour la calibration est de 101 — number of samples after the elimination of calibration outliers, 101 samples were choosen for calibration ; N2 : nombre d'échantillons utilisés pour la validation externe — number of samples used for the external validation ; SE<sub>CV</sub>, SEP : respectivement erreurs standard de la validation croisée et de la prédiction — standard errors of cross validation and prediction, respectively ; R<sup>2</sup><sub>CV</sub>, R<sup>2</sup><sub>val</sub> : coefficients de détermination de la validation croisée et de la validation externe — coefficients of determination of cross validation and external validation ; RPD<sub>CV</sub>, RPD<sub>val</sub> : ratios de déviation de la performance de la validation croisée et de la validation externe — ratios of performance to deviation for cross validation and external validation ; *** : significatif à p < 0,01 — significant at p < 0.01.

Les valeurs de Da et les teneurs en carbone et en azote ont été utilisées pour calculer les stocks de carbone et d'azote par masse de terre équivalente (MTE) pour les premiers 650 kg·m<sup>-2</sup> de sol. Cette valeur est proche de la masse totale de terre fine calculée sur la profondeur de 50 cm sous andain ou dans l'interligne. Le calcul de la masse totale de terre fine et des stocks de carbone et d'azote a été opéré selon la procédure suivante tirée des travaux de Ellert et al. (1995) et Ellert et al. (2002) :

- on a déterminé que l'on pouvait atteindre une MTE de 650 kg·m<sup>-2</sup> de sol pour toutes les fosses en utilisant les données des cinq couches de sol ;

- pour chaque fosse, on a calculé l'épaisseur exacte de la couche la plus profonde pour atteindre exactement la MTE, puis on a calculé l'épaisseur totale exacte (Ec) correspondante ;

- enfin, on a calculé le stock de C pour l'épaisseur Ec de sol en cumulant les stocks C<sub>i</sub> de chaque couche.

Le stock de C (QC) par MTE a été obtenu par la formule suivante :

$$QC = \sum_{i=1}^{i=5} C_i \times Da_i \times E_i \times TF_i$$

avec QC en Mg·ha<sup>-1</sup> ; teneurs pondérales en carbone C<sub>i</sub> en g·100 g<sup>-1</sup> ; densité apparente Da<sub>i</sub> en g·cm<sup>-3</sup> (ou kg·dm<sup>-3</sup>) ; épaisseurs E<sub>i</sub> en cm pour la couche i ; taux de terre fine TF<sub>i</sub> en g·g<sup>-1</sup>.

**Analyse statistique.** Les données de la première phase de l'étude (stocks de carbone des six parcelles de deux âges) ont permis d'étudier les variations de stock de carbone et d'azote par MTE pour RT et pour RN des palmeraies des deux classes d'âges. Les données de la deuxième phase ont permis de comparer les propriétés physico-chimiques des sols pour RT et RN dans les parcelles adultes.

Dans les deux phases de l'étude, l'effet sur les propriétés du sol des facteurs étudiés (âge, mode de gestion des feuilles et profondeur), les différences entre plantations d'une même classe d'âge et la variabilité aléatoire liée à chaque niveau d'observation (arbre, fosse et prélèvement), ont été évalués et testés en appliquant un Modèle Linéaire Mixte (PROC MIXED, SAS Institute 2011).

Pour la première phase ( étude du stock de C et N en fonction de l'âge et du mode de gestion), le modèle a intégré les variabilités aléatoires liées à l'effet de chaque arbre dans la parcelle de plantation et à l'interaction « parcelle x mode de restitution ». Les résidus calculés pour chacun des deux stocks mesurés autour d'un même arbre (l'un dans l'andain et l'autre
dans l'interligne) ont été considérés comme des erreurs
aléatoires indépendantes, conditionnellement à l'effet
de l'arbre.

Pour la deuxième phase, l'effet aléatoire de chaque
arbre a été inclus dans le modèle et les résidus calculés
pour chacun des prélèvements dans une fosse ont été
modélisés par un processus autorégressif d'ordre 1 qui
introduit une corrélation entre les mesures faites à des
niveaux contigus (Littell et al., 2006).

Les interactions entre facteurs étudiés ont été
détaillées en testant l’effet de chaque facteur par
tranches de combinaison des niveaux des autres
facteurs (tests connus sous le nom de « tests of simple
effects » [SAS Institute, 2011]). Des comparaisons
multiples de moyennes par le test de Tukey ont
permis d’identifier les différences significatives entre
paires. Les comparaisons ont été faites entre andain et
interligne et par profondeur.

3. RÉSULTATS

**3.1. Stocks de carbone et d'azote dans le sol**

Les stocks de carbone et d'azote pour 650 kg·m⁻² de terre (environ 50 cm de profondeur de sol) sont présentés pour chaque modalité de restitution et classe d'âge (Figure 2). Les stocks de carbone (58,0 Mg·ha⁻¹) et d'azote (4,3 Mg·ha⁻¹) obtenus avec RT en palmeraies adultes sont significativement (p < 0,001) supérieurs aux stocks obtenus pour toutes les autres situations dont les valeurs varient autour d'une moyenne de 34,6 Mg·ha⁻¹ et 2,9 Mg·ha⁻¹ respectivement pour le carbone et l'azote. Il n'y a pas de différence significative au seuil de 5 % entre les stocks sous RT et RN à l'âge

pré-adulte, ni entre les stocks sous RN des deux classes
d'âges.

**3.2. Propriétés chimiques du sol des palmeraies adultes**

Pour l'ensemble des variables examinées, les teneurs et propriétés chimiques des sols décroissent avec la profondeur (Annexe 1). Au sein d'une même zone, les teneurs des horizons 0-5 cm, 5-10 cm et 10-20 cm sont fréquemment significativement différentes entre elles et presque toujours significativement supérieures à celles des horizons plus profonds. Le pH de l'horizon 0-5 cm n'est supérieur à celui des horizons plus profonds que sous RT et il ne présente aucune variation significative avec la profondeur sous RN.

La comparaison entre zones pour un même horizon
montre aussi que les teneurs en C, N et les cations sont
significativement plus élevées sous RT jusqu'à une
profondeur de 20 cm. On observe les mêmes effets
statistiques jusqu'à 20 cm pour la somme des cations
changeables (S) et la capacité d'échange cationique
(CEC).

La figure 3 illustre les variations observées avec la profondeur et par zone, pour quatre variables importantes de la fertilité des sols sous palmeraie : C total, N total, pH et K. Elle montre une décroissance des propriétés du sol avec la profondeur et confirme qu'au-delà de 20 cm, les écarts entre RT et RN sont minimes.

**3.3. Densité apparente et porosité du sol dans les palmeraies adultes**

La densité apparente augmente avec la profondeur du
sol, que ce soit avec RT ou avec RN. Pour les deux

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

**Figure 2.** Stocks de carbone et d'azote pour 650 kg·m⁻² de masse de terre équivalente (~ 0-50 cm) pour les traitements restitution totale (RT) et restitution nulle (RN) des palmeraies de 7-12 ans et 13-24 ans — *Carbon and nitrogen stocks for 650 kg·m⁻² of equivalent soil mass (~ 0-50 cm) for the total recycling (TR) and no recycling (NR) treatments of oil palm plantation of 7-12 and 13-24 year-old.*

Les moyennes sont présentées avec les barres d'erreur standard — *the mean values are presented with standard error bars* ; les lettres a et b indiquent le classement des moyennes selon le test de Tukey (*p* < 5 %) — *the letters a and b indicate the classification of the mean values according to Tukey test* (*p* < 5%).

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Illustration de l'effet des modalités de recyclage des feuilles d'élagage (RT, restitution totale ; RN, restitution nulle) sur les paramètres physico-chimiques du sol — *Illustration of the effect of the recycling modes of pruned oil palm fronds (RT, total recycling; RN, no recycling) on the soil chemical parameters.*</figcaption></figure>

Chaque moyenne est présentée avec l'intervalle de confiance à 95 % — each mean is presented with the confidence interval to 95%.

modalités, les écarts entre deux horizons de surface sont significatifs entre eux et par rapport aux trois horizons profonds qui restent confondus (Figure 4). Les valeurs de densité apparente obtenues pour RT sont systématiquement inférieures à celles de RN et les écarts sont significatifs sur toute la profondeur étudiée.

La porosité totale du sol diminue avec la profondeur du sol. On obtient évidemment les mêmes différences entre horizons que pour la densité du sol (Tableau 3). Tout au long du profil, la porosité du sol sous RT est significativement supérieure à celle de RN. Les différences sont importantes sur les deux premiers horizons (respectivement 17 et 9 %), puis elles diminuent en profondeur où elles varient peu, autour de 4 % en moyenne.

## 4. DISCUSSION

Dans cette étude qui vise à évaluer l'effet du recyclage des feuilles élaguées sur la matière organique et sur

la fertilité du sol dans les palmeraies villageoises, on constate qu'après dix années de recyclage, les stocks de C et N sont significativement plus élevés sous RT par rapport à RN. Ils ne le sont pas pour les palmeraies pré-adultes où les feuilles sont recyclées pendant seulement une moyenne de trois années. On constate aussi que les stocks de C et N ont augmenté respectivement de 65 % et 47 % entre les deux classes d'âge pour RT, mais seulement de 16 % et 11 % entre les deux classes d'âge pour RN. On peut donc penser que c'est bien la durée d'accumulation des feuilles au sol qui conduit à augmenter progressivement le stock de matière organique du sol (MOS) jusqu'à ce que la différence entre RT et RN devienne significative. Ceci est cohérent avec les observations de Frazão et al. (2013) qui n'avaient pas trouvé de différence significative entre les stocks de carbone sous RT (27,4 Mg C·ha⁻¹) et sous RN (27,0 Mg C·ha⁻¹) sur une profondeur de 30 cm de sol dans des palmeraies brésiliennes de 8 ans, encore trop jeunes. Nos résultats confirment l'existence d'une variabilité à la fois spatiale et temporelle des stocks de
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Évolution de la densité apparente du sol sous restitution totale (RT) et restitution nulle (RN) jusqu'à 50 cm de profondeur du sol dans les palmeraies adultes (13-24 ans) — Evolution of the soil bulk density under total recycling (RT) and no recycling (RN) for 0-50 cm soil layer in adult palm (13-24 years).</figcaption></figure>

Les barres d'erreurs sont placées au niveau de chaque point pour indiquer l'intervalle de confiance à 95% — Each mean is presented with the confidence interval to 95%; les lettres indiquent la comparaison des moyennes par le test de Tukey pour les modalités RT (lettres romaines) et RN (lettres grecques) — the letters indicate the comparison of means by Tukey test for TR (Roman letters) and NR (Greek letters) modes.

MOS, qui avait été signalée par d'autres auteurs (Djègui, 1992; Law et al., 2009). Frazão et al. (2013) ont aussi montré des différences de stocks de carbone entre RT (29,6 Mg·ha⁻¹) et RN (25,4 Mg·ha⁻¹) dans des palmeraies de 25 ans au Brésil. Ils ont aussi constaté en 2014 que le stock de carbone sous andains de feuilles d'élagage d'une palmeraie de 34 ans était supérieur de 78 % à celui d'une plantation de 23 ans (Frazão et al., 2014).

En considérant que RT ne représente que 10 % de la surface agricole, on peut calculer à partir de la masse de terre équivalente (MTE) de 650 kg·m⁻², le stock moyen de C d'une palmeraie adulte où les feuilles d'élagage sont recyclées. On obtient une valeur de 39 Mg·ha⁻¹ pour une profondeur de sol d'environ 50 cm. Le gain est de

<table><caption>**Tableau 3.** Porosité moyenne du sol (%) par profondeur sous restitution totale (RT) et restitution nulle (RN) dans les palmeraies adultes — Soil porosity (%) for the soil layers under total recycling (RT) and no recycling (RN) in adult oil palm plantations.</caption><thead><tr><th>Profondeur (cm)</th><th>Porosité sol RT</th><th>Porosité sol RN</th></tr></thead><tbody><tr><td>0-5</td><td>66 ± 4<sup>a, α</sup></td><td>49 ± 4<sup>a, β</sup></td></tr><tr><td>5-10</td><td>52 ± 7<sup>b, α</sup></td><td>43 ± 4<sup>b, β</sup></td></tr><tr><td>10-20</td><td>44 ± 3<sup>c, α</sup></td><td>40 ± 2<sup>c, β</sup></td></tr><tr><td>20-30</td><td>44 ± 2<sup>c, α</sup></td><td>41 ± 2<sup>c, β</sup></td></tr><tr><td>30-50</td><td>45 ± 2<sup>c, α</sup></td><td>41 ± 2<sup>c, β</sup></td></tr></tbody></table>

Les valeurs moyennes de porosité sont présentées ± leur écart-type — mean values of soil porosity averages are presented ± the standard deviation; les lettres en exposant de chaque valeur indiquent la comparaison des moyennes selon le test de Tukey — superscript letters indicate the comparison of the means by the Tukey test; les lettres romaines comparent les moyennes des profondeurs dans un même mode de gestion des feuilles d'élagage, alors que les lettres grecques comparent les moyennes d'une même profondeur entre les deux modes de gestion — Roman letters compare the mean values of the same mode of pruned fronds management, while the Greek letters compare the same soil layer between the two modes.

2,1 Mg·ha⁻¹, soit 5,7 % par rapport à une palmeraie du même âge où l'on ne recyclerait pas les feuilles. Un calcul analogue pour l'azote conduit à un stock total de 3,21 Mg·ha⁻¹ dans une palmeraie RT, soit un gain de 0,12 Mg·ha⁻¹, qui représente 4 % par rapport à une palmeraie sans recyclage des feuilles. Notre stock de carbone est cohérent avec les 34,3 Mg·ha⁻¹ indiqués par Frazão et al. (2014) dans les palmeraies industrielles de 23 ans au Brésil, pour une profondeur de 30 cm du sol.

Nos résultats confirment que les palmeraies villageoises du Bénin seraient, à l'instar d'autres cultures, un « puits de carbone » que le recyclage de la biomasse peut alimenter (Arrouays et al., 2002; Bernoux et al., 2011). L'augmentation des stocks de carbone dans les palmeraies a été identifiée par Carlson et al. (2012) comme un moyen de mitigation des émissions de gaz à effet de serre dans certaines conditions de précédents et de pratiques culturaux.

Notre estimation pour une palmeraie adulte est proche du stock de carbone (42 Mg·ha⁻¹) trouvé par Barthès et al. (2004) sur 0-40 cm de profondeur du sol dans le système agricole maïs-mucuna installé sur une terre de barre du sud-Bénin, avec restitution des résidus du mucuna pendant 10 ans. Ce qui confirme que le potentiel de stockage des palmeraies est aussi important que d'autres cultures. Notre estimation est toutefois inférieure à celles mesurées dans des forêts secondaires par Djegui et al. (1992) au Bénin (48 Mg·C·ha⁻¹ sur 0-35 cm) et par Lardy Chapuis et al. (2002) au Brésil (53 Mg·C·ha⁻¹ sur 0-20 cm).
Cette réflexion indique que l'intérêt des palmeraies pour la séquestration du carbone doit donc être modulé en fonction du précédent cultural. Reinhardt et al. (2007) ont montré que lorsque les palmeraies ne sont pas directement plantées après une forêt primaire, le stock de carbone du sol augmente avec le temps en raison du recyclage de la biomasse aérienne. Les palmeraies béninoises dans leur grande majorité sont implantées sur d'anciennes cultures annuelles (Aholoukpè et al., 2013) qui ont contribué à l'appauvrissement en carbone du sol (Amadji et al., 2008). Elles peuvent donc être très efficaces pour restaurer les stocks de C lorsque les feuilles ne sont pas exportées.

L'augmentation des stocks de matière organique du sol avec RT en palmeraies adultes a eu des effets positifs sur les propriétés physico-chimiques du sol. Sur les 20 premiers cm, les écarts entre RT et RN pour la porosité du sol, la CEC, les cations et les teneurs en N, sont significatifs. On a noté aussi une amélioration de 0,7 unité de pH sur l'horizon 0-5 cm. Ces résultats confirment les effets positifs signalés par Goh et al. (2003) en Asie du Sud-Est et par Nelson et al. (2011) en Papouasie-Nouvelle-Guinée avec le recyclage des rafles en replantation : amélioration de la CEC, des cations échangeables et du pH du sol.

Ces changements dus au recyclage des feuilles sont d'importance puisqu'ils touchent des propriétés qui conditionnent directement les facteurs de productivité du palmier à huile :

- N et K sont deux éléments essentiels dans les processus métaboliques qui contrôlent la production de matière sèche (N : assimilation chlorophyllienne, K : effet sur la photosynthèse *via* le contrôle de l'ouverture stomatique, transport des carbohydrates vers les régimes) (Nodichao, 2008; Dubos et al., 2011);

- les réserves hydriques du sol sont améliorées par l'augmentation de la porosité concomitante à l'augmentation des teneurs en MOS, celles de Ca et Mg ainsi que la CEC qui favorisent l'obtention d'un complexe argilo-humique stable conférant au sol une meilleure structure.

Ces résultats montrent que le choix que nous avons fait pour une MTE de 650 kg·m⁻<fnref n="2" /> est donc justifié puisqu'en travaillant sur une profondeur qui n'excédait pas 50 cm, nous avons largement couvert la couche de sol où se produisent les changements majeurs induits par le recyclage des feuilles.

Ces effets intéressent certes un volume de sol limité. Cependant, sous l'andain, à mi-distance entre deux palmiers, l'essentiel de la biomasse des racines est située dans les 30 premiers cm du sol et en particulier les racines fines absorbantes sont situées dans les 10 premiers cm (Nodichao et al., 2008). La surface traitée ne représente que 10 % de la surface colonisée

par les racines d'une palmeraie adulte, mais le progrès mérite d'être souligné si on compare les propriétés du sol RN aux normes qui prévalent dans la profession pour l'horizon de surface : Goh et al.<fnref n="2" /> (1997), cités par Paramananthan (2003), considèrent en effet que des limitations du rendement sont attendues lorsque les seuils suivants ne sont pas atteints : C : 15 g·kg⁻¹, N : 1,5 g·kg⁻¹, K et Mg : 0,25 meq·100 g⁻¹. Les teneurs en C, N et K des 20 premiers centimètres du sol de la modalité RN sont donc très faibles en comparaison avec ces normes. La CEC est aussi très faible si on compare avec les normes optimales données par Paramananthan (2003) (18 meq·100 g⁻¹) ou Nelson et al. (2010).

Compte tenu de ces considérations, dans le contexte de plantations dans des conditions de faible fertilité, il est donc prudent de ne pas se priver d'une surface, même réduite comme celle de recyclage des feuilles, où la fertilité du sol est beaucoup plus proche des conditions minimales recommandées pour la culture. Nous n'avons toutefois pas pu corroborer par des observations de production, réalisées en nombre insuffisant (données non présentées), que le recyclage des feuilles aurait un impact direct sur les rendements des plantations de notre étude.

## 5. CONCLUSIONS

Le recyclage des feuilles d'élagage induit au fil des années des modifications des propriétés des sols sous les zones de restitution. Nous avons pu en mesurer les effets significatifs sur les stocks de C et N après 10 ans de recyclage dans les palmeraies adultes. L'augmentation de la MOS améliore aussi les propriétés physico-chimiques du sol par rapport aux zones non traitées. Les effets positifs sur les teneurs en éléments nutritifs (N, K, Ca, Mg), la porosité du sol et sa CEC sont significatifs jusqu'à une profondeur de 20 cm, or cette couche de sol est très colonisée par les racines. Même si la surface couverte par les andains ne représente qu'environ 10% de la surface cultivée, nous pensons qu'elle peut jouer un rôle favorable pour la culture, ne serait-ce qu'en raison des conditions de fertilité très médiocres que nous avons trouvées pour le sol hors recyclage. On peut conclure qu'il est recommandé de ne pas exporter les feuilles d'élagage, d'autant plus que l'espace amendé sera utilisé pour planter la prochaine génération de palmiers après 25 ans environ.

## Bibliographie

Adjadi E., 2012. *Le développement du palmier à huile sélectionné au Bénin*. Rapport d'étude. Actualisation des

<footnote n="2">² Goh K.J. et al., 1997. Developing decision support systems for plantations: a new role for agronomists. *The Planter*, **73**, 687-693.</footnote>
superficies plantées. Pobè, Bénin : Centre de Recherches Agricoles Plantes Pérennes (CRAPP).

Adjè A.I., 2004. *Contribution à l'analyse du micro-climat, de la fertilité des sols et de la productivité des systèmes de culture à base du palmier à huile* (Elaeis guineensis Jacq) au Sud-Bénin : cas des sites expérimentaux de Pobè et d'Obèkè-Ouèrè. Thèse de doctorat: Université Félix Houphouët Boigny, Abidjan-Cocody (Côte d'Ivoire).

Aho N. & Kossou D.K., 1997. *Précis d'agriculture tropicale*. Cotonou, Bénin : Les Éditions du Flamboyant.

Aholoukpè N.S.H., 2013. *Matière organique du sol et développement du palmier à huile sous différents modes de gestion des feuilles d'élagage*. *Cas des palmeraiès villageoises du département du Plateau au Bénin*. Thèse de doctorat : Université d'Abomey-Calavi (Bénin).

Aholoukpè H. et al., 2013. Typologie des plantations villageoises de palmier à huile (Elaeis guineensis Jacq.) dans le département du Plateau au Bénin. *Int. J. Biol. Chem. Sci.*, **7**, 978-999.

Ahoyo N.R., 2008. *Potentialités des filières anacarde et palmier à huile au Bénin*. *Rapport d'étude*. Cotonou, Bénin : SNV Bénin.

Akpo E. et al., 2012. A participatory diagnostic study of oil palm (Elaeis guineensis) seed system in Benin. *Wageningen J. Life Sci.*, **60-63**, 15-27.

Amadji G.L. & Aholoukpè H.N.S., 2008. Impact du niébé (*Vigna unguiculata*) et de la fumure minérale sur les propriétés chimiques de la terre de barre du Bénin. *Étud. Gestion Sols*, **15**(3), 147-160.

Arrouays D. et al., éds, 2002. *Contribution à la lutte contre l'effet de serre. Stocker du carbone dans les sols agricoles de France ? Expertise Scientifique Collective*. *Synthèse du rapport*. INRA.

Azontondé H.A., 1991. Propriétés physiques et hydrauliques des sols au Bénin. *IAHS Publ.*, **199**, 253-256.

Barcelos E., 1998. *Étude de la diversité génétique du genre Elaeis* (E. oleifera Kunth. Cortès et E. guineensis Jacq.) *par marqueur moléculaires (RFLP et AFLP)*. Thèse de doctorat : École Nationale Supérieure Agronomique de Montpellier (France).

Barthes B. et al., 2004. Effect of a legume cover crop (*Mucuna pruriens* var. utilis) on soil carbon in an Ultisol under maize cultivation in southern Benin. *Soil Use Manage.*, **20**, 231-239.

Barthès B.G. et al., 2006. Determination of total carbon and nitrogen content in a range of tropical soils using near infrared spectroscopy: influence of replication and sample grinding and drying. *J. Near Infrared Spectrosc.*, **14**, 341-348.

Bernoux M. et al., 2011. Le programme GESSOL 2 : impact des pratiques agricoles sur les matières organiques et les fonctions des sols. *Étud. Gestion Sols*, **18**(3), 137-145.

Brunet D., Barthès B., Chotte J.-L. & Feller C., 2007. Determination of carbon and nitrogen contents in Alfisols, Oxisols and Ultisols from Africa and Brazil

using NIRS analysis: effects of sample grinding and set heterogeneity. *Geoderma*, **139**(1-2), 106-117.

Carlson K.M. et al., 2013. Carbon emissions from forest conversion by Kalimantan oil palm plantations. *Nature Climate Change*, **3**, 283-287.

Chang C.-W., Laird D.A., Mausbach M.J. & Hurburgh C.R. Jr., 2001. Near-infrared reflectance spectroscopy-principal components regression analyses of soil properties. *Soil Sci. Soc. Am. J.*, **65**, 480-490.

Corley R.H.V. & Tinker P.B., eds, 2003. *The oil palm*. 4th ed. Oxford, UK: Blackwell Science.

Dassou O.S., 2011. *Impact du choix de la méthode d'échantillonnage racinaire dans la détermination des paramètres de développement racinaire chez le palmier à huile* (Elaeis guineensis Jacq.) au Bénin. *Mémoire d'ingénieur agronome* : Université d'Abomey-Calavi (Bénin).

Djegui N., 1992. *Matière organique et azote dans les sols cultivés sur terres de barre (Bénin)*. Thèse de doctorat : École Nationale Supérieure d'Agronomie de Toulouse (France).

Djegui N., 1995. Les stocks organiques dans les sols cultivés sous palmeraiè et cultures vivrières dans le sud du Bénin. *In : Actes du séminaire, Fertilité du milieu et stratégies paysannes sous les tropiques humides, 13-17 novembre 1995, CIRAD, Montpellier, France*, 189-193.

Djègui N., de Boissezon P. & Gavinelli E., 1992. Statut organique d'un sol ferrallitique du Sud-Bénin sous forêt et différents systèmes de cultures. *Cah. Orstom Sér. Pédol.*, **17**(1), 5-22.

Dubos B., Alarcon H.W., Lopez J.E. & Ollivier J., 2011. Potassium uptake and storage in oil palm organs: the role of chlorine and the influence of soil characteristics in the Magdalena valley, Colombia. *Nutr. Cycl. Agroecosyst.*, **89**, 219-227.

Ellert B.H. & Bettany J.R., 1995. Calculation of organic matter and nutrients stored in soils under contrasting management regimes. *Can. J. Sci.*, **75**, 529-538.

Ellert B.H., Janzen H.H. & Entz T., 2002. Assessment of a method to measure temporal change in soil carbon storage. *Soil Sci. Soc. Am. J.*, **66**, 1687-1695.

Frazão L.A., Paustian K., Pellegrino Cerri C.E. & Cerri C.C., 2013. Soil carbon stocks and changes after oil palm introduction in the Brazilian Amazon. *Global Change Biol. Bioenergy*, **5**, 384-390.

Frazão L.A., Paustian K., Pellegrino Cerri C.E. & Cerri C.C., 2014. Soil carbon stocks under oil palm plantations in Bahia State, Brazil. *Biomass Bioenergy*, **62**, 1-7.

Goh K.J. & Härdter R., 2003. General oil palm nutrition. *In: Fairhurst T.H. & Härdter R.*, eds. *Oil palm management for large and sustainable yields*. Singapore: Potash & Phosphate Institute (PPI), Potash & Phosphate Institute of Canada (PPIC) and International Potash Institute (IPI, Basel), 191-230.

Khalid H., Zin Z.Z. & Anderson J.M., 1999. Quantification of oil palm biomass and nutrient value in a mature
plantation. 1. Above-ground biomass. *J. Oil Palm Res.*, **11**(1), 23-32.

Khalid H., Zin Z.Z. & Anderson J.M., 2000a. Soil nutrient dynamics and palm growth performance in relation to residue management practices following replanting of oil palm plantation. *J. Oil Palm Res.*, **12**(1), 25-45.

Khalid H., Zin Z.Z. & Anderson J.M., 2000b. Decomposition processes and nutrients release patterns of oil palm residues. *J. Oil Palm Res.*, **12**(1), 46-63.

Lardy Chapuis L., Brossard M., Assad Lopes M.L. & Laurent J.Y., 2002. Carbon and phosphorus stocks of clayey Ferralsols in Cerrado native and agroecosystems, Brazil. *Agric. Ecosyst. Environ.*, **92**, 147-158.

Law M.C. et al., 2009. Spatial variability of soil organic carbon in oil palm: a comparison between young and mature stands. *Int. J. Agric. Res.*, **4**(12), 402-417.

Littell R.C. et al., 2006. *SAS for mixed models*. 2nd ed. Cary, NC, USA: SAS Institute Inc.

Nelson N.P. & Su N., 2010. Soil pH buffering capacity: a descriptive function and its application to some acidic tropical soils. *Aust. J. Soil Res.*, **48**, 201-207.

Nelson N.P. et al., 2011. Soil acidification under oil palm: rates and effects on yield. *Better Crops*, **95**(4), 22-25.

Nodichao L., 2008. *Biodiversité racinaire, absorption potassique et résistance à la sécheresse chez le palmier à huile au Bénin*. Thèse de doctorat : Université Félix Houphouët Boigny, Abidjan-Cocody (Côte d'Ivoire).

Nodichao L., Aké S. & Jourdan C., 2008. Développement du système racinaire chez le palmier à huile selon

l'origine génétique et le régime hydropotassique du sol. *Agron. Afr.*, **20**(3), 277-289.

Nouy B., Baudouin L., Djégui N. & Omoré A., 1999. Le palmier à huile en conditions hydriques limitantes. *Plant. Rech. Dev.*, **6**(1), 31-45.

Pansu M. & Gautheyrou J., 2006. *Handbook of soil analysis. Mineralogical, organic and inorganic methods*. Berlin; Heidelberg, Germany: Springer-Verlag.

Paramananthan, 2003. Land selection for oil palm. In: Fairhurst T. & Härdter R., eds. *Oil palm management for large and sustainable yields*. Singapore: Potash & Phosphate Institute (PPI), Potash & Phosphate Institute of Canada (PPIC) and International Potash Institute (IPI, Basel), 27-57.

Reinhardt G., Retteinmaier N. & von Falkenstein E., 2007. *Conclusive evaluation of studies. Assessing the environmental impact of the use of palm oil as a bioenergy carrier. Final Report*. Heidelberg, Germany: Institut für Energie- und Umweltforschung (IFEU).

SAS Institute Inc, 2011. *SAS/STAT® 9.3 User's Guide*. Cary, NC, USA: SAS Institute Inc.

Volkoff B. & Willaime P., 1976. *Carte pédologique de reconnaissance de la République Populaire du Bénin. Feuille de Porto-Novo. Notice explicative n°66(1)*. Paris : ORSTOM, 39, carte h.t.

(44 réf.)

<table><caption>**Annexe 1.** Paramètres physico-chimiques du sol sous RT (andain) et RN (interligne) des palmeraies adultes — *Soil chemical parameters under TR (pruned fronds pile) and NR (interrow) of the adult oil palm plantations*.</caption><thead><tr><th rowspan="2">Zone</th><th rowspan="2">Profondeur (cm)</th><th colspan="3">g·kg<sup>-1</sup></th><th rowspan="2">pH H<sub>2</sub>O</th><th colspan="5">meq·100 g<sup>-1</sup></th></tr><tr><th>C</th><th>N</th><th>C/N</th><th>Ca</th><th>Mg</th><th>K</th><th>S</th><th>CEC</th></tr></thead><tbody><tr><td rowspan="5">RT (Andain)</td><td>0-5</td><td>31,28<sup>a, α</sup></td><td>2,07<sup>a, α</sup></td><td>15,22<sup>a, α</sup></td><td>6,2<sup>a, α</sup></td><td>7,32<sup>a, α</sup></td><td>3,42<sup>a, α</sup></td><td>0,17<sup>a, α</sup></td><td>10,90<sup>a, α</sup></td><td>9,00<sup>a, α</sup></td></tr><tr><td>5-10</td><td>19,29<sup>b, α</sup></td><td>1,26<sup>b, α</sup></td><td>15,49<sup>a, α</sup></td><td>6,0<sup>b, α</sup></td><td>4,42<sup>b, α</sup></td><td>1,70<sup>b, α</sup></td><td>0,09<sup>b, α</sup></td><td>6,21<sup>b, α</sup></td><td>5,85<sup>b, α</sup></td></tr><tr><td>10-20</td><td>10,12<sup>c, α</sup></td><td>0,69<sup>c, α</sup></td><td>14,62<sup>a, α</sup></td><td>5,9<sup>b, α</sup></td><td>2,84<sup>c, α</sup></td><td>0,98<sup>c, α</sup></td><td>0,05<sup>c, α</sup></td><td>3,87<sup>c, α</sup></td><td>4,05<sup>c, α</sup></td></tr><tr><td>20-30</td><td>4,98<sup>d, α</sup></td><td>0,43<sup>d, α</sup></td><td>11,49<sup>b, α</sup></td><td>5,8<sup>b, α</sup></td><td>2,01<sup>d, α</sup></td><td>0,74<sup>c, α</sup></td><td>0,03<sup>d, α</sup></td><td>2,78<sup>d, α</sup></td><td>3,40<sup>c, α</sup></td></tr><tr><td>30-50</td><td>3,55<sup>d, α</sup></td><td>0,36<sup>d, α</sup></td><td>9,65<sup>c, α</sup></td><td>5,7<sup>b, α</sup></td><td>2,03<sup>d, α</sup></td><td>0,76<sup>c, α</sup></td><td>0,02<sup>d, α</sup></td><td>2,82<sup>d, α</sup></td><td>3,80<sup>c, β</sup></td></tr><tr><td rowspan="5">RN (Interligne)</td><td>0-5</td><td>12,08<sup>a, β</sup></td><td>0,84<sup>a, β</sup></td><td>14,48<sup>a, α</sup></td><td>5,5<sup>a, β</sup></td><td>2,68<sup>a, β</sup></td><td>0,87<sup>a, β</sup></td><td>0,06<sup>a, β</sup></td><td>3,62<sup>a, β</sup></td><td>4,71<sup>a, β</sup></td></tr><tr><td>5-10</td><td>9,04<sup>b, β</sup></td><td>0,64<sup>b, β</sup></td><td>13,92<sup>a, β</sup></td><td>5,6<sup>a, β</sup></td><td>2,40<sup>b, β</sup></td><td>0,70<sup>b, β</sup></td><td>0,04<sup>b, β</sup></td><td>3,14<sup>b, β</sup></td><td>3,81<sup>b, β</sup></td></tr><tr><td>10-20</td><td>6,19<sup>c, β</sup></td><td>0,49<sup>c, β</sup></td><td>12,45<sup>b, β</sup></td><td>5,7<sup>a, β</sup></td><td>2,20<sup>b, β</sup></td><td>0,64<sup>b, β</sup></td><td>0,03<sup>c, β</sup></td><td>2,86<sup>b, β</sup></td><td>3,50<sup>c, β</sup></td></tr><tr><td>20-30</td><td>3,68<sup>d, α</sup></td><td>0,36<sup>d, α</sup></td><td>9,89<sup>c, β</sup></td><td>5,7<sup>a, α</sup></td><td>1,84<sup>c, α</sup></td><td>0,68<sup>b, α</sup></td><td>0,02<sup>d, α</sup></td><td>2,54<sup>c, α</sup></td><td>3,35<sup>c, α</sup></td></tr><tr><td>30-50</td><td>3,41<sup>d, α</sup></td><td>0,37<sup>d, α</sup></td><td>8,63<sup>c, β</sup></td><td>5,7<sup>a, α</sup></td><td>2,22<sup>b, α</sup></td><td>0,86<sup>a, α</sup></td><td>0,02<sup>c, d, α</sup></td><td>3,10<sup>b, α</sup></td><td>4,38<sup>d, α</sup></td></tr></tbody></table>

S : somme cations — *som of cations*; CEC : capacité d'échange cationique — *cationic exchange capacity*; les lettres en exposant de chaque valeur indiquent la comparaison des moyennes selon le test de Tukey — *superscript letters indicate the comparison of the means by the Tukey test*; les lettres romaines comparent les moyennes des profondeurs dans le même mode de gestion des feuilles d'élageage, alors que les lettres grecques comparent les moyennes d'une même profondeur entre les deux modes — *Roman letters compare the soil layers of the same mode of pruned fronds management, while the Greek letters compare the same soil layer between the two modes.*