The jovian anticyclone BA: I. Motions and interaction with the GRS from observations and non-linear simulations

A study of the dynamics of the second largest anticyclone in Jupiter, Oval BA, and its red colour change that occurred in late 2005 is presented in a three part study. The first part, this paper, deals with its long-term kinematical and dynamical behaviour monitored since its formation in 2000 to September 2008 using ground-based observations archived at the public International Outer Planet Watch (IOPW) database. The vortex changed its zonal drift velocity from 1.8 m s⁻¹ in the period 2000-2002 to 0.8 m s⁻¹ in 2002-2003, and to 2.5 m s⁻¹ since late 2003. It also migrated southwards by 1.0 ± 0.5° in latitude between 2000 and 2004, remaining afterwards at an almost fixed latitude position. During the period 2000-2007, the oval also changed its triangular-like shape to a more symmetrical one. No latitudinal change was found in the months before the development of a red annulus in its interior. The colour change took place in less than 5 months in 2005-2006 and no red colour feature was observed to have been present or entrained by BA months before the annulus development. After detailed examination of the four encounters between BA and GRS that took place during this 9 year period, we did not detect any noticeable change in its drift rate or in apparent structure associated with the encounters at cloud level. Also, the area of BA did not significantly change in this period. Additionally, we found that BA displays a long-term oscillation of ∼160 days in its longitude position with peak to peak amplitude of 1.2°. Numerical experiments using the global circulation model EPIC reproduce accurately the shape, connecting it to its latitude migration, and morphology of the oval and confirm that no strong interaction between BA and the GRS is possible at least in the current situation.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Analysis of Jupiter’s Oval BA: A streamlined approach

We present a novel method of constructing streamlines to derive wind speeds within jovian vortices and demonstrate its application to Oval BA for 2001 pre-reddened Cassini flyby data, 2007 post-reddened New Horizons flyby data, and 1998 Galileo data of precursor Oval DE. Our method, while automated, attempts to combine the advantages of both automated and manual cloud tracking methods. The southern maximum wind speed of Oval BA does not show significant changes between these data sets to within our measurement uncertainty. The northern maximum does appear to have increased in strength during this time interval, which likely correlates with the oval’s return to a symmetric shape. We demonstrate how the use of closed streamlines can provide measurements of vorticity averaged over the encircled area with no a priori assumptions concerning oval shape. We find increased averaged interior vorticity between pre- and post-reddened Oval BA, with the precursor Oval DE occupying a middle value of vorticity between these two.     

Vertical structure of Jupiter's troposphere from nonlinear simulations of long-lived vortices

Numerical simulations of jovian vortices at tropical and temperate latitudes, under different atmospheric conditions, have been performed using the EPIC code [Dowling, T.E., Fisher, A.S., Gierasch, P.J., Harrington, J., LeBeau, R.P., Santori, C.M., 1998. Icarus 132, 221-238] to simulate the high-resolution observations of motions and of the lifetimes presented in a previous work [Legarreta, J., Sánchez-Lavega, A., 2005. Icarus 174, 178-191] and infer the vertical structure of Jupiter's troposphere. We first find that in order to reproduce the longevity and drift rate of the vortices, the Brunt-Väisälä frequency of the atmosphere in the upper troposphere (pressures P∼1 to 7 bar) should have a lower limit value of 5×¹⁰⁻³ s⁻¹, increasing upward up to 1.25×¹⁰⁻² s⁻¹ at pressures P∼0.5 bar (latitudes between 15° and 45° in both hemispheres). Second, the vortices drift also depend on the vertical structure of the zonal wind speed in the same range of altitudes. Simulations of the slowly drifting Southern hemisphere vortices (GRS, White Ovals and anticyclones at 40° S) require a vertically-constant zonal-wind with depth, but Northern hemisphere vortices (cyclonic “barges” and anticyclones at 19, 41 and 45° N) require decreasing winds at a rate of ∼5 m s⁻¹ per scale height. However vortices drifting at a high speed, close to or in the peak of East or West jets and in both hemispheres, require the wind speed slightly increasing with depth, as is the case for the anticyclones at 20° S and at 34° N. We deduce that the maximum absolute vertical shear of the zonal wind from P∼1 bar up to P∼7 bar in these jets is ∼15 m s⁻¹ per scale height. Intense vortices with tangential velocity at their periphery ∼100 m s⁻¹ tend to decay asymptotically to velocities ∼40 to 60 m s⁻¹ with a characteristic time that depends on the vortex intensity and static stability of the atmosphere. The vortices adjust their tangential velocity to the averaged peak to peak velocity of the opposed eastward and westward jets at their boundary. We show through our simulations that large-scale and long-lived vortices whose maximum tangential velocity is ∼100 m s⁻¹ can survive by absorbing smaller intense vortices.     

In situ measurements of particle load and transport in dust devils

In situ (mobile) sampling of 33 natural dust devil vortices reveals very high total suspended particle (TSP) mean values of 296 mg m⁻³ and fine dust loadings (PM10) mean values ranging from 15.1 to 43.8 mg m⁻³ (milligrams per cubic meter). Concurrent three-dimensional wind profiles show mean tangential rotation of 12.3 m s⁻¹ and vertical uplift of 2.7 m s⁻¹ driving mean vertical TSP flux of 1689 mg m⁻³ s⁻¹ and fine particle flux of ∼1.0 to ∼50 mg m⁻³ s⁻¹. Peak PM10 dust loading and flux within the dust column are three times greater than mean values, suggesting previous estimates of dust devil flux might be too high. We find that deflation rates caused by dust devil erosion are ∼2.5-50 μm per year in dust devil active zones on Earth. Similar values are expected for Mars, and may be more significant there where competing erosional mechanisms are less likely.     

Vertical structure of Jupiter’s Oval BA before and after it reddened: What changed?

To constrain the properties of Oval BA before and after it reddened, we use Hubble methane band images from 1994 to 2009 to find that the distribution of upper tropospheric haze atop the oval and its progenitors remained unchanged, with reflectivity variations of less than 10% over this time span. We quantify measurement uncertainties and short-term fluctuations in velocity fields extracted from Cassini and Hubble data, and show that there were no significant changes in the horizontal velocity field of Oval BA in 2000, 2006, and 2009. Based on models of the oval’s dynamics, the static stability of the oval’s surroundings was also unchanged.The vertical extent of the oval did not change, based on the unchanged haze reflectivity and unchanged stratification. Published vortex models require Brunt-Väisälä frequencies of about 0.08 s⁻¹ at the base of the vortex, and we combine this value with a review of prior constraints on the vertically variable static stability in Jupiter’s troposphere to show that the vortex must extend down to the condensation level of water in supersolar abundance.The only observable change was an increase in short-wavelength optical absorption that appeared not at the core of the oval, but in a red annulus. The secondary circulation in the vortex keeps this red annulus warmer than the vortex core. Although the underlying cause of the color change cannot be proven, we explore the idea that the new chromophores in the red annulus may be related to a global or hemispheric temperature change.     

Velocity and vorticity measurements of Jupiter's Great Red Spot using automated cloud feature tracking

We have produced mosaics of the Great Red Spot (GRS) using images taken by the Galileo spacecraft in May 2000, and have measured the winds of the GRS using an automated algorithm that does not require manual cloud tracking. Our technique yields a high-density, regular grid of wind velocity vectors that is advantageous over a limited number of scattered wind vectors that result from manual cloud tracking. The high-velocity collar of the GRS is clearly seen from our velocity vector map, and highest wind velocities are measured to be around 170 m s⁻¹. The high resolution of the mosaics has also enabled us to map turbulent eddies inside the chaotic central region of the GRS, similar to those mapped by Sada et al. [Sada, P.V., Beebe, R.F., Conrath, B.J., 1996. Icarus 119, 311-335]. Using the wind velocity measurements, we computed particle trajectories around the GRS as well as maps of relative and absolute vorticities. We have discovered a narrow ring of cyclonic vorticity that surrounds the main anti-cyclonic high-velocity collar. This narrow ring appears to correspond to a ring surrounding the GRS that is bright in 5 μm [Terrile, R.J., Beebe, R.F., 1979. Science 204, 948-951]. It appears that this cyclonic ring is not a transient feature of the GRS, as we have discovered it in a re-analysis of Galileo data taken in 1996 first analyzed by Vasavada et al. [Vasavada, A.R., and 13 colleagues, 1998. Icarus 135, 265-275]. We also calculate how absolute vorticity changes as a function of latitude along a trajectory around the GRS and compare these measurements to similar ones performed by Dowling and Ingersoll [Dowling, T.E., Ingersoll, A.P., 1988. J. Atmos. Sci. 45, 1380-1396] using Voyager data. We show no dramatic evolution in the structure of the GRS since the Voyager era except for additional evidence for a counter-rotating GRS core, an increase in velocity in the main velocity collar, and an overall decrease in the length of the GRS.     

Dynamics of Jupiter’s equatorial region at cloud top level from Cassini and HST images

We present a study of the equatorial region of Jupiter, between latitudes ∼15°S and ∼15°N, based on Cassini ISS images obtained during the Jupiter flyby at the end of 2000, and HST images acquired in May and July 2008. We examine the structure of the zonal wind profile and report the detection of significant longitudinal variations in the intensity of the 6°N eastward jet, up to 60 m s⁻¹ in Cassini and HST observations. These longitudinal variations are, in the HST case, associated with different cloud morphology. Photometric and radiative transfer analysis of the cloud features used as tracers in HST images show that at most there is only a small height difference, no larger than ∼0.5-1 scale heights, between the slow (∼100 m s⁻¹) and fast (∼150 m s⁻¹) moving features. This suggests that speed variability at 6°N is not dominated by vertical wind shears but instead we propose that Rossby wave activity is the responsible for the zonal variability. Removing this variability, we find that Jupiter’s equatorial jet is actually symmetric relative to equator with two peaks of ∼140-150 m s⁻¹ located at latitudes 6°N and 6°S and at a similar pressure level. We also study the local dynamics of particular equatorial features such as several dark projections associated with 5 μm hot spots and a large, long-lived feature called the White Spot (WS) located at 6°S. Convergent flow at the dark projections appears to be a characteristic which depends on the particular morphology and has only been detected in some cases. The internal flow field in the White Spot indicates that it is a weakly rotating quasi-equatorial anticyclone relative to the ambient meridionally sheared flow.     

Effects of a large convective storm on Saturn's equatorial jet

Hubble Space Telescope observations revealed that Saturn's equatorial jet at the cloud level blows at ∼275 m s⁻¹ today, approximately half the ∼470 m s⁻¹ wind during the Voyager flybys in 1980-1981. Radiative transfer calculations estimate the clouds to be significantly higher today than in 1980. The higher clouds make it difficult to observationally isolate any true slowdown from the vertical wind shear because Voyager and Cassini observations show that the winds become slower with altitude. Here, we test the hypothesis that the large equatorial storm in 1990 called the Great White Spot (GWS) decelerated the equatorial jet. We first use order of magnitude estimates to show: (1) if the GWS triggers vertical momentum redistribution, a minor speed change in the troposphere can lead to a substantial stratospheric wind speed change; (2) storm-triggered turbulent mixing slows a prograde equatorial jet; and (3) a prograde equatorial jet inhibits turbulent mixing in latitude. To test whether a GWS-like large storm decelerates the equatorial jet, we perform numerical experiments using the Explicit Planetary Isentropic Coordinate (EPIC) atmosphere model. Our simulation results are consistent with our order of magnitude predictions. We show that the storm excites waves, and the waves transport westward momentum from the troposphere to the stratosphere and decelerate the equatorial jet by as much as ∼40 m s⁻¹ at the 10-mbar level. However, our results show that the storm's effect is too weak at the cloud levels to halve the jet's speed from ∼470 m s⁻¹. Our results suggest that a combination of higher clouds and a true slowdown is necessary to explain the apparent equatorial jet slowdown. We also analyze the effect of waves on the apparent cloud motions, and show that waves can influence cloud-tracking wind speed measurements.     

Jupiter’s shrinking Great Red Spot and steady Oval BA: Velocity measurements with the ‘Advection Corrected Correlation Image Velocimetry’ automated cloud-tracking method

We show that between 1996 and 2006, the area circumscribed by the high-speed collar of the Great Red Spot (GRS) shrunk by 15%, while the peak velocities within its collar remained constant. This shrinkage indicates a dynamical change in the GRS because the region circumscribed by the collar is nearly coincident with the location of the potential vorticity anomaly of the GRS. It was previously observed that the area of the clouds associated with the GRS has been shrinking. However, the cloud cover of the GRS is not coincident with the location of its potential vorticity anomaly or any other of its known dynamical features. We show that the peak velocities of the Oval BA were nearly the same in 2000, when the Oval was white, and in 2006, when it was red, as were all of the other features of the two velocities fields. To measure temporal changes in the GRS and Oval, we extracted velocities from images taken with Galileo, Cassini, and the Hubble Space Telescope using a new iterative method called Advection Corrected Correlation Image Velocimetry (ACCIV). ACCIV finds correlations over image pairs with 10-h time separations when other automated velocity-extraction methods are limited to time separations of 2 h or less. Typically, ACCIV velocities produced from images separated by 10 h had errors that are 3-6 times smaller than similar velocities extracted from images separated by 2 h or less. ACCIV produces velocity fields containing hundreds of thousands of independent correlation vectors (tie-points). Dense velocity fields are needed to locate the loci of peak velocities and other features.     

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.     

Revisiting the thermal inertia of Iapetus: Clues to the thickness of the dark material

The energy balance at the surface of an airless planetary body is strongly influenced by the bolometric Bond albedo and the surface thermal inertia. Both of these values may be calculated through the application of a thermal model to measured surface temperatures. The accuracy of either, though, increases if the value of the other is better constrained. In this study, we used the improved global bolometric Bond albedo map of Iapetus derived from Cassini VIMS and ISS and Voyager ISS data in conjunction with Cassini CIRS temperature data to reevaluate surface thermal inertia across Iapetus. Results showed the thermal inertia of the dark terrain varies between 11 and 14.8 J m⁻² K⁻¹ s^−1/2 while the light material varies between 15 and 25 J m⁻² K⁻¹ s^−1/2. Using an approximation to the thermal properties of the dark overburden derived from our thermal inertia results, we can implement our thermal model to provide estimates on the dark material thickness, which was found to lie between 7 cm and 16 cm. In order to develop an accurate global thermal model, a weighted function that approximates the surface thermal inertia across Iapetus was developed and verified via our measurements. The global bolometric Bond albedo map, surface thermal inertia map, and the thermal model are then used to synthesize global temperature maps that may be used to study the stability of volatiles.     

Jupiter's cyclones and anticyclones vorticity from Voyager and Galileo images

We have measured the internal velocity field in jovian synoptic-scale cyclones and anticyclones by tracking cloud elements in very high spatial resolution images obtained by the Voyager 1 and 2 (in 1979) and Galileo (1996-2000) spacecrafts. In total we have studied 24 different closed vortices (6 cyclones, 18 anticyclones) spanning a latitude range from ∼60° N to 60° S and with East-West sizes larger than ∼2000 km. The tangential component of the velocity as a function of the distance to the vortex center and position angle is used to retrieve the vorticity field. We find that the velocity increases in all the vortices from a nearly quiescent center to a maximum at the vortex periphery, with a record of about 180 m s⁻¹ for the GRS. The vorticity of cyclones and anticyclones increases in general toward their periphery with absolute values in the range from ∼2-14×¹⁰⁻⁵ s⁻¹. There is a marked tendency to increase the vortices vorticity with their latitude location. However the vorticity does not depend on the vortex size, circulation sense, or ambient background meridional wind shear. The vortex Rossby number ranges from ∼0.2 to 0.5. A study of the interaction between the Great Red Spot with other vortices show that the GRS does not change its vorticity upon their absorption. The two White Ovals mergers showed contradictory results, with greater vorticity in the case of BE, but lower vorticity in the case of BA, although data are poorer for this last case. We present the case of a short lived but large coherent cyclone at −59° that was embedded in a weakly anticyclone wind shear domain. We show that jovian vortices do not follow the simple Kida vortex relationship between vorticity and aspect ratio as it has been previously suggested.     

Dust devil sediment flux on Earth and Mars: Laboratory simulations

Laboratory simulations using the Arizona State University Vortex Generator (ASUVG) were run to simulate sediment flux in dust devils in terrestrial ambient and Mars-analog conditions. The objective of this study was to measure vortex sediment flux in the laboratory to yield estimations of natural dust devils on Earth and Mars, where all parameters may not be measured. These tests used particles ranging from 2 to 2000 μm in diameter and 1300 to 4800 kg m⁻³ in density, and the results were compared with data from natural dust devils on Earth and Mars. Typically, the cores of dust devils (regardless of planetary environment) have a pressure decrease of ∼0.1-1.5% of ambient atmospheric pressure, which enhances the lifting of particles from the surface. Core pressure decreases in our experiments ranged from ∼0.01% to 5.00% of ambient pressure (10 mbar Mars cases and 1000 mbar for Earth cases) corresponding to a few tenths of a millibar for Mars cases and a few millibars for Earth cases. Sediment flux experiments were run at vortex tangential wind velocities of 1-45 m s⁻¹, which typically correspond to ∼30-70% above vortex threshold values for the test particle sizes and densities. Sediment flux was determined by time-averaged measurements of mass loss for a given vortex size. Sediment fluxes of ∼10⁻⁶-10⁰ kg m⁻² s⁻¹ were obtained, similar to estimates and measurements for fluxes in dust devils on Earth and Mars. Sediment flux is closely related to the vortex intensity, which depends on the strength of the pressure decrease in the core (ΔP). This study found vortex size is less important for lifting materials because many different diameters can have the same ΔP. This finding is critical in scaling the laboratory results to natural dust devils that can be several orders of magnitude larger than the laboratory counterparts.     

EPIC simulations of the merger of Jupiter's White Ovals BE and FA: altitude-dependent behavior

Observations of the merger of Jupiter's White Ovals BE and FA show altitude-dependent behavior that we seek to capture in numerical simulations. In particular, it was observed that the upper portions of the vortices orbited each other before merging, but the lower portions translated into each other without orbiting, a phenomenon we term the pair-orbit vertical dichotomy. To reproduce this dichotomy in the EPIC model, it is sufficient to have (i) a decrease with altitude of the background zonal winds above the cloud level with a scale height of ∼2.4, (ii) a height scale of the winds inside the vortices that is the same or larger, and (iii) a maximum tangential velocity in each vortex of ∼100 m s⁻¹ or larger. Condition (i) is expected from thermal-wind analyses, (ii) is consistent with thermal-wind and Shoemaker-Levy 9 debris-trajectory analyses, and (iii) is consistent with cloud-top wind tracking. The model generally does not reproduce the dichotomy for vortices with smaller vertical extent or weaker circulations. Our simulated mergers correctly reproduce the observed ∼9° separation at which vortices start to orbit in the upper layers before they merge and the ∼70% area ratio of the final vortex BA to the sum of BE and FA.     

Microscale simulations of Venus’ convective adjustment and mixing near the surface: Thermal and material transport processes

Heat and material transport processes caused by convective adjustment and mixing are important in modeling of Venus’ atmosphere. In the present study, microscale atmospheric simulations near the venusian surface were conducted using a Weather Research and Forecasting model to elucidate the thermal and material transport processes of convective adjustment and mixing. When convective adjustment occurs, the heat and passive tracer are rapidly mixed into the upper stable layer with convective penetration. The convective adjustment produces large eddy diffusions of heat and passive tracer, which may explain the large eddy diffusions estimated in the radiative-convective equilibrium model.For values of surface heat flux Q greater than a threshold (=0.064 K m s⁻¹ in the present study), the convectively mixed layer with high eddy diffusion coefficients grows with time. In contrast, the mixed layer decays with time for Q values smaller than the threshold. The thermal structure near the surface is controlled not only by extremely long-term radiative processes, but also by microscale dynamics with time scales of several hours. A mixed layer with high eddy diffusion coefficients may be maintained or grow with time if the surface heat flux is high in the volcanic hotspot and adjacent areas.     

In situ observation of penetration process in silica aerogel: Deceleration mechanism of hard spherical projectiles

A large number of cometary dust particles were captured with low-density silica aerogels by NASA’s Stardust Mission. Knowledge of the details of the capture mechanism of hypervelocity particles in silica aerogel is needed in order to correctly derive the original particle features from impact tracks. However, the mechanism has not been fully understood yet. We shot hard spherical projectiles of several different materials into silica aerogel of density 60 mg cm⁻³ and observed their penetration processes using an image converter or a high-speed video camera. In order to observe the deceleration of projectiles clearly, we carried out impact experiments at two velocity ranges; ∼4 km s⁻¹ and ∼200 m s⁻¹. From the movies we took, it was indicated that the projectiles were decelerated by hydrodynamic force which was proportional to v² (v: projectile velocity) during the faster penetration process (∼4 km s⁻¹) and they were merely overcoming the aerogel crushing strength during the slower penetration process (∼200 m s⁻¹). We applied these deceleration mechanisms for whole capture process to calculate the track length. Our model well explains the track length in the experimental data set by Burchell et al. (Burchell, M.J., Creighton, J.A., Cole, M.J., Mann, J., Kearsley, A.T. [2001]. Meteorit. Planet. Sci. 36, 209-221).     

Jovian chromophore characteristics from multispectral HST images

The chromophores responsible for coloring the jovian atmosphere are embedded within Jupiter’s vertical aerosol structure. Sunlight propagates through this vertical distribution of aerosol particles, whose colors are defined by ?₀(λ), and we remotely observe the culmination of the radiative transfer as I/F(λ). In this study, we employed a radiative transfer code to retrieve ?₀(λ) for particles in Jupiter’s tropospheric haze at seven wavelengths in the near-UV and visible regimes. The data consisted of images of the 2008 passage of Oval BA to the south of the Great Red Spot obtained by the Wide Field Planetary Camera 2 on-board the Hubble Space Telescope. We present derived particle colors for locations that were selected from 14 weather regions, which spanned a large range of observed colors. All ?₀(λ) curves were absorbing in the blue, and ?₀(λ) increased monotonically to approximately unity as wavelength increased. We found accurate fits to all ?₀(λ) curves using an empirically derived functional form: ?₀(λ) = 1 − A exp(−Bλ). The best-fit parameters for the mean ?₀(λ) curve were A = 25.4 and B = 0.0149 for λ in units of nm. We performed a principal component analysis (PCA) on our ?₀(λ) results and found that one or two independent chromophores were sufficient to produce the variations in ?₀(λ). A PCA of I/F(λ) for the same jovian locations resulted in principal components (PCs) with roughly the same variances as the ?₀(λ) PCA, but they did not result in a one-to-one mapping of PC amplitudes between the ?₀(λ) PCA and I/F(λ) PCA. We suggest that statistical analyses performed on I/F(λ) image cubes have limited applicability to the characterization of chromophores in the jovian atmosphere due to the sensitivity of I/F(λ) to horizontal variations in the vertical aerosol distribution.     

Vertical profile of H2SO4 vapor at 70-110 km on Venus and some related problems

The vertical profile of H₂SO₄ vapor is calculated using current atmospheric and thermodynamic data. The atmospheric data include the H₂O profiles observed at 70-112 km by the SOIR solar occultations, the SPICAV-UV profiles of the haze extinction at 220 nm, the VeRa temperature profiles, and a typical profile of eddy diffusion. The thermodynamic data are the saturated vapor pressures of H₂O and H₂SO₄ and chemical potentials of these species in sulfuric acid solutions. The calculated concentration of sulfuric acid in the cloud droplets varies from 85% at 70 km to a minimum of 70% at 90 km and then gradually increasing to 90-100% at 110 km. The H₂SO₄ vapor mixing ratio is ∼10⁻¹² at 70 and 110 km with a deep minimum of 3 × 10⁻¹⁸ at 88 km. The H₂O-H₂SO₄ system matches the local thermodynamic equilibrium conditions up to 87 km. The column photolysis rate of H₂SO₄ is 1.6 × 10⁵ cm⁻² s⁻¹ at 70 km and 23 cm⁻² s⁻¹ at 90 km. The calculated abundance of H₂SO₄ vapor at 90-110 km and its photolysis rate are smaller than those presented in the recent model by Zhang et al. (Zhang, X., Liang, M.C., Montmessin, F., Bertaux, J.L., Parkinson, C., Yung, Y.L. [2010]. Nat. Geosci. 3, 834-837) by factors of 10⁶ and 10⁹, respectively. Assumptions of 100% sulfuric acid, local thermodynamic equilibrium, too warm atmosphere, supersaturation of H₂SO₄ (impossible for a source of SO_X), and cross sections for H₂SO₄·H₂O (impossible above the pure H₂SO₄) are the main reasons of this huge difference. Significant differences and contradictions between the SPICAV-UV, SOIR, and ground-based submillimeter observations of SO_X at 70-110 km are briefly discussed and some weaknesses are outlined. The possible source of high altitude SO_X on Venus remains unclear and probably does not exist.     

Titan’s atomic hydrogen corona

Based on measurements performed by the Hydrogen Deuterium Absorption Cell (HDAC) aboard the Cassini orbiter, Titan’s atomic hydrogen exosphere is investigated. Data obtained during the T9 encounter are used to infer the distribution of atomic hydrogen throughout Titan’s exosphere, as well as the exospheric temperature.The measurements performed during the flyby are modeled by performing Monte Carlo radiative transfer calculations of solar Lyman-α radiation, which is resonantly scattered on atomic hydrogen in Titan’s exosphere. Two different atomic hydrogen distribution models are applied to determine the best fitting density profile. One model is a static model that uses the Chamberlain formalism to calculate the distribution of atomic hydrogen throughout the exosphere, whereas the second model is a Particle model, which can also be applied to non-Maxwellian velocity distributions.The density distributions provided by both models are able to fit the measurements although both models differ at the exobase: best fitting exobase atomic hydrogen densities of n_H = (1.5 ± 0.5) × 10⁴ cm⁻³ and n_H = (7 ± 1) × 10⁴ cm⁻³ were found using the density distribution provided by both models, respectively. This is based on the fact that during the encounter, HDAC was sensitive to altitudes above about 3000 km, hence well above the exobase at about 1500 km. Above 3000 km, both models produce densities which are comparable, when taking into account the measurement uncertainty.The inferred exobase density using the Chamberlain profile is a factor of about 2.6 lower than the density obtained from Voyager 1 measurements and much lower than the values inferred from current photochemical models. However, when taking into account the higher solar activity during the Voyager flyby, this is consistent with the Voyager measurements. When using the density profile provided by the particle model, the best fitting exobase density is in perfect agreement with the densities inferred by current photochemical models.Furthermore, a best fitting exospheric temperature of atomic hydrogen in the range of T_H = (150-175) ± 25 K was obtained when assuming an isothermal exosphere for the calculations. The required exospheric temperature depends on the density distribution chosen. This result is within the temperature range determined by different instruments aboard Cassini. The inferred temperature is close to the critical temperature for atomic hydrogen, above which it can escape hydrodynamically after it diffused through the heavier background gas.