Changes in Jupiter’s zonal velocity between 1979 and 2008

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

Adsorptive fractionation of HDO on JSC MARS-1 during sublimation with implications for the regolith of Mars

A chamber was constructed to simulate the boundary between the ice table, regolith and atmosphere of Mars and to examine fractionation between H₂O and HDO during sublimation under realistic martian conditions of temperature and pressure. Thirteen experimental runs were conducted with regolith overlying the ice. The thickness and characteristic grain size of the regolith layer as well as the temperature of the underlying ice was varied. From these runs, values for the effective diffusivity, taking into account the effects of adsorption, of the regolith were derived. These effective diffusivities ranged from 1.8 × 10⁻⁴ m² s⁻¹ to 2.2 × 10⁻³ m² s⁻¹ for bare ice and from 2.4 × 10⁻¹¹ m² s⁻¹ to 2.0 × 10⁻⁹ m² s⁻¹ with an adsorptive layer present. From these, latent heats of adsorption of 8.6 ± 2.6 kJ mol⁻¹ and 9.3 ± 2.8 kJ mol⁻¹ were derived at ice-surface temperatures above 223 ± 8 K and 96 ± 28 kJ mol⁻¹ and 104 ± 31 kJ mol⁻¹ respectively for H₂O and HDO were derived at colder temperatures. For temperatures below 223 K, the effective diffusivity of HDO was found to be lower than the diffusivity of H₂O by 40% on average, suggesting that the regolith was adsorptively fractionating the sublimating gas with a fractionation factor of 1.96 ± 0.74. Applying these values to Mars predicts that adsorbed water on the regolith is enriched in HDO compared to the atmosphere, particularly where the regolith is colder. Based on current observations, the D/H ratio of the regolith may be as high as 21 ± 8 times VSMOW at 12°S and L_S = 357° if the regolith is hydrated primarily by the atmosphere, neglecting any hydration from subsurface ice.     

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

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

Deep Impact, Stardust-NExT and the behavior of Comet 9P/Tempel 1 from 1997 to 2010

We present observational data for Comet 9P/Tempel 1 taken from 1997 through 2010 in an international collaboration in support of the Deep Impact and Stardust-NExT missions. The data were obtained to characterize the nucleus prior to the Deep Impact 2005 encounter, and to enable us to understand the rotation state in order to make a time of arrival adjustment in February 2010 that would allow us to image at least 25% of the nucleus seen by the Deep Impact spacecraft to better than 80 m/pixel, and to image the crater made during the encounter, if possible. In total, ∼500 whole or partial nights were allocated to this project at 14 observatories worldwide, utilizing 25 telescopes. Seventy percent of these nights yielded useful data. The data were used to determine the linear phase coefficient for the comet in the R-band to be 0.045 ± 0.001 mag deg⁻¹ from 1° to 16°. Cometary activity was observed to begin inbound near r ∼ 4.0 AU and the activity ended near r ∼ 4.6 AU as seen from the heliocentric secular light curves, water-sublimation models and from dust dynamical modeling. The light curve exhibits a significant pre- and post-perihelion brightness and activity asymmetry. There was a secular decrease in activity between the 2000 and 2005 perihelion passages of ∼20%. The post-perihelion light curve cannot be easily explained by a simple decrease in solar insolation or observing geometry. CN emission was detected in the comet at 2.43 AU pre-perihelion, and by r = 2.24 AU emission from C₂ and C₃ were evident. In December 2004 the production rate of CN increased from 1.8 × 10²³ mol s⁻¹ to Q_CN = 2.75 × 10²³ mol s⁻¹ in early January 2005 and 9.3 × 10²⁴ mol s⁻¹ on June 6, 2005 at r = 1.53 AU.     

Dust devil sediment flux on Earth and Mars: Laboratory simulations

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

Venus night airglow: Ground-based detection of OH, observations of O2 emissions, and photochemical model

Venus nightglow was observed at NASA IRTF using a high-resolution long-slit spectrograph CSHELL at LT = 21:30 and 4:00 on Venus. Variations of the O₂ airglow at 1.27 μm and its rotational temperature are extracted from the observed spectra. The mean O₂ nightglow is 0.57 MR at 21:30 at 35°S-35°N, and the temperature increases from 171 K near the equator to ∼200 K at ±35°. We have found a narrow window that covers the OH (1-0) P1(4.5) and (2-1) Q1(1.5) airglow lines. The detected line intensities are converted into the (1-0) and (2-1) band intensities of 7.2 ± 1.8 kR and <1.4 kR at 21:30 and 15.5 ± 2 kR and 4.7 ± 1 kR at 4:00. The f-component of the (1-0) P1(4.5) line has not been detected in either observation, possibly because of resonance quenching in CO₂. The observed Earth’s OH (1-0) and (2-1) bands were 400 and 90 kR at 19:30 and 250 and 65 kR at 9:40, respectively. A photochemical model for the nighttime atmosphere at 80-130 km has been made. The model involves 61 reactions of 24 species, including odd hydrogen and chlorine chemistries, with fluxes of O, N, and H at 130 km as input parameters. To fit the OH vibrational distribution observed by VEX, quenching of OH (v > 3) in CO₂ only to v ? 2 is assumed. According to the model, the nightside-mean O₂ emission of 0.52 MR from the VEX and our observations requires an O flux of 2.9 × 10¹² cm⁻² s⁻¹ which is 45% of the dayside production above 80 km. This makes questionable the nightside-mean O₂ intensities of ∼1 MR from some observations. Bright nightglow patches are not ruled out; however, the mean nightglow is ∼0.5 MR as observed by VEX and supported by the model. The NO nightglow of 425 R needs an N flux of 1.2 × 10⁹ cm⁻² s⁻¹, which is close to that from VTGCM at solar minimum. However, the dayside supply of N at solar maximum is half that required to explain the NO nightglow in the PV observations. The limited data on the OH nightglow variations from the VEX and our observations are in reasonable agreement with the model. The calculated intensities and peak altitudes of the O₂, NO, and OH nightglow agree with the observations. Relationships for the nightglow intensities as functions of the O, N, and H fluxes are derived.     

The near-IR spectrum of Titan modeled with an improved methane line list

We have obtained spatially resolved spectra of Titan in the near-infrared J, H and K bands at a resolving power of ∼5000 using the near-infrared integral field spectrometer (NIFS) on the Gemini North 8 m telescope. Using recent data from the Cassini/Huygens mission on the atmospheric composition and surface and aerosol properties, we develop a multiple-scattering radiative transfer model for the Titan atmosphere. The Titan spectrum at these wavelengths is dominated by absorption due to methane with a series of strong absorption band systems separated by window regions where the surface of Titan can be seen. We use a line-by-line approach to derive the methane absorption coefficients. The methane spectrum is only accurately represented in standard line lists down to ∼2.1 μm. However, by making use of recent laboratory data and modeling of the methane spectrum we are able to construct a new line list that can be used down to 1.3 μm. The new line list allows us to generate spectra that are a good match to the observations at all wavelengths longer than 1.3 μm and allow us to model regions, such as the 1.55 μm window that could not be studied usefully with previous line lists such as HITRAN 2008. We point out the importance of the far-wing line shape of strong methane lines in determining the shape of the methane windows. Line shapes with Lorentzian, and sub-Lorentzian regions are needed to match the shape of the windows, but different shape parameters are needed for the 1.55 μm and 2 μm windows. After the methane lines are modeled our observations are sensitive to additional absorptions, and we use the data in the 1.55 μm region to determine a D/H ratio of 1.77 ± 0.20 × 10⁻⁴, and a CO mixing ratio of 50 ± 11 ppmv. In the 2 μm window we detect absorption features that can be identified with the ν₅ + 3ν₆ and 2ν₃ + 2ν₆ bands of CH₃D.     

Methane evolution from UV-irradiated spacecraft materials under simulated martian conditions: Implications for the Mars Science Laboratory (MSL) mission

Fifteen organic and three inorganic compounds were tested for methane (CH₄) evolution under simulated martian conditions of 6.9 mbar; UVC (200-280 nm) flux of 4 W m⁻²; 20 °C; simulated optical depth of 0.1; and a Mars gas composition of CO₂ (95.3%), N₂ (2.7%), Ar (1.7%), O₂ (0.13%), and water vapor (0.03%). All three inorganic compounds (i.e., NaCl, CaCO₃, graphite) failed to evolve methane at the minimum detection level 0.5 ppm, or above. In contrast, all organic compounds evolved methane when exposed to UV irradiation under simulated martian conditions. The polycyclic aromatic hydrocarbon, pyrene, released the most methane per unit of time at 0.175 nmol CH₄ g⁻¹ h⁻¹, and a spectral reflectance target material used for the MER rovers and Phoenix lander released the least methane at 0.00065 nmol CH₄ cm⁻² h⁻¹. Methane was also released from UV-killed bacterial endospores of Bacillus subtilis. Although all organic compounds evolved methane when irradiated with UV photons under martian conditions, the concentrations of residual organics, biogenic signature molecules, and dead microbial cells should be relatively low on the exterior surfaces of the MSL rover, and, thus, not significant sources of methane contamination. In contrast, kapton tape was found to evolve methane at the rate of 0.00165 nmol CH₄ cm⁻² h⁻¹ (16.5 nmol m⁻² h⁻¹) under the UV and martian conditions tested. Although the evolution of methane from kapton tape was found to decline over time, the large amount of kapton tape used on the MSL rover (lower bound estimated at 3 m²) is likely to create a significant source of terrestrial methane contamination during the early part of the mission.     

The outgassing and composition of Comet 19P/Borrelly from radio observations

Radio spectroscopic observations of Comet 19P/Borrelly were performed during the 1994 apparition and at, and near, the time of the Deep Space 1 flyby in 2001. HCN, CS, CH₃OH, and H₂CO were detected using the 30-m telescope of the Institut de Radioastronomie Millimétrique and the James Clerk Maxwell Telescope, and their production rates relative to water are estimated to be 0.06-0.11, 0.07, 1.7, and 0.4%, respectively. Only upper limits are derived for H₂S and CO. The upper limit for CO/H₂O (<15%) is not very constraining, while the upper limit for the H₂S/H₂O ratio of 0.45% is near the bottom of the range of values measured for other comets. Observations of the OH radical at the Nançay radio telescope provide water production rates a few weeks before the 1994 and 2001 perihelia. Observations of the 1₁₀-1₀₁ water line at 557 GHz with the Odin satellite yield a water production rate of (2.5±0.5)×10²⁸ s⁻¹ on September 22, 2001, at the time of the Deep Space 1 encounter, and (3.3±0.6)×10²⁸ s⁻¹ averaged over the September 22-24, 2001 period. The line shapes are asymmetric and blueshifted by V₀∼−0.18 km s⁻¹ for the best observed HCN lines recorded one week after perihelion. The HCN line shapes, and the similar OH and HCN velocity shifts over the September-November 1994 and August-September 2001 periods, favor anisotropic outgassing towards the Sun. Strong outgassing directed along the primary dust jet seen on visible images is not excluded by the HCN line shapes, but unrealistically high gas expansion velocities are required to explain the line shapes in that case.     

In situ measurements of particle load and transport in dust devils

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

Kinematics of stellar populations with RAVE data

We study the kinematics of the Galactic thin and thick disk populations using stars from the RAVE survey’s second data release together with distance estimates from Breddels et al. (2010). The velocity distribution exhibits the expected moving groups present in the solar neighborhood. We separate thick and thin disk stars by applying the X (stellar-population) criterion of Schuster et al. (1993), which takes into account both kinematic and metallicity information. For 1906 thin disk and 110 thick disk stars classified in this way, we find a vertical velocity dispersion, mean rotational velocity and mean orbital eccentricity of (σ_W, 〈V_Φ〉, 〈e〉)_thin = (18 ± 0.3 km s⁻¹, 223 ± 0.4 km s⁻¹, 0.07 ± 0.07) and (σ_W, 〈V_Φ〉, 〈e〉)_thick = (35 ± 2 km s⁻¹, 163 ± 3 km s⁻¹, 0.31 ± 0.16), respectively. From the radial Jeans equation, we derive a thick disk scale length in the range 1.5-2.2 kpc, whose greatest uncertainty lies in the adopted form of the underlying potential. The shape of the orbital eccentricity distribution indicates that the thick disk stars in our sample most likely formed in situ with minor gas-rich mergers and/or radial migration being the most likely cause for their orbits. We further obtain mean metal abundances of 〈[M/H]〉_thin = +0.03 ± 0.17, and 〈[M/H]〉_thick = −0.51 ± 0.23, in good agreement with previous estimates. We estimate a radial metallicity gradient in the thin disk of −0.07 dex kpc⁻¹, which is larger than predicted by chemical evolution models where the disk grows inside-out from infalling gas. It is, however, consistent with models where significant migration of stars shapes the chemical signature of the disk, implying that radial migration might play at least part of a role in the thick disk’s formation.     

Numerical modelling of heating in porous planetesimal collisions

Collisions between planetesimals at speeds of several kilometres per second were common during the early evolution of our Solar System. However, the collateral effects of these collisions are not well understood. In this paper, we quantify the efficiency of heating during high-velocity collisions between planetesimals using hydrocode modelling. We conducted a series of simulations to test the effect on shock heating of the initial porosity and temperature of the planetesimals, the relative velocity of the collision and the relative size of the two colliding bodies. Our results show that while heating is minor in collisions between non-porous planetesimals at impact velocities below 10 km s⁻¹, in agreement with previous work, much higher temperatures are reached in collisions between porous planetesimals. For example, collisions between nearly equal-sized, porous planetesimals can melt all, or nearly all, of the mass of the bodies at collision velocities below 7 km s⁻¹. For collisions of small bodies into larger ones, such as those with an impactor-to-target mass ratio below 0.1, significant localised heating occurs in the target body. At impact velocities as low as 5 km s⁻¹, the mass of melt will be nearly double the mass of the impactor, and the mass of material shock heated by 100 K will be nearly 10 times the mass of the impactor. We present a first-order estimate of the cumulative effects of impact heating on a porous planetesimal parent body by simulating the impact of a population of small bodies until a disruptive event occurs. Before disruption, impact heating is volumetrically minor and highly localised; in no case was more than about 3% of the parent body heated by more than 100 K. However, heating during the final disruptive collision can be significant; in about 10% of cases, almost all of the parent body is heated to 700 K (from an initial temperature of ∼300 K) and more than a tenth of the parent body mass is melted. Hence, energetic collisions between planetesimals could have had important effects on the thermal evolution of primitive materials in the early Solar System.     

New laboratory measurements of CH4 in Titan's conditions and a reanalysis of the DISR near-surface spectra at the Huygens landing site

Laboratory spectra of methane-nitrogen mixtures have been recorded in the near-infrared range (1.0-1.65 μm) in conditions similar to Titan's near surface, to facilitate the interpretation of the DISR/DLIS (DISR—Descent Imager/Spectral Radiometer) spectra taken during the last phase of the descent of the Huygens Probe, when the surface was illuminated by a surface-science lamp. We used a 0.03 cm⁻¹ spectral resolution, adequate to resolve the lines at high pressure (p_N2∼1.5 bar). By comparing the laboratory spectra with synthetic calculations in the well-studied ν₂+2ν₃ band (7515-7620 cm⁻¹), we determine a methane absorption column density of 178±20 cm atm and a temperature of 118±10 K in our experiment. From this, we derive the methane absorption coefficients over 1.0-1.65 μm with a 0.03 cm⁻¹ sampling, allowing for the extrapolation of the results to any other methane column density under the relevant pressure and temperature conditions. We then revisit the calibration and analysis of the Titan “lamp-on” DLIS spectra. We infer a 5.1±0.8% methane-mixing ratio in the first 25 m of Titan's atmosphere. The CH₄ mixing ratio measured 90 s after landing from a distance of 45 cm is found to be 0.92±0.25 times this value, thus showing no post-landing outgassing of methane in excess of ∼20%. Finally, we determine the surface reflectivity as seen between 25 m and 45 cm and find that the 1500 nm absorption band is deeper in the post-landing spectrum as compared to pre-landing.     

A spatially resolved high spectral resolution study of Neptune’s stratosphere

Using TEXES, the Texas Echelon cross Echelle Spectrograph, mounted on the Gemini North 8-m telescope we have mapped the spatial variation of H₂, CH₄, C₂H₂ and C₂H₆ thermal-infrared emission of Neptune. These high-spectral-resolution, spatially resolved, thermal-infrared observations of Neptune offer a unique glimpse into the state of Neptune’s stratosphere in October 2007, L_S = 275.4° just past Neptune’s southern summer solstice (L_S = 270°). We use observations of the S(1) pure rotational line of molecular hydrogen and a portion of the ν₄ band of methane to retrieve detailed information on Neptune’s stratospheric vertical and meridional thermal structure. We find global-average temperatures of 163.8 ± 0.8, 155.0 ± 0.9, and 123.8 ± 0.8 K at the 7.0 × 10⁻³-, 0.12-, and 2.1-mbar levels with no meridional variations within the errors. We then use the inferred temperatures to model the emission of C₂H₂ and C₂H₆ in order to derive stratospheric volume mixing ratios (hence forth, VMR) as a function of pressure and latitude. There is a subtle meridional variation of the C₂H₂ VMR at the 0.5-mbar level with the peak abundance found at −28° latitude, falling off to the north and south. However, the observations are consistent within error to a meridionally constant C₂H₂ VMR of at 0.5 mbar. We find that the VMR of C₂H₆ at 1-mbar peaks at the equator and falls by a factor of 1.6 at −70° latitude. However, a meridionally constant VMR of at the 1-mbar level for C₂H₆ is also statistically consistent with the retrievals. Temperature predictions from a radiative-seasonal climate model of Neptune that assumes the hydrocarbon abundances inferred in this paper are lower than the measured temperatures by 40 K at 7 × 10⁻³ mbar, 30 K at 0.12 mbar and 25 K at 2.1 mbar. The radiative-seasonal model also predicts meridional temperature variations on the order of 10 K from equator to pole, which are not observed. Assuming higher stratospheric CH₄ abundance at the equator relative to the south pole would bring the meridional trends of the inferred temperatures and radiative-seasonal model into closer agreement.We have also retrieved observations of C₂H₄ emission from Neptune’s stratosphere using TEXES on the NASA Infrared Telescope Facility (IRTF) in June 2003, L_S = 266°. Using the observations from the middle of the planet and an average of the middle three latitude temperature profiles from the 2007 observations (9.5° of L_S later, the seasonal equivalent of 9.5 Earth days within Earth’s seasonal cycle), we infer a C₂H₄ VMR of at 1.5 × 10⁻³ mbar, a value that is 3.25 times that predicted by global-average photochemical models.     

Titan's aerosols: Comparison between our model and DISR findings

Our model [Dimitrov, V., Bar-Nun, A., 1999. A model of energy dependent agglomeration of hydrocarbon aerosol particles and implication to Titan's aerosol. J. Aerosol. Sci. 30(1), 35-49] describes the experimentally found polymerization of C₂H₂ and HCN to form aerosol embryos, their growth and adherence to form various aerosols objects [Bar-Nun, A., Kleinfeld, I., Ganor, E., 1988. Shape and optical properties of aerosols formed by photolysis of C₂H₂, C₂H₄ and HCN. J. Geophys. Res. 93, 8383-8387]. These loose fractal objects describe well the findings of DISR on the Huygens probe [Tomasko, M.G., Bézard, B., Doose, L., Engel, S., Karkoschka, E., 2008. Measurements of methane absorption by the descent imager/spectral radiometer (DISR) during its descent through Titan's atmosphere. Planet. Space Sci., this issue, doi:10.1016/j.pss.2007]. These include (1) various regular objects of R=(0.035-0.064)×10⁻⁶ m, as compared with DISR's 0.05×10⁻⁶ m; (2) diverse low and high fractal structures composed of random combinations of various regular and irregular objects; (3) the number density of fractal particles is 6.9×10⁶ m⁻³ at Z=100 km, as compared with DISR's finding of 5.0×10⁶ m⁻³ at Z=80 km; (4) the number of structural units per higher fractals in the atmosphere at Z∼100 km is (2400-2700), as compared with DISR's 3000, and their size being of R=(5.4-6.4)×10⁻⁶ m will satisfy this value and (5) condensation of CH₄ on the highly fractal structures could begin at the altitude where thin methane clouds were observed, filling somewhat the new open fractal structures.     

Nitrile gas chemistry in Titan’s atmosphere

This work presents the first study of the gaseous products resulting from the partial dissociation of methane and nitrogen in the PAMPRE experimental setup simulating Titan’s atmospheric chemistry.Using cryogenic trapping, the gaseous products generated from the chemical reactions occurring in the reactor have been trapped. Analyses of these products by gas chromatography coupled to mass spectrometry have allowed the detection and identification of more than 30 reaction products. Most of them are identified as nitrile species, accompanied by aliphatic hydrocarbons and a few aromatics compounds. The observed species are in agreement with the data from the recent Cassini-Huygens mission as well as from other laboratory setups capable of dissociating nitrogen and methane. This work emphasizes the probable importance of nitrogen-bearing compounds in the chemistry taking place in Titan’s atmosphere.Furthermore, a quantification of mono-nitriles with saturated alkyl chains has been performed relatively to hydrogen cyanide and shows a power law dependence in their concentration. This dependence is consistent with the Cassini-INMS data and Titan’s photochemical models.An empirical relationship has been extracted from our experimental data: [C_xH_2x−1N] = 100x⁻⁵, where x is the number of carbon atoms in the nitrile molecule. This relationship can be directly used in order to foretell the concentration of heavier nitriles induced by chemistry in Titan’s atmosphere.     

Kinetic studies at room temperature of the cyanide anion CN with cyanoacetylene (HC3N) reaction

Rate coefficient of the cyanide anion (CN⁻) with cyanoacetylene (HC₃N) reaction, has been studied in gas phase at room temperature using a Flowing Afterglow Langmuir Probe - Mass Spectrometer (FALP-MS) apparatus. The rate constant for the CN⁻ + HC₃N reaction is k = 4.8 × 10⁻⁹ cm³/s with an uncertainty of 30%.     

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in three primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations (Hubickyj, O., Bodenheimer, P., Lissauer, J.J. [2005]. Icarus 179, 415-431; Lissauer, J.J., Hubickyj, O., D’Angelo, G., Bodenheimer, P. [2009]. Icarus 199, 338-350). These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 M_⊕, for solid surface densities, σ, of 10, 6, and 4 g cm⁻², respectively. For σ = 10 and σ = 6 g cm⁻², respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.     

Calculations of the thermochemistry of six reactions leading to ammonia formation in Titan's atmosphere

The thermochemical properties of the six reactions: (1) N₂+hν (solar EUV) → N⁺ + N(⁴S) + e⁻, (2) N⁺ + H₂ → NH⁺ + H, (3) NH⁺ + H₂ → NH⁺₂ + H, (4) NH⁺₂ + H₂ → NH⁺₃ + H, (5) NH⁺₃ + H₂ → NH⁺₄ + H, and (6) NH⁺₄ + e⁻ → NH₃ + H, were theoretically proposed by Atreya in 1986 and were cited in 2003 by Bernard who assumed that this chain reaction would lead to ammonia formation in Titan's atmosphere. The thermochemical properties of these six reactions have been calculated by means of the coupled cluster singles and doubles (CCSD) at the CCSD/cc-pvdz level, and the CCSD/6-311++g(3df,3pd) level, and G2 method. The geometries of the reactants and products of reactions have been optimized, the energies of reactions have been computed. The analysis of the results shows that: (I) The free energies of four reactions among these six reactions are negative. It means that these reactions, namely reactions (1)-(6) except reaction (2), can react spontaneously in Titan's low temperature environment. The converted temperatures of reactions (3) and (5) are 11881.7 and 4596.9 K, respectively. (II) Reaction (2) is an endothermic reaction, its converted temperature is 1797.6 K. When T<1797.6 K, reaction (2) cannot react forward spontaneously. The barrier of reaction (2) is 26.154 kcal mol⁻¹, which is probably too high to allow it to occur in the atmosphere of Titan. The rate for this reaction at 300 K has been calculated, and the value is k=4.16×¹⁰⁻⁷ s⁻¹. (III) The results of the three methods are more or less the same. So it is concluded that this chain reaction cannot be a pathway to lead to ammonia (gas phase) formation in Titan's atmosphere.