Several measurement campaigns took place at the ALOMAR observatory at Andøya, Northern Norway during July-November 2001 to validate ionospheric electron density and dry temperature profiles in the troposphere and lower stratosphere derived from radio occultation measurements of the low earth orbiter satellite CHAMP. For temperature sounding, three balloons are released around GPS satellite occultation events that occurred inbetween a distance of 200 km around Andøya. At altitudes of 7–20 km the CHAMP profile shows a positive mean deviation increasing with height by about 1.5-2 Kelvin/ 10 km overlayed with variations of ±1 K when compared to the radiosonde. Taking into account the previous and following radiosonde ascents the mean deviation seems to be of systematic nature due to the occultation principle or the retrieval algorithm and the variations are related to geographical variations of temperature and to the horizontal averaging by the radio occultation technique. During the period from mid July to mid August, four occultations for ionospheric soundings occurred. The values of the F2 layer calculated from the CHAMP derived electron density profiles are compared to the readings of the Alomar and Tromsø ionosondes for these times. Comparison shows that using the radio occultation technique electron densities of the maximum value layer are calculated inbetween the same order of magnitude as the ionosondes measurements, however, they overestimate it in the cases discussed here. / Zur Validierung von Elekronendichte- und Temperaturprofilen, abgeleitet aus Radiookkultationsmessungen von CHAMP wurden im Juli-November 2001 mehrere Messkampagnen am ALOMAR Institut auf Andøya, Nordnorwegen durchgeführt. Zur Temperatursondierung wurden drei aufeinanderfolgende Radiosonden um den Zeitpunkt von Okkultationsereignissen im Umkreis von weniger als 200 km gestartet. Das hier diskutierte Temperaturprofil von CHAMP zeigt im Höhenbereich 7-20 km eine mit der Höhe zunehmende positive Abweichung von ca. 1,5-2 K/10 km mit Variationen um ±1 K verglichen mit dem Temperaturprofil der zum Okkultationszeitpunkt fliegenden Radiosonde. Der Vergleich mit den vorhergehenden und nachfolgenden Sondierungen lässt darauf schließen, dass die mittlere Abweichung durch systematische Fehler des Okkultationsverfahrens oder den Retrieval-Algorithmus bedingt sind, die Variationen jedoch durch die örtliche Abweichung und die horizontale Mittelung des Messverfahrens. Während des Zeitraumes von Mitte Juli bis Mitte August ereigneten sich vier Okkultationen zur Sondierung der Ionosphäre. Von den abgeleiteten Elektronendichteprofilen werden jeweils die Werte der F2-Schicht mit den zur gleichen Zeit gemessenen Elektronendichten der Ionosonden auf Andøya und bei Tromsø verglichen. Der Vergleich zeigt, dass mit Hilfe der Radiookkultaktionstechnik die Elektronendichtewerte der F2-Schicht in der gleichen Größenordnung berechnet, in diesen konkreten Fällen jedoch überschätzt werden.

In the height range 10–30 km atmospheric gravity waves lead to periodic perturbations of the background temperature field in the order of 2-3 K, that are resolved in temperature profiles derived from radio occultation measurements. Due to the spherical symmetry assumption in the retrieval algorithm and the low horizontal resolution of the measurement damping in the amplitude and phase shift of the waves occurs leading to remarkable errors in the retrieved temperatures. The influence of the geometric wave parameters and the measurement geometry on plane gravity waves in the range 100-1000 km horizontal and 1-10 km vertical wavelength is investigated with a 2D model ranging ±1000 km around the tangent point and 10-50 km in height. The investigation shows, that with radio occultation measurements more than 90 % of the simulated waves can be resolved and more than 50% with amplitudes above 90%. But the geometrical parameters cannot be identified, since one signal can be attributed to different combinations of wave parameters and view angle. Even short waves with horizontal wavelengths below 200 km can be derived correctly in amplitude and phase if the vertical tilt is small or the view angle of the receiver satellite is in direction of the wave crests. / Atmosphärische Schwerewellen führen im Höhenbereich 10-30 km zu periodischen Störungendes Hintergrundtemperaturfeldes in der Größenordnung von 2-3 K, die in Temperaturprofilen aus Radiookkultationsmessungen aufgelöst werden. Aufgrund der sphärischen Symmetrieannahme im Retrievalverfahren und durch die niedrige horizontale Auflösung des Messverfahrens werden Phasenverschiebungen und Dämpfung der Amplitude verursacht, die zu beachtlichen Fehlern bei den abgeleiteten Temperaturen führen. Der Einfluss der geometrischen Wellenparameter und der Messgeometrie auf ebene Schwerewellen im Bereich 100-1000 km horizontale und 1-10 km vertikale Wellenlänge wird untersucht mit einem 2D-Modell, dass sich auf ein Gebiet von ±1000 km um den Tangentenpunkt und von 10-50 km in der Höhe erstreckt. Die Untersuchung zeigt, dass mit Radiookkultationsmessungen mehr als 90% der simulierten Wellen aufgelöst werden und mehr als 50% mit Amplituden oberhalb von 90% der ursprünglichen. Die geometrischen Parameter können jedoch nicht aus Einzelmessungen abgeleitet werden, da ein Signal zu verschiedenen Kombinationen von Wellenparametern und Sichtwinkel zugeordnet werden kann. Auch relativ kurze Wellen mit horizontalen Wellenlängen unterhalb von 200 km können korrekt in der Amplitude und Phase aufgelöst werden, falls die Neigung des Wellenvektors gegen die vertikale gering ist oder der Sichtwinkel des Empfängersatelliten in Richtung der Wellenberge ist.

La thèse, dont le sujet est "Modélisation de la propagation électromagnétique en milieux à gradient d'indice basée sur les faisceaux gaussiens - Application à la propagation en atmosphère réaliste et à la radio-occultation entre satellites" a été commencée le 2 décembre 2013, au Département ÉlectroMagnétisme et Radar (DEMR) de l'Onera de Toulouse et avec le laboratoire LAPLACE de l'Université Paul Sabatier. Elle est co-financée par l'ONERA et par la Région Midi-Pyrénées. L'encadrement a été assuré par Jérôme Sokoloff (Laplace/UPS, directeur de thèse), Alexandre Chabory (ENAC, co-directeur) et Vincent Fabbro (ONERA). L'École Doctorale est l' "École Doctorale Génie Électrique, Électronique, Télécommunications : du système au nanosystème". Le faisceau gaussien a été principalement utilisé dans la recherche scientifique afin d'étudier les systèmes optiques tels que les lasers. Des études plus rares et plus récentes ont proposé de l'utiliser pour modéliser la propagation des ondes sismiques. Ses propriétés spatiales et spectrales ont amené certains auteurs à étudier son utilisation dans des modèles de propagation atmosphériques. Cette thèse a consisté à développer un modèle, appelé GBAR (Gaussian Beam for Atmospheric Refraction), de propagation troposphérique réaliste et déterministe en utilisant le formalisme des faisceaux gaussiens. La démarche adoptée a consisté à reprendre les équations fondamentales introduites par Cerveny et Popov décrivant de façon itérative la propagation d'un faisceau gaussien en milieu inhomogène, sous hypothèse de haute fréquence (modèle asymptotique). De nouvelles équations ont été développées à partir d'elles pour obtenir une description analytique de la propagation d'un faisceau gaussien dans un milieu troposphérique décrit par les variations spatiales de l'indice de réfraction. L'hypothèse de base pour l'obtention de la formulation analytique est que le gradient de l'indice de réfraction peut être considéré vertical et constant au voisinage du faisceau. Les équations analytiques pour la description de la propagation d'un seul faisceau ont ensuite été étendues à la modélisation d'un champ quelconque dans un milieu troposphérique pouvant contenir de fortes variations du gradient d'indice, y compris des inversions de gradient. Ceci a été réalisé en couplant les équations analytiques avec la procédure de décomposition multi-faisceaux développée dans sa thèse pas Alexandre Chabory. Le modèle GBAR a été validé dans des milieux troposphériques réalistes issus de simulations du modèle météo méso-échelle WRF (Weather Research and Forecasting). Dans un troisième temps, le modèle a été utilisé pour simuler des inversions de données de radio-occultation. Des outils existent pour fournir un modèle d'interprétation de ces données pour estimer les propriétés physiques de l'atmosphère à partir des mesures en phase, amplitude, Doppler et délai des signaux GNSS transmis entre satellites en orbite autour de la Terre / The subject of this PhD thesis is " Electromagnetic propagation modeling in inhomogeneous media with refractive index gradients based on Gaussian beams - Application to realistic atmospheric propagation and radio occultation between satellites ". The study started on december 2nd, 2013 at the DEMR (Département Électromagnétisme et Radar) department of the ONERA research laboratory, in Toulouse, France. It was funded both by the ONERA and Région Midi-Pyrénées. It was supervised by Jérôme Sokoloff (LAPLACE/UPS, thesis director), Alexandre Chabory (ENAC, thesis co-director) and Vincent Fabbro (ONERA). The doctoral school was "École Doctorale Génie Électrique, Électronique, Télécommunications : du système au nanosystème ". The Gaussian beam was mostly used in scientific investigations to study optical systems such as lasers. Rarer and more recent works suggested the use of the Gaussian beam formalism in order to model the propagation of seismic waves. The properties of the Gaussian beam also led some authors to develop models for atmospheric propagation. In this thesis a model based on Gaussian beams called GBAR (Gaussian Beam for Atmospheric Refraction) was developped for tropospheric propagation in realistic and deterministic conditions. The scientific approach consisted in rewritting the fundamental equations introduced by Cerveny and Popov describing iteratively the propagation of a Gaussian beam in inhomogeneous media, under the high-frequency assumption (asymptotic model). New equations were derived from them in order to get analytical equations of the propagation of a Gaussian beam in inhomogeneous media described by the variations of the refractive index. The basic assumption under to get the analytical equations is to consider that the refractive index gradient is vertical and constant around the beam axis. The analytical equations that describe the propagation of a Gaussian beam were extended to model the propagation of an arbitrary field in a tropospheric medium with strong variations and inversions of the refractive index. This was done by coupling the analytical equations with the multibeam expansion procedure developped by Alexandre Chabory in his PhD thesis. The GBAR model was validated in tropospheric conditions, using refractive index grids from the WRF (Weather Research and Forecasting) mesoscale meteorological model. In the third and final phase, the GBAR model was used to simulate Radio Occultation data inversions. Tools exist to allow for interpretations of Radio Occultation data in order to estimate the physical properties of the atmosphere from measured phased, amplitude, Doppler shift and delay of GNSS signals transmitted between satellites orbiting around the Earth

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第21579号 / 理博第4486号 / 新制||理||1644(附属図書館) / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)教授 橋口 浩之, 教授 塩谷 雅人, 教授 秋友 和典 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

The GPS radio occultation technique is used to study sporadic E (Es) layer plasma irregularities of the Earth’s ionosphere on a global scale using COSMIC/FORMOSAT-3 satellite constellation GPS signal-to-noise ratio (SNR) profiles. The maximum deviation from the mean SNR is attributed to the height of the Es layer. Es are produced by ion convergence due to vertical wind shear in the presence of a horizontal component of the Earth magnetic field, while the wind shear is provided mainly by solar tides. Indeed, close correlation between Es and wind shear phases have already been found for the semidiurnal and terdiurnal tidal components. Here, we present the global distribution of quarterdiurnal (QDT) signatures in Es occurrence rates. We find that, in accordance with upward energy flux, negative vertical phase gradients of QDT Es signatures are observed. The maximum signal of QDT Es is found at altitudes above 100 km. In the southern hemisphere, maximum QDT Es occurrence rates are found in winter and during equinoxes. In the northern hemisphere, however, at altitudes above 100 km strong QDT activity is also visible. / Die GPS-Radiookkultationstechnik wird verwendet, um ionosphärische Plasmairregularitäten in Verbindung mit sporadischen E- (Es) Schichten auf globaler Skala zu untersuchen. Verwendet werden Signal-Rauschverhältnis- (SNR-) Profile. Die maximale Abweichung vom mittleren SNR wird der Höhe der Es-Schicht zugeordnet. Es werden durch Ionenkonvergenz aufgrund von vertikaler Windscherung in Anwesenheit einer horizontalen Komponente des Erdmagnetfeldes hervorgerufen, wobei die Scherung hauptsächlich durch solare Gezeiten verursacht wird. Tatsächlich wurden schon früher deutliche Übereinstimmungen zwischen dem Auftreten von Es und den Phasen der halb- und dritteltägigen Gezeiten gefunden. Hier stellen wir die globale Verteilung der vierteltägigen (QDT) Signaturen in Es- Auftretensraten vor. Es zeigt sich dass, in Übereinstimmung mit einem aufwärts gerichteten Energietransport, negative vertikale Phasengradienten der QDT in Es auftreten. Die maximale Auftretenswahrscheinlichkeit liegt bei Höhen oberhalb von 100 km. Auf der Südhemisphäre fällt das Maximum der QDT in Es-Auftretensraten in den Winter, während auf der Nordhemisphäre oberhalb von 100 km auch im Sommer starke QDT-Aktivität zu verzeichnen ist.

The molecular compositions of the atmospheres of the giant planets (Jupiter, Saturn, Uranus and Neptune) are fundamental to understanding the processes which formed these planets and the solar system as a whole. Microwave observations of these planets probe regions in their atmospheres from approximately 0.1 to several bars, a process otherwise unachievable by visible and infrared means. Many gases and various cloud layers influence the millimeter wave spectra of the outer planets; however phosphine and ammonia are the main microwave absorbers at Saturn at pressures less than two bars. Understanding the pressure induced absorption of both constituents at observational frequencies is therefore vital to the analysis of any observational data. Laboratory measurements have been conducted to measure the microwave absorptivity and refractivity of phosphine and ammonia at Ka-band (32-40 GHz) and W-band (94 GHz), under conditions characteristic of the atmosphere of Saturn. The results were used to verify the accuracy of the phosphine formalism created by Hoffman et. al (2001) for use at millimeter wavelengths. Based on the laboratory measurements conducted, new formalisms were also created to express the opacity of ammonia at the measured frequencies. An important method for the study of planetary atmospheres is the radio occultation experiment ??method that uses radio links between Earth, and the spacecraft which passes behind the planet. The Cassini mission to Saturn, which will be conducting such experiments at Ka-band as well as S (2.3 GHz) and X (8.4 GHz) bands, has prompted the development of a radio occultation simulator used to calculate excess Doppler shifts and attenuation profiles for Saturn, utilizing the newest formalisms for phosphine and ammonia. The results indicate that there will be unambiguous detection and profiling of phosphine and ammonia, and predictions are made for the pressures at which loss of signal is anticipated.

This thesis is concerned with the role of internal gravity waves (IGWs) in the stratospheric dynamics and variability demonstrating the effect of spatiotemporal distribution of their activity on the stratospheric dynamics and transport. The first part introduces a theoretical overview of the most recent as well as classical approaches used for description of the wave-mean interaction in the middle atmosphere. Methodology for an IGW analysis from the GPS radio occultation density data is described in the next chapter and the advantages of utilization of density data are listed. The third chapter presents results describing the peculiar dynamics and anomalous IGW activity in the Eastern Asia/Northern Pacific region. An important part is dedicated to a discussion of accuracy limits and usability of different IGW activity proxies. The possible impact of the localized IGW activity is investigated using a mechanistic middle and upper atmosphere model in the last chapter. Sensitivity simulations are used to demonstrate an important role of the spatial distribution of IGW activity for a formation of planetary waves and for the longitudinal variability of the Brewer-Dobson circulation. Implications for the middle atmospheric and climate change research are discussed along with consequences for parameterizations of...

Abstract This thesis studies the ionospheres of the Earth and Mars by combining the observations of versatile instruments providing information on various aspects of the planetary environments. The work on the terrestrial ionosphere focuses particularly on solar wind–magnetosphere–ionosphere coupling, while the work on the Martian ionosphere is based on the development of a new approach to analyse radio-occultation data to retrieve the atmospheric and ionospheric profiles. In the Earth's ionosphere, two papers study the effects of solar wind high-speed streams on the ionospheric F-region peak electron density and on cosmic noise absorption resulting from the precipitation of energetic (>30 keV) electrons into the D region. For the first paper, a modified version of the superposed epoch analysis method, called phase-locked superposed epoch analysis method, was developed. The main finding is that a depletion near the F-region peak takes place in the afternoon and evening sectors during high-speed stream events. This could be explained by an increase in the electron loss rate as a consequence of ion-neutral frictional heating, which enhances the ion temperature and leads to neutral atmosphere expansion. In addition, dayside and post-midnight F-peak electron density increases are observed, probably related to soft particle precipitation. The second study reveals that cosmic noise absorption occurs during up to 4 days after the arrival of a high-speed stream, as substorm-injected energetic electrons precipitate in the midnight to early-afternoon ionosphere, principally at auroral latitudes. A third study reports for the first time observations of a modulation of cosmic noise absorption at periods near 10 s, associated with pulsating aurora. This suggests that the energetic component of the precipitating ux is modulated consistently with the auroral (1–10 keV) energies. At Mars, radio-occultation experiments have been performed by the Mars Express spacecraft since 2004. In this thesis, a new data analysis approach is developed, based on the numerical simulation of radio wave propagation through modelled Martian atmosphere and ionosphere. This approach enables one to overcome limitations inherent in the classical inversion method which has been in use so far. It also gives access to new parameters such as ion density profiles. The new method is tested by analysing the data from two radio-occultation experiments. / Tiivistelmä Tämä väitöskirja tutkii Maapallon ja Marsin ionosfäärejä yhdistämällä useiden eri instrumenttien havaintoja, joilla saadaan tietoa planeettojen ympäristöistä. Maapallon ionosfääriä koskeva työ tutkii aurinkotuuli–magnetosfääri–ionosfäärikytkentää, kun taas Marsin ionosfääriä koskevan työn tavoite on uuden radio-okkultaatiomittauksen data-analyysimenetelmän kehittäminen, joka tuottaa ilmakehän ja ionosfäärin profiileja. Maan ionosfäärin tapauksessa yhdessä julkaisussa tutkitaan nopeiden aurinkotuulivirtausten vaikutuksia F-kerroksen elektronitiheyteen ja toisessa julkaisussa tutkitaan energeettisten (>30 keV) elektronien sateesta johtuvaa kosmisen radiokohinan absorptiota D-kerroksessa. Ensimmäisessä julkaisussa on kehitetty uusi versio data-analyysimenetelmästä, jota kutsutaan vaihelukituksi epookkien superpositiomenetelmäksi. Julkaisun päätulos on, että nopeiden aurinkotuulivirtausten aikana F-kerroksen maksimielektronitiheys pienenee iltapäivän ja illan sektoreilla. Tämä voidaan selittää johtuvan siitä, että ioni-neutraalitörmäysten synnyttämä kitkalämpö kasvattaa ionilämpötilaa ja aiheuttaa lisäksi ilmakehän laajenemisen. Molemmat prosessit kasvattavat elektronien häviönopeutta. F-kerroksen elektronitiheysmaksimi puolestaan kasvaa sektorilla, joka ulottuu keskiyöstä aamun kautta keskipäivään, ja tämä johtuu matalaenergeettisestä elektronisateesta. Toisessa julkaisussa havaitaan, että lisääntynyt kosmisen radiokohinan absorptio kestää jopa neljä päivää nopean aurinkotuulivirtauksen saavuttua Maan kohdalle. Tämä johtuu siitä, että alimyrskyitse injektoidut energeettiset elektronit satavat keskiyön ja aamun ionosfääriin, pääasiassa revontuliovaalin alueella. Kolmas julkaisu raportoi ensimmäistä kertaa havainnon sykkiviin revontuliin liittyvästä kosmisen radiokohinan absorptiosta n. 10 s jaksollisuudella. Tämä osoittaa, että elektronivuon energeettinen komponentti on moduloitu samalla jaksollisuudella kuin revontulielektronien energiat (1–10 keV). Marsissa on tehty radio-okkultaatiomittauksia vuodesta 2004 saakka Mars Express -luotaimen avulla. Vaitoskirjassa on kehitetty uusi datan analyysimenetelmä, joka perustuu numeeriseen simulointiin radioaaltojen etenemisestä Marsin ilmakehässä ja ionosfäärissä. Tämän lähestymistavan avulla vältetään tähän asti käytetyn klassisen inversiomenetelmän rajoitukset. Lisäksi menetelmä tuottaa uusia parametrejä kuten ionitiheysprofiileja. Uutta menetelmää testattiin tulkiten kahden radio-okkultaatiomittauksen aineistoa. / Résumé Le travail présenté dans ce manuscrit de thèse s'articule autour de l'étude des ionosphères terrestre et martienne. Une approche multi-instrumentale est adoptée afin de combiner des observations permettant de mettre en perspective des manifestations de phénomènes physiques de natures différentes mais intervenant dans un même contexte global. Le travail doctoral comporte également un volet modélisation. Le manuscrit de thèse consiste en une partie introductrice à laquelle sont adossées cinq publications dans des revues scientifiques à comité de lecture. La partie introductrice de ce manuscrit de thèse a pour objectif de présenter le contexte scientifique sur lequel est basé le travail doctoral. Un premier chapitre passe en revue les principaux aspects théoriques dans lesquels s'inscrivent les études dont les résultats sont publiés dans les cinq articles. Les atmosphères et ionosphères de la Terre et de Mars y sont succinctement décrites, de même que les interactions entre ces planètes et le vent solaire, comprenant notamment la formation de magnétosphères. Les deux chapitres suivants présentent les instruments dont sont issues les données utilisées dans ce travail doctoral ainsi que les méthodes d'analyse des données. Le quatrième chapitre résume les principaux résultats obtenus autour des trois grandes thématiques abordées au cours de cette thèse. Enfin, des pistes quant à la continuation potentielle du travail présenté dans ce manuscrit de thèse sont évoquées en conclusion. Le premier article porte sur une étude statistique des effets des courants de vent solaire rapide sur la région F de l'ionosphère aurorale. Il s'appuie sur des données mesurées par l'ionosonde de Sodankylä entre 2006 et 2008. Au cours de cette période, 95 événements associés à des courants de vent solaire rapide ont été sélectionnés, et la réponse de l'ionosphère au-dessus de Sodankylä a été étudiée à partir des fréquences critiques des régions E et F de l'ionosphère, qui donnent la valeur du pic de concentration électronique dans ces deux régions. Pour cela, une version modifiée de la méthode des époques superposées a été développée, appelée “méthode des époques superposées avec verrouillage de phase”. Une augmentation du pic de concentration des régions E et F est observée du côté nuit et le matin, en lien avec une activité aurorale accrue, tandis qu'une déplétion de la région F est révélée aux temps magnétiques locaux situés entre 12 h et 23 h. Une estimation des effets d'une possible modification de l'équilibre photo-chimique résultant d'un accroissement du chauffage issu de la friction entre les ions et les éléments neutres est proposée. Le deuxième article s'intéresse aux précipitations énergétiques dans l'ionosphère aurorale durant ces mêmes 95 événements, en étudiant l'absorption du bruit cosmique qui en résulte. Il apparaît que les événements au cours desquels le vent solaire demeure rapide pendant plusieurs jours produisent davantage de précipitations énergétiques, qui peuvent atteindre les latitudes subaurorales. Par ailleurs, trois types de précipitations énergétiques sont étudiés séparément, selon qu'elles sont associées avec des signatures de sous-orage magnétique, avec des pulsations géomagnétiques, ou ni l'un ni l'autre. Les deux premiers types de précipitations semblent liés. En effet, l'analyse des données suggère que les électrons énergétiques sont injectés dans la magnétosphère interne durant les sous-orages. Tandis qu'une partie d'entre eux précipitent immédiatement du côté nuit, d'autres dérivent vers le côté matin, où ils subissent des interactions avec des ondes électromagnétiques de type siffleur (whistler en anglais), qui peuvent être modulées par des pulsations géomagnétiques, menant à leur précipitation. Le troisième article présente pour la première fois l'observation de signatures d'aurore pulsante dans des données d'absorption du bruit cosmique. Ces signatures sont consistantes avec les pulsations observables dans l'émission aurorale, et semblent indiquer une modulation cohérente des composantes aurorale (1–10 keV) et énergétique (> 30 keV) du spectre des précipitations électroniques au cours d'une aurore pulsante. Le quatrième article introduit une nouvelle méthode proposée pour analyser les données de radio-occultation mesurées par la sonde Mars Express. Cette approche vise à contourner des difficultés posées par les hypothèses fortes nécessaires à la mise en œuvre de la méthode classique d'inversion, notamment celle d'un environnement martien à symétrie sphérique — qui n'est pas acceptable lors de sondages proches du terminateur jour-nuit. La nouvelle méthode est basée sur la modélisation de l'atmosphère et de l'ionosphère de Mars, et sur la simulation de la propagation des ondes radio entre la station sol sur Terre et Mars Express lors d'une expérience de radio-occultation. Les paramètres libres contrôlant les profils atmosphériques et ionosphériques sont ajustés afin que la simulation reproduise le plus fidèlement possible les mesures. Le cinquième article est une réponse à un commentaire sur l'article précédent. Il vise d'une part à répondre aux critiques émises sur la méthode développée, montrant que celles-ci n'en remettent en cause ni la validité ni la pertinence, et d'autre part à y apporter quelques améliorations.