Arabia Terra is a region of Mars where signs of past-water occurrence are recorded in several landforms. Broad and local scale geomorphological, compositional and hydrological analyses point towards pervasive fluid circulation through time. In this work we focus on mound fields located in the interior of three casters larger than 40 km (Firsoff, Kotido and unnamed crater 20 km to the east) and showing strong morphological and textural resemblance to terrestrial mud volcanoes and spring-related features. We infer that these landforms likely testify the presence of a pressurized fluid reservoir at depth and past fluid upwelling. We have performed morphometric analyses to characterize the mound morphologies and consequently retrieve an accurate automated mapping of the mounds within the craters for spatial distribution and fractal clustering analysis. The outcome of the fractal clustering yields information about the possible extent of the percolating fracture network at depth below the craters. We have been able to constrain the depth of the pressurized fluid reservoir between ~2.5 and 3.2 km of depth and hence, we propose that mounds and mounds alignments are most likely associated to the presence of fissure ridges and fluid outflow. Their process of formation is genetically linked to the formation of large intra-crater bulges previously interpreted as large scale spring deposits. The overburden removal caused by the impact crater formation is the inferred triggering mechanism for fluid pressurization and upwelling, that through time led to the formation of the intra-crater bulges and, after compaction and sealing, to the widespread mound fields in their surroundings.

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.

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.

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.

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

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.

The European Space Agency's Rosetta spacecraft passed by the main belt asteroid (21) Lutetia the 10th July 2010. With its ~100km size, Lutetia is one of the largest asteroids ever imaged by a spacecraft. During the flyby, the on-board OSIRIS imaging system acquired spectacular images of Lutetia's northern hemisphere revealing a complex surface scarred by numerous impact craters, reaching the maximum dimension of about 55km. In this paper, we assess the cratering history of the asteroid. For this purpose, we apply 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 us to interpret the observed crater size-frequency distribution (SFD) and constrain the cratering history. Thanks to this approach, we derive the crater retention age of several regions on Lutetia, namely the time lapsed since their formation or global surface reset. We also investigate the influence of various factors -like Lutetia's bulk structure and crater obliteration- on the observed crater SFDs and the estimated surface ages. From our analysis, it emerges that Lutetia underwent a complex collisional evolution, involving major local resurfacing events till recent times. The difference in crater density between the youngest and oldest recognized units implies a difference in age of more than a factor of 10. The youngest unit (Beatica) has an estimated age of tens to hundreds of Myr, while the oldest one (Achaia) formed during a period when the bombardment of asteroids was more intense than the current one, presumably around 3.6Gyr ago or older.

In this paper we present a crater age determination of several terrains associated with the Raditladi and Rachmaninoff basins. These basins were discovered during the first and third MESSENGER flybys of Mercury, respectively. One of the most interesting features of both basins is their relatively fresh appearance. The young age of both basins is confirmed by our analysis on the basis of age determination via crater chronology. The derived Rachmaninoff and Raditladi basin model ages are about 3.6 Ga and 1.1 Ga, respectively. Moreover, we also constrain the age of the smooth plains within the basins' floors. This analysis shows that Mercury had volcanic activity until recent time, possibly to about 1 Ga or less. We find that some of the crater size-frequency distributions investigated suggest the presence of a layered target. Therefore, within this work we address the importance of considering terrain parameters, as geo-mechanical properties and layering, into the process of age determination. We also comment on the likelihood of the availability of impactors able to form basins with the sizes of Rachmaninoff and Raditladi in relatively recent times.

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.

Model Production Function chronology uses dynamic models of the Main Belt Asteroids (MBAs) and Near Earth Objects (NEOs) to derive the impactor flux to a target body. This is converted into the crater size-frequency-distribution for a specific planetary surface, and calibrated using the radiometric ages of different regions of the Moon's surface. This new approach has been applied to the crater counts on Mariner 10 images of the highlands and of several large impact basins on Mercury. MPF estimates for the plains show younger ages than those of previous chronologies. Assuming a variable uppermost layering of the Hermean crust, the age of the Caloris interior plains may be as young as 3.59 Ga, in agreement with MESSENGER results that imply that long-term volcanism overcame contractional tectonics. The MPF chronology also suggests a variable projectile flux through time, coherent with the MBAs for ancient periods and then gradually comparable also to the NEOs.

In this paper we present a new method for dating the surface of the Moon, obtained by modeling the incoming flux of impactors and converting it into a size distribution of resulting craters. We compare the results from this model with the standard chronology for the Moon showing their similarities and discrepancies. In particular, we find indications of a non-constant impactor flux in the last 500 Myr and also discuss the implications of our findings for the Late Heavy Bombardment hypothesis. We also show the potential of our model for accurate dating of other inner Solar System bodies, by applying it to Mercury.