-
A quantitative description of comet 67P's dust and gas production remains enigmatic
Authors:
Dorothea Bischoff,
Christian Schuckart,
Nicholas Attree,
Bastian Gundlach,
Jürgen Blum
Abstract:
The mechanism of dust emission from a cometary nucleus is still an open question and thermophysical models have problems reproducing outgassing and dust productions rates simultaneously. In this study, we investigate the capabilities of a rather simple thermophysical model to match observations from Rosetta instruments at comet 67P/Churyumov-Gerasimenko and the influence of model variations. We as…
▽ More
The mechanism of dust emission from a cometary nucleus is still an open question and thermophysical models have problems reproducing outgassing and dust productions rates simultaneously. In this study, we investigate the capabilities of a rather simple thermophysical model to match observations from Rosetta instruments at comet 67P/Churyumov-Gerasimenko and the influence of model variations. We assume a macro-porous surface structure composed of pebbles and investigate the influence of different model assumptions. Besides the scenario in which dust layers are ejected when the vapour pressure overcomes the tensile strength, we use artificial ejection mechanisms, depending on ice-depletion of layers. We find that dust activity following the pressure criterion is only possible for reduced tensile strength values or reduced gas diffusivity and is inconsistent with observed outgassing rates, because activity is driven by CO$_2$. Only when we assume that dust activity is triggered when the layer is completely depleted in H$_2$O, the ratio of CO$_2$ to H$_2$O outgassing rates is in the expected order of magnitude. However, the dust-to-H$_2$O ratio is never reproduced. Only with decreased gas diffusivity, the slope of the H$_2$O outgassing rate is matched, however absolute values are too low. To investigate maximum reachable pressures, we adapted our model equivalent to a gas-impermeable dust structure. Here, pressures exceeding the tensile strength by orders of magnitude are possible. Maximum activity distances of $3.1 \,\mathrm{au}$, $8.2 \,\mathrm{au}$, and $74 \,\mathrm{au}$ were estimated for H$_2$O-, CO$_2$-, and CO-driven activity of $1 \,\mathrm{cm}$-sized dust, respectively. In conclusion, the mechanism behind dust emission remains unclear.
△ Less
Submitted 12 June, 2023;
originally announced June 2023.
-
Activity distribution of comet 67P/Churyumov-Gerasimenko from combined measurements of non-gravitational forces and torques
Authors:
Nicholas Attree,
Laurent Jorda,
Olivier Groussin,
Jessica Agarwal,
Riccardo Lasagni Manghi,
Paolo Tortora,
Marco Zannoni,
Raphael Marschall
Abstract:
Aims. Understanding the activity is vital for deciphering the structure, formation, and evolution of comets. We investigate models of cometary activity by comparing them to the dynamics of 67P/Churyumov-Gerasimenko. Methods. We matched simple thermal models of water activity to the combined Rosetta datasets by fitting to the total outgassing rate and four components of the outgassing induced non-g…
▽ More
Aims. Understanding the activity is vital for deciphering the structure, formation, and evolution of comets. We investigate models of cometary activity by comparing them to the dynamics of 67P/Churyumov-Gerasimenko. Methods. We matched simple thermal models of water activity to the combined Rosetta datasets by fitting to the total outgassing rate and four components of the outgassing induced non-gravitational force and torque, with a final manual adjustment of the model parameters to additionally match the other two torque components. We parametrised the thermal model in terms of a distribution of relative activity over the surface of the comet, and attempted to link this to different terrain types. We also tested a more advanced thermal model based on a pebble structure. Results. We confirm a hemispherical dichotomy and non-linear water outgassing response to insolation. The southern hemisphere of the comet and consolidated terrain show enhanced activity relative to the northern hemisphere and dust-covered, unconsolidated terrain types, especially at perihelion. We further find that the non-gravitational torque is especially sensitive to the activity distribution, and to fit the pole-axis orientation in particular, activity must be concentrated (in excess of the already high activity in the southern hemisphere and consolidated terrain) around the south pole and on the body and neck of the comet over its head. This is the case for both the simple thermal model and the pebble-based model. Overall, our results show that water activity cannot be matched by a simple model of sublimating surface ice driven by the insolation alone, regardless of the surface distribution, and that both local spatial and temporal variations are needed to fit the data.
△ Less
Submitted 12 January, 2023;
originally announced January 2023.
-
Gas flow in Martian spider formation
Authors:
Nicholas Attree,
Erkia Kaufmann,
Axel Hagermann
Abstract:
Martian araneiform terrain, located in the Southern polar regions, consists of features with central pits and radial troughs which are thought to be associated with the solid state greenhouse effect under a CO$_{2}$ ice sheet. Sublimation at the base of this ice leads to gas buildup, fracturing of the ice and the flow of gas and entrained regolith out of vents and onto the surface. There are two p…
▽ More
Martian araneiform terrain, located in the Southern polar regions, consists of features with central pits and radial troughs which are thought to be associated with the solid state greenhouse effect under a CO$_{2}$ ice sheet. Sublimation at the base of this ice leads to gas buildup, fracturing of the ice and the flow of gas and entrained regolith out of vents and onto the surface. There are two possible pathways for the gas: through the gap between the ice slab and the underlying regolith, as proposed by Kieffer et al (2007), or through the pores of a permeable regolith layer, which would imply that regolith properties can control the spacing between adjacent spiders, as suggested by Hao et al. We test this hypothesis quantitatively in order to place constraints on the regolith properties. Based on previously estimated flow rates and thermophysical arguments, we suggest that there is insufficient depth of porous regolith to support the full gas flow through the regolith. By contrast, free gas flow through a regolith--ice gap is capable of supplying the likely flow rates for gap sizes on the order of a centimetre. This size of gap can be opened in the centre of a spider feature by gas pressure bending the overlying ice slab upwards, or by levitating it entirely as suggested in the original Kieffer et al (2007) model. Our calculations therefore support at least some of the gas flowing through a gap opened between the regolith and ice. Regolith properties most likely still play a role in the evolution of spider morphology, by regolith cohesion controlling the erosion of the central pit and troughs, for example.
△ Less
Submitted 29 January, 2021;
originally announced January 2021.
-
Potential Effects of Atmospheric Collapse on Martian Heat Flow and Application to the InSight Measurements
Authors:
Nicholas Attree,
Narissa Patel,
Axel Hagermann,
Matthias Grott,
Tilman Spohn,
Matthew Siegler
Abstract:
Heat flow is an important constraint on planetary formation and evolution. It has been suggested that Martian obliquity cycles might cause periodic collapses in atmospheric pressure, leading to corresponding decreases in regolith thermal conductivity (which is controlled by gas in the pore spaces). Geothermal heat would then build up in the subsurface, potentially affecting present-day heat flow -…
▽ More
Heat flow is an important constraint on planetary formation and evolution. It has been suggested that Martian obliquity cycles might cause periodic collapses in atmospheric pressure, leading to corresponding decreases in regolith thermal conductivity (which is controlled by gas in the pore spaces). Geothermal heat would then build up in the subsurface, potentially affecting present-day heat flow - and thus the measurements made by a heat-flow probe such as the InSight HP$^{3}$ instrument. To gauge the order of magnitude of this effect, we model the diffusion of a putative heat pulse caused by thermal conductivity changes with a simple numerical scheme and compare it to the heat-flow perturbations caused by other effects. We find that an atmospheric collapse to 300 Pa in the last 40 kyr would lead to a present-day heat flow that is up to $2-8\%$ larger than the average geothermal background. Considering the InSight mission with expected $5-15\%$ error bars on the HP$^{3}$ measurement, this perturbation would only be significant in the best-case scenario of full instrument deployment, completed measurement campaign, and a well-modelled surface configuration. The prospects for detecting long-term climate perturbations via spacecraft heat-flow experiments remain challenging.
△ Less
Submitted 10 October, 2019;
originally announced October 2019.
-
The thermal, mechanical, structural, and dielectric properties of cometary nuclei after Rosetta
Authors:
O. Groussin,
N. Attree,
Y. Brouet,
V. Ciarletti,
B. Davidsson,
G. Filacchione,
H. H. Fischer,
B. Gundlach,
M. Knapmeyer,
J. Knollenberg,
R. Kokotanekova,
E. Kührt,
C. Leyrat,
D. Marshall,
I. Pelivan,
Y. Skorov,
C. Snodgrass,
T. Spohn,
F. Tosi
Abstract:
The physical properties of cometary nuclei observed today relate to their complex history and help to constrain their formation and evolution. In this article, we review some of the main physical properties of cometary nuclei and focus in particular on the thermal, mechanical, structural and dielectric properties, emphasizing the progress made during the Rosetta mission. Comets have a low density…
▽ More
The physical properties of cometary nuclei observed today relate to their complex history and help to constrain their formation and evolution. In this article, we review some of the main physical properties of cometary nuclei and focus in particular on the thermal, mechanical, structural and dielectric properties, emphasizing the progress made during the Rosetta mission. Comets have a low density of 480 $\pm$ 220 kg m-3 and a low permittivity of 1.9 - 2.0, consistent with a high porosity of 70 - 80 %, are weak with a very low global tensile strength $<$100 Pa, and have a low bulk thermal inertia of 0 - 60 J K-1 m-2 s-1/2 that allowed them to preserve highly volatiles species (e.g. CO, CO2, CH4, N2) into their interior since their formation. As revealed by 67P/Churyumov-Gerasimenko, the above physical properties vary across the nucleus, spatially at its surface but also with depth. The broad picture is that the bulk of the nucleus consists of a weakly bonded, rather homogeneous material that preserved primordial properties under a thin shell of processed material, and possibly covered by a granular material; this cover might in places reach a thickness of several meters. The properties of the top layer (the first meter) are not representative of that of the bulk nucleus. More globally, strong nucleus heterogeneities at a scale of a few meters are ruled out on 67P small lobe.
△ Less
Submitted 3 May, 2019;
originally announced May 2019.
-
Constraining models of activity on comet 67P/Churyumov-Gerasimenko with Rosetta trajectory, rotation, and water production measurements
Authors:
Nicholas Attree,
Laurent Jorda,
Olivier Groussin,
Stefano Mottola,
Nick Thomas,
Yann Brouet,
Ekkehard Kührt,
Martin Knapmeyer,
Frank Preusker,
Frank Scholten,
Jorg Knollenberg,
Stubbe Hviid,
Paul Hartogh,
Rafael Rodrigo
Abstract:
Aims. We use four observational data sets, mainly from the Rosetta mission, to constrain the activity pattern of the nucleus of comet 67P/Churyumov-Gerasimenko. Methods. We develop a numerical model that computes the production rate and non-gravitational acceleration of the nucleus of comet 67P as a function of time, taking into account its complex shape with a shape model reconstructed from OSIRI…
▽ More
Aims. We use four observational data sets, mainly from the Rosetta mission, to constrain the activity pattern of the nucleus of comet 67P/Churyumov-Gerasimenko. Methods. We develop a numerical model that computes the production rate and non-gravitational acceleration of the nucleus of comet 67P as a function of time, taking into account its complex shape with a shape model reconstructed from OSIRIS imagery. We use this model to fit three observational data sets: the trajectory data from flight dynamics; the rotation state, as reconstructed from OSIRIS imagery; and the water production measurements from ROSINA, of 67P. The two key parameters of our model, adjusted to fit the three data sets all together, are the activity pattern and the momentum transfer efficiency (i.e., the so-called "$η$ parameter" of the non-gravitational forces). Results. We find an activity pattern able to successfully reproduce the three data sets simultaneously. The fitted activity pattern exhibits two main features: a higher effective active fraction in two southern super-regions ($\sim 10$~\%) outside perihelion compared to the northern ones ($< 4$~\%), and a drastic rise of the effective active fraction of the southern regions ($\sim 25-35$~\%) around perihelion. We interpret the time-varying southern effective active fraction by cyclic formation and removal of a dust mantle in these regions. Our analysis supports moderate values of the momentum transfer coefficient $η$ in the range $0.6-0.7$; values $η\leq0.5$ or $η\geq0.8$ degrade significantly the fit to the three data sets. Our conclusions reinforce the idea that seasonal effects linked to the orientation of the spin axis play a key role in the formation and evolution of dust mantles, and in turn largely control the temporal variations of the gas flux.
△ Less
Submitted 9 January, 2019;
originally announced January 2019.
-
Tensile Strength of 67P/Churyumov-Gerasimenko Nucleus Material from Overhangs
Authors:
N. Attree,
O. Groussin,
L. Jorda,
D. Nébouy,
N. Thomas,
Y. Brouet,
E. Kührt,
F. Preusker,
F. Scholten,
J. Knollenberg,
P. Hartogh,
H. Sierks,
C. Barbieri,
P. Lamy,
R. Rodrigo,
D. Koschny,
H. Rickman,
H. U. Keller,
M. F. A'Hearn,
A. -T. Auger,
M. A. Barucci,
J. -L. Bertaux,
I. Bertini,
D. Bodewits,
S. Boudreault
, et al. (30 additional authors not shown)
Abstract:
We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of mos…
▽ More
We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.
△ Less
Submitted 20 December, 2017;
originally announced December 2017.
-
Thermal fracturing on comets. Applications to 67P/Churyumov-Gerasimenko
Authors:
N. Attree,
O. Groussin,
L. Jorda,
S. Rodionov,
A-T. Auger,
N. Thomas,
Y. Brouet,
O. Poch,
E. Kührt,
M. Knapmeyer,
F. Preusker,
F. Scholten,
J. Knollenberg,
S. Hviid,
P. Hartogh
Abstract:
We simulate the stresses induced by temperature changes in a putative hard layer near the surface of comet 67P/Churyumov--Gerasimenko with a thermo-viscoelastic model. Such a layer could be formed by the recondensation or sintering of water ice (and dust grains), as suggested by laboratory experiments and computer simulations, and would explain the high compressive strength encountered by experime…
▽ More
We simulate the stresses induced by temperature changes in a putative hard layer near the surface of comet 67P/Churyumov--Gerasimenko with a thermo-viscoelastic model. Such a layer could be formed by the recondensation or sintering of water ice (and dust grains), as suggested by laboratory experiments and computer simulations, and would explain the high compressive strength encountered by experiments on board the Philae lander. Changes in temperature from seasonal insolation variation penetrate into the comet's surface to depths controlled by the thermal inertia, causing the material to expand and contract. Modelling this with a Maxwellian viscoelastic response on a spherical nucleus, we show that a hard, icy layer with similar properties to Martian permafrost will experience high stresses: up to tens of MPa, which exceed its material strength (a few MPa), down to depths of centimetres to a metre. The stress distribution with latitude is confirmed qualitatively when taking into account the comet's complex shape but neglecting thermal inertia. Stress is found to be comparable to the material strength everywhere for sufficient thermal inertia ($\gtrsim50$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) and ice content ($\gtrsim 45\%$ at the equator). In this case, stresses penetrate to a typical depth of $\sim0.25$ m, consistent with the detection of metre-scale thermal contraction crack polygons all over the comet. Thermal fracturing may be an important erosion process on cometary surfaces which breaks down material and weakens cliffs.
△ Less
Submitted 27 November, 2017;
originally announced November 2017.
-
A Survey of Low-Velocity Collisional Features in Saturn's F Ring
Authors:
Nicholas O. Attree,
Carl D. Murray,
Gareth A. Williams,
Nicholas J. Cooper
Abstract:
Small (~50km scale), irregular features seen in Cassini images to be emanating from Saturn's F ring have been termed mini-jets by Attree et al. (2012). One particular mini-jet was tracked over half an orbital period, revealing its evolution with time and suggesting a collision with a local moonlet as its origin. In addition to these data we present here a much more detailed analysis of the full ca…
▽ More
Small (~50km scale), irregular features seen in Cassini images to be emanating from Saturn's F ring have been termed mini-jets by Attree et al. (2012). One particular mini-jet was tracked over half an orbital period, revealing its evolution with time and suggesting a collision with a local moonlet as its origin. In addition to these data we present here a much more detailed analysis of the full catalogue of over 800 F ring mini-jets, examining their distribution, morphology and lifetimes in order to place constraints on the underlying moonlet population. We find mini-jets randomly located in longitude around the ring, with little correlation to the moon Prometheus, and randomly distributed in time, over the full Cassini tour to date. They have a tendency to cluster together, forming complicated `multiple' structures, and have typical lifetimes of ~1d. Repeated observations of some features show significant evolution, including the creation of new mini-jets, implying repeated collisions by the same object. This suggests a population of <~1km radius objects with some internal strength and orbits spread over 100km in semi-major axis relative to the F ring but with the majority within 20km. These objects likely formed in the ring under, and were subsequently scattered onto differing orbits by, the perturbing action of Prometheus. This reinforces the idea of the F ring as a region with a complex balance between collisions, disruption and accretion.
△ Less
Submitted 12 September, 2013;
originally announced September 2013.