Calibration Techniques for Studying Venus and Mars Atmospheres

The European Space Agency Venus Express mission (VEX) was sent to Venus in 2005 to unveil the unsolved mysteries regarding its atmosphere, the plasma environment and its temperatures. Radio occultation experiments performed by VeRa radio science instrument probed the planet’s atmosphere by studying the frequency shift on the radio signal sent by the spacecraft to Earth-based ground stations. This study carries out the calibration of the radio frequencies within a radio occultation experiment in order to correct the main sources of error as: thermal noise, spacecraft clock, spacecraft trajectory, and plasma noise. Any uncalibrated effects will bias the retrieval of atmospheric properties. A comparison of the occultation experiments between Venus and Mars is presented, both from the engineering and scientific point of view, through the analysis of Venus Express and Mars Global Surveyor (MGS) occultations data, highlighting stronger calibrations required for VEX, the extreme, hostile, thick Venus’ atmosphere, and a friendly, thin Mars’ atmosphere. This investigation analyzes Venus Express data recorded by the NASA Deep Space Network in 2014, and the results are compatible to previous studies of Venus atmosphere with VEX between 2006 and 2009.


Introduction
Venus has been the prime target in the early interplanetary exploration since it is the planet with the closest approach to Earth. The first radio occultation experiment on Venus was performed in 1967 during a flyby of Mariner V [1][2][3], and a strong contributor to the understanding of this planet was made by Pioneer Venus Orbiter, Venera, Mariner missions and Magellan, which mapped the entire surface of the planet and its gravity up to the degree and order 180 [4]. Venus revealed similarities in size, density, mass, volume, orbital radius and bulk composition to Earth. However, the missions discovered an atmosphere characterized by extremely high temperatures, pressures, and composition which renders it uninhabitable, pointing out two planets that had evolved very differently. In 2005, the European Space Agency sent the Venus Express mission (VEX) to solve the unanswered questions regarding the past, present and future of Venus. One of the main objectives was the study of the atmosphere, ionosphere and gravity of Venus through the VeRa radio science instrument, using S-and X-band radio links (2.3 and 8.4 GHz, respectively) between the spacecraft and Earthbased ground stations [5]. The first VEX scientific results have been presented by Pätzold [6], showing vertical profiles of Venus ionosphere between 100 and 500 km altitude and a neutral atmosphere between 40 and 90 km altitude. Several studies have been published by Tellmann, Piccialli, Lee, Ando [7][8][9][10][11] regarding also the dynamics of Venus atmosphere as well as its thermal structure. Nowadays Akatsuki spacecraft from the Japanese Space Agency (JAXA) is continuing Venus' exploration and more missions will be sent there in the near future.
The exploration of Mars started in 1965 with the first flyby of Mariner IV [12], followed by MarinerVI and Mari-nerVII from NASA. The Soviet Union participated in the early exploration of Mars with Mars 2 and Mars 3 orbiters. Then, in the 1970s, NASA launched the Viking missions, followed by Mars Global Surveyor in 1996 [13]. MGS made global observations of Mars through its scientific payload, and it sounded Mars neutral atmosphere and ionosphere too. MGS and previous missions discovered a very thin Mars atmosphere, which starts at only 40 km altitude, and it is also characterized by temperatures and pressures which are not 1 3 as high as the ones on Venus [14,15]. Mars Express, Mars Odissey and Mars Reconnaissance Orbiter continued MGS observations to further expand Mars' atmosphere, climate and weather knowledge [16][17][18]. These missions highlighted a more habitable environment compared to Venus and, as a result, dozens of spacecrafts as orbiters, landers and rovers (which can survive on Mars but not on Venus due to the high surface temperatures and pressures) have been sent to study the Red Planet. In the next years, two new missions will explore Mars, Mars2020 from NASA and Exo Mars 2022 by ESA-Roscosmos.
The goal of this investigation is to show the main differences and similarities between Venus and Mars atmospheres, both from the engineering and scientific point of view, through radio occultation experiments, with a focus on the calibrations performed on the frequency residuals. Section 2 reports the structure of a radio occultation experiment, the adopted method and the developed atmosphere algorithm, the main parameters, as well as the main assumptions and limitations of this investigation. Section 3 introduces how the radio signal sent by the spacecraft is acquired by the Deep Space Network, and how the time-domain signal is converted in the frequency domain, to obtain the frequency residuals. Section 4 presents the frequency residuals calibration process, and the differences between VEX and MGS calibrations. To conclude, Section 5 shows the main scientific results and the comparison between Venus and Mars atmospheres.

Radio Occultation Method
Atmospheric radio occultation investigations rely on the detection of the frequency change of a radio signal as it travels through the atmosphere of a Solar System object. This is caused by the refraction and bending of the radio signal in the neutral gas and ionospheric plasma around the target object (see Fig. 1), which modify the frequency of the radio signal, making it different with respect to what expected by the receiver if no occultation would have occurred. In this way the frequency residuals can be evaluated, which are defined as the difference between the modified frequency during occultation and the expected direct frequency without refraction and occultation.
However, the measurements do not contain only the atmosphere information. The propagation of the radio signal is affected by interplanetary plasma along its path (which can not be removed in this investigation, since multiple frequency experiments would be required [19]), as well as Earth's ionosphere and troposphere, and errors related to the spacecraft clock, trajectory, and thermal noise [20]. The calibration process is a crucial step of this work, since it aims to compensate all the noises and errors, in order to obtain measurements characterized only by the target's atmosphere information.
Planets' atmospheres are studied through radio occultation data, by implementing the method of Withers [21] within an atmosphere retrieval algorithm, developed in MATLAB environment. In the frame of this investigation, the non-relativistic solution is the one considered, which does not take into account the relativity effects on the radio signal sent by the spacecraft. For both planets, the temperature differences between the relativistic and non-relativistic solution are always lower than 0.5 K, so the relativity effects are negligible. The main assumptions and limitations of this method are: 1. One-way link: the radio signal is transmitted from the spacecraft and received at the ground station; 2. Single frequency: only X-band signals are used in these experiments; 3. Stable frequency source: the source of the transmitted frequency must be stable (use of ultra-stable oscillators by the transmitter is key to the success of one-way link occultations); 4. Spherically symmetric atmosphere and ionosphere at the target: assumption needed to keep the radio signal within the r-z plane defined by Withers (2014) [21], otherwise gradients of refractivity will exists perpendicular to this plane and the radio signal will travel outside the plane; 5. Well-mixed atmosphere: assumption needed to use the hydrostatic equilibrium equation; 6. Ideal gas behavior: the neutral atmosphere is assumed to behave as an ideal gas.
This investigation is based on the atmospheric profiles retrieved from the X-band of both VEX and MGS. The Ray bending in the atmosphere of Venus. r 0 is the ray path closest approach distance; is the bending angle; a is the impact parameter; n is the index of refraction. The radii of the planet and the atmosphere are not to scale. [5] VEX data are composed by 25

Signal Processing
Radio occultation experiments rely on their ability to obtain scientific information of a planetary body through the radio frequency of the signal sent by the spacecraft to Earth. Open-loop receivers in the Deep Space Network record the raw antenna voltages of the downconverted signal, relying on frequency predicts to remain tuned to the incoming signal [22]. To get the radio frequencies of the recorded time-domain signal (also called sky frequencies), a spectral fast Fourier transform algorithm (FFT) is adopted, which estimates the main component frequency of the signal for each time step. Then, the sky frequencies are used to evaluate the frequency residuals, which are defined as: where f o is the sky frequency while f c is the transmitted frequency from the spacecraft corrected for the relative Doppler effects due to the known relative motion of the Earth antenna with respect to the spacecraft, so that only the atmospheric Doppler shifts are measured. These frequency residuals are the input for the atmosphere retrieval algorithm, which converts them in atmospheric parameters of the target body. The first difference between Venus and Mars, can be seen by looking at the frequency residuals of Figs. 2, 3. The thick and dense Venus' atmosphere has a stronger effect on the refraction of the radio signal, which result in frequency residuals as high as 5000 Hz, on the other hand the small bending caused by the thin Mars' atmosphere is reflected in residuals of only 1-2 Hz. The frequency residuals are characterized by three main regions: the baseline (between − 2400 s and − 200 s of Fig. 2) is the region where the signal is travelling outside the atmosphere of the planet, so the frequency residuals are expected to be flat with a zero mean value; the ionosphere signature (clearly visible at − 100 s of Fig. 3), whose charged particle advance the Doppler phase of the signal resulting in a positive frequency shift; the atmosphere signature (the slope of the frequency residuals), which represents the signal travelling inside the atmosphere of the planet. (

Calibration
The accuracy of a one-way radio occultation experiment depends on the stability of the transmitted signal frequency. Both VEX and MGS were equipped with an ultra-stable oscillator (USO), an instrument which guarantees a high precise control of the transmitted radio frequency in terms of its Allan Standard Deviation over a wide range of integration times. Both VEX and MGS USO had a frequency stability with an Allan deviation of 10 −13 at integration times 1-100 seconds [5,14]. However, several noise sources and errors affect the observed frequency  1 3 (i.e. thermal noise, media propagation noise, Earth's ionosphere and troposphere, spacecraft clock and estimated trajectory), which result in fluctuations in the evaluated frequency residuals and, as a consequence, errors in the retrieved atmospheric parameters. The occultation geometry is directly linked to the frequency residuals (see Equation 13 in [21]), which is also analyzed within the Abel transform to compute the refractivity index of the planet's atmosphere, so as a result the retrieved atmospheric parameters will be affected too.
In this work, the order of magnitude of the frequency residual uncertainties is estimated to be: 10 −3 Hz for the thermal noise, between 10 −3 and 10 −4 Hz for Earth' ionosphere and troposphere (the DSN provides calibrations to remove the delay caused by the Earth's ionosphere to a precision of approximately 2-5 TECU [24]), 10 −3 Hz for interplanetary plasma, 10 −3 Hz for the spacecraft USO [25,26]. The uncertainty due to the estimated spacecraft trajectory, as well as the one of the residual Earth' ionosphere and troposphere particles after calibration are not quantified herein.
For these reasons, the calibration process is crucial to compensate the observed frequency for the noises present in the recorded signal, providing reliable and accurate scientific results. Furthermore, these errors can be easily seen by looking at the baseline of the frequency residuals (see Fig. 4), which in an ideal experiment should be flat with zero-mean value. A zero-mean value baseline indicates that the frequency residuals are fully calibrated for all non-atmospheric sources, so that only the atmospheric refraction signature remains.
The first calibration is performed to correct the frequency residuals for the Earth's ionosphere and troposphere, which cause a frequency delay on the Doppler observables. These corrections are in the form of polynomials whose coefficients are provided by the Tracking System Analytic Calibration (TSAC) Group at JPL, based on Global Positioning Systems (GPS) data analysis. In particular, the troposphere polynomial is subtracted from the frequency residuals, while the ionosphere one is added, since the charged particle effect advances the Doppler phase.
Then, a second calibration is required to compensate for the noises errors left as spacecraft clock, estimated trajectory, plasma noise and thermal noise. This is made through a polynomial fitting on the baseline, with the goal to better estimate the transmitted frequency from the unperturbed signal travelling outside the planet's atmosphere. The polynomial coefficients of order n are retrieved from the baseline and they are used to generate a new calibration polynomial for the entire experiment time-span. This polynomial is then subtracted from the frequency residuals to correct the measurements, so that the final calibrated frequency residuals are obtained.
Equation 2 summarizes the entire calibration process: where t is the time related to the frequency residuals, T and I are the troposphere and ionosphere calibration polynomials, respectively, and the polynomial of order n corrects for the spacecraft clock, trajectory, thermal noise and plasma noise. The effect of the calibration on the baseline can be seen by comparing Figs. 4, 5. Even small errors in the knowledge of the trajectory could propagate because of the high dynamics (as gravity) when the spacecraft approaches the planet. Regarding VEX, this is especially pronounced in the ingress region (see Figs. 4, 5, 6) and as a result, higher order polynomials are required for VEX ingress occultations, while VEX egress cases are calibrated with first and second order

Venus and Mars Atmosphere Comparison
Venus and Mars atmospheres vertical profiles, retrieved within this work, are presented in Figs. 7, 8, 9. The results clearly show an extreme, hostile environment on Venus, caused by its thick and dense atmosphere, which starts at 100 km altitude. In addition, the results are limited to 38 km altitude, due to the critical refraction of Venus atmosphere: the strong bending, refraction, absorption and defocussing of the radio signal at this altitude affects the signal, which becomes so weak that reliable detection is no longer possible [27]. At this altitude temperatures and pressures are as high as 450 K and 4 bar (see Figs. 8,9). On the other hand, Mars' results present a friendly environment, characterized by a thin atmosphere which starts at 40 km altitude, temperatures between 150 and 300 K, and extremely low pressure as 0.003 bar at the surface (see Fig. 7). Thanks to its thin and sparse atmosphere, the critical refraction does not occur so that the whole atmosphere, until the Red Planet's surface, can be studied through radio occultation experiments.
Moreover, the results on VEX data are consistent with the ones presented by Tellmann (2009) and Pätzold (2007) [6,7]. To conclude, a temperature-latitude dependence is also discovered, resulting in higher temperatures at Venus' equator than at the poles, see Fig. 9.

Discussion
This study highlights the calibration of the radio frequencies residuals within radio occultation investigations, a process needed to correct the major sources of noises and errors which affect the frequency residuals measurements. The main process' requirement is related to the polynomial fitting, which must be performed only in the baseline of the frequency residuals, the region where the signal is not perturbed by the target's atmosphere. The calibration process is a crucial part of a radio occultation experiment, and it influences the reliability and the accuracy of the retrieved scientific results. However, to increase the accuracy of this work, the VEX estimated trajectory should be corrected, since the major errors found are attributable to this effect, especially for the VEX ingress cases.
Furthermore, the major differences and similarities between Venus and Mars are presented, both from the engineering and scientific point of view. The similarities are related to the relativity effects, which are negligible for both planets, and the atmosphere retrieval algorithm which behaves in the same way. However, major differences than similarities are discovered. First, in the frequencies calibration process, where higher order polynomials for the VEX residuals are adopted, especially for the ingress case due to the errors in the estimated trajectory and spacecraft clock. Second, all the results enlighten two completely opposite atmospheres. The Venus one, very high, thick, dense, with high pressures and extremely hot, and the Mars one, smaller in height, characterized by lower pressures and temperatures. As a matter of fact, planning a space exploration on Venus is extremely challenging; therefore, this is one of the reasons which explains why in the last decades Mars has been explored more than Venus.
Moreover, Venus throughout the 2014 always shows the same general trend on the retrieved temperature profiles and, in addition, these results are comparable to the ones presented by Pätzold (2007) and Tellmann (2009) [6,7]. This temperature trend is linked to Venus' high thermal inertia caused by its robust greenhouse effect, which traps the heat inside the atmosphere, but also to its zonal winds which transfer the heat in the mid-low altitudes in an efficient way. Together, these effects, tend to keep Venus' temperature trend unchanged throughout the years. On the other hand, Mars does not experience these phenomena, and the temperature trend undergoes more variations depending on several factors, as the relative position between Mars and the Sun and the day/night condition. To conclude, the analysis on VEX data highlight a temperature-latitude dependence, so that Venus' temperatures at the equator are higher than at the poles.