Vertical cloud structure of Uranus from UKIRT/UIST observations and changes seen during Uranus’ northern spring equinox from 2006 to 2008

Long-slit spectroscopy observations of Uranus by the United Kingdom Infrared Telescope UIST instrument in 2006, 2007 and 2008 have been used to monitor the change in Uranus’ vertical and latitudinal cloud structure through the planet’s northern spring equinox in December 2007.The observed reflectance spectra in the Long J (1.17-1.31 μm) and H (1.45-1.65 μm) bands, obtained with the slit aligned along Uranus’ central meridian, have been fitted with an optimal estimation retrieval model to determine the vertical cloud profile from 0.1 to 6-8 bar over a wide range of latitudes. Context images in a number of spectral bands were used to discriminate general zonal cloud structural changes from passing discrete clouds. From 2006 to 2007 reflection from deep clouds at pressures between 2 and 6-8 bar increased at all latitudes, although there is some systematic uncertainty in the absolute pressure levels resulting from extrapolating the methane coefficients of Irwin et al. (Irwin, P.G.J., Sromovsky, L.A., Strong, E.K., Sihra, K., Teanby, N.A., Bowles, N., Calcutt, S.B., Remedios, J.J. [2006] Icarus, 181, 309-319) at pressures greater than 1 bar, as noted by Tomasko et al. and Karkoschka and Tomasko (Tomasko, M.G., Bezard, B., Doose, L., Engel, S., Karkoschka, E. [2008] Planet. Space Sci., 56, 624-647; Karkoschka, E., Tomasko, M. [2009] Icarus). However, from 2007 to 2008 reflection from these clouds throughout the southern hemisphere and from both northern and southern mid-latitudes (30° N,S) diminished. As a result, the southern polar collar at 45°S has diminished in brightness relative to mid-latitudes, a similar collar at 45°N has become more prominent (e.g. Rages, K.A., Hammel, H.B., Sromovsky, L. [2007] Bull. Am. Astron. Soc., 39, 425; Sromovsky, L.A., Fry, P.M., Ahue, W.M., Hammel, H.B., de Pater, I., Rages, K.A., Showalter, M.R., van Dam, M.A. [2008] vol. 40 of AAS/Division for Planetary Sciences Meeting Abstracts, pp. 488-489; Sromovsky, L.A., Ahue, W.K.M., Fry, P.M., Hammel, H.B., de Pater, I., Rages, K.A., Showalter, M.R. [2009] Icarus), and the lowering reflectivity from mid-latitudes has left a noticeable brighter cloud zone at the equator (e.g. Sromovsky, L.A., Fry, P.M. [2007] Icarus, 192, 527-557;Karkoschka, E., Tomasko, M. [2009] Icarus). For such substantial cloud changes to have occurred in just two years suggests that the circulation of Uranus’ atmosphere is much more vigorous and/or efficient than is commonly thought. The composition of the main observed cloud decks between 2 and 6-8 bar is unclear, but the absence of the expected methane cloud at 1.2-1.3 bar (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987] J. Geophys. Res., 92, 14987-15001) is striking (as previously noted by, among others, Sromovsky, L.A., Irwin, P.G.J., Fry, P.M. [2006] Icarus, 182, 577-593; Sromovsky, L.A., Fry, P.M. [2007] Icarus, 192, 527-557; Sromovsky, L.A., Fry, P.M. [2008] Icarus, 193, 252-266; Karkoschka, E., Tomasko, M. [2009] Icarus) and suggests that cloud particles may be considerably different from pure condensates and may be linked with stratospheric haze particles drizzling down from above, or that tropospheric hazes are generated near the methane condensation level and then drizzle down to deep pressures as suggested by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2009] Icarus).The retrieved cloud structures were also tested for different assumptions of the deep methane mole fraction, which Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2009] Icarus) find may vary from ∼1-2% in polar regions to perhaps as much as 4% equatorwards of 45°N,S. We found that such variations did not significantly affect our conclusions.     

Revised vertical cloud structure of Uranus from UKIRT/UIST observations and changes seen during Uranus’ Northern Spring Equinox from 2006 to 2008: Application of new methane absorption data and comparison with Neptune

Long-slit spectroscopy observations of Uranus by the United Kingdom InfraRed Telescope UIST instrument in 2006, 2007 and 2008 have been used to monitor the change in Uranus’ vertical and latitudinal cloud structure through the planet’s Northern Spring Equinox in December 2007.These spectra were analysed and presented by Irwin et al. (Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2009]. Icarus 203, 287-302), but since publication, a new set of methane absorption data has become available (Karkoschka, E., Tomasko, M. [2010]. Methane absorption coefficients for the jovian planets from laboratory, Huygens, and HST data. Icarus 205, 674-694.), which appears to be more reliable at the cold temperatures and high pressures of Uranus’ deep atmosphere. We have fitted k-coefficients to these new methane absorption data and we find that although the latitudinal variation and inter-annual changes reported by Irwin et al. (2009) stand, the new k-data place the main cloud deck at lower pressures (2-3 bars) than derived previously in the H-band of ∼3-4 bars and ∼3 bars compared with ∼6 bars in the J-band. Indeed, we find that using the new k-data it is possible to reproduce satisfactorily the entire observed centre-of-disc Uranus spectrum from 1 to 1.75 μm with a single cloud at 2-3 bars provided that we make the particles more back-scattering at wavelengths less than 1.2 μm by, for example, increasing the assumed single-scattering albedo from 0.75 (assumed in the J and H-bands) to near 1.0. In addition, we find that using a deep methane mole fraction of 4% in combination with the associated warm ‘F’ temperature profile of Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92, 14987-15001), the retrieved cloud deck using the new (Karkoschka and Tomasko, 2010) methane absorption data moves to between 1 and 2 bars.The same methane absorption data and retrieval algorithm were applied to observations of Neptune made during the same programme and we find that we can again fit the entire 1-1.75 μm centre-of-disc spectrum with a single cloud model, providing that we make the stratospheric haze particles (of much greater opacity than for Uranus) conservatively scattering (i.e. ω = 1) and we also make the deeper cloud particles, again at around the 2 bar level more reflective for wavelengths less than 1.2 μm. Hence, apart from the increased opacity of stratospheric hazes in Neptune’s atmosphere, the deeper cloud structure and cloud composition of Uranus and Neptune would appear to be very similar.     

Jupiter's ammonia clouds—localized or ubiquitous?

From an analysis of the Galileo Near Infrared Imaging Spectrometer (NIMS) data, Baines et al. (Icarus 159 (2002) 74) have reported that spectrally identifiable ammonia clouds (SIACs) cover less than 1% of Jupiter. Localized ammonia clouds have been identified also in the Cassini Composite Infrared Spectrometer (CIRS) observations (Planet. Space Sci. 52 (2004a) 385). Yet, ground-based, satellite and spacecraft observations show that clouds exist everywhere on Jupiter. Thermochemical models also predict that Jupiter must be covered with clouds, with the top layer made up of ammonia ice. For a solar composition atmosphere, models predict the base of the ammonia clouds to be at 720 mb, at 1000 mb if N/H were 4×solar, and at 0.5 bar for depleted ammonia of 10⁻²×solar (Planet. Space Sci. 47 (1999) 1243). Thus, the above NIMS and CIRS findings are seemingly at odds with other observations and cloud physics models. We suggest that the clouds of ammonia ice are ubiquitous on Jupiter, but that spectral identification of all but the freshest of the ammonia clouds and high altitude ammonia haze is inhibited by a combination of (i) dusting, starting with hydrocarbon haze particles falling from Jupiter's stratosphere and combining with an even much larger source—the hydrazine haze; (ii) cloud properties, including ammonia aerosol particle size effects. In this paper, we investigate the role of photochemical haze and find that a substantial amount of haze material can deposit on the upper cloud layer of Jupiter, possibly enough to mask its spectral signature. The stratospheric haze particles result from condensation of polycyclic aromatic hydrocarbons (PAHs), whereas hydrazine ice is formed from ammonia photochemistry. We anticipate similar conditions to prevail on Saturn.     

Limits on the size of aerosols from measurements of linear polarization in Titan’s atmosphere

The Descent Imager/Spectral Radiometer (DISR) instrument on the Huygens probe into the atmosphere of Titan yielded information on the size, shape, optical properties, and vertical distribution of haze aerosols in the atmosphere of Titan [Tomasko, M.G., Doose, L., Engel, S., Dafoe, L.E., West, R., Lemmon, M., Karkoschka, E., 2008. Planet. Space Sci. 56, 669-707] from photometric and spectroscopic measurements of sunlight in Titan’s atmosphere. This instrument also made measurements of the degree of linear polarization of sunlight in two spectral bands centered at 491 and 934 nm. Here we present the calibration and reduction of the polarization measurements and compare the polarization observations to models using fractal aggregate particles which have different sizes for the small dimension (monomer size) of which the aggregates are composed. We find that the Titan aerosols produce very large polarizations perpendicular to the scattering plane for scattering near 90° scattering angle. The size of the monomers is tightly constrained by the measurements to a radius of 0.04 ± 0.01 μm at altitudes from 150 km to the surface. The decrease in polarization with decreasing altitude observed in red and blue light is as expected by increasing dilution due to multiple scattering at decreasing altitudes. There is no indication of particles that produce small amounts of linear polarization at low altitudes.     

Methane on Uranus: The case for a compact CH4 cloud layer at low latitudes and a severe CH4 depletion at high-latitudes based on re-analysis of Voyager occultation measurements and STIS spectroscopy

Lindal et al. (Lindal, G.F., Lyons, J.R., Sweetnam, D.N., Eshleman, V.R., Hinson, D.P. [1987]. J. Geophys. Res. 92 (11), 14987-15001) presented a range of temperature and methane profiles for Uranus that were consistent with 1986 Voyager radio occultation measurements of refractivity versus altitude. A localized refractivity slope variation near 1.2 bars was interpreted to be the result of a condensed methane cloud layer. However, models fit to near-IR spectra found particle concentrations much deeper in the atmosphere, in the 1.5-3 bar range (Sromovsky, L.A., Irwin, P.G.J., Fry, P.M. [2006]. Icarus 182, 577-593; Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229; Irwin, P.G.J., Teanby, N.A., Davis, G.R. [2010]. Icarus 208, 913-926), and a recent analysis of STIS spectra argued for a model in which aerosol particles formed diffusely distributed hazes, with no compact condensation layer (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309). To try to reconcile these results, we reanalyzed the occultation observations with the He volume mixing ratio reduced from 0.15 to 0.116, which is near the edge of the 0.033 uncertainty range given by Conrath et al. (Conrath, B., Hanel, R., Gautier, D., Marten, A., Lindal, G. [1987]. J. Geophys. Res. 92 (11), 15003-15010). This allowed us to obtain saturated mixing ratios within the putative cloud layer and to reach above-cloud and deep methane mixing ratios compatible with STIS spectral constraints. Using a 5-layer vertical aerosol model with two compact cloud layers in the 1-3 bar region, we find that the best fit pressure for the upper compact layer is virtually identical to the pressure range inferred from the occultation analysis for a methane mixing ratio near 4% at 5°S. This strongly argues that Uranus does indeed have a compact methane cloud layer. In addition, our cloud model can fit the latitudinal variations in spectra between 30°S and 20°N, using the same profiles of temperature and methane mixing ratio. But closer to the pole, the model fails to provide accurate fits without introducing an increasingly strong upper tropospheric depletion of methane at increased latitudes, in rough agreement with the trend identified by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M. [2009]. Icarus 202, 287-309).     

Determination of the complex refractive indices of Titan haze analogs using photothermal deflection spectroscopy

The spectrometers of the Cassini mission to the Saturn system have detected haze layers reaching up to 800 km in Titan’s atmosphere. Knowledge of the complex refractive index (k) of the haze is important for modeling the surface and atmosphere of Titan and retrieving some information about the functional groups present in the aerosols. Plasma discharges or ultraviolet radiation are commonly used to drive the formation of solid organics assumed to be good analogs of the Titan aerosols. [Tran, B.N., Ferris, J.P., Chera, J.J., 2003a. The photochemical formation of a Titan haze analog. Structural analysis by X-ray photoelectron and infrared spectroscopy. Icarus 162, 114-124; Tran, B.N., Force, M., Briggs, R., Ferris J.P., Persans, P., Chera, J.J., 2008. Photochemical processes on Titan: Irradiation of mixtures of gases that simulate Titan’s atmosphere. Icarus 177, 106-115] reported the index of refraction of analogs synthesized by far ultraviolet irradiation of various gas mixtures. k was determined in the 200-800 nm wavelength range from transmission and reflection spectroscopy. However, this technique is limited by (i) uncertainties in the absorption values because of the small amounts of organics available, (ii) light scattering by the surface roughness and particulates in the sample. These limitations prompted us to perform new measurements using photothermal deflection spectroscopy (PDS), a technique based on the conversion of absorbed light into heat in the material of interest. By combining traditional spectroscopy (λ < 500 nm) and PDS (λ > 500 nm), we determined values of k over the 375-1550 nm range. k values as low as 10⁻⁴ above 1000 nm were determined. This is one order of magnitude lower than the measurements generally used as a reference for Titan’s aerosols analogs [Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callicott, T.A., Williams, M.W., 1984. Optical-constants of organic Tholins produced in a simulated Titanian atmosphere—from soft-X-ray to microwave-frequencies. Icarus 60(1), 127-137]. We recommend that these results were used in models to describe the optical properties of the aerosols produced in Titan’s stratosphere.     

A physical model of Titan's aerosols

Microphysical simulations of Titan's stratospheric haze show that aerosol microphysics is linked to organized dynamical processes. The detached haze layer may be a manifestation of 1 cm sec-1 vertical velocities at altitudes above 300 km. The hemispherical asymmetry in the visible albedo may be caused by 0.05 cm sec-1 vertical velocities at altitudes of 150 to 200 km, we predict contrast reversal beyond 0.6 micrometer. Tomasko and Smith's (1982, Icarus 51, 65-95) model, in which a layer of large particles above 220 km altitude is responsible for the high forward scattering observed by Rages and Pollack (1983, Icarus 55, 50-62), is a natural outcome of the detached haze layer being produced by rising motions if aerosol mass production occurs primarily below the detached haze layer. The aerosol's electrical charge is critical for the particle size and optical depth of the haze. The geometric albedo, particularly in the ultraviolet and near infrared, requires that the particle size be near 0.15 micrometer down to altitudes below 100 km, which is consistent with polarization observations (Tomasko and Smith 1982, West and Smith 1991, Icarus 90, 330-333). Above about 400 km and below about 150 km Yung et al.'s (1984, Astrophys. J. Suppl. Ser. 55, 465-506) diffusion coefficients are too small. Dynamical processes control the haze particles below about 150 km. The relatively large eddy diffusion coefficients in the lower stratosphere result in a vertically extensive region with nonuniform mixing ratios of condensable gases, so that most hydrocarbons may condense very near the tropopause rather than tens of kilometers above it. The optical depths of hydrocarbon clouds are probably less than one, requiring that abundant gases such as ethane condense on a subset of the haze particles to create relatively large, rapidly removed particles. The wavelength dependence of the optical radius is calculated for use in analyzing observations of the geometric albedo. The lower atmosphere and surface should be visible outside of regions of methane absorption in the near infrared. Limb scans at 2.0 micrometers wavelength should be possible down to about 75 km altitude.     

Deciphering the origin of the regular satellites of gaseous giants - Iapetus: The Rosetta ice-moon

Ever since their discovery the regular satellites of Jupiter and Saturn have held out the promise of providing an independent set of observations with which to test theories of planet formation. Yet elucidating their origins has proven elusive. Here we show that Iapetus can serve to discriminate between satellite formation models. Its accretion history can be understood in terms of a two-component gaseous subnebula, with a relatively dense inner region, and an extended tail out to the location of the irregular satellites, as in the SEMM model of Mosqueira and Estrada (2003a,b) (Mosqueira, I., Estrada, P.R. [2003a]. Icarus 163, 198-231; Mosqueira, I., Estrada, P.R. [2003b]. Icarus 163, 232-255). Following giant planet formation, planetesimals in the feeding zone of Jupiter and Saturn become dynamically excited, and undergo a collisional cascade. Ablation and capture of planetesimal fragments crossing the gaseous circumplanetary disks delivers enough collisional rubble to account for the mass budgets of the regular satellites of Jupiter and Saturn. This process can result in rock/ice fractionation as long as the make up of the population of disk crossers is non-homogeneous, thus offering a natural explanation for the marked compositional differences between outer solar nebula objects and those that accreted in the subnebulae of the giant planets. For a given size, icy objects are easier to capture and to ablate, likely resulting in an overall enrichment of ice in the subnebula. Furthermore, capture and ablation of rocky fragments become inefficient far from the planet for two reasons: the gas surface density of the subnebula is taken to drop outside the centrifugal radius, and the velocity of interlopers decreases with distance from the planet. Thus, rocky objects crossing the outer disks of Jupiter and Saturn never reach a temperature high enough to ablate either due to melting or vaporization, and capture is also greatly diminished there. In contrast, icy objects crossing the outer disks of each planet ablate due to the melting and vaporization of water-ice. Consequently, our model leads to an enhancement of the ice content of Iapetus, and to a lesser degree those of Titan, Callisto and Ganymede, and accounts for the (non-stochastic) compositions of these large, low-porosity outer regular satellites of Jupiter and Saturn. For this to work, the primordial population of planetesimals in the Jupiter-Saturn region must be partially differentiated, so that the ensuing collisional cascade produces an icy population of ?1 m size fragments to be ablated during subnebula crossing. We argue this is likely because the first generation of solar nebula ∼10 km planetesimals in the Jupiter-Saturn region incorporated significant quantities of ²⁶Al. This is the first study successfully to provide a direct connection between nebula planetesimals and subnebulae mixtures with quantifiable and observable consequences for the bulk properties of the regular satellites of Jupiter and Saturn, and the only explanation presently available for Iapetus’ low density and ice-rich composition.     

The source of widespread 3-μm absorption in Jupiter’s clouds: Constraints from 2000 Cassini VIMS observations

The Cassini flyby of Jupiter in 2000 provided spatially resolved spectra of Jupiter’s atmosphere using the Visual and Infrared Mapping Spectrometer (VIMS). A prominent characteristic of these spectra is the presence of a strong absorption at wavelengths from about 2.9 μm to 3.1 μm, previously noticed in a 3-μm spectrum obtained by the Infrared Space Observatory (ISO) in 1996. While Brooke et al. (Brooke, T.Y., Knacke, R.F., Encrenaz, T., Drossart, P., Crisp, D., Feuchtgruber, H. [1998]. Icarus 136, 1-13) were able to fit the ISO spectrum very well using ammonia ice as the sole source of particulate absorption, Sromovsky and Fry (Sromovsky, L.A., Fry, P.M. [2010]. Icarus 210, 211-229), using significantly revised NH₃ gas absorption models, showed that ammonium hydrosulfide (NH₄SH) provided a better fit to the ISO spectrum than NH₃, but that the best fit was obtained when both NH₃ and NH₄SH were present in the clouds. Although the large FOV of the ISO instrument precluded identification of the spatial distribution of these two components, the VIMS spectra at low and intermediate phase angles show that 3-μm absorption is present in zones and belts, in every region investigated, and both low- and high-opacity samples are best fit with a combination of NH₄SH and NH₃ particles at all locations. The best fits are obtained with a layer of small ammonia-coated particles (r ∼ 0.3 μm) overlying but often close to an optically thicker but still modest layer of much larger NH₄SH particles (r ∼ 10 μm), with a deeper optically thicker layer, which might also be composed of NH₄SH. Although these fits put NH₃ ice at pressures less than 500 mb, this is not inconsistent with the lack of prominent NH₃ features in Jupiter’s longwave spectrum because the reflectivity of the core particles strongly suppresses the NH₃ absorption features, at both near-IR and thermal wavelengths. Unlike Jupiter, Saturn lacks the broad 3-μm absorption feature, but does exhibit a small absorption near 2.965 μm, which resembles a similar jovian feature and suggests that both planets contain upper tropospheric clouds of sub-micron particles containing ammonia as a minor fraction.     

Phosphine on Jupiter and Saturn from Cassini/CIRS

The global distribution of phosphine (PH₃) on Jupiter and Saturn is derived using 2.5 cm⁻¹ spectral resolution Cassini/CIRS observations. We extend the preliminary PH₃ analyses on the gas giants [Irwin, P.G.J., and 6 colleagues, 2004. Icarus 172, 37-49; Fletcher, L.N., and 9 colleagues, 2007a. Icarus 188, 72-88] by (a) incorporating a wider range of Cassini/CIRS datasets and by considering a broader spectral range; (b) direct incorporation of thermal infrared opacities due to tropospheric aerosols and (c) using a common retrieval algorithm and spectroscopic line database to allow direct comparison between these two gas giants.The results suggest striking similarities between the tropospheric dynamics in the 100-1000 mbar regions of the giant planets: both demonstrate enhanced PH₃ at the equator, depletion over neighbouring equatorial belts and mid-latitude belt/zone structures. Saturn's polar PH₃ shows depletion within the hot cyclonic polar vortices. Jovian aerosol distributions are consistent with previous independent studies, and on Saturn we demonstrate that CIRS spectra are most consistent with a haze in the 100-400 mbar range with a mean optical depth of 0.1 at 10 μm. Unlike Jupiter, Saturn's tropospheric haze shows a hemispherical asymmetry, being more opaque in the southern summer hemisphere than in the north. Thermal-IR haze opacity is not enhanced at Saturn's equator as it is on Jupiter.Small-scale perturbations to the mean PH₃ abundance are discussed both in terms of a model of meridional overturning and parameterisation as eddy mixing. The large-scale structure of the PH₃ distributions is likely to be related to changes in the photochemical lifetimes and the shielding due to aerosol opacities. On Saturn, the enhanced summer opacity results in shielding and extended photochemical lifetimes for PH₃, permitting elevated PH₃ levels over Saturn's summer hemisphere.     

Characteristics of Saturn’s FUV airglow from limb-viewing spectra obtained with Cassini-UVIS

This study reports the analysis of far ultraviolet (FUV) limb spectra of the airglow of Saturn in the 1150-1850 Å spectral window, obtained with the Ultraviolet Imaging Spectrograph (UVIS) onboard Cassini, spanning altitudes from −1200 to 4000 km. The FUV limb emission consists of three main contributions: (1) H Ly-α peaking at 1100 km with a brightness of 0.8 kilo-Rayleighs (kR), (2) reflected sunlight longward of 1550 Å which maximizes at −950 km with 16.5 kR and (3) H₂ bands in the 1150-1650 Å bandwidth, peaking at 1050 km reaching a maximum of 3.9 kR.A vertical profile of the local H₂ volume emission rate has been derived using the hydrocarbon density profiles from a model of the Saturn equatorial atmosphere. It is well matched by a Chapman function, characterized by a maximum value of 3.5 photons cm⁻³ s⁻¹ in the 800-1650 Å UV bandwidth, peaking at 1020 km.Comparisons between the observed spectra and a first-order synthetic airglow H₂ model in the 1150-1650 Å bandwidth show that the spectral shape of the H₂ bands is accounted for by solar fluorescence and photoelectron excitation. The best fits are obtained with a combination of H₂ fluorescence lines and 20 eV electron impact spectra, the latter contributing ∼68% of the total H₂ airglow emission.     

Accurate and approximate calculations of Raman scattering in the atmosphere of Neptune

Raman scattering by H₂ in Neptune's atmosphere has significant effects on its reflectivity for λ<0.5 μm, producing baseline decreases of ∼20% in a clear atmosphere and ∼10% in a hazy atmosphere. However, few accurate Raman calculations are carried out because of their complexity and computational costs. Here we present the first radiation transfer algorithm that includes both polarization and Raman scattering and facilitates computation of spatially resolved spectra. New calculations show that Cochran and Trafton's (1978, Astrophys. J. 219, 756-762) suggestion that light reflected in the deep CH₄ bands is mainly Raman scattered is not valid for current estimates of the CH₄ vertical distribution, which implies only a 4% Raman contribution. Comparisons with IUE, HST, and groundbased observations confirm that high altitude haze absorption is reducing Neptune's geometric albedo by ∼6% in the 0.22-0.26 μm range and by ∼13% in the 0.35-0.45 μm range. A sample haze model with 0.2 optical depths of 0.2-μm radius particles between 0.1 and 0.8 bars fits reasonably well, but is not a unique solution. We used accurate calculations to evaluate several approximations of Raman scattering. The Karkoschka (1994, Icarus 111, 174-192) method of applying Raman corrections to calculated spectra and removing Raman effects from observed spectra is shown to have limited applicability and to undercorrect the depths of weak CH₄ absorption bands. The relatively large Q-branch contribution observed by Karkoschka is shown to be consistent with current estimates of Raman cross-sections. The Wallace (1972, Astrophys. J. 176, 249-257) approximation, produces geometric albedo ∼5% low as originally proposed, but can be made much more accurate by including a scattering contribution from the vibrational transition. The original Pollack et al. (1986, Icarus 65, 442-466) approximation is inaccurate and unstable, but can be greatly improved by several simple modifications. A new approximation based on spectral tuning of the effective molecular single scattering albedo provides low errors for zenith angles below 70° in a clear atmosphere, although intermediate clouds present problems at longer wavelengths.     

Optical properties of Titan and early Earth haze laboratory analogs in the mid-visible

Scattering and absorption of sunlight by aerosols are integral to understanding the radiative balance of any planetary atmosphere covered in a haze, such as Titan and possibly the early Earth. One key optical parameter of an aerosol is its refractive index. We have simulated both Titan and early Earth organic haze aerosols in the laboratory and measured the real and imaginary portion of their refractive index at λ = 532 nm using cavity ringdown aerosol extinction spectroscopy. This novel technique allows analysis on freely-floating particles minutes after formation. For our Titan analog particles, we find a real refractive index of n = 1.35 ± 0.01 and an imaginary refractive index k = 0.023 ± 0.007, and for the early Earth analog particles we find n = 1.81 ± 0.02 and k = 0.055 ± 0.020. The Titan analog refractive index has a smaller real and similar imaginary refractive index compared to most previous laboratory measurements of Titan analog films, including values from Khare et al. (Khare, B.N., Sagan, C., Arakawa, E.T., Suits, F., Callcott, T.A., Williams, M.W. [1984]. Icarus 60, 127-137). These newly measured Titan analog values have implications for spacecraft retrievals of aerosol properties on Titan. The early Earth analog has a significantly higher real and imaginary refractive index than Titan analogs reported in the literature. These differences suggest that, for a given amount of aerosol, the early Earth analog would act as a stronger anti-greenhouse agent than the Titan analog.     

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.     

Measurements of methane absorption by the descent imager/spectral radiometer (DISR) during its descent through Titan's atmosphere

New low-temperature methane absorption coefficients pertinent to the Titan environment are presented as derived from the Huygens DISR spectral measurements combined with the in-situ measurements of the methane gas abundance profile measured by the Huygens Gas Chromatograph/Mass Spectrometer (GCMS). The visible and near-infrared spectrometers of the descent imager/spectral radiometer (DISR) instrument on the Huygens probe looked upward and downward covering wavelengths from 480 to 1620 nm at altitudes from 150 km to the surface during the descent to Titan's surface. The measurements at continuum wavelengths were used to determine the vertical distribution, single-scattering albedos, and phase functions of the aerosols. The gas chromatograph/mass spectrometer (GCMS) instrument on the probe measured the methane mixing ratio throughout the descent. The DISR measurements are the first direct measurements of the absorbing properties of methane gas made in the atmosphere of Titan at the pathlengths, pressures, and temperatures that occur there. Here we use the DISR spectral measurements to determine the relative methane absorptions at different wavelengths along the path from the probe to the sun throughout the descent. These transmissions as functions of methane path length are fit by exponential sums and used in a haze radiative transfer model to compare the results to the spectra measured by DISR. We also compare the recent laboratory measurements of methane absorption at low temperatures [Irwin et al., 2006. Improved near-infrared methane band models and k-distribution parameters from 2000 to 9500 cm⁻¹ and implications for interpretation of outer planet spectra. Icarus 181, 309-319] with the DISR measurements. We find that the strong bands formed at low pressures on Titan act as if they have roughly half the absorption predicted by the laboratory measurements, while the weak absorption regions absorb considerably more than suggested by some extrapolations of warm measurements to the cold Titan temperatures. We give factors as a function of wavelength that can be used with the published methane coefficients between 830 and 1620 nm to give agreement with the DISR measurements. We also give exponential sum coefficients for methane absorptions that fit the DISR observations. We find the DISR observations of the weaker methane bands shortward of 830 nm agree with the methane coefficients given by Karkoschka [1994. Spectrophotometry of the jovian planets and Titan at 300- to 1000-nm wavelength: the methane spectrum. Icarus 111, 174-192]. Finally, we discuss the implications of our results for computations of methane absorption in the atmospheres of the outer planets.     

Latitudinal variation of Saturn photochemistry deduced from spatially-resolved ultraviolet spectra

We obtained spatially-resolved ultraviolet spectra of Saturn in 1994 with the Faint Object Spectrometer and Goddard High Resolution Spectrograph of the Hubble Space Telescope. We observed four areas on the planet at 15° N, 33° S, 41° S, and 52° S, with a field-of-view of less than 2 × 2 arcsec², compared to the 16-arcsec planet diameter. The wavelength range, 1550-2300 Å, encompasses absorption from major hydrocarbons (C₂H₆, C₂H₄, C₂H₂, CH₃C₂H, C₄H₂) and water. We find global hydrocarbon abundances and a C₂H₂ vertical distribution compatible with infrared observations, in contrast with previous analyses of ultraviolet spectra. The stratospheric haze opacity decreases from polar region to the equator. Saturn mid-latitudes are photochemically distinct from the rest of the planet. At 33° S, the spectrum requires either (1) a distinctly different C₂H₂ vertical distribution or (2) a locally enhanced water abundance. At 41° S, the hydrocarbon abundance exhibits a local minimum, within a global trend of increasing abundance from equator to pole. This global trend may result from an increased abundance of short-lived hydrocarbons such as C₄H₂. Photochemical models predict a depletion of hydrocarbon molecules in the presence of stratospheric water [Moses et al., 2000. Icarus 143, 166-202]. These results are consistent with a localized influx of water, in the form of high charge to mass ratio particles, flowing into Saturn's atmosphere at latitudes magnetically linked to the rings.     

Saturn: UBV photoelectric pinhole scans of the disk

UBV pinhole scans of the Saturn disk have been made with a photoelectric area-scanning photometer. Limb profiles, spaced parallel to the equator, were obtained over the entire southern hemisphere of the planet. Saturn was found to exhibit strong limb brightening in the ultraviolet, moderate limb brightening at blue wavelengths, and strong limb darkening in the visual region of the spectrum. Latitudinal variations in the disk profiles were found. In general, the degree of limb brightening decreases towards the polar region. Pronounced asymmetry is apparent in the disk profiles in each color. The sunward limb is significantly brighter than the opposite limb. This asymmetry depends on phase angle; approaching zero at opposition, it reaches a maximum near quadrature. Our observations are interpreted using an elementary radiative transfer model. The Saturn atmosphere is approximated by a finite homogeneous layer of isotropically scattering particles overlying a Lambert scattering haze or cloud layer. The reflectivity of the haze or clouds is a strongly dependent function of wavelength. Our best-fitting model consists of a clear H₂ layer of column density ～31 km-am above the haze or clouds; the maximum permitted H₂ column density is ～46 km-am. The H₂ column density above the equatorial region appears to be less than at temperate latitudes. The phase-dependent asymmetry in the disk profiles is a natural consequence of the scattering geometry. Our results are consistent with current knowledge of the Saturn atmosphere.     

Jupiter and Saturn rotation periods

Anderson and Schubert [2007. Saturn's Gravitational field, internal rotation, and interior structure. Science 317, 1384-1387 (paper I)] proposed that Saturn's rotation period can be ascertained by minimizing the dynamic heights of the 100 mbar isosurface with respect to the geoid; they derived a rotation period of 10 h 32 m 35 s. We investigate the same approach for Jupiter to see if the Jovian rotation period is predicted by minimizing the dynamical heights of its isobaric (1 bar pressure level) surface using zonal wind data. A rotation period of 9 h 54 m 29.7 s is found. Further, we investigate the minimization method by fitting Pioneer and Voyager occultation radii for both Jupiter and Saturn. Rotation periods of 9 h 55 m 30 s and 10 h 32 m 35 s are found to minimize the dynamical heights for Jupiter and Saturn, respectively. Though there is no dynamical principle requiring the minimization of the dynamical heights of an isobaric surface, the successful application of the method to Jupiter lends support to its relevance for Saturn.We derive Jupiter and Saturn rotation periods using equilibrium theory to explain the difference between equatorial and polar radii. Rotation periods of 9 h 55 m 20 s and 10 h 31 m 49 s are found for Jupiter and Saturn, respectively. We show that both Jupiter's and Saturn's shapes can be derived using solid-body rotation, suggesting that zonal winds have a minor effect on the planetary shape for both planets.The agreement in the values of Saturn's rotation period predicted by the different approaches supports the conclusion that the planet's period of rotation is about 10 h 32 m.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets - Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model”. We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size.