Descent motions of the Huygens probe as measured by the Surface Science Package (SSP): Turbulent evidence for a cloud layer

The Huygens probe underwent vigorous short-period motions during its parachute descent through the atmosphere of Saturn's moon Titan in January 2005, at least some of which were excited by the Titan environment. Several sensors in the Huygens Surface Science Package (SSP) detect these motions, indicating the transition to the smaller stabilizer parachute, the changing probe spin rate, aerodynamic buffeting, and pendulum motions. Notably, in an altitude range of about 20–30 km where methane drops will freeze, the frequency content and statistical kurtosis of the tilt data indicate excitation by turbulent air motions like those observed in freezing clouds on Earth, supporting the suggestion of Tokano et al. [Tokano, T., McKay, C.P., Neubauer, F.M., Atreya, S.K., Ferri, F., Fulchignoni, M., Niemann, H.B. (2006a). Methane drizzle on Titan. Nature 442, 432–435] that the probe passed through such a cloud layer. Motions are weak below 20 km, suggesting a quiescent lower atmosphere with turbulent fluctuations of nominally <0.15 m/s (to within a factor of ∼2) but more violent motions in the upper troposphere may have been excited by turbulent winds with amplitudes of 1–2 m/s. Descent in part of the stratosphere (150–120 km) was smooth despite strong ambient wind (∼100 m/s), but known anomalies in the probe spin prevent investigation of turbulence in the known wind-shear layer from 60 to 100 km.     

Titan global climate model: A new 3-dimensional version of the IPSL Titan GCM

We have developed a new 3-dimensional climate model for Titan’s atmosphere, using the physics of the IPSL Titan 2-dimensional climate model with the current version of the LMDZ General Circulation Model dynamical core. Microphysics and photochemistry are still computed as zonal averages. This GCM covers altitudes from surface to 500 km altitude, with barotropic waves now being resolved and the diurnal cycle included. The boundary layer scheme has been changed, yielding a strong improvement in the tropospheric zonal wind profile modeled at Huygens descent position and season. The potential temperature profile is fairly consistent with Huygens observations in the lowest 10 km. The latitudinal profile of the near-surface temperature is close to observed values. The minimum of zonal wind observed by the Huygens probe just above the tropopause is also present in these simulations, and its origin is discussed by comparing solar heating and dynamical transport of energy. The stratospheric temperature and wind fields are consistent with our previous works. Compared to observations, the zonal wind peak is too weak (around 120 m/s) and too low (around 200 km). The temperature structures appear to be compressed in altitude, and depart strongly from observations in the upper stratosphere. These discrepancies are correlated, and most probably related to the altitude of the haze production. The model produces a detached haze layer located more than 150 km lower than observed by the Cassini instruments. This low production altitude is due to the current position of the GCM upper boundary. However, the temporal behaviour of the detached haze layer in the model may explain the seasonal differences observed between Cassini and Voyager 1. The waves present in the GCM are analyzed, together with their respective roles in the angular momentum budget. Though the role of the mean meridional circulation in momentum transport is similar to previous work, and the transport by barotropic waves is clearly seen in the stratosphere, a significant part of the transport at high latitudes is done all year long through low-frequency tropospheric waves that may be baroclinic waves.     

Analytic theory of Titan’s Schumann resonance: Constraints on ionospheric conductivity and buried water ocean

This study presents an approximate model for the atypical Schumann resonance in Titan’s atmosphere that accounts for the observations of electromagnetic waves and the measurements of atmospheric conductivity performed with the Huygens Atmospheric Structure and Permittivity, Wave and Altimetry (HASI–PWA) instrumentation during the descent of the Huygens Probe through Titan’s atmosphere in January 2005. After many years of thorough analyses of the collected data, several arguments enable us to claim that the Extremely Low Frequency (ELF) wave observed at around 36 Hz displays all the characteristics of the second harmonic of a Schumann resonance. On Earth, this phenomenon is well known to be triggered by lightning activity. Given the lack of evidence of any thunderstorm activity on Titan, we proposed in early works a model based on an alternative powering mechanism involving the electric current sheets induced in Titan’s ionosphere by the Saturn’s magnetospheric plasma flow. The present study is a further step in improving the initial model and corroborating our preliminary assessments. We first develop an analytic theory of the guided modes that appear to be the most suitable for sustaining Schumann resonances in Titan’s atmosphere. We then introduce the characteristics of the Huygens electric field measurements in the equations, in order to constrain the physical parameters of the resonating cavity. The latter is assumed to be made of different structures distributed between an upper boundary, presumably made of a succession of thin ionized layers of stratospheric aerosols spread up to 150 km and a lower quasi-perfect conductive surface hidden beneath the non-conductive ground. The inner reflecting boundary is proposed to be a buried water–ammonia ocean lying at a likely depth of 55–80 km below a dielectric icy crust. Such estimate is found to comply with models suggesting that the internal heat could be transferred upwards by thermal conduction of the crust, while convective processes cannot be ruled out.     

A method to determine the atmospheric temperature profile from in situ pressure data: Application to Titan

We have developed a new method to analyse in situ observations of atmospheric variables of state: the reconstruction of the vertical temperature profile from pressure measurements accompanied by rough knowledge of the atmospheric composition and the aerodynamical response properties of the descent vehicle. We can use the method to construct the temperature profile when no direct measurements are available, as well as to analyse the consistency between data from different instruments. We applied the method to the Huygens measurements of Titan's atmosphere, determining the aerodynamical drag properties from radar altimeter data. We discovered that the temperature profile computed in this manner differs from the profile from the temperature sensor (TEM) of the probe by up to 5% in the altitude range of 0–60 km, and up to 10% at higher altitudes due to increased noise. The method gives a tropopause altitude of about 50 km and a surface temperature of about 98 K, in contrast to the TEM temperature measurements. Our error analysis shows that these differences are caused by the known discrepancy in the Huygens altimeter data, with the estimates made by the reconstruction algorithm contributing only 1–2% of error.     

Radiative transfer analyses of Titan’s tropical atmosphere

Titan’s optical and near-IR spectra result primarily from the scattering of sunlight by haze and its absorption by methane. With a column abundance of 92 km amagat (11 times that of Earth), Titan’s atmosphere is optically thick and only ～10% of the incident solar radiation reaches the surface, compared to 57% on Earth. Such a formidable atmosphere obstructs investigations of the moon’s lower troposphere and surface, which are highly sensitive to the radiative transfer treatment of methane absorption and haze scattering. The absorption and scattering characteristics of Titan’s atmosphere have been constrained by the Huygens Probe Descent Imager/Spectral Radiometer (DISR) experiment for conditions at the probe landing site (Tomasko, M.G., Bézard, B., Doose, L., Engel, S., Karkoschka, E. [2008a]. Planet. Space Sci. 56, 624–247; Tomasko, M.G. et al. [2008b]. Planet. Space Sci. 56, 669–707). Cassini’s Visual and Infrared Mapping Spectrometer (VIMS) data indicate that the rest of the atmosphere (except for the polar regions) can be understood with small perturbations in the high haze structure determined at the landing site (Penteado, P.F., Griffith, C.A., Tomasko, M.G., Engel, S., See, C., Doose, L., Baines, K.H., Brown, R.H., Buratti, B.J., Clark, R., Nicholson, P., Sotin, C. [2010]. Icarus 206, 352–365). However the in situ measurements were analyzed with a doubling and adding radiative transfer calculation that differs considerably from the discrete ordinates codes used to interpret remote data from Cassini and ground-based measurements. In addition, the calibration of the VIMS data with respect to the DISR data has not yet been tested. Here, VIMS data of the probe landing site are analyzed with the DISR radiative transfer method and the faster discrete ordinates radiative transfer calculation; both models are consistent (to within 0.3%) and reproduce the scattering and absorption characteristics derived from in situ measurements. Constraints on the atmospheric opacity at wavelengths outside those measured by DISR, that is from 1.6 to 5.0 μm, are derived using clouds as diffuse reflectors in order to derive Titan’s surface albedo to within a few percent error and cloud altitudes to within 5 km error. VIMS spectra of Titan at 2.6–3.2 μm indicate not only spectral features due to CH₄ and CH₃D (Rannou, P., Cours, T., Le Mouélic, S., Rodriguez, S., Sotin, C., Drossart, P., Brown, R. [2010]. Icarus 208, 850–867), but also a fairly uniform absorption of unknown source, equivalent to the effects of a darkening of the haze to a single scattering albedo of 0.63 ± 0.05. Titan’s 4.8 μm spectrum point to a haze optical depth of 0.2 at that wavelength. Cloud spectra at 2 μm indicate that the far wings of the Voigt profile extend 460 cm^?1 from methane line centers. This paper releases the doubling and adding radiative transfer code developed by the DISR team, so that future studies of Titan’s atmosphere and surface are consistent with the findings by the Huygens Probe. We derive the surface albedo at eight spectral regions of the 8 × 12 km² area surrounding the Huygens landing site. Within the 0.4–1.6 μm spectral region our surface albedos match DISR measurements, indicating that DISR and VIMS measurements are consistently calibrated. These values together with albedos at longer 1.9–5.0 μm wavelengths, not sampled by DISR, resemble a dark version of the spectrum of Ganymede’s icy leading hemisphere. The eight surface albedos of the landing site are consistent with, but not deterministic of, exposed water ice with dark impurities.     

Dynamical properties of the Venus mesosphere from the radio-occultation experiment VeRa onboard Venus Express

The dynamics of Venus’ mesosphere (60–100 km altitude) was investigated using data acquired by the radio-occultation experiment VeRa on board Venus Express. VeRa provides vertical profiles of density, temperature and pressure between 40 and 90 km of altitude with a vertical resolution of few hundred meters of both the Northern and Southern hemisphere. Pressure and temperature vertical profiles were used to derive zonal winds by applying an approximation of the Navier–Stokes equation, the cyclostrophic balance, which applies well on slowly rotating planets with fast zonal winds, like Venus and Titan. The main features of the retrieved winds are a midlatitude jet with a maximum speed up to 140 ± 15 m s^?1 which extends between 20°S and 50°S latitude at 70 km altitude and a decrease of wind speed with increasing height above the jet. Cyclostrophic winds show satisfactory agreement with the cloud-tracked winds derived from the Venus Monitoring Camera (VMC/VEx) UV images, although a disagreement is observed at the equator and near the pole due to the breakdown of the cyclostrophic approximation. Knowledge of both temperature and wind fields allowed us to study the stability of the atmosphere with respect to convection and turbulence. The Richardson number Ri was evaluated from zonal field of measured temperatures and thermal winds. The atmosphere is characterised by a low value of Richardson number from ～45 km up to ～60 km altitude at all latitudes that corresponds to the lower and middle cloud layer indicating an almost adiabatic atmosphere. A high value of Richardson number was found in the region of the midlatitude jet indicating a highly stable atmosphere. The necessary condition for barotropic instability was verified: it is satisfied on the poleward side of the midlatitude jet, indicating the possible presence of wave instability.     

Atomic oxygen on the Venus nightside: Global distribution deduced from airglow mapping

The Visible and Infra-Red Thermal Imaging Spectrometer (VIRTIS) instrument on board the Venus Express spacecraft has measured the O₂(a¹Δ) nightglow distribution at 1.27 μm in the Venus mesosphere for more than two years. Nadir observations have been used to create a statistical map of the emission on Venus nightside. It appears that the statistical 1.6 MR maximum of the emission is located around the antisolar point. Limb observations provide information on the altitude and on the shape of the emission layer. We combine nadir observations essentially covering the southern hemisphere, corrected for the thermal emission of the lower atmosphere, with limb profiles of the northern hemisphere to generate a global map of the Venus nightside emission at 1.27 μm. Given all the O₂(a¹Δ) intensity profiles, O₂(a¹Δ) and O density profiles have been calculated and three-dimensional maps of metastable molecular and atomic oxygen densities have been generated. This global O density nightside distribution improves that available from the VTS3 model, which was based on measurements made above 145 km. The O₂(a¹Δ) hemispheric average density is 2.1 × 10⁹ cm^?3, with a maximum value of 6.5 × 10⁹ cm^?3 at 99.2 km. The O density profiles have been derived from the nightglow data using CO₂ profiles from the empirical VTS3 model or from SPICAV stellar occultations. The O hemispheric average density is 1.9 × 10¹¹ cm^?3 in both cases, with a mean altitude of the peak located at 106.1 km and 103.4 km, respectively. These results tend to confirm the modeled values of 2.8 × 10¹¹ cm^?3 at 104 km and 2.0 × 10¹¹ cm^?3 at 110 km obtained by Brecht et al. [Brecht, A., Bougher, S.W., Gérard, J.-C., Parkinson, C.D., Rafkin, S., Foster, B., 2011a. J. Geophys. Res., in press] and Krasnopolsky [Krasnopolsky, V.A., 2010. Icarus 207, 17–27], respectively. Comparing the oxygen density map derived from the O₂(a¹Δ) nightglow observations, it appears that the morphology is very different and that the densities obtained in this study are about three times higher than those predicted by the VTS3 model.     

The abundance and vertical distribution of the unknown ultraviolet absorber in the venusian atmosphere from analysis of Venus Monitoring Camera images

Observations of Venus using the ultraviolet filter of the Venus Monitoring Camera (VMC) on ESA’s Venus Express Spacecraft (VEX) provide the best opportunity for study of the spatial and temporal distribution of the venusian unknown ultraviolet absorber since the Pioneer Venus (PV) mission. We compare the results of two sets of 125 radiative transfer models of the upper atmosphere of Venus to each pixel in a subset of VMC UV channel images. We use a quantitative best fit criterion based upon the notion that the distribution of the unknown absorber should be independent of the illumination and observing geometry. We use the product of the cosines of the incidence and emission angles and search for absorber distributions that are uncorrelated with this geometric parameter, finding that two models can describe the vertical distribution of the unknown absorber. One model is a well-mixed vertical profile above a pressure level of roughly 120 mb (～63 km). This is consistent with the altitude of photochemical formation of sulfuric acid. The second model describes it as a thin layer of pure UV absorber at a pressure level roughly around 24 mb (～71 km) and this altitude is consistent with the top of upper cloud deck. We find that the average abundance of unknown absorber in the equatorial region is 0.21 ± 0.04 optical depth and it decreases in the polar region to 0.08 ± 0.05 optical depth at 365 nm.     

Titan atmosphere profiles from Huygens engineering (temperature and acceleration) sensors

Data from several Huygens probe housekeeping sensors (an engineering accelerometer and housekeeping temperature sensors) are studied to determine how effectively such nonideal instruments may characterize the density or temperature structure of the atmosphere. While only confirming the results of the dedicated atmospheric structure instrument, this exercise is of relevance to possible future missions to various bodies which might not be equipped with such science-grade sensors able to accurately profile the atmosphere top-to-bottom. It is found that for typical engineering accelerometers with 8-bit resolution, the atmosphere density for ∼4 scale heights above the peak deceleration altitude may be recovered. If, as with Huygens, the peak deceleration exceeds the range of the accelerometers, recovery of an additional scale height or so below the peak is still possible, but relies on accurate total velocity knowledge. Engineering temperature sensors can, with care, be analyzed to recover the temperature structure of at least the lowest ∼30 km of Titan's atmosphere. Fortunately, in the case of Huygens, data from the surface after landing were available to constrain models of heat leaks which offset the observed temperature from that of the ambient air during descent; data from before and during the entry phase on other missions would be similarly useful. When corrections are made for the estimated heat transfer processes, the atmospheric temperature can be recovered to within about 3 K.     

Vertical profiling of SO2 and SO above Venus’ clouds by SPICAV/SOIR solar occultations

New measurements of sulfur dioxide (SO₂) and monoxide (SO) in the atmosphere of Venus by SPICAV/SOIR instrument onboard Venus Express orbiter provide ample statistics to study the behavior of these gases above Venus’ clouds. The instrument (a set of three spectrometers) is capable to sound atmospheric structure above the clouds in several observation modes (nadir, solar and stellar occultations) either in the UV or in the near IR spectral ranges. We present the results from solar occultations in the absorption ranges of SO₂ (190–230 nm, and at 4 μm) and SO (190–230 nm). The dioxide was detected by the SOIR spectrometer at the altitudes of 65–80 km in the IR and by the SPICAV spectrometer at 85–105 km in the UV. The monoxide’s absorption was measured only by SPICAV at 85–105 km. We analyzed 39 sessions of solar occultation, where boresights of both spectrometers are oriented identically, to provide complete vertical profiling of SO₂ of the Venus’ mesosphere (65–105 km). Here we report the first firm detection and measurements of two SO₂ layers. In the lower layer SO₂ mixing ratio is within 0.02–0.5 ppmv. The upper layer, also conceivable from microwave measurements by Sandor et al. (Sandor, B.J., Todd Clancy, R., Moriarty-Schieven, G., Mills, F.P. [2010]. Icarus 208, 49–60) is characterized by SO₂ increasing with the altitude from 0.05 to 2 ppmv, and the [SO₂]/[SO] ratio varying from 1 to 5. The presence of the high-altitude SO_x species could be explained by H₂SO₄ photodissociation under somewhat warmer temperature conditions in Venus mesosphere. At 90–100 km the content of the sulfur dioxide correlates with temperature increasing from 0.1 ppmv at 165–170 K to 0.5–1 ppmv at 190–192 K. It supports the hypothesis of SO₂ production by the evaporation of H₂SO₄ from droplets and its subsequent photolysis at around 100 km.     

Sulfur chemistry in the middle atmosphere of Venus

Venus Express measurements of the vertical profiles of SO and SO₂ in the middle atmosphere of Venus provide an opportunity to revisit the sulfur chemistry above the middle cloud tops (～58 km). A one dimensional photochemistry-diffusion model is used to simulate the behavior of the whole chemical system including oxygen-, hydrogen-, chlorine-, sulfur-, and nitrogen-bearing species. A sulfur source is required to explain the SO₂ inversion layer above 80 km. The evaporation of the aerosols composed of sulfuric acid (model A) or polysulfur (model B) above 90 km could provide the sulfur source. Measurements of SO₃ and SO (a¹Δ → X³Σ^-) emission at 1.7 μm may be the key to distinguish between the two models.     

Direct wind measurements from November 2007 in Venus’ upper atmosphere using ground-based heterodyne spectroscopy of CO2 at 10 μm wavelength

Between November 23 and 28, 2007, the Cologne Tuneable Heterodyne Infrared Spectrometer THIS was installed at the McMath-Pierce Solar Telescope (Kitt Peak, Arizona, USA) to determine zonal wind velocities and to estimate the subsolar-to-antisolar flow. We investigate dynamics in the upper atmosphere of Venus by measuring the Doppler shift of fully-resolved non-LTE CO₂ emission lines at 959.3917 cm^?1 (10.423 μm), which probe a narrow altitude region in Venus’ atmosphere around 110 ± 10 km (～1 μbar). The results show no significant zonal wind velocity at the equator. An increase with latitude up to 43 ± 13 m/s at a latitude of 33°N was observed. This confirms the deduction of a minor influence of Venus superrotation at an altitude of 110 km from previous measurements in May 2007 (Sornig et al., 2008). The specific observing geometry enables estimating the maximum cross terminator velocity of the subsolar-to-antisolar flow at 72 ± 47 m/s.     

The planetary fourier spectrometer (PFS) onboard the European Venus Express mission

The planetary fourier spectrometer (PFS) for the Venus Express mission is an infrared spectrometer optimized for atmospheric studies. This instrument has a short wavelength (SW) channel that covers the spectral range from 1700 to 11400 cm⁻¹ (0.9–5.5 μm) and a long wavelength (LW) channel that covers 250–1700 cm⁻¹ (5.5–45 μm). Both channels have a uniform spectral resolution of 1.3 cm⁻¹. The instrument field of view FOV is about 1.6 ° (FWHM) for the short wavelength channel and 2.8 ° for the LW channel which corresponds to a spatial resolution of 7 and 12 km when Venus is observed from an altitude of 250 km. PFS can provide unique data necessary to improve our knowledge not only of the atmospheric properties but also surface properties (temperature) and the surface-atmosphere interaction (volcanic activity).PFS works primarily around the pericentre of the orbit, only occasionally observing Venus from larger distances. Each measurements takes 4.5 s, with a repetition time of 11.5 s. By working roughly 1.5 h around pericentre, a total of 460 measurements per orbit will be acquired plus 60 for calibrations. PFS is able to take measurements at all local times, enabling the retrieval of atmospheric vertical temperature profiles on both the day and the night side.The PFS measures a host of atmospheric and surface phenomena on Venus. These include the:(1) thermal surface flux at several wavelengths near 1 μm, with concurrent constraints on surface temperature and emissivity (indicative of composition); (2) the abundances of several highly-diagnostic trace molecular species; (3) atmospheric temperatures from 55 to 100 km altitude; (4) cloud opacities and cloud-tracked winds in the lower-level cloud layers near 50-km altitudes; (5) cloud top pressures of the uppermost haze/cloud region near 70–80 km altitude; and (6) oxygen airglow near the 100 km level. All of these will be observed repeatedly during the 500-day nominal mission of Venus Express to yield an increased understanding of meteorological, dynamical, photochemical, and thermo-chemical processes in the Venus atmosphere. Additionally, PFS will search for and characterize current volcanic activity through spatial and temporal anomalies in both the surface thermal flux and the abundances of volcanic trace species in the lower atmosphere.Measurement of the 15 μm CO₂ band is very important. Its profile gives, by means of a complex temperature profile retrieval technique, the vertical pressure-temperature relation, basis of the global atmospheric study.PFS is made of four modules called O, E, P and S being, respectively, the interferometer and proximity electronics, the digital control unit, the power supply and the pointing device.     

The distribution of carbon monoxide in the lower atmosphere of Venus

We have obtained spatially resolved near-infrared spectroscopy of the Venus nightside on 15 nights over three observing seasons. We use the depth of the CO absorption band at 2.3 μm to map the two-dimensional distribution of CO across both hemispheres. Radiative transfer models are used to relate the measured CO band depth to the volume mixing ratio of CO. The results confirm previous investigations in showing a general trend of increased CO abundances at around 60° latitude north and south as compared with the equatorial regions. Observations taken over a few nights generally show very similar CO distributions, but significant changes are apparent over longer periods. In past studies it has been assumed that the CO latitudinal variation occurs near 35 km altitude, at which K-band sensitivity to CO is greatest. By modeling the detailed spectrum of the excess CO at high latitudes we show that it occurs at altitudes around 45 km, much higher than has previously been assumed, and that there cannot be significant contribution from levels of 36 km or lower. We suggest that this is most likely due to downwelling of CO-rich gas from the upper atmosphere at these latitudes, with the CO being removed by around 40 km through chemical processes such as the reaction with SO₃.     

A photochemical model for the Venus atmosphere at 47–112 km

The model is intended to respond to the recent findings in the Venus atmosphere from the Venus Express and ground-based submillimeter and infrared observations. It extends down to 47 km for comparison with the kinetic model for the lower atmosphere (Krasnopolsky, V.A. [2007]. Icarus 191, 25–37) and to use its results as the boundary conditions. The model numerical accuracy is significantly improved by reduction of the altitude step from 2 km in the previous models to 0.5 km. Effects of the NUV absorber are approximated using the detailed photometric observations at 365 nm from Venera 14. The H₂O profile is not fixed but calculated in the model. The model involves odd nitrogen and OCS chemistries based on the detected NO and OCS abundances. The number of the reactions is significantly reduced by removing of unimportant processes. Column rates for all reactions are given, and balances of production and loss may be analyzed in detail for each species.The calculated vertical profiles of CO, H₂O, HCl, SO₂, SO, OCS and of the O₂ dayglow at 1.27 μm generally agree with the existing observational data; some differences are briefly discussed. The OH dayglow is ～30 kR, brighter than the OH nightglow by a factor of 4. The H + O₃ process dominates in the nightglow excitation and O + HO₂ in the dayglow, because of the reduction of ozone by photolysis. A key feature of Venus’ photochemistry is the formation of sulfuric acid in a narrow layer near the cloud tops that greatly reduces abundances of SO₂ and H₂O above the clouds. Delivery of SO₂ and H₂O through this bottleneck determines the chemistry and its variations above the clouds. Small variations of eddy diffusion near 60 km result in variations of SO₂, SO, and OCS at and above 70 km within a factor of ～30. Variations of the SO₂/H₂O ratio at the lower boundary have similar but weaker effect: the variations within a factor of ～4 are induced by changes of SO₂/H₂O by ±5%. Therefore the observed variations of the mesospheric composition originate from minor variations of the atmospheric dynamics near the cloud layer and do not require volcanism. NO cycles are responsible for production of a quarter of O₂, SO₂, and Cl₂ in the atmosphere. A net effect of photochemistry in the middle atmosphere is the consumption of CO₂, SO₂, and HCl from and return of CO, H₂SO₄, and SO₂Cl₂ to the lower atmosphere. These processes may be balanced by thermochemistry in the lower atmosphere even without outgassing from the interior, though the latter is not ruled out by our models. Some differences between the model and observations and the previous models are briefly discussed.     

First ever in situ observations of Venus’ polar upper atmosphere density using the tracking data of the Venus Express Atmospheric Drag Experiment (VExADE)

On its highly elliptical 24 h orbit around Venus, the Venus Express (VEX) spacecraft briefly reaches a periapsis altitude of nominally 250 km. Recently, however, dedicated and intense radio tracking campaigns have taken place in August 2008, October 2009, February and April 2010, for which the periapsis altitude was lowered to the 186–176 km altitude range in order to be able to probe the upper atmosphere of Venus above the North Pole for the first time ever in situ. As the spacecraft experiences atmospheric drag, its trajectory is measurably perturbed during the periapsis pass, allowing us to infer total atmospheric mass density at the periapsis altitude. A Precise Orbit Determination (POD) of the VEX motion is performed through an iterative least-squares fitting process to the Doppler tracking data, acquired by the VEX radioscience experiment (VeRa). The drag acceleration is modelled using an initial atmospheric density model (VTS3 model, Hedin, A.E., Niemann, H.B., Kasprzak, W.T., Seiff, A. [1983]. J. Geophys. Res. 88, 73–83). A scale factor of the drag acceleration is estimated for each periapsis pass, which scales Hedin’s density model in order to best fit the radio tracking data. Reliable density scale factors have been obtained for 10 passes mainly from the second (October 2009) and third (April 2010) VExADE campaigns, which indicate a lower density by a factor of about 1.8 than Hedin’s model predicts. These first ever in situ polar density measurements at solar minimum have allowed us to construct a diffusive equilibrium density model for Venus’ thermosphere, constrained in the lower thermosphere primarily by SPICAV-SOIR measurements and above 175 km by the VExADE drag measurements (Müller-Wodarg et al., in preparation). The preliminary results of the VExADE campaigns show that it is possible to obtain with the POD technique reliable estimates of Venus’ upper atmosphere densities at an altitude of around 175 km. Future VExADE campaigns will benefit from the planned further lowering of VEX pericenter altitude to below 170 km.     

Small-scale temperature fluctuations seen by the VeRa Radio Science Experiment on Venus Express

The Venus Express Radio Science Experiment VeRa retrieves atmospheric profiles in the mesosphere and troposphere of Venus in the approximate altitude range of 40–90 km. A data set of more than 500 profiles was retrieved between the orbit insertion of Venus Express in 2006 and the end of occultation season No. 11 in July 2011. The atmospheric profiles cover a wide range of latitudes and local times, enabling us to study the dependence of vertical small-scale temperature perturbations on local time and latitude.Temperature fluctuations with vertical wavelengths of 4 km or less are extracted from the measured temperature profiles in order to study small-scale gravity waves. Significant wave amplitudes are found in the stable atmosphere above the tropopause at roughly 60 km as compared with the only shallow temperature perturbations in the nearly adiabatic region of the adjacent middle cloud layer, below.Gravity wave activity shows a strong latitudinal dependence with the smallest wave amplitudes located in the low-latitude range, and an increase of wave activity with increasing latitude in both hemispheres; the greatest wave activity is found in the high-northern latitude range in the vicinity of Ishtar Terra, the highest topographical feature on Venus.We find evidence for a local time dependence of gravity wave activity in the low latitude range within ±30° of the equator. Gravity wave amplitudes are at their maximum beginning at noon and continuing into the early afternoon, indicating that convection in the lower atmosphere is a possible wave source.The comparison of the measured vertical wave structures with standard linear-wave theory allows us to derive rough estimates of the wave intrinsic frequency and horizontal wavelengths, assuming that the observed wave structures are the result of pure internal gravity waves. Horizontal wavelengths of the waves at 65 km altitude are on the order of ≈300–450 km with horizontal phase speeds of roughly 5–10 m/s.     

The Huygens Probe Descent Trajectory Working Group: Organizational framework,goals, and implementation

Cassini/Huygens, a flagship mission to explore the rings, atmosphere, magnetic field, and moons that make up the Saturn system, is a joint endeavor of the National Aeronautics and Space Administration, the European Space Agency, and Agenzia Spaziale Italiana. Comprising two spacecraft—a Saturn orbiter built by NASA and a Titan entry/descent probe built by the European Space Agency—Cassini/Huygens was launched in October 1997. The Huygens probe parachuted to the surface of Titan in January 2005. During the descent, six science instruments provided in situ measurements of Titan's atmosphere, clouds, and winds, and photographed Titan's surface. To correctly interpret and correlate results from the probe science experiments, and to provide a reference set of data for ground-truth calibration of orbiter remote sensing measurements, an accurate reconstruction of the probe entry and descent trajectory and surface landing location is necessary. The Huygens Descent Trajectory Working Group was chartered in 1996 as a subgroup of the Huygens Science Working Team to develop and implement an organizational framework and retrieval methodologies for the probe descent trajectory reconstruction from the entry altitude of 1270 km to the surface using navigation data, and engineering and science data acquired by the instruments on the Huygens Probe. This paper presents an overview of the Descent Trajectory Working Group, including the history, rationale, goals and objectives, organizational framework, rules and procedures, and implementation.     

Reconstruction of the trajectory of the Huygens probe using the Huygens Atmospheric Structure Instrument (HASI)

The Huygens probe returned scientific measurements from the atmosphere and surface of Titan on 14 January 2005. Knowledge of the trajectory of Huygens is necessary for scientific analysis of those measurements. We use measurements from the Huygens Atmospheric Structure Instrument (HASI) to reconstruct the trajectory of Huygens during its mission. The HASI Accelerometer subsystem measured the axial acceleration of the probe with errors of 3E−6 m s⁻², the most accurate measurements ever made by an atmospheric structure instrument on another planetary body. The atmosphere was detected at an altitude of 1498 km. Measurements of the normal acceleration of the probe, which are important for determining the probe's attitude during hypersonic entry, were significantly less accurate and limited by transverse sensitivity of the piezo sensors. Peak acceleration of 121.2 m s⁻² occurred at 234.9 km altitude. The parachute deployment sequence started at 157.1 km and a speed of 342.1 m s⁻¹. Direct measurements of pressure and temperature began shortly afterwards. The measured accelerations and equations of motion have been used to reconstruct the trajectory prior to parachute deployment. Measured pressures and temperatures, together with the equation of hydrostatic equilibrium and the equation of state, have been used to reconstruct the trajectory after parachute deployment. Uncertainties in the entry state of Huygens at the top of the atmosphere are significant, but can be reduced by requiring that the trajectory and atmospheric properties be continuous at parachute deployment.     

The O2 nightglow in the martian atmosphere by SPICAM onboard of Mars-Express

We present observations of the O₂(a¹Δ_g) nightglow at 1.27 μm on Mars using the SPICAM IR spectrometer onboard of the Mars Express orbiter. In contrast to the O₂(a¹Δ_g) dayglow that results from the ozone photodissociation, the O₂(a¹Δ_g) nightglow is a product of the recombination of O atoms formed by CO₂ photolysis on the dayside at altitudes higher than 80 km and transported downward above the winter pole by the Hadley circulation. The first detections of the O₂(a¹Δ_g) nightglow in 2010 indicate that it is about two order of magnitude less intense than the dayglow (Bertaux, J.-L., Gondet, B., Bibring, J.-P., Montmessin, F., Lefèvre, F. [2010]. Bull. Am. Astron. Soc. 42, 1040; Clancy et al. [2010]. Bull. Am. Astron. Soc. 42, 1041). SPICAM IR sounds the martian atmosphere in the near-IR range (1–1.7 μm) with the spectral resolution of 3.5 cm^?1 in nadir, limb and solar occultation modes. In 2010 the vertical profiles of the O₂(a¹Δ_g) nightside emission have been obtained near the South Pole at latitudes of 82–83°S for two sequences of observations: L_s = 111–120° and L_s = 152–165°. The altitude of the emission maximum varied from 45 km on L_s = 111–120° to 38–49 km on L_s = 152–165°. Averaged vertically integrated intensity of the emission at these latitudes has shown an increase from 0.22 to 0.35 MR. Those values of total vertical emission rate are consistent with the OMEGA observations on Mars-Express in 2010. The estimated density of oxygen atoms at altitudes from 50 to 65 km varies from 1.5 × 10¹¹ to 2.5 × 10¹¹ cm^?3. Comparison with the LMD general circulation model with photochemistry (Lefèvre, F., Lebonnois, S., Montmessin, F., Forget, F. [2004]. J. Geophys. Res. 109, E07004; Lefèvre et al. [2008]. Nature 454, 971–975) shows that the model reproduces fairly well the O₂(a¹Δ_g) emission layer observed by SPICAM when the large field of view (>20 km on the limb) of the instrument is taken into account.