Jupiter's zone-belt structure at radio wavelengths: II. Comparison of observations with model atmosphere calculations

Model atmosphere calculations are presented which simulate high-resolution maps of Jupiter's radio emission. They are compared with observations recently obtained at the Very Large Array at 1.3, 2.0, 6.1, and 20.5 cm with resolutions ranging from 0.075 to 0.218 Jovian radii (I. de Pater and J. R. Dickel (1986). Jupiter's zone-belt structure at radio wavelengths. I. Observations. Astrophys. J., in press). The models indicate that ammonia gas is strongly depleted in the upper atmosphere with respect to the solar value both in zones and belts. At very high levels in the atmosphere (P < 0.3−0.5 bar) the gas is undersaturated and distributed uniformly over the planet. In the cloud formation region (0.5 < P < 2 bar), the ammonia depletion is largest in the belts, where it extends down to depths corresponding to 1.8–2 bar. In the zones, the lower ammonia abundances are found down to pressures of 1 bar. Deeper into the Jovian atmosphere, at pressures ≥2.2 bar, the gas is overabundant relative to the solar value by nearly a factor of 2 in both zones and belts. The altitude distribution of the ammonia gas is explained in terms of chemistry, cloud physics, and atmospheric dynamics. The undersaturation at high levels in the atmosphere is attributed to photodissociation of ammonia gas under influence of solar UV photons, coupled with Jupiter's meteorology (up- and downward drafts in the atmosphere). The general depletion of this gas throughout Jupiter's upper atmosphere may be caused by trapping of the gas in a layer of NH₄SH particles, and/or in an aqueous ammonia cloud. The cloud deck responsible for trapping ammonia gas is thicker above zones than belts. If the observed depletion of ammonia gas is entirely due to trapping in an NH₄SH cloud, the difference in thickness of this cloud between zones and belts gives rise to a temperature difference of 3–4°K between the two regions. This temperature difference may trigger the zonal wind motions in Jupiter's atmosphere near the cloud tops.     

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.     

The cloud structure of the jovian atmosphere as seen by the Cassini/CIRS experiment

We analyze the thermal infrared spectra of Jupiter obtained by the Cassini-CIRS instrument during the 2000 flyby to infer temperature and cloud density in the jovian stratosphere and upper troposphere. We use an inversion technique to derive zonal mean vertical profiles of cloud absorption coefficient and optical thickness from a narrow spectral window centered at 1392 cm⁻¹ (7.18 μm). At this wavenumber atmospheric absorption due to ammonia gas is very weak and uncertainties in the ammonia abundance do not impact the cloud retrieval results. For cloud-free conditions the atmospheric transmission is limited by the absorption of molecular hydrogen and methane. The gaseous optical depth of the atmosphere is of order unity at about 1200 mbar. This allows us to probe the structure of the atmosphere through a layer where ammonia cloud formation is expected. The results are presented as height vs latitude cross-sections of the zonal mean cloud optical depth and cloud absorption coefficient. The cloud optical depth and the cloud base pressure exhibit a significant variability with latitude. In regions with thin cloud cover (cloud optical depth less than 2), the cloud absorption coefficient peaks at 1.1±0.05 bar, whereas in regions with thick clouds the peak cloud absorption coefficient occurs in the vicinity of 900±50 mbar. If the cloud optical depth is too large the location of the cloud peak cannot be identified. Based on theoretical expectations for the ammonia condensation pressure we conclude that the detected clouds are probably a system of two different cloud layers: a top ammonia ice layer at about 900 mbar covering only limited latitudes and a second, deeper layer at 1100 mbar, possibly made of ammonium hydrosulfide.     

Jupiter's 24° N highest speed jet: Vertical structure deduced from nonlinear simulations of a large-amplitude natural disturbance

The evolution of a large-amplitude disturbance at cloud level in Jupiter's 24° N jet stream in 1990 is used to constrain the vertical structure of a realistic atmospheric model down to the 6 bar pressure level. We use the EPIC model (Dowling et al., 1998, The explicit planetary isentropic-coordinate (EPIC) atmospheric model, Icarus 132, 221-238) to perform long-term, three-dimensional, nonlinear simulations with a series of systematic variations in vertical structure and find that the details of the 1990 disturbance combine with the characteristics of the 24° N jet, the fastest on Jupiter, to yield a tight constraint on the solution space. The most important free parameters are the vertical dependence of the zonal-wind profile, and the thermal structure, below the cloud tops (p>0.7 bar) at the jet's central latitude. The temporal evolution of the disturbed cloud patterns, which spans more than 2 years, can be reproduced if the jet peak reaches ∼180 ms⁻¹ at the cloud level and increases to ∼210 ms⁻¹ at 1 bar and up to ∼240 ms⁻¹ at 6 bar; the observations were not reproduced for other configurations investigated. This trend is consistent with that measured by the Galileo Probe at 7° N; the implication is that this jovian jet extends well below the solar radiation penetration level situated near the 2 bar level.     

The 2009–2010 fade of Jupiter’s South Equatorial Belt: Vertical cloud structure models and zonal winds from visible imaging

The South Equatorial Belt (SEB) of Jupiter is known to alternate its appearance at visible wavelengths from a classical belt-like band most of the time to a short-lived zone-like aspect which is called a “fade” of the belt, hereafter SEBF. The albedo change of the SEB is due to a change in the structure and properties of the clouds and upper hazes. Recent works based on infrared observations of the last SEBF have shown that the aerosol density below 1 bar increased in parallel with the reflectivity change. However, the nature of the change in the upper clouds and hazes that produces the visible reflectivity change and whether or not this reflectivity change is accompanied by a change in the winds at the upper cloud level remained unknown. In this paper we focus in the near ultraviolet to near infrared reflected sunlight (255–953 nm) to address these two issues. We characterize the vertical cloud structure above the ammonia condensation level from Hubble Space Telescope images, and the zonal wind velocities from long-term high-quality images coming from the International Outer Planet Watch database, both during the SEB and SEBF phases. We show that reflectivity changes do not happen simultaneously in this wavelength range, but they start earlier in the most deep-sensing filters and end in 2010 with just minor changes in those sensing the highest particle layers. Our models require a substantial increase of the optical thickness of the cloud deck at 1.0 ± 0.4 bar from τ_cloud = 6 ± 2 in July 2009 (SEB phase) to semiinfinite at visual wavelengths in 2010 (SEBF). Upper tropospheric particles (～240–1400 mbar) are also required to become substantially reflectant and their single scattering albedo in the blue changes from ?₀ = 0.905 ± 0.005 in November 2009 to ?₀ = 0.95 ± 0.01 in June 2010. No significant changes were found in the cloud top heights or in the particle density of the tropospheric haze. The disturbance travels from the levels below ～3 bar to a level about 400 ± 100 mbar. We derive an upward velocity of 0.15 ± 0.05 cm/s, in agreement with a diffusive process in Jupiter’s upper troposphere requiring a mean eddy coefficient K ～ 8 × 10⁵ cm² s^?1. On the other hand, cloud tracking on the IOPW imaging showed no significant changes in the zonal wind profile between the SEB and SEBF stages. As in other visually huge changes in Jupiter’s cloud morphology and structure, the wind profile remains robust, possibly indicating a deeply rooted dynamical regime.     

A new model of the hydrogen and helium-broadened microwave opacity of ammonia based on extensive laboratory measurements

Close to 2000 laboratory measurements of the microwave opacity and refractivity of gaseous NH₃ in an H₂/He atmosphere have been conducted in the 1.1-20 cm wavelength range (1.5-27 GHz) at pressures from 30 mbar to 12 bar and at temperatures from 184 to 450 K. The mole fraction of NH₃ ranged from 0.06 to 6% with some additional measurements of pure NH₃. The high accuracy of these results have enabled development of a new model for the opacity of NH₃ in a H₂/He atmosphere under jovian conditions. The model employs the Ben-Reuven lineshape applied to the published inversion line center frequencies and intensities of NH₃ (JPL Catalog—[Pickett, H.M., Poynter, R.L., Cohen, E.A., Delitsky, M.L., Pearson, J.C., Müller, H.S.P., 1998. J. Quant. Spectrosc. Radiat. Trans. 60, 883-890]) with empirically-fitted line parameters for H₂ and He broadening, and for the self-broadening of some previously unmeasured ammonia inversion lines. The new model for ammonia opacity will provide reliable results for temperatures from 150 to 500 K, at pressures up to 50 bar and at frequencies up to 40 GHz. These results directly impact the retrieval of jovian atmospheric constituent abundances from the Galileo Probe radio signal absorption measurements, from microwave emission measurements conducted with Earth-based radio telescopes and with the future NASA Juno mission, and studies of Saturn's atmosphere conducted with the Cassini Radio Science Experiment and the Cassini RADAR 2.1 cm passive radiometer.     

Statistical Characteristics on SEPs,Radio-Loud CMEs,Low Frequency Type II and Type III Radio Bursts Associated with Impulsive and Gradual Flares

We have statistically analyzed a set of 115 low frequency (Deca-Hectometer wavelengths range) type II and type III bursts associated with major Solar Energetic Particle (SEP: E_p?>?10 MeV) events and their solar causes such as solar flares and coronal mass ejections (CMEs) observed from 1997 to 2014. We classified them into two sets of events based on the duration of the associated solar flares:75 impulsive flares (duration <?60 min) and 40 gradual flares (duration >?60 min).On an average, the peak flux (integrated flux) of impulsive flares?×?2.9 (0.32 J m^?2) is stronger than that of gradual flares M6.8 (0.24 J m^?2). We found that impulsive flare-associated CMEs are highly decelerated with larger initial acceleration and they achieved their peak speed at lower heights (??27.66 m s^?2 and 14.23 R_o) than the gradual flare-associated CMEs (6.26 m s^?2 and 15.30 R_o), even though both sets of events have similar sky-plane speed (space speed) within LASCO field of view. The impulsive flare-associated SEP events (R_t?=?989.23 min: 2.86 days) are short lived and they quickly reach their peak intensity (shorter rise time) when compared with gradual flares associated events (R_t?=?1275.45 min: 3.34 days). We found a good correlation between the logarithmic peak intensity of all SEPs and properties of CMEs (space speed: cc?=?0.52, SE_cc?=?0.083), and solar flares (log integrated flux: cc?=?0.44, SE_cc?=?0.083). This particular result gives no clear cut distinction between flare-related and CME-related SEP events for this set of major SEP events. We derived the peak intensity, integrated intensity, duration and slope of these bursts from the radio dynamic spectra observed by Wind/WAVES. Most of the properties (peak intensity, integrated intensity and starting frequency) of DH type II bursts associated with impulsive and gradual flare events are found to be similar in magnitudes. Interestingly, we found that impulsive flare-associated DH type III bursts are longer, stronger and faster (31.30 min, 6.43 sfu and 22.49 MHz h^?1) than the gradual flare- associated DH type III bursts (25.08 min, 5.85 sfu and 17.84 MHz h^?1). In addition, we also found a significant correlation between the properties of SEPs and key parameters of DH type III bursts. This result shows a closer association of peak intensity of the SEPs with the properties of DH type III radio bursts than with the properties DH type II radio bursts, atleast for this set of 115 major SEP events.

     

Star forming activities in W75N and DR21 – Are the two regions in collision?

We have observed the massive star formation region W75N in ¹²CO J = 3 ? 2 with KOSMA. The profile of ¹²CO J = 3 ? 2 indicated that besides the 9 km s^?1 component, there is another component of ?3 km s^?1, which is associated with another star formation region, DR21N, located to the north of DR21. We derived the physical and dynamical parameters of the core and high velocity gas associated with the two components separately. Star forming activities were investigated, including outflows and infall analysis. The two regions overlap in space and are not connected in velocity. We found that the cloud–cloud collision scenario may not apply for the DR21/W75N case.     

Vertical structure modeling of Saturn's equatorial region using high spectral resolution imaging

A series of narrow-band images of Saturn was acquired on 7-11 February 2002 with an acousto-optic imaging spectrometer (AImS) at about 160 wavelengths between 500 and 950 nm. Our unique data set with high spectral agility and wide spectral coverage enabled us to extensively study the cloud structure and aerosol properties of Saturn's equatorial region at −10° latitude. Theoretical center-limb profiles based on twelve cloud models were fit to the observations at 23 wavelengths across the 619-, 727-, and 890-nm methane bands. A simultaneous multiwavelength multivariable fitting algorithm was adopted in varying up to 9 free parameters to efficiently explore the vast multidimensional parameter space, and a total of ∼12,000 initial conditions were tested. From the acceptable ranges of the model parameters, we obtained the following major conclusions: (1) the brightening of Saturn's equatorial region observed near 890 nm in February 2002 (I/F∼0.25 at the central meridian) results from high altitudes of a stratospheric haze layer (τ?∼0.05 above ∼0.04-bar level) and an upper tropospheric cloud (τ∼6 above ∼0.25-bar level), (2) if the upper tropospheric cloud is composed of ammonia ice particles and the Mie theory is applied, the mean particle size is larger than about 0.5 μm, (3) an optically thick cloud layer exists at a level of 0.5-2.2 bar below the upper cloud deck in Saturn's equatorial region. The ongoing observations by the Cassini spacecraft over wider spectral range and from various phase angles will further constrain Saturn's cloud structure and aerosol properties.     

Long-term monitoring of the long-period variable Y Cassiopeiae in the 1.35-cm water-vapor line

Observations of the circumstellar maser emission from the long-period variable star Y Cas in the 1.35-cm water-vapor line are presented. The observations were performed with the RT-22 radio telescope at the Pushchino Radio Astronomy Observatory (Astrospace Center, Lebedev Physical Institute, Russian Academy of Sciences) in the period 1982–2005. The variations in the integrated flux F_int in the H₂O line correlate with the visual light curve of the star. The phase delay Δ? between the F_int variations and the light curve is 0.2–0.4P (P is the period of the star). The H₂O maser Y Cas belongs to transient sources: peaks of high maser activity alternate with intervals of a low emission level when the H₂O-line flux does not exceed (0.1–0.5) × 10^?20 W m^?2. A “superperiod” of ～5.7 yr was found in the occurrence of activity peaks. A particularly strong maximum of maser radio emission took place at the end of 1997, when the flux F_int reached 15.6 × 10^?20 W m^?2. A model for the H₂O maser variability in Y Cas is discussed. The variability is caused by a periodic action of shock waves driven by stellar pulsations. The H₂O maser flares may be associated with short-lived episodes of enhanced mass loss by the star or with the propagation of a particularly strong shock wave when a planet orbiting the star passes through its periastron.     

Nuclear spin temperature of ammonia in Comet 9P/Tempel 1 before and after the Deep Impact event

The Deep Impact mission succeeded in excavating inner materials from the nucleus of Comet 9P/Tempel 1 on 2005 July 04 (at 05:52 UT). Comet 9P/Tempel 1 is one of Jupiter family short period comets, which might originate in the Kuiper belt region in the solar nebula. In order to characterize the comet and to support the mission from the ground-based observatory, optical high-dispersion spectroscopic observations were carried out with the echelle spectrograph (UVES) mounted on the 8-m telescope VLT (UT2) before and after the Deep Impact event. Ortho-to-para abundance ratios (OPRs) of cometary ammonia were determined from the NH₂ emission spectra. The OPRs of ammonia on July 3.996 UT and 4.997 UT were derived to be 1.28±0.07 (nuclear spin temperature: Tspin=24±2 K) and 1.26±0.08 (Tspin=25±2 K), respectively. There is no significant change between before and after the impact. Actually, most materials ejected from the impact site could have moved away from the nucleus on July 4.997 UT, about 17 h after the impact. However, a small fraction of the ejected materials might remain in the slit of UVES instrument at that time because an excess of about 20% in the NH₂ emission flux is observed above the normal activity level was found [Manfroid, J., Hutsemékers, D., Jehin, E., Cochran, A.L., Arpigny, C., Jackson, W.M., Meech, K.J., Schulz, R., Zucconi, J.-M., 2007. Icarus. This issue]. If the excess of NH₂ on July 04.997 UT was produced from icy materials excavated by the Deep Impact, then an upper-limit of the ammonia OPR would be 1.75 (Tspin>17 K) for those materials. On the other hand, the OPR of ammonia produced from the quiescent sources was similar to that of the Oort cloud comets observed so far. This fact may imply that physical conditions where cometary ices formed were similar between Comet 9P/Tempel 1 and the Oort cloud comets.     

Solar radiation scattered in the Venus atmosphere: The Venera 11,12 data

Results of the scattered solar radiation spectrum measurements made deep in the Venus atmosphere by the Venera 11 and 12 descent probes are presented. The instrument had two channels: spectrometric (to measure downward radiation in the range 0.45 < γ < 1.17 μm) and photometric (four filters and circular angle scanning in an almost vertical plane). Spectra and angular scans were made in the height range from 63 km above the planet surface. The integral flux of solar radiation is 90 ± 12 W m^?2 measured on the surface at the subsolar point. The mean value of surface absorbed radiation flux per planetary unit area is 17.5 ± 2.3 W m^?2. For Venera 11 and 12 landing sites the atmospheric absorbed radiation flux is ～15 W m^?2 for H >; 43 km and ～45 W m^?2 for H < 48 km in the range 0.45 to 1.55 μm. At the landing sites of the two probes the investigated portion of the cloud layer has almost the same structure: it consists of three parts with boundaries between them at about 51 and 57 km. The base of clouds is near 48 km above the surface. The optical depth of the cloud layer (below 63 km) in the range 0.5 to 1 μm does not depend on the wavelength and is ～29 and ～38 for the Venera 11 and 12 landing sites, respectively. The single-scattering albedo, ω₀, in the clouds is very close to 1 outside the absorption bands. Below 58 km the parameter (1 ? ω₀) is <10^?3 for 0.49 and 0.7 μm. The parameter (1 ? ω₀) obviously increases above 60 km. Below 48 km some aerosol is present. The optical depth here is a strong function of wavelength. It varies from 1.5 to 3 at λ = 0.49 μm and from 0.13 to 0.4 at 1.0 μm. The mean size of particles below the cloud deck is about 0.1 μm. Below 35 km true absorption was found at λ < 0.55 μm with the (1 ? ω₀) maximum at H ≈ 15 km. The wavelength and height dependence of the absorption coefficient are compatible with the assumption that sulfur with a mixing ratio ～2 × 10^?8 normalized to S₂ molecules is the absorber. The upper limits of the mixing ratio for Cl₂, Br₂, and NO₂ are 4 × 10^?8, 2 × 10^?11, and 4 × 10^?10, respectively. The CO₂ and H₂O bands are confidently identified in the observed spectra. The mean value of the H₂O mixing ratio is 3 × 10^?5 < F_H2O < 10^?4 in the undercloud atmosphere. The H₂O mixing ratio evidently varies with height. The most probable profile is characterized by a gradual increase from F_H2O = 2 × 10^?5 near the surface to a 10 to 20 times higher value in the clouds.     

Coupling of acoustic waves to clouds in the jovian troposphere

Seismology is the best tool for investigating the interior structure of stars and giant planets. This paper deals with a photometric study of jovian global oscillations. The propagation of acoustic waves in the jovian troposphere is revisited in order to estimate their effects on the planetary albedo. According to the standard model of the jovian cloud structure there are three major ice cloud layers (e.g., [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]). We consider only the highest layers, composed of ammonia ice, in the region where acoustic waves are trapped in Jupiter's atmosphere. For a vertical wave propagating in a plane parallel atmosphere with an ammonia ice cloud layer, we calculate first the relative variations of the reflected solar flux due to the smooth oscillations at about the ppm level. We then determine the phase transitions induced by the seismic waves in the clouds. These phase changes, linked to ice particle growth, are limited by kinetics. A Mie model [Mishchenko et al., 2002. Scattering, Absorption, and Emission of Light by Small Particles. Cambridge Univ. Press, Cambridge, pp. 158-190] coupled with a simple radiation transfer model allows us to estimate that the albedo fluctuations of the cloud perturbed by a seismic wave reach relative variations of 70 ppm for a 3-mHz wave. This albedo fluctuation is amplified by a factor of ∼70 relative to the previously published estimates that exclude the effect of the wave on cloud properties. Our computed amplifications imply that jovian oscillations can be detected with very precise photometry, as proposed by the microsatellite JOVIS project, which is dedicated to photometric seismology [Mosser et al., 2004. JOVIS: A microsatellite dedicated to the seismic analysis of Jupiter. In: Combes, F., Barret, D., Contini, T., Meynadier, F., Pagani, L. (Eds.), SF2A-2004, Semaine de l'Astrophysique Francaise, Les Ulis. In: EdP-Sciences Conference Series, pp. 257-258].     

Fresh Ammonia Ice Clouds in Jupiter: I. Spectroscopic Identification, Spatial Distribution, and Dynamical Implications

We report the first spectroscopic detection of discrete ammonia ice clouds in the atmosphere of Jupiter, as discovered utilizing the Galileo Near-Infrared Mapping Spectrometer (NIMS). Spectrally identifiable ammonia clouds (SIACs) cover less than 1% of the globe, as measured in complete global imagery obtained in September 1996 during Galileo's second orbit. More than half of the most spectrally prominent SIACs reside within a small latitudinal band, extending from 2° to 7° N latitude, just south of the 5-μm hot spots. The most prominent of these are spatially correlated with nearby 5-μm-bright hot spots lying 1.5°-3.0° of latitude to the north: they reside over a small range of relative longitudes on the eastward side of hot spots, about 37% of the longitudinal distance to the next hot spot to the east. This strong correlation between the positions of hot spots and the most prominent equatorial SIACs suggests that they are linked by a common planetary wave. Good agreement is demonstrated between regions of condensation predicted by the Rossby wave model of A. J. Friedson and G. S. Orton (1999, Bull. Am. Astron. Assoc31, 1155-1156) and the observed longitudinal positions of fresh ammonia clouds relative to 5-μm hot spots. Consistency is also demonstrated between (1) the lifetime of particles as determined by the wave phase speed and cloud width and (2) the sedimentation time for 10-μm radius particles consistent with previously reported ammonia particle size by T. Y. Brooke et al. (1998, Icarus136, 1-13). A young age (<two days) for most SIAC cloud particles is indicated. To the south, the most prominent SIACs are located to the northwest of the Great Red Spot, in a region where a westward flow of jovian air, diverted approximately 10° of latitude northward by the Great Red Spot, encounters a large eastward flow. SIACs have been observed repeatedly by NIMS at this location during Galileo's first four years in Jupiter orbit. It is speculated that due to the three-dimensional interactions of these flows, relatively large amounts of ammonia gas are steadily transported from the sub-cloud troposphere (below the ∼600-mbar level) to the high troposphere, nearly continuously forming fresh ammonia ice clouds to the northwest of the Great Red Spot.     

Physical properties of asteroid dust bands and their sources

Disruptive collisions in the main belt can liberate fragments from parent bodies ranging in size from several micrometers to tens of kilometers in diameter. These debris bodies group at initially similar orbital locations. Most asteroid-sized fragments remain at these locations and are presently observed as asteroid families. Small debris particles are quickly removed by Poynting-Robertson drag or comminution but their populations are replenished in the source locations by collisional cascade. Observations from the Infrared Astronomical Satellite (IRAS) showed that particles from particular families have thermal radiation signatures that appear as band pairs of infrared emission at roughly constant latitudes both above and below the Solar System plane. Here we apply a new physical model capable of linking the IRAS dust bands to families with characteristic inclinations. We use our results to constrain the physical properties of IRAS dust bands and their source families. Our results indicate that two prominent IRAS bands at inclinations ≈2.1° and ≈9.3° are byproducts of recent asteroid disruption events. The former is associated with a disruption of a ≈30-km asteroid occurring 5.8 Myr ago; this event gave birth to the Karin family. The latter came from the breakup of a large >100-km-diameter asteroid 8.3 Myr ago that produced the Veritas family. Using an N-body code, we tracked the dynamical evolution of ≈10⁶ particles, 1 μm to 1 cm in diameter, from both families. We then used these results in a Monte Carlo code to determine how small particles from each population undergo collisional evolution. By computing the thermal emission of particles, we were able to compare our results with IRAS observations. Our best-fit model results suggest the Karin and Veritas family particles contribute by 5-9% in 10-60-μm wavelengths to the zodiacal cloud's brightness within 50° latitudes around the ecliptic, and by 9-15% within 10° latitudes. The high brightness of the zodiacal cloud at large latitudes suggests that it is mainly produced by particles with higher inclinations than what would be expected for asteroidal particles produced by sources in the main belt. From these results, we infer that asteroidal dust represents a smaller fraction of the zodiacal cloud than previously thought. We estimate that the total mass accreted by the Earth in Karin and Veritas particles with diameters 20-400 μm is ≈15,000-20,000 tons per year (assuming 2 g cm⁻³ particles density). This is ≈30-50% of the terrestrial accretion rate of cosmic material measured by the Long Duration Exposure Facility. We hypothesize that up to ≈50% of our collected interplanetary dust particles and micrometeorites may be made up of particle species from the Veritas and Karin families. The Karin family IDPs should be about as abundant as Veritas family IDPs though this ratio may change if the contribution of third, near-ecliptic source is significant. Other sources of dust and/or large impact speeds must be invoked to explain the remaining ≈50-70%. The disproportional contribution of Karin/Veritas particles to the zodiacal cloud (only 5-9%) and to the terrestrial accretion rate (30-50%) suggests that the effects of gravitational focusing by the Earth enhance the accretion rate of Karin/Veritas particles relative to those in the background zodiacal cloud. From this result and from the latitudinal brightness of the zodiacal cloud, we infer that the zodiacal cloud emission may be dominated by high-speed cometary particles, while the terrestrial impactor flux contains a major contribution from asteroidal sources. Collisions and Poynting-Robertson drift produce the size-frequency distribution (SFD) of Karin and Veritas particles that becomes increasingly steeper closer to the Sun. At 1 AU, the SFD is relatively shallow for small particle diameters D (differential slope exponent of particles with D?100 μm is ≈2.2-2.5) and steep for D?100 μm. Most of the mass at 1 AU, as well as most of the cross-sectional area, is contributed by particles with D≈100-200 μm. Similar result has been found previously for the SFD of the zodiacal cloud particles at 1 AU. The fact that the SFD of Karin/Veritas particles is similar to that of the zodiacal cloud suggests that similar processes shaped these particle populations. We estimate that there are ≈5×¹⁰²⁴ Karin and ≈10²⁵ Veritas family particles with D>30 μm in the Solar System today. The IRAS observation of the dust bands may be satisfactorily modeled using ‘averaged’ SFDs that are constant with semimajor axis. These SFDs are best described by a broken power-law function with differential power index α≈2.1-2.4 for D?100 μm and by α?3.5 for 100 μm?D?1 cm. The total cross-sectional surface area of Veritas particles is a factor of ≈2 larger than the surface area of the particles producing the inner dust bands. The total volumes in Karin and Veritas family particles with 1 μm<D<1 cm correspond to D=11 km and D=14 km asteroids with equivalent masses ≈1.5×¹⁰¹⁸ g and ≈3.0×¹⁰¹⁸ g, respectively (assuming 2 g cm⁻³ bulk density). If the size-frequency and radial distribution of particles in the zodiacal cloud were similar to those in the asteroid dust bands, we estimate that the zodiacal cloud represents ∼3×¹⁰¹⁹ g of material (in particles with 1 μm<D<1 cm) at ±10° around the ecliptic and perhaps as much as ∼¹⁰²⁰ g in total. The later number corresponds to about a 23-km-radius sphere with 2 g cm⁻³ density.     

On the structure of Mars' atmosphere at 120–220 km

Altitude dependences of [CO₂] and [CO₂⁺] are deduced from Mariner 6 and 7 CO₂⁺ airglow measurements. CO₂ densities are also obtained from n_e radio occultation measurements. Both [CO₂] profiles are similar and correspond to the model atmosphere of Barth et al. (1972) at 120 km, but at higher altitudes they diverge and at 200–220 km the obtained [CO₂] values are three times less the model. Both the airglow and radio occultation observations show that a correction factor of 2.5 should be included into the values for solar ionization flux given by Hinteregger (1970). The ratio of [CO₂⁺]/n_e is 0.15–0.2 and, hence, [O]/[CO₂] is ～3% at 135 km. An atmospheric and ionospheric model is developed for 120–220 km. The calculated temperature profile is characterized by a value of T ≈ 370°K at h ? 220 km, a steep gradient (～2°/km) at 200-160 km, a bend in the profile at 160 km, a small gradient (～0.7°/km) below and a value of T ≈ 250°K at 120 km. The upper point agrees well with the results of the Lyman-α measurements; the steep gradient may be explained by molecular viscosity dissipation of gravity and acoustical waves (the corresponding energy flux is 4 × 10^?2 erg cm^?2sec^?1 at 180 km). The bend at 160 km may be caused by a sharp decrease of the eddy diffusion coefficient and defines K ≈ 2 × 10⁸cm²sec^?1; and the low gradient gives an estimate of the efficiency of the atmosphere heating by the solar radiation as ? ≈ 0.1.     

Chandra X-ray observations of two unusual BAL quasars

We report sensitive Chandra X-ray non-detections of two unusual, luminous Iron Low-Ionization Broad Absorption Line Quasars (FeLoBALs). The observations do detect a non-BAL, wide-binary companion quasar to one of the FeLoBAL quasars. We combine X-ray-derived column density lower limits (assuming solar metallicity) with column densities measured from ultraviolet spectra and CLOUDY photoionization simulations to explore whether constant-density slabs at broad-line region densities can match the physical parameters of these two BAL outflows, and find that they cannot. In the “overlapping-trough” object SDSS J0300+0048, we measure the column density of the X-ray absorbing gas to be N_H ? 1.8 × 10²⁴ cm^?2. From the presence of Fe ii UV78 absorption but lack of Fe ii UV195/UV196 absorption, we infer the density in that part of the absorbing region to be n_e ? 10⁶ cm^?3. We do find that a slab of gas at that density might be able to explain this object’s absorption. In the Fe iii-dominant object SDSS J2215–0045, the X-ray absorbing column density of N_H ? 3.4 × 10²⁴ cm^?2 is consistent with the Fe iii-derived N_H ? 2 × 10²² cm^?2 provided the ionization parameter is log U > 1.0 for both the n_e = 10¹¹ cm^?3 and n_e = 10¹² cm^?3 scenarios considered (such densities are required to produce Fe iii absorption without Fe iiabsorption). However, the velocity width of the absorption rules out its being concentrated in a single slab at these densities. Instead, this object’s spectrum can be explained by a low density, high ionization and high temperature disk wind that encounters and ablates higher density, lower ionization Fe iii-emitting clumps.     

Neptune's far-infrared spectrum from the ISO long-wavelength and short-wavelength spectrometers

Neptune was observed by the Infrared Space Observatory (ISO) Long-Wavelength Spectrometer (LWS) between 46 and 185 μm. At wavelengths between 50 and 110 μm the accuracy of these measurements is ?0.3 K. Observations of this planet made by the ISO Short-Wavelength Spectrometer between 28 and 44 μm were combined with the LWS data to determine a disk-averaged temperature profile and derive several physical quantities. The combined spectra are matched best by a He/(H₂+He) mass ratio of 26.4^+2.6_−3.5%, reflecting a He molar fraction of 14.9^+1.7_−2.2%, assuming the molar fraction of CH₄ to be 2% in the troposphere. This He abundance is consistent with one derived from analysis of joint Voyager-2 IRIS and radio occultation experiment data, a technique whose accuracy has recently been called into question. For a disk average, the para-H₂ fraction is found to be no more than ∼1.5% different from its equilibrium value, and the N₂ mixing ratio is probably less than 0.7%. The composite spectrum is best fit by invoking a CH₄ ice condensate cloud. Using a Mie approximation to particle scattering and absorption, best-fit particle sizes lie between 15 and 40 μm. The composite spectra are relatively insensitive to the vertical distribution of the cloud, but the particle scale height must be greater than 5% of the gas scale height. The best models are consistent with an effective temperature for Neptune that is 59.5±0.6 K, a value slightly lower than derived by the Voyager IRIS experiment—possibly Neptune's mid- and far-infrared emission has changed during the seven years that lie between its encounter with Voyager 2 and the first spectra taken of this planet with ISO. The model spectra are also ostensibly lower than ground-based observations in the spectral range of 17-24 μm, but this discrepancy can be relieved by perturbing the temperature of the lower stratosphere where the LWS spectrum is not particularly sensitive, combined with the uncertainty in the absolute calibration of the ground-based measurements.     

Assessing the long-term variability of Venus winds at cloud level from VIRTIS–Venus Express

The Venus Express (VEX) mission has been in orbit to Venus for more than 4 years now. The Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) instrument onboard VEX observes Venus in two channels (visible and infrared) obtaining spectra and multi-wavelength images of the planet that can be used to sample the atmosphere at different altitudes. Day-side images in the ultraviolet range (380 nm) are used to study the dynamics of the upper cloud at 66–72 km while night-side images in the near infrared (1.74 μm) map the opacity of the lower cloud deck at 44–48 km. Here we present a long-term analysis of the global atmospheric dynamics at these levels using a large selection of orbits from the VIRTIS-M dataset covering 860 Earth days that extends our previous work (Sánchez-Lavega, A. et al. [2008]. Geophys. Res. Lett. 35, L13204) and allows studying the variability of the global circulation at the two altitude levels. The atmospheric superrotation is evident with equatorial to mid-latitudes westward velocities of 100 and 60 m s^?1 in the upper and lower cloud layers. These zonal velocities are almost constant in latitude from the equator to 50°S. From 50°S to 90°S the zonal winds at both cloud layers decrease steadily to zero at the pole. Individual cloud tracked winds have errors of 3–10 m s^?1 with a mean of 5 m s^?1 and the standard deviations for a given latitude of our zonal and meridional winds are 9 m s^?1. The zonal winds in the upper cloud change with the local time in a way that can be interpreted in terms of a solar tide. The zonal winds in the lower cloud are stable at mid-latitudes to the tropics and present variability at subpolar latitudes apparently linked to the activity of the South polar vortex. While the upper cloud presents a net meridional motion consistent with the upper branch of a Hadley cell with peak velocity v = 10 m s^?1 at 50°S, the lower cloud meridional motions are less organized with some cloud features moving with intense northwards and southwards motions up to v = ±15 m s^?1 but, on average, with almost null global meridional motions at all latitudes. We also examine the long-term behavior of the winds at these two vertical layers by comparing our extended wind tracked data with results from previous missions.     

Densities Probed by Coronal Type III Radio Burst Imaging

We present coronal density profiles derived from low-frequency (80?–?240 MHz) imaging of three Type III solar radio bursts observed at the limb by the Murchison Widefield Array (MWA). Each event is associated with a white-light streamer at larger heights and is plausibly associated with thin extreme-ultraviolet rays at lower heights. Assuming harmonic plasma emission, we find average electron densities of 1.8\(\times10^{8}\) cm^?3 down to 0.20\(\times10^{8}\) cm^?3 at heights of 1.3 to 1.9 R_⊙. These values represent approximately 2.4?–?5.4× enhancements over canonical background levels and are comparable to the highest streamer densities obtained from data at other wavelengths. Assuming fundamental emission instead would increase the densities by a factor of four. High densities inferred from Type III source heights can be explained by assuming that the exciting electron beams travel along overdense fibers or by radio propagation effects that may cause a source to appear at a larger height than the true emission site. We review the arguments for both scenarios in light of recent results. We compare the extent of the quiescent corona to model predictions to estimate the impact of propagation effects, which we conclude can only partially explain the apparent density enhancements. Finally, we use the time- and frequency-varying source positions to estimate electron beam speeds of between 0.24 and 0.60 c.