Europa Thermal Emission Imaging System (E-THEMIS) cruise observations of Mars

Copernicus Publications (2025)

Authors:

Philip Christensen, John Spencer, Sylvain Piqueux, Greg Mehall, Saadat Anwar, Oleg Abramov, Paul Hayne, Carly Howett, Michael Mellon, Francis Nimmo, Julie Rathbun, Bonnie Buratti, Robert Pappalardo

Abstract:

The Europa Thermal Emission Imaging System (E-THEMIS) on the Europa Clipper spacecraft will investigate the temperature and physical properties of Europa using thermal infrared images in three wavelength bands at 7-14 µm, 14-28 µm and 28-70 μm [Christensen et al., 2024]. The specific objectives of the investigation are to 1) understand the formation of surface features, including sites of recent or current activity, in order to understand regional and global processes and evolution and 2) to identify safe sites for future landed missions. The E-THEMIS radiometric calibration includes removing the thermal emission from the instrument housing, optical elements, and filters using observations of space and an internal calibration flag [Christensen et al., 2024]. On February 28, 2025, the Clipper spacecraft performed a close flyby of Mars for a trajectory gravity assist. Twenty four hours prior to closest approach the spacecraft pointed the E-THEMIS instrument at Mars and performed a sequence that scanned E-THEMIS across the planet at a slew rate of 100 micro-radians per second. This rate is the same as what be used to image Europa during each flyby [Pappalardo et al., 2024]. This activity accomplished two primary objectives: 1) collect images of a well-characterized target (Mars) to validate the E-THEMIS calibration methodology and software prior to the first observations of Europa; and 2) rehearse the data collection procedure that will be used to obtain global observations of Europa.Mars makes an excellent thermal calibration target because it has been extensively studied and characterized by numerous thermal infrared instruments. The E-THEMIS observations were simulated using modeled surface temperatures generated using global maps of thermal inertia albedo made from the MGS TES data [Christensen et al., 2001], together with the krc thermal model [Kieffer, 2013]. The wavelength-dependent atmospheric absorption and emission was modeled using data from the UAE Emirates Mars Mission EMIRS thermal infrared spectrometer [Amiri et al., 2022; Edwards et al., 2021]. EMIRS collects global scans of hyperspectral data from 6-100 µm at 5 and 10 cm-1 spectral sampling at ~200 km spatial resolution [Edwards et al., 2021]. These spectra were resampled from wavenumber to wavelength and weighted by the three E-THEMIS spectral bandpasses to produce 3-band simulated E-THEMIS global images. EMIRS data were not collected simultaneously with the E-THEMIS imaging, but global observations were acquired at the same season and within 5° of latitude, 10° of longitude, and 0.3 H local time of the E-THEMIS data. Fig. 1 shows an example of the nearest EMIRS observation to the E-THEMIS observing conditions of sub-spacecraft latitude=20.3° N, longitude=163.0° E, local time=11.48 H, and Ls=50.5°. A transfer function from the krc-generated surface temperatures and the bandpass-weighted EMIRS data was created by averaging the ratio of forty-five EMIRS observations to the krc-generated surface temperatures. Simulated E-THEMIS observations were produced using the average of these ratios and the krc surface temperatures. The results are given in Fig. 2. The E-THEMIS data could not be transmitted to Earth until Clipper was more than 2 AU from Sun due to spacecraft thermal constraints. As a result the data were received on Earth on May 7, 2025, and the results and an assessment of the E-THEMIS calibration will be discussed.Fig. 1. Measured Mars temperatures. Comparison of temperature globes for surface temperature (krc model) and E-THEMIS-bandpass-weighted EMIRS data for Bands 1, 2, and 3. The EMIRS observations were acquired on Feb. 17, 2025, at a sub-spacecraft viewing geometry of 16.0° N latitude, 174.1° E longitude, 11.60 H local time, and 45.5° Ls. Fig. 2. Simulated E-THEMIS temperature images. The data for each E-THEMIS band were created using the krc model surface temperatures transferred to E-THEMIS wavelength bands using a transfer function derived from EMIRS observations. ReferencesAmiri, H., D. Brain, O. Sharaf, P. Withnell, M. McGrath, M. Alloghani, M. Al Awadhi, S. Al Dhafri, O. Al Hamadi, and H. Al Matroushi (2022), The emirates Mars mission, Space Science Reviews, 218(1), 4.Christensen, P. R., et al. (2001), The Mars Global Surveyor Thermal Emission Spectrometer experiment: Investigation description and surface science results, J. Geophys. Res., 106, 23,823-823,871.Christensen, P. R., J. R. Spencer, G. L. Mehall, M. Patel, S. Anwar, M. Brick, H. Bowles, Z. Farkas, T. Fisher, and D. Gjellum (2024), The Europa Thermal Emission Imaging System (E-THEMIS) Investigation for the Europa Clipper Mission, Space Science Reviews, 220(4), 1-65.Edwards, C. S., P. R. Christensen, G. L. Mehall, S. Anwar, E. A. Tunaiji, K. Badri, H. Bowles, S. Chase, Z. Farkas, and T. Fisher (2021), The Emirates Mars Mission (EMM) Emirates Mars InfraRed Spectrometer (EMIRS) Instrument, Space science reviews, 217, 1-50.Kieffer, H. H. (2013), Thermal model for analysis of Mars infrared mapping, J. Geophys. Res, 116, 451-470.Pappalardo, R. T., B. J. Buratti, H. Korth, D. A. Senske, D. L. Blaney, D. D. Blankenship, J. L. Burch, P. R. Christensen, S. Kempf, and M. G. Kivelson (2024), Science Overview of the Europa Clipper Mission, Space Science Reviews, 220(4), 1-58.

Independent constraint of Enceladus’ ice shell thickness using thermal observations

Copernicus Publications (2025)

Authors:

Georgina Miles, Carly Howett, Francis Nimmo, Douglas Hemingway

Abstract:

Enceladus maintains its global, unconsolidated ocean around its rocky, porous core by tidal dissipation with Saturn and torque from its resonance with Dione [1].  The active South Polar Terrain (SPT) region is associated with intense concentrations of endogenic heat, but it is the significantly lower-power conductive heat flow that dominates global heat loss as it occurs over the entire surface.  If Enceladus’ global ocean is to be sustained over a significant fraction of its existence, heating rates would have to be balanced endogenic heat loss. Estimates of heating rates from models vary from 1.5-150 GW [2].  The large range results from uncertainty in both the structure of the bodies’ interiors and their evolution.  Ice shell thickness/shape models, which interpret gravity, libration and topographic data, produce global conductive heat loss estimates of around 18-35 GW [3,4,5].Endogenic heat loss from the SPT has been estimated using thermal observations from Cassini Composite Infrared Spectrometer (CIRS) to be between 5-19 GW [6,7,8], resulting in a conventional, combined heat loss estimate of around 50 GW [9].Detecting endogenic heat loss using thermal observations presents a significant challenge, principally relating to limited data coverage and uncertainty about the surface thermal properties but is possible under some circumstances [10].We use thermal observations CIRS to identify endogenic heat at the north pole of Enceladus in the form of conductive heat flow.   From this estimate we can infer global average heat loss.  We are then able to invoke the same mechanisms used to estimate the global average heat loss from ice shell thickness models [5, 9] to characterize the first north polar and global average ice shell thicknesses independently derived from thermal observations.Acknowledgments: This work was made possible through NASA’s support of Cassini Data Analysis Program Grant Number 80NSSC20K0477. References[1] Nimmo, F., Barr, A.C., Behounková, M. and McKinnon, W.B., 2018. The thermal and orbital evolution of Enceladus: observational constraints and models. Enceladus and the icy moons of Saturn, 475, pp.79-94.[2] Lainey, V., Casajus, L.G., Fuller, J., Zannoni, M., Tortora, P., Cooper, N., Murray, C., Modenini, D., Park, R.S., Robert, V. and Zhang, Q., 2020. Resonance locking in giant planets indicated by the rapid orbital expansion of Titan. Nature Astronomy, 4(11), pp.1053-1058.[3] Thomas, P.C., Tajeddine, R., Tiscareno, M.S., Burns, J.A., Joseph, J., Loredo, T.J., Helfenstein, P. and Porco, C., 2016. Enceladus’s measured physical libration requires a global subsurface ocean. Icarus, 264, pp.37-47.[4] Čadek, O., Tobie, G., Van Hoolst, T., Massé, M., Choblet, G., Lefèvre, A., Mitri, G., Baland, R.M., Běhounková, M., Bourgeois, O. and Trinh, A., 2016. Enceladus's internal ocean and ice shell constrained from Cassini gravity, shape, and libration data. Geophysical Research Letters, 43(11), pp.5653-5660.[5] Hemingway, D.J. and Mittal, T., 2019. Enceladus's ice shell structure as a window on internal heat production. Icarus, 332, pp.111-131.[6] Spencer, J.R., Pearl, J.C., Segura, M., Flasar, F.M., Mamoutkine, A., Romani, P., Buratti, B.J., Hendrix, A.R., Spilker, L.J. and Lopes, R.M.C., 2006. Cassini encounters Enceladus: Background and the discovery of a south polar hot spot. science, 311(5766), pp.1401-1405.[7] Howett, C.J.A., Spencer, J.R., Pearl, J. and Segura, M., 2011. High heat flow from Enceladus' south polar region measured using 10–600 cm− 1 Cassini/CIRS data. Journal of Geophysical Research: Planets, 116(E3).[8] Spencer, J.R., Nimmo, F., Ingersoll, A.P., Hurford, T.A., Kite, E.S., Rhoden, A.R., Schmidt, J. and Howett, C.J., 2018. Plume origins and plumbing: from ocean to surface. Enceladus and the icy moons of Saturn, 163.[9] Nimmo, F., Neveu, M. and Howett, C., 2023. Origin and evolution of Enceladus’s tidal dissipation. Space Science Reviews, 219(7), p.57[10] Miles, G., Howett., C., Spencer J., Vol. 16, EPSC2022-1190, 2022, https://doi.org/10.5194/epsc2022-119

Investigating Phobos' Origin using X-ray Diffraction and Reflectance Spectroscopy of Meteorites.

(2025)

Authors:

Emelia Branagan-Harris, Neil E Bowles, Ashley J King, Katherine A Shirley, Helena C Bates, Sara S Russell

Abstract:

Introduction: The origins of Mars' moons, Phobos and Deimos, remain uncertain, with two main hypotheses under consideration: formation from debris following a high-energy impact between Mars and an asteroid [1], or capture of primitive asteroids [2]. To address this, JAXA's Martian Moons eXploration (MMX) mission aims to return samples from Phobos by 2031 [3]. The characterisation of these samples will determine the origin of Phobos.To ground-truth remote observations of Phobos, we have used X-ray diffraction (XRD), and Fourier transform infrared (FTIR) reflectance spectroscopy to characterise the bulk mineralogy and IR spectral properties of ureilites, carbonaceous and ordinary chondrites, the composition of which could be indicative of a captured asteroid [4], and Martian meteorites that could represent a collisional formation. By acquiring XRD and IR data from the same material, mineral abundances can be directly correlated with features in reflectance spectra [5]. When MMX reaches Phobos, meteorite data collected in the laboratory will play a crucial role towards interpreting the mineralogy and composition of materials on its surface.Methods: We have characterised the mineralogy and spectral properties of six CM (Mighei-like) carbonaceous chondrites, Tarda (C2-ung), the CO (Ornans-like) chondrite Kainsaz, a range of shock darkened ordinary chondrites (mostly falls) including L4-6, and H5-6, four CR2 chondrites, four ureilites, Martian meteorites Nakhla and Tissint, and a Tagish Lake (C2-ung) based simulant created by the University of Tokyo, known as UTPS-TB [6]. For the meteorites, chips of approximately 200 mg were ground to produce powders with grain sizes of less than 40 microns. The UTPS-TB sample came in a powdered form which was ground to the same grain size as the meteorites.Diffuse reflectance spectra (1.7 - 50 μm) were collected using a Bruker VERTEX 70V FTIR spectrometer at the University of Oxford Planetary Spectroscopy Facility. Spectra were calibrated at the start of each measurement day and between measurements of samples using a gold standard. The powdered sample was measured under a vacuum to reduce terrestrial atmospheric contributions.XRD patterns of the same powders were collected using an INEL X-ray diffractometer with a position-sensitive detector at the Natural History Museum, London. Around 50 mg of powdered sample was measured for 16 hours to achieve good signal-to-noise. Measurements of well-characterised standard minerals were collected for 30 minutes and compared with meteorite patterns to identify minerals and quantify their abundance in the sample [e.g. 7].Results & Discussion: The mineralogical and spectral characteristics of meteorites in this investigation are compared the reflectance spectra of Phobos’ surface. The CR chondrites are primitive, containing both anhydrous silicates (e.g. olivine and pyroxene) and aqueous alteration phases such as phyllosilicates, carbonates, magnetite, and sulfides. Their albedo is ~3-5% reflectance with a weak red slope in the visible to near-infrared (VNIR). The CRs have a 3 μm hydration band, due to partial aqueous alteration. Their low VNIR reflectance, red-sloped continuum, and weak 3 μm spectral absorption feature is like that of Phobos, supporting the captured asteroid origin theory. The CM chondrites share similar spectral features but have a lower albedo and a stronger μm hydration band, corresponding to a higher phyllosilicate composition.   The ureilites are achondritic ultramafic meteorites containing olivine, pyroxene and carbon phases. These samples have a low albedo (~6-15% in VNIR) due to their opaque carbonaceous composition. However, their VNIR spectra are blue-sloped, inconsistent with Phobos’ red-sloped spectra. Ureilites are also anhydrous and therefore lack the 3 μm hydration band seen in Phobos spectra. Their low reflectance and feature-poor spectra could resemble Phobos, however there is a significant difference in spectral slope and hydration features. Therefore, Phobos were composed of ureilitic material, its surface would need to be significantly modified by space weathering.Martian meteorites Nakhla (a nakhlite) and Tissint (a shergottite) have mineralogical and spectral features consistent with their basaltic origin. XRD measurements of these meteorites are dominated by pyroxene (augite, pigeonite), and olivine, consistent with their origin in the Martian crust. Their reflectance spectra have relatively high albedo, mafic absorption bands at ~1 and 2 μm, and a lack of hydration features. These features are inconsistent with the spectra of Phobos, which lack 1 or 2 μm bands and show significantly lower reflectance.CR and CM chondrites are the closest spectral match to Phobos from the samples studied. Their low albedo, red-sloped, hydrated spectra are consistent with surface measurements of Phobos. Ureilites share low reflectance but differ significantly in slope and hydration, while Martian meteorites differ in more spectral characteristics. These results support the interpretation that Phobos is composed of primitive, carbon-rich material, likely of outer solar system origin, and favour a capture scenario over a collisional formation from Martian ejecta. The similarities between the carbonaceous chondrites and Phobos indicates that the Martian moons may be captured asteroids and further demonstrates the importance of the MMX mission sample return for solving the mystery of their origin definitively.References: [1] R. Citron et al. (2015) Icarus 252:334-338. [2] M. Pajola et al. (2013) The Astrophysical Journal 777:127. [3] K. Kuramoto et al. (2022) Earth, Planets and Space 74:12. [4] K. D. Pang et al. (1978) Science 199(4324):64-66. [5] H. C. Bates et al. (2023) Meteoritics & Planetary Science 1-23. [6] H. Miyamoto et al. (2021) Earth, Planets and Space 73:1-17 [7] G. Cressey et al. (1996) Powder Diffraction 11:35-39.

Investigating the Vertical Variability of Titan’s 14N/15N in HCN

(2025)

Authors:

Alexander Thelen, Katherine de Kleer, Nicholas Teanby, Amy Hofmann, Martin Cordiner, Conor Nixon, Jonathon Nosowitz, Patrick Irwin

Abstract:

Titan’s substantial atmosphere is primarily composed of molecular nitrogen (N2) and methane (CH4), which are dissociated by solar UV photons and subsequently generate a vast chemical network of trace gases. The composition of Titan’s atmosphere is markedly different than that of Saturn, including both the complex molecular inventory and the hitherto measured isotopic ratios – including that of nitrogen (14N/15N). Atmospheric and interior evolution models (e.g., Mandt et al., 2014) indicate that the atmospheres of Saturn and Titan did not form in the same manner or from the same constituents, and that Titan’s atmospheric N2 may have originated from its interior as NH3. The evolution of 14N/15N in Titan’s atmosphere over time does not result in a value comparable to that measured on Saturn and instead is closer to cometary values; this indicates that the origin of Titan’s atmosphere appears to be from protosolar planetesimals enriched in ammonia and not from the sub-Saturnian nebula. However, selective isotopic fractionation of molecular species in Titan’s atmosphere complicates this picture, as the isotopic ratios may vary as a function of altitude (Figure 1). To further constrain the evolution of Titan’s atmosphere – and indeed, its origin – isotopic ratios must be measured throughout its atmosphere, instead of being interpreted from bulk values likely only representative of the stratosphere.While the measurement of Titan’s 14N/15N in N2 (167.7; Niemann et al. 2010) places it firmly below the lower limit derived for Saturn (~350; Fletcher et al., 2014), Titan’s atmospheric nitriles (e.g., HCN, HC3N, CH3CN) are further enriched in 15N, resulting in ratios closer to 70 (Molter et al., 2016; Cordiner et al., 2018; Nosowitz et al., 2025). The variation in nitrogen isotopic ratios between the nitriles and N2 is thought to be the result of higher photolytic efficiency of 15N14N compared to N2 in the upper atmosphere (~900 km), resulting in increased 15N incorporated into nitrogen-bearing species (Liang et al., 2007; Dobrijevic & Loison, 2018; Vuitton et al., 2019). As these species are advected to lower altitudes, the nitrogen isotope ratio may vary vertically (Figure 1, red and black profiles), but previous measurements have only presented bulk atmospheric isotope ratios primarily representing Titan’s stratosphere (Figure 1, blue lines).Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) have allowed for the derivation of vertical abundance profiles of Titan’s trace atmospheric species and measurements of N, D, and O-bearing isotopologues (Molter et al., 2016; Serigano et al., 2016; Cordiner et al., 2018; Thelen et al., 2019; Nosowitz et al., 2025). However, vertical isotopic ratio profiles have yet to be derived. Here, we utilize observations acquired with ALMA in July 2022 containing high sensitivity measurements of the HC15N J=4–3 transition at 344.2 GHz (~ 0.87 mm) to investigate vertical variations in the 14N/15N of Titan’s HCN. We compare the results of the vertical 14N/15N profile to those predicted by photochemical models to determine the impact of the isotopic-selective photodissociation of nitrogen-bearing molecular species in Titan’s atmosphere, and the impact of the Saturnian and space environments that vary between model implementations.Figure 1. 14N/15N profile for HCN predicted by photochemical models from Vuitton et al. (2019; black line) and Dobrijevic & Loison (2018; red line). Blue colored bars in the lower atmosphere represent previous HCN nitrogen isotope ratios from Cassini, Herschel, and ground-based (sub)millimeter observations (see Molter et al., 2016, and references therein). Measurements are offset vertically for clarity, and all refer to HC14N/HC15N measurements for the bulk stratosphere.References:Cordiner et al., 2018, The Astrophysical Journal Letters, 859, L15.Dobrijevic & Loison, 2018, Icarus, 307, 371.Fletcher et al., 2014, Icarus, 238, 170.Liang et al., 2007, The Astrophysical Journal Letters, 644, L115.Mandt et al. 2014, The Astrophysical Journal Letters, 788, L24.Molter et al., 2016, The Astronomical Journal, 152, 42.Niemann et al., 2010, Journal of Geophysical Research, 115, E12006.Nosowitz et al., 2025, The Planetary Science Journal, 6, 107.Serigano et al., 2016, The Astrophysical Journal Letters, 821, L8.Thelen et al., 2019, The Astronomical Journal, 157, 219.Vuitton et al., 2019, Icarus, 324, 120.

Jovian chromophore and upper hazes from CARMENES spectra

(2025)

Authors:

José Ribeiro, Pedro Machado, Santiago Pérez-Hoyos, Asier Anguiano-Arteaga, Patrick Irwin

Abstract:

The nature of the red colouration of Jupiter’s belts and some of its major anticyclones is still debated to this day. Sromovsky et al. (2017) proposed the existence of an “universal chromophore” by fitting Cassini/VIMS-V observations. Baines et al. (2019) concluded that this chromophore should be located in a thin layer above the ammonia clouds, giving rise to the so called “Crème Brûlée” model. Both of these works had as a basis the red compound that formed through the reaction of photolyzed ammonia with acetylene as obtained in the laboratory by Carlson et al. (2016).However,  both Pérez-Hoyos et al. (2020) and Braude et al. (2020) found that a less blue and more vertically extended chromophore layer would fit better their HST/ WFC3 North Temperate Belt disturbance observations for the former and latitudinal swath from MUSE/VLT observations for the later, without fully discarding the possible existence of an “universal chromophore”. Recently, analysis of HST/WFC3 images of Jupiter’s Great Red Spot, its surroundings, and, Oval BA by Anguiano-Arteaga et al. (2021,2023) suggest the presence of two distinct colouring aerosols. The first being similar to the “universal chromophore” and the second one being a new UV-absorbing species below the main chromophore layer at tropospheric levels. This highlights the uncertainties on the vertical distribution of aerosols, their properties and their variability.To address this uncertainty, we used new Jupiter spectra obtained with CARMENES (The Calar Alto High-Resolution search for M dwarfs with Exoearths with Near-infrared and optical Échelle Spectrographs) in 2019. This instrument consists of two separated spectrographs with spectral resolutions R = 80,000-100,000, covering the wavelength ranges of 0.52 to 0.96 μm and of 0.96 to 1.71 μm. The original purpose of these observations was to measure winds through the Doppler velocimetry method. We used a downgraded resolution version (R = 173-570) so the observations match the available spectral data for methane, as this resolution is enough for constraining aerosol properties. Due to the original nature of the observations, no calibration star was recorded. In order to achieve flux calibration, we used  2017 observations of Saturn with CARMENES. We employed Saturn’s B ring to obtain the response function of the instrument, since no other sources of calibration are available at the desired resolution or epoch.We used the reflectivity (I/F) spectrum obtained with Cassini/VIMS (Cuzzi et al., 2009) at phase angles less than 3º. We applied the response function to the centre of disc spectrum of Saturn and compared the obtained reflectivity spectrum with results from Clark and McCord (1979) and Mendikoa, et al. (2017). Lastly, we applied the flux calibration to the Jupiter observations and compared them results from Mendikoa, et al., (2017) and Irwin et al. (2018) (Figure 1). All calibrations agree within 10% with MUSE calibration.We were able to perform a Minnaert Limb-darkening approximation and produce 2 synthetic spectra (zenith angle = 0º/61.45º) for five distinct sample areas (EZ (Figure 2), SEB, NEB, transition region from EZ to SEB, and from NEB to NTrZ). We performed retrievals using the same a priori atmospheric parameterization as presented in Braude et al. (2020), Pérez-Hoyos et al. (2020) and Anguiano-Arteaga et al. (2021), comparing the retrieved results of each in order to constrain the uncertainties in the Jovian aerosol scheme. To achieve this, we used the NEMESIS (Nonlinear Optimal Estimator for MultivariatE Spectral analySIS) radiative transfer suite (Irwin et al., 2008). We present here the results of this analysis.Figure 1: Comparison of centre of disk Jupiter spectrum after flux calibration with EZ spectrum from Irwin et al. (2018) and 0º latitude spectrum from Mendikoa et al. (2017).Figure 2: Observation spectra compared to the obtained synthetic spectra after retrieving the atmospheric parameters for the EZ using Braude et al. (2020) model. Top row corresponds to nadir (incidence and emission angle = 0º) and bottom row to limb (incidence and emission angle = 61.45º). Figure 3: Comparison between the a priori aerosol vertical profiles and the retrieved profiles for every region for the model from Braude et al. (2020).  References:Carlson, R. W., et al. (2016). Chromophores from photolyzed ammonia reacting with acetylene: Application to Jupiter's Great Red Spot. Icarus, 274, 106–115. Sromovsky, L. A., et al. (2017). A possibly universal red chromophore for modeling color variations on Jupiter. Icarus, 291, 232–244. Baines, K. H., et al. (2019). The visual spectrum of Jupiter's Great Red Spot accurately modelled with aerosols produced by photolyzed ammonia reacting with acetylene. Icarus, 330, 217–229. Pérez-Hoyos, S., et al. (2020). Color and aerosol changes in Jupiter after a North temperate belt disturbance. Icarus, 132, 114021. Braude, A. S., et al. (2020). Colour and tropospheric cloud structure of Jupiter from MUSE/VLT: Retrieving a universal chromophore. Icarus, 338, 113589. Anguiano-Arteaga, A., et al. (2021). Vertical Distribution of Aerosols and Hazes Over Jupiter's Great Red Spot and Its Surroundings in 2016 From HST/WFC3 Imaging. Journal of Geophysical Research: Planets, 126, e2021JE006996. Anguiano-Arteaga, A., et al. (2023). Temporal variations in vertical cloud structure of Jupiter's Great Red Spot, its surroundings and Oval BA from HST/WFC3 imaging. Journal of Geophysical Research: Planets, 128, e2022JE007427. Karkoschka, E. (1994). Spectrophotometry of the Jovian Planets and Titan at 300- to 1000-nm Wavelength: The Methane Spectrum. Icarus, 111, 1, 174–192. Irwin, P., et al. (2008). The NEMESIS planetary atmosphere radiative transfer and retrieval tool. J. Quant. Spectrosc. Radiat. Transf., 109, 1136–1150. Rodgers CD. (2000). Inverse methods for atmospheric sounding: theory and practice. Singapore: World Scientific. Cuzzi, J., et al., 2009. Ring Particle Composition and Size Distribution. Springer Netherlands, Dordrecht. pp. 459–509. Clark, R.N., McCord, T.B., 1979. Jupiter and Saturn: Near-infrared spectral albedos. Icarus 40, 180–188. Mendikoa, I., et al., 2017. Temporal and spatial variations of the absolute reflectivity of Jupiter and Saturn from 0.38 to 1.7 𝜇m with planetcam-upv/ehu. A&A 607, A72. Irwin, P.G., et al., 2018. Analysis of gaseous ammonia (NH3) absorption in the visible spectrum of Jupiter. Icarus 302, 426–436