Amandine Liénard, Laurent Bock & Gilles Colinet

Intérêt des cartes des sols pour l'élaboration d'une stratégie d'échantillonnage en sols contaminés par retombées atmosphériques : application à l'étude de l'effet sol sur le devenir des éléments traces métalliques

Document(s) associé(s)

Annexes

• pdf (2,6M)

Résumé

Les sites contaminés par des retombées d'éléments traces métalliques (ETM) tels que le cadmium, le plomb et le zinc ne sont pas rares en Wallonie (Belgique). Ces milieux enrichis en ETM sont dénommés « sites calaminaires » et sont le résultat du passé industriel de la région. L'objectif principal de cette étude est la caractérisation de la distribution multi-échelle des ETM dans ces sols et paysages contaminés. La stratégie d'échantillonnage développée dans cet article est basée sur l'étude de l'influence du type de sol et de son occupation sur la distribution des ETM, ceci en vue d'évaluer les risques encourus par des changements environnementaux, par exemple en cas d'exploitation agricole de zones boisées.

Abstract

Interest of soil maps to elaborate a sampling strategy in soils contaminated by atmospheric fallouts: application to the study of soil effect on the fate of metallic trace elements. As a result of pollution from atmospheric dusts, some sites of Wallonia (Belgium) are contaminated by metallic trace elements (MTE) such as cadmium, lead and zinc. These metalliferous soils are called ‘calaminary sites'. The major aim of this study concerns the multi-scalar characterization of MTE distribution in these contaminated soils and landscapes. The strategy of sampling described in this paper is based on the influence of soil type and its land use on the fate of MTE. After analysis of collected samples, it will be possible to have an idea of the distribution of MTE among soil components, and therefore, to evaluate the risks caused by environmental changes such as an agricultural use of afforested areas.

Table des matières

1. Introduction

1Au cours des deux siècles passés, des usines présentes en Wallonie (Belgique) se sont spécialisées dans la transformation de minerais riches en métaux. Suite à leurs activités, elles ont émis des poussières chargées en éléments traces métalliques (ETM), principalement du cadmium, du plomb, du zinc et des fumées riches en oxydes acides (SO₂, SO₃, NO₂, etc.) (Graitson, 2005). Ces retombées ont favorisé la formation de sites calaminaires¹, dont le sol est enrichi en ETM induisant une forte toxicité du milieu (Graitson, 2005 ; Graitson et al., 2005 ; Rosengarten, 2010). La majorité des sites calaminaires wallons se situent en province de Liège dans les bassins de la Gueule, de la Vesdre et de la Basse Ourthe ainsi que dans la vallée de la Meuse. Leur surface est estimée à 360 ha (Graitson, 2005) (Figure 1). Les sites sont classés en trois grands types selon l'origine de la contamination :

2– les sites primaires, qui sont les affleurements et les filons de minerais,

3– les sites secondaires, qui sont soit les déblais résultant de l'exploitation des mines et des excavations, soit les cendrées et scories sous-produits des usines produisant le métal,

4– les sites tertiaires qui résultent des retombées atmosphériques de particules métalliques (Duvigneaud et al., 1993 ; Saintenoy-Simon et al., 1996).

5Depuis longtemps, les naturalistes portent un grand intérêt à ces milieux calaminaires. La diversité de leur faciès ainsi que leurs origines diverses ont permis le développement d'une flore spécifique, la flore calaminaire (Faucon et al., 2011). Cette flore est capable de coloniser ces milieux ouverts, pauvres en matières organiques et en nutriments mais riches en éléments toxiques comme les ETM (Bizoux et al., 2004 ; Faucon et al., 2011). Les taxons rares et endémiques de ces milieux sont classés en deux groupes (Lambinon et al., 1963 ; Graitson, 2005) : les métallophytes (taxons présents uniquement sur les sols calaminaires dans une région déterminée) et les pseudométallophytes (taxons également présents sur les sols non calaminaires). Grâce à cette diversité de milieux, les sites calaminaires attirent également une entomofaune formée pour la plupart d'espèces thermophiles. Ces insectes apprécient la végétation rase et éparse, ainsi que le microclimat chaud et sec des haldes² (Ertz et al., 2001 ; Graitson et al., 2005 ; Quoilin et al., 2007 ; Rosengarten, 2010).

6Sans la présence des usines de traitements de minerais, la très grande majorité des sites calaminaires n'existerait pas. En effet, seuls quatre sites sur les trente-sept que compte la Wallonie sont de type primaire (Graitson, 2005). De fait, c'est en partie suite à la production de poussières dans l'environnement proche des usines et par leur altération minérale dans les cours d'eau voisins que l'étendue des sols contaminés a dépassé largement les limites des sites calaminaires (Tordoff et al., 2000). Cette pollution diffuse peut avoir un impact sur la qualité des terres arables environnantes et engendrer des risques en matière de santé publique. Dans la région de Prayon (Province de Liège, Vallée de la Vesdre), dès 1964, Ramaut avait constaté que la contamination en zinc et les dégâts engendrés par les fumées chargées en SO₂ provenant des usines de la Société Métallurgique de Prayon (Nouvelle Montagne) s'étaient étendus au-delà de la vallée. L'ampleur de la dégradation des écosystèmes dépend de la nature et des quantités des retombées polluantes, de la direction des vents dominants, des précipitations atmosphériques, du relief, du type de sol et de la végétation originelle (Denaeyer-De Smet et al., 1974).

7Outre leur action toxique sur la végétation, les poussières provoquent une contamination importante de l'horizon de surface des sols (Ramaut, 1964). Cette contamination de surface peut engendrer un risque pour la durabilité des écosystèmes, la qualité des productions végétales et la santé humaine et ceci même quand la contamination est moins lourde (Cambier et al., 2009). Dans ce cas, le risque n'est pas seulement dû à la quantité totale d'éléments traces présents dans le sol, mais bien à la capacité du sol à retenir un ETM et à empêcher son transport, soit à l'intérieur de la matrice solide, soit à l'extérieur de celle-ci (Blanchard, 2000). Les phases réactives des sols que sont les argiles, les carbonates, les oxydes métalliques, les silicates et la matière organique influencent la mobilité des ETM dans les sols. Ces divers constituants jouent en effet un rôle direct dans les relations d'échange (sorption, complexation et précipitation) entre phase aqueuse et solide (Blanchard, 2000 ; Citeau, 2004 ; Lamy et al., 2006).

8Le devenir des ETM apportés aux sols dépend des constituants et des propriétés des sols, mais également de la végétation qui s'y développe ou de l'usage qui est fait du sol (Fernandez-Cornudet, 2006). L'objectif de nos recherches s'inscrit dans l'évaluation des risques liés à la mobilité environnementale des constituants métalliques. Les objectifs spécifiques sont d'évaluer les contributions respectives des facteurs sol, occupation et direction des vents sur l'ampleur de la contamination et sur l'évolution des contaminants dans le sol. Ce papier présente et discute la méthodologie d'échantillonnage qui a été conçue et appliquée à trois paysages contaminés, chacun d'eux entourant un site calaminaire.

2. Matériel et méthodes

2.1. Contexte historique des sites calaminaires

9Les zones d'étude sont toutes trois centrées sur un site calaminaire de type secondaire et/ou tertiaire. Le choix des sites calaminaires s'est porté sur :

10– le site de Sclaigneaux repris dans le site Natura 2000 « Vallée de la Meuse de Marche-les-Dames à Andenne » (Code BE35006), dont l'entière superficie est une réserve naturelle (Maingeot, 2003 ; Michaux, 2010),

11– le site calaminaire de Prayon (La Rochette) qui fait partie du site Natura 2000 « Basse vallée de la Vesdre » (Code BE33016),

12– l'ancien site minier de Plombières (Bleiberg) intégré dans le site Natura 2000 « Vallée de la Gueule en aval de Kelmis » (Code BE33006) dont une partie est érigée en réserve naturelle (Graitson, 2005).

13Le contexte historique de formation des trois sites calaminaires (Sclaigneaux, Prayon et Plombières) est décrit dans les paragraphes suivants.

14Dès le Moyen Âge, des veines de zinc et de plomb ont été découvertes à proximité du site de Sclaigneaux, au lieu-dit « Haies-Monet » (Dejonghe, 1998 ; Maingeot, 2003). Entre 1848 et 1873, une concession a été ouverte pour des mines de plomb, de zinc et de pyrite sur les anciennes communes de Landenne et Seilles (Dejonghe, 2010). Cette concession a fourni 40 000 tonnes de minerai concentré, soit environ 75 000 tonnes de minerai brut extrait (Dejonghe, 2010). Parallèlement, dès 1856, la société Dumont installe ses fours à zinc au pied des falaises de Sclaigneaux (Michaux, 2010). Elle achète alors un ravin entre le ruisseau de Loysse et les rochers de Sclaigneaux, ravin qu'elle va combler à l'aide des déchets du traitement des minerais (De Nul, 2010). Vers 1880, de hautes cheminées sont construites sur la crête, elles sont reliées à l'usine et rejettent ainsi, jusqu'en 1940, des fumées fortement chargées en métaux lourds (plomb, zinc, cuivre, cadmium, arsenic, etc.) et oxydes gazeux de soufre et d'azote (Michaux, 2010). L'usine cesse toute activité en 1971 et les cheminées sont abattues au début des années 1980 (Michaux, 2010).

15Les minerais de plomb du site de Prayon, plus précisément du lieu-dit « La Rochette », n'ont été découverts qu'au cours du 16^e siècle et seulement exploités à partir de 1730. De 1902 à 1925, la S.A. Vieille Montagne a extrait du plomb, du zinc et de la pyrite de la concession (Dejonghe et al., 1993). Auparavant vers 1829, l'industrie du traitement des minerais avait installé les premiers fours à zinc à proximité de Prayon. Ces derniers ont été transformés en une fonderie à zinc et plomb en 1851 (Jadot, 2007). Durant toute cette période et jusqu'en 1976, huit cheminées ont émis des fumées chargées en poussières métalliques à partir du complexe de traitement des minerais posté dans la vallée (Duvigneaud et al., 1987 ; Meerts et al., 2001 ; Jadot, 2007).

16L'ancien site minier de Plombières apparait dans la littérature pour la première fois en 1365. Cependant, son exploitation n'a démarré réellement qu'en 1828 (Dejonghe et al., 1993). En 1882, les gisements sont abandonnés par manque de rentabilité et suite aux problèmes d'exhaure dus à la Gueule toute proche. Cette fermeture a eu lieu après que quelque 115 000 tonnes de galène (sulfure de plomb) et 110 000 tonnes de blende (sulfure de zinc) aient été extraites (Dejonghe et al., 1993 ; Saintenoy-Simon et al., 1996). Malgré la fermeture de la mine, la fonderie a fonctionné avec cinq fours jusqu'en 1922 grâce aux tonnes de minerai importées par la société. Toutes les activités ont cessé entièrement en 1969 (Saintenoy-Simon et al., 1996 ; Dispas, 2009).

17Le développement d'usines de traitements de minerais est à l'origine des trois sites. Dans un premier temps, les usines ont exploité les gisements existants : ceux de Plombières et Sclaigneaux qui sont dans les calcaires ou la dolomie du Dinantien et ceux de Prayon dans des calcaires plus anciens du Frasnien et du Givetien (Dejonghe, 1998). Chacune des concessions a été fermée par manque de rentabilité, l'extraction devenant trop contraignante. Les usines se sont alors tournées vers des minerais étrangers, plus riches en plomb et en zinc (Dejonghe et al., 1993).

18Au niveau floristique, les trois sites se distinguent les uns des autres. L'ancien site minier de Plombières (20 ha) présente la plus grande diversité d'espèces calaminaires avec cinq espèces métallophytes, le site de Prayon (75 ha) compte cinq espèces métallophytes dont deux introduites récemment, tandis qu'aucune flore métallophyte n'a été observée au sein de la réserve de Sclaigneaux (39 ha) (Graitson, 2005).

2.2. Délimitation des zones d'étude

19Sur les sites calaminaires, la présence de teneurs élevées en ETM se remarque par une végétation rase et spécifique (Graitson, 2005 ; Rosengarten, 2010). Cependant, un amendement organique apporté aux sols peut détruire cette flore en donnant l'avantage absolu aux plantes classiques de la prairie pâturée (Hermanns et al., 1996). Dans ce cas, il est difficile de différencier sur le terrain un sol sain d'un sol contaminé sans en mesurer les teneurs. La tâche est d'autant plus complexe que les sols contaminés peuvent être urbains, agricoles ou forestiers et peuvent s'étendre sur une distance de plusieurs kilomètres (Cambier et al., 2009).

20Actuellement en Europe, de nombreuses études ont été réalisées afin de diagnostiquer et modéliser l'ampleur et les effets toxicologiques des contaminations atmosphériques en ETM (Hough et al., 2004 ; Verma et al., 2007 ; Descamps et al., 2009 ; Gillet et al., 2009 ; Mench et al., 2009 ; van Oort et al., 2009). Un échantillonnage de qualité doit prendre en compte la complexité du milieu contaminé et présenter une échelle de perception adaptée à la question posée (van Oort et al., 2009). Cependant, il n'existe pas d'instructions ou de directives européennes en matière d'échantillonnage de sols contaminés (Theocharopoulos et al., 2001). Chaque pays a sa propre démarche et ses propres critères de sélection. Il n'est donc pas aisé de choisir un plan et un périmètre d'étude en se basant uniquement sur la littérature. De plus, les principaux critères qui orientent les campagnes d'échantillonnage sont généralement d'ordre économique (Back, 2007). Néanmoins, la caractérisation du milieu biophysique (relations entre géologie, relief, sols, utilisation des sols et historique des activités humaines) sur base des documents existants semble indispensable pour orienter l'échantillonnage dès qu'on souhaite évaluer les redistributions d'éléments dans les paysages en dehors des sites identifiés comme contaminés.

21Dans le cadre de cette étude, le périmètre défini est une maille de 3 km de rayon centrée sur l'emplacement historique des cheminées des usines de traitements de minerai. Cette échelle de perception est suffisante pour à la fois prendre en considération le facteur de distance et rencontrer une importante diversité d'unités paysagères. De plus, la superficie d'environ 2 800 ha de chaque zone d'étude est entièrement circonscrite dans un seul bassin versant : la Meuse aval, la Vesdre et la Gueule respectivement pour les zones de Sclaigneaux, Prayon et Plombières.

2.3. Méthodologie d'échantillonnage

22La répartition des poussières contaminantes au sein d'un paysage est influencée par :

23– la distance à la source de contamination,

24– la direction des vents qui crée des zones plus exposées à la contamination,

25– la nature des sols et leur position dans le relief qui conditionnent les redistributions d'éléments suite aux phénomènes d'érosion et aux migrations verticales, de même que leur vocation agricole ou forestière (Meshalkina et al., 1996 ; Fernandez-Cornudet, 2006 ; van Oort et al., 2009).

26Par le plan d'échantillonnage développé ci-après, nous entendons déterminer l'importance relative des facteurs influençant la distribution spatiale des teneurs en ETM dans le paysage et leur devenir dans le sol.

27L'échantillonnage a été stratifié en fonction des trois critères : position par rapport à la source, nature du sol et occupation. La même procédure a été appliquée aux trois zones d'étude. Les données cartographiques utilisées pour orienter l'échantillonnage sont :

28– la Carte des Principaux Types de Sols de Wallonie à 1/250 000 (Bock et al., 2006), basée sur la Carte des Sols de la Belgique éditée à 1/20 000 par l’IRSIA de 1947 à 1991 (Legrain et al., 2011),

29– la Carte d'Occupation des Sols de Wallonie (COSW) de 2005 à 1/10 000. Les données cartographiques sont exploitées à l'aide du logiciel ArcGis 9.3.

30La procédure de sélection des plages de sols échantillonnées comporte six étapes (Figure 2) :

311. La CNSW250 et la COSW sont projetées et découpées à la dimension de la zone étudiée, soit un buffer de 3 km autour de l'emplacement des cheminées.

322. Parmi les principaux types de sols (PTS) de la CNSW250 répartis sur la zone, certains sont abandonnés faute de caractérisation suffisante. Ce sont les sols non cartographiés ou artificiels et les regroupements de complexes de sols de textures différentes ou sur fortes pentes et de sols de fonds de vallons limoneux ou rocailleux.

333. La COSW est utilisée pour sélectionner les occupations étudiées :

34– les terres arables, les cultures permanentes et les friches agricoles,

35– les surfaces enherbées,

36– les forêts et milieux semi-naturels respectivement dénommés dans cette étude culture, prairie et forêt.

374. Les deux nouvelles couches cartographiques ainsi créées sont croisées pour donner des entités « Sol-Occupation ». Suite à cette étape, seule la superficie de sol conservée est celle couverte par les trois occupations choisies.

385. Le troisième facteur pris en considération pour l'échantillonnage est la direction des vents. Pour le représenter, les zones d'étude sont divisées en quatre quadrants chacun orienté dans une direction, soit le Nord, le Sud, l'Est et l'Ouest. Seules les entités « Sol-Occupation » présentes dans tous les quadrants sont retenues pour les prélèvements. Par zone d'étude, il ne reste que trois principaux types de sols qui répondent aux critères précités.

396. La dernière étape consiste à identifier des points d'échantillonnage au sein des plages cartographiques retenues de manière à couvrir un maximum de gammes de distances au centre de la zone d'étude. Seules les plages de dimension pertinente, soit minimum 300 m², sont reprises dans le dispositif. Les combinaisons « Sol-Occupation » sont échantillonnées en fonction de leur importance surfacique. Cette approche doit permettre d'équilibrer le positionnement des échantillons sur l'ensemble du quadrant.

2.4. Analyse spatiale

40La dispersion de la distribution spatiale des lieux de sondages est évaluée à l'aide de l'outil Average Nearest Distance du logiciel ArcGis 9.3. Cet outil mesure, pour chaque point, la distance comprise entre ce point et le point voisin le plus proche dans une zone d'étude de dimension fixe. Le traitement des points est réalisé séparément pour chaque type de sol et chaque occupation des zones d'étude. Le mode de calcul des distances se base sur la distance euclidienne, distance en ligne droite à vol d'oiseau. Les coordonnées métriques X et Y correspondent au système de projection belge Lambert 72.

41La méthode repose sur la comparaison entre la moyenne de toutes les distances au plus proche voisin (D) et une moyenne calculée pour une distribution aléatoire hypothétique (d). Si D est significativement supérieure à d, les points sont considérés comme dispersés, si D est significativement inférieure à d, les points tendent à être groupés. Les p-values et les z-scores indiquent le niveau de signification statistique du test de l'hypothèse nulle du caractère aléatoire de la distribution des points dans l'espace, à savoir ni agrégée ni dispersée.

42Une analyse variographique a ensuite été réalisée pour les trois principaux types de sols retenus et les 3 occupations du sol dans les différentes zones d'étude. Six variables binaires ont été créées pour chaque point de prélèvement, correspondant à une fonction d'appartenance à la classe considérée (PTS ou occupation). Chaque point forestier s'est ainsi vu attribuer la valeur 1 à la variable « Forêt » et 0 aux variables « Culture » et « Prairie », et ainsi de suite pour chaque variable. Un variogramme est alors modélisé pour chacune de ces variables. Celui-ci exprime la continuité spatiale des différents facteurs et leur degré d'autocorrélation spatiale.

3. Résultats et discussion

43Le dispositif d'échantillonnage obtenu par cette méthodologie (Figure 3) est de type stratifié et raisonné. Stratifié car les points sont classés suivant l'importance de chaque facteur : type de sol, occupation et orientation. Raisonné car le positionnement des lieux de sondage a été choisi et pointé manuellement sur écran de manière équilibrée sur l'ensemble des quadrants et de la zone. Le nombre d'échantillons récoltés par campagne est variable, soit 196, 231 et 256 respectivement pour Prayon, Plombières et Sclaigneaux (Tableau 1). Cette différence quantitative est due à la fréquence de rencontre des entités « Sol-Occupation » par quadrant. Leur nombre a permis dans certains cas de répéter les sondages jusqu'à 14 fois. L'échantillonnage est donc proportionnel plutôt qu'équilibré, chaque entité étant échantillonnée en fonction de son importance surfacique. Lorsque des répétitions sont possibles, plusieurs positionnements sont prévus afin de répartir l'échantillonnage sur différentes gammes de distance à la source (Figure 4 et Tableau 2). Néanmoins, la répartition des sondages selon la distance n'est pas comparable entre les zones car elle dépend fortement des modèles de distribution des sols et de leurs occupations.

44Un des sols retenus est commun aux trois zones étudiées, il s'agit du sol limoneux à drainage naturel modéré ou imparfait (code 4020) (Tableau 1). Ce type de sol permettra une comparaison entre les zones. Les sols uniquement présents sur une zone reflètent les caractéristiques intrinsèques de celle-ci.

3.1. Relation entre type de sol et occupation du sol

45Les types de sols étudiés sur la zone de Sclaigneaux sont des sols limoneux dont le drainage naturel est « favorable » (code 4010) ou « modéré ou imparfait » (code 4020) et des sols limono-caillouteux à charge de silexite ou de gravier ou de conglomérat et à drainage naturel principalement favorable (code 7610) (Tableau 1). Ces trois types de sols représentent 67 % de la superficie de sols pouvant être choisis lors de la procédure de sélection (Tableau 3). C'est dans la partie centrale de la zone d'étude traversée par la Meuse que se rencontrent les sols à charge de silexite et de graviers issus des dépôts de terrasse et/ou d'anciennes alluvions de la Meuse. Ces dépôts sont très caillouteux et constitués de fragments rocheux aux arêtes émoussées (cailloux roulés) de nature très diverse (Oldenhove de Guertechin, 1973). Au nord du sillon, les plateaux sont principalement recouverts de sols limoneux à drainage variable selon leur disposition dans le relief. Cette partie nord de la zone correspond à la bordure sud de la région limoneuse également dénommée « Hesbaye ». Celle-ci est exploitée par de grandes cultures annuelles et des prairies en fonds humides pour laisser les quelques boisements sur les sols moins profonds et les sols de versants (Oldenhove de Guertechin, 1973). Par contre, la partie sud de la zone marque la bordure nord du plateau condruzien et l'apparition de sols gréseux à drainage variable entrecoupés de sols schisteux et de placages limoneux. Le mode d'exploitation agricole qui y est pratiqué se rapproche d'un type mixte agriculture-élevage (Oldenhove de Guertechin, 1973). Les pratiques agricoles de la zone sont centrées sur la culture et la prairie qui se partagent équitablement les surfaces arables (Feltz et al., 2010). Le relief mouvementé conditionne la présence de nombreux bois sur les pentes des ruisseaux affluents de la Meuse (Feltz et al., 2010). Cette répartition des occupations dans les paysages de la zone d'étude se traduit par un nombre de prélèvements élevé en cultures et en prairies avec, respectivement, 43 % et 35,5 % des prélèvements (Tableau 4). Les forêts occupent souvent les pentes escarpées des ruisseaux ou certaines crêtes sur des sols non représentatifs de la zone et ne représentent que 21,5 % des échantillons prélevés sur les trois sols sélectionnés.

46La zone de Prayon se situe entre le bocage occidental du pays de Herve (au nord), le massif de la Vesdre (au centre) et l'Ardenne liégeoise (au sud) (Pahaut, 1963). Les sols y sont très différents et très diversement distribués. De fait, la partie du bocage occidental regroupe des sols limoneux à drainage modéré ou imparfait et des sols à charges psammitiques et schisteuses. Le massif de la Vesdre comprend les sols de versants sur calcaire qui alternent avec des sols psammitiques, le tout entrecoupé de placages limoneux. Tandis que la partie située en Ardenne liégeoise est majoritairement gréseuse avec, sur ses plateaux, des sols limoneux peu caillouteux à drainage naturel imparfait à pauvre (Pahaut, 1963). Il existe une grande diversité pédologique sur cette zone mais les exigences de l'échantillonnage excluent les matériaux mal répartis. Le choix s'est orienté vers les sols limono-caillouteux à charge calcaire ou contenant du calcaire et à drainage naturel quasi exclusivement favorable (code 7510), les sols limono-caillouteux à charge psammitique ou schisto-psammitique et à drainage naturel principalement favorable (code 7410) et les sols de type limoneux à drainage naturel modéré ou imparfait (code 4020) (Tableau 1). Au niveau surfacique, ils représentent 47 % des sols pouvant faire partie de l'échantillonnage (Tableau 3). Le sol 4020 ne couvre que 10 % de la surface, mais est un des rares PTS à être présent dans tous les quadrants de manière homogène. D'un point de vue paysager, la zone de Prayon se situe dans la Basse Vesdre, partie la plus encaissée de la vallée avec un relief entaillant des versants abrupts et fermée par des fonds étroits (Godart et al., 2007a). Les occupations du sol sont influencées par ce relief et les propriétés des sols. On rencontre principalement des boisements, occupation dominante dans le paysage mais confinée sur les versants de la Vesdre et de ses affluents (Godart et al., 2007a). Ils n'occupent pas toujours les trois sols étudiés sur la zone et ne représentent que 31,2 % des sondages prévus (Tableau 4). L'occupation la plus échantillonnée est la prairie, soit 41,3 % des prélèvements. Le reste est échantillonné dans le peu de cultures mises en labour sur les anciennes couronnes de prairies entourant les villages. Les cultures ne pouvant s'accommoder des sols de faible aptitude agricole, acides et peu profonds des versants (Pahaut, 1963 ; Godart et al., 2007a).

47Les sols de Plombières sont des sols limoneux peu caillouteux relativement épais. La minorité de sols de type limono-caillouteux que la zone présente se situe sur les plateaux. Dans ce cas, le manteau mixte comprend des éléments du substrat et une couverture limoneuse plus ou moins épaisse (Pahaut, 1965). De fait, les trois sols sélectionnés par la méthodologie d'échantillonnage sont tous des sols limoneux dont le drainage naturel varie de « favorable » (code 4010) à « assez pauvre à très pauvre » (code 4030) en passant par « modéré ou imparfait » (code 4020) (Tableau 1). À eux trois, ces types de sols couvrent 65 % de la zone d'étude (Tableau 3). La zone est située dans l'aire paysagère de la Vallée de la Gueule herbagère (Godart et al., 2007b). On y rencontre de grands espaces dédiés au pâturage et entourés de haies, sans pour autant parler de bocage strict (Godart et al., 2007b). La prairie est l'occupation la plus représentée dans l'échantillonnage avec quasiment la moitié des prélèvements prévus (46,5 %) (Tableau 4). Entre les parcelles dédiées au pâturage, on retrouve des cultures qui seront échantillonnées à raison de 32,5 %. Les boisements sont souvent regroupés sur les sols des crêtes moins profonds (Godart et al., 2007b), d'où leur plus faible représentation dans la campagne d'échantillonnage (21,2 %).

3.2. Analyse spatiale

48Le traitement des points à l'aide de l'outil Average Nearest Distance montre que les distributions des lieux de prélèvement des trois zones d'étude ne peuvent pas être considérées comme aléatoires (Tableau 5). Le z-score mesure l'écart entre moyenne observée et moyenne attendue rapporté à l'écart-type. La p-value est la probabilité de se tromper en rejetant l'hypothèse du caractère aléatoire de la distribution spatiale des points. Sur base du seuil de 5 %, 7 tests sur les 18 réalisés sont compatibles avec une distribution aléatoire. On y retrouve les trois occupations du sol de Prayon, les prairies de Plombières et les PTS limoneux de Sclaigneaux et Plombières, à l'exception des sols limoneux à drainage naturel favorable de Plombières. Les autres modalités peuvent être considérées comme présentant des distributions spatiales en grappes, sans risque élevé de se tromper. C'est particulièrement le cas pour les cultures (Sclaigneaux et Plombières), les forêts de Plombières, les sols limoneux à drainage naturel favorable (Plombières) ou modéré ou imparfait (Prayon) et les sols limono-caillouteux à charge calcaire de Prayon.

49Les résultats de l'analyse variographique des variables liées à la nature des sols et leur occupation, exprimées sous forme d'un indice présence/absence du caractère étudié, sont donnés au tableau 6. Le rapport C/(C0+C) exprime la part de la variance qui présente une continuité spatiale et la portée, la distance au-delà de laquelle deux points ne sont plus corrélés spatialement.

50Deux cas de figure se présentent. Premièrement, la moitié des variables étudiées ne présente pas de continuité spatiale (effet pépite pur). Les PTS à Sclaigneaux et à Plombières et l'occupation à Prayon. Cela signifie qu'il y a autant de différences (de PTS ou d'occupation) entre deux points distants de 300 m qu'entre deux points distants de 3 000 m pour ces modalités. Cela témoigne d'un morcellement important de ces variables dans le plan d'échantillonnage. Deuxièmement, l'autre moitié des variables présente une structure de continuité spatiale modélisable par un modèle sphérique. La portée pour les PTS de Prayon est comprise entre 870 et 1 940 m. Les portées des cultures et prairies sont inférieures à celles des forêts dans les deux autres zones. La longueur des portées traduit la dimension des ensembles spatiaux caractérisés par une même modalité. Toutefois, cela peut correspondre à des réalités différentes en termes de dimensions des plages cartographiques, comme on peut le voir à la figure 3 pour les PTS de Prayon. Les sols limoneux et les sols à charge calcaire présentent des portées relativement similaires, mais là où les seconds sont présents sous forme de plages continues relativement importantes, les premiers présentent un morcellement marqué. L'existence d'une continuité spatiale dans ce cas est liée à l'absence des PTS voisins du 4020 dans le plan d'échantillonnage.

51L'analyse spatiale des motifs et de l'autocorrélation entre les points présente des éléments de convergence. Les points, quelle que soit l'occupation du sol, à Prayon et sur les cultures et les prairies de Plombières, apparaissent peu ou pas groupés et peuvent être considérés comme spatialement indépendants à très faible distance (moins de 500 m). Cela s'explique en partie par le fait que sur la zone d'étude de Prayon, les cultures et prairies sont relativement dispersées dans le paysage, contrairement aux forêts qui, malgré une grande superficie, se concentrent davantage dans les vallées. À Plombières, la prairie est le mode d'exploitation le plus courant et le mieux distribué sur l'ensemble de la zone, d'où un mode de répartition ni groupé ni dispersé. En pratique, cela signifie que pour l'interpolation future des résultats des teneurs en ETM, le facteur occupation du sol ne devrait jouer qu'un rôle limité pour les faibles gammes de distance pour les cultures et les prairies. Pour les forêts et les PTS à Prayon, l'influence des facteurs, si elle s'avérait significative, pourrait porter sur des distances plus grandes, de 1 500 à 2 000 m.

4. Conclusion et perspectives

52Les trois dispositifs d'échantillonnage obtenus par la stratégie de sélection des plages de sols développée dans ce papier permettront l'étude de l'influence des facteurs sol, occupation et direction des vents sur la distribution spatiale des ETM et leur répartition au sein des phases porteuses du sol, après analyse des différentes teneurs en laboratoire.

53Le schéma opéré repose sur un équilibre statistique de l'échantillonnage et sur une représentativité des facteurs sur l'ensemble des zones étudiées. Ce ne sont donc pas des critères comme la relation cout-efficacité ou les exigences des autorités en charge du projet qui ont dicté la stratégie, comme dans de nombreux programmes d'échantillonnages à large échelle en sols contaminés (Back, 2007).

54Pour conclure, la compilation des données de terrain et de laboratoire permettra d'évaluer les risques encourus par la valorisation agricole et forestière actuelle de ces terres plus ou moins contaminées, ce en vue de proposer des solutions adaptées selon le mode d'exploitation sans pour autant compromettre leur valorisation économique future.

Bibliographie

Notes

Pour citer cet article

A propos de : Amandine Liénard

A propos de : Laurent Bock

A propos de : Gilles Colinet