Open Access
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

© The Authors 2023

Licence Creative CommonsOpen 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).

thumbnail 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.

thumbnail 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 2 1.8 + 12 $ 2^{+12}_{-1.8} $ mbar level and HCN lines probe the 0 . 4 0.3 + 0.5 $ 0.4^{+0.5}_{-0.3} $ 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:

W = n ( a n × ϕ n ) , $$ \begin{aligned} W = \sum _n (a_n \times \phi ^n) ,\end{aligned} $$(1)

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:

χ 2 = i = 1 N p ( W test W ALMA rms ALMA ) 2 , $$ \begin{aligned} \chi ^2 = \sum _{i=1}^{N_{\rm p}}\left(\frac{W_{\rm test}-W_{\rm ALMA}}{\mathrm{rms_{ALMA}}} \right)^2 ,\end{aligned} $$(2)

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:

χ 2 χ min 2 C × χ min 2 N × f oversampling . $$ \begin{aligned} \chi ^2 - \chi ^2_{\rm min} \le C \times \frac{\chi ^2_{\rm min}}{N} \times f_{\rm oversampling} .\end{aligned} $$(3)

Here, χ min 2 $ \chi^2_{\rm min} $ is the minimum χ2 value among all the samples in the run, and χ min 2 /N $ \chi^2_{\rm min}/N $ is the reduced value of χ min 2 $ \chi^2_{\rm min} $, 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).

thumbnail 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).

Table 1.

Polynomial coefficients of the best fits to the CO and HCN zonal wind measurements.

At the equator, we retrieved retrograde (westward) zonal winds of 180 60 + 70 $ -180 ^{+70}_{-60} $ m s−1 from CO measurements and 190 70 + 90 $ -190 ^{+90}_{-70} $ m s−1 from HCN. We find a decrease in the wind intensity towards mid-latitudes. At 20°S, wind speeds are 170 40 + 50 $ -170 ^{+50}_{-40} $ m s−1 for CO ( 140 60 + 60 $ -140 ^{+60}_{-60} $ m s−1 for HCN) and at 40°S, they are 90 60 + 50 $ -90 ^{+50}_{-60} $ m s−1 for CO ( 40 80 + 60 $ -40 ^{+60}_{-80} $ 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 ( 180 110 + 130 $ 180 ^{+130}_{-110} $ m s−1 for CO, and 180 140 + 110 $ 180 ^{+110}_{-140} $ 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 150 80 + 30 $ -150 ^{+30}_{-80} $ m s−1 for CO and 170 50 + 80 $ -170 ^{+80}_{-50} $ 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 200 80 + 100 $ 200^{+100}_{-80} $ m s−1 less intense than the Voyager reference. Assuming nominal probed levels of ∼1 mbar and ∼1 bar, respectively, this indicates a + 70 20 + 30 $ +70^{+30}_{-20} $ 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 70 170 + 180 $ 70^{+180}_{-170} $ 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

  1. Bardet, D., Spiga, A., & Guerlet, S. 2022, Nat. Astron., 6, 804 [NASA ADS] [CrossRef] [Google Scholar]
  2. Benmahi, B., Cavalié, T., Fouchet, T., et al. 2022, A&A, 666, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Cavalié, T., Billebaud, F., Encrenaz, T., et al. 2008, A&A, 489, 795 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Cavalié, T., Benmahi, B., Hue, V., et al. 2021, A&A, 647, L8 [EDP Sciences] [Google Scholar]
  5. Conrath, B., Flasar, F. M., Hanel, R., et al. 1989, Science, 246, 1454 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cordiner, M. A., Garcia-Berrios, E., Cosentino, R. G., et al. 2020, ApJ, 904, L12 [NASA ADS] [CrossRef] [Google Scholar]
  7. Cosentino, R. G., Morales-Juberías, R., Greathouse, T., et al. 2017, J. Geophys. Res. (Planets), 122, 2719 [NASA ADS] [CrossRef] [Google Scholar]
  8. de Pater, I., Fletcher, L. N., Luszcz-Cook, S., et al. 2014, Icarus, 237, 211 [Google Scholar]
  9. Fitzpatrick, P. J., de Pater, I., Luszcz-Cook, S., Wong, M. H., & Hammel, H. B. 2014, Ap&SS, 350, 65 [Google Scholar]
  10. Fletcher, L. N., de Pater, I., Orton, G. S., et al. 2014, Icarus, 231, 146 [NASA ADS] [CrossRef] [Google Scholar]
  11. Fletcher, L. N., Kaspi, Y., Guillot, T., & Showman, A. P. 2020, Space Sci. Rev., 216, 30 [NASA ADS] [CrossRef] [Google Scholar]
  12. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  13. French, R. G., McGhee, C. A., & Sicardy, B. 1998, Icarus, 136, 27 [CrossRef] [Google Scholar]
  14. Goodman, J., & Weare, J. 2010, Commun. Appl. Math. Comput. Sci., 5, 65 [Google Scholar]
  15. Hammel, H. B., Beebe, R. F., de Jong, E. M., et al. 1989, Science, 245, 1367 [NASA ADS] [CrossRef] [Google Scholar]
  16. Helled, R., Anderson, J. D., & Schubert, G. 2010, Icarus, 210, 446 [NASA ADS] [CrossRef] [Google Scholar]
  17. Högbom, J. A. 1974, A&AS, 15, 417 [Google Scholar]
  18. Hueso, R., de Pater, I., Simon, A., et al. 2017, Icarus, 295, 89 [NASA ADS] [CrossRef] [Google Scholar]
  19. Ingersoll, A. P., Atreya, S., Bolton, S. J., et al. 2021, Geophys. Res. Lett., 48, e2021GL095756 [CrossRef] [Google Scholar]
  20. Irwin, P. G. J., Teanby, N. A., Davis, G. R., et al. 2011, Icarus, 216, 141 [NASA ADS] [CrossRef] [Google Scholar]
  21. Irwin, P. G. J., Fletcher, L. N., Tice, D., et al. 2016, Icarus, 271, 418 [NASA ADS] [CrossRef] [Google Scholar]
  22. Karkoschka, E. 2011, Icarus, 215, 439 [NASA ADS] [CrossRef] [Google Scholar]
  23. Lecacheux, A., Zarka, P., Desch, M. D., & Evans, D. R. 1993, Geophys. Res. Lett., 20, 2711 [NASA ADS] [CrossRef] [Google Scholar]
  24. Lellouch, E., Goldstein, J. J., Bougher, S. W., Paubert, G., & Rosenqvist, J. 1991, ApJ, 383, 401 [Google Scholar]
  25. Lellouch, E., Gurwell, M. A., Moreno, R., et al. 2019, Nat. Astron., 3, 614 [NASA ADS] [CrossRef] [Google Scholar]
  26. Lellouch, E., Paubert, G., Moreno, R., & Moullet, A. 2008, Planet. Space Sci., 56, 1355 [NASA ADS] [CrossRef] [Google Scholar]
  27. Limaye, S. S., & Sromovsky, L. A. 1991, J. Geophys. Res., 96, 18941P [NASA ADS] [CrossRef] [Google Scholar]
  28. Liu, J., & Schneider, T. 2010, J. Atmos. Sci., 67, 3652 [NASA ADS] [CrossRef] [Google Scholar]
  29. Luszcz-Cook, S. H., de Kleer, K., de Pater, I., Adamkovics, M., & Hammel, H. B. 2016, Icarus, 276, 52 [NASA ADS] [CrossRef] [Google Scholar]
  30. Milcareck, G., Guerlet, S., Spiga, A., et al. 2021, European Planetary Science Congress, EPSC2021-54. [Google Scholar]
  31. Moreno, R., Marten, A., & Hidayat, T. 2005, A&A, 437, 319 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Moreno, R., Lellouch, E., Forget, F., et al. 2009, Icarus, 201, 549 [NASA ADS] [CrossRef] [Google Scholar]
  33. Moreno, R., Lellouch, E., Cavalié, T., & Moullet, A. 2017, A&A, 608, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Orton, G. S., Encrenaz, T., Leyrat, C., Puetter, R., & Friedson, A. J. 2007, A&A, 473, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. 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]
  36. Smith, B. A., Soderblom, L. A., Banfield, D., et al. 1989, Science, 246, 1422 [NASA ADS] [CrossRef] [Google Scholar]
  37. Sromovsky, L. A., Limaye, S. S., & Fry, P. M. 1993, Icarus, 105, 110 [NASA ADS] [CrossRef] [Google Scholar]
  38. Sromovsky, L. A., Fry, P. M., Dowling, T. E., Baines, K. H., & Limaye, S. S. 2001a, Icarus, 149, 459 [NASA ADS] [CrossRef] [Google Scholar]
  39. Sromovsky, L. A., Fry, P. M., Dowling, T. E., Baines, K. H., & Limaye, S. S. 2001b, Icarus, 150, 244 [CrossRef] [Google Scholar]
  40. Tollefson, J., de Pater, I., Marcus, P. S., et al. 2018, Icarus, 311, 317 [NASA ADS] [CrossRef] [Google Scholar]
  41. 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

thumbnail 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

thumbnail 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 χ min 2 /N $ \chi^2_{min}/N $ is shown in the legend. We also show the residual maps of CO (middle row) and HCN (bottom row) after subtracting the best-fitting wind profile to the measurements.

All Tables

Table 1.

Polynomial coefficients of the best fits to the CO and HCN zonal wind measurements.

All Figures

thumbnail 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
thumbnail 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
thumbnail 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
thumbnail 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
thumbnail 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 χ min 2 /N $ \chi^2_{min}/N $ is shown in the legend. We also show the residual maps of CO (middle row) and HCN (bottom row) after subtracting the best-fitting wind profile to the measurements.

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.