Characterizing Io’s Pele,Tvashtar and Pillan plumes: Lessons learned from Hubble

Hubble Space Telescope/Wide Field and Planetary Camera 2 (HST/WFPC2) images of Io obtained between 1995 and 2007 between 0.24 and 0.42 μm led to the detection of the Pele plume in reflected sunlight in 1995 and 1999; imaging of the Pele plume via absorption of jovian light in 1996 and 1999; detection of the Prometheus-type Pillan plume in reflected sunlight in 1997; and detection of the 2007 Pele-type Tvashtar plume eruption in reflected sunlight and via absorption of jovian light. Based on a detailed analysis of these observations we characterize and compare the gas and dust properties of each of the detected plumes. In each case, the brightness of the plumes in reflected sunlight is less at 0.26 μm than at 0.33 μm. Mie scattering analysis of the wavelength dependence of each plume’s reflectance signature suggests that range of particle sizes within the plumes is quite narrow. Assuming a normal distribution of particle sizes, the range of mean particle sizes is ～0.035–0.12 μm for the 1997 Pillan eruption, ～0.05–0.08 μm for the 1999 Pele and 2007 Tvasthar plumes, and ～0.05–0.11 μm for the 1995 Pele plume, and in each case the standard deviation in the particle size distribution is <15%. The Mie analysis also suggests that the 2007 Tvashtar eruption released ～10⁹ g of sulfur dust, the 1999 Pele eruption released ～10⁹ g of SO₂ dust, the 1997 Pillan eruption released ～10¹⁰ g of SO₂ dust, and the 1995 Pele plume may have released ～10¹⁰ g of SO₂ dust. Analysis of the plume absorption signatures recorded in the F255W filter bandpass (0.24–0.28 μm) indicates that the opacity of the 2007 Tvashtar plume was 2× that of the 1996 and 1999 Pele plume eruptions. While the sulfur dust density estimated for the Tvashtar from the reflected sunlight data could have produced 61% of the observed plume opacity, <10% of the 1999 Pele F255W plume opacity could have resulted from the SO₂ dust detected in the eruption. Accounting for the remaining F255W opacity level of the Pele and Tvasthar plumes based on SO₂ and S₂ gas absorption, the SO₂ and S₂ gas density inferred for each plume is almost equivalent corresponding to ～2–6 × 10¹⁶ cm^?2 and 3–5 × 10¹⁵ cm^?2, respectively, producing SO₂ and S₂ gas resurfacing rates ～0.04–0.2 cm yr^?1 and 0.007–0.01 cm yr^?1; and SO₂ and S₂ gas masses ～1–4 × 10¹⁰ g and ～2–3 × 10⁹ g; for a total dust to gas ratio in the plumes ～10^?1–10^?2. The 2007 Tvashtar plume was detected by HST at ～380 ± 40 km in both reflected sunlight and absorbed jovian light; in 1999, the detected Pele plume altitude was 500 km in absorbed jovian light, but in reflected sunlight the detected height was ～2× lower. Thus, for the 1999 Pele plume, similar to the 1979 Voyager Pele plume observations, the most efficient dust reflections occurred in the region closest to the plume vent. The 0.33–0.42 μm brightness of the 1997 Pillan plume was 10–20× greater than the Pele or Tvashtar plumes, exceeding by a factor of 3 the average brightness levels observed within 200 km of 1979 Loki eruption vent. But, the 0.26 μm brightness of the 1997 Pillan plume in reflected sunlight was significantly lower than would be predicted by the dust scattering model. Presuming that the 0.26 μm brightness of the 1997 Pillan plume was attenuated by the eruption plume’s gas component, then an SO₂ gas density ～3–6 × 10¹⁸ cm^?2 is inferred from the data (for S₂/SO₂ ratios ?4%), comparable to the 0.3–2 × 10¹⁸ cm^?2 SO₂ density detected at Loki in 1979 (Pearl, J.C. et al. [1979]. Nature 280, 755; Lellouch et al., 1992), and producing an SO₂ gas mass ～3–8 × 10¹¹ g and an SO₂ resurfacing rate ～8–23 cm yr^?1. These results confirm the connection between high (?10¹⁷ cm^?2) SO₂ gas content and plumes that scatter strongly at nearly blue wavelengths, and it validates the occurrence of high density SO₂ gas eruptions on Io. Noting that the SO₂ gas content inferred from a spectrum of the 2003 Pillan plume was significantly lower ～2 × 10¹⁶ cm^?2 (Jessup, K.L., Spencer, J., Yelle, R. [2007]. Icarus 192, 24–40); and that the Pillan caldera was flooded with fresh SO₂ frost/slush just prior to the 1997 Pillan plume eruption (Geissler, P., McEwen, A., Phillips, C., Keszthelyi, L., Spencer, J. [2004a]. Icarus 169, 29–64; Phillips, C.B. [2000]. Voyager and Galileo SSI Views of Volcanic Resurfacing on Io and the Search for Geologic Activity at Europa. Ph.D. Thesis, Univ. of Ariz., Tucson); we propose that the density of SO₂ gas released by this volcano is directly linked to the local SO₂ frost abundance at the time of eruption.     

Trapped Ar isotopes in meteorite ALH 84001 indicate Mars did not have a thick ancient atmosphere

Water is not currently stable in liquid form on the martian surface due to the present mean atmospheric pressure of ～7 mbar and mean global temperature of ～220 K. However, geomorphic features and hydrated mineral assemblages suggest that Mars’ climate was once warmer and liquid water flowed on the surface. These observations may indicate a substantially more massive atmosphere in the past, but there have been few observational constraints on paleoatmospheric pressures. Here we show how the ⁴⁰Ar/³⁶Ar ratios of trapped gases within martian meteorite ALH 84001 constrain paleoatmospheric pressure on Mars during the Noachian era [～4.56–3.8 billion years (Ga)]. Our model indicates that atmospheric pressures did not exceed ～1.5 bar during the first 400 million years (Ma) of the Noachian era, and were <400 mbar by 4.16 Ga. Such pressures of CO₂ are only sufficient to stabilize liquid water on Mars’ surface at low latitudes during seasonally warm periods. Other greenhouse gases like SO₂ and water vapor may have played an important role in intermittently stabilizing liquid water at higher latitudes following major volcanic eruptions or impact events.     

The impact of a weak south pole on thermal convection in Enceladus’ ice shell

Enceladus exhibits a strong hemispheric dichotomy of tectonism and heat flux, with geologically young, heavily tectonized terrains and a high heat flux in the South Polar Terrain (SPT) and relatively ancient terrains with presumably lower heat fluxes over the rest of the satellite. To understand the convective pattern and its relationship with surface tectonics, we present three-dimensional numerical models of convection in Enceladus’ ice shell including basal heating and tidal heating. Our thermal boundary conditions exhibit no north–south asymmetries, but because the tectonism at the SPT may weaken the ice there, we impose a mechanically weak lithosphere within the SPT. The weakening is parameterized by adopting a reduced viscosity contrast within the SPT. Without such a weak zone, convection (if any) resides in stagnant-lid mode and exhibits no hemispheric dichotomy. In the presence of such an SPT weak zone, however, we find vigorous convection in the ice underneath the SPT, with convective plumes rising close to the surface. In contrast, only stagnant lid convection, or no convection at all, occurs elsewhere over the satellite. Away from the SPT, the heat flux in our models is small (5–10 mW m^?2) and the surface strains are small enough to imply surface ages >10⁹ years. Within the SPT, however, our models yield peak heat fluxes of ～70–200 mW m^?2, implying heat flows integrated across the SPT of up to 5 GW, similar to that inferred from Cassini thermal observations. The surface strains in our models are high enough near the south pole to cause intense tectonism and imply surface ages of ～10⁶–10⁷ years, consistent with age estimates of the SPT.     

Constraints on mantle plumes on Venus: Implications for volatile history

The analysis of Venus’ gravity field and topography suggests the presence of a small number of deep mantle plumes (～9). This study predicts the number of plumes formed at the core–mantle boundary, their characteristics, and the production of partial melt from adiabatic decompression. Numerical simulations are performed using a 3D spherical code that includes large viscosity variations and internal heating. This study investigates the effect of several parameters including the core–mantle boundary temperature, the amount of internal heating, and the mantle viscosity. The smallest number of plumes is achieved when no internal heating is present. However, scaling Earth’s radiogenic heating to Venus suggests a value of ～16 TW. Cases with internal heating produce more realistic lid thickness and partial melting, but produce either too many plumes or no plumes if a high mantle temperature precludes the formation of a hot thermal boundary layer. Mantle viscosity must be reduced to at least 10²⁰ Pa s in order to include significant internal heating and still produce hot plumes. In all cases that predict melting, melting occurs throughout the upper mantle. Only cases with high core temperature (>1700 K) produce dry melting. Over time the upper mantle may have lost significant volatiles. Depending on the water content of the lower mantle, deep plumes may contribute to present-day atmospheric water via volcanic outgassing. Assuming 50 ppm water in mantle, 10 plumes with a buoyancy flux of 500 kg/s continuously erupting for 4 myr will outgas an amount of water on the order of that in the lower atmosphere. A higher level of internal heating than achieved to date, as well as relatively low mantle viscosity, may be required to achieve simulations with ～10 plumes and a thinner lid. Alternatively, if the mantle is heating up due to the stagnant lid, the effect is equivalent to having lower rates of internal heating. A temperature increase of 110 K/byr is equivalent to ?13 TW. This value along with the internal heating of 3 TW used in this study may represent the approximate heat budget of Venus’ mantle.     

A whole-moon thermal history model of Europa: Impact of hydrothermal circulation and salt transport

A whole-moon numerical model of Europa is developed to simulate its thermal history. The thermal evolution covers three phases: (i) an initial, roughly 0.5 Gyr-long period of radiogenic heating and differentiation, (ii) a long period from 0.5 Gyr to 4 Gyr with continuing radiogenic heating but no tidal dissipative heating (TDH), and (iii) a final period covering the last 0.5 Gyr until the present, during which TDH is active. Hydrothermal plumes develop after the initial period of heating and differentiation and transport heat and salt from Europa’s silicate mantle to its ice shell. We find that, even without TDH, vigorous hydrothermal convection in the rocky mantle can sustain flow in an ocean layer throughout Europa’s history. When TDH becomes active, the ice shell melts quickly to a thickness of about 20 km, leaving an ocean 80 km or more deep. Parameterized convection in the ice shell is non-uniform spatially, changes over time, and is tied to the deeper ocean–mantle dynamics. We also find that the dynamics are affected by salt concentrations. An initially non-uniform salt distribution retards plume penetration, but is homogenized over time by turbulent diffusion and time-dependent flow driven by initial thermal gradients. After homogenization, the uniformly distributed salt concentrations are no longer a major factor in controlling plume transport. Salt transport leads to the formation of a heterogeneous brine layer and salt inclusions at the bottom of the ice shell; the presence of salt in the ice shell could strongly influence convection in that layer.     

Reduced sulfur–carbon–water systems on Mars may yield shallow methane hydrate reservoirs

Methane clathrate hydrate reservoirs capped by overlying permafrost have been proposed as potential sources of atmospheric methane plumes on Mars. However, the surface flux of methane from hydrate dissociation is limited by the diffusion rate of methane through the overlying ice. Assuming hydrates underlay the entire plume footprint, the maximum diffusion path length is expected to be less than 15 m, depths too shallow to stabilize pure methane hydrates under Mars geothermal and lithostatic conditions at low to mid latitudes. Therefore, pure methane hydrates confined within permafrost could not produce methane surface fluxes of the magnitude observed near the equator. However, the addition of 10% H₂S, a secondary gas commonly associated with methane production on Earth, expands the hydrate stability field, with clathrates expected within 10 m of the surface at the equator and at depths less than 1 m at higher latitudes. This indicates that H₂S would also be expected to be released as well as methane if the plumes have a confined hydrate reservoir source.     

Cassini RPWS observations of dust in Saturn's E Ring

The Cassini radio and plasma wave science (RPWS) instrument is sensitive to few-micron dust grains impacting on the spacecraft at relative speeds of order 10 km/s. Through the first year or so of operations in orbit at Saturn, the RPWS has made a number of both inclined and equatorial crossings of the E ring, particularly near the orbit of Enceladus. Assuming water ice grains, the typical size particle detected by the RPWS has a radius of a few microns. Peak impact rates of about 50 s⁻¹ are found near the orbit of Enceladus corresponding to densities of order 5×10⁻⁴ m⁻³. The variation of dust fluxes as a function of height above or below the equator is well described by a Gaussian distribution with a scale height of about 2800 km although there is usually some non-Gaussian variation near the peak fluxes suggesting some structure in the core of the ring. Offsets of the peak number densities are typically of the order of a few hundred km from the geometric equator. A near-equatorial radial profile through the orbit of Enceladus shows a sharply peaked distribution at the orbit of the moon. A size distribution averaged over several passes through the orbit of Enceladus is determined which varies as m^−2.80. The peak in dust number density at the orbit of Enceladus is consistent with previous optical measurements and strongly supports the suggestion that Enceladus is a primary source for E ring particles.     

The Cassini Enceladus encounters 2005–2010 in the view of energetic electron measurements

The moon Enceladus, embedded in Saturn’s radiation belts, is the main internal source of neutral and charged particles in the Kronian magnetosphere. A plume of water ice molecules and dust released through geysers on the south polar region provides enough material to feed the E-ring and also the neutral torus of Saturn and the entire magnetosphere. In the time period 2005–2010 the Cassini spacecraft flew close by the moon 14 times, sometimes as low as 25 km above the surface and directly through the plume. For the very first time measurements of plasma and energetic particles inside the plume and its immediate vicinity could be obtained. In this work we summarize the results of energetic electron measurements in the energy range 27 keV to 21 MeV taken by the Low Energy Magnetospheric Measurement System (LEMMS), part of the Magnetospheric Imaging Instrument (MIMI) onboard Cassini in the vicinity of the moon in combination with measurements of the magnetometer instrument MAG and the Electron Spectrometer ELS of the plasma instrument CAPS onboard the spacecraft. Features in the data can be interpreted as that the spacecraft was connected to the plume material along field lines well before entering the high density region of the plume. Sharp absorption signatures as the result of losses of energetic electrons bouncing along those field lines, through the emitted gas and dust clouds, clearly depend on flyby geometry as well as on measured pitch angle/look direction of the instrument. We found that the depletion signatures during some of the flybys show “ramp-like” features where only a partial depletion has been observed further away from the moon followed by nearly full absorption of electrons closer in. We interpret this as partially/fully connected to the flux tube connecting the moon with Cassini. During at least two of the flybys (with some evidence of one additional encounter) MIMI/LEMMS data are consistent with the presence of dust in energetic electron data when Cassini flew directly through the south polar plume. In addition we found gradients in the magnetic field components which are frequently found to be associated with changes in the MIMI/LEMMS particles intensities. This indicates that complex electron drifts in the vicinity of Enceladus could form forbidden regions for electrons which may appear as intensity drop-outs.     

UV spectrum of Enceladus

We present a far ultraviolet (FUV) spectrum of Saturn’s moon Enceladus from the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). We have put upper limits on emission from C, N, and O lines in Enceladus’ atmosphere and column densities for the C lines assuming solar resonance scattering. We find these upper limits to be relatively low—on the order of tens to thousands of Rayleighs and with C column densities on the order of 10⁸–10¹⁵ cm^?2, depending on the assumed source size. We also present a segment of a reflectance spectrum in the FUV from ～1900–2130 Å. This region was sensitive to the different ice mixtures in the model spectra reported by Hendrix et al. (Hendrix, A.R., Hansen, C.J., Holsclaw, G.M. [2010]. Icarus, 206, 608). We find the spectrum brightens quickly longward of ～1900 Å, constraining the absorption band observed by Hendrix et al. from ～170 to 190 nm. We find our data is consistent with the suggestion of Hendrix et al. of the presence of ammonia ice (or ammonia hydrate) to darken that region, and also possibly tholins to darken the mid-UV, as reported by Verbiscer et al. (Verbiscer, A.J., French, R.G., McGhee, C.A. [2005]. Icarus, 173, 66).     

Pressure dependent trace gas trapping in amorphous water ice at 77 K: Implications for determining conditions of comet formation

The nature of cometary volatile materials is subject to debate. Theoretical models of cometary nuclei and laboratory studies suggest that these objects could be made of amorphous water ice in addition to other volatile molecules and refractory grains. This water ice structure has the ability to encapsulate the gases of surrounding environment, reflecting the physical and chemical conditions during their deposition. Therefore, the knowledge of the chemical composition of volatile molecules trapped in amorphous water ice provides a tool for probing the formation environment of cometary ice grains. Experimental studies of gas trapping efficiency in amorphous water ice have been previously conducted mostly under kinetic conditions, where dynamic pumping and temperature gradients prevented rigorous calibrations. In this work, we investigated the trapping efficiencies of Ar, CO, CH₄, Kr and N₂ by depositing water vapor as ice in the presence of trace gases in a volume submerged in liquid nitrogen at 77 K. The gas trapping efficiencies were determined simply by monitoring the pressure difference of the trace gases before and after the deposition of a known amount of water molecules as amorphous ice.Our results show that the trapped gas to water molecule ratio in amorphous ice is controlled primarily by the partial pressure of the gas during water ice deposition, and is independent of the ice deposition rate as well as the gas to water ratio in the vapor phase. The trapping efficiencies of gases decrease in the order of Kr > CH₄ > CO > Ar > N₂ in accordance with previous studies. Assuming that the water ice structure of comets is at least partially amorphous water ice at the time of their formation, these results suggest that the total pressure and composition of the surrounding environment of amorphous ice formation are significant controlling factors of trace gas concentrations in cometary ice. This further indicates that the evolution of the solar nebula and timing of cometary ice condensation can also be important parameters in linking the volatile contents of comets and their formation process.     

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.     

The formation heritage of Jupiter Family Comet 10P/Tempel 2 as revealed by infrared spectroscopy

We present spectral and spatial information for major volatile species in Comet 10P/Tempel 2, based on high-dispersion infrared spectra acquired on UT 2010 July 26 (heliocentric distance R_h = 1.44 AU) and September 18 (R_h = 1.62 AU), following the comet’s perihelion passage on UT 2010 July 04. The total production rate for water on July 26 was (1.90 ± 0.12) × 10²⁸ molecules s^?1, and abundances of six trace gases (relative to water) were: CH₃OH (1.58% ± 0.23%), C₂H₆ (0.39% ± 0.04%), NH₃ (0.83% ± 0.20%), and HCN (0.13% ± 0.02%). A detailed analysis of intensities for water emission lines provided a rotational temperature of 35 ± 3 K. The mean OPR is consistent with nuclear spin populations in statistical equilibrium (OPR = 3.01 ± 0.18), and the (1σ) lower bound corresponds to a spin temperature >38 K. Our measurements were contemporaneous with a jet-like feature observed at optical wavelengths. The spatial profiles of four primary volatiles display strong enhancements in the jet direction, which favors release from a localized vent on the nucleus. The measured IR continuum is much more sharply peaked and is consistent with a dominant contribution from the nucleus itself. The peak intensities for H₂O, CH₃OH, and C₂H₆ are offset by ～200 km in the jet direction, suggesting the possible existence of a distributed source, such as the release of icy grains that subsequently sublimed in the coma. On UT September 18, no obvious emission lines were present in our spectra, nevertheless we obtained a 3σ upper limit Q(H₂O) < 2.86 × 10²⁷ molecules s^?1.     

Io: Volcanic thermal sources and global heat flow

We have examined thermal emission from 240 active or recently-active volcanic features on Io and quantified the magnitude and distribution of their volcanic heat flow during the Galileo epoch. We use spacecraft data and a geological map of Io to derive an estimate of the maximum possible contribution from small dark areas not detected as thermally active but which nevertheless appear to be sites of recent volcanic activity. We utilize a trend analysis to extrapolate from the smallest detectable volcanic heat sources to these smallest mapped dark areas. Including the additional heat from estimates for “outburst” eruptions and for a multitude of very small (“myriad”) hot spots, we account for ～62 × 10¹² W (～59 ± 7% of Io’s total thermal emission). Loki Patera contributes, on average, 9.6 × 10¹² W (～9.1 ± 1%). All dark paterae contribute 45.3 × 10¹² W (～43 ± 5%). Although dark flow fields cover a much larger area than dark paterae, they contribute only 5.6 × 10¹² W (～5.3 ± 0.6%). Bright paterae contribute ～2.6 × 10¹² W (～2.5 ± 0.3%). Outburst eruption phases and very small hot spots contribute no more than ～4% of Io’s total thermal emission: this is probably a maximum value. About 50% of Io’s volcanic heat flow emanates from only 1.2% of Io’s surface. Of Io’s heat flow, 41 ± 7.0% remains unaccounted for in terms of identified sources. Globally, volcanic heat flow is not uniformly distributed. Power output per unit surface area is slightly biased towards mid-latitudes, although there is a stronger bias toward the northern hemisphere when Loki Patera is included. There is a slight favoring of the northern hemisphere for outbursts where locations were well constrained. Globally, we find peaks in thermal emission at ～315°W and ～105°W (using 30° bins). There is a minimum in thermal emission at around 200°W (almost at the anti-jovian longitude) which is a significant regional difference. These peaks and troughs suggest a shift to the east from predicted global heat flow patterns resulting from tidal heating in an asthenosphere. Global volcanic heat flow is dominated by thermal emission from paterae, especially from Loki Patera (312°W, 12°N). Thermal emission from dark flows maximises between 165°W and 225°W. Finally, it is possible that a multitude of very small hot spots, smaller than the present angular resolution detection limits, and/or cooler, secondary volcanic processes involving sulphurous compounds, may be responsible for at least part of the heat flow that is not associated with known sources. Such activity should be sought out during the next mission to Io.     

The outcome of the protoplanetary disk of very massive stars

We suggest that planets, brown dwarfs, and even low mass stars can be formed by fragmentation of protoplanetary disks around very massive stars (M ? 100 M_⊙). We discuss how fragmentation conditions make the formation of very massive planetary systems around very massive stars favorable. Such planetary systems are likely to be composed of brown dwarfs and low mass stars of ～0.1–0.3 M_⊙, at orbital separations of ～ few × 100–10⁴ AU. In particular, scaling from solar-like stars suggests that hundreds of Mercury-like planets might orbit very massive stars at ～10³ AU where conditions might favor liquid water. Such fragmentation objects can be excellent targets for the James Webb Space Telescope and other large telescopes working in the IR bands. We predict that deep observations of very massive stars would reveal these fragmentation objects, orbiting in the same orbital plane in cases where there are more than one object.     

Formation of spherules in impact produced vapor plumes

We have constructed a numerical model of spherule formation in an impact produced vapor plume. This model tracks the expansion of the vapor plume using a one-dimensional Lagrangian hydrocode coupled with the ANEOS equation of state for silica. We then include the equations for nucleation and growth as described by homogeneous nucleation theory to describe the process of spherule formation. We use this model to determine the number and size of the spherules that an impact creates. We also explore when and where spherules are formed in the vapor plume, and how this affects the size of the spherules. In general we find that smaller spherules form in the outer, faster moving, portions of the vapor plume at earlier times. This work also explores the effect of impactor size and impact velocity on the resultant spherule size. We report a simple linear dependence on impactor size and a complex dependence on impact velocity. We find that a 10 km diameter asteroid impacting at a velocity of ～21 km/s creates spherules that are ～250 μm in diameter which is comparable to the spherules found in the K/Pg boundary layer.     

Sulfur chemistry in the middle atmosphere of Venus

Venus Express measurements of the vertical profiles of SO and SO₂ in the middle atmosphere of Venus provide an opportunity to revisit the sulfur chemistry above the middle cloud tops (～58 km). A one dimensional photochemistry-diffusion model is used to simulate the behavior of the whole chemical system including oxygen-, hydrogen-, chlorine-, sulfur-, and nitrogen-bearing species. A sulfur source is required to explain the SO₂ inversion layer above 80 km. The evaporation of the aerosols composed of sulfuric acid (model A) or polysulfur (model B) above 90 km could provide the sulfur source. Measurements of SO₃ and SO (a¹Δ → X³Σ^-) emission at 1.7 μm may be the key to distinguish between the two models.     

Water vapor near the cloud tops of Venus from Venus Express/VIRTIS dayside data

Observations of the dayside of Venus performed by the high spectral resolution channel (–H) of the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) on board the ESA Venus Express mission have been used to measure the altitude of the cloud tops and the water vapor abundance around this level with a spatial resolution ranging from 100 to 10 km. CO₂ and H₂O bands between 2.48 and 2.60 μm are analyzed to determine the cloud top altitude and water vapor abundance near this level. At low latitudes (±40°) mean water vapor abundance is equal to 3 ± 1 ppm and the corresponding cloud top altitude at 2.5 μm is equal to 69.5 ± 2 km. Poleward from middle latitudes the cloud top altitude gradually decreases down to 64 km, while the average H₂O abundance reaches its maximum of 5 ppm at 80° of latitude with a large scatter from 1 to 15 ppm. The calculated mass percentage of the sulfuric acid solution in cloud droplets of mode 2 (～1 μm) particles is in the range 75–83%, being in even more narrow interval of 80–83% in low latitudes. No systematic correlation of the dark UV markings with the cloud top altitude or water vapor has been observed.     

Episodes of floods in Mangala Valles,Mars, from the analysis of HRSC,MOC and THEMIS images

The Mangala Valles is a 900-km long outflow channel system in the highlands adjacent to the south-eastern flank of the Tharsis bulge. This work was intended to answer the following two questions unresolved in previous studies: (1) Was there only one source of water (Mangala Fossa at the valley head which is one of the Medusae Fossae troughs or graben) or were other sources also involved in the valley-carving water supply, and (2) Was there only one episode of flooding (maybe with phases) or were there several episodes significantly separated in time. The geologic analysis of HRSC image 0286 and mapping supported by analysis of MOC and THEMIS images show that Mangala Valles was carved by water released from several sources. The major source was Mangala Fossa, which probably formed in response to magmatic dike intrusion. The graben cracked the cryosphere and permitted the release of groundwater held under hydrostatic pressure. This major source was augmented by a few smaller-scale sources at localities in (1) two mapped heads of magmatic dikes, (2) heads of two clusters of sinuous channels, and (3) probably several large knob terrain locals. The analysis of results of crater counts at more than 60 localities showed that the first episode of formation of Mangala Valles occurred ～3.5 Ga ago and was followed by three more episodes, one occurred ～1 Ga ago, another one ～0.5 Ga ago, and the last one ～0.2 Ga ago. East of the mapped area there are extended and thick lava flows whose source may be the eastern continuation of the Mangala source graben. Crater counts in 10 localities on these lava flows correlate with those taken on the Mangala valley elements supporting the idea that the valley head graben was caused by dike intrusions. Our observations suggest that the waning stage of the latest flooding episode (～0.2 Ga ago) led to the formation at the valley head of meander-like features sharing some characteristics with meanders of terrestrial rivers. If this analogy is correct this could suggest a short episode of global warming in Late Amazonian time.     

Sea level budget over 2003–2008: A reevaluation from GRACE space gravimetry,satellite altimetry and Argo

From the IPCC 4th Assessment Report published in 2007, ocean thermal expansion contributed by ～ 50% to the 3.1 mm/yr observed global mean sea level rise during the 1993–2003 decade, the remaining rate of rise being essentially explained by shrinking of land ice. Recently published results suggest that since about 2003, ocean thermal expansion change, based on the newly deployed Argo system, is showing a plateau while sea level is still rising, although at a reduced rate (～ 2.5 mm/yr). Using space gravimetry observations from GRACE, we show that recent years sea level rise can be mostly explained by an increase of the mass of the oceans. Estimating GRACE-based ice sheet mass balance and using published estimates for glaciers melting, we further show that ocean mass increase since 2003 results by about half from an enhanced contribution of the polar ice sheets – compared to the previous decade – and half from mountain glaciers melting. Taking also into account the small GRACE-based contribution from continental waters (< 0.2 mm/yr), we find a total ocean mass contribution of ～ 2 mm/yr over 2003–2008. Such a value represents ～ 80% of the altimetry-based rate of sea level rise over that period. We next estimate the steric sea level (i.e., ocean thermal expansion plus salinity effects) contribution from: (1) the difference between altimetry-based sea level and ocean mass change and (2) Argo data. Inferred steric sea level rate from (1) (～ 0.3 mm/yr over 2003–2008) agrees well with the Argo-based value also estimated here (0.37 mm/yr over 2004–2008). Furthermore, the sea level budget approach presented in this study allows us to constrain independent estimates of the Glacial Isostatic Adjustment (GIA) correction applied to GRACE-based ocean and ice sheet mass changes, as well as of glaciers melting. Values for the GIA correction and glacier contribution needed to close the sea level budget and explain GRACE-based mass estimates over the recent years agree well with totally independent determinations.     

Search for methane and upper limits to ethane and SO2 on Mars

The IRTF/CSHELL observations in February 2006 at L_S = 10° and 63–93°W show ～10 ppb of methane at 45°S to 7°N and ～3 ppb outside this region that covers the deepest canyon Valles Marineris. Observations in December 2009 at L_S = 20° and 0–30°W included spectra of the Moon at a similar airmass as a telluric calibrator. A technique for extraction of the martian methane line from a combination of the Mars and Moon spectra has been developed. The observations reveal no methane with an upper limit of 8 ppb. The results of both sessions agree with the observations by Mumma et al. (Mumma, M.J. et al. [2009]. Science 323, 1041–1045) at the same season in February 2006 and are smaller than those in the PFS and TES maps. Production and removal of the biological methane on Mars do not significantly change the redox state of the atmosphere and the balance of hydrogen. A search for ethane at 2977 cm^?1 results in an upper limit of 0.2 ppb. However, this limit does not help to establish the origin of methane on Mars. Reanalysis of our search for SO₂ using TEXES confirms the recently established upper limit of 0.3 ppb. Along with the lack of hot spots and gas vents with endogenic heat sources in the THEMIS observations, the very low upper limit to SO₂ on Mars does not favor geological methane that is less abundant than SO₂ in the outgassing from the terrestrial planets.