The growth of wind-waves in Titan’s hydrocarbon seas

Ocean wave growth on Titan is considered. The classic Sverdrup–Munk theory for terrestrial wave growth is applied to Titan, and is compared with a simple energy balance model that exposes the effect of Titan’s environmental parameters (air density, gravity, and fluid density). These approaches are compared with the only previously-published (semi-empirical) model (Ghafoor, N.A.-L., Zarnecki, J.C., Challenor, P., Srokosz, M.A. [2000] J. Geophys. Res. 105, 12,077–12,091, hereafter G2k), and allow the impact of various parameters such as atmospheric density to be transparently explored.Our model, like G2k, suggests fully-developed significant wave heights on Titan H_s = 0.2 U², where U is the windspeed (SI units): in dimensionless terms this is rather close to H_s = 0.2 U²/g, a rule of thumb previously noted for terrestrial waves (we find various datasets where the prefactor varies by ～2). It is noted that liquid and air densities affect the growth rate of waves, but not their fully-developed height: for 1 m/s winds wave amplitude reaches 0.15 m (75% of fully-developed) with a fetch of only 1 km, rather faster than predicted by G2k. Liquid viscosity has no major effect on gravity wave growth, but does influence the threshold windspeed at which gravity–capillary waves form in the first place.The model is used to develop predicted ranges for wave height to guide the design of the Titan Mare Explorer (TiME), a proposed Discovery-class mission to float a capsule on Ligeia Mare in 2023. For the expected maximum 1 m/s winds, a significant wave height of 0.2 m and wavelength of ～4 m can be expected. Assuming that wave heights follow Rayleigh statistics as they do on Earth, then given the wave period of ～4 s, individual waves of ～0.6 m might be encountered over a 3 month period.For predicted Titan winds at Kraken Mare, significant wave heights may reach ～0.6 m in the peak of summer but do not exceed the tidal amplitude at its northern end, consistent with the area around Mayda Insula being a tidal flat, while elsewhere on Kraken and Ligeia and at Ontario Lacus, shorelines may be wave- or tidally-dominated, depending on the specific location.     

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.     

The mass,orbit, and tidal evolution of the Quaoar–Weywot system

Here we present new adaptive optics observations of the Quaoar–Weywot system. With these new observations we determine an improved system orbit. Due to a 0.39 day alias that exists in available observations, four possible orbital solutions are available with periods of ～11.6, ～12.0, ～12.4, and ～12.8 days. From the possible orbital solutions, system masses of 1.3–1.5 ± 0.1 × 10²¹ kg are found. These observations provide an updated density for Quaoar of 2.7–5.0 g cm^?3. In all cases, Weywot’s orbit is eccentric, with possible values ～0.13–0.16. We present a reanalysis of the tidal orbital evolution of the Quaoar–Weywot system. We have found that Weywot has probably evolved to a state of synchronous rotation, and has likely preserved its initial inclination over the age of the Solar System. We find that for plausible values of the effective tidal dissipation factor tides produce a very slow evolution of Weywot’s eccentricity and semi-major axis. Accordingly, it appears that Weywot’s eccentricity likely did not tidally evolve to its current value from an initially circular orbit. Rather, it seems that some other mechanism has raised its eccentricity post-formation, or Weywot formed with a non-negligible eccentricity.     

Descent motions of the Huygens probe as measured by the Surface Science Package (SSP): Turbulent evidence for a cloud layer

The Huygens probe underwent vigorous short-period motions during its parachute descent through the atmosphere of Saturn's moon Titan in January 2005, at least some of which were excited by the Titan environment. Several sensors in the Huygens Surface Science Package (SSP) detect these motions, indicating the transition to the smaller stabilizer parachute, the changing probe spin rate, aerodynamic buffeting, and pendulum motions. Notably, in an altitude range of about 20–30 km where methane drops will freeze, the frequency content and statistical kurtosis of the tilt data indicate excitation by turbulent air motions like those observed in freezing clouds on Earth, supporting the suggestion of Tokano et al. [Tokano, T., McKay, C.P., Neubauer, F.M., Atreya, S.K., Ferri, F., Fulchignoni, M., Niemann, H.B. (2006a). Methane drizzle on Titan. Nature 442, 432–435] that the probe passed through such a cloud layer. Motions are weak below 20 km, suggesting a quiescent lower atmosphere with turbulent fluctuations of nominally <0.15 m/s (to within a factor of ∼2) but more violent motions in the upper troposphere may have been excited by turbulent winds with amplitudes of 1–2 m/s. Descent in part of the stratosphere (150–120 km) was smooth despite strong ambient wind (∼100 m/s), but known anomalies in the probe spin prevent investigation of turbulence in the known wind-shear layer from 60 to 100 km.     

A multi-wavelength study of the 2009 impact on Jupiter: Comparison of high resolution images from Gemini,Keck and HST

Within several days of A. Wesley’s announcement that Jupiter was hit by an object on UT 19 July 2009, we observed the impact site with (1) the Hubble Space Telescope (HST) at UV through visible (225–924 nm) wavelengths, (2) the 10-m W.M. Keck II telescope in the near-infrared (1–5 μm), and (3) the 8-m Gemini-North telescope in the mid-infrared (7.7–18 μm). All observations reported here were obtained between 22 and 25 July 2009. Observations at visible and near-infrared wavelengths show that large (～0.75-μm radius) dark (imaginary index of refraction m_i ～ 0.01–0.1) particulates were deposited at atmospheric pressures between 10 and 200–300 mbar; analysis of HST-UV data reveals that in addition smaller-sized (～0.1 μm radius) material must have been deposited at the highest altitudes (～10 mbar). Differences in morphology between the UV and visible/near-IR images suggest three-dimensional variations in particle size and density across the impact site, which probably were induced during the explosion and associated events. At mid-infrared wavelengths the brightness temperature increased due to both an enhancement in the stratospheric NH₃ gas abundance and the physical temperature of the atmosphere. This high brightness temperature coincides with the center part of the impact site as seen with HST. This observation, combined with (published) numerical simulations of the Shoemaker-Levy 9 impacts on Jupiter and the Tunguska airburst on Earth, suggests that the downward jet from the terminal explosion probably penetrated down to the ～700-mbar level.     

Comparisons of Cassini flybys of the Titan magnetospheric interaction with an MHD model: Evidence for organized behavior at high altitudes

Recent papers suggest the significant variability of conditions in Saturn’s magnetosphere at the orbit of Titan. Because of this variability, it was expected that models would generally have a difficult time regularly comparing to data from the Titan flybys. However, we find that in contrast to this expectation, it appears that there is underlying organization of the interaction features roughly above ～1800 km (1.7 Rt) altitude by the average external field due to Saturn’s dipole moment. In this study, we analyze Cassini’s plasma and magnetic field data collected at 9 Titan encounters during which the external field is close to the ideal southward direction and compare these observations to the results from a 2-fluid (1 ion, 1 electron) 7-species MHD model simulations obtained under noon SLT conditions. Our comparative analysis shows that under noon SLT conditions the Titan plasma interaction can be viewed in two layers: an outer layer between 6400 and 1800 km where interaction features observed in the magnetic field are in basic agreement with a purely southward external field interaction and an inner layer below 1800 km where the magnetic field measurements show strong variations and deviate from the model predictions. Thus the basic features inferred from the Voyager 1 flyby seem to be generally present above ～1800 km in spite of the ongoing external variations from SLT excursions, time variability and magnetospheric current systems as long as a significant southward external field component is present. At around ～1800 km kinetic effects (such as mass loading and heavy ion pickup) and below 1800 km ionospheric effects (such as drag of ionospheric plasma due to coupling with neutral winds and/or magnetic memory of Titan’s ionosphere) complicate what is observed.     

Direct wind measurements from November 2007 in Venus’ upper atmosphere using ground-based heterodyne spectroscopy of CO2 at 10 μm wavelength

Between November 23 and 28, 2007, the Cologne Tuneable Heterodyne Infrared Spectrometer THIS was installed at the McMath-Pierce Solar Telescope (Kitt Peak, Arizona, USA) to determine zonal wind velocities and to estimate the subsolar-to-antisolar flow. We investigate dynamics in the upper atmosphere of Venus by measuring the Doppler shift of fully-resolved non-LTE CO₂ emission lines at 959.3917 cm^?1 (10.423 μm), which probe a narrow altitude region in Venus’ atmosphere around 110 ± 10 km (～1 μbar). The results show no significant zonal wind velocity at the equator. An increase with latitude up to 43 ± 13 m/s at a latitude of 33°N was observed. This confirms the deduction of a minor influence of Venus superrotation at an altitude of 110 km from previous measurements in May 2007 (Sornig et al., 2008). The specific observing geometry enables estimating the maximum cross terminator velocity of the subsolar-to-antisolar flow at 72 ± 47 m/s.     

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.     

The layered structure of lunar maria: Identification of the HF-radar reflector in Mare Serenitatis using multiband optical images

Comparison of the Lunar Radar Sounder (LRS) data to the Multiband Imager (MI) data is performed to identify the subsurface reflectors in Mare Serenitatis. The LRS is FM-CW radar (4–6 MHz) and the 2 MHz bandwidth leads to the range resolution of 75 m in a vacuum, whereas the sampling interval in the flight direction is about 75 m when an altitude of the spacecraft with polar orbit is nominal (100 km). Horizontally continuous reflectors were clearly detected by LRS in limited areas that consist of about 9% of the whole maria. The typical depth of the reflectors is estimated to be a few hundred meters. Layered structures of mare basalts are also discernible on some crater walls in the MI data of the visible bands (VIS). The VIS range has nine wavelengths of 415, 750, 900, 950, and 1000 nm, and their spatial resolution is 20 m/pixel at a nominal altitude. The stratigraphies around Bessel and Bessel-H craters in Mare Serenitatis are examined in this paper. It was revealed that the subsurface reflectors lie on the boundaries between basalt units with different chemical compositions. In addition, model calculations using the simplified radar equation indicate that the subsurface reflectors are not compositional interfaces but layer boundaries with a high-porosity contrast. These results suggest that the detected reflectors in Mare Serenitatis are regolith accumulated during so long hiatus of mare volcanisms enough for chemical composition of magma to change, not instantaneously. Therefore combination of the LRS and MI data has a potential to reveal characteristics of a series of magmatism forming each lithostratigraphic unit in Mare Serenitatis and other maria.     

Water vapor near the cloud tops of Venus from Venus Express/VIRTIS dayside data

Observations of the dayside of Venus performed by the high spectral resolution channel (–H) of the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) on board the ESA Venus Express mission have been used to measure the altitude of the cloud tops and the water vapor abundance around this level with a spatial resolution ranging from 100 to 10 km. CO₂ and H₂O bands between 2.48 and 2.60 μm are analyzed to determine the cloud top altitude and water vapor abundance near this level. At low latitudes (±40°) mean water vapor abundance is equal to 3 ± 1 ppm and the corresponding cloud top altitude at 2.5 μm is equal to 69.5 ± 2 km. Poleward from middle latitudes the cloud top altitude gradually decreases down to 64 km, while the average H₂O abundance reaches its maximum of 5 ppm at 80° of latitude with a large scatter from 1 to 15 ppm. The calculated mass percentage of the sulfuric acid solution in cloud droplets of mode 2 (～1 μm) particles is in the range 75–83%, being in even more narrow interval of 80–83% in low latitudes. No systematic correlation of the dark UV markings with the cloud top altitude or water vapor has been observed.     

Hydrogen density in the dayside venusian exosphere derived from Lyman-α observations by SPICAV on Venus Express

A series of observations of the venusian hydrogen corona made by SPICAV on Venus Express are analyzed to estimate the amount of hydrogen in the exosphere of Venus. These observations were made between November 2006 and July 2007 at altitudes from 1000 km to 8000 km on the dayside. The Lyman-α brightness profiles derived are reproduced by the sum of a cold hydrogen population dominant below ～2000 km and a hot hydrogen population dominant above ～4000 km. The temperature (～300 K) and hydrogen density at 250 km (～10⁵ cm^?3) derived for the cold populations, near noon, are in good agreement with previous observations. Strong dawn–dusk exospheric asymmetry is observed from this set of observations, with a larger exobase density on the dawn side than on the dusk side, consistent with asymmetry previously observed in the venusian thermosphere, but with a lower dawn/dusk contrast. The hot hydrogen density derived is very sensitive to the sky background estimate, but is well constrained near 5000 km. The density of the hot population is reproduced by the exospheric model from Hodges (Hodges, R.R. [1999]. J. Geophys. Res. 104, 8463–8471) in which the hot population is produced by neutral–ions interactions in the thermosphere of Venus.     

Circulation of the Venus upper mesosphere/lower thermosphere: Doppler wind measurements from 2001–2009 inferior conjunction,sub-millimeter CO absorption line observations

Sub-millimeter ¹²CO (346 GHz) and ¹³CO (330 GHz) line absorptions, formed within the mesospheric to lower thermospheric altitude (70–120 km) region of the Venus atmosphere, have been mapped across the nightside disk of Venus 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 described in a companion paper (Clancy et al., 2012). Here, we consider the analysis of the sharp line absorption cores of these CO spectra in terms of accurate Doppler wind profile measurements at 95–115 km altitudes versus local time (～8 pm–4 am) and latitude (～60N–60S). These Doppler wind measurements support determinations of the nightside zonal and subsolar-to-antisolar (SSAS) circulation components over a variety of timescales. The average behavior fitted from 21 retrieved maps of ¹²CO Doppler winds (obtained over hourly, daily, weekly, and interannual intervals) indicates stronger average zonal (85 m/s retrograde) versus SSAS (65 m/s) circulation at the 1 μbar pressure (108–110 km altitude) level. However, the absolute and relative magnitudes of these circulation components exhibit extreme variability over daily to weekly timescales. Furthermore, the individual Doppler wind measurements within each nightside mapping observation generally show significant deviations (20–50 m/s, averaged over 5000 km horizontal scales) from the simple zonal/SSAS solution, with distinct local time and latitudinal characters that are also time variable. These large scale residual circulations contribute 30–70% of the observed nightside Doppler winds at any given time, and may be most responsible for global variations in nightside lower thermospheric trace composition and temperatures, as coincidentally retrieved CO abundance and temperature distributions do not correlate with solution retrograde zonal and SSAS winds (see companion paper, Clancy et al., 2012). Limited comparisons of these nightside submillimeter results with dayside infrared Doppler wind measurements suggest distinct dayside versus nightside circulations, in terms of zonal winds in particular. Combined ¹²CO and ¹³CO Doppler wind mapping observations obtained since 2004 indicate that the average zonal and SSAS wind components increase by 50–100% between altitudes of 100 and 115 km. If gravity waves originating from the cloud levels are responsible for the extension of zonal winds into the thermosphere (Alexander, M.J. [1992]. Geophys. Res. Lett. 19, 2207–2210), such waves deposit substantial momentum (i.e., break) in the lower nightside thermosphere.     

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.     

Ground-based near-infrared observations of water vapour in the Venus troposphere

We present a study of water vapour in the Venus troposphere obtained by modelling specific water vapour absorption bands within the 1.18 μm window. We compare the results with the normal technique of obtaining the abundance by matching the peak of the 1.18 μm window. Ground-based infrared imaging spectroscopy of the night side of Venus was obtained with the Anglo-Australian Telescope and IRIS2 instrument with a spectral resolving power of R ～ 2400. The spectra have been fitted with modelled spectra simulated using the radiative transfer model VSTAR. We find a best fit abundance of 31 ppmv (?6 +9 ppmv), which is in agreement with recent results by Bézard et al. (Bézard, B., Fedorova, A., Bertaux, J.-L., Rodin, A., Korablev, O. [2011]. Icarus, 216, 173–183) using VEX/SPICAV (R ～ 1700) and contrary to prior results by Bézard et al. (Bézard, B., de Bergh, C., Crisp, D., Maillard, J.P. [1990]. Nature, 345, 508–511) of 44 ppmv (±9 ppmv) using VEX/VIRTIS-M (R ～ 200) data analyses. Comparison studies are made between water vapour abundances determined from the peak of the 1.18 μm window and abundances determined from different water vapour absorption features within the near infrared window. We find that water vapour abundances determined over the peak of the 1. 18 μm window results in plots with less scatter than those of the individual water vapour features and that analyses conducted over some individual water vapour features are more sensitive to variation in water vapour than those over the peak of the 1. 18 μm window. No evidence for horizontal spatial variations across the night side of the disk are found within the limits of our data with the exception of a possible small decrease in water vapour from the equator to the north pole. We present spectral ratios that show water vapour absorption from within the lowest 4 km of the Venus atmosphere only, and discuss the possible existence of a decreasing water vapour concentration towards the surface.     

Sodium atoms in the lunar exotail: Observed velocity and spatial distributions

The lunar sodium tail extends long distances due to radiation pressure on sodium atoms in the lunar exosphere. Our earlier observations measured the average radial velocity of sodium atoms moving down the lunar tail beyond Earth (i.e., near the anti-lunar point) to be ～12.5 km/s. Here we use the Wisconsin H-alpha Mapper to obtain the first kinematically resolved maps of the intensity and velocity distribution of this emission over a 15° × 15 ° region on the sky near the anti-lunar point. We present both spatially and spectrally resolved observations obtained over four nights bracketing new Moon in October 2007. The spatial distribution of the sodium atoms is elongated along the ecliptic with the location of the peak intensity drifting 3° east along the ecliptic per night. Preliminary modeling results suggest the spatial and velocity distributions in the sodium exotail are sensitive to the near surface lunar sodium velocity distribution. Future observations of this sort along with detailed modeling offer new opportunities to describe the time history of lunar surface sputtering over several days.     

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.     

An investigation of the field-aligned currents associated with a large-scale ULF wave in the morning sector

Previous work by Scoffield, H.C., Yeoman, T.K., Wright, D.M., Milan, S.E., Wright, A.N., Strangeway, R.J. [2005. An investigation of the field aligned currents associated with a large scale ULF wave using data from CUTLASS and FAST. Ann. Geophys. 23, 487–498) investigated a large-scale ULF wave, occurring in the dusk sector (∼1900 MLT). The wave had a period of ∼800 s (corresponding to 1.2 mHz frequency), an azimuthal wave number of ∼7 and a full-width at half-maximum (FWHM) across the resonance of 350 km. IMAGE ground magnetometer and SuperDARN radar observations of the wave's spatial and temporal characteristics were used to parameterise a simple, two-dimensional field line resonance (FLR) model. The model-calculated field-aligned current (FAC) was compared with FACs derived from the FAST energetic particle spectra and magnetic field measurement. Here the authors use the same method to investigate the FAC structure of a second large-scale ULF wave, with a period of ∼450 s, occurring the dawn sector (∼0500 MLT) with an opposite sense background region 1–region 2 current system. This wave has a much larger longitudinal scale (m∼4.5) and a smaller latitude scale (FWHM=150 km). Unlike the dusk sector wave, which was dominated by upward FAC, FAST observations of the dawn sector wave show an interval of large-scale downward FAC of ∼1.5 μA m⁻². Downgoing magnetospheric electrons with energies of a few keV were observed, which are associated with upward FACs of ∼1 μA m⁻². For both wave studies, downward currents appear to be carried partially by upgoing electrons below the FAST energy detection threshold (5 eV), but also consist of a mixture of hotter downgoing magnetospheric electrons and upgoing ionospheric electrons of energies 30 eV–1 keV. Strong intervals of upward current show that small-scale structuring of scale ∼50 km has been imposed on the current carriers. In general, this study confirms the findings of Scoffield, H.C., Yeoman, T.K., Wright, D.M., Milan, S.E., Wright, A.N., Strangeway, R.J. [2005. An investigation of the FACs associated with a large-scale ULF wave using data from CUTLASS and FAST. Ann. Geophys. 23, 487–498).     

The abundance and vertical distribution of the unknown ultraviolet absorber in the venusian atmosphere from analysis of Venus Monitoring Camera images

Observations of Venus using the ultraviolet filter of the Venus Monitoring Camera (VMC) on ESA’s Venus Express Spacecraft (VEX) provide the best opportunity for study of the spatial and temporal distribution of the venusian unknown ultraviolet absorber since the Pioneer Venus (PV) mission. We compare the results of two sets of 125 radiative transfer models of the upper atmosphere of Venus to each pixel in a subset of VMC UV channel images. We use a quantitative best fit criterion based upon the notion that the distribution of the unknown absorber should be independent of the illumination and observing geometry. We use the product of the cosines of the incidence and emission angles and search for absorber distributions that are uncorrelated with this geometric parameter, finding that two models can describe the vertical distribution of the unknown absorber. One model is a well-mixed vertical profile above a pressure level of roughly 120 mb (～63 km). This is consistent with the altitude of photochemical formation of sulfuric acid. The second model describes it as a thin layer of pure UV absorber at a pressure level roughly around 24 mb (～71 km) and this altitude is consistent with the top of upper cloud deck. We find that the average abundance of unknown absorber in the equatorial region is 0.21 ± 0.04 optical depth and it decreases in the polar region to 0.08 ± 0.05 optical depth at 365 nm.     

Interpretation of combined infrared,submillimeter, and millimeter thermal flux data obtained during the Rosetta fly-by of Asteroid (21) Lutetia

The European Space Agency’s Rosetta spacecraft is the first Solar System mission to include instrumentation capable of measuring planetary thermal fluxes at both near-IR (VIRTIS) and submillimeter–millimeter (smm–mm, MIRO) wavelengths. Its primary mission is a 1 year reconnaissance of Comet 67P/Churyumov–Gerasimenko beginning in 2014. During a 2010 close fly-by of Asteroid 21 Lutetia, the VIRTIS and MIRO instruments provided complementary data that have been analyzed to produce a consistent model of Lutetia’s surface layer thermal and electrical properties, including a physical model of self-heating. VIRTIS dayside measurements provided highly resolved 1 K accuracy surface temperatures that required a low thermal inertia, I < 30 J/(K m² s^0.5). MIRO smm and mm measurements of polar night thermal fluxes produced constraints on Lutetia’s subsurface thermal properties to depths comparable to the seasonal thermal wave, yielding a model of I < 20 J/(K m² s^0.5) in the upper few centimeters, increasing with depth in a manner very similar to that of Earth’s Moon. Subsequent MIRO-based model predictions of the dayside surface temperatures reveal negative offsets of ～5–30 K from the higher VIRTIS-measurements. By adding surface roughness in the form of 50% fractional coverage of hemispherical mini-craters to the MIRO-based thermal model, sufficient self-heating is produced to largely remove the offsets relative to the VIRTIS measurements and also reproduce the thermal limb brightening features (relative to a smooth surface model) seen by VIRTIS. The Lutetia physical property constraints provided by the VIRTIS and MIRO data sets demonstrate the unique diagnostic capabilities of combined infrared and submillimeter/millimeter thermal flux measurements.     

Spatial and temporal variations in Titan's surface temperatures from Cassini CIRS observations

We report a wide-ranging study of Titan's surface temperatures by analysis of the Moon's outgoing radiance through a spectral window in the thermal infrared at 19 μm (530 cm^?1) characterized by lower atmospheric opacity. We begin by modeling Cassini Composite Infrared Spectrometer (CIRS) far infrared spectra collected in the period 2004–2010, using a radiative transfer forward model combined with a non-linear optimal estimation inversion method. At low-latitudes, we agree with the HASI near-surface temperature of about 94 K at 10°S (Fulchignoni et al., 2005). We find a systematic decrease from the equator toward the poles, hemispherically asymmetric, of ～1 K at 60° south and ～3 K at 60° north, in general agreement with a previous analysis of CIRS data (Jennings et al., 2009), and with Voyager results from the previous northern winter. Subdividing the available database, corresponding to about one Titan season, into 3 consecutive periods, small seasonal changes of up to 2 K at 60°N became noticeable in the results. In addition, clear evidence of diurnal variations of the surface temperatures near the equator are observed for the first time: we find a trend of slowly increasing temperature from the morning to the early afternoon and a faster decrease during the night. The diurnal change is ～1.5 K, in agreement with model predictions for a surface with a thermal inertia between 300 and 600 J m^?2 s^?1/2 K^?1. These results provide important constraints on coupled surface–atmosphere models of Titan's meteorology and atmospheric dynamic.