Non-LTE polarized radiative transfer in intermediate magnetic fields: Numerical problems and results

This paper presents some numerical results relative to a solution, based on the density matrix formalism, of the non-LTE, polarized radiative transfer problem for a two-level atom. The results concern the atomic upper level population and alignment, and the emergent radiation Stokes profiles, for a plane-parallel, static, isothermal atmosphere embedded in a magnetic field of intermediate strength, such that the Zeeman splitting has to be taken into account in the line profile. Zeeman coherences are neglected, whereas magneto-optical effects are taken into account, resulting in a full 4×4 absorption matrix. Induced emission is neglected and complete frequency redistribution, in the rest and laboratory frames, is assumed. Pure Doppler absorption profile (gaussian shape) has also been assumed. The presentation of the results is preceded by a brief discussion of their accuracy and of the numerical difficulties that were met in the solution of the problem.On leave from the Dipartimento di Astronomia e Scienza dello Spazio, Università di Firenze, Largo E. Fermi 5, I-50125 Firenze, Italia     

The importance of continuum radiation for the stellar wind hydrodynamics of hot stars

The contribution of bound-free and free-free processes to the outward acceleration of Pup is studied and is found to be negligible.     

Stokes profiles inversion techniques

Inversion techniques of the radiative transfer equation for polarized light are presented as one of the best current procedures to infer the vector magnetic field, as well as other quantities governing the physical state of the atmospheric layers that photons are coming from. Several characteristics of the various available inversion procedures are pointed out. They are mostly based on the diagnostic contents of the spectral lines as well as on the main hypotheses assumed in these procedures. In particular, the role of gradients in the atmospheric quantities is emphasized as of paramount importance in any diagnostic analysis and, hence, in any interpretation of inversion results.     

Numerical methods in polarized radiative transfer

This paper looks at three aspects of numerical methods for solving polarized radiative transfer problems associated with spectral line formation in the presence of a magnetic field. First we prove Murphy's law for Stokes evolution operators which is the basis of the efficient algorithm used in the SPSR software package to compute the Stokes line depression contribution functions. Then we use a two-stream model to explain the efficacy of the field-free method in which the non-LTE line source function in a uniform magnetic field is approximated by the source function neglecting the magnetic field. Finally we introduce a totally new and computationally efficient approach to solving non-LTE problems based on a method of sparsely representing integral operators using wavelets. As an illustration, the wavelet method is used to solve the source function integral equation for a two-level atomic model in a finite atmosphere with coherent scattering, ignoring polarization.     

Particle size and abundance of HC3N ice in Titan’s lower stratosphere at high northern latitudes

Up to now, there has been no corroboration from Cassini CIRS of the Voyager IRIS-discovery of cyanoacetylene (HC₃N) ice in Titan’s thermal infrared spectrum. We report the first compelling spectral evidence from CIRS for the ν₆ HC₃N ice feature at 506 cm⁻¹ at latitudes 62°N and 70°N, from which we derive particle sizes and column abundances in Titan’s lower stratosphere. We find mean particle radii of 3.0 μm and 2.3 μm for condensed HC₃N at 62°N and 70°N, respectively, and corresponding ice phase molecular column abundances in the range 1-10 × 10¹⁶ mol cm⁻². Only upper limits for cloud abundances can be established at latitudes of 85°N, 55°N, 30°N, 10°N, and 15°S. Under the assumption that cloud tops coincide with the uppermost levels at which HC₃N vapor saturates, we infer geometric thicknesses for the clouds equivalent to 10-20 km or so, with tops at 165 km and 150 km at 70°N and 62°N, respectively.     

Methane absorption in the atmosphere of Jupiter from 1800 to 9500 cm and implications for vertical cloud structure

New measurements of the low-temperature near-infrared absorption of methane (Sihra, 1998, Laboratory measurements of near-infrared methane bands for remote sensing of the jovian atmosphere, Ph.D. thesis, University of Oxford) have been combined with existing, longer path-length, higher-temperature data of Strong et al. (1993, Spectral parameters of self- and hydrogen-broadened methane from 2000 to 9500 cm⁻¹ for remote sounding of the atmosphere of Jupiter, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) and fitted with band models. The combined data set is found to be more consistent with previous low-temperature methane absorption measurements than that of Strong et al. (1993, J. Quant. Spectrosc. Radiat. Trans. 50, 309-325) but covers the same wider wavelength range and accounts for both self- and hydrogen-broadening conditions. These data have been fitted with k-coefficients in the manner described by Irwin et al. (1996, Calculated k-distribution coefficients for hydrogen- and self-broadened methane in the range 2000-9500 cm⁻¹ from exponential sum fitting to band modelled spectra, J. Geophys. Res. 101, 26,137-26,154) and have been used in multiple-scattering radiative transfer models to assess their impact on our previous estimates of the jovian cloud structure obtained from Galileo Near-Infrared Mapping Spectrometer (NIMS) observations (Irwin et al., 1998, Cloud structure and atmospheric composition of Jupiter retrieved from Galileo NIMS real-time spectra, J. Geophys. Res. 103, 23,001-23,021; Irwin et al., 2001, The origin of belt/zone contrasts in the atmosphere of Jupiter and their correlation with 5-μm opacity, Icarus 149, 397-415; Irwin and Dyudina, 2002, The retrieval of cloud structure maps in the equatorial region of Jupiter using a principal component analysis of Galileo/NIMS data, Icarus 156, 52-63). Although significant differences in methane opacity are found at cooler temperatures, the difference in the optical depth of the atmosphere due to methane is found to diminish rapidly with increasing pressure and temperature and thus has negligible effect on the cloud structure inferred at deeper levels. Hence the main cloud opacity variation is still found to peak at around 1-2 bar using our previous analytical approach, and is thus still in disagreement with Galileo Solid State Imager (SSI) determinations (Banfield et al., 1998, Jupiter's cloud structure from Galileo imaging data, Icarus 135, 230-250; Simon-Miller et al., 2001, Color and the vertical structure in Jupiter's belts, zones and weather systems, Icarus 154, 459-474) which place the main cloud deck near 0.9 bar. Further analysis of our retrievals reveals that this discrepancy is probably due to the different assumptions of the two analyses. Our retrievals use a smooth vertically extended cloud profile while the SSI determinations assume a thin NH₃ cloud below an extended haze. When the main opacity in our model is similarly assumed to be due to a thin cloud below an extended haze, we find the main level of cloud opacity variation to be near the 1 bar level—close to that determined by SSI and moderately close to the expected condensation level of ammonia ice of 0.85 bar, assuming that the abundance of ammonia on Jupiter is (7±1)×¹⁰⁻⁴ (Folkner et al., 1998, Ammonia abundance in Jupiter's atmosphere derived from the attenuation of the Galileo probe's radio signal, J. Geophys. Res. 103, 22,847-22,855; Atreya et al., 1999, A comparison of the atmospheres of Jupiter and Saturn: deep atmospheric composition, cloud structure, vertical mixing, and origin, Planet. Space Sci. 47, 1243-1262). However our data in the 1-2.5 μm range have good height discrimination and our lowest estimate of the cloud base pressure of 1 bar is still too great to be consistent with the most recent estimates of the ammonia abundance of 3.5 × solar. Furthermore the observed limited spatial distribution of ammonia ice absorption features on Jupiter suggests that pure ammonia ice is only present in regions of localised vigorous uplift (Baines et al., 2002, Fresh ammonia ice clouds in Jupiter: spectroscopic identification, spatial distribution, and dynamical implications, Icarus 159, 74-94) and is subsequently rapidly modified in some way which masks its pure absorption features. Hence we conclude that the main cloud deck on Jupiter is unlikely to be composed of pure ammonia ice and instead find that it must be composed of either NH₄SH or some other unknown combination of ammonia, water, and hydrogen sulphide and exists at pressures of between 1 and 2 bar.     

Colliding Blast Waves Driven by the Interaction of a Short-Pulse Laser with a Gas of Atomic Clusters

Collisions between shocks are commonly found in many astrophysical objects, however robust numerical models or laboratory analogues of these complex systems remain challenging to implement. We report on the development of scaled laboratory experiments which employ new techniques for launching and diagnosing colliding shocks and high Mach number blast waves, scalable to a limited subset of astrophysically-relevant regimes. Use of an extended medium of atomic clusters enables efficient (>80%) coupling of 700 fs, 1 J, 1054 nm laser pulses to a “cluster” gas with an average density of ≈10¹⁹ particles cm⁻³, producing an initial energy density >10⁵ J cm⁻³, equivalent to ≈5×10⁹ J/g. Multiple laser foci are used to tailor the spatial profile of energy deposition, or to launch pairs of counter-propagating cylindrical shocks which then collide. By probing the collision interferometrically at multiple view angles in 5^{^} increments and applying an inverse Radon transform to the resulting phase projections we have been able to tomographicall reconstruct the full three-dimensional, time-framed electron density profile of the system.     

Re-analysis of previous laboratory phase curves: 1. Variations of the opposition effect morphology with the textural properties,and an application to planetary surfaces

Typical variations in the opposition effect morphology of laboratory samples at optical wavelengths are investigated to probe the role of the textural properties of the surface (roughness, porosity and grain size). A previously published dataset of 34 laboratory phase curves is re-analyzed and fit with several morphological models. The retrieved morphological parameters that characterize the opposition surge, amplitude, width and slope (A, HWHM and S respectively) are correlated to the single scattering albedo, the roughness, the porosity and the grain size of the samples. To test the universality of the laboratory samples’ trends, we use previously published phase curves of planetary surfaces, including the Moon, satellites and rings of the giant planets. The morphological parameters of the surge (A and HWHM) for planetary surfaces are found to have a non-monotonic variation with the single scattering albedo, similar to that observed in asteroids (Belskaya, I.N., Shevchenko, V.G. [2000]. Icarus 147, 94–105), which is unexplained so far. The morphological parameters of the surge (A and HWHM) for laboratory samples seem to exhibit the same non-monotonic variation with single scattering albedo. While the non-monotonic variation with albedo was already observed by Nelson et al. (Nelson, R.M., Hapke, B.W., Smythe, W.D., Hale, A.S., Piatek, J.L. [2004]. Planetary regolith microstructure: An unexpected opposition effect result. In: Mackwell, S., Stansbery, E. (Eds.), Proc. Lunar Sci. Conf. 35, p. 1089), we report here the same variation for the angular width.     

An investigation of Titan's aerosols using microwave analog measurements and radiative transfer modeling

A combination of laboratory experiments, theoretical modeling, and spacecraft observations is employed to characterize the aerosols in the atmosphere of Titan. The scattering properties of model aerosols were measured using the Microwave Analog Light Scattering Facility at the University of Florida and complemented with theoretical modeling of single scattering characteristics and radiative transfer in Titan's atmosphere. This study compares these modeling results with photopolarimetric observations made over a range of phase angles by the Pioneer 11 and Voyagers 1 and 2 spacecraft. Important results of this work include a survey of the scattering properties of different particle morphologies and compositions necessary to accurately interpret these observations without introducing non-physical assumptions about the particles or requiring additional free parameters to the radiative transfer models. Previous studies use calculation methods which, due to computing memory and processing time requirements, a priori exclude much of the parameter space that the microwave analog laboratory is ideal for exploring. The goal of the present work, to directly constrain aerosol physical characteristics, is addressed by studying in a consistent manner how a variety of particle morphologies and refractive indices affect the polarization and intensity reflected by Titan's atmosphere. Based on comparisons of model results to spacecraft observations, many model morphologies are excluded from further consideration. The most plausible physical particle models suggest that a combination of Rayleigh-like single particles and aggregates that are larger than those previously suggested and investigated [West, R.A., Smith, P.H., 1991. Evidence for aggregate particles in the atmospheres of Titan and Jupiter. Icarus 90, 330-333; Rannou, P., Cabane, M., Botet, R., Chassefière, E., 1997. A new interpretation of scattered light measurements at Titan's limb. J. Geophys. Res. 102, 10997-11013] provide the best fit to the existing data. Additional laboratory experiments and more refined modeling awaits the results of the new rich observational dataset from the Cassini/Huygens encounter with Titan.     

Examining the wake structure in Saturn's rings from microwave observations over varying ring opening angles and wavelengths

Over the last 15 to 20 years several high quality, high resolution data have been taken with the very large array (VLA). These data exhibit a wide range of ring opening angles (|B|=0 to 26°) and wavelengths (λ=0.7 to 20 cm). At these wavelengths the primary flux from the rings is scattered saturnian thermal emission, with a small contribution coming from the ring particles' own thermal emission. Much of the data do show signs of asymmetries due to wakes either on the ansae or the portion of the rings which occult the planet. As in previous work, we use our Monte Carlo radiative transfer code including idealized wakes [Dunn, D.E., Molnar, L.A., Fix, J.D., 2002. Icarus 160, 132-160; Dunn, D.E., Molnar, L.A., Niehof, J.T., de Pater, I., Lissauer, J.L., 2004. Icarus 171, 183-198] to model the relative contributions of the scattered and thermal radiation emanating from the rings and compare the results to that seen in the data. Although the models do give satisfactory fits to all of our data, we find that no single model can simulate the data at all different |B| and λ. We find that one model works best for moderate and low |B| and another one at higher |B|. The main difference between these models is the ratio of the wake width to their separation. We similarly find that the 2 cm data require higher density wakes than the longer wavelength data, perhaps caused by a preponderance of somewhat smaller ring material in the wakes. We further find evidence for an increase in the physical temperature of the rings with increasing |B|. Continuous observations are required to determine whether the above results regarding variations in wake parameters with |B| and λ are indeed caused by these parameters, or instead by changes over time.     

Testing the Hapke photometric model: Improved inversion and the porosity correction

In conjunction with a companion paper (Shepard, M.K., Helfenstein, P. [2011]. Icarus, submitted for publication), we derive, test, and apply a detailed approach for visualizing the phase angle dependence of light scattering in particulate soils from both whole-disk and disk-resolved observations. To reduce the number of model parameters and provide stronger constraints on model fits, we combine Hapke’s (Hapke, B. [2008]. Icarus 195, 918-926) recent correction for effects of porosity with his (Hapke, B. [1986]. Icarus 67, 264-280) model of the shadow hiding opposition effect. We further develop our method as a tool for least-squares fitting of Hapke’s model to photometric data. Finally, we present an improved method for estimating uncertainties in retrieved values of Hapke model parameters. We perform a preliminary test of the model on spectrogoniometric measurements from three selected laboratory samples from Shepard and Helfenstein (Shepard, M.K., Helfenstein, P. [2007]. J. Geophys. Res. 112 (E03001), 17). Our preliminary suite of test samples is too small and selective to permit the drawing of general conclusions. However, our results suggest that Hapke’s porosity correction improves the fidelity of fits to samples composed of low- and moderate-albedo particles and may allow for more reliable retrieval of porosity estimates in these materials. However, we find preliminary evidence that in high-albedo surfaces, the effects of porosity may be difficult to detect.     

Effects of surface roughness and graybody emissivity on martian thermal infrared spectra

Slopes are present in martian apparent surface emissivity observations collected by the Thermal Emission Spectrometer (TES) and the Thermal Emission Imaging System (THEMIS). These slopes are attributed to misrepresenting the surface temperature, either through incorrect assumptions about the maximum emissivity of surface materials or the presumption of a uniform surface temperature within the field of view. These incorrect assumptions leave distinct characteristics in the resulting apparent emissivity data that can be used to gain a better understanding of the surface properties. Surfaces with steep slopes typically have a variable surface temperatures within the field of view that cause distinct and highly variable slopes in apparent emissivity spectra based on the observing conditions. These properties are documented on the southwestern flank of Apollinaris Patera and can be reasonably approximated by modeled data. This spectral behavior is associated with extremely rough martian surfaces and includes surfaces south of Arsia Mons and near Warrego Valles that also appear to have high slopes in high resolution images. Surfaces with low maximum values of emissivity have apparent emissivity spectra with more consistent spectral slopes that do not vary greatly based on observing conditions. This spectral surface type is documented in Terra Serenum and is consistent with associated high resolution images that do not indicate the presence of a surface significantly rougher that the surrounding terrain.     

Interpreting photometry of regolith-like surfaces with different topographies: shadowing and multiple scattering

The effects of various types of topography on the shadow-hiding effect and multiple scattering in particulate surfaces are studied. Two bounding cases were examined: (1) the characteristic scale of the topography is much larger than the surface particle size, and (2) the characteristic scale of the topography is comparable to the surface particle size. A Monte Carlo ray-tracing method (i.e., geometric optics approximation) was used to simulate light scattering. The computer modeling shows that rocky topographies generated by randomly distributed stones over a flat surface reveal much steeper phase curves than surface with random topography generated from Gaussian statistics of heights and slopes. This is because rocks may have surface slopes greater than 90°. Consideration of rocky topography is important for interpreting rover observations. We show the roughness parameter in the Hapke model to be slightly underestimated for bright planetary surfaces, as the model neglects multiple scattering on large-scale topographies. The multiple scattering effect also explains the weak spectral dependences of the roughness parameter in Hapke's model found by some authors. Multiple scattering between different parts of a rough surface suppresses the effect of shadowing, thus the effects produced by increases in albedo on the photometric behavior of a surface can be compensated for with the proper decreases in surface roughness. This defines an effective (photometric) roughness for a surface. The interchangeability of albedo and roughness is shown to be possible with fairly high accuracy for large-scale random topography. For planetary surfaces that have a hierarchically arranged large-scale random topography, predictions made with the Hapke model can significantly differ from real values of roughness. Particulate media with surface borders complicated by Gaussian or clumpy random topographies with characteristic scale comparable to the particle size reveal different photometric behaviors in comparison with particulate surfaces that are flat or the scale of their topographies is much larger than the particle size.     

Hydrodynamical and radiative transfer modeling of meteoroid impacts into Saturn's rings

In a small hypervelocity impact, superheated gas and particles glow brightly with thermal emission for a brief time interval at short wavelengths; this phenomenon is referred to as an impact flash. Over the past decade, impact flashes have been observed on the Moon and in the laboratory in both the IR and visible portions of the spectrum. These phenomena have been used to constrain impactor parameters, such as impact size, velocity and composition. With the arrival of the Cassini spacecraft at Saturn, we embarked on a study of impact flashes in Saturn's rings. We present results on the feasibility of observing impact flashes and therefore estimating the flux of meteoroids impacting Saturn's rings using Cassini's Ultraviolet Imaging Spectrograph (UVIS). Our modeling effort is two-fold. We start by simulating impacts using the CTH hydrodynamical code. Impacts involve an icy ring particle and a serpentine meteoroid, modeled with the ANEOS equation of state. The objects are centimeters to meters in diameter and collide at 30 to 50 km s⁻¹. We then use the resulting temperatures and densities of the impact plumes in a radiative transfer calculation. We calculate bound-free, free-free, electron scattering and negative ion opacities along a line-of-sight through the center of each impact plume. Our model has shown that impact flashes will not be seen with the UVIS because (1) the plumes are optically thick when their central temperatures are high, with photosphere temperatures too cool to emit observable UV flux and (2) when the plumes become optically thin, even the hottest region of the plume is too cool to observe in the UV. This corroborates the lack of UVIS impact flash detections to date. Impact flashes are not likely to be seen by other Cassini instruments because of the short lifetimes of the plumes.     

The Pioneer anomaly and new physics

The Pioneer anomaly is one of the most important problems in modern physics. The observed blueshift of the Doppler signals coming back from the space probes Pioneer 10 and 11 is interpreted as being due to an anomalous acceleration a_p = (8.74 ± 1.33) × 10^–8 cm s^–2 towards the Sun. In this paper the blueshift is explained by the frequency shifts of the receivers. These frequency shifts result from an increase in elementary particle masses in time, the rate of increase being tied up with the present‐day Hubble parameter H₀. The result is that the seeming acceleration a_p is the product of H₀ and the velocity of light. Taking new physics into consideration, this paper presents a new explanation of the Pioneer anomaly based on the assumption that the Universe is eternal and infinite without expansion or contraction (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The effects and characteristics of atmospheric dust during martian global dust storm 2001A

We present retrieved trends in dust optical depth, dust effective radius and surface temperature from our analysis of Mars Global Surveyor Thermal Emission Spectrometer daytime data from global dust storm 2001A, and describe their significance for the martian dust cycle. The dust optical depth becomes correlated with surface pressure during southern spring and summer in years both with and without a global dust storm, indicating that global dust mixing processes are important at those seasons. The correlation is low at other times of the year. We found that the observed decay of optical depths at the later stages of the dust storm match, to first-order, theoretical values of clearing from Stokes–Cunningham fallout of the dust. Zonally averaged effective radius is constant within standard deviation of results (between 1.2 and 2.0 μm, with a global mean for all seasons of 1.7 μm), at all latitudes and seasons except at southern latitudes of 35° and higher around equinoxes in both martian years, where it is larger than average (2–3 μm). The emergence and disappearance of these larger particles correlates with observations of polar cap edge storms at those latitudes. Northern latitude observations under similar conditions did not yield a similar trend of larger average effective radii during the equinoxes. We also report on a linear correlation between daytime surface temperature drop and rise in optical depth during the global dust storm. Global dust storm 2001A produced a significant optical depth and surface temperature change.     

A multilayer model for thermal infrared emission of Saturn's rings: Basic formulation and implications for Earth-based observations

We present our new model for the thermal infrared emission of Saturn's rings based on a multilayer approximation. In our model, (1) the equation of classical radiative transfer is solved directly for both visible and infrared light, (2) the vertical heterogeneity of spin frequencies of ring particles is taken into account, and (3) the heat transport due to particles motion in the vertical and azimuthal directions is taken into account. We adopt a bimodal size distribution, in which rapidly spinning small particles (whose spin periods are shorter than the thermal relaxation time) with large orbital inclinations have spherically symmetric temperatures, whereas non-spinning large particles (conventionally called slow rotators) with small orbital inclinations are heated up only on their illuminated sides. The most important physical parameters, which control ring temperatures, are the albedo in visible light, the fraction of fast rotators (ffast) in the optical depth, and the thermal inertia. In the present paper, we apply the model to Earth-based observations. Our model can well reproduce the observed temperature for all the main rings (A, B, and C rings), although we cannot determine exact values of the physical parameters due to degeneracy among them. Nevertheless, the range of the estimated albedo is limited to 0-0.52±0.05, 0.55±0.07-0.74±0.03, and 0.51±0.07-0.74±0.06 for the C, B, and A rings, respectively. These lower and upper limits are obtained assuming all ring particles to be either fast and slow rotators, respectively. For the C ring, at least some fraction of slow rotators is necessary (ffast?0.9) in order for the fitted albedo to be positive. For the A and B rings, non-zero fraction of fast rotators (ffast?0.1-0.2) is favorable, since the increase of the brightness temperature with increasing solar elevation angle is enhanced with some fraction of fast rotators.     

Visible/near-infrared spectrogoniometric observations and modeling of dust-coated rocks

Interpretations of visible/near-infrared reflectance spectra of Mars are often complicated by the effects of dust coatings that obscure the underlying materials of interest. The ability to separate the spectral reflectance signatures of coatings and substrates requires an understanding of how their individual and combined reflectance properties vary with phase angle. Toward this end, laboratory multispectral observations of rocks coated with different amounts of Mars analog dust were acquired under variable illumination and viewing geometries using the Bloomsburg University Goniometer (BUG). These bidirectional reflectance distribution function (BRDF) data were fit with a two-layer radiative transfer model, which replicated BUG observations of dust-coated basaltic andesite substrates relatively well. Derived single scattering albedo and phase function parameters for the dust were useful in testing the model's ability to derive the spectrum of a “blind” substrate (unknown to the modeler) coated with dust. Subsequent tests were run using subsets of the BUG data restricted by goniometric or coating thickness coverage. Using the entire data set provided the best constraints on model parameters, although some reductions in goniometric coverage could be tolerated without substantial degradation. Predictably, the most thinly coated samples provided the best information on the substrate, whereas the thickest coatings best replicated the dust. Dust zenith optical thickness values ∼0.6-0.8 best constrain the substrate and coating simultaneously, particularly for large ranges of incidence or emission angles. The lack of sufficient angles can be offset by having a greater number and range of coatings thicknesses. Given few angles and thicknesses, few constraints can be placed concurrently on the spectral properties of the coating and substrate.     

The jovian anticyclone BA: III. Aerosol properties and color change

A study of the vertical cloud structure of oval BA and its red color change is presented in this third part of our complete analysis. A large interest in Jupiter’s anticyclone BA was created by its reddening that occurred between 2005 and 2006. In this work we quantify the color change in oval BA by using images taken with the Advanced Camera for Surveys (ACS) onboard the Hubble Space Telescope (HST) in six filters from the near ultraviolet (F250W) to the deep methane band in the near infrared (F892N). Reflectivity changes are noteworthy in nadir viewing geometry at the ultraviolet and blue wavelengths (F250W, F330W and F435W filters) but almost undetectable or inside error bars in the rest of filters (F550M, F658N and F892N). The observed reflectivity variations are discussed in terms of a commonly accepted vertical cloud structure model for jovian anticyclones in order to explore some causes for the color alteration. Our models of the observed reflectivity variation show that the vortex clouds did not change its vertical extension (top pressure) or its optical depth. We find that a change occurred in the absorbing properties of the particles populating the upper aerosols (single scattering albedo and imaginary refractive index). A discussion on the thermo-physical and dynamical properties of the vortex that could be in the origin of the color change is also presented.