Télécharger le fichier original (Mémoire de fin d’études)
Le Système Global d’Observation des Océans
L’océan, difficilement observable par l’homme, est resté longtemps méconnu. Ce n’est qu’à partir du XIX me siècle que les premières campagnes d’explorations maritimes autant europèennes qu’américaines débutent notamment avec la cam-pagne océanographique anglaise du Challenger (1872 − 1876). D’ambitieux pro-grammes internationaux d’observation de l’océan ont été menés dans les années 80 et 90 (TOGA, WOCE) et ont conduit à des avancées majeures dans la compréhen-sion de l’océan et de son rôle sur le climat. La mise en place depuis une dizaine d’années des moyens nouveaux d’observations temps réel (satellites, instruments autonomes in-situ dont le réseau global de flotteurs profilants Argo) ainsi que les progrès sur la modélisation numérique ont permis plus récemment le développe-ment de l’océanographie opérationnelle et des services océaniques associés.
Un système global d’observation des océans performant est indispensable pour observer, comprendre et prévoir l’océan et son rôle sur le climat. Le Système Global d’Observation des Océans actuel (GOOS, http ://www.ioc-goos.org/, Dexter and Summerhayes [2010]) résulte de la juxtaposition et de l’intégration progressive de systèmes de mesures variés et complémentaires (mesures in-situ et satellites, voir la figure 1.1). Ces systèmes d’observation sont indispensables aux centres d’anal-yses et de prévisions océaniques. Le centre français Mercator Ocean assimile ainsi en temps réel l’ensemble de ces données dans des modèles océaniques globaux et régionaux.
Cette thèse s’inscrit dans les activités de l’Ifremer 1 en collaboration avec Mer-cator Océan autour du programme Argo et ses évolutions. Dès le début des années 2000, le programme Argo s’amorçe. L’année 2007 le réseau contient plus de 3000 flotteurs, fournissant des mesures globales de profils de température et de salinité sur les 2000 premiers mètres de l’océan. Le déploiement d’instruments sous-marins autonomes et programmés à dériver sur toutes les mers du globe est sans précé-dent dans l’histoire des mesures in situ. Les observations Argo possèdent, en outre une forte complémentarité avec les observations satellitaires. C’est une des raisons pour laquelle le réseau Argo est une composante essentielle du système mondial d’observation des océans et du climat permettant de :
– détecter la variabilité climatique aux échelles saisonnières et décennales, et observer à long terme le changement climatique dans les océans ;
– fournir des données essentielles pour contraindre des modèles d’analyse et de prévision océanique (Service Marin de Copernicus 2 , Mercator Océan), afin d’initialiser des modèles couplés océan-atmosphère de prévision saisonnière et décennale et de valider les modèles climatiques ;
– fournir l’information nécessaire à l’étalonnage des capteurs et la validation des données satellitaires.
Apport des données Argo
L’observation opérationnelle de l’océan est ainsi devenue une réalité. La ca-pacité d’observer l’océan en temps réel via les mesures des satellites et via les réseaux in-situ comme Argo ont permis la mise en place de modèles d’analyse et de prévision de tout l’océan analogues aux modèles métérologiques. Ainsi le centre Mercator Océan fournit maintenant en routine des analyses et prévisions océaniques sur l’ensemble des océans. Des cartes de prévisions à 5, 7 ou 14 jours des courants, température, salinité sont actualisées de façon quotidienne ou hebdomadaire. A ce jour le réseau Argo représente une des plus importantes sources de données in-situ en assurant une couverture globale spatiale et temporelle. Le nombre d’études utilisant les données Argo grandit et des résultats notables ont été obtenus 3 . Les données Argo sont également utilisées de manière systématique au sein des systèmes d’océanographie opérationnelle. Il est donc nécessaire de con-tinuer acquérir une bonne connaissance de “la donnée Argo” et de ses limitations pour une utilisation optimale et une pérennisation du réseau.
Les flotteurs Argo ont permis de combler les manques significatifs présents dans le réseau global d’échantillonnage de l’océan (Roemmich and the Argo Steer-ing Team [2009]). Ainsi, ils fournissent des informations très utiles quant à l’esti-mation de l’état en surface des océans et de la structure 3D de la température et de la salinité. Les cartes globales du niveau de la mer obtenues par altimétrie satel-litaire, à des échelles hebdomadaires comme annuelles, peuvent maintenant être interprétées grâce à une meilleure connaissance complète de la structure tridimen-sionnelle des champs de température et de salinité. Dès lors les modèles océaniques globaux et climatiques peuvent être plus précisément initialisés, testés, et con-traints. (Guinehut et al. [2012]. De plus Juza et al. [2011b,a] ont montré, en par-ticulier, que le réseau Argo est un système d’observations in-situ bien adapté à l’étude d’indices climatiques. Notamment il capture une grande partie de la vari-abilité saisonnière et interannuelle du contenu thermique de l’océan global.
Depuis le déploiement du réseau Argo la communauté océanographique a dé-montré les apports essentiels des données Argo dans les systèmes d’assimilation de données et de modélisation numérique, ainsi que leur cohérence avec les données al-timétriques et les autres types de données in-situ (documentés dans Roemmich and the Argo Steering Team [2009]). Les études conduites par la communauté au sein de GODAE 4 détaillées dans Oke et al. [2009] exposent les différentes approches utilisées pour quantifier l’impact de ces observations dans des systèmes avec assim-ilation de données. Il est usuel de quantifier l’impact des observations au sein d’un système de prévisions/analyse sous la forme d’OSEs 5 (Oke and Schiller [2007], Smith and Haines [2009] et Oke et al. [2014]). Mercator Océan met régulièrement en œuvre ces techniques pour quantifier l’apport des données Argo. Les récentes études menées (Turpin et al. [2015]) mettent en évidence l’apport très important des données Argo au niveau global et régional. Une réduction de la couverture du réseau Argo dégrade systématiquement et de façon marquée la qualité des analyses et prévisions océaniques. D’autre part les résultats d’une étude réalisée au CEP-MMT 6 , basée sur la même approche que celle décrite pour le système de Mercator Océan montrent l’importance des observations Argo pour initialiser des modèles couplés océan/atmosphère pour la prévision saisonnière (Oke et al. [2014]). Enfin l’OSE utilisant le système MOVE-G d’assimilation de données japonais souligne que l’impact des données Argo est plus important en salinité qu’en température. Ce résultat est cohérent avec les conclusions des études passées (Oke and Schiller [2007]).
Objectifs de la thèse
Les systèmes d’observations océaniques fournissent aux océanographes un grand nombre d’estimations de diverses quantités géophysiques (température, salinité, niveau de surface de la mer, vitesses des courants . . . ) à des résolutions spatiales et temporelles hétéroclites. De telles observations peuvent être utilisées pour créer des références communes sous forme d’atlas climatologiques construits sur de longues périodes (World Atlas Ocean documenté dans Levitus et al. [1998, 2005, 2013]). La description de l’état complet de l’océan et la prévision de son évolution nécessitent cependant d’intégrer ces observations dans des systèmes d’assimilation. On com-bine alors les observations avec un modèle océanique. Par conséquent caractériser les erreurs du modèle est alors primordial. Celles-ci peuvent être évaluées par com-paraison entre la quantité simulée et la quantité observée. On attend de la part des modèles océaniques qu’ils parviennent à représenter certaines caractéristiques décisives de l’état et de la variabilité physique de l’océan. Grâce aux données Argo il est maintenant possible d’évaluer ces caractéristiques à différentes profondeurs en température et en salinité. Dans cette thèse nous confrontons les observations Argo et une simulation numérique libre dans le but d’évaluer la capacité d’un mod-èle de circulation générale au 1/4° dont elle est issue, à reproduire correctement les caractéristiques principales de la circulation océanique.
Les études Kamenkovich et al. [2009, 2011], Vecchi and Harrison [2005], Guine-hut et al. [2012] et Oke et al. [2014] s’intéressent aux capacités de reconstructions de divers champs océaniques à l’aide des observations Argo. Certains facteurs tels que les mouvements des flotteurs, leurs cycles d’échantillonnage, les tourbillons à méso-échelles peuvent influencer la qualité et la précision des résultats obtenus. Analyser de tels facteurs est essentiel à la compréhension des limitations du sys-tème d’observations Argo. Certaines étapes préliminaires aux estimations statis-tiques réalisées dans cette thèse ont été effectuées dans le but d’évaluer l’impact de facteurs tels que le nombre et les mouvements des flotteurs selon les régions observées. Nos choix en ce qui concerne l’approche d’estimation en tiendront compte.
Plusieurs études témoignent des nouvelles analyses scientifiques rendues pos-sibles grâce au réseau global Argo. Cependant ce potentiel implique un besoin de validation et de contrôle de qualité au sein de ce large jeu de données. En effet, Gaillard et al. [2009] et Roemmich and Gilson [2010] soulignent que con-naître les variations des échelles spatiales en température et salinité est nécessaire pour améliorer les contrôles de qualité mais aussi pour caractériser au mieux la dynamique océanique et spécifier les schémas d’assimilation de données. L’un des objectifs de cette thèse concerne la caractérisation des échelles spatiales et tem-porelles en température et salinité de la surface à 1500 m, et il a été motivé par plusieurs faits :
– L’assimilation de données permet en pratique de déterminer la meilleure com-binaison entre les observations et un des états à priori de l’océan représenté par un modèle numérique. Une part importante en assimilation de données consiste à déterminer dans quelle mesure le modèle est réaliste. Il s’agit de lui affecter un poids approprié dans le but d’approcher la réalité et d’améliorer les capacités d’analyse et de prévision du modèle. Ce poids doit prendre en compte les échelles spatio-temporelles des erreurs modèles (qui sont souvent approximées par les échelles du signal océanique lui même).
– Le système d’assimilation de Mercator Ocean basé sur la formulation SEEK utilise des informations sur les échelles de corrélation. En effet une carac-téristique importante au sein de cette formulation est la spécification de la matrice de covariance d’erreurs de prévision du modèle. Sa définition repose sur une technique de localisation impliquant l’utilisation d’échelles de cor-rélation spatiales. Parce que la matrice de covariance d’erreurs de prévision du modèle est construite à partir d’un nombre fini d’anomalies de l’état de l’océan, on choisit de mettre à zéro les covariances à partir d’une certaine distance définie comme étant deux fois l’échelle de corrélation spatiale. Ce choix est légitime du fait que l’ensemble fini des anomalies n’est plus signi-ficatif à partir d’une certaine distance du point d’analyse. Ces échelles ont été estimées jusqu’à présent en surface via les anomalies des hauteurs d’eau observées par satellites (SLA, SSALTO/DUACS), et à faible profondeur (100 et 300 m) à partir de champs 2D en température issus de la ré-analyse glob-ale au 1/4◦ GLORYS2V1 (GLobal Ocean ReanalYsis ans Simulation, Ferry et al. [2012]; Lellouche et al. [2013]). Cependant Kessler et al. [1996], Jacobs et al. [2001] ont souligné des dépendances en fonction des latitudes et de la profondeur au sein des échelles spatiales. Il s’agit donc de poursuivre au mieux la caractérisation des échelles spatiales en température et salinité de la surface à 1500 mètres, ainsi que leur variation en fonction de la profondeur, et des latitudes et longitudes.
– Par ailleurs ces échelles permettent aussi de définir une “bulle d’influence” au voisinage du point d’analyse et d’y sélectionner les données utilisées (Lel-louche et al. [2013]). Les observations utilisées dans le cadre de l’assimila-tion de données possèdent une contribution différentes selon leurs nature et qualité (Turpin et al. [2015]). Il est donc intéressant de savoir mesurer objectivement la quantité d’information apportée par les données Argo par comparaison aux observations de natures différentes composants les jeux de données assimilés.
On dispose également d’outils permettant d’évaluer la contribution relative des observations dans n’importe quel système d’assimilation de données : l’un d’eux est le DFS pour Degrees of Freedom. La notion de DFS permet d’estimer la contribution relative des observations mais il est aussi intéressant de comprendre le lien entre les échelles de corrélation spatiales et la contribution relative des observations. Cette dépendance est illustrée par la suite avec un exemple simplifié pour ensuite tester différentes méthodes de calcul approximatif du DFS applicables facilement à des systèmes d’assimilation tel que celui opéré à Mercator Océan.
Structure du manuscrit
Dans cette thèse j’ai tout d’abord exploité les observations hydrographiques issues du réseau Argo, pour caractériser les erreurs d’un modèle global au 1/4° n’ayant pas assimilé de données Argo (chapitre 3). Dans ce contexte, il a fallu utiliser des méthodes d’interpolation afin d’extraire de la simulation numérique choisie pour l’étude, le pendant exact (en temps et en espace) des observations Argo. Cette étape correspond à la colocalisation au point d’observation des sorties modèles en température et salinité. Dans un deuxième temps j’ai effectué un traite-ment statistique en temps et en espace dans le but d’obtenir les erreurs modèles moyennes ainsi que les variances d’erreurs en température et salinité en surface et en profondeur. L’analyse de ces résultats permet d’évaluer le comportement de la simulation numérique au regard des observations Argo.
La variabilité océanique possède une large gamme d’échelles spatiales et tem-porelles. Afin d’étudier plus finement la manière dont cette simulation numérique reconstruit les champs océaniques de température et de salinité, nous avons car-actérisé les échelles spatiales des erreurs modèles et du signal océanique lui même (chapitre 4). J’ai développé une méthode pour l’estimation de ces échelles spatiales à partir des observations Argo. Une analyse de l’impact de l’échantillonnage du réseau Argo sur l’estimation de ces échelles grâce à des simulations et un calcul d’erreur formelle ont été menés. Les échelles spatiales de corrélation en température et en salinité de la surface à 1500 m ont été, enfin, calculées à partir d’anomalies par rapport à la climatologie et par rapport aux sorties modèles sur une période de plusieurs années. Les résultats ont été analysés afin de mettre en évidence les vari-ations géographiques de ces échelles, les variations en profondeur et les variations entre température et salinité.
Dans un troisième temps, j’ai analysé les méthodes basées sur le contenu en information. Ces méthodes permettent de mesurer directement l’impact ou la con-tribution relative des observations dans un système d’assimilation comme celui mis en œuvre à Mercator Ocean. J’ai mené une étude de sensibilité (bruit de mesure, longueur de corrélation) et testé différentes méthodes permettant d’approximer le contenu en information dans un contexte d’analyse objective. Ces tests ont permis de déterminer la formulation qui conviendrait le mieux à une implémentation dans le système d’assimilation de données de Mercator Ocean.
Table des matières
1 Introduction générale
1.1 Le Système Global d’Observation des Océans
1.2 Apport des données Argo
1.3 Objectifs de la thèse
1.4 Structure du manuscrit
2 Le programme international Argo et le centre d’océanographie opérationnelle Mercator Ocean
2.1 Le programme international Argo
2.1.1 Les objectifs
2.1.2 Cycle de plongée et flux de données Argo
2.1.3 Résultats prometteurs du réseau Argo
2.1.4 Participation très active de la France
2.1.5 Les nouveaux défis pour Argo
2.2 Mercator Ocean, Centre francais d’analyse et de prévision océanique
2.2.1 Histoire et projets européens
2.2.2 Modèle océanique et systèmes d’assimilation
3 Caractérisation des écarts modèle observations pour une simulation sans assimilation
3.1 Méthodes statistiques
3.2 Caractéristiques propres du modèle
3.3 Présentation des observations in-situ Argo
3.4 Evaluation des écarts entre la simulation et les observations in-situ Argo
3.4.1 Séries temporelles des écarts moyens
3.4.2 Moyennes et variances en température et salinité sur la période 2006-2009
3.4.3 Moyennes saisonnières en température et salinité sur la période 2006-2009
4 Caractérisation des échelles spatiales des erreurs modèles et de la variabilité océanique pour la température et la salinité
4.1 Introduction
4.2 Méthodologie pour l’estimation des échelles spatiales
4.2.1 Discrétisation des covariances
4.2.2 Ajustement d’une fonction analytique de corrélation
4.2.3 Validation de la méthodologie à partir de données simulées
4.2.4 Conclusion
4.3 Caractérisation des échelles spatiales des erreurs modèles et du signal à partir des données Argo
4.3.1 Structures et caractéristiques des régions examinées
4.3.2 Utilisation d’une nouvelle formulation analytique du modèle de corrélation
4.3.3 Analyse globale des échelles spatiales des erreurs modèles en température et salinité
4.3.4 Comparaison des échelles spatiales des erreurs modèles par rapport aux échelles de la variabilité océanique
4.4 Estimation des échelles spatiales des variations de température et salinité à partir des données Argo
4.4.1 Article paru dans le journal Ocean Science
4.4.2 Erreur formelle sur l’estimation de la covariance
4.5 Conclusion
5 Analyse des méthodes basées sur le contenu en information et étude de sensibilité
5.1 Introduction
5.2 Mesure d’impact des observations sur l’analyse
5.2.1 Définition analytique
5.2.2 Définition dans le cadre du BLUE
5.2.3 Definition dans le cadre de l’analyse objective
5.2.4 Méthodes alternatives pour l’approximation du DFS
5.3 Etude de sensibilité du contenu en information
5.3.1 Simulation de distributions gaussiennes et estimateur de Gauss- Markov
5.3.2 Tests d’influence sur le DFS
5.4 Comparaison des estimations du contenu en information
5.4.1 Illustration à partir d’une application au système d’assimilation de Mercator Océan
5.4.2 Mise en pratique dans le cadre d’analyse objective
5.5 Conclusion
Conclusion et perspectives
Annexes
A Catalogue des biais annuels et saisonniers modèle-observations sur la période 2006-2009
B Catalogue des échelles spatiales de corr