# Photometry and high-resolution spectroscopy of comet 21P/Giacobini-Zinner during its 2018 apparition

Y. Moulane, E. Jehin, P. Rousselot, J. Manfroid, Y. Shinnaka, F. Pozuelos, D.

Hutsemékers, C. Opitom, Bin Yang, Z. Benkhaldoun

► To cite this version:

Y. Moulane, E. Jehin, P. Rousselot, J. Manfroid, Y. Shinnaka, et al., Photometry and high-resolution spectroscopy of comet 21P/Giacobini-Zinner during its 2018 apparition. Astronomy and Astrophysics
- A&A, 2020, 640, pp.A54. 10.1051/0004-6361/202037997. hal-02962955

HAL Id: hal-02962955

https://hal.science/hal-02962955

Submitted on 11 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Photometry and high-resolution spectroscopy of comet 21P/Giacobini-Zinner during its 2018 apparition*

Y. Moulane<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, E. Jehin<fnref n="2" />, P. Rousselot<fnref n="4" />, J. Manfroid<fnref n="2" />, Y. Shinnaka<fnref n="5" />, F. J. Pozuelos<fnref n="2" />, D. Hutsemékers<fnref n="2" />, C. Opitom<fnref n="1" />,<fnref n="6" />, B. Yang<fnref n="1" />, and Z. Benkhaldoun<fnref n="3" />

Received 20 March 2020 / Accepted 8 June 2020

## ABSTRACT

We report on photometry and high-resolution spectroscopy of the chemically peculiar Jupiter-family comet (hereafter JFC) 21P/Giacobini-Zinner. Comet 21P is a well-known member of the carbon-chain-depleted family, but also displays a depletion of amines. We continuously monitored the comet over more than seven months with the two TRAPPIST telescopes (TN and TS), covering a wide heliocentric distance range from 1.60 au inbound to 2.10 au outbound with a perihelion at 1.01 au on September 10, 2018. We computed and followed the evolution of the dust- (represented by $A\phi\rho$) and gas-production rates of the daughter species OH, NH, CN, $C_3$, and $C_2$ and their relative abundances to OH and to CN over the cometary orbit. We compared them to those measured in the previous apparitions. The activity of the comet and its water production rate reached a maximum of $(3.72 \pm 0.07) \times 10^{28} \text{ mol s}^{-1}$ on August 17, 2018 ($r_h = 1.07$ au), 24 days before perihelion. The peak value of $A(0)\phi\rho$ was reached on the same date (1646 ± 13) cm in the red filter. Using a sublimation model for the nucleus, we constrained the active surface of the nucleus using the slow-rotator model. The abundance ratios of the various species are remarkably constant over a wide range of heliocentric distances before and after perihelion, showing a high level of homogeneity of the ices in the surface of the nucleus. The behaviour and level of the activity of the comet is also remarkably similar over the last five orbits. In the coma dust colour, 21P shows reflective gradients similar to JFCs. We obtained a high-resolution spectrum of 21P with UVES at ESO Very Large Telescope one week after perihelion. Using the CN B-X (0, 0) violet band, we measured $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N isotopic ratios of $100 \pm 10$ and $145 \pm 10$, respectively, both in very good agreement with the ratios commonly found in comets. We measured an ortho-para abundance ratio of $NH_3$ of $1.16 \pm 0.02$, corresponding to a nuclear spin temperature of $T_{\text{spin}} = 27 \pm 1$ K, which is similar to other comets. While the abundance ratios in the gaseous coma reveal a peculiar composition, the isotopic and ortho-to-para ratios appear entirely normal. We performed a dynamical simulation of 21P and found that it is likely a young member of the JFC population. We favour a pristine composition scenario to explain the chemical peculiarities of this comet.

**Key words.** comets: general – comets: individual: 21P/Giacobini-Zinner – techniques: photometric – techniques: spectroscopic

## 1. Introduction

Comet 21P/Giacobini-Zinner (hereafter 21P) is a Jupiter-family comet (JFC) with a short period of 6.5 yr. 21P was discovered in 1900 by Michel Giacobini and was rediscovered by Ernst Zinner in 1913¹. After its discovery, 21P was observed in most of its apparitions, and many photometric and spectroscopic measurements were reported. In September 1985, 21P was the first comet visited by the International Cometary Exporer (ICE) spacecraft to study the interaction between the solar wind and the cometary atmosphere (Von Rosenvinge et al. 1986; Scarf et al. 1986). 21P is also known as the parent body of the Draconids meteor shower (Beech 1986; Egal et al. 2019). Many spectroscopic-photometric studies at various wavelength ranges have been performed since its discovery (Schleicher et al. 1987; Cochran & Barker 1987;

Fink & Hicks 1996; Weaver et al. 1999; Lara et al. 2003; Combi et al. 2011), including production rate measurements and atomic and molecular abundances. 21P is the prototype of depleted comets in $C_2$ and $C_3$ with respect to CN and to OH (Schleicher et al. 1987; A'Hearn et al. 1995). The relative abundances of $C_2$ and $C_3$ are about five and ten times lower than those measured in typical comets (A'Hearn et al. 1995). 21P was found to also be depleted in NH (Schleicher et al. 1987; Kiselev et al. 2000) and $NH_2$ (Konno & Wyckoff 1989; Beaver et al. 1990; Fink 2009).

The physical properties of 21P were investigated in the previous apparitions. Its nucleus size is not well determined; the average estimate of its radius is about 1–2 km (Tancredi et al. 2000; Królikowska et al. 2001; Pittichová et al. 2008). Its rotation period is not well constrained; a wide range from 9.5 to 19 h was estimated (Leibowitz & Brosch 1986). The activity of 21P has shown an asymmetric light curve with respect to perihelion in the previous apparitions. In 1985, the production rates pre-perihelion were twice higher than post-perihelion at heliocentric distances of 1.0–1.5 au (Schleicher et al. 1987). The highest gas and dust

<footnote n="1">¹ European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile</footnote>

<footnote n="2">² Space sciences, Technologies & Astrophysics Research (STAR) Institute, University of Liège, Liège, Belgium
e-mail: youssef.moulane@doct.uliege.be</footnote>

<footnote n="3">³ Oukaimeden Observatory, High Energy Physics and Astrophysics Laboratory, Cadi Ayyad University, Marrakech, Morocco</footnote>

<footnote n="4">⁴ Institut UTINAM UMR 6213, CNRS, Univ. Bourgogne Franche-Comté, OSU THETA, BP 1615, 25010 Besançon Cedex, France</footnote>

<footnote n="5">⁵ Koyama Astronomical Observatory, Kyoto Sangyo University, Motoyama, Kamigamo, Kita-ku, Kyoto 603-8555, Japan</footnote>

<footnote n="6">⁶ Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK</footnote>

<footnote>* Based on observations collected at the European Southern Observatory under ESO program 2101.C-5051.</footnote>

<footnote n="1">¹ https://ssd.jpl.nasa.gov/sbdb.cgi?sstr=21P;old=0; orb=0;cov=0;log=0;cad=0#discovery</footnote>
production was observed about one month before perihelion for the previous apparitions (Schleicher et al. 1987; Hanner et al. 1992; Lara et al. 2003). Both its unusual composition and the behaviour of its activity during multiple apparitions make 21P an object of great interest. In addition, as it is the parent body of the Draconids, a study of its dust properties might be valuable. The 2018 apparition was very favourable for ground-based observations because the comet was close to both the Sun and Earth at the same time and reached high elevation. Many observed the comet again for the 2018 return using various advanced instrumentation (IR and optical spectrographs on large telescopes) in order to better understand these peculiarities.

This work is organised as follows: after the introduction and historical background given in Sect. 1, we describe the observing circumstances and the reduction process of images obtained with the Transiting Planets and Planetesimals Small Telescope (TRAPPIST) and spectra obtained with UVES at the Very Large Telescope (VLT) in Sect. 2. In Sect. 3 we compute the production rates and discuss the gas and dust activity pre- and post-perihelion as well as the properties of the dust. The relative molecular abundances and their evolution with respect to the heliocentric distance are discussed in Sect. 4 and compared to the IR abundances of the mother species. In Sect. 5 we present the nitrogen and carbon isotopic ratios and the NH₂ (and NH₃) ortho-para ratio derived from the UVES high-resolution spectrum. In Sect. 6 we investigate the dynamical evolution of comet 21P within the last 10⁵ yr. Discussion of the chemical composition of 21P and the possible scenarios of its depletion in carbon species are given in Sect. 7. The summary and conclusions are given in Sect. 8.

## 2. Observation and data reduction

### 2.1. Photometry (TRAPPIST)

We started monitoring 21P with TRAPPIST-North (hereafter TN) at the beginning of June 2018 when the comet was at 1.55 au from the Sun. The comet was then observed from the Southern Hemisphere with TRAPPIST-South (hereafter TS) from the beginning of September 2018. The pair of TRAPPIST telescopes (Jehin et al. 2011) is very useful in this case as it allowed a continuous monitoring of the comet before and after perihelion. We collected images with the cometary HB narrow-band filters (Farnham et al. 2000) to measure the production rates of the radicals OH, NH, CN, C₃, and C₂. We also acquired images with the dust-continuum filters BC, GC, and RC for blue, green, and red (Farnham et al. 2000). We used the broad-band filters B, V, Rc, and Ic (Bessell 1990) to compute the Afp parameter, which is a proxy of the dust production rate (A'Hearn et al. 1984), and to derive the dust colours.

Throughout the passage of the comet, we made a high-cadence monitoring of 21P with images taken about twice a week. On photometric nights, we also obtained long series of observations with the gas narrow-band filters, especially CN and C₂ filters, to measure the variations in production rates during the same night that are caused by the rotation of the nucleus. We chose the exposure time of the different filters depending on the brightness of the comet. We used exposure times between 60 and 240 s for the broad-band filters and between 600 and 1500 s for the narrow-band filters. Observational circumstances and the number of sets of each filter are summarised in Table A.1. We started to collect data three months before perihelion to four months after perihelion. The comet reached perihelion on September 10, 2018, at a heliocentric distance of 1.01 au and a

<table><caption>**Table 1.** Scale lengths and fluorescence efficiency of different molecules at 1 au scaled by $r_h^{-2}$.</caption><thead><tr><th>Molecules</th><th>Parent (km)</th><th>Daughter (km)</th><th>g-factors (erg s<sup>-1</sup> mol<sup>-1</sup>)</th></tr></thead><tbody><tr><td>OH(0,0)</td><td>2.4 × 10<sup>4</sup></td><td>1.6 × 10<sup>5</sup></td><td>1.49 × 10<sup>-15</sup></td></tr><tr><td>NH(0,0)</td><td>5.0 × 10<sup>4</sup></td><td>1.5 × 10<sup>5</sup></td><td>6.27 × 10<sup>-14</sup></td></tr><tr><td>CN(Δv = 0)</td><td>1.3 × 10<sup>4</sup></td><td>2.1 × 10<sup>5</sup></td><td>2.62 × 10<sup>-13</sup></td></tr><tr><td>C<sub>3</sub>(λ = 4050 Å)</td><td>2.8 × 10<sup>3</sup></td><td>2.7 × 10<sup>5</sup></td><td>1.00 × 10<sup>-12</sup></td></tr><tr><td>C<sub>2</sub>(Δv = 0)</td><td>2.2 × 10<sup>4</sup></td><td>6.6 × 10<sup>4</sup></td><td>4.50 × 10<sup>-13</sup></td></tr></tbody></table>

**Notes.** The scale lengths are equivalent to the lifetimes of molecules as we used a constant radial velocity of 1 km s⁻¹ (A'Hearn et al. 1995). The fluorescence efficiency is taken from Schleicher's website (https://asteroid.lowell.edu/comet/gfactor.html).

geocentric distance of 0.39 au. In total, the comet was observed on 50 different nights with both telescopes, 13 nights before perihelion and 37 nights after. We used the same procedures as described in our previous papers (e.g. Opitom et al. 2015a; Moulane et al. 2018 and references therein) to reduce the data and to perform the flux calibration. To compute the production rates, we converted the flux of different gas species into column densities and adjusted their profiles with a Haser model (Haser 1957). This simple model, which is widely used, is based on a number of assumptions. Outgassing is assumed to be isotropic, and the gas has a constant radial velocity of 1 km s⁻¹. Parent molecules emitted by the nucleus decay by photodissociation to produce the observed daughter molecules. The model adjustment was performed at a physical distance of 10 000 km from the nucleus. Table 1 shows the scale lengths and g-factors of different molecules at 1 au scaled by $r_h^{-2}$. More details about the Haser model and its parameters are given in our previous works (see Moulane et al. 2018 and references therein). We would like to point out that we used the same parameters as Schleicher & Knight (2018) for the previous apparitions of 21P. We derived the Afp parameter, a proxy for the dust production (A'Hearn et al. 1984), from the dust profiles in the cometary dust continuum BC, GC, and RC filters and the broad-band Rc and Ic filters. It was computed at 10 000 km from the nucleus and corrected for the phase-angle effect according to the phase function normalised at $\theta=0^\circ$ derived by D. Schleicher<fnref n="2" />.

### 2.2. Spectroscopy (UVES/VLT)

We obtained one spectrum of comet 21P with the Ultraviolet-Visual Echelle Spectrograph (UVES) mounted on the Unit 2 telescope (UT2) at ESO's VLT on September 18, 2018 (a week after perihelion, $r_h = 1.01$ au and $\Delta = 0.40$ au) under Director's Discretionary Time. We used the UVES standard settings DIC#1 346+580 covering the range 3030 to 3880 Å in the blue and 4760 to 6840 Å in the red. We used a 0.44" wide slit, providing a resolving power $R \sim 80 000$. We obtained one single exposure of 3000 s at 8h35 UT with a mean air-mass of 1.7. This exposure provided two different spectra, which covered one of the above mentioned spectral ranges each.

The ESO UVES pipeline was used to reduce the spectra in the extended-object mode, keeping the spatial information. The spectra were corrected for the extinction and flux calibrated using the UVES master response curve provided by ESO.

<footnote n="2">² http://asteroid.lowell.edu/comet/dustphase.html</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Logarithm of the production rates (in mol s⁻¹) of each observed species and of the A(0)fρ parameter (in cm) of comet 21P during its 2018 return (this work and Schleicher & Knight 2018) compared with two previous apparitions in 1985 and 1998 (Schleicher & Knight 2018) as a function of heliocentric distance. The dashed vertical line represents the perihelion distance on September 10, 2018. The maximum of the gas and dust activity was reached at 1.07 au from the Sun on August 17, 2018, 24 days before perihelion.</figcaption></figure>

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 Earth. More details about the UVES data reduction are given in the UVES manual<fnref n="3" />. The dust-reflected sunlight was finally removed using a reference solar spectrum BASS2000⁴. A more detailed description of the steps for UVES data reduction and the solar spectrum subtraction is given in Manfroid et al. (2009) and references therein.

**3. Activity and composition**

**3.1. Gas-production rates**

Along with our monitoring of comet 21P with the TRAPPIST telescopes, we derived the OH, NH, CN, C₃, and C₂ production rates. They are summarised in Table A.2, and their evolution as a function of the heliocentric distance is compared with two previous passages in Fig. 1.

We started to detect most of the radicals in the coma by the end of June 2018 (except for NH, which was detected one month later). The various production rates and the dust activity

slowly increased as the comet approached the Sun (from 1.52 to 1.07 au). The highest activity was reached at 1.07 au from the Sun, on August 17, 24 days before perihelion. It then started to decrease rapidly after perihelion. CN was detected in our data until the end of 2018 at 1.66 au, while OH, C₃, and C₂ were no longer detected after the beginning of November at 1.4 au and NH in early October at 1.2 au. We found that like for the previous apparitions (Schleicher et al. 1987), the production rates pre-perihelion are higher by more than a factor two than post-perihelion. It is clear from Fig. 1 that the asymmetric activity is seen for all species, and this behaviour does not change over the various apparitions (Schleicher et al. 1987; Combi et al. 2011), as shown also in Fig. 2 for the water-production rate. The same behaviour has been reported for the parent molecule (H₂O, CO, CH₄, C₂H₂, C₂H₆, NH₃, and CH₃OH) production rates derived at IR wavelengths during the 2018 passage (Faggi et al. 2019; Roth et al. 2020). This might be due to the shape of the nucleus and its spin-axis orientation. This effect has been observed in several comets, such as 9P/Tempel 1 (Schleicher 2007), 81P/Wild 2 (Farnham & Schleicher 2005), and also for comet 67P/Churyumov-Gerasimenko (Schleicher 2006; Opitom et al. 2017). It has been shown very clearly by the Rosetta mission that the highest activity of 67P was clearly associated with the illumination of the most southern regions, which received the highest solar flux, and were subject to intense erosion

<footnote n="3">³ ftp://ftp.eso.org/pub/dfs/pipelines/uves/uves-pipeline-manual-22.17.pdf</footnote>

<footnote n="4">⁴ http://bass2000.obspm.fr/solarspect.php</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** H₂O production rates of comet 21P as a function of days to perihelion in 2018 compared to the previous apparitions (1985, 1998, and 2005). The data are described in more detail in Sect. 3.2.</figcaption></figure>

(Lai et al. 2019). Recently, Marshall et al. (2019) have shown for those three comets that the nucleus shape, the spin-axis orientation, and the distribution of activity on the cometary surface can explain the light curve of the water-production rate as a function of the heliocentric distance.

Around the activity maximum, the production rates are almost the same as those measured in the previous apparitions, showing that the activity level of 21P did not decrease over the last five orbits. Our production rates usually agree very well with those derived by Schleicher & Knight (2018), who used the same technique, while we note a discrepancy at large heliocentric distance and post-perihelion in the 1985 and 1998 apparition data (Schleicher et al. 1987; Schleicher & Knight 2018). This could be due to a sensitivity problem in their data as the production rates seem to level off on either side of perihelion, while the distance increased to a higher activity level of the comet after perihelion in the past. It has been also found that there is no significant change in the production rates of hypervolatile molecules (CO, CH₄, and C₂H₆) in comet 21P over the three different apparitions, 1998 (Weaver et al. 1999; Mumma et al. 2000), 2005 (DiSanti et al. 2012), and 2018 (Faggi et al. 2019; Roth et al. 2020).

## 3.2. H₂O production rate

The water-production rate is the most significant indicator of the activity of a comet. It can be measured directly from near-infrared observations or derived from OH emission at 3090 Å and radio wavelengths or from H Lyman-α emission at 1216 Å (Combi et al. 1986) assuming that both OH and H arise from the dissociation of H₂O. In this work, we computed the vector-equivalent water-production rates according to an empirical procedure based on a comparison of OH and water-production rates derived from the mean lifetimes, velocities, and scale lengths given by Cochran & Schleicher (1993). Schleicher et al. (1998) built an empirical relationship $Q(H_2O) = 1.361r_h^{-0.5}Q(OH)$ based on an $r_h^{-0.5}$ dependence of the H₂O outflow velocity, a photodissociation branching ratio for water to OH of 90%, and the heliocentric distance. Figure 2 shows the

water production we derived compared to previous apparitions (with different techniques) as a function of days to perihelion. We used the formula given above to convert $Q(OH)$ into water production rates for the Schleicher & Knight (2018) data. Combi et al. (2011) derived the production rates from the H Ly-α emission observed by the Solar Wind ANisotropy (SWAN) instrument on board the Solar and Heliospheric Observatory (SOHO) in 1998 and in 2005. The Combi & Feldman (1992) values are derived from H Ly-α emission observed by the IUE mission for the 1985 apparition. From the Pioneer Venus Orbiter ultraviolet system (UVS) instrument, McFadden et al. (1987) derived the water-production rates from OH (3090 Å) emission. Fink & Hicks (1996) derived the water-production rates from the [OI]($^1$D) forbidden-line doublet using the correlation between the water-production rates and the total photon luminosity. Faggi et al. (2019) and Roth et al. (2020) measured the water-production rates directly from near-infrared spectra. TRAPPIST and UVES data points are from this work (see Sects. 3.2 and 5.1). The maximum in the last four apparitions was reached about one month before perihelion and does not change over all apparitions, but we observe a clear systematic difference between the narrow-band and spectroscopic methods in the optical on one hand and the measurements made from the space observations of the H Ly-α emission in the UV on the other hand. The maximum of the water production we measured was on August 17, 24 days before perihelion, and it reached (3.72 ± 0.07) × 10²⁸ mol s⁻¹, in good agreement with the Schleicher & Knight (2018) measurement of 4.20 × 10²⁸ mol s⁻¹ at the heliocentric distance of 1.07 au. Using the same technique for the 1985 apparition, Schleicher et al. (1987) reported $Q(H_2O) = 4.85 × 10²⁸ mol s⁻¹$ when the comet was at 1.05 au from the Sun. Using high-resolution infrared spectroscopy, Weaver et al. (1999) measured ~2–3 × 10²⁸ mol s⁻¹ at $r_h = 1.10$ au in 1998. For the 2005 apparition, Combi et al. (2011) found a value of 5.80 × 10²⁸ mol s⁻¹ from the H Ly-α emission observed by the SWAN/SOHO at $r_h = 1.08$ au. Comparing these data, we found that the water-production rates measured by Combi & Feldman (1992) in 1985 and Combi et al. (2011) in 2005 are systematically higher by a factor of about two than our results in
<table><caption>**Table 2.** Active area (km²) and active fraction of the surface (%) for 21P using the slow-rotator model at some interesting heliocentric distances.</caption><thead><tr><th>Date UT</th><th>r<sub>h</sub><br>(au)</th><th>Active area<br>(km²)</th><th>Active fraction<br>(%)</th></tr></thead><tbody><tr><td>(a) 2018 Jun. 22</td><td>-1.49</td><td>4.9 ± 0.1</td><td>17.5 ± 11.6</td></tr><tr><td>(b) 2018 Aug. 17</td><td>-1.07</td><td>12.0 ± 0.2</td><td>42.5 ± 28.0</td></tr><tr><td>(c) 2018 Sep. 09</td><td>-1.01</td><td>9.0 ± 0.2</td><td>32.0 ± 21.3</td></tr><tr><td>(d) 2018 Sep. 15,16,17</td><td>+1.02</td><td>6.9 ± 0.1</td><td>21.5 ± 14.0</td></tr><tr><td>(e) 2018 Nov. 09</td><td>+1.31</td><td>0.90 ± 0.04</td><td>3.3 ± 2.2</td></tr></tbody></table>

**Notes.** (a) The first measurement during our monitoring campaign, (b) the maximum activity during our monitoring campaign, (c) the last measurement before perihelion passage, (d) the mean of the first measurements after perihelion passage with similar heliocentric distances of 1.02 au, and (e) the last and minimum measurements during our monitoring campaign. Large errors in the active fractions of the surface come from the large uncertainties in the radius of the nucleus, which we adopted as 1.5 ± 0.5 km (see Sect. 3.3 for details).

2018. This offset between various techniques has been reported in previous studies and as early as in studies for comet 1P/Halley (Schleicher et al. 1998). The origin of this discrepancy is not clear, but it is obvious that the agreement is good when the same technique is used. This indicates that the level of activity of 21P was the same over the past four decades and did not decrease, like that of comet 41P/Tuttle-Giacobini-Kresak, which lost as much as 30 to 40% of its activity from one orbit to the next (Moulane et al. 2018).

### 3.3. Active area of the nucleus

To estimate the active area of the nucleus surface, we modelled the water production using the sublimation model of Cowan & A'Hearn (1979). Because of the low thermal inertia of cometary nuclei (Gulkis et al. 2015), the slow-rotator approach was adopted in a number of cases as the most appropriate way for computing the cometary out-gassing (see e.g. Bodewits et al. 2014; Lis et al. 2019). The slow-rotator model assumes that every facet of the nucleus is in equilibrium with the solar radiation incident upon it, with the rotational pole pointed at the Sun. As mentioned previously, the size of the 21P nucleus, which is necessary to convert the active area into the active fraction of the whole surface, is not well constrained so far, with a radius ranging from 1 to 2 km. To estimate the active fraction of the surface, we therefore assumed a radius of 1.5 ± 0.5 km. Moreover, we assumed a bond albedo of 5 and a 100% infrared emissivity (see e.g. A'Hearn et al. 1989; McKay et al. 2018, 2019). We found that the active area of 21P during our monitoring campaign varied from ~5 km² at 1.49 au pre-perihelion, reached a maximum of ~12 km² at 1.07 au pre-perihelion, and decreased to ~1 km² at 1.31 au post-perihelion. Table 2 shows the lowest and highest active areas and active fraction for 21P using the slow-rotator model at some interesting heliocentric distances. We obtained different values than previous estimates given by Combi et al. (2019). The reason is twofold. First, the already mentioned discrepancy in the water-production rates found through different observational techniques (see Sect. 3.2), and second, the model used by the authors (fast rotator), which is less appropriate to describe the cometary outgassing.

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** A(0) *f*ρ parameter measurements, computed at 10 000 km from the nucleus and corrected for the phase-angle effect, for the broad band (Rc and Ic) and narrow-band cometary filters (RC, GC, and BC) as a function of days to perihelion.</figcaption></figure>

### 3.4. Dust properties

We computed the A(0) *f*ρ parameter at 10 000 km, as defined by A'Hearn et al. (1984), using broad-band (Rc and Ic) and narrow-band dust continuum filters (RC, GC, BC) (see Table A.2). Figure 3 shows its evolution as a function of time to perihelion. Our results agree very well with those reported by Ehlert et al. (2019), as shown in the bottom right panel of Fig. 1. Like for the gas, the maximum was reached on August 17 with a value of (1646.1 ± 12.8) cm in the red narrow-band RC filter. About the same value was reported on previous apparitions at the same heliocentric distance (Schleicher et al. 1987; Lara et al. 2003; Pittichová et al. 2008). For a detailed description of the dust environment and its evolution, a more sophisticated model should be used, such as the Monte Carlo model presented by Moreno et al. (2012), which was used successfully in a number of cases (see e.g. Pozuelos et al. 2015, 2018; Moreno et al. 2016a,b, 2017). This study will be presented in a separate paper.

We used the A(0) *f*ρ values obtained with the narrow-band filters, which are not contaminated by the gas emission, to derive the dust colours. The normalised reflectivity gradients between wavelength λ₁ and λ₂ are defined as (A'Hearn et al. 1984; Jewitt & Meech 1986)

$$S_v(\%/1000 \AA) = \frac{Af\rho_1 - Af\rho_2}{Af\rho_1 + Af\rho_2} \times \frac{2000}{\lambda_1 - \lambda_2} \quad (1)$$

λ₁ and λ₂ are the effective wavelengths of the filters: BC[4450 Å], GC[5260 Å], and RC[7128 Å].

We found that the RC-GC, RC-BC, and GC-BC colours are redder than those of the Sun, with mean values of (14.8 ± 3.3), (13.2 ± 2.6), and (12.4 ± 7.5)%/1000 Å, respectively (see Fig. 4). This result agrees with previous apparitions, with values of S<sub>v</sub> = 15%/1000 Å in 1985 (Schleicher et al. 1987) and S<sub>v</sub> = 13%/1000 Å in 1998 (Lara et al. 2003). These values are consistent with the colour of the nucleus of 21P (12.8 ± 2.7)%/1000 Å measured at a heliocentric distance of 3.5 au in 1991 (Luu 1993). They fall within the range observed for most JFCs (Lamy & Toth 2009; Solontoi et al. 2012;
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Normalised reflectivity gradients $S_v$ (% per 1000 Å) of comet 21P for different colour indices as a function of days to perihelion.</figcaption></figure>

Jewitt 2015). During our long monitoring, we did not detect any significant variation in the colour of the dust in the coma (or any outburst).

## 4. Abundance ratios

The molecular abundances and their ratios with respect to the distance to the Sun give information about the homogeneity of a cometary nucleus and the chemical processes involved in the coma. Based on the relative abundance of 41 comets, A'Hearn et al. (1995) classified comets into two groups based on their C₂-to-CN ratio. Typical comets are defined as those having a $log[Q(C₂)/Q(CN)] ≥ -0.18$, while the carbon-chain-depleted comets are those below that value. This classification was confirmed later by other photometric and spectroscopic studies of large data sets (Schleicher 2008; Fink 2009; Langland-Shula & Smith 2011; Cochran et al. 2012) and must reflect some differences between the formation conditions (the pristine scenario) or a change of relative composition with time (several perihelion passages) of these comets (the evolutionary scenario). Figure 5 shows the evolution of the 21P abundance ratios of the various radicals with respect to OH (a proxy of water) and CN as a function of heliocentric distance. It is clear that the 21P abundance ratios in the 2018 return agree with the mean values of depleted comets given in A'Hearn et al. (1995). Table 3 summarizes the relative abundances in 2018 compared to 1985 and 1998 data using the same technique and the same Haser model parameters (Schleicher & Knight 2018). Our 2018 ratios are the mean values for all the data obtained (see Table A.2). Like for the activity level over the past passages, the relative abundances did not change over the last five orbits. We note that the *Afρ* values derived in 1985 and in 1998 by Schleicher & Knight (2018) were computed for the narrow-band GC[5260 Å] filter, while we used the BC[4450 Å] filter. After correcting their *Afρ* values for the phase-angle effect using the same function as for the TRAP-PIST data (see Sect. 2.1), the two data sets agree. This indicates that there is no evidence of changes in the chemical composition in the coma of the comet at different heliocentric distances (in the range 1.0 to 1.5 au) and over the five orbits, which is an argument for rejecting the evolutionary origin of the carbon-chain depletion in this comet.

The bottom panel of Fig. 5 shows that there is no evidence either that the dust-to-gas ratio represented by $A(0)\rho/Q(CN)$

<table><caption>**Table 3.** Relative molecular abundances of comet 21P over the last passages compared to the mean values for typical comets.</caption><thead><tr><th rowspan="2"></th><th colspan="4">Log production rate ratio</th></tr><tr><th>1985 <sup>(a)</sup></th><th>1998 <sup>(a)</sup></th><th>2018 <sup>(b)</sup></th><th>Typical comets <sup>(c)</sup></th></tr></thead><tbody><tr><td>C<sub>2</sub>/CN</td><td>-0.64</td><td>-0.50</td><td>-0.52 ± 0.10</td><td>0.06 ± 0.10</td></tr><tr><td>C<sub>3</sub>/CN</td><td>-1.42</td><td>-1.30</td><td>-1.39 ± 0.12</td><td>-1.09 ± 0.11</td></tr><tr><td>CN/OH</td><td>-2.59</td><td>-2.67</td><td>-2.62 ± 0.08</td><td>-2.50 ± 0.18</td></tr><tr><td>C<sub>2</sub>/OH</td><td>-3.23</td><td>-3.17</td><td>-3.16 ± 0.21</td><td>-2.44 ± 0.20</td></tr><tr><td>C<sub>3</sub>/OH</td><td>-4.02</td><td>-3.98</td><td>-4.03 ± 0.16</td><td>-3.59 ± 0.29</td></tr><tr><td>NH/OH</td><td>-2.66</td><td>-2.87</td><td>-2.68 ± 0.14</td><td>-2.37 ± 0.27</td></tr><tr><td>A(0)ρ/CN</td><td>-22.74</td><td>-22.73</td><td>-22.70 ± 0.04</td><td>-23.30 ± 0.32</td></tr><tr><td>A(0)ρ/OH</td><td>-25.33</td><td>-25.42</td><td>-25.32 ± 0.04</td><td>-25.82 ± 0.40</td></tr></tbody></table>

**References.** <sup>(a)</sup>Schleicher & Knight (2018), <sup>(b)</sup>this work, <sup>(c)</sup>A'Hearn et al. (1995).

and $A(0)\rho/Q(OH)$ depends on the heliocentric distance. We found that this ratio in 21P is consistent with the average value of depleted comets and higher than the mean value of the typical comets as defined in A'Hearn et al. (1995) (see Table 3). Lara et al. (2003) obtained a value of $log[A(0)\rho/Q(CN)] = -22.91 \pm 0.10$ in 1998 that agrees with our measurement. Like for the gas relative abundances, we conclude that the coma of 21P does not show significant variation in the dust-to-gas ratio over the previous apparitions and that it has a similar ratio as the depleted comets defined by A'Hearn et al. (1995).

The comparison with abundances of parent molecules derived from IR data (Faggi et al. 2019; Roth et al. 2020) allows us to investigate the origin of the radicals observed in the atmosphere of 21P. Table 4 shows the comparison between the production rates of daughter molecules observed in the optical and the possible parent molecules observed at IR. Using high-resolution infrared spectra obtained in 1998, Weaver et al. (1999) reported upper limits for different species relative to H₂O such as C₂H₆ (2–3%), HCN(0.3–0.4%), and C₂H₂(0.5–0.8%) assuming that all species are parent molecules. C₂H₂ has been found to be depleted with respect to HCN by a factor five compared to other comets such as Hyakutake and Hale-Bopp. This result has been confirmed at this apparition by Faggi et al. (2019).

We derived a $Q(C₂) = 1.26 \times 10^{25} \text{ mol s}^{-1}$, which is consistent with the upper limit of $Q(C₂H₂) < 4.52 \times 10^{25} \text{ mol s}^{-1}$ reported by Faggi et al. (2019) and $<1.80 \times 10^{25} \text{ mol s}^{-1}$ reported by Roth et al. (2018) at 1.18 au from the Sun. This agreement indicates that C₂ might be a daughter species of C₂H₂. C₂ may also come from C₂H₆ and HC₂N (Helbert et al. 2005; Weiler 2012; Hölscher 2015) or be released from organic-rich grains (Combi & Fink 1997), but a detailed chemical model of the coma is required for more details. We also found a good match between $Q(CN) = 4.40 \times 10^{25}$ and $Q(HCN) = 4.30 \times 10^{25} \text{ mol s}^{-1}$ (Faggi et al. 2019) at 1.01 au, showing that HCN might be the main parent molecule of CN in 21P. This result is known for several comets using different methods, including a comparison between the HCN and CN production rates (Rauer et al. 2003; Opitom et al. 2015b), coma morphology (Woodney et al. 2002), and also carbon and nitrogen isotopic ratios in both species (Manfroid et al. 2009; Bockelée-Morvan et al. 2015). We note, however, that in some cases, the two abundances do not agree, and other sources, for instance extended sources, have been claimed for the CN origin (Fray et al. 2005).

Some molecules such as C₄H₂, CH₂C₂H₂, and CH₃C₂H are proposed to be the parent molecules of C₃ (Helbert et al. 2005; Mumma & Charnley 2011; Hölscher 2015), but these complex
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Evolution of the logarithmic production rates ratios of each species with respect to OH and to CN as a function of heliocentric distance. The dashed red line represents the mean value of typical comets as defined in A'Hearn et al. (1995), while the blue line represents the mean value of the depleted group. *Bottom panels:* dust-to-gas ratio represented by $A(0)f\rho$-to-OH and $A(0)f\rho$-to-CN. The vertical dashed line shows the perihelion distance on September 10, 2018, at $r_h = 1.01$ au.</figcaption></figure>

<table><caption>**Table 4.** Comparison of daughter molecules and possible parent molecule production rates derived from optical and infrared data of comet 21P in the 2018 passage.</caption><thead><tr><th rowspan="2">UT date</th><th rowspan="2">r<sub>h</sub><br>(au)</th><th rowspan="2">Δ<br>(au)</th><th colspan="8">Production rates (10<sup>25</sup> mol s<sup>-1</sup>)</th><th rowspan="2">Reference</th></tr><tr><th>Q(OH)</th><th>Q(H<sub>2</sub>O)</th><th>Q(CN)</th><th>Q(HCN)</th><th>Q(C<sub>2</sub>)</th><th>Q(C<sub>2</sub>H<sub>2</sub>)</th><th>Q(C<sub>2</sub>H<sub>6</sub>)</th><th>Q(NH)</th><th>Q(NH<sub>3</sub>)</th></tr></thead><tbody><tr><td>2018 Jul 30</td><td>1.17</td><td>0.61</td><td>1810 ± 28</td><td></td><td>4.55 ± 0.07</td><td>—</td><td>1.26 ± 0.07</td><td>—</td><td>—</td><td>3.13 ± 0.50</td><td>—</td><td>This work</td></tr><tr><td>2018 Jul 30</td><td>1.17</td><td>0.61</td><td>—</td><td>2401 ± 394</td><td>—</td><td>&lt;3.20</td><td>—</td><td>&lt;4.52</td><td>4.49 ± 1.45</td><td>—</td><td>&lt;63.72</td><td>Faggi et al. (2019)</td></tr><tr><td>2018 Jul 31</td><td>1.16</td><td>0.59</td><td>—</td><td>2503 ± 385</td><td>—</td><td>6.16 ± 0.12<sup>(a)</sup></td><td>—</td><td>&lt;1.80<sup>(a)</sup></td><td>6.05 ± 0.77</td><td>—</td><td>&lt;16.19<sup>(a)</sup></td><td>Roth et al. (2020)</td></tr><tr><td>2018 Sep 07</td><td>1.01</td><td>0.39</td><td>3036 ± 357</td><td>3206 ± 112</td><td>—</td><td>—</td><td>—</td><td>—</td><td>10.60 ± 1.10</td><td>—</td><td>—</td><td>Roth et al. (2020)</td></tr><tr><td>2018 Sep 09</td><td>1.01</td><td>0.39</td><td>—</td><td>2623 ± 586</td><td>—</td><td>4.30 ± 0.32</td><td>—</td><td>&lt;0.62</td><td>8.30 ± 1.38</td><td>—</td><td>&lt;12.59</td><td>Faggi et al. (2019)</td></tr><tr><td>2018 Sep 09</td><td>1.01</td><td>0.39</td><td>2360 ± 33</td><td>—</td><td>4.39 ± 0.07</td><td>—</td><td>1.67 ± 0.06</td><td>—</td><td>—</td><td>5.66 ± 0.38</td><td>—</td><td>This work</td></tr><tr><td>2018 Oct 07</td><td>1.07</td><td>0.49</td><td>—</td><td>2583 ± 864</td><td>—</td><td>—</td><td>—</td><td>—</td><td>4.55 ± 1.44</td><td>—</td><td>—</td><td>Faggi et al. (2019)</td></tr><tr><td>2018 Oct 08</td><td>1.08</td><td>0.49</td><td>834 ± 25</td><td>—</td><td>1.68 ± 0.05</td><td>—</td><td>0.40 ± 0.06</td><td>—</td><td>—</td><td>1.07 ± 0.35</td><td>—</td><td>This work</td></tr><tr><td>2018 Oct 10</td><td>1.10</td><td>0.51</td><td>—</td><td>2028 ± 255</td><td>—</td><td>—</td><td>—</td><td>—</td><td>2.92 ± 0.39</td><td>—</td><td>—</td><td>Roth et al. (2020)</td></tr></tbody></table>

Notes. <sup>(a)</sup>From Roth et al. (2018) measured on July 29, 2018. Upper limits are $3\sigma$ for the results of Roth et al. (2020) and Faggi et al. (2019).

species were not observed at infrared or at radio wavelengths. NH and NH<sub>2</sub> were found to be depleted in 21P in the previous apparitions (A'Hearn et al. 1995; Fink 2009). New infrared observations in 2018 show very low NH<sub>3</sub> in 21P, with an upper limit ratio of $Q(NH_3)/Q(H_2O) < 0.6\%$ (Faggi et al. 2019). In this work, we derive $Q(NH)/Q(OH) = 0.2\%$, which is consistent with $Q(NH_3)/Q(H_2O)$.

# 5. Optical high-resolution spectrum

## 5.1. Water-production rate

The UVES spectrum offered the possibility of independently computing the water-production rate at the time of observation. We first measured the overall flux for the OH (0,0) band near 309 nm, integrated over the whole slit. We found
<figure><img src="image_7.png" /><figcaption>Fig. 6. Observed and synthetic CN spectra of the *R* branch of the B-X (0, 0) violet band in comet 21P.</figcaption></figure>

1.47 × 10⁻¹³ erg s⁻¹ cm⁻² arcsec⁻². The fluorescence efficiency computed for this band and the heliocentric distance and velocity at the time of observation was 2.62 × 10⁻⁴ s⁻¹ (or 1.71 × 10⁻¹⁵ erg s⁻¹ mol⁻¹ when scaled to 1 au; see details on the fluorescence model in Rousselot et al. 2019). From these values and a Monte Carlo simulation of the water molecules creating OH radicals in the inner coma (model based on equations given by Combi & Delsemme 1980) it is possible to compute the corresponding water-production rate for the number of OH radicals observed in the slit (0.44 × 9.5 arcsec) centred on the nucleus. Using the parameters of H₂O radial velocity, OH and H₂O lifetimes given in Cochran & Schleicher (1993) and assuming that 91.8% of water molecules dissociate to OH (Crovisier 1989), we found Q(H₂O) = 1.7 × 10²⁸ mol s⁻¹. This result is in excellent agreement with the water-production rates computed from TRAPPIST observations in the same period (see Fig. 2). It must, nevertheless, be pointed out that it depends of the different parameters and can change slightly with them, especially with the water lifetime.

## 5.2. ¹²C/¹³C and ¹⁴N/¹⁵N isotopic ratios

The study of the isotopic ratios in comets has attracted considerable attention as it contains information about the conditions that prevailed at the time of formation of these objects in the early Solar System (Jehin et al. 2009; Hyodo et al. 2013). The carbon ¹²C/¹³C ratio has been determined for several comets from the analysis of the C₂ Swan band and CN B-X system in the optical (Manfroid et al. 2009; Bockelée-Morvan et al. 2015, and references therein). Some in situ measurements have also been obtained in comet 67P by the ROSINA mass spectrometer on board the Rosetta spacecraft for C₂H₄, C₂H₅, CO (Rubin et al. 2017), and CO₂ molecules (Hässig et al. 2017). All derived values are compatible with the terrestrial ratio of 89, except for CO, which might be slightly enriched in ¹³C. The

nitrogen ¹⁴N/¹⁵N isotopic ratio was measured for the first time from high-resolution spectra of the CN violet band in comets C/2000 WM1 (LINEAR) and C/1995 O1 (Hale-Bopp) and was found to be enriched by a factor of two in ¹⁵N with respect to the Earth value (Arpigny 2003). The same ratio was found later from sub-millimeter observations of HCN in comet 17P/Holmes during its outburst and archival data of C/1995 O1 (Hale-Bopp) (Bockelée-Morvan et al. 2008). It has also been possible recently to measure the ¹⁴N/¹⁵N ratio in ammonia via the NH₂ radical (Rousselot et al. 2014). The values obtained are similar to the value found in HCN and CN, which was confirmed by subsequent works (Shinnaka et al. 2014, 2016; Rousselot et al. 2015; Shinnaka & Kawakita 2016; Yang et al. 2018). Recent measurements performed by the ROSINA mass spectrometer in comet 67P provided a ratio ¹⁴N/¹⁵N = 118 ± 25 for NH₃ and 130 ± 30 for N₂ molecules (Altwegg et al. 2019).

We used the ¹²C¹⁴N B-X (0,0) band to estimate the ¹²C/¹³C and ¹⁴N/¹⁵N isotopic ratios of 21P. We used a CN fluorescence model to create synthetic spectra of ¹³C¹⁴N, ¹²C¹⁵N, and ¹²C¹⁴N. More details of the model are given in Manfroid et al. (2009). Figure 6 shows the observed CN spectrum compared to the synthetic spectrum made under the same observational conditions. The ratios found for ¹²C/¹³C and ¹⁴N/¹⁵N are 100 ± 10 and 145 ± 10, 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 ¹²C/¹³C and ¹⁴N/¹⁵N, respectively (Manfroid et al. 2009; Bockelée-Morvan et al. 2015).

## 5.3. NH₂ and NH₃ ortho-to-para ratios

We measured the ortho-to-para abundance ratio (OPR) of NH₂ from the three rovibronic emissions bands (0,7,0), (0,8,0), and (0,9,0), see Fig. 7, following the method described in Shinnaka et al. (2011). The derived OPRs of NH₂ and of its parent molecule NH₃ are listed for each band in Table 5, and they have
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Comparison between the observed and modelled spectra of the NH₂(0,9,0), (0,8,0), and (0,7,0) bands. The modelled spectrum of C₂ is also plotted in the NH₂ (0,9,0) band panel, but because of the depleted nature of 21P, the C₂ lines do not affect the NH₂ spectrum. The ortho- and para-lines of NH₂ are labelled in these modelled spectra. The two strong emission lines at 6300 and 6364 Å are the forbidden oxygen lines in the NH₂ (0,8,0) band panel. We note that the intensity ratio among bands is not correct because we scaled intensity for each plot independently.</figcaption></figure>

average values of 3.38 ± 0.06 and 1.19 ± 0.03, respectively. The latter is in very good agreement with the Subaru/HDS determination (NH₃ OPR = 1.16 ± 0.02; Shinnaka et al. 2020). A nuclear spin temperature ($T_{\text{spin}}$) for ammonia of 27 ± 1 K was derived. The 21P value is consistent with typical values measured in comets (see Fig. 8). 21P therefore cannot be distinguished from other comets based on its NH₂ OPR (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 demonstrate that the OPR of water does not retain the memory of its formation temperature (Hama et al. 2011, 2016; Hama & Watanabe 2013). It is likely that this is also the case for ammonia.

<table><caption>**Table 5.** Derived NH₂ and NH₃ OPRs in comet 21P.</caption><thead><tr><th>NH₂ band</th><th>NH₂ OPR</th><th>NH₃ OPR</th><th>T<sub>spin</sub> (K)</th></tr></thead><tbody><tr><td>(0,7,0)</td><td>3.30 ± 0.13</td><td>1.15 ± 0.07</td><td>28<sup>+5</sup>/<sub>-3</sub></td></tr><tr><td>(0,8,0)</td><td>3.55 ± 0.08</td><td>1.28 ± 0.04</td><td>23<sup>+2</sup>/<sub>-1</sub></td></tr><tr><td>(0,9,0)</td><td>3.15 ± 0.10</td><td>1.08 ± 0.05</td><td>34<sup>+8</sup>/<sub>-4</sub></td></tr><tr><td>Average</td><td>3.38 ± 0.06</td><td>1.19 ± 0.03</td><td>27 ± 1</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Summary of NH₃ OPRs in comets. The orange and green crosses are the NH₃ OPR of 21P by VLT/UVES (this work) and by Subaru/HDS (Shinnaka et al. 2020), respectively.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Summary of the NH₃ $T_{\text{spin}}$ in 28 comets of various origins. The orange and green crosses are the NH₃ $T_{\text{spin}}$ of 21P by VLT/UVES (this work) and by Subaru/HDS (Shinnaka et al. 2020), respectively.</figcaption></figure>

The OPRs of cometary volatiles might have been modified by an ortho-to-para conversion process in the inner coma or other catalyst activities of dust crust surfaces of the nucleus rather than be reflected by a formation temperature in the solar nebula 4.6 Gyr ago. OPRs might be a diagnostic of the physico-chemical conditions in the innermost coma or beneath the surface.

## 6. Dynamical evolution

In this section we analyse the dynamical evolution of the comet within the last $10^5$ yr. JFCs are highly chaotic objects, whose
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Orbital evolution of 21P and its 200 clones for 3000 yr backward in time from January 1, 2020. *From top to the bottom:* closest approaches with Jupiter, semi-major axis, eccentricity, inclination, perihelion, and aphelion distance. In all panels, the grey lines correspond to the evolution of each clone, the black line is the mean value of the clones, and the red line is that of the nominal comet 21P. The vertical blue line corresponds to the time of the closest encounter with Jupiter. The initial orbital elements were taken from the JPL Small-Body Data Browser (orbital solution JPL K182/3).</figcaption></figure>

dynamic evolution must be studied in terms of statistics (Levison & Duncan 1994). With this in mind, we analysed the evolution of the original object, that is, comet 21P, by considering the nominal values of its orbital parameters as they are defined in JPL-HORIZONS (orbital solution JPL K182/3). In the analysis, 200 clones were generated following the covariance matrix of its orbital parameters<fnref n="5" />. We performed the integration with the numerical package MERCURY (Chambers 1999), using the integration algorithm of Bulirsch-Stoer (Bulirsch & Stoer 1964) with a time step of 8 d, and we included the Sun, all planets, and Pluto in the simulation. In addition, we also included non-gravitational forces. The results of the simulations are displayed in Fig. 10.

We find that the orbits of all the clones in the simulation were very compact for a period of ~1650 yr. After that period, the orbits started to scatter, which was provoked by a close encounter with Jupiter, at a mean distance of 0.1 au. Because the nature of

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Logarithm of C₂-to-CN ratio of 110 comets as a function of the Tisserand invariant parameter with respect to Jupiter (T<sub>j</sub>). Filled symbols present typical comets, while the open symbols present the carbon-chain-depleted comets. Our measurement of comet 21P is represented by a blue diamond. The vertical dashed line at T<sub>j</sub>=2 separates the families of JFCs and Oort cloud comets.</figcaption></figure>

JFCs is chaotic, a comparison of results from different authors who applied different methods (e.g. different integration algorithm, different number of clones, or different manner in which their clones were generated) is difficult to perform, and any superficial comparison might yield incorrect conclusions. Only one analysis identical to that performed here has been carried out for comet 66P/du Toit by Yang et al. (2019). The authors found that the comet belongs to the Jupiter family for at least ~60 × 10³ yr, and the stable nature of its orbit was evident. This result indicates that 21P likely is a young member of the Jupiter family that has crossed its perihelion ~230 times with similar distances of *q* ~ 1.013 au. The youth of 21P could explain its unusual composition. However, the lack of a statistic sample prevents us from robustly confirming this hypothesis.

## 7. Discussion

As mentioned above, C₂ and C₃ have been found to be depleted compared to CN in 21P more than 50 yr ago (Mianes et al. 1960; Herbig 1976; Schleicher et al. 1987). In the data set of A'Hearn et al. (1995), 21P was classified as the prototype of the group that is depleted in carbon-chain molecules. Figure 11 shows our C₂-to-CN ratio compared to 120 comets (90 comets from Schleicher 2008 and 30 comets from TRAPPIST database Opitom 2016) as a function of the Tisserand invariant parameter with respect to Jupiter (T<sub>j</sub>). About 30% of the comets analysed were found to be depleted in carbon-chain elements by varying amounts, including different dynamical types of comets, two-thirds are JFCs, and one-third are LPCs (A'Hearn et al. 1995; Schleicher 2008; Fink 2009; Cochran et al. 2012). 21P was found to be also depleted in NH with respect to OH (A'Hearn et al. 1995). This result was confirmed by its depletion in NH₂ using spectro-photometric observations by Konno & Wyckoff (1989) in the 1985 apparition and later by Fink & Hicks (1996) in the 1998 passage. This depletion in both NH and NH₂ indicates that 21P is likely depleted in the parent molecule NH₃, which was recently confirmed by Faggi et al. (2019). 21P is not a unique case of a comet depleted in both carbon-chain and ammonia daughter species. A few others have been found, but with a lesser degree of depletion, such

<footnote n="5">⁵ Both sets of the orbital parameters and the covariance matrix of the orbit for 21P are published together in the NASA/JPL small-body browser: https://ssd.jpl.nasa.gov/sbdb.cgi?sstr=21P;old=0;orb=0;cov=1;log=0;cad=0#elem</footnote>
as 43P/Wolf-Harrington and the split comet 73P/Schwassmann-Wachmann 3 (A'Hearn et al. 1995; Fink 2009; Cochran et al. 2012). This indicates that there might be a small group of similar comets that formed under similar conditions and is different from other comets. According to taxonomy studies, however, no clear grouping associated with an abundance in NH has been identified.

Our long monitoring of the abundance ratios combined with previous studies (Schleicher & Knight 2018; Combi et al. 2011) rules out the evolutionary scenario (peculiar composition due to repeated passages to perihelion). Our observations show remarkably constant abundance ratios of the different species, especially the depleted C₂ and NH, before and after perihelion, and over months. The fact that these ratios are still the same after five orbits is in favour of a pristine composition rather than compositional changes due to repeated passages of the comet at perihelion.

It was argued that this peculiar composition might be linked to a higher formation temperature, closer to that of the Sun (Schleicher et al. 1987), or in a local disc around Jupiter, as was proposed for comet 73P (Shinnaka et al. 2011). We obtained high-resolution optical spectra with high S/N in order to investigate the C and the N isotopic ratios, as well as the NH₂ OPR. 21P appears to have a normal ¹⁴N/¹⁵N ratio and a normal NH₃ OPR, similar to other comets. This is in contrast with comet 73P, which has both peculiar ¹⁴N/¹⁵N and OPR (see Fig. 9 and Shinnaka et al. 2011). The two comets therefore are apparently not related, and the peculiar composition of 21P still needs to be explained. These peculiarities are clearly linked to the ice composition of the nucleus because the IR studies of the mother molecules also show the same kind of depletion, with an obvious link to the daughter species.

## 8. Summary and conclusion

We performed an extensive monitoring of comet 21P on either side of perihelion with TRAPPIST. The gas species production rates as well as the dust proxy, A(0)fρ parameter, were computed until the detection limit. We derived the water-production rates for this apparition, and we compared it, as well as the various abundance ratios, to previous passages. Using a sublimation model for the nucleus and the water-production rates, we constrained the active area of the nucleus surface using the slow-rotator approach. An accurate determination of the 21P nucleus parameters is needed to better constrain the active area fraction. Comet 21P shows an asymmetric activity with respect to perihelion, which might be due to the nucleus shape, the spin axis orientation, and the distribution of activity on the cometary surface. The maximum of the gas and dust activity was about 24 days before perihelion, similar to the previous apparitions. According to the molecular abundance relative to CN and OH, we confirm that 21P is depleted in C₂, C₃, and NH with respect to CN and to OH. A very good agreement between the abundance of the potential mother molecules measured in the IR (HCN, C₂H₂, and NH₃) and the daughter species from our optical observations has been found. We obtained a high-resolution UVES spectrum of 21P a week after perihelion, and we derived ¹²C/¹³C and ¹⁴N/¹⁵N isotopic ratios of 100 ± 10 and 145 ± 10 from the CN R-branch of the B-X (0, 0) violet band. The ammonia OPR was found to be equal to 1.19 ± 0.03, corresponding to a spin temperature of 27 ± 1 K. All these values agree with those found for several comets of different dynamical types and origins and do not show any peculiarity that could be related to the low-carbon-chain species and ammonia abundances. Our observations favour

a pristine origin for this composition, rather than heterogeneity or evolutionary scenarios of the surface composition.

*Acknowledgments.* The research leading to these results has received funding from the ARC grant for Concerted Research Actions, financed by the Wallonia-Brussels Federation. TRAPPIST-South is a project funded by the Belgian Fonds (National) de la Recherche Scientifique (F.R.S.-FNRS) under grant FRFC 2.5.594.09.F. TRAPPIST-North is a project 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 observations made with ESO Telescopes at the La Silla Paranal Observatory under program DDT proposal 2101.C-5051.

## References

* A'Hearn, M. F., Schleicher, D. G., Millis, R. L., Feldman, P. D., & Thompson, D. T. 1984, AJ, 89, 579

* A'Hearn, M. F., Campins, H., Schleicher, D. G., & Millis, R. L. 1989, ApJ, 347, 1155

* A'Hearn, M. F., Millis, R. C., Schleicher, D. O., Osip, D. J., & Birch, P. V. 1995, Icarus, 118, 223

* Altwegg, K., Balsiger, H., & Fuselier, S. A. 2019, ARA&A, 57, 113

* Arpigny, C. 2003, Science, 301, 1522

* Beaver, J. E., Wagner, R. M., Schleicher, D. G., & Lutz, B. L. 1990, ApJ, 360, 696

* Beech, M. 1986, AJ, 91, 159

* Bessell, M. S. 1990, PASP, 102, 1181

* Bockelée-Morvan, D., Biver, N., Jehin, E., et al. 2008, ApJ, 679, L49

* Bockelée-Morvan, D., Calmonte, U., Charnley, S., et al. 2015, Space Sci. Rev., 197, 47

* Bodewits, D., Farnham, T. L., A'Hearn, M. F., et al. 2014, ApJ, 786, 48

* Bulirsch, R., & Stoer, J. 1964, Numer. Math., 6, 413

* Chambers, J. E. 1999, MNRAS, 304, 793

* Cochran, A. L., & Barker, E. S. 1987, AJ, 93, 239

* Cochran, A. L., & Schleicher, D. G. 1993, Icarus, 105, 235

* Cochran, A., Barker, E., & Gray, C. 2012, Icarus, 218, 144

* Combi, M. R., & Delsemme, A. H. 1980, ApJ, 237, 633

* Combi, M. R., & Feldman, P. D. 1992, Icarus, 97, 260

* Combi, M. R., & Fink, U. 1997, ApJ, 484, 879

* Combi, M. R., Stewart, A. I. F., & Smyth, W. H. 1986, Geophys. Res. Lett., 13, 385

* Combi, M. R., Bertaux, J.-L., Quémerais, E., Ferron, S., & Mäkinen, J. T. T. 2011, ApJ, 734, L6

* Combi, M. R., Mäkinen, T. T., Bertaux, J. L., Quémerais, E., & Ferron, S. 2019, Icarus, 317, 610

* Cowan, J. J., & A'Hearn, M. F. 1979, Moon Planets, 21, 155

* Crovisier, J. 1989, A&A, 213, 459

* DiSanti, M. A., Bonev, B. P., Villanueva, G. L., & Mumma, M. J. 2012, ApJ, 763, 1

* Egal, A., Wiegert, P., Brown, P. G., et al. 2019, Icarus, 330, 123

* Ehlert, S., Moticska, N., & Egal, A. 2019, AJ, 158, 7

* Faggi, S., Mumma, M. J., Villanueva, G. L., Paganini, L., & Lippi, M. 2019, AJ, 158, 254

* Farnham, T., & Schleicher, D. 2005, Icarus, 173, 533

* Farnham, T. L., Schleicher, D. G., & A'Hearn, M. F. 2000, Icarus, 147, 180

* Fink, U. 2009, Icarus, 201, 311

* Fink, U., & Hicks, M. D. 1996, ApJ, 459, 729

* Fray, N., Bénilan, Y., Cottin, H., Gazeau, M.-C., & Crovisier, J. 2005, Planet. Space Sci., 53, 1243

* Gulkis, S., Allen, M., von Allmen, P., et al. 2015, Science, 347, aaa0709

* Hama, T., & Watanabe, N. 2013, Chem. Rev., 113, 8783

* Hama, T., Watanabe, N., Kouchi, A., & Yokoyama, M. 2011, ApJ, 738, L15

* Hama, T., Kouchi, A., & Watanabe, N. 2016, Science, 351, 65

* Hanner, M., Veeder, G., & Tokunaga, A. 1992, AJ, 104, 386

* Haser, L. 1957, Bull. Soc. R. Sci. Liege, 43, 740

* Hässig, M., Altwegg, K., Balsiger, H., et al. 2017, A&A, 605, A50

* Helbert, J., Rauer, H., Boice, D. C., & Huebner, W. F. 2005, A&A, 442, 1107

* Herbig, G. H. 1976, Rev. of cometary spectra, 393, 136

* Hölscher, A. 2015, PhD thesis, Technische Universität, Berlin

* Hyodo, F., Kusaka, S., Wardle, D. A., & Nilsson, M.-C. 2013, Plant Soil, 364, 315

* Jehin, E., Manfroid, J., Hutsemékers, D., Arpigny, C., & Zucconi, J.-M. 2009, Earth Moon Planets, 105, 167
Jehin, E., Gillon, M., Queloz, D., et al. 2011, The Messenger, 145, 2

Jewitt, D. 2015, AJ, 150, 201

Jewitt, D., & Meech, K. J. 1986, ApJ, 310, 937

Kiselev, N. N., Jockers, K., Rosenbush, V. K., et al. 2000, Planet. Space Sci., 48, 1005

Konno, I., & Wyckoff, S. 1989, Adv. Space Res., 9, 163

Królikowska, M., Sitarski, G., & Szutowicz, S. 2001, A&A, 368, 676

Lai, I.-L., Ip, W.-H., Lee, J.-C., et al. 2019, A&A, 630, A17

Lamy, P., & Toth, I. 2009, Icarus, 201, 674

Langland-Shula, L. E., & Smith, G. H. 2011, Icarus, 213, 280

Lara, L.-M., Licandro, J., Oscoz, A., & Motta, V. 2003, A&A, 399, 763

Leibowitz, E. M., & Brosch, N. 1986, Icarus, 68, 430

Levison, H. F., & Duncan, M. J. 1994, Icarus, 108, 18

Lis, D. C., Bockelée-Morvan, D., Güsten, R., et al. 2019, A&A, 625, L5

Luu, J. X. 1993, Icarus, 104, 138

Manfroid, J., Jehin, E., Hutsemékers, D., et al. 2009, A&A, 503, 613

Marshall, D., Rezac, L., Hartogh, P., Zhao, Y., & Attree, N. 2019, A&A, 623, A120

McFadden, L. A., A'Hearn, M. F., Feldman, P. D., et al. 1987, Icarus, 69, 329

McKay, A. J., Cochran, A. L., DiSanti, M. A., et al. 2018, Icarus, 309, 1

McKay, A. J., DiSanti, M. A., Kelley, M. S. P., et al. 2019, AJ, 158, 128

Mianes, P., Grudzinska, S., & Stawikowski, A. 1960, Ann. Astrophys., 23, 788

Moreno, F., Pozuelos, F., Aceituno, F., et al. 2012, ApJ, 752, 136

Moreno, F., Licandro, J., Cabrera-Lavers, A., & Pozuelos, F. J. 2016a, ApJ, 826, 137

Moreno, F., Licandro, J., Cabrera-Lavers, A., & Pozuelos, F. J. 2016b, ApJ, 826, L22

Moreno, F., Pozuelos, F. J., Novaković, B., et al. 2017, ApJ, 837, L3

Moulane, Y., Jehin, E., Opitom, C., et al. 2018, A&A, 619, A156

Mumma, M. J., & Charnley, S. B. 2011, ARA&A, 49, 471

Mumma, M. J., DiSanti, M. A., Russo, N. D., Magee-Sauer, K., & Rettig, T. W. 2000, ApJ, 531, L155

Opitom, C. 2016, Monitoring of the chemical composition of comets in the framework of the TRAPPIST survey, PhD thesis, Université de Liège, Liège, Belgium

Opitom, C., Jehin, E., Manfroid, J., et al. 2015a, A&A, 574, A38

Opitom, C., Jehin, E., Manfroid, J., et al. 2015b, A&A, 584, A121

Opitom, C., Snodgrass, C., Fitzsimmons, A., et al. 2017, MNRAS, 469, S222

Pittichová, J., Woodward, C. E., Kelley, M. S., & Reach, W. T. 2008, AJ, 136, 1127

Pozuelos, F. J., Cabrera-Lavers, A., Licandro, J., & Moreno, F. 2015, ApJ, 806, 102

Pozuelos, F. J., Jehin, E., Moulane, Y., et al. 2018, A&A, 615, A154

Rauer, H., Helbert, J., Arpigny, C., et al. 2003, A&A, 397, 1109

Roth, N., Gibb, E., Dello Russo, N., et al. 2018, AAS/Division for Planetary Sciences Meeting Abstracts, 50, 210.11

Roth, N. X., Gibb, E. L., Bonev, B. P., et al. 2020, AJ, 159, 42

Rousselot, P., Pirali, O., Jehin, E., et al. 2014, ApJ, 780, L17

Rousselot, P., Decock, A., Korsun, P. P., et al. 2015, A&A, 580, A3

Rousselot, P., Opitom, C., Jehin, E., et al. 2019, A&A, 628, A22

Rubin, M., Altwegg, K., Balsiger, H., et al. 2017, A&A, 601, A123

Scarf, F. L., Coroniti, F. V., Kennel, C. F., et al. 1986, Science, 232, 377

Schleicher, D. G. 2006, Icarus, 181, 442

Schleicher, D. 2007, Icarus, 190, 406

Schleicher, D. G. 2008, AJ, 136, 2204

Schleicher, D., & Knight, M. 2018, AAS/Division for Planetary Sciences Meeting Abstracts, 210.12

Schleicher, D. G., Millis, R. L., & Birch, P. V. 1987, A&A, 187, 531

Schleicher, D. G., Millis, R. L., & Birch, P. V. 1998, Icarus, 132, 397

Shinnaka, Y., & Kawakita, H. 2016, AJ, 152, 145

Shinnaka, Y., Kawakita, H., Kobayashi, H., et al. 2011, ApJ, 729, 81

Shinnaka, Y., Kawakita, H., Kobayashi, H., Nagashima, M., & Boice, D. C. 2014, ApJ, 782, L16

Shinnaka, Y., Kawakita, H., Jehin, E., et al. 2016, MNRAS, 462, S195

Shinnaka, Y., Kawakita, H., & Tajitsu, A. 2020, AJ, 159, 203

Solontoi, M., Ivezić, Ž., Jurić, M., et al. 2012, Icarus, 218, 571

Tancredi, G., Fernández, J. A., Rickman, H., & Licandro, J. 2000, A&AS, 146, 73

Von Rosenvinge, T. T., Brandt, J. C., & Farquhar, R. W. 1986, Science, 232, 353

Weaver, H., Chin, G., Bockelee-Morvan, D., et al. 1999, Icarus, 142, 482

Weiler, M. 2012, A&A, 538, A149

Woodney, L., A'Hearn, M., Schleicher, D. G., et al. 2002, Icarus, 157, 193

Yang, B., Hutsemékers, D., Shinnaka, Y., et al. 2018, A&A, 609, L4

Yang, B., Jehin, E., Pozuelos, F. J., et al. 2019, A&A, 631, A168
## Appendix A: Observational circumstances and production rates of comet 21P with TRAPPIST telescopes

<table><caption>**Table A.1.** Observational circumstances of comet 21P with TRAPPIST telescopes.</caption><thead><tr><th rowspan="2">UT date</th><th rowspan="2">r<sub>h</sub><br>(au)</th><th rowspan="2">Δ<br>(au)</th><th rowspan="2">ΔT<br>(Days)</th><th rowspan="2">PA<br>(°)</th><th colspan="4">Gas filters</th><th colspan="4">Dust filters</th><th rowspan="2">Telescope<br>TN/TS</th></tr><tr><th>OH</th><th>NH</th><th>CN</th><th>C<sub>2</sub></th><th>C<sub>3</sub></th><th>BC</th><th>RC</th><th>GC</th><th>Rc</th><th>Ic</th></tr></thead><tbody><tr><td>2018 Jun. 09</td><td>1.61</td><td>1.07</td><td>-93.20</td><td>38.01</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>6</td><td>1</td><td>TN</td></tr><tr><td>2018 Jun. 19</td><td>1.52</td><td>0.94</td><td>-83.25</td><td>40.98</td><td>1</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>1</td><td>1</td><td>TN</td></tr><tr><td>2018 Jun. 22</td><td>1.49</td><td>0.90</td><td>-80.16</td><td>41.98</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>5</td><td>1</td><td>TN</td></tr><tr><td>2018 Jun. 28</td><td>1.44</td><td>0.88</td><td>-74.20</td><td>44.13</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>6</td><td>2</td><td>TN</td></tr><tr><td>2018 Jul. 09</td><td>1.34</td><td>0.78</td><td>-63.30</td><td>48.62</td><td></td><td></td><td>2</td><td>1</td><td></td><td></td><td>1</td><td></td><td>3</td><td>1</td><td>TN</td></tr><tr><td>2018 Jul. 26</td><td>1.21</td><td>0.65</td><td>-46.25</td><td>56.96</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>3</td><td>1</td><td>TN</td></tr><tr><td>2018 Jul. 30</td><td>1.18</td><td>0.62</td><td>-42.25</td><td>59.13</td><td>1</td><td></td><td>2</td><td>2</td><td></td><td></td><td></td><td></td><td>4</td><td></td><td>TN</td></tr><tr><td>2018 Aug. 17</td><td>1.07</td><td>0.48</td><td>-24.25</td><td>69.44</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>6</td><td>1</td><td>TN</td></tr><tr><td>2018 Aug. 18</td><td>1.07</td><td>0.49</td><td>-23.08</td><td>69.99</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>5</td><td>1</td><td>TN</td></tr><tr><td>2018 Aug. 23</td><td>1.04</td><td>0.45</td><td>-18.13</td><td>72.65</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>5</td><td>1</td><td>TN</td></tr><tr><td>2018 Aug. 29</td><td>1.03</td><td>0.42</td><td>-12.11</td><td>75.38</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>5</td><td>1</td><td>TN</td></tr><tr><td>2018 Sep. 05</td><td>1.01</td><td>0.39</td><td>-05.10</td><td>77.67</td><td>1</td><td></td><td>3</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>5</td><td>1</td><td>TN</td></tr><tr><td>2018 Sep. 09</td><td>1.01</td><td>0.39</td><td>-01.10</td><td>78.01</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>8</td><td>4</td><td>TN</td></tr><tr><td>2018 Sep. 15</td><td>1.01</td><td>0.39</td><td>+05.09</td><td>77.40</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>5</td><td>1</td><td>TS</td></tr><tr><td>2018 Sep. 17</td><td>1.01</td><td>0.40</td><td>+06.90</td><td>76.83</td><td>1</td><td>1</td><td>3</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>11</td><td>2</td><td>TN</td></tr><tr><td>2018 Sep. 17</td><td>1.01</td><td>0.40</td><td>+06.90</td><td>76.83</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Sep. 18</td><td>1.01</td><td>0.40</td><td>+08.13</td><td>76.53</td><td></td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Sep. 20</td><td>1.02</td><td>0.40</td><td>+10.12</td><td>75.84</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Sep. 21</td><td>1.02</td><td>0.41</td><td>+11.12</td><td>75.45</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>3</td><td>1</td><td>TN</td></tr><tr><td>2018 Sep. 21</td><td>1.02</td><td>0.41</td><td>+11.12</td><td>75.45</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Sep. 22</td><td>1.02</td><td>0.41</td><td>+12.12</td><td>75.04</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Sep. 23</td><td>1.02</td><td>0.41</td><td>+12.95</td><td>74.61</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>3</td><td>1</td><td>TN</td></tr><tr><td>2018 Sep. 25</td><td>1.03</td><td>0.42</td><td>+15.10</td><td>73.63</td><td>1</td><td>1</td><td>2</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Sep. 29</td><td>1.04</td><td>0.44</td><td>+18.90</td><td>71.61</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>6</td><td>2</td><td>TN</td></tr><tr><td>2018 Oct. 01</td><td>1.05</td><td>0.45</td><td>+21.10</td><td>70.49</td><td>2</td><td>1</td><td>2</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>3</td><td>TS</td></tr><tr><td>2018 Oct. 02</td><td>1.06</td><td>0.46</td><td>+22.10</td><td>69.92</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td></td><td>TS</td></tr><tr><td>2018 Oct. 04</td><td>1.07</td><td>0.47</td><td>+23.90</td><td>68.69</td><td>1</td><td>1</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>6</td><td>2</td><td>TN</td></tr><tr><td>2018 Oct. 07</td><td>1.08</td><td>0.49</td><td>+27.13</td><td>66.92</td><td></td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>1</td><td></td><td>TS</td></tr><tr><td>2018 Oct. 08</td><td>1.08</td><td>0.49</td><td>+27.90</td><td>66.37</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>11</td><td>1</td><td>TN</td></tr><tr><td>2018 Oct. 12</td><td>1.10</td><td>0.51</td><td>+31.90</td><td>64.55</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>2</td><td>1</td><td>10</td><td>2</td><td>TN</td></tr><tr><td>2018 Oct. 14</td><td>1.12</td><td>0.53</td><td>+34.10</td><td>62.73</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>1</td><td>TS</td></tr><tr><td>2018 Oct. 15</td><td>1.12</td><td>0.54</td><td>+35.10</td><td>62.12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>5</td><td>1</td><td>TS</td></tr><tr><td>2018 Oct. 18</td><td>1.14</td><td>0.56</td><td>+38.11</td><td>60.33</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>3</td><td>2</td><td>TS</td></tr><tr><td>2018 Oct. 18</td><td>1.14</td><td>0.56</td><td>+37.90</td><td>60.33</td><td></td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>1</td><td></td><td>TN</td></tr><tr><td>2018 Oct. 21</td><td>1.16</td><td>0.58</td><td>+40.95</td><td>58.55</td><td></td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>3</td><td></td><td>TN</td></tr><tr><td>2018 Oct. 22</td><td>1.17</td><td>0.59</td><td>+42.06</td><td>57.97</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>5</td><td>2</td><td>TS</td></tr><tr><td>2018 Oct. 25</td><td>1.19</td><td>0.61</td><td>+45.10</td><td>56.25</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td><td>1</td><td>TS</td></tr><tr><td>2018 Oct. 28</td><td>1.21</td><td>0.63</td><td>+48.09</td><td>54.55</td><td></td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>1</td><td>1</td><td>TS</td></tr><tr><td>2018 Nov. 01</td><td>1.24</td><td>0.66</td><td>+52.04</td><td>52.37</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>4</td><td>1</td><td>TS</td></tr><tr><td>2018 Nov. 05</td><td>1.27</td><td>0.68</td><td>+56.05</td><td>50.80</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>4</td><td>1</td><td>TS</td></tr><tr><td>2018 Nov. 09</td><td>1.31</td><td>0.71</td><td>+59.98</td><td>48.27</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>4</td><td>1</td><td>TS</td></tr><tr><td>2018 Nov. 14</td><td>1.35</td><td>0.74</td><td>+65.10</td><td>45.84</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>5</td><td>1</td><td>TS</td></tr><tr><td>2018 Nov. 30</td><td>1.49</td><td>0.85</td><td>+81.08</td><td>38.84</td><td></td><td></td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td>3</td><td>1</td><td>TS</td></tr><tr><td>2018 Dec. 07</td><td>1.55</td><td>0.90</td><td>+88.12</td><td>36.12</td><td></td><td></td><td>1</td><td></td><td></td><td></td><td></td><td></td><td>3</td><td>1</td><td>TS</td></tr><tr><td>2018 Dec. 10</td><td>1.58</td><td>0.92</td><td>+91.10</td><td>35.01</td><td></td><td></td><td>1</td><td></td><td></td><td></td><td></td><td></td><td>3</td><td>1</td><td>TS</td></tr><tr><td>2018 Dec. 19</td><td>1.66</td><td>0.98</td><td>+100.05</td><td>32.04</td><td></td><td></td><td>1</td><td></td><td></td><td></td><td></td><td></td><td>4</td><td>1</td><td>TS</td></tr><tr><td>2018 Dec. 29</td><td>1.75</td><td>1.06</td><td>+110.12</td><td>29.28</td><td></td><td></td><td>1</td><td></td><td></td><td></td><td></td><td></td><td>1</td><td>1</td><td>TS</td></tr><tr><td>2019 Jan. 14</td><td>1.90</td><td>1.20</td><td>+126.13</td><td>26.20</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>2</td><td></td><td>TS</td></tr><tr><td>2019 Jan. 15</td><td>1.91</td><td>1.21</td><td>+127.10</td><td>26.08</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>3</td><td>1</td><td>TS</td></tr><tr><td>2019 Jan. 29</td><td>2.04</td><td>1.37</td><td>+141.12</td><td>24.74</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>4</td><td>1</td><td>TS</td></tr><tr><td>2019 Feb. 02</td><td>2.08</td><td>1.42</td><td>+145.08</td><td>24.52</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>4</td><td>1</td><td>TS</td></tr><tr><td>2019 Feb. 04</td><td>2.10</td><td>1.44</td><td>+147.11</td><td>24.43</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>4</td><td>4</td><td>TS</td></tr></tbody></table>

Notes. $r_h$ and $\Delta$ are the heliocentric and geocentric distances, $\Delta T$ is the time to perihelion in days, (-) for pre-perihelion and (+) for post-perihelion. PA is the solar phase angle.
<table><caption>**Table A.2.** OH, NH, CN, C₂, and C₃ production rates and A(θ=0)fρ measurements for comet 21P with the TN and TS telescopes.</caption><thead><tr><th rowspan="2">UT date</th><th rowspan="2">r<sub>h</sub><br>(au)</th><th colspan="5">Production rates (× 10<sup>24</sup> mol s<sup>-1</sup>)</th><th colspan="5">A(θ=0)fρ</th><th rowspan="2">TN/TS</th></tr><tr><th>OH</th><th>NH</th><th>CN</th><th>C<sub>2</sub></th><th>C<sub>3</sub></th><th>BC</th><th>RC</th><th>GC</th><th>Rc</th><th>Ic</th></tr></thead><tbody><tr><td>2018 Jun. 09</td><td>1.61</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>236.5 ± 14.0</td><td>265.5 ± 18.7</td><td>TN</td></tr><tr><td>2018 Jun. 19</td><td>1.52</td><td>—</td><td>—</td><td>25.20 ± 0.52</td><td>4.30 ± 0.65</td><td>—</td><td>234.8 ± 13.6</td><td>351.1 ± 17.5</td><td>—</td><td>—</td><td>326.9 ± 15.6</td><td>TN</td></tr><tr><td>2018 Jun. 22</td><td>1.49</td><td>6200 ± 253</td><td>—</td><td>28.00 ± 0.53</td><td>—</td><td>1.50 ± 0.20</td><td>—</td><td>—</td><td>—</td><td>333.6 ± 16.8</td><td>367.5 ± 18.4</td><td>TN</td></tr><tr><td>2018 Jun. 28</td><td>1.44</td><td>9560 ± 365</td><td>—</td><td>32.00 ± 0.64</td><td>—</td><td>—</td><td>—</td><td>673.7 ± 21.4</td><td>—</td><td>411.0 ± 13.7</td><td>475.7 ± 12.9</td><td>TN</td></tr><tr><td>2018 Jul. 26</td><td>1.21</td><td>15 200 ± 322</td><td>31.30 ± 4.90</td><td>37.80 ± 0.54</td><td>8.40 ± 0.66</td><td>1.08 ± 0.20</td><td>618.6 ± 36.4</td><td>801.5 ± 12.9</td><td>—</td><td>687.8 ± 11.5</td><td>758.2 ± 15.9</td><td>TN</td></tr><tr><td>2018 Jul. 30</td><td>1.18</td><td>18 100 ± 278</td><td>—</td><td>45.50 ± 0.63</td><td>12.60 ± 0.65</td><td>—</td><td>—</td><td>—</td><td>—</td><td>959.8 ± 15.8</td><td>—</td><td>TN</td></tr><tr><td>2018 Aug. 17</td><td>1.07</td><td>28 300 ± 558</td><td>54.10 ± 4.59</td><td>52.00 ± 0.65</td><td>17.90 ± 0.64</td><td>2.52 ± 0.19</td><td>1087.6 ± 13.7</td><td>1646.1 ± 12.7</td><td>—</td><td>1502.7 ± 12.1</td><td>1609.3 ± 12.7</td><td>TN</td></tr><tr><td>2018 Aug. 18</td><td>1.07</td><td>25 600 ± 276</td><td>—</td><td>50.80 ± 0.64</td><td>18.00 ± 0.62</td><td>2.47 ± 0.17</td><td>—</td><td>—</td><td>—</td><td>1522.1 ± 11.6</td><td>1648.9 ± 10.3</td><td>TN</td></tr><tr><td>2018 Aug. 23</td><td>1.04</td><td>27 300 ± 304</td><td>50.20 ± 3.60</td><td>55.20 ± 0.67</td><td>19.20 ± 0.64</td><td>2.55 ± 0.17</td><td>1028.9 ± 12.5</td><td>1570.7 ± 11.8</td><td>—</td><td>1463.7 ± 11.7</td><td>1570.2 ± 10.0</td><td>TN</td></tr><tr><td>2018 Aug. 29</td><td>1.03</td><td>23 000 ± 286</td><td>43.50 ± 4.07</td><td>42.50 ± 0.62</td><td>13.70 ± 0.62</td><td>1.66 ± 0.18</td><td>887.6 ± 29.6</td><td>1229.9 ± 12.5</td><td>—</td><td>1100.0 ± 12.7</td><td>1176.3 ± 16.1</td><td>TN</td></tr><tr><td>2018 Sep. 05</td><td>1.01</td><td>20 800 ± 312</td><td>—</td><td>38.20 ± 0.64</td><td>14.00 ± 0.63</td><td>1.95 ± 0.17</td><td>651.2 ± 12.4</td><td>—</td><td>—</td><td>933.1 ± 10.5</td><td>984.0 ± 14.8</td><td>TN</td></tr><tr><td>2018 Sep. 09</td><td>1.01</td><td>23 600 ± 328</td><td>56.60 ± 3.86</td><td>43.90 ± 0.68</td><td>16.70 ± 0.64</td><td>2.26 ± 0.17</td><td>723.3 ± 12.4</td><td>1103.1 ± 15.7</td><td>815.5 ± 12.4</td><td>1062.9 ± 12.7</td><td>1123.3 ± 12.1</td><td>TN</td></tr><tr><td>2018 Sep. 17</td><td>1.01</td><td>17 100 ± 300</td><td>39.80 ± 3.77</td><td>37.10 ± 0.66</td><td>11.60 ± 0.64</td><td>1.91 ± 0.18</td><td>574.6 ± 14.1</td><td>861.6 ± 12.4</td><td>642.0 ± 10.7</td><td>799.0 ± 15.1</td><td>833.7 ± 14.0</td><td>TN</td></tr><tr><td>2018 Sep. 17</td><td>1.01</td><td>14 700 ± 326</td><td>38.30 ± 3.84</td><td>37.60 ± 0.53</td><td>12.30 ± 0.63</td><td>2.02 ± 0.22</td><td>571.6 ± 19.6</td><td>851.4 ± 17.2</td><td>599.9 ± 19.2</td><td>767.0 ± 16.2</td><td>832.1 ± 14.8</td><td>TS</td></tr><tr><td>2018 Sep. 18</td><td>1.01</td><td>—</td><td>—</td><td>30.00 ± 0.52</td><td>11.00 ± 0.65</td><td>2.07 ± 0.28</td><td>459.5 ± 34.2</td><td>872.2 ± 22.2</td><td>595.9 ± 14.8</td><td>733.3 ± 14.8</td><td>849.6 ± 20.4</td><td>TS</td></tr><tr><td>2018 Sep. 21</td><td>1.02</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>604.7 ± 16.6</td><td>652.2 ± 14.9</td><td>TN</td></tr><tr><td>2018 Sep. 21</td><td>1.02</td><td>1280 ± 212</td><td>28.50 ± 3.79</td><td>33.00 ± 0.52</td><td>10.90 ± 0.63</td><td>1.72 ± 0.22</td><td>542.3 ± 17.4</td><td>808.3 ± 12.5</td><td>606.4 ± 19.9</td><td>658.3 ± 14.9</td><td>779.5 ± 21.4</td><td>TS</td></tr><tr><td>2018 Sep. 22</td><td>1.02</td><td>13 600 ± 216</td><td>29.40 ± 3.84</td><td>32.80 ± 0.51</td><td>10.80 ± 0.65</td><td>1.64 ± 0.23</td><td>510.9 ± 17.4</td><td>744.0 ± 15.0</td><td>569.2 ± 12.5</td><td>649.8 ± 12.5</td><td>740.3 ± 14.7</td><td>TS</td></tr><tr><td>2018 Sep. 23</td><td>1.02</td><td>12 800 ± 249</td><td>—</td><td>28.60 ± 0.61</td><td>9.48 ± 0.62</td><td>—</td><td>394.5 ± 12.5</td><td>596.1 ± 12.5</td><td>—</td><td>—</td><td>629.3 ± 15.2</td><td>TN</td></tr><tr><td>2018 Sep. 25</td><td>1.03</td><td>10 880 ± 216</td><td>29.00 ± 4.15</td><td>27.30 ± 0.54</td><td>9.62 ± 0.76</td><td>1.24 ± 0.31</td><td>406.8 ± 22.4</td><td>611.7 ± 19.9</td><td>460.0 ± 20.1</td><td>540.5 ± 14.8</td><td>611.6 ± 15.3</td><td>TS</td></tr><tr><td>2018 Sep. 29</td><td>1.04</td><td>11 300 ± 254</td><td>23.80 ± 4.11</td><td>24.50 ± 0.58</td><td>7.21 ± 0.67</td><td>0.79 ± 0.19</td><td>368.1 ± 20.4</td><td>517.1 ± 15.1</td><td>—</td><td>530.7 ± 17.8</td><td>535.4 ± 15.2</td><td>TN</td></tr><tr><td>2018 Oct. 01</td><td>1.05</td><td>10 220 ± 217</td><td>18.40 ± 3.68</td><td>42.40 ± 0.52</td><td>7.84 ± 0.64</td><td>1.13 ± 0.23</td><td>381.4 ± 23.2</td><td>490.7 ± 20.6</td><td>368.1 ± 21.9</td><td>418.3 ± 15.3</td><td>478.8 ± 15.7</td><td>TS</td></tr><tr><td>2018 Oct. 02</td><td>1.06</td><td>9860 ± 182</td><td>19.50 ± 3.27</td><td>21.70 ± 0.52</td><td>6.55 ± 0.71</td><td>1.04 ± 0.21</td><td>350.8 ± 15.3</td><td>—</td><td>—</td><td>—</td><td>—</td><td>TS</td></tr><tr><td>2018 Oct. 04</td><td>1.07</td><td>9460 ± 268</td><td>—</td><td>20.20 ± 0.55</td><td>5.95 ± 0.63</td><td>0.70 ± 0.18</td><td>282.6 ± 14.0</td><td>412.3 ± 12.2</td><td>—</td><td>411.2 ± 12.2</td><td>437.4 ± 15.4</td><td>TN</td></tr><tr><td>2018 Oct. 08</td><td>1.08</td><td>8340 ± 248</td><td>10.70 ± 3.53</td><td>16.80 ± 0.53</td><td>4.03 ± 0.61</td><td>0.36 ± 0.17</td><td>—</td><td>318.5 ± 19.9</td><td>—</td><td>350.5 ± 12.2</td><td>350.5 ± 18.3</td><td>TN</td></tr><tr><td>2018 Oct. 12</td><td>1.10</td><td>5950 ± 248</td><td>—</td><td>15.80 ± 0.52</td><td>4.24 ± 0.63</td><td>—</td><td>199.7 ± 12.3</td><td>280.1 ± 19.1</td><td>—</td><td>294.6 ± 12.8</td><td>296.5 ± 20.4</td><td>TN</td></tr><tr><td>2018 Oct. 14</td><td>1.12</td><td>4730 ± 247</td><td>11.00 ± 4.16</td><td>13.30 ± 0.52</td><td>3.90 ± 0.67</td><td>—</td><td>190.9 ± 25.7</td><td>—</td><td>—</td><td>234.0 ± 17.3</td><td>269.8 ± 20.5</td><td>TS</td></tr><tr><td>2018 Oct. 15</td><td>1.12</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>166.9 ± 18.3</td><td>257.9 ± 18.6</td><td>192.1 ± 21.2</td><td>226.1 ± 13.3</td><td>—</td><td>TS</td></tr><tr><td>2018 Oct. 18</td><td>1.14</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>151.0 ± 17.2</td><td>—</td><td>—</td><td>208.2 ± 12.4</td><td>230.8 ± 11.7</td><td>TS</td></tr><tr><td>2018 Oct. 18</td><td>1.14</td><td>—</td><td>—</td><td>10.80 ± 0.54</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>229.7 ± 14.3</td><td>—</td><td>TN</td></tr><tr><td>2018 Oct. 21</td><td>1.16</td><td>—</td><td>—</td><td>9.87 ± 0.51</td><td>3.02 ± 0.65</td><td>—</td><td>—</td><td>—</td><td>—</td><td>198.2 ± 19.4</td><td>—</td><td>TN</td></tr><tr><td>2018 Oct. 22</td><td>1.17</td><td>3500 ± 188</td><td>—</td><td>9.04 ± 0.50</td><td>2.30 ± 0.65</td><td>0.34 ± 0.21</td><td>137.6 ± 21.1</td><td>149.0 ± 26.1</td><td>—</td><td>114.4 ± 13.3</td><td>127.9 ± 14.4</td><td>TS</td></tr><tr><td>2018 Oct. 25</td><td>1.19</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>111.5 ± 20.0</td><td>—</td><td>TS</td></tr><tr><td>2018 Oct. 28</td><td>1.21</td><td>—</td><td>—</td><td>6.00 ± 0.55</td><td>—</td><td>—</td><td>—</td><td>151.2 ± 13.8</td><td>117.2 ± 16.5</td><td>137.6 ± 12.1</td><td>144.7 ± 14.1</td><td>TS</td></tr><tr><td>2018 Nov. 01</td><td>1.24</td><td>2320 ± 234</td><td>—</td><td>6.88 ± 0.56</td><td>—</td><td>—</td><td>—</td><td>108.1 ± 18.6</td><td>85.7 ± 23.6</td><td>109.6 ± 15.7</td><td>123.2 ± 12.4</td><td>TS</td></tr><tr><td>2018 Nov. 05</td><td>1.27</td><td>2010 ± 191</td><td>—</td><td>5.62 ± 0.52</td><td>1.16 ± 0.66</td><td>—</td><td>95.7 ± 20.4</td><td>120.5 ± 21.8</td><td>50.1 ± 16.2</td><td>121.1 ± 15.4</td><td>139.0 ± 15.0</td><td>TS</td></tr><tr><td>2018 Nov. 09</td><td>1.31</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>88.8 ± 13.4</td><td>129.3 ± 15.6</td><td>TS</td></tr><tr><td>2018 Nov. 14</td><td>1.35</td><td>—</td><td>—</td><td>4.39 ± 0.55</td><td>1.27 ± 0.69</td><td>—</td><td>85.5 ± 13.2</td><td>88.7 ± 18.7</td><td>87.2 ± 14.6</td><td>86.9 ± 13.6</td><td>91.3 ± 12.0</td><td>TS</td></tr><tr><td>2018 Dec. 07</td><td>1.55</td><td>—</td><td>—</td><td>3.30 ± 0.52</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>42.9 ± 10.5</td><td>71.7 ± 13.2</td><td>TS</td></tr><tr><td>2018 Dec. 10</td><td>1.58</td><td>—</td><td>—</td><td>3.00 ± 0.56</td><td>—</td><td>—</td><td>—</td><td>75.9 ± 14.2</td><td>—</td><td>28.1 ± 13.9</td><td>—</td><td>TS</td></tr><tr><td>2018 Dec. 19</td><td>1.66</td><td>—</td><td>—</td><td>3.30 ± 0.54</td><td>—</td><td>—</td><td>—</td><td>80.6 ± 12.5</td><td>—</td><td>27.2 ± 11.5</td><td>172.3 ± 18.2</td><td>TS</td></tr><tr><td>2018 Dec. 29</td><td>1.75</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>57.9 ± 18.1</td><td>TS</td></tr><tr><td>2019 Jan. 14</td><td>1.90</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>60.7 ± 11.1</td><td>—</td><td>TS</td></tr><tr><td>2019 Jan. 15</td><td>1.91</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>57.2 ± 11.0</td><td>60.4 ± 11.6</td><td>TS</td></tr><tr><td>2019 Jan. 29</td><td>2.04</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>85.5 ± 11.4</td><td>93.5 ± 11.5</td><td>TS</td></tr><tr><td>2019 Feb. 02</td><td>2.08</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>83.7 ± 12.5</td><td>91.6 ± 11.3</td><td>TS</td></tr><tr><td>2019 Feb. 04</td><td>2.10</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>71.5 ± 22.0</td><td>96.3 ± 17.1</td><td>TS</td></tr></tbody></table>

Notes. $r_h$ is the heliocentric distance. The A(0)fρ values are printed at 10 000 km from the nucleus and corrected for the phase-angle effect.