Open Access
Issue
A&A
Volume 653, September 2021
Article Number A132
Number of page(s) 14
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202141014
Published online 22 September 2021

© S. Fornasier et al. 2021

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.

1 Introduction

The Rosetta mission of the European Space Agency carried out an in-depth investigation of the short-period comet 67P/Churyumov-Gerasimenko (hereafter 67P) from 1 August 2014 until 30 September 2016, when it ended its operations landing on the comet surface. The Philae lander on board Rosetta was delivered on 12 November 2014 and, after the unexpected rebound on the Agilkia landing site (a second landing at 30 m away from the final stop that was recently identified; see O’Rourke et al. 2020), it finally landed on Abydos, where most of the foreseen in situ scientific measurements were performed.

The 67P comet has a peculiar bilobate shape with a surface characterized by extensive layering and showing a variety of terrains including both consolidated and smooth areas, pits, depressions, boulders, and talus (Sierks et al. 2015; Thomas et al. 2015; Vincent et al. 2015; Massironi et al. 2015), of which 26 regions have been defined by structural or morphological boundaries (El-Maarry et al. 2015, 2016).

Spectroscopy of comet 67P from the Visible, InfraRed, and Thermal Imaging Spectrometer (VIRTIS) reveals a red spectrum featureless until ~ 2.9 μm, then showing a broad absorption in the 2.9–3.6 μm region, indicating a composition dominated by refractory materials rich in organics (Capaccioni et al. 2015; Quirico et al. 2016). Recent recalibration of the VIRTIS instrument permits the identification of some faint absorption bands inside the broad 2.9–3.6 μm feature, as well as bands that reveal the presence of aliphatic organics (Raponi et al. 2020), and of ammonium salts (Poch et al. 2020).

The southern hemisphere regions of the comet became observable by Rosetta only since February 2015, that is about six months before the perihelion passage of the comet. For safety reasons, the spacecraft was kept at relatively high altitude from the nucleus (>100 km), permitting only low-medium spatial resolution (i.e. >2 m px−1) investigationof the southern hemisphere regions. These regions were observed at high spatial resolution (sub-meter) only during the Rosetta extended phase in 2016, when the spacecraft approached again to the nucleus as done for the northern hemisphere in 2014 during the Philae landing site selection and characterization (Sierks et al. 2015; Thomas et al. 2015).

The data acquired in 2016 permitted a detailed investigation of the physical properties and of the surface evolution due to activity for some southern hemisphere regions. The detailed analysis of Anhur and Khonsu, two regions of the big lobe of the comet, were reported in Fornasier et al. (2017, 2019a) and Hasselmann et al. (2019), respectively, while that of the Abydos landing site and surroundings was recently published in Hoang et al. (2020).

In this paper, we aim at investigating Wosret, a southern hemisphere region located in the small lobe of the nucleus (Fig. 1). As the other southern hemisphere regions, Wosret experiences intense thermal changes during the orbit because it isilluminated close to the perihelion, when the incoming solar heating on the comet reaches its peak (Marshall et al. 2017).

Wosret is a region with peculiar surface morphology: it is dominated by consolidated outcropping material showing lineaments and long fractures, and it looks flattened-out probably because of the high erosion rate associated with the strong thermal effects (El-Maarry et al. 2016). Moreover, it is an active region source of several jets and it includes the second and the last (Abydos) landing site of Philae. Therefore, the investigation of Wosret provides the context of the physical properties measured for the landing site. Finally, covering most of the southern part of the small lobe and being highly eroded, Wosret exposes the inner layers of the small lobe (Penasa et al. 2017). Thus, the comparison of its spectrophotometric properties with those of southern hemisphere regions in the big lobe affected by the same high erosion rate allows for the investigation of surface colors and possible compositional variation between the two lobes.

The article is organized as follows: Sect. 2 summarizes the OSIRIS observations of Wosret and the methodology applied in our analysis. Section 3 describes the geomorphology of the region and the surface morphological changes. In Sect. 4, we present the analysis of Wosret spectrophotometric properties and constraints on the water ice content in the localized bright spots identified. Section 5 focuses on the cometary activity events reported for Wosret. Finally, in Sect. 6 we discuss our findings and we compare them to the published results on comet 67P. This comparison allows us to highlights some differences in the physical properties between the small and big lobe of comet 67P.

thumbnail Fig. 1

Wosret region on comet 67P. Left: 3D view of the southern hemisphere of 67P (from an image acquired on 2 May 2015 at 10h42) and overlay regional boundaries to facilitate locating the Wosret region on the nucleus. The final landing site of Philae, Abydos, is indicated by the yellow arrow. Center: 3D view of the small lobe (always from the same May 2015 data shown on the left side) with the geomorphological map of Wosret from Lee et al. (2106) superimposed. Right: Sub-units on top identified by Thomas et al. (2018) in Wosret (Wr a, Wr b, and Wr c) and other regions (Ne, So, An, Gb, and Be refer to Neith, Sobek, Anhur, Geb, and Bes regions, respectively), with the associated comet view from an image acquired on 2 January 2016 at 17h23 (bottom-right).

2 Observations and data reduction

For the purposes of this study, we used data from the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) on board Rosetta. This imaging system is composed of a Narrow Angle Camera (NAC), which was mostly devoted to the investigation of the 67P nucleus, and of a Wide Angle Camera (WAC) which investigated the wide-field coma (Keller et al. 2007).

Wosret was observed only since February 2015 and until the end of the Rosetta mission (Table 1), except for small areas close to the Maftet, Hatmehit and Bastet boundaries that were visible also before. Because of the high activity phase of the comet and of the associated higher altitude orbits of Rosetta during most of the 2015 observations, Wosret was observed at sub-meter spatial resolution only in 2016, during the extended phase of the Rosetta mission.

The images analyzed here come from the NAC camera, which was equipped with 11 filters covering the 250–1000 nm range. We searchedin the OSIRIS archive, which includes more than 75 000 images of the comet, for all the NAC spectrophotometric sequences covering the Wosret region and having at least three filters.

We used radiance factor (also known as IF) data produced by the OSIRIS standard pipeline after the bias, flat-field, geometric distortion, and absolute flux calibration steps corrections described in Tubiana et al. (2015): (1)

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

To perform the spectrophotometric analysis, the images of a given sequence were first coregister using the F22 NAC filter (centered at 649.2 nm) as a reference and applying devoted python scripts, as done previously in similar studies of 67P comet (Hasselmann et al. 2019; Fornasier et al. 2019a). Then the Lommel–Seeliger disk function (D(i, e)) was applied to correct the illumination conditions: (2)

where μi and μe are the cosine of the solar incidence and emission angles, respectively. To retrieve the illumination conditions we used the 3D stereophotoclinometric (SPC, SHAP8 version) shape model of the 67P nucleus (Jorda et al. 2016). The RGB images were generated using the filters centered at 882 nm (R), 649 nm (G), and 480 nm (B), and the STIFF code (Bertin 2012).

To compute the spectral slope, we first normalized the data at 535 nm (filter F23), then we computed it in the 882–535 nm wavelength range to be consistent with the methodology applied for 67P spectrophotometry in previous publications. When the F23 filter was not in the observing sequence, we produce an artificial image at 535 nm through linear interpolation of the images acquired at 480 nm and 649 nm.

The average spectral slope reported in Table 1 was evaluated fitting the slope distribution with a Gaussian function, where the average spectral slope is derived from the peak value of the Gaussian and the associated error from the standard deviation of the fit. For the medium-low resolution data covering also other regions, the slope maps were first cut selecting the area covering Wosret before applying the Gaussian fit.

Details on the observing conditions and the Wosret spectral slope values are reported in Table 1.

To estimate the height h of a given surface feature, we considered its projected shadow and the associated geometric conditions, determined by ancillary images produced using an image simulator (Hasselmann et al. 2019) supplied with pre-calculated NAIF SPICE Kernels1 and the SPC shape model (Jorda et al. 2016). Any misalignment between the ancillary images and the observed one was corrected through translation and rotation transformations calculated from uniquely identified tie landmarks in both images. We thus used the following equation to determine the height (h) of a given surface feature: (3)

where Lsha is the length of the shadow of the feature, and i the incident angle. The incidence angle here considered is the angle between the vector centered in the tip of the shadow, connecting the top of the structure, and the normal vector to the surface of the feature (see Fig. 2 from Cambianica et al. 2020).

This method has been tested and applied in many analyses of meter-sized features on the OSIRIS images of the nucleus to estimate local surface changes (El-Maarry et al. 2017; Hasselmann et al. 2019; Cambianica et al. 2020).

Table 1

Observing conditions of the Wosret region for the NAC images.

thumbnail Fig. 2

Comparison between images acquired on May 2015 (on the left, spatial resolution (res) of 2.3 m px−1) and on July 2016 (on the right, res = 0.3 m px−1) covering part the Wosret region (area a following Thomas et al. 2018 sub-regions definition). A new cluster of outcrops (named A and indicated by the white arrow) is clearly visible on the July 2016 images. Some boulders are indicated for reference (and numbered 1, 2, and 3). The ellipse named B includes a region where the dust coating may also have been thinned during the perihelion passage, as observed close to the cluster of outcrops named A and to the goosebump features.

thumbnail Fig. 3

Comparison between images acquired on March 2015 (left side, res = 1.9 m px−1) and June 2016 (right side, res = 0.5 m px−1) showing theformation of a new cavity inside the ellipse named C. Some boulders are indicated for reference by white bars and numbers, following the same numbering scheme of Fig. 2.

3 Geomorphology and surface evolution

Wosret is surrounded by Maftet and Bastet in the West and East sides, respectively, by the Hatmehit depression at North, and by the Neith region at South (Fig. 1). Wosret is dominated by consolidated outcropping material showing ubiquitous fractures up to ~ 300 m in length, lineaments (El-Maarry et al. 2016), and some pits and niches observed to be periodically active during the perihelion passage (Vincent et al. 2016; Fornasier et al. 2019b). The region has a surface of about 2.36 km2 covering most of the small lobe’s southern hemisphere, and it has been divided by Thomas et al. (2018) into three areas (Wr a, Wr b, and Wr c), showing topographic and textural difference (Fig. 1). Area Wr a is relatively flat and smooth, Wr b is dominated by fractured terrains, while area Wr c shows a higher roughness texture and it includes almost circular depressions and ridges (Thomas et al. 2018). As shown in Fig. 1, Wosret lacks of wide-spread dust coatings, except for few pits filled by smooth materials (El-Maarry et al. 2016).

Wosret has a peculiar flattened-out aspect (El-Maarry et al. 2016), probably because of the high erosion rate associated with the strong thermal effects. In fact, Keller et al. (2015, 2017) estimated an erosion rate up to four times higher in the southern hemisphere than in the northern one. Concerning Wosret, Lai et al. (2016) estimated a dust loss up to 1.8 m in depth per orbit. Wosret is also the region which has the highest water production rate during perihelion, according to the Microwave Instrument for the Rosetta Orbiter (MIRO) measurements (Marschall et al. 2017).

Figures 2 and 3 show part of the Wosret region (sub-unit Wr a following Thomas et al. 2018 classification) observed in 2015, before the perihelion passage and in 2016. The images taken on March and May 2015 have not sufficient spatial resolution (2.3 m px−1) nor good illumination conditions (i ~ 30°) to accurately measure the typical depth of surface features from their shadows. Therefore, it is not possible to precisely quantify the possible evolution of the dust cover in this area through the perihelion passage.

Conversely, the NAC images acquired on June and July 2016, have a high spatial resolution (0.27–0.52 mpx−1), which highlights some interesting features such as highly fractured terrains showing polygonal block patterns (F2 structure in El-Maarry et al. 2016). The texture reminds the surface features referred to as “goosebumps” (see Fig. 2) or clods, first noticed on the 67P comet inside the Seth pits (Sierks et al. 2015), but later observed in other regions. The Wosret blocks appear consolidated, highly irregular in shape and height. We estimate the size of 251 blocks from the 2 July 2016 image, finding an equivalent diameter ranging from 2 to ~ 12 m, with an average value of 4.74 ± 1.51 m (Fig. 4). Wosret clods are larger than the average diameter values (⟨ D⟩) reported forsimilar structures observed in the big lobe of the comet, such as Seth (⟨ D⟩ = 2.2 m, Sierks et al. 2015), Imhotep and Anubis (⟨D⟩ = 3.2 m and 2.5 m, respectively, Davidsson et al. 2016), or Atum, where the average diameter is 2.2 m but, in its southern part, some clods are larger, with size ranging from 3.2 to 4.6 m (Davidsson et al. 2016).

In comparing the pre- and post-perihelion images, we also identify a couple of clear morphological changes: a cluster of bright outcrops at lat = −15.15°, and lon = 343.2° (feature A in Fig. 2); a new cavity at lat = −19.7°, and lon = 335.9° (inside the ellipse called C in Fig. 3).

thumbnail Fig. 4

Goosebumps equivalent diameter from the 2 July 2016 image (see Fig. 2).

thumbnail Fig. 5

Spectral slope evolution with heliocentric distance comparing data acquired on 2 May 2015 (rh = 1.73 AU inbound) and 2 July 2016 (rh = 3.32 AU outbound) images. The July 2016 data were phase reddening corrected from 93° to 62° using the phase reddening coefficient of 0.055 × 10−4 nm−1 deg−1 determined in Sect. 4.1 for Wosret in the post-perihelion 2016 observations. Comparing to the original data, we applied a zoom factor of 2 to the May 2015 data, while the July 2016 image was binned by a factor of 2.

3.1 Outcrops cluster

The cluster of outcrops is clearly visible in the July 2016 image (Fig. 2), and it includes 15 outcrops having 2–5 m in size. In the 2 July 2016 image acquired at 15h22, the outcrops look brighter than the surroundings. This is not, however, an intrinsic effect related to a different composition, but simply an illumination effect. In fact, OSIRIS also acquired an image of the same area one hour before the one presented in Fig. 2, and there the brightness of the cluster is similar to that of the average terrain. In the 2015 images, only one block of the cluster is visible, while most elements of the cluster are buried beneath the dust. As said before, the resolution of the Wosret pre-perihelion images is too low and our attempts to estimate the features height were unsuccessful, thus, we cannot precisely quantify the dust coating variation. However, we were able to measure the height of the outcrops with shadows clearly ending on the background terrain in the 2016 images, finding an average height of 1.6 ± 0.3 m. We can thus infer that the thickness of the dust coating in this area diminished by ~ 1 m from pre-perihelion to post-perihelion to reveal the 15 outcrops of the cluster. We visually inspected and compared the pre-perihelion images and the post-perihelion ones binned by a factor of 8 to have a comparable spatial resolution between the two datasets. The pre-perihelion images look smoother in several areas, including the one around the outcrops, showing the goosebumps features, and the area at the northern side of Wosret, indicated by the ellipse named B in Fig. 2. We thus may tentatively deduce that the dust coating was thinned in the aforementioned areas too.

Dust movements in this region are also supported by the observed colors variation reported in Fig. 5, where the spectral slope of the same pre- and post-perihelion images presented in Fig. 2 are shown. As these images were acquired at different heliocentric distances and phase angles, and considering that comet 67P has a strong spectral phase reddening effect, the July 2016 image was phase reddening corrected from phase angle 93° to 62° using the Wosret phase reddening coefficient here determined from 2016 observations (see Sect. 4.1), in order to match the viewing geometry of the May 2015 image. As indicated in Fig. 5, the region included inside the A ellipse was relatively spectrally bluer (i.e., it has a lower spectral slope) than that included in the B ellipse in pre-perihelion data. Conversely, in July 2016 data the region included inside the A ellipse looks redder or as red as the one inside the B ellipse, indicating local surface evolution related to the composition or to different grain size or roughness of the material.

3.2 New cavity

The cavity that formed between March 2015 and June 2016 has apparently expanded from another small depression visible on 25 March 2015 image. It is ~ 30 m long and 11–15 m large, having an area of 492 ± 32 m2. We traced seven shadow lengths throughout the cavity, from which we derive a depth of 6.5 ± 0.8 m. Assuming for the loss material the density of the bulk nucleus (537.8 ± 2 kg m−3 (Pätzold et al. 2016; Jorda et al. 2016; Preusker et al. 2017), we estimate that (1.2 ± 0.2) × 106 kg of cometary material was lost in this cavity. For comparison, this mass is similar to that lost in a couple of new scarps (or cliffs) observed in the Anhur region (see Figs. 2 and 10 in Fornasier et al. 2019a). However, Anhur, as well as the Khonsu region or the Aswan site, also shows some drastic surface changes with an estimated mass loss at least ten times larger than that of the cavity here reported (Fornasier et al. 2019a; Hasselmann et al. 2019).

thumbnail Fig. 6

Spectral slope of the Wosret region, evaluated in the 535–882 nm range, from different images in the April 2015–July 2016 time frame. The Wosret region is indicated by white rectangles in the lower resolution images. The high resolution image acquired on July 2016 was binned by a factor of 2.

4 Spectrophotometry and bright spots

Comet 67P shows compositional heterogeneities which have been observed both with the VIRTIS spectrometer and the OSIRIS cameras (Fornasier et al. 2016; Barucci et al. 2016; Filacchione et al. 2016b; Deshapriya et al. 2016; Oklay et al. 2017). From OSIRIS pre-perihelion data Fornasier et al. (2015) reported spectral slope variations from ~ 11 to 21%/(100 nm) in the 535–882 nm range and at phase ~ 50o, distinguishing three kind of terrains across the nucleus, from the spectrally bluer and water ice enriched terrains observed in Hapi, to the redder ones, associated mostly with dusty regions. A similar color variation has been reported later for the southern hemisphere (Fornasier et al. 2016). Looking at Table 1 and Fig. 6, Wosret is shown to have an average slope of 15.5%/(100 nm) at α = 61° (pre-perihelion), belonging thus to the intermediate color group # 2, according to Fornasier et al. (2015) spectral classification of 67P nucleus terrains.

4.1 Spectral phase reddening

The variation of the spectral slope with phase angle, knows as the spectral phase reddening effect, is the results of small-scale surface roughness and multiple scattering in the surface medium at extreme geometries (high phase angles) and it is commonly observed on several Solar System bodies. This has been reported for asteroids (see Fornasier et al. 2020 and references therein), comets (Fornasier et al. 2015; Ciarniello et al. 2015; Longobardo et al. 2017), as well as planets and their satellites (Gehrels et al. 1964; Warell & Bergfors 2008; Nelson et al. 1987; Cuzzi et al. 2002; Filacchione et al. 2012). In case of a low-albedo surface, such as that of cometary nuclei, the spectral phase reddening effect has been interpreted as due to the presence of fines, namely, particles of ~ micron size, or to the irregular surface structure of larger grains, having micron-scale surface roughness (Li et al. 2019; Schröder et al. 2014; Pilorget et al. 2016).

Comet 67P has a strong spectral phase reddening effect (Fornasier et al. 2015, 2017; Ciarniello et al. 2015; Longobardo et al. 2017), which is, moreover not fixed over time but it shows seasonal variations depending on the cometary solar distance and activity (Fornasier et al. 2016). In fact, the spectral slope value was found to be smaller close to the perihelion compared to the one measured pre- and post-perihelion at similar phase angles. This less red color was interpreted as a consequence of the cometary activity peak during perihelion, with the partial removal of the dust mantle and the consequent exposure of the volatile enriched material underneath. The 67P phase reddening coefficient is monotonic and wavelength dependent (Table 2), with values of 0.04–0.1 × 10−4 nm−1 deg−1 in the 535–882 nm range (Fornasier et al. 2015, 2016), and lower values, 0.015–0.018 × 10−4 nm−1 deg−1 in the 1–2 μm range (Ciarniello et al. 2015; Longobardo et al. 2017).

As the whole nucleus, also the Wosret region shows changes in the spectral slope with the phase angle (Figs. 6 and 7, and Table 1). From the reported spectral slopes at different phase angles, we computed the phase reddening coefficient (γ) and the spectral slope at 0° phase angle (Y0) applying a linear fit to the data. Because of the different heliocentric distances and cometary activity levels, we evaluated the spectral reddening coefficient for data acquired in 2016 and for those acquired close to the perihelion passage during the May–October 2015 timeframe (Fig. 7). Unfortunately, the Wosret region was not observable during 2014, except for tiny fractions, and no data are available at a low phase angle.

For the Wosret region, we observe a much shallower phase reddening coefficient compared to the results reported in Fornasier et al. (2015). For the post-perihelion observations of 2016, we find for Wosret a γ2016 = (0.05462 ± 0.0042) × 10−4 nm−1 deg−1, which is about a factor of 2 lower than that determined for the northern hemisphere regions in inbound orbits (Fornasier et al. 2015). For the 2015 observations acquired close to the perihelion passage the phase angle coverage is very limited, and the associated phase reddening coefficient (γperih) is (0.0396±0.0067)×10−4 nm−1 deg−1, a value veryclose to the one found globally for the comet during the perihelion passage (Table 2).

Similar trend with progressive lower phase reddening and lower spectral slope values approaching perihelion were also observed elsewhere on comet 67P: in the northern hemisphere area (at 5–14° latitude) across the Imhotep-Ash regions (γ = 0.06 × 10−4 nm−1 deg−1) during the February 2015 Rosetta flyby (Feller et al. 2016); in the southern hemisphere Anhur region (Fornasier et al. 2017, 2019b); and in the Abydos final landing site of Philae (Hoang et al. 2020), which is located in Wosret but close to the boundaries with the Bastet and Hatmehit regions. For Abydos, Hoang et al. (2020) found a phase reddening coefficient value (Table 2) very close to the one here determined for the whole Wosret region.

The phase reddening coefficient determined for the Wosret region is comparable to the one found for other dark and presumed organic rich bodied such as D-type asteroids (0.05 ± 0.03 × 10−4 nm−1 deg−1 in the 0.45–2.45 μm range, Lantz et al. 2017), and the dwarf planet Ceres (γ = 0.046 × 10−4 nm−1 deg−1, Ciarniello et al. 2017, 2020). It is however two to four times higher than the value reported for low albedo and carbonaceous rich bodies such as the near-Earth asteroids: (101955) Bennu and (162173) Ryugu (Table 2), that were recently visited by the OSIRIS-REx and Hayabusa 2 missions, respectively.

Considering that the surfaces of Bennu and Ryugu are dominated by boulders and display a lack of extended fine dust layers (Lauretta et al. 2019; Jaumann et al. 2019), as well as the fact that the 67P phase reddening decreases approaching perihelion, when the activity lifts-up part of the dust mantle, and that Wosret has a lower phase reddening coefficient than the global 67P nucleus, we may deduce that the surface structure of Wosret has an intermediate roughness level between that of dust poor surfaces, such as those of Bennu and Ryugu, and the dust-rich ones such as most of the 67P northern hemisphere regions. The lack of widespread dust coating observed in the southern hemisphere regions of the comet (El-Maarry et al. 2016; Thomas et al. 2018), including Wosret, is probably responsible of its reduced surface micro-roughness relative to most of the northern hemisphere regions of the comet.

This is corroborated by the results of Keller et al. (2015, 2017), who found that the dust particles ejected from the southern hemisphere of comet 67P during the peak of activity partially fall back and are then deposited in the northern hemisphere. Thus, the southern hemisphere regions become progressively depleted of dust.

Table 2

Phase reddening coefficients evaluated for Wosret during different observations and those reported in the literature for comet 67P and for other dark small bodies.

thumbnail Fig. 7

Phase reddening in spectral slopes evaluated in the 535–882 nm range. The black circles represent data acquired in 2016, redcircles those acquired in May–October 2015, that is three months before and after the perihelion passage, while the remaining 2015 data are represented with blue circles. The black and red lines are the linear fit of the 2016 and near perihelion data, respectively.

thumbnail Fig. 8

RGB images acquired on 17 December 2015 taken one hour apart. The right image acquired at UT23h56 shows a relatively blue region (indicated by the white rectangle) at dawn. The same area was observed in a sequence acquired one hour later (on 18 December 2015, UT 00:56) and it shows a color indistinguishable from that of the average Wosret region. The bluing at dawn is periodic and observed in other sequences in December 2015 but also in June 2015, as shown in the top insert of the image on the right side.

4.2 Water ice exposures

In our spectral analysis of the Wosret region, we notice some local compositional heterogeneity: (i) a few areas that look bluer in colors than their surroundings; (ii) some localized bright spots showing evidence of exposed volatiles.

Figure 8 shows RGB images acquired on 17 December 2015 separated by 1 h. On the right side of Fig. 8 an area (indicated by the white rectangle) emerging from the shadows shows a bluer color than its surroundings. Further images acquired 1–2 h later, when the region is fully illuminated by the Sun, reveal colors and slope as red as the surroundings. This diurnal color variation – namely, a relatively bluer slope at dawn or for an area just emerging from shadows followed later by a color as red as the surrounding terrains when the area is illuminated by the Sun – is periodic. The diurnal color variation is observed in different data acquired in December 2015, but also before and after the perihelion passage, for example on June 2015, as shown in the insert in Fig. 8, and on June 2016 (Fig. 9). Another example of bluer area is shown in Fig. 10 from June 2016 data (the region inside the blue rectangle in the left side image).

This diurnal color variation has been noted at dawn on different areas on both lobes of the comet (Fornasier et al. 2016), and associated with a surface enrichment of water ice-frost condensed during the previous night, which results in a relatively blue color at dawn. During the comet morning the increasing temperature produces the sublimation of the volatiles and, consequently, the surface gets redder in color because of depleted in water ice-frost. Previous studies made with OSIRIS and VIRTIS have in fact demonstrate that the relatively blue colors of a given area of comet 67P are associated with a local enrichment in volatiles (de Sanctis et al. 2015; Barucci et al. 2016; Fornasier et al. 2015, 2016; Filacchione et al. 2016a,b).

Beside local bluing, we also observe the presence of tiny bright spots (BS), having a surface of 1–2 m2 usually located close to shadows. For these regions of interest (ROIs), we computed the relative spectrophotometry and radiance by integrating them over a box of 3 × 3 pixels, as reported in Fig. 10. These features, besides being brighter than the comet dark terrain, have a distinct spectrophotometric behavior characterized by a low to neutral slope, indicating localized exposures of water ice. For the bright spots we attempt to estimate their water-ice content following the method described in Fornasier et al. (2019a), namely, using a simple areal mixing model with two components: the cometary dark terrain (represented by the circle in Fig. 10) and water ice (4)

where R is the reflectance of the bright patches, Rice and RDT are the reflectance of the water ice and of the cometary dark terrain, respectively, and p is the relative surface fraction of water ice. The water-ice spectrum was derived from the synthetic reflectance from Hapke modeling starting from optical constants published in Warren & Brandt (2008) and adopting a grain size of 30–100 μm, as is typical for ice grains on cometary nuclei (Sunshine et al. 2006; Capaccioni et al. 2015; Filacchione et al. 2016a), but also testing larger grain size of 1000 μm. Before applying the linear mixing model, the images were photometrically corrected using the Hapke model parameters determined by Fornasier et al. (2015, see their Table 4) to obtain an estimation of the normal albedo of the different ROIs. In doing so, we assume that the phase function, determined by Fornasier et al. (2015) from the comet surface at 649 nm also applies to the other wavelengths and to the volatile rich spots. This introduces some uncertainties in the water-ice abundance estimation, but because of the short visibility of the bright spots and of the limited phase angle coverage, it is impossible to determine Hapke photometric parameters for the Wosret region or for the BS on comet 67P.

We report in Table 3 the coordinated of the ROIs investigated, the associated spectral slope, and the estimated water ice content that was determined for BS, having an incidence and emission angle lower than 75o, in line with the illumination conditions for which the photometric corrections are reliable (Hasselmann et al. 2017).

For the 13 June 2016 data, the ROI indicated by the blue asterisk (Fig. 10) is a typical example of a relatively blue terrain just emerging from shadows, where the spectral slope is lower (~ 12 (%/100 nm)) than that of the average dark terrain, but far from the null value. We estimate a local enrichment of ~ 7% in water ice in this ROI, and the best fit is given by small water ice particles. The ROIs indicated by the red star and the green triangle have a flat spectrophotometric behavior indicating some exposure of water ice. They are close to shadows and very tiny, their bright surface being only of ~ 1 m2, that is smaller than the box (2.3 m2) used to integrate the signal. Therefore, the water-ice estimation (around 14% for the ROI indicated by the green triangle in Table 3) should be taken with caution because the integrated signal is partially affected by the shadowed regions. We show these ROIs as example of tiny BS seen on Wosret.

Of greater interest are the BS observed on 2 July 2016 at a spatial resolution that is about the double than in June 2016 (Fig. 10, right image). These BS are indicated by the magenta square and the blue asterisk, they also have a small size (~ 1.2 and 1.5 m2, respectively), but they are larger than the box where the signal was integrated and they are extremely bright, by a factor of 12 and 5, respectively, than the comet dark terrain. Their spectrophotometry is peculiar, showing a moderate negative slope in the 535–882 nm range. These ROIs are inside shadows and represent the top of morphological features illuminated only in these areas. Our best-fit compositional models for these two ROIs are shown in Fig. 11. The BS represented by the blue asterisk is best fit by a model including 29 ± 3% water ice with large grain size (1000 μm), which better mimics the decrease of reflectance at 880 nm than the model with smaller grain size. For the ROI indicated by the magenta square the estimated normal reflectance is very high, about 60–70%, as it is the associated water-ice abundance: 65.5 ± 6.0% for water ice with a grain size of 100 μm. However, these models have a higher chi-squared value than those reproducing the blue asterisk ROIs normal albedo and do not fully match the observed spectrophotometry of the bright spot. In this case, we favor the model including medium size water ice grains (100 μm) because the one with larger grains produces a strong decrease in reflectance in the NIR range, which is not observed in the data.

Although the models are not optimal, we have clear indications of a local water ice abundance beyond 50%, which is among the highest reported in the literature for bright spots on comet 67P, but in very small areas. Other tiny bright spots have been observed in Wosret in the area surrounding the Abydos landing site, including part of the Wosret, Bastet, and Hatmehit regions, which were recently investigated by Hoang et al. (2020). They found a number of bright spots with size ranging from 0.1 m2 to 27 m2, having a distinct flat spectrum compared to the comet average dark and red terrain, and with an estimated local ice abundance reaching 50–80% for 3 of them, thus indicating a fresh exposure of volatiles. A high content of water ice (~ 46%) was also determined for the boulder onto which Philae stamped two minutes during the second landing site located 30 m away from Abydos, unveiling the primordial water ice inside it (O’Rourke et al. 2020).

Other ROIs exposing fresh water ice with an estimated ice content higher than ~ 25% were found in several areas in the Anhur, Bes, Khonsu, and Imhotep regions (Fornasier et al. 2016, 2017, 2019a; Deshapriya et al. 2016; Oklay et al. 2017; Hasselmann et al. 2019), as well as in the Aswan site (Pajola et al. 2017). It should be noted that these last bright and water-enriched areas are usually much more frequently observed and much larger than the ones detected in Wosret.

thumbnail Fig. 9

RGB image taken on 12 June 2016 at 12h48 showing a blue patch at latitude − 21° and longitude 318° in an area observed to be active by Vincent et al. (2016, jet numbered 11 in their paper).

thumbnail Fig. 10

Wosret spectrophotometry. Top: examples of blue patches and bright spots from 13 June 2016 (left) and 2 July 2016 images (right). Inserts show a zoom of some ROIs. Middle and bottom panels: corresponding relative reflectance and the I over F value at the given phase angle (70° and 93° for the June and July images, respectively). The dark terrain used in the mixing model is represented by a circle in white for clarity in the RGB images on the top and in black in the plots.

Table 3

Water-ice content estimates for the bright features observed in the 13 June and 2 July 2016 images.

thumbnail Fig. 11

Normal albedo of the bright spots (indicated by the magenta squares and the blue asterisks, see Fig. 10) observed on 2 July 2016 UT15h22. The symbols correspond to the bright spots shown on the right panel of Fig. 10. The black circles represent the mean spectrum of the comet from a region close to the bright patches. Continuous and dashed lines show the best fit spectral models to the bright patches (in gray the associated uncertainty), produced by the linear mixture of the comet dark terrain enriched with 29 ± 3% of water ice with 1000 μm grain size (continuous line) for the ROI indicated by the blue asterisk, and with 65.5 ± 6.0% of water ice with 100 μm grain size (dashed line) for that represented by the magenta square. Models including water ice of medium and largegrain size are represented for comparison: dashed lines indicate models with water ice of 100 μm grain size, and continuous lines those with 1000 μm grain size.

5 Activity

Wosret is also one of the most active regions: it has the highest estimated water production rate, which peaks at ~ 1028 molecules s−1 during the perihelion passage (Marshall et al. 2017), and it is the source of about 40 activity events, including one of the brightest outburst caught by the Rosetta observations (Vincent et al. 2016). In fact, 33 distinct source regions were reported by Fornasier et al. (2019b) and six more were identified by Vincent et al. (2016). Their position is shown in Fig. 12, with black squares representing four cavities showing repeated activity close to the perihelion passage (see Fig. 18 in Fornasier et al. 2019b). These cavities were showing activity events in 16 to 60 distinct observations. The jets departing from these cavities are often relatively faint, sometimes with peculiar morphology.

On Wosret, the activity mechanism is triggered by changes in local insolation linked to the diurnal and seasonal cycles of water on the nucleus. In fact, most of the activity sources, notably the periodically active cavities, are in or close to shadowed areas, as similarly observed for the majority of active sources of comet 67P during the perihelion passage (Fornasier et al. 2019b). In areas casting shadows, subsurface water ice may easily recondense during the cometary night, and it may survive at the cometary surface as frost or ice until the area is sufficiently illuminated and heated by the Sun to produce its sublimation.

Fornasier et al. (2019b) reported the presence of tiny water ice patches in one of the active cavities, while for the others, unfortunately, no observations at high spatial resolution are available. Some images caught however few areas enriched in volatiles, as discussed in the previous section. An interesting example is reported in Fig. 9, which displays the southern part of Wosret at high resolution (310 < lon < 330°, −30, lat < −20°), including a blue patch located at latitude −21° and longitude 318°, and few tiny bright spots. This patch has a moderate spectral slope (12.7%/(100 nm)) about 50% lower than that of the average dark terrain (18.1%/(100 nm)), and it is located close to the source region of a summer jet identified by Vincent et al. (2016, jet #11). Similar processes associated with the condensation and sublimation cycle of volatiles are expected to be ongoing in the other cavities of Wosret that show periodic activity.

thumbnail Fig. 12

Map of the comet with superposed the sources of jets identified in the Wosret region during the comet southern hemisphere summer. The location of some nearby jets are averaged for clarity, and some notable jets are represented with larger symbols. Blue circles represent jets reported in Vincent et al. (2016), green triangles those reported in Fornasier et al. (2019a), and black squares represent cavities showing repeated activity in several data sets (Fornasier et al. 2019a).

6 Discussion and conclusions

Among the 26 regions of comet 67P, Wosret is shown to have unique geomorphological features. It is one of the regions receiving the highest solar flux, that is, between 400 and 600 W m−2 at perihelion (Marshall et al. 2017). This high insolation is responsible for both of the high erosion and activity level in Wosret. The high erosion level produces the flattened-out aspect, the pervasive fracturing observed in the consolidated areas in its southern part, and, together with the activity, the overall lack of dust deposits compared to other regions (El-Maarry et al. 2016). Consequently, Wosret exposes part of the inner layers of the small lobe (Penasa et al. 2017).

The final landing site of Philae, Abydos, is located in the northern part of Wosret, close to the boundary with the Hatmehit depression and Bastet region (Fig. 1). Its spectrophotometric properties have been studied by Hoang et al. (2020), who investigated an area of about 5 × 5° in longitudeand latitude surrounding the Philae final landing site. This terrain is as dark as the overall nucleus albedo and it shows relatively red colors, with the presence of some localized bright spots of 0.1 m2 to 27 m2 size rich in water ice. They also identified few possible morphological changes in this area with an estimated total mass loss of 4.7–7.0 × 105 kg.

Globally, the surface properties of Abydos are nearly indistinguishable from those of Wosret, thus, the results derived from the Philae in situ measurements may be considered representative of the whole Wosret region. The ubiquitous fractures observed in Wosret should extend also at sub-meter scale, as derived for the surrounding of Philae landing site (Bibring et al. 2015; Poulet et al. 2016). The fact that the lander instrument detected no dust impacts (Krüger et al. 2015) may support the overall paucity of dust in Wosret, even if this could also be related to shielding effects by the cliffs close to Philae combined with the relatively low level of activity in November 2014. Regarding the compressive strength of Abydos, very different numbers ranging from a few tenths of Pa to 2 MPa were reported in the literature (Spohn et al. 2015; Biele et al. 2015; Knapmeyer et al. 2018; Heinisch et al. 2019), although the measurements derived from the penetrator sensors, which give the higher strengths values, should be taken with caution because they are affected by deployment uncertainties.

A recent reconstruction of the final Philae trajectory allowed for the identification of an additional landing site touched during about 2 min before Philae stopped in Abydos. This side is located 30 m apart from Abydos and it is also located in Wosret (O’Rourke et al. 2020). During this landing, Philae collided with a boulder, producing some morphological changes, lifting up the dust covering the surface and uncovering primitive buried water ice in a 3.5 m2 size spot, six to eight times brighter than the average comet. This confirms previous findings showing that the water ice is abundant beneath the dark dust layer and drives the observed seasonal and diurnal color variations (Fornasier et al. 2016). In fact, cliff collapses observed in Seth and Anhur regions also exposed underneath layers that are rich in water ice (Pajola et al. 2017; Fornasier et al. 2019a).

O’Rourke et al. (2020) were also able to measure the local porosity of the boulder into which Philae stamped, finding high values (75%) consistent with those determined previously for the overall nucleus (Kofman et al. 2015; Herique et al. 2019; Pätzold et al. 2016). They also found an extremely low (< 12 Pascals) value for the compressive strength, indicating that the mixture of dust and ice is extremely soft (O’Rourke et al. 2020) even inside the boulders.

As shown in the previous section, Wosret spectrophotometric properties are similar to those observed in other regions of the 67P nucleus. Here, we report the presence of local bright spots and relatively blue areas, located close to shadows and associated with the diurnal cycle of water ice. These last areas are however very small and less frequent than the water ice enriched zones observed in Anhur, a region in the big lobe that is subject to the same strong thermal heating during the cometary summer and that is also highly active (Fornasier et al. 2017, 2019a).

Even if no appreciable differences between the two lobes have been reported in the literature on the global surface composition (Capaccioni et al. 2015; Fornasier et al. 2015), and on the deuterium to hydrogen (DH) ratio (Schroeder et al. 2019), our study on Wosret highlights some differences in the physical and mechanical properties among the two lobes, which are summarized in Table 4:

  1. First, the water-ice enriched regions directly exposed at the surface of Wosret are less frequent and smaller in size than those observed in Anhur, where water-ice rich areas of sizes ranging from a few to 1600 m2 were observed (Fornasier et al. 2016);

  2. We report the presence of spectrally bluer area enriched in frost located close to shadows and related to the diurnal cycle of water. However, these areas are less frequently observed in Wosret than in other regions. For example, frost is repeatedly observed in Anhur inside shadowed regions (Fornasier et al. 2019b);

  3. The comparison of Wosret surface evolution with that of southern hemisphere regions such as Anhur and Khonsu, which experience the same high incoming solar flux than Wosret at perihelion (see Fig. 6 of Marshall et al. 2017), indicates important differences. In fact, all these regions are highly active, but the activity results in different surface re-shaping. Indeed important morphological changes are observed in Anhur: the formation of two new scarps exposing the subsurface water ice layer; the local removal of the dust coating up to ~14 ± 2 m in depth within a canyon-like structure; vanishing structures of several tenth of meters in length; and boulder displacement and fragmentation, for a total mass loss estimation higher than 50 × 106 kg (Fornasier et al. 2017, 2019a). For Khonsu, another southern hemisphere region of the big lobe, Hasselmann et al. (2019) reported a number of morphological changes (boulder displacements, cavities formation, dust bank removal), producing an estimated mass loss of ~2 × 108 kg. Conversely, no major morphological changes are observed in Wosret, except for the formation of the new cavity, and the dust coating removal tentatively estimated in ~1 m depth, locally. This difference cannot be related to observational biases because these regions have been observed under similar spatial resolution and illumination conditions pre- and post-perihelion;

  4. Polygonal block in Wosret are on average two times larger than the goosebumps (or clods) observed in different regions on the big lobe. These structures have been interpreted as being representative of the original cometesimals forming, by aggregation, cometary nuclei, or as the result of fracturing processes produced by the seasonal and diurnal thermal gradients in the material (Sierks et al. 2015; El-Maarry et al. 2015; Auger et al. 2015).

El-Maarry et al. (2016), based on the morphological analysis of the 67P regions, deduced that the material composing the two lobes has different physical and mechanical properties. We can support this conclusion by the fact that the clods-goosebumps features in Wosret are larger than those observed in the big lobe. Moreover, the limited morphological changes observed in Wosret as compared to Anhur, shown in Table 4 (while both regions were among the most active ones and the most heatedby the Sun during the perihelion passage), indicate that the material on Wosret should be less fragile and more consolidatedthan the material present on the big lobe southern regions. Furthermore, the low number of exposed volatiles at the surface of Wosret, as compared to Anhur, might indicate that the small lobe has a lower volatile content, at least in its top layers, than the biglobe.

All this evidence support the hypothesis formulated by Massironi et al. (2015) based on the analysis of the layering of the two lobes stating that comet 67P is composed of two distinct bodies that merged during a low-velocity collision in the early Solar System.

Table 4

Summary of the observed behavior and physical properties for the small and big lobe of comet 67P.

Acknowledgements

OSIRIS was built by a consortium led by the Max-Planck-Institut für Sonnensystemforschung, Goettingen, Germany, in collaboration with CISAS, University of Padova, Italy, the Labo ratoire 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 (programme Classy, ANR-17-CE31-0004).

References

  1. Auger, A. T., Groussin, O., Jorda, L., et al. 2015, A&A, 583, A35 [NASA ADS] [CrossRef] [EDP Sciences] [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 Conf. Ser., 461, eds. P. Ballester, D. Egret, & N. P. F. Lorente, 263 [Google Scholar]
  4. Bibring J. P., Langevin, Y., Carter, J., et al. 2015, Science, 349, aab0671 [CrossRef] [Google Scholar]
  5. Biele J., Ulamec, S., Maibaum, M., et al. 2015, Science, 349, aaa9816 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cambianica, P., Fulle, M., Cremonese, G., et al. 2020, A&A, 636, A91 [EDP Sciences] [Google Scholar]
  7. Capaccioni F., Coradini, A., Filacchione, G., et al. 2015, Science, 347, a0628 [NASA ADS] [CrossRef] [Google Scholar]
  8. Ciarniello, M., Capaccioni, F., Filacchione, G., et al. 2015, A&A, 583, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Ciarniello, M., De Sanctis, M. C., Ammannito, E., et al. 2017, A&A, 598, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Ciarniello, M., De Sanctis, M. C., Raponi, A., et al. 2020, A&A, 634, A39 [CrossRef] [EDP Sciences] [Google Scholar]
  11. Cuzzi, J. N., French, R. G., & Dones, L. 2002, Icarus, 158, 199 [NASA ADS] [CrossRef] [Google Scholar]
  12. Davidsson, B. J. R., Sierks, H., Güttler, C., et al. 2016, A&A, 592, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. de Sanctis, M.C., Capaccioni, F., Ciarniello, M., et al. 2015, Nature, 525, 500 [NASA ADS] [CrossRef] [Google Scholar]
  14. Deshapriya, J.D. P., Barucci, M. A., Fornasier, S., et al. 2016, MNRAS, 462, S274 [NASA ADS] [CrossRef] [Google Scholar]
  15. El-Maarry, M.R., Thomas, N., Giacomini, L., et al. 2015, A&A, 583, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. El-Maarry, M.R., Thomas, N., Gracia-Berna, A., et al. 2016, A&A, 593, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. El-Maarry, R., Groussin, O., Thomas, N., et al. 2017, Science, 355, 1392 [NASA ADS] [CrossRef] [Google Scholar]
  18. Feller, C., Fornasier, S., Hasselmann, P. H., et al. 2016, MNRAS, 462, S287 [NASA ADS] [CrossRef] [Google Scholar]
  19. Filacchione, G., Capaccioni, F., Ciarniello, M., et al. 2012, Icarus, 220, 1064 [NASA ADS] [CrossRef] [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 [NASA ADS] [CrossRef] [Google Scholar]
  22. Fornasier, S., Hasselmann, P. H., Barucci, M. A., et al. 2015, A&A, 583, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Fornasier, S., Mottola, S., Keller, H. U., et al. 2016, Science, 354, 1566 [NASA ADS] [CrossRef] [Google Scholar]
  24. Fornasier, S., Feller, C., Lee, J.-C., et al. 2017, MNRAS, 469, S93 [CrossRef] [Google Scholar]
  25. Fornasier, S., Feller, C., Hasselmann, P. H., et al. 2019a, A&A, 630, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Fornasier, S., Hoang, V. H., Hasselmann, P. H., et al. 2019b, A&A, 630, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Fornasier, S., Hasselmann, P. H., Deshapriya, J. D. P., et al. 2020, A&A, 644, A142 [CrossRef] [EDP Sciences] [Google Scholar]
  28. Gehrels, T., Coffeen, D., Owings, D. 1964, AJ, 69, 826 [NASA ADS] [CrossRef] [Google Scholar]
  29. Hasselmann, P. H., Barucci, M. A., Fornasier, S., et al. 2017, MNRAS, 469, S550 [CrossRef] [Google Scholar]
  30. Hasselmann P. H., Barucci, M. A., Fornasier, S., et al. 2019, A&A, 630, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Heinisch P., Auster, H. U., Gundlach, B., et al. 2019, A&A, 630, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Herique, A., Kofman, W., Zine, S., et al. 2019, A&A, 630, A6 [CrossRef] [EDP Sciences] [Google Scholar]
  33. Hoang, H. V., Fornasier, S., Quirico, E., et al. 2020, MNRAS, 498, 1221 [NASA ADS] [CrossRef] [Google Scholar]
  34. Jaumann, R., Schmitz, N., Ho T.-M., et al. 2019, Science, 365, 817 [NASA ADS] [CrossRef] [Google Scholar]
  35. Jorda, L., Gaskell, R., Capanna, C., et al. 2016, Icarus, 277, 257 [NASA ADS] [CrossRef] [Google Scholar]
  36. Keller, H. U., Barbieri, C., Lamy, P., et al. 2007, Space Sci. Rev., 128, 433 [NASA ADS] [CrossRef] [Google Scholar]
  37. Keller, H. U., Mottola, S., Davidsson, B., et al. 2015, A&A, 583, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Keller H. U., Mottola, S., Hviid, S. F., et al. 2017, MNRAS, 469, S357 [CrossRef] [Google Scholar]
  39. Knapmeyer, M., Fischer, H. H., Knollenberg, J., et al. 2018, Icarus, 310, 165 [NASA ADS] [CrossRef] [Google Scholar]
  40. Kofman, W., Herique, A., Barbin, Y., et al. 2015, Science, 349, aab0639 [Google Scholar]
  41. Krüger, H., Seidensticker, K. J., Fischer, H. H., et al. 2015, A&A, 583, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Lai, I.-L., Ip, W. H., Su, C.-C, et al. 2016, MNRAS, 462, S533 [CrossRef] [Google Scholar]
  43. Lantz, C., Brunetto, R., Barucci, M. A., et al. 2017, Icarus, 285, 43 [NASA ADS] [CrossRef] [Google Scholar]
  44. Lauretta D. S., Dellagiustina, D. N., Bennett, C. A., et al. 2019, Nature, 568, 55 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  45. Lee, J.-C., Massironi, M., Ip, W-H., et al. 2016, MNRAS, 462, S573 [Google Scholar]
  46. Li, J.-Y., Helfenstein, P., Buratti, B., Takir, D., & Clark, B. E., 2015, in Asteroids IV, eds. P. Michel, F. E. DeMeo, & W. F. Bottke (Tucson: University of Arizona Press), 129 [Google Scholar]
  47. Li, J.-Y., Schröder, S. E., Mottola, S., et al. 2019, Icarus, 322, 144 [NASA ADS] [CrossRef] [Google Scholar]
  48. Longobardo, A., Palomba, E., Capaccioni, F., et al. 2017, MNRAS, 469, S346 [CrossRef] [Google Scholar]
  49. Massironi, M., Simioni, E., Marzari, F., et al. 2015, Nature, 526, 402 [NASA ADS] [CrossRef] [Google Scholar]
  50. Marshall, D. W., Hartogh, P., Rezac, L., et al. 2017, A&A, 603, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Nelson, R. M., Buratti, B. J., & Wallis B. D. 1987, J. Geoph. Res., 92, 14905 [NASA ADS] [CrossRef] [Google Scholar]
  52. Oklay, N., Mottola, S., Vincent, J.-B., et al. 2017, MNRAS, 469, S582 [CrossRef] [Google Scholar]
  53. O’Rourke, L., Heinisch, P., Blum, J., et al. 2020, Nature, 586, 697 [CrossRef] [Google Scholar]
  54. Pätzold, M., Andert, T., Hahn, M., et al. 2016, Nature, 530, 63 [NASA ADS] [CrossRef] [Google Scholar]
  55. Pajola, M., Hofner, S., Vincent, J. B., et al. 2017, Nat. Astron. 1, 0092 [CrossRef] [Google Scholar]
  56. Perna, D., Barucci, M. A., Fulchignoni, M., et al. 2018, Plan. Space Sci., 157, 82 [NASA ADS] [CrossRef] [Google Scholar]
  57. Penasa, L., Massironi, M., Naletto, G., et al. 2017, MNRAS, 469, S741 [CrossRef] [Google Scholar]
  58. Pilorget, C., Fernando, J., Ehlmann, B. L., Schmidt, F., & Hiroi, T. 2016, Icarus, 267, 296 [NASA ADS] [CrossRef] [Google Scholar]
  59. Poch, O., Istiqomah, I., Quirico, E., et al. 2020, Science, 367, aaw7462 [NASA ADS] [CrossRef] [Google Scholar]
  60. Poulet, F., Lucchetti, A., Bibring, J. -P., et al. 2016, MNRAS, 462, S23 [CrossRef] [Google Scholar]
  61. Quirico, E., Moroz, L. V., Schmitt, B., et al. 2016, Icarus, 272, 32 [NASA ADS] [CrossRef] [Google Scholar]
  62. Raponi, A., Ciarniello, M., Capaccioni, F., et al. 2020, Nat. Astron., 4, 500 [CrossRef] [Google Scholar]
  63. Schröder, S. E., Grynko, Y., Pommerol, A., et al. 2014, Icarus, 239, 201 [NASA ADS] [CrossRef] [Google Scholar]
  64. Schroeder I. R. H. G., Altwegg, K., Balsiger, H., et al. 2019, MNRAS 489, 4734 [NASA ADS] [CrossRef] [Google Scholar]
  65. Sierks, H., Barbieri, C., Lamy, P. L., et al. 2015, Science, 347, aaa1044 [Google Scholar]
  66. Spohn T., Knollenberg, J., Ball, A. J., et al. 2015, Science, 349, aab0464 [NASA ADS] [CrossRef] [Google Scholar]
  67. Sunshine, J. M., A’Hearn, M. F., Groussin, O., et al. 2006, Science, 311, 1453 [NASA ADS] [CrossRef] [Google Scholar]
  68. Tatsumi, E., Domingue, D., Schröder, S., et al. 2020, A&A, 639, A83 [CrossRef] [EDP Sciences] [Google Scholar]
  69. Thomas, N., Sierks, H., Barbieri, C., et al. 2015, Science, 337, a0440 [Google Scholar]
  70. Thomas, N., El Maarry, M. R., Theologou, P., et al. 2018, Plan. Space Sci., 164, 19 [NASA ADS] [CrossRef] [Google Scholar]
  71. Tubiana, C., Güttler, C., Kovacs, G., et al. 2015, A&A, 583, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Vincent, J.-B., Bodewits, D., Besse, S., et al. 2015, Nature, 523, 63 [NASA ADS] [CrossRef] [Google Scholar]
  73. Vincent, J.-B., A’Hearn, M. F., Lin, Z.-Y., et al. 2016, MNRAS, 462, S184 [CrossRef] [Google Scholar]
  74. Warell, J., & Bergfors, C. 2008, Plan. Space Sci., 56, 1939 [NASA ADS] [CrossRef] [Google Scholar]
  75. Warren, S. G., & Brandt, R. E., 2008, J. Geophys. Res. 113, D14220 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Observing conditions of the Wosret region for the NAC images.

Table 2

Phase reddening coefficients evaluated for Wosret during different observations and those reported in the literature for comet 67P and for other dark small bodies.

Table 3

Water-ice content estimates for the bright features observed in the 13 June and 2 July 2016 images.

Table 4

Summary of the observed behavior and physical properties for the small and big lobe of comet 67P.

All Figures

thumbnail Fig. 1

Wosret region on comet 67P. Left: 3D view of the southern hemisphere of 67P (from an image acquired on 2 May 2015 at 10h42) and overlay regional boundaries to facilitate locating the Wosret region on the nucleus. The final landing site of Philae, Abydos, is indicated by the yellow arrow. Center: 3D view of the small lobe (always from the same May 2015 data shown on the left side) with the geomorphological map of Wosret from Lee et al. (2106) superimposed. Right: Sub-units on top identified by Thomas et al. (2018) in Wosret (Wr a, Wr b, and Wr c) and other regions (Ne, So, An, Gb, and Be refer to Neith, Sobek, Anhur, Geb, and Bes regions, respectively), with the associated comet view from an image acquired on 2 January 2016 at 17h23 (bottom-right).

In the text
thumbnail Fig. 2

Comparison between images acquired on May 2015 (on the left, spatial resolution (res) of 2.3 m px−1) and on July 2016 (on the right, res = 0.3 m px−1) covering part the Wosret region (area a following Thomas et al. 2018 sub-regions definition). A new cluster of outcrops (named A and indicated by the white arrow) is clearly visible on the July 2016 images. Some boulders are indicated for reference (and numbered 1, 2, and 3). The ellipse named B includes a region where the dust coating may also have been thinned during the perihelion passage, as observed close to the cluster of outcrops named A and to the goosebump features.

In the text
thumbnail Fig. 3

Comparison between images acquired on March 2015 (left side, res = 1.9 m px−1) and June 2016 (right side, res = 0.5 m px−1) showing theformation of a new cavity inside the ellipse named C. Some boulders are indicated for reference by white bars and numbers, following the same numbering scheme of Fig. 2.

In the text
thumbnail Fig. 4

Goosebumps equivalent diameter from the 2 July 2016 image (see Fig. 2).

In the text
thumbnail Fig. 5

Spectral slope evolution with heliocentric distance comparing data acquired on 2 May 2015 (rh = 1.73 AU inbound) and 2 July 2016 (rh = 3.32 AU outbound) images. The July 2016 data were phase reddening corrected from 93° to 62° using the phase reddening coefficient of 0.055 × 10−4 nm−1 deg−1 determined in Sect. 4.1 for Wosret in the post-perihelion 2016 observations. Comparing to the original data, we applied a zoom factor of 2 to the May 2015 data, while the July 2016 image was binned by a factor of 2.

In the text
thumbnail Fig. 6

Spectral slope of the Wosret region, evaluated in the 535–882 nm range, from different images in the April 2015–July 2016 time frame. The Wosret region is indicated by white rectangles in the lower resolution images. The high resolution image acquired on July 2016 was binned by a factor of 2.

In the text
thumbnail Fig. 7

Phase reddening in spectral slopes evaluated in the 535–882 nm range. The black circles represent data acquired in 2016, redcircles those acquired in May–October 2015, that is three months before and after the perihelion passage, while the remaining 2015 data are represented with blue circles. The black and red lines are the linear fit of the 2016 and near perihelion data, respectively.

In the text
thumbnail Fig. 8

RGB images acquired on 17 December 2015 taken one hour apart. The right image acquired at UT23h56 shows a relatively blue region (indicated by the white rectangle) at dawn. The same area was observed in a sequence acquired one hour later (on 18 December 2015, UT 00:56) and it shows a color indistinguishable from that of the average Wosret region. The bluing at dawn is periodic and observed in other sequences in December 2015 but also in June 2015, as shown in the top insert of the image on the right side.

In the text
thumbnail Fig. 9

RGB image taken on 12 June 2016 at 12h48 showing a blue patch at latitude − 21° and longitude 318° in an area observed to be active by Vincent et al. (2016, jet numbered 11 in their paper).

In the text
thumbnail Fig. 10

Wosret spectrophotometry. Top: examples of blue patches and bright spots from 13 June 2016 (left) and 2 July 2016 images (right). Inserts show a zoom of some ROIs. Middle and bottom panels: corresponding relative reflectance and the I over F value at the given phase angle (70° and 93° for the June and July images, respectively). The dark terrain used in the mixing model is represented by a circle in white for clarity in the RGB images on the top and in black in the plots.

In the text
thumbnail Fig. 11

Normal albedo of the bright spots (indicated by the magenta squares and the blue asterisks, see Fig. 10) observed on 2 July 2016 UT15h22. The symbols correspond to the bright spots shown on the right panel of Fig. 10. The black circles represent the mean spectrum of the comet from a region close to the bright patches. Continuous and dashed lines show the best fit spectral models to the bright patches (in gray the associated uncertainty), produced by the linear mixture of the comet dark terrain enriched with 29 ± 3% of water ice with 1000 μm grain size (continuous line) for the ROI indicated by the blue asterisk, and with 65.5 ± 6.0% of water ice with 100 μm grain size (dashed line) for that represented by the magenta square. Models including water ice of medium and largegrain size are represented for comparison: dashed lines indicate models with water ice of 100 μm grain size, and continuous lines those with 1000 μm grain size.

In the text
thumbnail Fig. 12

Map of the comet with superposed the sources of jets identified in the Wosret region during the comet southern hemisphere summer. The location of some nearby jets are averaged for clarity, and some notable jets are represented with larger symbols. Blue circles represent jets reported in Vincent et al. (2016), green triangles those reported in Fornasier et al. (2019a), and black squares represent cavities showing repeated activity in several data sets (Fornasier et al. 2019a).

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.