Phosphine on Jupiter and Saturn from Cassini/CIRS

The global distribution of phosphine (PH₃) on Jupiter and Saturn is derived using 2.5 cm⁻¹ spectral resolution Cassini/CIRS observations. We extend the preliminary PH₃ analyses on the gas giants [Irwin, P.G.J., and 6 colleagues, 2004. Icarus 172, 37-49; Fletcher, L.N., and 9 colleagues, 2007a. Icarus 188, 72-88] by (a) incorporating a wider range of Cassini/CIRS datasets and by considering a broader spectral range; (b) direct incorporation of thermal infrared opacities due to tropospheric aerosols and (c) using a common retrieval algorithm and spectroscopic line database to allow direct comparison between these two gas giants.The results suggest striking similarities between the tropospheric dynamics in the 100-1000 mbar regions of the giant planets: both demonstrate enhanced PH₃ at the equator, depletion over neighbouring equatorial belts and mid-latitude belt/zone structures. Saturn's polar PH₃ shows depletion within the hot cyclonic polar vortices. Jovian aerosol distributions are consistent with previous independent studies, and on Saturn we demonstrate that CIRS spectra are most consistent with a haze in the 100-400 mbar range with a mean optical depth of 0.1 at 10 μm. Unlike Jupiter, Saturn's tropospheric haze shows a hemispherical asymmetry, being more opaque in the southern summer hemisphere than in the north. Thermal-IR haze opacity is not enhanced at Saturn's equator as it is on Jupiter.Small-scale perturbations to the mean PH₃ abundance are discussed both in terms of a model of meridional overturning and parameterisation as eddy mixing. The large-scale structure of the PH₃ distributions is likely to be related to changes in the photochemical lifetimes and the shielding due to aerosol opacities. On Saturn, the enhanced summer opacity results in shielding and extended photochemical lifetimes for PH₃, permitting elevated PH₃ levels over Saturn's summer hemisphere.     

Stability of jets on Jupiter and Saturn

The linear stability of a zonal jet that decays with depth is investigated under the assumption that the thermal stratification is very small. A westerly cosine jet is found to be more stable than it is in a thin fluid shell with two-dimensional flow. This is in agreement with observations of Jupiter and Saturn, where jet curvature exceeds the barotropic stability criterion. This result constitutes an alternative hypothesis to that of Ingersoll and Pollard [Icarus 52 (1982) 62], who showed that deep jets extending through the interior are also more stable than thin shell jets. The flow regime assumed in the present work requires that a small stratification can exist and persist even in the presence of horizontal temperature gradients. Further work will be needed to test whether this is realistic.     

Interaction of moist convection with zonal jets on Jupiter and Saturn

Observations suggest that moist convection plays an important role in the large-scale dynamics of Jupiter's and Saturn's atmospheres. Here we use a reduced-gravity quasigeostrophic model, with a parameterization of moist convection that is based on observations, to study the interaction between moist convection and zonal jets on Jupiter and Saturn. Stable jets with approximately the same width and strength as observations are generated in the model. The observed zonal jets violate the barotropic stability criterion but the modeled jets do so only if the flow in the deep underlying layer is westward. The model results suggest that a length scale and a velocity scale associated with moist convection control the width and strength of the jets. The length scale and velocity scale offer a possible explanation of why the jets of Saturn are stronger and wider than those of Jupiter.     

Turbulent convection in rapidly rotating spherical shells: A model for equatorial and high latitude jets on Jupiter and Saturn

The origin of zonal jets on the jovian planets has long been a topic of scientific debate. In this paper we show that deep convection in a spherical shell can generate zonal flow comparable to that observed on Jupiter and Saturn, including a broad prograde equatorial jet and multiple alternating jets at higher latitudes. We present fully turbulent, 3D spherical numerical simulations of rapidly rotating convection with different spherical shell geometries. The resulting global flow fields tend to be segregated into three regions (north, equatorial, and south), bounded by the tangent cylinder that circumscribes the inner boundary equator. In all of our simulations a strong prograde equatorial jet forms outside the tangent cylinder, whereas multiple jets form in the northern and southern hemispheres, inside the tangent cylinder. The jet scaling of our numerical models and of Jupiter and Saturn is consistent with the theory of geostrophic turbulence, which we extend to include the effect of spherical shell geometry. Zonal flow in a spherical shell is distinguished from that in a full sphere or a shallow layer by the effect of the tangent cylinder, which marks a reversal in the sign of the planetary β-parameter and a jump in the Rhines length. This jump is manifest in the numerical simulations as a sharp equatorward increase in jet widths—a transition that is also observed on Jupiter and Saturn. The location of this transition gives an estimate of the depth of zonal flow, which seems to be consistent with current models of the jovian and saturnian interiors.     

Long slit spectropolarimetry of Jupiter and Saturn

We present ground-based limb polarization measurements of Jupiter and Saturn consisting of full disk imaging polarimetry for the wavelength 7300 Å and spatially resolved (long-slit) spectropolarimetry covering the wavelength range 5200-9350 Å.For the polar region of Jupiter we find for λ = 6000 Å a very strong radial (perpendicular to the limb) fractional polarization with a seeing corrected maximum of about +11.5% in the South and +10.0% in the North. This indicates that the polarizing haze layer is thicker at the South pole. The polar haze layers extend down to 58° in latitude. The derived polarization values are much higher than reported in previous studies because of the better spatial resolution of our data and an appropriate consideration of the atmospheric seeing. Model calculations demonstrate that the high limb polarization can be explained by strongly polarizing (p ≈ 1.0), high albedo (ω ≈ 0.98) haze particles with a scattering asymmetry parameter of g ≈ 0.6 as expected for aggregate particles of the type described by West and Smith (West, R.A., Smith, P.H. [1991]. Icarus 90, 330-333). The deduced particle parameters are distinctively different when compared to lower latitude regions.The spectropolarimetry of Jupiter shows a decrease in the polar limb polarization towards longer wavelengths and a significantly enhanced polarization in strong methane bands when compared to the adjacent continuum. This is a natural outcome for a highly polarizing haze layer above an atmosphere where multiple scatterings are suppressed in absorption bands. For lower latitudes the fractional polarization is small, negative, and it depends only little on wavelength except for the strong CH₄-band at 8870 Å.The South pole of Saturn shows a lower polarization (p ≈ 1.0-1.5%) than the poles of Jupiter. The spectropolarimetric signal for Saturn decrease rapidly with wavelength and shows no significant enhancements in the fractional polarization in the absorption bands. These properties can be explained by a vertically extended stratospheric haze region composed of small particles <100 nm as suggested previously by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2005]. Icarus 179, 195-221).In addition we find in the V- and R-band a previously not observed strong polarization feature (p = 1.5-2.0%) near the equator of Saturn. The origin of this polarization signal is unclear but it could be related to a seasonal effect.Finally we discuss the potential of ground-based limb polarization measurements for the investigation of the scattering particles in the atmospheres of Jupiter and Saturn.     

Kinetic models for the exospheres of Jupiter and Saturn

The exospheric theory based on the Kappa velocity distribution function (VDF) is used to model the exosphere of the giant planets Jupiter and Saturn. Such Kappa velocity distribution functions with an enhanced population of suprathermal particles are indeed often observed in space plasmas and in the space environment of the planets. The suprathermal particles have significant effects on the escape flux, density and temperature profiles of the particles in the exosphere of the giant planets. The polar wind flux becomes several orders larger when suprathermal electrons are considered, so that the planetary ionosphere becomes then a significant source for their inner magnetosphere. Moreover, the number density of the particles decreases slower as a function of the altitude when a Kappa distribution is considered instead of a Maxwellian one. Two-dimensional maps of density are calculated for typical values of the temperatures. The exospheric formalism is also applied to study the escape flux from the exospheres of Io and Titan, respectively, moons of Jupiter and Saturn.     

Simulations of high-latitude spots on Jupiter: Constraints on vortex strength and the deep wind

We combine high-resolution observations of the dynamical behavior of small vortices (diameters ?5000 km) located at latitude 60°N on Jupiter with forward modeling, using the EPIC atmospheric model, to address two open questions: the dependence of the zonal winds with depth, and the strength of vortices that are too small to apply cloud tracking to their internal structure. The observed drift rates of the vortices can only be reproduced in the model when the zonal winds increase slightly with depth below the cloud tops, with a vertical shear that is less than was measured at 7°N at the southern rim of a 5-μm hotspot by the Galileo Probe Doppler Wind Experiment (DWE). This supports the idea that Jupiter's vertical shear may vary significantly with latitude. Our simulations suggest that the morphology of the mergers between vortices mainly depends on their maximum tangential velocities, the best results occurring when the tangential velocity is close to the velocity difference of the alternating jets constraining the zone in which the vortices are embedded. We use this correlation, together with the high-resolution data available for the White Ovals, to derive an empirical relationship between the maximum tangential velocity of a jovian vortex and its size, normalized by the strength and size of the encompassing shear zone. The Great Red Spot stands out as a significant anomaly to this relationship, but interestingly it is becoming less so with time.     

The Structure of Jupiter, Saturn, and Exoplanets: Key Questions for High-Pressure Experiments

We give an overview of our current understanding of the structure of gas giant planets, from Jupiter and Saturn to extrasolar giant planets. We focus on addressing what high-pressure laboratory experiments on hydrogen and helium can help to elucidate about the structure of these planets.     

Phase separation in giant planets: inhomogeneous evolution of Saturn

We present the first models of Jupiter and Saturn to couple their evolution to both a radiative-atmosphere grid and to high-pressure phase diagrams of hydrogen with helium and other admixtures. We find that prior calculated phase diagrams in which Saturn's interior reaches a region of predicted helium immiscibility do not allow enough energy release to prolong Saturn's cooling to its known age and effective temperature. We explore modifications to published phase diagrams that would lead to greater energy release, and propose a modified H-He phase diagram that is physically reasonable, leads to the correct extension of Saturn's cooling, and predicts an atmospheric helium mass fraction Y_atmos=0.185, in agreement with recent estimates. We also explore the possibility of internal separation of elements heavier than helium, and find that, alternatively, such separation could prolong Saturn's cooling to its known age and effective temperature under a realistic phase diagram and heavy element abundance (in which case Saturn's Y_atmos would be solar but heavier elements would be depleted). In none of these scenarios does Jupiter's interior evolve to any region of helium or heavy-element immiscibility: Jupiter evolves homogeneously to the present day. We discuss the implications of our calculations for Saturn's primordial core mass.     

Renewed acceleration of the 24° N jet on Jupiter

Jupiter's eastward jet at 24° N, which formerly had the fastest winds on the planet, has maintained a less extreme speed of ∼135 m/s since 1991, carrying a series of long-lived vortices at 125 m/s. In 2002-2003, as the albedo of the adjacent North Temperate Belt increased, the tracks of the vortices accelerated slightly, and they had disappeared by 2005. In 2005, small tracers had a mean speed of 146.4 (±0.9) m/s, significantly faster than the previous mean speed of the jet, suggesting that the jet peak itself has accelerated at cloud-top level, and that the jet is beginning to return to the super-fast state. These changes may resemble the even greater transformations occurring in the equatorial jet of Saturn.     

Time-dependent flexure of the lithospheres on the icy satellites of Jupiter and Saturn

Previous analyses into flexural deformation on the icy satellites of Jupiter and Saturn have assumed static, elastic lithospheres. Viscous creep within the lithosphere, however, can cause evolution over time. Here, we apply a finite-element model that employs a time-dependent elastic–viscous-plastic rheology in order to investigate flexure on icy satellites. Factors that affect this time-dependent response are those that control creep rates; surface temperature, heat flow, and grain size. Our results show that surface temperature is by far the dominant factor. At higher surface temperatures (100–130 K), the evolution of the deformation is such that the thickness of a modeled elastic lithosphere could vary by up to an order of magnitude, depending on the time scale over which the deformation occurred. Because the flexure observed on icy satellites generally indicates transient high heat flow events, our results indicate that the duration of the heat pulse is an important factor. For the icy worlds of Jupiter and Saturn, static models of lithospheric flexure should be used with caution.     

The 10 October 1999 HIP 9369 occultation by the northern polar region of Jupiter: ingress and egress lightcurves analysis

The occultation of bright star HIP9369 by the northern polar region of Jupiter was observed from four locations in North and South America, providing four data sets for ingress and egress. The inversion of the eight occultation lightcurves provides temperature profiles at different latitudes ranging from 55°N to 73.2°N. We estimate the errors on the profiles due to the uncertainties of the inversion method and compare the value of the temperature at the deepest level probed (∼ 50 μbar) with previous observations. The shape of the temperature gradient profile is found similar to previous investigations of planetary atmospheres with propagating and breaking gravity waves. We analyze the small scale structures in both lightcurves and temperature profiles using the continuous wavelet transform. The calculated power spectra of localized fluctuations in the temperature profiles show slopes close to −3 for all eight profiles. We also isolate and reconstruct the three-dimensional geometry of a single wave mode with vertical and horizontal wavelengths of respectively 3 and 70 km. The identified wave is consistent with the gravity wave regime, with a horizontal phase speed nearly parallel to the planetary meridian. Nevertheless, the dissipation of the corresponding wave in Jupiter’s stratosphere should preclude its detection at the observed levels and an acoustic wave cannot be ruled out.     

Emergence of polar-jet polygons from jet instabilities in a Saturn model

Voyager flybys of Saturn in 1980-1981 revealed a circumpolar wave at ≈78° north planetographic latitude. The feature had a dominant wavenumber 6 mode, and has been termed the Hexagon from its geometric appearance in polar-projected mosaics. It was also noted for being stationary with respect to Saturn’s Kilometric Radiation (SKR) rotation rate. The Hexagon has persisted for over 30 years since the Voyager observations until now. It has been observed from ground based telescopes, Hubble Space Telescope and multiple instruments onboard Cassini in orbit around Saturn. Measurements of cloud motions in the region reveal the presence of a jet stream whose path closely follows the Hexagon’s outline. Why the jet stream takes the characteristic six-sided shape and how it is stably maintained across multiple saturnian seasons are yet to be explained. We present numerical simulations of the 78.3°N jet using the Explicit Planetary Isentropic-Coordinate (EPIC) model and demonstrate that a stable hexagonal structure can emerge without forcing when dynamic instabilities in the zonal jet nonlinearly equilibrate. For a given amplitude of the jet, the dominant zonal wavenumber is most strongly dependent on the peak curvature of the jet, i.e., the second north-south spatial derivative of the zonal wind profile at the center of the jet. The stable polygonal shape of the jet in our simulations is formed by a vortex street with cyclonic and anticyclonic vortices lining up towards the polar and equatorial side of the jet, respectively. Our result is analogous to laboratory experiments of fluid motions in rotating tanks that develop polygonal flows out of vortex streets. However, our results also show that a vortex street model of the Hexagon cannot reproduce the observed propagation speed unless the zonal jet’s speed is modified beyond the uncertainties in the observed zonal wind speed, which suggests that a vortex street model of the Hexagon and the observed zonal wind profile may not be mutually compatible.     

UV and IR auroral emission model for the outer planets: Jupiter and Saturn comparison

Planetary aurora display the dynamic behavior of the plasma gas surrounding a planet. The outer planetary aurora are most often observed in the ultraviolet (UV) and the infrared (IR) wavelengths. How the emissions in these different wavelengths are connected with the background physical conditions are not yet well understood. Here we investigate the sensitivity of UV and IR emissions to the incident precipitating auroral electrons and the background atmospheric temperature, and compare the results obtained for Jupiter and Saturn. We develop a model which estimates UV and IR emission rates accounting for UV absorption by hydrocarbons, ion chemistry, and non-LTE effects. Parameterization equations are applied to estimate the ionization and excitation profiles in the H₂ atmosphere caused by auroral electron precipitation. The dependences of UV and IR emissions on electron flux are found to be similar at Jupiter and Saturn. However, the dependences of the emissions on electron energy are different at the two planets, especially for low energy (<10 keV) electrons; the UV and IR emissions both decrease with decreasing electron energy, but this effect in the IR is less at Saturn than at Jupiter. The temperature sensitivity of the IR emission is also greater at Saturn than at Jupiter. These dependences are interpreted as results of non-LTE effects on the atmospheric temperature and density profiles. The different dependences of the UV and IR emissions on temperature and electron energy at Saturn may explain the different appearance of polar emissions observed at UV and IR wavelengths, and the differences from those observed at Jupiter. These results lead to the prediction that the differences between the IR and UV aurora at Saturn may be more significant than those at Jupiter. We consider in particular the occurrence of bright polar infrared emissions at Saturn and quantitatively estimate the conditions for such IR-only emissions to appear.     

Scaling laws for convection and jet speeds in the giant planets

Three dimensional studies of convection in deep spherical shells have been used to test the hypothesis that the strong jet streams on Jupiter, Saturn, Uranus, and Neptune result from convection throughout the molecular envelopes. Due to computational limitations, these simulations must be performed at parameter settings far from jovian values and generally adopt heat fluxes 5-10 orders of magnitude larger than the planetary values. Several numerical investigations have identified trends for how the mean jet speed varies with heat flux and viscosity in these models, but no previous theories have been advanced to explain these trends. Here, we show using simple arguments that if convective release of potential energy pumps the jets and viscosity damps them, the mean jet speeds split into two regimes. When the convection is weakly nonlinear, the equilibrated jet speeds should scale approximately with F/ν, where F is the convective heat flux and ν is the viscosity. When the convection is strongly nonlinear, the jet speeds are faster and should scale approximately as (F/ν)^1/2. We demonstrate how this regime shift can naturally result from a shift in the behavior of the jet-pumping efficiency with heat flux and viscosity. Moreover, both Boussinesq and anelastic simulations hint at the existence of a third regime where, at sufficiently high heat fluxes or sufficiently small viscosities, the jet speed becomes independent of the viscosity. We show based on mixing-length estimates that if such a regime exists, mean jet speeds should scale as heat flux to the 1/4 power. Our scalings provide a good match to the mean jet speeds obtained in previous Boussinesq and anelastic, three-dimensional simulations of convection within giant planets over a broad range of parameters. When extrapolated to the real heat fluxes, these scalings suggest that the mass-weighted jet speeds in the molecular envelopes of the giant planets are much weaker—by an order of magnitude or more—than the speeds measured at cloud level.     

Potential for stratospheric Doppler windspeed measurements of Jupiter by sub-millimetre spectroscopy

The sub-millimetre/microwave range of the spectrum has been exploited in the field of Earth observation by many instruments over the years and has provided a plethora of information on atmospheric chemistry and dynamics—however, this spectral range has not been fully explored in planetary science, having been exclusively employed to carry out ground-based measurements. To this end, a sub-millimetre instrument, the Orbiter Terahertz Infrared Spectrometer (ORTIS), is studied by the University of Oxford and the Rutherford Appleton Laboratory, to meet the requirements of the European Space Agency's Cosmic Visions 2015-2025 programme—in particular, the Europa Jupiter System Mission (EJSM), which has the European Space Agency and the National Aeronautics and Space Administration as partners. ORTIS is designed to measure atmospheric temperature, the abundance of stratospheric water vapour and other jovian gases, and is intended to be capable of retrieving vertical profiles of horizontal windspeed in the stratosphere for the first time, from Doppler-shifted emission lines measured at high spectral resolution. In this work, a preliminary study and implementation of the estimation of windspeed profiles on simulated spectra representative of Jupiter is presented, detailing the development of the retrieval algorithm, showing that a sub-millimetre instrument such as ORTIS should be able to retrieve windspeed profiles to an accuracy of about 15 m/s between 70 and 200 km/0.1-10 mb using a single near-limb measurement, for expected noise amplitudes.     

Compressible convection in the deep atmospheres of giant planets

Fast rotating giant planets such as Jupiter and Saturn possess alternate prograde and retrograde zonal winds which are stable over long periods of time. We consider a compressible model of convection in a spherical shell with rapid rotation, using the anelastic approximation, to explore the parameter range for which such zonal flows can be produced.We consider models with a large variation in density across the layer. Our models are based only on the molecular H/He region above the metallic hydrogen transition at about 2 Mbar, and we do not include the hydromagnetic effects which may be important if the electrical conductivity is significant. We find that the convective velocities are significantly higher in the low density regions of the shell, but the zonal flow is almost independent of the z-coordinate parallel to the rotation axis. We analyse how this behaviour is consistent with the Proudman-Taylor theorem.We find that deep prograde zonal flow near the equator is a very robust feature of our models. Prograde and retrograde jets alternating in latitude can occur inside the tangent cylinder in compressible as well as Boussinesq models, particularly at lower Prandtl numbers. However, the zonal jets inside the tangent cylinder are suppressed if a no-slip condition is imposed at the inner boundary. This suggests that deep high latitude jets may be suppressed if there is significant magnetic dissipation.Our compressible calculations include the viscous dissipation in the entropy equation, and we find this is comparable to, and in some cases exceeds, the total heat flux emerging from the surface. For numerical reasons, these simulations cannot reach the extremely low Ekman number found in giant planets, and they necessarily also have a much larger heat flux than planets. We therefore discuss how our results might scale down to give solutions with lower dissipation and lower heat flux.     

Hadley circulations and Kelvin wave-driven equatorial jets in the atmospheres of Jupiter and Saturn

We propose a dynamical mechanism that can plausibly explain the origin of the broad prograde equatorial winds observed on Jupiter and Saturn, and examine the feasibility of this mechanism using two- (2D) and three-dimensional (3D) numerical simulation models. The idea is based on combining a narrow Gaussian jet peaking at the equator, which is induced by the momentum transfer from an upward propagating equatorial Kelvin-wave, and a pair of off-equatorial jets due to a meridional-vertical circulation similar to the tropical Hadley circulation on Earth. We employ for this feasibility study a 2D mechanistic mean-flow model which incorporates the influence of prescribed waves, and a 3D general circulation model, based on the generalised primitive equations of atmospheric motion. We then confirm that the dynamical models of both kinds can successfully reproduce theoretically expected flows of a reasonable magnitude, and that when two mechanisms are combined, a broad super-rotating jet is produced with off-equatorial maxima in zonal velocity for both Jupiter and Saturn, approximately in accordance with observations.     

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.     

Long-term evolution of the aerosol debris cloud produced by the 2009 impact on Jupiter

We present a study of the long-term evolution of the cloud of aerosols produced in the atmosphere of Jupiter by the impact of an object on 19 July 2009 (Sánchez-Lavega, A. et al. [2010]. Astrophys. J. 715, L155-L159). The work is based on images obtained during 5 months from the impact to 31 December 2009 taken in visible continuum wavelengths and from 20 July 2009 to 28 May 2010 taken in near-infrared deep hydrogen-methane absorption bands at 2.1-2.3 μm. The impact cloud expanded zonally from ∼5000 km (July 19) to 225,000 km (29 October, about 180° in longitude), remaining meridionally localized within a latitude band from 53.5°S to 61.5°S planetographic latitude. During the first two months after its formation the site showed heterogeneous structure with 500-1000 km sized embedded spots. Later the reflectivity of the debris field became more homogeneous due to clump mergers. The cloud was mainly dispersed in longitude by the dominant zonal winds and their meridional shear, during the initial stages, localized motions may have been induced by thermal perturbation caused by the impact’s energy deposition. The tracking of individual spots within the impact cloud shows that the westward jet at 56.5°S latitude increases its eastward velocity with altitude above the tropopause by 5-10 m s⁻¹. The corresponding vertical wind shear is low, about 1 m s⁻¹ per scale height in agreement with previous thermal wind estimations. We found evidence for discrete localized meridional motions with speeds of 1-2 m s⁻¹. Two numerical models are used to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the winds and their shears. The other uses the EPIC code, a nonlinear calculation of the evolution of the potential vorticity field generated by a heat pulse that simulates the impact. Both models reproduce the observed global structure of the cloud and the dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. The reflectivity of the impact cloud decreased exponentially with a characteristic timescale of 15 days; we can explain this behavior with a radiative transfer model of the cloud optical depth coupled to an advection model of the cloud dispersion by the wind shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) for the small aerosol particles forming the cloud is 45-200 days, thus aerosols were removed vertically over the long term following their zonal dispersion. No evidence of the cloud was detected 10 months after the impact.