A re-analysis of the 1971 Beta Scorpii occultation by Jupiter: study of temperature fluctuations and detection of wave activity

Data from the 13 May 1971 β Scorpii occultation by the southern polar region of Jupiter (Vapillon et al., 1973, Astron. Astrophys. 29, 135-149) are re-analyzed with current methods. We correct the previous results for an inacurrate background estimation and calculate new temperature profiles, that are now consistent with the results of other observers of this occultation, as well as with the current knowledge of the jovian atmosphere. The characteristics of the profiles of temperature gradient and the spectral behavior of the temperature fluctuations are found to be similar to the results of previous investigations of planetary atmospheres and in agreement with the presence of atmospheric propagating gravity waves in the jovian atmosphere. We use a wavelet analysis of the temperature profiles to identify the dominant modes of wave activity and compare the reconstructed temperature fluctuations to model-generated gravity waves.     

The structure of Titan’s atmosphere from Cassini radio occultations

We present results from the two radio occultations of the Cassini spacecraft by Titan in 2006, which probed mid-southern latitudes. Three of the ingress and egress soundings occurred within a narrow latitude range, 31-34°S near the surface, and the fourth at 52.8°S. Temperature-altitude profiles for all four occultation soundings are presented, and compared with the results of the Voyager 1 radio occultation (Lindal, G.F., Wood, G.E., Hotz, H.B., Sweetnam, D.N., Eshleman, V.R., Tyler, G.L. [1983]. Icarus 53, 348-363), the HASI instrument on the Huygens descent probe (Fulchignoni, M. et al. [2005]. Nature 438, 785-791), and Cassini CIRS results (Flasar, F.M. et al. [2005]. Science 308, 975-978; Achterberg, R.K., Conrath, B.J., Gierasch, P.J., Flasar, F.M., Nixon, C.A. [2008b]. Icarus 194, 263-277). Sources of error in the retrieved temperature-altitude profiles are also discussed, and a major contribution is from spacecraft velocity errors in the reconstructed ephemeris. These can be reduced by using CIRS data at 300 km to make along-track adjustments of the spacecraft timing. The occultation soundings indicate that the temperatures just above the surface at 31-34°S are about 93 K, while that at 53°S is about 1 K colder. At the tropopause, the temperatures at the lower latitudes are all about 70 K, while the 53°S profile is again 1 K colder. The temperature lapse rate in the lowest 2 km for the two ingress (dawn) profiles at 31 and 33°S lie along a dry adiabat except within ∼200 m of the surface, where a small stable inversion occurs. This could be explained by turbulent mixing with low viscosity near the surface. The egress profile near 34°S shows a more complex structure in the lowest 2 km, while the egress profile at 53°S is more stable.     

A reanalysis of Venus winds at two cloud levels from Galileo SSI images

The global circulation of the Venus atmosphere is characterized at cloud level by a zonal super rotation studied over the years with data from a battery of spacecrafts: orbiters, balloons and probes. Among them, the Galileo spacecraft monitored the Venus atmosphere in a flyby in February 1990 in its route toward Jupiter. Since the flyby was almost equatorial, published analysis of zonal winds obtained from displacements of cloud elements on images obtained by the SSI camera [Belton, M.J.S., and 20 colleagues, 1991. Science 253, 1531-1536] stop at latitudes 50° north and south. In this paper we present new results on Venus winds based on a reanalysis of an extended set of images obtained at two wavelengths, 418 nm (violet) and 986 nm (near infrared), that sense different altitude levels in the upper cloud. Our main result is that we have been able to extend the zonal wind profile up to the polar latitudes: 70° N and 70° S at 418 nm and 70° N at 986 nm. Binned and smoothed profiles are given in tabular form. We show that the zonal winds drop in their velocity poleward of latitudes 45° N and 50° S where an intense meridional wind shear develops at the two cloud levels. Our data confirm the magnitude of this shear, retrieved previously from radio occultation data, but disagrees with it in the latitudinal location of the sheared region. The new wind data can be used to recalibrate the zonal winds retrieved from the previous measurements of the temperature field and the cyclostrophic balance assumption. The meridional profiles of the zonal winds at the two cloud levels are used to assess the vertical wind shear in the upper cloud layer as a function of latitude and locate the most unstable region.     

Radio occultation measurements and MGCM simulations of Kelvin waves on Mars

We have derived new results concerning thermal tides on Mars from a combination of radio occultation measurements and numerical simulations by a Mars General Circulation Model (MGCM). This investigation exploits a set of concurrent observations by Mars Express (MEX) and Mars Global Surveyor (MGS) in mid-2004, when the season on Mars was midspring in the northern hemisphere. The MEX occultations sampled the atmosphere near the evening terminator at latitudes ranging from 54° N to 15° S. The MGS occultations provided complementary coverage near the morning terminator at latitudes of 35° N and 71° S. The geopotential field derived from these measurements contains distinctive modulation caused by solar-asynchronous thermal tides. Through careful analysis of the combined observations, we characterized two prominent wave modes, obtaining direct solutions for some properties, such as the amplitude and phase, as well as constraints on others, such as the period, zonal wave number, and meridional structure. We supplemented these observations with MGCM simulations. After evaluating the performance of the MGCM against the measurements, we used the validated simulation to deduce the identity of the two tidal modes and to explore their behavior. One mode is a semidiurnal Kelvin wave with a zonal wave number of 2 (SK2), while the other is a diurnal Kelvin wave with a zonal wave number of 1 (DK1). Both modes are known to be close to resonance in the martian atmosphere. Our observations of the SK2 are more complete and less ambiguous than any previous measurement. The well-known DK1 is the dominant solar-asynchronous tide in the martian atmosphere, and our results confirm and extend previous observations by diverse instruments.     

Strong turbulence and atmospheric waves in stellar occultations

Many of the problems of stellar occultation observations stem from the difficulty of determining the effects of realistic atmospheric structure on the lightcurves. General techniques for producing model lightcurves for a variety of realistic atmospheric irregularities, including turbulence and inertia-gravity waves, are presented and applied. Using numerical simulations which model the propagation of a wave through a phase-changing screen, the limit of strong scintillations for one-dimensional, Kolmogorov-like turbulence, both for a point source and for extended sources, is investigated in some detail, and significant departures from the behavior in the weak scintillation regime are found. The results are compared with published analytical results and recent occultation data. The effects of large-scale atmospheric waves with realistic horizontal structure are examined, and the reliability of the numerical inversion method of retrieving the true atmospheric vertical structure under circumstances of strong ray crossing and horizontal inhomogeneities is assessed. The simulations confirm that large-scale layered features of the atmosphere are accurately recovered; horizontally inhomogeneous structures (including turbulence) with coherence scale $\text{L ? (2πRH)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ (where R = planetary radius and H = scale height) have little effect on the derived temperature profiles. It is concluded that analysis of occultations may eventually allow us to determine both the quasiglobal atmospheric structure and the statistical characteristics of small-scale refractivity variations.     

Propagation of tropospheric gravity waves into the upper atmosphere of Mars

We study the propagation of gravity waves in the martian atmosphere using a linearized one-dimensional full-wave model. Calculations are carried out for atmospheric parameters characteristic of Mars Orbiter Laser Altimeter (on Mars Global Surveyor MGS) observations of apparent gravity waves in high latitude clouds and MGS radio occultation measurements of temperature variations with height suggestive of gravity wave activity. Waves that reach the thermosphere produce fluctuations in density comparable in amplitude with the density variations detected in Mars Odyssey aerobraking data. Gravity waves of modest amplitude are found to deposit momentum and generate significant heating and cooling in the martian atmosphere. The largest heating and cooling effects occur in the thermosphere, at altitudes between about 130 and 150 km, with heating occurring at the lower altitudes and cooling taking place above.     

Methane on Uranus: The case for a compact CH4 cloud layer at low latitudes and a severe CH4 depletion at high-latitudes based on re-analysis of Voyager occultation measurements and STIS spectroscopy

Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92 (11), 14987-15001) presented a range of temperature and methane profiles for Uranus that were consistent with 1986 Voyager radio occultation measurements of refractivity versus altitude. A localized refractivity slope variation near 1.2 bars was interpreted to be the result of a condensed methane cloud layer. However, models fit to near-IR spectra found particle concentrations much deeper in the atmosphere, in the 1.5-3 bar range (Sromovsky, L.A., Irwin, P.G.J., Fry, P.M. [2006]. Icarus 182, 577-593; Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229; Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2010]. Icarus 208, 913-926), and a recent analysis of STIS spectra argued for a model in which aerosol particles formed diffusely distributed hazes, with no compact condensation layer (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309). To try to reconcile these results, we reanalyzed the occultation observations with the He volume mixing ratio reduced from 0.15 to 0.116, which is near the edge of the 0.033 uncertainty range given by Conrath et al. (Conrath, B., Hanel, R., Gautier, D., Marten, A., Lindal, G. [1987]. J. Geophys. Res. 92 (11), 15003-15010). This allowed us to obtain saturated mixing ratios within the putative cloud layer and to reach above-cloud and deep methane mixing ratios compatible with STIS spectral constraints. Using a 5-layer vertical aerosol model with two compact cloud layers in the 1-3 bar region, we find that the best fit pressure for the upper compact layer is virtually identical to the pressure range inferred from the occultation analysis for a methane mixing ratio near 4% at 5°S. This strongly argues that Uranus does indeed have a compact methane cloud layer. In addition, our cloud model can fit the latitudinal variations in spectra between 30°S and 20°N, using the same profiles of temperature and methane mixing ratio. But closer to the pole, the model fails to provide accurate fits without introducing an increasingly strong upper tropospheric depletion of methane at increased latitudes, in rough agreement with the trend identified by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309).     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

Impact of atmospheric gravity waves on the jovian ionosphere

We investigate the effects of atmospheric gravity waves on the vertical and horizontal structure of the ionosphere of Jupiter. The presented non-linear, two-dimensional model of the jovian ionosphere allows for spatially and temporally varying neutral wind and temperature fields and tracks the time evolution of six ionospheric species, , and . An analytical approach is used to validate the model results for linear, small-amplitude waves and to elucidate the mechanisms that leads to perturbations in the density of the main ion species, H⁺ and . We demonstrate that the long-lived H⁺ ions are perturbed directly by wave dynamics whereas short-lived ions such as are perturbed by chemical interactions with other perturbed ion species. The model is then applied using larger gravity wave amplitudes consistent with observations. Atmospheric gravity waves propagating at high altitudes create layers of enhanced electron density similar to the system of layers observed during the J0-ingress radio occultation of the Galileo spacecraft. Our best fit to the J0-ingress observation is achieved using an 82 min period forcing wave with horizontal and vertical wavelengths of 500 km and 60 km respectively, and peaks at 510 km above the 1 bar pressure level. We further investigate the effects of the wave-induced ion flux on the background ionospheric structure and demonstrate that in the presence of a gravity wave the background density profiles of the H⁺ and ions are significantly modified. We also find that the column density of has variations that can exceed 10% as the wave propagates.     

Estimates of gravity wave drag on Mars: Indication of a possible lower thermospheric wind reversal

Spectral gravity wave parameterization suitable for planetary thermospheres applied to wind and temperature from Mars Climate Database reveals enormously strong drag incompatible with the wind distribution. It points out to a possible wind reversal in the 110-140 km layer similar to the one in the Earth’s lower thermosphere.     

Numerical simulation of atmospheric bore waves on Mars

The Viking Orbiters imaged early morning, long, linear wave clouds along the flanks of the Tharsis volcanoes during late northern spring and early summer. These clouds are believed to be a product of either an atmospheric bore wave or a hydraulic jump generated by nightly katabatic winds. The Mars Regional Atmospheric Modeling System was used to study the interaction of the katabatic flows with the surrounding atmosphere to determine what mechanism is responsible for the clouds. Simulations at Ls=90°, 100°, 142°, 180°, 270°, and 358° were conducted focusing on the eastern flank of Olympus Mons. Model results compare well with Viking observations and closely approximate theoretical treatments of atmospheric bores. Strong downslope flows are simulated during the night, with a bore wave forming on and behind a well-defined katabatic front. The observed seasonality of the clouds was reproduced in the simulations; the bore was deeper and faster during northern summer and weakest during the winter. When the bore was strong, it was undular in form, and generated vertically propagating gravity waves in the atmosphere above. During the winter, the atmospheric structure was such that any gravity waves generated damped with height. Less atmospheric water vapor abundance during northern winter, as compared to the summer, is also a factor in the seasonality of the wave clouds. This study concludes that bore waves are the most likely mechanism for the generation of the observed linear wave clouds.     

Radar and optical observations and physical modeling of near-Earth Asteroid 10115 (1992 SK)

We estimate Asteroid 1992 SK's physical properties from delay-Doppler images and Doppler-only echo spectra obtained during March 22-27, 1999, at Goldstone and from optical lightcurves obtained during February-March 1999 at Ond?ejov Observatory. The images span only about 15° of sky motion and are not strong, but they place up to twenty 40 m by 160 m pixels on the asteroid and have complete rotational phase coverage. Our analysis establishes that the radar observations are confined to subradar latitudes between −20° and −40°. The echo spectra and optical lightcurves span ∼80° of sky motion, which provides important geometric leverage on the pole direction. The lightcurves are essential for accurate estimation of the asteroid's shape and spin state. We estimate the asteroid's period to be 7.3182±0.0003 h and its pole direction to be at ecliptic longitude, latitude=(99°±5°,−3°±5°). The asteroid is about 1.4 km in maximum extent and mildly asymmetric, with an elongation of about 1.5 and relatively subdued topography. The OC radar albedo is 0.11±0.02 and the SC/OC ratio is 0.34±0.05. The current orbital solution permits accurate identification of planetary close approaches during 826-2690. We use our model to predict salient characteristics of radar images and optical lightcurves obtainable during the asteroid's March 2006 approach.     

The atmospheric temperatures over Olympus Mons on Mars: An atmospheric hot ring

We study the thermal fields over Olympus Mons separating seasons (northern spring and summer against southern spring and summer) and local time observations (day side versus night side). Temperature vertical profiles retrieved from Planetary Fourier Spectrometer on board Mars Express (PFS-MEX) data have been used. In many orbits (running north to south along a meridian) passing over the top of the volcano there is evidence of a hot air on top of the volcano, of two enhancement of the air temperature both north and south of it and in between a collar of air that is colder than nearby at low altitudes, and warmer than nearby at high altitudes. Mapping together many orbits passing over or near the volcano we find that the hot air has the tendency to form an hot ring around it. This hot structure occurs mostly between LT = 10.00 and 15.00 and during the northern summer. Distance of the hot structure from the top of the volcano is about 600 km (10° of latitude). The hot atmospheric region is 300-420 km (5-7°) wide. Hot ring temperature contrasts of about 40 K occur at 2 km above the surface and decrease to 20 K at 5 km and to 10 K at 10 km. The atmospheric circulation over an area of 40° × 40° (latitudes and longitudes) is affected by the topography and the presence of Olympus Mons (−133°W, 18°N). We discuss also the thermal stability of the atmosphere over the selected area using the potential temperatures. The temperature field over the top of the volcano shows unstable atmosphere within 10 km from the surface. Finally, we interpret the hot temperatures around volcano as an adiabatic compression of down-welling branch coming from over the top of volcano.Different air temperature profiles are observed in the same seasons during the night, or in different seasons. In northern spring-summer during the night the isothermal contours do not show the presence of the volcano until we reach close to the surface very much, where a thermal inversion is observed. The surface temperature shows higher values (by 10 K) in correspondence of the scarp (an abrupt altimetry variation of roughly 5 km) on south (6°N) and north (30°N) sides of volcano. During the southern spring-summer, on the contrary the isothermal curves run parallel to the surface even on top the volcano, just like the GCM have predicted.     

The 2003 November 14 occultation by Titan of TYC 1343-1865-1: II. Analysis of light curves

We observed a stellar occultation by Titan on 2003 November 14 from La Palma Observatory using ULTRACAM with three Sloan filters: u^′, g^′, and i^′ (358, 487, and 758 nm, respectively). The occultation probed latitudes 2° S and 1° N during immersion and emersion, respectively. A prominent central flash was present in only the i^′ filter, indicating wavelength-dependent atmospheric extinction. We inverted the light curves to obtain six lower-limit temperature profiles between 335 and 485 km (0.04 and 0.003 mb) altitude. The i^′ profiles agreed with the temperature measured by the Huygens Atmospheric Structure Instrument [Fulchignoni, M., and 43 colleagues, 2005. Nature 438, 785–791] above 415 km (0.01 mb). The profiles obtained from different wavelength filters systematically diverge as altitude decreases, which implies significant extinction in the light curves. Applying an extinction model [Elliot, J.L., Young, L.A., 1992. Astron. J. 103, 991–1015] gave the altitudes of line of sight optical depth equal to unity: 396±7 and 401±20 km (u^′ immersion and emersion); 354±7 and 387±7 km (g^′ immersion and emersion); and 336±5 and 318±4 km (i^′ immersion and emersion). Further analysis showed that the optical depth follows a power law in wavelength with index 1.3±0.2. We present a new method for determining temperature from scintillation spikes in the occulting body's atmosphere. Temperatures derived with this method are equal to or warmer than those measured by the Huygens Atmospheric Structure Instrument. Using the highly structured, three-peaked central flash, we confirmed the shape of Titan's middle atmosphere using a model originally derived for a previous Titan occultation [Hubbard, W.B., and 45 colleagues, 1993. Astron. Astrophys. 269, 541–563].     

Titan's middle-atmospheric temperatures and dynamics observed by the Cassini Composite Infrared Spectrometer

The Composite Infrared Radiometer-Spectrometer (CIRS) instrument, on the NASA Cassini Saturn orbiter, has been acquiring thermal emission spectra from the atmosphere of Titan since orbit insertion in 2004. Observation sequences for measuring stratospheric temperatures have been obtained using both a nadir mapping mode and a limb viewing mode. The limb observations give better vertical resolution, and give information from higher altitudes, while the nadir observations provide more complete longitude coverage. Because the scale height of Titan's atmosphere is large enough so that emission from a grazing ray is influenced by horizontal temperature variations in the atmosphere, we have developed a two-dimensional temperature retrieval algorithm for reducing the limb spectra, which solves simultaneously for meridional and vertical temperature variations. The analyzed nadir mapping data have sampled nearly all longitudes at latitudes from about 90° S to 60° N, providing temperatures between pressure levels of about 5 to 0.2 mbar. The limb data covers latitudes between about 75° S and 85° N, and yields temperatures between about 1 and 0.005 mbar, at a small number of longitudes. The retrieved temperatures are consistent with early results from nadir observations [Flasar, F.M., and 44 colleagues, 2005. Science 308, 975-978] between 0.5 and 5 mbar where both results are valid, with the warmest temperatures at the equator, and much stronger meridional temperature gradients in the northern (winter) hemisphere than in the southern. At higher altitudes not probed by nadir viewing, the limb data reveal that the stratopause is nearly 20 K warmer in the northern polar regions than at the equator and southern hemisphere, and that the altitude of the stratopause shifts from ≈0.1 mbar (300 km) near the equator to 0.01 mbar (400 km) poleward of about 40° N. When the gradient wind equation is used to construct a zonal mean wind, the reversal in sign of the temperature leads to capping of the winter westerly flow. The core of the resulting jet is about 190 m s⁻¹ in magnitude, spans between 30° N and 60° N, and peaks near 0.1 mbar. Estimates of the radiative heating associated with the radiative disequilibrium lead to a meridional overturning timescale of about three Earth years.     

On the long-term variability of Jupiter’s winds and brightness as observed from Hubble

Hubble Space Telescope Wide Field Planetary Camera 2 imaging data of Jupiter were combined with wind profiles from Voyager and Cassini data to study long-term variability in Jupiter’s winds and cloud brightness. Searches for evidence of wind velocity periodicity yielded a few latitudes with potential variability; the most significant periods were found nearly symmetrically about the equator at 0°, 10-12°N, and 14-18°S planetographic latitude. The low to mid-latitude signals have components consistent with the measured stratospheric temperature Quasi-Quadrennial Oscillation (QQO) period of 4-5 years, while the equatorial signal is approximately seasonal and could be tied to mesoscale wave formation. Robustness tests indicate that a constant or continuously varying periodic signal near 4.5 years would appear with high significance in the data periodograms as long as uncertainties or noise in the data are not of greater magnitude. However, the lack of a consistent signal over many latitudes makes it difficult to interpret as a QQO-related change. In addition, further analyses of calibrated 410-nm and 953-nm brightness scans found few corresponding changes in troposphere haze and cloud structure on QQO timescales. However, stratospheric haze reflectance at 255-nm did appear to vary on seasonal timescales, though the data do not have enough temporal coverage or photometric accuracy to be conclusive. Sufficient temporal coverage and spacing, as well as data quality, are critical to this type of search.     

The depth of the convective boundary layer on Mars

We are using observations obtained with Mars Express to explore the structure and dynamics of the martian lower atmosphere. We consider a series of radio occultation experiments conducted in May-August 2004, when the season on Mars was midspring of the northern hemisphere. The measurements are widely distributed in latitude and longitude, but the local time remained within a narrow range, 17.0-17.2 h. Most of the atmospheric profiles retrieved from these data contain a distinct, well-mixed convective boundary layer (CBL). We have accurately determined the depth of the CBL and its spatial variations at fixed local time through analysis of these profiles. The CBL extends to a height of 3-10 km above the surface at the season and locations of these measurements. Its depth at fixed local time is clearly correlated with variations in surface elevation on planetary scales, with a weaker dependence on spatial variations in surface temperature. In general, the CBL is deep (8-10 km) where the surface elevation is high, as in Tharsis Montes and Syrtis Major, and shallow (4-6 km) where the surface elevation is low, as in Amazonis and Utopia. This variability results from the combined effects of conditions near the surface and in the atmosphere above the CBL. Convection arises from solar heating of the ground, and the impact of this heat source on thermal structure is largest where the surface pressure and atmospheric density are smallest, at high surface elevations. The vertical extent of the CBL is in turn constrained by the static stability of the overlying atmosphere. These results greatly reduce the long-standing uncertainty concerning the depth of the CBL.     

Validation of martian meteorological data assimilation for MGS/TES using radio occultation measurements

We describe an assimilation of thermal profiles below about 40 km altitude and total dust opacities into a general circulation model (GCM) of the martian atmosphere. The data were provided by the Thermal Emission Spectrometer (TES) on board the Mars Global Surveyor (MGS) spacecraft. The results of the assimilation are verified against an independent source of contemporaneous data represented by radio occultation measurements with an ultra-stable radio oscillator, also aboard MGS. This paper describes a comparison between temperature profiles retrieved by the radio occultation experiments and the corresponding profiles given by both an independent, carefully tuned GCM simulation and by an assimilation of TES observations performed over the period of time from middle, northern summer in martian year 24, corresponding to May 1999, until late, northern spring in martian year 27, corresponding to August 2004. This study shows that the assimilation of TES measurements improves the overall agreement between radio occultation observations and the GCM analysis, in particular below 20 km altitude, where the radio occultation measurements are known to be most accurate. Discrepancies still remain, mostly during the global dust storm of year 2001 and at latitudes around 60° N in northern winter-early spring. These are the periods of time and locations, however, for which discrepancies between TES and radio occultation profiles are also shown to be the largest. Finally, a further direct validation is performed, comparing stationary waves at selected latitudes and time of year. Apart from biases at high latitudes in winter time, data assimilation is able to represent the correct wave behaviour, which is one major objective for martian assimilation.     

Vertical profiles of HCN, HC3N, and C2H2 in Titan's atmosphere derived from Cassini/CIRS data

Mid-infrared limb spectra in the range 600-1400 cm⁻¹ taken with the Composite InfraRed Spectrometer (CIRS) on-board the Cassini spacecraft were used to determine vertical profiles of HCN, HC₃N, C₂H₂, and temperature in Titan's atmosphere. Both high (0.5 cm⁻¹) and low (13.5 cm⁻¹) spectral resolution data were used. The 0.5 cm⁻¹ data gave profiles at four latitudes and the 13.5 cm⁻¹ data gave almost complete latitudinal coverage of the atmosphere. Both datasets were found to be consistent with each other. High temperatures in the upper stratosphere and mesosphere were observed at Titan's northern winter pole and were attributed to adiabatic heating in the subsiding branch of a meridional circulation cell. On the other hand, the lower stratosphere was much colder in the north than at the equator, which can be explained by the lack of solar radiation and increased IR emission from volatile enriched air. HC₃N had a vertical profile consistent with previous ground based observations at southern and equatorial latitudes, but was massively enriched near the north pole. This can also be explained in terms of subsidence at the winter pole. A boundary observed at 60° N between enriched and un-enriched air is consistent with a confining polar vortex at 60° N and HC₃N's short lifetime. In the far north, layers were observed in the HC₃N profile that were reminiscent of haze layers observed by Cassini's imaging cameras. HCN was also enriched over the north pole, which gives further evidence for subsidence. However, the atmospheric cross section obtained from 13.5 cm⁻¹ data indicated a HCN enriched layer at 200-250 km, extending into the southern hemisphere. This could be interpreted as advection of polar enriched air towards the south by a meridional circulation cell. This is observed for HCN but not for HC₃N due to HCN's longer photochemical lifetime. C₂H₂ appears to have a uniform abundance with altitude and is not significantly enriched in the north. This is consistent with observations from previous CIRS analysis that show increased abundances of nitriles and hydrocarbons but not C₂H₂ towards the north pole.     

Analysis of stellar occultation data: Effects of photon noise and initial conditions

A new inversion technique for obtaining temperature, pressure, and number density profiles of a planetary atmosphere from an occultation light curve is described. This technique employs an improved boundary condition to begin the numerical inversion and permits the computation of errors in the profiles caused by photon noise in the light curve. We present our assumptions about the atmosphere, optics, and noise and develop the equations for temperature, pressure, and number density and their associated errors. By inverting in equal increments of altitude, Δh, rather than in equal increments of time, Δt, the inversion need not be halted at the first negative point on the light curve as required by previous methods. The importance of the boundary condition is stressed, and a new initial condition is given. Numerical results are presented for the special case of inversion of a noisy isothermal light curve. From these results, simple relations are developed which can be used to predict the noise quality of an occultation. It is found that fractional errors in temperature profiles are comparable to those of pressure and number density profiles. An example of the inversion method is shown. Finally, we discuss the validity of our assumptions. In an appendix we demonstrate that minimum fractional errors in scale height determined from the inversion are comparable to those from an isothermal fit to a noisy isothermal light curve.