Deriving thermophysical properties such as thermal inertia from thermal infrared observations provides useful insights into the structure of the surface material on planetary bodies. The estimation of these properties is usually done by fitting temperature variations calculated by thermophysical models to infrared observations. For multiple free model parameters, traditional methods such as Least-Squares fitting or Markov-Chain Monte-Carlo methods become computationally too expensive. Consequently, the simultaneous estimation of several thermophysical parameters together with their corresponding uncertainties and correlations is often not computationally feasible and the analysis is usually reduced to fitting one or two parameters. Data assimilation methods have been shown to be robust while sufficiently accurate and computationally affordable even for a large number of parameters. This paper will introduce a standard sequential data assimilation method, the Ensemble Square Root Filter, to thermophysical modelling of asteroid surfaces. This method is used to re-analyse infrared observations of the MARA instrument, which measured the diurnal temperature variation of a single boulder on the surface of near-Earth asteroid (162173) Ryugu. The thermal inertia is estimated to be $295 \pm 18$ $\mathrm{J\,m^{-2}\,K^{-1}\,s^{-1/2}}$, while all five free parameters of the initial analysis are varied and estimated simultaneously. Based on this thermal inertia estimate the thermal conductivity of the boulder is estimated to be between 0.07 and 0.12 $\mathrm{W\,m^{-1}\,K^{-1}}$ and the porosity to be between 0.30 and 0.52. For the first time in thermophysical parameter derivation, correlations and uncertainties of all free model parameters are incorporated in the estimation procedure which is more than 5000 times more efficient than a comparable parameter sweep.

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.

After coming to rest on the night side of comet 67P/Churyumov-Gerasimenko, the ROLIS camera on-board Rosetta's Philae lander acquired five images of the surface below the lander, four of which were with the aid of LED illumination of different colors. The images confirm that Philae was perched on a sloped surface. A local horizon is visible in one corner of the image, beyond which we can see the coma. Having spent a full day on the surface Philae was commanded to lift and rotate, after which a final, sixth, LED image was acquired. The change in perspective allowed us to construct a shape model of the surface. The distance to the foreground was about 80 cm, much larger than the nominal 30 cm. This caused stray light, rather than directly reflected LED light, to dominate the image signal, complicating the analysis. The images show a lumpy surface with a roughness of apparently fractal nature. Its appearance is completely different from that of the first landing site, which was characterized by centimeter to meter-sized debris (Mottola et al., 2015). We recognize neither particles nor pores at the image resolution of 0.8 mm per pixel and large color variations are absent. The surface has a bi-modal brightness distribution that can be interpreted in terms of the degree of consolidation, a hypothesis that we support with experimental evidence. We propose the surface below the lander to consist of smooth, cracked plates with unconsolidated edges, similar to terrain seen in CIVA images.