On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

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.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough information to clearly isolate contribution from Hapi and Inhotep independently, and compare it to the nominal case of activity from the entire illuminated surface. We also demonstrate that spectral line shapes differ from the 1D model for different viewing geometries and coma conditions relevant to this study. Specifically, line shapes are somewhat sensitive to the location of the terminator in the coma. At last, fitting the MIRO observations we show that the Imhotep region (possible extended source of H$_{2}$O due to CO$_{2}$ activity) contributes only a small fraction of the total number of water molecules into MIRO beam in the early activity. On the other hand, a strong enhancement of water activity from the Hapi region seems required to fit the MIRO line shapes. This is consistent with earlier Rosetta results. Nevertheless, within the assumption of our coma and surface boundary conditions, we cannot get a reasonable fit to all MIRO mapping observations in July 2014, which may illustrate that a more sophisticated coma model or more accurate temperature/velocity distribution is needed.

In April 2016, the Rosetta spacecraft performed a low-altitude low-phase-angle flyby over the Imhotep-Khepry transition of 67P/Churyumov-Gerasimenko's nucleus. The OSIRIS/Narrow-Angle-Camera (NAC) acquired 112 images with mainly 3 broadband filters in the visible at a resolution of up to 0.53 m/px and for phase angles between 0.095{\deg} and 62{\deg}. Using those images, we have investigated the morphological and spectrophotometrical properties of this area. We assembled the images into coregistered color cubes. Using a 3D shape model, we produced the illumination conditions and georeference for each image. We projected the observations on a map to investigate its geomorphology. Observations were photometrically corrected using the Lommel-Seeliger disk law. Spectrophotometric analyses were performed on the coregistered color cubes. These data were used to estimate the local phase reddening. This region of the nucleus hosts numerous and varied types of terrains and features. We observe an association between a feature's nature, its reflectance, and its spectral slope. Fine material deposits exhibit an average reflectance and spectral slope, while terrains with diamictons, consolidated material, degraded outcrops, or features such as somber boulders, present a lower-than-average reflectance and higher-than-average spectral slope. Bright surfaces present here a spectral behavior consistent with terrains enriched in water-ice. We find a phase-reddening slope of 0.064{\pm}0.001{\%}/100nm/{\deg} at 2.7 au outbound, similarly to the one obtained at 2.3 au inbound during the February 2015 flyby. Identified as the source region of multiple jets and a host of water-ice material, the Imhotep-Khepry transition appeared in April 2016, close to the frost line, to further harbor several potential locations with exposed water-ice material among its numerous different morphological terrain units.

(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.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

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.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

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 Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

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.

The Descent Imager/Spectral Radiometer aboard the Huygens probe successfully acquired images and spectra of the surface of Titan. To counter the effects of haze and atmospheric methane absorption it carried a Surface Science Lamp to illuminate the surface just before landing. We reconstruct the reflectance spectrum of the landing site in the 500-1500 nm range from Downward Looking Visual and Infrared Spectrometer data that show evidence of lamp light. Our reconstruction is a followup to the analysis by Tomasko et al. (2005), who scaled their result to the ratio of the up- and down flux measured just before landing. Instead, we use the lamp flux from the calibration experiment, and find a significantly higher overall reflectance. We attribute this to a phase angle dependance, possibly representing the opposition surge commonly encountered on solar system bodies. The reconstruction in the visible wavelength range is greatly improved. Here, the reflectance spectrum features a red slope, consistent with the presence of organic material. We confirm the blue slope in the near-IR, featureless apart from a single shallow absorption feature at 1500 nm. We agree with Tomasko et al. that the evidence for water ice is inconclusive. By modeling of absorption bands we find a methane mixing ratio of 4.5% +/- 0.5% just above the surface. There is no evidence for the presence of liquid methane, but the data do not rule out a wet soil at a depth of several centimeters.

During its year-long orbital mission, the Dawn spacecraft has mapped the surface of main-belt asteroid Vesta multiple times at different spatial resolutions and illumination and viewing angles. The onboard Framing Camera has acquired thousands of clear filter and narrow band images, which, with the availability of high-resolution global shape models, allows for a photometric characterization of the surface in unprecedented detail. We analyze clear filter images to retrieve the photometric properties of the regolith. In the first part of the paper we evaluate different photometric models for the global average. In the second part we use these results to study variations in albedo and steepness of the phase curve over the surface. Maps of these two photometric parameters show large scale albedo variations, which appear to be associated with compositional differences. They also reveal the location of photometrically extreme terrains, where the phase curve is unusually shallow or steep. We find that shallow phase curves are associated with steep slopes on crater walls and faults, as calculated from a shape model. On the other hand, the phase curve of ejecta associated with young impact craters is steep. We interpret these variations in phase curve slope in terms of physical roughness of the regolith. The lack of rough ejecta around older craters suggests that initially rough ejecta associated with impact craters on Vesta are smoothed over a relatively short time of several tens of Myr. We propose that this process is the result of impact gardening, and as such represents a previously unrecognized aspect of Vesta space weathering (Pieters et al., 2012). If this type of space weathering is common, we may expect to encounter this photometric phenomenon on other main belt asteroids.

The OSIRIS camera onboard Rosetta successfully acquired images of asteroid 2867 Steins through a variety of color filters during the flyby on 5 September 2008. The best images of this 5 km diameter asteroid have a resolution of 78 meters per pixel. We process the images by deconvolving with the point spread function and enlarging through the Mitchell-Netravali filter. The enhanced set is analyzed by means of various techniques (PCA, band ratios, stereo anaglyphs) to study surface morphology and search for variegation. We identify a landslide, which supports a YORP origin for Steins' unusual diamond shape. In addition, we find that the interior of one of two large craters on the south pole is bluer than the rest of the body.

We present a method for calibrating images acquired by the Dawn Framing Camera (FC) that is based on the results of an in-flight calibration campaign performed during the cruise from Earth to Vesta. We describe this campaign and the data analysis in full. Both the primary camera FC2 and the backup camera FC1 are radiometrically and geometrically calibrated through observations of standard stars, star fields, and solar system objects. The calibration in each spectral filter is accurate to within a few percent for point sources. Geometric distortion, small by design, is characterized with high accuracy. Dark current, monitored on a regular basis, is very low at flight operational temperatures. Out-of-field stray light was characterized using the Sun as a stray light source. In-field stray light is confirmed in narrow-band filter images of Vesta. Its magnitude and distribution are scene-dependent, and expected to contribute significantly to images of extended objects. Description of a method for in-field stray light correction is deferred to a follow-up paper, as is a discussion of the closely related topic of flat-fielding.

The visible reflectance spectrum of many solar system bodies changes with changing viewing geometry for reasons not fully understood. It is often observed to redden (increasing spectral slope) with increasing solar phase angle, an effect known as phase reddening. Only once, in an observation of the Martian surface by the Viking 1 lander, was reddening observed up to a certain phase angle with bluing beyond, making the reflectance ratio as a function of phase angle shaped like an arch. However, in laboratory experiments this arch-shape is frequently encountered. To investigate this, we measured the bidirectional reflectance of particulate samples of several common rock types in the 400-1000 nm wavelength range and performed ray-tracing simulations. We confirm the occurrence of the arch for surfaces that are forward scattering, i.e. are composed of semi-transparent particles and are smooth on the scale of the particles, and for which the reflectance increases from the lower to the higher wavelength in the reflectance ratio. The arch shape is reproduced by the simulations, which assume a smooth surface. However, surface roughness on the scale of the particles, such as the Hapke and van Horn (1963) fairy castles that can spontaneously form when sprinkling a fine powder, leads to monotonic reddening. A further consequence of this form of microscopic roughness (being indistinct without the use of a microscope) is a flattening of the disk function at visible wavelengths, i.e. Lommel-Seeliger-type scattering. The experiments further reveal monotonic reddening for reflectance ratios at near-IR wavelengths. The simulations fail to reproduce this particular reddening, and we suspect that it results from roughness on the surface of the particles. Given that the regolith of atmosphereless solar system bodies is composed of small particles, our results indicate that the prevalence of monotonic (...)

During the flyby in 2010, the OSIRIS camera on-board Rosetta acquired hundreds of high-resolution images of asteroid Lutetia's surface through a range of narrow-band filters. While Lutetia appears very bland in the visible wavelength range, Magrin et al. (2012) tentatively identified UV color variations in the Baetica cluster, a group of relatively young craters close to the north pole. As Lutetia remains a poorly understood asteroid, such color variations may provide clues to the nature of its surface. We take the color analysis one step further. First we orthorectify the images using a shape model and improved camera pointing, then apply a variety of techniques (photometric correction, principal component analysis) to the resulting color cubes. We characterize variegation in the Baetica crater cluster at high spatial resolution, identifying crater rays and small, fresh impact craters. We argue that at least some of the color variation is due to space weathering, which makes Lutetia's regolith redder and brighter.

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.

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.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

Since the orbital insertion of the Rosetta spacecraft, comet 67P/Churyumov-Gerasimenko (67P/C-G) has been mapped by OSIRIS camera and VIRTIS spectro-imager, producing a huge quantity of images and spectra of the comet's nucleus. The aim of this work is to search for the presence of H$_2$O on the nucleus which, in general, appears very dark and rich in dehydrated organic material. After selecting images of the bright spots which could be good candidates to search for H$_2$O ice, taken at high resolution by OSIRIS, we check for spectral cubes of the selected coordinates to identify these spots observed by VIRTIS. The selected OSIRIS images were processed with the OSIRIS standard pipeline and corrected for the illumination conditions for each pixel using the Lommel-Seeliger disk law. The spots with higher I/F were selected and then analysed spectrophotometrically and compared with the surrounding area. We selected 13 spots as good targets to be analysed by VIRTIS to search for the 2 micron absorption band of water ice in the VIRTIS spectral cubes. Out of the 13 selected bright spots, eight of them present positive H$_2$O ice detection on the VIRTIS data. A spectral analysis was performed and the approximate temperature of each spot was computed. The H$_2$O ice content was confirmed by modeling the spectra with mixing (areal and intimate) of H$_2$O ice and dark terrain, using Hapke's radiative transfer modeling. We also present a detailed analysis of the detected spots.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

We present OSIRIS/NAC observations of decimetre-sized, likely ice-containing aggregates ejected from a confined region on the surface of comet 67P/Churyumov-Gerasimenko. The images were obtained in January 2016 when the comet was at 2 AU from the Sun out-bound from perihelion. We measure the acceleration of individual aggregates through a two-hour image series. Approximately 50% of the aggregates are accelerated away from the nucleus, and 50% towards it, and likewise towards either horizontal direction. The accelerations are up to one order of magnitude stronger than local gravity, and are most simply explained by the combined effect of gas drag accelerating all aggregates upwards, and the recoil force from asymmetric outgassing, either from rotating aggregates with randomly oriented spin axes and sufficient thermal inertia to shift the temperature maximum away from an aggregate's subsolar region, or from aggregates with variable ice content. At least 10% of the aggregates will escape the gravity field of the nucleus and feed the comet's debris trail, while others may fall back to the surface and contribute to the deposits covering parts of the northern hemisphere. The rocket force plays a crucial role in pushing these aggregates back towards the surface. Our observations show the future back fall material in the process of ejection, and provide the first direct measurement of the acceleration of aggregates in the innermost coma (<2km) of a comet, where gas drag is still significant.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

The ESA Rosetta spacecraft, currently orbiting around comet 67P, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. A Monte Carlo dust tail code has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is $\alpha$=--3, instead of $\alpha$=--2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 au, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 au. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

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 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.

We use the gravitational instability formation scenario of cometesimals to derive the aggregate size that can be released by the gas pressure from the nucleus of comet 67P/Churyumov-Gerasimenko for different heliocentric distances and different volatile ices. To derive the ejected aggregate sizes, we developed a model based on the assumption that the entire heat absorbed by the surface is consumed by the sublimation process of one volatile species. The calculations were performed for the three most prominent volatile materials in comets, namely, H_20 ice, CO_2 ice, and CO ice. We find that the size range of the dust aggregates able to escape from the nucleus into space widens when the comet approaches the Sun and narrows with increasing heliocentric distance, because the tensile strength of the aggregates decreases with increasing aggregate size. The activity of CO ice in comparison to H_20 ice is capable to detach aggregates smaller by approximately one order of magnitude from the surface. As a result of the higher sublimation rate of CO ice, larger aggregates are additionally able to escape from the gravity field of the nucleus. Our model can explain the large grains (ranging from 2 cm to 1 m in radius) in the inner coma of comet 67P/Churyumov-Gerasimenko that have been observed by the OSIRIS camera at heliocentric distances between 3.4 AU and 3.7 AU. Furthermore, the model predicts the release of decimeter-sized aggregates (trail particles) close to the heliocentric distance at which the gas-driven dust activity vanishes. However, the gas-driven dust activity cannot explain the presence of particles smaller than ~1 mm in the coma because the high tensile strength required to detach these particles from the surface cannot be provided by evaporation of volatile ices. These smaller particles can be produced for instance by spin-up and centrifugal mass loss of ejected larger aggregates.

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.

One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 67P/Churyumov-Gerasimenko at close range. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general and in the framework of current solar system formation models. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the early epochs of the solar system. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution.

In this work, a novel approach to explain the survival of sungrazing comets within the Roche limit is presented. It is shown that in the case of low tensile strength of the cometary nucleus, tidal splitting of the nucleus can be prevented by the reaction force caused by the sublimation of the icy constituents. The survival of Comet C/2011 W3 (Lovejoy) within the Roche limit of the Sun is, thus, the result of high tensile strength of the nucleus, or the result of the reaction force caused by the strong outgassing of the icy constituents near the Sun.

The gas transport through non-volatile random porous media is investigated numerically. We extend our previous research of the transport of molecules inside the uppermost layer of a cometary surface (Skorov and Rickmann, 1995; Skorov et al. 2001). We assess the validity of the simplified capillary model and its assumptions to simulate the gas flux trough the porous dust mantle as it has been applied in cometary physics. A new microphysical computational model for molecular transport in random porous media formed by packed spheres is presented. The main transport characteristics such as the mean free path distribution and the permeability are calculated for a wide range of model parameters and compared with those obtained by more idealized models. The focus in this comparison is on limitations inherent in the capillary model. Finally a practical way is suggested to adjust the algebraic Clausing formula taking into consideration the nonlinear dependence of permeability on layer porosity. The retrieved dependence allows us to accurately calculate the permeability of layers whose thickness and porosity vary in the range of values expected for the near-surface regions of a cometary nucleus.

The cratering history of main belt asteroid (2867) Steins has been investigated using OSIRIS imagery acquired during the Rosetta flyby that took place on the 5th of September 2008. For this purpose, we applied current models describing the formation and evolution of main belt asteroids, that provide the rate and velocity distributions of impactors. These models coupled with appropriate crater scaling laws, allow the cratering history to be estimated. Hence, we derive Steins' cratering retention age, namely the time lapsed since its formation or global surface reset. We also investigate the influence of various factors -like bulk structure and crater erasing- on the estimated age, which spans from a few hundred Myrs to more than 1Gyr, depending on the adopted scaling law and asteroid physical parameters. Moreover, a marked lack of craters smaller than about 0.6km has been found and interpreted as a result of a peculiar evolution of Steins cratering record, possibly related either to the formation of the 2.1km wide impact crater near the south pole or to YORP reshaping.

We compute the distribution of velocities of the particles ejected by the impact of the projectile released from NASA Deep Impact spacecraft on the nucleus of comet 9P/Tempel 1 on the successive 20 hours following the collision. This is performed by the development and use of an ill-conditioned inverse problem approach, whose main ingredients are a set of observations taken by the Narrow Angle Camera (NAC) of OSIRIS onboard the Rosetta spacecraft, and a set of simple models of the expansion of the dust ejecta plume for different velocities. Terminal velocities are derived using a maximum likelihood estimator. We compare our results with published estimates of the expansion velocity of the dust cloud. Our approach and models reproduce well the velocity distribution of the ejected particles. We consider these successful comparisons of the velocities as an evidence for the appropriateness of the approach. This analysis provides a more thorough understanding of the properties of the Deep Impact dust cloud.

Context: In 2004 asteroid (2867) Steins has been selected as a flyby target for the Rosetta mission. Determination of its spin period and the orientation of its rotation axis are essential for optimization of the flyby planning. Aim: Measurement of the rotation period and light curve of asteroid (2867) Steins at a phase angle larger than achievable from ground based observations, providing a high quality data set to contribute to the determination of the orientation of the spin axis and of the pole direction. Methods: On March 11, 2006, asteroid (2867) Steins was observed continuously for 24 hours with the scientific camera system OSIRIS onboard Rosetta. The phase angle was 41.7 degrees, larger than the maximum phase angle of 30 degrees when Steins is observed from Earth. A total of 238 images, covering four rotation periods without interruption, were acquired. Results: The light curve of (2867) Steins is double peaked with an amplitude of $\approx$ 0.23 mag. The rotation period is 6.052 $\pm$ 0.007 hours. The continuous observations over four rotation periods exclude the possibility of period ambiguities. There is no indication of deviation from a principal axis rotation state. Assuming a slope parameter of G = 0.15, the absolute visual magnitude of Steins is 13.05 $\pm$ 0.03.

Comet Hale-Bopp was observed five times with ISOPHOT, the photometer on board the Infrared Space Observatory (ISO). Each time broadband photometry was performed using 4 different detectors, 5 apertures and, 10 filters covering the range between 3.6 to 175 $\mu$m. Calibration measurements using the internal Fine Calibration Source were done together with each individual measurement. Background observations were performed with identical instrument settings at the same positions on the sky several days after the comet observations. The observation strategy and the initial data reduction steps are described in some detail and the resulting in-band power values of the Hale-Bopp observations and their uncertainties are derived. Initial reduction of these measurements was performed in 3 steps: (1) processing of raw data by removing instrumental and energetic particle effects, (2) averaging of the individual signals, and (3) determination of the detector responsivities and their uncertainties. The detector signal is determined by two different methods and the best value is chosen. At the present level of processing uncertainties range from 10% to a factor of 3 for the low power levels at short wavelengths. The in-band power levels at different wavelengths varied over 3 orders of magnitude.