Open Access
Issue
A&A
Volume 672, April 2023
Article Number A136
Number of page(s) 18
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202245614
Published online 14 April 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Comet 67P/Churyumov-Gerasimenko (hereafter 67P) was the main target of the Rosetta mission of the European Space Agency. Launched in 2004, Rosetta took ten years to reach the comet before orbiting around it for ~25 months, from July 2014 to September 2016, permitting an in-depth investigation of the 67P nucleus morphology, physical properties, and composition, and of the cometary activity and the dust–gas interaction on the nucleus surface and the inner coma at different heliocentric distances. For the first time in space exploration Rosetta delivered a lander, Philae, on a cometary surface on 12 November 2014. Even if Philae rebounded from the original selected landing site, after an adventurous trajectory and a second rebound (O’Rourke et al. 2020), it reached the Abydos site where most of the foreseen in situ measurements were successfully achieved.

Rosetta revealed a complex morphology of the nucleus, with different kinds of terrains (Thomas et al. 2015), including layers, boulders, cliffs, and pits, sometime active (Vincent et al. 2015). The cometary surface shows pervasive fractures ranging from millimeters (Bibring et al. 2015) to several tens of meters long produced by thermal insolation weathering (El-Maarry et al. 2015), as well as goosebumps or clod features on a scale of a few meters (Sierks et al. 2015; Davidsson et al. 2016; Fornasier et al. 2021) interpreted as remnants of the original pebbles or results of fracturing processes. Twenty-six regions, named after Egyptian deities, were identified based on the surface geomor-phological properties (see El-Maarry et al. 2015, 2016 for the cometary regions definition and location). The bilobate shape of the nucleus, which shows extensive layering but with different centers of gravity between the large and small lobes, is associated with a binary structure resulting from the collision at low speed of two distinct bodies in the early Solar System (Massironi et al. 2015). The binary structure interpretation is also supported by the different mechanical and physical properties reported for the two lobes (El-Maarry et al. 2016; Fornasier et al. 2021).

The comet is dark with a geometric albedo of 6.5±0.2% at 649 nm (Fornasier et al. 2015). The nucleus composition is dominated by refractory material mixed with opaque minerals and organics. The spectrum is red (i.e., the reflectance increases in a steep way with the wavelength) and is characterized by a wide absorption band in the 2.8–3.6 µm region indicating the presence of a complex mixture of organics (Capaccioni et al. 2015; Quirico et al. 2016). The latest recalibration of the Visible, InfraRed, and Thermal Imaging Spectrometer (VIRTIS) gave evidence of different structures in the broad band, attributed to ammonium salts (Poch et al. 2020) and aliphatic organics (Raponi et al. 2020), with a possible contribution from hydroxylated amorphous silicates to the overall absorption (Mennella et al. 2020).

The nucleus shows compositional heterogeneities on several spatial scales, resulting in different spectral slopes and albedo in regional and local areas. On the dark and red average cometary terrain, exposures of volatiles stand out because they are very bright and with a bluer spectrum (i.e., less steep). Two volatile species were detected as exposed ice on comet 67P, mainly crystalline water ice (De Sanctis et al. 2015; Barucci et al. 2016; Filacchione et al. 2016a) and also carbon dioxide, the latter found for the very first time exposed at a comet surface (Filacchione et al. 2016b).

Joint observations of the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) and VIRTIS spectrometers have proven that the bright and spectrally bluer features observed with the cameras display the typical water ice bands in the infrared spectra. Based on this correlation, a number of bluer and bright features detected with OSIRIS have been attributed to exposure of water ice (Barucci et al. 2016). Pommerol et al. (2015) reported the first OSIRIS detection of volatile exposures on comet 67P with features being 5–10 times brighter than the cometary dark terrain. Deshapriya et al. (2018) generated the first catalog of volatile exposures including 57 entries. Other studies highlight the presence of bright spots associated with water ice exposures in the northern hemisphere (Pommerol et al. 2015; Fornasier et al. 2015; Raponi et al. 2016; Barucci et al. 2016; Filacchione et al. 2016a; Lucchetti et al. 2017; La Forgia et al. 2015; De Sanctis et al. 2015; Oklay et al. 2017) and in the southern hemisphere (Fornasier et al. 2016, 2019a, 2021; Deshapriya et al. 2016, 2018; Hasselmann et al. 2019; Hoang et al. 2020), sometimes freshly exposed on the surface after cliff collapses or outbursts (Pajola et al. 2017a; Agarwal et al. 2017; Filacchione et al. 2016a) or due to the mechanical action of Philae (O’Rourke et al. 2020).

The estimated water ice abundance varies from a few percent (Barucci et al. 2016; Filacchione et al. 2016a,b; Raponi et al. 2016; De Sanctis et al. 2015; Ciarniello et al. 2016) to more than 20–30% in several bright areas observed in the Imhotep, Seth, Khonsu, Bes, Anhur, and Wosret regions (Deshapriya et al. 2016, 2018; Oklay et al. 2017; Pajola et al. 2017a; Fornasier et al. 2016, 2017, 2019a, 2021; Hasselmann et al. 2019; Hoang et al. 2020), with peaks up to ~50–80% in few localized tiny bright spots (Oklay et al. 2017; Hoang et al. 2020; O’Rourke et al. 2020; Fornasier et al. 2021), indicating fresh exposures of volatiles.

In this paper we present the most extended catalog of exposures of volatiles of comet 67P built upon a systematic analysis of the color sequences acquired with the OSIRIS cameras. We investigate their distribution in the different cometary regions and morphological terrains, their spectral slope evolution, their size distribution, and their duration with the aim of understanding volatile properties in comets, and of constraining cometary models.

2 Observations and methodology

The analysis is based on data from the OSIRIS imaging system of the Rosetta mission. OSIRIS included two cameras, the Narrow Angle Camera (NAC) for the high-resolution study of the nucleus, and the Wide Angle Camera (WAC) for the coma investigation (Keller et al. 2007).

We analyzed medium- and high-resolution images acquired with the NAC camera with several filters in the 250–1000 nm range, investigating 109 different color sequences taken between August 2014 and September 2016, and covering a spatial resolution from a few m px−1 up to 0.1 m px−1. We searched in the OSIRIS archive all the NAC spectrophotometric sequences pointing to the 67P nucleus and having at least three filters. Exposures of volatiles are usually brighter than the comet dark terrain, and are characterized by a neutral to moderate spectral slope in the visible range, which has been proven to be associated with a local enrichment in water ice thanks to joint observations carried out with the OSIRIS cameras and the VIRTIS spectrometer (Barucci et al. 2016; Filacchione et al. 2016a). Thus, to determine whether a bright feature on the surface is ice dominated and not simply brighter because of illumination conditions, information on the reflectance value and on the spectral slope of a region of interest (ROI) are needed. This implies that some bright spots observed only in the NAC orange filter centered at 649 nm, the filter most frequently used to generate the nucleus shape model or to investigate the comet morphology, are not included in our analysis. With only one filter available we cannot determine if the higher brightness is due to a geometric effect, to the presence of a bright mineral, or to a real exposure of ice.

Therefore, we applied the following methodology to identify exposures of volatiles on the 67P nucleus. Bright spots exposing volatiles should be both brighter (by at least 50%) than the comet dark terrain, and should have neutral to moderate spectral slope values in the visible range (535–882 nm). The adopted upper limit in the spectral slope value is 11%/(100nm), but usually bright spots have spectral slopes much lower than 8%/(100nm), and often close to zero. Moreover, we considered only the bright spots larger than 3 pixels. Smaller features are indeed difficult to characterize because of residuals in the images co-registration process.

We used the NAC images generated by the instrument pipeline (Tubiana et al. 2015) corrected by bias, flat field, geometric distortion, absolutely calibrated in radiance, and finally converted in radiance factor (also named I/F) Radiancefactor(i,e,α,λ)=πI(i,e,α,λ)Fsunλ,$ {\rm{Radiancefactor}}\left( {i,e,\alpha ,\lambda } \right) = {{\pi {\bf{I}}\left( {i,e,\alpha ,\lambda } \right)} \over {F\,su{n_\lambda }}}, $(1)

where I is the scattered radiance at a given incidence (i), emission (e), phase (α) angles and wavelength (λ), and F sunλ is the incoming irradiance of the Sun at the heliocentric distance of the comet and at a given wavelength (λ).

As done in previous studies of the 67P nucleus, the NAC images of a given sequence were first co-registered using the F22 NAC filter (centered at 649.2 nm) as a reference, then corrected by the illumination conditions using the Lommel–Seeliger disk function and the 3D stereophotoclinometric shape model of the 67P nucleus (Jorda et al. 2016), adopting the same methodology already presented in Hasselmann et al. (2019) and Fornasier et al. (2017, 2019a).

We also created RGB images with the STIFF software, which converts scientific FITS images to TIFF (Bertin 2012), mostly using the filters centered at 882 nm (R), 649 nm (G), and 480 nm (B). These RGB images are very helpful in identifying volatile exposures since they look bright and blue compared to the dark and red cometary terrain.

For each bright feature, the spectral slope (S l) was computed in the 535–882 nm range as: sl=R882R535R535×(882  nm535  nm),$ sl = {{{R_{882}} - {R_{535}}} \over {{R_{535}} \times \left( {882\,\,{\rm{nm}} - 535\,\,{\rm{nm}}} \right)}}, $(2)

where R882 and R535 are the radiance factors in the filters centered at 882 nm and 535 nm, respectively.

Details on the observing conditions are reported in Table A.1.

thumbnail Fig. 1

Maps of ice exposure on comet 67P. The color-coding is as follows: red, pre-perihelion (August 2014–May 2015); cyan, perihelion (June–October 2015); green, post-perihelion (November 2015–September 2016). The spectrally blue spots, those having negative spectral slope (<−3%/(100 nm)) in the 535–882nm range, are shown in light green (post-perihelion) and orange (pre-perihelion). The symbol size represents three ranges of volatile exposure area; they are not in scale compared to the nucleus surface (51.74 km2 in total; Thomas et al. 2018), but are enlarged for clarity.

3 Catalog of volatiles exposures

We identified and characterized 603 bright spots (hereafter BS) having a spectral slope much lower than the typical value of the cometary dark terrain (Table A.1), thus indicating local exposures of volatiles, very likely water ice on the 67P nucleus. This is the most complete catalog of volatile exposures on 67P published to date, increasing by a factor of ~10 the number of identified bright spots on the comet compared to data already published in the literature. However, this catalog does not include the totality of the volatile exposures for the following reasons: (i) some BS might have been present on the surface but not captured by OSIRIS observations because Rosetta was pointing elsewhere or because they fully sublimated between two consecutive OSIRIS sequences covering a given region; (ii) in this catalog we included only the BS observed within color sequences, thus we do not consider those captured by a single filter where the spectrophotometric analysis cannot be performed; (iii) we considered only the bright spots larger than 3 pixels in size; (iv) in the case of clusters of icy exposure, not all the individual points (often smaller than 3 pixels in size) were counted.

For each bright feature we computed the surface, the coordinates, the spectral slope, and the minimum duration, when possible, and for a few showing negative spectral slope values we also estimated the water ice abundance using geographical mixtures of the comet dark terrain and water ice. The full catalog of bright spots is reported in Table A.1, and their distribution across the nucleus is shown in Fig. 1.

4 Bright feature distribution and type

Bright patches and spots could be found isolated on the nucleus surface or grouped in clusters, usually at the bottom of cliffs. In Tables 1 and A.1 we report the BS type following the Deshapriya et al. (2018) classification scheme: type (1) isolated BS on smooth terrains; type (2) isolated BS close to irregular structures; type (3) BS on top of boulders; type (4) clusters of bright patches and BS.

Examples of the different types of volatiles exposures are reported in Figs. 24, while in Table 1 we summarize the BS identification per type. The majority of them are of type 2, and are thus identified close to irregular structures. This is quite expected because of the complex geomorphology of the comet.

The largest icy exposure belongs to type 2 and was observed on Imhotep shortly after the perihelion passage, on 23 August 2015 (BS 188 in Table A.1). This bright patch occupied a vast surface of ~5260 m2, and was repeatedly observed for 4 h by NAC sequences capturing that region, and was still observed one week later, even if part of it sublimated during this time lapse. It is worth mentioning that this area was brighter and spectrally bluer than the comet dark terrain, but its spectral slope was relatively high compared to other BS (around 10%/(100 nm)), indicating a local surface enrichment of volatiles, but highly mixed with the cometary dust. In addition, the spatial resolution was relatively low (about 6 m px−1), impeding an accurate study of this BS.

More than one-third of the type 2 BS are smaller than 1 m2, and the average size is of 86 m2, or 72 m2 when excluding the largest patch previously described. For 136 out of 373 BS of type 2, we estimated their minimum duration (i.e., the time between the first and last sequences capturing a BS), and its average value is 41 days.

Smooth terrains (type 1, Fig. 2) host few BS, but they tend to be larger (~370 m2) than the isolated ones observed close to irregular structures or on boulders, and their average duration is the shortest (about 18 days, Table 1). This may be attributed to mixing processes with the surrounding dust, and/or to a longer illumination time and/or intensity compared to rough terrains where mutual shadows favor a longer ice survival.

Conversely, BS on boulders (type 3, Fig. 3) are found to be smaller in size but with the longest duration (72 days). This may be associated with the presence of fractures and small cavities on boulders, which slow down the volatile sublimation. Bright spots on boulders may also be fed by internal reservoirs of volatiles. During the second Philae touch down, the lander imprinted in a boulder, revealing a 3.5 m2 bright area containing the primordial water ice embedded inside it (O’Rourke et al. 2020). They estimated a water ice fraction of 46% and a dust-to-ice mass ratio of 2.3 in this boulder.

Clusters of bright spots (type 4, Fig. 4) are located at the bases of cliffs and likely formed as a result of cliff collapses, such as the large clustered features (CFs) named CF1, CF2, and CF3 described in Oklay et al. (2017). Some of them, especially in post-perihelion images, look clearly associated with frost recondensation, like the ones in the Ash and Babi regions shown in Fig. 4. We identified 145 BS in clusters, located mainly in the Anuket, Ash, Aten, Babi, Geb, Hatmehit, and Seth regions. The individual BS in clusters are relatively small (23 m2), and some of them are long-lived features; for example, the blue enriched areas in the Seth alcove (bottom right part of Fig. 4) was observable for more than two years. Tiny spots (i.e., under 1 m2) were frequently identified inside the Hatmehit rim in high-resolution post-perihelion images (Hoang et al. 2020), but mostly with a duration of a few minutes or 1 day, indicating the presence of frost, with the notable exception of a few BS that survived for 6–8 days.

Volatile exposures are found at different latitudes post-perihelion. During the pre-perihelion observations they are more concentrated in the equator and mid-latitudes, between −40° and 40° (Fig. 5).

thumbnail Fig. 2

Example of isolated bright features on smooth terrains (top) and close to irregular structures (bottom), feature types 1 and 2, respectively, following the Deshapriya et al. (2018) classification scheme. The bright spot numbers correspond to those listed in Table A.1.

Table 1

Volatile exposures types from the catalog here presented (Table A.1) following the Deshapriya et al. (2018) classification scheme.

thumbnail Fig. 3

Example of bright features resting on boulders (BS 32 in Anuket, BS 374 in Seth, and BS 542 and BS 544 in Anhur), type 3 following Deshapriya et al. (2018) classification scheme. The bright spot numbers correspond to those listed in Table A.1. Several BS show blue colors, indicating a small or negative spectral slope value (see Table A.1).

thumbnail Fig. 4

Example of clusters of bright spots, type 4 following the Deshapriya et al. (2018) classification scheme. The bright spot numbers correspond to those listed in Table A.1. Several BS show blue colors, indicating a small or negative spectral slope value (see Table A.1).

thumbnail Fig. 5

Frequency of the bright spots vs. latitude. The histogram in black represents the pre-perihelion BS.

thumbnail Fig. 6

Histograms showing the bright spot spectral slopes, evaluated in the 535–882 nm wavelength range, pre-perihelion (August 2014–May 2015), during perihelion (June–October 2015), and post-perihelion (November 2015–September 2016).

5 Spectral slope distribution of the bright spots

We investigated the spectral slope distribution, evaluated in the 535–882 nm range, of the BS during the different comet orbital periods (Fig. 6), which we defined as follows: pre-perihelion from August 2014 to the end of May 2015; during perihelion from June to October 2015; post-perihelion from November 2015 to the end of the Rosetta mission on 30 September 2016.

The volatile exposures show a distinct spectral slope distribution in the post-perihelion period, with the presence of several bright spots having negative spectral slope values, while pre-perihelion and during perihelion the BS spectral slope values were close to zero or moderately positive, with very few BS showing negative values (Fig. 6). It should be noted that the observations at perihelion were acquired at relatively high distances and low spatial resolution (3–10 m px−1), preventing the identification of square-meter-sized bright spots. The observing conditions explain the limited number of bright spots observed close to perihelion. The average spectral slope values in the 535–882 nm range of pre-perihelion, perihelion, and post-perihelion bright spots are 5.0±0.3 %/(100nm), 6.8±0.4%/(100nm), and 0.9±0.2%/(100 nm), respectively. To exclude that the lower spectral slope values of the BS in the post-perihelion images are related to spatial resolution effects, we investigated the BS spectral slope distribution for pre- and post-perihelion data acquired at similar high resolution, between 0.33 m px−1 (the highest pre-perihelion resolution available) and 1 m px−1. The histograms shown in Fig. 7 confirm the trend observed in Fig. 6, and thus the decrease in the BS spectral slope in post-perihelion images. The average spectral slope values of the BS investigated at similar high resolution is of 4.9±0.3%/(100nm) for the pre-perihelion BS, and 1.3±0.3%/(100nm) for the post-perihelion BS.

Moreover, we found that 57 spots have unusually negative slope values (i.e., below −3%/(100nm)) in the 535–882 nm range, which we call blue spots. All except one were observed after perihelion. The only pre-perihelion one was detected on 5 September 2014 in Imhotep (BS 42 in Table A.1). This region also hosts the first blue sloped post-perihelion BS, which was detected at the end of November 2015 (BS 205 in Table A.1). In the same period a blue sloped BS was also observed in Anuket (BS 199).

Twenty-two of these blue spots were observed in the Anhur region (a few examples are shown in Fig. 8), close to the canyon structure (see Fornasier et al. 2017 for the Anhur morphological description) that also hosted one of the brightest outbursts reported for comet 67P during the Rosetta observations, called the perihelion outburst, which took place on 12 August 2015 (Fornasier et al. 2019a). Another region showing blue spots (a total of 13) is Khonsu. These features are located in the low bank area (i.e., a flat area between −20° and −30° latitude) defined and investigated by Hasselmann et al. (2019), where they report important morphological changes. This area was also the source of several activity events during the perihelion passage, including an outburst on 14 September 2015 (Vincent et al. 2016; Hasselmann et al. 2019).

The presence of blue spots is not restricted to the big lobe or to the southern hemisphere of the comet (where Anhur and Khonsu are located). They were also detected in the small lobe (e.g., four BS in Wosret) and in the northern hemisphere (e.g., four BS each in the Babi and Seth regions), even if they are less frequent.

The evolution of the BS spectral slope toward smaller values in post-perihelion images is also confirmed in the analysis of the BS in individual regions, as shown in Table 2, where the slope systematically decreases in outbound orbits compared to the pre-perihelion images for regions observed during both periods and showing at least a few BS. It is worth noting that in Table 2 the spectral slope of the BS increased during the perihelion passage, while several studies reported a bluing of the comet colors (i.e., a decrease in the spectral slope) during the perihelion passage due to the seasonal cycle of water (Fornasier et al. 2016; Ciarniello et al. 2016, 2022). The higher spectral slope value of BS at perihelion is very likely an observational bias, due to the low spatial resolution of the data. In fact, BS are usually meter sized, thus during perihelion their spectral behavior is highly affected by the cometary red dark terrain, because BS are not resolved.

Most of the blue spots are located in shadowed areas, as shown in Fig. 8 or in Fornasier et al. (2021, their Fig. 10), and some others were found within frost fronts at the bottom of cliffs. The majority of the blue spots are only a few m2 in size or even smaller, revealed thus in high-resolution images, while some other spots are a few tens of m2 and the largest blue spot, with a surface of about 50 m2, is located in the Babi region (BS 395 in Table A.1). The estimated lifetime of the blue spots ranges from at least 20 minutes (BS 599) to nearly two weeks for BS 489 (Fig. 8).

The first bright blue spots we found in our analysis, and which motivated a deeper study, are located in the Anhur region and were observed in high-resolution images acquired on 16 July 2016 when the comet was 3.4 au outbound, at a resolution of 0.19 m px−1 (Fig. 8). BS 486 and 487 in Table A.1 display a sharp decrease in the reflectance after 650 nm, with spectral slope value (in the 535–882nm range) of −8 and −4.3%/(100nm), respectively. These BS are sub-meter sized, with a dimension of 0.95 m × 0.76 m for BS 486 and 0.76 m × 0.56 m for the BS 487. The decrease in reflectance is not associated with image saturation, nor to local fast sublimation between the sequence color images (that last 230 s) because the order of acquisition of the filters was not in order of increasing wavelength. The sequence order was 649–535–480–360–700–882–989 nm. Moreover, some BS survived several days, as shown in Fig. 8, excluding fast sublimation.

This peculiar spectral behavior was never reported in the existing literature in pre-perihelion BS, which were usually spectrally flat in the visible range. We observed only one pre-perihelion bright spot having negative slope, located in the Imhotep region (BS 42 in Table A.1), very likely related to frost because of its short lifetime (a few hours).

We attempted to estimate the water ice content of the Anhur blue bright patches using linear mixing model of the cometary dark terrain and water ice with different grain sizes, using the method detailed in Fornasier et al. (2016, 2019a): R=p×Rice+(1p)×RDT.$ R = p \times {R_{{\rm{ice}}}} + \left( {1 - p} \right) \times {R_{{\rm{DT}}}} .$(3)

Here R is the reflectance of the bright patches, Rice and RDT are respectively the reflectance of the water ice and of the cometary dark terrain, and p is the relative surface fraction of water ice.

Water ice spectra were produced using the Warren and Brandt (2008) optical constants and the Hapke radiative transfer model (Hapke 2002). To produce the absolute reflectance of the regions of interest, we applied the Hapke (2002) photometric model correction with parameters derived from Fornasier et al. (2016). However, we neglect the disk function correction, which is set to unity, because none of the shape models produced by the OSIRIS team has a spatial resolution high enough to correctly reproduce these tiny features. Considering this limitation, the estimated water ice abundance is about 20% for spot BS 489, and 75% for the brightest spot (BS 486). Conversely to most of the pre-perihelion BS observed with OSIRIS, which were matched by linear mixtures of the cometary dark terrain and water ice with grain sizes in the 30–100 µm range (Fornasier et al. 2016, 2021, 2019a; Oklay et al. 2017; O’Rourke et al. 2020), the ones analyzed here cannot be reproduced by water ice with relatively small grains. Even the BS showing a positive slope (BS 489, magenta square in Fig. 8) is better fitted by a model including water ice with large grains (1000 µm) than by the model including water ice with 100 µm grains. BS 487 is about five times brighter than the comet dark terrain, and its reflectance is reproduced by a mixture including 46% of water ice with large grains (3000 µm). Finally, BS 486 is 7-8 times brighter than the comet dark terrain (at 650 nm), and is dominated by water ice (75%) with very large grains (5000 µm, Fig. 8).

The fact that the BS spectrophotometry is reproduced by water ice with large grains in areal mixing with the cometary dark terrain is not unique, and was already reported for the modeling of some bright patches in comet 67P observed with VIRTIS. In the infrared the spectra of the exposures of volatiles show the water ice absorption bands, permitting a deeper analysis of the composition and allows us to better constrain the components’ physical parameters (grains size, temperature). Both Raponi et al. (2016) and Filacchione et al. (2016a) found that the compositional models reproducing the analyzed BS require the presence of both areal and intimate mixtures of the cometary dark terrain and water ice, this last having different grain sizes: a few tens of microns for the intimate mixture, and millimeter-sized grains for the areal mixture, as found in this paper for the tiny blue BS in Anhur (Fig. 8). The BS analyzed by Raponi et al. (2016) and Filacchione et al. (2016a) were detected pre-perihelion; they showed a positive spectral slope in the visible range, survived for 2–4 months with variable abundance of the water ice related both to the seasonal and diurnal cycle, and were observed in different regions. Unfortunately, infrared spectra of the negative sloped BS reported here are not available because the VIRTIS infrared channel has not been operational since May 2015. The presence of millimeter-sized grains in some water ice exposures of comet 67P may be due to the sintering of smaller water ice grains or to the growth of ice crystals produced by vapor diffused from the colder and water ice enriched material beneath the surface (Filacchione et al. 2016a).

However, water ice in comets is usually characterized by having small to medium grain size values: very fine grains 2 µm) were used to model the frost fronts in the Hapi region (De Sanctis et al. 2015), to model the overall decrease in the 67P nucleus spectral slope with the increasing activity approaching perihelion (Filacchione et al. 2016c; Ciarniello et al. 2016), to model the material excavated by the Deep Impact impactor on 9P/Tempel 1 (Sunshine et al. 2007), and were also found in the 103P/Hartley coma (Protopapa et al. 2014) and 17P/Holmes outbursts (Yang et al. 2009); coarser grains (30–100 µm) were used to model several BS of comet 67P (Filacchione et al. 2016a; Barucci et al. 2016; Pommerol et al. 2015; Fornasier et al. 2016; Oklay et al. 2017), and also on ice patches in comet 9P/Tempel 1 (Sunshine et al. 2006).

The high water ice abundance is not exclusive for the blue sloped BS reported in Fig. 8. In the literature high water ice abundances were reported for a few other BS (see Table A.1 for the BS numbers): over 80% in BS 480 (Hoang et al. 2020), 64–69.5% in BS 476 (Fornasier et al. 2021), and ~46% (O’Rourke et al. 2020) in the boulder onto which Philae impacted (BS 385), exposing primordial ice. All these BS, located in the Wosret region, should be fresh exposures of volatiles. Other BS in the Anhur region also have a large water ice fraction (30–40%), for example BS 461 (magenta symbol in Fig. 13, bottom panel, from Fornasier et al. 2019a), or the few large BS (1500 m2 each) observed pre-perihelion in the Anhur-Bes boundary (Fornasier et al. 2016, BS 143 and BS 144). Oklay et al. (2017) reported long-lived BS in the Imhotep region, with estimated water ice content up to 48% in the feature called IF2.

Some of the post-perihelion blue spots may also be associated with frosts, especially those having short lifetimes. Frost was repeatedly observed in areas exiting from the cometary night or shadows, as shown in Fig. 9, and was already reported in the Anhur region (Fornasier et al. 2017, 2019a), near the final landing site Abydos (Hoang et al. 2020), and in Seth (Lucchetti et al. 2017), while it was less frequently observed pre-perihelion, with the notable exception of the Hapi region (De Sanctis et al. 2015).

thumbnail Fig. 7

Histograms showing the bright spots spectral slopes, evaluated in the 535–882 nm wavelength range, in images having similar high resolution (in the 0.33–1 m px−1 range), in the pre-perihelion (August 2014–May 2015), and post-perihelion (November 2015–September 2016) periods.

thumbnail Fig. 8

Duration, spectrophotometry and compositional modeling of some blue bright spots in the Anhur region. Top: RGB images showing the blue BS observed in boulders named A, B, and C on 16 July 2016 (top central panel). This area is located close to the canyon-like structure of Anhur (Fornasier et al. 2017). Frosts and tiny blue spots were already present on 25 July 2016 images (top left panel) near the boulders named A and C, and a few close to boulder A were still visible in the 30 July 2016 images (top right panel). Bottom: relative reflectance and I/F of four selected blue BS located in the aforementioned A and B boulders: BS 486 (red star, slope = −8.0%/(100nm)), BS 487 (blue asterisk, slope = −4.3%/(100nm)), BS 488 (green triangle, slope = −1.8%/(100 nm)), and BS 489 (magenta square, slope = 4.6%/(100nm)) in Table A.1. DT indicates the cometary dark terrain. Bottom right: linear mixing models of the cometary dark terrain and water ice with different grain sizes reproducing the spectral behavior of the blue bright features. The estimated water ice abundance is indicated for the different models.

Table 2

Bright spots per region and spectral slope values.

thumbnail Fig. 9

Examples of RGB images showing frosts during the outbound orbit of comet 67P, near 3 au.

6 Volatile exposures and cometary activity

Volatile exposures are directly linked to cometary activity. They may either be sources of jets or appear as the result of cometary activity that generated self-cleaning and erosion of the nucleus, or they may be produced by morphological changes that expose the underlying volatile-rich materials.

Several relatively faint jets were observed directly departing from BS, as reported in Fornasier et al. (2019b), see, e.g., their Fig. 6), and some other BS appeared after cometary activity departing from their location or surroundings (Barucci et al. 2016; Deshapriya et al. 2016; Vincent et al. 2016; Fornasier et al. 2017).

Examples of self-cleaning of the comet are the two large bright patches reported by Fornasier et al. (2016, 2017) close to the Anhur-Bes regions boundary (BS 143 and 144 in Table A.1) that suddenly appeared in April–May 2015 in a smooth area. The ice survived exposed for about 7–10 days, and then these BS fully sublimated leaving a surface spectrally indistinguishable from the surrounding terrain. Notably, in the patch located within Anhur, the VIRTIS spectrometer detected in March 2015 for the very first time the exposure of carbon dioxide ice on a comet (Filacchione et al. 2016b). The discovery of carbon dioxide and water ices in this region demonstrates that different layers of volatiles are present within the nucleus, and points to compositional heterogeneity on large scales, on the order of tens of meters, on comet 67P (Fornasier et al. 2016).

Examples of ices exposed after activity and morphological changes are the following: the area rich in fresh water ice six times brighter than the surrounding dark terrain that appeared in the Aswan site after the cliff collapse observed on 10 July 2015, generated by an outburst, and that stayed exposed until at least the end of 2015 (Pajola et al. 2017a, BS168 in Table A.1); a 15×5m2 bright patch formed inside a circular basin in Imhotep after an outburst on 3 July 2016 and that lasted for at least seven weeks (Agarwal et al. 2017, BS 485); a new scarp, 140 m long and 10 m high, formed in the Bes region, near the Anhur-Bes boundary and very close to BS 144 (where CO2 ice was detected), exposing in the talus volatiles with water ice abundance estimated at 17±2% (Fornasier et al. 2017); a second 10 m high scarp, located in the Anhur region, formed in January 2016 and showing at its base a surface of about 500 m2, brighter and bluer than the cometary dark terrain, exposing water ice for at least 5 months with estimated abundances of 26–30% (Fornasier et al. 2019a, BS 261); the appearance of a BS on a boulder in the Khonsu region at the beginning of January 2016, surviving for more than 6 months and associated with a cometary jet close to the perihelion passage (Deshapriya et al. 2018); again in the Khonsu region, a number of important morphological changes with dust removal up to 10 m in height and exposing several ice enriched patches (Hasselmann et al. 2019); water ice exposure in front of large expanding structures in Imhotep (Groussin et al. 2015, Deshapriya et al. 2018), observed shortly before the perihelion passage; the clustered features named CF2 and CF3 in Imhotep (Oklay et al. 2017), related to activity events and exposing water ice with abundances of ~25%; the mechanical action of Philae that impacted into a boulder in the Wosret region, exposing fresh ice with abundance estimated at 46% (O’Rourke et al. 2020).

The comparison between Fig. 1, showing the BS distribution on the nucleus, and Fig. 10, displaying the sources of jets and outburst observed close to the perihelion passage, clearly indicates a correlation between activity and water ice exposures, in particular for the Anhur, Bes, Khepry, Imhotep, Khonsu, Atum, Anubis, and Anuket regions. As already noted in Fornasier et al. (2019b), activity events are correlated with local compositional heteregeneities (i.e., with local exposure of volatiles) and activated by solar illumination.

thumbnail Fig. 10

Map with the sources of activity identified on the 67P nucleus in Vincent et al. (2016), blue asterisks, and Fornasier et al. (2019b).

7 Size distribution of the volatiles exposures

We investigated the surface distribution of the volatile exposures. The area of the spots ranges from sub-m2, for example the spots observed on the Hatmehit rim in May 2016 (Hoang et al. 2020) or in Anhur (Fig. 8), to a few thousand square meters, like the two 1500 m2 BS in the Anhur & Bes regions, observed in April–May 2015 (Fornasier et al. 2016), up to ~5300 m2 for the spot observed in Imhotep shortly after the perihelion passage (BS 188 in Table A.1).

Even if numerous bright spots are detected, the total surface of exposed water ice reported here (Table A.1) is 46 714 m2, which is 0.1% of the total 67P nucleus surface. This confirms that the surface of comet 67P is dominated by refractory dark terrains, while volatiles may occupy only a tiny areal fraction.

The majority of the volatile exposures are found in the regions of the big lobe. Their total area is about 42 800 m2, which correspond to 0.13% of the large lobe surface (31.66 km2, Thomas et al. 2018). The percentage of volatile exposures is six times smaller for the small lobe and the neck area, corresponding to the Hapi and Sobek regions, with values of 0.02% (3450 m2 over 17.26 km2) and 0.016% (460 m2 over 2.82 km2), respectively.

The histogram showing the frequency of the BS diameter, estimated from the BS area assuming a square shape, is reported in Fig. 11. Bright spots are predominantly small; about one-third of them are smaller than 3 m2, and 170 are smaller than 1 m2.

Figure 12 shows the area distribution of the bright spots identified in the pre- and post-perihelion images (Table A.1), observed at resolutions higher than 1 m px−1, and having a surface lower than 10 m2; because of the low resolution, we did not consider the observations close to the perihelion passage. Even if larger bright spots were present, we chose this limit in the analysis and to compute their surface mean values because in the pre-and post-perihelion data the great majority of these bright features have areas of a few square meters. With these criteria we count 98 pre-perihelion and 299 post-perihelion bright spots.

In addition to the different spectral slope distribution previously discussed, there is also a clear distinction in the areal distribution of the bright spots pre- and post-perihelion (Fig. 12). While pre-perihelion BS have a median surface of about 4.5 m2, in the post-perihelion images most of them are smaller than 1–2 m2, with an average value of 0.7 m2, clearly indicating that high spatial resolution is mandatory to identify ice exposures on cometary surfaces.

The largest number of post-perihelion BS is clearly associated with an observational bias; approaching the end of the mission Rosetta went closer to the cometary surface achieving very high spatial resolution. More than one-third of the BS reported in Table A.1 (216 BS) were observed at a spatial resolution better than 0.33 m px−1 in post-perihelion images. When comparing the BS distribution pre- and post-perihelion observed at similar high resolution (i.e., between 0.33 m and 1 m px−1; Fig. 13), the average value of post-perihelion BS surface increases to 3.9 m2, and the distribution is more compatible with the one observed pre-perihelion, even if about one-fourth of the post-perihelion BS are relatively small (<1.5 m2). It should be noted that small BS are still not observed at a similar resolution in pre-perihelion images.

We also investigated the area distribution for the different cometary regions showing more frequently exposures of volatiles. The regions having the higher number of BS (80–60) are located in the large lobe: Anhur, Imhotep, Ash, Khonsu, and Seth (Table 2), while Wosret and Hatmehit, with ~30 BS observed each, are the regions with more BS observed in the small lobe.

Our results on the bright spot dimensions support the findings of Ciarniello et al. (2022) and Fulle et al. (2020), who deduced that the bright spots on comets are exposure of the water ice enriched blocks (WEBs) forming, together with the refractory matrix, cometary nuclei, and whose dominant size is on the order of 0.5–1 m. WEBs should be formed of water ice-rich pebbles mixed with drier material, and exposed to the surface when the nucleus is eroded by the cometary activity. The fact that the majority of the bright spots are sub-meter sized is thus in agreement with these predictions and with the radar measurements the 67P comet provided by the Comet Nucleus Sounding Experiment by Radiowave Transmission (CONSERT), which indicate that the nucleus is homogeneous up to scales of a few meters (Ciarletti et al. 2017).

thumbnail Fig. 11

Histogram density of the bright spot diameters identified on the 67P nucleus and reported in this study. Only those having diameter smaller than 80 m are represented, for clarity. The cumulative probability is shown in blue.

thumbnail Fig. 12

Frequency of bright spot areas identified on the 67P nucleus for the pre-perihelion (top) and the post-perihelion (bottom) periods, observed at spatial resolution lower than 1 m px−1.

thumbnail Fig. 13

Frequency of the bright spot areas identified on the 67P nucleus for the pre-perihelion (top) and the post-perihelion (bottom) periods, and observed at comparable high resolution (0.33–1 m px−1).

8 Duration and evolution of the bright spots

The lifetime of bright spots was estimated as the time in which a bright spot remains visible in different observing sequences, and thus it should be considered a minimum duration (Table A.1). This estimation is biased by the observing frequency-conditions; therefore, the real lifetime is usually longer than that reported here.

Volatiles may survive exposed at the surface for a period varying from a few hours, in which case they are very likely frosts, to a few days and, for some of them, even to several months. A longer duration is usually found for water ice exposed after cliff collapse or the formation of new scarps, which may survive at the surface for several months, as detailed in the previous sections.

The evolution and water ice content of some clusters of BS was already reported in the literature, and it varies locally and with time: the CF1, CF2, and CF3 features in Imhotep contained spots with water ice fractions from 6.5% to 24.5% (Oklay et al. 2017); the water ice content of individual spots in the Hatmehit rim was estimated to be below 15% in late November and early December 2015, and exceeding 50% in late December 2015 and early 2016 (Hoang et al. 2020); joint OSIRIS and VIRTIS data analysis of BS indicates water ice abundances ranging from 0.1 to 7.2% (Barucci et al. 2016; Raponi et al. 2016), with seasonal variability during the BS lifetime (on the order of 2–4 months).

A few examples of BS duration are reported in Figs. 1417. For pre-perihelion images, BS 32 and BS 16 in Anuket survive for at least 2 and 3 months, respectively (Figs. 14, and 15), and BS 35 and BS 36 in Neith for about 2 months (Fig. 16). An example of post-perihelion duration is reported for Atum, where BS 253 survived exposed at the surface for at least 84 days, and BS 204 for more than 4 months (Fig. 17).

Figure 18 shows the spectral slope value versus the duration for all the bright spots here reported, and for which the minimum duration was evaluated. Black points are post-perihelion BS, while red and blue points indicate the pre-perihelion observations, for the August 2014–January 2015, and February–May 2015 periods, respectively.

We separated the pre-perihelion observations into two datasets to investigate if any correlation exists between the blue spot color evolution and the overall evolution of the comet color and/or composition, due to the seasonal cycle of water, already observed and reported in the literature (Fornasier et al. 2016; Filacchione et al. 2016b; Ciarniello et al. 2022). The seasonal cycle of water produces on comet 67P, and very likely on comets in general, a bluing of the colors approaching perihelion, when the intense cometary activity erodes the nucleus exposing the underlying water ice enriched material. Models predict that between 20 and 70% of dust removed from the southern hemisphere (Keller et al. 2017; Hu et al. 2017; Fulle et al. 2021), which is illuminated for a short time but with high intensity during the perihelion passage, falls back in the northern hemisphere mainly in decimeter-sized dehydrated aggregate, relatively poor in water ice (Keller et al. 2017). For this reason, when the cometary activity decreases, the colors of the comet go back to red values because of the dehydrated dust blanket covering the nucleus. Already at 2.2 au outbound colors were reported to be as red as in pre-perihelion observations (Fornasier et al. 2016). A substantial bluing of comet 67P surface started since February 2015, according to Ciarniello et al. (2022). These authors predict that WEBs start to be exposed mostly from this period. Figure 18 shows that very few BS were observed in the February–May 2015 period. Most have short lifetimes of a few minutes to hours, and are very likely frosts. Instead, the two 1500 m2 bright patches observed close to the Anhur–Bes boundaries (BS 143 and BS 144 in Table A.1) have a spectral slope close to zero, a water ice abundance up to 30%, and a lifetime of 1–2 weeks, and thus are very probably primitive WEB exposures. Long-lived BS show variable spectral slope values, with negative values observed predominantly in the post-perihelion period, as already described in Sect. 5.

To further constrain whether the BS exposed after February 2015 (inbound phase) and up to the post-perihelion phase can be interpreted in terms of primordial WEB exposures, the measured BS color (spectral slope) can be directly compared with the expected value for WEBs, as computed by spectral modeling following an approach similar to Ciarniello et al. (2022), based on Hapke’s theory (Hapke 2012)1. WEBs are modeled as intimate mixtures of water ice and 67P dark terrain. For the dark terrain, we assume the optical properties (single scattering albedo and single particle phase function) estimated by Ciarniello et al. (2015) from VIRTIS observations (see also Ciarniello et al. 2022 for further details), while the single scattering albedo of water ice is computed from optical constants. To constrain the color, we compute the spectral slope from the simulated spectra as a function of the ratio of the reflectance factor at 882 nm and 535 nm, following the definition given in Sect. 2. VIRTIS data are affected by a calibration residue in the 800–1000 nm spectral interval (Ciarniello et al. 2015; Filacchione et al. 2016a), thus the estimated reflectance of the dark terrain at those wavelengths is likely underestimated. This would result in a slight underestimation of the modeled spectral slope, as it involves the computation of the reflectance factor at 882 nm. To bracket the correct spectral slope value, we then perform the computation also in the reduced interval 535–800 nm, where the dark terrain reflectance from VIRTIS observations is more reliable. Given the overall shape of the 67P spectrum, with a steeper slope shortward of 800 nm, this spectral slope value likely overestimates the expected value over the 535–882 nm interval.

In Table 3, we report the estimated spectral slope of different intimate mixtures of water ice and dark terrain with variable abundances. For completeness, we run simulations varying the water ice grain size (10, 100, 1000 µm), assuming this is the same as the dark material. In the model of Ciarniello et al. (2022) (hereafter WEB model), BS from WEB exposures are assumed to be characterized by a dust-to-ice mass ratio δ ≈ 2 (O’Rourke et al. 2020), and thus their predicted spectral slope would be on the order of ≈1.5–2.4%/(100 nm) (Table 3). This is consistent with the average BS color during the post-perihelion phase and with the spectral slope of Anhur-Bes bright patches (≈0.7–1.6%/(100nm)) observed in April 2015, supporting the idea that such BS are exposures of original WEBs. The average BS color during the pre-perihelion phase, from observations carried out mostly before February 2015 (≈4.9 ± 0.3%/(100nm)), is consistent with BS having δ ≈ 8, and thus likely does not represent the exposure of WEBs. The pre-perihelion BS color could also be possibly explained by BS undergoing partial self-cleaning (see Sect. 9).

It is difficult, however, to interpret the spectral slope in terms of dust-to-ice ratio because it depends on the spatial resolution, on the intrinsic water ice abundance in mixtures, on the waterice grain size, and also on the spectral phase reddening. This last (i.e., the phenomenon related to the increase in colors and spectral slope values with increasing phase angle) is very important and is well characterized for the dark terrain (Fornasier et al. 2015), but has not been determined for the bright spots because of their limited phase angle coverage coupled with their relatively short visibility.

thumbnail Fig. 14

Duration of BS 32 in Anuket region.

thumbnail Fig. 15

Duration of BS 16 in Anuket region.

thumbnail Fig. 16

Duration of BS 36 in Neith region.

thumbnail Fig. 17

Duration of BS 204 in Atum region.

thumbnail Fig. 18

Spectral slope value, evaluated in the 535–882 nm wavelength range, vs. BS duration.

Table 3

Bright spot spectral slope (535–882 nm) as a function of the volumetric abundance and grain size of water ice mixed with 67P dark terrain; delta represents the dust-to-ice mass ratio.

9 Discussion

Figure 19 shows the BS duration versus their surface (Table A.1). It clearly demonstrates that the BS lifetimes are inconsistent with an explanation of all the observed BS in terms of diurnal or seasonal frost. The BS lifetimes do not depend on the heliocentric distance rh: the largest range of lifetimes approaches six orders of magnitudes for BS all observed at rh > 3 au, thus evidencing that the main factor driving the BS lifetime is its thickness. Since water-ice frost can be observed if its thickness is T ≥ 1 µm, the BS with the longest lifetimes would require T ≈ 1 m, inconsistent with frost build-up by diurnal or seasonal thermal inversion below the nucleus surface (De Sanctis et al. 2015; Fornasier et al. 2016) or due to the condensation on the nucleus of water vapor from the coma.

Bright features due to diurnal frost have a lifetime shorter than the nucleus rotation period (12.4 h, Mottola et al. 2014) and form the lower cluster in Fig. 19, which is independent of the BS area, because diurnal frost occurs in sunrise nucleus shadows along the terminator, which may have any area. For instance, Fornasier et al. (2016) analyze color sequences close to the perihelion passage, and find frost fronts moving with the cometary morning shadows that have a permanence time of a few minutes and a solid ice equivalent thickness of 1.5–2.7 mm.

BS lasting more than 2 days form the upper cluster in Fig. 19, and are consistent with the WEB model (Ciarniello et al. 2022) (i.e., the proposed nucleus structural model is consistent with the observed seasonal evolution of the 67P nucleus color), which is based on the only available activity model consistent with dust ejection (Fulle et al. 2020). WEBs are blocks of sizes ranging from 0.5 to 1 m and composed of water-rich pebbles, embedded in a matrix of water-poor pebbles (Ciarniello et al. 2022). Water-rich pebbles have a refractory-to-water-ice mass ratio δ ≈ 2 (O’Rourke et al. 2020), whereas water-poor pebbles have δ ≈ 50 (Fulle et al. 2021). Water-rich pebbles probably formed at the water-snow line of the protoplanetary disk due to the recondensation of water vapor (Fulle et al. 2021), so that their water ice is crystalline and may reach millimeter sizes, consistent with the extreme blue color of some BS (Fig. 8). Such millimeter-sized water ice is always embedded inside the dust particles composing all the centimeter-sized pebbles structuring every comet nucleus (Blum et al. 2017). Such a pebble structure is a necessary condition driving cometary activity (Fulle et al. 2020).

At rh = 3 au, the WEB model predicts a δ-independent nucleus erosion E = 5 cm day−1 (Fulle et al. 2020), best fitting the observed erosion during August–November 2014 in the Hapi deposits (Cambianica et al. 2020), which are composed of water-poor pebbles (Fulle et al. 2021). Therefore, the model predicts the same erosion also for all the BS lasting more than a few days listed in Table A.1, which are composed of water-rich pebbles. The line in Fig. 19 shows the observed BS lifetimes longer than 2 days compared to the BS lifetime predicted by the WEB model according to the computed erosion E, assuming that each BS is a WEB exposed by the ongoing nucleus erosion and of thickness T=A$T = \sqrt A $, where A is the measured BS area (Table A.1). The predicted BS lifetimes best fit the available data for 0.1 ≤ A ≤ 10 m2, with observed lifetimes longer and shorter by up of a factor ten than predicted, due to the fact that WEBs may have aspect ratios ranging from 0.1 to 10. WEBs of T<A$T &gt; \sqrt A $ have lifetimes shorter than predicted, the opposite if T>A$T > \sqrt A $. On average, BS observed at rh < 3 au show a shorter lifetime (Fig. 19), as predicted by the WEB model: E = 13 cm day−1 at 67P perihelion (Fulle et al. 2020).

For A > 10 m2, most BS have lifetimes shorter than predicted (Fig. 19); we recall that the duration reported here is a lower limit. This suggests that BS of A > 10 m2 are in fact clusters of BS of A ≤ 10 m2 that appear as a single larger BS, due to the low resolution of the OSIRIS images. Most BS of A > 10 m2 have been observed at rh < 3 au (Fig. 19), when Rosetta was farther from the nucleus than for rh > 3 au. We can conclude that the data shown in Fig. 19 constrain the average WEB cross section to A ≤ 10 m2. The cross section distribution of all BS of A ≤ 10 m2 shown in Fig. 12 constrains the median WEB cross section to 0.7 m2, perfectly matching the WEB model (Ciarniello et al. 2022). A total BS area covering 0.1% of the nucleus surface is lower than the uncertainty affecting the BS aerial fraction consistent with the seasonal evolution of the nucleus color (Ciarniello et al. 2022). This is evidence that the observed BS are a few percent of the exposed WEBs driving the seasonal evolution of the nucleus color. This may be related to the fact that most OSIRIS images have a resolution insufficient to resolve exposed sub-meter WEBs, and/or related to the criteria here adopted to identify BS (i.e., having lower spectral slope values than the cometary dark terrain and higher reflectance). The average cometary dark terrain is a strong darkening agent in mixtures with ice, decreasing rapidly the BS reflectance. Thus, blue areas that are water ice enriched are often not as bright as the BS here reported (i.e., 50% brighter than the cometary dark terrain), and are therefore not included in this study.

Bright spots with lifetimes between 0.5 and 2 days form another area-independent cluster, which cannot be diurnal frost of lifetime longer than the nucleus rotation period. The most probable explanation of this cluster is the uniform fallout lasting at least the first three months of 2016, explaining the evolution of the dust backscattering observed in the 67P coma (Bertini et al. 2019), and most post-perihelion nucleus reddening (Ciarniello et al. 2022). Such uniform fallout is composed of water-poor decimeter-sized chunks eroded by the CO2-driven activity from the southern hemisphere during three months centered on the 67P perihelion (Fulle et al. 2021). If these chunks fall back on BS where the self-cleaning is occurring (Pajola et al. 2017b), then they are eroded into sub-chunk dust in about 2 days at the computed erosion rate of 5 cm day−1 (Fulle et al. 2020).

The self-cleaning BS are never buried by the fallout because the southern area is eroded at a perihelion rate of about 0.1 m day−1 and because the fallout is distributed over an area five times larger than the southern eroded area (Fulle et al. 2021), so that the first chunk monolayer would occur in about 5 days. On the other hand, most sub-chunk dust eroded by the self-cleaning falls back on non-self-cleaning BS, forming a dust monolayer in about 1 day (thus matching the average cluster lifetime shown in Fig. 19) because self-cleaning and non-self-cleaning areas are similar (Fig. 19) and because the largest falling back dust is half the size of the parent chunks (Fulle et al. 2020). Close observations of dust deposits show that the self-cleaning areas are uniformly distributed among the non-self-cleaning areas, and that the non-self-cleaning deposits are much richer in sub-chunk dust than the self-cleaning deposits (Pajola et al. 2017b). Non-self-cleaning deposits may be composed of chunks of δ > 104 (Ciarniello et al. 2022), thus preventing any water-driven erosion of the chunks themselves (Fulle et al. 2020). All of this explains why exposed WEBs form two different clusters in Fig. 19: the area-independent cluster which is fallout-driven and having a shorter lifetime than the area-dependent cluster, which is erosion driven.

thumbnail Fig. 19

BS duration vs. area distribution. The heliocentric distance of the 67P comet (rh) is indicated using symbols with different colors. The black line represents the expected lifetime of bright spots predicted by the WEB model, assuming an erosion rate of 5 cm day−1.

10 Conclusions

We built the most extensive catalog of exposures of volatiles on comet 67P, based on OSIRIS-NAC color images, including 603 individual entries. We investigate the type, the spectral slope, the area distribution, and their minimum duration in a homogeneous way. Our main findings are the following:

  1. Bright spots are found isolated or in clusters, with lifetimes ranging from a few minutes-hours, in which case they are very likely frost, to several days-months, in which case they should be considered exposure of original water ice enriched blocks.

  2. Bright spots are more often observed post-perihelion and have typical sub-meter sizes, with a median value of 0.7 m2, indicating that high resolution is mandatory to observe icy exposures in comets.

  3. The BS spectral slope evolved toward negative values in post-perihelion observations, indicating the presence of frost, for the short-lived ones, and of water ice with large grain sizes (≥1000 µm) for those having longer duration.

  4. The BS lifetimes form three clusters (Fig. 19): (a) the area-independent cluster lasting less than 0.5 days, best explained by diurnal frost; (b) the area-independent cluster persisting from 0.5 to 2 days, best explained by a seasonal fallout lasting many months (Bertini et al. 2019, Ciarniello et al. 2022); and (c) the area-dependent cluster lasting more than 2 days, best explained by water-driven erosion of WEBs (Fulle et al. 2020, Ciarniello et al. 2022).

  5. The observed erosion of BS lasting more than 2 days and of δ ~ 2 matches the one observed in Hapi’s deposits of δ ~ 100 (Cambianica et al. 2020). The erosion is therefore independent of the refractory-to-ice ratio, as predicted by the water-driven activity model (Fulle et al. 2020).

  6. The observed BS lifetimes longer than 2 days are consistent with the predictions of the WEB model (Ciarniello et al. 2022).

  7. The observed median of the BS cross section (0.7 m2) fits that of the WEBs constrained by the seasonal evolution of the nucleus color (Ciarniello et al. 2022).

  8. The observed evolution of the BS spectral slope matches the predictions of the WEB model. After February 2015, when the exposition of WEBs starts to make the average nucleus color bluer (Ciarniello et al. 2022; Fornasier et al. 2016), it has a value consistent with WEBs δ ≈ 2. Before this date it shows a redder color, probably due to partial self-cleaning of the few exposed WEBs.

  9. The total integrated surface of volatile exposure is less than 50000 m2, which is less than 0.1% of the 67P nucleus surface, indicating that the top layer of cometary nucleus composition (and the whole nucleus, according to the WEB model) is dominated by refractory material.

  10. The majority of the BSs are observed in the large lobe of the comet, where they occupied 0.13% of the surface, while the small lobe has only 0.02% of volatile exposures. This finding supports the hypothesis of Massironi et al. (2015) on the binary structure of the comet, and the findings of El-Maarry et al. (2016) and Fornasier et al. (2021), on the fact that the small and large lobes of 67P comet have different mechanical and physical properties. Our study on volatile exposures highlights that the small lobe is pauperized in volatile abundance, at least in its upper layer, compared to the large lobe, even though it shows a seasonal evolution of colors (Ciarniello et al. 2022), driven by the exposure of WEBs, similar to the large lobe.

Acknowledgements

OSIRIS was built by a consortium led by the Max-PlanckInstitut für Sonnensystemforschung, Goettingen, Germany, in collaboration with CISAS, University of Padova, Italy, the Laboratoire d’Astrophysique de Marseille, France, the Instituto de Astrofísica de Andalucia, CSIC, Granada, Spain, the Scientific Support Office of the European Space Agency, Noordwijk, The Netherlands, the Instituto Nacional de Técnica Aeroespacial, Madrid, Spain, the Universidad Politéchnica de Madrid, Spain, the Department of Physics and Astronomy of Uppsala University, Sweden, and the Institut für Datentechnik und Kommunikationsnetze der Technischen Universitat Braunschweig, Germany. The support of the national funding agencies of Germany (DLR), France (CNES), Italy (ASI), Spain (MEC), Sweden (SNSB), and the ESA Technical Directorate is gratefully acknowledged. We thank the Rosetta Science Ground Segment at ESAC, the Rosetta Mission Operations Centre at ESOC and the Rosetta Project at ESTEC for their outstanding work enabling the science return of the Rosetta Mission. We acknowledge the financial support from the France Agence Nationale de la Recherche (program Classy, ANR-17-CE31-0004). We thank the anonymous reviewer for his/her comments and suggestions which helped us to improve this article.

Appendix A Supplementary material: Table

Table A.1

List of the volatile exposure positions, types, and characteristics (slope, surface, duration) identified in this paper and in the literature.

References

  1. Agarwal, J., Della Corte, V., Feldman, P. D., et al. 2017, MNRAS, 469, S606 [Google Scholar]
  2. Barucci, M. A., Filacchione, G., Fornasier, S., et al. 2016, A&A, 595, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bertin, E. 2012, ASP Conference Series, 461, 263 [NASA ADS] [Google Scholar]
  4. Bertini, I., La Forgia, F., Fulle, M., et al. 2019, MNRAS, 482, 2924 [Google Scholar]
  5. Bibring, J. P., Langevin, Y., Carter, J., et al. 2015, Science, 349, aab0671 [CrossRef] [Google Scholar]
  6. Biele, J., Ulamec, S., Maibaum, M., et al. 2015, Science, 349, aaa9816 [Google Scholar]
  7. Blum, J., Gundlach, B., Krause, M., et al. 2017, MNRAS, 469, S755 [Google Scholar]
  8. Cambianica, P., Fulle, M., Cremonese, G., et al. 2020, A&A, 636, A91 [EDP Sciences] [Google Scholar]
  9. Capaccioni, F., Coradini, A., Filacchione, G., et al. 2015, Science, 347, a0628 [NASA ADS] [CrossRef] [Google Scholar]
  10. Ciarletti, V., Herique, A., Lasue, J., et al. 2017, MNRAS, 469, S805 [NASA ADS] [CrossRef] [Google Scholar]
  11. Ciarniello, M., Capaccioni, F., Filacchione, G., et al. 2015, A&A, 583, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Ciarniello, M., Raponi, A., Capaccioni, F., et al. 2016, MNRAS, 462, S443 [Google Scholar]
  13. Ciarniello, M., Fulle, M., Raponi, A., et al. 2022, Nat. Astron., 6, 546 [NASA ADS] [CrossRef] [Google Scholar]
  14. Davidsson, B.J.R., Sierks, H., Güttler, C., et al. 2016, A&A, 592, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. De Sanctis, M.C., Capaccioni, F., Ciarniello, M., et al. 2015, Nature, 525, 500 [NASA ADS] [CrossRef] [Google Scholar]
  16. Deshapriya, J. D. P., Barucci, M. A., Fornasier, S., et al. 2016, MNRAS, 462, S274 [NASA ADS] [CrossRef] [Google Scholar]
  17. Deshapriya, J. D. P., Barucci, M. A., Fornasier, S., et al. 2018, A&A, 613, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. El-Maarry, M. R., Thomas, N., Giacomini, L., et al. 2015, A&A, 583, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. El-Maarry, M. R., Thomas, N., Gracia-Berna, A., et al. 2016, A&A, 593, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Filacchione, G., De Sanctis, M.C., Capaccioni, F., et al. 2016a, Nature, 529, 368 [NASA ADS] [CrossRef] [Google Scholar]
  21. Filacchione, G., Raponi, A., Capaccioni, F., et al. 2016b, Science, 354, 1563 [Google Scholar]
  22. Filacchione, G., Capaccioni, F., Ciarniello, M., et al. 2016c, Icarus, 274, 334 [NASA ADS] [CrossRef] [Google Scholar]
  23. Fornasier, S., Hasselmann, P. H., Barucci, M. A., et al. 2015, A&A, 583, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Fornasier, S., Mottola, S., Keller, H. U., et al. 2016, Science, 354, 1566 [NASA ADS] [CrossRef] [Google Scholar]
  25. Fornasier, S., Feller, C., Lee, J.-C., et al. 2017, MNRAS, 469, S93 [NASA ADS] [CrossRef] [Google Scholar]
  26. Fornasier, S., Feller, C., Hasselmann, P. H., et al. 2019a, A&A, 630, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Fornasier, S., Hoang, V.H., Hasselmann, P. H., et al. 2019b, A&A, 630, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Fornasier, S., Bourdelle de Micas, J., Hasselmann, P. H., et al. 2021, A&A, 653, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Fulle, M. 2021, MNRAS, 505, 3112 [Google Scholar]
  30. Fulle, M., Blum, J., Rotundi, A., et al. 2020, MNRAS, 493, 4039 [CrossRef] [Google Scholar]
  31. Groussin, O., Sierks, H., Barbieri, C., et al. 2015, A&A 583, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Hapke, B. 2002, Icarus, 157, 523 [NASA ADS] [CrossRef] [Google Scholar]
  33. Hapke, B. 2012, Icarus, 221, 1079 [NASA ADS] [CrossRef] [Google Scholar]
  34. Hasselmann, P. H., Barucci, M. A., Fornasier, S., et al. 2019, A&A, 630, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Hoang, H.V., Fornasier, S., Quirico, E., et al. 2020, MNRAS, 498, 1221 [NASA ADS] [CrossRef] [Google Scholar]
  36. Hu, X., Shi, X., Sierks, H., et al. 2017, MNRAS, 469, S295 [NASA ADS] [CrossRef] [Google Scholar]
  37. Jorda, L., Gaskell, R., Capanna, C., et al. 2016, Icarus, 277, 257 [Google Scholar]
  38. Keller, H. U., Barbieri, C., Lamy, P., et al. 2007, Space Sci. Rev., 128, 433 [Google Scholar]
  39. Keller H. U., Mottola, S., Hviid, S. F., et al. 2017, MNRAS, 469, S357 [NASA ADS] [CrossRef] [Google Scholar]
  40. La Forgia, F., Giacomini, L., Lazzarin, M., et al. 2015, A&A, 583, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Lucchetti A., Pajola, M., Fornasier, S., et al. 2017, MNRAS, 469, S238 [NASA ADS] [CrossRef] [Google Scholar]
  42. Massironi, M., Simioni, E., Marzari, F., et al. 2015, Nature, 526, 402 [NASA ADS] [CrossRef] [Google Scholar]
  43. Mennella, A., Ciarniello, M., Raponi, A., et al. 2020, ApJ, 897, L37 [NASA ADS] [CrossRef] [Google Scholar]
  44. Mottola, S., Lowry, S., Snodgrass, C., et al. 2014, A&A, 569, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Oklay, N., Mottola, S., Vincent, J.-B., et al. 2017, MNRAS, 469, S582 [NASA ADS] [CrossRef] [Google Scholar]
  46. O’Rourke, L., Heinisch, P., Blum, J., et al. 2020, Nature, 586, 697 [Google Scholar]
  47. Pajola, M., Hofner, S., Vincent, J. B., et al. 2017a, Nat. Astron., 1, 0092 [NASA ADS] [CrossRef] [Google Scholar]
  48. Pajola, M., Lucchetti, A., Fulle, M., et al. 2017b, MNRAS, 469, S636 [Google Scholar]
  49. Poch, O., Istiqomah, I., Quirico, E., et al. 2020, Science, 367, aaw7462 [NASA ADS] [CrossRef] [Google Scholar]
  50. Pommerol, A., Thomas, N., El-Maarry, M. R., et al. 2015, A&A, 583, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Protopapa, S., Sunshine, J. M., Feaga, L. M., et al. 2014, Icarus, 238, 191 [NASA ADS] [CrossRef] [Google Scholar]
  52. Quirico, E., Moroz, L. V., Schmitt, B., et al. 2016, Icarus, 272, 32 [CrossRef] [Google Scholar]
  53. Raponi, A., Ciarniello, M., Capaccioni, F., et al. 2016, MNRAS, 462, S476 [Google Scholar]
  54. Raponi, A., Ciarniello, M., Capaccioni, F., et al. 2020, Nat. Astron., 4, 500 [NASA ADS] [CrossRef] [Google Scholar]
  55. Sierks, H., Barbieri, C., Lamy, P. L., et al. 2015, Science, 347, aaa1044 [Google Scholar]
  56. Sunshine, J. M., A’Hearn, M. F., Groussin, O., et al. 2006, Science, 311, 1453 [NASA ADS] [CrossRef] [Google Scholar]
  57. Sunshine, J. M., Groussin, O., Schultz, P. H., et al. 2007, Icarus, 190, 284 [NASA ADS] [CrossRef] [Google Scholar]
  58. Thomas, N., Sierks, H., Barbieri, C., et al. 2015, Science, 337, a0440 [Google Scholar]
  59. Thomas, N., El Maarry, M. R., Theologou, P., et al. 2018, Planet. Space Sci., 164, 19 [NASA ADS] [CrossRef] [Google Scholar]
  60. Tubiana, C., Güttler, C., Kovacs, G., et al. 2015, A&A, 583, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Vincent, J.-B., Bodewits, D., Besse, S., et al. 2015, Nature, 523, 63 [Google Scholar]
  62. Vincent, J.-B., A’Hearn, M. F., Lin, Z.-Y., et al. 2016, MNRAS, 462, S184 [Google Scholar]
  63. Warren, S. G., & Brandt, R. E. 2008, J. Geophys. Res., 113, D14220 [NASA ADS] [CrossRef] [Google Scholar]
  64. Yang, B., Jewitt, D., & Bus, S. J. 2009, AJ, 137, 4538 [NASA ADS] [CrossRef] [Google Scholar]

1

In Ciarniello et al. (2022) the modeled quantity from which the color is inferred is the “effective single scattering albedo,” a proxy of the reflectance factor, while in the present work our computations refer to the reflectance factor itself, to be directly compared with OSIRIS observations.

All Tables

Table 1

Volatile exposures types from the catalog here presented (Table A.1) following the Deshapriya et al. (2018) classification scheme.

Table 2

Bright spots per region and spectral slope values.

Table 3

Bright spot spectral slope (535–882 nm) as a function of the volumetric abundance and grain size of water ice mixed with 67P dark terrain; delta represents the dust-to-ice mass ratio.

Table A.1

List of the volatile exposure positions, types, and characteristics (slope, surface, duration) identified in this paper and in the literature.

All Figures

thumbnail Fig. 1

Maps of ice exposure on comet 67P. The color-coding is as follows: red, pre-perihelion (August 2014–May 2015); cyan, perihelion (June–October 2015); green, post-perihelion (November 2015–September 2016). The spectrally blue spots, those having negative spectral slope (<−3%/(100 nm)) in the 535–882nm range, are shown in light green (post-perihelion) and orange (pre-perihelion). The symbol size represents three ranges of volatile exposure area; they are not in scale compared to the nucleus surface (51.74 km2 in total; Thomas et al. 2018), but are enlarged for clarity.

In the text
thumbnail Fig. 2

Example of isolated bright features on smooth terrains (top) and close to irregular structures (bottom), feature types 1 and 2, respectively, following the Deshapriya et al. (2018) classification scheme. The bright spot numbers correspond to those listed in Table A.1.

In the text
thumbnail Fig. 3

Example of bright features resting on boulders (BS 32 in Anuket, BS 374 in Seth, and BS 542 and BS 544 in Anhur), type 3 following Deshapriya et al. (2018) classification scheme. The bright spot numbers correspond to those listed in Table A.1. Several BS show blue colors, indicating a small or negative spectral slope value (see Table A.1).

In the text
thumbnail Fig. 4

Example of clusters of bright spots, type 4 following the Deshapriya et al. (2018) classification scheme. The bright spot numbers correspond to those listed in Table A.1. Several BS show blue colors, indicating a small or negative spectral slope value (see Table A.1).

In the text
thumbnail Fig. 5

Frequency of the bright spots vs. latitude. The histogram in black represents the pre-perihelion BS.

In the text
thumbnail Fig. 6

Histograms showing the bright spot spectral slopes, evaluated in the 535–882 nm wavelength range, pre-perihelion (August 2014–May 2015), during perihelion (June–October 2015), and post-perihelion (November 2015–September 2016).

In the text
thumbnail Fig. 7

Histograms showing the bright spots spectral slopes, evaluated in the 535–882 nm wavelength range, in images having similar high resolution (in the 0.33–1 m px−1 range), in the pre-perihelion (August 2014–May 2015), and post-perihelion (November 2015–September 2016) periods.

In the text
thumbnail Fig. 8

Duration, spectrophotometry and compositional modeling of some blue bright spots in the Anhur region. Top: RGB images showing the blue BS observed in boulders named A, B, and C on 16 July 2016 (top central panel). This area is located close to the canyon-like structure of Anhur (Fornasier et al. 2017). Frosts and tiny blue spots were already present on 25 July 2016 images (top left panel) near the boulders named A and C, and a few close to boulder A were still visible in the 30 July 2016 images (top right panel). Bottom: relative reflectance and I/F of four selected blue BS located in the aforementioned A and B boulders: BS 486 (red star, slope = −8.0%/(100nm)), BS 487 (blue asterisk, slope = −4.3%/(100nm)), BS 488 (green triangle, slope = −1.8%/(100 nm)), and BS 489 (magenta square, slope = 4.6%/(100nm)) in Table A.1. DT indicates the cometary dark terrain. Bottom right: linear mixing models of the cometary dark terrain and water ice with different grain sizes reproducing the spectral behavior of the blue bright features. The estimated water ice abundance is indicated for the different models.

In the text
thumbnail Fig. 9

Examples of RGB images showing frosts during the outbound orbit of comet 67P, near 3 au.

In the text
thumbnail Fig. 10

Map with the sources of activity identified on the 67P nucleus in Vincent et al. (2016), blue asterisks, and Fornasier et al. (2019b).

In the text
thumbnail Fig. 11

Histogram density of the bright spot diameters identified on the 67P nucleus and reported in this study. Only those having diameter smaller than 80 m are represented, for clarity. The cumulative probability is shown in blue.

In the text
thumbnail Fig. 12

Frequency of bright spot areas identified on the 67P nucleus for the pre-perihelion (top) and the post-perihelion (bottom) periods, observed at spatial resolution lower than 1 m px−1.

In the text
thumbnail Fig. 13

Frequency of the bright spot areas identified on the 67P nucleus for the pre-perihelion (top) and the post-perihelion (bottom) periods, and observed at comparable high resolution (0.33–1 m px−1).

In the text
thumbnail Fig. 14

Duration of BS 32 in Anuket region.

In the text
thumbnail Fig. 15

Duration of BS 16 in Anuket region.

In the text
thumbnail Fig. 16

Duration of BS 36 in Neith region.

In the text
thumbnail Fig. 17

Duration of BS 204 in Atum region.

In the text
thumbnail Fig. 18

Spectral slope value, evaluated in the 535–882 nm wavelength range, vs. BS duration.

In the text
thumbnail Fig. 19

BS duration vs. area distribution. The heliocentric distance of the 67P comet (rh) is indicated using symbols with different colors. The black line represents the expected lifetime of bright spots predicted by the WEB model, assuming an erosion rate of 5 cm day−1.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.