Comet 49P/Arend-Rigaux is a well known low-activity Jupiter Family comet. Despite the low activity, we have witnessed outgassing activity in 1992, 2004, and 2012. In 2012 a broad tail-like feature (PA$\sim270^\circ, \sim2.3\times10^5$ km) and a narrow jet-like feature (PA$\sim180^\circ, \sim9.3\times10^4$ km) were seen simultaneously. Using Finson-Probstein (FP) dust dynamical models we determine: grain sizes released in each event; duration of activity; when activity peaked; and velocity of the dust particles, allowing us to make comparisons between the events. We find that the tail feature in 2012 is similar to the tail in 1992 with large grains (40-4000 $\mu$m) peaking in activity near perihelion with a long outgassing duration greater than 150 days. The jet feature from 2012, however, is more similar to the 2004 event which we model with small grains (1-8 $\mu$m) with a short duration of activity ($\sim$1 month). The main difference between these two features is that the 2004 event occurs prior to perihelion, while the 2012 event is post-perihelion. We use the grain sizes from the FP models to constrain ice sublimation models. Between 1985 and 2018 we cover 6 apparitions with 26 nights of our own observations plus data from the literature and the Minor Planet Center, which together, allow us to model the heliocentric light curve. We find that the models are consistent with H$_2$O ice sublimation as the volatile responsible for driving activity over most of the active phases and a combination of H$_2$O and CO$_2$ ices are responsible for driving activity near perihelion. We measure the fractional active area over time for H$_2$O and discover that the activity decreases from an average active area of $\sim3\%$ to $\sim0.2\%$. This secular decrease in activity implies that the comet is becoming depleted of volatiles and is in the process of transitioning to a dormant or dead state.

Aims. We investigate the influence of the Yarkovsky force on the long-term orbital evolution of Jupiter Trojan asteroids. Methods. Clones of the observed population with different sizes and different thermal properties were numerically integrated for 1 Gyr with and without the Yarkovsky effect. The escape rate of these objects from the Trojan region as well as changes in the libration amplitude, eccentricity, and inclination were used as a metric of the strength of the Yarkovsky effect on the Trojan orbits. Results. Objects with radii $R\leq$1 km are significantly influenced by the Yarkovsky force. The effect causes a depletion of these objects over timescales of a few hundred million years. As a consequence, we expect the size-frequency distribution of small Trojans to show a shallower slope than that of the currently observable population ($R$ $\gtrsim$ 1 km), with a turning point between $R$ = 100 m and $R$ = 1 km. The effect of the Yarkovsky acceleration on the orbits of Trojans depends on the sense of rotation in a complex way. The libration amplitude of prograde rotators decreases with time while the eccentricity increases. Retrograde rotators experience the opposite effect, which results in retrograde rotators being ejected faster from the 1:1 resonance region. Furthermore, for objects affected by the Yarkovsky force, we find indications that the effect tends to smooth out the differences in the orbital distribution between the two clouds.

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.

Comets are made of volatile and refractory material and naturally experience various degrees of sublimation as they orbit around the Sun. This gas release, accompanied by dust, represents what is traditionally described as activity. Although the basic principles are well established, most details remain elusive, especially regarding the mechanisms by which dust is detached from the surface and subsequently accelerated by the gas flows surrounding the nucleus. During its 2 years rendez-vous with comet 67P/Churyumov-Gerasimenko, ESA's Rosetta has observed cometary activity with unprecedented details, in both the inbound and outbound legs of the comet's orbit. This trove of data provides a solid ground on which new models of activity can be built. In this chapter, we review how activity manifests at close distance from the surface, establish a nomenclature for the different types of observed features, discuss how activity is at the same time transforming and being shaped by the topography, and finally address several potential mechanisms.

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.

(Context) The change of the rotation period and the orientation of the rotation axis of comet 67P/Churyumov-Gerasimenko (67P/C-G) is deducible from images taken by the scientific imaging instruments on-board the Rosetta mission with high precision. Non gravitational forces are a natural explanation for these data. (Aims) We describe observed changes for the orientation of the rotation axis and the rotation period of 67P/C-G. For these changes we give an explanation based on a sublimation model with a best-fit for the surface active fraction (model P). Torque effects of periodically changing gas emissions on the surface are considered. (Methods) We solve the equation of state for the angular momentum in the inertial and the body-fixed frames and provide an analytic theory of the rotation changes in terms of Fourier coefficients, generally applicable to periodically forced rigid body dynamics. (Results) The torque induced changes of the rotation state constrain the physical properties of the surface, the sublimation rate and the local active fraction of the surface. (Conclusions) We determine a distribution of the local surface active fraction in agreement with the rotation properties, period and orientation, of 67P/C-G. The torque movement confirms that the sublimation increases faster than the insolation towards perihelion. The derived relatively uniform activity pattern is discussed in terms of related surface features.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

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.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

Infrared observations of the coma of 67P/Churyumov-Gerasimenko were carried out from July to September 2015, i.e., around perihelion (13 August 2015), with the high-resolution channel of the VIRTIS instrument onboard Rosetta. We present the analysis of fluorescence emission lines of H$_2$O, CO$_2$, $^{13}$CO$_2$, OCS, and CH$_4$ detected in limb sounding with the field of view at 2.7-5 km from the comet centre. Measurements are sampling outgassing from the illuminated southern hemisphere, as revealed by H$_2$O and CO$_2$ raster maps, which show anisotropic distributions, aligned along the projected rotation axis. An abrupt increase of water production is observed six days after perihelion. In the mean time, CO$_2$, CH$_4$, and OCS abundances relative to water increased by a factor of 2 to reach mean values of 32%, 0.47%, and 0.18%, respectively, averaging post-perihelion data. We interpret these changes as resulting from the erosion of volatile-poor surface layers. Sustained dust ablation due to the sublimation of water ice maintained volatile-rich layers near the surface until at least the end of the considered period, as expected for low thermal inertia surface layers. The large abundance measured for CO$_2$ should be representative of the 67P nucleus original composition, and indicates that 67P is a CO$_2$-rich comet. Comparison with abundance ratios measured in the northern hemisphere shows that seasons play an important role in comet outgassing. The low CO$_2$/H$_2$O values measured above the illuminated northern hemisphere are not original, but the result of the devolatilization of the uppermost layers.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

We present THERMAP, a mid-infrared (8-16 {\mu}m) spectro-imager for space missions to small bodies in the inner solar system, developed in the framework of the MarcoPolo-R asteroid sample return mission. THERMAP is very well suited to characterize the surface thermal environment of a NEO and to map its surface composition. The instrument has two channels, one for imaging and one for spectroscopy: it is both a thermal camera with full 2D imaging capabilities and a slit spectrometer. THERMAP takes advantage of the recent technological developments of uncooled microbolometers detectors, sensitive in the mid-infrared spectral range. THERMAP can acquire thermal images (8-18 {\mu}m) of the surface and perform absolute temperature measurements with a precision better than 3.5 K above 200 K. THERMAP can acquire mid-infrared spectra (8-16 {\mu}m) of the surface with a spectral resolution {\Delta}{\lambda} of 0.3 {\mu}m. For surface temperatures above 350 K, spectra have a signal-to-noise ratio >60 in the spectral range 9-13 {\mu}m where most emission features occur.

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. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) 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.

The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

Understanding the concept of habitability is related to an evolutionary knowledge of the particular planet-in-question. Additional indications so-called "systemic aspects" of the planetary system as a whole governs a particular planet's claim on habitability. Here we focus on such systemic aspects and discuss their relevance to the formation of an 'Earth-like' habitable planet. We summarize our results obtained by lunar sample work and numerical models within the framework of the Research Alliance "Planetary Evolution and Life". We consider various scenarios which simulate the dynamical evolution of the Solar System and discuss the likelihood of forming an Earth-like world orbiting another star. Our model approach is constrained by observations of the modern Solar System and the knowledge of its history. Results suggest that the long-term presence of terrestrial planets is jeopardized due to gravitational interactions if giant planets are present. But habitability of inner rocky planets may be supported in those planetary systems hosting giant planets. Gravitational interactions within a complex multiple-body structure including giant planets may supply terrestrial planets with materials which formed in the colder region of the proto-planetary disk. During these processes, water, the prime requisite for habitability, is delivered to the inner system. This may occur either during the main accretion phase of terrestrial planets or via impacts during a post-accretion bombardment. Results for both processes are summarized and discussed with reference to the lunar crater record. Starting from a scenario involving migration of the giant planets this contribution discusses the delivery of water to Earth, the modification of atmospheres by impacts in a planetary system context and the likelihood of the existence of extrasolar Earth-like habitable worlds.