---
title: Comet Evolutionary Diagram (CED)
url: https://www.emergentmind.com/topics/comet-evolutionary-diagram-ced
type: topic
---

# Comet Evolutionary Diagram (CED)

Searching arXiv for the specified papers to ground the encyclopedia entry in the cited literature.
Search query: 2005.09698 OR 2604.09941 OR 2510.26549 OR 2210.02741 OR 1305.2621 OR 1306.5010
The Comet Evolutionary Diagram (CED) denotes a family of diagrammatic frameworks used to represent comet evolution, but it is not a single standardized construct across the literature. In one major usage, the CED is an empirical, comparative plane whose axes are a mass-loss-based age proxy and the number of remaining returns, intended to organize Oort-cloud comets, Jupiter-family comets, asteroidal-belt comets, disintegrators, and suffocated objects. In another usage, developed for the CODE catalogue of long-period comets, the CED is a five-node dynamical track linking a comet’s previous perihelion, original barycentric orbit, present osculating orbit, future barycentric orbit, and next perihelion or escape. More recent work also uses the term for process-resolved or stage-wise schemata that encode volatile drivers, reservoir depths, and seasonal or orbital evolution [1306.5010] [2005.09698] [2510.26549].

## 1. Terminological scope and principal variants

In the cited literature, the term “Comet Evolutionary Diagram” refers to several technically distinct representations rather than to one universally fixed diagram. The main variants can be summarized as follows.

| Variant | Coordinates or nodes | Primary purpose |
|---|---|---|
| Empirical RR–ML-Age / RR–WB-AGE diagram | Horizontal: ML-Age or WB-AGE; vertical: RR | Comparative evolutionary classification |
| CODE five-stage dynamical diagram | Previous → original → osculating → future → next | Orbital-energy and perihelion-distance evolution |
| Model-derived activity/reservoir schema | $r$, $Q_i(r)$, depth, epoch, dominant driver | Thermal-compositional interpretation |
| Stage-wise mission synthesis | Orbital stage A–D with materials and observables | Seasonal and morphological evolution |

The empirical RR–ML-Age form is described as a two-parameter diagnostic that places objects according to how much volatile mass they lose per apparition and how many returns they have left before either sublimating away or being choked by a dust mantle. In that framework, the CED is explicitly compared to the Hertzsprung–Russell diagram in the sense that it organizes comets into physically meaningful families and regimes [2604.09941] [1306.5010].

The CODE form is different in both ontology and data product. It is a compact visualization of a long-period comet’s five standardized dynamical snapshots: a previous perihelion, the inbound original barycentric orbit at 250 au, the observed osculating orbit near the current perihelion, the outbound future barycentric orbit at 250 au, and the next perihelion or escape. Its emphasis is not comparative placement on a population plane, but the evolution of orbital energy and perihelion distance under non-gravitational outgassing, planetary scattering, and Galactic and stellar perturbations [2005.09698].

A third usage appears in the coupled orbital–thermal model of Oort Cloud comets. That paper does not present a named CED, but it explicitly states that its outputs provide all the ingredients to construct one, organized by activity drivers versus heliocentric distance and epoch, volatile retention and depletion by depth and size, and production rates with detectability thresholds and active-area adjustments [2510.26549]. A related stage-wise synthesis for 67P/Churyumov–Gerasimenko organizes composition, morphology, frost, erosion, and seasonal forcing into a structured CED anchored by Rosetta results [2210.02741].

## 2. The empirical RR–age framework

The empirical CED is built on two scalar quantities: an age proxy derived from mass loss, and the number of remaining returns implied by the current erosion rate. In the mass-loss formulation, the horizontal axis is Mass-Loss Age, ML-Age, defined as
$$
\mathrm{ML\text{-}Age} = \frac{k}{M_{\mathrm{loss,tot}}},
$$
with $k = 3.58 \times 10^{11}\ \mathrm{kg\cdot cy}$ and $M_{\mathrm{loss,tot}}$ the integrated mass lost over the active arc. In the water-budget formulation, the horizontal axis is Water-Budget Age,
$$
\mathrm{WB\text{-}AGE}\,[\mathrm{cy}] = \frac{3.58 \times 10^{11}}{\mathrm{WB}\,[\mathrm{kg}]},
$$
with the normalization chosen so that comet 28P/Neujmin 1 corresponds to $100$ comet years. Larger ML-Age or WB-AGE means smaller mass loss per orbit and therefore a more evolved object in this comparative sense [2604.09941] [1305.2621] [1306.5010].

The vertical axis is Remaining Returns, RR. Two equivalent forms are used:
$$
\mathrm{RR} = \frac{r}{\Delta r},
$$
and
$$
\mathrm{RR} = \frac{3\,M_\mathrm{nuc}}{\Delta M}.
$$
Here $r$ is nucleus radius, $\Delta r$ is the shell thickness removed in one active passage, $M_\mathrm{nuc}$ is nucleus mass, and $\Delta M$ is the total mass loss per return. The shell-thickness relation is
$$
\Delta r = \frac{\Delta M}{4\pi \rho r^2},
$$
or, in the water-budget formulation with dust-to-gas ratio $\delta$,
$$
\Delta r = \frac{(1+\delta)\,\mathrm{WB}}{4\pi r^2 \rho}.
$$
RR therefore estimates how many comparable future returns remain before the nucleus either reaches the disintegration limit, $\mathrm{RR}\to 1$, or evolves toward a suffocated state in which $\Delta r \to 0$ and $\mathrm{RR}\to\infty$ [2604.09941] [1305.2621] [1306.5010].

In this framework, the CED is plotted on logarithmic axes spanning many orders of magnitude. The left side contains young Oort-cloud comets with large mass budgets and small ML-Age. The middle contains Jupiter-family comets. The upper right contains old, dust-mantled, or “graveyard” objects, including asteroidal-belt comets and active asteroids. The lower boundary at $\mathrm{RR}=1$ marks disintegrators such as C/2002 O4 Hönig and C/2012 S1 ISON. One paper further defines a “graveyard of comets” by $\mathrm{WB\text{-}AGE}>1000\ \mathrm{cy}$, and another places main-belt or asteroidal-belt comets in the upper-right graveyard at large age and large RR [1305.2621] [1306.5010].

The empirical CED is closely tied to Secular Light Curve methodology. Photometry is reduced to absolute magnitude, often using
$$
m_V(1,1,0) = m_V(\Delta, R, \alpha) - 5\log_{10}(\Delta R) - \beta \alpha,
$$
with the important caveat that for optically thick comae the nucleus is not seen and $\beta$ is effectively undefined, so an envelope method is used rather than a phase correction to the coma-dominated flux. The envelope of the data is taken as the “truth,” because practical observing limitations bias magnitudes downward. Production-rate integration over the active interval then yields the mass budget. For dust, one widely used empirical conversion is Af$\rho = 1$ cm corresponding to $1\ \mathrm{kg\ s^{-1}}$, while the paper on ISON writes the equivalent scaling as $1000\ \mathrm{cm} = 1000\ \mathrm{kg}$ [2604.09941] [1306.5010].

Several interpretive borders are integral to this version of the CED. A heuristic separation is drawn between sublimating-away comets, which lose ice with little residue and evolve downward toward $\mathrm{RR}=1$, and suffocating comets, where dust fallback quenches activity and drives migration up and to the right. The “suffocation–sublimation border” is estimated in one study at $\mathrm{RR(SSB)} = (6\pm5)\times 10^4$, and a lower-right “desert” is predicted because very old, nearly exhausted comets should become undetectably small [1306.5010] [2604.09941].

## 3. The CODE dynamical CED for long-period comets

The CODE catalogue defines a different CED centered on orbital dynamics. It records each long-period comet in five stages that cover three successive perihelion passages and form the five nodes of a polyline. Snapshot 1 is the observed osculating orbit near the current perihelion and is heliocentric. Snapshot 2 is the original barycentric orbit at 250 au inbound. Snapshot 3 is the future barycentric orbit at 250 au outbound. Snapshot 4 is the previous perihelion, obtained by long-term backward integration under external perturbations. Snapshot 5 is the next perihelion or the escape state at 120,000 au [2005.09698].

The original and future barycentric states isolate the net effect of planetary perturbations during the observed perihelion passage. The previous and next perihelia connect that short-timescale planetary scattering to long-timescale Galactic and stellar forcing. In practice, the nominal orbit is cloned into a swarm of 5001 virtual comets using a Monte Carlo method, and each swarm is integrated to 250 au for the original and future states. For the previous and next states, the integrations include the Galactic potential and a curated list of 643 stars or stellar systems that may approach within 4 pc over approximately $\pm 10$ Myr. Integrations stop at the previous or next perihelion for returning virtual comets, or at an escape limit of 120,000 au for escaping ones [2005.09698].

The most common vertical axis in this CED is inverse semimajor axis, $1/a$, plotted in units of $10^{-6}\ \mathrm{AU}^{-1}$, because it is a proxy for specific orbital energy:
$$
E = -\frac{\mu}{2a}.
$$
Higher $1/a$ corresponds to more tightly bound orbits; for hyperbolas, $a<0$ and $1/a<0$. A complementary axis is perihelion distance,
$$
q = a(1-e),
$$
used either as the horizontal axis or as an annotation at the five nodes to show how energy and perihelion evolve together [2005.09698].

This formulation explicitly incorporates non-gravitational acceleration. CODE uses the Marsden–Sekanina–Yeomans formalism,
$$
\mathbf{a}_{\mathrm{NG}} = g(r)\,[A_1\hat{\mathbf{r}} + A_2\hat{\boldsymbol{\theta}} + A_3\hat{\mathbf{w}}],
$$
with
$$
F_i = A_i\,g(r), \qquad i=1,2,3,
$$
and
$$
g(r) = \alpha\,(r/r_0)^m\,[1 + (r/r_0)^n]^k,
$$
normalized so that $g(1\ \mathrm{au})=1$. CODE lists standard parameter sets for water-ice sublimation, typically for $q\leq 3.1$ au, and for CO-driven activity, often for $q\geq 3.1$ au. Asymmetric variants are used when warranted by residuals, but symmetric models are preferred unless an asymmetric model significantly improves the fit [2005.09698].

The dynamical classifications attached to this CED also differ from those in the RR–age literature. CODE’s sample centers on Oort-spike comets defined by purely gravitational original semimajor axis $a_{\mathrm{orig}} > 10{,}000$ au, or equivalently
$$
1/a_{\mathrm{orig}} < 100 \times 10^{-6}\ \mathrm{AU}^{-1}.
$$
However, dynamical new versus old is classified primarily by the previous perihelion distance obtained from the virtual-comet swarm: $q_{\mathrm{prev}} > 20$ au implies dynamically new, $q_{\mathrm{prev}} < 10$ au dynamically old, and $10 \leq q_{\mathrm{prev}} \leq 20$ au uncertain. For hyperbolic cases, the five-stage scheme is adapted so that the next snapshot is taken at the 120,000 au escape limit, often with $1/a_{\mathrm{fut}} < 0$ [2005.09698].

Uncertainty treatment is intrinsic to the CODE CED. At snapshots 2 and 3, distributions at 250 au are often approximately Gaussian, so mean and standard deviation are reported. At snapshots 4 and 5, non-Gaussian outcomes are common, so CODE provides counts of returning, escaping, and hyperbolic virtual comets and summarizes distributions by deciles. On the diagram, this permits error bars at the original and future stages and percentile boxes or ranges at the previous and next stages [2005.09698].

## 4. Thermophysical and mission-based extensions

A model-driven extension of the CED concept appears in the coupled orbital–thermal evolution of Oort Cloud comets. There the orbit $r(t)$ is fed continuously into a thermophysical solver over $4.6$ Gyr, with four dynamical phases: early planetary evolution near Neptune with $q\approx 30$ au for about $30$ Myr, outward migration over about $3$ Myr, Oort Cloud residence with $a\approx 10{,}000$ au for about $4.5$ Gyr, and a final return inward with $q\approx 5$ au. The nuclei are porous, spherical dust–ice mixtures with radii of 2, 10, and 50 km, and the initial volatile inventory includes CO and CO$_2$ ices, amorphous water ice with trapped CO and CO$_2$, and dust [2510.26549].

The governing physics includes porous heat conduction, sublimation and refreezing within pores, Knudsen-regime gas transport, crystallization of amorphous ice, and a test of long-lived radiogenic heating. The surface energy balance is written as
$$
(1-A)\,\frac{S_0}{r^2}\,\cos\zeta
=
\epsilon\sigma T_s^4
+
\sum_i L_i J_i
+
k\left.\frac{\partial T}{\partial x}\right|_{x=0},
$$
and the crystallization kinetics are
$$
\lambda(T) = c\exp\!\left(-\frac{b}{T}\right),
\qquad
c = 1.05\times 10^{13}\ \mathrm{s}^{-1},
\qquad
b = 5370\ \mathrm{K}.
$$
The resulting CED content is not a single point but a structured map of dominant activity drivers, volatile reservoir accessibility, production rates, and detectability thresholds [2510.26549].

Several quantitative thresholds emerge from that model. CO sublimation is the dominant large-distance driver and first exceeds detectability at $r\approx 50$ au for $R=10$ km and $r\approx 150$ au for $R=50$ km. CO$_2$ becomes detectable around $r\approx 13$ au in the global-average implementation, and under subsolar or limited-active-area assumptions can be detectable out to about $20$–$25$ au. Amorphous ice begins crystallizing around $r\approx 8$ au, and the release of trapped gases becomes the main source of outgassing around $r\approx 6.5$–$7$ au, with peak total production reaching about $10^{28}\ \mathrm{molecules\ s^{-1}}$ for the 10 km model and about $10^{29}\ \mathrm{molecules\ s^{-1}}$ for the 50 km model. The model also finds that CO ice is completely depleted in the 2 km nucleus but preserved from depths of about 500 m inward in the larger nuclei, while CO$_2$ and amorphous H$_2$O are entirely preserved; radiogenic heating by long-lived isotopes is negligible [2510.26549].

A mission-anchored extension appears in the structured CED for 67P/Churyumov–Gerasimenko. That synthesis uses Rosetta observations to define four stages: aphelion dormancy, inbound activation, perihelion peak activity, and outbound relaxation. The nucleus is described as very dark, red-sloped, and covered by a variable-thickness dust mantle, with geometric albedo $5.9\pm0.3\%$ at 535 nm and $6.2\pm0.2\%$ at 550 nm, porosity around $85$–$87\%$, patchy H$_2$O and transient CO$_2$ exposures, ammoniated salts, refractory organics with aliphatic and aromatic components, and hydroxylated Mg-amorphous silicates contributing to the broad 3.2 $\mu$m absorption [2210.02741].

The stage-wise 67P CED is explicitly tied to heliocentric forcing and energy balance. Frost lines are given as H$_2$O at about 3 AU, CO$_2$ at about 14 AU, CO at about 28 AU, and N$_2$ at about 32 AU. The instantaneous energy balance is written as
$$
(1-A)F_0 r^{-2}\cos\theta = \epsilon \sigma T^4 + L\dot{m} + k\frac{\partial T}{\partial z},
$$
with sublimation flux represented by the Hertz–Knudsen relation. Observationally, the structured diagram encodes dawn frost cycles in Hapi, CO$_2$ ice in Anhur, perihelion outbursts from cliffs and pits, net erosion of about 4 m near perihelion, and post-perihelion fallback that blankets northern terrains and restores redder colors. This suggests a CED usage in which composition, morphology, seasonal illumination asymmetry, and observable spectral or photometric changes are treated as a single evolutionary system rather than as scalar summary indices [2210.02741].

## 5. Representative placements and case studies

The empirical RR–age literature uses the CED to place specific comets within a comparative evolutionary taxonomy. For exocomet 3I/ATLAS, the Secular Light Curve and integrated production rates yield
$$
\mathrm{ML\text{-}Age} \approx 0.16\ \mathrm{cy},
\qquad
\mathrm{RR} \approx 24.
$$
The total mass loss is reported as $2.274\times 10^{12}$ kg, composed of dust $3.9\%$, H$_2$O $2.66\%$, CO$_2$ $93.23\%$, and CO $0.20\%$, making the object strongly CO$_2$-dominated. On the log–log CED it lies in the extreme lower-left domain characteristic of young, icier comets with minimal dust mantling, colocated with the Oort-cloud group. The same work reports a photometric anomaly from $-120$ to $-45$ days before perihelion, interpreted as an eclipse consistent with a possible binary nucleus, and a maximum reduced absolute magnitude of $m_V(1,1,\alpha)=6.8\pm0.1$ about 15 days after perihelion [2604.09941].

The same empirical framework places C/2012 S1 ISON at the disintegration limit. Its Secular Light Curve showed a Slope Discontinuity Event at about $-4.1$ AU, followed by a near-standstill lasting more than 132 days. The paper reports $\mathrm{WB}=3.94\times10^{10}$ kg, dust budget $2.35\times10^{11}$ kg, CO budget $2.5\times10^9$ kg, and total $\mathrm{ML\text{-}Budget}=3.05\times10^{11}$ kg, implying
$$
\mathrm{ML\text{-}AGE} = 1.2\ \mathrm{cy},
\qquad
\mathrm{RR}=1.
$$
With $\rho = 530\ \mathrm{kg\ m^{-3}}$, the corresponding diameter is $1030\pm70$ m. The paper treated the SDE+standstill signature as a predictor of disintegration and reports that the comet indeed fragmented at $R\approx -0.67$ AU. By contrast, C/2011 L4 Panstarrs, with an SDE at $-4.97\pm0.03$ AU and post-SDE brightening to $m(1,q)=-1.2\pm0.2$, is placed in the Oort-cloud region on the left side of the CED [1306.5010].

Asteroidal-belt or asteroidal-belt-comet placements occupy the opposite corner. For 133P/Elst–Pizarro, the paper gives
$$
\mathrm{WB} \approx 1.81\times10^8\ \mathrm{kg},
\qquad
\mathrm{WB\text{-}AGE} \approx 2.0\times10^3\ \mathrm{cy},
$$
with $\mathrm{RR}\approx 2.2\times10^5$ for $\delta=1$ and $\approx 4.0\times10^5$ for $\delta=0.1$. For 107P/Wilson–Harrington, $\mathrm{WB\text{-}AGE}\approx1.8\times10^4$ cy and $\mathrm{RR}\approx 7\times10^5$ to $1.3\times10^6$. For 300163 = 2006 VW139, $\mathrm{RR}(\delta=1)\approx 9.7\times10^6$. These values place the objects firmly in the upper-right graveyard, where activity is weak, mantled, and episodically rejuvenated by decreases in perihelion distance. The same paper gives dust-mantle thicknesses of $2.0\pm0.5$ m for 107P, $4.7\pm1.2$ m for 133P, and $1.9\pm0.5$ m for a sample of nine comets, derived from post-perihelion activation lags and a thermal-wave model [1305.2621].

The CODE dynamical CED yields a different type of case study. For C/2002 T7 LINEAR, the preferred NG solution “d6” gives observed osculating perihelion $q = 0.615$ au and
$$
1/a_{\mathrm{orig}} = 21.65 \pm 0.86 \times 10^{-6}\ \mathrm{AU}^{-1}
$$
at 250 au inbound. Across water-ice NG variants, $1/a_{\mathrm{orig}}$ ranges from $17.75$ to $28.93 \times 10^{-6}\ \mathrm{AU}^{-1}$; a CO-driven model yields $53.21 \pm 3.63 \times 10^{-6}\ \mathrm{AU}^{-1}$ but is judged implausible for this small-$q$ comet. The previous perihelion is at approximately 150 au, making the comet dynamically new by the $q_{\mathrm{prev}}$ criterion, while the future orbit is hyperbolic, with the next snapshot recorded at the 120,000 au escape limit. In the CODE CED, the track therefore moves from small positive $1/a$ at the original stage to negative $1/a$ at the future and next stages, visually encoding ejection by strong planetary scattering [2005.09698].

## 6. Interpretation, misconceptions, and unresolved issues

A recurrent misconception is to treat the CED as a single standard diagram. The literature does not support that usage. At minimum, the term denotes a comparative RR–age plane, a five-snapshot dynamical trajectory for long-period comets, and model-derived or stage-wise schemata that encode activity drivers and material evolution. These variants are related by their concern with comet evolution, but their variables, timescales, and inferential targets differ [2005.09698] [2604.09941] [2510.26549].

A second misconception is to read ML-Age or WB-AGE as chronological age. In the empirical CED, these are explicitly proxies based on inverse mass loss, normalized in comet years by a calibration constant. They express comparative evolutionary state under current or inferred activity, not time since formation. Likewise, RR is not a deterministic lifetime; it is an estimate under the assumption that the per-apparition mass loss remains roughly constant [2604.09941] [1305.2621].

In the CODE framework, dynamical new versus old is also frequently misunderstood. CODE does not classify long-period comets primarily by a hard threshold in $1/a_{\mathrm{orig}}$. The operative criterion is the previous perihelion distance from the virtual-comet swarm, with $q_{\mathrm{prev}} > 20$ au marking dynamically new, $q_{\mathrm{prev}} < 10$ au dynamically old, and the intermediate range treated as uncertain. The Tisserand parameter with respect to Jupiter can aid interpretation of planetary-scattering signatures during the osculating stage,
$$
T_J = \frac{a_J}{a} + 2\sqrt{\frac{a}{a_J}(1-e^2)}\cos i,
$$
but CODE does not classify long-period comets by $T_J$, and for LPCs it is best used only as a qualitative indicator because it is not conserved across deep planetary encounters [2005.09698].

The limitations of the various CEDs are explicit in the source literature. In RR–age studies, photometric apertures can under-extract coma flux, water datasets may require scale corrections by factors of about 5 to 12, and the Af$\rho\rightarrow\mathrm{kg\ s^{-1}}$ conversion is empirical and order-of-magnitude in nature. Derived quantities depend on adopted density, albedo, size, and production-rate integration, although the authors emphasize that the log–log diagram is “forgiving” because factor-of-two uncertainties move points little [2604.09941] [1306.5010].

In CODE, previous and next perihelion outcomes depend on the completeness and accuracy of the stellar perturber list, on the adopted non-gravitational model, and on whether asymmetric $g(r)$ or data-subset solutions are needed to avoid biased osculating fits. Uncertainties become non-Gaussian at long timescales, which is why deciles rather than means and standard deviations are often reported for the previous and next stages [2005.09698].

The thermophysical extensions add another layer of uncertainty. The coupled orbital–thermal model is one-dimensional and spherically symmetric, anchored to a single representative dynamical history; porosity, pore size, tortuosity, thermal conductivity, dust fraction, and initial volatile abundances all shift depletion depths, refreezing profiles, onset distances, and predicted production rates. The 67P synthesis, although exceptionally rich observationally, remains a stage-wise abstraction of a nucleus whose activity is strongly controlled by local shading, morphology, and seasonal asymmetry [2510.26549] [2210.02741].

Taken together, the literature presents the CED not as one diagram but as a compact language for encoding comet evolution at different levels: comparative population state, dynamical transfer through the planetary system and the Oort cloud, and thermophysical or seasonal reorganization of surface and subsurface volatiles. A plausible implication is that future convergence among these variants will require models that connect mass-loss-based age proxies, orbit-resolved energy changes, and internally resolved volatile evolution within a single formalism.

Source: https://www.emergentmind.com/topics/comet-evolutionary-diagram-ced