Temporal variations of Titan’s middle-atmospheric temperatures from 2004 to 2009 observed by Cassini/CIRS

We use five and one-half years of limb- and nadir-viewing temperature mapping observations by the Composite Infrared Radiometer-Spectrometer (CIRS) on the Cassini Saturn orbiter, taken between July 2004 and December 2009 (L_S from 293° to 4°; northern mid-winter to just after northern spring equinox), to monitor temperature changes in the upper stratosphere and lower mesosphere of Titan. The largest changes are in the northern (winter) polar stratopause, which has declined in temperature by over 20 K between 2005 and 2009. Throughout the rest of the mid to upper stratosphere and lower mesosphere, temperature changes are less than 5 K. In the southern hemisphere, temperatures in the middle stratosphere near 1 mbar increased by 1-2 K from 2004 through early 2007, then declined by 2-4 K throughout 2008 and 2009, with the changes being larger at more polar latitudes. Middle stratospheric temperatures at mid-northern latitudes show a small 1-2 K increase from 2005 through 2009. At north polar latitudes within the polar vortex, temperatures in the middle stratosphere show a ∼4 K increase during 2007, followed by a comparable decrease in temperatures in 2008 and into early 2009. The observed temperature changes in the north polar region are consistent with a weakening of the subsidence within the descending branch of the middle atmosphere meridional circulation.     

Formation of the regular satellites of giant planets in an extended gaseous nebula I: subnebula model and accretion of satellites

We model the subnebulae of Jupiter and Saturn wherein satellite accretion took place. We expect each giant planet subnebula to be composed of an optically thick (given gaseous opacity) inner region inside of the planet’s centrifugal radius (where the specific angular momentum of the collapsing giant planet gaseous envelope achieves centrifugal balance, located at rCJ ∼ 15R_J for Jupiter and rCS ∼ 22R_S for Saturn) and an optically thin, extended outer disk out to a fraction of the planet’s Roche-lobe (R_H), which we choose to be ∼R_H/5 (located at ∼150 R_J near the inner irregular satellites for Jupiter, and ∼200R_S near Phoebe for Saturn). This places Titan and Ganymede in the inner disk, Callisto and Iapetus in the outer disk, and Hyperion in the transition region. The inner disk is the leftover of the gas accreted by the protoplanet. The outer disk may result from the nebula gas flowing into the protoplanet during the time of giant planet gap-opening (or cessation of gas accretion). For the sake of specificity, we use a solar composition “minimum mass” model to constrain the gas densities of the inner and outer disks of Jupiter and Saturn (and also Uranus). Our model has Ganymede at a subnebula temperature of ∼250 K and Titan at ∼100 K. The outer disks of Jupiter and Saturn have constant temperatures of 130 and 90 K, respectively.Our model has Callisto forming in a time scale ∼10⁶ years, Iapetus in 10⁶-10⁷ years, Ganymede in 10³-10⁴ years, and Titan in 10⁴-10⁵ years. Callisto takes much longer to form than Ganymede because it draws materials from the extended, low density portion of the disk; its accretion time scale is set by the inward drift times of satellitesimals with sizes 300-500 km from distances ∼100R_J. This accretion history may be consistent with a partially differentiated Callisto with a ∼300-km clean ice outer shell overlying a mixed ice and rock-metal interior as suggested by Anderson et al. (2001), which may explain the Ganymede-Callisto dichotomy without resorting to fine-tuning poorly known model parameters. It is also possible that particulate matter coupled to the high specific angular momentum gas flowing through the gap after giant planet gap-opening, capture of heliocentric planetesimals by the extended gas disk, or ablation of planetesimals passing through the disk contributes to the solid content of the disk and lengthens the time scale for Callisto’s formation. Furthermore, this model has Hyperion forming just outside Saturn’s centrifugal radius, captured into resonance by proto-Titan in the presence of a strong gas density gradient as proposed by Lee and Peale (2000). While Titan may have taken significantly longer to form than Ganymede, it still formed fast enough that we would expect it to be fully differentiated. In this sense, it is more like Ganymede than like Callisto (Saturn’s analog of Callisto, we expect, is Iapetus). An alternative starved disk model whose satellite accretion time scale for all the regular satellites is set by the feeding of planetesimals or gas from the planet’s Roche-lobe after gap-opening is likely to imply a long accretion time scale for Titan with small quantities of NH₃ present, leading to a partially differentiated (Callisto-like) Titan. The Cassini mission may resolve this issue conclusively. We briefly discuss the retention of elements more volatile than H₂O as well as other issues that may help to test our model.     

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.     

Titan impacts and escape

We report on hydrodynamic calculations of impacts of large (multi-kilometer) objects on Saturn’s moon Titan. We assess escape from Titan, and evaluate the hypothesis that escaping ejecta blackened the leading hemisphere of Iapetus and peppered the surface of Hyperion.We carried out two- and three-dimensional simulations of impactors ranging in size from 4 to 100 km diameter, impact velocities between 7 and 15 km s⁻¹, and impact angles from 0° to 75° from the vertical. We used the ZEUSMP2 hydrocode for the calculations. Simulations were made using three different geometries: three-dimensional Cartesian, two-dimensional axisymmetric spherical polar, and two-dimensional plane polar. Three-dimensional Cartesian geometry calculations were carried out over a limited domain (e.g. 240 km on a side for an impactor of size d_i = 10 km), and the results compared to ones with the same parameters done by Artemieva and Lunine (2005); in general the comparison was good. Being computationally less demanding, two-dimensional calculations were possible for much larger domains, covering global regions of the satellite (from 800 km below Titan’s surface to the exobase altitude 1700 km above the surface). Axisymmetric spherical polar calculations were carried out for vertical impacts. Two-dimensional plane-polar geometry calculations were made for both vertical and oblique impacts. In general, calculations among all three geometries gave consistent results.Our basic result is that the amount of escaping material is less than or approximately equal to the impactor mass even for the most favorable cases. Amounts of escaping material scaled most strongly as a function of velocity, with high-velocity impacts generating the largest amount, as expected. Dependence of the relative amount of escaping mass f_esc = m_esc/M_i on impactor diameter d_i was weak. Oblique impacts (impact angle θ_i > 45°) were more effective than vertical or near-vertical impacts; ratios of m_esc/M_i ∼ 1-2 were found in the simulations.     

Trapping and dynamical evolution of interplanetary dust particles in Earth’s quasi-satellite resonance

We used numerical simulations to model the orbital evolution of interplanetary dust particles (IDPs) evolving inward past Earth’s orbit under the influence of radiation pressure, Poynting–Robertson light drag (PR drag), solar wind drag, and gravitational perturbations from the planets. A series of β values (where β is the ratio of the force from radiation pressure to that of central gravity) were used ranging from 0.0025 up to 0.02. Assuming a composition consistent with astronomical silicate and a particle density of 2.5 g cm⁻³ these β values correspond to dust particle diameters ranging from 200 μm down to 25 μm. As the dust particle orbits decay past 1 AU between 4% (for β = 0.02, or 25 μm) and 40% (for β = 0.0025, or 200 μm) of the population became trapped in 1:1 co-orbital resonance with Earth. In addition to traditional horseshoe type co-orbitals, we found about a quarter of the co-orbital IDPs became trapped as so-called quasi-satellites. Quasi-satellite IDPs always remain relatively near to Earth (within 0.1–0.3 AU, or 10–30 Hill radii, R_H) and undergo two close-encounters with Earth each year. While resonant perturbations from Earth halt the decay in semi-major axis of quasi-satellite IDPs their orbital eccentricities continue to decrease under the influence of PR drag and solar wind drag, forcing the IDPs onto more Earth-like orbits. This has dramatic consequences for the relative velocity and distance of closest approach between Earth and the quasi-satellite IDPs. After 10⁴–10⁵ years in the quasi-satellite resonance dust particles are typically less than 10R_H from Earth and consistently coming within about 3R_H. In the late stages of evolution, as the dust particles are escaping the 1:1 resonance, quasi-satellite IDPs can have deep close-encounters with Earth significantly below R_H. Removing the effects of Earth’s gravitational acceleration reveals that encounter velocities (i.e., velocities “at infinity”) between quasi-satellite IDPs and Earth during these close-encounters are just a few hundred meters per second or slower, well below the average values of 2–4 km s⁻¹ for non-resonant Earth-crossing IDPs with similar initial orbits. These low encounter velocities lead to a factor of 10–100 increase in Earth’s gravitationally enhanced impact cross-section (σ_grav) for quasi-satellite IDPs compared to similar non-resonant IDPs. The enhancement in σ_grav between quasi-satellite IDPs and cometary Earth-crossing IDPs is even more pronounced, favoring accretion of quasi-satellite dust particles by a factor of 100–3000 over the cometary IDPs. This suggests that quasi-satellite dust particles may dominate the flux of large (25–200 μm) IDPs entering Earth’s atmosphere. Furthermore, because quasi-satellite trapping is known to be directly correlated with the host planet’s orbital eccentricity the accretion of quasi-satellite dust likely ebbs and flows on 10⁵ year time scales synchronized with Earth’s orbital evolution.     

Ionospheric photoelectrons: Comparing Venus, Earth, Mars and Titan

The sunlit portion of planetary ionospheres is sustained by photoionization. This was first confirmed using measurements and modelling at Earth, but recently the Mars Express, Venus Express and Cassini-Huygens missions have revealed the importance of this process at Mars, Venus and Titan, respectively. The primary neutral atmospheric constituents involved (O and CO₂ in the case of Venus and Mars, O and N₂ in the case of Earth and N₂ in the case of Titan) are ionized at each object by EUV solar photons. This process produces photoelectrons with particular spectral characteristics. The electron spectrometers on Venus Express and Mars Express (part of ASPERA-3 and 4, respectively) were designed with excellent energy resolution (ΔE/E=8%) specifically in order to examine the photoelectron spectrum. In addition, the Cassini CAPS electron spectrometer at Saturn also has adequate resolution (ΔE/E=16.7%) to study this population at Titan. At Earth, photoelectrons are well established by in situ measurements, and are even seen in the magnetosphere at up to 7R_E. At Mars, photoelectrons are seen in situ in the ionosphere, but also in the tail at distances out to the Mars Express apoapsis (∼3R_M). At both Venus and Titan, photoelectrons are seen in situ in the ionosphere and in the tail (at up to 1.45R_V and 6.8R_T, respectively). Here, we compare photoelectron measurements at Earth, Venus, Mars and Titan, and in particular show examples of their observation at remote locations from their production point in the dayside ionosphere. This process is found to be common between magnetized and unmagnetized objects. We discuss the role of photoelectrons as tracers of the magnetic connection to the dayside ionosphere, and their possible role in enhancing ion escape.     

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.     

The present-day thermal state of Mars

The present-day thermal state of the martian interior is a very important issue for understanding the internal evolution of the planet. Here, in order to obtain an improved upper limit for the heat flow at the north polar region, we use the lower limit of the effective elastic thickness of the lithosphere loaded by the north polar cap, crustal heat-producing elements (HPE) abundances based on martian geochemistry, and a temperature-dependent thermal conductivity for the upper mantle. We also perform similar calculations for the south polar region, although uncertainties in lithospheric flexure make the results less robust. Our results show that the present-day surface and sublithospheric heat flows cannot be higher than 19 and 12 mW m⁻², respectively, in the north polar region, and similar values might be representative of the south polar region (although with a somewhat higher surface heat flow due to the radioactive contribution from a thicker crust). These values, if representative of martian averages, do not necessarily imply sub-chondritic HPE bulk abundances for Mars (as previously suggested), since (1) chondritic composition models produce a present-day total heat power equivalent to an average surface heat flow of 14-22 mW m⁻² and (2) some convective models obtain similar heat flows for the present time. Regions of low heat flow may even have existed during the last billions of years, in accordance with several surface heat flow estimates of ∼20 mW m⁻² or less for terrains loaded during Hesperian or Amazonian times. On the other hand, there are some evidences suggesting the current existence of regions of enhanced heat flow, and therefore average heat flows could be higher than those obtained for the north (and maybe the south) polar region.     

The transition from complex crater to peak-ring basin on the Moon: New observations from the Lunar Orbiter Laser Altimeter (LOLA) instrument

Impact craters on planetary bodies transition with increasing size from simple, to complex, to peak-ring basins and finally to multi-ring basins. Important to understanding the relationship between complex craters with central peaks and multi-ring basins is the analysis of protobasins (exhibiting a rim crest and interior ring plus a central peak) and peak-ring basins (exhibiting a rim crest and an interior ring). New data have permitted improved portrayal and classification of these transitional features on the Moon. We used new 128 pixel/degree gridded topographic data from the Lunar Orbiter Laser Altimeter (LOLA) instrument onboard the Lunar Reconnaissance Orbiter, combined with image mosaics, to conduct a survey of craters >50 km in diameter on the Moon and to update the existing catalogs of lunar peak-ring basins and protobasins. Our updated catalog includes 17 peak-ring basins (rim-crest diameters range from 207 km to 582 km, geometric mean = 343 km) and 3 protobasins (137-170 km, geometric mean = 157 km). Several basins inferred to be multi-ring basins in prior studies (Apollo, Moscoviense, Grimaldi, Freundlich-Sharonov, Coulomb-Sarton, and Korolev) are now classified as peak-ring basins due to their similarities with lunar peak-ring basin morphologies and absence of definitive topographic ring structures greater than two in number. We also include in our catalog 23 craters exhibiting small ring-like clusters of peaks (50-205 km, geometric mean = 81 km); one (Humboldt) exhibits a rim-crest diameter and an interior morphology that may be uniquely transitional to the process of forming peak rings. A power-law fit to ring diameters (D_ring) and rim-crest diameters (D_r) of peak-ring basins on the Moon [D_ring = 0.14 ± 0.10(D_r)^1.21±0.13] reveals a trend that is very similar to a power-law fit to peak-ring basin diameters on Mercury [D_ring = 0.25 ± 0.14(D_rim)^1.13±0.10] [Baker, D.M.H. et al. [2011]. Planet. Space Sci., in press]. Plots of ring/rim-crest ratios versus rim-crest diameters for peak-ring basins and protobasins on the Moon also reveal a continuous, nonlinear trend that is similar to trends observed for Mercury and Venus and suggest that protobasins and peak-ring basins are parts of a continuum of basin morphologies. The surface density of peak-ring basins on the Moon (4.5 × 10⁻⁷ per km²) is a factor of two less than Mercury (9.9 × 10⁻⁷ per km²), which may be a function of their widely different mean impact velocities (19.4 km/s and 42.5 km/s, respectively) and differences in peak-ring basin onset diameters. New calculations of the onset diameter for peak-ring basins on the Moon and the terrestrial planets re-affirm previous analyses that the Moon has the largest onset diameter for peak-ring basins in the inner Solar System. Comparisons of the predictions of models for the formation of peak-ring basins with the characteristics of the new basin catalog for the Moon suggest that formation and modification of an interior melt cavity and nonlinear scaling of impact melt volume with crater diameter provide important controls on the development of peak rings. In particular, a power-law model of growth of an interior melt cavity with increasing crater diameter is consistent with power-law fits to the peak-ring basin data for the Moon and Mercury. We suggest that the relationship between the depth of melting and depth of the transient cavity offers a plausible control on the onset diameter and subsequent development of peak-ring basins and also multi-ring basins, which is consistent with both planetary gravitational acceleration and mean impact velocity being important in determining the onset of basin morphological forms on the terrestrial planets.     

MARCI and MOC observations of the atmosphere and surface cap in the north polar region of Mars

We used MGS-MOC and MRO-MARCI daily mapping images of the North Polar Region of Mars from 16 August 2005 (L_s = 270°) to 21 May 2009 (L_s = 270°), covering portions of three consecutive martian years (MY 27-MY 29), to observe the seasonal behavior of the polar ice cap and atmospheric phenomena. The rate of cap regression was similar in MY 28 and MY 29, but was advanced by 3.5° of L_s (∼7-8 sols) in MY 29. The spatial and temporal behaviors of dust and condensate clouds were similar in the two years and generally in accord with prior years. Dust storms (>100 km²) were observed in all seasons, with peak activity occurring at L_s = 10-20° from 50°N to 70°N and at L_s = 135-140° from 70°N to 90°N. The most active quadrant was 0-90°W in MY 28, shifting to 180-270°W in MY 29. The majority of regional storms in both years developed in longitudes from 10°W to 60°W. During late summer the larger storms obscure the North Polar Region in a cloud of dust that transitions to north polar hood condensate clouds around autumnal equinox.Changes in the distribution of perennial ice deposits, especially in Olympia Planum, were observed between the 2 years, with the MY 29 ice distribution being the most extensive observed to date. Modeling suggests that the small, bright ice patches on the residual cap are not the result of slope or elevation effects. Rather we suggest that they are the result of local meteorological effects on ice deposition. The annual darkening and brightening of peripheral areas of the residual cap around summer solstice can be explained by the sublimation of a brighter frost layer revealing an underlying darker, ice rich layer that itself either sublimes to reveal brighter material below or acts as a cold trap, attracting condensation of water vapor that brightens the surface. An alternative explanation invokes transport and deposition of dust on the surface from the cap interior, and later removal of that dust. The decrease in cap albedo and accompanying increase in near surface atmospheric stability may be related to the annual minimum of polar storm activity near northern summer solstice.     

Do planetary encounters reset surfaces of near Earth asteroids?

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these ‘space weathering’ processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r^∗, at which planetary encounters can reset the surface. We found that t_sw ∼ 10⁶ yr (at 1 AU) and r^∗ ∼ 5R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 10⁵ yr would require that r^∗ > 20R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 10⁷ yr. We found that t_sw ∝ q², where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.     

Nonhydrostatic effects and the determination of icy satellites’ moment of inertia

We compare the moment of inertia (MOI) of a simple hydrostatic, two layer body as determined by the Radau–Darwin Approximation (RDA) to its exact hydrostatic MOI calculated to first order in the parameter q = Ω²R³/GM, where Ω, R, and M are the spin angular velocity, radius, and mass of the body, and G is the gravitational constant. We show that the RDA is in error by less than 1% for many configurations of core sizes and layer densities congruent with those of solid bodies in the Solar System. We then determine the error in the MOI of icy satellites calculated with the RDA due to nonhydrostatic effects by using a simple model in which the core and outer shell have slight degree 2 distortions away from their expected hydrostatic shapes. Since the hydrostatic shape has an associated stress of order ρΩ²R² (where ρ is density) it follows that the importance of nonhydrostatic effects scales with the dimensionless number σ/ρΩ²R², where σ is the nonhydrostatic stress. This highlights the likely importance of this error for slowly rotating bodies (e.g., Titan and Callisto) and small bodies (e.g., Saturn moons other than Titan). We apply this model to Titan, Callisto, and Enceladus and find that the RDA-derived MOI can be 10% greater than the actual MOI for nonhydrostatic stresses as small as ∼0.1 bars at the surface or ∼1 bar at the core–mantle boundary, but the actual nonhydrostatic stresses for a given shape change depends on the specifics of the interior model. When we apply this model to Ganymede we find that the stresses necessary to produce the same MOI errors as on Titan, Callisto, and Enceladus are an order of magnitude greater due to its faster rotation, so Ganymede may be the only instance where RDA is reliable. We argue that if satellites can reorient to the lowest energy state then RDA will always give an overestimate of the true MOI. Observations have shown that small nonhydrostatic gravity anomalies exist on Ganymede and Titan, at least at degree 3 and presumably higher. If these anomalies are indicative of the nonhydrostatic anomalies at degree 2 then these imply only a small correction to the MOI, even for Titan, but it is possible that the physical origin of nonhydrostatic degree 2 effects is different from the higher order terms. We conclude that nonhydrostatic effects could be present to an extent that allows Callisto and Titan to be fully differentiated.     

Minimum estimates of the amount and timing of gases released into the martian atmosphere from volcanic eruptions

Volcanism has been a major process during most of the geologic history of Mars. Based on data collected from terrestrial basaltic eruptions, we assume that the volatile content of martian lavas was typically ∼0.5 wt.% water, ∼0.7 wt.% carbon dioxide, ∼0.14 wt.% sulfur dioxide, and contained several other important volatile constituents. From the geologic record of volcanism on Mars we find that during the late Noachian and through the Amazonian volcanic degassing contributed ∼0.8 bar to the martian atmosphere. Because most of the outgassing consisted of greenhouse gases (i.e., CO₂ and SO₂) warmer surface temperatures resulting from volcanic eruptions may have been possible. Our estimates suggest that ∼1.1 × 10²¹ g (∼8 ± 1 m m⁻²) of juvenile water were released by volcanism; slightly more than half the amount contained in the north polar cap and atmosphere. Estimates for released CO₂ (1.6 × 10²¹ g) suggests that a large reservoir of carbon dioxide is adsorbed in the martian regolith or alternatively ∼300 cm cm⁻² of carbonates may have formed, although these materials would not occur readily in the presence of excess SO₂. Up to ∼120 cm cm⁻² (2.2 × 10²⁰ g) of acid rain (H₂SO₄) may have precipitated onto the martian surface as the result of SO₂ degassing. The hydrogen flux resulting from volcanic outgassing may help explain the martian atmospheric D/H ratio. The amount of outgassed nitrogen (∼1.3 mbar) may also be capable of explaining the martian atmospheric ¹⁵N/¹⁴N ratio. Minor gas constituents (HF, HCl, and H₂S) could have formed hydroxyl salts on the surface resulting in the physical weathering of geologic materials. The amount of hydrogen fluoride emitted (1.82 × 10¹⁸ g) could be capable of dissolving a global layer of quartz sand ∼5 mm thick, possibly explaining why this mineral has not been positively identified in spectral observations. The estimates of volcanic outgassing presented here will be useful in understanding how the martian atmosphere evolved over time.     

Impact-induced N2 production from ammonium sulfate: Implications for the origin and evolution of N2 in Titan’s atmosphere

Chemical reactions and volatile supply through hypervelocity impacts may have played a key role for the origin and evolution of both planetary and satellite atmospheres. In this study, we evaluate the role of impact-induced N₂ production from reduced nitrogen-bearing solids proposed to be contained in Titan’s crust, ammonium sulfate ((NH₄)₂SO₄), for the replenishment of N₂ to the atmosphere in Titan’s history. To investigate the conversion of (NH₄)₂SO₄ into N₂ by hypervelocity impacts, we measured gases released from (NH₄)₂SO₄ that was exposed to hypervelocity impacts created by a laser gun. The sensitivity and accuracy of the measurements were enhanced by using an isotope labeling technique for the target. We obtained the efficiency of N₂ production from (NH₄)₂SO₄ as a function of peak shock pressure ranging from ∼8 to ∼45 GPa. Our results indicate that the initial and complete shock pressures for N₂ degassing from (NH₄)₂SO₄ are ∼10 and ∼25 GPa, respectively. These results suggest that cometary impacts on Titan (i.e., impact velocity v_i > ∼8 km/s) produce N₂ efficiently; whereas satellitesimal impacts during the accretion (i.e., v_i < 4 km/s) produce N₂ only inefficiently. Even when using the proposed small amount of (NH₄)₂SO₄ content in the crust (∼4 wt.%) (Fortes, A.D. et al., 2007. Icarus 188, 139-153), the total amount of N₂ provided through cometary impacts over 4.5 Ga reaches ∼2-6 times the present atmospheric N₂ (i.e., ∼7 × 10²⁰-2 × 10²¹ [mol]) based on the measured production efficiency and results of a hydrodynamic simulation of cometary impacts onto Titan. This implies that cometary impacts onto Titan’s crust have the potential to account for a large part of the present N₂ through the atmospheric replenishment after the accretion.     

Mass-radius curve for extrasolar Earth-like planets and ocean planets

By comparison with the Earth-like planets and the large icy satellites of the Solar System, one can model the internal structure of extrasolar planets. The input parameters are the composition of the star (Fe/Si and Mg/Si), the Mg content of the mantle (Mg# = Mg/[Mg + Fe]), the amount of H₂O and the total mass of the planet. Equation of State (EoS) of the different materials that are likely to be present within such planets have been obtained thanks to recent progress in high-pressure experiments. They are used to compute the planetary radius as a function of the total mass. Based on accretion models and data on planetary differentiation, the internal structure is likely to consist of an iron-rich core, a silicate mantle and an outer silicate crust resulting from magma formation in the mantle. The amount of H₂O and the surface temperature control the possibility for these planets to harbor an ocean. In preparation to the interpretation of the forthcoming data from the CNES led CoRoT (Convection Rotation and Transit) mission and from ground-based observations, this paper investigates the relationship between radius and mass. If H₂O is not an important component (less than 0.1%) of the total mass of the planet, then a relation (R/REarth)=ab(M/MEarth) is calculated with (a,b)=(1,0.306) and (a,b)=(1,0.274) for ¹⁰⁻²MEarth<M<MEarth and MEarth<M<10MEarth, respectively. Calculations for a planet that contains 50% H₂O suggest that the radius would be more than 25% larger than that based on the Earth-like model, with (a,b)=(1.258,0.302) for ¹⁰⁻²MEarth<M<MEarth and (a,b)=(1.262,0.275) for MEarth<M<10MEarth, respectively. For a surface temperature of 300 K, the thickness of the ocean varies from 150 to 50 km for planets 1 to 10 times the Earth's mass, respectively. Application of this algorithm to bodies of the Solar System provides not only a good fit to most terrestrial planets and large icy satellites, but also insights for discussing future observations of exoplanets.     

Infrared spectroscopy of crystalline and amorphous diacetylene (C4H2) and implications for Titan's atmospheric composition

New laboratory spectra of crystalline and amorphous diacetylene ice have been recorded in the range of 7000-500 cm⁻¹ (1.4-20 μm) to aid in the identification of solid diacetylene on Saturn's moon Titan. We have established that amorphous diacetylene ice is stable only at temperatures less than 70±1 K. With respect to observations on Titan, the best approach would be to utilize future space-based telescopes to search for the ν₄ (3277/3271 cm⁻¹) in absorption against the reflected light from the sun and the slightly weaker ν₈ absorption bands (676/661 cm⁻¹) in absorption against the continuum emission.     

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.     

High spectral resolution infrared studies of Titan: Winds, temperature, and composition

The Cassini Huygens mission provides a unique opportunity to combine ground-based and spacecraft investigations to increase our understanding of chemical and dynamical processes in Titan’s atmosphere. Spectroscopic measurements from both vantage points enable retrieving global wind structure, temperature structure, and atmospheric composition. An updated analysis of Titan data obtained with the NASA Goddard Space Flight Center’s Infrared Heterodyne Spectrometer (IRHS) and Heterodyne Instrument for Planetary Wind and Composition (HIPWAC) prior to and during the Cassini Huygens mission is compared to retrievals from measurements with the Cassini Composite Infrared Spectrometer (CIRS). IRHS/HIPWAC results include the first direct stratospheric wind measurements on Titan, constraints on stratospheric temperature, and the study of atmospheric molecular composition. These results are compared to CIRS retrievals of wind and temperature profile from thermal mapping data and ethane abundance at 10-15° South latitude, near the equatorial region. IRHS/HIPWAC wind results are combined with other direct techniques, stellar occultation measurements, and CIRS results to explore seasonal variability over nearly one Titan year and to provide an empirical altitude profile of stratospheric winds, varying from ∼50 to 210 m/s prograde. The advantage of fully resolved line spectra in species abundance measurements is illustrated by comparing the possible effect on retrieved ethane abundance by blended spectral features of other molecular constituents, e.g., acetylene (C₂H₂), ethylene (C₂H₄), allene (C₃H₄), and propane (C₃H₈), which overlap the ν₉ band of ethane, and are not resolved at lower spectral resolution. IR heterodyne spectral resolution can discriminate weak spectral features that overlap the ν₉ band of ethane, enabling ethane lines alone to be used to retrieve abundance. Titan’s stratospheric mean ethane mole fraction (8.6±3 ppmv) retrieved from IRHS/HIPWAC emission line profiles (resolving power λ?Δλ∼10⁶) is compared to past values obtained from lower resolution spectra and from CIRS measurements (resolving power λ?Δλ∼2×10³) and more compatible recent analysis. Results illustrate how high spectral resolution ground-based studies complement the spectral and spatial coverage and resolution of moderate spectral resolution space-borne spectrometers.     

Meridional transport of HCN from SL9 impacts on Jupiter

In July 1994, the Shoemaker-Levy 9 (SL9) impacts introduced hydrogen cyanide (HCN) to Jupiter at a well confined latitude band around −44°, over a range of specific longitudes corresponding to each of the 21 fragments (Bézard et al. 1997, Icarus 125, 94-120). This newcomer to Jupiter's stratosphere traces jovian dynamics. HCN rapidly mixed with longitude, so that observations recorded later than several months after impact witnessed primarily the meridional transport of HCN north and south of the impact latitude band. We report spatially resolved spectroscopy of HCN emission 10 months and 6 years following the impacts. We detect a total mass of HCN in Jupiter's stratosphere of 1.5±0.7×10¹³ g in 1995 and 1.7±0.4×10¹³ g in 2000, comparable to that observed several days following the impacts (Bézard et al. 1997, Icarus 125, 94-120). In 1995, 10 months after impact, HCN spread to −30° and −65° latitude (half column masses), consistent with a horizontal eddy diffusion coefficient of K_yy=2-3×10¹⁰ cm² s⁻¹. Six years following impact HCN is observed in the northern hemisphere, while still being concentrated at 44° south latitude. Our meridional distribution of HCN suggests that mixing occurred rapidly north of the equator, with K_yy=2-5×10¹¹ cm² s⁻¹, consistent with the findings of Moreno et al. (2003, Planet. Space Sci. 51, 591-611) and Lellouch et al. (2002, Icarus 159, 112-131). These inferred eddy diffusion coefficients for Jupiter's stratosphere at 0.1-0.5 mbar generally exceed those that characterize transport on Earth. The low abundance of HCN detected at high latitudes suggests that, like on Earth, polar regions are dynamically isolated from lower latitudes.     

Saturn's magnetospheric interaction with Titan as defined by Cassini encounters T9 and T18: New results

We present new results of Cassini's T9 flyby with complementary observations from T18. Based on Cassini plasma spectrometer (CAPS) and Cassini magnetometer (MAG), compositional evidence shows the upstream flow for both T9 and T18 appears composed of light ions (H⁺ and H₂⁺), with external pressures ∼30 times lower than that for the earlier TA flyby where heavy ions dominated the magnetospheric plasma. When describing the plasma heating and sputtering of Titan's atmosphere, T9 and T18 can be considered interactions of low magnetospheric energy input. On the other hand, T5, when heavy ion fluxes are observed to be higher than typical (i.e., TA), represents the limiting case of high magnetospheric energy input to Titan's upper atmosphere. Anisotropy estimates of the upstream flow are 1<T_⊥/T_‖<3 and the flow is perpendicular to B, indicative of local picked up ions from Titan's H and H₂ coronae extending to Titan's Hill sphere radius. Beyond this distance the corona forms a neutral torus that surrounds Saturn. The T9 flyby unexpectedly resulted in observation of two “wake” crossings referred to as Events 1 and 2. Event 2 was evidently caused by draped magnetosphere field lines, which are scavenging pickup ions from Titan's induced magnetopause boundary with outward flux ∼2×10⁶ ions/cm²/s. The composition of this out flow is dominated by H₂⁺ and H⁺ ions. Ionospheric flow away from Titan with ion flux ∼7×10⁶ ion/cm²/s is observed for Event 1. In between Events 1 and 2 are high energy field aligned flows of magnetosphere protons that may have been accelerated by the convective electric field across Titan's topside ionosphere. T18 observations are much closer to Titan than T9, allowing one to probe this type of interaction down to altitudes ∼950 km. Comparisons with previously reported hybrid simulations are made.