Issue |
A&A
Volume 674, June 2023
|
|
---|---|---|
Article Number | L3 | |
Number of page(s) | 8 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202346621 | |
Published online | 02 June 2023 |
Letter to the Editor
Doppler wind measurements in Neptune’s stratosphere with ALMA
1
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
e-mail: oscar.carrion.gonzalez@gmail.com, oscar.carrion@obspm.fr
2
Laboratoire d’Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, Allée Geoffroy Saint-Hilaire, 33615 Pessac, France
3
Laboratoire de Météorologie Dynamique/Institut Pierre-Simon Laplace (LMD/IPSL), Sorbonne Université, CNRS, École Polytechnique, Institut Polytechnique de Paris, École Normale Supérieure (ENS), PSL Research University, Paris, France
Received:
7
April
2023
Accepted:
10
May
2023
Context. Neptune’s tropospheric winds are among the most intense in the Solar System, but the dynamical mechanisms that produce them are still unclear. Measuring wind speeds at different pressure levels may help shed light on the atmospheric dynamics of the planet.
Aims. The goal of this work is to directly measure winds in Neptune’s stratosphere with ALMA Doppler spectroscopy.
Methods. We derived the Doppler lineshift maps of Neptune at the CO(3-2) and HCN(4-3) lines at 345.8 GHz (λ ∼ 0.87 mm) and 354.5 GHz (0.85 mm), respectively. For that purpose, we used spectra obtained with ALMA in 2016 and recorded with a spatial resolution of ∼0.37″ on Neptune’s 2.24″ disk. After subtracting the planet’s solid rotation, we inferred the contribution of zonal winds to the measured Doppler lineshifts at the CO and HCN lines. We developed an MCMC-based retrieval methodology to constrain the latitudinal distribution of wind speeds.
Results. We find that CO(3-2) and HCN(4-3) lines probe the stratosphere of Neptune at pressures of 2−1.8+12 mbar and 0.4−0.3+0.5 mbar, respectively. The zonal winds at these altitudes are less intense than the tropospheric winds based on cloud tracking from Voyager observations. We find equatorial retrograde (westward) winds of −180−60+70 m s−1 for CO, and −190−70+90 m s−1 for HCN. Wind intensity decreases towards mid-latitudes and wind speeds at 40°S are −90−60+50 m s−1 for CO and −40−80+60 m s−1 for HCN. Wind speeds become 0 m s−1 at about 50°S. We find that the circulation reverses to a prograde jet southwards of 60°S. Overall, our direct stratospheric wind measurements match previous estimates from stellar occultation profiles and expectations based on thermal wind equilibrium.
Conclusions. These are the first direct Doppler wind measurements performed on the Icy Giants, opening up a new method for the study and monitoring of their stratospheric dynamics.
Key words: planets and satellites: atmospheres / planets and satellites: gaseous planets / radiative transfer
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
Voyager 2 measurements of Neptune’s atmosphere have revealed some of the most intense zonal winds ever measured in the Solar-System planets (Hammel et al. 1989; Smith et al. 1989). This was done by tracking cloud motions in images of Voyager’s ISS instrument, thereby constraining Neptune’s winds to be −400 m s−1 (retrograde) at the equator, with a prograde jet of about +250 m s−1 at latitude 70°S (Limaye & Sromovsky 1991; Sromovsky et al. 1993). Subsequent cloud-tracking measurements with the Hubble Space Telescope (Sromovsky et al. 2001a,b) and ground-based observatories with adaptive optics (Fitzpatrick et al. 2014; Tollefson & Pater 2018) have confirmed this general wind pattern (Fletcher et al. 2020).
However, cloud-tracking methods generally lack an accurate constraint of the atmospheric altitudes probed by the observations. Although clouds are mostly expected to be located in the upper troposphere (∼100 mbar–1 bar), determining the exact cloud-top pressure levels requires multiple-scattering radiative-transfer computations. These are highly dependent on the assumed optical properties of the atmospheric aerosols and on the assumed vertical distribution of gaseous species and aerosols (e.g. Luszcz-Cook et al. 2016). Furthermore, cloud pressure levels have been found to vary with latitude by about an order of magnitude (Irwin et al. 2011, 2016; Hueso et al. 2017). Even clouds at similar latitudes have been found to be located at pressure levels across such a varied range of pressure levels as 30–300 mbar (de Pater et al. 2014).
Above the cloud level, wind information has been derived primarily from the thermal wind equation that relates the latitudinal temperature gradient to the vertical wind shear. Initially, temperature fields from Voyager/IRIS spectra were used for this (Conrath et al. 1989). Later, additional ground-based thermal observations became available (Fletcher et al. 2014). At tropospheric levels, these data consistently indicate a warm equator, cool mid-latitudes, and a warmer south pole. The thermal wind equation (which cannot be applied at the equator) implies decaying winds with altitude at low latitudes (−30° to 30°) and essentially altitude-independent winds at mid-to-high southern latitudes. Relatively small changes were found in the predicted circulation pattern from 1989 to 2003 (southern summer solstice), except near the south pole, which was found to be anomalously warm (by ∼5 K) in 2003 (Orton et al. 2007).
Additional wind measurements with other techniques were performed for a reduced number of latitudes. French et al. (1998) measured stratospheric winds from three stellar occultations over 1985–1989, constraining the wind-distorted shape of the planet. They inferred that wind speeds at the 0.38 mbar level are 0.6±0.2 times those of Voyager cloud-tracking from Sromovsky et al. (1993). Also, they found tentative evidence of winds at 0.7 μbar further decaying to ∼0.17 times the tropospheric values. This was the first direct evidence for the decay of winds above the tropopause, with an estimated wind shear in agreement with that inferred from Voyager temperature fields at deeper levels.
In this work, we aim to directly measure Neptune’s stratospheric winds based on Doppler lineshift measurements of CO and HCN lines with the Atacama Large Millimeter/submillimeter Array (ALMA). Millimetric and submillimetric measurements have already been used to directly measure the winds in the atmospheres of Venus (e.g., Lellouch et al. 2008, and references therein), Mars (Lellouch et al. 1991; Cavalié et al. 2008; Moreno et al. 2009), Jupiter (Cavalié et al. 2021), Saturn (Benmahi et al. 2022), and Titan (Moreno et al. 2005; Lellouch et al. 2019; Cordiner et al. 2020). No such analysis has been available for the Icy Giants to date. Besides probing the stratospheres of these planets, which is otherwise challenging, an advantage of the technique is that it is insensitive to aerosols, reducing considerably the uncertainties on the probed levels.
2. Observations
We observed Neptune on April 30, 2016, over 20 min on-source, with 41 antennas of the ALMA interferometer in the C36-2/3 hybrid configuration, yielding an angular resolution of about 0.37″. The spectral setup included the CO(3-2) line at 345.7959899 GHz and the HCN(4-3) line at 354.5054773 GHz, as well as the CS(7-6) line at 342.883 GHz (whose detection was reported in Moreno et al. 2017), with a spectral resolution of 1 MHz. The bandpass, amplitude, and phase calibration procedure of the CS visibilities was described in Moreno et al. (2017), and we applied the same calibration procedure to the CO and HCN data using the ALMA/CASA data reduction software.
The resulting calibrated visibilities were then exported into the GILDAS package with the following aims: (i) to apply a self-calibration technique using Neptune’s continuum to improve the image quality (ii) to perform the imaging and deconvolution using the Högbom algorithm (Högbom 1974). We obtained a synthetic elliptical beam (robust weighting 0.5) of 0.39″×0.34″ (polar angle (PA) of −47°) for CO and 0.37″×0.35″ (PA = −82°) for HCN. The planet’s angular diameter was 2.24″. This yields a spatial resolution of about 20° at the equator. The resulting spectral maps are shown Fig. A.1 with signal-to-noise ratio (S/N) at line peak at the limbs of about 150 and 50, respectively, for CO and HCN. The final clean images were built with a sampling of 0.05″ over −1.4″ and +1.4″ (3136 points).
An example of the observed lines is shown in Fig. 1. The detailed analysis of the lineshapes in terms of spatial and vertical distribution of temperature, as well as CO and HCN, is left for future work. For the purpose of measuring winds from line central positions, we performed Gaussian fits: we used three-Gaussian fits for CO (with initial FWHM of 4, 20, and 80 km s−1) and two-Gaussian fits for HCN (with initial FWHM of 3 and 16 km s−1) as shown in Fig. 1. We retained the narrowest component to measure the Doppler lineshift. The high S/N in our maps allowed us to derive the Doppler lineshifts with an averaged velocity accuracy of 25 and 37 m s−1, respectively, for CO and HCN (Fig. 2, right column).
![]() |
Fig. 1. Example of measured spectra (left column) of CO(3-2) and HCN(4-3) lines (black) at the position of the sky western equatorial limb, at an offset from Neptune centre (−1.1″, 0.0″). Individual Gaussian fit components with FWHMs of about 4, 20, 80 km s−1 are shown in dashed blue and their sum in red. Also shown: the beam-convolved solid body velocity (black dotted line) and the fit velocity (red dotted line). The difference of these two velocities allows us to derive the Doppler winds. Right: normalized wind contribution functions (WCF) at the limb for each molecule. |
![]() |
Fig. 2. ALMA Doppler measurements for the CO(3-2) line (top row) and the HCN(4-3) line (bottom row). Left: measured Doppler lineshift. Middle: line-of-sight winds, after subtracting the solid-body rotation. Right: root mean square (rms) of the measured winds. Top-right ellipse in each subplot shows the synthetic beam. Latitudes are shown in 20° steps, with the equator displayed as the thicker line. |
3. Model
3.1. Radiative transfer
We used the same Neptune radiative transfer model described in Moreno et al. (2017), along with their thermal profile and their CO and HCN vertical distributions, to model the spectral lines of CO and HCN, as well as to compute the wind weighting function shown in Fig. 1. Wind weighting functions account for the spectrally-dependent wind information content of each channel within lines (see e.g., Lellouch et al. 2019), and were convolved by the beam. At the limb (where most of the wind information comes from) CO lines probe the mbar level and HCN lines probe the
mbar level.
3.2. Wind retrieval methodology
In order to interpret the measured Doppler lineshifts, we developed a retrieval framework for the wind profiles based on the Markov chain Monte Carlo (MCMC) emcee sampler (Foreman-Mackey et al. 2013). The observed lineshifts correspond to the sum of the line-of-sight Doppler displacement due to the planetary rotation and the winds. The planet rotation was modelled as solid-body rotation at the altitude of the 1 mbar level above the local planet radius. Although somewhat different rotation periods have been proposed (Helled et al. 2010; Karkoschka 2011), a period of 16.11 h (Warwick et al. 1989; Lecacheux et al. 1993) is adopted here to enable comparisons with the Voyager winds.
We parameterized the wind profiles (W), assumed to be purely zonal, as a polynomial function depending on the latitude (ϕ, given in degrees) that takes the form:
and we used emcee to explore the parameter space described by the coefficients an, given in m s−1. The box priors which set the limits of the parameter space were [−500, 500] for a0 and [−1, 1] for the other coefficients. We tested polynomial orders from 0 to 5 for the wind parameterization, as discussed in Sect. 4.
The MCMC sampler tests points of the n-dimensional parameter space of an coefficients. For each test point, we computed the wind profile at latitudes ϕ ∈ [ − 90° ,90° ] and the resulting wind map, with a pixel stepsize of 0.05″. To simulate the line-of-sight Doppler lineshifts at each point of the map, we projected the zonal winds onto the planet geometry, accounting for the sub-observer latitude of 26.2°S. We then weighted the modelled lineshifts by the local intensity of the CO (resp. HCN) line and convolved by the ALMA beam, following a similar procedure to Lellouch et al. (2019). After adding the contribution of the solid body rotation, the test wind map (Wtest) was compared with the one measured by ALMA (WALMA) for the molecule under study (Fig. 2, middle) by means of the χ2 figure of merit:
where rmsALMA is the 1-σ error of the ALMA measurement for the molecule under study (Fig. 2, rightmost column) and Np is the number of pixels in the wind map.
We used 50 chains (or “walkers”) to simultaneously explore the parameter space independently in order to avoid possible χ2 local minima. The 50 walkers ran for up to 5×104 steps, with a convergence criterion to stop the run at a number of steps larger than 50 times the autocorrelation time (τ, see Foreman-Mackey et al. 2013, for details). This convergence criterion ensures that the sampling of the parameter space has been completed and the sampled test points are effectively independent (Goodman & Weare 2010). For the analysis of each retrieval run, we discarded the samples of the first τ steps (“burn-in phase”).
For each retrieval run, we defined the ensemble of good fits as those sampled wind profiles with a χ2 value in the 68.3% (i.e., 1-σ) confidence level. That is as follows, with χ2 verifying:
Here, is the minimum χ2 value among all the samples in the run, and
is the reduced value of
, with N equal to the degrees of freedom of the retrieval (N = Np − n). The factor foversampling = Areabeam/Areapixel accounts for the fact that with a spatial sampling step of 0.05″ and a beam of 0.37″, the measurements are considerably oversampled, and the number of independent measurements is Np/foversampling. The coefficient C denotes the 1-σ confidence region in the n-dimensional parameter space, where n is the number of an coefficients in the polynomial fit. For polynomial orders 0–5, C is 1.0, 2.3, 3.53, 4.72, 5.89, and 7.04, respectively (see Ch. 15.6 in Press et al. 2007).
4. Results
With the method above, we carried out wind retrievals for the CO and HCN measurements. We ran retrievals for wind parameterizations with polynomials of orders between 0 and 5. Previous cloud-tracking studies had generally assumed latitudinally symmetric wind profiles and thereby omitted odd polynomial orders (Sromovsky et al. 1993; Fitzpatrick et al. 2014; Tollefson & Pater 2018). However, the spatial resolution achieved in our data is potentially sensitive to latitudinal wind asymmetries. Therefore, we kept odd polynomial orders in our retrievals.
The retrieval results for CO and HCN with polynomial fits of orders 0–5 are shown in Fig. B.1. The sixth-order polynomial fit to Voyager’s cloud-tracking measurements (Sromovsky et al. 1993) is shown for comparison. For each tested polynomial order, we include in Fig. B.1 the map of residuals between the observed Doppler lineshifts and the best-fitting model (solid rotation plus wind) to assess the variation of fit quality as a function of polynomial degree. For both CO and HCN, we find that the value of χ2/N is practically the same for retrievals of orders 3, 4, and 5. This implies that these polynomials are similarly able to fit the measurements at the latitudes we are sensitive to. Indeed, the three parameterizations retrieve similar values of the zonal winds at latitudes over 20°N to 70°S.
Given Neptune’s sub-observer latitude of 26.2°S and projection effects, the wind information is restricted to latitudes southward of 40°N. In addition, the beam size of 0.37″ prevents the retrieval of detailed information at southern polar latitudes. Angular speed considerations indicate that the wind speed should theoretically be zero at the pole. We find that this constraint is met within the error bars for the fourth- and fifth-order polynomial fits. Therefore, the order n = 4 represents a good compromise between the model’s mathematical complexity and physical realism, and we adopted it as the reference wind parameterization.
Figure 3 shows the retrieved wind profiles for the reference n = 4 parameterization. Table 1 shows the polynomial coefficients for CO and HCN best-fit wind profile, as well as the coefficients that parameterize the envelope of good solutions – namely, those verifying Eq. (3) – using a fourth-order polynomial fit. Figure 3 also shows our retrieved wind measurements at 0°, 20°, and 70°S as a function of pressure. For reference, the plot includes previously reported zonal winds (see discussion in Sect. 5).
![]() |
Fig. 3. Retrieval results. Left: retrieved best-fitting wind profiles for CO (red line) and HCN (blue line) measurements, using a fourth-order polynomial. Red and blue shadowed regions contain the ensemble of good fits according to the χ2 criterion from Eq. (3). The semi-transparent grey rectangle indicates the unobservable northern latitudes. The solid black line shows the sixth-order fit to Voyager’s cloud-tracking winds (Sromovsky et al. 1993). Right: wind variations with altitude at the equator, 20°S and 70°S, both for our measurements and for a set of references in the literature (see Sect. 5 for details). Cloud-tracking winds from Fitzpatrick et al. (2014) and Tollefson & Pater (2018) are not shown at 70°S as they have uncertainties of about 1000 m s−1. Winds from Fletcher et al. (2014) do not correspond to a direct measurement, but to the computed thermal wind equation applied to a reanalysis of the 1989 IRIS/Voyager data (dashed grey lines) and to 2003 Keck data (dashed cyan lines). |
Polynomial coefficients of the best fits to the CO and HCN zonal wind measurements.
At the equator, we retrieved retrograde (westward) zonal winds of m s−1 from CO measurements and
m s−1 from HCN. We find a decrease in the wind intensity towards mid-latitudes. At 20°S, wind speeds are
m s−1 for CO (
m s−1 for HCN) and at 40°S, they are
m s−1 for CO (
m s−1 for HCN). Winds continue along this trend towards southern latitudes, becoming 0 m s−1 at about 50°S. At 70°S, winds are prograde (eastward) although uncertainties are larger (
m s−1 for CO, and
m s−1 for HCN). Further south than 70°S, the winds remain unconstrained due to the limited spatial resolution. Wind speeds in the observable northern-hemisphere latitudes are compatible, within errors, with a symmetric wind profile. At 20°N, zonal winds are
m s−1 for CO and
m s−1 for HCN. Additional measurements will be needed to confirm this behaviour, eventually when regions further north become observable.
5. Discussion
Our results offer a new method for probing Neptune’s stratospheric winds and wind shear, while also yielding latitudinal information. In Fig. 3, we plot the representative values of wind speeds at 0°, 20°, and 70°S from CO and HCN in the context of previous measurements. Within the errors, our values agree with the thermal wind values calculated by Fletcher et al. (2014), both with respect to their reanalysis of the Voyager/IRIS data and for their 2003 Keck data. Our retrieved winds also agree with the wind speed at 0.38 mbar derived by French et al. (1998) from stellar occultations. A comparison with the Voyager 2 cloud-tracking winds (Sromovsky et al. 1993) also confirms a decay of the wind intensity with latitude, as well as a lower value of the wind shear at high latitudes. Although the various observations pertain to different epochs, our results validate predictions of the thermal wind equation and argue for a preservation of the general circulation pattern over the 28 year interval (∼60 deg in heliocentric longitude) spanned by the data around the 2005 southern summer solstice.
Our equatorial wind is about m s−1 less intense than the Voyager reference. Assuming nominal probed levels of ∼1 mbar and ∼1 bar, respectively, this indicates a
m s−1 wind shear per pressure decade (or ∂u/∂z = +30 ± 10 m s−1 per scale height), where the positive sign is related to the retrograde wind direction. At 70°S, our winds are about
m s−1 less intense than Voyager’s. We find a much smaller wind shear at 70°S, although the uncertainties are larger than for equatorial winds: −20 ± 60 m s−1 per pressure decade (or ∂u/∂z = −9 ± 25 m s−1 per scale height). Our results compare well with the estimates from French et al. (1998, their Fig. 11b), who studied the wind-shear between Voyager’s cloud-tracking winds (which they assumed at 100 mbar) and their occultation data at 0.38 mbar. French et al. (1998) determined a wind shear of about +30 m s−1 per scale height at the equator and −15 m s−1 per scale height at 70°S.
In contrast, cloud-tracking measurements from Tollefson & Pater (2018) appear somewhat at odds with our estimated wind shear, as their H-band measurements – assumed by the authors to probe deeper levels – indicate less intense winds than the K′ band. Tollefson & Pater (2018) assumed that the H-band (resp. K′) winds probe the 1–2 bar (resp. 10–100 mbar) level. This led them to an inverted wind shear, with the equatorial winds becoming more intense with increasing altitude. The authors attempted to explain this behavior by invoking a thermal-compositional wind equation that accounts for density changes associated to latitudinal variations of the methane abundance. We find that such an approach is not warranted according to our results. Furthermore, the absolute sounded pressures are uncertain and highly model-dependent and both bands might not be probing such different pressure levels (Tollefson & Pater 2018, Fig. 16 therein). Similarly, Fitzpatrick et al. (2014) found differences among their cloud-tracking H- and K′-band winds. The pressure levels of the observed clouds are also uncertain in this case, with both H- and K’-band clouds spanning pressure levels between 0.1 and 0.6 bar (Fitzpatrick et al. 2014, Fig. 11 therein).
In itself, the consistency of our direct wind measurements with thermal wind calculations does not highlight a particular mechanism responsible for the wind decay with altitude; namely, the thermal wind equation simply states a balance between vertical wind shears and temperature meridional gradients. The wind decay reported here indicates a drag source, which could be the propagation and breaking of gravity and/or planetary waves (common in planetary stratospheres), although this has to be tested in dynamical simulations. On Saturn and Jupiter, interactions between vertically-propagating waves and the mean zonal flow drive the strong acceleration and deceleration of the stratospheric equatorial zonal flow (e.g., Cosentino et al. 2017; Bardet et al. 2022). Wave-breaking as a source of friction was also hypothesized by Ingersoll et al. (2021) to maintain the stacked circulation cells in Jupiter’s upper troposphere.
Our measurements open up a new window on the study of Neptune’s stratospheric dynamics. In addition, our findings provide useful information for general-circulation modelling studies (Liu & Schneider 2010; Milcareck et al. 2021), which require observations to compare with the numerical simulations. Nevertheless, our wind measurements remain modest in precision, as a result of combined limited integration time and low spatial resolution. Future dedicated observations, possibly combined with long-term monitoring (given the duration of Neptune seasons), are expected to yield further insights into the topic.
Acknowledgments
This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.01471.S. (PI: R. Moreno). ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The authors acknowledge the support of the French Agence Nationale de la Recherche (ANR), under grant ANR-20-CE49-0009 (project SOUND). T. Cavalié acknowledges funding from CNES.
References
- Bardet, D., Spiga, A., & Guerlet, S. 2022, Nat. Astron., 6, 804 [NASA ADS] [CrossRef] [Google Scholar]
- Benmahi, B., Cavalié, T., Fouchet, T., et al. 2022, A&A, 666, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cavalié, T., Billebaud, F., Encrenaz, T., et al. 2008, A&A, 489, 795 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cavalié, T., Benmahi, B., Hue, V., et al. 2021, A&A, 647, L8 [EDP Sciences] [Google Scholar]
- Conrath, B., Flasar, F. M., Hanel, R., et al. 1989, Science, 246, 1454 [NASA ADS] [CrossRef] [Google Scholar]
- Cordiner, M. A., Garcia-Berrios, E., Cosentino, R. G., et al. 2020, ApJ, 904, L12 [NASA ADS] [CrossRef] [Google Scholar]
- Cosentino, R. G., Morales-Juberías, R., Greathouse, T., et al. 2017, J. Geophys. Res. (Planets), 122, 2719 [NASA ADS] [CrossRef] [Google Scholar]
- de Pater, I., Fletcher, L. N., Luszcz-Cook, S., et al. 2014, Icarus, 237, 211 [Google Scholar]
- Fitzpatrick, P. J., de Pater, I., Luszcz-Cook, S., Wong, M. H., & Hammel, H. B. 2014, Ap&SS, 350, 65 [Google Scholar]
- Fletcher, L. N., de Pater, I., Orton, G. S., et al. 2014, Icarus, 231, 146 [NASA ADS] [CrossRef] [Google Scholar]
- Fletcher, L. N., Kaspi, Y., Guillot, T., & Showman, A. P. 2020, Space Sci. Rev., 216, 30 [NASA ADS] [CrossRef] [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
- French, R. G., McGhee, C. A., & Sicardy, B. 1998, Icarus, 136, 27 [CrossRef] [Google Scholar]
- Goodman, J., & Weare, J. 2010, Commun. Appl. Math. Comput. Sci., 5, 65 [Google Scholar]
- Hammel, H. B., Beebe, R. F., de Jong, E. M., et al. 1989, Science, 245, 1367 [NASA ADS] [CrossRef] [Google Scholar]
- Helled, R., Anderson, J. D., & Schubert, G. 2010, Icarus, 210, 446 [NASA ADS] [CrossRef] [Google Scholar]
- Högbom, J. A. 1974, A&AS, 15, 417 [Google Scholar]
- Hueso, R., de Pater, I., Simon, A., et al. 2017, Icarus, 295, 89 [NASA ADS] [CrossRef] [Google Scholar]
- Ingersoll, A. P., Atreya, S., Bolton, S. J., et al. 2021, Geophys. Res. Lett., 48, e2021GL095756 [CrossRef] [Google Scholar]
- Irwin, P. G. J., Teanby, N. A., Davis, G. R., et al. 2011, Icarus, 216, 141 [NASA ADS] [CrossRef] [Google Scholar]
- Irwin, P. G. J., Fletcher, L. N., Tice, D., et al. 2016, Icarus, 271, 418 [NASA ADS] [CrossRef] [Google Scholar]
- Karkoschka, E. 2011, Icarus, 215, 439 [NASA ADS] [CrossRef] [Google Scholar]
- Lecacheux, A., Zarka, P., Desch, M. D., & Evans, D. R. 1993, Geophys. Res. Lett., 20, 2711 [NASA ADS] [CrossRef] [Google Scholar]
- Lellouch, E., Goldstein, J. J., Bougher, S. W., Paubert, G., & Rosenqvist, J. 1991, ApJ, 383, 401 [Google Scholar]
- Lellouch, E., Gurwell, M. A., Moreno, R., et al. 2019, Nat. Astron., 3, 614 [NASA ADS] [CrossRef] [Google Scholar]
- Lellouch, E., Paubert, G., Moreno, R., & Moullet, A. 2008, Planet. Space Sci., 56, 1355 [NASA ADS] [CrossRef] [Google Scholar]
- Limaye, S. S., & Sromovsky, L. A. 1991, J. Geophys. Res., 96, 18941P [NASA ADS] [CrossRef] [Google Scholar]
- Liu, J., & Schneider, T. 2010, J. Atmos. Sci., 67, 3652 [NASA ADS] [CrossRef] [Google Scholar]
- Luszcz-Cook, S. H., de Kleer, K., de Pater, I., Adamkovics, M., & Hammel, H. B. 2016, Icarus, 276, 52 [NASA ADS] [CrossRef] [Google Scholar]
- Milcareck, G., Guerlet, S., Spiga, A., et al. 2021, European Planetary Science Congress, EPSC2021-54. [Google Scholar]
- Moreno, R., Marten, A., & Hidayat, T. 2005, A&A, 437, 319 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Moreno, R., Lellouch, E., Forget, F., et al. 2009, Icarus, 201, 549 [NASA ADS] [CrossRef] [Google Scholar]
- Moreno, R., Lellouch, E., Cavalié, T., & Moullet, A. 2017, A&A, 608, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Orton, G. S., Encrenaz, T., Leyrat, C., Puetter, R., & Friedson, A. J. 2007, A&A, 473, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 2007, Numerical Recipes: The Art of Scientific Computing, 3rd edn. (USA: Cambridge University Press) [Google Scholar]
- Smith, B. A., Soderblom, L. A., Banfield, D., et al. 1989, Science, 246, 1422 [NASA ADS] [CrossRef] [Google Scholar]
- Sromovsky, L. A., Limaye, S. S., & Fry, P. M. 1993, Icarus, 105, 110 [NASA ADS] [CrossRef] [Google Scholar]
- Sromovsky, L. A., Fry, P. M., Dowling, T. E., Baines, K. H., & Limaye, S. S. 2001a, Icarus, 149, 459 [NASA ADS] [CrossRef] [Google Scholar]
- Sromovsky, L. A., Fry, P. M., Dowling, T. E., Baines, K. H., & Limaye, S. S. 2001b, Icarus, 150, 244 [CrossRef] [Google Scholar]
- Tollefson, J., de Pater, I., Marcus, P. S., et al. 2018, Icarus, 311, 317 [NASA ADS] [CrossRef] [Google Scholar]
- Warwick, J. W., Evans, D. R., Peltzer, G. R., et al. 1989, Science, 246, 1498 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Measured spectral maps
![]() |
Fig. A.1. Measured spectral maps. Top: Spectral map of the CO(3–2) transition at 345.796 GHz on Neptune observed on April 30, 2016 with ALMA (top row). The spectral resolution is 1 MHz. The planet’s angular diameter is 2.24" and is shown by a blue circle. The synthesized beam is indicated with a red ellipse. The right box represent the scales of the spectra expressed in flux (Jy/beam) against velocity (km s−1). Bottom: Spectral map of the HCN(4–3) transition at 354.505 GHz, measured simultaneously with that of CO. |
Appendix B: Parameterizing wind profiles with different polynomial orders
![]() |
Fig. B.1. As in Fig. 3, top row shows the retrieved wind profiles for CO and HCN, but using different polynomial orders for the wind parameterization. The value of |
All Tables
Polynomial coefficients of the best fits to the CO and HCN zonal wind measurements.
All Figures
![]() |
Fig. 1. Example of measured spectra (left column) of CO(3-2) and HCN(4-3) lines (black) at the position of the sky western equatorial limb, at an offset from Neptune centre (−1.1″, 0.0″). Individual Gaussian fit components with FWHMs of about 4, 20, 80 km s−1 are shown in dashed blue and their sum in red. Also shown: the beam-convolved solid body velocity (black dotted line) and the fit velocity (red dotted line). The difference of these two velocities allows us to derive the Doppler winds. Right: normalized wind contribution functions (WCF) at the limb for each molecule. |
In the text |
![]() |
Fig. 2. ALMA Doppler measurements for the CO(3-2) line (top row) and the HCN(4-3) line (bottom row). Left: measured Doppler lineshift. Middle: line-of-sight winds, after subtracting the solid-body rotation. Right: root mean square (rms) of the measured winds. Top-right ellipse in each subplot shows the synthetic beam. Latitudes are shown in 20° steps, with the equator displayed as the thicker line. |
In the text |
![]() |
Fig. 3. Retrieval results. Left: retrieved best-fitting wind profiles for CO (red line) and HCN (blue line) measurements, using a fourth-order polynomial. Red and blue shadowed regions contain the ensemble of good fits according to the χ2 criterion from Eq. (3). The semi-transparent grey rectangle indicates the unobservable northern latitudes. The solid black line shows the sixth-order fit to Voyager’s cloud-tracking winds (Sromovsky et al. 1993). Right: wind variations with altitude at the equator, 20°S and 70°S, both for our measurements and for a set of references in the literature (see Sect. 5 for details). Cloud-tracking winds from Fitzpatrick et al. (2014) and Tollefson & Pater (2018) are not shown at 70°S as they have uncertainties of about 1000 m s−1. Winds from Fletcher et al. (2014) do not correspond to a direct measurement, but to the computed thermal wind equation applied to a reanalysis of the 1989 IRIS/Voyager data (dashed grey lines) and to 2003 Keck data (dashed cyan lines). |
In the text |
![]() |
Fig. A.1. Measured spectral maps. Top: Spectral map of the CO(3–2) transition at 345.796 GHz on Neptune observed on April 30, 2016 with ALMA (top row). The spectral resolution is 1 MHz. The planet’s angular diameter is 2.24" and is shown by a blue circle. The synthesized beam is indicated with a red ellipse. The right box represent the scales of the spectra expressed in flux (Jy/beam) against velocity (km s−1). Bottom: Spectral map of the HCN(4–3) transition at 354.505 GHz, measured simultaneously with that of CO. |
In the text |
![]() |
Fig. B.1. As in Fig. 3, top row shows the retrieved wind profiles for CO and HCN, but using different polynomial orders for the wind parameterization. The value of |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.