vendredi, 24 avril, 2026
Isotopie CSIA
La CSIA (Compound-Specific Isotope Analysis) mesure la composition isotopique d'un composé individuel séparé d'une matrice complexe, au lieu d'un signal global bulk. Cette approche permet de suivre la transformation d'un polluant spécifique à travers sa dégradation, grâce au fractionnement isotopique qui enrichit le composé résiduel à mesure qu'il disparaît. Dans IsoFind, la CSIA est exposée au moteur de simulation et au module de rapports via le bridge molecule_csia_bridge qui interroge 56 fractionnements documentés sur 28 molécules dans 9 familles. Cette page présente les principes, les données tabulées et les modes d'accès programmatiques.
Principe du fractionnement isotopique
Lorsqu'une molécule polluante se dégrade, les liaisons chimiques ciblées par la réaction ont des cinétiques très légèrement différentes selon que l'atome impliqué est un isotope léger ou lourd. Pour un composé donné, la liaison C-¹²C se rompt un peu plus facilement que la liaison C-¹³C. Le résultat est que la fraction non encore dégradée s'enrichit progressivement en isotope lourd, tandis que le produit de dégradation part plus léger. Cet enrichissement se mesure et se quantifie, ce qui transforme une signature isotopique en preuve de dégradation.
δ(t) = δ₀ + ε · ln(f)
L'équation de Rayleigh décrit cette évolution : δ(t) est la signature isotopique du résiduel à un instant t, δ₀ la signature initiale de la source, ε le facteur d'enrichissement isotopique (en ‰), et f la fraction du composé encore présente par rapport à la masse initiale (f = C/C₀ entre 0 et 1). Cette équation est le cœur du diagnostic CSIA : quand on mesure δ(t) sur un échantillon et qu'on connaît δ₀ et ε, on peut calculer f et donc quantifier le taux de dégradation, même sans connaître l'âge exact du panache.
Tous les 56 fractionnements tabulés dans IsoFind sont en mode rayleigh. Les modes equilibrium et kinetic sont supportés par le bridge mais aucune entrée n'est actuellement en mode équilibre dans le catalogue moléculaire, ce qui est cohérent avec la réalité des dégradations terrain : les processus biotiques et abiotiques dominants sont quasi-exclusivement irréversibles et suivent un modèle de Rayleigh.
Le facteur epsilon et le AKIE
Le facteur d'enrichissement ε est la pente qui lie δ à ln(f). Il est exprimé en ‰ (pour mille) et se négocie entre 0 et −70 ‰ environ selon les cas. Une valeur proche de zéro signale un fractionnement faible (la dégradation ne distingue presque pas les isotopes) ; une valeur très négative signale au contraire un fractionnement marqué, donc un signal CSIA diagnostique puissant.
Le AKIE (Apparent Kinetic Isotope Effect) est une autre expression du même phénomène, utilisée dans la littérature mécaniste. Il s'exprime comme un ratio des constantes cinétiques des isotopes léger et lourd, normalement supérieur à 1. Pour un fractionnement intramoléculaire typique où seul un atome est impliqué dans la réaction, AKIE ≈ 1 − ε/1000. Les valeurs extrêmes du catalogue vont de 1,005 (TCA, fractionnement modéré) à 1,045 (DCM par Methylobacterium, le plus grand fractionnement biotique reporté dans la littérature).
Les quatre éléments supportés
Le bridge CSIA d'IsoFind supporte quatre éléments isotopiques principaux. La distribution actuelle des 56 fractionnements tabulés dans la base reflète l'état de l'art de la littérature.
| Élément | Système | Nombre d'entrées | Usage principal |
|---|---|---|---|
| Carbone | ¹³C / ¹²C (VPDB) | 36 | Élément universel, présent dans toutes les molécules organiques |
| Azote | ¹⁵N / ¹⁴N (AIR) | 9 | Pesticides triazines, médicaments, explosifs nitrés |
| Chlore | ³⁷Cl / ³⁵Cl (SMOC) | 9 | Solvants chlorés (diagnostique dual C/Cl) |
| Hydrogène | ²H / ¹H (VSMOW) | 1 | HAP (δ²H très diagnostique pour naphtalène) |
| Oxygène (cas particulier) | ¹⁸O / ¹⁶O | 1 | Perchlorate ClO₄⁻ |
Le bridge resolve_csia() accepte un paramètre preferred_element (défaut 'C') qui choisit l'élément privilégié pour la résolution. Si l'élément demandé n'est pas disponible pour la voie choisie, le bridge retombe automatiquement sur le carbone.
Répartition par famille chimique
Les 28 molécules avec CSIA dans le catalogue se répartissent en neuf familles, avec une concentration marquée sur les solvants chlorés et les pesticides.
| Famille | Molécules CSIA | Fractionnements |
|---|---|---|
| Solvants chlorés | 9 | 21 |
| Pesticides | 5 | 13 |
| Médicaments / PE | 3 | 6 |
| PFAS | 5 | 5 |
| HAP | 2 | 4 |
| Explosifs | 2 | 4 |
| Perchlorates | 1 | 2 |
| PCB | 1 | 1 |
Le dual isotope plot et le paramètre Lambda
Pour les molécules mesurées sur deux éléments isotopiques différents pour la même voie de dégradation, on peut construire un dual isotope plot : Δδ de l'élément secondaire en fonction de Δδ de l'élément primaire. La pente de ce plot, notée Λ (lambda), est un discriminant mécanistique bien plus puissant que le fractionnement sur un seul élément. Le bridge expose ce calcul via resolve_dual_isotope().
Λ = ε_secondaire / ε_primaire
Interprétation Λ pour C/Cl (solvants chlorés)
Le cas C/Cl est le mieux codifié, c'est la raison pour laquelle le bridge inclut des seuils interprétatifs explicites pour ce couple.
| Valeur absolue de Λ | Interprétation |
|---|---|
| Λ < 0,4 | Mécanisme biotique (Dehalococcoides, co-métabolisme) |
| 0,4 ≤ Λ < 0,7 | Voie biotique variable ou mixte |
| Λ ≥ 0,7 | Mécanisme abiotique (Fe-ZVI, oxydation chimique) |
Interprétation Λ pour C/N (pesticides et médicaments)
Pour les pesticides triazines, les médicaments et les explosifs, le couple pertinent est C/N avec des seuils différents.
| Valeur absolue de Λ | Interprétation |
|---|---|
| Λ ≤ 1,5 | Les deux atomes participent à la rupture |
| 1,5 < Λ ≤ 3 | Cible préférentielle azote |
| Λ > 3 | Azote périphérique (amine secondaire, impact minoritaire) |
Les couples dual isotope disponibles
Dix-huit couples de mesures sur deux éléments pour une même voie sont tabulés, ce qui permet le diagnostic dual isotope pour un large panel.
| Famille | Molécule | Couple |
|---|---|---|
| Solvants chlorés | PCE, TCE, cDCE, DCM, Chloroforme, 1,2-DCA | C / Cl |
| Pesticides | Atrazine (2 voies), Diuron, Glyphosate | C / N |
| Pesticides | S-Métolachlore | C / Cl |
| Médicaments / PE | Carbamazépine, Diclofénac, Sulfaméthoxazole | C / N |
| Explosifs | RDX, TNT | C / N |
| HAP | Naphtalène | C / H |
| Perchlorates | Perchlorate (ClO₄⁻) | O / Cl |
Sélection automatique de la voie de dégradation
Une molécule peut avoir plusieurs voies de dégradation tabulées pour des conditions différentes (aérobie vs anaérobie, biotique vs abiotique). Le bridge CSIA utilise une règle de sélection explicite à trois niveaux pour choisir la voie pertinente.
| Priorité | Critère | Comportement |
|---|---|---|
| 1 | Nom de voie explicite (pathway_name) | La voie demandée est utilisée directement, avec matching partiel en fallback |
| 2 | Conditions redox (Eh, O₂) | Si Eh < 50 mV ou O₂ < 0,5 mg/L : voie anaérobie. Sinon : aérobie stricte |
| 3 | Défaut sub-surface | Voie naturelle la plus rapide (demi-vie la plus courte), hors photolytique et PRB |
Le contexte par défaut est aquifère sub-surface, ce qui exclut les voies photolytiques (pas d'UV en profondeur) et les voies de barrière réactive artificielle (Fe-ZVI). Ces voies restent accessibles en les demandant explicitement via pathway_name ou en désactivant le flag sub-surface.
La structure de la réponse resolve_csia
L'appel resolve_csia(molecule_name, pathway_name, preferred_element, conditions) retourne un dictionnaire riche qui alimente à la fois le moteur de simulation et les blocs de rapport.
| Champ | Type | Contenu |
|---|---|---|
| molecule, molecule_id | str, int | Molécule résolue (tolérant aux abréviations) |
| famille, masse_molaire | str, float | Contexte taxinomique |
| pathway, pathway_category, environment | str | Voie choisie et ses conditions |
| element | str | Élément retenu (C, N, Cl, H, O) |
| epsilon, epsilon_min, epsilon_max, uncertainty | float | Fractionnement et son incertitude |
| mode | str | rayleigh | equilibrium | kinetic |
| akie, lambda | float | AKIE et couplage multi-élément |
| k_deg, half_life_days | float | Constante cinétique (j⁻¹) et demi-vie (jours) |
| metabolites | list | Liste [{name, formula, mass, molar_yield_max}] |
| reference, n_records | str, int | Référence bibliographique et nombre d'entrées trouvées |
Chaîne d'utilisation dans IsoFind
La CSIA n'est pas un module isolé : elle alimente activement plusieurs autres composants du logiciel. La même résolution sert à plusieurs destinations.
| Destination | Usage du résultat CSIA |
|---|---|
| Moteur de simulation 3D | Remplace l'epsilon par défaut, fournit k_deg pour la cinétique, déclenche la cascade des métabolites |
| Interface échantillon (onglet Molécules) | Suggestions de mesures isotopiques à réaliser pour un échantillon donné |
| Module Nexus Quick Match | Diagnostic mécanistique (biotique vs abiotique) à partir de signatures terrain |
| Blocs de rapport | Blocs csia_diagnosis, dual_isotope_plot, degradation_fingerprint |
Accès via l'API
Quatre endpoints exposent directement le bridge CSIA aux applications clientes.
| Endpoint | Usage |
|---|---|
| GET /api/molecules/csia/list | Liste des molécules avec CSIA disponibles (28 molécules, avec compteurs) |
| GET /api/molecules/csia/{molecule_name}/pathways | Voies de dégradation tabulées pour une molécule |
| POST /api/molecules/csia/resolve | Résolution CSIA complète avec conditions optionnelles (Eh, O₂, pH, T) |
| POST /api/molecules/csia/dual | Dual isotope avec calcul de Λ et interprétation codée C/Cl ou C/N |
Limites et bonnes pratiques
- La CSIA demande des concentrations suffisantes pour l'analyse : en dessous de la LOQ CSIA (généralement supérieure à la LOQ classique), l'incertitude sur δ devient prohibitive. Il est fréquent qu'un échantillon conforme à un seuil réglementaire soit en dessous de la LOQ CSIA.
- L'epsilon tabulé est une valeur représentative issue d'études en conditions contrôlées. La variabilité inter-souches peut être importante ; le catalogue IsoFind préserve les plages min-max pour encadrer cette incertitude.
- Une valeur δ enrichie seule ne prouve pas une voie de dégradation spécifique. Seul le dual isotope plot avec la valeur de Λ permet une attribution mécanistique robuste.
- La signature δ₀ de la source est rarement mesurable directement. Elle est souvent estimée à partir de fiches matières premières, de témoins non dégradés sur site ou de la valeur la moins enrichie du panache.
- Les modèles de Rayleigh supposent un système fermé. Sur un panache ouvert avec apports continus, l'interprétation doit intégrer le mélange, ce que le moteur de simulation IsoFind fait automatiquement via son prior géochimique.
Pour aller plus loin
- Solvants chlorés : la famille de référence pour la CSIA, 21 fractionnements tabulés.
- Pesticides : CSIA sur C/N pour les triazines et les phosphonates.
- Voies de dégradation : les 50 voies tabulées qui sous-tendent les fractionnements.
- Spéciation et propagation : comment le moteur de simulation utilise le résultat CSIA.