Neptune's Latitudinal H2S Distribution: Reconciling Near-Infrared and Microwave Observations

Copernicus Publications (2025)

Authors:

Joseph Penn, Patrick Irwin, Jack Dobinson

Abstract:

In 2018, analysis of Gemini-NIFS near-infrared observations revealed the probable presence of H2S above the main cloud deck on Neptune [1]. The spectral signature of the gas was found to be much stronger at Neptune's south pole compared to regions nearer the equator.Conversely, analysis of Neptune's microwave emission with ALMA suggested strongly enhanced H2S abundances at midlatitudes [2], with much less at the south pole. Determining the true variation of H2S with latitude is crucial for understanding the tropospheric circulation of Neptune.We present our analysis of observations of Neptune from VLT-SINFONI in 2018. Using a limb-darkening approximation, we are able to fit the reflected solar radiance from multiple zenith angles, which allows us to discriminate between gas and aerosol opacity. Despite the lower spectral resolution of this instrument compared to Gemini-NIFS, we are able to detect the H2S spectral signature. With our radiative transfer retrieval code, archNEMESIS [3], we use nested sampling to fit a parameterised cloud model (similar to that of [4]) to these observations over a range of latitudes. We prescribe a latitudinally varying deep methane abundance derived from recent VLT-MUSE observations [5], which enables us to constrain the depth of the cloud top.Our retrieved results are in agreement with the results derived from ALMA [2] - we find a significant enhancement of deep H2S at Neptune's southern midlatitudes, decreasing towards the equator and the pole. Our results show a much deeper cloud top towards the pole, resulting in the increased cloud top column abundance of H2S observed here in the previous near-infrared analysis [1].Figure 1: A comparison of fits to a spectrum extracted from the 50°S to 60°S latitude band, with a model including H2S (blue) and a model without H2S (red). Note the significant discrepancy around 1.58 microns. The models are fitted to spectra at two zenith angles simultaneously.[1] Irwin, P. G., Toledo, D., Garland, R., Teanby, N. A., Fletcher, L. N., Orton, G. S., & Bézard, B. (2019). Probable detection of hydrogen sulphide (H2S) in Neptune’s atmosphere. Icarus, 321, 550-563.[2] Tollefson, J., de Pater, I., Luszcz-Cook, S., & DeBoer, D. (2019). Neptune's latitudinal variations as viewed with ALMA. The Astronomical Journal, 157(6), 251.[3] Alday, J., Penn, J., Irwin, P. G., Mason, J. P., & Yang, J. (2025). archNEMESIS: an open-source Python package for analysis of planetary atmospheric spectra. arXiv preprint arXiv:2501.16452.[4] Irwin, P. G., Teanby, N. A., Fletcher, L. N., Toledo, D., Orton, G. S., Wong, M. H., ... & Dobinson, J. (2022). Hazy blue worlds: a holistic aerosol model for Uranus and Neptune, including dark spots. Journal of Geophysical Research: Planets, 127(6), e2022JE007189.[5] Irwin, P. G., Dobinson, J., James, A., Wong, M. H., Fletcher, L. N., Roman, M. T., ... & de Pater, I. (2023). Latitudinal variations in methane abundance, aerosol opacity and aerosol scattering efficiency in Neptune's atmosphere determined from VLT/MUSE. Journal of Geophysical Research: Planets, 128(11), e2023JE007980.

Optically Observed Ammonia in the Northern Equatorial Zone

(2025)

Authors:

Steven M Hill, Patrick Irwin, John Rogers, Leigh Fletcher

Abstract:

IntroductionJupiter’s northern Equatorial Zone (EZn) and southern North Equatorial Belt (NEBs) are dominated by three features: five-micron hotspots (seen as North Equatorial Dark Features, NEDFs, in the optical), white cloud plumes, and complex local circulation. These features are influenced by the NEBs jet, which is modulated by a meridionally trapped Rossby wave, in conjunction with the high concentration of ammonia in the EZ and the ammonia depletion in the NEB. Numerous measurements have been made of the temperature, aerosol, and ammonia distributions in this region (c.f. Fletcher et al., 2020). And a number of models have been partially successful at explaining the interrelationships between the observed features (c.f Showman & Dowling, 2000). Here we explore the ammonia and cloud height distribution during 2024-25, when NEDFs and five-micron hotspots were prominent, using the optical band-average technique (Hill et al., 2024, Irwin et al., 2025). We show that while many sensing methods highlight the ammonia and aerosol depletion in five-micron hotspots, this band average method highlights enhancements in ammonia to the south of the hotspots.ObservationsMultiple observations on 2025-01-06 were made allowing coverage of a wide range of longitudes and coverage of a given longitude at several zenith angles. Figure 1 shows maps constructed using the method of Hill et al. (2024). An empirical limb correction is applied in addition to a weighted averaging scheme for overlapping observations. The data clearly show that enhanced ammonia regions lie to the south of NEDFs (labeled 1-4 in order of ascending longitude). For the ammonia enhancements we observe a planetary wave number of nine, within the range of hotspot and NEDF wavenumbers typically observed.DiscussionThe NEBs jet speed peaks at about  7° N, which in fact marks the boundary between the NEDFs and the ammonia enhancements. Anticyclonic gyres are a known feature seen in the same location as we show ammonia enhancements (c.f. Choi et al., 2013). We hypothesize that these gyres are regions of uplift and outflow, bringing up ammonia rich air from deeper levels of the atmosphere. The NEDFs are thought to be areas of subsidence, with cyclonic flow, where dryer air descends from above and results in a clearing of aerosols. Figure 1D shows this schematically with upwelling occurring at the gyres, horizontal winds carrying condensates from the upwelling source to the east and northeast as the visible cloud plumes, and descending clear air in the NEDFs.To further support this hypothesis, we analyze the ammonia mole fraction and cloud pressure at the NEDFs, gyres, and in the plumes through a regions-of-interest (ROI) approach. Figure 2 shows a longitudinal subset of the data in Figure 1, focusing on ammonia regions 3 and 4. Rectangles outline the ROIs which are analyzed for three observation times in Figure 2A. Figure 2B shows a time series of average values at each observation time for cloud pressure and ammonia mole fraction along with statistical errors. Finally, 2C shows scatter plots of the average cloud pressure versus the ammonia abundance. Note the very clear clustering of points where the NEB sample provides a consistent reference with relatively high pressure and very low ammonia abundance. Following the upwelling ammonia, eastward advection of plume aerosols, and NEDF subsidence from Figure 1, we can trace an ammonia cycle between its gaseous source and sink, with an intermediary aerosol state.Future WorkHundreds of observations of NEDFs and ammonia enhancements in the EZn have been made in 2024-25 using the Hill et al. (2024) technique. This data set will be analyzed and assessed for the statistical consistency of the results presented here. In addition, this data set will be compared to complementary multispectral observations to help discriminate why the optical method seems to so clearly detect ammonia enhancements at the 1-2 bar pressure level and why these enhancements appear broad enough to overlap NEDFs.Figure 1. Ammonia mole fraction, cloud pressure, and visual context maps created from observations on 2025-01-06 using an 11 inch Schmidt-Cassegrain telescope. A) Ammonia mole fraction (ppm) with enhanced areas labeled 1-4 in order of ascending longitude. The black circle at left shows the approximate spatial resolution of the data. B) Cloud pressure (mbar). C) Visual context image with selected contour overlays to show enhanced ammonia mole fraction and lowest pressure (highest) clouds. D) Same as C), but with arrows indicating presumed upwelling (black ⊙), downwelling (white ⦻), and horizontal flow (red arrows). Figure 2. Two ammonia enhancements (4 & 3 from Figure 1), associated plumes, and NEDFs are analyzed for cloud pressure and ammonia abundance. Three observations are assessed with the targets near nadir viewing. A) Ammonia mole fraction, cloud pressure, and visual context image with overlaid rectangles indicating regions-of-interest (ROIs). B) Time series of cloud pressure (left) and ammonia mole fraction (right) over the three observations. C) Scatter plot of all ammonia and cloud measurements in each ROI (left) and of the averages over the three observations. Note that the NEB data are provided as a stable reference.ReferencesChoi, D. S. et al. 2013. Icarus, 223, 832. Hill, S. M. et al. 2024. Earth and Space Science, 11(8), e2024EA003562.Fletcher, L. N. et al. 2020. Journal of Geophysical Research (Planets), 125, e06399. Irwin, P. G. J. et al. 2025. Journal of Geophysical Research: Planets, 130(1), e2024JE008622. Showman, A. P., & Dowling, T. E. 2000. Science, 289, 1737-1740. 

Photochemistry versus Escape in the Trappist-1 planets.

(2025)

Authors:

Sarah Blumenthal, Richard Chatterjee, Harrison Nicholls, Louis Amard, Shang-Min Tsai, Tad Komacek, Raymond Pierrehumbert

Abstract:

Survive or not survive, that is the question of the 500-hour JWST Rocky Worlds DDT Program. Whether a terrestrial planets’ atmosphere can suffer under the intense XUV of its host, or if it completely escapes, these are the questions we explore. Zahnle & Catling (2017) defined the Cosmic Shoreline, but recent observations from JWST reveal airless worlds around M-stars, calling for a refinement of this “receding” shoreline (Pass et al. 2025). M-stars spend a longer time in pre-main sequence, subjecting their orbiting worlds to some higher intensity XUV activity. This complicates our present understanding of this shoreline. Investigating chemical effects of planet-star interactions could be the key to a more complete picture of this shoreline.  We investigate the interplay between photochemistry, mixing, and escape of carbon dioxide atmospheres under intense and mild XUV fluxes as follow on work to both Johnstone et al. (2018) and Nakayama et al. (2022). We expand on this work by adopting thermal structure models from Nakayama et al. (2022) and apply them to identify key chemical pathways for escape. We create a reduced C-O chemical network including neutral and ionic species to identify these pathways. As photochemistry simulations take into account many reactions, these 1D calculations are too computationally expensive to be done in 3D. Although rudimentary at best, the mixing parameter– eddy diffusion term, K_zz, comprises the dynamical element of 1D photochemical simulations. Here, we consider the mixing of photochemical products in competition with escape to explore the chemical pathways of retention and loss. We compare the photochemical model results for active and inactive cases for the Trappist-1 system planets. Then, using the resulting composition-dependent heating and cooling rates for Trappist-1 planets, we assess their propensity for efficient atomic line cooling versus escape. We follow the work of Chatterjee & Pierrehumbert (2024) in this assessment.  Finally, using our pathway analysis, we find an analytical formula for calculating an energy-limited escape boundary for these planets based on composition.  It is important here to note the limitations of 1D work. First, there exists an exchange of rigor between modelling chemistry and dynamics. Insights from this work are ripe for implementation into 3D GCMs, especially in response to incorporating UV-driven processes for thermospheric modelling mentioned in Ding and Wordsworth (2019). Second, interaction with the interior is important in the early phase of planetary formation, i.e., the magma ocean phase. Due to exchange between atmosphere and magma early in the planet’s formation, incorporation with an interior-atmosphere model would better constrain higher pressure chemical abundances. Although this work focuses on the upper atmosphere, extrapolation to the surface environment is a key goal for understanding a planet.  Considering planet-star interaction is imperative for the selection of targets for observation. However, it is also important when considering anomalous detections of atmospheres around planets predicted to not have an atmosphere. This could be a first step in determining an atmosphere as non-primary and/or distinguishing between an airless planet and one with high altitude haze. 

Quantifying Thin Dust Layer Effects on Thermal-IR Spectra of Bennu-Like Regolith: FTIR Experiments with CI Asteroid Simulant 

(2025)

Authors:

Emma Belhadfa, Neil Bowles, Katherine Shirley

Abstract:

Introduction: The surfaces of airless bodies, such as asteroid (101955) Bennu, are typically composed of a regolith mixture containing both coarse and fine particulates. Observations from NASA’s Origins, Spectral Interpretation, Resource Identification, Security, Regolith Explorer (OSIRIS-REx) mission demonstrated a discontinuity between the remote sensing derived thermophysical properties and thermal spectroscopy results, indicating that a fine layer of dust may be coating the large boulders and coarse regolith surface [1]. To better understand the impact of such a coating on the thermal infrared spectra measured at Bennu, this work developed experimental methods for simulating dust coverings using Space Resource Technology’s CI simulant, based on the bulk composition of the Orgueil meteorite [2].    Figure 1: FTIR Reflectance Spectra of Control Samples of CI simulant. Figure 2: Figure 2: Microscope Camera Images of Sample Surfaces (7%, 10%, 15%, 20%, 25%, 50% Fines wt) of CI simulant Methods: The CI simulant was first sieved into seven size fractions: 1000 mm. An unsieved sample was used as a control. The spectra of the eight samples were measured using a Bruker Vertex 70v Fourier Transform Infrared Reflectance (FTIR) spectrometer, normalized using a gold standard prior to and after measurements, in the range of 1000-650 cm-1 (Figure 1). The dust coating was simulated by placing increasing mass fractions of fine particulates (10% change in the spectral slope).  Implications for OSIRIS-Rex Findings: From the data returned by the OSIRIS-Rex Thermal Emission Spectrometer (OTES) [3],  thermal inertia modelling imply that the surface is porous;, however, the spectral findings indicate that the surface is composed of non-porous? rocks with thin dust coatings [4]. Our experiments find that as little as ~7 wt % of 7% wt) was sufficient to overwhelm and dominate mid-infrared emissivity spectra. The results indicate that the discontinuity in OTES data could be linked back to dust coating on the larger rocks and boulders.   References: [1] Tinker C. et al. (2023) RAS Techniques and Instruments (Vol. 2, Issue 1). [2] Landsman Z. et al. (2020) EPSC 2020. [3] Christensen P. R. et al. (2018) Space Science Reviews (Vol. 214, Issue 5). [4] Rozitis B. et al. (2022) JGR: Planets (Vol. 127, Issue 6). [5] Rivera-Hernandez F. et al. (2015) Icarus (Vol. 262). 

Saturn’s Local and Seasonal Aerosol Variations Inferred from Cassini Combined UV, Visual, and Near-IR Observations  

(2025)

Authors:

James Sinclair, Emma Dahl, Kevin Baines, Tom Momary, Lawrence Sromovsky, Pat Fry, Patrick Irwin

Abstract:

Clouds are the manifestations of atmospheric dynamics, chemistry, thermal evolution, and orbital characteristics; thus, understanding their physical and spectral properties and their spatial and temporal variability is critical to understanding the planet as a whole.  Observations of Saturn by the Hubble Space Telescope since 1994 and by Cassini from 2004 - 2017 have spanned almost one Saturn year (29.5 Earth years).  Despite the wealth of data, a self-consistent picture of the seasonal variations in Saturn’s haze and cloud structure remain elusive. In this work, we present a radiative transfer analysis of Cassini-VIMS (Visible and Infrared Mapping Spectrometer) spectra in order to derive the vertical structure and color properties of Saturn’s clouds and their latitudinal and seasonal variability.  VIMS records spectra over visible (0.3 to 1.05 micron) and infrared (0.85 to 5.1 micron) channels at spectral resolutions of 7 and 16 nm, respectively.  After a review of the VIMS dataset, we have identified dayside spectra that capture unique cloud features in a given latitude circle at multiple emission angles, allowing for improved vertical discrimination of cloud models. Data are additionally available over multiple epochs, allowing us to analyze any seasonal evolution.  Using the NEMESIS radiative transfer code (Irwin et al., 2008, JQSRT 109, 1136-1150), we invert the VIMS spectra to derive the vertical profiles of phosphine (PH3), ammonia (NH3) and the vertical structure of 4 haze/cloud layers (using the cloud model and cloud/gas parameters shown in Figure 1).  In preliminary findings, in adopting the chromophore optical constants derived by Sromovsky et al., 2021 (Icarus 362, 114409) for a north polar cloud observed in 2016, we find we can adequately fit the spectra for a subset of clouds observed in September 2014.  At other locations/times, the chromophore optical constants derived by Sromovsly et al., 2021, need to be varied in order to fit the spectra within uncertainty, which indicates seasonal evolution of Saturn’s chromophore.  In this work, we present derived cloud properties and the optical constants of the derived chromophore as a function of latitude and season in order to shed light on the complex interplay between cloud structure, color, chemistry, and orbital characteristics.