Soutenance de thèse de ZHANG Yutong
Titre de thèse
Vers une modélisation parcimonieuse de la biogéochimie océanique basée sur des indicateurs clés pour des projections décennales robustes en Méditerranée
Towards parsimonious ocean biogeochemical modelling based on key indicators for robust decadal projections in the Mediterranean
Résumé de la thèse
Les modèles biogéochimiques (BGC) marins sont devenus des outils incontournables pour étudier le fonctionnement des écosystèmes marins, notamment son évolution future dans le cadre de projections climatiques. Les coûts de calcul associés, généralement élevés, et très dépendants de la complexité du modèle, peuvent cependant constituer un frein à l'exploration d'une diversité de scénarios, ce qui plaide pour des modèles simples et d'exécution rapide. A contrario, la confiance dans les prévisions dépend à la fois du réalisme des sorties du modèle, qui nécessite de conserver un certain niveau de complexité, et de la quantification des incertitudes associées.
C'est dans ce contexte que s'inscrit cette thèse qui s'est d'abord focalisée sur une simplification du modèle BGC marin Eco3M-MED ajustée à l'objectif visé, puis sur la quantification de ses incertitudes paramétriques et de celles des versions simplifiées obtenues. Une approche de simplification semi-automatique basée sur une analyse de sensibilité a été mise au point pour identifier les paramètres les moins influents et simplifier le modèle en conséquence, selon des critères objectifs. Cette analyse a bénéficié d'une nouvelle stratégie de paramétrisation distinguant les paramètres indépendants des paramètres dépendants, ce qui permet d'appliquer l'analyse de sensibilité à un modèle comportant centaines de paramètres avec un coût de calcul limité. Un autre avantage de l'approche SAIS est sa flexibilité: le niveau de simplification peut être progressivement augmenté en fonction de la question scientifique et des sorties du modèle que l'on souhaite conserver. Trois versions simplifiées, correspondant à des objectifs de moins en moins contraignants, ont été élaborées à l'aide de la configuration 1DV au site DYFAMED en Méditerranée Nord-Occidentale. Les résultats montrent que la simplification peut réduire le coût de calcul jusqu'à environ 30% tout en conservant une bonne fidélité des sorties ciblées, et que le modèle Eco3M-MED ne présente pas de sur-complexité évidente si conserver ses bases mécanistes et son adaptabilité aux différentes bio-régions de la Méditerranée reste une priorité. La version la plus simplifiée a également été implémentée en configuration 3D à l'échelle du bassin Méditerranéen. Les écarts à la simulation de référence restent assez modérés, y compris ailleurs qu'en Méditerranée nord-occidentale, et varient en fonction des grandeurs considérées.
Nous avons ensuite développé et automatisé un cadre pour quantifier les incertitudes paramétriques du modèle complet et étudier leur évolution avec la simplification du modèle. Malgré l'exploration de larges plages de valeurs plausibles des paramètres, le modèle reste globalement stable, avec seulement environ 2-3% des simulations pour lesquelles la structure initiale du modèle n'est plus conservée (disparition d'un ou plusieurs compartiments). Les résultats montrent surtout qu'il n'existe pas de relation directe entre la simplification du modèle et son niveau d'incertitude. La suppression de processus peut modifier la propagation des incertitudes et produire des effets de compensation, de sorte qu'une simplification ne conduit pas nécessairement à une diminution ou à une augmentation de l'incertitude. Cet effet dépend également des variables, des périodes et des profondeurs considérées, particulièrement pendant le bloom et la stratification, dans la couche près de la nutricline et du maximum profond de chlorophylle. Enfin, les indicateurs écosystémiques intégrés peuvent masquer certaines différences temporelles ou verticales d'incertitude des sorties du modèle, ce qui souligne l'importance de choisir les indicateurs en fonction de la question scientifique étudiée.
Thesis resume
Marine biogeochemical (BGC) models have become essential tools for studying the marine ecosystems functioning, particularly their future evolution in the context of climate projections. However, the associated computational costs, which are generally high and strongly dependent on model complexity, can limit the exploration of a wide range of scenarios, which argues in favor of simple and computationally efficient models. Conversely, confidence in projections depends both on the realism of model outputs, which requires maintaining a certain level of complexity, and on the quantification of the associated uncertainties.
It is within this context that this thesis was conducted. It first focused on simplifying the marine BGC model Eco3M-MED according to the targeted objective, and then on quantifying its parametric uncertainties and those of the obtained simplified versions. A semi-automated simplification approach based on sensitivity analysis was developed to identify the least influential parameters and simplify the model accordingly, using objective criteria. This analysis benefited from a new parameterization strategy that distinguishes independent from dependent parameters, making it possible to apply sensitivity analysis to a model containing hundreds of parameters at a limited computational cost. Another advantage of the SAIS approach is its flexibility: the level of simplification can be progressively increased depending on the scientific question and on the model outputs that need to be preserved. Three simplified versions, corresponding to progressively less restrictive objectives, were developed using the 1DV configuration at the DYFAMED site in the northwestern Mediterranean. The results show that simplification can reduce computational cost by up to about 30% while maintaining good fidelity for the targeted outputs, and that the Eco3M-MED model does not exhibit obvious over-complexity if preserving its mechanistic basis and its adaptability to the different Mediterranean bio-regions remains a priority. The most simplified version was also implemented in a 3D configuration at the Mediterranean basin scale. Differences relative to the reference simulation remain relatively moderate, including outside the northwestern Mediterranean, and vary depending on the quantities considered.
We then developed and automated a framework to quantify the parameter-induced uncertainties of the full model and to investigate how they evolve with model simplification. Despite exploring broad ranges of plausible parameter values, the model remains globally stable, with only about 2-3% of the simulations in which the initial model structure is no longer preserved (disappearance of one or more compartments). The results mainly show that there is no direct relationship between model simplification and its level of uncertainty. The removal of processes can modify the propagation of uncertainties and generate compensation effects, meaning that simplification does not necessarily lead to either a decrease or an increase in uncertainty. This effect also depends on the state variables, periods, and depths considered, particularly during the bloom and stratification periods, near the nutricline and the deep chlorophyll maximum. Finally, integrated ecosystem indicators can mask some spatio-temporal differences in the model outputs uncertainty, highlighting the importance of selecting indicators according to the scientific question being addressed.