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.

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.

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 "$\eta$ 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 $\eta$ in the range $0.6-0.7$; values $\eta\leq0.5$ or $\eta\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.

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.

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.

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.