Investigation of the force balance in the Titan ionosphere: Cassini T5 flyby model/data comparisons

Cassini’s Titan flyby on 16 April, 2005 (T5) is the only encounter when the two main ionizing sources of the moon’s atmosphere, solar radiation and corotating plasma, align almost anti-parallel. In this paper a single-fluid multi-species 3D MHD model of the magnetospheric plasma interaction for T5 conditions is analyzed. Model results are compared to observations to investigate the ionospheric dynamics at Titan as well as to understand the deviations from a typical solar wind interaction, such as Venus’ interaction with the solar wind. Model results suggest that for the T5 interaction configuration, corotating plasma is the dominant driver determining the global interaction features at high altitudes. In the lower ionosphere below ～1500 km altitude – where the control of the ionospheric composition transfers from dynamic to chemical processes – magnetic and thermal pressure gradients oppose each other locally, complicating the ionospheric dynamics. Model results also imply that the nightside ionosphere – produced only by the impact ionization in the model – does not provide enough thermal pressure to balance the incident plasma dynamic pressure. As a result, the induced magnetic barrier penetrates into the ionosphere by plasma convection down to ～1000 km altitude and by magnetic diffusion below this altitude. Moreover, strong horizontal drag forces due to ion-neutral collisions and comparable drag forces estimated from possible neutral winds in the lower ionosphere below ～1400 km altitude oppose over local regions, implying that the Titan interaction must be treated as a 3D problem. Ion and electron densities calculated from the model generally agree with the Cassini Ion Neutral Mass Spectrometer and Langmuir probe measurements; however, there are significant differences between the calculated and measured magnetic fields. We discuss possible explanations for the discrepancy in the magnetic field predictions.     

Analytic theory of Titan’s Schumann resonance: Constraints on ionospheric conductivity and buried water ocean

This study presents an approximate model for the atypical Schumann resonance in Titan’s atmosphere that accounts for the observations of electromagnetic waves and the measurements of atmospheric conductivity performed with the Huygens Atmospheric Structure and Permittivity, Wave and Altimetry (HASI–PWA) instrumentation during the descent of the Huygens Probe through Titan’s atmosphere in January 2005. After many years of thorough analyses of the collected data, several arguments enable us to claim that the Extremely Low Frequency (ELF) wave observed at around 36 Hz displays all the characteristics of the second harmonic of a Schumann resonance. On Earth, this phenomenon is well known to be triggered by lightning activity. Given the lack of evidence of any thunderstorm activity on Titan, we proposed in early works a model based on an alternative powering mechanism involving the electric current sheets induced in Titan’s ionosphere by the Saturn’s magnetospheric plasma flow. The present study is a further step in improving the initial model and corroborating our preliminary assessments. We first develop an analytic theory of the guided modes that appear to be the most suitable for sustaining Schumann resonances in Titan’s atmosphere. We then introduce the characteristics of the Huygens electric field measurements in the equations, in order to constrain the physical parameters of the resonating cavity. The latter is assumed to be made of different structures distributed between an upper boundary, presumably made of a succession of thin ionized layers of stratospheric aerosols spread up to 150 km and a lower quasi-perfect conductive surface hidden beneath the non-conductive ground. The inner reflecting boundary is proposed to be a buried water–ammonia ocean lying at a likely depth of 55–80 km below a dielectric icy crust. Such estimate is found to comply with models suggesting that the internal heat could be transferred upwards by thermal conduction of the crust, while convective processes cannot be ruled out.     

Titan’s aerosol and stratospheric ice opacities between 18 and 500 μm: Vertical and spectral characteristics from Cassini CIRS

Vertical distributions and spectral characteristics of Titan’s photochemical aerosol and stratospheric ices are determined between 20 and 560 cm^?1 (500–18 μm) from the Cassini Composite Infrared Spectrometer (CIRS). Results are obtained for latitudes of 15°N, 15°S, and 58°S, where accurate temperature profiles can be independently determined.In addition, estimates of aerosol and ice abundances at 62°N relative to those at 15°S are derived. Aerosol abundances are comparable at the two latitudes, but stratospheric ices are ～3 times more abundant at 62°N than at 15°S. Generally, nitrile ice clouds (probably HCN and HC₃N), as inferred from a composite emission feature at ～160 cm^?1, appear to be located over a narrow altitude range in the stratosphere centered at ～90 km. Although most abundant at high northern latitudes, these nitrile ice clouds extend down through low latitudes and into mid southern latitudes, at least as far as 58°S.There is some evidence of a second ice cloud layer at ～60 km altitude at 58°S associated with an emission feature at ～80 cm^?1. We speculate that the identify of this cloud may be due to C₂H₆ ice, which in the vapor phase is the most abundant hydrocarbon (next to CH₄) in the stratosphere of Titan.Unlike the highly restricted range of altitudes (50–100 km) associated with organic condensate clouds, Titan’s photochemical aerosol appears to be well-mixed from the surface to the top of the stratosphere near an altitude of 300 km, and the spectral shape does not appear to change between 15°N and 58°S latitude. The ratio of aerosol-to-gas scale heights range from 1.3–2.4 at about 160 km to 1.1–1.4 at 300 km, although there is considerable variability with latitude. The aerosol exhibits a very broad emission feature peaking at ～140 cm^?1. Due to its extreme breadth and low wavenumber, we speculate that this feature may be caused by low-energy vibrations of two-dimensional lattice structures of large molecules. Examples of such molecules include polycyclic aromatic hydrocarbons (PAHs) and nitrogenated aromatics.Finally, volume extinction coefficients Nχ_E derived from 15°S CIRS data at a wavelength of λ = 62.5 μm are compared with those derived from the 10°S Huygens Descent Imager/Spectral Radiometer (DISR) data at 1.583 μm. This comparison yields volume extinction coefficient ratios Nχ_E(1.583 μm)/Nχ_E(62.5 μm) of roughly 70 and 20, respectively, for Titan’s aerosol and stratospheric ices. The inferred particle cross-section ratios χ_E(1.583 μm)/χ_E(62.5 μm) appear to be consistent with sub-micron size aerosol particles, and effective radii of only a few microns for stratospheric ice cloud particles.     

The solar wind interaction with Venus through the eyes of the Pioneer Venus Orbiter

Venus has no internal magnetic dynamo and thus its ionosphere and hot oxygen exosphere dominate the interaction with the solar wind. The solar wind at 0.72 AU has a dynamic pressure that ranges from 4.5 nPa (at solar max) to 6.6 nPa (at solar min), and its flow past the planet produces a shock of typical magnetosonic Mach number 5 at the subsolar point. At solar maximum the pressure in the ionospheric plasma is sufficient to hold off the solar wind at an altitude of 400 km above the surface at the subsolar point, and 1000 km above the terminators. The deflection of the solar wind occurs through the formation of a magnetic barrier on the inner edge of the magnetosheath, or shocked solar wind. Under typical solar wind conditions the time scale for diffusion of the magnetic field into the ionosphere is so long that the ionosphere remains field free and the barrier deflects almost all the incoming solar wind. Any neutral atoms of the hot oxygen exosphere that reach the altitude of the magnetosheath are accelerated by the electric field of the flowing magnetized plasma and swept along cycloidal paths in the antisolar direction. This pickup process, while important for the loss of the Venus atmosphere, plays a minor role in the deceleration and deflection of the solar wind. Like at magnetized planets, the Venus shock and magnetosheath generate hot electrons and ions that flow back along magnetic field lines into the solar wind to form a foreshock. A magnetic tail is created by the magnetic flux that is slowed in the interaction and becomes mass-loaded with thermal ions.The structure of the ionosphere is very much dependent on solar activity and the dynamic pressure of the solar wind. At solar maximum under typical solar wind conditions, the ionosphere is unmagnetized except for the presence of thin magnetic flux ropes. The ionospheric plasma flows freely to the nightside forming a well-developed night ionosphere. When the solar wind pressure dominates over the ionospheric pressure the ionosphere becomes completely magnetized, the flow to the nightside diminishes, and the night ionosphere weakens. Even at solar maximum the night ionosphere has a very irregular density structure. The electromagnetic environment of Venus has not been well surveyed. At ELF and VLF frequencies there is noise generated in the foreshock and shock. At low altitude in the night ionosphere noise, presumably generated by lightning, can be detected. This paper reviews the plasma environment at Venus and the physics of the solar wind interaction on the threshold of a new series of Venus exploration missions.     

The growth of wind-waves in Titan’s hydrocarbon seas

Ocean wave growth on Titan is considered. The classic Sverdrup–Munk theory for terrestrial wave growth is applied to Titan, and is compared with a simple energy balance model that exposes the effect of Titan’s environmental parameters (air density, gravity, and fluid density). These approaches are compared with the only previously-published (semi-empirical) model (Ghafoor, N.A.-L., Zarnecki, J.C., Challenor, P., Srokosz, M.A. [2000] J. Geophys. Res. 105, 12,077–12,091, hereafter G2k), and allow the impact of various parameters such as atmospheric density to be transparently explored.Our model, like G2k, suggests fully-developed significant wave heights on Titan H_s = 0.2 U², where U is the windspeed (SI units): in dimensionless terms this is rather close to H_s = 0.2 U²/g, a rule of thumb previously noted for terrestrial waves (we find various datasets where the prefactor varies by ～2). It is noted that liquid and air densities affect the growth rate of waves, but not their fully-developed height: for 1 m/s winds wave amplitude reaches 0.15 m (75% of fully-developed) with a fetch of only 1 km, rather faster than predicted by G2k. Liquid viscosity has no major effect on gravity wave growth, but does influence the threshold windspeed at which gravity–capillary waves form in the first place.The model is used to develop predicted ranges for wave height to guide the design of the Titan Mare Explorer (TiME), a proposed Discovery-class mission to float a capsule on Ligeia Mare in 2023. For the expected maximum 1 m/s winds, a significant wave height of 0.2 m and wavelength of ～4 m can be expected. Assuming that wave heights follow Rayleigh statistics as they do on Earth, then given the wave period of ～4 s, individual waves of ～0.6 m might be encountered over a 3 month period.For predicted Titan winds at Kraken Mare, significant wave heights may reach ～0.6 m in the peak of summer but do not exceed the tidal amplitude at its northern end, consistent with the area around Mayda Insula being a tidal flat, while elsewhere on Kraken and Ligeia and at Ontario Lacus, shorelines may be wave- or tidally-dominated, depending on the specific location.     

Winds and tides of Ligeia Mare,with application to the drift of the proposed time TiME (Titan Mare Explorer) capsule

We use two independent General Circulation Models (GCMs) to estimate surface winds at Titan’s Ligeia Mare (78° N, 250° W), motivated by a proposed mission to land a floating capsule in this ～500 km hydrocarbon sea. The models agree on the overall magnitude (～0.5–1 m/s) and seasonal variation (strongest in summer) of windspeeds, but details of seasonal and diurnal variation of windspeed and direction differ somewhat, with the role of surface exchanges being more significant than that of gravitational tides in the atmosphere. We also investigate the tidal dynamics in the sea using a numerical ocean dynamics model: assuming a rigid lithosphere, the tidal amplitude is up to ～0.8 m. Tidal currents are overall proportional to the reciprocal of depth—with an assumed central depth of 300 m, the characteristic tidal currents are ～1 cm/s, with notable motions being a slosh between Ligeia’s eastern and western lobes, and a clockwise flow pattern.We find that a capsule will drift at approximately one tenth of the windspeed, unless measures are adopted to augment the drag areas above or below the waterline. Thus motion of a floating capsule is dominated by the wind, and is likely to be several km per Earth day, a rate that will be readily measured from Earth by radio navigation methods. In some instances, the wind vector rotates diurnally such that the drift trajectory is epicyclic.     

Is the primordial crust of Mars magnetized?

A meteorite impact capable of creating a 200 km diameter crater can demagnetize the entire crust beneath, and produce an appreciable magnetic anomaly at satellite altitudes of ～400 km in case the pre-existing crust is magnetized. In this study we examine the magnetic field over all of the craters and impact-related Quasi-Circular Depressions (QCDs) with diameters larger than 200 km that are located on the highlands of Mars, excluding the Tharsis bulge, in order to estimate the mean magnetization of the highland crust. Using the surface topography and the gravity of Mars we first identify those QCDs that are likely produced by impacts. The magnetic map of a given crater or impact-related QCD is derived using the Mars Global Surveyor high-altitude nighttime radial magnetic data. Two extended ancient areas are identified on the highlands, the South Province and the Tempe Terra, which have large number of craters and impact-related QCDs but none of them has an appreciable magnetic signature. The primordial crust of these areas is not magnetized, or is very weakly magnetized at most. We examine some plausible scenarios to explain the weak magnetization of these areas, and conclude that no strong dynamo existed in the first ～100 Myr of Mars’ history when the newly formed primordial crust was cooling below the magnetic blocking temperatures of its minerals.     

Titan’s global crater population: A new assessment

We report a revised crater population for Titan using Cassini RADAR data through January 2010 (flyby T65), and make a size-dependent correction for the incomplete coverage (～33%) using a Monte-Carlo model. Qualitatively, Titan’s landscape is more heavily cratered than Earth, but much less than Mars or Ganymede: the area fraction covered by craters is in fact comparable with that of Venus. Quantitative efforts to interpret crater densities for Titan as surface age have been confounded by widely divergent crater production rates proposed in the literature. We elucidate the specific model assumptions that lead to these differences (assumed projectile density, scaling function for simple crater diameter, and complex crater size exponent) and suggest these are reasonable bounding models, with the Korycansky and Zahnle (2005) model representing a crater retention age of ～1 Ga, and the Artemieva and Lunine (2005) model representing a crater retention age of ～200 Ma. These estimates are consistent with models of Titan’s evolution that predict a thickening of its crust 0.3–1.2 Gyr ago.     

Hydrogen density in the dayside venusian exosphere derived from Lyman-α observations by SPICAV on Venus Express

A series of observations of the venusian hydrogen corona made by SPICAV on Venus Express are analyzed to estimate the amount of hydrogen in the exosphere of Venus. These observations were made between November 2006 and July 2007 at altitudes from 1000 km to 8000 km on the dayside. The Lyman-α brightness profiles derived are reproduced by the sum of a cold hydrogen population dominant below ～2000 km and a hot hydrogen population dominant above ～4000 km. The temperature (～300 K) and hydrogen density at 250 km (～10⁵ cm^?3) derived for the cold populations, near noon, are in good agreement with previous observations. Strong dawn–dusk exospheric asymmetry is observed from this set of observations, with a larger exobase density on the dawn side than on the dusk side, consistent with asymmetry previously observed in the venusian thermosphere, but with a lower dawn/dusk contrast. The hot hydrogen density derived is very sensitive to the sky background estimate, but is well constrained near 5000 km. The density of the hot population is reproduced by the exospheric model from Hodges (Hodges, R.R. [1999]. J. Geophys. Res. 104, 8463–8471) in which the hot population is produced by neutral–ions interactions in the thermosphere of Venus.     

Dynamical properties of the Venus mesosphere from the radio-occultation experiment VeRa onboard Venus Express

The dynamics of Venus’ mesosphere (60–100 km altitude) was investigated using data acquired by the radio-occultation experiment VeRa on board Venus Express. VeRa provides vertical profiles of density, temperature and pressure between 40 and 90 km of altitude with a vertical resolution of few hundred meters of both the Northern and Southern hemisphere. Pressure and temperature vertical profiles were used to derive zonal winds by applying an approximation of the Navier–Stokes equation, the cyclostrophic balance, which applies well on slowly rotating planets with fast zonal winds, like Venus and Titan. The main features of the retrieved winds are a midlatitude jet with a maximum speed up to 140 ± 15 m s^?1 which extends between 20°S and 50°S latitude at 70 km altitude and a decrease of wind speed with increasing height above the jet. Cyclostrophic winds show satisfactory agreement with the cloud-tracked winds derived from the Venus Monitoring Camera (VMC/VEx) UV images, although a disagreement is observed at the equator and near the pole due to the breakdown of the cyclostrophic approximation. Knowledge of both temperature and wind fields allowed us to study the stability of the atmosphere with respect to convection and turbulence. The Richardson number Ri was evaluated from zonal field of measured temperatures and thermal winds. The atmosphere is characterised by a low value of Richardson number from ～45 km up to ～60 km altitude at all latitudes that corresponds to the lower and middle cloud layer indicating an almost adiabatic atmosphere. A high value of Richardson number was found in the region of the midlatitude jet indicating a highly stable atmosphere. The necessary condition for barotropic instability was verified: it is satisfied on the poleward side of the midlatitude jet, indicating the possible presence of wave instability.     

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.     

Thermospheric/mesospheric temperatures on Venus: Results from ground-based high-resolution spectroscopy of CO2 in 1990/1991 and comparison to results from 2009 and between other techniques

We report temperatures in Venus’ upper mesosphere/lower thermosphere, deduced from reanalyzing very high resolution infrared spectroscopy of CO₂ emission lines acquired in 1990 and 1991. Kinetic temperatures at ～110 km altitude (0.15 Pa) are derived from the Doppler width of fully-resolved single line profiles measured near 10.4 μm wavelength using the NASA GSFC Infrared Heterodyne Spectrometer (IRHS) at the NASA IRTF on Mauna Kea, HI, close to Venus inferior conjunction and two Venus solstices. Measured temperatures range from ～200 to 240 K with uncertainty typically less than 10 K. Temperatures retrieved from similar measurement in 2009 using the Cologne Tuneable Heterodyne Infrared Spectrometer (THIS) at the NOAO McMath Telescope at Kitt Peak, AZ are 10–20 K lower. Temperatures retrieved more recently from the SOIR instrument on Venus EXpress are consistent with these results when the geometry of observation is accounted for. It is difficult to compare ground-based sub-mm retrievals extrapolated to 110 km due to their much larger field of view, which includes the night side regions not accessible to infrared heterodyne observations. Temperature variability appears to be high on day-to-day as well as longer timescales. Observed short term and long term variability may be attributed to atmospheric dynamics, diurnal variability and changes over solar activity and seasons. The Venus International Reference Atmosphere (VIRA) model predicts cooler temperatures at the sampled altitudes in the lower thermosphere/upper mesosphere and is not consistent with these measurements.     

Small-scale variability in Saturn’s lower ionosphere

We perform and present a wavelet analysis on all 31 Cassini electron density profiles published to date (Nagy, A.F. et al. [2006]. J. Geophys. Res. 111 (A6), CiteID A06310; Kliore, A.J. et al. [2009]. J. Geophys. Res. 114 (A4), CiteID A04315). We detect several discrete scales of variability present in the observations. Small-scale variability (S < 700 km) is observed in almost all data sets at different latitudes, both at dawn and dusk conditions. The most typical scale of variability is 300 km with scales between 200 km and 450 km being commonly present in the vast majority of the profiles. A low latitude dawn/dusk asymmetry is noted in the prevalent scales with the spectrum peaking sharply at the 300 km scale at dusk conditions and being broader at dawn conditions. Compared to dawn conditions the dusk ionosphere also shows more significant variability at the 100 km scale. The 300 km vertical scale is also present in the few available profiles from the northern hemisphere. Early observations from 2005 show a dominant scale at 350 km whereas later in 2007–2008 the spectrum shifts to the shorter scales with the most prominent scale being 300 km. The performed wavelet analysis and the obtained results are independent of assumptions about the nature of the layers and do not require a definition for a “background” electron density profile.In the second part of the paper we present a gravity wave propagation/dissipation model for Saturn’s upper atmosphere and compare the wave properties to the characteristics of the observed electron density variability at different scales. The general features observed in the data are consistent with gravity waves being present in the lower ionosphere and causing layering of the ions and the electrons. The wave-driving mechanism provides a simultaneous explanation for several of the properties of the observed variability: (i) lack of variability in the electron density above the predicted region of wave dissipation; (ii) in most cases the peak amplitude of variability occurs within the altitude range for dissipation of gravity waves or below; (iii) shorter scales have smaller amplitudes than the longer scales; (iv) shorter scales are present at lower altitudes whereas longer scales persist to higher altitudes; and (v) several layers often form a system of equally spaced maxima and minima that can be traced over a large altitude range.     

Optical constants of Titan tholins at mid-infrared wavelengths (2.5–25 μm) and the possible chemical nature of Titan’s haze particles

Complex organic materials may exist as haze layers in the atmosphere of Titan and as dark coloring agents on icy satellite surfaces. Laboratory measurements of optical constants of plausible complex organic materials are necessary for quantitative evaluation from remote sensing observations, and to document the existence of complex organic materials in the extraterrestrial environments. The recent Cassini VIMS and CIRS observations provide new constraints on Titan’s haze properties in the mid-infrared wavelength region. Here, we present the optical constants (2.5–25 μm) of Titan tholins generated with cold plasma irradiation of a N₂/CH₄ (90/10) gas mixture at pressures of 0.26 mbar, 1.6 mbar, and 23 mbar. Our new optical constants of three types of Titan tholins suggest that no single Titan tholin in this study fulfills all the observational constraints of the Titan haze material. The discrepancy remains a challenge for future modeling and laboratory efforts that aim toward a better understanding of Titan’s haze material.     

Titan’s internal structure and the evolutionary consequences

Titan’s moment of inertia (MoI), estimated from the quadrupole gravity field measured by the Cassini spacecraft, is 0.342, which has been interpreted as evidence of a partially differentiated internal mass distribution. It is shown here that the observed MoI is equally consistent with a fully differentiated internal structure comprising a shell of water ice overlying a low-density silicate core; depending on the chemistry of Titan’s subsurface ocean, the core radius is between 1980 and 2120 km, and its uncompressed density is 2570–2460 kg m^?3, suggestive of a hydrated CI carbonaceous chondrite mineralogy. Both the partially differentiated and fully differentiated hydrated core models constrain the deep interior to be several hundred degrees cooler than previously thought. I propose that Titan has a warm wet core below, or buffered at, the high-pressure dehydration temperature of its hydrous constituents, and that many of the gases evolved by thermochemical and radiogenic processes in the core (such as CH₄ and ⁴⁰Ar, respectively) diffuse into the icy mantle to form clathrate hydrates, which in turn may provide a comparatively impermeable barrier to further diffusion. Hence we should not necessarily expect to see a strong isotopic signature of serpentinization in Titan’s atmosphere.     

Optical constants of Titan’s stratospheric aerosols in the 70–1500 cm−1 spectral range constrained by Cassini/CIRS observations

We utilized aerosol extinction coefficient inferred from Cassini/CIRS spectra in the far and mid infrared region to derive the extinction cross-section near an altitude of 190 km at 15°S (from far-IR) and 20°S (from mid-IR). By comparing the extinction cross section that are derived from observations with theoretical calculations for a fractal aggregate of 3000 monomers, each having a radius of 0.05 μm, and a fractal dimension of 2, we are able to constrain the refractive index of Titan’s aerosol between 70 and 1500 cm^?1 (143 and 6.7 μm). As the real and imaginary parts of the refractive index are related by the Kramers–Kronig equation, we apply an iterative process to determine the optical constants in the thermal infrared. The resulting spectral dependence of the imaginary index displays several spectral signatures, some of which are also seen for some Titan’s aerosol analogues (tholins) produced in laboratory experiments. We find that Titan’s aerosols are less absorbent than tholins in the thermal infrared. The most prominent emission bands observed in the mid-infrared are due to CH bending vibrations in methyl and methylene groups. It appears that Titan’s aerosols predominantly display vibrations implying carbon and hydrogen atoms and perhaps marginally nitrogen. In the mid infrared, all the aerosol spectral signatures are observed at three additional latitudes (56°S, 5°N and 30°N) and in the 193–274 km altitude range, which implies that Titan’s aerosols exhibit the same chemical composition in all investigated latitude and altitude regions.     

Optical constants from 370 nm to 900 nm of Titan tholins produced in a low pressure RF plasma discharge

Determining the optical constants of Titan aerosol analogues, or tholins, has been a major concern for the last three decades because they are essential to constrain the numerical models used to analyze Titan’s observational data (albedo, radiative transfer, haze vertical profile, surface contribution, etc.). Here we present the optical constant characterization of tholins produced with an RF plasma discharge in a (95%N₂–5%CH₄) gas mixture simulating Titan’s main atmospheric composition, and deposited as a thin film on an Al–SiO₂ substrate. The real and imaginary parts, n and k, of the tholin complex refractive index have been determined from 370 nm to 900 nm wavelength using spectroscopic ellipsometry. The values of n decrease from n = 1.64 (at 370 nm) to n = 1.57 (at 900 nm) as well as the values of k which feature two behaviors: an exponential decay from 370 nm to 500 nm, with k = 12.4 × e^?0.018λ (where λ is expressed in nm), followed by a plateau, with k = (1.8 ± 0.2) × 10^?3. The trends observed for the PAMPRE tholins optical constants are compared to those determined for other Titan tholins, as well as to the optical constants of Titan’s aerosols retrieved from observational data.     

Ionization by cosmic rays of the atmosphere of Titan

Cosmic ray radiation is the main mechanism for ionizing the lower atmosphere of Titan. Their higher penetration power, in comparison with solar photons, allows cosmic rays to penetrate deep into the atmosphere of Titan, ionizing the neutral molecules and generating an ionosphere with an electron density peak, placed at around 90 km, similar in magnitude to the ionospheric peak produced by solar radiation in the upper atmosphere. In the lower atmosphere, the electron density profile, in the absence of a magnetic field, depends mainly on the modulation of cosmic rays by the solar wind and on the nature of the ionizable particles. We present here the first results of a new numerical model developed to calculate the concentration of electrons and most abundant ions in the Titan lower atmosphere. The present knowledge of Titan’s atmosphere permits us to include new neutral and ionic species, such as oxygen derivates, in a more detailed ion-chemistry calculation than previous lower ionospheric models of Titan. The electron density peaks at 90 km with a magnitude of 2150 cm⁻³. The ion distribution obtained predicts that cluster cations and hydrocarbon cations are the most abundant ions below and above the electron density peak, respectively. We also discuss the effect of solar activity at the distance of the Saturn orbit on the spectrum of the cosmic particles. We obtain that from solar minimum to solar maximum the ionization rate at the energy deposition peak changes by a factor of 1.2 at 70 km, and by a factor of 2.6 at altitudes as high as 400 km. The electron density at the concentration peak changes by a factor of 1.1 at 90 km, and by a factor of 1.6 at 400 km.     

Radiative transfer analyses of Titan’s tropical atmosphere

Titan’s optical and near-IR spectra result primarily from the scattering of sunlight by haze and its absorption by methane. With a column abundance of 92 km amagat (11 times that of Earth), Titan’s atmosphere is optically thick and only ～10% of the incident solar radiation reaches the surface, compared to 57% on Earth. Such a formidable atmosphere obstructs investigations of the moon’s lower troposphere and surface, which are highly sensitive to the radiative transfer treatment of methane absorption and haze scattering. The absorption and scattering characteristics of Titan’s atmosphere have been constrained by the Huygens Probe Descent Imager/Spectral Radiometer (DISR) experiment for conditions at the probe landing site (Tomasko, M.G., Bézard, B., Doose, L., Engel, S., Karkoschka, E. [2008a]. Planet. Space Sci. 56, 624–247; Tomasko, M.G. et al. [2008b]. Planet. Space Sci. 56, 669–707). Cassini’s Visual and Infrared Mapping Spectrometer (VIMS) data indicate that the rest of the atmosphere (except for the polar regions) can be understood with small perturbations in the high haze structure determined at the landing site (Penteado, P.F., Griffith, C.A., Tomasko, M.G., Engel, S., See, C., Doose, L., Baines, K.H., Brown, R.H., Buratti, B.J., Clark, R., Nicholson, P., Sotin, C. [2010]. Icarus 206, 352–365). However the in situ measurements were analyzed with a doubling and adding radiative transfer calculation that differs considerably from the discrete ordinates codes used to interpret remote data from Cassini and ground-based measurements. In addition, the calibration of the VIMS data with respect to the DISR data has not yet been tested. Here, VIMS data of the probe landing site are analyzed with the DISR radiative transfer method and the faster discrete ordinates radiative transfer calculation; both models are consistent (to within 0.3%) and reproduce the scattering and absorption characteristics derived from in situ measurements. Constraints on the atmospheric opacity at wavelengths outside those measured by DISR, that is from 1.6 to 5.0 μm, are derived using clouds as diffuse reflectors in order to derive Titan’s surface albedo to within a few percent error and cloud altitudes to within 5 km error. VIMS spectra of Titan at 2.6–3.2 μm indicate not only spectral features due to CH₄ and CH₃D (Rannou, P., Cours, T., Le Mouélic, S., Rodriguez, S., Sotin, C., Drossart, P., Brown, R. [2010]. Icarus 208, 850–867), but also a fairly uniform absorption of unknown source, equivalent to the effects of a darkening of the haze to a single scattering albedo of 0.63 ± 0.05. Titan’s 4.8 μm spectrum point to a haze optical depth of 0.2 at that wavelength. Cloud spectra at 2 μm indicate that the far wings of the Voigt profile extend 460 cm^?1 from methane line centers. This paper releases the doubling and adding radiative transfer code developed by the DISR team, so that future studies of Titan’s atmosphere and surface are consistent with the findings by the Huygens Probe. We derive the surface albedo at eight spectral regions of the 8 × 12 km² area surrounding the Huygens landing site. Within the 0.4–1.6 μm spectral region our surface albedos match DISR measurements, indicating that DISR and VIMS measurements are consistently calibrated. These values together with albedos at longer 1.9–5.0 μm wavelengths, not sampled by DISR, resemble a dark version of the spectrum of Ganymede’s icy leading hemisphere. The eight surface albedos of the landing site are consistent with, but not deterministic of, exposed water ice with dark impurities.     

Non-thermal water loss of the early Mars: 3D multi-ion hybrid simulations

In this study we analyze the non-thermal loss rates of O⁺, O₂⁺ and CO₂⁺ ions over the last 4.5 billion years (Gyr) in the Martian history by using a 3D hybrid model. For this reason we derived the past solar wind conditions in detail. We take into account the intensified particle flux of the early Sun as well as an Martian atmosphere, which was exposed to a sun's extreme ultraviolet (EUV) radiation flux 4.5 Gyr ago that was 100 times stronger than today. Furthermore, we model the evolution of the interplanetary magnetic field by a Weber & Davis solar wind model. The ‘external’ influences of the Sun's radiation flux and solar wind flux lead to the formation of an ionospheric obstacle by photoionization, charge exchange and electron impact. For the early Martian conditions we could show that charge exchange was the dominant ionization mechanism. Several hybrid simulations for different stages in the evolution of the Martian atmosphere, at 1, 2, 5, 10, 30 and 100 EUV, were performed to analyze the non-thermal escape processes by ion pick-up, momentum transfer from the solar wind to the ionosphere and detached ionospheric plasma clouds. Our results show a non-linear evolution of the loss rates. Using mean solar wind parameters the simulations result in an oxygen loss equivalent to the depth of a global Martian ocean of about 2.6 m over the last 4.5 Gyr. The induced magnetic field strength could be increased up to about 2000 nT. A simulation run with high solar wind density results in an oxygen loss of a Martian ocean up to 205 m depth during 150 million years after the sun reached the zero age mean sequence (ZAMS).