GCM simulation of balloon trajectories on Titan

Possible trajectories of passive balloons in Titan's troposphere are simulated with the instantaneous wind field predicted by a GCM (general circulation model). In most areas the basic motion of a balloon is a predominantly eastward or westward drift, depending on altitude, latitude and season of the balloon release point. Some meridional oscillation is always superposed on this basic motion, resulting in a wavy trajectory, with a maximum extent (of 40^°) at high latitudes of the winter hemisphere. As a general rule, the meridional oscillation can be maximised if a balloon is deployed in altitudes and latitudes where the mean zonal wind is eastward and smaller than the phase speed of Saturn's gravitational tide on Titan. A balloon's groundtrack does not repeat as it makes successive circuits around Titan, but rather makes a spiral or braided pattern. The summer pole is rather difficult to access for a balloon not directly introduced there because of small meridional oscillation, while the winter pole can be readily accessed and left several times. A preferred zonal drift direction can be achieved by choosing a proper altitude and hemisphere, but choosing a preferred meridional direction is not possible.     

Near-surface winds at the Huygens site on Titan: Interpretation by means of a general circulation model

This study aims at interpreting the zonal and meridional wind in Titan's troposphere measured by the Huygens probe by means of a general circulation model. The numerical simulation elucidates the relative importance of the seasonal variation in the Hadley circulation and Saturn's gravitational tide in affecting the actual wind profile. The observed reversal of the zonal wind at two altitudes in the lower troposphere can be reproduced with this model only if the near-surface temperature profile is asymmetric about the equator and substantial seasonal redistribution of angular momentum by the variable Hadley circulation takes place. The meridional wind near the surface is mainly caused by the meridional pressure gradient and is thus a manifestation of the Hadley circulation. Southward meridional wind in the PBL (planetary boundary layer) is consistent with the near-surface temperature at the equator being lower than at mid southern latitudes. Even small changes in the radiative heating profile in the troposphere can substantially affect the mean zonal and meridional wind including their direction. Saturn's gravitational tide is rather weak at the Huygens site due to the proximity to the equator, and does not clearly manifest itself in the instantaneous vertical profile of wind. Nevertheless, the simulated descent trajectory is more consistent with the observation if the tide is present. Because of a different force balance in Titan's atmosphere from terrestrial conditions, PBL-specific wind systems like on Earth are unlikely to exist on Titan.     

A tidal explanation for the Titan haze layers

Though Titan is in synchronous rotation around Saturn, it experiences gravitational tides as a consequence of its eccentric orbit. It is proposed that the vertical transport of aerosols by these tides produces the haze layers in Titan's upper atmosphere. Analysis shows that the zonal winds in Titan's superrotating atmosphere have a profound influence on which tidal components are effective in establishing the multiple detached-haze layers. If the Huygens Doppler winds are representative of the equatorial global superrotation, then the westward propagating s=2 mode is the responsible tidal component even though its forcing is significantly weaker than that of the s=0 and eastward s=2 components. The eastward s=2 tidal mode is eliminated by critical levels while the s=0 mode is viscously damped in the strong high altitude winds. At polar latitudes, however, the gravest s=0 mode is the one most likely to produce layering. It is also suggested that the atmospheric gravitational tides could be responsible for decelerating the superrotating atmosphere as seen in the Huygens Doppler wind velocity profile at about 80 km altitude.     

Thermal structure of Saturn from Voyager infrared measurements: Implications for atmospheric dynamics

Saturn atmospheric temperatures at the 150-mbar level retrieved from Voyager IRIS measurements indicate the presence of small-scale meridional gradients which are approximately symmetric with respect to the equator, but are superposed on a large-scale hemispheric thermal asymmetry. Under the assumption that the retrieved values at this atmospheric level represent kinetic temperatures on a constant pressure surface, it is suggested that the small-scale structure is produced by a meriodional circulation associated with the dissipative decay of the zonal winds with height, while the hemispheric asymmetry represents a thermal response to the seasonally varying insolation. The small-scale gradients are correlated with zonal winds derived from Voyager images at mid and high latitudes through the thermal wind relation; the calculated thermal wind shears suggest a decay with height of the jet system toward a state of uniform eatward flow. The existence of the approximately symmetric zonal winds and associated temperature gradients in the presence of a large-scale seasonal thermal response suggests that the jet system is driven at depths substantially below the levels where seasonally modulated insolation is important ($\text{p?0.5}\text{bar}$).     

Westward rotation of the atmospheric angular momentum vector of Titan by thermal tides

The orientation of the atmospheric angular momentum vector of Titan and its temporal variation predicted by a general circulation model are analysed and interpreted. The atmospheric angular momentum vector is tilted by a few degrees from the polar axis and the vector rotates (precesses) westward with a constant period of 1 Titan day. The fast westward rotation is likely to be caused by migrating diurnal thermal tides. The tilt is almost cancelled out in the troposphere by the wavenumber 2 pattern of Saturn's gravitational tide, but is more pronounced in the stratosphere, where thermal tides are significant. The predicted tilt angle and the equatorial angular momentum vary with season and maximize when the hemispheric asymmetry of the axial angular momentum or superrotation attains its peak.     

On the maintenance of thermal wind balance and equatorial superrotation in Titan's stratosphere

Titan's atmospheric winds, like those on Venus, exhibit superrotation at high altitudes. Titan general circulation models have yielded conflicting results on whether prograde winds in excess of 100 m s⁻¹ at the 1 mbar level are possible based on known physical processes that drive wind systems. A comprehensive two-dimensional (2D) model for Titan's stratosphere was constructed to systematically explore the physical mechanisms that produce and maintain stratospheric wind systems. To ensure conservation of angular momentum in the limit of no net exchange of atmospheric angular momentum with the solid satellite and no external sources and sinks, the zonal momentum equation was solved in flux form for total angular momentum. The relationships among thermal wind balance, meridional circulation, and zonal wind were examined with numerical experiments over a range of values for fundamental input parameters, including planetary rotation rate, radius, internal friction due to wave stresses, and net radiative drive. The magnitude of mid-latitude jets is most sensitive to a single parameter, the planetary rotation rate and results from the conversion of planetary angular momentum to relative angular momentum by the meridional circulation, whereas the strength of meridional circulation is mainly determined by the magnitude of the radiative drive. For Titan's slowly rotating atmosphere, the meridional temperature gradient is vanishingly small, even when the radiative drive is enhanced beyond reasonable magnitudes, and can be inferred from zonal winds in gradient/thermal wind balance. In our 2D model large equatorial superrotation in Titan's stratosphere can be only produced through internal drag forcing by eddy momentum fluxes, which redistribute angular momentum within the atmosphere, while still conserving the total angular momentum of the atmosphere with time. We cannot identify any waves, such as gravitational or thermal tides, that are sufficiently capable of generating the required eddy forcing of >50 m s⁻¹ Titan-day⁻¹ to maintain peak prograde winds in excess of 100 m s⁻¹ at the 1 mbar level.     

Constraints on Saturn’s tropospheric general circulation from Cassini ISS images

An automated cloud tracking algorithm is applied to Cassini Imaging Science Subsystem high-resolution apoapsis images of Saturn from 2005 and 2007 and moderate resolution images from 2011 and 2012 to define the near-global distribution of zonal winds and eddy momentum fluxes at the middle troposphere cloud level and in the upper troposphere haze. Improvements in the tracking algorithm combined with the greater feature contrast in the northern hemisphere during the approach to spring equinox allow for better rejection of erroneous wind vectors, a more objective assessment at any latitude of the quality of the mean zonal wind, and a population of winds comparable in size to that available for the much higher contrast atmosphere of Jupiter. Zonal winds at cloud level changed little between 2005 and 2007 at all latitudes sampled. Upper troposphere zonal winds derived from methane band images are ～10 m s^?1 weaker than cloud level winds in the cores of eastward jets and ～5 m s^?1 stronger on either side of the jet core, i.e., eastward jets appear to broaden with increasing altitude. In westward jet regions winds are approximately the same at both altitudes. Lateral eddy momentum fluxes are directed into eastward jet cores, including the strong equatorial jet, and away from westward jet cores and weaken with increasing altitude on the flanks of the eastward jets, consistent with the upward broadening of these jets. The conversion rate of eddy to mean zonal kinetic energy at the visible cloud level is larger in eastward jet regions (5.2 × 10^?5 m² s^?3) and smaller in westward jet regions (1.6 × 10^?5 m² s^?3) than the global mean value (4.1 × 10^?5 m² s^?3). Overall the results are consistent with theories that suggest that the jets and the overturning meridional circulation at cloud level on Saturn are maintained at least in part by eddies due to instabilities of the large-scale flow near and/or below the cloud level.     

Dune-forming winds on Titan and the influence of topography

Numerous extended dunes on Saturn's moon Titan detected by the Cassini RADAR constrain the long-term pattern of surface winds. We analyse the statistics of surface wind speed and direction and their spatial and temporal variability predicted by a general circulation model (GCM) in order to constrain surface wind predictions of this GCM by dune observations. The model shows that modern winds are sufficient for saltation and dune formation at low latitudes, in agreement with the presence of dunes there. The best condition for the dune-forming wind occurs with a threshold friction speed of 0.02 m s⁻¹ or slightly less. The equatorial region is conducive to longitudinal dunes because of a combined effect of a high sand drift potential and obtuse bimodal wind pattern oblique to the equator caused by the seasonal reversal of the Hadley circulation. The cross-equatorial wind is steady, and is characterised by a high Weibull shape parameter (k∼4). The wind pattern at higher latitudes is more complex and gusty, and neither longitudinal nor transversal dunes would be able to form. Putative large-scale topography is found to have a profound influence on the near-surface wind pattern. Generally mountains cause a convergence and speeding up on the flank, while basins weaken the wind and cause a divergent flow. Longitudinal dunes can be deflected on the foot of mountains by up to 90°. If Xanadu is a hypothetical large mountain, a wind pattern converging in Xanadu that entirely disagrees with the dune observations is predicted. If instead Xanadu is a large basin, the wind arcs clockwise north of Xanadu and anti-clockwise west and southwest of Xanadu, in agreement with the dune orientations in the vicinity of Xanadu. The albedo pattern has comparatively little influence on the wind field. Isolated mountain chains cause only local-scale change in the wind pattern. However, the persistent surface easterlies in Belet, which are in conflict with the dune orientations, do not disappear by any combination of large-scale topography.     

Wind variations in Jupiter's equatorial atmosphere: A QQO counterpart?

Jupiter's equatorial atmosphere, much like the Earth's, is known to show quasi-periodic variations in temperature, particularly in the stratosphere, but variations in other jovian atmospheric tracers have not been studied for any correlations to these oscillations. Data taken at NASA's Infrared Telescope Facility (IRTF) from 1979 to 2000 were used to obtain temperatures at two levels in the atmosphere, corresponding to the upper troposphere (250 mbar) and to the stratosphere (20 mbar). We find that the data show periodic signals at latitudes corresponding to the troposphere zonal wind jets, with periods ranging from 4.4 (stratosphere, 95% confidence at 4° S planetographic latitude) to 7.7 years (troposphere, 97% confidence at 6° N). We also discuss evidence that at some latitudes the troposphere temperature variations are out of phase from the stratosphere variations, even where no periodicity is evident. Hubble Space Telescope images were used, in conjunction with Voyager and Cassini data, to track small changes in the troposphere zonal winds from 20° N to 20° S latitude over the 1994-2000 time period. Oscillations with a period of 4.5 years are found near 7°-8° S, with 80-85% significance. Further, the strongest evidence for a QQO-induced tropospheric wind change tied to stratospheric temperature change occurs near these latitudes, though tropospheric temperatures show little periodicity here. Comparison of thermal winds and measured zonal winds for three dates indicate that cloud features at other latitudes are likely tracked at pressures that can vary by up to a few hundred millibar, but the cloud altitude change required is too large to explain the wind changes measured at 7° S.     

Numerical simulation of the general circulation of the atmosphere of Titan

The atmospheric circulation of Titan is investigated with a general circulation model. The representation of the large-scale dynamics is based on a grid point model developed and used at Laboratoire de Météorologie Dynamique for climate studies. The code also includes an accurate representation of radiative heating and cooling by molecular gases and haze as well as a parametrization of the vertical turbulent mixing of momentum and potential temperature. Long-term simulations of the atmospheric circulation are presented. Starting from a state of rest, the model spontaneously produces a strong superrotation with prograde equatorial winds (i.e., in the same sense as the assumed rotation of the solid body) increasing from the surface to reach 100 m sec-1 near the 1-mbar pressure level. Those equatorial winds are in very good agreement with some indirect observations, especially those of the 1989 occultation of Star 28-Sgr by Titan. On the other hand, the model simulates latitudinal temperature contrasts in the stratosphere that are significantly weaker than those observed by Voyager 1 which, we suggest, may be partly due to the nonrepresentation of the spatial and temporal variations of the abundances of molecular species and haze. We present diagnostics of the simulated atmospheric circulation underlying the importance of the seasonal cycle and a tentative explanation for the creation and maintenance of the atmospheric superrotation based on a careful angular momentum budget.     

A problem in resonance

Saturn's satellite Hyperion experiences large perturbations by Titan, the largest of Saturn's satellites, because of the closeness of Hyperion's orbital period to three-quarters of that of Titan.The motion of Hyperion is a superposition of periodic fluctuations (both free and forced) onto a motion which is periodic in a suitable uniformly rotating frame of reference, and in which Hyperion would be at a maximum distance from Saturn at each conjunction with Titan.Successive attempts to determine the mass of Titan from observed perturbations of Hyperion have suffered from omissions in the theory of terms subsequently found to be significant.An attempt is in progress which, it is believed, comprises all long-period changes in the osculating elements which are of second degree in the mass of Titan, and of third degree in the eccentricity of its orbit. Results so far obtained indicate that the period of the free motion of the orbit plane of Hyperion is better determined by Woltjer's reduction of the observations than by his theoretical calculations.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.     

Saturn’s zonal wind profile in 2004-2009 from Cassini ISS images and its long-term variability

Five years of Cassini Imaging Science Subsystem images, from 2004 to 2009, are analyzed in this work to retrieve global zonal wind profiles of Saturn’s northern and southern hemispheres in the methane absorbing bands at 890 and 727 nm and in their respective adjacent continuum wavelengths of 939 and 752 nm. A complete view of Saturn’s global circulation, including the equator, at two pressure levels, in the tropopause (60 mbar to 250 mbar with the MT filters) and in the upper troposphere (from ∼350 mbar to ∼500 mbar with the CB filter set), is presented. Both zonal wind profiles (available at the Supplementary Material Section), show the same structure but with significant differences in the peak of the eastward jets and the equatorial region, including a region of positive vertical shear symmetrically located around the equator between the 10° < |φ_c| < 25° where zonal velocities close to the tropopause are higher than at 500 mbar. A comparison of previously published zonal wind sets obtained by Voyager 1 and 2 (1980-1981), Hubble Space Telescope, and ground-based telescopes (1990-2004) with the present Cassini profiles (2004-2009) covering a full Saturn year shows that the shape of the zonal wind profile and intensity of the jets has remained almost unchanged except at the equator, despite the seasonal insolation cycle and the variability of Saturn’s emitted power. The major wind changes occurred at equatorial latitudes, perhaps following the Great White Spot eruption in 1990. It is not evident from our study if the seasonal insolation cycle and its associated ring shadowing influence the equatorial circulation at cloud level.     

Zonal winds at high latitudes on Venus: An improved application of cyclostrophic balance to Venus Express observations

Recent retrievals of zonal thermal winds obtained in a cyclostrophic regime on Venus are generally consistent with cloud tracking measurements at mid-latitudes, but become unphysical in polar regions where the values obtained above the clouds are often less than or close to zero. Using a global atmospheric model, we show that the main source of errors that appear in the polar regions when retrieving the zonal thermal winds is most likely due to uncertainties in the zonal wind intensity in the choice of the lower boundary condition.Here we suggest a new and robust method to better estimate the lower boundary condition for high latitudes, thereby improving the retrieved zonal thermal winds throughout the high latitudes middle atmosphere. This new method is applied to temperature fields derived from Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) data on board the Venus Express spacecraft. We obtain a zonal thermal wind field that is in better agreement with other, more direct methods based on either retrieving the zonal winds from cloud tracking or from direct measurements of the meridional slope of pressure surfaces.     

The influence of the Martian polar caps on the diurnal tide

The large horizontal heating gradients that exist near the edge of the Martian polar caps during spring are shown to be capable of exciting large oscillations in the diurnal tide. To a lesser extent, the daily mass cycling between cap and atmosphere can also contribute. The calculations which demonstrate this are based on classical tidal theory as applied to the cylindrical coordinate system. This is done to facilitate the representation of the heating function. Results are presented for the horizontal surface winds only. They indicate a circulation at the cap edge somewhat analogous to the smaller scale terrestrial sea breeze. The amplitude of the zonal component is largest and is increased from 1 to 10 m sec^?1 by the modeled influence of the polar cap. When coupled with the basic flow these cap-edge tides can produce strong surface winds during spring. Such a mechanism may contribute to the ability of the south polar cap winds to generate the local dust storms observed near the cap edge at this season.     

Zonal winds near Venus' cloud top level: A model study of the interaction between the zonal mean circulation and the semidiurnal tide

A primitive equation wave-mean flow interaction model, designed by J. R. Holton and used originally to study Earth's middle atmosphere, has been adapted to Venus in order to clarify our understanding of the interaction between the semidiurnal tide and the thermally driven mean meridional circulation near the cloud top level. With or without the tide the model produces midlatitude jets whose structure is insensitive to vertical shear of the background angular velocity above and below the cloud top level, but it is sensitive to background angular velocity at the cloud top level. When this background angular velocity is close to that of Venus, the latitudes and speeds of these jets are similar to the latitudes and speeds of jets at the Venus cloud top level as inferred from observed temperatures and the cyclostrophic balance condition. In agreement with the hypothesis of Fels and Lindzen, the model tide accelerates the equatorial zonal wind near the cloud top level and decelerates it at higher levels. The tidal vertical wavelength, maximum amplitude, dissipative decay length, and zonal flow accelerations are sensitive functions of background angular velocity, in agreement with elementary gravity wave theory. In the equatorial cloud top region, tidal acceleration is comparable in magnitude to the decelerative effects of vertical advection and the model's Rayleigh friction damping. For sufficiently rapid initial zonal flow near the cloud top level, the area-weighted global mean cloud top level zonal wind increases with time over a 50-day model run as a result of tidal acceleration. Agreement between the model tide and the observed tide, or the tide determined in the more detailed calculations of Pechmann and Ingersoll, is best when the background angular velocity at the jet level is about 30% larger than that observed.     

The effect of gravitational and pressure torques on Titan's length-of-day variations

Cassini radar observations show that Titan's spin is slightly faster than synchronous spin. Angular momentum exchange between Titan's surface and the atmosphere over seasonal time scales corresponding to Saturn's orbital period of 29.5 year is the most likely cause of the observed non-synchronous rotation. We study the effect of Saturn's gravitational torque and torques between internal layers on the length-of-day (LOD) variations driven by the atmosphere. Because static tides deform Titan into an ellipsoid with the long axis approximately in the direction to Saturn, non-zero gravitational and pressure torques exist that can change the rotation rate of Titan. For the torque calculation, we estimate the flattening of Titan and its interior layers under the assumption of hydrostatic equilibrium. The gravitational forcing by Saturn, due to misalignment of the long axis of Titan with the line joining the mass centers of Titan and Saturn, reduces the LOD variations with respect to those for a spherical Titan by an order of magnitude. Internal gravitational and pressure coupling between the ice shell and the interior beneath a putative ocean tends to reduce any differential rotation between shell and interior and reduces further the LOD variations by a few times. For the current estimate of the atmospheric torque, we obtain LOD variations of a hydrostatic Titan that are more than 100 times smaller than the observations indicate when Titan has no ocean as well as when a subsurface ocean exists. Moreover, Saturn's torque causes the rotation to be slower than synchronous in contrast to the Cassini observations. The calculated LOD variations could be increased if the atmospheric torque is larger than predicted and or if fast viscous relaxation of the ice shell could reduce the gravitational coupling, but it remains to be studied if a two order of magnitude increase is possible and if these effects can explain the phase difference of the predicted rotation variations. Alternatively, the large differences with the observations may suggest that non-hydrostatic effects in Titan are important. In particular, we show that the amplitude and phase of the calculated rotation variations are similar to the observed values if non-hydrostatic effects could strongly reduce the equatorial flattening of the ice shell above an internal ocean.     

Another look at atmospheric dynamics on Titan and some of its general consequences

Mean wind velocities, U, and horizontal temperature differences, δT, are estimated for the Titan atmosphere using the similarity theory of the author. It is found that U is of order 1 m/sec and δT ～ 0.1 K. The last value agrees with its estimate by Leovy and Pollack (1973, Icarus19, 195–201); however the values of U are an order of magnitude less. While analyzing the causes of the difference it is found that the circulation models developed in I overestimate considerably the atmospheric efficiency in transformation of solar energy into the kinetic energy of motions. Possible reasons for such an overestimate are discussed. We conclude that the efficiency coefficient is a very sensitive characteristic of a circulation regime and that its determination is an efficient means for checking the correspondence of various circulation models with reality. Arguments are presented stressing the conclusion of I that the Titan atmospheric circulation is in the Hadley symmetric regime, which is strongly influenced by the satellite's own rotation. At the same time a thermal tide should be a noticeable feature of the circulation. In the upper part of the Titan atmosphere something like the phenomenon of the 4-day Venus circulation may be developed. It is noted that the analogy between the Titan and Venus atmospheric circulations might be a very close one.     

Remote sensing of mesospheric winds with the high-resolution doppler imager

Observations of the winds in the upper atmosphere obtained with the High-Resolution Doppler Imager (HRDI) on the Upper Atmosphere Research Satellite (UARS) are discussed. This instrument is a very stable high-resolution triple-etalon Fabry-Perot interferometer, which is used to observe the slight Doppler shifts of absorption and emission lines in the O₂ Atmospheric bands induced by atmospheric motions. Preliminary observations indicate that the winds in the mesosphere and lower thermosphere are a mixture of migrating and non-migrating tides, and planetary-scale waves. The mean meridional winds are dominated by the 1,1 diurnal tide which is easily extracted from the daily zonal means of the satellite observations. The daily mean zonal winds are a mixture of the diurnal tide and a zonal flow which is consistent with theoretical expectations.     

Zonal mean properties of Jupiter's upper troposphere from voyager infrared observations

The highest spatial resolution Voyager IRIS spectra are used to produce zonal averages of the temperature at the 150- and 270-mb pressure levels, of the para-hydrogen fraction at 270 mb, of the ammonia abundance near the 680-mb level, and of two infrared cloud optical depths, one near 5 μm and one near 45 μm wavelength. There are two cloud components, one uniformly distributed and only apparent at 5 μm, and another that correlates strongly with the ammonia abundance and that is apparent at both 5 and 45 μm. From the ratio of optical depths at the two wavelengths, the particles in the variable cloud are between 3 and 10 μm in radius. This cloud is located near the ammonia condensation level. The other particles are either smaller or deeper. The cloud and ammonia distribution is consistent with concentration by upward vertical motion at the equatorward edges of prograde atmospheric jets. The temperature field is also consistent with such vertical motion, with radiative heating balancing adiabatic expansional cooling. The para-hydrogen distribution also appears consistent, but noise levels are high. The thermal wind shear indicates decay of the jets with height within the upper troposphere, with a vertical scale of two or three scale heights. The entire set of upper troposphere data is consistent with a simple axisymmetric dynamical model with Coriolis acceleration of the zonal wind balanced by a linear drag. The meridional residual mean circulation in the model, if interpreted also as a Lagrangian mean circulation, would explain nicely the distribution of ammonia and para-hydrogen. The circulation is a response to a deeper tropospheric flow of unknown origin. However, the horizontal scale of jets is on the order of the deformation radius based on a scale height at the base of the upper troposphere. It is conjectured that the physics of the flow may require this to be true, and may also require that the relative vorticity gradient be of the same order as the planetary vorticity gradient, thereby fixing both the dimensions and amplitudes of the jets.     

Characterization of zonal winds in the stratosphere of Titan with UVES

Titan has been observed with UVES, the UV-Visual Echelle Spectrograph at the Very Large Telescope, with the aim of characterizing the zonal wind flow. We use a retrieval scheme originally developed for absolute stellar accelerometry [Connes, P., 1985. Astrophys. Space Sci., 110, 211-255] to extract the velocity signal by simultaneously taking into account all the lines present in the spectrum. The method allows to measure the Doppler shift induced at a given point by the zonal wind flow, with high precision. The short-wavelength channel (4200-5200 Å) probes one scale height higher than the long-wavelength one (5200-6200 Å), and we observe statistically significant evidence for stronger winds at higher altitudes. The results show a high dispersion. Globally, we detect prograde zonal winds, with lower limits of 62 and 50 m s⁻¹ at the regions centered at 200 and 170 km altitude, but approximately a quarter of the measurements indicates null or retrograde winds.