A&A 430, L37-L40 (2005)
DOI: 10.1051/0004-6361:200400127
E. Lellouch 1 - R. Moreno 1 - G. Paubert 2
1 - Laboratoire d'Études Spatiales et d'Instrumentation en Astrophysique (LESIA),
Observatoire de Paris, 92195 Meudon, France
2 -
Institut de Radio-Astronomie Millimétrique, 18080 Granada, Spain
Received 5 November 2004 / Accepted 14 December 2004
Abstract
Heterodyne observations of Neptune have provided a measurement
of the CO(2-1) line profile with a total bandpass of almost
8 GHz and a resolution of 4 MHz. The lineshape indicates that the CO mole
fraction in Neptune's atmosphere is not uniform, but increases by a factor of
2 from the troposphere/lower stratosphere (0.5 ppm at p>20 mbar)
to the upper stratosphere (1 ppm at p<20 mbar). This indicates the existence
of both external and internal sources of CO. The equivalent flux associated
with the external source is
cm-2 s-1. We
propose that the stratospheric CO results from a large (2 km) cometary
impact that occurred
200 years ago, although there remains problems
with this hypothesis.
Key words: planets and satellites: Neptune - radio-lines: solar system
The case for Neptune's CO being of internal origin is primarily made from
the claim (Marten et al. 1993; Guilloteau et al. 1993; Naylor et al.
1994) that the CO stratospheric mixing ratio also reflects its
abundance in Neptune's troposphere. However, for this, the first two
papers used
broad-band measurements in the vicinity of the CO(3-2) and (1-0) lines.
Naylor et al. (1994) reported a clear detection of the CO(3-2)
line in absorption from FTS/JCMT observations, but the quality of their
spectra
was limited by channel fringing.
From similar measurements in the CO(2-1) line, Encrenaz et al. (1996) rather inferred
a 3-
upper limit of 1 ppm for CO in Neptune's troposphere. Thus none of the earlier observations
have allowed a satisfactory determination of the CO line profiles; in addition,
the presence of an absorption feature does not per se
imply the presence of CO in Neptune's troposphere.
From the theoretical standpoint, models that attempt to explain a CO abundance of 1 ppm in the bulk of Neptune's atmosphere have to invoke an oxygen enhancement in Neptune's interior O/H = 440 times solar (Lodders & Fegley 1994), an extraordinarily high value that implies that water constitutes 60% of the gas phase in Neptune's interior.
The last decade has brought a wealth of observations that have shed new light
on the origin of oxygen compounds in Outer Planets. First, the Shoemaker-Levy
9 collision with Jupiter has revealed that cometary impacts may supply large
amounts of CO and H2O to planetary atmospheres (see e.g. review
in Lellouch 1995). The discovery by ISO of H2O and CO2
in the atmospheres of the Giant Planets and Titan (Feuchtgruber et al.
1997; Coustenis et al. 1998), has provided a second proof for an
external supply of oxygen to the Outer Planets in the form or water and probably
additional species (CO, CO2, CH3OH..., see Moses et al. 2000), and the likely role
of interplanetary micrometeorites, planetary environments (satellites, rings)
and cometary impacts has been demonstrated (Lellouch et al. 2002). In Jupiter's
case, the existence of two distinct sources of carbon monoxide (besides the
recent SL9 event), resulting respectively from internal transport and (probably)
episodic cometary deposition,
has been clearly established (Bézard et al. 2002). Finally,
CO was recently detected in Uranus with a
mixing ratio of
at 0.1-1 bar, and an external origin
was tentatively favored (Encrenaz et al. 2004).
In this context, it appears necessary to reassess the origin of carbon monoxide on Neptune. Here, we present new observations of the CO(2-1) line, that provide a better characterization of its absorption component and from which the vertical distribution of CO is inferred.
Table 1: Observational details.
![]() |
Figure 1: An overview of all observations. Data are plotted directly in antenna temperatures, i.e. are not corrected for forward and main beam efficiencies, or Neptune's filling factor. Black: Aug. 15, 2003. Blue (dashed line): Sep. 1, 2004. Red: May 12, 2004. Green (dashed line): May 13, 2004. The absorption near 231 280 MHz is due to telluric ozone. |
Open with DEXTER |
Figure 1 shows all observations, in antenna temperature scale.
Due to calibration uncertainties, pointing uncertainties,
and residual sky fluctuations, the individual pieces do not, in general, line up "naturally''.
A notable exception is provided by the August 15, observations, which
cover the 226.572-231.432 GHz
range, and clearly show the CO absorption without any need for piecewise adjustment. The observations
of August 15, May 12 and May 13 encompass the CO emission core. After calibration, the measured
antenna temperature minimum at the foot of the emission (230.470 GHz) is 56.1, 47.5 and 64.4 K for the
three observations respectively. Radiative transfer models
predict a value of 64 K at this frequency for a
uniform
mixing
ratio. Allowing for a 15% calibration uncertainty, the May 12 observations thus show a
significant flux deficit compared to expectations. We therefore treat these observations with caution,
but given that they detect the CO emission with appropriate line-to-continuum
contrast, we do not discard them. The procedure to restore the full shape of the CO line was as
follows. All individual pieces were first averaged
according to their central frequency.
Then, the two pieces corresponding to setup 1 were rescaled individually by a constant factor so that their
common part (230.412-230.664 GHz) matches the averaged value in this range (measured 6 times).
Pieces corresponding to other setups were then rescaled sequentially, to ensure proper overlap with adjacent
bands. The resulting line shape is shown in Figs. 2 and 3, in which the precise vertical scale is unimportant.
Our procedure thus applies multiplicative factors to the various pieces of the spectrum, rather than
shifting them vertically, as would be appropriate if (additive) sky fluctuations were the dominant cause of
flux dispersion. Justification is provided by the fact that the constrast
of the CO emission
scales with the local continuum rather than being constant
from one day to another. We note, anyway, that given the restricted frequency coverage of the individual pieces,
applying vertical shifts would produce very small changes in the CO profile. The
final CO line profile shows asymmetry in the far-wings (>2 GHz
from line center), and given the caveats above on the May 12 observations, we regard the high-frequency
wing as the most reliable.
![]() |
Figure 2: Models of the total CO line with uniform CO mole fraction. Data: solid black lines. Models ( from top to bottom): red: CO = 1.0 ppm; green: 0.8 ppm; dark blue: 0.5 ppm; light blue: 0.3 ppm. |
Open with DEXTER |
Figure 2 shows a comparison of the data with models in which CO is assumed to be well mixed
throughout Neptune's atmosphere, with
mole fractions ranging from 0.3 to 1.0 ppm.
A value of 0.8 ppm matches the central emission,
consistent with the initial findings of Marten et al. (1993) and Rosenqvist et al.
(1992). However, while abundances of 0.8-1 ppm also provide a satisfactory match of the absorption
up to
1.5 GHz from line center, they do not allow to fit the more distant wings, which
are more flat, and rather suggest a
mole fraction of
ppm.
The same conclusion is reached even if only the best quality August
15 observations, which extend up to 2.4 GHz from line center,
are considered. We tested its robustness to temperature uncertainties
by using an alternate thermal profile, based on the results of Bézard et al.
(1999) at p>10 mbar, Marten et al. (2005) at p< 1 mbar, and interpolation
between. In the troposphere, this
profile is warmer than ours by
2 K at 2 bar and
4 K at 0.1 bar.
With this profile, the far wings and central emission are fit respectively with
CO abundances of 0.65 and 0.9 ppm, however, the close wings (0.2-0.8 GHz from line
center) are significantly too broad. We conclude that the data
indicate a non-uniformity of the CO mixing ratio in Neptune's stratosphere,
and proceed to test vertically-varying models.
The presence of CO in Neptune's stratosphere is directly demonstrated by the existence of the
emission core. In contrast, the absorption feature does not, in itself, prove that CO is
present in the troposphere, since the absorption could conceivably be formed in the lower
stratosphere (4 mbar < p < 100 mbar) where temperatures are colder than
the tropospheric
continuum near 230 GHz (93 K). However, we found models in which CO is restricted
to pressures less than 100 mbar produce too narrow an absorption (Fig. 3). We then tested
"two-level'' models, in which CO was characterized by a deep value (q1) below a given
pressure ps,
and a high altitude value (q2) at pressures less than ps. The best fit to the grand
average line was achieved for q1 = 0.5 ppm, q2 = 1.0 ppm, and ps = 20 mbar (Fig. 3). The
accuracy on these mole fractions is about 15%, and the ps level is determined to within
a factor of 2.
![]() |
Figure 3:
Models with altitude-varying CO mixing profiles.
Light blue (label S): CO restricted to p< 100 mbar with 1.0 ppm
mole fraction. Green (label 2L): two-level model, with CO = 1.0 ppm
at p< 20 mbar and CO = 0.5 ppm at p> 20 mbar.
Red (label Pa): physical model, with CO = 0.5 ppm in deep troposphere,
an external flux of
![]() |
Open with DEXTER |
A higher mixing ratio in Neptune's middle and upper stratosphere compared to the troposphere/lower
stratosphere value implies the existence of an external source of CO along with the internal
source that is responsible for the deep abundance. The excess of column
density that must be maintained
by this external source is
cm-2,
to within a factor of 2.
Physical profiles of CO were generated by solving the vertical transport equation,
accounting for eddy and molecular transport. Because CO is chemically stable, its downward
flux
is constant and the CO mole fraction q(z) is given by
(omitting here molecular
transport for simplicity). CO profiles were generated for a series of values of
,
with the boundary condition q(
) =
.
Two K(z) profiles were used,
namely the "A'' and "B'' profiles of the Romani et al. (1993)
photochemical model. These profiles have minimum values of 600 and 2000 cm2 s-1, respectively, at the tropopause. Below this level, we
assumed that K(z) is proportional to n(z). These physical models
match the data for a CO flux of (
)
108 cm-2 s-1.
However, CO profiles resulting from the eddy K "A"
profile tend to produce somewhat too narrow an emission, and profile "B'', characterized
by very rapid mixing down to 0.5 mbar, allows a slightly better match to the data. This is
consistent with the findings of Romani et al. (1993) who also favored
the eddy K profile "B'' on the basis of the C2H2 and C2H6 abundance.
We find a deep CO abundance (1.3-2) times lower than the previous estimates of Rosenqvist et al. (1992) and Marten et al. (1993). In the framework of the Lodders & Fegley (1994) thermochemical equilibrium model, this, however, does not alleviate much the need to invoke a huge oxygen abundance, since the CO mixing ratio at observable levels is an extremely sensitive function of the oxygen abundance, increasing typically by 5 orders of magnitude for a factor-of-10 increase in the oxygen enrichment. Bézard et al. (2002) note, however, that the chemical scheme used by the authors is kinetically too ambitious, which may require a reassessment of their work.
Our primary new result is that of an additional, external source of CO,
with a rate of
cm-2 s-1. This is 10-500 times
the magnitude of the external flux of H2O into Neptune, as estimated
by Feuchtgruber et al. (1997) for the Romani et al. (1993) "B''
and "A''
models respectively.
The situation is similar to Jupiter, where - excluding the recent
input due to comet Shoemaker Levy 9 - the ratio of the CO
((1.5-10)
106 cm-2 s-1) to H2O (<8
104 cm-2 s-1) deposition rates is larger than 20 and may be as large as 250 (Lellouch et al. 2002; Bézard et al. 2002). Large CO/H2O
ratios seem to be inconsistent with an oxygen supply by interplanetary grains.
At Jupiter, Bézard
et al. (2002) favored the case for a stratospheric CO resulting from
the downward diffusion of material delivered by the impact of
(sub)kilometer-size comets at sub-millibar levels. We suggest that the same
mechanism takes place on Neptune. The timescale for diffusion
down to
20 mbar is
/
.
Taking H = 30 km and
cm2 s-1, this gives
years.
Integrated over this timescale and the surface of Neptune, the measured
flux gives an input of
37 CO molecules, or
15 CO grams.
This is typically 3 times the SL9 delivery of CO to Jupiter, and may
be produced by a 2 km diameter comet with density 0.5 and 50% CO yield
at impact.
We thus tentatively suggest that the stratospheric CO we currently observe
was produced by a 2 km-comet impact 200 years ago. We note that smaller and more
recent impacts could provide the required equivalent flux,
but not the presence of external CO down to the 20-mbar level. Additional aspects
supporting this scenario are: (i) the CO resulting profile after 200 years would mimick
the two-level profile of Fig. 3; (ii) the external CO/HCN ratio in Neptune's
stratosphere above 20 mbar is about 500, but the condensation of HCN near 3 mbar removes 85% of its column; the corrected ratio
is
75, similar to what was measured in the SL9 debris,
suggesting that the hypothesized event may have produced HCN as well.
Possible problems with the comet scenario, however, are that (i)
CS was detected after the SL9 collision but not
in Neptune (Moreno 1998), and (ii) collisions of 2 km-comets with Neptune are
expected to occur only once every 8000 years according to Bottke
et al. (2002), though, based on cratering rates on the Galilean satellites,
these estimates might be pessimistic by a factor of 10 (Zahnle et al. 1998).