A fully 3-dimensional thermal model of a comet nucleus

A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the internal thermal structure was obtained as function of depth, longitude and hour angle. The effect of the spin axis tilt on the surface temperature distribution was studied in detail. It was found that for small tilt angles, relatively low temperatures may prevail on near-pole areas, despite lateral heat conduction. A high-resolution run for a comet model of 67P/Churyumov-Gerasimenko with low tilt angle, allowing for crystallization of amorphous ice, showed that the amorphous/crystalline ice boundary varies significantly with depth as a function of cometary latitude.     

Tracking the edge of the south seasonal polar cap of Mars

The advance and retreat of the polar caps were one of the first observations that indicated Mars had seasons. Because a large portion of the atmosphere is cycled in and out of the seasonal caps during the year, the frost deposits play a significant role in regional and global atmospheric circulation. Understanding the nature of the seasonal polar caps is imperative if we are to understand the current Martian climate. In this study, we track the seasonal cap edges as a function of season and longitude for the fall and winter seasons (MY27), using data from the Planetary Fourier Spectrometer (PFS) onboard the Mars Express (MEX) ESA mission. Making use of the rapid rise (decrease) in surface temperature that occurs when CO₂ ice is removed (deposited), in a first approach, we defined the advancing cap edge to be where the surface temperature drops below 150 K, and the retreating cap edge where the surface temperature rises above 160 K. In this case, starting from Ls∼50°, the edge progression speed start to be longitude dependent. In the hemisphere that extends form the eastern limit of the Hellas basin to the western limit of the Argyrae basin (and containing the two) the edges progression speed is about a half than that of the other hemisphere; the cap is thus asymmetric and, unexpectedly, no CO₂ ice seems to be present inside the basins. This is because the above mentioned surface temperatures used in this approach to detect the cap edges are not adequate (too low) for the high-pressure regions inside the basins where, following the Clausius–Clapeyron's law, the CO₂ condensation temperature can be several degrees higher than that of the adjacent lower-pressure regions. In the second, final approach, special attention has been given to this aspect and the advancing and retreating cap edges are defined where, respectively, the surface temperatures drop below and rise above the CO₂ condensation temperature for the actual surface pressure values. Now, the results show an opposite situation than the previous one, with the progression speed being higher and the cap more extended (up to −30° latitude) in the hemisphere containing the two major Martian basins. During the fall season, up to Ls∼50° the South Martian polar cap consists of CO₂ frost deposits that advance towards lower latitudes at a constant speed of 10° of latitude per 15 degrees of Ls. The maximum extension (−40° latitude) of the South polar cap occurs somewhere in the 80°–90° Ls range. At the winter solstice, when the edges of the polar night start moving poleward, the cap recession has already started, in response to seasonal changes in insolation. The CO₂ ice South polar cap will recede with a constant speed of ∼5° of latitude every 25° degrees of Ls during the whole winter. The longitudinal asymmetries reduce during the cap retreat and completely disappear around Ls=145°.     

The water cycle in the general circulation model of the martian atmosphere

Within the numerical general-circulation model of the Martian atmosphere MAOAM (Martian Atmosphere: Observation and Modeling), we have developed the water cycle block, which is an essential component of modern general circulation models of the Martian atmosphere. The MAOAM model has a spectral dynamic core and successfully predicts the temperature regime on Mars through the use of physical parameterizations typical of both terrestrial and Martian models. We have achieved stable computation for three Martian years, while maintaining a conservative advection scheme taking into account the water–ice phase transitions, water exchange between the atmosphere and surface, and corrections for the vertical velocities of ice particles due to sedimentation. The studies show a strong dependence of the amount of water that is actively involved in the water cycle on the initial data, model temperatures, and the mechanism of water exchange between the atmosphere and the surface. The general pattern and seasonal asymmetry of the water cycle depends on the size of ice particles, the albedo, and the thermal inertia of the planet’s surface. One of the modeling tasks, which results from a comparison of the model data with those of the TES experiment on board Mars Global Surveyor, is the increase in the total mass of water vapor in the model in the aphelion season and decrease in the mass of water ice clouds at the poles. The surface evaporation scheme, which takes into account the turbulent rise of water vapor, on the one hand, leads to the most complete evaporation of ice from the surface in the summer season in the northern hemisphere and, on the other hand, supersaturates the atmosphere with ice due to the vigorous evaporation, which leads to worse consistency between the amount of the precipitated atmospheric ice and the experimental data. The full evaporation of ice from the surface increases the model sensitivity to the size of the polar cap; therefore, the increase in the latter leads to better results. The use of a more accurate dust scenario changes the model temperatures, which also strongly affects the water cycle.     

Early Mars climate near the Noachian–Hesperian boundary: Independent evidence for cold conditions from basal melting of the south polar ice sheet (Dorsa Argentea Formation) and implications for valley network formation

Currently, and throughout much of the Amazonian, the mean annual surface temperatures of Mars are so cold that basal melting does not occur in ice sheets and glaciers and they are cold-based. The documented evidence for extensive and well-developed eskers (sediment-filled former sub-glacial meltwater channels) in the south circumpolar Dorsa Argentea Formation is an indication that basal melting and wet-based glaciation occurred at the South Pole near the Noachian–Hesperian boundary. We employ glacial accumulation and ice-flow models to distinguish between basal melting from bottom-up heat sources (elevated geothermal fluxes) and top-down induced basal melting (elevated atmospheric temperatures warming the ice). We show that under mean annual south polar atmospheric temperatures (?100 °C) simulated in typical Amazonian climate experiments and typical Noachian–Hesperian geothermal heat fluxes (45–65 mW/m²), south polar ice accumulations remain cold-based. In order to produce significant basal melting with these typical geothermal heat fluxes, the mean annual south polar atmospheric temperatures must be raised from today’s temperature at the surface (?100 °C) to the range of ?50 to ?75 °C. This mean annual polar surface atmospheric temperature range implies lower latitude mean annual temperatures that are likely to be below the melting point of water, and thus does not favor a “warm and wet” early Mars. Seasonal temperatures at lower latitudes, however, could range above the melting point of water, perhaps explaining the concurrent development of valley networks and open basin lakes in these areas. This treatment provides an independent estimate of the polar (and non-polar) surface temperatures near the Noachian–Hesperian boundary of Mars history and implies a cold and relatively dry Mars climate, similar to the Antarctic Dry Valleys, where seasonal melting forms transient streams and permanent ice-covered lakes in an otherwise hyperarid, hypothermal climate.     

A historical search for habitable ice at the Phoenix landing site

A time-resolved energy balance model in the latitude range targeted by Phoenix, and extending back in time over the past 10 Ma, has been developed and used to predict the time-varying temperature field in ground ice over scales ranging from minutes to millions of years. The temperature history is compared to the population doubling times of terrestrial psychrophiles as a function of temperature, and the lifetime of analog microbe spores against de-activation by galactic cosmic rays (GCR), in order to assess the habitability of ground ice and surrounding materials that may be sampled by Phoenix. Metrics are derived to quantify “habitability” and compare different model configurations, including total and maximum continuous time, per year, that ground ice temperatures exceed various thresholds, maximum and average dormancy periods, and maximum and average consecutive growing seasons. The key unknowns in assessing the position, and hence the temperature, of the ground ice table at high northern latitude is the fate of the perennial north polar cap at high obliquity. If enough H₂O ice can persist at polar latitudes to buffer at least the high-latitude atmosphere at all orbital configurations, ground ice is found to be relatively shallow over much of the past 10 Ma, and regularly achieves temperatures in excess of those required for the growth of terrestrial psychrophiles. The dry overburden expected at the landing site can easily be sampled by Phoenix, and includes the “sweet spot” that is characterized by the optimal habitability metrics over the past 10 Ma. If the atmosphere is buffered only by low-latitude ice deposits at obliquities greater than about 30°, the frequency and duration of habitable ice is considerably diminished, and the intervening dormancy periods, during which cosmic ray damage accumulates, are correspondingly longer. In all cases, the maximum dormancy period that must be survived by putative martian psychrophiles is at least an order of magnitude greater than the amount of time required to reduce terrestrial psychrophile spore viability by 10⁻⁶ (∼7×¹⁰⁴ years). Depending on the fate of high-obliquity polar ice, the maximum dormancy period can exceed 4×¹⁰⁶ years, a factor of 60 longer than terrestrial psychrophile spore lifetimes. Habitability of martian ground ice is therefore dependent on putative martian psychrophiles developing robustness against GCR deactivation at least an order of magnitude greater than their terrestrial counterparts. Simulations of ground ice throughout the 65° N-72° N latitude range accessible to Phoenix suggest that higher-latitude ground ice has better habitability metrics, although the discrepancy is less than an order of magnitude for all metrics and across the entire latitude range.     

Water clouds and dust aerosols observations with PFS MEX at Mars

Observations of water ice clouds and dust are among the main scientific goals of the Planetary Fourier Spectrometer (PFS), a payload instrument of the European Mars Express mission. We report some results, obtained in three orbits: 37, 41 and 68. The temperature profile, and dust and water ice cloud opacities are retrieved from the thermal infrared (long-wavelength channel of PFS) in a self-consistent way using the same spectrum. Orographic ice clouds are identified above Olympus (orbit 37) and Ascraeus Mons (orbit 68). Both volcanoes were observed near noon at Ls=337° and 342°, respectively. The effective radius of ice particles is preliminary estimated as 1-3 μm, changing along the flanks. The corresponding visual opacity changes in the interval 0.2-0.4 above Olympus and 0.1-0.6 above Ascraeus Mons. In the case of Ascraeus Mons, the ice clouds were observed mainly above the Southern flank of the volcano with maximum opacity near the summit. In the case of Olympus, the clouds were found above both sides of the top. A different type of ice cloud is observed at latitudes above 50°N (orbit 68) in the polar hood: the effective particle radius is estimated to be 4 μm. Below the 1 mb level an inversion in the temperature profiles is found with maximum temperature at around 0.6 mb. Along orbit 68 it appears above Alba Patera, then it increases to the north and decreases above the CO₂ polar cap. Beginning from latitude 20°S above Tharsis (orbit 68), the ice clouds and dust contribute equally to the spectral shape. Further on, the ice clouds are found everywhere along orbit 68 up to the Northern polar cap, except the areas between the Northern flank of Ascraeus Mons (below 10 km) and the edge of Alba Patera. Orbit 41 is shifted from the orbit 68 by roughly 180° longitude and passes through Hellas. Ice clouds are not visible in this orbit at latitudes below 80°S. The dust opacity is anticorrelated with the surface altitude. From 70°S to 25°N latitude the vertical dust distribution follows an exponential law with a scale height of 11.5±0.5 km, which corresponds to the gaseous scale height near noon and indicates a well-mixed condition. The 9 μm dust opacity, reduced to zero surface altitude, is found to be 0.25±0.05, which corresponds to a visual opacity of 0.5-0.7 (depending on the particle size).     

Martian Seasonal CO2 Ice in Polygonal Troughs in Southern Polar Region: Role of the Distribution of Subsurface H2O Ice

The Mars Orbiter Camera onboard the Mars Global Surveyor has obtained several images of polygonal features in the southern polar region. In images taken during the end of the southern spring, when the surrounding surface is free of the seasonal frost, CO₂ ice still appears to be present within the polygonal troughs. In Earth's polar regions, polygons such as these are indicative of water ice in the ground below. We analyzed the seasonal evolution of the thermal state and the CO₂ content of these features. Our 2-D model includes condensation and sublimation of the CO₂ ice, a self consistent treatment of the variations of the thermal properties of the regolith, and the seasonal variations of the local atmospheric pressure which we take from the results of a general circulation model. We find that the residence time of seasonal CO₂ ice in troughs depends not only on atmospheric opacity and albedo of the CO₂ ice, but also and most significantly on the distribution of water ice in the regolith. Optical properties of the atmosphere and surface CO₂ ice can be independently obtained from observations. To date this is not true about the distribution of water ice below the surface. Our analysis quantifies the dependence of the seasonal cycle of the CO₂ ice within the troughs on the assumed distribution of the water ice below the surface. We show that presence of water ice in the ground at a depth smaller than the depth of the troughs reduces winter condensation rate of CO₂ ice. This is due to higher heat flux conducted from the water ice rich regolith toward the facets of the troughs.     

Surface temperatures and retention of H2O frost on Ganymede and Callisto

Surface temperatures and ice evaporation rates are calculated for Ganymede and Callisto as a function of latitude, time of day, and albedo. The model uses surface thermal properties determined by eclipse radiometry (Morrison and Cruikshank, 1973Icarus18 224–236) and albedos determined from photometrically decalibrated Voyager images. Daytime temperatures on Callisto are roughly 8°K warmer than those in Ganymede's cratered terrain and 11°K warmer than those in Ganymede's grooved terrain. Diurnal mean ice evaporation rates are high enough on both bodies that the surface material probably consists of a very low density lag deposit of primarily silicate dust overlying a denser regolith of silicates and ice. The difference in temperature between Ganymede and Callisto is not great enough to account for the lack of bright polar caps on Callisto. This lack seems instead to reflect a real deficiency in the amount of available H₂O frost relative to Ganymede. The temperature difference between Ganymede's grooved and cratered terrains also cannot account for the strong concentration of bright ray craters in grooved terrain. This concentration suggests instead that an internal geologic process has enriched the grooved terrain in ice relative to the cratered terrain.     

Spatial inhomogeneity of the martian subsurface water distribution: implication from a global water cycle model

In an effort to test and to understand the global hydrogen distribution in the shallow subsurface of Mars retrieved by the Mars Odyssey gamma-ray spectrometer, the present state and movement of water are investigated by a coupled global subsurface-atmosphere water cycle model. It was found that the observed global subsurface hydrogen distribution is largely consistent with the modeled global water cycle, so a large fraction of hydrogen is likely to exist as water, at low and mid latitudes in the form of adsorbed water. Under the present climate the water content in the shallow subsurface becomes higher in the northern hemisphere than in the southern hemisphere as a result of global water cycle, regardless of the initial water distribution in the soil or adsorptive capacity. The higher annual maximum soil temperature in the south, stronger net northward transport of atmospheric water vapor, and the emission of vapor from the northern residual polar cap in northern summer contribute to this hemispheric asymmetry. The generally higher adsorptive capacity of clay minerals in the northern plains may further increase this bias. The longitudinal inhomogeneity is caused by several factors, such as thermal inertia, adsorptive capacity, and atmospheric surface pressure. The water abundance is locally high in low thermal inertia regions (e.g., Arabia Terra) and at deep places where the surface pressure is high (e.g., Hellas); it is low in soil with a low adsorptive capacity (e.g., Tharsis) and high thermal inertia regions (e.g., Solis Planum). Most of the soil humidity near the surface at low and mid latitudes may originate from the atmosphere. The model implies that the upper soil layer should be largely ice-free because otherwise an excessive sublimation and vapor emission into the atmosphere in warm seasons would violate the observational constraints. Moreover, the more uniform latitudinal variation of the observed hydrogen abundance near the surface compared to that of deeper layers is indicative of the presence of adsorbed water instead of ground ice because the adsorbed water content does not as steeply depend on latitude as the ground ice stability. Concerning the regolith mineralogy, montmorillonite can much better account for the observed water cycle than palagonite. While the presence of permanent ground ice appears likely in the polar region below a thin layer, large seasonal cycle of phase change between pore ice and adsorbed water may be possible. Regolith adsorption/desorption is neither negligible nor crucial for the seasonal atmospheric water cycle, but the surface-atmosphere coupling is a major prerequisite for the long-term evolution of subsurface water distribution.     

Polar warming in the middle atmosphere of Mars

We report ground-based laser heterodyne spectroscopy of non-thermal emission in the cores of the 10.33-μmR(8) and 10.72-μmP(32) lines of ¹²C¹⁶O₂, obtained at 23 locations on the disk of Mars during the 1984 opposition, at L_s = 130°. The data were obtained at a sub-Doppler spectral resolution, and the temperature of the middle Martian atmosphere (50–85 km) is derived from the frequency width and intensity of the R(8) emission, and from the total intensity of the P(32) emission. We find that the temperature of the middle Martian atmosphere varies with latitude. Near the subsolar latitude, the average 50- to 85-km temperature is close to the radiative equilibrium value for a CO₂ atmosphere. However, at high latitudes in both the northern (summer) and southern (winter) hemispheres the 50- to 85-km temperature exceeds the CO₂ radiative equilibrium value; a meridional gradient in the range of 0.4 – 0.9°K per degree of latitude is indicated by our data. The highest temperatures are seen at high latitudes in the winter hemisphere, reminiscent of the seasonal effects seen at the Earth's mesopause. As in the terrestrial case, this winter polar warming in the Martian middle atmosphere necessitates departures from radiative equilibrium; dynamical heating of order 4 × 10² ergs g⁻¹ sec⁻¹ is required at the edge of the winter polar night. A comparison with 2-D circulation models shows that the presence of atmospheric dust may enhance this dynamical heating at high winter latitudes, and may also account for heating at high latitudes in the summer hemisphere.     

Numerical modeling of endogenic thermal anomalies on Europa

A variety of recent resurfacing features have been observed on Europa, which may produce thermal anomalies detectable by a future mission. However, the likelihood of such a detection depends on their size and lifetimes. The results of this numerical study suggest that the lifetime of a thermal anomaly associated with the emplacement of 100 m of water onto the surface of Europa is several hundred years, and ∼10 years for 10 m of water. If warm ice is emplaced on the surface instead of liquid water, these lifetimes decrease by up to a factor of two. Exploration of model parameters indicates that a thin insulating surface layer can double thermal anomaly lifetimes, anomalies emplaced at a latitude of 80° can remain detectable nearly a factor of two longer than those at equatorial latitudes, and anomalies on the night side can remain detectable for up to ∼20% longer than those on the day side. High temperatures are very short-lived as the surface ice cools very rapidly to below 200 K due to sublimation cooling. Assuming steady-state resurfacing, the number of detectable thermal anomalies associated with the emplacement of 100 m of water would be on the order of 10 if the typical resurfacing area is 15 km². If recent resurfacing is dominated by chaos regions with typical areas of 100 to 1000 km² and lifetimes of 1000 to 4000 years, the number of detectable thermal anomalies would be on the order of 1 to 10.     

Global distribution and migration of subsurface ice on mars

A thermal/diffusive model of H₂O kinetics and equilibrium was developed to investigate the long-term evolution and depth distribution of subsurface ice on Mars. The model quantitatively takes into account (1) obliquity variations; (2) eccentricity variations; (3) long-term changes in the solar luminosity; (4) variations in the argument of subsolar meridian (in planetocentric equatorial coordinates); (5) albedo changes at higher latitudes due to seasonal phase changes of CO₂ and the varying extent of CO₂ ice cover; (6) planetary internal heat flow; (7) temperature variations in the regolith as a function of depth, time, and latitude due to the above factors; (8) atmospheric pressure variations over a 10⁴-year time scale; (9) the effects of factors (1) through (5) on seasonal polar cap temperatures; and (10) Knudsen and molecular diffusion of H₂O through the regolith. The migration of H₂O into or out of the regolith is determined by two boundary conditions, the H₂O vapor pressure at the subsurface ice boundary and the annual average H₂O concentration at the base of the atmosphere. These are controlled respectively by the annual average regolith temperature at the given depth and seasonal temperatures at the polar cap. Starting from an arbitrary initial uniform depth distribution of subsurface ice, H₂O fluxes into or out of the regolith are calculated for 100 selected obliquity cycles, each representing a different epoch in Mars' history. The H₂O fluxes are translated into ice thicknesses and extrapolated over time to give the subsurface ice depth as a function of latitude and time. The results show that obliquity variations influence annual average regolith temperatures in varying degrees, depending on latitude, with the greatest effect at the poles and almost no effect at 40° lat. Insolation changes at the pole, due to obliquity, argument of subsolar meridian, and eccentricity variations can produce enormous atmospheric H₂O concentration variations of ≈6 orders of magnitude over an obliquity cycle. Superimposed on these cyclic variations is a slow, monotonic change due to the increasing solar luminosity. Albedo changes at the polar cap due to seasonal phase changes of CO₂ and the varying thickness of the CO₂ ice cover are critically important in determining annual average atmospheric H₂O concentrations. Despite the strongly oscillating character of the boundary conditions, only small amounts of H₂O are exchanged between the regolith and the atmosphere per obliquity cycle (<10 g/cm²). The net result of H₂O migration is that the regolith below 30–40° lat is depleted of subsurface ice, while the regolith above 30–40° lat contains permanent ice due to the depth of penetration of the annual thermal wave. This result is supported by recent morphological studies. The rate of migration of H₂O is strongly dependent on average pore/capillary radius for which we have assumed values of 1 and 10 μm. We estimate that the H₂O ice removed from the regolith would produce a permanent ice cap with a volume between 2 × 10⁶ and 6 × 10⁶ km³. This generally agrees with estimates deduced from deflationary features at lower latitudes, depositional features at higher latitudes, and the mass of the polar caps.     

MGS TES observations of the water vapor above the seasonal and perennial ice caps during northern spring and summer

We report on new retrievals of water vapor column abundances from the Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) data. The new retrievals are from the TES nadir data taken above the ‘cold’ surface areas in the North polar region (T_surf < 220 K, including seasonal frost and permanent ice cap) during spring and summer seasons, where retrievals were not performed initially. Retrievals are possible (with some modifications to the original algorithm) over cold surfaces overlaid by sufficiently warm atmosphere. The retrieved water vapor column abundances are compared to the column abundances observed by other spacecrafts in the Northern polar region during spring and summer and good agreement is found. We detect an annulus of water vapor growing above the edge of the retreating seasonal cap during spring. The formation of the vapor annulus is consistent with the previously proposed mechanism for water cycling in the polar region, according to which vapor released by frost sublimation during spring re-condenses on the retreating seasonal CO₂ cap. The source of the vapor in the vapor annulus, according to this model, is the water frost on the surface of the CO₂ at the retreating edge of the cap and the frost on the ground that is exposed by the retreating cap. Small contribution from regolith sources is possible too, but cannot be quantified based on the TES vapor data alone. Water vapor annulus exhibits interannual variability, which we attribute to variations in the atmospheric temperature. We propose that during spring and summer the water ice sublimation is retarded by high relative humidity of the local atmosphere, and that higher atmospheric temperatures lead to higher vapor column abundances by increasing the water holding capacity of the atmosphere. Since the atmospheric temperatures are strongly influenced by the atmospheric dust content, local dust storms may be controlling the release of vapor into the polar atmosphere. Water vapor abundances above the residual polar cap also exhibit noticeable interannual variability. In some years abundances above the cap are lower than the abundances outside of the cap, consistent with previous observations, while in the other years the abundances above the cap are higher or similar to abundances outside of the cap. We speculate that the differences may be due to weaker off-cap transport in the latter case, keeping more vapor closer to the source at the surface of the residual cap. Despite the large observed variability in water vapor column abundances in the Northern polar region during spring and summer, the latitudinal distribution of the vapor mass in the atmosphere is very similar during the summer season. If the variability in vapor abundances is caused by the variability of vapor sources across the residual cap then this would mean that they annually contribute relatively little vapor mass to significantly affect the vapor mass budget. Alternatively this may suggest that the vapor variability is caused by the variability of the polar atmospheric circulation. The new water vapor retrievals should be useful in tuning the Global Circulation Models of the martian water cycle.     

Mars and Earth: Comparison of cold-climate features

On Earth, glacial and periglacial features are common in areas of cold climate. On Mars, the temperature of the present-day surface is appropriate for permafrost, and the presence of water is suspected from data relating to the outgassing of the planet, from remote-sensing measurements over the polar caps and elsewhere on the Martian surface, and from recognition of fluvial morphological features such as channels. These observations and the possibility that ice could be in equilibrium with the atmosphere in the high latitudes north and south of ±40° latitude suggest that glacial and periglacial features should exist on the planet. Morphological studies based mainly on Viking pictures indicate many features that can be attributed to the action of ice. Among these features are extensive talus aprons; debris avalanches; flows that resemble glaciers or rock glaciers; ridges that look like moraines; various types of patterned ground, scalloped scarps, and chaotically collapsed terrain that could be attributed to thermokarst processes; and landforms that may reflect the interaction of volcanism and ice.     

Radiative habitable zones in martian polar environments

The biologically damaging solar ultraviolet (UV) radiation (quantified by the DNA-weighted dose) reaches the martian surface in extremely high levels. Searching for potentially habitable UV-protected environments on Mars, we considered the polar ice caps that consist of a seasonally varying CO2 ice cover and a permanent H2O ice layer. It was found that, though the CO2 ice is insufficient by itself to screen the UV radiation, at approximately 1 m depth within the perennial H2O ice the DNA-weighted dose is reduced to terrestrial levels. This depth depends strongly on the optical properties of the H2O ice layers (for instance snow-like layers). The Earth-like DNA-weighted dose and Photosynthetically Active Radiation (PAR) requirements were used to define the upper and lower limits of the northern and southern polar Radiative Habitable Zone (RHZ) for which a temporal and spatial mapping was performed. Based on these studies we conclude that photosynthetic life might be possible within the ice layers of the polar regions. The thickness varies along each martian polar spring and summer between approximately 1.5 and 2.4 m for H2O ice-like layers, and a few centimeters for snow-like covers. These martian Earth-like radiative habitable environments may be primary targets for future martian astrobiological missions. Special attention should be paid to planetary protection, since the polar RHZ may also be subject to terrestrial contamination by probes.     

Results from the Greenland Search for Meteorites expedition

Abstract— Following discoveries of blue ice areas in Greenland resembling meteorite‐bearing blue ice fields in Antarctica, a surface search of several of the most promising sites was carried out in August 2003. The ice fields are located in Kong Christian X Land, in northeastern Greenland around 74°N at elevations between 2100 and 2400 m. No meteorites were found in any of the localities that were searched. Evidence of occasional significant melting (filled crevasses and melt sheets) suggest that summer temperatures are sometimes high enough that dark rocks, like meteorites, can melt through the upper layers of ice. Small terrestrial rocks and cryogenite were found down to 50 cm below the ice surface. Meter‐sized terrestrial rocks were found on top of the ice downstream from nunataks. These rocks shade the ice below, and since they were apparently too massive to warm up during warm days, they remained at the surface as the surrounding ice ablated away. Our findings strongly suggest that Greenland is currently unlikely to harbor significant meteorite concentrations on blue ice fields.     

An idealized short-period comet model: Surface insolation,H2O flux,dust flux,and mantle evolution

A model of cometary activity is developed which integrates the feedback processes involving heat, gas, and dust transport, and dust mantle development. The model includes the effects of latitude, rotation, and spin axis orientation. Results are obtained for various grain size distributions, dust-to-ice ratios, and spin axis orientations. Attention is focused on the development, change of structure and distribution of dust mantles and their mutual interaction with ice surface temperature and gas and dust production. In this model the dust mantle controls the mechanism of gas transport not onlu by its effect on the temperature but, more importantly, by its own dynamic stability. Results suggest that an initially homogeneous short-period comet with a “cosmic” dust-to-water ice ratio, typical orbit, rotation rate, and grain size distribution would develop at most only a thin (<1 mm) cyclic mantle at all points on the nucleus. Such a fully developed temporary mantle would exist throughout the diurnal cycle only beyond ～4AU. Thus, cyclic behavior would be expected for such an idealized comet, at least for most of its lifetime. Long-term irreversible mantle development on comets with typical rotation rates was not found except regionally on Encke and also on objects with perihelia ?1.5 AU. Even in these cases, free silicate exists, after a few cycles, only as relatively rare large grains and agglomerates with radii ～1 cm scattered over exposed ice. Full mantle development would require hundreds to thousands of cycles. In the case of an initially homogeneous comet Encke, this slow incipient mantle development is shown to be the direct result of its peculiar axial orientation. High obliquity appears required for long-term mantle development for typical rotation rates and perihelia ?1.5 AU. Heat conduction into the nucleus for an incompletely mantled or bald comet has been found to be very important in maintaining relatively higher ice surface temperatures, and hence fluxes, during those portions of the diurnal and orbital cycles which would otherwise be cooler. It is also shown to be at least one cause of post perihelion brightness asymmetries, especially in lower obliquity comets. Maximum heliocentric distances at which 1-μm dust, sand, pebbles, cobbles, and boulders can be permanently ejected from the subsolar point by H₂O (CO₂) are (in AU): 6.9 (16.8), 5.2 (11.5), 1.8 (3.0), 0.21 (0.34) and 0.07 (0.11), respectively. A detailed anatomy of temperature, gas and dust fluxes vs latitude and longitude for a homogeneous rotating comet with fixed axis is given for comparison with future observations. Most H₂O flux histories deduced from brightness data are found to be in reasonable agreement with the model, allowing for uncertainty in radius and albedo. A clear exception is Encke. It is shown that the large discrepancy between Encke's observed and model predicted fluxes, based on radar cross section, can be used to evaluate the extent of exposed ice (<10%). The model is then used to place an active area so as to explain a reported sharp drop in flux on approach to the Sun at 0.78 AU. An active area or areas, <10% of the comet's surface, centered near 65°N latitude appears indicated. Although cyclic mantles are generally indicated for the set of parameters we used, our results show that a global mantle only 1 to 3 cm thick (depending on the orbit) consisting of a full range of grain sizes can cause irresversible evolution to a noncometary body. We investigated the long-term evolution of such a postulated initially thinly mantled cometary object. It was found that after the first few passes and until the end of its dynamic lifetime the object averaged <3 × 10^?12 g cm^?1 sec^?1 H₂O flux. Therefore, if cometary objects evolve into Apollo asteroids, ice should always be accessible within 10 m of the surface despite numerous close perihelion passages. The possible impact of factors not included in the model, such as initial inhomogeneities, coma scattering of radiation, and global redistribution of ejected silicate around the nucleus, are discussed.     

Endogenic heat from Enceladus' south polar fractures: New observations, and models of conductive surface heating

Linear features dubbed “tiger stripes” in the south polar region of Enceladus have anomalously high heat fluxes and are the apparent source of the observed plume. Several explanations for the observed activity have been proposed, including venting from a subsurface reservoir of liquid water, sublimation of surface ice, dissociation of clathrates, and shear heating. Thermal modeling presented in this work, coupled with observations from the Cassini Composite Infrared Spectrometer (CIRS) instrument, seeks to elucidate the underlying physical mechanism by constraining vent temperatures and thermal emission sources, using a model in which the observed thermal signature results primarily from conductive heating of the surface by warm subsurface fractures. The fractures feed surface vents, which may themselves contribute to the observed thermal emission. Model variables include vent temperature, presence of a surface insulating layer, vent width, time-variable heat input, and heat sources other than the central vent. Results indicate that CIRS spectra are best fitted with a model in which the surface is heated by narrow vents at temperatures as high as 223 K. Although equally good fits can be obtained for vent temperatures in the range of 130 to 155 K if the vents are wider (180 m and 22 m respectively) and dominate the emission spectrum, these models are probably less realistic because vents with these temperatures and widths cannot supply the observed H₂O vapor flux. The lack of emission angle dependence of the thermal emission when July 2005 and November 2006 CIRS observations are compared also argues against thermal emission being dominated by the vents themselves. Thus, results favor high-temperature models, possibly venting from a subsurface liquid water reservoir. However, a fracture filled with liquid water near the surface would produce significantly higher radiances than were detected unless masked by a thermally insulating surface layer. Models that best match the CIRS data are characterized by small fractions of the surface at high temperatures, which strengthens the case for the vents and/or their conductively-heated margins being the primary heat source. Models where the thermal emission is dominated by conductive heating of the surface from below by a laterally-extensive buried heat source cannot reproduce the observed spectrum. Models with a 10 cm thick upper insulating layer produce a poor match to the CIRS spectra, suggesting high thermal inertias near the tiger stripes. Finally, tiger stripe thermal emission measured by CIRS varied by less than 15% over the 16 month period from July 2005 to November 2006.     

Evaporation of ice in planetary atmospheres: Ice-covered rivers on mars

The evaporation rate of water ice on the surface of a planet with an atmosphere involves an equilibrium between solar heating and radiative and evaporative cooling of the ice layer. The thickness of the ice is governed principally by the solar flux which penetrates the ice layer and then is conducted back to the surface. These calculations differ from those of Lingenfelter et al. [(1968) Science161, 266–269] for putative lunar channels in including the effect of the atmosphere. Evaporation from the surface is governed by two physical phenomena: wind and free convection. In the former case, water vapor diffuses from the surface of the ice through a lamonar boundary layer and then is carried away by eddy diffusion above, provided by the wind. The latter case, in the absence of wind, is similar, except that the eddy diffusion is caused by the lower density of water vapor than the Martian atmosphere. For mean Martian insolations the evaporation rate above the ice is ～ 10^?8 g cm^?2 sec^?1. Thus, even under present Martian conditions a flowing channel of liquid water will be covered with ice which evaporates sufficiently slowly that the water below can flow for hundreds of kilometers even with quite modest discharges. Evaporation rates are calculated for a wide range of frictional velocities, atmospheric pressures, and insolations and it seems clear that at least some subset of observed Martian channels may have formed as ice-choked rivers. Typical equilibrium thicknesses of such ice covers are ～ 10 to 30 m; typical surface temperatures are 210 to 235°K. Ice-covered channels or lakes on Mars today may be of substantial biological interest. Ice is a sufficiently poor conductor of heat that sunlight which penetrates it can cause melting to a depth of several meters or more. Because the obliquity of Mars can vary up to some 35°, the increased polar heating at such times seems able to cause subsurface melting of the ice caps to a depth which corresponds to the observed lamina thickness and may be responsible for the morphology of these polar features.     

The effects of reduced land fraction and solar forcing on the general circulation: results from the NCAR CCM

Land fraction and the solar energy at the top of the atmosphere (solar constant) may have been significantly lower early in Earth's history. It is likely that both of these factors played some important role in the climate of the early earth. The climate changes associated with a global ocean(i.e. no continents) and reduced solar constant are examined with a general circulation model and compared with the present-day climate simulation. The general circulation model used in the study is the NCAR CCM with a swamp ocean surface. First, all land points are removed in the model and then the solar constant is reduced by 10% for this global ocean case.Results indicate that a 4 K increase in air temperature occurs with global ocean simulation compared to the control. When solar constant is reduced by 10% under global ocean conditions a 23 K decrease in air temperature is noted. The global ocean warms much of the troposphere and stratosphere, while a reduction in the solar constant cools the troposphere and stratosphere. The largest cooling occurs near the surface with the lower solar constant.Global mean values of evaporation, water vapor amounts, absorbed solar radiation and the downward longwave radiation are increased under global ocean conditions, while all are reduced when the solar constant is lowered. The global ocean simulation produces sea ice only in the highest latitudes. A frozen planet does not occur when the solar constant is reduced—rather, the ice line settles near 30° of latitude. It is near this latitude that transient eddies transport large amounts of sensible heat across the ice line acting as a negative feedback under lower solar constant conditions keeping sea ice from migrating to even lower latitudes.Clouds, under lower solar forcing, also act as a negative feedback because they are reduced in higher latitudes with colder atmospheric temperatures allowing additional solar radiation to reach the surface. The overall effect of clouds in the global ocean is to act as a positive feedback because they are slightly reduced thereby allowing additional solar radiation to reach the surface and increase the warming caused by the removal of land. The relevance of the results to the “Faint-Young Sun Paradox” indicates that reduced land fraction and solar forcing affect dynamics, heat transport, and clouds. Therefore the associated feedbacks should be taken into account in order to understand their roles in resolving the “Faint-Young Sun Paradox”.