Decomposing Titan’s far-infrared Haze B feature with PCA and NMF analysis
(2026)
Authors:
Joshua S Ford, Nicholas A Teanby, Patrick GJ Irwin, Conor A Nixon, Lucy Wright
Abstract:
During its flyby of Titan in 1980, Voyager 1 unveiled an atmosphere thick with an opaque, orange haze that completely obscured the surface, and provide the first close-up images of its vertical layering (Smith et al. 1981, Hanel et al 1981). Subsequent observations from the Cassini-Huygens mission (Flasar et al. 2004) expanded this view, revealing a world dominated by photochemical hydrocarbons and nitriles with lakes of methane (Stofan et al. 2007, Mastrogiuseppe et al. 2019), tholin-like organic sand (Lorenz et al. 2006), and seasonally evolving ice clouds (Jennings et al. 2012, West et al. 2016) as seen in Figure 1. While recent observations and laboratory studies (Anderson et al. 2011, Chatain et al. 2020) have provided greater insights into the physical properties of these aerosols and how they form, their composition and complexity remain a mystery.Figure 1: Image taken by Cassini Imaging Science Subsystem (Credit: NASA/JPL-Caltech/Space Science Institute) showing a large ice cloud possibly made of HCN at 300km in Titan’s south pole, 2012 (West et. 2016, Vinatier et al. 2018)In the upper atmosphere, nitrogen and methane are photodissociated by UV radiation and energetic particles to produce ions which act as embryos for the growth of large organic molecules (Vuitton et al. 2024). As these molecules sink through the atmosphere, they recombine, coagulate, and accumulate into larger aerosol particles, onto which trace species can condense, forming both photochemical hazes and stratospheric ice clouds (Vuitton et al. 2024). Cassini CIRS (Composite Infrared Spectrometer) far infra-red spectra (FP1) exhibit four types of hazes: Haze 0, Haze A, Haze B (otherwise known as the “Haystack”) and Haze C (de Kok et al. 2007). While Haze 0 is present throughout the CIRS spectral range and affects the spectral continuum, Haze A, B and C appear as broad features in the far-infrared range. Relatively little is known about these hazes, making their features difficult to fit and estimate. In addition, their signatures overlap with H2O, C4H2, C2N2 and CH4 rotational lines, complicating retrievals of these gases (Sylvestre et al.2017). Among these hazes, Haze B exhibits the strongest spectral signature. Anderson et al. 2018 suggested that the composition must be a mixture of more than one chemical compound due to its magnitude, breadth and opacity. Haze B is observed to only be present at the winter poles, forming and dissipating with the changing seasons, most likely caused by reduced sunlight and temperature (Jennings et al. 2012a,b , Anderson et al 2018).To enable accurate fitting of water features, Ford et al. (in review) retrieved an effective spectral cross-section of Cassini CIRS FP1 spectra between 147-257cm-1 of 156 FIRNADCMP 0.5cm-1 observations (Figure 2) across different latitudes and times (see also Ford et al. 2025). This was achieved by scaling gaussian basis functions to fit Haze A, B and C (and any other unknown aerosols) in the spectra baselines using the NEMESIS radiative transfer code (Irwin et al. 2008). Previous aerosol cross-sections did not account for latitude or time variation and therefore this technique empirically and agnostically modelled the changing baseline. In this study, we decomposed those 156 effective spectral cross-sections and extrapolated unique Haze B cross-sections for latitudes of -89° to 88° from June 2004 to April 2017. We used principal component analysis to determine the number of complete hazes present in the spectral range, finding that Haze B accounts for 98% of the observed variation. We then applied non-negative matrix factorisation to seperate the spectra into two unique and stable components: one representing Haze B, and a second representing a mixture of other haze contributions, along with numerical and spectral noise.Figure 2: Effective spectral cross-section of all 156 Cassini CIRS FIRNADCMP 0.5cm-1 observations before decomposition. Each colour represents a different observation. The plot shows the variation extent of the Haze B feature at ~220cm-1. Using the maximum Haze B cross-section of each observation as proxy, we find a large increase at the winter poles with the increase during northern spring being nearly twice that of northern winter, consistent with previous investigations. We also find Haze B extends to 60°, different to 70° proposed by Anderson et al. 2018. The results also give an insight into the time of formation/dispersion of Haze B at the poles. We see that Haze B at the south pole forms around 2013, reaching its peak during the summer solstice. The final dissipation of Haze B at the north pole cannot be determined due to limited data. ReferencesAnderson, C.M. et al. (2011), Icarus, 212.2, 762-778. DOI: 10.1016/j.pss.2010.10.009Anderson, C.M. et al. (2018), Organic Ices in Titan’s Stratosphere in Space Science Reviews, 214.8, 125 DOI: 10.1007/s11214-018-0559-5Chatain, A. et al. (2020), Icarus, 345, 113741. DOI: 10.1016/j.icarus/2020.113741Coustenis, A. et al. (1999), Planet Space Science, 47, 1305-1329. DOI: 10.1016/S0032-0633(99)00053-7de Kok, R. et al. (2007), Icarus, 191, 223. DOI:10.1016/j.icarus.2007.04.003Flasar, F.M. et al. (2004), Space Science Reviews, 115, 169–297. DOI: 10.1007/s11214-004-1454-9Ford, J.S. et al. (2025), EGU General Assembly 2025, Vienna, Austria, EGU25-3741. DOI: 10.5194/egusphere-egu25-3741Ford, J.S. et al. In review. Titan’s Stratospheric Water: Latitudinal and Seasonal Variation from Cassini CIRS Data and Implications for External Oxygen Sources. PSJHanel, R.A. et al. (1981), Science, 212, 192–200. DOI: 10.1126/science.212.4491.192Irwin, P.G.J et al. (2008), JQSRT, 109, 1136. DOI: 10.1016/j.jqsrt.2007.11.006Jennings, D.E. et al. (2012), The Astrophysical Journal Letters, 754, L3. DOI: 10.1088/2041-8205/754/1/L3Khanna, R.K. (2005), Icarus, 178,165-170. doi:10.1016/j.icarus.2005.03.011Lorenz, R. et al. (2006), Science, 312, 724–727. DOI: 10.1126/science.1123257Mastrogiuseppe, M. et al. (2019), Nature Astronomy, 3, 535, DOI: 10.1038/s41550-019-0714-2Samuelson, R.E. Clouds and aerosols of Titan’s atmosphere, in The Atmospheres of Saturn and Titan. Proc.Int. Workshop, vol. ESA SP-241 (1985)Smith, B.A. et al. (1981), Science, 212, 163–191. DOI: 10.1126/science.212.4491.163Stofan, E.R. et al. (2007), Nature, 445, 61–64. DOI: 10.1038/nature05438Sylvestre, M. et al. (2017), A&A, 609:A64. DOI: 10.1051/0004-6361/201630255Vinatier, S. et al. (2018), Icarus, 310, 89-104. DOI: 10.1016/j.icarus.2017.12.040Vuitton, V. et al. (2024), Chapter 6:Titan’s Atmospheric Structure, Composition, Haze, and Dynamics in Titan after Cassini–Huygens, COSPAR Scientific Symposium Series.West, R.A. et al, (2016), Icarus, 270, 399-408. DOI: 10.1016/j.icarus.2014.11.038
