Résumé

Cet article présente une revue bibliographique relative au développement de la spectroscopie proche infrarouge appliquée à l’analyse de terre en vue de contribuer à l’évaluation de la fertilité chimique des sols. Cette technique permet d’estimer la composition chimique des échantillons de terre grâce à leurs propriétés d’absorption de la lumière. Il s’agit donc d’une méthode de mesure indirecte nécessitant une phase d’étalonnage préalable à la prédiction des propriétés étudiées. Elle présente de nombreux atouts en comparaison aux analyses réalisées selon les méthodes de référence : rapidité, facilité de mise en œuvre sans consommation d’extractifs, etc. Très courante dans certains domaines, comme l’agro-alimentaire, elle a cependant tardé à être utilisée en analyse de terre en raison de la complexité de la matrice sol. Depuis une dizaine d’années, grâce notamment au développement des techniques chimiométriques, de très nombreuses recherches ont été effectuées pour évaluer son potentiel. La majorité d’entre elles ont conclu à la faisabilité de la prédiction pour un certain nombre de propriétés de sol, mais cette technique reste toutefois peu utilisée en routine à l’heure actuelle. Après une présentation de la technique et des méthodes d’analyses statistiques s’y rapportant, les caractéristiques du sol qui peuvent être prédites par cette technique sont présentées. Enfin, quelques exemples de résultats obtenus en analyse de terre sont détaillés, principalement pour la détermination du taux d’argile, du taux de carbone organique et de la capacité d’échange cationique.

Abstract

Use of near-infrared reflectance spectroscopy in soil analysis. A review. This paper presents a literature review on the development of near infrared reflectance spectroscopy for soil analysis and the contribution of this technique to the evaluation of soil fertility analysis. This technique is used to estimate the chemical composition of soil samples on the basis of their absorption properties. It is therefore an indirect method of measurement, which requires a calibration phase for the prediction of these properties. NIR spectroscopy offers many advantages compared to reference analysis: it is known to be a physical, non-destructive, rapid, reproducible and low cost method. Often employed in other analytical domains, such as agro-food, NIR spectroscopy has, however, seldom been used in soil characterization, due to the complexity of the soil matrix. Thanks to the development of chemometric methods, numerous studies have recently been conducted to evaluate the feasibility of the application of the technique in soil analysis. Most authors conclude that NIR spectroscopy is promising; however, to date, use of the technique has not spread to routine laboratories. The paper is organized as follows. Firstly, we provide an overview of the NIR spectroscopy technique and related chemometric methods. Secondly, we describe the soil characteristics that can be predicted using this technique. Finally, we detail examples of results that have been obtained through the use of the technique, mainly in the determination of clay and organic carbon content, and of cation exchange capacity.

1. Introduction

1La fertilité d’un sol peut être définie comme son aptitude à produire de la biomasse végétale, le sol étant alors le support de la production. Il est fréquemment visualisé comme une boite noire dans laquelle on introduit des éléments fertilisants et dont on extrait des produits. Mais le sol est le siège d’une série de processus intervenant au sein de systèmes (sol – eau – plante – atmosphère) dynamiques et dont la compréhension est indispensable pour contribuer à une gestion qui concilie objectifs de production et durabilité environnementale (Genot et al., 2009). L’évaluation de la fertilité chimique passe par la réalisation d’une analyse en laboratoire, domaine dans lequel la spectroscopie proche infrarouge (SPIR ou NIRS pour Near InfraRed Spectroscopy) offre des perspectives intéressantes d’autant que, vu sa large utilisation dans le domaine agro-alimentaire, les laboratoires d’analyse de terre disposent très souvent d’un spectromètre. En effet, cette technique ne nécessite qu’une préparation minimale de l’échantillon, un faible besoin en matériel, une formation réduite des manipulateurs et offre une grande rapidité de mesure ainsi que la détermination simultanée de plusieurs composants (Russel et al., 2002). De plus, aucun extractif chimique n’est utilisé et des mesures in situ de propriétés physico-chimiques des sols sont également envisageables par cette technique (Viscarra Rossel et al., 2006), ce qui ouvre des perspectives en matière d’agriculture de précision. Enfin, la mesure de certaines propriétés chimiques ou biologiques des sols qui ne se réalise pas en routine actuellement, car les méthodes sont trop couteuses malgré leur intérêt dans l’évaluation de la fertilité du sol, pourrait être remplacée par une prédiction sur base de la spectroscopie proche infrarouge (Stenberg et al., 2010).

2La SPIR s’est tout d’abord très largement diffusée dans le monde de l’agro-alimentaire. Hart et al. (1962) et Massie et al. (1965) ont été parmi les pionniers, en utilisant la technique pour des analyses sur les céréales et les fruits. En parallèle, le développement de méthodes chimiométriques et de processeurs de plus en plus performants ont également grandement contribué à son succès. Elle a cependant tardé à être utilisée en analyse de terre compte tenu du caractère complexe du sol. Islam et al. (2003) décrivent le sol comme un système non idéal chimiquement et minéralogiquement, qui est plus complexe que les systèmes purs fréquemment étudiés dans les laboratoires. Bellon Maurel et al. (2010) justifient cette complexité par le fait que les échantillons de terre, contrairement aux produits issus du domaine agro-alimentaire, n’ont pas de code génétique avec, pour conséquence, un comportement beaucoup plus variable, des distributions asymétriques, etc. Les mécanismes du sol sont ainsi malaisés à comprendre complètement et les liens fondamentaux entre les mesures chimiques et les propriétés du sol sont complexes. Par conséquent, il est difficile de mettre directement en relation une caractéristique physique, chimique ou physico-chimique du sol avec une longueur d’onde particulière. Shepherd et al. (2002) indiquent ainsi que, dans l’objectif de réaliser un étalonnage des propriétés des sols, il est préférable d’utiliser l’information du spectre en entier, plutôt que de tenter d’interpréter les bandes d’absorption individuellement.

3Les premières publications sur les potentialités de la SPIR en analyse de terre apparaissent au début des années 1990 avec les travaux de Sudduth et al. (1991) et de Ben-Dor et al. (1995). Elles ont ensuite connu une croissance exponentielle au cours des années 2000. Comme l’ont montré de nombreux auteurs (Dalal et al., 1986 ; Morra et al., 1991 ; Ben-Dor et al., 1995 ; Chang et al., 2001 ; Hummel et al., 2001 ; Chang et al., 2002 ; Fystro, 2002 ; Ludwig et al., 2002 ; Malley et al., 2002 ; Mc Carty et al., 2002 ; Reeves et al., 2002 ; Shepherd et al., 2002 ; Islam et al., 2003 ; Udelhoven et al., 2003 ; Van Groeningen et al., 2003 ; He et al., 2005 ; Brown et al., 2006 ; Nanni et al., 2006 ; Viscarra Rossel et al., 2006 ; Brunet et al., 2007 ; Terhoeven-Urselmans et al., 2008 ; Genot et al., 2011a ; etc.), la technique SPIR appliquée à l’analyse de terre présente ainsi de sérieux atouts en vue d’une utilisation en routine, en remplacement ou en parallèle aux protocoles classiques. Très présente dans le monde de la recherche, elle a cependant encore très peu émergé en routine dans les laboratoires d’analyse de terre ou en technique embarquée pour des mesures in situ. Actuellement, les laboratoires d’analyse de terre membres du réseau REQUASUD en Belgique (Genot et al., 2011b) utilisent cette technique pour prédire certaines propriétés de sol comme le carbone organique, la CEC, l’azote total, le taux d’humidité et le taux d’argile d’échantillons de terre en remplacement ou en complément des méthodes analytiques de référence.

4L’objectif de cette revue bibliographique est de présenter un état des recherches menées sur l’application de la SPIR pour la prédiction de propriétés des sols indicatrices de leur fertilité et répondre ainsi aux exigences d’une agriculture performante et respectueuse de l’environnement. Dans un premier temps, la technique SPIR et les méthodes chimiométriques associées seront brièvement présentées. Ensuite, nous présenterons les principales propriétés du sol qui influencent le spectre infrarouge. Nous présenterons également les principaux résultats obtenus pour la prédiction du taux de carbone organique et de la CEC. Enfin, nous discuterons des paramètres importants à considérer pour évaluer les potentialités de la technique en analyses de terre et envisager une utilisation en routine.

2. Fondements de la SPIR et analyse statistique multi-variée

2.1. SPIR et analyse spectrale

5L’analyse par spectroscopie dans le proche infrarouge utilise les propriétés du rayonnement proche infrarouge et permet d’établir une relation entre la réponse de l’échantillon à ce rayonnement et les caractéristiques de l’échantillon qu’on souhaite étudier.

6Le rayonnement infrarouge est une radiation de nature électromagnétique de longueur d’onde supérieure à celle de la lumière visible et plus courte que celle des micro-ondes (Figure 1). Ce rayonnement est un type d’énergie qui se propage à travers l’espace et la matière et dont la plupart des propriétés peuvent être décrites à l’aide d’un modèle dans lequel il est assimilé à une onde sinusoïdale. Le rayonnement infrarouge est traditionnellement subdivisé en trois régions : les domaines du proche, du moyen et du lointain infrarouge. En proche infrarouge, la gamme de longueurs d’onde retenue varie généralement entre 700 à 780 et 2 500 nm, une distinction entre le très proche infrarouge et le proche infrarouge est parfois faite, la limite entre ces deux domaines se situant entre 900 et 1 000 nm.

7La technique consiste à soumettre un échantillon à un rayonnement situé dans le domaine du proche infrarouge, l’échantillon est ainsi scanné. Pour chaque longueur d’onde, la part du rayonnement réfléchi (réflectance, R) par l’échantillon est mesurée à l’aide de détecteurs et convertie en absorbance (Log1/R). L’ensemble de ces absorbances constituent le spectre qui peut être considéré comme une empreinte globale, reflétant la composition chimique de l’objet analysé.

8Lorsqu’un rayonnement atteint une couche de matière, certaines de ses fréquences peuvent être sélectivement absorbées, l’énergie électromagnétique correspondante est transférée aux atomes, ions ou molécules constituant la matière. Ce transfert d’énergie s’accomplit par absorption de photons et les atomes ou molécules passent d’un état fondamental à un état énergétique supérieur qualifié d’état excité (Lachenal, 2006). Quelle que soit sa fréquence, un rayonnement proche infrarouge n’est pas suffisamment énergétique pour provoquer des modifications des états d’énergie des électrons et des particules. Il induit uniquement des changements d’état d’énergie vibratoire ou rotatoire. Les bandes d’absorption observées dans le proche infrarouge sont ainsi des harmoniques ou des combinaisons de bandes fondamentales observées dans l’infrarouge moyen (Prévot, 2004). Elles sont plus larges et contiennent moins d’informations que les spectres obtenus dans le moyen infrarouge et peuvent donc être plus difficiles à interpréter. Néanmoins, elles contiennent des informations utiles sur les constituants organiques et inorganiques caractéristiques du sol. Comme illustré par la figure 2, les bandes d’absorption observées dans la région du proche infrarouge résultent des harmoniques des groupements OH, SO₄ et CO₃, ainsi que des combinaisons des bandes fondamentales de H₂O et CO₂ (Clark, 1999). Cette figure (©ASD Inc.) se base sur les recherches menées par Ellis (1928), Kaye (1954), Goddu (1960), Goddu et al. (1960), Weyer et al. (2002) et Workman (2000). Les minéraux argileux présentent également des bandes d’absorption dans la région du proche infrarouge en relation avec des liaisons OH-métal et des harmoniques des liaisons O-H (Viscarra Rossel et al., 2006). Les carbonates ont également des pics d'absorption faible dans le domaine du proche infrarouge (Hunt et al., 1970). Enfin, l'eau influence fortement le spectre avec des bandes d'absorption dominantes autour de 1 450 et 1 940 nm et des bandes plus faibles dans d'autres parties du domaine proche infrarouge (Liu et al., 2002).

2.2. Traitement du spectre

9Compte tenu de la complexité du spectre, il est toutefois difficile de tenter d’interpréter les bandes une à une pour prédire un paramètre quantitatif. Des approches statistiques d’analyses multi-variées ou chimiométriques sont développées pour progresser dans l’exploitation des spectres SPIR. La chimiométrie peut être définie comme l’ensemble des méthodes statistiques, graphiques ou symboliques permettant de mieux comprendre les informations obtenues dans le domaine de la chimie (Bertrand et al., 2006). Les spectres sont le résultat d’une interaction de la lumière avec la matière qui ne peut être complètement décrite sur le plan théorique. Les lacunes de la connaissance théorique rendent indispensable l’approche expérimentale. Pour tenir compte des variations incontrôlées inhérentes à l’expérimentation, il est nécessaire d’accumuler les spectres de nombreux échantillons et de traiter la collection spectrale ainsi obtenue par des méthodes chimiométriques (Bertrand et al., 2006). C’est l’étape d’étalonnage des données qui est indispensable avant une utilisation en routine de la SPIR. Les analyses multi-variées s’appliquent à une collection homogène de spectres qui sert à l’établissement d’un modèle. La constitution, au sein de la population ou du territoire étudié, d’une collection d’échantillons représentative de la diversité de cette population ou de ce territoire est une étape clé de la réussite de la phase d’étalonnage des données. Un modèle prédictif basé sur un lot d’échantillons non représentatifs ne permettra pas de prédire de nouveaux échantillons parce qu’ils seront trop éloignés de la collection ayant servi à l’étalonnage. Au contraire, des modèles élaborés sur base d’une collection représentative de la diversité rencontrée permettront de prédire correctement des échantillons faisant partie de la même population ou du même territoire. Un spectre est, par nature, un signal continu, mais il est transformé en une suite de nombres traitables mathématiquement et statistiquement. Cette transformation, appelée la numérisation du spectre, aboutit à la constitution d’une séquence de nombres présentée sous forme d’un vecteur. Par exemple, le spectromètre FOSS NIR-Systems 5000 peut enregistrer les absorbances d’un échantillon entre 1 100 et 2 498 nm avec un intervalle de mesure constant de 2 nm, le vecteur du spectre numérisé comprendra alors 700 valeurs d’absorbance (adapté de Bertrand et al., 2006). L’ensemble de ces vecteurs ou échantillons forme une matrice sur laquelle seront appliqués les traitements de statistiques multi-variées, dont les techniques de régression multiple.

10Prétraitements de l’information spectrale. Avant d’être exploitées, les données spectrales brutes font l’objet de divers prétraitements dont les objectifs sont :

11– une atténuation ou une élimination de la non-linéarité présente entre les variables dépendantes et explicatives,

12– une élimination des interférences,

13– une élimination ou une atténuation du bruit aléatoire lié aux conditions expérimentales et au bruit électronique de l’appareil de mesure,

14– une réduction des corrélations entre les variables explicatives pour permettre l’application des techniques de régression linéaire multiple.

15La stratégie la plus courante a pour principe de soumettre les données brutes à une ou plusieurs transformation(s) mathématique(s) destinée(s) à les rendre aptes à une modélisation linéaire (Prévot, 2004). Le choix de la transformation est fonction du type de non-linéarité affectant les données et donc de la qualité des modèles prédictifs élaborés pour chaque prétraitement.

16Développement des modèles d’étalonnage. La méthode d’étalonnage la plus classique consiste à établir une relation entre la variable à prédire (y) et les variables explicatives (x₁ à x_p correspondant aux différentes longueurs d’onde) en faisant appel aux méthodes de régression linéaire multiple. La régression selon les moindres carrés partiels (PLS1) est actuellement une méthode parmi les plus employées dans le domaine de la spectroscopie infrarouge. Cette méthode sur la régression de la variable à prédire à partir de variables t₁ à t_p qui sont des variables latentes ou des combinaisons linéaires de x₁ à x_p. À la différence de la régression par composantes principales (PCA), autre méthode utilisée, les variables latentes tiennent compte à la fois des x et de y, le vecteur des facteurs est ainsi calculé de manière à maximiser le carré de la covariance entre x, le vecteur et y (Palm et al., 1995).

17Une méthode dérivée consiste à prédire les échantillons en local et en direct (Genot et al., 2011a). Dans ce cas, il n’y a pas d’équation ou de modèle à la fin de l’étape d’étalonnage, seule une base de données spectrales et les paramètres utiles à la régression PLS sont nécessaires pour prédire de nouveaux échantillons en routine, chaque spectre étant associé à une valeur du paramètre y. En prédiction, lorsqu’un nouveau spectre est saisi, il est comparé à la base spectrale et les spectres dont il est le plus proche servent à la réalisation d’une régression PLS. Cette méthode permet de prendre en compte la non-linéarité de la relation entre le paramètre étudié et la réponse spectrale. Lorsque la matrice est très hétérogène, comme c’est le cas pour les échantillons de terre, son principal inconvénient réside alors dans le nombre important de spectres et de valeurs de référence que doit contenir la base de données pour que chaque nouvel individu soit associé à un nombre suffisant de spectres proches.

18D’autres méthodes non linéaires existent également et sont de plus en plus utilisées, comme la régression localement pondérée (LWR, locally weighted regression), la régression multivariée par fonction spline adaptative (MARS, Multivariate Adaptive Regression Splines ; Shepherd et al., 2002), les arbres de régression stimulés (BRT, Boosted Regression Trees ; Brown et al., 2006) et les réseaux neuronaux (artificial neural network ; Fidêncio et al., 2002). Viscarra Rossel et al. (2010) ont comparé plusieurs méthodes d’étalonnage pour prédire le carbone organique, le taux d’argile et le pH à l’eau : la régression selon les moindres carrés partiels, les techniques MARS, forêt d’arbre décisionnel (RF, Random Forest), BRT, machines à support de vecteurs (SVM, Support Vector Machine), la décomposition en ondelettes (DWT, Discrete Wavelet Transform) et les réseaux neuronaux. Dans leur cas, la technique combinée de DWT et réseaux neuronaux proposait les meilleurs résultats en termes de prédiction. Il n’y a cependant pas de méthode universelle et selon les conditions de l’expérience (population, échantillons, paramètre à prédire, etc.), des méthodes différentes peuvent s’avérer les meilleures.

19Indicateurs de performance des modèles de prédiction. Il est indispensable de disposer de paramètres permettant de valider l’ajustement d’un modèle et/ou de comparer plusieurs modèles afin de choisir celui qui propose les meilleurs résultats en termes de prédiction. Chang et al. (2002) proposent ainsi trois groupes de paramètres permettant de tester l’ajustement du modèle pour le ou les paramètres étudiés. Il s’agit du r², du RMSEP2 ou du RMSECV3 et du RPD4 .

20Le r² correspond au coefficient de détermination entre la variable dépendante y et l’estimation de cette variable ŷ par le modèle et donc à la part de la variabilité de y expliquée par le modèle. En pratique, une valeur de r² supérieure à 0,9 indique une bonne réponse de la SPIR pour le paramètre. Entre 0,9 et 0,7, la réponse est moyenne et doit être améliorée par l’ajout d’échantillons ou par la modification de certains paramètres du modèle. En-dessous de 0,7, la réponse est considérée comme mauvaise. Si le nombre d’échantillons est faible, il faut être prudent quant à l’interprétation de ce coefficient qui peut être surévalué. Il est toujours préférable de considérer le r² et l’erreur (RMSEP ou RMSECV).

21La MSE5 est un des paramètres de mesure de l’erreur de prédiction qui correspond à la moyenne de la différence entre la valeur mesurée et la valeur prédite élevée au carré. De manière à travailler dans les mêmes unités que la variable dépendante, on l’exprime plus souvent par sa racine carrée, la RMSE. Cette dernière n’est pas connue et doit donc être estimée par une validation externe ou une validation interne. Dans le premier cas, en plus du lot d’échantillons destiné à l’étalonnage, on dispose d’un lot pour la validation. On lui applique le modèle et on définit ainsi l’erreur quadratique de prédiction (MSEP) ou sa racine carrée, la RMSEP, qui correspond à l’erreur résiduelle de prédiction. Celle-ci permet d’estimer correctement la RMSE, mais elle est cependant rarement calculée en raison du cout élevé lié à l’obtention de ce nouveau lot de données destiné à la validation. C’est pourquoi la validation interne est plus fréquemment suivie. Elle consiste à utiliser le jeu de données d’étalonnage pour réaliser la validation. Il est possible d’estimer ainsi l’erreur quadratique moyenne d’étalonnage (MSEC)6 qui correspond à l’erreur résiduelle de la régression linéaire multiple. Elle est un mauvais estimateur de la MSE qu’elle sous-évalue en ne prenant pas en compte les variances liées aux coefficients de régression. C’est pourquoi il est recommandé de travailler sur base de la validation croisée qui consiste à enlever successivement un échantillon ou un lot d’échantillons pour obtenir un nouveau lot de données indépendant auquel le modèle peut être appliqué. L’erreur quadratique moyenne de la validation croisée (MSECV) est alors obtenue. La racine carrée correspond au RMSECV qui est le paramètre le plus couramment utilisé pour valider un modèle. Il donne un résultat plus optimiste et donc inférieur à l’erreur résiduelle de prédiction mais comme il varie généralement dans le même sens que celle-ci, il est utile pour valider la méthode de modélisation choisie.

22Le RPD est le rapport entre la déviation standard du lot de données d’étalonnage pour le paramètre (sur base de la mesure obtenue selon la méthode de référence) et le RMSEP ou le RMSECV. Un RPD élevé garantit une bonne qualité de prédiction. Chang et al. (2002) estiment qu’une valeur de RPD supérieure à 2 est synonyme d’une bonne qualité de prédiction. Entre 1,4 et 2, l’ajustement est moyen et il est considéré comme mauvais pour une valeur inférieure à 1,4. Ce paramètre complète utilement le RMSEP ou le RMSECV en intégrant la variabilité de la population pour le paramètre prédit et en évitant ainsi une surestimation de la qualité de la prédiction pour des populations peu diversifiées ou représentant des territoires limités géographiquement (van Groenigen et al., 2003).

3. Prédiction des propriétés du sol

3.1. Signatures spectrales

23Comme explicité dans le paragraphe précédent, le sol est une matrice minérale et organique complexe et il est donc particulièrement délicat d’interpréter les bandes d’absorption qui y sont observées (Ben-Dor et al., 1995 ; Fidêncio et al., 2002). Outre cet aspect chimique, le spectre infrarouge est également influencé par les caractéristiques physiques du sol. La taille et la forme des particules, l’espace entre celles-ci et leur arrangement affectent la longueur de la transmission lumineuse atteignant l’échantillon et donc son absorbance. Il est donc délicat d’interpréter directement le spectre d’un échantillon de sol, d’autant plus que si les vibrations sont bien distinctes dans l’infrarouge moyen, ce n’est pas le cas dans le proche infrarouge où les harmoniques se superposent rendant très délicate l’identification de pics et donc l’utilisation de bandes spécifiques, comme illustrées à la figure 2. Pour ces raisons, les équations d’étalonnage se basent prioritairement sur une approche multi-variée en considérant la totalité du spectre, plutôt que sur une recherche de longueurs d’onde caractéristiques pour l’élément considéré. Néanmoins, des recherches fondamentales ont été conduites en vue de définir les propriétés spectrales des sols. L’objectif sous-jacent est alors la possibilité d’éliminer l’effet de ces éléments sur le spectre de manière à augmenter la sensibilité à la détection d’autres paramètres. Bogrekci et al. (2004) notent ainsi que la propriété photo-physique d’un échantillon de sol est très utile pour déterminer la distribution de la taille des particules et sa granulométrie, mais que ce même fait n’est pas désirable pour évaluer sa concentration en phosphore. Il est également avantageux de connaitre les conditions environnementales affectant les propriétés d’absorbance du sol. Si les relations entre celle-ci et les propriétés physiques et chimiques peuvent être établies quantitativement pour des conditions environnementales données, alors ces spectres pourront être utilisés pour extraire des informations utiles sur les sols (Baumgardner et al., 1985).

3.2. Humidité

24He et al. (2005), ainsi que d’autres auteurs dont Ben-Dor et al. (1995), Chang et al. (2002), Fidêncio et al. (2002), caractérisent le spectre par trois pics d’absorption de forte intensité dans le domaine du proche infrarouge (1 400, 1 900 et 2 200 nm). Viscarra Rossel et al. (1998) placent ces mêmes pics à 1 450, 1 940 et 2 200 nm. Mouazen et al. (2007) se basent sur l’intensité du pic à 1 450 nm pour évaluer le degré d’humidité du sol. Bertrand al. (2006) identifient trois bandes à 1 190, 1 450 et 1 940 nm. Les propriétés générales d’absorbance sont également influencées par le degré d’humidité. Les sols apparaissent plus foncés lorsqu’ils sont plus frais, avec pour conséquence une augmentation de l’absorbance.

25Cette propriété est un des principaux freins à l’utilisation de la SPIR pour prédire d’autres propriétés du sol par le biais de mesures in situ. Minasny et al. (2011) ont notamment étudié l’opportunité d’éliminer l’effet de l’humidité du sol pour améliorer la prédiction du carbone organique in situ.

3.3. Matières organiques

26Principes. La réponse du spectre à la quantité et à la nature de la matière organique, dans le domaine du proche infrarouge, est un fait connu et une des premières et des plus courantes applications de la SPIR en analyse de terre. Néanmoins, la relation entre ce spectre et la matière organique est très complexe, en lien avec la superposition des harmoniques dans le proche infrarouge.

27Il n’empêche que les perspectives liées à la SPIR sont très intéressantes :

28– une quantification du taux de matière organique (par l’intermédiaire de la détermination du pourcentage de carbone organique),

29– une prédiction du taux de minéralisation de l’azote et du carbone en vue d’évaluer la dynamique de la matière organique et la disponibilité des éléments pour les plantes (Fystro, 2002 ; Russel et al., 2002 ; Mutuo et al., 2006),

30– une évaluation du niveau de décomposition de la matière organique dans les sols et dans diverses matières amendantes (Coûteaux et al., 2003 ; Cozzolino et al., 2006),

31– une évaluation de la distribution du taux de matière organique au sein des différentes fractions du sol (Barthès et al., 2008),

32– une quantification physique, chimique et biologique de la matière organique (Niemeyer et al., 1992 ; Ben-Dor et al., 1995 ; Haberhauer et al., 1998 ; Hartmann et al., 2006). Certains auteurs (Ertlen et al., 2010) ont également étudié l’opportunité d’utiliser la SPIR pour identifier l’origine végétale de la matière organique.

33Performances de prédiction. L’évaluation du taux de carbone organique (CO) par la SPIR a fait l’objet de nombreuses études comme détaillé dans le tableau 1 et dans la revue bibliographique effectuée par Stenberg et al. (2010). Lorsque le taux de CO varie dans la gamme des valeurs habituellement rencontrées pour la couche arable des sols cultivés, les modèles prédisent avec une précision de l’ordre de 0,2 à 0,5 g·100 g^-1, à l’exception de celui d’Udelhoven et al. (2003) qui prédit avec une erreur d’environ 1,5 g·100 g^-1. Les auteurs ayant travaillé sur un plus large éventail de taux de CO (entre 0 et plus de 20 g·100 g^-1 de CO) obtiennent de bons résultats en termes de RMSECV ou RMSEP, variant entre 0,6 et 0,8 g·100 g^-1 et des valeurs de RPD supérieures ou proches de 2. Ces résultats montrent l’acuité du dosage du taux de CO par SPIR. À l’exception de l’étude d’Udelhoven et al. (2003), le taux de CO peut être correctement prédit par la SPIR. Les méthodes classiques, détaillées dans le tableau 1, fournissent des résultats avec une variation de l’ordre de 0,2 g·100 g^-1 de CO. Pour le domaine des sols sous culture, avec un taux de CO variant entre 0,5 et 5g·100 g^-1, les prédictions de la SPIR approchent les résultats obtenus par les méthodes de référence. Ces écarts de performances, malgré le nombre important de bandes liées à la matière organique, peuvent être dus au fait que le signal dans la région proche infrarouge est faible, particulièrement dans des sols contenant peu de matière organique et caractérisée par une matrice minérale variable (Viscarra Rossel et al., 1998). Même si nous n’avons comparé que des études portant sur un grand nombre d’échantillons, la diversité des dispositifs expérimentaux, la taille et l’homogénéité de la population d’échantillons peuvent également expliquer ces différences.

34Certaines recherches, Morra et al. (1991) et Cozzolino et al. (2006), se basent sur une détermination du CO en séparant au préalable les différentes fractions granulométriques. Cozzolino et al. (2006) ont, de plus, réalisé pour chaque fraction un modèle sur base des spectres bruts et de leur dérivée seconde. Sur base de 180 échantillons, ils observent que les taux de CO varient en fonction de la fraction considérée : pour celle des sables, les taux s’échelonnent de 0,02 à 1,1 % et de 0,5 à 4,8 % pour les fractions limoneuses et argileuses. Les résultats du modèle obtenu en effectuant la dérivée seconde des spectres sont légèrement moins bons, en termes de prédiction. Les valeurs du r² sont de l’ordre de 0,90, celles du RMSECV varient aux alentours de 0,09 %, celles du RMSEP varient entre 0,09 % (fraction sableuse) et 0,25 % pour les fractions limoneuse et argileuse. En comparaison aux approches classiques menées sur l’échantillon de terre, les résultats s’ajustent mieux lorsqu’une séparation granulométrique est effectuée, chaque fraction étant irradiée par la lumière infrarouge. Il parait néanmoins délicat de suivre cette procédure dans le cadre d’un travail de routine, le temps passé à séparer ces fractions étant trop important.

3.4. Taille des particules

35Comme pour les matières organiques, la réponse du spectre à la taille des particules est en même temps très intéressante pour déterminer la granulométrie de l’échantillon et, au contraire, peut masquer l’effet de certains autres constituants du sol sur le spectre. Les études menées dans ce domaine ont permis d’identifier certains pics liés aux groupements OH et Al-OH en relation avec le taux et la minéralogie des argiles. Les bandes d’absorption relatives aux groupements OH sont également retenues, l’essentiel des autres pics identifiés s’observant au-delà de 2 000 nm.

36Selon Orlov7 (1966), Kojima8 (1958) et Zwerman9 (1940) cités par Baumgardner et al. (1985), l’effet de la taille des particules sur l’intensité du spectre est une fonction inverse du diamètre de la particule. Les terres les plus sableuses, présentant un diamètre plus important, sont logiquement celles pour lesquelles les bandes d’absorption de l’eau sont les moins marquées.

37Bogrekci et al. (2004) ont étudié l’effet de la taille des particules d’échantillons sableux sur la détermination spectrale de la concentration en phosphore disponible. Ils concluent :

38– qu’il existe une différence spectrale significative entre les fractions sableuses fines, moyennes et larges, les particules les plus grandes absorbant davantage la lumière que les particules moyennes et fines,

39– que de meilleurs résultats sont obtenus, en termes d’étalonnage, lorsque le phosphore est évalué pour chacune des fractions plutôt que sur l’échantillon complet duquel l’effet granulométrique a été ou non soustrait.

40Dans le même sens, Brunet et al. (2007) ont également étudié l’importance de la préparation de l’échantillon pour prédire les taux de carbone et d’azote. Ils ont comparé la qualité de la prédiction d’échantillons émottés et tamisés à 2 mm ou broyés à 200 µm sur base d’un set d’échantillons représentatif de situations pédologiques variées. Ils ont souligné l’importance de travailler sur des échantillons broyés et de même nature.

3.5. Capacité d’échange cationique

41Baumgardner et al. (1985) indiquent que la CEC présente fréquemment des corrélations plus élevées avec l’absorbance que la matière organique ou le pourcentage d’argile considéré individuellement. Ceci est à mettre en relation avec les constituants du sol en termes de quantité et de nature, ainsi qu’avec les phénomènes de complexations organo-minérales, non considérés lors des dosages individuels de l’argile et de la matière organique.

42Comme le précisent Chang et al. (2001), la CEC est une propriété secondaire du spectre. La CEC est effectivement fonction des propriétés organiques, minérales et organo-minérales du sol. Ces propriétés étant des propriétés primaires du spectre, la CEC peut être prédite par la SPIR. Les résultats de cette prédiction sont plus précis que ceux obtenus via une équation liant les taux d’argile et de matière organique.

43Chang et al. (2001), Shepherd et al. (2002), Islam et al. (2003) et van Groenigen et al. (2003) ont notamment testé l’opportunité d’appliquer la SPIR pour évaluer la CEC de l’échantillon de terre analysé (Tableau 2).

44van Groenigen et al. (2003) obtiennent de très bonnes prédictions, avec une précision de l’ordre de 2 cmol⁺·kg^-1 de TS et un RPD supérieur à 2, ceci pour une gamme de valeurs de CEC assez large, puisque variant environ entre 2 et 40 cmol⁺·kg^-1 de TS. Ce fait peut être mis en relation avec la méthode utilisée qui fournit, ici, une estimation de la CEC effective du sol, obtenue après sommation des cations échangeables et de l’acidité d’échange. Il est possible que le spectre obtenu soit plus directement lié à une valeur effective de la CEC plutôt qu’à une valeur potentielle obtenue en tamponnant le sol à une valeur de pH donné (pH 7,0 ou 8,2 selon les méthodes). Les autres études montrent également une bonne réponse du sol pour l’évaluation de la CEC. La précision de la valeur prédite est d’environ 4 cmol⁺·kg^-1 de TS.

3.6. Composition minéralogique

45Le spectre du sol est influencé par sa composition minérale (Baumgardner et al., 1995). Ce fait rend plus difficile l’évaluation de la granulométrie du sol ou nécessite de travailler par grande famille minéralogique pour mieux évaluer la texture. Cette propriété présente, par contre, l’avantage de pouvoir envisager une première évaluation de la composition minéralogique du sol par l’étude de l’empreinte spectrale. La présence d’ions ferriques, en faible quantité dans les kaolinites par opposition aux montmorillonites, se marque également au niveau du spectre par un pic d’absorption d’intensité variable, aux alentours de 1 100 nm. Les illites se distinguent également par une absorbance plus importante pour les longueurs d’onde inférieures à 1 700 nm. De même, les bandes d’absorption caractéristiques de l’eau et des hydroxydes sont d’une intensité plus faible pour les illites, par comparaison aux kaolinites et aux montmorillonites (Baumgardner et al., 1995).

46En ciblant les ions ferriques et les hydroxydes, il est également possible de distinguer des minéraux comme l’hématite, la limonite, la goethite. La gibbsite, en relation avec les hydroxydes d’aluminium, est également reconnaissable sur le spectre infrarouge (Baumgardner et al., 1985).

47Demattê et al. (2004) ont conclu que le quartz, la magnétite, la kaolinite et la smectite étaient les minéraux qui influençaient principalement le spectre et permettaient de discriminer les sols.

48La couleur plus blanche des sols riches en carbonate de calcium influence les propriétés spectrales du sol. Chang et al. (2002) notent qu’un sol plus blanc entraine une plus forte réflectance et donc une plus faible absorbance du sol. Ils indiquent, en outre, deux pics de forte intensité à 2 340 et 2 500 nm et trois plus petits respectivement à 1 870, 2 000 et 2 160 nm, dont l’intensité permet de doser la teneur en carbonate de calcium du sol. Ils attribuent ces pics aux harmoniques et aux combinaisons des bandes fondamentales observées dans l’infrarouge moyen en relation avec la vibration des liaisons de la molécule CO₃. Gomez et al. (2012) ont également montré l’opportunité de prédire le taux de carbonate de calcium et le taux de fer par SPIR.

4. Discussion relative aux principaux résultats en analyse de terre

49Les différentes études de faisabilité liées à l’application de la SPIR pour évaluer les paramètres physiques, chimiques et physico-chimiques s’inscrivent dans trois cadres principaux :

501. Des approches à l’échelle régionale. L’objectif est, alors, de tester les paramètres répondant positivement à l’analyse et de développer des équations d’étalonnage pour ceux-ci à travers le développement de librairies spectrales à l’échelle d’une région, d’un ou de plusieurs pays. De nombreuses études sont réalisées dans ce sens. Leur principal point faible, à l’exception de certaines études comme celles menées par Chang et al. (2001), Shepherd et al. (2002) et Brown et al. (2006), est le nombre restreint d’échantillons utilisés. Seules les perspectives de la SPIR sont alors envisagées et une application en routine n’est pas directement possible. De plus, le risque principal est une surestimation de la valeur des résultats de la prédiction (van Groenigen et al., 2003).

512. Des approches à l’échelle de la parcelle en vue de répondre aux besoins de l’agriculture de précision (He et al., 2005). Cette voie demande de connaitre précisément la variabilité spatiale des paramètres physiques et chimiques d’une parcelle. Un nombre important d’échantillons est soumis à l’analyse et la technique de la SPIR permet donc une économie considérable. Cette approche voit également se développer des appareillages permettant des mesures in situ (Mouazen et al., 2007 ; Maleki et al., 2008).

523. Des approches plus fines, permettant de caractériser plus précisément certains paramètres du sol et leur évolution. L’objectif n’est plus alors de remplacer des techniques existantes, mais de les compléter. Malgré leur intérêt dans le conseil de fumure, certaines propriétés de sol ne sont pas déterminées parce qu’elles font appel à des techniques trop longues ou trop couteuses. Une approche plus fine de la composition des matières organiques (Coûteaux et al., 2003), mais aussi de la composition minérale du sol, pourraient apporter une réelle plus-value au conseil de fumure (Vendrame et al., 2012).

53Les paramètres physiques, chimiques et physico-chimiques susceptibles d’être déterminés par la SPIR sont l’humidité résiduelle, le pH, la matière organique (carbone organique, azote total et autres composantes), les fractions granulométriques (taux d’argile, de limon et de sable et caractérisation minéralogique), la capacité d’échange cationique. En outre, des recherches ont été menées pour prédire certains éléments disponibles, certains éléments totaux majeurs ou traces. C’est ainsi que Mouazen et al. (2007) ont évalué l’opportunité, notamment, de doser le phosphore disponible sur des échantillons séchés, émottés et tamisés.

54Les principales sources de variation entre les études sont :

55– le nombre d’échantillons et la diversité pédologique ; certaines études ne concernent que des échantillons de sols très homogènes, alors que d’autres tentent de couvrir une large diversité de situations pédologiques (Chang et al., 2001 ; Shepherd et al., 2002 ; Brown et al., 2006 ; Barthès et al., 2008 ; Viscarra Rossel, 2009 ; Genot et al., 2011a) ;

56– l’application d’analyses de références variables en fonction de la région considérée, les méthodes ne sont pas toutes normalisées et comparables, cette particularité rendant parfois difficile la comparaison des résultats des modèles prédictifs (Viscarra Rossel, 2009) ;

57– la gamme de longueurs d’onde ainsi que le type de spectrophotomètre ;

58– le choix des méthodes statistiques qui diffère selon les objectifs de l’étude, selon l’état des connaissances en termes de puissance des ordinateurs et des techniques statistiques et mathématiques existants au moment de la réalisation des études ;

59– les prétraitements réalisés sur les spectres avant la phase d’étalonnage, ainsi que les méthodes de détection des valeurs aberrantes.

60La SPIR fait l’objet de nombreux travaux tant son potentiel d’utilisation est large. Cependant, pour envisager une utilisation de la technique à une échelle plus large que celle de la recherche, il apparait important de bien tenir compte d’un certain nombre d’éléments, ceci afin d’éviter que les conclusions ne valent que pour la parcelle, le cas d’étude, etc. et qu’elles ne surestiment les capacités réelles de la technique à prédire les propriétés des échantillons. L’étude préliminaire du territoire pour lequel on souhaite établir des modèles de prédiction est une étape essentielle pour s’assurer que les modèles puissent être utilisés pour l’ensemble de celui-ci. Un modèle d’étalonnage développé sur base d’échantillons prélevés sur une même parcelle ne pourra être utilisé que pour prédire d’autres échantillons provenant de cette parcelle ou des échantillons provenant d’une parcelle présentant les mêmes propriétés. Développer des modèles à l’échelle d’une région demande d’appréhender la variabilité de cette région de manière à ce que le lot d’échantillons d’étalonnage couvre cette variabilité. En fonction de la diversité rencontrée, le nombre d’échantillons à prélever sera plus ou moins important, de même que le nombre d’analyses de référence à réaliser. Il est également essentiel d’apporter un soin particulier aux analyses de référence qui seront réalisées sur les échantillons de cette base de données d’étalonnage, la qualité des prédictions y sera étroitement liée puisque l’incertitude liée à la méthode de référence sera cumulée à celle de la technique infrarouge. En outre, il est indispensable que les échantillons constitutifs du lot d’étalonnage soient scannés selon la même procédure et que tout nouveau spectre soit également saisi selon cette même procédure. Enfin, comme l’ont montré Hunt et al. (197310, 197411) cités par Baumgardner et al. (1985) ou Barthès et al. (2008), il peut être préférable de grouper les échantillons de sol selon le type de matériaux parentaux, de manière à obtenir une plus grande acuité des résultats pour les autres paramètres à déterminer. D’autres études ont été réalisées dans ce sens : Léone et al. (2000) utilisent les propriétés du spectre infrarouge pour grouper les échantillons de sol et comparer ces groupements aux principaux types de sol rencontrés. Demattê et al. (2004) qui, en mettant en évidence les principaux paramètres influençant le spectre et permettant une discrimination des sols, concluent que ceux-ci sont la matière organique, le taux de fer total, de limon et de sable et la composition minéralogique du sol. Viscarra Rossel et al. (2011) ont également utilisé la spectroscopie proche infrarouge pour classifier les grands types de sols australiens et les représenter cartographiquement. Si l’utilisation en routine de la spectroscopie proche infrarouge est peu onéreuse une fois le spectromètre acquis, rapide, simple d’utilisation et non consommatrice d’extractifs chimiques, la phase d’étalonnage préalable à son utilisation peut être longue et onéreuse selon les propriétés étudiées. C’est pourquoi il est tentant d’établir des librairies spectrales associant les valeurs de références pour différentes propriétés de sol et les spectres proche infrarouge. Ces librairies spectrales pourraient alors être utilisées par toutes les structures de recherche et les laboratoires de routine travaillant sur un même territoire. Dans le même sens, il est intéressant de pouvoir profiter de pédothèques que détiennent souvent les universités et centre de recherche. Cette mise en commun de librairies spectrales a été tentée au niveau mondial par la mise en place du « Soil Spectroscopic Group » (Viscarra Rossel, 2009). Plusieurs universités ont collaboré pour scanner leurs échantillons sur base d’un protocole bien arrêté et mettre en commun les spectres, les valeurs de référence et les caractéristiques pédologiques des échantillons. La principale faiblesse de ces librairies spectrales à l’échelle mondiale est l’absence de maitrise des analyses de référence. Les laboratoires ayant contribué à l’étude ne travaillent pas tous selon les mêmes protocoles. Une part importante de l’erreur de prédiction obtenue peut être liée à cette absence d’harmonisation des protocoles analytiques. Au niveau d’un territoire plus restreint regroupant des laboratoires travaillant selon des protocoles normalisés, cette approche est certainement celle à réaliser pour envisager à court ou à moyen terme une utilisation en routine de la SPIR en analyse de terre.

5. Perspectives

61Dans le cadre spécifique d’une contribution de la technique à l’évaluation de la fertilité du sol, l’interprétation du spectre est utile d’un point de vue quantitatif :

62– en remplacement de méthodes existantes couteuses en temps, matériel ou extractifs chimiques nuisibles pour l’environnement,

63– en apportant des informations spatialement plus précises sur des paramètres pour lesquels des moyennes régionales sont utilisées par défaut,

64– en complément du menu classique d’analyse, de manière à apporter d’autres informations utiles à l’évaluation de la fertilité du sol. Une approche qualitative peut également être intéressante pour reconnaitre l’appartenance du sol à un groupe de sols pour lesquels des référentiels régionaux existent ou reconnaitre la nature de la matière organique.

65Les perspectives d’utilisation de la spectroscopie proche infrarouge sont nombreuses et variées, depuis une utilisation classique en remplacement ou en complément des analyses classiquement effectuées en laboratoire, jusqu’à une utilisation in situ ou par le biais des techniques de télédétection (Stevens et al., 2006). Son champ d’application est très large :

66– Elle contribue à mieux appréhender la variabilité intraparcellaire des parcelles agricoles en vue d’une bonne gestion économique et environnementale des intrants. Par sa facilité d’utilisation et sa rapidité, cette technique utilisée au laboratoire en remplacement des techniques classiques permet d’envisager l’analyse d’un plus grand nombre d’échantillons et ainsi de mieux caractériser la parcelle et le conseil de fumure lié (Genot et al., 2011b). En l’utilisant sur du matériel embarqué (Mouazen et al., 2007), cette technique s’intègre dans les outils nécessaires à l’agriculture de précision en permettant d’adapter les intrants en fonction des caractéristiques physico-chimiques de la parcelle.

67– La technique est également utilisée pour cartographier sur un territoire plus ou moins large des propriétés physico-chimiques des sols, en complément des méthodes géostatistiques. À l’échelle parcellaire, cette technique permet de cartographier la variabilité d’une parcelle pour les propriétés de sol étudiées et de l’utiliser dans le cadre de l’agriculture de précision (Cressie et al., 2008). À l’échelle d’un territoire, la technique SPIR couplée à la géostatistique permet de cartographier les différentes propriétés physico-chimiques du sol (Viscarra Rossel et al., 2011). Selon ces auteurs, il est alors possible de construire sur cette base une carte des sols utiles pour de nombreuses applications agro-environnementales sur base de leurs caractéristiques morphologiques mais aussi physiques, chimiques et physico-chimiques telles que mesurées en laboratoire.

68– La technique présente également des atouts majeurs en vue du monitoring de la qualité des sols agricoles ou forestiers (Cécillon et al., 2009). L’intérêt de suivre au cours du temps la qualité environnementale ou agronomique des sols n’est plus à démontrer ; néanmoins, le cout des analyses est souvent un frein au maintien de ces projets. Dans ce cadre, la technique peut nettement permettre de diminuer ce cout, une fois l’étape d’étalonnage effectuée. Un second problème inhérent à ces projets est l’évolution des techniques et des objectifs au cours du temps. Les échantillons de terre ne sont pas toujours gardés et sont susceptibles d’évoluer dans le temps, au contraire des spectres qui peuvent être conservés à l’infini. Grâce à ceux-ci, il sera donc toujours possible de prédire de nouvelles propriétés de sol, même sans les échantillons. Vu la facilité de mise en œuvre de la technique, tout échantillon prélevé devrait donc être scanné, même si l’intérêt immédiat ne le justifie pas.

69La SPIR est une technique d’avenir pour l’analyse de terre tant au laboratoire qu’in situ. Les nombreuses recherches montrent toutes ses atouts et son potentiel dans ce domaine. Il lui reste à s’inscrire dans la routine des laboratoires d’analyse en remplacement ou en complément des techniques analytiques actuelles.

Bibliographie

Notes

Para citar este artículo