ALMA Spectral Imaging of Titan Contemporaneous with Cassini's Grand Finale
M. A. Cordiner, N. A. Teanby, C. A. Nixon, V. Vuitton, A. E. Thelen,, S. B. Charnley

TL;DR
This study uses ALMA to perform spectral imaging of Titan's atmosphere, complementing Cassini data, revealing detailed spatial and vertical distributions of various molecules and their seasonal and dynamical variations.
Contribution
First detailed ALMA measurements of Titan's atmospheric molecular distributions, providing new insights into chemical lifetimes and spatial asymmetries, complementing Cassini observations.
Findings
ALMA maps of HCN, HNC, HC3N, C2H3CN, C2H5CN, C3H8 show seasonal and dynamical effects.
HNC is confined above 730 km with east-west hemispheric asymmetry.
Chemical lifetimes of C2H3CN and C2H5CN are a few years, matching chemical models.
Abstract
The Cassini mission performed 127 targeted flybys of Titan during its 13-year mission to Saturn, culminating in the Grand Finale between April-September 2017. Here we demonstrate the use of the Atacama Large Millimeter/submillimeter Array (ALMA) to continue Cassini's legacy for chemical and climatological studies of Titan's atmosphere. Whole-hemisphere, interferometric spectral maps of HCN, HNC, HC3N, CH3CN, C2H3CN, C2H5CN and C3H8 were obtained using ALMA in May 2017 at moderate (~0.2'', or 1300 km) spatial resolution, revealing the effects of seasonally-variable chemistry and dynamics on the distribution of each species. The ALMA sub-mm observations of HCN and HC3N are consistent with Cassini infrared data on these species, obtained in the same month. Chemical/dynamical lifetimes of a few years are inferred for C2H3CN and C2H5CN, in reasonably close agreement with the latest chemical…
| Species | Rest Freq. | Transition | Spec. Res. | |
|---|---|---|---|---|
| (MHz) | (kHz) | (K) | ||
| CO | 345796 | 977 | 33.2 | |
| HCN | 354505 | 977 | 42.5 | |
| HCN | 354460 | 977 | 1067 | |
| HNC | 362630 | 244 | 42.5 | |
| H13CN | 345340 | 977 | 41.4 | |
| H13CN | 345239 | 977 | 1057 | |
| HC15N | 344200 | 977 | 41.3 | |
| HC3N | 345609 | 977 | 324 | |
| HC3N | 354697 | 977 | 341 | |
| HC3N | 355278 | 977 | 1059 | |
| HC3N | 355557 | 977 | 1059 | |
| HC3N | 355566 | 977 | 662 | |
| HC3N | 356072 | 977 | 662 | |
| HC3N | 363785 | 244 | 358 | |
| H13CCCN | 343739 | 977 | 330 | |
| H13CCCN | 361353 | 488 | 364 | |
| HC13CCN | 344143 | 977 | 322 | |
| HC315N | 344385 | 977 | 331 | |
| CH3CN | 349025 | 488 | 624 | |
| CH3CN | 349125 | 488 | 517 | |
| CH3CN | 349212 | 488 | 425 | |
| CH3CN | 349286 | 488 | 346 | |
| CH3CN | 349346 | 488 | 282 | |
| CH3CN | 349393 | 488 | 232 | |
| CH3CN | 349427 | 488 | 196 | |
| CH3CN | 349447 | 488 | 175 | |
| CH3CN | 349454 | 488 | 168 | |
| CH313CN | 349221 | 488 | 232 | |
| CH313CN | 349254 | 488 | 196 | |
| CH313CN | 349274 | 488 | 175 | |
| CH313CN | 349280 | 488 | 168 | |
| C2H3CN | 360876 | 488 | 476 | |
| C2H3CN | 360876 | 488 | 476 | |
| C2H3CN | 360878 | 488 | 512 | |
| C2H3CN | 360878 | 488 | 512 | |
| C2H3CN | 361652 | 488 | 373 | |
| C2H5CN | 343195 | 977 | 347 | |
| C2H5CN | 354477 | 977 | 361 | |
| C2H5CN | 355756 | 977 | 351 | |
| C2H5CN | 349443 | 488 | 446 | |
| C2H5CN | 349443 | 488 | 446 | |
| C2H5CN | 349547 | 488 | 425 | |
| C2H5CN | 349547 | 488 | 425 | |
| C2H5CN | 349731 | 488 | 407 | |
| C2H5CN | 349731 | 488 | 407 | |
| C2H5CN | 349796 | 488 | 354 | |
| C3H8 | 360978 | 488 | 64.2 |
| Species | (km) | (km) | (yr)a |
| HCN | 2670 | 530 | 350 |
| HNC | 3140 | 970 | 0.6 |
| HC3N | 2930 | 760 | 2.6 |
| C2H3CN | — | — | 1.0 |
| C2H5CN | 2670 | 610 | 3.2b |
| CH3CN | 2590 | 270 | 91 |
| C3H8 | 2400 | 130 | 13 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
ALMA Spectral Imaging of Titan Contemporaneous with Cassini’s Grand Finale
M. A. Cordiner11affiliation: NASA Goddard Space Flight Center, 8800 Greenbelt Road, MD 20771, USA. 22affiliation: Institute for Astrophysics and Computational Sciences, The Catholic University of America, Washington, DC 20064, USA. , N. A. Teanby33affiliation: School of Earth Sciences, University of Bristol, Wills Memorial Building, Queens Road, Bristol, BS8 1RJ, UK. , C. A. Nixon11affiliation: NASA Goddard Space Flight Center, 8800 Greenbelt Road, MD 20771, USA. , V. Vuitton44affiliation: Université Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France , A. E. Thelen11affiliation: NASA Goddard Space Flight Center, 8800 Greenbelt Road, MD 20771, USA. 22affiliation: Institute for Astrophysics and Computational Sciences, The Catholic University of America, Washington, DC 20064, USA. , S. B. Charnley11affiliation: NASA Goddard Space Flight Center, 8800 Greenbelt Road, MD 20771, USA.
(AJ; submitted March 2019; accepted for publication June 2019)
Abstract
The Cassini mission performed 127 targeted flybys of Titan during its 13-year mission to Saturn, culminating in the Grand Finale between April-September 2017. Here we demonstrate the use of the Atacama Large Millimeter/submillimeter Array (ALMA) to continue Cassini’s legacy for chemical and climatological studies of Titan’s atmosphere. Whole-hemisphere, interferometric spectral maps of HCN, HNC, HC3N, CH3CN, C2H3CN, C2H5CN and C3H8 were obtained using ALMA in May 2017 at moderate (, or km) spatial resolution, revealing the effects of seasonally-variable chemistry and dynamics on the distribution of each species. The ALMA sub-mm observations of HCN and HC3N are consistent with Cassini infrared data on these species, obtained in the same month. Chemical/dynamical lifetimes of a few years are inferred for C2H3CN and C2H5CN, in reasonably close agreement with the latest chemical models incorporating sticking of C2H5CN to stratospheric aerosol particles. ALMA radial limb flux profiles provide column density information as a function of altitude, revealing maximum abundances in the thermosphere (above 600 km) for HCN, HNC, HC3N and C2H5CN. This constitutes the first detailed measurement of the spatial distribution of HNC, which is found to be confined predominantly to altitudes above km. The HNC emission map shows an east-west hemispheric asymmetry of %. These results are consistent with very rapid production (and loss) of HNC in Titan’s uppermost atmosphere, making this molecule an effective probe of short-timescale (diurnal) ionospheric processes.
Subject headings:
planets and satellites: individual (Titan) — planets and satellites: atmospheres — techniques: interferometric — techniques: imaging spectroscopy — submillimeter: planetary systems
1. Introduction
Saturn’s largest moon Titan has a thick (1.45 bar) atmosphere composed primarily of molecular nitrogen (%) and methane (%). Remote and in-situ measurements have identified a wealth of other species therein (see Bézard et al., 2014), the abundances, distributions and temperatures of which are linked to seasonal variations in Titan’s combined photochemistry and atmospheric circulation patterns (e.g. Teanby et al. 2009; Vinatier et al. 2010; Achterberg et al. 2011; Cordiner et al. 2015; Thelen et al. 2019). Hydrocarbons and nitriles spanning a broad range in size and complexity have been detected, from acetylene (C2H2) and hydrogen cyanide (HCN) to benzene (C6H6), propionitrile (CH2CH3CN) and organic aerosols (Hörst, 2017), which are theorized to originate from photochemistry initiated primarily by N2 and CH4 photolysis in the upper atmosphere (Vuitton et al., 2019). However, despite 13 years of dedicated study by the Cassini spacecraft during over 120 scientific flybys (from 2004-2017), major gaps remain in our knowledge of Titan’s atmospheric chemistry and dynamics (Hörst, 2017; Nixon et al., 2018), particularly regarding the sources and sinks of complex organic molecules.
Strong latitudinal anisotropies in the abundances of Titan’s smaller hydrocarbons and nitriles (such as HCN, HC3N and C2H2) were first observed by Voyager 1 (Coustenis & Bézard, 1995). Continued monitoring of these, and other species, by the Cassini Composite Infrared Spectrometer (CIRS) has revealed temporal variations in the abundance patterns of photochemically-produced gases over time periods as short as (Earth) years or even months (Teanby et al., 2012; Coustenis et al., 2013; Vinatier et al., 2015; Teanby et al., 2017). The observed changes in molecular distributions are unique among the bodies of the Solar System, and are believed to be driven by seasonal insolation variations, mediated by a global meridional circulation cell that changes direction during Titan’s 29.5 yr seasonal cycle. To provide new insights into Titan’s atmospheric chemistry and global circulation system, we undertook a program to perform comprehensive mapping of short-lived photochemical products using the Atacama Large Millimeter/submillimeter Array (ALMA) and Cassini CIRS during the transition to Titan’s northern summer solstice in May 2017. These observations capitalize on (1) the recent completion of the full array of m ALMA antennas, enabling the highest-resolution studies of Titan to-date at mm/sub-mm wavelengths (see also Lellouch et al. 2019), and (2) the final opportunity for contemporaneous observations using ALMA and CIRS. The aim of this study was to combine the complementary diagnostic power of the sub-mm and infrared spectral regions to provide the most complete picture to-date of Titan’s rapidly evolving chemistry and dynamics.
Cassini CIRS spectroscopy benefits from high sensitivity and spatial resolution for measuring rovibrational molecular emission lines, but provides an incomplete picture of Titan, hindered by a relatively low spectral resolution () and reduced sensitivity at far-infrared wavelengths. Moreover, CIRS is only sensitive to a subset of Titan’s gases (including CO, HCN, H2O, HC3N, C2N2 and various CnHm hydrocarbons), and crucially, lacks the ability to detect several important and abundant nitrogen-bearing species. Polar nitriles (such as HNC, HC3N, C2H5CN and C2H3CN) are among the species with the shortest predicted chemical lifetimes (less than a few years in the stratosphere; see Wilson & Atreya 2004; Krasnopolsky 2009; Vuitton et al. 2019), and are therefore extremely useful probes of Titan’s short-timescale atmospheric variability. The mm/sub-mm provides optimal sensitivity for rotational emission from polar nitriles, and with its combined high sensitivity and spatial resolution, ALMA can provide detailed maps of the emission from these species in Titan’s atmosphere. Through radiative transfer analysis, the combined Cassini plus ALMA dataset enables molecular abundances to be derived as a function of spatial coordinate, from which a comprehensive evaluation of our understanding of the dominant nitrile formation and destruction pathways in nitrogen and carbon-rich planetary atmospheres can be obtained. Such information is crucial for constraining theories regarding the formation of aerosols and complex nitrogen-bearing organics (including possible pre-biotic species) in primitive (exo-)planetary atmospheres throughout the Galaxy.
The observations presented here provide an overview of the complete ALMA dataset obtained in synergy with the final Titan flybys by Cassini. Molecular maps for HCN, HNC, HC3N, CH3CN, C2H3CN, C2H5CN and C3H8 demonstrate ALMA’s unique, high-sensitivity interferometric imaging capabilities, resulting in spatially resolved distributions for individual molecules, from which we derive preliminary abundances and chemical/dynamical decay timescales. By comparison with results from the latest atmospheric chemical models, new insights are obtained into the relationship between Titan’s climate, chemistry, and global atmospheric circulation system.
2. ALMA Observations
Observations of Titan were carried out using the ALMA Band 7 receiver on 2017-05-08 and 2017-05-16, as part of the director’s discretionary time program 2016.A.00014.S. The observations were targeted to be close in time to Cassini’s May 24 (solstice) flyby of Titan, for synergistic observations with CIRS. The ALMA correlator was configured using two different spectral settings to observe multiple frequency windows covering our molecules of interest in the range (1) 342-356 GHz (877-842 m) and (2) 349-364 GHz (859-824 m), at a spectral resolution of 244-977 kHz (0.2-0.8 km ). For more details of the targeted spectral lines, see Table 1.
The ALMA configuration was moderately extended, with 46 active antennas providing baselines in the range 15-1124 m, resulting in a spatial resolution of (using natural visibility weighting). The total on-source observing time was 18 min in setting 1 and 138 min in setting 2, leading to RMS noise levels in the range 2-4 mJy beam*-1* MHz*-1*. Weather conditions were very good, with zenith precipitable water vapor (PWV) in the range 0.74-0.83 mm, resulting in low atmospheric opacity and good phase stability.
After the initial calibration of the antenna phase delays, the observations consisted of bandpass and flux calibration scans, followed by an interleaved sequence of three visits each to Titan and the phase calibrator J1751-1950. The phase-center was updated in real-time to track Titan’s moving position on the sky. The flux of the calibration quasar J1733-1304 was measured on 2017-05-04 and 2017-05-17, with an accuracy of %, implying a similar accuracy in the flux scale of our calibrated spectra.
Data were flagged and calibrated in CASA 5.1 (Jaeger, 2008) using the automated pipeline scripts supplied by the Joint ALMA Observatory (see Shinnaga et al., 2015). The data from May 8 and May 16 were combined after regridding and Doppler-shifting to Titan’s rest frame. A small offset in declination (of unknown origin) was identified and corrected in Titan’s position with respect to the ALMA phase center on May 8 ( in setting 1 and in setting 2). Titan’s continuum flux was subtracted using low-order polynomial fits to the (line-free) spectral regions adjacent to our lines of interest. Imaging and deconvolution were performed using the Clark clean algorithm with a pixel size of , a flux threshold of twice the expected RMS noise and a mask diameter of (8700 km at Titan’s distance of 9.26 AU), encircling the entirety of the detected flux from Titan’s atmosphere.
The coordinate scales of the cleaned images were transformed to physical distances with respect to the center of Titan. Titan’s north pole was oriented counter-clockwise from celestial north, and tilted towards the observer by . This is close to the maximum polar tilt due to the proximity of our observations to Titan’s southern winter solstice on 2017-05-24.
3. Cassini Observations
Cassini’s Composite Infra-Red Spectrometer (CIRS) is a Fourier transfer spectrometer covering 10-1500 cm*-1* over three focal planes: FP1 (10-600 cm*-1*); FP3 (600-1100 cm*-1*); and FP4 (1100-1500 cm*-1*). For additional details regarding the CIRS instrument see Flasar et al. (2004) and Jennings et al. (2017). Here, we use data obtained using the FP3 and FP4 focal planes, which comprise linear arrays of 10 pixels each with a size of 0.27 mrad (Nixon et al., 2009b). These focal planes permit high spatial resolution nadir mapping of Titan’s disk, as well as limb profiling of the atmosphere. The data presented here were obtained on orbits #261 (2017-02-17) and #275 (2017-05-24), the latter being Cassini’s final flyby of Titan on which whole-hemisphere CIRS (nadir) mapping was performed.
The CIRS nadir observations were taken with moderate (2.5 cm*-1*) spectral resolution, which allowed rapid mapping of an entire hemisphere of Titan by scanning the FP3/4 pixel arrays in a ‘pushbroom’ pointing sequence. Observations were obtained on orbit #275 at a distance from Titan of 300 000 km, which resulted in a pixel size of 100 km, or 3*∘* of latitude when projected to Titan’s location. FP3 contained the prominent rovibrational emission bands of HCN, HC3N, C2H2, C2H6, C3H4 and C4H2 in the range 628-814 cm*-1* (see Teanby et al. 2009 for example spectra). The interval 650-660 cm*-1* is free from trace gas emission features and was used to determine the continuum level. In the present study, we focus primarily on the HCN and HC3N data for comparison with ALMA observations of these species.
CIRS limb observations were obtained on orbits #261 and #275 with the FP3 and FP4 pixel arrays orientated perpendicular to the limb to give spectra covering tangent altitudes 0-600 km. Observations were in ‘sit-and-stare’ mode at 0.5 cm*-1* spectral resolution, resulting in high signal-to-noise, high-resolution spectra at a single latitude. Orbit #261 covered the latitude (S) with a field-of-view size of 45 km and #275 covered latitude (S) with a field-of-view size of 60 km. On orbit 261 only every other pixel was measured, so this observation had an equivalent vertical resolution of 120 km.
4. Results and Analysis
4.1. Spectra and Maps
ALMA spectra for our main species of interest are shown in the left panels of Figures 1 and 2. These were obtained by integrating the reduced data cubes inside a circular aperture with radius 4300 km, centered on Titan (compared with the solid-body radius of 2575 km). This was sufficiently large to include all detected flux from the most extended (HNC) emission map. The resulting disk-integrated spectrum was too noisy to permit detection of C2H3CN, so the spectrum for this species was extracted from a beam centered at the emission peak for this molecule (200 km outside Titan’s southern polar limb; see Figure 1).
Spectral features greater than above the noise were identified and integrated to produce the emission maps shown in the right panels of Figures 1 and 2. For C2H3CN, only three of the expected nine spectral lines were detected with greater than confidence (as shown in Table 1; two of these lines are blends of two closely-spaced rotational transitions). The unexpected weakness of the C2H3CN emission could be explained by temporal variability with respect to the earlier studies of this molecule by Palmer et al. (2017) and Lai et al. (2017), but may also be due to relatively poor coupling between our ALMA beam and the intrinsic C2H3CN distribution — the earlier studies used a larger beam size (), which was more sensitive to emission spread across Titan’s entire disk plus atmosphere — whereas our beam was only able to recover the compact C2H3CN emission peak near the south pole. To improve the signal-to-noise (S/N) ratio of the C2H3CN map in Figure 1, integration was performed over the frequency ranges corresponding to the nine expected strongest lines in one of our ALMA spectral windows (360.8-361.7 GHz), using the model of Palmer et al. (2017) as a guide for the line strengths, and adopting a 1.2 km integration width for each spectral line.
The resulting ALMA maps provide the highest-resolution (instantaneous) views of HNC, CH3CN, C2H3CN, and C2H5CN in Titan’s entire sunward-facing-hemisphere published to-date using any instrument. Characteristic emission patterns are apparent for each species, and provide for the first time, a detailed view of the global column density distributions for these gases. Since most of the detected emission lines are relatively weak (and unsaturated), the majority of the observed molecular emission is expected to be optically thin (a notable exception is HC3N, which becomes optically thick in the vicinity of the south pole; Cordiner et al. 2018). Titan’s sub-mm continuum forms near the tropopause (at altitudes km; Thelen et al. 2018), where the temperature is close to 70 K. Our species of interest are synthesized primarily at higher altitudes (above a few hundred km; e.g. Vuitton et al. 2019), where -180 K, and are subsequently transported to the lower stratosphere where they condense out (at km; ). Across Titan’s disk (and limb), the targeted spectral lines are observed in emission because these gases are much warmer than the 70 K background continuum source of Titan (and the 2.73 K cosmic microwave background). Furthermore, while Titan’s atmospheric temperatures are known to vary by up to K with latitude (Achterberg et al., 2011), our observed lines are only weakly sensitive to temperature variations in this range (see e.g. Lai et al., 2017; Thelen et al., 2019). These ALMA emission maps therefore provide a good indication of the intrinsic column density distribution of each species.
The HC3N, CH3CN, C2H3CN, C2H5CN, HCN and HNC maps obtained using ALMA all show strong variations in column density across Titan’s disk. Limb brightening dominates the HCN, HNC and CH3CN maps, characteristic of a relatively uniform latitudinal distribution, whereas the HC3N, C2H3CN and C2H5CN maps are dominated by polar emission, with the southern (winter) pole being strongest for the latter three species.
The HCN and HC3N radiance maps obtained using Cassini CIRS (on orbit #275) are shown in Figure 3. While the spatial resolution is significantly better than obtained by ALMA, the basic morphology of the resulting maps is similar, with a dominant south-polar emission peak for HC3N (at latitudes ) and a more uniform, limb-brightened distribution for HCN. The intensities of these vibronic emission bands observed by CIRS have a stronger dependence on temperature than the pure rotational transitions observed by ALMA. As a result, the relationship between CIRS observed radiances and molecular column densities is non-trivial, making the interpretation of these maps more difficult. For example, it is likely that the drop in HCN radiance with decreasing latitude (south of the equator) is a consequence of falling stratospheric temperatures towards the winter pole (see e.g. Teanby et al., 2017, 2019).
A unique aspect of our ALMA observations compared with typical (ground or space-based) infrared measurements is the extremely high spectral resolution (up to ), which enables the detection of emission lines from higher in the atmosphere ( km). This is due to the relative lack of pressure broadening above km, resulting in the concentration of emission into narrow (thermally-broadened) spectral line peaks. This benefit comes with the tradeoff that we are less sensitive to low-altitude emission from the pressure-broadened line wings (below km), although for the stronger lines (such as the pure rotational transitions of HC3N and CH3CN), we maintain good sensitivity down to the lower stratosphere (around km; see e.g. Thelen et al. 2019). While HCN and HC3N have previously been extensively mapped at higher resolution by Cassini CIRS, instantaneous, whole-hemisphere limb coverage has not been possible at high resolution using that instrument; instead, CIRS limb maps may be built up by assembling data from multiple Cassini flybys over a period of months or years (during which time the atmosphere may have changed). Furthermore, CIRS nadir mapping is most sensitive to emission from altitudes around 150 km; our ALMA maps thus provide unprecedented spatial information at higher altitudes, for the 2017 May (solstice) epoch.
While the main focus of our present study is on the nitrile observations, we also obtained a serendipitous detection of the C3H8 multiplet at 360.978 GHz (comprised of a blend of eight transitions with ). This is the first time propane has been definitively detected and mapped in any extraterrestrial source at radio wavelengths (following a tentative ALMA detection on Titan by Lai et al. 2017), and is therefore worthy of further investigation. While C3H8 has an extremely weak radio spectrum due to its small dipole moment and large partition function, the large abundance of this molecule in Titan’s atmosphere enables its detection using ALMA. Vibrational emission from C3H8 was previously detected on Titan by Voyager 1 (Maguire et al., 1981), ISO (Coustenis et al., 2003) and Cassini (Nixon et al., 2009a). Our ALMA map shows a relatively uniform distribution for C3H8 across Titan’s disk, with little evidence for any latitudinal variability above the (admittedly high) noise level.
Maps and spectra of the HC3N vibrationally excited lines and isotopologues have already been presented by Cordiner et al. (2018), and are consistent with the overall distribution of the HC3N line in Figure 1, but provide a more accurate representation of the HC3N column density distribution (albeit at lower S/N), due to their reduced optical depth (see Section 4.2). We also obtained high-sensitivity CO and HCN spectra (including the main line of HCN as well as its 13C and 15N isotopologues). Unfortunately, these lines were too optically thick to provide useful emission maps, and require detailed radiative transfer modeling for their analysis, which will be presented in a future article. We use the vibrationally excited HCN () line as an alternative to mapping the ground-state () line of this molecule, the core of which is completely optically thick (self-absorbed) across much of Titan’s disk. The HCN line is predominantly optically thin, but sits atop the strong, pressure-broadened wing of the (ground state) HCN line wing, which was subtracted using a polynomial fit, thus allowing the line to be mapped.
4.2. Latitudinal Profiles
The latitudinal variability observed in Figures 1 and 2 can be quantified by plotting the flux as a function of (sky-projected) angle about Titan’s disk-center. To eliminate the impact of the differing (elliptical) beam dimensions as a function of angle, it was first necessary to degrade the resolution of our images so as to produce a circular PSF (with diameter equal to that of the long-axis of the original beam ellipse). The average flux was then obtained for each species inside radial wedges (of angular size ), originating from the center of Titan’s disk, and plotted as a function of angle from the perpendicular bisector of Titan’s polar axis in Figure 4. A wedge size was used for all molecules apart from C2H3CN and C3H8, for which was used, to improve the S/N for these species. The angular coordinate system and binning scheme are shown in the top-right panel for CH3CN. The relationship between the angular range covered by each wedge and the corresponding latitude on Titan is complicated due to the finite wedge size and the tilt of Titan’s pole. This results in complete obscuration of the south polar region, while at the north pole, a relatively large range of latitudes are included in each wedge. The weighted mean latitude was taken within each wedge, using tangent latitudes for pixels falling outside Titan’s limb. These weighted mean latitudes are shown on the upper abscissa axes in Figure 4.
The isomeric partners HNC and HCN are shown on the same plot (with HNC multiplied by a factor of 5 for display). The profile for vibrationally-excited HC3N () is shown in addition to the ground-state () transition for that molecule (multiplied by 3 for display). Cordiner et al. (2018) determined that the ground-state HC3N rotational lines in ALMA Band 7 become highly saturated at the south pole, which makes the transition unreliable as a measure of the HC3N column density at that location. The transition, by contrast, is reasonably optically thin and thus provides an improved measure of the true HC3N distribution.
Significant latitudinal variations are apparent for HCN, HC3N, CH3CN and C2H5CN. The strongest variations are for HC3N and C2H5CN, which peak in the vicinity of the south pole and fall sharply towards the equator (over a distance smaller than the spatial resolution element). This behaviour is consistent with high concentrations of these gases within regions km from the pole. HC3N and C2H5CN show a weaker peak over the north pole. By contrast, CH3CN peaks most strongly in the north, whereas the HCN northern and southern abundance peaks are of similar strength. C3H8 shows no evidence for significant latitudinal variation.
4.3. Longitudinal Variability of HNC
The only species to show significant longitudinal (east-west) variability is HNC, manifested by a relative weakness in the HNC emission on the western (dawn) limb. Figure 5 (top panel) shows the average flux (in radial wedges), as a function of angle () about the disk-center, measured from the sky-projected north polar axis. The HNC emission is strongest at the north pole, relatively uniform in the eastern hemisphere, then begins to fall towards the south pole (from around ), reaching a trough in the west (), before rising back to the northern peak. The longitudinal asymmetry of HNC is demonstrated more clearly in the bottom panel of Figure 5, where the east-west flux difference ((East-West)) is plotted as a function of angle from the bisector of Titan’s sky-projected polar axis. Corresponding weighted-mean latitudes (shown on the upper abscissa), were calculated as in Section 4.2. The east-west flux difference is seen to reach a peak close to (or slightly south of) the equator. The southerly bias could be a result of viewing geometry, due to the tilt of Titan’s polar axis.
Such longitudinal asymmetry is surprising given the presence of fast, superrotating zonal winds (Lellouch et al., 2019), which traverse Titan’s globe in a period of less than a few Earth days, and are therefore expected to quickly even out atmospheric inhomogeneities as a function of longitude. While Cordiner et al. (2014) claimed the first evidence for longitudinal variations in Titan’s HNC, their measurements were from relatively low resolution, low-sensitivity flux calibration data, and the observed east-west asymmetries could plausibly have resulted from the combined effects of noise and smearing by a strongly elliptical telescope beam (M. Gurwell, private communication 2016). The HNC asymmetries presented here are clearer and more robust than those identified by Cordiner et al. (2014). We measure a total HNC flux of mJy km from Titan’s eastern hemisphere and mJy km from the western hemisphere, which amounts to a ratio %. Thus, a conclusive, difference is identified in the HNC flux between the hemispheres. Use of a circularized beam renders this result free from any geometrical artifacts related to (elliptical) beam smearing.
4.4. ALMA Limb Abundance Retrievals
The variation in integrated line flux () as a function of sky-projected radial distance from Titan’s limb () contains information on how the gas density () varies as a function of altitude (). The vertical abundance profile can thus be retrieved for each gas by fitting an atmospheric model to the ALMA observations. Despite coarse radial sampling due to limited angular resolution, our measurements can nevertheless provide unique information on the high-altitude abundances for our species of interest. To generate at the highest resolution (at the expense of some sensitivity), Briggs weighting of the interferometric data was used (with a parameter of zero), leading to a resolution of ( km) at 362 GHz.
We generated profiles by remapping the image pixels onto a polar grid with origin at the center of Titan’s disk. Fluxes as a function of radius were extracted and averaged within angular regions -135*∘* and -315*∘*. These ranges were chosen in order to capture a sufficiently large range of equatorial latitudes for maximum S/N while excluding the more complex polar regions (seen for HC3N, CH3CN, C2H3CN and C2H5CN in Figure 1). Abundance enhancements at the winter pole are associated with subsidence within a highly localized region inside Titan’s polar vortex (Teanby et al., 2008; Vinatier et al., 2015; Teanby et al., 2017). This leads to downwelling of photochemically-enriched gases, resulting in a distinctly different abundance (and temperature) profile at the winter pole compared with equatorial regions. Our simplistic, 1D modeling procedure provides an estimate of the mean abundance profile at equatorial latitudes (where the abundances are dominated by high-altitude chemistry), to the exclusion of the polar regions, which are strongly influenced by seasonally-variable gas dynamics.
The peak-normalized profiles are shown in Figure 6 (top panel). Clear separations are present between the peak emission radii () for each species, which were derived from the centroids of the profiles, and are given in Table 2. While the C3H8 and CH3CN fluxes peak at km and 2590 km, respectively (close to the edge of Titan’s solid km disk), C2H5CN, HCN, HC3N, and HNC peak well outside this radius, indicating flux contributions from high altitudes — the very high HNC peak emission radius of 3140 km is particularly surprising, and is consistent with a dominant thermospheric/ionospheric flux contribution for this gas. Unfortunately, an profile (and corresponding value) could not be generated for C2H3CN due to a lack of flux from this molecule near the equator.
To further interpret the observed limb profiles, it is necessary to account for the line-of-sight averaging of different altitudes within each ALMA beam, as well as the atmospheric density (and temperature) variation with altitude. This was done by constructing a simple, spherically-symmetric molecular emission model for Titan. The density profile was taken from Krasnopolsky (2009), and temperatures were retrieved from CIRS CH4 limb observations at latitude (observed during rev#275), using the method of Teanby et al. 2007. Molecular abundances as a function of altitude () followed the simplest polygonal scheme possible while preserving a good fit to the observations: a logarithmic abundance gradient was adopted between altitudes and , with zero abundance below and constant abundance above . Spectral line fluxes were generated using parameters from the Cologne Database for Molecular Spectroscopy (Müller et al., 2001), under the assumption of optically thin emission. To reproduce the ALMA spatial response, the resulting images were passed through the ALMA simulator (using the Simobserve and Clean tasks), with the same observational and imaging parameters as described in Section 2. Synthetic image cubes were integrated over the same spectral ranges as the observations (Figures 1 and 2), and radial flux profiles were extracted for comparison with observation.
The best-fitting , and abundance gradient parameters were obtained through nonlinear least-squares optimization using the MPFIT routine (Markwardt, 2012), and parameter errors were determined from the covariance matrix. The corresponding optimized limb flux models are shown for comparison with observations in Appendix A (Figure 9), and retrieved abundance profiles are shown in Figure 6 (lower panel). Given the simplicity of our abundance parameterization scheme, the model limb flux profiles provide a surprisingly good fit for all species, with respect to the observational noise errors.
To determine the altitude sensitivity of our abundance models, contribution functions (total flux as a function of altitude) were calculated within an ALMA beam centered at the respective value for each species. The retrieved abundances are considered to be most reliable around the flux-weighted mean altitude () of the contribution function, and these values are given in Table 2. Pressure-broadened line wings become significant in models for HC3N and C2H5CN sub-mm emission at altitudes below km (Cordiner et al., 2015; Thelen et al., 2019), but such line wings are absent (or below the noise threshold) in our ALMA limb spectra, so the retrievals for these molecules contain little useful information below this altitude. We therefore truncate the abundance profiles for these species below 300 km in Figure 6. We also truncate the HCN profile below 300 km due to the difficulty in reliably separating the wings of the line from those of the ground-state HCN () line. This region of the HCN () line wing originates predominantly from altitudes 100-200 km (Molter et al., 2016), and is strong enough to call into question the validity of the optically thin approximation for this molecule, which may therefore have also resulted in underestimating the overall HCN abundances. By virtue of its intrinsic line strength and high abundance at lower altitudes, pressure-broadened wings were clearly detected for CH3CN, allowing the abundance for this species to be retrieved down to km (determined using the CH3CN contribution function of Thelen et al. 2019).
Our ALMA limb observations provide the first well-constrained abundance retrievals for HNC and C2H5CN in the upper atmosphere. Prior (disk-averaged) observations of these molecules by Moreno et al. (2011) and Cordiner et al. (2015) were only able to place lower limits (of 400 km and 200 km, respectively), on the altitudes from which the emission from these species originated. Our limb profile modeling shows that the detected HNC emission is confined to thermospheric altitudes above km, whereas C2H5CN follows a steep (logarithmic) abundance gradient from at 340 km to at 1200 km. While the present article was undergoing peer review, the study of Lellouch et al. (2019) was published, also confirming that HNC is restricted to the thermosphere.
5. Comparison with Cassini CIRS
Accounting for variations in atmospheric temperature and instrumental resolution, the primary morphology for HCN and HC3N observed by CIRS is consistent with that revealed by our ALMA observations. Both Cassini and ALMA reveal that most of the HCN emission is concentrated in an atmospheric band (relatively uniform with latitude), at altitudes between 50-400 km. The HC3N emission is dominated by a higher-altitude, southern polar peak, confined to latitudes south of . The northern HC3N peak detected by ALMA is barely visible in the Cassini nadir observations, which is consistent with the north-to-south column density ratio of determined from the ALMA lines. As explained in Section 4.1, CIRS HCN radiances have a significant temperature dependency, which can explain the lack of a southern polar emission peak — in fact, the CIRS HCN flux reaches a minimum at the south pole, where the ALMA flux reaches a maximum.
To enable a more direct comparison of the ALMA and Cassini observations, vertical profiles of temperature, HCN, and HC3N abundance were retrieved from the CIRS data by modeling the limb radiances. The NEMESIS spectral inversion tool (Irwin et al., 2008) was used to fit the limb observations from orbit #261, following the method of Teanby et al. (2007). Prior to inverting for composition, a cubic b-spline curve was fitted to the measured radiance profiles at each wavenumber using a knot spacing of 100 km (see Teanby 2007 for further details). This improved the signal to noise by smoothing the radiance through multiple pixel measurements at similar but not identical tangent altitudes. The profiles obtained are reliable for 180-400 km; below 180 km the line opacity was too high to allow an accurate determination, and above 400 km the signal-to-noise was too low.
Retrieved equatorial limb abundances for HCN and HC3N are shown in Figure 7. By comparison with Figure 6, these are seen to be in reasonable agreement with the ALMA limb abundance profiles, especially when considering the error bars on the ALMA and CIRS retrievals, as well as the (unquantified) uncertainties in the ALMA limb abundances due to the low spatial resolution and simplified (polygonal) parameterization scheme used for those data.
6. Discussion
6.1. Chemistry and dynamics of the thermosphere/ionosphere probed by HNC
Our HNC observations reveal that this molecule is confined predominantly to the thermosphere, with an abundance of above 1030 km. The retrieved abundance profile (Figure 6) is remarkably similar to the photochemical model results of Vuitton et al. (2019) (shown in Figure 8). The HNC model abundance profile rises steeply in the range 600-1000 km, and reaches a peak abundance in the range 1050-1100 km. The very good agreement between model and observation for this species provides a strong validation for the chemical reaction scheme of Vuitton et al. (2019), in which HNC is produced by electron recombination of the abundant HCNH+ ion, and is rapidly converted back to the lower-energy HCN isomer (on a timescale of about a day), primarily by collisions with ionospheric H atoms. The short chemical lifetime means that HNC is destroyed before it can diffuse down through the atmosphere from its high-altitude production site, and is therefore confined to the thermosphere, with negligible abundance at lower altitudes. The rapid destruction of HNC makes it a potentially unique probe of very short-timescale dynamical processes in Titan’s upper atmosphere.
The observed asymmetry in HNC emission between Titan’s eastern and western hemispheres can be explained as a result of the combined effects of photochemistry and atmospheric dynamics. The presence of mesospheric zonal winds (in the altitude range 250-450 km) was initially inferred on Titan based on the atmospheric structure derived from stellar occultation observations by Hubbard et al. (1993). High-altitude (450-1000 km) superrotating winds were subsequently detected by Moreno et al. (2005) using ground-based interferometry. Cui et al. (2009) measured a strong (order of magnitude) diurnal variation in Titan’s ionospheric HCNH+ abundance using Cassini mass spectrometry, with the ion strongly depleted in the night hemisphere. They constructed a chemical/dynamical model for the ionosphere, which showed that east-west asymmetry in the distribution of HCNH+ (and other ions) can be produced by a zonal wind that sweeps ion-rich material from the day side towards the eastern (dusk) terminator. The rapid conversion of HCNH+ to HNC via electron recombination would then lead to an enhanced HNC abundance at the eastern limb. As it passes across the night side, collisional isomerization converts HNC back to HCN, resulting in subsequent emergence of HNC-depleted air in the west. The detection by Lellouch et al. (2019) of superrotating zonal HNC winds extending up to the ionosphere lends support to this theory, so in situ production of an asymmetric HCNH+ distribution by electrons from Saturn’s magnetosphere (e.g. Cravens et al., 2009) may not be required.
The HNC line flux has a relatively weak dependence on temperature, so a relatively large, 30-40 K hemispheric difference in rotational temperature would be required to explain the observed asymmetry. The atmospheric circulation model of Müller-Wodarg et al. (2000) predicted a diurnal temperature variation of 10-20 K (the precise value depends on solar activity), at an altitude of 1300 km, corresponding to the top of Titan’s thermosphere. However, this variation was found to fall with decreasing altitude, becoming negligible at the km altitude to which our HNC observations are most sensitive. Combined with the relatively low Solar activity around the time of our observations, diurnal temperature variations seem an unlikely explanation for the observed HNC asymmetry.
6.2. High altitude nitrile chemistry
6.2.1 Inferred molecular lifetimes
The vertical abundance profiles retrieved for the remaining nitriles (HCN, CH3CN, HC3N and C2H5CN) are broadly consistent (within about an order of magnitude) with the results of recent chemical models (e.g. Krasnopolsky, 2009; Loison et al., 2015; Willacy et al., 2016; Dobrijevic et al., 2016; Vuitton et al., 2019), confirming the efficient production of organics through photochemistry in the thermosphere (above 600 km), initiated by Solar ultraviolet radiation. An abundance gradient develops as photochemical products mix downward towards the lower stratosphere, where the temperature becomes low enough for them to precipitate out. The steepness of the abundance gradient is affected by the chemical lifetime of each species, so from Figure 6, we may infer that after HNC, HC3N has the next shortest lifetime, followed by HCN, C2H5CN then CH3CN. The overall trend is compatible with the lifetimes predicted at 300 km by Vuitton et al. (2019) (see Table 2), although it should be noted that their base model overestimated the stratospheric C2H5CN abundance, so the lifetime for this species is taken from their alternative model where C2H5CN is lost from the gas-phase via sticking to stratospheric aerosol particles. The steep abundance profile for HCN also seems to be inconsistent with a long (350 yr) lifetime at 300 km. However, the HCN model lifetime varies strongly with altitude, and takes a value of only a few years in the lower thermosphere.
As a result of the coupling between the observed gas spatial distributions and Titan’s seasonally varying global circulation, approximate lifetimes may also be derived from our ALMA emission maps, based on the rate of decay of abundance at the summer pole since the time of the last seasonal reversal. We begin with the assumption that sufficient time has elapsed since the last seasonal reversal (following the 2009 equinox) that the molecular abundances at the southern (winter) pole are in a quasi steady-state, maintained by a balance between the source (Titan’s global meridional circulation system) and sink (photochemical destruction). Then, neglecting temperature effects (which are expected to have only a small impact on the measured fluxes; see Section 4.1), the peak of the spectrally integrated flux in the southern hemisphere () for a given species at the present epoch is taken to be representative of the peak flux that would have been present in the north at the time of last circulation reversal. This is supported by Teanby et al. (2019) who suggest that north and south polar stratospheric abundances are closely comparable at similar seasonal phases, based on trends observed in CIRS nadir data across the entire Cassini mission. Using temporally-resolved molecular observations from CIRS, Vinatier et al. (2015) determined that the last circulation reversal occurred in 2011, so the peak north polar flux () has decayed to reach its present value over a period of 6 years.
Assuming exponential decay of the form , the following lifetimes were derived — HC3N: 3.4 yr, C2H3CN: yr, C2H5CN: 9.3 yr. These values are in reasonable agreement with the chemical/dynamical lifetimes () from the model of Vuitton et al. (2019) (see Table 2). We measure for CH3CN, which implies that the south polar abundance for this species has not yet reached steady state, so its lifetime cannot be derived using our method. However, the fact that the north and south polar fluxes differ indicates that the lifetime for this species is probably less than (or similar to) the 29.5 yr Saturnian year; if it were much longer, would be expected as a result of the repeated replenishment of chemically-enriched polar gases following successive winters. For HCN we find , which implies a relatively long lifetime for the observed (polar) gases, consistent with the larger values of 60 yr and 350 yr at 300 km calculated by Dobrijevic et al. (2016) and Vuitton et al. (2019), respectively.
These values represent an averaged measure of the molecular lifetime due to the combined influx/loss rates of molecules to/from the north polar region as a result of photolysis, chemical reactions, precipitation, diffusion, advection and winds, and may therefore not be considered as directly comparable with values from 1-dimensional chemical models. Additional caution is required for their interpretation due to the gross assumptions made in their derivation. In particular, the south polar flux at the present epoch may not be precisely representative of the northern value at the time of the 2011 circulation reversal, which may be the case if steady state was not reached in time, or in the case of seasonal climate asymmetry (i.e., if the atmosphere in Titan’s northern winter is not a mirror image of the southern winter). Furthermore, the tilt of Titan’s polar axis means that the south polar region was partially obscured from view at the time of our observations, so that provides an underestimate of the total flux. However, the majority of our observed flux is from altitudes above 200 km and latitudes , all of which is within view, so geometrical obscuration is not expected to be a dominant source of uncertainty.
6.2.2 Comparisons with prior observations and models
Our ALMA limb abundance profiles for HCN, HC3N and CH3CN are in general, closely consistent with the prior equatorial retrievals from Thelen et al. (2019), obtained by modeling ALMA spectra at lower spatial resolution. Our C2H5CN profile between 300 and 600 km is quite consistent with the disk-averaged gradient profile from Cordiner et al. (2015). The thermospheric HCN profile also matches reasonably well with the Cassini VIMS limb observations of Adriani et al. (2011). It should be noted that the vertical abundance retrievals in our present study (based on spectrally-integrated limb fluxes) contain only limited information due to the relatively coarse ( km) spatial resolution of our data, and correspondingly simplified (polygonal) parameterization scheme for the abundances as a function of altitude. As such, our ALMA retrievals should be considered only approximate; full retrievals taking into account the detailed (pressure broadened) spectral line profiles will be presented in a future article. Nevertheless, some useful comparisons between our results and prior observations and models can still be made, particularly at higher altitudes (above 500 km) where previous (IR and radio) retrievals contain little useful information.
The retrieved HCN, HNC, HC3N, C2H5CN and C3H8 abundances at an altitude of 1000 km are close (within errors) to those predicted by the models of Loison et al. (2015), Willacy et al. (2016) and Vuitton et al. (2019). Predicted vertical abundance profiles from Vuitton et al. (2019) are shown for comparison with our observations in Figure 8. For CH3CN, these models over-estimate the ionospheric abundance by about an order of magnitude, which may indicate the presence of unaccounted-for (or underestimated) destruction pathway(s) for this molecule. It should be noted that although our HCN profile matches the model of Vuitton et al. (2019) reasonably well in the mesosphere and above, our retrieved HCN abundances are more uncertain towards lower altitudes due to the difficulty in separating the wings of the line from those of the main () line.
ALMA equatorial limb data show that HC3N maintains a high abundance to lower altitudes than predicted by the Loison et al. (2015) and Dobrijevic et al. (2016) models, which fall too rapidly in the lower thermosphere (between 800-600 km). The models of Willacy et al. (2016) and Vuitton et al. (2019) provide a better overall match for our HC3N observations. Note, however, there is still some considerable uncertainty regarding the rate of loss of HC3N via sticking to aerosol particles in the Vuitton et al. (2019) model, so additional studies are needed to fully understand the chemistry of this species.
In summary, the reasonably good agreement between our nitrile observations and photochemical models for Titan’s mesosphere and thermosphere indicates that the fundamental processes involved in the high-altitude production of nitriles are quite well understood. However, towards the stratosphere, the base model of Vuitton et al. (2019) does a relatively poorer job of matching the HC3N and C2H5CN observational data, which are over-predicted by at least an order of magnitude (a similar problem occurs in the lower stratosphere for the Loison et al. 2015 and Dobrijevic et al. 2016 models), indicating that important loss mechanisms for these molecules may be missing from chemical networks. This issue may be addressed with the inclusion of sticking of these molecules to aerosol particles as an additional loss process in the stratosphere.
7. Conclusions
The ALMA observations presented here constitute a comprehensive, moderate-resolution mapping of Titan’s most abundant atmospheric nitriles around the time of the 2017 solstice, coinciding with the end of the Cassini mission. This is the first time detailed, whole-hemisphere maps have been published for the complete set of nitriles accessible in the microwave/sub-mm band, and constitutes a unique resource for investigating their global distributions. Each species shows a characteristic spatial distribution, resulting from the interplay between photochemical production, destruction and atmospheric dynamics. ALMA’s unique capability to probe narrow molecular emission lines from high altitudes (-1000 km) in the mesosphere and thermosphere, in addition to pressure-broadened emission from lower, stratospheric altitudes (-300 km) provides complimentary data to that obtainable using infrared instruments such as Cassini CIRS and IRTF TEXES.
Our global maps for HC3N, C2H3CN and C2H5CN reveal strong, compact emission peaks over the southern (winter) pole, indicative of rapid photochemical production (and destruction) at high altitudes, combined with transport towards southern latitudes by Titan’s global meridional circulation system. Approximate lifetimes for these species are on the order of a few years, which matches the stratospheric chemical/dynamical lifetimes (at km) predicted by the latest chemical models for Titan’s atmosphere, indicating that their dominant production and loss mechanisms are reasonably well understood. The HCN and CH3CN maps show polar peaks with reduced contrast compared with the other species, implying that they are relatively more well-mixed in latitude, consistent with longer lifetimes for these gases (between about half and a few Saturnian years).
Limb flux profiles confirm that CH3CN tends to be one of the more abundant nitriles at lower (stratospheric) altitudes, whereas the HC3N abundance in the stratosphere is low, and rises sharply with altitude throughout the mesosphere. The HNC vertical abundance profile is unique among the gases measured, as it does not become detectable until altitudes km in Titan’s thermosphere. Thus, HNC provides a rare opportunity to probe chemical and physical processes occurring in Titan’s thermosphere/ionosphere, using ground-based observations.
The ALMA HNC map shows significant longitudinal asymmetry. This is explained as a result of the combined effects of zonal winds and diurnal variations in the ionospheric HCNH+ abundance, which is theorized to be lower on the night side due to reduced solar insolation, giving rise to a reduced HNC abundance at the western (dawn) limb. The observed east-west asymmetry confirms the very rapid theorized production and loss rates for this species (on the order of a day at 1000 km altitude), making HNC a unique probe of short-timescale ionospheric processes.
This work was supported by the National Science Foundation under Grant No. AST-1616306, the NASA Astrobiology Institute (NAI), the NASA Solar System Observations program, and the UK Science and Technology Facilities Council (STFC). It makes use of ALMA data set ADS/JAO.ALMA#2016.A.00014.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.
Appendix A Limb flux profile fits
Least-squares fits to the equatorial fluxes as a function of (sky-projected) radius are shown in Figure 9, using a simple, four-parameter model (two parameters for C3H8) for the abundance of each species () as a function of altitude (). See Section 4.4 for details of the modeling procedure. The corresponding optimized profiles are shown in Figure 6.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Achterberg et al. (2011) Achterberg, R. K., Gierasch, P. J., Conrath, B. J., Michael Flasar, F., & Nixon, C. A. 2011, Icarus, 211, 686
- 2Adriani et al. (2011) Adriani, A., Dinelli, B. M., López-Puertas, M., et al. 2011, Icarus, 214, 584
- 3Bézard et al. (2014) Bézard, B., Yelle, R., & Nixon, C. 2014, Titan: Surface, Atmosphere and Magnetosphere, Ed.s I. Muller-Wodarg, C. Griffith, E. Lellouch and T. Cravens, Chap. 6, Cambridge University Press, 158
- 4Boucher et al. (1980) Boucher, D., Burie, J., Bauer, A., Dubrulle, A., & Demaison, J. 1980, JPCRD, 9, 659
- 5Brauer et al. (2009) Brauer, C. S., Pearson, J. C., Drouin, B. J., & Yu, S. 2009, Ap JS, 184, 133
- 6Cordiner et al. (2018) Cordiner, M. A., Nixon, C. A., Charnley, S. B., et al. 2018, Ap J, 859, L 15
- 7Cordiner et al. (2014) Cordiner, M. A., Nixon, C. A., Teanby, N. A., et al. 2014, Ap J, 795, L 30
- 8Cordiner et al. (2015) Cordiner, M. A., Palmer, M. Y., Nixon, C. A., et al. 2015, Ap J, 800, L 14
