Monitoring Volatile Evolution in Disrupting Comet D/2021 A1 (Leonard) with NOEMA and APEX
Pith reviewed 2026-05-13 19:04 UTC · model grok-4.3
The pith
Disruption and fragmentation in comet D/2021 A1 drove changes in volatile mixing ratios beyond what solar distance alone predicts.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Comet D/2021 A1 exhibited a clear rise in CS mixing ratios from 0.02% to 0.10% as heliocentric distance decreased, while HCN remained near 0.07% with added variability once cross-facility data were included; line kinematics showed HCN released near the nucleus and CS from a distributed source, with both species displaying changes tied to mid-December outbursts and fragmentation.
What carries the argument
Mixing ratios of HCN and CS relative to H2O, derived from line intensities and spatial constraints, that separate nucleus versus distributed sources and track changes across epochs.
If this is right
- CS production increases steadily inward because the molecule forms in the coma from a distributed parent.
- HCN shows no robust single trend with distance once all data are merged, indicating multiple release mechanisms.
- Outbursts and fragmentation cause short-term spikes in both species that solar insolation models alone cannot explain.
- Multi-facility, multi-epoch coverage is required to resolve species-dependent responses on timescales of days.
Where Pith is reading between the lines
- Similar monitoring of other fragmenting comets could reveal whether buried ices are preferentially released during breakup events.
- The observed CS rise may imply that distributed-source chemistry becomes more efficient closer to the Sun in active comets.
- Accounting for disruption effects could improve predictions of volatile delivery to inner planets during cometary passages.
Load-bearing premise
That HCN measurements from different telescopes can be combined directly without calibration offsets large enough to erase the apparent lack of monotonic trend with distance.
What would settle it
A new campaign using a single calibrated instrument across the same distance range that finds HCN mixing ratios either constant or strictly decreasing with no outburst-linked spikes would undermine the claim that disruption drives the observed variability.
Figures
read the original abstract
We report a pre-perihelion survey of volatile emissions from comet D/2021 A1 (Leonard) with the Northern Extended Millimeter Array (NOEMA; UT 2021 Nov. 5, 21, and Dec. 1) and the Atacama Pathfinder Experiment (APEX; UT 2021 Dec. 9-10), spanning heliocentric distances ($r_H$) from 1.3 to 0.80 au. We securely detected HCN and CS and place 3$\sigma$ upper limits on CH$_3$OH, H$_2$CO, and CO abundances. Line kinematics and NOEMA spatial constraints indicate that HCN was released at or near the nucleus (parent scale length $<300$ km), while CS showed higher gas expansion velocities and mixing ratios that increased with decreasing $r_H$ $-$ consistent with production from a distributed source. Across our campaign, CS mixing ratios relative to H$_2$O increased by a factor of $\sim$5, from $0.02 \pm 0.01\%$ at $r_H$ = 1.3 au to $0.10\pm0.02\%$ by $r_H$ = 0.80 au. HCN mixing ratios in our data rose modestly, from $0.04 \pm 0.02\%$ at $r_H$ = 1.3 au to $0.07 \pm 0.02\%$ by $r_H$ = 0.81 au. However, contemporaneous measurements from other facilities placed HCN consistently at a higher absolute level ($\sim\!0.08\%$) with additional variability. Once cross-facility measurements were included, the HCN abundance showed no statistically robust monotonic dependence on $r_H$. Variability in both species during the mid-December outbursts and fragmentation suggests that D/2021 A1's volatile evolution reflected not only solar insolation but also disruption processes, underscoring the value of multi-epoch, multi-instrument monitoring to capture rapid, species-dependent changes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports pre-perihelion millimeter observations of comet D/2021 A1 (Leonard) with NOEMA (Nov 5, 21 and Dec 1 2021) and APEX (Dec 9-10 2021) at heliocentric distances 1.3–0.80 au. Secure detections of HCN and CS are presented together with 3σ upper limits on CH₃OH, H₂CO and CO. Line kinematics and NOEMA spatial constraints indicate HCN release at or near the nucleus (parent scale length <300 km) while CS exhibits higher expansion velocities and increasing mixing ratios with decreasing r_H, consistent with a distributed source. CS/H₂O rises by a factor of ~5 (0.02±0.01% to 0.10±0.02%); the authors’ HCN/H₂O values rise modestly (0.04±0.02% to 0.07±0.02%), but inclusion of contemporaneous external data (~0.08%) removes any monotonic r_H trend and reveals additional variability during mid-December outbursts and fragmentation. The paper concludes that volatile evolution reflects both solar insolation and disruption processes.
Significance. If the cross-facility HCN data can be placed on a common absolute scale, the work provides concrete evidence that fragmentation events drive species-dependent volatile changes on short timescales, beyond the usual r_H dependence. The kinematic and spatial arguments for parent versus distributed sources are a clear strength, and the multi-epoch, multi-instrument approach demonstrates the scientific return from coordinated monitoring of disrupting comets.
major comments (1)
- [Abstract] Abstract: the interpretation that mid-December variability in HCN (and CS) indicates disruption-driven evolution beyond insolation rests on combining the authors’ NOEMA/APEX mixing ratios (0.04–0.07%) with external measurements (~0.08%). No section quantifies beam-filling factors, absolute flux calibration offsets, or differences in production-rate modeling assumptions between facilities; if these systematics are comparable to the reported factor-of-2 swings, the physical conclusion is not yet secured.
Simulated Author's Rebuttal
We thank the referee for their constructive review and positive assessment of the manuscript's significance. We address the single major comment below and will incorporate revisions to strengthen the cross-facility comparison.
read point-by-point responses
-
Referee: [Abstract] Abstract: the interpretation that mid-December variability in HCN (and CS) indicates disruption-driven evolution beyond insolation rests on combining the authors’ NOEMA/APEX mixing ratios (0.04–0.07%) with external measurements (~0.08%). No section quantifies beam-filling factors, absolute flux calibration offsets, or differences in production-rate modeling assumptions between facilities; if these systematics are comparable to the reported factor-of-2 swings, the physical conclusion is not yet secured.
Authors: We agree that the abstract's interpretation would be strengthened by explicit quantification of cross-facility systematics. In the revised manuscript we will add a new subsection (Section 4.2) that (i) estimates beam-filling factors using the NOEMA synthesized beam and the comet's known angular size at each epoch, (ii) propagates absolute flux calibration uncertainties (typically 10–15% for NOEMA and 20% for APEX) into the mixing ratios, and (iii) compares production-rate modeling assumptions (Haser vs. vectorial models, expansion velocity prescriptions) with those adopted in the external data sets. Preliminary estimates show that the combined uncertainties are ~25–30%, smaller than the observed factor-of-~2 variability in HCN. We will also revise the abstract and conclusions to state that the disruption-driven interpretation is supported once these systematics are accounted for, while noting that contemporaneous multi-facility observations would further reduce residual offsets. These changes directly address the concern and secure the physical conclusion. revision: yes
Circularity Check
No significant circularity; claims rest on direct empirical measurements
full rationale
The paper presents observational results from NOEMA and APEX telescope data, reporting detected mixing ratios for HCN and CS, upper limits for other species, and trends with heliocentric distance. These are derived from line detections, kinematics, and spatial constraints without any equations, fitted parameters, or self-citations that reduce a claimed prediction back to the input data by construction. Cross-facility comparisons are included as external data points to contextualize variability, but the paper does not define its own results in terms of those or invoke uniqueness theorems. The central claims about disruption-driven evolution are empirical interpretations of measured abundances, not forced by internal definitions or self-referential loops. This is a standard observational astronomy paper with no load-bearing circular steps.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard assumptions about cometary gas expansion velocities and parent scale lengths for interpreting line profiles and spatial constraints.
Reference graph
Works this paper leans on
-
[1]
Altwegg, K., Balsiger, H., H¨ anni, N., et al. 2020, Evidence of ammonium salts in comet 67P as explanation for the nitrogen depletion in cometary comae, NatAs, 4, 533
work page 2020
-
[2]
Bergman, P., Lerner, M., Olofsson, A., et al. 2022, Emission from hcn and ch3oh in comets-onsala 20-m observations and radiative transfer modelling, Astronomy & Astrophysics, 660, A118
work page 2022
-
[3]
Biver, N., Bockel´ ee-Morvan, D., Colom, P., et al. 2011, Molecular investigations of comets C/2002 X5 (Kudo-Fujikawa), C/2002 V1 (NEAT), and C/2006 P1 (McNaught) at small heliocentric distances, Astronomy & Astrophysics, 528, A142
work page 2011
-
[4]
Biver, N., Boissier, J., Bockel´ ee-Morvan, D., et al. 2022, Observations of comet C/2020 F3 (NEOWISE) with IRAM telescopes, Astronomy & Astrophysics, 668, A171
work page 2022
-
[5]
1999, Spectroscopic Monitoring of Comet C/1996 B2 (Hyakutake) with the JCMT and IRAM Radio
Biver, N., Bockel´ ee-Morvan, D., Crovisier, J., et al. 1999, Spectroscopic Monitoring of Comet C/1996 B2 (Hyakutake) with the JCMT and IRAM Radio
work page 1999
-
[6]
Biver, N., Bockel´ ee-Morvan, D., Handzlik, B., et al. 2024, Chemical composition of comets C/2021 A1 (Leonard) and C/2022 E3 (ZTF) from radio spectroscopy and the abundance of HCOOH and HNCO in comets, Astronomy & Astrophysics, 690, A271 Bockel´ ee-Morvan, D., & Biver, N. 2017, The composition of cometary ices, PTRSA, 375, 20160252, doi: 10.1098/rsta.201...
-
[7]
Weaver, H. A. 2004, in Comets II, ed. H. U. Keller & H. A. Weaver (University of Arizona Press), 391 Bockel´ ee-Morvan, D., Lis, D. C., Wink, J. E., et al. 2000, New molecules found in comet C/1995 O1 (Hale-Bopp). Investigating the link between cometary and interstellar material, A&A, 353, 1101 Bockel´ ee-Morvan, D., Boissier, J., Biver, N., & Crovisier, ...
-
[8]
Boissier, J., Bockel´ ee-Morvan, D., Biver, N., et al. 2007, Interferometric imaging of the sulfur-bearing molecules H2S, SO, and CS in comet C/1995 O1 (Hale-Bopp), A&A, 475, 1131
work page 2007
-
[9]
P., Dello Russo, N., DiSanti, M
Bonev, B. P., Dello Russo, N., DiSanti, M. A., et al. 2021, First Comet Observations with NIRSPEC-2 at Keck: Outgassing Sources of Parent Volatiles and Abundances Based on Alternative Taxonomic Compositional Baselines in 46P/Wirtanen, PSJ, 2, 45
work page 2021
-
[10]
Bregman, J. D., Witteborn, F. C., Allamandola, L. J., et al. 1987, Airborn and groundbased spectrophotometry of comet P/Halley from 5-13 micrometers, Astronomy and Astrophysics, 187, 616
work page 1987
-
[11]
Brooke, T., Weaver, H. A., Chin, G., et al. 2003, Spectroscopy of Comet Hale–Bopp in the infrared, Icarus, 166, 167
work page 2003
-
[12]
Ferron, S. 2023, Water production rates from SOHO/SWAN observations of comets C/2020S3 (Erasmus), C/2021 A1 (Leonard) and C/2021 O3 (PanSTARRS), Icarus, 398, 115543, doi: https://doi.org/10.1016/j.icarus.2023.115543
-
[13]
Cordiner, M. A., Remijan, A. J., Boissier, J., et al. 2014, Mapping the Release of Volatiles in the Inner Comae of Comets C/2012 F6 (Lemmon) and C/2012 S1 (ISON) Using the Atacama Large Millimeter/Submillimeter
work page 2014
-
[14]
Cordiner, M. A., Boissier, J., Charnley, S. B., et al. 2017b, ALMA Mapping of Rapid Gas and Dust Variations in Comet C/2012 S1 (ISON):New Insights into the Origin of Cometary HNC, ApJ, 838, 147
work page 2012
-
[15]
Cordiner, M. A., Coulson, I. M., Garcia-Berrios, E., et al. 2022, A SUBLIME 3D Model for Cometary Coma Emission: The Hypervolatile-rich Comet C/2016 R2 (PanSTARRS), The Astrophysical Journal, 929, 38, doi: 10.3847/1538-4357/ac5893
-
[16]
Cordiner, M. A., Roth, N. X., Milam, S. N., et al. 2023, Gas Sources from the Coma and Nucleus of Comet 46P/Wirtanen Observed Using ALMA, https://arxiv.org/abs/2305.04822
-
[17]
2008, Distributed Sources in Comets, Space Science Reviews, 138, 179
Cottin, H., & Fray, N. 2008, Distributed Sources in Comets, Space Science Reviews, 138, 179
work page 2008
-
[18]
2021, Central Bureau Electronic Telegrams, 5087, de Val-Borro, M., Milam, S., Cordiner, M., et al
Crovisier, J., Biver, N., & Bockel´ ee-Morvan, D. 2021, Central Bureau Electronic Telegrams, 5087, de Val-Borro, M., Milam, S., Cordiner, M., et al. 2018, Measuring molecular abundances in comet C/2014 Q2 (Lovejoy) using the APEX telescope, Monthly Notices of the Royal Astronomical Society, 474, 1099 Dello Russo, N., Kawakita, H., Vervack, R. J. J., & Wea...
work page 2021
-
[19]
Magee-Sauer, K., & Griep, D. M. 2003, Evidence for a dominant native source of carbon monoxide in Comet C/1996 B2 (Hyakutake), Journal of Geophysical Research: Planets, 108
work page 2003
-
[20]
Magee-Sauer, K. 2001, Carbon monoxide production and excitation in Comet C/1995 O1 (Hale-Bopp): Isolation of native and distributed CO sources, Icarus, 153, 361
work page 2001
-
[21]
DiSanti, M. A., Mumma, M. J., Russo, N. D., et al. 1999, Identification of two sources of carbon monoxide in Comet Hale–Bopp, Nature, 399, 662
work page 1999
-
[22]
DiSanti, M. A., Bonev, B. P., Gibb, E. L., et al. 2016, En Route to Destruction: The Evolution in Composition of Ices in Comet D/2012 S1 (ISON) between 1.2 and 0.34 AU from the Sun as Revealed at Infrared Wavelengths, ApJ, 820, 20
work page 2016
-
[23]
Dubernet, M., & Quintas-S´ anchez, E. 2019, First quantum study of the rotational excitation of HCN by para-H2O: Convergence of quantum results, influence of the potential energy surface, and approximate rate coefficients of interest for cometary atmospheres., Molecular Astrophysics, 16, 100046, doi: https://doi.org/10.1016/j.molap.2019.100046
-
[24]
1988, in Exploration of Halley’s Comet (Springer), 481–484
Eberhardt, P., Krankowsky, D., Schulte, W., et al. 1988, in Exploration of Halley’s Comet (Springer), 481–484
work page 1988
-
[25]
Faggi, S., Lippi, M., Mumma, M. J., & Villanueva, G. L. 2023, Strongly Depleted Methanol and Hypervolatiles in Comet C/2021 A1 (Leonard): Signatures of Interstellar Chemistry? The Planetary Science Journal, 4, 8, doi: 10.3847/PSJ/aca64c
-
[26]
Faure, A., & Josselin, E. 2010, VizieR Online Data Catalog: Collisional excitation of water in warm media (Faure+, 2008), VizieR Online Data Catalog, J/A+A/492/257
work page 2010
-
[27]
Feldman, P. D., Cochran, A. L., & Combi, M. R. 2004, in Comets II, ed. H. U. Keller & H. A. Weaver (University of Arizona Press), 425
work page 2004
-
[28]
Fray, N., B´ enilan, Y., Biver, N., et al. 2006a, Heliocentric evolution of the degradation of polyoxymethylene: Application to the origin of the formaldehyde (H2CO) extended source in Comet C/1995 O1 (Hale–Bopp), Icarus, 184, 239
work page 1995
-
[29]
Fray, N., B´ enilan, Y., Biver, N., et al. 2006b, Heliocentric evolution of the degradation of polyoxymethylene: Application to the origin of the formaldehyde (H2CO) extended source in Comet C/1995 O1 (Hale-Bopp), Icarus, 184, 239
work page 1995
-
[30]
F., Tsiganis, K., & Morbidelli, A
Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, Origin of the cataclysmic Late Heavy Bombardment period of the terrestrial planets, Nature, 435, 466
work page 2005
-
[31]
1957, Distribution d’intensit´ e dans la tˆ ete d’une com` ete, BSRSL, 43, 740
Haser, L. 1957, Distribution d’intensit´ e dans la tˆ ete d’une com` ete, BSRSL, 43, 740
work page 1957
-
[32]
Huebner, W. F., & Mukherjee, J. 2015, Photoionization and photodissociation rates in solar and blackbody radiation fields, P&SS, 106, 11, doi: 10.1016/j.pss.2014.11.022
-
[33]
Irvine, W., Bergman, P., Lowe, T., et al. 2004, Hcn and Hnc in Comets C/2000 Wm1 (Linear) and C/2002 C1 (Ikeya-Zhang), Origins of life and evolution of the biosphere : the journal of the International Society for the Study of the Origin of Life, 33, 609, doi: 10.1023/A:1025791101127
-
[34]
Irvine, W. M., Bergin, E. A., Dickens, J. E., et al. 1998, Chemical processing in the coma as the source of cometary HNC, Nature, 393, 547
work page 1998
-
[35]
2022, Outburst of comet C/2021 A1 (Leonard) on January 7, 2021, The Astronomer’s Telegram, 15189, 1
Jehin, E., Moulane, Y., Manfroid, J., et al. 2022, Outburst of comet C/2021 A1 (Leonard) on January 7, 2021, The Astronomer’s Telegram, 15189, 1
work page 2022
-
[36]
Jewitt, D., Kim, Y., Mattiazzo, M., et al. 2023, Disintegration of Long-period Comet C/2021 A1 (Leonard), The Astronomical Journal, 165, 122, doi: 10.3847/1538-3881/acb53b
-
[37]
F., Morbidelli, A., Tsiganis, K., Nesvorny, D., & Gomes, R
Levison, H. F., Morbidelli, A., Tsiganis, K., Nesvorny, D., & Gomes, R. 2011, Late Orbital Instabilities in the Outer Planets Induced by Interaction with a Self-gravitating Planetesimal Disk, AJ, 142, 152
work page 2011
-
[38]
Lique, F., Spielfiedel, A., & Cernicharo, J. 2006, Rotational excitation of carbon monosulfide by collisions with helium, A&A, 451, 1125, doi: 10.1051/0004-6361:20054363
-
[39]
C., Bockel´ ee-Morvan, D., Boissier, J., et al
Lis, D. C., Bockel´ ee-Morvan, D., Boissier, J., et al. 2008, Hydrogen Isocyanide in Comet 73P/Schwassmann-Wachmann (Fragment B), ApJ, 675, 931
work page 2008
-
[40]
Loomis, R. A., ¨Oberg, K. I., Andrews, S. M., et al. 2018, Detecting Weak Spectral Lines in Interferometric Data through Matched Filtering, The Astronomical Journal, 155, 182, doi: 10.3847/1538-3881/aab604
-
[41]
Raymond, S. N., & Zolensky, M. 2024, Comets III (University of Arizona Press)
work page 2024
-
[42]
Meier, R., Eberhardt, P., Krankowsky, D., & Hodges, R. R. 1993, The extended formaldehyde source in comet P/Halley, A&A, 277, 677 21
work page 1993
-
[43]
Milam, S. N., Remijan, A. J., Womack, M., et al. 2006, Formaldehyde in Comets C/1995 O1 (Hale-Bopp), C/2002 T7 (LINEAR), and C/2001 Q4 (NEAT): Investigating the Cometary Origin of H2CO, The Astrophysical Journal, 649, 1169, doi: 10.1086/506501
-
[44]
Morbidelli, A., Levison, H. F., Tsiganis, K., & Gomes, R. 2005, Chaotic capture of Jupiter’s Trojan asteroids in the early Solar System, Nature, 435, 462
work page 2005
-
[45]
Mumma, M. J., & Charnley, S. B. 2011, The Chemical Composition of Comets - Emerging Taxonomies and Natal Heritage, ARA&A, 49, 471 Nesvorn´ y, D., Vokrouhlick´ y, D., Dones, L., et al. 2017, Origin and Evolution of Short-period Comets, ApJ, 845, 27, doi: 10.3847/1538-4357/aa7cf6
- [46]
-
[47]
2005, in SF2A-2005: Semaine de l’Astrophysique Francaise, ed
Pety, J. 2005, in SF2A-2005: Semaine de l’Astrophysique Francaise, ed. F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 721
work page 2005
-
[48]
Rabli, D., & Flower, D. R. 2010, The rotational excitation of methanol by molecular hydrogen, MNRAS, 406, 95, doi: 10.1111/j.1365-2966.2010.16671.x
-
[49]
Rodgers, S. D., & Charnley, S. B. 2005, Suprathermal chemical reactions driven by fast hydrogen atoms in cometary comae, MNRAS, 356, 1542, doi: 10.1111/j.1365-2966.2004.08606.x
-
[50]
Roth, N. X., Milam, S. N., Cordiner, M. A., et al. 2021a, Leveraging the ALMA Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide, The Astrophysical Journal, 921, 14, doi: 10.3847/1538-4357/ac0441
-
[51]
Black, J. H. 2005, An atomic and molecular database for analysis of submillimetre line observations, A&A, 432, 369, doi: 10.1051/0004-6361:20041729
-
[52]
2013, GILDAS: Grenoble Image and Line Data Analysis Software, Astrophysics Source Code Library, ascl
Team, G. 2013, GILDAS: Grenoble Image and Line Data Analysis Software, Astrophysics Source Code Library, ascl
work page 2013
-
[53]
Wiesenfeld, L., & Faure, A. 2013, Rotational quenching of H2CO by molecular hydrogen: cross-sections, rates and pressure broadening, MNRAS, 432, 2573, doi: 10.1093/mnras/stt616
-
[54]
Willacy, K., Turner, N., Bonev, B., et al. 2022, Comets in Context: Comparing Comet Compositions with Protosolar Nebula Models, The Astrophysical Journal, 931, 164, doi: 10.3847/1538-4357/ac67e3
-
[55]
Farnham, T. L. 2021, Preview of Comet C/2021 A1 (Leonard) and Its Encounter with Venus, The Astronomical Journal, 162, 194, doi: 10.3847/1538-3881/ac19ba
-
[56]
Zolensky, M. E., Zega, T. J., Yano, H., et al. 2006, Mineralogy and Petrology of Comet 81P/Wild 2 Nucleus
work page 2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.