Spéciation redox

La spéciation redox décrit la répartition d'un élément entre ses différents états d'oxydation en solution ou dans une matrice solide. Pour les éléments redox actifs (Cr, Fe, Se, As, Sb, Mo, Mn, U), cette répartition contrôle la mobilité, la toxicité, l'adsorption et le fractionnement isotopique. IsoFind résout la spéciation à partir des conditions physico-chimiques (Eh, pH, O₂ dissous) via un modèle couplé entre estimation thermodynamique (diagramme de Pourbaix) et prédiction par apprentissage machine du bridge Nexus. Cette page présente le modèle interne, les seuils codés et les systèmes supportés.

Les deux approches résolues par IsoFind

Le moteur combine deux niveaux de prédiction. La prédiction par apprentissage machine, quand elle est disponible via le bridge Nexus, utilise des modèles entraînés sur des bases de spéciation thermodynamique étendues. Le fallback thermodynamique reste actif quand Nexus est indisponible ou quand les conditions fournies sont insuffisantes.

Approche	Moteur	Entrées	Source dans le retour
ML Nexus	Bridge /api/nexus/predict_speciation	Élément, concentration, pH, Eh, DO, matrice	source = 'ml_speciation'
Thermodynamique	Fonction _estimate_cr_speciation_from_eh	pH, Eh uniquement	source = 'thermodynamic'
Fallback Rayleigh	Fonction fraction_reduite_depuis_redox	Eh, DO, profondeur (estimations automatiques)	Fraction réduite et epsilon selon le régime

Modèle thermodynamique pour Cr(VI) / Cr(III)

Le cas du chrome est le mieux codifié dans IsoFind. Le Cr(VI), oxydé et mobile (chromate CrO₄²⁻), est le contaminant d'intérêt, tandis que le Cr(III) réduit est nettement moins mobile et moins toxique. La conversion Cr(VI) → Cr(III) est donc un objectif de remédiation naturelle ou active. Le modèle thermodynamique d'IsoFind repose sur le diagramme pe-pH simplifié du chrome.

            pe = Eh / 59,16 (conversion mV → pe à 25 °C)

            pe_frontière = 23,0 − 1,33 × pH

            f_Cr(VI) = 1 / (1 + exp(−0,8 × (pe − pe_frontière)))

Cette formulation donne une courbe sigmoïde de la fraction Cr(VI) centrée sur la frontière thermodynamique Cr(VI)/Cr(III). La pente de 0,8 dans la sigmoïde représente l'incertitude intrinsèque autour de la frontière : aucune transition nette n'est réaliste dans un système naturel où coexistent plusieurs phases minérales et espèces complexantes.

Les valeurs numériques (pe° = 23,0 ; pente de pH −1,33 ; facteur de sigmoïde 0,8) sont codées en dur dans la fonction _estimate_cr_speciation_from_eh de prediction_routes.py. Elles correspondent à une approximation du diagramme de Pourbaix du Cr en conditions de dilution modérée, à 25 °C, sans ajustement pour la force ionique. Pour un travail critique, la prédiction Nexus ML est préférable parce qu'elle intègre plus de variables.

Seuils codés pour la réduction

Indépendamment du modèle thermodynamique pur, IsoFind utilise une cascade de seuils Eh et O₂ pour estimer la fraction réduite et sélectionner l'epsilon isotopique approprié. Ces seuils sont codés comme constantes et guident le modèle de Rayleigh.

Constante	Valeur	Signification
EH_SEUIL_OXIQUE	+350 mV	Au-dessus : conditions oxiques, pas de réduction significative
EH_SEUIL_REDUCTION_PARTIEL	+200 mV	Au-dessous : réduction partielle possible
EH_SEUIL_REDUCTION_FORT	−50 mV	Au-dessous : réduction complète probable
DO_SEUIL_ANOXIQUE	0,5 mg/L	Au-dessous : conditions anoxiques strictes

Les trois régimes de réduction

La fonction fraction_reduite_depuis_redox classe chaque échantillon dans un des trois régimes selon Eh et DO, puis applique l'epsilon approprié. Le choix de l'epsilon traduit le fait que la voie biologique (bactéries sulfato-réductrices) et la voie abiotique (Fe(II) aqueux, surfaces minérales) fractionnent différemment le chrome.

Régime	Conditions	ε ⁵³Cr (‰)	Référence bibliographique
Biologique	DO < 0,5 mg/L et Eh < +200 mV	−1,5	Bain & Bullen 2005, Zink et al. 2010
Abiotique	Eh < +200 mV mais DO non anoxique	−3,5	Zink et al. 2010 (réduction par Fe(II), surfaces minérales)
Mixte / oxique	Eh > +200 mV ou transition	−2,5 (valeur par défaut)	Valeur consensus si données insuffisantes

Inhibition par l'oxygène dissous

L'O₂ dissous inhibe la réduction du Cr(VI) même quand le potentiel Eh est modérément favorable. Le facteur d'inhibition est codé comme suit dans la fonction fraction_reduite_depuis_redox.

            do_factor = max(0 ; 1 − DO / 3)
        

Interprétation pratique : DO = 0 mg/L donne un facteur 1 (pas d'inhibition) ; DO = 1 mg/L donne 0,67 ; DO = 2 mg/L donne 0,33 ; DO = 3 mg/L donne 0 (inhibition totale). Cette correction empirique compense le fait que l'Eh mesuré à la sonde ne capte pas toujours correctement l'O₂ résiduel, et que la réduction du Cr(VI) est cinétiquement ralentie en présence d'oxygène même si la thermodynamique la permettrait.

Fractionnement isotopique du résiduel

Quand une fraction du Cr(VI) est réduite en Cr(III), la fraction restante de Cr(VI) est enrichie en ⁵³Cr selon l'équation de Rayleigh. C'est le fondement de l'utilisation du δ⁵³Cr comme preuve de dégradation réductrice active.

            δ⁵³Cr(t) = δ₀ + ε × ln(f)

            avec f = fraction de Cr(VI) non encore réduite

Pour un epsilon biotique de −1,5 ‰, une réduction de 50 % du Cr(VI) enrichit le résiduel de +1 ‰ ; une réduction de 90 % l'enrichit de +3,5 ‰. Ces ordres de grandeur sont détectables en routine par MC-ICP-MS. Le moteur calcule aussi le δ⁵³Cr du Cr(III) cumulé qui intègre l'ensemble du produit réduit : δ_produit_cumulé = δ₀ − ε × (1−f) / f × ln(f).

Autres systèmes redox importants

Au-delà du chrome, plusieurs systèmes redox méritent une attention particulière. Les principes sont analogues mais les seuils thermodynamiques, les epsilons et les contextes d'usage diffèrent.

Fe(II) / Fe(III)

Le fer est le plus abondant des éléments redox actifs dans les aquifères. Sa spéciation contrôle le cycle redox local et les couples secondaires (Cr, As, Mn). Fe(II) est soluble, Fe(III) précipite sous forme d'oxyhydroxydes peu mobiles à pH neutre.

Paramètre	Valeur
ε ⁵⁶Fe abiotique (réduction)	−1,5 ‰
ε ⁵⁶Fe biologique (réduction dissimilatoire)	−2,0 à −3,0 ‰
Frontière Pourbaix à pH 7	Eh ≈ +220 mV (précipitation Fe(OH)₃)
Anomalie Fe(II) en eau oxique	Indique apport anthropique ou couche réductrice déconnectée

As(III) / As(V)

L'arsenic est présent sous deux formes principales qui portent des toxicités différentes. L'arsénite As(III) est plus mobile et plus toxique que l'arséniate As(V). La frontière redox est proche de celle du Cr, autour de Eh ≈ +100 à +200 mV à pH 7.

Paramètre	Valeur
ε ⁷⁵As abiotique	−2,0 ‰
ε ⁷⁵As biologique	−3,0 ‰
Diagnostic par δ⁷⁵As	Réduction As(V) → As(III) en zone anoxique
Méthylation (As(III) → MMA, DMA)	Bactéries sulfato-réductrices, fractionnement additionnel

Se(VI), Se(IV), Se(0), Se(−II)

Le sélénium a quatre états d'oxydation stables, ce qui en fait un des systèmes les plus riches pour la forensique redox. Se(VI) et Se(IV) sont solubles ; Se(0) natif précipite ; Se(−II) forme des séléniures de métaux.

État	Nom	Mobilité	Indicateur
Se(VI)	Sélénate SeO₄²⁻	Très mobile	Eaux oxiques, irrigation
Se(IV)	Sélénite SeO₃²⁻	Mobile avec adsorption	Sols neutres à acides
Se(0)	Sélénium natif	Peu mobile (précipité)	Zone réductrice, produit bactérien
Se(−II)	Séléniure	Immobile (séléniure métallique)	Conditions très réductrices sulfurées

Epsilon biologique de Se (−5 ‰) est le plus grand de la table _EPSILON_TABLE pour les éléments redox, ce qui fait du δ⁸²Se un marqueur particulièrement sensible de la réduction microbienne.

Éléments sans redox en solution

Rappel : onze éléments du catalogue IsoFind restent sous une forme unique en solution aqueuse et ne supportent pas de spéciation dynamique : Pb, Zn, Cu, Cd, Ni, Co, Sr, Ca, Na, K, Mg. Pour ces éléments, la vue Géochimie rapporte uniquement la concentration totale. L'adsorption et le transport restent modélisés mais sans branche de spéciation.

Articulation avec l'adsorption et le fractionnement

La spéciation redox d'un élément est la porte d'entrée de plusieurs comportements environnementaux. Le moteur IsoFind propage la fraction réduite vers ces modules secondaires.

Couplage	Mécanisme
Adsorption différentielle	Cr(III) adsorbé préférentiellement sur oxydes Fe(III), Cr(VI) reste mobile
Correction epsilon par adsorption	Correction = favorability × f_reduced × ε_adsorption
Précipitation	Fe(III) précipite en Fe(OH)₃, Se(VI)→Se(IV)→Se(0) avec précipitation terminale
Couplage entre éléments	Fe(II) réduit As(V) et Cr(VI), cascade redox secondaire

Visualisation dans IsoFind

Les résultats de spéciation sont intégrés à plusieurs visualisations d'IsoFind. L'utilisateur retrouve la spéciation selon le contexte de consultation.

Vue	Affichage de la spéciation
Onglet Géochimie d'un échantillon	Badge Cr(VI) / Cr(III) avec pourcentages si pH et Eh renseignés
Simulation 3D	Mode de rendu « Simulation » colorie par fraction réduite
Nexus Quick Match Process	Estimation de f_reduced à partir des conditions saisies
Bloc de rapport Spéciation	Tableau des pourcentages par élément et par échantillon

Accès via l'API

Endpoint	Usage
POST /api/predict/speciation	Spéciation redox d'un élément à partir de pH, Eh, DO
POST /api/predict/infer	Inférence complète avec spéciation + fractionnement + adsorption
POST /api/nexus/predict_speciation	Appel direct au bridge Nexus ML (si disponible)

Limites et bonnes pratiques

La spéciation prédite est une estimation thermodynamique ou ML. Elle ne remplace pas une mesure directe de spéciation par HPLC-ICP-MS, méthode de référence quand le dossier l'exige.
Les valeurs d'Eh mesurées en puits sont souvent peu fiables (contamination par O₂, cinétique lente des électrodes). Prioriser les indicateurs chimiques (Fe²⁺ mesuré, sulfure, NO₃⁻ disparu) quand disponibles.
Les seuils codés (+350 oxique, +200 partiel, −50 fort) sont des repères convenables en aquifères neutres. Sur des milieux très acides (sols sulfatés) ou très alcalins (aquifères carbonatés), une calibration locale améliore la prédiction.
La cinétique n'est pas prise en compte dans le modèle thermodynamique : une eau peut rester transitoirement dans un état hors équilibre (par exemple Cr(VI) résiduel en zone réductrice pendant des mois). Le bridge Nexus ML intègre partiellement cet aspect via ses données d'entraînement.
Pour les dossiers contentieux, conserver une trace écrite des conditions physico-chimiques ayant servi à la prédiction. Le bloc de rapport Spéciation d'IsoFind porte automatiquement les valeurs de pH, Eh et DO utilisées.

Pour aller plus loin

Majeurs et traces : les concentrations totales qui alimentent le modèle de spéciation.
Ratios diagnostiques : compléments de diagnostic par combinaisons élémentaires.
Spéciation et propagation : usage de la spéciation dans le moteur de simulation 3D.
Cas panache Cr(VI) : exemple complet de simulation avec spéciation redox active.