Amélioration de la modélisation des effets écotoxicologiques terrestres des métaux en évaluation des impacts du cycle de vie

L'analyse du cycle de vie (ACV) est un outil d'aide à la décision qui permet de dresser l'empreinte environnementale d'un produit ou d'un service tout le long de son cycle de vie. Les substances émises ou les ressources extraites du produit analysé sont traduites en impacts pour différentes catégories (réchauffement climatique, écotoxicologie terrestre et aquatique) via un facteur de caractérisation (CF). Pour les métaux, le CF de la catégorie écotoxicité est calculé comme le produit de trois facteurs : un facteur de devenir, (FF), un facteur de biodisponibilité (FB) et un facteur d'effet (FE). Le FF évalue la fraction de la substance qui est transportée du compartiment environnemental d'émission au compartiment de réception et son temps de résidence. Le BF décrit la fraction de substance biodisponible et le FE modélise l'effet sur la qualité et la biodiversité des écosystèmes, c'est à dire la fraction d'espèces potentiellement affectées, par concentration d'exposition, qui est calculé à partir des données de tests écotoxicologiques. Dans le cas des organismes terrestres, les FEs ne tiennent actuellement pas compte des propriétés des métaux. En effet, le FE d'un métal n'est pas habituellement calculé à partir des données d'effet écotoxicologiques terrestres, qui sont rares et ne couvrent pas le nombre de données recommandées par la méthode AMI (Assessment of the Mean Impact) utilisée dans USEtox. Le FE est extrapolé à partir des données des organismes aquatiques. Cette extrapolation, appelée l'équation de partition à l'équilibre (EqP) est basée sur le coefficient de partition sol-eau Kd et la densité du sol. Un unique Kd est présentement utilisé pour représenter la répartition du métal entre les phases solides et aqueuses du sol, quel que soit le sol dans lequel le métal est émis. Or dans le cas des substances métalliques, la toxicité d'un métal dépend de sa spéciation, et donc de sa répartition dans le sol. La spéciation d'un métal est fonction des propriétés physicochimiques des sols (pH, taux de matière organique, texture, etc.). Ainsi, utiliser une seule valeur de Kd quel que soit le sol de réception signifie que l'EqP considère que les propriétés des sols ne varient pas spatialement et donc que la toxicité des métaux ne varie pas selon ces propriétés. De plus, l'EqP suppose que la sensibilité des écosystèmes aquatiques et terrestres est équivalente. Malgré les approximations de cette équation, il est intéressant d'utiliser les données écotoxiques aquatiques, qui sont bien plus nombreuses que les données terrestres, et de les extrapoler pour obtenir des données d'effet sur les organismes terrestres. Néanmoins, il est nécessaire de tenir compte dans l'EqP de la variation spatiale des propriétés des sols, et donc de la spéciation des métaux. La plupart des travaux ayant pour but d'intégrer la spéciation des métaux dans le calcul du CF portent sur l'amélioration du FF, et peu sur le FE (Plouffe, Bulle, & Deschênes, 2015d). Owsianiak et al (2013) ont inclus la spéciation du métal dans le calcul du FE, en utilisant la méthode du Terrestrial Biotic Ligand Model (TBLM) qui considère que seule la fraction du métal étant lié aux ligands biotique (site récepteur des organismes) aura un effet toxique. Il tient compte de la compétition entre les cations présents dans les sols et les cations métallique pour ces sites récepteurs. Ce modèle est intéressant pour intégrer la spéciation du métal dans le calcul du FE, mais il est pour l'instant uniquement applicable au cuivre et au nickel, en suivant les recommandations de méthode AMI. Dans les méthodes actuelles d'évaluation des impacts, des FEs non-régionalisés sont actuellement utilisés pour les métaux. Ces limites nous amènent à définir les objectifs de recherche suivants : Déterminer si la forme du métal, et donc si sa spéciation, a une influence significative (de plus d'un ordre de grandeur) sur la valeur du FE. Inclure la spéciation du métal, pour une grande variété de sols, dans l'EqP et vérifier si ces nouveaux FE régionalisés peuvent être utilisés pour quantifier l'effet des métaux en les comparant aux autres approches de la littérature, à partir des données écotoxicologiques terrestres et le TBLM. Estimer la contribution relative des FEs, FBs et FFs régionalisés à la variabilité spatiale des FCs Afin de réaliser ces objectifs, des FEs ont été calculés en utilisant l'EqP régionalisée via des Kd qui tiennent compte de la spéciation des métaux pour le cuivre, le nickel et le zinc. Ces métaux ont été choisis car le cuivre et le zinc sont les principaux contributeurs aux impacts écotoxicologiques terrestres en ACV au Canada, USA et en Europe, tandis que le nickel et le zinc sont les deux métaux dont la contribution aux impacts varient significativement selon la méthode d'EICV appliquée. Puis ces FEs ont été comparés aux FEs utilisés dans les méthodes d'EICV et développés dans la littérature. Dans chaque cas, la méthode AMI a été utilisée, pour laquelle FE=0,5/〖HC50〗_EC50 où la HC50EC50 est la concentration de risque affectant 50% des espèces présentes dans l'écosystème à leur EC50. La EC50 est la concentration ayant un effet sur 50% de la population d'une espèce particulière. Ainsi les HC50EC50 ont été déterminée de quatre façons différentes : a- HC50EC50 terrestres, extrapolées à partir des HC50EC50 aquatiques via l'EqP régionalisée. Des Kd régionalisés ont été utilisés pour obtenir des FEs pour le métal total (FEEqPregionalise). Les Kd régionalisés sont déterminés avec WHAM 7, un modèle de spéciation géochimique pour 5213 sols sélectionnés à partir de la HSWD (Harmonized World Soil Database). Aucune étude prouve actuellement que le modèle de spéciation géochimique tient compte des mécanismes de précipitations, l'étude s'est donc principalement portée sur 2908 sols non-calcaires. Des FEs sont obtenus pour les sols calcaires afin de servir de valeur intérimaires mais doivent être utilisées avec précaution. ; b- HC50EC50 calculées à partir des EC50 provenant de données écotoxicologiques terrestres – lorsque celles-ci sont disponibles – et séparées en fonction de la forme du métal (total, soluble, labile) (FEterrnon-regionalise). c- HC50EC50 terrestres, extrapolées à partir des HC50EC50 aquatiques via l'EqP non-régionalisée, telle que utilisée habituellement dans IMPACT 2002+, pour déterminer un FEs pour le métal total (FEEqPnon-regionalise ) d- HC50EC50 calculés à partir des EC50 terrestres pour la fraction ionique du métal obtenus via le modèle de ligand biotique terrestre (TBLM), uniquement pour la cuivre et le nickel et la modélisation de la spéciation d'après WHAM 7. Des FEs régionalisés sont calculés pour les 5213 sols utilisés en (a) (FETBLMregionalise) Les FEs déterminés à partir de l'EqP régionalisées (a) ont été comparés aux FEs provenant de ces différentes sources de données écotoxicologiques (b, c et d). La corrélation des FEs régionalisés calculés à partir de l'EqP et du TBLM a également été étudiée. Afin de déterminer le facteur contribuant le plus à la variabilité spatiale du FC, de nouveaux FCs ont été calculés pour le zinc, en utilisant les FEEqPregionalise et les FF et BF régionalisés développés par Plouffe (2014). De plus, une analyse de l'influence des propriétés physicochimiques du sol sur la valeur des FEs a été réalisée, ainsi que le calcul des FCs pour le zinc pour deux sols différentiés, avec un FE régionalisé et non-régionalisé. Ces différentes méthodes ont abouti à l'obtention de quatre types de FEs, dont deux sont régionalisés, pour l'écotoxicologie terrestre. Des variations 2 à 3 ordres de grandeurs en fonction de la forme du métal des FEterrnon-regionalise indiquent que la spéciation du métal influe la valeur des FEs et donc l'estimation de l'effet potentiel. Les FEs régionalisés selon l'approche EqP (FEEqPregionalise) varient de plusieurs ordres de grandeurs en fonction du type de sols (entre 6 et 13 ordres de grandeur) ce qui montre le nécessité d'inclure un Kd régionalisé lors de l'application de l'EqP. De façon plus globale, l'importante variabilité des EFs régionalisés, quelle que soit la méthode, indique que la spéciation du métal, et donc les propriétés des sols, contribuent fortement à la variation spatiale des FEs. Comparés aux FEs régionalisés calculés à partir de la méthode TBLM (FETBLMregionalise) , les FEsEqPregionalise couvrent une plus faible plage de valeurs. La médiane de ces deux méthodes est proche de 2 ordres de grandeurs maximum de la valeur du FE calculé à partir des données écotoxicologiques terrestres pour le métal total, ce qui indique que les valeurs de FE calculées sont cohérentes. Néanmoins, ces deux méthodes ne sont pas corrélées. La différence entre les FEsEqPregionalise et les FETBLMregionalise peut provenir de la fraction biodisponible considérée (fraction soluble pour l'EqP et la fraction d'ion métallique libre pour le TBLM) et de la différence de sensibilité entre les organismes aquatiques et terrestres, considérés égale dans l'EqP. Les valeurs des CFs régionalisés obtenues pour le zinc à partir des FEs, FBs et FFs spécifiques aux sols dépendent fortement de la valeur du FEEqPregionalise, celui-ci étant le facteur le plus corrélé. Ce mémoire a permis d'obtenir des FEs spécifiques aux paramètres physicochimiques des sols à partir des l'EqP et de données écotoxiques aquatiques. Ainsi l'effet potentiel du cuivre, du nickel et du zinc sur les organismes terrestres a été modélisé dans de nombreux sols aux propriétés différentes. De plus, l'influence de la forme du métal sur la valeur du FE a été quantifiée ainsi que la contribution des FEs à la variabilité spatiale des FCs. Il a été possible de conclure que tenir compte de la spéciation dans le calcul des FEs du zinc, du nickel et du cuivre en écotoxicité terrestre, permet de mieux différencier les effets potentiels spatialement et dans le cas du zinc, les impacts. La faible corrélation entre l'approche EqP régionalisée et le TBLM a monté qu'extrapoler les données aquatiques selon la méthode EqP, même en intégrant la spéciation du métal, n'estime pas de façon robuste la toxicité du métal par rapport à la modélisation du TBLM.