Les chondrites à enstatite partagent un réservoir commun de l'isotope de l'oxygène avec la Terre et la Lune. Cette caractéristique soulève la question d'une possible participation de ce type de météorites à la construction de la Terre. Egalement, l'étude des chondrites à enstatite permet d'apporter de nouvelles contraintes pour comprendre plus en détail l'histoire de la différenciation de petits corps planétaires soumis à des conditions réductrices. La première étude expérimentale systématique à haute pression et haute température d'une chondrite à enstatite a ainsi été mise en oeuvre, entre 1 et 25 GPa, et entre 1200°C et 2500°C, couvrant ainsi les conditions de pression et de température du manteau supérieur terrestre. La météorite qui a été choisie pour cette étude s'appelle Indarch. Il s'agit d'une chondrite à enstatite de type EH4. A 1 GPa, l'effet de la fO2 sur les relations de phases de la météorite, ainsi que sur les coefficients de partage entre métal et silicate liquide des éléments S, Si, Cr, Mn, Ni, et Mo est étudié spécifiquement. Les phases silicatées et métalliques subissent de grands changements lorsque la fO2 varie de IW-1.5 à IW-4.5. Des monosulfures contenant les éléments (Fe, Mg, Mn, Ca, Cr) sont présents aux fO2 les plus réduites. La phase métallique de nos échantillons est toujours liquide et comprend deux phases immiscibles riches en fer : l'une pauvre en S et riche en C, et l'autre riche en S et pauvre en C. La phase riche en C est également riche en Si aux fO2 les plus réduites. Les relations de phases d'Indarch ont été déterminées entre 3 et 25 GPa, et entre 1500°C et 2500°C, nous permettant de proposer le premier diagramme de phases en pression et température pour une chondrite à enstatite. La phase métallique est toujours liquide et il s'agit d'un alliage Fe-Ni-S à ces plus hautes pressions. La formation d'un noyau planétaire à basse pression (0 à 5 GPa) sous des conditions très réductrices ne peut expliquer les appauvrissements observés de S dans les manteaux silicatés alors que plusieurs wt% de Si pourront être piégés dans le noyau. En revanche, à plus haute pression, l'incorporation de S dans le noyau sera facilitée sous des conditions plus oxydantes. Finalement, notre étude permet de proposer un modèle d'évolution d'un matériau chondritique en cours d'accrétion sous des conditions rédox variables

Les océans Pacifique et Indien tropicaux se partagent la plus grande étendue d'eau chaude et de convection profonde de la planète. Cette région est le siège de la branche ascendante de la circulation de Walker, circulation atmosphérique d'échelle planétaire parfois décrite comme la " machine thermique " de la Terre. Cette région, dont les répercussions sur le climat sont importantes, est aussi source de variabilité océanique et atmosphérique aux échelles intrasaisonnières et interannuelles . En effet, la variabilité interannuelle associée à El Niño, dans l'océan Pacifique, et - dans une moindre mesure - au dipôle de l'Océan Indien (DOI) ont des conséquences climatiques marquées sur les pourtours de ces bassins et à l'échelle du globe. La variabilité intrasaisonnière liée à l'oscillation de Madden-Julian (OMJ) a également des conséquences climatiques marquées : modulations des moussons indiennes et australiennes et un rôle potentiellement important dans le déclenchement d'ENSO. Dans ce mémoire, je vais décrire mes travaux de recherche sur la variabilité océanique et atmosphérique aux échelles intrasaisonnière et interannuelle dans les Océans Indien et Pacifique. El Niño ou le DOI sont des modes couplés : c'est la rétroaction positive découlant des interactions océan-atmosphère qui est source de variabilité (le " Bjerknes feedback "). À l'échelle intrasaisonnière, le rôle du couplage océan-atmosphère semble moins primordial, et modifie seulement des modes de variabilité essentiellement atmosphériques ou océaniques. Par exemple, les ondes d'instabilité dans le Pacifique Est sont le résultat d'une instabilité interne océanique. Cependant, elles affectent la stabilité atmosphérique, les vents de surface, et cela tend à réduire légèrement leur activité. À l'inverse, l'OMJ est un phénomène dont la source est atmosphérique, naissant du couplage entre dynamique et convection dans les tropiques. Toutefois, nous verrons que ce phénomène a une réponse océanique forte dans l'Océan Indien, à la fois en termes de dynamique et de thermodynamique. Le degré d'influence du couplage dans les propriétés de l'OMJ reste toutefois une question largement ouverte. Nous nous intéresserons aussi à la question des interactions entre ces différents modes de variabilité. Nous verrons par exemple comment la variabilité intrasaisonnière atmosphérique peut déclencher un El Niño, comment El Niño peut supprimer l'activité des ondes tropicales d'instabilité et l'effet retour, comment le DOI module l'activité de l'OMJ et enfin, quelles sont les interactions entre DOI et El Niño. Je présenterai alors une région de l'Océan Indien assez emblématique de ces interactions d'échelle, et dans laquelle j'ai développé une activité d'observations (campagnes océanographiques Cirene de 2005 à 2008 et projet de campagne TRIO). La bande 5°S-10°S dans l'Océan Indien est une région très particulière. En raison de la structure des vents, la thermocline y est proche de la surface et la couche de mélange est peu profonde, ce qui induit une forte réactivité de la température de surface aux sollicitations de l'atmosphère. De plus, la température de surface en hiver boréal est proche du seuil de convection, impliquant une sensibilité accrue de l'atmosphère à de petites variations de température. Ces deux facteurs augmentent le couplage océan atmosphère dans cette région qui a une variabilité très marquée aux échelles synoptiques (cyclones), intrasaisonnières (OMJ) et interannuelle (réponse à El Niño, mais aussi au DOI). Cette région a enfin des conséquences climatiques marquées (sur l'intensité des pluies de la mousson suivante, sur le nombre de cyclones dans le secteur La Réunion-Madagascar, sur la convection au-dessus du continent maritime, et même sur l'Amérique du Nord). Pour conclure, je présenterai ma réflexion sur mes axes de recherches futurs, ainsi que mes projets en termes de campagnes et réseaux d'observations.

Les apports atmosphériques sont suspectés de jouer un rôle important dans les processus biogéochimiques des régions océaniques HNLC. Des échantillons de dépôt atmosphérique total et d'aérosols ont été collectés en continu entre 2009 et 2010 dans l'océan austral indien, aux îles Kerguelen et de Crozet. Ils ont été mesurés par HR-ICP-MS. Les concentrations en surface des aérosols terrigènes sont basses et ne peuvent rendre compte des flux de dépôt de poussières qui sont plus élevés que ceux précédemment observés pour la région mais qui restent en adéquation avec les modèles atmosphériques. Le calcul des coefficients d'abattement suggère que l'aérosol mesuré en surface n'est pas représentatif de celui qui va générer le dépôt et qui est présent plus haut dans la colonne troposphérique, ce qui interdit ici l'estimation du flux à partir de ces mesures. Les flux de dépôt ont été établis pour une large gamme d'éléments identifiés comme venant de source terrigène ou marine. Pour certains d'entre eux (Pb, As, ...), une contribution anthropique a été observée pendant l'hiver austral. Les flux des éléments crustaux Al, Fe, Mn et Si sont similaires pour Kerguelen et Crozet ; ceux des autres éléments non marins sont différents d'un facteur deux à cinq avec un gradient décroissant de Crozet aux Kerguelen. Les rapports isotopiques de plomb montrent que le dépôt reçu aux Kerguelen provient principalement d'Amérique du sud ; le dépôt à Crozet semble fortement influencé par les apports africains. Les profils de REE pourraient confirmer ces observations. Les concentrations de métaux sont très basses dans les pluies avec une solubilité très élevée.

Le but principal de cette thèse était de réaliser une simulation du climat du 21e siècle, en tenant compte non seulement des émissions de gaz à effet de serre et d'aérosols produites par l'activité humaine, mais aussi des changements de végétation résultant du développement démographique et économique et des rétroactions climatiques. Pour réaliser cette simulation, il a été décidé d'utiliser le modèle d'impact intégré IMAGE2.2 (développé au RIVM, en Hollande) qui simule l'évolution de la concentration des gaz à effet de serre et de la végétation, et de le coupler au modèle climatique de circulation générale ARPEGE/OPA, disponible au CNRM. Avant de réaliser la simulation, des expériences de sensibilité de chacun des deux modèles au forçage produit par l'autre ont été réalisées. Finalement, la simulation couplée a montré qu'à l'échelle du siècle les interactions entre le climat et la végétation sont relativement peu importantes.

Depuis une vingtaine d'années, la qualité de la prévision de la trajectoire des cyclones tropicaux a fortement progressé mais peu d'améliorations ont été apportées à la prévision de l'intensité. Cette thèse s'intéresse aux changements d'intensité des cyclones tropicaux du sud-ouest de l'océan Indien sous l'angle original de l'activité électrique. Une première étude climatologique s'appuyant sur les données du réseau de détection d'éclairs World Wide Lightning Location Network (WWLLN) a permis de valider la qualité de ces données et de montrer que les cyclones tropicaux pouvaient être une source importante de production d'éclairs à l'échelle du bassin. À partir de ce résultat, une étude ciblée a été menée en considérant la spécificité du bassin et trois régions ont été définies (océan ouvert, région à proximité de la côte est de Madagascar et Canal du Mozambique). L'activité électrique des cyclones tropicaux dépend de la région où ils se trouvent ainsi que de leur stade d'intensité. Les éclairs semblent être un marqueur des phases d'intensification et d'affaiblissement dans certains situations. La troisième étape de cette thèse a consisté à simuler de manière idéalisée un cyclone tropical mature et a tenté d'expliquer les processus physiques à l'origine de l'activité électrique. Le modèle reproduit bien le comportement sporadique des éclairs habituellement observé. Des bilans montrent que ce cyclone présente trois phases caractéristiques, chacune associée à un comportement dynamique, microphysique et électrique différent. Enfin, une deuxième étude purement numérique a tenté d'identifier des estimateurs de l'activité électrique. Pour celà deux approches ont été adoptées : une étude globale et une étude par cellule. Dans les deux cas, les meilleurs estimateurs sont la masse totale de graupel, le volume d'updraft et le produit des flux des masses de glace précipitante et non précipitante.

Le but de cette thèse est de montrer quels peuvent être les apports à l'étude de la dynamique océanique des données de télédétection spatiale dans l'infrarouge proche (1000 nm) et moyen 10000 nm). On présente tout d'abord le rayonnement observé (réflexion, émission). On examine les influences de l'état de la mer et de la présence de l'atmosphère entre l'océan et le capteur ou le soleil sur ce rayonnement. Le chapitre 3 aborde les relations entre les distributions spatiales et superficielles de température et de courant. Le courant peut être estimé à l'aide d'un couple de thermographies, soit par suivi de structures particulières, soit par résolution de l'équation de conservation de la chaleur. Des exemples sont donnés pour le courant Ligure, les tourbillons anticycloniques issus du Gulf Stream et les eaux au large de la côte est des Etats-Unis. Enfin, les densités de variance de température calculées en fonction de l'échelle à partir des thermographies, sont comparées aux spectres d'énergie cinétique résultants de théories de turbulence bidimensionnelle. Le chapitre 4 traite des images d'énergie réfléchie par la surface de la mer à partir desquelles sont dérivés l'état d'agitation de la mer et le champ de vitesse du vent. Grâce à un modèle de distribution statistique des pentes des facettes de vague, on examine le signal parvenant au radiomètre en fonction de l'état de surface et des géométries d'éclairement et d'observation, en s'attachant aux cas des satellites SPOT et Landsat. Par ailleurs, on montre que le champ du module de la vitesse du vent à la surface de la mer peut être calculé, à partir des images AVHRR à 1000 nm et qu'il se compare favorablement au champ obtenu d'après une carte d'isobares de surface. Le chapitre 5 est une application des résultats des chapitres précédents, à la connaissance de la dynamique de la couche superficielle de la mer Ligure. La circulation générale dans ce bassin et ses causes présumées y sont discutées. Les variations, les instabilités *du courant et les dissymétries saisonnières nord-sud de la répartition spatiale de ces instabilités et de leur amplitude sont mises en relief. La conclusion est consacrée à la synthèse des apports de la télédétection à la connaissance du milieu marin mis en évidence par les travaux présentés, synthèse particulièrement illustrée par l'exemple de la mer Ligure.

Le problème de la paramétrisation des mécanismes turbulents dans les couches marines superficielles est étudié par comparaison des résultats obtenus avec différents modèles de turbulence dans diverses situations types. On montre que la description fine des mécanismes turbulents est profondément affectée par les hypothèses de modélisation retenues mais que la structure thermique prévue reste globalement identique dans les situations considérées. Plusieurs simulations examinant l'influence de la variabilité temporelle des échanges air-mer ont mis en évidence qu'une représentation fidèle de cette variabilité est le facteur critique dans ce problème.

Les travaux ont été effectués sur le golfe du Lion. Ils présentent des résultats combinés de l'interprétation d'images de satellites, obtenues dans l'infrarouge thermique (capteur AVHRR) et dans le domaine visible (capteur CZCS) avec des connaissances historiques et des mesures recueillies simultanément en mer. Ils permettent, d'une part, d'estimer la forte variabilité spatio-temporelle de la zone de dilution rhodanienne, en particulier sous l'influence du vent local et du débit du fleuve, et, d'autre part, d'observer la formation hivernale d'eau profonde dans le golfe du Lion. L'évolution saisonnière des zones de production dans le golfe du Lion est suivie à l'aide d'images de la concentration relative en pigments chlorophylliens et de la température de la mer. Le traitement d'images du spectre visible du satellite METEOSAT II permet l'estimation mensuelle de l'éclairement solaire global au niveau de la mer.

" ... Notre objectif n'est pas le satellite, mais ce que nous pouvons atteindre à travers lui. Ce qui nous importe est de savoir en quoi, comment et sous quelles conditions, il peut nous aider dans notre tâche. ... " (Mr PIKETTY, Allocution d'ouverture des journées nationales d'études : Utilisation pour l'océanographie des satellites d'observation de la Terre. Brest, 6-8 février 1978). Cette thèse constitue une réponse partielle à cette question. En effet, ce travail a toujours été guidé par le problème de l'utilisation des données VHRR des satellites NOAA4 et NOAAS, qui permettent de tracer des cartes de températures de surface, à des fins d'océanographie dynamique. De la lecture de cette thèse, peut résulter une certaine impression de dispersion et d'inachèvement. Ayant développé deux axes de recherche (variations saisonnières, variations spatiales) et ne les ayant traités qu'à l'aide des données satellite, il était fatal que les spécialistes de chaque sujet n'y trouvent pas leur compte. Si, pour les variations saisonnières, j'ai pu élaborer minutieusement une méthode d'étude, il n'en est pas de même pour le second axe de recherche. Il n'est pas habituel en océanographie de traiter des champs bi-dimensionnels. J'ai donc utilisé plus ou moins les méthodes usuelles d'étude des séries temporelles, en les alliant avec celles issues de la géostatique minière que j'ai apprises à l'Ecole des Mines. J'ai conscience des imperfections et des problèmes statistiques attachés à mon étude, mais, étant d'abord un océanographe, j'ai préféré ne pas m'y attarder et continuer, afin de pouvoir évaluer la richesse potentielle de cette étude spatiale. J'ai également essayé de calculer numériquement les courants d'advection thermique à partir des thermographies, à l'aide de l'équation de continuité de la quantité de chaleur. Mais un calcul d'incertitude m'a vite montré que le bruit du VHRR était trop important pour que les résultats aient une quelconque signification. L'utilisation prochaine des mesures moins bruitées du satellite TIROS-N autorisera une telle étude. L'utilisation de ces mesures de température par satellite pour l'océanographie dynamique s'avère très fructueuse. Elle a permis de compléter notre connaissance des phénomènes de surface dans le bassin liguro-provençal. Elle peut également apporter des éléments de solution à l'heure où la communauté océanographique se pose des questions sur la turbulence océanique à moyenne et grande échelle. Elle devrait aussi autoriser le calcul en routine des courants d'advection thermique. En conclusion, je pense que les satellites fournissent actuellement des mesures fiables, qui sont souvent des bons traceurs de mouvement, et que l'avènement des satellites plus performants du type TIROS-N va autoriser des analyses très fines, qui vont considérablement augmenter notre compréhension des phénomènes marins.

Nous avons mis en oeuvre une méthode statistique déjà existante (Cano, 1982) d'évaluation du rayonnement global à 1' aide des mesures du satellite Météosat. Cette méthode a été testée sur environ 1800 données : 11 mois de l'année 1983, 11 mois de l'année 1984 et 7 mois de l'année 1985. Les grands points de cette méthode statistique sont les suivants : 1) la construction d'une carte d'albédo planétaire par ciel clair d'après une série temporelle d'images sous conditions ciel clair, 2) l'évaluation d'un indice d'ennuagement, n, représentant le pourcentage du pixel recouvert par les nuages, par comparaison avec la carte d'albédo planétaire (albédo référence), 3) l'évaluation de la relation statistique entre la transmittance atmosphérique horaire (K), définie comme le rapport entre le rayonnement global et le rayonnement hors atmosphère (G0h), et n (l'indice d'ennuagement). Cette relation est du type : K = -a n + b ou Gh / G0h = -a n + b, où a et b > 0. Les coefficients de régression de cette relation ont été établis à partir des données du réseau de la Météorologie Nationale (Gh) et des données quotidiennes du canal visible de Météosat (n). En conclusion, nous pouvons dire que cette méthode est satisfaisante, elle permet 1) d'affiner ou de compléter (océans) les atlas du rayonnement déjà existants, en mettant en évidence des micro-climats (précision à 5 km près) peu connus ou inconnus, 2) l'évaluation du rayonnement global horaire moyen mensuel avec une précision de l'ordre de 0,1 kWh/m2, 3) un grand nombre d'application dites solaires (précision de rentabilité, méthode de conception architecturale assistée par micro-ordinateur, dimensionnement ...).