Open Access
Issue
A&A
Volume 648, April 2021
Article Number L8
Number of page(s) 4
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202140837
Published online 26 April 2021

© E. Marcq et al. 2021

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.

1. Introduction

Although sulphur dioxide (SO2) is known to vary by several orders of magnitude at the cloud top of Venus (Esposito 1984; Esposito et al. 1988; Marcq et al. 2013, 2020; Vandaele et al. 2017a,b), such a variability could never be evidenced in the troposphere below the clouds, mostly since remote measurements rely on the sole spectroscopic analysis of a narrow spectral feature in the night-side thermal emission of Venus near 2.46 μm (Bézard et al. 1993). Various spectroscopic analyses of this night-side emission have led to the commonly accepted range for the SO2 mixing ratio of 130 ± 50 ppmv (Marcq et al. 2008; Arney et al. 2014; Vandaele et al. 2017b). Yet, in order to understand SO2 variability above the clouds, and possibly relate it to volcanic outgassing (Esposito 1984), its behaviour in the troposphere (which acts as a proximal reservoir for cloud-top SO2) must be better characterised. Investigations using high-spectral-resolution observations provide a means to achieve this goal since they allow for a better separation between the various minor gaseous species absorbing in the 2.3 − 2.5 μm near-infrared window, namely: CO, OCS, H2O, HDO, and HF (Taylor et al. 1997). To address this, we performed new high-spectral-resolution observations, which are described in Sect. 2. We compare them to our radiative transfer calculations in Sect. 3. We present our co-located SO2 and CO measurements and discuss their implications in Sect. 4, before summarising our findings in the last section, Sect. 5.

2. Observations

2.1. Acquisition

We performed our observations in January 2019 and August 2020, during western elongations of Venus (another observation campaign was scheduled in late April 2020 during an eastern elongation but could not be performed since the Mauna Kea observatory was closed due to the COVID-19 pandemic). We used the infrared high-resolution iSHELL echelle spectrograph (Rayner et al. 2016) mounted at the NASA Infrared Telescope Facility (IRTF) in Hawaii. We chose to use the 5″ × 1.5″ slit as a compromise between minimising acquisition time on the night side of Venus and limiting stray-light contamination from the day side of Venus. In the K3 grating mode of iSHELL, the spectral interval ranges from 2.26 μm to 2.55 μm, covering the longer wavelength portion of the 2.2 − 2.5 μm thermal infrared window of Venus, and extending beyond at longer wavelengths.

Since the angular diameter of Venus during maximal elongations far exceeds 5″, we also recorded images from the guider at 3.46 μm to avoid saturation and properly locate the day side of Venus, including the terminator. These guider images were then compared to ephemeris computations from the Institut de mécanique céleste et de calcul des éphémérides (IMCCE) in order to properly locate the slit geometry projected on the Venusian disk (latitude, local time, emission angle). The slit was always kept parallel to the terminator.

During each session, we recorded (1) spectra from standard A stars used for intensity calibration by the SpeXTool suite, (2) spectra recorded on the night side of Venus, (3) spectra recorded on the day side of Venus, and (4) sky images for background subtraction. Stellar spectra were recorded early during each session, when the sky was dark enough and Venus was too low on the horizon to observe, and at a similar airmass to later Venus observations. The typical required integration time is close to 30 min for a slit position on the night side, 2 min on the day side, and 10 min for calibration stars (see Table 1).

Table 1.

Observation summary.

2.2. Processing

As previously mentioned, we performed the first stages of our data processing using the SpeXTool data reduction package (Cushing et al. 2004; Vacca et al. 2003). In order to increase the signal-to-noise ratio (S/N) and considering the typical seeing extent (∼1″), we divided each slit position on the night side into two separate 2″ × 1.5″ locations. The SpeXTool package provides calibrated spectra, in both wavelength and spectral radiance, as well as error intervals for each of the 27 dispersion orders of the echelle grating.

Following our past IRTF observations of Venus using SpeX (Marcq et al. 2005, 2006), we cleaned the night-side thermal emission spectra from the day-side solar stray light by subtracting a day-side spectrum acquired near the terminator. The day-side spectrum was interpolated using several day-side observations in order to match the same terrestrial airmass as the night-side observation. The amount of stray light to subtract was determined by taking advantage of the saturation of spectrally resolved CO lines in the 2.32 − 2.35 μm interval (Fig. 1) since telluric absorption is too strong to make use of wavelengths longer than 2.5 μm. Data were then filtered by removing outliers (resulting from the faulty correction of some telluric absorption lines by SpeXTool) and then spectrally smoothed in order to increase the S/N while keeping our effective spectral resolution close to λλ ≃ 20 000.

thumbnail Fig. 1.

Raw night-side spectrum (blue), interpolated day-side spectrum (×0.05, green), and resulting cleaned night-side spectrum (red) acquired on January 3, 2019.

3. Fitting

3.1. Forward radiative transfer model

We employed an updated version of the radiative transfer model previously used by Marcq et al. (2005, 2006, 2008). The two main updates are the following: (1) An eight-stream radiative transfer solver based on DISORT (Stamnes et al. 1988) was used, which allowed us to take the effect of a variation in emission angle into account; and (2) the CO2 line database was updated using the ab initio list from Huang et al. (2014).

3.1.1. Limb darkening

We found that our model yields a remarkably simple wavelength-independent limb darkening parametrisation: I(μ) = I(μ = 1) × (0.74μ + 0.26), where μ stands for the cosine of the emission angle. This is easily understood considering the strong multiple scattering in the optically thick clouds above the emission region and the very small opacity of the atmosphere above the clouds.

3.1.2. Vertical profiles of gaseous species

Our updated model uses the same method of specifying the vertical profiles of minor species as in Marcq et al. (2008): Volume mixing ratios (qi) are specified for a limited number of pressure levels (Pi), and the mixing ratios are linearly interpolated (in log q vs. log P coordinates) in between. The nominal qi(Pi) values are given in Table 2 for CO and SO2.

Table 2.

Nominal vertical profiles for CO and SO2.

3.2. Fitting algorithm

The fitting algorithm proceeds in two stages; a Levenberg-Marquardt loop (Newville et al. 2014) interpolating from look-up tables computed by the forward radiative transfer model described in Sect. 3.1 is used in both. In the first stage, lower cloud opacity and a carbon monoxide scaling factor are retrieved using grating orders 13 to 19 (2.325 to 2.402 μm) and assuming the same CO vertical profile shape as Marcq et al. (2008). In the second stage, assuming the previously retrieved CO and lower cloud opacity values and using grating orders 7 to 10 (2.422 to 2.470 μm), OCS mixing ratios at 30 and 37 km are retrieved, as is an SO2 scaling factor. A typical fit is shown in Fig. 2.

thumbnail Fig. 2.

Spectrum observed on January 3, 2019 (colour coding for the different, partially overlapping grating orders) and its best fit (in black). The grey shaded area near 2.46 μm shows the sensitivity to the SO2 mixing ratio (shaded between ×2 and ×0.5 the best fit value).

4. Results

4.1. Carbon monoxide (CO)

Our retrievals for CO are shown in Fig. 3. These results are in agreement with previous measurements (Marcq et al. 2018): Carbon monoxide increases with increasing latitude, although the latitude of minimum CO is shifted to circa 15°S. Such a shift was also seen by previous observers (Marcq et al. 2008; Arney et al. 2014) and interpreted as the influence of topography (especially in Aphrodite Terra) on the large-scale circulation, which is responsible for creating this latitudinal gradient – CO being produced photochemically at higher altitudes and carried below the cloud through the downwelling branch of the Hadley cell (Tsang et al. 2008).

thumbnail Fig. 3.

CO and SO2 volume mixing ratio measurements. Error bars stand for ±3 ⋅ σ standard deviations. Systematic uncertainties of the absolute SO2 mixing ratio are not displayed.

Additionally, this shift makes our observed CO distribution less correlated with the emission angle than for a hemispherically symmetric CO distribution. It is therefore less likely that it originates from a faulty limb darkening correction, which further increases our confidence in these results.

4.2. Sulphur dioxide (SO2)

Our co-located retrievals for SO2 are also shown in Fig. 3. Previous measurements suffered very large error bars due to the narrowness of SO2 spectral features near 2.46 μm, and they agreed on a latitudinally constant value of 130 ± 50 ppmv below the clouds (Vandaele et al. 2017a; Marcq et al. 2018). Very interestingly, Arney et al. (2014, their Fig. 26) observed a quantitatively similar pattern in their SO2 maps but dismissed it on the basis of the large error bars and other ‘ghosting effects’ inherent to their low-to-medium spectral resolution. However, our high-resolution iSHELL observations circumvent this issue, such that an increase in SO2 with increasing latitude is unambiguous in our data. Interestingly, the latitudinal offset is the same as for CO, hinting at a common mechanism behind these variations, which is further discussed in Sect. 4.3.

We also note that our latitudinally averaged value is closer to 180 ppmv, 40% higher than the commonly accepted value. However, we found that this average value depends on several factors, such as the assumed CO2-CO2 continuum value, the parametrisation of cloud opacity, and the exact amount of stray light that is subtracted. On the other hand, relative variations are much better constrained, and so we focus on this new result hereafter.

4.3. Discussion

The aforementioned correlation (R2 ≃ 0.694) between SO2 and CO can be seen in Fig. 4 and hints at a common origin for their horizontal variability. For CO, this tropospheric variability comes from the competition between the positive vertical gradient (itself originating from the relative altitude of its photochemical source, which is located above the sounded altitude, and a lower atmospheric sink, which is located below) and the general Hadley-cell circulation (Tsang et al. 2008). We therefore infer that SO2 should also exhibit a positive vertical gradient below the clouds in order to account for its observed latitudinal variability. We can even estimate the typical value of the required vertical gradient based on the already known CO gradient value, assuming there are no sources or sinks for either species and that the probed altitudes are similar: Since the typical gradient for CO is close to +1 ppm km−1 (Marcq et al. 2018), we estimate a SO2 vertical gradient of the order of +10 ppm km−1 in the probed altitude range.

thumbnail Fig. 4.

Correlation plot between CO and SO2 measurements. Error bars stand for ±3 ⋅ σ standard deviations, and colour is coded for the measurement latitudes. Systematic uncertainties of the absolute SO2 mixing ratio are not displayed.

Coincidentally, this gradient value is close to the value found by in situ ISAV measurements on board VEGA descent probes (Bertaux et al. 1996), though ISAV measured substantially lower mixing ratios for SO2 than what we found. However, currently available thermochemical models of the lower atmosphere (Krasnopolsky 2007) do not predict such a vertical gradient (even suggesting a moderate decrease in the 30 − 40 km altitude range). Further modelling work is definitely required in order to account for this newly discovered SO2 variability.

5. Conclusion

Our iSHELL observations of the night-side thermal spectra of Venus showed, for the first time, an increase in SO2 with increasing latitude below the clouds of Venus, which we tentatively ascribe to the general circulation acting on an increasing-with-height SO2 mixing ratio. This study also highlights the interest of high-spectral-resolution infrared spectra for accurate measurements of trace gases in Venus’s atmosphere, either from long-term-monitoring ground-based instruments, such as iSHELL, or dedicated space-borne instruments, such as the VenSpec-H channel (Robert et al., in prep.) on board ESA M5 candidate EnVision.

Acknowledgments

E. Marcq wishes to thank L. Bercot for his logistic support during remote IRTF observations in August 2020, as well as A. Boogert, G. Osterman, G. Engh, and J. Rayner for their support during iSHELL operations in both 2019 and 2020 observations campaigns. The authors also thank C. F. Wilson for his valuable comments and remarks.

References

  1. Arney, G., Meadows, V., Crisp, D., et al. 2014, J. Geophys. Res.: Planets, 119, 1860 [Google Scholar]
  2. Bertaux, J., Widemann, T., Hauchecorne, A., Moroz, V. I., & Ekonomov, A. P. 1996, J. Geophys. Res, 101, 12709 [Google Scholar]
  3. Bézard, B., de Bergh, C., Fegley, B., et al. 1993, Geophys. Res. Lett., 20, 1587 [Google Scholar]
  4. Cushing, M. C., Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 362 [Google Scholar]
  5. Esposito, L. W. 1984, Science, 223, 1072 [Google Scholar]
  6. Esposito, L. W., Copley, M., Eckert, R., et al. 1988, J. Geophys. Res., 93, 5267 [Google Scholar]
  7. Huang, X., Gamache, R. R., Freedman, R. S., Schwenke, D. W., & Lee, T. J. 2014, J. Quant. Spectr. Rad. Transf., 147, 134 [Google Scholar]
  8. Krasnopolsky, V. A. 2007, Icarus, 191, 25 [Google Scholar]
  9. Marcq, E., Bézard, B., Encrenaz, T., & Birlan, M. 2005, Icarus, 179, 375 [Google Scholar]
  10. Marcq, E., Encrenaz, T., Bézard, B., & Birlan, M. 2006, Planet. Space Sci., 54, 1360 [Google Scholar]
  11. Marcq, E., Bézard, B., Drossart, P., et al. 2008, J. Geophys. Res.: Planets, 113, E00B07 [Google Scholar]
  12. Marcq, E., Bertaux, J.-L., Montmessin, F., & Belyaev, D. 2013, Nature Geosci., 6, 25 [Google Scholar]
  13. Marcq, E., Mills, F. P., Parkinson, C. D., & Vandaele, A. C. 2018, Space Sci. Rev., 214, 10 [Google Scholar]
  14. Marcq, E., Lea Jessup, K., Baggio, L., et al. 2020, Icarus, 335, 113368 [Google Scholar]
  15. Newville, M., Stensitzki, T., Allen, D. B., & Ingargiola, A. 2014, https://doi.org/10.5281/zenodo.11813 [Google Scholar]
  16. Rayner, J. T., Tokunaga, A., Jaffe, D., et al. 2016, SPIE Conf. Ser., 9908E, 84R [Google Scholar]
  17. Stamnes, K., Tsay, S., Jayaweera, K., & Wiscombe, W. 1988, Appl. Opt., 27, 2502 [Google Scholar]
  18. Taylor, F. W., Crisp, D., & Bézard, B. 1997, in Venus II, eds. D. M. Hunten, L. Glin, T. M. Donahue, & V. I. Moroz (Tucson: University of Arizona Press), 325 [Google Scholar]
  19. Tsang, C. C. C., Irwin, P. G. J., Wilson, C. F., et al. 2008, J. Geophys. Res.: Planets, 113, E00B08 [Google Scholar]
  20. Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389 [Google Scholar]
  21. Vandaele, A. C., Korablev, O., Belyaev, D., et al. 2017a, Icarus, 295, 16 [Google Scholar]
  22. Vandaele, A. C., Korablev, O., Belyaev, D., et al. 2017b, Icarus, 295, 1 [Google Scholar]

All Tables

Table 1.

Observation summary.

Table 2.

Nominal vertical profiles for CO and SO2.

All Figures

thumbnail Fig. 1.

Raw night-side spectrum (blue), interpolated day-side spectrum (×0.05, green), and resulting cleaned night-side spectrum (red) acquired on January 3, 2019.

In the text
thumbnail Fig. 2.

Spectrum observed on January 3, 2019 (colour coding for the different, partially overlapping grating orders) and its best fit (in black). The grey shaded area near 2.46 μm shows the sensitivity to the SO2 mixing ratio (shaded between ×2 and ×0.5 the best fit value).

In the text
thumbnail Fig. 3.

CO and SO2 volume mixing ratio measurements. Error bars stand for ±3 ⋅ σ standard deviations. Systematic uncertainties of the absolute SO2 mixing ratio are not displayed.

In the text
thumbnail Fig. 4.

Correlation plot between CO and SO2 measurements. Error bars stand for ±3 ⋅ σ standard deviations, and colour is coded for the measurement latitudes. Systematic uncertainties of the absolute SO2 mixing ratio are not displayed.

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.