Evaporation effects on the formation of martian gullies

In order to investigate the formation of martian gullies and the stability of fluids on Mars, we examined about 120 gully images. Twelve HiRISE images contained a sufficient number of Transverse Aeolian Ridges (TARs) associated with the gullies to make the following measurements: overall gully length, length of the alcove, channel and apron, and we also measured the frequency of nearby TARs. Six of the 12 images examined showed a statistically significant negative correlation between overall gully length (alcove, channel and apron length) and TAR frequency. Previous experimental work from our group has shown that at temperatures below ∼200 K, evaporation rate increases by about an order of magnitude as wind speed increases from 0 to ∼15 m/s. Thus the negative correlations we observe between gully length and dune frequency can be explained by formation at temperatures below ∼200 K where wind speed/evaporation is a factor governing gully length. In these cases evaporation of the fluid carving the gully was a constraint on their dimensions. Cases where there is no correlation between gully length and TAR frequency, can be explained by formation at temperatures >200 K. The temperatures are consistent with Global Circulation Model and Thermal Emission Spectrometer (TES) data for these latitudes. The temperatures suggested by these trends are consistent with the fluid responsible for gully formation being a strong brine, such as Fe₂(SO₄)₃ which has a eutectic temperature of ∼200 K. We also find that formation timescales for gullies are 10⁵-10⁶ years.     

Structure of the Earth’s circumsolar dust ring

Interplanetary dust particles from comets and asteroids pervade the Solar System and become temporarily trapped into orbital resonances with Earth, leading to a circumsolar dust ring. Using the unique vantage point of the Spitzer Space Telescope from its Earth-trailing solar orbit, we have measured for the first time the azimuthal structure of the Earth’s resonant dust ring. There is a relative paucity of particles within 0.1 AU of the Earth, followed by an enhancement in a cloud that is centered 0.2 AU behind Earth with a width of 0.08 AU along the Earth’s orbit. The North ecliptic pole is ∼3% brighter at 8 μm wavelength when viewed from inside the enhancement. The presence of azimuthal asymmetries in debris disks around other stars is considered strong evidence for planets. By measuring the properties of the Earth’s resonant ring, we can provide “ground truth” to models for interactions of planets and debris disks, possibly leading to improved predictions for detectability of life-bearing planets. The low amplitude of the azimuthal asymmetry in the Earth’s circumsolar ring suggests significant contributions to the zodiacal light from particles that are large (>30 μm) or have large orbital eccentricity that makes capture into mean motion resonances inefficient.     

A new analysis of the ESO Very Large Telescope (VLT) observations of Titan at 2 μm

We have performed an analysis of ESO Very Large Telescope (VLT) observations of Titan at 2 μm. The data were acquired with the Nasmyth Adaptative Optics System Near-Infrared Imager and Spectrograph (NAOS/CONICA), on the 16th of January 2005, that is 2 days after the landing of the Huygens probe (Hirtzig et al., 2007). The data consist in 21 spectra taken along two diameters of Titan’s disk at wavelengths between 2.03 and 2.5 μm. This range covers a part of the 2 μm methane window and the adjacent band. The data received a preliminary analysis in a recent paper (Negrão et al., 2007), essentially focused on the surface albedo near Huygens landing site. In this work, we perform an in-depth analysis to retrieve information about several aspects: the latitude haze distribution in the stratosphere and in the low atmosphere, the latitudinal variation of the surface albedo and its spectral behaviour. Also, this analysis allowed us to make sensitivity tests on the influence of the scatterer profiles on the retrieved surface albedo and its spectral slope. The news analysis confirms that, as was the case with VIMS observations at the same epoch, the Northern (currently winter) Hemisphere contains more haze than the southern one (Summer Hemisphere). The sensitivity tests show that the scatterer profiles have just a little impact on the surface albedo and its spectral slope. The analysis seems to confirm the presence of H₂O and CH₄ ices.     

Determination of asteroid masses from their close encounters with Mars

An obstacle to the asteroid mass determination lies in the difficulty in isolating the gravitational perturbation exerted by a single asteroid on the planets, being strongly correlated and mixed up with those of many other asteroids. This hindrance may be avoided by the method of analysis presented here: an asteroid mass is estimated in correspondence with its close encounters with Mars where the acceleration it induces on the planet can be sufficiently disentangled from those generated by the remaining asteroid masses to calculate. We test this technique in the analysis of range observations to Mars Global Surveyor and Mars Express performed from 1999 to 2007. For this purpose, we adopt the dynamical model of the planetary ephemeris INPOP06 (Fienga et al., 2008), which includes the gravitational influences of the 300 most perturbing asteroids of the Martian orbit. We obtain the solutions of 10 asteroid masses that have the largest effects on this orbit over the period examined: they are generally in good agreement with determinations recently published.     

The light curve of asteroid 2867 Steins measured by VIRTIS-M during the Rosetta fly-by

On 5 September 2008, the Rosetta spacecraft encountered the asteroid 2867 Steins on its way to the comet 67P/Churyumov-Gerasimenko. This was the first of two planned asteroid fly-bys performed by the probe, the second fly-by being with the much larger asteroid 21 Lutetia in July 2010. The VIRTIS imaging spectrometer (IFOV 0.250 mrad, overall spectral range 0.25-5.1 μm) onboard Rosetta acquired data of Steins already before the closest approach phase, when the target was spatially unresolved, in order to obtain a light curve of the asteroid in the infrared spectral range extending up to 5 μm, that was never explored before. The VIRTIS light curve campaign started at 11:30 UTC onboard time, when the spacecraft was about 221,377 km away from the target, and ended at 17:58 UTC, at a distance of 20,741 km away from Steins. During this timeframe, the solar phase angle of the asteroid was roughly constant, ranging from 38.2° to 36.3°.Assuming the most recent value derived for the rotational period of Steins (Lamy et al., 2008), the VIRTIS observations covered slightly more than one rotation of the asteroid. In this interval, VIRTIS collected 8 hyperspectral cubes where Steins was captured 119 times, both in the visual and in the infrared range. Given the low signal and the unresolved appearance of the source, for which the instrument was not designed, only a small subset of wavelengths turned out to be suitable to sample the light curve. Nevertheless, in both the VIS and NIR ranges we find a similar trend, with two different maxima and minima during one rotational period, and amplitudes consistent with the results in the visual range obtained in previous works, including the data set acquired by the OSIRIS camera onboard Rosetta. We also report the presence of a new broad feature centered at approximately 0.81-0.82 μm, which is seen in the visual data throughout the rotation of the asteroid.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

The opposition and tilt effects of Saturn’s rings from HST observations

The two major factors contributing to the opposition brightening of Saturn’s rings are (i) the intrinsic brightening of particles due to coherent backscattering and/or shadow hiding on their surfaces, and (ii) the reduced interparticle shadowing when the solar phase angle α → 0°. We utilize the extensive set of Hubble Space Telescope observations (Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 199–223) for different elevation angles B and wavelengths λ to disentangle these contributions. We assume that the intrinsic contribution is independent of B, so that any B dependence of the phase curves is due to interparticle shadowing, which must also act similarly for all λ’s. Our study complements that of Poulet et al. (Poulet, F., Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 224), who used a subset of data for a single B ∼ 10°, and the French et al. (French, R.G., Verbiscer, A., Salo, H., McGhee, C.A., Dones, L. [2007b] PASP 119, 623–642) study for the B ∼ 23° data set that included exact opposition. We construct a grid of dynamical/photometric simulation models, with the method of Salo and Karjalainen (Salo and Karjalainen [2003]. Icarus 164, 428–460), and use these simulations to fit the elevation-dependent part of opposition brightening. Eliminating the modeled interparticle component yields the intrinsic contribution to the opposition effect: for the B and A rings it is almost entirely due to coherent backscattering; for the C ring, an intraparticle shadow hiding contribution may also be present.Based on our simulations, the width of the interparticle shadowing effect is roughly proportional to B. This follows from the observation that as B decreases, the scattering is primarily from the rarefied low filling factor upper ring layers, whereas at larger B’s the dense inner parts are visible. Vertical segregation of particle sizes further enhances this effect. The elevation angle dependence of interparticle shadowing also explains most of the B ring tilt effect (the increase of brightness with elevation). From comparison of the magnitude of the tilt effect at different filters, we show that multiple scattering can account for at most a 10% brightness increase as B → 26°, whereas the remaining 20% brightening is due to a variable degree of interparticle shadowing. The negative tilt effect of the middle A ring is well explained by the the same self-gravity wake models that account for the observed A ring azimuthal brightness asymmetry (Salo, H., Karjalainen, R., French, R.G. [2004]. Icarus 170, 70–90; French, R.G., Salo, H., McGhee, C.A., Dones, L. [2007]. Icarus 189, 493–522).     

The effects and characteristics of atmospheric dust during martian global dust storm 2001A

We present retrieved trends in dust optical depth, dust effective radius and surface temperature from our analysis of Mars Global Surveyor Thermal Emission Spectrometer daytime data from global dust storm 2001A, and describe their significance for the martian dust cycle. The dust optical depth becomes correlated with surface pressure during southern spring and summer in years both with and without a global dust storm, indicating that global dust mixing processes are important at those seasons. The correlation is low at other times of the year. We found that the observed decay of optical depths at the later stages of the dust storm match, to first-order, theoretical values of clearing from Stokes–Cunningham fallout of the dust. Zonally averaged effective radius is constant within standard deviation of results (between 1.2 and 2.0 μm, with a global mean for all seasons of 1.7 μm), at all latitudes and seasons except at southern latitudes of 35° and higher around equinoxes in both martian years, where it is larger than average (2–3 μm). The emergence and disappearance of these larger particles correlates with observations of polar cap edge storms at those latitudes. Northern latitude observations under similar conditions did not yield a similar trend of larger average effective radii during the equinoxes. We also report on a linear correlation between daytime surface temperature drop and rise in optical depth during the global dust storm. Global dust storm 2001A produced a significant optical depth and surface temperature change.     

Retrieval algorithm for atmospheric dust properties from Mars Global Surveyor Thermal Emission Spectrometer data during global dust storm 2001A

We present a new parameter retrieval algorithm for Mars Global Surveyor Thermal Emission Spectrometer data. The algorithm uses Newtonian first-order sensitivity functions of the infrared spectrum in response to variations in physical parameters to fit a model spectrum to the data at 499, 1099, and 1301 cm⁻¹. The algorithm iteratively fits the model spectrum to data to simultaneously retrieve dust extinction optical depth, effective radius, and surface temperature. There are several sources of uncertainty in the results. The assumed dust vertical distribution can introduce errors in retrieved optical depth of a few tens of percent. The assumed dust optical constants can introduce errors in both optical depth and effective radius, although the systematic nature of these errors will not affect retrieval of trends in these parameters. The algorithm does not include the spectral signature of water ice, and hence data needs to be filtered against this parameter before the algorithm is applied. The algorithm also needs sufficient dust spectral signature, and hence surface-to-atmosphere temperature contrast, to successfully retrieve the parameters. After the application of data filters the algorithm is both relatively accurate and very fast, successfully retrieving parameters, as well as meaningful parameter variability and trends from tens of thousands of individual spectra on a global scale (Elteto, A., Toon, O.B. [2010]. Icarus, this issue). Our results for optical depth compare well with TES archive values when corrected by the single scattering albedo. Our results are on average 1–4 K higher in surface temperatures from the TES archive values, with greater differences at higher optical depths. Our retrieval of dust effective radii compare well with the retrievals of Wolff and Clancy (Wolff, M.J., Clancy, R.T. [2003]. J. Geophys. Res. 108 (E9), 5097) for the corresponding data selections from the same orbits.     

Formation of an eroded lava channel within an Elysium Planitia impact crater: Distinguishing between a mechanical and thermal origin

A lava channel identified on the wall of an Elysium Planitia impact crater is investigated to identify the dominant erosion mechanism, mechanical vs. thermal, acting during channel formation. Observations of channel morphology are used to supplement analytical models of lava channel formation in order to calculate the duration of channel formation, the velocity of the lava flowing through the channel, and the erosion rate in each erosion regime considered. Results demonstrate that the channel observed in the Elysium Planitia impact crater formed primarily due to mechanical erosion. In a more general sense, results of this study suggest that lava channels can form primarily due to thermal erosion in the presence of more gradual slopes and more consolidated substrates whereas lava channels can form primarily due to mechanical erosion in the presence of more energetic flows on steeper slopes and more poorly consolidated substrates. Therefore, both erosion regimes must be considered when analyzing origins of eroded lava channels that cut through strata of different strengths.