Polygonal cracks in the seasonal semi-translucent CO2 ice layer in Martian polar areas

In this paper, we use morphological and numerical methods to test the hypothesis that seasonally formed fracture patterns in the Martian polar regions result from the brittle failure of seasonal CO2 slab ice. The observations by the High Resolution Imaging Science Experiment (HiRISE) of polar regions of Mars show very narrow dark elongated linear patterns that are observed during some periods of time in spring, disappear in summer and re-appear again in the following spring. They are repeatedly formed in the same areas but they do not repeat the exact pattern from year to year. This leads to the conclusion that they are cracks formed in the seasonal ice layer. Some of models of seasonal surface processes rely on the existence of a transparent form of CO2 ice, so-called slab ice. For the creation of the observed cracks the ice is required to be a continuous media, not an agglomeration of relatively separate particles like a firn. The best explanation for our observations is a slab ice with relatively high transparency in the visible wavelength range. This transparency allows a solid state green-house effect to act underneath the ice sheet raising the pressure by sublimation from below. The trapped gas creates overpressure and the ice sheet breaks at some point creating the observed cracks. We show that the times when the cracks appear are in agreement with the model calculation, providing one more piece of evidence that CO2 slab ice covers polar areas in spring.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

The Aeolian Flux of Calcium, Chloride and Nitrate to the McMurdo Dry Valleys Landscape: Evidence from Snow Pit Analysis

We have determined the flux of calcium, chloride and nitrate to the McMurdo Dry Valleys region by analysing snow pits for their chemical composition and their snow accumulation using multiple records spanning up to 48 years. The fluxes demonstrate patterns related to elevation and proximity to the ocean. In general, there is a strong relationship between the nitrate flux and snow accumulation, indicating that precipitation rates may have a great influence over the nitrogen concentrations in the soils of the valleys. Aeolian dust transport plays an important role in the deposition of some elements (e.g. C(2+)) into the McMurdo Dry Valleys' soils. Because of the antiquity of some of the soil surfaces in the McMurdo Dry Valleys regions, the accumulated atmospheric flux of salts to the soils has important ecological consequences. Although precipitation may be an important mechanism of salt deposition to the McMurdo Dry Valley surfaces, it is poorly understood because of difficulties in measurement and high losses from sublimation.

The rotation state of 67P/Churyumov-Gerasimenko from approach observations with the OSIRIS cameras on Rosetta

Aims. Approach observations with the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard Rosetta are used to determine the rotation period, the direction of the spin axis, and the state of rotation of comet 67P’s nucleus. Methods. Photometric time series of 67P have been acquired by OSIRIS since the post wake-up commissioning of the payload in March 2014. Fourier analysis and convex shape inversion methods have been applied to the Rosetta data as well to the available ground-based observations. Results. Evidence is found that the rotation rate of 67P has significantly changed near the time of its 2009 perihelion passage, probably due to sublimation-induced torque. We find that the sidereal rotation periods P1 = 12.76129 ± 0.00005 h and P2 = 12.4043 ± 0.0007 h for the apparitions before and after the 2009 perihelion, respectively, provide the best fit to the observations. No signs of multiple periodicity are found in the light curves down to the noise level, which implies that the comet is presently in a simple rotation state around its axis of largest moment of inertia. We derive a prograde rotation model with spin vector J2000 ecliptic coordinates λ = 65° ± 15°, β = + 59° ± 15°, corresponding to equatorial coordinates RA = 22°, Dec = + 76°. However, we find that the mirror solution, also prograde, at λ = 275° ± 15°, β = + 50° ± 15° (or RA = 274°, Dec = + 27°), is also possible at the same confidence level, due to the intrinsic ambiguity of the photometric problem for observations performed close to the ecliptic plane.

Chemical compositions of solid particles present in the Greenland NEEM ice core over the last 110,000 years

This study reports the chemical composition of particles present along Greenland’s North Greenland Eemian Ice Drilling (NEEM) ice core, back to 110,000 years before present. Insoluble and soluble particles larger than 0.45 μm were extracted from the ice core by ice sublimation, and their chemical composition was analyzed using scanning electron microscope and energy dispersive X-ray spectroscopy and micro-Raman spectroscopy. We show that the dominant insoluble components are silicates, whereas NaCl, Na₂SO₄, CaSO ₄, and CaCO₃ represent major soluble salts. For the first time, particles of CaMg(CO₃)₂ and Ca(NO₃)₂ 4H₂O are identified in a Greenland ice core. The chemical speciation of salts varies with past climatic conditions. Whereas the fraction of Na salts (NaCl + Na₂SO₄) exceeds that of Ca salts (CaSO₄+ CaCO₃) during the Holocene (0.6–11.7 kyr B.P.), the two fractions are similar during the Bølling-Allerød period (12.9–14.6 kyr B.P.). During cold climate such as over the Younger Dryas (12.0–12.6 kyr B.P.) and the Last Glacial Maximum (15.0–26.9 kyr B.P.), the fraction of Ca salts exceeds that of Na salts, showing that the most abundant ion generally controls the salt budget in each period. High-resolution analyses reveal changing particle compositions: those in Holocene ice show seasonal changes, and those in LGM ice show a difference between cloudy bands and clear layers, which again can be largely explained by the availability of ionic components in the atmospheric aerosol body of air masses reaching Greenland.

An ice-rich flow origin for the banded terrain in the Hellas basin, Mars

The interior of Hellas Basin displays a complex landscape and a variety of geomorphological domains. One of these domains, the enigmatic banded terrain covers much of the northwestern part of the basin. We use high-resolution (CTX and HiRISE) Digital Terrain Models to show that most of the complex viscous flowing behavior exhibited by the banded terrain is controlled by topography and flow-like interactions between neighboring banded terrain. Furthermore, the interior of the basin hosts several landforms suggestive of the presence of near-surface ice, which include polygonal patterns with elongated pits, scalloped depressions, isolated mounds and collapse structures. We suggest that thermal contraction cracking and sublimation of near-surface ice are responsible for the formation and the development of most of the ice-related landforms documented in Hellas. The relatively pristine form, lack of superposed craters, and strong association with the banded terrain, suggest an Amazonian (<3 Ga) age of formation for these landforms. Finally, relatively high surface pressures (above the triple point of water) expected in Hellas and summer-time temperatures often exceeding the melting point of water ice suggest that the basin may have recorded relatively “temperate” climatic conditions compared to other places on Mars. Therefore, the potentially ice-rich banded terrain may have deformed with lower viscosity and stresses compared to other locations on Mars, which may account for its unique morphology.

Understanding measured water rotational temperatures and column densities in the very innermost coma of Comet 73P/Schwassmann–Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J., 685, 659−677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J., 732) to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann–Wachmann 3 obtained with a very high spatial resolution of ∼30 km using IRCS/Subaru in May 2006 (Bonev, B.P., Mumma, M.J., Kawakita, H., Kobayashi, H., Villanueva, G.L. [2008]. Icarus, 196, 241−248). For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10–20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.

Inventory of the volatiles on comet 67P/Churyumov-Gerasimenko from Rosetta/ROSINA

Context. The ESA Rosetta spacecraft (S/C) is tracking comet 67P/Churyumov-Gerasimenko in close vicinity. This prolonged encounter enables studying the evolution of the volatile coma composition. Aims. Our work aims at comparing the diversity of the coma of 67P/Churyumov-Gerasimenko at large heliocentric distance to study the evolution of the comet during its passage around the Sun and at trying to classify it relative to other comets. Methods. We used the Double Focussing Mass Spectrometer (DFMS) of the ROSINA experiment on ESA’s Rosetta mission to determine relative abundances of major and minor volatile species. This study is restricted to species that have previously been detected elsewhere. Results. We detect almost all species currently known to be present in cometary coma with ROSINA DFMS. As DFMS measured the composition locally, we cannot derive a global abundance, but we compare measurements from the summer and the winter hemisphere with known abundances from other comets. Differences between relative abundances between summer and winter hemispheres are large, which points to a possible evolution of the cometary surface. This comet appears to be very rich in CO2 and ethane. Heavy oxygenated compounds such as ethylene glycol are underabundant at 3 AU, probably due to their high sublimation temperatures, but nevertheless, their presence proves that Kuiper belt comets also contain complex organic molecules.

Regional surface morphology of comet 67P/Churyumov-Gerasimenko from Rosetta/OSIRIS images

Aims. The OSIRIS camera onboard the Rosetta spacecraft has been acquiring images of the comet 67P/Churyumov-Gerasimenko (67P)'s nucleus at spatial resolutions down to similar to 0.17 m/px ever since Aug. 2014. These images have yielded unprecedented insight into the morphological diversity of the comet's surface. This paper presents an overview of the regional morphology of comet 67P. Methods. We used the images that were acquired at orbits similar to 20-30 km from the center of the comet to distinguish different regions on the surface and introduce the basic regional nomenclature adopted by all papers in this Rosetta special feature that address the comet's morphology and surface processes. We used anaglyphs to detect subtle regional and topographical boundaries and images from close orbit (similar to 10 km from the comet's center) to investigate the fine texture of the surface. Results. Nineteen regions have currently been defined on the nucleus based on morphological and/or structural boundaries, and they can be grouped into distinctive region types. Consolidated, fractured regions are the most common region type. Some of these regions enclose smooth units that appear to settle in gravitational sinks or topographically low areas. Both comet lobes have a significant portion of their surface covered by a dusty coating that appears to be recently placed and shows signs of mobilization by aeolian-like processes. The dusty coatings cover most of the regions on the surface but are notably absent from a couple of irregular large depressions that show sharp contacts with their surroundings and talus-like deposits in their interiors, which suggests that short-term explosive activity may play a significant role in shaping the comet's surface in addition to long-term sublimation loss. Finally, the presence of layered brittle units showing signs of mechanical failure predominantly in one of the comet's lobes can indicate a compositional heterogeneity between the two lobes.

Temporal morphological changes in the Imhotep region of comet 67P/Churyumov-Gerasimenko

Aims. We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko in the smooth terrains of the Imhotep region. Methods. We used images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. Results. The morphological changes observed on the surface are visible in the form of roundish features that are growing in size from a given location in a preferential direction at a rate of 5.6-8.1 x 10(-5) m s(-1) during the observational period. The location where the changes started and the contours of the expanding features are bluer than the surroundings, which suggests that ices (H2O and/or CO2) are exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.

Insolation, erosion, and morphology of comet 67P/Churyumov-Gerasimenko

Context. The complex shape of comet 67P and its oblique rotation axis cause pronounced seasonal effects. Irradiation and hence activity vary strongly. Aims. We investigate the insolation of the cometary surface in order to predict the sublimation of water ice. The strongly varying erosion levels are correlated with the topography and morphology of the present cometary surface and its evolution. Methods. The insolation as a function of heliocentric distance and diurnal (spin dependent) variation is calculated using >10(5) facets of a detailed digital terrain model. Shading, but also illumination and thermal radiation by facets in the field of view of a specific facet are iteratively taken into account. We use a two-layer model of a thin porous dust cover above an icy surface to calculate the water sublimation, presuming steady state and a uniform surface. Our second model, which includes the history of warming and cooling due to thermal inertia, is restricted to a much simpler shape model but allows us to test various distributions of active areas. Results. Sublimation from a dirty ice surface yields maximum erosion. A thin dust cover of 50 pm yields similar rates at perihelion. Only about 6% of the surface needs to be active to match the observed water production rates at perihelion. A dust layer of 1 mm thickness suppresses the activity by a factor of 4 to 5. Erosion on the south side can reach more than 10 m per orbit at active spots. The energy input to the concave neck area (Hapi) during northern summer is enhanced by about 50% owing to self-illumination. Here surface temperatures reach maximum values along the foot of the Hathor wall. Integrated over the whole orbit this area receives the least energy input. Based on the detailed shape model, the simulations identify "hot spots" in depressions and larger pits in good correlation with observed dust activity. Three-quarters of the total sublimation is produced while the sub-solar latitude is south, resulting in a distinct dichotomy in activity and morphology. Conclusions. The northern areas display a much rougher morphology than what is seen on Imhotep, an area at the equator that will be fully illuminated when 67P is closer to the Sun. Self-illumination in concave regions enhance the energy input and hence erosion. This explains the early activity observed at Hapi. Cliffs are more prone to erosion than horizontal, often dust covered, areas, which leads to surface planation. Local activity can only persist if the forming cliff walls are eroding. Comet 67P has two lobes and also two distinct sides. Transport of material from the south to the north is probable. The morphology of the Imhotep plain should be typical for the terrains of the yet unseen southern hemisphere.

Morphology and dynamics of the jets of comet 67P/Churyumov-Gerasimenko: Early-phase development

Aims. The OSIRIS camera onboard the Rosetta spacecraft obtained close-up views of the dust coma of comet 67P. The jet structures can be used to trace their source regions and to examine the possible effect of gas-surface interaction. Methods. We analyzed the wide-angle images obtained in the special dust observation sequences between August and September 2014. The jet features detected in different images were compared to study their time variability. The locations of the potential source regions of some of the jets are identified by ray tracing. We used a ring-masking technique to calculate the brightness distribution of dust jets along the projected distance. Results. The jets detected between August and September 2014 mostly originated in the Hapi region. Morphological changes appeared over a timescale of several days in September. The brightness slope of the dust jets is much steeper than the background coma. This might be related to the sublimation or fragmentation of the emitted dust grains. Interaction of the expanding gas flow with the cliff walls on both sides of Hapi could lead to erosion and material down-fall to the nucleus surface.

Pits formation from volatile outgassing on 67P/Churyumov-Gerasimenko

We investigate the thermal evolution of comet 67P/Churyumov-Gerasimenko's subsurface in the Seth_01 region, where active pits have been observed by the ESA/Rosetta mission. Our simulations show that clathrate destabilization and amorphous ice crystallization can occur at depths corresponding to those of the observed pits in a timescale shorter than 67P/Churyumov-Gerasimenko's lifetime in the comet's activity zone in the inner solar system. Sublimation of crystalline ice down to such depths is possible only in the absence of a dust mantle, which requires the presence of dust grains in the matrix small enough to be dragged out by gas from the pores. Our results are consistent with both pits formation via sinkholes or subsequent to outbursts, the dominant process depending on the status of the subsurface porosity. A sealed dust mantle would favor episodic and disruptive outgassing as a result of increasing gas pressure in the pores, while high porosity should allow the formation of large voids in the subsurface due to the continuous escape of volatiles. We finally conclude that the subsurface of 67P/Churyumov-Gerasimenko is not uniform at a spatial scale of similar to 100-200 m.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.