Atomic oxygen on the Venus nightside: Global distribution deduced from airglow mapping

The Visible and Infra-Red Thermal Imaging Spectrometer (VIRTIS) instrument on board the Venus Express spacecraft has measured the O₂(a¹Δ) nightglow distribution at 1.27 μm in the Venus mesosphere for more than two years. Nadir observations have been used to create a statistical map of the emission on Venus nightside. It appears that the statistical 1.6 MR maximum of the emission is located around the antisolar point. Limb observations provide information on the altitude and on the shape of the emission layer. We combine nadir observations essentially covering the southern hemisphere, corrected for the thermal emission of the lower atmosphere, with limb profiles of the northern hemisphere to generate a global map of the Venus nightside emission at 1.27 μm. Given all the O₂(a¹Δ) intensity profiles, O₂(a¹Δ) and O density profiles have been calculated and three-dimensional maps of metastable molecular and atomic oxygen densities have been generated. This global O density nightside distribution improves that available from the VTS3 model, which was based on measurements made above 145 km. The O₂(a¹Δ) hemispheric average density is 2.1 × 10⁹ cm^?3, with a maximum value of 6.5 × 10⁹ cm^?3 at 99.2 km. The O density profiles have been derived from the nightglow data using CO₂ profiles from the empirical VTS3 model or from SPICAV stellar occultations. The O hemispheric average density is 1.9 × 10¹¹ cm^?3 in both cases, with a mean altitude of the peak located at 106.1 km and 103.4 km, respectively. These results tend to confirm the modeled values of 2.8 × 10¹¹ cm^?3 at 104 km and 2.0 × 10¹¹ cm^?3 at 110 km obtained by Brecht et al. [Brecht, A., Bougher, S.W., Gérard, J.-C., Parkinson, C.D., Rafkin, S., Foster, B., 2011a. J. Geophys. Res., in press] and Krasnopolsky [Krasnopolsky, V.A., 2010. Icarus 207, 17–27], respectively. Comparing the oxygen density map derived from the O₂(a¹Δ) nightglow observations, it appears that the morphology is very different and that the densities obtained in this study are about three times higher than those predicted by the VTS3 model.     

Atomic oxygen distributions in the Venus thermosphere: Comparisons between Venus Express observations and global model simulations

Nightglow emissions provide insight into the global thermospheric circulation, specifically in the transition region (～70–120 km). The O₂ IR nightglow statistical map created from Venus Express (VEx) Visible and InfraRed Thermal Imaging Spectrometer (VIRTIS) observations has been used to deduce a three-dimensional atomic oxygen density map. In this study, the National Center of Atmospheric Research (NCAR) Venus Thermospheric General Circulation Model (VTGCM) is utilized to provide a self-consistent global view of the atomic oxygen density distribution. More specifically, the VTGCM reproduces a 2D nightside atomic oxygen density map and vertical profiles across the nightside, which are compared to the VEx atomic oxygen density map. Both the simulated map and vertical profiles are in close agreement with VEx observations within a ～30° contour of the anti-solar point. The quality of agreement decreases past ～30°. This discrepancy implies the employment of Rayleigh friction within the VTGCM may be an over-simplification for representing wave drag effects on the local time variation of global winds. Nevertheless, the simulated atomic oxygen vertical profiles are comparable with the VEx profiles above 90 km, which is consistent with similar O₂ (¹Δ) IR nightglow intensities. The VTGCM simulations demonstrate the importance of low altitude trace species as a loss for atomic oxygen below 95 km. The agreement between simulations and observations provides confidence in the validity of the simulated mean global thermospheric circulation pattern in the lower thermosphere.     

Spatial correlation of OH Meinel and O2 infrared atmospheric nightglow emissions observed with VIRTIS-M on board Venus Express

We present the two-dimensional distribution of the O₂ a¹Δ–X³Σ (0–0) band at 1.27 μm and the OH Δv = 1 Meinel airglow measured simultaneously with the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) on board Venus Express. We show that the two emissions present very similar spatial structures. A cross-correlation analysis indicates that the highest level of correlation is reached with only very small relative shifts of the pairs of images. In spite of the strong spatial correlation between the morphology of the bright spots in the two emissions, we also show that their relative intensity is not constant, in agreement with earlier statistical studies of their limb profiles. We conclude that the two emissions have a common precursor that controls the production of both excited species. We argue that atomic oxygen, which produces O₂ (¹Δ) molecules by three-body recombination and is the precursor of ozone formation, also governs to a large extent the OH airglow morphology through the H + O₃ → OH^* + O₂ reaction.     

Characteristics of planetary-scale waves simulated by a new venusian mesosphere and thermosphere general circulation model

We have developed a new general circulation model (GCM) for the venusian mesosphere and thermosphere (80-about 180 km). Our GCM simulations show that winds in the subsolar-to-antisolar direction (SS–AS) are predominant above about 90 km. A weak return flow of the SS–AS is seen below about 90 km. We performed GCM simulations imposing the planetary-scale waves (thermal tides, Rossby wave, and Kelvin wave) at the lower boundary. Although the diurnal and semidiurnal tides are damped below 95 km, the Rossby wave propagates up to around 130 km. However, the amplitude of the Rossby wave is too small (<1 m/s) to affect the general circulation. On the other hand, the Kelvin wave propagates up to about 130 km with a maximum zonal wind fluctuation of approximately 5.9 m/s on average. The amplitude of the Kelvin wave sometimes exceeds 10 m/s around the terminator. The Kelvin wave causes a temporal variation in the wind velocity at the altitude of the O₂-1.27 μm nightglow emission (about 95 km). Using a newly developed 1-D nightglow model and the composition distribution calculated from our GCM, we investigated the impact of the Kelvin wave on the nightglow distribution. Our results suggest that the Kelvin wave would cause temporal variations in the nightglow emission in the 23:50–00:20 LT region with an intensity of 1.1–1.3 MR and a period of approximately 4 days.     

Sulfur chemistry in the middle atmosphere of Venus

Venus Express measurements of the vertical profiles of SO and SO₂ in the middle atmosphere of Venus provide an opportunity to revisit the sulfur chemistry above the middle cloud tops (～58 km). A one dimensional photochemistry-diffusion model is used to simulate the behavior of the whole chemical system including oxygen-, hydrogen-, chlorine-, sulfur-, and nitrogen-bearing species. A sulfur source is required to explain the SO₂ inversion layer above 80 km. The evaporation of the aerosols composed of sulfuric acid (model A) or polysulfur (model B) above 90 km could provide the sulfur source. Measurements of SO₃ and SO (a¹Δ → X³Σ^-) emission at 1.7 μm may be the key to distinguish between the two models.     

Galactic longitude dependent galactic model parameters

We present the Galactic model parameters for thin disc estimated by Sloan Digital Sky Survey (SDSS) data of 14 940 stars with apparent magnitudes 16 < g₀ ⩽ 21 in six intermediate latitude fields in the first Galactic quadrant. Star/galaxy separation was performed by using the SDSS photometric pipeline and the isodensity contours in the (g − r)₀ − (r − i)₀ two colour diagram. The separation of thin disc stars is carried out by the bimodal distribution of stars in the (g − r)₀ histogram, and the absolute magnitudes were evaluated by a procedure presented in the literature (Bilir, S., Karaali, S., Tunçel, S. 2005. AN 326, 321). Exponential density law fits better to the derived density functions for the absolute magnitude intervals 8 < M(g) ⩽ 9 and 11 < M(g) ⩽ 12, whereas sech/sech² laws are more appropriate for absolute magnitude intervals 9 < M(g) ⩽ 10 and 10 < M(g) ⩽ 11. We showed that the scaleheight and scalelength are Galactic longitude dependent. The average values and ranges of the scaleheight and the scalelength are 〈H〉 = 220 pc (196 ⩽ H ⩽ 234 pc) and 〈H〉 = 1900 pc (1561 ⩽ h ⩽ 2280 pc) respectively. This result would be useful to explain different numerical values claimed for those parameters obtained by different authors for the fields in different directions of the Galaxy.     

Titan’s aerosol and stratospheric ice opacities between 18 and 500 μm: Vertical and spectral characteristics from Cassini CIRS

Vertical distributions and spectral characteristics of Titan’s photochemical aerosol and stratospheric ices are determined between 20 and 560 cm^?1 (500–18 μm) from the Cassini Composite Infrared Spectrometer (CIRS). Results are obtained for latitudes of 15°N, 15°S, and 58°S, where accurate temperature profiles can be independently determined.In addition, estimates of aerosol and ice abundances at 62°N relative to those at 15°S are derived. Aerosol abundances are comparable at the two latitudes, but stratospheric ices are ～3 times more abundant at 62°N than at 15°S. Generally, nitrile ice clouds (probably HCN and HC₃N), as inferred from a composite emission feature at ～160 cm^?1, appear to be located over a narrow altitude range in the stratosphere centered at ～90 km. Although most abundant at high northern latitudes, these nitrile ice clouds extend down through low latitudes and into mid southern latitudes, at least as far as 58°S.There is some evidence of a second ice cloud layer at ～60 km altitude at 58°S associated with an emission feature at ～80 cm^?1. We speculate that the identify of this cloud may be due to C₂H₆ ice, which in the vapor phase is the most abundant hydrocarbon (next to CH₄) in the stratosphere of Titan.Unlike the highly restricted range of altitudes (50–100 km) associated with organic condensate clouds, Titan’s photochemical aerosol appears to be well-mixed from the surface to the top of the stratosphere near an altitude of 300 km, and the spectral shape does not appear to change between 15°N and 58°S latitude. The ratio of aerosol-to-gas scale heights range from 1.3–2.4 at about 160 km to 1.1–1.4 at 300 km, although there is considerable variability with latitude. The aerosol exhibits a very broad emission feature peaking at ～140 cm^?1. Due to its extreme breadth and low wavenumber, we speculate that this feature may be caused by low-energy vibrations of two-dimensional lattice structures of large molecules. Examples of such molecules include polycyclic aromatic hydrocarbons (PAHs) and nitrogenated aromatics.Finally, volume extinction coefficients Nχ_E derived from 15°S CIRS data at a wavelength of λ = 62.5 μm are compared with those derived from the 10°S Huygens Descent Imager/Spectral Radiometer (DISR) data at 1.583 μm. This comparison yields volume extinction coefficient ratios Nχ_E(1.583 μm)/Nχ_E(62.5 μm) of roughly 70 and 20, respectively, for Titan’s aerosol and stratospheric ices. The inferred particle cross-section ratios χ_E(1.583 μm)/χ_E(62.5 μm) appear to be consistent with sub-micron size aerosol particles, and effective radii of only a few microns for stratospheric ice cloud particles.     

Rotation periods of binary asteroids with large separations – Confronting the Escaping Ejecta Binaries model with observations

Durda et al. (Durda, D.D., Bottke, W.F., Enke, B.L., Merline, W.J., Asphaug, E., Richardson, D.C., Leinhardt, Z.M. [2004]. Icarus 170, 243–257), using numerical models, suggested that binary asteroids with large separation, called Escaping Ejecta Binaries (EEBs), can be created by fragments ejected from a disruptive impact event. It is thought that six binary asteroids recently discovered might be EEBs because of the high separation between their components (～100 > a/R_p > ～20).However, the rotation periods of four out of the six objects measured by our group and others and presented here show that these suspected EEBs have fast rotation rates of 2.5–4 h. Because of the small size of the components of these binary asteroids, linked with this fast spinning, we conclude that the rotational-fission mechanism, which is a result of the thermal YORP effect, is the most likely formation scenario. Moreover, scaling the YORP effect for these objects shows that its timescale is shorter than the estimated ages of the three relevant Hirayama families hosting these binary asteroids. Therefore, only the largest (D ～ 19 km) suspected asteroid, (317) Roxane, could be, in fact, the only known EEB.In addition, our results confirm the triple nature of (3749) Balam by measuring mutual events on its lightcurve that match the orbital period of a nearby satellite in addition to its distant companion. Measurements of (1509) Esclangona at different apparitions show a unique shape of the lightcurve that might be explained by color variations.     

New wind measurements in Venus’ lower mesosphere from visible spectroscopy

The dynamics of Venus’ mesosphere (70–110 km) is characterized by the superposition of two different wind regimes: (1) Venus’ retrograde superrotation; (2) a sub-solar to anti-solar (SS–AS) flow pattern, driven by solar EUV heating on the sunlit hemisphere. Here, we report on new ground-based velocity measurements in the lower part of the mesosphere. We took advantage of two essentially symmetric Venus elongations in 2001 and 2002 to perform high-resolution Doppler spectroscopy (R=120,000) in ¹²C¹⁶O₂ visible lines of the 5ν₃ band and in a few solar Fraunhofer lines near 8700 Å. These measurements, mapped over several points on Venus’ illuminated hemisphere, probe the region of cloud tops. More precisely, the solar Fraunhofer lines sample levels a few kilometers below the UV features (i.e. near ∼67 km), while the CO₂ lines probe an altitude higher by about 7 km. The wind field over Venus’ disk is retrieved with an rms uncertainty of 15–25 m s⁻¹ on individual measurements. Kinematical fit to a one- or two-component circulation model indicates the dominance of the zonal retrograde flow with a mean equatorial velocity of ∼75 m s⁻¹, exhibiting very strong day-to-day variations (±65 m s⁻¹). Results are very consistent for the two kinds of lines, suggesting a negligible vertical wind shear over 67–74 km. The SS–AS flow is not detected in single-day observations, but combining the results from all data suggests that this component may invade the lower mesosphere with a ∼40 m s⁻¹ velocity.     

Identification of Zirconium Oxide molecular lines in the sunspot umbral spectra

The ZrO molecule has been detected in sunspot umbrae through the identification of following laboratory molecular transitions: ¹Σ⁺ ? X¹Σ⁺ (0, 0), A³Φ₂ ? X²Δ₁ (0, 0), A³Φ₃ ? X²Δ₂ (0, 0), A³Φ₄ ? X²Δ₃ (0, 0), B³Π₂ ? X³Δ₃ (0, 0), B³Π₁ ? X³Δ₂ (0, 0) and B³Π₀ ? X³Δ₁ (0, 0) in red – infrared region using high resolution, visible range Fourier Transform Spectrum of sunspot umbra observed at the National Solar Observatory in Kitt Peak (NSO/KP). Much new identification has been made in the searched spectral wavenumber region from 16650 cm^?1 to 18007 cm^?1 of sunspot spectrum. Equivalent widths of well resolved lines, versus rotational quantum number J have been used to determine the effective rotational temperature for seven bands of the ZrO molecule. This result agrees well with the temperatures derived for other molecules’ presence in sunspot umbrae. It is evident that ZrO molecular lines are formed in higher layers of the atmosphere of relatively “cold” sunspots.     

A photochemical model for the Venus atmosphere at 47–112 km

The model is intended to respond to the recent findings in the Venus atmosphere from the Venus Express and ground-based submillimeter and infrared observations. It extends down to 47 km for comparison with the kinetic model for the lower atmosphere (Krasnopolsky, V.A. [2007]. Icarus 191, 25–37) and to use its results as the boundary conditions. The model numerical accuracy is significantly improved by reduction of the altitude step from 2 km in the previous models to 0.5 km. Effects of the NUV absorber are approximated using the detailed photometric observations at 365 nm from Venera 14. The H₂O profile is not fixed but calculated in the model. The model involves odd nitrogen and OCS chemistries based on the detected NO and OCS abundances. The number of the reactions is significantly reduced by removing of unimportant processes. Column rates for all reactions are given, and balances of production and loss may be analyzed in detail for each species.The calculated vertical profiles of CO, H₂O, HCl, SO₂, SO, OCS and of the O₂ dayglow at 1.27 μm generally agree with the existing observational data; some differences are briefly discussed. The OH dayglow is ～30 kR, brighter than the OH nightglow by a factor of 4. The H + O₃ process dominates in the nightglow excitation and O + HO₂ in the dayglow, because of the reduction of ozone by photolysis. A key feature of Venus’ photochemistry is the formation of sulfuric acid in a narrow layer near the cloud tops that greatly reduces abundances of SO₂ and H₂O above the clouds. Delivery of SO₂ and H₂O through this bottleneck determines the chemistry and its variations above the clouds. Small variations of eddy diffusion near 60 km result in variations of SO₂, SO, and OCS at and above 70 km within a factor of ～30. Variations of the SO₂/H₂O ratio at the lower boundary have similar but weaker effect: the variations within a factor of ～4 are induced by changes of SO₂/H₂O by ±5%. Therefore the observed variations of the mesospheric composition originate from minor variations of the atmospheric dynamics near the cloud layer and do not require volcanism. NO cycles are responsible for production of a quarter of O₂, SO₂, and Cl₂ in the atmosphere. A net effect of photochemistry in the middle atmosphere is the consumption of CO₂, SO₂, and HCl from and return of CO, H₂SO₄, and SO₂Cl₂ to the lower atmosphere. These processes may be balanced by thermochemistry in the lower atmosphere even without outgassing from the interior, though the latter is not ruled out by our models. Some differences between the model and observations and the previous models are briefly discussed.     

Vertical profiling of SO2 and SO above Venus’ clouds by SPICAV/SOIR solar occultations

New measurements of sulfur dioxide (SO₂) and monoxide (SO) in the atmosphere of Venus by SPICAV/SOIR instrument onboard Venus Express orbiter provide ample statistics to study the behavior of these gases above Venus’ clouds. The instrument (a set of three spectrometers) is capable to sound atmospheric structure above the clouds in several observation modes (nadir, solar and stellar occultations) either in the UV or in the near IR spectral ranges. We present the results from solar occultations in the absorption ranges of SO₂ (190–230 nm, and at 4 μm) and SO (190–230 nm). The dioxide was detected by the SOIR spectrometer at the altitudes of 65–80 km in the IR and by the SPICAV spectrometer at 85–105 km in the UV. The monoxide’s absorption was measured only by SPICAV at 85–105 km. We analyzed 39 sessions of solar occultation, where boresights of both spectrometers are oriented identically, to provide complete vertical profiling of SO₂ of the Venus’ mesosphere (65–105 km). Here we report the first firm detection and measurements of two SO₂ layers. In the lower layer SO₂ mixing ratio is within 0.02–0.5 ppmv. The upper layer, also conceivable from microwave measurements by Sandor et al. (Sandor, B.J., Todd Clancy, R., Moriarty-Schieven, G., Mills, F.P. [2010]. Icarus 208, 49–60) is characterized by SO₂ increasing with the altitude from 0.05 to 2 ppmv, and the [SO₂]/[SO] ratio varying from 1 to 5. The presence of the high-altitude SO_x species could be explained by H₂SO₄ photodissociation under somewhat warmer temperature conditions in Venus mesosphere. At 90–100 km the content of the sulfur dioxide correlates with temperature increasing from 0.1 ppmv at 165–170 K to 0.5–1 ppmv at 190–192 K. It supports the hypothesis of SO₂ production by the evaporation of H₂SO₄ from droplets and its subsequent photolysis at around 100 km.     

O2 1.27 μm emission maps as derived from OMEGA/MEx data

We present a method to derive the 2D maps of the O₂ (a¹Δ_g) airglow emission at 1.27 μm from the OMEGA/MEx nadir observations. The OMEGA imaging capabilities allow monitoring the 2D distribution, daily and seasonal variation of the O₂ emission intensities with a detection limit of 4 MR. The highest values, of the order of ∼31 MR, are found on the south pole for 11 h < LT < 13 h, during the early spring (186° < Ls < 192°) of martian year (MY) 27, according to the Mars Year numbering scheme of Clancy et al. [Clancy, R.T., Wolff, M.J., Christensen, P.R., 2003. Mars aerosol studies with the MGS TES emission phase function observations: Optical depths, particle sizes, and ice cloud types versus latitude and solar longitude. J. Geophys. Res. 108. doi: 10.1029/2003JE002058]. In the polar regions the day-by-day variability, associated with polar vortex turbulences, is obtained of the order of 30-50% as predicted by the model [Lefévre, F., Lebonnois, S., Montmessin, F., Forget, F., 2004. Three-dimensional modeling of ozone on Mars. J. Geophys. Res. 109, E07004. doi: 10.1029/2004JE002268] and found by SPICAM [Perrier, S., Bertaux, J.-L., Lebonnois, S., Korablev, O., Fedorova, A., 2006. Global distribution of total ozone on Mars from SPICAM/MEX UV measurements. J. Geophys. Res. 111, E09S06. doi: 10.1029/2006JE002681]. In the considered set of data a maximum of the O₂ emission is observed between 11 h and 15 h LT in the latitude range 70-85° during early spring on both hemispheres, while for the southern autumn-winter season a maximum is found between 50° and 60° in the southern hemisphere for MY28. Increase of intensity of the O₂ emission observed from Ls 130° to 160° at southern high latitudes may be explained by increase of solar illumination conditions in the maps acquired during the considered period.Atmospheric waves crossing the terminator on the southern polar regions are observed for the first time during the MY28 early spring. The spatial scale of the waves ranges from 100 to 130 km, and the intensity fluctuations are of the order of 4MR.This study confirms the high potentiality of O₂ (a¹Δ_g) day glow as a passive tracer of the martian atmosphere dynamics at high latitudes.     

Thermal structure and CO distribution for the Venus mesosphere/lower thermosphere: 2001–2009 inferior conjunction sub-millimeter CO absorption line observations

Sub-millimeter ¹²CO (346 GHz) and ¹³CO (330 GHz) line absorptions, formed in the mesosphere and lower thermosphere of Venus (70–120 km), have been mapped across the nightside Venus disk during 2001–2009 inferior conjunctions, employing the James Clerk Maxwell Telescope (JCMT). Radiative transfer analysis of these thermal line absorptions supports temperature and CO mixing profile retrievals, as well as Doppler wind fields (described in the companion paper, Clancy et al., 2012). Temporal sampling over the hourly, daily, weekly and interannual timescales was obtained over 2001–2009. On timescales inferred as several weeks, we observe changes between very distinctive CO and temperature nightside distributions. Retrieved nightside CO, temperature distributions for January 2006 and August 2007 observations display strong local time, latitudinal gradients consistent with early morning (2–3 am), low-to-mid latitude (0–40NS) peaks of 100–200% in CO and 20–30 K in temperature. The temperature increases are most pronounced above 100 km altitudes, whereas CO variations extend from 105 km (top altitude of retrieval) down to below 80 km in the mesosphere. In contrast, the 2004 and 2009 periods of observation display modest temperature (5–10 K) and CO (30–60%) increases, that are centered on antisolar (midnight) local times and equatorial latitudes. Doppler wind derived global (zonal and should be SSAS) circulations from the same data do not exhibit variations correlated with these CO, temperature short-term variations. However, large-scale residual wind fields not fit by the zonal, SSAS circulations are observed in concert with the strong temperature, CO gradients observed in 2006 and 2007 (Clancy et al., 2010). These short term variations in nightside CO, temperature distributions may also be related to observed nightside variations in O₂ airglow (Hueso, H., Sánchez-Lavega, A., Piccioni, G., Drossart, P., Gérard, J.C., Khatuntsev, I., Zasova, L., Migliorini, A. [2008]. J. Geophys. Res. 113, E00B02. doi:10.1029/2008JE003081) and upper mesospheric SO and SO₂ layers (Sandor, B.J., Clancy, R.T., Moriarty-Schieven, G.H., Mills, F.P. [2010]. Icarus 208, 49–60).The retrieved temperature profiles also exhibit 20 K long-term (2001–2009) variations in nightside (whole disk) average mesospheric (80–95 km) temperatures, similar to 1982–1991 variations identified in previous millimeter CO line observations (Clancy et al., 1991). Global average diurnal variations in lower thermospheric temperatures and mesospheric CO abundances decreased by a factor-of-two between 2000–2002 versus 2007–2009 periods of combined dayside and nightside observations. The infrequency and still limited temporal extent of the observations make it difficult to assign specific timescales to such longer term variations, which may be associated with longer term variations observed for cloud top SO₂ (Esposito, L.W., Bertaux, J.-L., Krasnopolsky, V., Moroz, V.I., Zasova, L.V. [1997]. Chemistry of lower atmosphere and clouds. In: Bougher, S.W., Hunten, D.M., Phillips, R.J. (Eds.), VENUS II, 1362pp) and mesospheric water vapor (Sandor, B.J., Clancy, R.T. [2005]. Icarus 177, 129–143) abundances.     

Characterizing Io’s Pele,Tvashtar and Pillan plumes: Lessons learned from Hubble

Hubble Space Telescope/Wide Field and Planetary Camera 2 (HST/WFPC2) images of Io obtained between 1995 and 2007 between 0.24 and 0.42 μm led to the detection of the Pele plume in reflected sunlight in 1995 and 1999; imaging of the Pele plume via absorption of jovian light in 1996 and 1999; detection of the Prometheus-type Pillan plume in reflected sunlight in 1997; and detection of the 2007 Pele-type Tvashtar plume eruption in reflected sunlight and via absorption of jovian light. Based on a detailed analysis of these observations we characterize and compare the gas and dust properties of each of the detected plumes. In each case, the brightness of the plumes in reflected sunlight is less at 0.26 μm than at 0.33 μm. Mie scattering analysis of the wavelength dependence of each plume’s reflectance signature suggests that range of particle sizes within the plumes is quite narrow. Assuming a normal distribution of particle sizes, the range of mean particle sizes is ～0.035–0.12 μm for the 1997 Pillan eruption, ～0.05–0.08 μm for the 1999 Pele and 2007 Tvasthar plumes, and ～0.05–0.11 μm for the 1995 Pele plume, and in each case the standard deviation in the particle size distribution is <15%. The Mie analysis also suggests that the 2007 Tvashtar eruption released ～10⁹ g of sulfur dust, the 1999 Pele eruption released ～10⁹ g of SO₂ dust, the 1997 Pillan eruption released ～10¹⁰ g of SO₂ dust, and the 1995 Pele plume may have released ～10¹⁰ g of SO₂ dust. Analysis of the plume absorption signatures recorded in the F255W filter bandpass (0.24–0.28 μm) indicates that the opacity of the 2007 Tvashtar plume was 2× that of the 1996 and 1999 Pele plume eruptions. While the sulfur dust density estimated for the Tvashtar from the reflected sunlight data could have produced 61% of the observed plume opacity, <10% of the 1999 Pele F255W plume opacity could have resulted from the SO₂ dust detected in the eruption. Accounting for the remaining F255W opacity level of the Pele and Tvasthar plumes based on SO₂ and S₂ gas absorption, the SO₂ and S₂ gas density inferred for each plume is almost equivalent corresponding to ～2–6 × 10¹⁶ cm^?2 and 3–5 × 10¹⁵ cm^?2, respectively, producing SO₂ and S₂ gas resurfacing rates ～0.04–0.2 cm yr^?1 and 0.007–0.01 cm yr^?1; and SO₂ and S₂ gas masses ～1–4 × 10¹⁰ g and ～2–3 × 10⁹ g; for a total dust to gas ratio in the plumes ～10^?1–10^?2. The 2007 Tvashtar plume was detected by HST at ～380 ± 40 km in both reflected sunlight and absorbed jovian light; in 1999, the detected Pele plume altitude was 500 km in absorbed jovian light, but in reflected sunlight the detected height was ～2× lower. Thus, for the 1999 Pele plume, similar to the 1979 Voyager Pele plume observations, the most efficient dust reflections occurred in the region closest to the plume vent. The 0.33–0.42 μm brightness of the 1997 Pillan plume was 10–20× greater than the Pele or Tvashtar plumes, exceeding by a factor of 3 the average brightness levels observed within 200 km of 1979 Loki eruption vent. But, the 0.26 μm brightness of the 1997 Pillan plume in reflected sunlight was significantly lower than would be predicted by the dust scattering model. Presuming that the 0.26 μm brightness of the 1997 Pillan plume was attenuated by the eruption plume’s gas component, then an SO₂ gas density ～3–6 × 10¹⁸ cm^?2 is inferred from the data (for S₂/SO₂ ratios ?4%), comparable to the 0.3–2 × 10¹⁸ cm^?2 SO₂ density detected at Loki in 1979 (Pearl, J.C. et al. [1979]. Nature 280, 755; Lellouch et al., 1992), and producing an SO₂ gas mass ～3–8 × 10¹¹ g and an SO₂ resurfacing rate ～8–23 cm yr^?1. These results confirm the connection between high (?10¹⁷ cm^?2) SO₂ gas content and plumes that scatter strongly at nearly blue wavelengths, and it validates the occurrence of high density SO₂ gas eruptions on Io. Noting that the SO₂ gas content inferred from a spectrum of the 2003 Pillan plume was significantly lower ～2 × 10¹⁶ cm^?2 (Jessup, K.L., Spencer, J., Yelle, R. [2007]. Icarus 192, 24–40); and that the Pillan caldera was flooded with fresh SO₂ frost/slush just prior to the 1997 Pillan plume eruption (Geissler, P., McEwen, A., Phillips, C., Keszthelyi, L., Spencer, J. [2004a]. Icarus 169, 29–64; Phillips, C.B. [2000]. Voyager and Galileo SSI Views of Volcanic Resurfacing on Io and the Search for Geologic Activity at Europa. Ph.D. Thesis, Univ. of Ariz., Tucson); we propose that the density of SO₂ gas released by this volcano is directly linked to the local SO₂ frost abundance at the time of eruption.     

SPICAV on Venus Express: Three spectrometers to study the global structure and composition of the Venus atmosphere

Spectroscopy for the investigation of the characteristics of the atmosphere of Venus (SPICAV) is a suite of three spectrometers in the UV and IR range with a total mass of 13.9 kg flying on the Venus Express (VEX) orbiter, dedicated to the study of the atmosphere of Venus from ground level to the outermost hydrogen corona at more than 40,000 km. It is derived from the SPICAM instrument already flying on board Mars Express (MEX) with great success, with the addition of a new IR high-resolution spectrometer, solar occultation IR (SOIR), working in the solar occultation mode. The instrument consists of three spectrometers and a simple data processing unit providing the interface of these channels with the spacecraft.A UV spectrometer (118–320 nm, resolution 1.5 nm) is identical to the MEX version. It is dedicated to nadir viewing, limb viewing and vertical profiling by stellar and solar occultation. In nadir orientation, SPICAV UV will analyse the albedo spectrum (solar light scattered back from the clouds) to retrieve SO₂, and the distribution of the UV-blue absorber (of still unknown origin) on the dayside with implications for cloud structure and atmospheric dynamics. On the nightside, γ and δ bands of NO will be studied, as well as emissions produced by electron precipitations. In the stellar occultation mode the UV sensor will measure the vertical profiles of CO₂, temperature, SO₂, SO, clouds and aerosols. The density/temperature profiles obtained with SPICAV will constrain and aid in the development of dynamical atmospheric models, from cloud top (∼60 km) to 160 km in the atmosphere. This is essential for future missions that would rely on aerocapture and aerobraking. UV observations of the upper atmosphere will allow studies of the ionosphere through the emissions of CO, CO⁺, and CO₂⁺, and its direct interaction with the solar wind. It will study the H corona, with its two different scale heights, and it will allow a better understanding of escape mechanisms and estimates of their magnitude, crucial for insight into the long-term evolution of the atmosphere.The SPICAV VIS-IR sensor (0.7–1.7 μm, resolution 0.5–1.2 nm) employs a pioneering technology: an acousto-optical tunable filter (AOTF). On the nightside, it will study the thermal emission peeping through the clouds, complementing the observations of both VIRTIS and Planetary Fourier Spectrometer (PFS) on VEX. In solar occultation mode this channel will study the vertical structure of H₂O, CO₂, and aerosols.The SOIR spectrometer is a new solar occultation IR spectrometer in the range λ=2.2–4.3 μm, with a spectral resolution λ/Δλ>15,000, the highest on board VEX. This new concept includes a combination of an echelle grating and an AOTF crystal to sort out one order at a time. The main objective is to measure HDO and H₂O in solar occultation, in order to characterize the escape of D atoms from the upper atmosphere and give more insight about the evolution of water on Venus. It will also study isotopes of CO₂ and minor species, and provides a sensitive search for new species in the upper atmosphere of Venus. It will attempt to measure also the nightside emission, which would allow a sensitive measurement of HDO in the lower atmosphere, to be compared to the ratio in the upper atmosphere, and possibly discover new minor atmospheric constituents.     

A photometric study of the detached binary WY Hydrae with twin components

We present a multicolor photometry for the eclipsing binary WY Hydrae, observed on four nights of 2008 December. From our new observations and Carr’s data, the photometric solutions were deduced by using the updated W–D program. The results show that WY Hya is a detached binary with a mass ratio of q = 0.970(±0.005).By analyzing the O–C curve, it is found that there exists either a continuous period increase or a cyclic variation. From Eq. (2), the orbital period of WY Hya secularly increases at a rate of dP/dt = +3.56(±0.37) × 10^?7 days/yr, which may be interpreted by some mass transfer for the near-contact configuration or tidal dissipation. From Eq. (3), the period and semi-amplitude of the periodic oscillation are P₃ = 95.4(±4.2) yr and A = 0^d.0087(±0^d.0003), respectively. This may be likely attributed by light-time effect via the presence of the assumed third body. Assumed in the coplanar orbit with the binary, the mass of the third body should be M₃ = 0.18 M_⊙. If the unseen additional companion exists, it will extract angular momentum from the binary system. Finally, WY Hya with high fill-out factors (i.e., f_1,2 > 80%), may evolve into a semi-detached configuration.     

Ensemble properties of comets in the Sloan Digital Sky Survey

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u, g, r, i, z) photometry is used to determine the comets’ colors, sizes, surface brightness profiles, and rates of dust production in terms of the Afρ formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(<H) ∝ 10^{(0.49±0.05)H} for H < 18, with evidence of a much shallower fit N(<H) ∝ 10^{(0.19±0.03)H} for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 ± 0.05 in g ? r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.     

A Study of Potential Earth-like Planets in the GJ 876 Planetary System

We perform numerical simulations to investigate potential Earth-like planets in the GJ 876 planetary system. We show that the secular resonances ν₁ and ν₂ (resulting respectively from the inner and outer giant planets) can excite the eccentricities of Earth-like planets with orbits 0.21 AU ≤ a < 0.50 AU and cause them to be ejected out of the system in a short time. However, in the dynamical sense, Earth-like planets potentially exist in the region 0.50 AU ≤ a ≤ 1.00 AU, in stable low-eccentricity orbits which may last up to 10⁵ yr.     

Optical constants from 370 nm to 900 nm of Titan tholins produced in a low pressure RF plasma discharge

Determining the optical constants of Titan aerosol analogues, or tholins, has been a major concern for the last three decades because they are essential to constrain the numerical models used to analyze Titan’s observational data (albedo, radiative transfer, haze vertical profile, surface contribution, etc.). Here we present the optical constant characterization of tholins produced with an RF plasma discharge in a (95%N₂–5%CH₄) gas mixture simulating Titan’s main atmospheric composition, and deposited as a thin film on an Al–SiO₂ substrate. The real and imaginary parts, n and k, of the tholin complex refractive index have been determined from 370 nm to 900 nm wavelength using spectroscopic ellipsometry. The values of n decrease from n = 1.64 (at 370 nm) to n = 1.57 (at 900 nm) as well as the values of k which feature two behaviors: an exponential decay from 370 nm to 500 nm, with k = 12.4 × e^?0.018λ (where λ is expressed in nm), followed by a plateau, with k = (1.8 ± 0.2) × 10^?3. The trends observed for the PAMPRE tholins optical constants are compared to those determined for other Titan tholins, as well as to the optical constants of Titan’s aerosols retrieved from observational data.