Open Access
Issue
A&A
Volume 670, February 2023
Article Number A159
Number of page(s) 10
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202244779
Published online 23 February 2023

© The Authors 2023

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

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

1 Introduction

46P/Wirtanen is a Jupiter-family comet (hereafter, JFC), with an orbital period of 5.5 yr, that was discovered on January 17, 1948 by Carl Wirtanen at the Lick Observatory (Jeffers 1948). It has been classified as an hyperactive comet, namely, it is part of a small group of comets whose activity levels are higher than expected based on the sizes of their nuclei. It was an initial target of the Rosetta mission and many studies have been performed to determine its nucleus size, shape, rotation period, albedo, and color, as well as its gas and dust activity (Farnham & Schleicher 1998; Lamy et al. 1998; Schulz et al. 1998; Bertaux et al. 1999; Crovisier et al. 2002). New radar observations have estimated the nucleus radius to be (0.70±0.05) km during its 2018 passage1.

The comet's return to perihelion was highly anticipated, as it made an unusually close approach to Earth (0.07 au) in December 2018 that was only about 30 times the distance to the moon, and with an excellent visibility from both hemispheres. This passage allowed for observations to be undertaken around the world and in great detail, with a broad set of ground-based and space telescopes at a wide range of wavelengths. This close approach provided a rare spatial resolution of about 50 km arcsec−1, which offered unique conditions to study the inner coma of a comet. Several results have been published on the characterization of 46P including measurements of the nucleus' rotation period (Farnham et al. 2021), the detection of mini-outbursts (Kelley et al. 2019, Kelley et al. 2021), measurement of the D/H ratio (Lis et al. 2019), abundances of several organic molecules observed at near-infrared wavelengths at different heliocentric distances (McKay et al. 2021; Khan et al. 2021; Roth et al. 2021), search for icy grains in the coma (Protopapa et al. 2021; Bonev et al. 2021), and coma features from optical imaging (Knight et al. 2021). Millimeter and submillimeter observations of 46P have reported on the detection of several complex organic molecules and showed that the comet is relatively rich in methanol, but relatively depleted in CO, CS, HNC, HNCO, and HCOOH (Biver et al. 2021; Roth et al. 2021b; Coulson et al. 2020; Bergman et al. 2022). Using both TRAPPIST-North (TN) and TRAPPIST-South (TS), we collected photometric observation over almost a year. 46P was visible for many hours during the night, allowing for long observing sequences for months during its return, passing near to the Sun.

In this work, we report on the activity, composition, and rotation period of comet 46P during its 2018 perihelion passage. An introduction and historical background of comet 46P is given in Sect. 1. We describe the observing circumstances and data reduction process in Sect. 2. In Sect. 3, we discuss the evolution of the gas and dust activity before and after perihelion and we compare the behavior of 46P to previous apparitions. The relative molecular abundances and their evolution with respect to the heliocentric distances are discussed in the same section. In Sect. 4, we measure the nucleus' rotation period and its variation around perihelion. In Sect. 5, we present the C and N isotopic ratios and the NH2 (and NH3) ortho-to-para ratio (OPR) and the green-to-red [OI] lines (G/R) ratio derived from high-resolution spectra. The summary and conclusions of this work are given in Sect. 6.

2 Observations and data analysis

2.1 TRAPPIST

We used both TN and TS telescopes (Jehin et al. 2011) to observe and follow comet 46P over almost a year. We started monitoring the comet at the beginning of August 2018 (rh = 1.88 au, pre-perihelion) until the end of March 2019 (rh = 1.70 au, post-perihelion). More than 2400 broad and narrow-band images of the comet were collected over 45 nights with TS and 40 nights with TN. 46P reached its perihelion on December 12, 2018 at 1.06 au from the Sun and only at 0.08 au from Earth. For the data reduction, we followed standard procedures using frequently updated master bias, flat, and dark frames. The sky contamination was removed and a flux calibration was performed using standard stars. In order to derive the production rates, we converted the flux of different gaseous species (OH, NH, CN, C3, and C2), measured through the HB narrow band cometary filters (Farnham et al. 2000), to the column densities using heliocentric distance and heliocentric velocity-dependent fluorescence efficiencies (Schleicher 2010, and references therein). To convert the column densities into production rates, we used a Haser model (Haser 1957) to fit the coma profile. The model adjustment was performed at a physical distance of 10 000 km from the nucleus. The scale lengths and g-factors of different molecules at 1 au were scaled by rh2$r_{\rm{h}}^{ - 2}$ (A'Hearn et al. 1995; Schleicher 2010). More details about the Haser model and its parameters are given in our previous works (Moulane et al. 2018, 2020). We derived the water-production rate from our the OH production rates using the empirical formula Q(H2O)=1.36rh0.5${\rm{Q}}\left( {{{\rm{H}}_{\rm{2}}}{\rm{O}}} \right) = 1.36\,r_{\rm{h}}^{ - 0.5}$ Q(OH), given in Cochran & Schleicher (1993) and Schleicher et al. (1998). We derived the Afρ parameter, a proxy for the dust production (A'Hearn et al. 1984), from the dust profiles using the HB cometary dust continuum BC, GC, and RC filters and the broadband R and I filters. We corrected the A(θ)fρ for the phase angle effect to obtain A(0)fρ using the phase function given in Schleicher (2007). The main uncertainties in the gas-production rates and Afρ come from the absolute flux calibration and the sky background subtraction. For the absolute flux calibration, we used an uncertainty of 5% on the extinction coefficients in all filters as seen in our long-term observations of the standard stars. This is almost negligible compared to the sky level uncertainty at lower airmass, but it becomes significant at high airmasses. We estimated the three-sigma-level uncertainty on the sky background value and use it to computed the error on production rates due to the sky subtraction. As a result, errors given in the following sections are a quadratic combination of sky background and extinction coefficient uncertainties.

2.2 UVES and ESPRESSO at 8-m ESO VLT

We obtained a high-resolution spectrum of comet 46P with the Ultraviolet-Visual Echelle Spectrograph (UVES) mounted on the Unit 2 telescope (UT2) of the Very Large Telescopes (VLT) at the European Southern Observatory (ESO) on December 9, 2018 (a week after perihelion, rh = 1.05 au and Δ = 0.09 au). We used the UVES standard setting DIC#1 346+580 covering the range 3030–3880 Å in the blue and 4760–6840 Å in the red. We used a 0.44″ wide slit, providing a resolving power of R ~ 80 000. We took two spectra at UT 0:59:52 (exposure time of 2300 s) and UT 1:19:54 (exposure time of 2300 s) as a start time.

The ESO UVES pipeline (Dekker et al. 2000) was used to reduce the spectra in the extended object mode, in which the spatial information is kept. The spectra were corrected for extinction and flux calibrated using the UVES master response curve provided by ESO. One-dimensional spectra were then extracted by averaging the 2D spectra with simultaneous cosmic ray rejection and then corrected for the Doppler shift due to the velocity of the comet with respect to the Earth. More details about the UVES data reduction are given in the UVES manual2. The dust-reflected sunlight was finally removed using a reference solar spectrum, BASS20003. A more detailed description of the steps for the UVES data reduction and the solar spectrum subtraction is given in Manfroid et al. (2009) and references therein.

A set of six spectra were also obtained with the new Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO) at the VLT (Pepe et al. 2021) over two nights on December 9 and 10, 2018. Exposure times ranged from 5400 s for the first spectrum to 6300 s for the following five spectra. We used the ultra-high-resolution mode of the instrument, providing a spectral resolution of 190 000 between 3800 and 7880 Å. In this mode, the spectrograph is fed by two 0.5″-size fibers: one centered on the object and another 7″ away usually used for simultaneous observation of the sky or a wavelength drift reference. However, in the case of the very extended coma of 46P, both fibers contained comet signal and provided an opportunity to probe different parts of the coma. Observations were executed in visitor mode, with good seeing (typically less than 1 arcsec) and the data were reduced using the ESPRESSO pipeline. which is publicly available from the ESO pipeline repository4.

thumbnail Fig. 1

OH, NH, CN, C2, and C3 logarithmic production rates of comet 46P as a function of the heliocentric distance and of time to perihelion. The vertical dashed line indicates the perihelion at 1.06 au on December 12, 2018. Open circles indicate pre-perihelion observations, filled symbols are used for post-perihelion observations. Negative signs in x-axis represents pre-perihelion data. The values of the production rates and their relative uncertainties are given in Table A.1.

3 Activity and composition

3.1 Gas-production rates and ratios

We followed the evolution of the production rates of different gaseous species over a year. The derived production rates for each species and the A(0)fρ parameter, as well as their uncertainties are given in Table A.1. Figure 1 shows the logarithmic evolution of OH, NH, CN, C2, and C3 production rates as a function of the distance to the Sun and time to perihelion. The gas activity of the comet was increasing slightly as the comet was approaching the Sun and decreased rapidly after perihelion, and more sharply for NH, OH and the dust. This fall-off could be due to the orientation of the active source regions on the surface of the comet toward the Earth before and after perihelion. The same behavior has been observed in the Lowell Discovery Telescope data (Knight et al. 2021). This explanation is supported by a study of the morphological features seen in 46P coma during this passage, reported by Knight et al. (2021). They show that 46P had large changes in both the sub-solar and sub-Earth latitudes during the apparition as well as changes in the pole orientation of the nucleus, assuming that there is no significant non-principal-axis rotation considered. It is also possible that the ice material present under the surface across the nucleus of the comet has different sublimation efficiency at the same heliocentric distances pre- versus post-perihelion. As mentioned above, we used narrow and broadband dust filters to estimate the dust production using the Afρ parameter introduced by A'Hearn et al. (1984), see Fig. 2. The behavior of the comet dust activity is similar to the gas on both sides of perihelion. This kind of asymmetry has been seen in many LPCs and JFCs.

3.2 Molecular abundances and dust-to-gas ratio

Using the production rates above, we derived the ratios with respect to CN and to OH as well as the dust-to-gas ratios. Figure 3 shows the evolution of the logarithm of the production rate ratios with respect to CN and to OH as well as ratios of A(0)fρ-to-Q(OH) and A(0)fρ-to-Q(CN) with time to perihelion. The C2/CN ratio decreases when the comet gets far away from the Sun (40 days post-perihelion). The abundances ratios did not show any significant variation on both sides of perihelion. This implies that 46P has an homogeneous composition along its orbit around the Sun. Comet 46P has been observed with many ground-based telescopes over the previous passages using narrow-band photometry, since it was an important target selected initially for the Rosetta mission in 1994. Table 1 summarizes the mean relative molecular abundances in 46P measured in 2018 (this work), in the 1991 apparition (A'Hearn et al. 1995), and in the 1997 apparition (Farnham & Schleicher 1998) compared to the database of comets given in Schleicher (2008). The measurements in the 1991, 1997 and 2018 apparitions show that 46P has a “typical” coma composition and a normal dust-to-gas ratio. Given the uncertainties on the ratios, no significant change is observed in the molecular abundances of 46P over the last four orbits. This shows that even an hyperactive comet, while outgassing a lot, does not change its main surface composition properties after several passages close to the Sun.

thumbnail Fig. 2

A(0)fρ parameter from the broad and narrow-band filters as a function of time to perihelion in days. Triangle symbols indicate data acquired with TN. Filled circles indicate TS.

thumbnail Fig. 3

Logarithm of the production rate and dust-to-gas ratios of comet 46P with respect to CN (Log10[Q(xx)/Q(CN)]) and to OH (Log10[Q(xx)/Q(OH)]) as a function of time to perihelion. Pre-perihelion values are represented with open circles and post-perihelion values with filled circles. The dust-to-gas ratio is expressed in cm second per molecule. The horizontal dashed line indicates the mean ratios of typical comets given in Schleicher (2008).

Table 1

Mean relative molecular abundances of comet 46P for the perihelion 1991–2018 compared to typical and depleted comets.

3.3 Parent and daughter molecules

In order to investigate the origin of the radicals observed in 46P's coma, we compared our daughter species abundances with parent molecular abundances obtained at Infrared wavelengths. Table 2 summarizes the abundances for some radicals observed in the optical compared to their possible parent abundances observed at NIR wavelengths (Roth et al. 2021a; Bonev et al. 2021; McKay et al. 2021; Khan et al. 2021), measured during the same period of time. We compared these abundances at two epochs, perihelion (rh = 1.06 au) and post-perihelion (rh = 1.12), when the infrared observations are available. Our CN/OH ratios are in agreement with the HCN/H2O ratios derived directly from IR observations around the same dates, namely, Dec. 16−23, 2018 (see Table 2). In addition, the HCN production rates (1.06±0.03) × 1025 molec s−1 obtained from millimeter (Biver et al. 2021) and (1.60±0.10) × 1025 molec s−1 from sub-millimeter observations (Bergman et al. 2022) show an agreement with our Q(CN) values obtained during the same nights of Dec. 15 and 20. This indicates that the observed CN abundances are consistent with HCN being the dominant parent. The C2 abundance is much higher than C2H2 which suggests that C2 could be dissociated from both C2H2 and C2H6. It has been shown that C2 is still linked to both molecules even at large heliocentric distances ≥3 au in comet Hale-Bopp (Helbert et al. 2005). We looked for a C3 parent but none has been identified in 46P at infrared wavelengths pointing to a C3 production coming from chemical reactions in the coma. Based on direct detections of NH3 lines, the mixing ratio of NH3/H2O ~ 0.50 in 46P (See Table 2) is lower than the NH/OH ~ 0.72 obtained with TRAPPIST. This result suggests that NH is not dissociated directly from NH3, but obviously via NH2 which has been demonstrated for several comets (Shinnaka et al. 2011, Shinnaka et al. 2016).

Table 2

Abundances ratios derived from the optical and infrared data of comet 46P in its 2018 apparition.

thumbnail Fig. 4

Water-production rate of comet 46P as a function of time from perihelion during the 2018 apparition. Our measurements and those from Knight et al. (2021) were derived from the OH production rates. Combi et al. (2020) measurements were derived from H Lyman-α spectra using the SWAN instrument on SOHO. iSHELL/IRTF and McKay et al. (2021) data were derived from NIR spectra. Lis et al. (2019) data were derived from far-infrared spectra using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA).

3.4 Water-production rates and active area

We derived the water-production rate from OH using an empirical relationship (see Sect. 2) proposed by Cochran & Schleicher (1993) and Schleicher et al. (1998). Figure 4 shows a comparison of our water-production rates for the 2018 apparition with other measurements derived at different wavelengths using various techniques. First, we have a very good agreement between the TRAPPIST data (red circles) and the Lowell Observatory data (black squares) on both sides of perihelion, with both using the same technique based on narrow-band filters and the same model to derive the water-production rates. We also have a good agreement with H2O measured directly from the NIR spectra from iSHELL at the NASA-IRTF and NIRSPEC-2 at Keck (McKay et al. 2021; Roth et al. 2021a; Khan et al. 2021; Bonev et al. 2021). The data from Combi et al. (2020) who derived the water-production rates from the hydrogen Lyman-α emission observed by the SWAN instrument on board SOHO show an offset on both sides of perihelion. This disagreement is probably the result of using different techniques and models. It is also tricky to compare these results as SWAN/SOHO instrument threshold only allows for Lyman-α emission to be obtained within a small range of heliocentric distances. SWAN's values have rather large error-bars and dispersion among their measurements during this passage and the same for the previous ones (see Fig. 5). We have seen this systematic difference also for comet 21P/Giacobini-Zinner (Moulane et al. 2020).

As 46P was the original target selected for Rosetta mission, many observations were performed in the previous passages in order to characterize the comet. Here, we compare the water-production rates obtained in the previous passages 1991 and 1997 with those obtained in 2008 and 2018 (see Fig. 5). Besides the offset between narrow-band and spectroscopic data, the comet s activity was similar in 1991 and 1997 but 30% lower in 2008–2018 apparitions. We did not include the 2002 observations in this comparison as they were dominated by large outbursts5 after perihelion (Combi et al. 2020). The maximum water-production rate in 2018 was measured ten days pre-perihelion with Q(H2O) = (8.00±0.30)×1027 molec s−1, while it was (7.40±0.20) × 1027 molec s−1 a month after perihelion in 2008, (1.23±0.10) × 1028 molec s−1 a week before perihelion in 1997, and 1.25 × 1028 molec s−1 three weeks after perihelion in 1991 (Farnham & Schleicher 1998; Fink et al. 1998; Bertaux et al. 1999; Crovisier et al. 2002; Kobayashi & Kawakita 2010). There is a decrease in the gas-production rates by a factor of two over the last two decades (from 1997–2018), taking into account the uncertainties and systematic offset between different techniques. Our results are similar to those reported by others who found a similar decrease rate of the water-production rate over several previous apparitions (Combi et al. 2020; Knight et al. 2021).

Comet 46P is considered a hyperactive comet as it is emitting more water than it would be expected from its nucleus radius. Based on the standard water vaporization model (Cowan & A'Hearn 1979)6, we estimate the evolution of the active area of 46P during our observational campaign in 2018 by modeling the water-production rate as a function of the heliocentric distances. We assumed a bond albedo of 5%, 100% infrared emissivity, and the rotational pole pointed at the Sun (see, e.g., A'Hearn et al. 1989; McKay et al. 2018, 2019). We found that the active area of 46P varied from ~0.9 km2 at 1.52 au pre-perihelion, reached a maximum of ~2.5 km2 at 1.06 au during the perihelion passage, and decreased to ~0.7 km2 at 1.53 au post-perihelion. Based on a radius of ~0.7 km, the active fraction of the nucleus was ~14, ~40 and ~11%, respectively. These values are much larger than the active fraction of less than 3% derived for the majority of JFCs that have radius measurements (A'Hearn et al. 1995). Since the comet had been observed in its previous passages, we compared the active fraction of its nucleus at perihelion in the past. Based on the water-production rates derived from the empirical conversion of Haser OH production rates (Schleicher et al. 1998; Farnham et al. 2007), we estimated the active fractions of its nucleus to be ~72% at rh = 1.11 au in 1991, ~65% at rh = 1.07 au in 1997, ~37% at rh = 1.05 au in 2008 and ~40% at rh = 1.06 au in 2018. It is clear that there is a decrease in the active fraction of the nucleus with time. Comet 103P/Hartley 2 was the first known “hyperactive” comet and it was visited by the Deep Impact extended Investigation mission, demonstrating that this hyperactivity was produced by sublimation of icy grains in the coma (A'Hearn et al. 2011; Protopapa et al. 2014). Several authors reported on the presence of an extended source of water ice sublimation in 46P coma in the form of small grains or large chunks. Bonev et al. (2021) provided indirect arguments for the presence of an extended source of water vapor in the coma of 46P at NIR spectra, but they did not provide information on the properties of the ice that produces the water vapor, such as particle size and ice-to-dust ratio, while Protopapa et al. (2021) and Kareta et al. (priv. comm.) did not detect any water ice absorption features in their NIR spectra.

thumbnail Fig. 5

Water-production rates of comet 46P as a function of time from perihelion during 1991 (orange), 1997 (black), 2008 (blue), and 2018 (in red) apparitions. Fink et al. (1998) derived the water-production rate from the [OI](1D) forbidden line doublet emission in the optical, Bertaux et al. (1999) values were derived from the H Ly-α emission observed by the SWAN instrument on board SOHO and Crovisier et al. (2002) values were derived from the OH emission observed at the radio wavelengths.

Table 3

Rotation period of 46P nucleus at different epochs.

4 Rotation period and CN coma morphology

We took the opportunity provided by this very close approach of 46P to investigate the rotation period of its nucleus, using long CN series of images collected during the same night and at different epochs on both sides of the perihelion. Thanks to the visibility of the comet in both hemispheres, we collected long series of CN images for many hours on several nights with both telescopes. These series were taken over 12 nights with TS and 8 nights with TN. We determined the rotation period from CN light curves. We tested different size apertures and we finally used 20″ which gave the largest amplitude in the light-curve. We phased multiple nights data to construct more extensive light curves and searched for the best alignment of the overlapping segments. We derived a rotation period of 9.12±0.05 h in late November and early December (pre-perihelion), 9.18±0.05 h between December 7 and 10 (perihelion time) and 9.00±0.04 h in mid January 2019 (post-perihelion). These results are in agreement (within the error bars) with those derived by Farnham et al. (2021) from the CN coma morphology and photometric series (see Table 3). Given the uncertainties on the measurements, we did not detect any significant change in the rotation period of the nucleus on both sides of the perihelion. Figure 6 shows the CN light curves at different epochs, showing two asymmetric maxima and minima, which could be due to two active areas on the surface of the nucleus.

We searched for CN coma features using the rotational filter technique as described in our previous work (Moulane et al. 2018). Figure 7 shows the detection and evolution of a CN jet at different epochs and on both sides of perihelion. We did not detect any clear rotation of this jet with time in our images for many hours during the same night. This could be due to an orientation of the jet towards the Earth. Around November 16, another small jet appeared in the CN images but it disappeared in the following nights. This could suggest that the nucleus might be in a non-principal-axis (NPA) rotational state, as was previously mentioned by Samarasinha et al. (1996). The orientation of the jet remains the same before perihelion, while it changes by ~95 degrees counterclockwise after perihelion because of the viewing geometry of the comet with respect to Earth. Farnham et al. (2021) used the 4.3-m Lowell Discovery Telescope to investigate the CN coma morphology of 46P. They were able to detect two spiral jets in the coma. The first one appears to have been active (at varying levels) throughout most of a rotation, corresponding to the bright jet detected in our data. The second seems to turn on and off with time and it could correspond to the small jet that appears around mid November in our images. Compared to the previous apparitions, Lamy et al. (1998) derived a rotation period of the nucleus of 6.0±0.3 using data from HST, while Meech et al. (1997) found a possible rotation of 7.6 h and a very small amplitude variation of 0.09 mag using ground-based telescope during the 1997 apparition. Based on these results, we conclude that the rotation of the 46P nucleus might have slowed down by 16–45% over the last four orbits. These changes might be due to some unusual activity in the previous apparitions. There is some evidence for potentially significant outburst activity during the 2002 apparition (Combi et al. 2020) that could explain these changes.

thumbnail Fig. 6

Phased CN light curves for the nucleus rotation period of comet 46P at different epochs. Time origin for the phase is JD 2458440 on Nov. 17, 2018. Data combined for (a) Nov. 23, 28 and Dec. 1, 2018 preperihelion, (b) Dec. 7–10, 2018 at perihelion and for (c) Jan. 12 and 15, 2019 post-perihelion. The red curve shows the best fit of the data.

5 Optical high-resolution spectra

5.1 The 12C13C and 14N/15N isotopic ratios

We obtained high-resolution observations with UVES/VLT (Dekker et al. 2000) of comet 46P on December 9, 2018 (rh = 1.05 au and Δ = 0.09 au). Using the 12C14N B-X(0,0) band, we estimated the 12C/13C and 14N/15N isotopic ratios of 46P. We used a CN fluorescence model to create synthetic spectra of 13C14N, 12C15N, and 12C14N. More details of the model are given in Manfroid et al. (2009). Figure 8 shows the observed CN spectrum compared to the synthetic one made under the same observing conditions. The ratios found for 12C/13C and 14N/15N are 100±20 and 150±30, respectively. These values are consistent with those of about 20 comets with different dynamical origins, 91.0±3.6 and 147.8±5.7 for 12C/13C and 14N/15N, respectively (Manfroid et al. 2009; Bockelée-Morvan et al. 2015).

thumbnail Fig. 7

Evolution of the CN jets of comet 46P along its orbit around the Sun. The images are oriented north (N) up and east (E) right.

thumbnail Fig. 8

Observed and synthetic CN isotopologues spectra of the R branch of the B-X (0, 0) violet band in comet 46P.

5.2 The NH2 and NH3 ortho-to-para ratios

From the UVES spectra, we measured the ortho-to-para (OPR, hereafter) abundance ratio of NH2 using the three robivronic emission bands (0,7,0), (0,8,0), and (0,9,0) following the method described in Shinnaka et al. (2011). The derived OPRs of NH2 and of its parent molecule NH3 are listed for each band in Table 4, and they have average values of 3.31 ± 0.03 and 1.19 ± 0.02, respectively. A nuclear spin temperature (Tspin) for ammonia of 28 ± 1 Κ was derived. The 46P value is consistent with typical values measured in comets (see Fig. 9), which is a possible cosmogonic indicator linked to the formation temperature of the molecule.

We would like to point out that recent laboratory experiments show that the OPR of water does not retain the memory of its formation temperature in the natal molecular cloud or in the solar nebula 4.6 Gyrs ago (Hama et al. 2011, 2016; Hama & Watanabe 2013). The OPR of water just after its sublimation from the solid phase is set to a statistical weight ratio (i.e., 3 for water). It is likely that this is also the case for ammonia (but the statistical weights ratio is unity for ammonia) although no experimental studies have been reported so far. The OPRs of cometary volatiles might have been modified by some kind of ortho-to-para conversion processes in the inner coma such as water clusters, or by other catalyst activities of dust crust surfaces of the nucleus (Shinnaka et al. 2016). Thus, the observed OPRs might be diagnostic of the physico-chemical conditions in the innermost coma or beneath the surface.

Table 4

Derived NH2 and NH3 OPRs in comet 46P.

thumbnail Fig. 9

NH3 Tspin in comets of various dynamical origin. The red filled circle is the value of 46P obtained with UVES at VLT (this work).

5.3 The [OI] G/R ratio

Using the very high-resolution ESPRESSO spectra (Pepe et al. 2021), we were able to resolve the telluric and cometary forbidden oxygen lines even with a Doppler shift as small as 5 km s−1. We detected the cometary [OI] lines at 557.73, 630.03, and 636.38 nm in all six spectra (see Fig. 10) and used them to measure the ratio between the green line (557.73 nm) and the sum of the two red lines (630.03 and 636.38 nm), commonly referred to as the G/R ratio. Due to the low geocentric velocity of the comet at the time of our observations (between −5.4 and −4.5 km s−1), the cometary and telluric lines were slightly blended, even at the very high resolution of ESPRESSO. We used two Gaussian to fit both lines simultaneously. We performed the measurements for the center and the sky fiber separately for the six spectra and subsequently averaged the values for each fiber. The uncertainty of the measurement is the standard deviation of the measurements made of six spectra. We measured a G/R ratio of 0.23±0.02 in the central fiber and of 0.05±0.01 only in the sky fiber (7″ away). Both measurements are in agreement with what has been measured in several comets at the same nucleo-centric projected distances (Decock et al. 2015). The difference of G/R ratio at different distances from the nucleus is mainly due to the collisional quenching of O(1S) and O(1D) by water molecules in the inner coma. This fits with the general behavior observed in several comets, where the G/R ratio decreases monotonically with the projected distance to the nucleus (Decock et al. 2015). The intensity of these lines can be used to derive the CO/CO2/H2O abundances ratios and, hence, the C/O ratio, which can provide a diagnostic of solar nebula chemistry in the comet-forming region (Öberg et al. 2011). The detail analysis and modeling of the [OI] lines are out of the scope of this paper.

thumbnail Fig. 10

Three forbidden oxygen [OI] lines in comet 46P, from ESPRESSO observations, at 557.73 nm (left) for the green line, and 630.03 nm (middle) and 636.38 nm (right) for the red doublet. The cometary oxygen lines are identified by tick marks and the telluric lines marked by dashed tick. The y-axis is in arbitrary units for all three panels. The red line is a fit for one of the six spectra taken in the centered fiber.

6 Summary and conclusions

We conducted an extensive monitoring of the Jupiter Family hyperactive comet 46P on both sides of perihelion with TRAPPIST telescopes. These observations led us to an overview of the comet activity along its journey around the Sun in 2018–2019. Using narrow and broadband filters with the TRAPPIST telescopes, we derived different quantities related to the activity of the comet such as the gas-production rates, the A(0)fρ parameter, the rotation period of the nucleus, and the morphology of the coma. Comet 46P shows an asymmetric activity with respect to perihelion, steeper decrease after perihelion for OH, NH, and dust (Afρ parameter). This might be due to the spin axis orientation or the distribution of activity on the comet's surface. The comet has exhibited a decrease in its activity from the 1997–2018 apparition by about 30%. This implies a decrease of the nucleus's active fraction to ~40% which confirms its hyperactivity compared to most of JFCs. The molecular abundances relative to CN and OH remained the same as a function of the heliocentric distance, showing that 46P has a typical composition similar to what has been found in previous apparitions. These derived abundances in the optical showed a good correlation with their possible mother molecules measured in the IR (HCN, C2H2, and NH3) and sub-millimeter (HCN), confirming the link between different molecules observed in several comets. Given the close approach of the comet, we obtained high spatial and high resolution spectra for the comet around the perihelion. Using UVES spectra, we measured typical 12C/13C and 14N/15N isotopic ratios of 100±20 and 150±30, and average NH2 OPR of 3.31±0.03 and 1.16±0.03 for NH3, with a nuclear spin temperature of 28±1 K. The ammonia OPR was found equal to 1.19±0.03, corresponding to a spin temperature of 27±1 K. Using very-high-resolution ESPRESSO spectra, we measured a forbidden oxygen lines ratio G/R of 0.23±0.02. These measurements are in agreement with those found for several comets of different dynamical types and origins and do not reveal any peculiarity with regard to the hyperactive comet 46P.

Acknowledgements

This publication makes use of data products from TRAPPIST project. TRAPPIST-South is funded by the Belgian National Fund for Scientific Research (F.R.S.-FNRS) under grant PDR T.0120.21. TRAPPIST-North is funded by the University of Liège, and performed in collaboration with Cadi Ayyad University of Marrakesh. E.J. and D.H. are FNRS Senior Research Associates. J.M. is Honorary Research Director of the FNRS. We thank NASA, David Schleicher and the Lowell Observatory for the loan of a set of HB comet filters. UVES and ESPRESSO observations were collected with VLTs telescopes at the European Southern Observatory Paranal in Chile under programs 0102.C-0438(A) and 0102.C-0438(B), respectively.

Appendix A Gas-production rates and A(0)fρ parameter of comet 46P/Wirtanen

Table A.1

Gas-production rates and A(0)fρ parameter of comet 46P/Wirtanen measured with TRAPPIST-South (TS) and -North (TN). The A(0)fρ values are computed at 10 000 km from the nucleus and corrected for the phase angle effect.

References

  1. A’Hearn, M. F., Schleicher, D. G., Millis, R. L., Feldman, P. D., & Thompson, D. T. 1984, AJ, 89, 579 [Google Scholar]
  2. A’Hearn, M. F., Campins, H., Schleicher, D. G., & Millis, R. L. 1989, ApJ, 347, 1155 [CrossRef] [Google Scholar]
  3. A’Hearn, M. F., Millis, R. C., Schleicher, D. O., Osip, D. J., & Birch, P.V. 1995, Icarus, 118, 223 [CrossRef] [Google Scholar]
  4. A’Hearn, M. F., Belton, M. J., Delamere, W. A., et al. 2011, Science, 332, 1396 [CrossRef] [Google Scholar]
  5. Bergman, P., Lerner, M. S., Olofsson, A. O. H., et al. 2022, A&A, 660, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bertaux, J., Costa, J., Mäkinen, T., et al. 1999, Planet. Space Sci., 47, 725 [NASA ADS] [CrossRef] [Google Scholar]
  7. Biver, N., Bockelée-Morvan, D., Boissier, J., et al. 2021, A&A, 648, A49 [EDP Sciences] [Google Scholar]
  8. Bockelée-Morvan, D., Calmonte, U., Charnley, S., et al. 2015, Space Sci. Rev., 197, 47 [Google Scholar]
  9. Bonev, B. P., Russo, N. D., DiSanti, M. A., et al. 2021, Planet. Sci. J., 2, 45 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cochran, A. L., & Schleicher, D. G. 1993, Icarus, 105, 235 [Google Scholar]
  11. Combi, M. R., Mäkinen, T., Bertaux, J. L., et al. 2020, Planet. Sci. J., 1, 72 [Google Scholar]
  12. Coulson, I. M., Liu, F.-C., Cordiner, M. A., et al. 2020, AJ, 160, 182 [Google Scholar]
  13. Cowan, J. J., & A’Hearn, M.F. 1979, Moon Planets, 21, 155 [NASA ADS] [CrossRef] [Google Scholar]
  14. Crovisier, J., Colom, P., Gérard, E., Bockelée-Morvan, D., & Bourgois, G. 2002, A&A, 393, 1053 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Decock, A., Jehin, E., Rousselot, P., et al. 2015, A&A, 573, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Dekker, H., D'Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, in SPIE Proceedings, ed. M. Iye, & A.F.M. Moorwood (SPIE) [Google Scholar]
  17. Farnham, T. L., & Schleicher, D. G. 1998, A&A, 335, L50 [NASA ADS] [Google Scholar]
  18. Farnham, T.L., Schleicher, D. G., & A’Hearn, M. F. 2000, Icarus, 147, 180 [NASA ADS] [CrossRef] [Google Scholar]
  19. Farnham, T. L., Samarasinha, N. H., Mueller, B. E. A., & Knight, M. M. 2007, AJ, 133, 2001 [NASA ADS] [CrossRef] [Google Scholar]
  20. Farnham, T.L., Knight, M.M., Schleicher, D.G., et al. 2021, Planet. Sci. J., 2, 7 [NASA ADS] [CrossRef] [Google Scholar]
  21. Fink, U., Hicks, M. D., Fevig, R. A., & Collins, J. 1998, A&A, 335, L37 [NASA ADS] [Google Scholar]
  22. Hama, T., & Watanabe, N. 2013, Chem. Rev., 113, 8783 [Google Scholar]
  23. Hama, T., Watanabe, N., Kouchi, A.M.Y., 2011, ApJ, 738, L15 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hama, T., Kouchi, A., & Watanabe, N. 2016, Science, 351, 65 [NASA ADS] [CrossRef] [Google Scholar]
  25. Haser, L. 1957, Bull. Soc. Roy. Sci. Liège, 43, 740 [Google Scholar]
  26. Helbert, J., Rauer, H., Boice, D. C., & Huebner, W. F. 2005, A&A, 442, 1107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Jeffers, H. M. 1948, Lick Observ. Bull., 19, 189 [Google Scholar]
  28. Jehin, E., Gillon, M., Queloz, D., et al. 2011, The Messenger, 145, 2 [NASA ADS] [Google Scholar]
  29. Kelley, M. S. P., Bodewits, D., Ye, Q., et al. 2019, Res. Notes Am. Astron. Soc., 3, 126 [Google Scholar]
  30. Kelley, M. S. P., Farnham, T. L., Li, J.-Y., et al. 2021, Planet. Sci. J., 2, 131 [NASA ADS] [CrossRef] [Google Scholar]
  31. Khan, Y., Gibb, E. L., Bonev, B. P., et al. 2021, Planet. Sci. J., 2, 20 [NASA ADS] [CrossRef] [Google Scholar]
  32. Knight, M. M., Schleicher, D. G., & Farnham, T. L. 2021, Planet. Sci. J., 2, 104 [NASA ADS] [CrossRef] [Google Scholar]
  33. Kobayashi, H., & Kawakita, H. 2010, PASJ, 62, 1025 [NASA ADS] [CrossRef] [Google Scholar]
  34. Lamy, P. L., Toth, I., Jorda, L., Weaver, H. A., & A’Hearn, M. 1998, A&A, 335, L25 [NASA ADS] [Google Scholar]
  35. Lis, D. C., Bockelée-Morvan, D., Güsten, R., et al. 2019, A&A, 625, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Manfroid, J., Jehin, E., Hutsemékers, D., et al. 2009, A&A, 503, 613 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. McKay, A. J., Cochran, A. L., DiSanti, M. A., et al. 2018, Icarus, 309, 1 [NASA ADS] [CrossRef] [Google Scholar]
  38. McKay, A. J., DiSanti, M. A., Kelley, M. S. P., et al. 2019, AJ, 158, 128 [NASA ADS] [CrossRef] [Google Scholar]
  39. McKay, A. J., DiSanti, M. A., Cochran, A. L., et al. 2021, Planet. Sci. J., 2, 21 [NASA ADS] [CrossRef] [Google Scholar]
  40. Meech, K. J., Bauer, J. M., & Hainaut, O. R. 1997, A&A, 326, 1268 [NASA ADS] [Google Scholar]
  41. Moulane, Y., Jehin, E., Opitom, C., et al. 2018, A&A, 619, A156 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Moulane, Y., Jehin, E., Rousselot, P., et al. 2020, A&A, 640, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011, ApJ, 743, L16 [Google Scholar]
  44. Pepe, F., Cristiani, S., Rebolo, R., et al. 2021, A&A, 645, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Protopapa, S., Sunshine, J. M., Feaga, L. M., et al. 2014, Icarus, 238, 191 [NASA ADS] [CrossRef] [Google Scholar]
  46. Protopapa, S., Kelley, M. S. P., Woodward, C. E., & Yang, B. 2021, Planet. Sci. J. 2, 176 [NASA ADS] [CrossRef] [Google Scholar]
  47. Roth, N. X., Bonev, B. P., DiSanti, M. A., et al. 2021a, Planet. Sci. J., 2, 54 [Google Scholar]
  48. Roth, N. X., Milam, S. N., Cordiner, M. A., et al. 2021b, Planet. Sci. J., 2, 55 [NASA ADS] [CrossRef] [Google Scholar]
  49. Samarasinha, N. H., Mueller, B. E., & Belton, M. J. 1996, Planet. Space Sci., 44, 275 [NASA ADS] [CrossRef] [Google Scholar]
  50. Schleicher, D. G. 2007, Icarus, 191, 322 [NASA ADS] [CrossRef] [Google Scholar]
  51. Schleicher, D. G. 2008, AJ, 136, 2204 [Google Scholar]
  52. Schleicher, D. G. 2010, AJ, 140, 973 [Google Scholar]
  53. Schleicher, D. G., & Osip, D. J. 2002, Icarus, 159, 210 [Google Scholar]
  54. Schleicher, D. G., Millis, R. L., & Birch, P.V. 1998, Icarus, 132, 397 [NASA ADS] [CrossRef] [Google Scholar]
  55. Schulz, R., Arpigny, C., Manfroid, J., et al. 1998, A&A, 335, L46 [NASA ADS] [Google Scholar]
  56. Shinnaka, Y., Kawakita, H., Kobayashi, H., et al. 2011, ApJ, 729, 81 [Google Scholar]
  57. Shinnaka, Y., Kawakita, H., Jehin, E., et al. 2016, MNRAS, 462, S195 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Mean relative molecular abundances of comet 46P for the perihelion 1991–2018 compared to typical and depleted comets.

Table 2

Abundances ratios derived from the optical and infrared data of comet 46P in its 2018 apparition.

Table 3

Rotation period of 46P nucleus at different epochs.

Table 4

Derived NH2 and NH3 OPRs in comet 46P.

Table A.1

Gas-production rates and A(0)fρ parameter of comet 46P/Wirtanen measured with TRAPPIST-South (TS) and -North (TN). The A(0)fρ values are computed at 10 000 km from the nucleus and corrected for the phase angle effect.

All Figures

thumbnail Fig. 1

OH, NH, CN, C2, and C3 logarithmic production rates of comet 46P as a function of the heliocentric distance and of time to perihelion. The vertical dashed line indicates the perihelion at 1.06 au on December 12, 2018. Open circles indicate pre-perihelion observations, filled symbols are used for post-perihelion observations. Negative signs in x-axis represents pre-perihelion data. The values of the production rates and their relative uncertainties are given in Table A.1.

In the text
thumbnail Fig. 2

A(0)fρ parameter from the broad and narrow-band filters as a function of time to perihelion in days. Triangle symbols indicate data acquired with TN. Filled circles indicate TS.

In the text
thumbnail Fig. 3

Logarithm of the production rate and dust-to-gas ratios of comet 46P with respect to CN (Log10[Q(xx)/Q(CN)]) and to OH (Log10[Q(xx)/Q(OH)]) as a function of time to perihelion. Pre-perihelion values are represented with open circles and post-perihelion values with filled circles. The dust-to-gas ratio is expressed in cm second per molecule. The horizontal dashed line indicates the mean ratios of typical comets given in Schleicher (2008).

In the text
thumbnail Fig. 4

Water-production rate of comet 46P as a function of time from perihelion during the 2018 apparition. Our measurements and those from Knight et al. (2021) were derived from the OH production rates. Combi et al. (2020) measurements were derived from H Lyman-α spectra using the SWAN instrument on SOHO. iSHELL/IRTF and McKay et al. (2021) data were derived from NIR spectra. Lis et al. (2019) data were derived from far-infrared spectra using the GREAT spectrometer aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA).

In the text
thumbnail Fig. 5

Water-production rates of comet 46P as a function of time from perihelion during 1991 (orange), 1997 (black), 2008 (blue), and 2018 (in red) apparitions. Fink et al. (1998) derived the water-production rate from the [OI](1D) forbidden line doublet emission in the optical, Bertaux et al. (1999) values were derived from the H Ly-α emission observed by the SWAN instrument on board SOHO and Crovisier et al. (2002) values were derived from the OH emission observed at the radio wavelengths.

In the text
thumbnail Fig. 6

Phased CN light curves for the nucleus rotation period of comet 46P at different epochs. Time origin for the phase is JD 2458440 on Nov. 17, 2018. Data combined for (a) Nov. 23, 28 and Dec. 1, 2018 preperihelion, (b) Dec. 7–10, 2018 at perihelion and for (c) Jan. 12 and 15, 2019 post-perihelion. The red curve shows the best fit of the data.

In the text
thumbnail Fig. 7

Evolution of the CN jets of comet 46P along its orbit around the Sun. The images are oriented north (N) up and east (E) right.

In the text
thumbnail Fig. 8

Observed and synthetic CN isotopologues spectra of the R branch of the B-X (0, 0) violet band in comet 46P.

In the text
thumbnail Fig. 9

NH3 Tspin in comets of various dynamical origin. The red filled circle is the value of 46P obtained with UVES at VLT (this work).

In the text
thumbnail Fig. 10

Three forbidden oxygen [OI] lines in comet 46P, from ESPRESSO observations, at 557.73 nm (left) for the green line, and 630.03 nm (middle) and 636.38 nm (right) for the red doublet. The cometary oxygen lines are identified by tick marks and the telluric lines marked by dashed tick. The y-axis is in arbitrary units for all three panels. The red line is a fit for one of the six spectra taken in the centered fiber.

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.