Cas : dégradation pesticides

Ce cas illustre la simulation d'une contamination de nappe par des pesticides et l'exploitation du couplage entre concentration, métabolites et signature isotopique CSIA pour distinguer ce qui relève de la dégradation active, de la dilution et de l'adsorption. La question pratique : pourquoi la concentration en atrazine baisse-t-elle plus vite que ne le suggère la décroissance attendue, et cette baisse reflète-t-elle une dégradation réelle ou un simple étalement du panache dans un volume plus grand.

Contexte

Nappe alluviale en zone agricole. L'atrazine, interdite d'usage en France depuis 2003, est encore régulièrement détectée dans plusieurs forages de contrôle du bassin, souvent accompagnée de ses métabolites DEA et DIA. Les concentrations en atrazine décroissent globalement année après année, mais à un rythme variable d'un piézomètre à l'autre. La question qui se pose est celle de la nature de cette décroissance : s'agit-il d'une véritable atténuation naturelle par dégradation microbienne, ou d'un simple transport et dilution qui déplacent le problème ailleurs sans réellement le résoudre.

Trois piézomètres du réseau sont équipés pour la CSIA sur atrazine (δ¹³C et δ¹⁵N) depuis trois ans. Les mesures isotopiques sont disponibles et peuvent être confrontées à la simulation.

Mise en place de la scène

Élément de scène	Configuration
Emprise horizontale	2000 × 800 m, zone agricole alluviale
Profondeur	12 m, référence niveau du sol
Grille de calcul	100 × 40 × 24 = 96 000 mailles
Couche superficielle (0 à 2 m)	Sol agricole, matière organique 1,5 à 3 %, pH 6,5 à 7,5
Nappe alluviale (2 à 10 m)	Alluvions graveleuses, porosité 0,20, K = 5e-4 m/s
Substratum (10 à 12 m)	Argile marneuse, porosité 0,05, K = 5e-9 m/s
Points d'échantillonnage	12 piézomètres, dont 3 équipés CSIA

Composés et paramètres

La simulation suit l'atrazine et ses deux principaux métabolites DEA et DIA. Trois voies de dégradation sont activées dans le moteur, chacune produisant un métabolite dominant différent et chacune ayant son epsilon carbone caractéristique.

Voie	Métabolite dominant	Epsilon δ¹³C (‰)	Conditions favorables
Désalkylation aérobie microbienne	DEA	-2 à -5	O₂ dissous > 1 mg/L, matière organique disponible
Photodégradation superficielle	DEA et DIA	proche de 0	Zones peu profondes, < 0,5 m sous la surface
Hydrolyse abiotique	Hydroxyatrazine	très faible, < 1	pH acide ou basique

Le taux de dégradation k varie maille par maille selon les conditions géochimiques interpolées. Dans les zones aérobies riches en matière organique, la voie microbienne domine avec un k de l'ordre de 0,003 /jour. Dans les zones anoxiques profondes, la dégradation est beaucoup plus lente avec un k de l'ordre de 0,0001 /jour. Le moteur propage simultanément la concentration et la signature isotopique selon l'equation de Rayleigh.

Contrairement à une simulation classique où un k global est imposé, ici le k effectif émerge du couplage avec les conditions géochimiques locales. Cela implique que si l'utilisateur décrit correctement la géochimie spatiale du site, le modèle produit naturellement une hétérogénéité de dégradation qui explique les différences observées entre piézomètres.

Source

La source est modélisée comme une source diffuse surfacique active de 1970 à 2003, reflétant les applications agricoles historiques de l'atrazine. Après 2003, l'apport direct s'arrête mais une désorption lente depuis le sol continue à alimenter la nappe pendant encore plusieurs années. L'intensité de la source réflète l'usage agricole typique de la région, calibré pour reproduire l'ordre de grandeur des concentrations mesurées au début du suivi.

Résultats

La simulation sur 55 ans produit un panache globalement homogène à l'échelle du site mais avec des variations locales marquées. Les concentrations d'atrazine simulées aux 12 piézomètres reproduisent les décroissances observées, y compris les différences entre piézomètres : certains montrent une décroissance rapide, d'autres plus lente. La différenciation est cohérente avec les conditions géochimiques locales interpolées sur la grille.

Le ratio atrazine/DEA simulé est particulièrement informatif. Dans les zones où la dégradation aérobie est active, le ratio descend à des valeurs inférieures à 1 avec un DEA majoritaire. Dans les zones peu dégradantes, le ratio reste proche ou supérieur à 1 même après plusieurs décennies. Cette différenciation spatiale est une empreinte directe de l'hétérogénéité de la dégradation.

Piézomètre	Zone géochimique	Atrazine mesurée (ng/L)	DEA mesuré (ng/L)	δ¹³C atrazine mesuré (‰)	δ¹³C simulé (‰)
PZ-03	Aérobie, MO élevée	80	220	-22,1	-22,5
PZ-07	Intermédiaire	180	150	-25,3	-25,0
PZ-11	Profond, peu oxygéné	320	40	-27,4	-27,2

Interprétation

La reproduction simultanée des concentrations, des ratios parent/métabolite et des signatures isotopiques aux trois piézomètres équipés en CSIA valide l'hypothèse centrale : la décroissance observée de l'atrazine dans la nappe n'est pas seulement due au transport et à la dilution, mais reflète une dégradation microbienne réellement active dans certaines zones du site. L'enrichissement isotopique mesuré ne peut pas être expliqué par un simple mélange ni par une adsorption préférentielle ; il est la signature directe d'un fractionnement cinétique associé à la dégradation.

Cette conclusion a une portée pratique : dans les zones aérobies du site, l'atténuation naturelle est effective et la concentration continuera de décroître. Dans les zones moins oxygénées (PZ-11 et sa zone géochimique), l'atrazine sera persistante à l'échelle de plusieurs décennies supplémentaires. Les zones qui demandent une surveillance renforcée sont donc clairement identifiables, et les zones qui s'auto-réhabilitent peuvent être laissées à leur évolution naturelle.

Sans la CSIA, la simulation aurait probablement ajusté un k global uniforme calé sur la baisse moyenne des concentrations. Le résultat aurait paru correct en moyenne mais aurait manqué l'hétérogénéité spatiale qui est précisément l'information la plus utile pour la gestion du site. L'intégration CSIA et simulation est donc bien plus qu'une vérification : c'est un levier d'interprétation qui oriente les décisions opérationnelles.

Prolongements

Variante « source persistante » : simuler le cas où la désorption depuis le sol est plus longue que prévu, pour évaluer la durée pendant laquelle la nappe restera contaminée même sans nouvel apport.
Ajout de simazine et terbuthylazine : étendre la simulation aux autres triazines du panel, pour voir si elles présentent les mêmes hétérogénéités ou si la persistance varie selon la molécule.
Intégration du δ¹⁵N sur atrazine : ajouter la dimension azote pour distinguer la désalkylation (fractionne surtout C) de l'hydrolyse (fractionne surtout N), utile si les conditions pH varient.
Scénario avec apport ponctuel : tester l'hypothèse que certaines concentrations élevées soient dues à des apports plus récents que l'arrêt officiel, plutôt qu'à la seule mémoire du sol.

Pour aller plus loin

Pesticides : famille triazines, métabolites, CSIA.
Voies de dégradation : chaînes métabolites atrazine.
Spéciation et propagation : cinétique variable spatialement.
Isotopie CSIA : principes de l'enrichissement résiduel.