Dissociative recombination of nitrile ions with implications for Titan's upper atmosphere

Nitrile ions are abundant in Titan's upper atmosphere and are expected to be lost mainly via dissociative recombination with free electrons. We review in this paper a series of experimental results on the dissociative recombination reactions of nitrile ions known/expected to be present in Titan's upper atmosphere. The experiments were all performed at the heavy ion storage ring CRYRING in Stockholm, Sweden, and the results presented here include information on rate coefficients at electron temperatures relevant for Titan's upper atmosphere as well as information on the product branching fractions of the reactions. We discuss implications of the results for Titan's atmosphere. As an example the presented results support a statement by Krasnopolsky (2009) that nitriles do not degrade to yield N₂ again in Titan's atmosphere, indicating that condensation and polymerization with precipitation to the surface are their ultimate fate.     

Photochemical kinetics uncertainties in modeling Titan's atmosphere: First consequences

Uncertainties carried by the different kinetic parameters included in photochemical models of planetary atmospheres have rarely been considered even if they are supposed to be contributing mostly to the inconsistencies between observations and computed predictions. In this paper, we report the first detailed analysis of the propagation of uncertainties carried by the reaction rate coefficients included in an up-to-date photochemical model of Titan's atmosphere. Monte Carlo calculations performed on these reaction rate coefficients have been used to introduce their uncertainties and to investigate their significance on the photochemical modeling of Titan's atmosphere. Crude approximations in the implemented physical processes have been adopted to limit the number of free parameters. This allows us to pinpoint specifically the importance of chemical processes uncertainties in Titan's photochemical models and to evaluate their chemical robustness. First implications of this preliminary study related to purely chemical rate coefficient uncertainties are discussed. They are important enough to question indeed any comparisons between theoretical models with observations as well as any potential conclusions subsequently inferred. Since the latest missions, such as Cassini–Huygens, are likely to induce an ever-increasing interest for such kind of comparing studies, our conclusions show that it is crucial to reform the way we think of, and use, current photochemical models to understand the processes occurring in the atmospheres of the outer Solar System.     

Radiocarbon on Titan

Abstract— We explore the likely production and fate of ¹⁴C in the thick nitrogen atmosphere of Saturn's moon Titan and investigate the constraints that measurements of ¹⁴C might place on Titan's photochemical, atmospheric transport and surface‐atmosphere interaction processes. Titan's atmosphere is thick enough that cosmic‐ray flux limits the production of ¹⁴C: absence of a strong magnetic field and the increased distance from the Sun suggest production rates of ?9 atom/cm²/s, ?4x higher than Earth. The fate and detectability of ¹⁴C depends on the chemical species into which it is incorporated: as methane it would be hopelessly diluted even in only the atmosphere. However, in the more likely case that the ¹⁴C attaches to the haze that rains out onto the surface (as tholin, HCN or acetylene and their polymers), haze in the atmosphere or recently deposited on the surface would be quite radioactive. Such radioactivity may lead to a significant enhancement in the electrical conductivity of the atmosphere which will be measured by the Huygens probe. Measurements with simple detectors on future missions could place useful constraints on the mass deposition rates of photochemical material on the surface and identify locations where surface deposits of such material are “freshest”.     

Carbon isotopic enrichment in Titan's tholins? Implications for Titan's aerosols

Since the discovery of the main composition of Titan's atmosphere, many laboratory experiments have been carried out to reproduce its chemical evolution, particularly the formation of organic haze particles found throughout this atmosphere. Some of these simulations have produced solid products—referred to as Titan's tholins—that are assumed to have properties similar to those of Titan's aerosols. In the present work, we focus on the possible isotopic fractionation of carbon during the processes involved in the formation of Titan's tholins. Initial ¹²C/¹³C isotopic ratios measured on tholins made in the laboratory using cold plasma discharges are presented. Measurements of isotopic enhancement in ¹³C (δ¹³C), both on tholins and on the initial gas mixture (N₂:CH₄ (98:2)) used to produce them do not show any clear deficit or enrichment in ¹³C relative to ¹²C in the lab-made tholins compared to the initial gas mixture. Preliminary data recovered from the Aerosol Collector Pyrolyzer (ACP) experiment of the Huygens probe suggests that Titan's aerosols may also be exempt of carbon isotopic enrichment. This observation creates possibilities for deeper analysis of ACP experiment data.     

The history of titan,of Saturn's rings and magnetic field,and the nature of short-period comets

Titan in the past, just as Ganymede, had a massive ice envelope subjected to volumetric electrolysis under the action of unipolar electric current generated through the interaction of the satellite with Saturn's magnetosphere. The electrolysis products concentration required to cause explosion could become accumulated only under conditions of an exponential decay of Saturn's magnetic field in time (with τ_? ≈ 0.55 Gyr) which implies a relict nature of the field and agrees with the present ideas on the planet's structure. The explosion of the electrolysis products contained in the ice envelope resulted in Titan's having lost ～13% of its mass in the form of gas (mainly of water vapor) and solid ice fragments, as well as in the appearance on Titan of an atmosphere (of volatile products from incomplete combustion of hydrogen and hydrocarbons) and a deep (～1000 km) ocean of liquid water. The presence of liquid water on Titan's surface is confirmed by an analysis of the available microwave measurements of brightness temperature. The condensation of the water vapor lost by Titan produced the visible inner rings of Saturn while large solid fragments of the ice envelope govern their dynamics. These are also located in the gap between Rhea and Titan (the G ring?). Most of the ice fragments were swept out from Saturn's system through perturbations by Titan. They made up a reservoir of cometary nuclei beyond Jupiter's orbit. Arguments are presented in favor of a recent (3–10 thousand years ago) explosion of Titan. Some implications from these concepts, lending themselves to observational testing, are pointed out.     

Titan's damp ground: Constraints on Titan surface thermal properties from the temperature evolution of the Huygens GCMS inlet

Abstract— A simple thermal model is developed to determine the temperature history of the inlet tube of the Huygens probe gas chromatograph mass spectrometer (GCMS) after its fortuitous emplacement on the surface of Saturn's moon Titan. The model parameters are adjusted to match the recorded temperature history of a nearby heater, taking into account heat losses by conduction to the rest of the probe and to Titan's cold atmosphere. The model suggests that after impact when forced convective cooling ceased, the inlet temperature rose from ?110 K to an asymptotic value of only ?145 K. This requires that the inlet was embedded in a surface that acted as an effective heat sink, most plausibly interpreted as wet or damp with liquid methane. The data appear inconsistent with a tar or dry, fine‐grained surface, and the inlet was not warm enough to devolatilize methane hydrate.     

The distribution of Titan's high-altitude (out to ∼50,000 km) exosphere from energetic neutral atom (ENA) measurements by Cassini/INCA

We report observations of Titan's high-altitude exosphere detected out to about 50,000 km altitude. The observations were made by the Ion Neutral Camera (INCA) on board the Cassini spacecraft. INCA detects energetic neutral atoms (ENA) that are formed when the ambient magnetospheric ions charge exchange with Titan's neutral atmosphere and exosphere. We find that Titan's exospheric H₂ distribution follows closely a full Chamberlain distribution including ballistic, escaping and satellite distributions. As expected, neutral densities are dominated by a satellite distribution above about 10,000 km. The maximum detectable extent of the exosphere (～50,000 km) coincides with the radius of the Hill sphere of gravitational influence from Saturn. While we find no direct indications of a neutral Titan torus with densities greater than about 1000 cm^?3, we observe interesting asymmetries in the distribution that warrants further investigation. Based on these findings we compute the average precipitating ENA flux to be about 5×10⁶ keV/(cm² s), or 8×10^?3 erg/(cm² s), which is directly comparable to that of precipitating energetic ions (Sittler, et al., 2009) and slightly higher than that of solar EUV (Tobiska, 2004). Thus, the energy deposited by precipitating ENAs must also be taken into consideration when studying the energy balance of Titan's thermosphere.     

A numerical assessment of simple airblast models of impact airbursts

Asteroids and comets 10–100 m in size that collide with Earth disrupt dramatically in the atmosphere with an explosive transfer of energy, caused by extreme air drag. Such airbursts produce a strong blastwave that radiates from the meteoroid's trajectory and can cause damage on the surface. An established technique for predicting airburst blastwave damage is to treat the airburst as a static source of energy and to extrapolate empirical results of nuclear explosion tests using an energy‐based scaling approach. Here we compare this approach to two more complex models using the iSALE shock physics code. We consider a moving‐source airburst model where the meteoroid's energy is partitioned as two‐thirds internal energy and one‐third kinetic energy at the burst altitude, and a model in which energy is deposited into the atmosphere along the meteoroid's trajectory based on the pancake model of meteoroid disruption. To justify use of the pancake model, we show that it provides a good fit to the inferred energy release of the 2013 Chelyabinsk fireball. Predicted overpressures from all three models are broadly consistent at radial distances from ground zero that exceed three times the burst height. At smaller radial distances, the moving‐source model predicts overpressures two times greater than the static‐source model, whereas the cylindrical line‐source model based on the pancake model predicts overpressures two times lower than the static‐source model. Given other uncertainties associated with airblast damage predictions, the static‐source approach provides an adequate approximation of the azimuthally averaged airblast for probabilistic hazard assessment.     

Titan atmosphere profiles from Huygens engineering (temperature and acceleration) sensors

Data from several Huygens probe housekeeping sensors (an engineering accelerometer and housekeeping temperature sensors) are studied to determine how effectively such nonideal instruments may characterize the density or temperature structure of the atmosphere. While only confirming the results of the dedicated atmospheric structure instrument, this exercise is of relevance to possible future missions to various bodies which might not be equipped with such science-grade sensors able to accurately profile the atmosphere top-to-bottom. It is found that for typical engineering accelerometers with 8-bit resolution, the atmosphere density for ∼4 scale heights above the peak deceleration altitude may be recovered. If, as with Huygens, the peak deceleration exceeds the range of the accelerometers, recovery of an additional scale height or so below the peak is still possible, but relies on accurate total velocity knowledge. Engineering temperature sensors can, with care, be analyzed to recover the temperature structure of at least the lowest ∼30 km of Titan's atmosphere. Fortunately, in the case of Huygens, data from the surface after landing were available to constrain models of heat leaks which offset the observed temperature from that of the ambient air during descent; data from before and during the entry phase on other missions would be similarly useful. When corrections are made for the estimated heat transfer processes, the atmospheric temperature can be recovered to within about 3 K.     

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

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

Near-surface winds at the Huygens site on Titan: Interpretation by means of a general circulation model

This study aims at interpreting the zonal and meridional wind in Titan's troposphere measured by the Huygens probe by means of a general circulation model. The numerical simulation elucidates the relative importance of the seasonal variation in the Hadley circulation and Saturn's gravitational tide in affecting the actual wind profile. The observed reversal of the zonal wind at two altitudes in the lower troposphere can be reproduced with this model only if the near-surface temperature profile is asymmetric about the equator and substantial seasonal redistribution of angular momentum by the variable Hadley circulation takes place. The meridional wind near the surface is mainly caused by the meridional pressure gradient and is thus a manifestation of the Hadley circulation. Southward meridional wind in the PBL (planetary boundary layer) is consistent with the near-surface temperature at the equator being lower than at mid southern latitudes. Even small changes in the radiative heating profile in the troposphere can substantially affect the mean zonal and meridional wind including their direction. Saturn's gravitational tide is rather weak at the Huygens site due to the proximity to the equator, and does not clearly manifest itself in the instantaneous vertical profile of wind. Nevertheless, the simulated descent trajectory is more consistent with the observation if the tide is present. Because of a different force balance in Titan's atmosphere from terrestrial conditions, PBL-specific wind systems like on Earth are unlikely to exist on Titan.     

“Comet‐tail” ejecta streaks: A predicted cratering landform unique to Titan

Abstract— A model for an impact ejecta landform peculiar to Saturn's moon Titan is presented. Expansion of the ejecta plume from moderate‐sized craters is constrained by Titan's thick atmosphere. Much of the plume is collimated along the incoming bolide's trajectory, as was observed for plumes from impacts on Jupiter of P/Shoemaker‐Levy‐9, but is retained as a linear, diagonal ejecta cloud, unlike on Venus where the plume “blows out.” On Titan, the blowout is suppressed because the vertically‐extended atmosphere requires a long wake to reach the vacuum of space, and the modest impact velocities mean plume expansion along the wake is slow enough to allow the wake to close off. Beyond the immediate ejecta blanket around the crater, distal ejecta is released into the atmosphere from an oblique line source: this material is winnowed by the zonal wind field to form streaks, with coarse radar‐bright particles transported less far than fine radar‐dark material. Thus, the ejecta form two distinct streaks faintly reminiscent of dual comet tails, a sharply W‐E radar‐dark one, and a less swept and sometimes comma‐shaped radar‐bright one.     

Huygens Probe descent dynamics inferred from Channel B signal level measurements

The signal strength of the Huygens Probe Channel B transmission to the Cassini Orbiter was monitored during the Probe descent through Titan's atmosphere on 14 January 2005. A model of the Probe motion during the mission was constructed to include Probe spin, coning motion and tilt caused by varying wind speeds. This simple model is sufficient to reproduce the most prominent features seen in the signal level measurements. It provides estimates of the coning and tilt angles as well as the direction of the Huygens coordinate axes over extended time intervals in the mission.     

The Huygens Probe Descent Trajectory Working Group: Organizational framework,goals, and implementation

Cassini/Huygens, a flagship mission to explore the rings, atmosphere, magnetic field, and moons that make up the Saturn system, is a joint endeavor of the National Aeronautics and Space Administration, the European Space Agency, and Agenzia Spaziale Italiana. Comprising two spacecraft—a Saturn orbiter built by NASA and a Titan entry/descent probe built by the European Space Agency—Cassini/Huygens was launched in October 1997. The Huygens probe parachuted to the surface of Titan in January 2005. During the descent, six science instruments provided in situ measurements of Titan's atmosphere, clouds, and winds, and photographed Titan's surface. To correctly interpret and correlate results from the probe science experiments, and to provide a reference set of data for ground-truth calibration of orbiter remote sensing measurements, an accurate reconstruction of the probe entry and descent trajectory and surface landing location is necessary. The Huygens Descent Trajectory Working Group was chartered in 1996 as a subgroup of the Huygens Science Working Team to develop and implement an organizational framework and retrieval methodologies for the probe descent trajectory reconstruction from the entry altitude of 1270 km to the surface using navigation data, and engineering and science data acquired by the instruments on the Huygens Probe. This paper presents an overview of the Descent Trajectory Working Group, including the history, rationale, goals and objectives, organizational framework, rules and procedures, and implementation.     

A method to determine the atmospheric temperature profile from in situ pressure data: Application to Titan

We have developed a new method to analyse in situ observations of atmospheric variables of state: the reconstruction of the vertical temperature profile from pressure measurements accompanied by rough knowledge of the atmospheric composition and the aerodynamical response properties of the descent vehicle. We can use the method to construct the temperature profile when no direct measurements are available, as well as to analyse the consistency between data from different instruments. We applied the method to the Huygens measurements of Titan's atmosphere, determining the aerodynamical drag properties from radar altimeter data. We discovered that the temperature profile computed in this manner differs from the profile from the temperature sensor (TEM) of the probe by up to 5% in the altitude range of 0–60 km, and up to 10% at higher altitudes due to increased noise. The method gives a tropopause altitude of about 50 km and a surface temperature of about 98 K, in contrast to the TEM temperature measurements. Our error analysis shows that these differences are caused by the known discrepancy in the Huygens altimeter data, with the estimates made by the reconstruction algorithm contributing only 1–2% of error.     

Estimated impact shock production of N2 and organic compounds on early Titan

Bombardment of Titan by Uranus-Neptune planetesimals and/or fragments of a disrupted Hyperion progenitor supplied more than enough energy to drive vigorous atmospheric shock chemistry. Chemical equilibrium modeling of the shock products in simulated atmospheres indicates that impact energy has produced large amounts of N₂ and organic compounds over Titan's history. The mole fraction of organic compounds in the shocked gas mixture (T = 1200−2500°K, P = 10⁻¹−10³bar) reaches a maximum of approximately 3% in a current Titan mixture and 12% in a primordial CH₄, NH₃-rich mixture. Atmospheric water mixing ratio controls the organic yield in shock reactions, but its limiting effect may have been reduced by cold-trapping of water in a cooling atmosphere. Kinetic inhibition of graphite formation in the shocked gas enhanced the yield of radicals and organic. The resulting mixture of carbonaceous soot and condensed hydrocarbons subsequently settled onto the surface; the depth of the generated layer was on the order of hundreds of meters. Impact shock energy was capable of converting massive amounts of NH₃ to N₂ early in Titan history—over twice the present atmospheric and 1.5 times the total ocean-atmospheric inventory of N₂. Shock conversion of NH₃ into N₂ bypasses the difficulties of other schemes of N₂ production and may have been of singular importance in Titan's atmospheric evolution.     

Atmospheric shock waves after impacts of cosmic bodies up to 1000 m in diameter

The results of numerical modeling of meteoroids' interaction with Earth's atmosphere are presented. We model the entry in two dimensions and then interpolate the results into a 3‐D model to calculate interaction of shock waves with the surface. Maximum shock pressures, wind speeds, and areas subjected to substantial overpressure are calculated for oblique impacts of asteroids and comets. We show that vertical impacts produce a smaller damage zone on the surface than oblique ones. Damage caused by shock waves covers an order of magnitude larger area than any other hazardous effects. The function of energy release in the atmosphere, which is traditionally used in meteoritics, has a limited application if cosmic bodies are larger than tens of meters in diameter: at each time moment energy is smoothed along a substantial length of the trajectory; both emitted radiation (routinely used for calibration of semi‐analytical models) and shock wave amplitude are complex functions of temperature–density distributions in atmosphere.     

On the influence of nonlinearities on the eigenfrequencies of five-minute oscillations of the Sun

Fitting the results of linear normal-mode analysis of the solar five-minute oscillations to the observed k - ω diagram selects a class of models of the Sun's envelope. It is a property of all the models in this class that their convection zones are too deep to permit substantial transmission of internal g modes of degree 20 or more. This is in apparent conflict with Hill and Caudell's (1979) claim to have detected such modes in the photosphere. A proposal to resolve the conflict was made by Rosenwald and Hill (1980). They pointed out that despite the impressive agreement between linearized theory and observation, nonlinear phenomena in the solar atmosphere might influence the eigenfrequencies considerably. In particular, they suggested that a correct nonlinear analysis could predict a shallow convection zone. This paper is an enquiry into whether their hypothesis is plausible. We construct k - ω diagrams assuming that the modes suffer local nonlinear distortions in the atmosphere that are insensitive to the amplitude of oscillation over the range of amplitudes that are observed. The effect of the nonlinearities on the eigenfrequencies is parameterized in a simple way. Taking a class of simple analytical models of the Sun's envelope, we compute the linear eigenfrequencies of one model and show that no other model can be found whose nonlinear eigenfrequencies agree with them. We show also that the nonlinear eigenfrequencies of a particular solar model with a shallow convective zone, computed with more realistic physics, cannot be made to agree with observation. We conclude, therefore, that the hypothesis of Rosenwald and Hill is unlikely to be correct.     

The rate of small impacts on Earth

Abstract— Asteroids tens to hundreds of meters in diameter constitute the most immediate impact hazard to human populations, yet the rate at which they arrive at Earth's surface is poorly known. Astronomic observations are still incomplete in this size range; impactors are subjected to disruption in Earth's atmosphere, and unlike the Moon, small craters on Earth are rapidly eroded. In this paper, we first model the atmospheric behavior of iron and stony bodies over the mass range 1–10¹² kg (size range 6 cm‐1 km) taking into account deceleration, ablation, and fragmentation. Previous models in meteoritics deal with rather small masses (<10⁵–10⁶ kg) with the aim of interpreting registered fireballs in atmosphere, or with substantially larger objects without taking into account asteroid disruption to model cratering processes. A few earlier attempts to model terrestrial crater strewn fields did not take into account possible cascade fragmentation. We have performed large numbers of simulations in a wide mass range, using both the earlier “pancake” models and also the separated fragments model to develop a statistical picture of atmosphere‐bolide interaction for both iron and stony impactors with initial diameters up to ?1 km. Second, using a compilation of data for the flux at the upper atmosphere, we have derived a cumulative size‐frequency distribution (SFD) for upper atmosphere impactors. This curve is a close fit to virtually all of the upper atmosphere data over 16 orders of magnitude. Third, we have applied our model results to scale the upper atmosphere curve to a flux at the Earth's surface, elucidating the impact rate of objects <1 km diameter on Earth. We find that iron meteorites >5 times 10⁴ kg (2.5 m) arrive at the Earth's surface approximately once every 50 years. Iron bodies a few meters in diameter (10⁵–10⁶ kg), which form craters ?100 m in diameter, will strike the Earth's land area every 500 years. Larger bodies will form craters 0.5 km in diameter every 20,000 years, and craters 1 km in diameter will be formed on the Earth's land area every 50,000 years. Tunguska events (low‐level atmospheric disruption of stony bolides >10⁸ kg) may occur every 500 years. Bodies capable of producing hazardous tsunami (?200 m diameter projectiles) should strike the Earth's surface every ?100,000 years. This data also allows us to assess the completeness of the terrestrial crater record for a given area over a given time interval.     

The reflectivity spectrum and opposition effect of Titan's surface observed by Huygens' DISR spectrometers

We determined Titan's reflectivity spectrum near the Huygens' landing site from observations taken with the Descent Imager/Spectral Radiometer below 500 m altitude, in particular the downward-looking photometer and spectrometers. We distinguish signal coming from illumination by sunlight and the lamp onboard Huygens based on their different spectral signatures. For the sunlight data before landing, we find that spatial variations of Titan's reflectivity were only ～0.8%, aside from the phase angle dependence, indicating that the probed area within ～100 m of the landing site was very homogeneous. Only the very last spectrum taken before landing gave a 3% brighter reflectivity, which probably was caused by one bright cobble inside its footprint. The contrast of the cobble was higher at 900 nm wavelength than at 600 nm.For the data from lamp illumination, we confirm that the phase function of Titan's surface displays a strong opposition effect as found by Schröder and Keller (2009. Planetary and Space Science 57, 1963–1974). We extend the phase function to even smaller phase angles (0.02°), which are among the smallest phase angles observed in the solar system. We also confirm the reflectivity spectrum of the dark terrain near the Huygens' landing site between 900 and 1600 nm wavelength by Schröder and Keller (2008. Planetary and Space Science 56, 753–769), but extend the spectrum down to 435 nm wavelength. The reflectivity at zero phase angle peaks at 0.45±0.06 around 750 nm wavelength and drops down to roughly 0.2 at both spectral ends. Our reflectivity of 0.45 is much higher than all previously reported values because our observations probe lower phase angles than others. The spectrum is very smooth except for a known absorption feature longward of 1350 nm. We did not detect any significant variation of the spectral shape along the slit for exposures after landing, probing a 25×4 cm² area. However, the recorded spectral shape was slightly different for exposures before and after landing. This difference is similar to the spectral differences seen on scales of kilometers (Keller et al., 2008. Planetary and Space Science 56, 728–752), indicating that most observations may probe spatially variable contributions from two basic materials, such as a dark soil partially covered by bright cobbles.We used the methane absorption features to constrain the methane mixing ratio near the surface to 5.0±0.3%, in agreement with the 4.92±0.24% value measured in situ by Niemann et al. (2005. Nature 438, 779–784), but smaller than their revised value of 5.65±0.18% (Niemann et al., 2010. Journal of Geophysical Research 115, E12006). Our results were made possible by an in depth review of the calibration of the spectroscopic and photometric data.