Introduction

Les charges critiques sont définies comme des seuils de contamination au-delà desquels des effets nocifs peuvent survenir sur les composantes de notre environnement, en particulier forestier [1]. Le respect de ces seuils a pour principal objectif de protéger l'équilibre à long terme des écosystèmes. Cette notion de charges critiques est la base d'évaluation des protocoles définis par la Commission économique européenne des Nations Unies (CEE-ONU) pour proposer des réductions d'émissions polluantes (SO₂ et NO_x, en particulier). Cette méthode a été mise en œuvre depuis la fin des années 80 dans les pays d'Europe du Nord.

En France, les nombreuses données nécessaires aux calculs de ces charges critiques ont été acquises depuis 1992 [2, 3] et sont progressivement mises en base de données [4-6]. La mise en œuvre depuis 1991 de réseaux d'observation des forêts [CEE, Réseau national de suivi à long terme des écosystèmes forestiers (RENECOFOR)] a aussi permis d'affiner et de valider les hypothèses cartographiques réalisées plus particulièrement vis-à-vis des charges critiques d'acidité [7].

Une fois les cartes de charges critiques réalisées, l'intérêt principal réside dans leur confrontation avec les dépôts actuels de polluants atmosphériques. Ceci permet ainsi de confirmer les espaces où des effets pourraient être constatés sur les peuplements forestiers.

Méthodes du calcul des charges critiques d'acidité

Le calcul des charges critiques d'acidité repose sur une estimation du flux de cations basiques libéré par l'altération des minéraux des sols et des roches et du drainage de protons et d'aluminium. Il est établi sur la base d'un bilan de masse à l'état d'équilibre (méthode SSMB,  Steady State Mass Balance [8]) selon l'équation suivante :

­

CL_AC = BC_W + Q.([H⁺]_crit + [Al⁺⁺⁺]_crit)

­

BC_W, correspond à l'altération en cations basiques ;

Q est le drainage ;

[H⁺]_crit  est la concentration critique de protons dans les eaux de drainage ; celle-ci correspond pour les sols forestiers français à 25 µeq/l, soit un pH de 4,5-4,6 [7] ;

[Al⁺⁺⁺]_crit = R_crit . (BC_d + BC_W - BC_u),

avec R_crit le rapport Al/BC = 1,2 (BC pour cations basiques, avec BC_d pour le dépôt atmosphérique de cations basiques, BC_u le prélèvement de cations basiques par les arbres forestiers [9]).

Cette équation suppose donc que les charges critiques d'acidité correspondent à la capacité de neutralisation de l'acidité par altération de la roche, de l'altérite et du sol (BC_W) augmentée d'un drainage acceptable de protons et d'aluminium. Ceci signifie qu'une acidification s'accompagnant d'une libération de protons et d'aluminium peut être tolérée tant qu'elle reste compatible avecl'équilibre écologique et l'état sanitaire des peuplements forestiers. Ce drainage acceptable est estimé d'une part avecles concentrations en [H⁺] critique et, d'autre part, selon le rapport Al/BC (voir ci-dessus) . Lorsque ces seuils sont dépassés, les conditions de la croissance racinaire des arbres forestiers deviennent défavorables. Les valeurs de ces seuils établis d'après de nombreuses données européennes [10] ont été vérifiées pour des peuplements forestiers vosgiens et à l'aide des données du réseau RENECOFOR (sous-réseau CATAENAT [Charge acide totale d'origine atmosphérique dans les écosystèmes naturels terrestres]) pour les solutions de sols.

Enfin les dépassements de ces charges critiques par les dépôts atmosphériques acides (EX pour « exceedance » en anglais) ont été calculés comme suit :

EX = S_d + N_d - BC_d - CL_AC

S_d correspondant au dépôt par voie atmosphérique des composés du soufre, N_d des composés de l'azote, BC_d des cations basiques, et CL_AC aux charges critiques d'acidité calculées.

Démarche pour la cartographie des charges critiques d'acidité

La cartographie a été réalisée en utilisant des modèles empiriques établis à partir des données de la carte des sols de France d'une part [11]. et sur la base des données acquises spécifiquement pour ce travail d'autre part [7]. Elles ont ensuite été validées à l'aide des données disponibles sur les réseaux existants (CEE et RENECOFOR) . Pour cela, plusieurs étapes ont été nécessaires :

1. Une carte « pédogéochimique » des sols de France (Figure 1, ci-contre) prenant en compte, outre la pédogénèse, la géochimie des matériaux constitutifs des sols, a d'abord été construite avecla base de données de la carte des sols de France au 1/1 000 000 [11].

Figure 1. Matériaux géologiques et sols du territoire français (carte construite d'après les données de la carte pédologique de France, INRA, 1998. (In : Party, 1999 [7]).
1 à 4 :sols issus de formations sédimentaires meubles et calcaires (ex.dunes, alluvions et conglomérats calcaires, loess).
5 à 16 : sols issus de formations sédimentaires meubles et acides (ex. sables des Landes, de Fontainebleau, formations à silex, lehms, alluvions anciennes ou récentes, formations glaciaires acides).
17 à 20 :sols issus de formations sédimentaires dures et calcaires (ex. calcaires, marnes, certains flyschs).
21 à 24 :sols issus de formations sédimentaires dures et acides (ex. schistes et grès).
25 à 29 :sols issus de formations métamorphiques (ex. micaschistes, gneiss, granitoïdes).
30 et 31 :sols issus de formations volcaniques (ex. rhyolites, tufs, basaltes).
Parent materials and soils of France (based on the pedological map of France, INRA, 1998. (ln :Party 1999 [7]).

2. Les valeurs d'altération ont ensuite été calculées sur 12 sites de référence représentatifs des grandes familles pédogéochimiques définies, à l'aide du modèle PROFILE [12] sur la base d'analyses chimiques totales et d'une minéralogie reconstituée à l'aide de deux méthodes normatives : UPPSALA model [13, 14] et la mésonorme de Barth (Barth [15], adaptée par Stussi [16, 17]).

3. Une relation observée entre les valeurs d'altération calculées, la granulométrie et la somme des cations basiques mesurés dans les sols des 12 sites de référence a été établie. Ce modèle empirique a été ensuite appliqué aux 102 placettes du réseau RENECOFOR.

4. Les valeurs d'altération calculées (5e centile) en utilisant la mésonorme de Barth, ont ensuite été extrapolées aux unités de la carte « pédogéochimique » de France. Ceci a permis d'obtenir une carte des valeurs d'altération pour l'ensemble du territoire français (Figure 2).

Figure 2. Classes d'altération des matériaux géologiques et des sols du territoire français, calculées selon le modèle PROFILE (Party, 1999 [7]).
Weathering rates for French parent materials and soils, calculated using the PROFILE model (Party, 1999 [7]) (keq.ha^-l.an^-1).

5. Le drainage Q, nécessaire au calcul des charges critiques d'acidité, a été cartographié à l'aide de la carte des pluies efficaces pour le territoire français [18], pondérée par le calcul des réserves utiles selon les profondeurs d'enracinement constatées pour les principales essences forestières (données du réseau RENECOFOR).

6. Les charges critiques d'acidité ont été calculées sur les 12 sites de référence, puis une relation établie entre les valeurs d'altération et les charges critiques calculées a été appliquée aux 102 placettes du réseau RENECOFOR.

7. Les charges critiques d'acidité pour les sols forestiers ont ensuite été cartographiées par modélisation de la relation établie au point 6 en utilisant les cartes de l'altération et du drainage (Figure 3, ci-contre). Pour apprécier l'influence des méthodes de reconstitution minéralogique qui sont à la base du calcul des valeurs d'altération, les charges critiques ont été cartographiées en utilisant les deux méthodes normatives citées plus haut au point 2.

Figure 3. Charges critiques d'acidité pour les sols forestiers du territoire français (Party, 1999 [7]).
Critical loads of acidity for French forest soils (Party, 1999 [7]) (keq.ha^-1.an^-1).

8. Enfin, les calculs menés sur les données des 102 placettes du réseau RENECOFOR ont permis de définir les sols, les essences forestières et les régions qui peuvent être concernées par les dépassements de charges critiques et de leur affecter un niveau de risque. Cette approche a été confirmée à l'aide des données des 540 placettes du réseau CEE [19] d'une part, et du résultat des calculs de charges critiques également réalisés pour les eaux de ruisseaux d'autre part [20]. Une étude récente concernant les dépôts atmosphériques à l'échelle de la France [21] permettra sans doute à court terme une spatialisation de ces dépassements à cette échelle.

Résultats obtenus

La carte pédogéochimique établie (Figure 1, p. 521) est fondamentale. En effet, elle permet de représenter les différentes familles du continuum roche-altérite-sol, et de leur associer un potentiel d'altération chiffré à l'aide d'un ou plusieurs profils types. C'est ce qui a été fait à l'aide des profils des 102 placettes du réseau RENECOFOR. Une carte des classes de valeurs d'altération (Figure 2, p. 522) peut ensuite en être dérivée. En ce sens, elle représente une amélioration quantifiée de la carte pédogéochimique présentée par Pedro et Sherer [22].

Quelle que soit la méthode de reconstitution minéralogique utilisée, les valeurs d'altération dépassent 2,0 keq/ha^-1.an^-1 pour environ 40 % du territoire français (roches carbonatées ou basiques) alors que 25 % du territoire présentent des valeurs inférieures à 0,5 keq.ha^-1.an^-1 (Figure 2, p. 522). Pour les charges critiques (Figure 3, p. 523), les valeurs supérieures à 2,0 keq.ha^-1.an^-1 concernent environ 60 % du territoire et moins de 10 % pour les valeurs inférieures à 0,5 keq.ha^-1.an^-1.

Les charges critiques d'acidité ont été calculées pour les forêts de production (environ 130 000 km² selon l'inventaire forestier national), les prairies permanentes principalement adjacentes aux espaces forestiers (environ 50 000 km²,en fait les communes dont le rapport surface toujours en herbe/surface agricole utilisée, STH/SAU est supérieur à 85 %), ainsi que les formations boisées non productives (environ 15 000 km² selon l'inventaire forestier national). Dans l'espace forestier de production, les charges critiques faibles, inférieures à 1,0 keq.ha^-1.an^-1, représentent un peu moins de 15 % des surfaces, soit environ 17 000 km². Les surfaces prairiales ne présentent pratiquement pas de surfaces avecdes charges critiques d'acidité < 1,0 keq.ha^-1.an^-1.

Les surfaces qui présentent les valeurs de charges critiques les plus faibles (inférieures à 0,5 keq.ha^-1.an^-1) correspondent surtout à cinq grands types pédogéochimiques :

• les sols bruns acides et podzolisés sur grès et granites acides, notamment dans les Vosges ;

• les sols bruns acides (plus rarement podzolisés) des granites les plus acides du Massif central ;

• les sols bruns acides (parfois podzolisés) des granites et surtout des schistes de l'ouest de la France (Vendée, Bretagne et Normandie) ;

• les sols bruns acides des schistes acides des Ardennes ;

• les sols bruns acides, lessivés et podzolisés des sables des Landes et des formations résiduelles à silex des vallées de la Seine et de la Loire, cet ensemble paraissant le plus sensible (avec le grès des Vosges) avecdes valeurs de charges critiques inférieures à 0,2 keq.ha^-1.an^-1.

Ceci est vrai quel que soit le modèle minéralogique utilisé. Dans les deux modélisations testées, les valeurs de charges critiques obtenues sont similaires pour 95 % des surfaces (± 0,2 keq.ha^-1.an^-1). Pour les 5 % restants (soit près de 9 000 km²), les différences observées sont plus importantes. Elles concernent pour 3 % des surfaces dont les charges critiques sont supérieures à 2,0 keq.ha^-1.an^-1 et pour 2 % la classe 1,0-2,0 keq.ha^-1.an^-1. Ce chiffre de 2 % de surfaces françaises concernées peut paraître faible, mais représente en fait un peu plus de 3 000 km², soit une surface équivalente, par exemple , à toutes les plantations de douglas en France ! Ces surfaces correspondent à des sols lessivés hydromorphes sur alluvions anciennes ou récentes acides et à des formations résiduelles à silex (unités 11, 13 et 16, Figure 1, p. 521).

Les calculs de dépassement des charges critiques d'acidité pour les 102 sites du réseau RENECOFOR montrent :

• 18 placettes pour lesquelles les dépôts atmosphériques acides dépassent actuellement les charges critiques ;

• 18 autres placettes pour lesquelles les dépôts atmosphériques acides sont très voisins des charges critiques calculées, 10 d'entre elles présentant un risque potentiel fort (EX < à 0,2 keq.ha^-1.an^-1), 8 autres présentant un risque potentiel modéré (EX compris entre 0,2 et 0,5 keq.ha^-1.an^-1).

On ne peut cependant déduire trop rapidement de ce décompte qu'un tiers des surfaces forestières françaises est concerné. En effet, les sols du réseau RENECOFOR comportent majoritairement des sols forestiers acides et s'écartent ainsi de la moyenne des sols forestiers français, qui en comportent proportionnellement deux fois moins environ [23]. Ceci est surtout lié à des régions non prises en compte par le réseau RENECOFOR.

Néanmoins, les 18 placettes présentant des dépassements de charges critiques permettent de situer les grands ensembles de matériaux concernés. Ce sont principalement :

• les sables acides (unités 5 à 10, Figure 1, p. 521) comportant surtout les sables du Tertiaire et les formations à silex à recouvrement sableux pour 7 placettes ;

• les schistes (1 placette) et les grès (6 placettes) (unités 22 et 23, Figure 1), soit les sols acides sur schistes du Cambrien, les schistes gréseux et siliceux, plus ponctuellement les grès argileux ;

• les granitoïdes (unités 26 à 28, Figure 1), soit des granites plus ou moins acides (3 placettes) et certains gneiss et migmatites (1 placette).

Ces placettes concernent 8 grands ensembles régionaux . Il s'agit dans l'ordre des dépassements ou des risques de dépassement de charges critiques acides des ensembles suivants :

• les Ardennes [24] ;

• les Vosges, l'Alsace et la Moselle ;

• le sud-ouest de l'Ile-de-France (Yvelines, Fontainebleau), plus ponctuellement la Champagne ;
  ces 3 ensembles présentant 10 des 18 cas de dépassement ;

• la vallée de la Loire, plus ponctuellement la vallée de la Seine ;

• le rebord Est du Massif central, plus ponctuellement le rebord Ouest (Limousin, Haute-Vienne) ;

• les Landes de Gascogne ;

• la Bretagne et la Normandie ;
  ces 4 unités présentant 7 cas de dépassement et 16 des 18 cas à risque de dépassement ;

• les Pyrénées, très ponctuellement (1 placette à risque modéré).

Enfin, par rapport aux peuplements forestiers, dans le réseau RENECOFOR, toutes les essences (sauf le mélèze) paraissent concernées, mais 3 d'entre elles le sont plus particulièrement :

• le pin sylvestre d'abord, avec 10 placettes sur les 14 existantes présentant un dépassement ou un risque de dépassement ; on peut aussi y associer les plantations de pin maritime et de pin Laricio, en dehors de leur aire écologique d'origine ;

• l'épicéa avec 6 placettes sur 10 existantes présentant un dépassement ou un risque de dépassement ;on peut lui associer aussi les plantations de douglas, mais à moindre titre, et plus rarement encore le sapin ;

• le chêne sessile, avec 6 placettes sur 21 existantes présentant un dépassement ou un risque de dépassement.

Deux autres ensembles présentent encore des risques de dépassement de charges critiques d'acidité :

• le pin maritime dans les Landes tout d'abord, avec 3 placettes à risque potentiel fort (sur 7 existantes, dans le réseau RENECOFOR) ;

• le hêtre en Normandie, dans l'est de la France ou le Massif central, avec 5 placettes à risque modéré (sur 21 existantes, dans le réseau).

Les valeurs de charges critiques d'acidité calculées indiquent que la plupart des écosystèmes ne sont donc pas tous en mesure de supporter les apports atmosphériques acides actuels. Les cartes présentées ne correspondent pas à une cartographie de cas avérés, mais sont le reflet de risques potentiels.

Conclusion et perspectives

Les charges critiques d'acidité modélisées à l'aide d'un système d'information géographique sur l'ensemble du territoire français à l'échelle du 1/1 000 000 montrent qu'un peu moins de 15 % des surfaces forestières de production (soit 17 000 km²) sont concernées par des valeurs de charges critiques inférieures à 1 keq.ha^-1.an^-1 (Figure 4). Ceci concerne surtout les plantations de pin sylvestre et d'épicéa sur limons et sables acides développés sur grès, granites et schistes du nord de la France [Vosges, Ardennes et plus ponctuellement sud-ouest de l'Ile-de-France (vers Fontainebleau), Massif central, Bretagne et Normandie].

Figure 4. Modélisation des charges critiques d'acidité à partir de 12 sites de référence représentatifs des grandes familles pédogéochimiques et extrapolation des résultats aux 102 placettes du réseau RENECOFOR.
Critical loads of acidity modelled for 12 reference sites and extended to the RENECOFOR Network (102 sites).

Ces résultats obtenus sur sites de référence représentatifs des écosystèmes forestiers français [7], ont été validés à l'aide des données du réseau RENECOFOR et extrapolés à l'ensemble de la France à l'aide des cartes numériques existantes (pédologique, géologique, pluies efficaces, occupation forestière et prairiale des sols). Ils sont également confirmés par les données des 540 placettes du réseau CEE qui ont été complétées en 1999 d'environ 60 placettes dont 20 dans des situations à risque de dépassement des charges critiques en Bretagne, Normandie et sur le rebord Est du Massif central, sous pin sylvestre, pin maritime, épicéa, douglas, chêne sessile ou hêtre, et sur matériaux sensibles, schistes et grès ou granitoïdes.

Enfin, les calculs de charges critiques réalisés pour les eaux de ruisseaux concordent également dans plus de 80 % des cas avec les charges critiques cartographiées pour les sols forestiers [7]. Une cartographie fine des dépassements de charges critiques constatés sur les 102 placettes du réseau RENECOFOR pourrait être prochainement effectuée à l'aide de la carte des dépôts atmosphériques réalisée récemment [21].

Ainsi, dans les sites identifiés comme sensibles, trois voies sont possibles pour conserver à terme l'intégrité de ces écosystèmes forestiers :

• continuer d'agir à court et moyen terme sur les sources d'émissions pour diminuer la part des composés acidifiants dans les dépôts atmosphériques ;

• agir sur les choix de gestion forestière à long terme, notamment dans les choix d'essence et de traitement forestier ;

• agir à court terme par des chaulages raisonnés et mesurés. Cette disposition ne doit être réservée qu'aux cas les plus sensibles dès leur origine et/ou ayant des dépassements de charges critiques sur de grandes surfaces (Vosges et Ardennes, principalement).

En particulier, les choix d'essences pourraient être réalisés non seulement sur la base des possibilités actuelles offertes par le milieu, ce que permet actuellement l'analyse des potentialités des stations forestières, mais aussi en fonction des conséquences à long terme de la production d'une essence forestière donnée sur le continuum roche-altérite-sol. En effet, une forte capacité de production forestière sur un site minéralogiquement pauvre provoque des pertes susceptibles d'épuiser en quelques dizaines d'années les réserves minérales d'un écosystème. C'est par exemple le cas des pessières sur grès des régions du nord-est de la France [25].

À terme, pour définir les priorités infrarégionales sur le sujet, il sera sans doute nécessaire de cartographier ces résultats à l'échelle du 1/250 000 [26-28], ce que devrait permettre des cartes de sols réalisées dans la plupart des régions dans le cadre du programme Inventaire Gestion et Conservation des Sols (IGCS) du ministère de l'Agriculture au cours de la prochaine décennie.

Enfin, une première évaluation des charges critiques azotées à l'échelle de la France complète les charges critiques d'acidité sur la base d'une cartographie numérique de la végétation.

Les travaux sur les charges critiques ont pour l'essentiel été financés par l'ADEME (Agence de l'environnement et de la maîtrise de l'énergie) et le ministère de l'Environnement.