Scenarios for the formation of Chasma Boreale, Mars

The martian polar caps feature large chasmata and smaller trough systems which have no counterpart in terrestrial ice sheets. Chasma Boreale cuts about 500 km into the western part of the north-polar cap, is up to 100 km wide and up to 2 km deep. One possible formation mechanism is by a temporary heat source under the ice due to tectono-thermal or volcanic activity, which melts the ice from below. It is demonstrated by model simulations that this process is feasible, a moderately increased heat flux of 0.5-1 W m⁻², sustained over at least tens of thousands of years, producing a topographic depression which resembles the real chasma. Associated meltwater discharge rates are small (), but can exceed 10 km³ a⁻¹ if a stronger heat flux of 10 W m⁻² is assumed. Local ice-flow velocities during the process of chasma formation can exceed 1 m a⁻¹ at the head and scarps of the chasma. However, if the thermal anomaly shuts down, glacial flow quickly decreases, so that the chasma can stay open for an indefinite amount of time without an ongoing, sustaining process under the climate conditions of the most recent millions of years.     

Importance of aeolian processes in the origin of the north polar chasmata, Mars

A detailed examination of the location and orientation of sand dunes and other aeolian features within the north polar chasmata indicates that steep scarps strongly influence the direction and intensity of prevailing winds. These steep scarps are present at the heads and along the margins of the north polar chasmata. Topographic profiles of the arcuate head scarps and equator-facing wall of Chasma Boreale reveal unusually steep polar slopes ranging from ∼6°-30°. The relatively steep-sloped (∼8°), sinuous scarp at the head of two smaller chasmata, located west of Chasma Boreale, exhibits an obvious resistant cap-forming unit. Scarp retreat is occurring in places where the cap unit is actively being undercut by descending slope winds. Low-albedo surfaces lacking sand dunes or dust mantles are present at the base of the polar scarps. A ∼100-300 m deep moat, located at the base of the scarps, corresponds with these surfaces and indicates an area of active aeolian scour from descending katabatic winds. Small local dust storms observed along the equator-facing wall of Chasma Boreale imply that slope wind velocities in Chasma Boreale are sufficient to mobilize dust and sand-sized particles in the Polar Layered Deposits (PLD). Two amphitheater forms, located above the cap-forming unit of the sinuous scarp and west of Chasma Boreale, may represent an early stage of polar scarp and chasma formation. These two forms are developing within a younger section of polar layered materials. The unusually steep scarps associated with the polar chasmata have developed where resistant layers are present in the PLD, offering resistance during the headward erosion and poleward retreat of the scarps. Steep slopes that formed under these circumstances enhance the flow of down-scarp katabatic winds. On the basis of these observations, we reject the fluvial outflood hypothesis for the origin of the north polar chasmata and embrace a wind erosion model for their long-term development. In the aeolian model, off-pole katabatic winds progressively remove materials from the steep slopes below chasmata scarps, undermining resistant layers at the tops of scarps and causing retreat by headward erosion. Assuming a minimum age for the onset of formation of Chasma Boreale (10⁵ yr) results in a maximum volumetric erosion rate of . Removal of this volume of material from the equator-facing wall and head scarps of chasma would require a rate for scarp retreat of .     

Light-toned layered deposits in Juventae Chasma, Mars

We examine hypotheses for the formation of light-toned layered deposits in Juventae Chasma using a combination of data from Mars Global Surveyor's Mars Orbiter Camera (MOC), Mars Orbiter Laser Altimeter (MOLA), and Thermal Emission Spectrometer (TES), as well as Mars Odyssey's Thermal Emission Imaging System (THEMIS). We divide Juventae Chasma into geomorphic units of (i) chasm wall rock, (ii) heavily cratered hummocky terrain, (iii) a mobile and largely crater-free sand sheet on the chasm floor, (iv) light-toned layered outcrop (LLO) material, and (v) chaotic terrain. Using surface temperatures derived from THEMIS infrared data and slopes from MOLA, we derive maps of thermal inertia, which are consistent with the geomorphic units that we identify. LLO thermal inertias range from ∼400 to 850 J m⁻² K⁻¹ s^−1/2. Light-toned layered outcrops are distributed over a remarkably wide elevation range () from the chasm floor to the adjacent plateau surface. Geomorphic features, the absence of small craters, and high thermal inertia show that the LLOs are composed of sedimentary rock that is eroding relatively rapidly in the present epoch. We also present evidence for exhumation of LLO material from the west wall of the chasm, within chaotic and hummocky terrains, and within a small depression in the adjacent plateau. The data imply that at least some of the LLO material was deposited long before the adjacent Hesperian plateau basalts, and that Juventae Chasma underwent, and may still be undergoing, enlargement along its west wall due to wall rock collapse, chaotic terrain evolution, and exposure and removal of LLO material. The new data allow us to reassess possible origins of the LLOs. Gypsum, one of the minerals reported elsewhere as found in Juventae Chasma LLO material, forms only at low temperatures () and thus excludes a volcanic origin. Instead, the data are consistent with either multiple occurrences of lacustrine or airfall deposition over an extended period of time prior to emplacement of Hesperian lava flows on the plateau above the chasm.     

The E ring in the vicinity of Enceladus: I. Spatial distribution and properties of the ring particles

Saturn's diffuse E ring is the largest ring of the Solar System and extends from about (Saturn radius RS=60,330 km) to at least encompassing the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea. After Cassini's insertion into her saturnian orbit in July 2004, the spacecraft performed a number of equatorial as well as steep traversals through the E ring inside the orbit of the icy moon Dione. Here, we report about dust impact data we obtained during 2 shallow and 6 steep crossings of the orbit of the dominant ring source—the ice moon Enceladus. Based on impact data of grains exceeding 0.9 μm we conclude that Enceladus feeds a torus populated by grains of at least this size along its orbit. The vertical ring structure at agrees well with a Gaussian with a full-width-half-maximum (FWHM) of ∼4200 km. We show that the FWHM at is due to three-body interactions of dust grains ejected by Enceladus' recently discovered ice volcanoes with the moon during their first orbit. We find that particles with initial speeds between 225 and 235 m s⁻¹ relative to the moon's surface dominate the vertical distribution of dust. Particles with initial velocities exceeding the moon's escape speed of 207 m s⁻¹ but slower than 225 m s⁻¹ re-collide with Enceladus and do not contribute to the ring particle population. We find the peak number density to range between 16×¹⁰⁻² m⁻³ and 21×¹⁰⁻² m⁻³ for grains larger 0.9 μm, and 2.1×¹⁰⁻² m⁻³ and 7.6×¹⁰⁻² m⁻³ for grains larger than 1.6 μm. Our data imply that the densest point is displaced outwards by at least with respect of the Enceladus orbit. This finding provides direct evidence for plume particles dragged outwards by the ambient plasma. The differential size distribution for grains >0.9 μm is described best by a power law with slopes between 4 and 5. We also obtained dust data during ring plane crossings in the vicinity of the orbits of Mimas and Tethys. The vertical distribution of grains >0.8 μm at Mimas orbit is also well described by Gaussian with a FWHM of ∼5400 km and displaced southwards by ∼1200 km with respect to the geometrical equator. The vertical distribution of ring particles in the vicinity of Tethys, however, does not match a Gaussian. We use the FWHM values obtained from the vertical crossings to establish a 2-dimensional model for the ring particle distribution which matches our observations during vertical and equatorial traversals through the E ring.     

Extension of atmospheric dust loading to high altitudes during the 2001 Mars dust storm: MGS TES limb observations

Mars Global Surveyor (MGS) visible (solarband bolometer) and thermal infrared (IR) spectral limb observations from the Thermal Emission Spectrometer (TES) support quantitative profile retrievals for dust opacity and particle sizes during the 2001 global dust event on Mars. The current analysis considers the behavior of dust lifted to altitudes above 30 km during the course of this storm; in terms of dust vertical mixing, particle sizes, and global distribution. TES global maps of visible (solarband) limb brightness at 60 km altitude indicate a global-scale, seasonally evolving (over 190-240° solar longitudes, L_S) longitudinal corridor of vertically extended dust loading (which may be associated with a retrograde propagating, wavenumber 1 Rossby wave). Spherical radiative transfer analysis of selected limb profiles for TES visible and thermal IR radiances provide quantitative vertical profiles of dust opacity, indicating regional conditions of altitude-increasing dust mixing ratios. Observed infrared spectral dependences and visible-to-infrared opacity ratios of dust scattering over 30-60 km altitudes indicate particle sizes characteristic of lower altitudes (cross-section weighted effective radius, ), during conditions of significant dust transport to these altitudes. Conditions of reduced dust loading at 30-60 km altitudes present smaller dust particle sizes . These observations suggest rapid meridional transport at 30-80 km altitudes, with substantial longitudinal variation, of dust lifted to these altitudes over southern hemisphere atmospheric regions characterized by extraordinary (m/s) vertical advection velocities. By L_S=230° dust loading above 50 km altitudes decreased markedly at southern latitudes, with a high altitude (60-80 km) haze of fine (likely) water ice particles appearing over 10°S-40°N latitudes.     

Stability of ice on Mars and the water vapor diurnal cycle: Experimental study of the sublimation of ice through a fine-grained basaltic regolith

In order to advance our understanding of the long-term stability of subsurface ice, the diurnal martian water cycle, and implications for liquid water, we determined diffusion coefficients and adsorption kinetics for the water vapor produced by the sublimation of ice buried beneath various layers of fine-grained (<63, 63-125, and 125-250 μm) basaltic powder under simulated martian conditions. Sublimation rates at shallower depths, <10 mm, were determined to be affected by mass transfer through the atmosphere in addition to the basalt layer. For greater depths, the measured diffusion coefficients for water vapor moving through basalt grains were 1.56±0.53×¹⁰⁻⁴, 2.05±0.82×¹⁰⁻⁴, and for the <63, 63-125, and 125-250 μm basaltic layers, respectively. Through the Brunauer, Emmett and Teller (BET) isotherm, which assumes multiple molecular layers of adsorbed water, we determined the adsorption constants of 52.6±8.3 at 270 K for <63 μm, 39.0±6.4 at 267 K for 63-125 μm, and 54.3±9.3 at 266 K for 125-250 μm, resulting in surface areas of 2.6±0.1×¹⁰⁴, 1.7±0.3×¹⁰⁴, , respectively. These results suggest that while diffusion is too rapid to explain the purported diurnal cycle in water content of the atmosphere, adsorption is efficient and rapid, and does provide an effective mechanism to explain such a cycle. The present diffusion data suggest that very thin, <50 pr μm, shallow, 10 mm, ice deposits would last for >10 h at ∼224 K, just above the freezing point of saturated CaCl₂. Temperatures can remain above ∼224 K over most of the planet, which means that water, even as saturated brine, will sublimate before the freezing point is reached and liquid could be formed. On the other hand, 1 m ice layers below 1 m of fine-grained basaltic regolith at 235 K and 10 Pa of atmospheric water could last 600 to 1300 years. At deeper depths and lower temperatures, ice could last since the last major obliquity change 400,000 years ago.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.     

A sensitive search for nitric oxide in the lower atmospheres of Venus and Mars: Detection on Venus and upper limit for Mars

High-resolution spectra of Venus and Mars at the NO fundamental band at 5.3 μm with resolving power ν/δν=76,000 were acquired using the TEXES spectrograph at NASA IRTF on Mauna Kea, Hawaii. The observed spectrum of Venus covered three NO lines of the P-branch. One of the lines is strongly contaminated, and the other two lines reveal NO in the lower atmosphere at a detection level of 9 sigma. A simple photochemical model for NO and N at 50-112 km was coupled with a radiative transfer code to simulate the observed equivalent widths of the NO and some CO₂ lines. The derived NO mixing ratio is 5.5±1.5 ppb below 60 km and its flux is . Predissociation of NO at the (0-0) 191 nm and (1-0) 183 nm bands of the δ-system and the reaction with N are the only important loss processes for NO in the lower atmosphere of Venus. The photochemical impact of the measured NO abundance is significant and should be taken into account in photochemical modeling of the Venus atmosphere. Lightning is the only known source of NO in the lower atmosphere of Venus, and the detection of NO is a convincing and independent proof of lightning on Venus. The required flux of NO is corrected for the production of NO and N by the cosmic ray ionization and corresponds to the lightning energy deposition of . For a flash energy on Venus similar to that on the Earth (∼¹⁰⁹ J), the global flashing rate is ∼90 s⁻¹ and ∼6 km⁻² y⁻¹ which is in reasonable agreement with the existing optical observations. The observed spectrum of Mars covered three NO lines of the R-branch. Two of these lines are contaminated by CO₂ lines, and the line at 1900.076 cm⁻¹ is clean and shows some excess over the continuum. Some photochemical reactions may result in a significant excitation of NO (v=1) in the lowest 20 km on Mars. However, quenching of NO (v=1) by CO₂ is very effective below 40 km. Excitation of NO (v=1) in the collisions with atomic oxygen is weak because of the low temperature in the martian atmosphere, and we do not see any explanation of a possible emission of NO at 5.3 μm. Therefore the data are treated as the lack of absorption with a 2 sigma upper limit of 1.7 ppb to the NO abundance in the lower atmosphere of Mars. This limit is above the predictions of photochemical models by a factor of 3.     

Heavy ion escape from Mars, influence from solar wind conditions and crustal magnetic fields

We have used more than 4 years of Mars Express ion data to estimate the escape of heavy ions ( and ) from Mars. To take the limited field of view of the instrument into account, the data has been binned into spatial bins and angular bins to create average distribution functions for different positions in the near Mars space. The net escape flux for the studied low solar activity period, between May 2007 and May 2011, is 2.0 ± 0.2 × 10²⁴ s⁻¹. The escape has been calculated independently for four different quadrants in the Y_MSO − Z_MSO plane, south, dusk, north and dawn. Escape is highest from the northern and dusk quadrants, 0.6 ± 0.1 × 10²⁴ s⁻¹, and smallest from the south and dawn quadrants, 0.4 ± 0.1 × 10²⁴ s⁻¹. The flux ratio of molecular ( and ) to O⁺ ions is 0.9 ± 0.1, averaged over all quadrants. The flux difference between the north and south quadrants is statistically significant, and is presumed to be due to the presence of significant crustal magnetic fields in the southern hemisphere, reducing the outflow. The difference between the dawn and dusk quadrants is likely due to the magnetic tension associated with the nominal Parker angle spiral, which should lead to higher average magnetic tension on the dusk side. The escape increases during periods of high solar wind flux and during times when co-rotating interaction regions (CIR) affect Mars. In the latter case the increase is a factor 2.4-2.9 as compared to average conditions.     

Modelling meteor ablation in the venusian atmosphere

A comparative study of meteor ablation in the atmospheres of the Earth and Venus is presented. The classical single body meteor ablation model is extended to incorporate a heat penetration depth estimate allowing the simulation of larger meteoroids, than would an isothermal model. The ablation of icy and rocky meteoroids, with densities of 1.0 and 3.4 g cm⁻³, respectively, and initial radii of up to for rock and for ice (equivalent to an initial mass of in both cases), was simulated in both atmospheres. In general venusian meteors are brighter than terrestrial equivalents. Large, slow, rocky objects may be up to 0.7 mag brighter on Venus, while small, icy particles with entry speeds in the range 30-60 km s⁻¹, are found to be upwards of 2.7 mag brighter than at the Earth. Venusian meteors reach maximum brightness at greater altitudes than would similar particles at the Earth. Rocky meteoroids have their points of maximum brightness some 15-35 km higher up at Venus, between 90 and 120 km, whereas, for icy particles this altitude difference is about 5-25 km higher up than at the Earth, in the range 100-125 km. These findings agree, for the most part, with recent analytical studies. Venusian meteors, which last from 100 ms to , tend to be shorter-lived than terrestrial meteors, with correspondingly shorter visible trails. Large (), slow () icy particles reach a maximum magnitude of ∼−2 at Venus and remain visible for about one second, with a large section of the smaller faster meteoroids simulated here remaining visible for several hundred milliseconds. In light of recent space-based meteor observations at the Earth [Jenniskens, P., Tedesco, E., Muthry, J., Laux, C.O., Price, S., 2002. Meteorit. Planet. Sci. 37, 1071-1078], such brightness, height and duration estimates as suggested in this work, may be used in developing future observational campaigns to be carried out from Venus orbit.     

The global shape of Europa: Constraints on lateral shell thickness variations

The global shape of Europa is controlled by tidal and rotational potentials and possibly by lateral variations in ice shell thickness. We use limb profiles from four Galileo images to determine the best-fit hydrostatic shape, yielding a mean radius of 1560.8±0.3 km and a radius difference a−c of 3.0±0.9 km, consistent with previous determinations and inferences from gravity observations. Adding long-wavelength topography due to proposed lateral variations in shell thickness results in poorer fits to the limb profiles. We conclude that lateral shell thickness variations and long-wavelength isostatically supported topographic variations do not exceed 7 and 0.7 km, respectively. For the range of rheologies investigated (basal viscosities from 10¹⁴ to ) the maximum permissible (conductive) shell thickness is 35 km. The relative uniformity of Europa's shell thickness is due to either a heat flux from the silicate interior, lateral ice flow at the base of the shell, or convection within the shell.     

Parametric analysis of modeled ion escape from Mars

We develop a parametric fit to the results of a detailed magnetohydrodynamic (MHD) study of the response of ion escape rates (O⁺, and ) to strongly varied solar forcing factors, as a way to efficiently extend the MHD results to different conditions. We then use this to develop a second, evolutionary model of solar forced ion escape. We treat the escape fluxes of ion species at Mars as proportional to the product of power laws of four factors - that of the EUV flux R_euv, the solar wind particle density R_ρ, its velocity (squared) R_v², and the interplanetary magnetic field pressure R_B², where forcing factors are expressed in units of the current epoch-averaged values. Our parametric model is: , where ?(i) is the escape flux of ion i. We base our study on the results of just six provided MHD model runs employing large forcing factor variations, and thus construct a successful, first-order parametric model of the MHD program. We perform a five-dimensional least squares fit of this power law model to the MHD results to derive the flux normalizations and the indices of the solar forcing factors. For O⁺, we obtain the values, 1.73 × 10²⁴ s⁻¹, 0.782, 0.251, 0.382, and 0.214, for ?₀, α, β, γ, and δ, respectively. For , the corresponding values are 1.68 × 10²⁴ s⁻¹, −0.393, 0.798, 0.967, and 0.533. For , they are 8.66 × 10²² s⁻¹, −0.427, 1.083, 1.214, and 0.690. The fit reproduces the MHD results to an average error of about 5%, suggesting that the power laws are broadly representative of the MHD model results. Our analysis of the MHD model shows that by itself an increase in R_EUV enhances O⁺ loss, but suppresses the escape of and , whereas increases in solar wind (i.e., in , and R_B², with R_euv constant) favors the escape of heavier ions more than light ions. The ratios of escaping ions detectable at Mars today can be predicted by this parametric fit as a function of the solar forcing factors. We also use the parametric model to compute escape rates over martian history. This second parametric model expresses ion escape functions of one variable (per ion), ?(i) = ?₀(i)(t/t₀)^−ξ(i). The ξ(i) are linear combinations of the epoch-averaged ion escape sensitivities, which are seen to increase with ion mass. We integrate the and oxygen ion escape rates over time, and find that in the last 3.85 Gyr, Mars would have lost about mbars of , and of water (from O⁺ and ) from ion escape.     

The Planetary Fourier Spectrometer (PFS) onboard the European Mars Express mission

The Planetary Fourier Spectrometer (PFS) for the Mars Express mission is an infrared spectrometer optimised for atmospheric studies. This instrument has a short wave (SW) channel that covers the spectral range from 1700 to (1.2-) and a long-wave (LW) channel that covers 250- (5.5-). Both channels have a uniform spectral resolution of . The instrument field of view FOV is about 1.6^° (FWHM) for the Short Wavelength channel (SW) and 2.8^° (FWHM) for the Long Wavelength channel (LW) which corresponds to a spatial resolution of 7 and 12 km when Mars is observed from an height of 250  km. PFS can provide unique data necessary to improve our knowledge not only of the atmosphere properties but also about mineralogical composition of the surface and the surface-atmosphere interaction.The SW channel uses a PbSe detector cooled to 200-220 K while the LW channel is based on a pyroelectric (L_iT_aO₃) detector working at room temperature. The intensity of the interferogram is measured every 150 nm of physical mirrors displacement, corresponding to 600 nm optical path difference, by using a laser diode monochromatic light interferogram (a sine wave), whose zero crossings control the double pendulum motion. PFS works primarily around the pericentre of the orbit, only occasionally observing Mars from large distances. Each measurements take 4 s, with a repetition time of 8.5 s. By working roughly 0.6 h around pericentre, a total of 330 measurements per orbit will be acquired 270 looking at Mars and 60 for calibrations. PFS is able to take measurements at all local times, facilitating the retrieval of surface temperatures and atmospheric vertical temperature profiles on both the day and the night side.     

Experimental study of the sublimation of ice through an unconsolidated clay layer: Implications for the stability of ice on Mars and the possible diurnal variations in atmospheric water

We have studied the sublimation of ice and water vapor transport through various thicknesses of clay (<63 μm grain size). We experimentally demonstrate that both adsorption and diffusion strongly affect the transport of water, and that the processes of diffusion and adsorption can be separately quantified once the system comes to a steady state. At shallow depths of clay, water vapor transport is determined by diffusion through both the atmosphere and the clay layer, whereas at greater depth the rate of sublimation of the ice is governed only by diffusion through the clay. Using two different models, we determine the diffusion coefficient for water vapor through unconsolidated clay layer to be 1.08±0.04×¹⁰⁻⁴ and . We also determined the adsorption isotherms for the clay layer, which follow the Langmuir theory at low water vapor pressure (<100 Pa, where a monolayer of water molecules forms on the surface of the clay) and the BET theory at higher pressure (where multiple water layers form). From our analysis of both types of isotherms we determined the adsorption constants to be and c=30±10, respectively, and specific surface areas of 1.10±0.2×¹⁰⁵ and , respectively. Finally, we report a theoretical kinetic model for the simultaneous diffusion and adsorption from which we determine adsorption kinetic constants according to the Langmuir theory of and . If the martian regolith possesses diffusive properties similar to those of the unconsolidated montmorillonite soil we investigated here, it would not represent a significant barrier to the sublimation of subsurface ice. However, at the low subsurface temperatures of high latitude (180 K on average), ice could survive from the last glaciation period (about 300 to 400,000 years ago). Higher subsurface temperatures in the equatorial regions would prevent long-timescale survival of ice in the shallow subsurface. In agreement with previous work, we show that adsorption of water by a clay regolith could provide a significant reservoir of subsurface water and it might account for the purported diurnal cycle in the water content of the atmosphere.     

New Horizons Alice ultraviolet observations of a stellar occultation by Jupiter’s atmosphere

The Alice ultraviolet spectrograph onboard the New Horizons spacecraft observed two occultations of the bright star χ Ophiucus by Jupiter’s atmosphere on February 22 and 23, 2007 during the approach phase of the Jupiter flyby. The ingress occultation probed the atmosphere at 32°N latitude near the dawn terminator, while egress probed 18°N latitude near the dusk terminator. A detailed analysis of both the ingress and egress occultations, including the effects of molecular hydrogen, methane, acetylene, ethylene, and ethane absorptions in the far ultraviolet (FUV), constrains the eddy diffusion coefficient at the homopause level to be  cm² s⁻¹, consistent with Voyager measurements and other analyses (Festou, M.C., Atreya, S.K., Donahue, T.M., Sandel, B.R., Shemansky, D.E., Broadfoot, A.L. [1981]. J. Geophys. Res. 86, 5717-5725; Vervack Jr., R.J., Sandel, B.R., Gladstone, G.R., McConnell, J.C., Parkinson, C.D. [1995]. Icarus 114, 163-173; Yelle, R.V., Young, L.A., Vervack Jr., R.J., Young, R., Pfister, L., Sandel, B.R. [1996]. J. Geophys. Res. 101 (E1), 2149-2162). However, the actual derived pressure level of the methane homopause for both occultations differs from that derived by [Festou et al., 1981] and [Yelle et al., 1996] from the Voyager ultraviolet occultations, suggesting possible changes in the strength of atmospheric mixing with time. We find that at 32°N latitude, the methane concentration is  cm⁻³ at 70,397 km, the methane concentration is  cm⁻³ at 70,383 km, the acetylene concentration is  cm⁻³ at 70,364 km, and the ethane concentration is  cm⁻³ at 70,360 km. At 18°N latitude, the methane concentration is  cm⁻³ at 71,345 km, the methane concentration is  cm⁻³ at 71,332 km, the acetylene concentration is cm⁻³ at 71,318 km, and the ethane concentration is  cm⁻³ at 71,315 km. We also find that the H₂ occultation light curve is best reproduced if the atmosphere remains cold in the microbar region such that the base of the thermosphere is located at a lower pressure level than that determined by in situ instruments aboard the Galileo probe (Seiff, A., Kirk, D.B., Knight, T.C.D., Young, R.E., Mihalov, J.D., Young, L.A., Milos, F.S., Schubert, G., Blanchard, R.C., Atkinson, D. [1998]. J. Geophys. Res. 103 (E10), 22857-22889) - the Sieff et al. temperature profile leads to too much absorption from H₂ at high altitudes. However, this result is highly model dependent and non-unique. The observations and analysis help constrain photochemical models of Jupiter’s atmosphere.     

High outflow channels on Mars indicate Hesperian recharge at low latitudes and the presence of Canyon Lakes

A number of martian outflow channels were carved by discharges from large dilational fault zones. These channels were sourced by groundwater, not surface water, and when observed on high-standing plateaus they provide indicators of elevated paleo-groundwater levels. We identify three outflow channels of Hesperian age that issued from a 750-km-long fault zone extending from Candor Chasma to Ganges Chasma. Two of these channels, Allegheny Vallis and Walla Walla Vallis, have sources >2500 m above the topographic datum, too high to be explained by discharge from a global aquifer that was recharged solely in the south polar region. The indicated groundwater levels likely required regional sources of recharge at low latitudes. The floodwaters that erupted from Ophir Cavus to form Allegheny Vallis encountered two ridges that restricted the flow, forming temporary lakes. The flow probably breached or overtopped these obstructions quickly, catastrophically draining the lakes and carving several scablands. After the last obstacle had been breached, a single main channel formed that captured all subsequent flow. We performed hydrologic analyses of this intermediate phase of the flooding, prior to incision of the channel to its present depth. Using floodwater depths of 30-60 m, we calculated flow velocities of 6-15 m s⁻¹ and discharges in the range of . Locally higher flow velocities and discharges likely occurred when the transient lakes were drained. Variable erosion at the channel and scabland crossing of MOLA pass 10644 suggests that the upper 25-30 m may consist of poorly consolidated surface materials underlain by more cohesive bedrock. We infer that an ice-covered lake with a surface elevation >2500 m probably existed in eastern Candor Chasma because this canyon is intersected by the Ophir Catenae fault system from which Allegheny Vallis and Walla Walla Vallis originated. We introduce a new hydrology concept for Mars in which the groundwater system was augmented by recharge from canyon lakes that were formed when water was released by catastrophic melting of former ice sheets in Tharsis by effusions of flood basalts. This model could help to reconcile the expected presence of a thick cryosphere during the Hesperian with the abundant evidence for groundwater as a source for some of the circum-Chryse outflow channels.     

Hubble Space Telescope observations of the nucleus fragment 73P/Schwassmann-Wachmann 3-C

The investigation of fragmented comets provides information on the physical properties and internal structure of cometary nuclei, as well as insights into the mechanisms responsible for cometary breakups. The Jupiter-family Comet 73P/Schwassmann-Wachmann 3 (73P/SW3) fragmented non-tidally into at least four components, and probably more, in the autumn of 1995. Fragment C was detected with the Wide Field Planetary Camera 2 (WFPC2) of the Hubble Space Telescope (HST) on 26 November 2001 when it was 3.26 AU from the Sun and 2.34 AU from the Earth. The high spatial resolution of the HST allowed us to separate the signal of the fragment from that of its coma, and to determine its R magnitude in the Johnson-Kron-Cousins photometric system from four images taken with the F675W filter. Assuming a spherical body with a geometric albedo of 0.04 and a linear phase coefficient of 0.04 mag deg⁻¹ for the R band, we derived an effective radius of . The pre-breakup radius of the original nucleus was estimated to be 1.1 km, which implies that the volume of fragment C is ∼25% of the total volume of the pre-breakup nucleus. The limited temporal coverage of our observations preclude deriving an accurate shape or rotational period; our measurements are consistent with a rather spherical body but an elongated shape cannot be excluded. Fragment C was very active despite its rather large heliocentric distance, with an estimated dust production rate of (∼130 metric tons day⁻¹). A very large fraction of the surface area of fragment C must have been sublimating to sustain such a high level of activity. Fragment C may be recovered at its next return in 2006, if it does not experience further fragmentation.     

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).