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 ± 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's small lobe.

Aims. Using data from the Rosetta mission to comet 67P/Churyumov–Gerasimenko, we evaluate the physical properties of the surface and subsurface of the nucleus and derive estimates for the thermal inertia (TI) and roughness in several regions on the largest lobe of the nucleus. Methods. We have developed a thermal model to compute the temperature on the surface and in the uppermost subsurface layers of the nucleus. The model takes heat conduction, self-heating, and shadowing effects into account. To reproduce the brightness temperatures measured by the MIRO instrument, the thermal model is coupled to a radiative transfer model to derive the TI. To reproduce the spatially resolved infrared measurements of the VIRTIS instrument, the thermal model is coupled to a radiance model to derive the TI and surface roughness. These methods are applied to Rosetta data from September 2014. Results. The resulting TI values from both instruments are broadly consistent with each other. From the millimetre channel on MIRO, we determine the TI in the subsurface to be <80 JK−1 m−2 s−0.5 for the Seth, Ash, and Aten regions. The submillimetre channel implies similar results but also suggests that higher values could be possible. A low TI is consistent with other MIRO measurements and in situ data from the MUPUS instrument at the final landing site of Philae. The VIRTIS results give a best-fitting value of 80 JK−1 m−2 s−0.5 and values in the range 40–160 JK−1 m−2 s−0.5 in the same areas. These observations also allow the subpixel scale surface roughness to be estimated and compared to images from the OSIRIS camera. The VIRTIS data imply that there is significant roughness on the infrared scale below the resolution of the available shape model and that, counter-intuitively, visually smooth terrain (centimetre scale) can be rough at small (micrometre–millimetre) scales, and visually rough terrain can be smooth at small scales.

Context. With hundreds of exoplanets detected, it is necessary to revisit giant planets accretion models to explain their mass distribution. In particular, formation of sub-jovian planets remains unclear, given the short timescale for the runaway accretion of massive atmospheres. However, gas needs to pass through a circum-planetary disc. If the latter has a low viscosity (as expected if planets form in "dead zones"), it might act as a bottleneck for gas accretion. Aims. We investigate what the minimum accretion rate is for a planet under the limit assumption that the circum-planetary disc is totally inviscid, and the transport of angular momentum occurs solely because of the gravitational perturbations from the star. Methods. To estimate the accretion rate, we present a steady-state model of an inviscid circum-planetary disc, with vertical gas inflow and external torque from the star. Hydrodynamical simulations of a circum-planetary disc were conducted in 2D, in a planetocentric frame, with the star as an external perturber in order to measure the torque exerted by the star on the disc. Results. The disc shows a two-armed spiral wave caused by stellar tides, propagating all the way in from the outer edge of the disc towards the planet. The stellar torque is small and corresponds to a doubling time for a Jupiter mass planet of the order of 5 Myr. Given the limit assumptions, this is clearly a lower bound of the real accretion rate. Conclusions. This result shows that gas accretion onto a giant planet can be regulated by circum-planetary discs. This suggests that the diversity of masses of extra-solar planets may be the result of different viscosities in these discs.

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 (∼ 10 %) outside perihelion compared to the northern ones (< 4 %), and a drastic rise of the effective active fraction of the southern regions (∼ 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 η ≤ 0.5 or η ≥ 0.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.

Aims. Permittivity measurements on porous samples of volcanic origin have been performed in the 0. GHz range under laboratory conditions in support of the Rosetta mission to comet 67P/Churyumov-Gerasimenko, specifically with the MIRO radiometric experiment and CONSERT radar experiment. Methods. The samples were split into several subsamples with different size ranges covering a few µm to 500 µm. Bulk densities of the subsamples were estimated to be in the 800 to 1500 kg/m 3 range. The porosities were in the range of 48% to 65%. From 50 MHz to 6 GHz and at 190 GHz, permittivity has been determined with a coaxial cell and with a quasi-optical bench, respectively. Results. Without taking into account the volume-scattering effect at 190 GHz, the real part of the permittivity, normalized by the bulk density, is in the range of 2.1 to 2.6. The results suggest that the real part of the permittivity of an ice-free dust mantle covering the nucleus is in the 1.5−2.2 range at 190 GHz. From these values, a lower limit for the absorption length for the millimeter receiver of MIRO has been estimated to be between 0.6 and 2 cm, in agreement with results obtained from MIRO in September 2014. At frequencies of interest for CONSERT experiment, the real part of the permittivity of a suspected ice-free dust mantle should be below 2.2. It may be in the range of 1.2 to 1.7 for the nucleus, in agreement with first CONSERT results, taking into account a mean temperature of 110 K and different values for the dust-to-ice volumetric ratio. Estimations of contributions of the different parameters to the permittivity variation may indicate that the porosity is the main parameter.

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 (≳50 J m−2 K−1 s−1∕2) and ice content (≳45% at the equator). In this case, stresses penetrate to a typical depth of ~0.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.

No abstract

We use a laboratory facility to study the sputtering properties of centimeterthick porous water ice subjected to the bombardment of ions and electrons to better understand the formation of exospheres of the icy moons of Jupiter. Our ice samples are as similar as possible to the expected moon surfaces but surface charging of the samples during ion irradiation may distort the experimental results. We therefore monitor the time scales for charging and discharging of the samples when subjected to a beam of ions. These experiments allow us to derive an electric conductivity of deep porous ice layers. The results imply that electron irradiation and sputtering play a non-negligible role for certain plasma conditions at the icy moons of Jupiter. The observed ion sputtering yields from our ice samples are similar to previous experiments where compact ice films were sputtered off a micro-balance.