Diurnal variation of dust and gas production in comet 67P/Churyumov-Gerasimenko at the inbound equinox as seen by OSIRIS and VIRTIS-M on board Rosetta

Context. On 27 April 2015, when comet 67P/Churyumov-Gerasimenko was at 1.76 au from the Sun and moving toward perihelion, the OSIRIS and VIRTIS-M instruments on board the Rosetta spacecraft simultaneously observed the evolving dust and gas coma during a complete rotation of the comet. Aims. We aim to characterize the spatial distribution of dust, H 2 O, and CO 2 gas in the inner coma. To do this, we performed a quantitative analysis of the release of dust and gas and compared the observed H 2 O production rate with the rate we calculated using a thermophysical model. Methods. For this study we selected OSIRIS WAC images at 612nm (dust) and VIRTIS-M image cubes at 612nm, 2700nm (H 2 O emission band), and 4200nm (CO 2 emission band). We measured the average signal in a circular annulus to study the spatial variation around the comet, and in a sector of the annulus to study temporal variation in the sunward direction with comet rotation, both at a ﬁxed distance of 3.1km from the comet center. Results. The spatial correlation between dust and water, both coming from the sunlit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO 2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0 ◦ subsolar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO 2 production. An excess of measured water production with respect to the value calculated using a simple thermophysical model is observed when the head lobe and regions of the southern hemisphere with strong seasonal variations are illuminated (subsolar longitude 270 ◦ –50 ◦ ). A drastic decrease in dust production when the water production (both measured and from the model) displays a maximum occurs when typical northern consolidated regions are illuminated and the southern hemisphere regions with strong seasonal variations are instead in shadow (subsolar longitude 50 ◦ –90 ◦ ). Possible explanations of these observations presented


Introduction
After a ten-year journey and 30 months of deep-space hibernation, the ESA Rosetta spacecraft woke up on 20 January 2014. Rosetta had the unique opportunity to stay in the vicinity of comet 67P/Churyumov-Gerasimenko (hereafter 67P) for 2.5 yr and to observe how the comet evolved while moving along its orbit.
One of the main goals of the Rosetta mission was to understand cometary activity, that is, the physical processes that generate the dust and gas coma from the nucleus. While the broad picture of the Whipple model, that is, that ices in the nucleus sublimate when heated by the Sun and the resulting gas outflow lifts dust (Whipple 1950), has been confirmed by observations, details of the processes involved remain the subject of debate. Observations of dust jets have been traced to certain areas of the surface , and various models have been developed to trace gas and dust flow in the inner coma of 67P, with varying degrees of complexity (e.g., Fougere et al. 2016;Kramer & Noack 2016;Kramer et al. 2017;Zakharov et al. 2018), but the models do not yet uniquely identify the surface features that are responsible for activity. Many of the models show that the bulk activity can be explained by more or less homogeneous activity from all illuminated surface facets (Keller et al. 2015) and that jets in the inner coma are controlled more by the complex shape of the nucleus than by anything special about their apparent source on the surface (Shi et al. 2018a). There are clear variations in activity with seasonal illumination of the comet, however, which appear to be related to the very different morphology of the northern and southern surfaces, and models that explain the early activity seen by Rosetta when 67P was far from the Sun do not reproduce the perihelion behavior (Shi et al. 2018b). Recent models attempt to reproduce the complexities of the changing activity (Attree et al. 2019;Marschall et al. 2019); they also include investigations of the relative contribution of different sublimating ices (water or CO 2 ) to the driving activity, that is, how the ices are related to each other (Gasc et al. 2017) and to dust release.
The aim of this work is to take advantage of the capabilities of two instruments on Rosetta to analyze the dust and gas coma behavior in the pre-perihelion phase, when the comet was at heliocentric distances of 1.76 au and close to the equinox between the changing seasons. This allows us to investigate the differences in observed dust and gas distributions and their longitudinal variations simultaneously over a full rotation of the nucleus. We investigate both the spatial variation in dust and gas around the comet and how it varies with time as different areas are illuminated throughout the cometary day.

Instruments
The Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) and the Visible InfraRed Thermal Imaging Spectrometer (VIRTIS) were 2 of the 12 scientific instruments on board the Rosetta orbiter (Glassmeier et al. 2007). OSIRIS (Keller et al. 2007) was the scientific camera system. It comprised a Narrow Angle Camera (NAC) and a Wide Angle Camera (WAC) with fields of view (FOV) of 2.20 • × 2.22 • and 11.35 • × 12.11 • , respectively. Both cameras used a 2048 × 2048 pixel backside-illuminated CCD detector with a UV optimized antireflection coating. The CCDs were equipped with lateral antiblooming that allowed overexposure of the nucleus without creating saturation artifacts, enabling the study of details in the faint coma structures next to the illuminated limb. The NAC was equipped with 11 filters covering the wavelength range 250-1000 nm, while the WAC had 14 filters covering the range 240-720 nm (Tubiana et al. 2015).
The VIRTIS spectrometer (Coradini et al. 2007) was composed of two spectral channels: VIRTIS-M and VIRTIS-H. VIRTIS-M was the visible (230-1000 nm, 432 bands) and infrared (1000-5000 nm, 432 bands) imaging spectrometer with a field of view of 3.6 • (along the slit axis) and an instantaneous field of view (IFOV) of 250 µrad. The instrument acquired hyperspectral cubes by scanning in time the target scene line by line. The duration of the acquisition (∆t in Table 1) is given by the number of lines (including periodic dark current frames) times the internal repetition time, where the repetition time is the time between two consecutive steps that are required to move the internal scan mirror by one IFOV. The integration time (t exp in Table 1) set for the VIS and IR channels was lower than the internal repetition time. The maximum 3.6 • × 3.6 • FOV was imaged by repeating acquisition on successive 256 scan mirror steps (lines). From a distance of 100 km, this corresponds to a 6.4 km × 6.4 km swath with a resolution of 25 m pix −1 . As an example, we show in Fig. 1 a VIRTIS-M hyperspectral cube with the line and time axes. VIRTIS-H was the infrared high-spectral resolution point spectrograph operating in the 1900-5000 nm spectral range with a λ/∆λ = 1300 ÷ 3000. The instrument observed in a single IFOV of 580 µrad × 1740 µrad, which corresponds to a resolution of 58 m × 174 m from a 100 km distance. Because the VIRTIS-H and VIRTIS-M boresights were coaligned, the point spectrograph was able to acquire the same area several times during the time required for a scan for the imaging channel.

Datasets
Throughout the entire mission, we regularly carried out socalled dust monitoring sequences, which were designed to observe the coma of unresolved dust particles. Typically, these sequences spanned at least 12 h (i.e., a full comet rotation) with hourly cadence. In the frame of this work, we have analyzed one dust monitoring sequence that was acquired on 27 April 2015 1 . At the time of the observation, 67P was at a heliocentric distance of 1.76 au, moving toward perihelion, and close to the equinox between the long but cool northern summer and the short intensely illuminated southern summer around perihelion. The Rosetta spacecraft was at a distance between 125 and 142 km from the comet. The Rosetta +Z-axis was pointing to the comet nucleus (IlluminatedPoint pointing), and because of the spacecraft-comet distance, the nucleus was entirely contained in all OSIRIS WAC images and in seven out of nine VIRTIS-M images.
The spacecraft was approximately in a terminator orbit (phase angle ∼73 • -75 • ) so that one side of the comet was illuminated by the Sun and the other side was in darkness. The observational details are summarized in Tables 1 and 2. The dataset of 27 April 2015 is one of the best acquired by VIRTIS-M for this purpose because it covers both the VIS and IR channels and spans more than 12 h. It was one of the final monitoring observations before the failure of the VIRTIS-M cryocooler disabled the IR channel. It was optimized for coma observations and achieved a high signal-to-noise ratio (S/N). The same dataset was previously analyzed by Rinaldi et al. (2016) and Fink et al. (2016). Rinaldi et al. (2016) focused their analysis on the comparison between the spatial distributions of dust and H 2 O and CO 2 gas, radial profiles, and azimuthal distributions to search for any correlation between them. Fink et al. (2016) focused their investigation on the emission intensity of CO 2 and H 2 O and provided an explanation for the large observed variations reported in the literature for the CO 2 -to-H 2 O ratio. In this work we study temporal variation in more detail, and can better study areas of the dust continuum that were heavily affected by stray light in the VIRTIS-M data by combining it with OSIRIS imaging.

OSIRIS dataset
The dust monitoring sequence STP053_DUST_MON_006 contained 45 WAC full-frame images. For this study we selected the 15 images acquired with the VIS610 filter (λ cent = 612.6 nm, ∆λ = 9.8 nm) and with exposure time optimized for dust coma studies. One image was acquired with the WAC door closed and it was excluded from the analysis. Thus, we used 14 OSIRIS images in total. We used OSIRIS level 3 (CODMAC Level 4) images, which are radiometric calibrated and geometric distortion corrected (for details, see a description of the OSIRIS calibration pipeline in Tubiana et al. 2015). The images, scaled to the same intensity levels, are shown in Fig. 3 (first column).

VIRTIS-M dataset
A set of nine image cubes was obtained. The image cubes in the two channels were taken at the same time, with the VIS exposure   Fig. 1 and Sect. 2.3); Col. 10: subsolar longitude of each observation at the time of the middle of the annulus (see Fig. 1 and Sect. 2.3). sequences starting about 2-4 s after the IR sequences. The first seven cubes have a FOV of about 9.0 km × 7.7 km (at the nucleus center distance), and the position of the cometary nucleus is roughly in the center of the image. The last two cubes have a FOV of about 9.0 km × 4 km, and the nucleus is only partially contained in the frame. Figure 3 shows the spatial distribution of dust (612 nm), water vapor (2700 nm), and carbon dioxide (4200 nm). Each map is a composite image where the comet nucleus, taken at 4200 nm, is superimposed on the maps of the dust continuum at 612 nm (second column), the water vapor (third column), and the CO 2 band intensities (fourth column). When the bright nucleus partially illuminates the instrument slit, a sizable fraction of the incoming photons is spread into the adjacent coma pixels. The measured dust continuum is contaminated by stray light caused by the nucleus and cannot be used. The data cubes used for the analysis were calibrated using the The dust continuum, in the range 200-3000 nm, is due to sunlight scattered by the dust particles in the coma. The IR spectrum shows the gas fluorescence emission of water vapor and CO 2 at 2700 and 4200 nm, respectively. VIRTIS reduction pipeline as described by Ammannito et al. (2006), Filacchione et al. (2006), and Rinaldi et al. (2016). A typical radiance spectrum of the 67P coma in the VIS and IR is shown in Fig. 2. It demonstrates the unique capability of the VIRTIS-M instrument to simultaneously measure the dust continuum in the range 200-3000 nm and the fluorescence emission of water vapor and CO 2 . The calculation of H 2 O and CO 2 band intensity was described by Migliorini et al. (2016) and Fink et al. (2016).
To measure the dust continuum intensity, we chose a 9.8 nm wide band centered at 612.6 nm, which corresponds to the OSIRIS WAC VIS610 filter. This allowed a comparison of the results obtained by the two instruments.

Aperture photometry
In each image, we measured the average signal in an annulus, or in a sector of the annulus, at fixed distance (in kilometers) from the center of the comet (Fig. 1). We call the selected area where the flux is measured "mask". We chose a mask width of 0.2 km in radial direction for the entire dataset. As distance we chose the maximum distance from the comet center for which the mask was fully included in the VIRTIS-M frame, which is smaller than the OSIRIS WAC FOV. For this dataset the selected distance is 3.1 km. The irregular shape of the comet means that a single-circle mask with a fixed distance from the center of the comet is not equidistant to the limb of the comet.
Owing to the VIRTIS-M operation mode (see Sect. 2.1), each image cube line was acquired at a different time and therefore has a different longitude of the subsolar point. For each image cube, we determined the mid-point inside each mask, as shown in Fig. 1, and calculated the subsolar longitude of the mid-point. The subsolar longitude is mask dependent because the number of considered lines is different in the two masks, therefore the subsolar longitude is different for a sector or for an annulus in the same image cube (see Table 1).
For the OSIRIS images, the statistical error associated with each measurement was determined using the sigma (or error) map (Tubiana et al. 2015). It contains the error associated with the intensity of each pixel, calculated using Poisson statistics and the readout noise error. This statistical error is very small, about ±0.3%. In addition to this statistical error, the images have a systematic error due to the radiometric calibration of ±1% (Tubiana et al. 2015). VIRTIS-M measurements have a statistical error given by the standard deviation of the average flux inside the mask and a systematic error due to the radiometric calibration. The total uncertainty on the measurement, calculated using error propagation, is about ±10% Coradini et al. 2007).

Azimuthal profiles
To determine the azimuthal profiles of dust and gas, we used the circular mask at 3.1 ± 0.2 km from the center of the comet. The selected angular step was 10 • . The profiles were measured clockwise. Zero degree is in the subsolar direction, as sketched in Fig. 1.

Longitudinal variation
The strongest dust signal is observed in the subsolar direction ( Fig. 3; see also Sect. 5). To study the overall dust coma, we selected a mask with 90 • angular size (±45 • ) in the subsolar direction, as shown in Fig. 1 (right panel). We chose this angular size to minimize the contribution of fine-scale structures (i.e., jets) in the images. Figure 4 shows the longitudinal variation of the dust flux. Each point of the figure represents the average dust flux inside the ±45 • mask measured in OSIRIS (blue triangles) and VIRTIS-M (orange squares) images. The measured average flux in the ±45 • mask is summarized in Table A.1.
The OSIRIS and VIRTIS-M measurements are in very good agreement at the wavelength used for the analysis (Fig. 4). In addition, the good agreement gives us the possibility to directly compare OSIRIS measurements with those of VIRTIS-M at different wavelengths, without having to consider possible instrumental effects. As described in Sect. 2.2.2, the analysis of the dust continuum in the VIRTIS-M data is limited by in-field stray light when the instrument slit is partially illuminated by the bright nucleus, as is the case in our observations. For this reason, we used only the OSIRIS images to study the dust in all subsequent sections. This avoids interpolating across stray-light regions, as was necessary in previous work . We used VIRTIS-M to measure the gas.

Dust brightness and Afρ
We calculated A f ρ, which is commonly used to quantify dust brightness in comets (A'Hearn et al. 1984) and is often used as a proxy of dust production; it is proportional to the dust-loss rate if the dust size distribution and velocity are constant. As we do not know the size distribution or velocity, we only used the observed flux of scattered light to quantify the dust in each image, and included the conversion into A f ρ only to allow convenient comparison with other observations, but did not attempt to derive any absolute dust production rate in kg s −1 . The average flux in a full annulus as a function of subsolar longitude is shown in Fig. 5 and summarized in Table A.1. About 50% of the total dust flux in the full annulus comes from the 90 • sector in subsolar direction, as shown in Table A.1.
To translate the observed scattered light intensity along the line of sight (LoS) into a local A f ρ, we used the method developed in Fink & Rubin (2012) and Fink & Rinaldi (2015). When the considered annulus is at a sufficiently large distance from the nucleus, in the collision-free flowing zone, and when no additional production or destruction of dust occurs, the calculated A f ρ can provide a global measure of A f ρ in the immediate vicinity of the nucleus, and it will miss only a small fraction of the total emitted dust. Rinaldi et al. (2016) found that at a distance closer than 4 km from the surface, the dust intensity decreases much faster than 1/ρ, which implies that the dust acceleration region is sampled by our measurements at 3.1 km from the cometary center. Gerig et al. (2018) determined that the average starting point of the 1/ρ behavior is (11.9 ± 2.8) km. Because the steady-state condition is not fulfilled here, the determined A f ρ cannot be directly compared with ground-based global measurements. Nevertheless, the result we obtain is similar to the ground-based value of A f ρ ∼ 0.9−1.0 m at this time (Snodgrass et al. 2016). We obtained an A f ρ between 1.0 m and 1.4 m, as shown in Fig. 5.

Gas production rate derived from VIRTIS-M data
To determine the gas production rates, we calculated the average emitted band intensity inside the 3.1 km annulus and translated it into a gas column density (n (ρ)) (Migliorini et al. 2016;Fink et al. 2016). For these calculations we discarded the last image cube (# 9) because of some radiometric problems at the wavelengths close to the H 2 O gas emission. In image cube # 8 (see Fig. 3) only the dayside part of the coma is observed, so that this image cube was also discarded. Only in the first seven image cubes (# 1-7) the full nucleus is inside the FOV, which is required for retrieving the complete azimuthal behavior of the gas in the coma.
For the combined H 2 O bands at 2660 nm and 2730 nm we used the fluorescence efficiency at 1 au g 0 = 2.745 × 10 23 W molec −1 (Bockelée- Morvan et al. 2015), and for the CO 2 band we used g 0 = 1.25 × 10 22 W molec −1 (Debout et al. 2016). The uncertainty on the H 2 O and CO 2 column density calculations is ∼10%. The measured average water and CO 2 column densities as functions of subsolar longitude are displayed in Fig. 5 (center and bottom panels) and summarized in Table A.1. At 3.1 km from the comet center, 78% of the water column density is contained within an angle of ±90 • and 52% within an angle of ±45 • in the subsolar direction. There is very little scatter in the percentages in the first seven observations despite the different configurations of nucleus, Sun, and spacecraft. The CO 2 distribution does not follow the direct solar illumination and has essentially no correlation with the water column density distribution (Fig. 5). The CO 2 molecules are emitted mostly from the southern hemisphere of the comet, as shown in Fig. 3 (fourth column). This is the reason why 50% of the CO 2 column density is contained within an angle of ±90 • in the subsolar direction (day side) and the same percentage is found on the night side, with little variability in these percentages (Table A.1). The large variability of the CO 2 column density, from 19 to 38%, in the ±45 • sector in subsolar direction is due to the orientation of the cometary spin axis during these observations.
To convert the gas column densities into production rates, we used the method described in Fink et al. (2016). The gas production rates are listed in Table A Fink et al. (2016) and the gas spatial distribution maps obtained by Migliorini et al. (2016). Figure 5 shows that there is little correlation between the measured productions of dust, water, and CO 2 . In addition, none of the patterns corresponds to the variation of the cross-section area of the illuminated nucleus ( Fig. 7 bottom panel). The lack of correlation between dust and CO 2 is not surprising, indicating that CO 2 outgassing does not play a dominant role in driving the global dust emission. On the other hand, the discrepancy between the patterns of dust and water measurements is harder to explain. Water outgassing is the dominant driver of dust activity from the northern hemisphere (De Sanctis et al. 2015;Shi et al. 2018a) before the observations we analyzed here (days before northern autumn equinox). The deviation of water production from the variation of the illuminated cross section, in particular, suggests that the nucleus surface is not homogeneously active, as presented by Marschall et al. (2019). In the following section, we perform a simple thermophysical analysis in order to shed some light on this discrepancy.

H 2 O production rate computed with a thermophysical model
We wish to compare the observed pattern of dust and gas release to a simplified model that describes what would be expected in the case of homogenous activity that is only controlled by illumination of each surface element. In order to do this, we employed a thermophysical model (Hu et al. 2017a) to estimate the total water production rate of the nucleus over a full rotation of 67P that encompasses all observations described above.

Model description
The shape of the nucleus is approximated by a model consisting of 1500 facets (Preusker et al. 2015). We assumed that the nucleus is covered by a desiccated dust layer, or dust mantle, of constant thickness that is composed of uniform spherical dust aggregates. Water ice is present underneath the dust mantle (Fig. 6). The sublimation flux is strongly influenced by the temperature of the ice front. The temperatures of the nucleus subsurface as a function of depth are estimated by solving the 1D heat equation, balancing the input energy from solar illumination of the surface with heat that is reradiated, conducted into the surface, or used in sublimating ice. However, the interior of the nucleus (beneath the mantle) is not composed of pure water ice, and sublimation can therefore not take place everywhere. For this reason, the factor f i ∈ (0, 1) is introduced. It formally measures the areal fraction of water ice and is approximately inverse to the dust-to-ice ratio of the subsurface (Crifo 1997). We assumed that heat flux vanishes beyond several (diurnal) skin depths. At any given epoch, the position vector of the Sun with respect to (the body-fixed frame of) the nucleus is obtained through the SPICE kernels 2 for 67P. It is subsequently transformed into local horizontal coordinates through a series of rotations of the coordinate system to yield the solar incidence angle. We made use of a "Landscape" database for 67P that delineates the skyline at each location on the nucleus (each facet of the shape model) in order to efficiently determine the local illumination (Hu et al. 2017a).
The 1D heat equation is solved with the Crank-Nicolson method. The solutions are diurnally equilibrated temperatures T and water production rates Z that repeat or coincide, if the seasonal cycle is neglected, with the rates exactly one comet rotation apart, for example, T (t ± t p ) = T (t) and Z(t ± t p ) = Z(t), where t P is the rotation period of 67P.

Choice of model parameters
A summary of the key parameters of the thermophysical model is given in Table 3. For simplicity, it is assumed that the nucleus subsurface is homogeneous. We also neglect variability of the parameters, for instance, changes in mantle thickness and loss of ice underneath. With these assumptions, the model parameters are treated as constants. Following the argument by Blum et al. (2017) that the dust aggregates, clearly nonuniform in size, are about a few millimeter in diameter, we adopted a diameter of the dust aggregates of d P = 1 mm. The OSIRIS observation of the dust activity that continued for about one hour after sunset indicates that water ice was present at some depth of less than 1 cm (Shi et al. 2016). The long-term (seasonal) evolution of the total water production of 67P throughout perihelion can also be modeled with a mantle thickness of 5 ≤ X i ≤ 10 mm. Overall, water ice is rarely exposed . When detected, it is usually present in small quantities, for example, a few percent (De Sanctis et al. 2015;Filacchione et al. 2016a;Barucci et al. 2016). There is observational evidence that the average water-ice abundance in the top about 1 m above the northern hemisphere cannot exceed 10% (Hu et al. 2017b). Therefore, the area fraction of water ice ( f i ) is probably similar. The scarcity of ice suggests that the thermophysical properties of the nucleus are dominated by those of the refractory component. We are aware that the thermal parameters, that is, the conductivity (κ) and specific heat capacity (c), are dependent on temperature. In particular, thermal radiation through the pores enhances the efficiency of heat transfer. In this case, the radiative component of the conductivity varies with T 3 . However, the radiation most significantly affects temperatures below the diurnal skin. In the case of millimeter-sized particles and when ice sublimation is assumed to occur mostly from above the diurnal skin depth, the enhancement of water production by radiation is not notable (Hu et al. 2019). In addition, the dependence of heat capacity on temperature is linear (Orosei et al. 1995), while the exact behavior of material on 67P is largely unknown. We therefore neglected the temperature dependence of the parameters and adopted κ d = κ i = 0.002 W m −1 K −1 (conductivity of the dry dust mantle and the underlying icy dust, respectively), c = 1000 J kg −1 K −1 and = 500 kg m −3 , which corresponds to a thermal inertia of 30 W m −2 K −1 s 1/2 , as measured by MIRO (Schloerb et al. 2015).

Model results
The water production rates at each observation time from the output of the thermal model and the VIRTIS-M observations described in Sect. 4.1 are shown in Fig. 7. Regulated by the mantle thickness and ice abundance, the modeled water production rate agrees well with the measurements in terms of overall magnitude. The variation largely follows that of the cross-section area of the illuminated nucleus, as expected for a model of homogeneous nucleus activity. A phase shift of about 20 • in subsolar longitude is evidently attributable to the presence of the dust mantle, which causes a thermal lag of about half an hour at the depth of the ice front (i.e., 5 mm). However, the modeled production rate is clearly underestimated compared with the measurements in the subsolar longitude range (270 • -50 • ). This indicates that the real gas production varies between areas and must depend on local (sub-)surface properties, not just illumination and topography of the nucleus; this is discussed further below.   Figure 8 shows the angular distribution of water (red squares) and CO 2 (green circles) band intensity, measured in the first seven VIRTIS-M image cubes, and the angular distribution of the dust brightness (blue triangles), measured in the OSIRIS images that are closest in the subsolar longitude to the VIRTIS-M data. In all observations, the absolute maximum for dust and water vapor is located in the subsolar direction. A secondary water vapor peak occurs at roughly 130 • and is related to the variable illumination of the neck area. CO 2 peaks between 180 • and 270 • , consistent with the southern hemisphere illuminated (Fig. 3).

Azimuthal gas and dust distribution
The spatial (or angular) correlation between dust and water, both coming from the subsolar side of the comet, which has previously been observed by Rinaldi et al. (2016), shows that water is the main driver of dust activity in this time period. This is also generally consistent with observations from the ground, showing that long-term variations in total water production rate correlate with the total dust brightness (Hansen et al. 2016).
The presence of CO 2 ice in these regions of the southern hemisphere and in the same period is not surprising because VIRTIS detected a CO 2 -ice rich area in the Anhur region at the A23, page 8 of 14 C. Tubiana et al.: Properties of dust and gas in the coma of 67P  (Filacchione et al. 2016b). In the southern hemisphere, the CO 2 ice is generally closer to the surface and therefore more easily accessible. When the comet approaches the Sun, the southern summer is very intense. In this strong sunlight, erosion rates are fast enough to expose fresh primordial layers of the interior, which are rich in CO 2 . Quantitative calculations of the comet erosion rates during its orbit are presented by Keller et al. (2015). The increase in CO 2 abundance in the coma during the perihelion time-frame is analyzed in Bockelée-Morvan et al. (2016). The authors found that during this period, the measured abundance ratios of CO 2 increased by a factor of 30 with respect to what was found above the illuminated northern hemisphere.

Discussion
To understand the illumination conditions across the surface of the nucleus of 67P at the time of each OSIRIS observation, we generated insolation maps. One example is shown in Fig. 9 (top panel), along with a map in which the different regions of the nucleus are identified (bottom panel). The distribution of solar irradiance is calculated on a polyhedral shape model representing the nucleus of 67P with 499 902 facets (Preusker et al. 2017). For each epoch, the Sun's position in the Cheops bodyfixed frame of the comet (as defined in Preusker et al. 2015) was derived from the reconstructed ephemeris and rotational status of 67P using the SPICE tool-kit (Acton 1996). The concave shape of 67P prevents the display of certain surface areas in maps with equidistant cylindrical projection. However, the overall illumination pattern shown in the figure is not affected. Figure 10 displays the insolation maps at the time of each OSIRIS observation. In the central inset we show the modeled and observed H 2 O production rates (gray and red curves, respectively), the observed CO 2 production rate (green curve), and the total dust brightness (blue curve), calculated inside the full 360 • annulus (see Table A.2). There is no strong temporal correlation between total dust brightness and water production rates (neither observed nor from the simple model). It should be emphasized (as noted in Sect. 5), however, that water is still the main driver of dust activity in this time period. Only the ratio between water production and dust activity changes.
The green box in Fig. 10 highlights the excess of water production (red curve) compared to the simple homogeneous model (gray curve) at subsolar longitudes between 270 • and 50 • . The green lines emphasize the corresponding insolation maps. This corresponds to epochs when the head lobe and regions of the southern hemisphere with strong seasonal variations (e.g., Bes, Geb, Anhur) are illuminated. While we caution against over interpreting the differences between our simple thermal model and the observed water production, it could be argued that these are most easily explained by a higher activity of these southern regions with respect to the northern regions.
The observed dust brightness (blue curve) shows a pronounced maximum around 0 • subsolar longitude, which is not pronounced in the water production (red curve). Previous studies have already shown that the Anhur and Bes regions, which are illuminated at this time, are highly active and sources of several jets Fornasier et al. 2017), thus in agreement with our findings.
At the same epoch (green region) there is also a maximum of CO 2 . The analysis of the azimuthal profiles in Sect. 5 (Fig. 8) showed that the dust is correlated with water and not CO 2 . The increased CO 2 production in this epoch is therefore not responsible for the peak in the dust activity.
The largest discrepancy between dust and gas production rates can be observed in the red box of Fig. 10. The dust brightness drastically decreases in this subsolar longitude range (50 • -90 • ), while instead the water production (measured and from the model) displays a maximum. These epochs correspond to when northern consolidated regions (e.g., Bastet, Aker, Khepry, Aten, Babi; Thomas et al. 2018;see Fig. 9, bottom panel) are illuminated and the southern hemisphere regions with strong seasons are in shadow.
This temporal noncorrelation (red box) can either be explained by regional variations of surface properties or regional variations in the scattering properties of the lifted dust particles. More specifically, we discuss here the effect of regional variations in (a) thickness of the desiccated layer, (b) surface cohesion, and the presence of large particles that affect (c) gas coupling, thus lifting, and (d) permeability of the dust layer. For the scattering properties of the dust particles, we discuss (e) the composition of lifted dust particles and (f) their sizes.
(a) We do not attribute the discrepancy of dust and gas production in the red box to the thickness of the desiccated layer on the nucleus surface because a thicker layer would also quench gas production. (b) A higher cohesion of surface material could quench dust activity (Bischoff et al. 2019), where the water vapor would escape without lifting dust. This is possible and likely in consolidated regions, which consist of centimeter-sized pebbles. (c) Even if large (decimeter-to meter-sized) particles were easily lifted against their cohesion, they still carry a high mass inertia. Their low size-to-mass ratio could prevent large enough particles from being carried into the coma from the gas drag. Fallback for particles in the considered size range was observed by Agarwal et al. (2016), and many more of these might not even be considerably lifted. (d) If the upper dust layer is dominated by large particles, the gas permeability would also be enhanced. For granular materials with macroscopic voids, Gundlach et al. (2011) have shown that the gas permeability increases with the size of the constituent particles. A layer of decimeter-sized dust A23, page 9 of 14 A&A 630, A23 (2019) Fig. 10. Insolation maps for the time of each OSIRIS observation. Inset: total dust brightness and gas production rates, measured in the full 3.1 km annulus, as functions of subsolar longitude. The green box encloses the insolation maps that correspond to the subsolar longitudes where an excess of measured water production was found, relative to a simple model. The red box frames the insolation maps that correspond to the subsolar longitudes where a minimum of dust and a maximum of water production were measured. particles (cf. Pajola et al. 2017) would thereafter have a ten times higher permeability than a layer of the same thickness of centimeter-sized particles. A higher gas permeability of the dust layer would result in a reduced pressure buildup and thus reduced dust production. The gas could simply escape through the large voids. Decimeter-to meter-sized particles, which are present in fallback regions (e.g., Ma'at), are almost entirely cleaned up in consolidated regions, which exposes the underlying consolidated material. This is consistent with the idea of "self-cleaning" of the northern hemisphere proposed by Fulle et al. (2019). (c) and (d) could quench dust activity in fallback regions, but do not play a role in consolidated areas.
For the coma-related effects (e and f), it is worth repeating that the observed dust brightness (or equivalently A f ρ) is proportional to the dust-loss rate only if the dust velocity, size distribution, and composition do not change. A variation in dust brightness might therefore in principle be interpreted as a change in dust-loss rate or as a change in any of the parameters above. (e) We rule out significant differences in the dust properties in different areas (i.e., at different times) due to composition variation, as none was observed in VIRTIS coma observations ). (f) If the coma were dominated by larger particles, they would tend to reduce the observed dust brightness for the same dust production rate in kg s −1 because the fewer but larger particles would constitute a smaller reflecting area. This would play a role in fallback regions, but not in the consolidated regions that we considered. Measurements of the dust size distribution from the in situ instruments on board Rosetta are difficult to separate by different surface areas as they were measured over extended periods, during which the sub-spacecraft location changed considerably, but there are hints of variation. For example, GIADA detected more compact particles from Hapi (a fallback region) and more fluffy aggregates elsewhere (Della Corte et al. 2015), suggesting variations in the size distribution; but again, this does not demonstrate differences between different consolidated terrains. Differences were observed in size distribution in the dust that was released during outbursts relative to the background coma, but this may be related to the (poorly understood) outburst process rather than regional differences (Bockelée-Morvan et al. 2017). In summary, the best explanation of our observations in the red box in Fig. 10 is a quenched dust activity due to high cohesion of surface material typical of consolidated regions.
The observations in the green box in Fig. 10, an increased water activity in the southern regions with respect to model expectations, are likely to be attributable to regional changes in volatile content or access to these. Our simple thermal model assumes the same thickness of the desiccated layer for the northern and southern hemisphere. The erosion rates in the south are higher, therefore this is likely not the case, and a shallower desiccated layer or a larger area fraction of ice (Fougere et al. 2016) would explain our observations. A23, page 10 of 14 C. Tubiana et al.: Properties of dust and gas in the coma of 67P A more thorough investigation in the near future requires applying more sophisticated thermophysical models to treat not only nonuniform properties of the nucleus subsurface, but also activity of multiple volatile species. It would be necessary to resolve the temperatures in layers deeper than the diurnal skin depth as considered in the current study. A more detailed characterization of the physical processes in the subsurface, such as phase change and mass transfer of volatiles as well as the resulting material loss, is also desired, as the phenomena already proved to influence the energy budget of the system strongly (de Sanctis et al. 1999;Capria et al. 2000;Prialnik et al. 2004;Gortsas et al. 2011).

Summary and conclusion
We have analyzed one OSIRIS and one VIRTIS-M dataset acquired on 27 April 2015, when the comet was at 1.76 au from the Sun in the inbound arc. No strong temporal correlation between total dust brightness and water production rates is found, although water is still the main driver of dust activity at this period in time. The observed increased water activity in the southern regions with strong seasonal variations with respect to model expectations is likely to be attributable to regional changes in volatile content or access to this. The best explanation for the drastic decrease in dust brightness when consolidated regions are illuminated is a quenched dust activity due to the high cohesion of surface material. These observations show that when 67P is approaching perihelion, the dust activity cannot be understood based on water-driven activity alone. This is in agreement with other modeling results on the seasonal evolution of the near-nucleus coma, which show that the correlation that was observed earlier in the mission between the observed dust coma and a modeled water coma from a homogeneously sublimating nucleus is significantly degraded (Shi et al. 2018b).