# The H 3 + ionosphere of Uranus: decades-long cooling and local-time morphology

Henrik Melin, L. Fletcher, T. Stallard, S. Miller, L. Trafton, L. Moore, J.

O'Donoghue, R. Vervack, N. Dello Russo, L. Lamy, et al.

► To cite this version:

Henrik Melin, L. Fletcher, T. Stallard, S. Miller, L. Trafton, et al. The H 3 + ionosphere of Uranus: decades-long cooling and local-time morphology. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2019, 377 (2154), pp.20180408. 10.1098/rsta.2018.0408. hal-03278908

HAL Id: hal-03278908

https://hal.science/hal-03278908v1

Submitted on 28 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research

<figure><img src="image_2.png" /></figure>

Check for
updates

**Cite this article:** Melin H *et al*. 2019 The H<sup>+</sup>
ionosphere of Uranus: decades-long cooling
and local-time morphology. *Phil. Trans. R.
Soc. A* **377**: 20180408.

http://dx.doi.org/10.1098/rsta.2018.0408

Accepted: 7 April 2019

One contribution of 18 to a discussion meeting
issue 'Advances in hydrogen molecular ions:
H<sup>+</sup>, H<sub>5</sub><sup>+</sup> and beyond'.

**Subject Areas:**

Solar System

**Keywords:**

Uranus, aeronomy, spectroscopy

**Author for correspondence:**

Henrik Melin
e-mail: henrik.melin@leicester.ac.uk

# The H<sup>+</sup> ionosphere of Uranus: decades-long cooling and local-time morphology

Henrik Melin<fnref n="1" />, L. N. Fletcher<fnref n="1" />, T. S. Stallard<fnref n="1" />, S. Miller<fnref n="2" />,
L. M. Trafton<fnref n="3" />, L. Moore<fnref n="4" />, J. O’Donoghue<fnref n="5" />,
R. J. Vervack Jr<fnref n="6" />, N. Dello Russo<fnref n="6" />, L. Lamy<fnref n="7" />, C. Tao<fnref n="8" />
and M. N. Chowdhury<fnref n="1" />

<figure><img src="image_3.png" /></figure>

HM, 0000-0001-5971-2633; LNF, 0000-0001-5834-9588;
TSS, 0000-0003-3990-670X; LM, 0000-0003-4481-9862;
MNC, 0000-0003-3577-9395

The upper atmosphere of Uranus has been observed to be slowly cooling between 1993 and 2011. New analysis of near-infrared observations of emission from H<sup>+</sup> obtained between 2012 and 2018 reveals that this cooling trend has continued, showing that the upper atmosphere has cooled for 27 years, longer than the length of a nominal season of 21 years. The new observations have offered greater spatial resolution and higher sensitivity than previous ones, enabling the characterization of the H<sup>+</sup> intensity as a function of local time. These profiles peak between 13 and 15 h local time, later than models suggest. The NASA Infrared Telescope Facility

<footnote n="1">¹Department of Physics & Astronomy, University of Leicester, Leicester, UK</footnote>

<footnote n="2">²Department of Physics & Astronomy, University College London, London, UK</footnote>

<footnote n="3">³Department of Astronomy, University of Texas, Austin, TX, USA</footnote>

<footnote n="4">⁴Center for Space Physics, Boston University, Boston, MA, USA</footnote>

<footnote n="5">⁵Goddard Space Flight Centre, Greenbelt, MD, USA</footnote>

<footnote n="6">⁶Johns Hopkins Applied Physics Laboratory, Laurel, MD, USA</footnote>

<footnote n="7">⁷LESIA, Observatoire de Paris, PSL, CNRS, Sorbonne Université, Meudon, France</footnote>

<footnote n="8">⁸National Institute of Information and Communications Technology, Tokyo, Japan</footnote>
iSHELL instrument also provides the detection of a bright $H_3^+$ signal on 16 October 2016, rotating into view from the dawn sector. This feature is consistent with an auroral signal, but is the only of its kind present in this comprehensive dataset.

This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: $H_3^+$, $H_5^+$ and beyond'.

# 1. Introduction

William Herschel's discovery of Uranus in 1781, announced in this journal [1], marked the first addition to the roster of planets in our Solar System since antiquity. With Herschel's discovery of two moons orbiting the planet 6 years later [2], Uranus was demonstrated to have a very large tilt, with the rotational axis almost aligned with the plane of the ecliptic, producing extreme seasons over its 84 year journey about the Sun. This subjects certain regions of the atmosphere to extreme contrasts, from being fully illuminated by the Sun to being entirely in darkness.

In 1986, the Voyager 2 spacecraft flew past Uranus, providing our first and only detailed close-up view of the planet. It revealed a strange magnetic field, offset approximately 60° from the rotational axis, with the dipole axis offset $0.3 R_U$ ($1 R_U = 25362$ km) from the centre of the planet [3]. With the rotational axis aligned with the ecliptic plane, the magnetic field configuration with respect to the interplanetary magnetic field (IMF) changes dramatically throughout each Uranian day ($17.24 \pm 0.01$ h; [4]), but also throughout the Uranian year ($P = 84$ yr). In addition, the angle of attack between the magnetic dipole and the solar-wind flow varies less at solstice than at equinox. During the flyby, Voyager 2 discovered weak auroral emissions in the ultraviolet, dotted about the magnetic poles [5], with the total emitted energy flux being approximately 50 times weaker than the aurora observed at Saturn (e.g. [5,6]). More recently, observations using the Hubble Space Telescope (HST) [7,8] have re-detected the ultraviolet aurora, imaging the dayside aurora before and after equinox in both hemispheres. Near-equinox aurora are weak, spot-like and intermittent but are observed regularly between 2011 and 2014, with the power, size and occurrence rate all increasing during that time. The HST observations obtained in 2014 were about as bright as those observed by Voyager, and the morphology was described well by model auroral ovals.

Owing to the large offset between the rotational axis and the magnetic dipole axis, auroral emission appears close to the rotational equator, about the magnetic poles, with the HST observations showing emissions at a maximum latitude of -50° [8]. Here, the Uranus Longitude System (ULS) is used [3], which defines the visible pole during the Voyager 2 encounter as the northern, while the International Astronomical Union defines it as the southern. These emissions are likely to be driven by impulsive dayside reconnection, but could also exhibit more stable field-aligned currents that close in the ionosphere. These processes can inject significant amounts of energy into the upper atmosphere by way of Joule heating. Theoretical considerations [9,10] have indicated that the geometry of Uranus at equinox in 2007 was not favourable for building up significant flux in the magnetotail, consistent with weak and intermittent auroral emissions. By contrast, during solstice conditions, reconnection is favourable at the magnetopause and magnetic flux can be accumulated in the tail, conditions favourable for generating aurora. This picture is consistent with a seasonal dependence on the intensity of the observed auroral emissions.

As highlighted throughout this special issue, the molecular ion $H_3^+$ is an important tracer of energy being injected into a system dominated by molecular hydrogen, and by analysing its spectrum one can determine the temperature of the region in which the ions are formed, along with the column integrated line-of-sight density. $H_3^+$ is a very efficient emitter in the infrared and can provide significant radiative cooling in hot environments (e.g. at Jupiter [11]). An increase in the $H_3^+$ temperature can be indicative of localized heating, while the density of $H_3^+$ is determined by the balance between $H_2$ ionization rates and the loss via recombination with electrons.

$H_3^+$ emissions from Uranus were discovered on the 1 April 1992 [12] using the United Kingdom Infrared Telescope (UKIRT). In the period between 1993 and 2009, a variety of
observing programmes, using a range of ground-based telescopes, recorded $H_3^+$ emissions from Uranus. Melin *et al.* [13] returned to these datasets and self-consistently re-analysed the observations, fitting temperatures and line-of-sight densities to the observed $H_3^+$ spectra. Owing to the relatively low signal to noise (S/N) of the data and small angular size of Uranus, these temperatures represent globally averaged temperatures, providing a unique view of the ionosphere as a unitary system. In this time series, they discovered that the upper atmosphere of Uranus had dramatically cooled from $715 \pm 47$ K in 1992 to $534 \pm 39$ K in 2008. This slow yet consistent cooling was initially interpreted to be related to seasonal solar irradiance, modulating the effectiveness of the Joule heating by changing ionospheric conductivity, where the solstice (1986) would appear hotter than the equinox (2007) since the illuminated area on the planet over the course of a day is greater by about a factor 2 at equinox. This hypothesis predicted that the cooling of the upper atmosphere would reverse to an epoch of heating at equinox, after some time-lag governed by the thermal insulation of this region. However, subsequent observations in 2011 [14] revealed that the upper atmosphere had continued to cool further, down to $520 \pm 32$ K, indicating either that the thermal lag is extremely long, or that the hypothesis of purely solar-driven seasons in the upper atmosphere is incomplete.

On a much larger perspective, one of the outstanding questions in giant planet aeronomy is why the upper atmosphere of all these planets is several hundreds of degrees Kelvin hotter than solar input alone can produce. Both heating by global re-distribution of auroral energy injected about the magnetic poles, and heating by breaking of acoustic or gravity waves generated by the turbulent lower atmosphere have been proposed as potential solutions to this 'energy crisis'. There is evidence that both of these could contribute significant heat under specific circumstances (for Jupiter, see e.g. [15,16]), but a cohesive picture that completes the energy budget at all the giant planets is still missing. Understanding how the temperature of Uranus' upper atmosphere evolves over time will add important constraints to our understanding of these processes.

This study extends the baseline of observations with ones obtained in 2012–2018, bringing the total period of $H_3^+$ observations of Uranus to 27 years, longer than an individual solar season at Uranus (21 years), adding data from the NASA Infrared Telescope Facility (IRTF), Keck, the Very Large Telescope (VLT) and Gemini. These facilities also enable the characterization of the $H_3^+$ emission across the disc of Uranus.

## 2. Observations

Because of the small angular size of Uranus in the sky (approx. 3.7") and large distance from the Sun (approx. 19 AU), the $H_3^+$ signal observed from the planet generally has low S/N. This, in combination with the large pixel scales of early instrumentation (e.g. [12] used a 3.1"/pixel instrument) means that the temperatures derived by Melin *et al.* [13,14] are effectively a global average over the period of observation, which is often a significant fraction of complete longitude coverage. More recent observations with smaller pixel scales, narrower slits and higher sensitivity, allow a more detailed view of the ionosphere of Uranus. However, in order to be able to compare like-for-like with historical observations, the more recent observations have initially been averaged over each night of observation to produce a similar global view of the ionosphere for each particular epoch. This works particularly well when the spectrograph slit is scanned across the disc of the planet, providing broader latitude coverage. Throughout this paper, the ULS [3] is employed, defining the north pole as the sunlit pole during the Voyager 2 encounter.

Here we analyse $H_3^+$ observations of Uranus obtained between 2012 and 2018. This includes observations from five different instruments on four different telescopes. Table 1 details the observations analysed in this study, listing the observations day-of-year (DOY) and mid-time (UTC), the length of time between the first and the last exposure, the calibration star used (all of spectral type A0), and the telescope-specific programme ID. Figure 2 shows an example spectrum from each instrument used in this study, and a brief description of each is provided below.
<table><caption>**Table 1.** The separate nights of mid-infrared observations of H<sub>3</sub><sup>+</sup> emission from Uranus analysed in this study. The mid-observation day-of-year (DOY) and UTC time are listed, along with the time between the first and the last exposure, the A0 calibration star used and the telescope specific programme ID.</caption><thead><tr><th>ID</th><th>facility and instrument</th><th>DOY mid-obs</th><th>t (h)</th><th>Cal. Star</th><th>program ID</th></tr></thead><tbody><tr><td>1</td><td>NASA IRTF SpeX</td><td>2012-229 13.15</td><td>4.8</td><td>HR 8911</td><td>2012B070</td></tr><tr><td>2</td><td>NASA IRTF SpeX</td><td>2012-321 07.04</td><td>4.6</td><td>HR 8911</td><td>2012B070</td></tr><tr><td>3</td><td>Gemini GNIRS</td><td>2012-348 06.36</td><td>3.2</td><td>HR 718</td><td>GN-2012B-Q-114</td></tr><tr><td>4</td><td>NASA IRTF SpeX</td><td>2013-215 13.12</td><td>4.6</td><td>HR 8573</td><td>2013B010</td></tr><tr><td>5</td><td>NASA IRTF SpeX</td><td>2013-216 13.24</td><td>4.2</td><td>HR 8573</td><td>2013B010</td></tr><tr><td>6</td><td>VLT CRIRES</td><td>2013-295 04.44</td><td>3.2</td><td>HR 8911</td><td>092.C-0077(A)</td></tr><tr><td>7</td><td>VLT CRIRES</td><td>2013-296 04.54</td><td>4.6</td><td>HR 8911</td><td>092.C-0077(A)</td></tr><tr><td>8</td><td>VLT CRIRES</td><td>2013-297 02.54</td><td>1.9</td><td>HR 125</td><td>092.C-0077(A)</td></tr><tr><td>9</td><td>VLT CRIRES</td><td>2013-300 04.08</td><td>3.8</td><td>HR 125</td><td>092.C-0077(A)</td></tr><tr><td>10</td><td>NASA IRTF SpeX</td><td>2013-319 09.32</td><td>4.2</td><td>HR 658</td><td>2013B010</td></tr><tr><td>11</td><td>Gemini North GNIRS</td><td>2013-340 08.08</td><td>2.5</td><td>HR 378</td><td>GN-2013B-Q-93</td></tr><tr><td>12</td><td>Gemini North GNIRS</td><td>2013-359 06.14</td><td>2.2</td><td>HR 8826</td><td>GN-2013B-Q-93</td></tr><tr><td>13</td><td>NASA IRTF SpeX</td><td>2014-255 12.04</td><td>5.2</td><td>HR 658</td><td>2014B034</td></tr><tr><td>14</td><td>NASA IRTF SpeX</td><td>2014-256 12.22</td><td>4.6</td><td>HR 658</td><td>2014B034</td></tr><tr><td>15</td><td>Keck II NIRSPEC</td><td>2014-285 11.29</td><td>0.9</td><td>HR 718</td><td>2014BN122NS</td></tr><tr><td>16</td><td>Keck II NIRSPEC</td><td>2014-286 11.15</td><td>1.3</td><td>HR 718</td><td>2014BN122NS</td></tr><tr><td>17</td><td>NASA IRTF SpeX</td><td>2014-328 06.28</td><td>3.0</td><td>HR 8518</td><td>2014B034</td></tr><tr><td>18</td><td>NASA IRTF SpeX</td><td>2014-337 07.26</td><td>5.2</td><td>HR 8518</td><td>2014B034</td></tr><tr><td>19</td><td>NASA IRTF SpeX</td><td>2015-325 08.16</td><td>7.1</td><td>HR 658</td><td>2015B031</td></tr><tr><td>20</td><td>NASA IRTF iSHELL</td><td>2016-284 09.32</td><td>4.1</td><td>HR 311</td><td>2016A041</td></tr><tr><td>21</td><td>NASA IRTF iSHELL</td><td>2016-285 10.15</td><td>5.5</td><td>HR 8518</td><td>2016A041</td></tr><tr><td>22</td><td>NASA IRTF iSHELL</td><td>2016-320 07.42</td><td>6.0</td><td>HR 8518</td><td>2016A041</td></tr><tr><td>23</td><td>NASA IRTF iSHELL</td><td>2017-230 14.58</td><td>1.2</td><td>HR 840</td><td>2017B077<sup>a</sup></td></tr><tr><td>24</td><td>NASA IRTF iSHELL</td><td>2017-244 13.30</td><td>1.2</td><td>HR 8865</td><td>2017A037</td></tr><tr><td>25</td><td>NASA IRTF iSHELL</td><td>2017-245 09.52</td><td>1.3</td><td>HR 8865</td><td>2017A037</td></tr><tr><td>26</td><td>NASA IRTF iSHELL</td><td>2017-258 12.56</td><td>5.1</td><td>HR 1718</td><td>2017A037</td></tr><tr><td>27</td><td>NASA IRTF iSHELL</td><td>2017-259 12.58</td><td>5.2</td><td>HR 1718</td><td>2017A037</td></tr><tr><td>28</td><td>NASA IRTF iSHELL</td><td>2017-260 12.56</td><td>5.1</td><td>HR 1718</td><td>2017A037</td></tr><tr><td>29</td><td>NASA IRTF iSHELL</td><td>2018-287 12.53</td><td>4.1</td><td>HR 1061</td><td>2018A037</td></tr><tr><td>30</td><td>NASA IRTF iSHELL</td><td>2018-288 11.33</td><td>6.8</td><td>HR 1061</td><td>2018A037</td></tr><tr><td>31</td><td>NASA IRTF iSHELL</td><td>2018-312 09.36</td><td>6.7</td><td>HR 7981</td><td>2018A037</td></tr><tr><td>32</td><td>NASA IRTF iSHELL</td><td>2018-314 08.54</td><td>8.1</td><td>HR 7891</td><td>2018A037</td></tr></tbody></table>

<footnote><sup>a</sup>This dataset was published as reference spectra in the search for H<sub>3</sub><sup>+</sup> at Neptune [17].</footnote>
(a) NASA IRTF iSHELL

The NASA IRTF iSHELL instrument [18] is a high-resolution spectrograph that became operational in August 2016, with Uranus being the first science target after commissioning (programme 2016A041). Using the Lp3 setting and the 0.375" × 15" slit produces a cross dispersed spectra on a 2048 × 2048 pixel detector, at a spectral resolution of R ~ 70 000. The spectral range is 3.83–4.14 µm, giving a near complete coverage of the H<sub>3</sub><sup>+</sup> Q branch about 4 µm. Once extracted, a complete spectrum contains over 28 000 spectral pixels.

During these observations, the spectrograph slit was scanned across the disc of the planet, using three spatial positions, each separated by 0.8", centred about the centre of the disc. An example iSHELL spectrum can be seen in figure 2a, with the dashed lines indicating regions where the spectra have been truncated in order to emphasize the individual H<sub>3</sub><sup>+</sup> lines.

(b) Keck NIRSPEC

The Near-Infrared Spectrograph (NIRSPEC; [19]) is a medium- to high-resolution spectrograph mounted on the Keck II telescope on the summit of Mauna Kea in Hawaii, USA. Using an echelle angle of 62.02° and a cross disperser angle of 33.65°, with the 0.432" × 24" slit, produces a cross dispersed spectrum covering the H<sub>3</sub><sup>+</sup> Q branch region about 4 µm at a spectral resolution of R ~ 25 000. A NIRSPEC H<sub>3</sub><sup>+</sup> spectrum is shown in figure 2b.

(c) Gemini GNIRS

The Gemini Near InfraRed Spectrograph (GNIRS; [20]) is a long-slit medium resolution spectrograph mounted on Gemini North on Mauna Kea, Hawaii. Using the 32 lines mm<sup>-1</sup> grating and the 0.1" wide slit in the L band produces a spectrum with a spectral resolution of R ~ 5400. These observations were designed to be acquired during 'Band 3' conditions, which includes periods of significant cloud-cover and high water vapour content in Earth's atmosphere, unsuitable for the majority of science programmes. An example Gemini GNIRS spectrum is shown in figure 2c.

(d) NASA IRTF SpeX

The NASA IRTF SpeX instrument [21] is a medium resolution spectrograph, here operated at the setting LongXD1.9 giving a spectral resolution of R ~ 2500 with the 0.5" × 15" slit. The wavelength coverage encompasses the H<sub>3</sub><sup>+</sup> Q branch spectrum. The instrument was upgraded to a Teledyne H2RG detector array in August 2014, increasing the sensitivity and pixel density. These observations employed a single slit-position across the centre of the disc of Uranus, as indicated in figure 1b. An example SpeX H<sub>3</sub><sup>+</sup> spectrum is shown in figure 2d.

(e) VLT CRIRES

The CRIRES instrument [22] mounted on the VLT Unit Telescope 3 (UT3) is a long-slit spectrograph, capable of very high spectral resolution. The 0.2" × 25" slit and a central wavelength of 3.965 µm produces a single-order spectrum distributed across four Raytheon detectors, at a resolution of R ~ 100 000, covering the Q branch of H<sub>3</sub><sup>+</sup>. CRIRES was removed from UT3 in mid-2014 for a component upgrade, due to return to the telescope in late 2019. The high-resolution H<sub>3</sub><sup>+</sup> spectrum observed by CRIRES is shown in figure 2e, with the dashed lines separating the three very narrow spectral regions.
<figure><img src="image_4.png" /></figure>

**Figure 1.** The geometry of Uranus as seen from Earth in 1992, 2007, 2012 and 2018. The last uranian equinox was in 2007 and the next solstice will be in 2028. The rings are not visible at L' (approx. 4.0 µm) in the near-infrared. The orientation of the spectrograph slit is shown in orange for 2018, with the same orientation used for all the data used in this study.

<figure><img src="image_5.png" /></figure>

**Figure 2.** (a–e) Example Uranus $H_3^+$ spectra, one from each of the facilities and instruments used in this study. The very high-resolution spectra shown in (a) and (e) show the spectral regions containing three $H_3^+$ Q branch transitions as a function of spectral pixel due to the high spectral dispersion, where the dashed lines indicate a discontinuity in the wavelength coverage. (f) The dependence of the $H_3^+$ $Q(1, 0^-)/Q(3, 0^-)$ ratio on temperature.

## 3. Analysis

### (a) Long-term evolution of the temperature of the ionosphere

The first-order data product obtained from each night of observations is a spectrum of intensity versus wavelength, obtained in the L telluric window, between 3.3 and 4.1 µm (e.g. figure 2). At
Uranus, this wavelength region contains almost exclusively emission from $H_3^+$, since methane in the stratosphere absorbs any incident sunlight or thermal emission generated below the homopause. (This is in contrast to Jupiter and, particularly, Saturn.) This greatly simplifies the analysis of the observed spectrum. The L band window contains emission principally from the R and Q branch of $H_3^+$, regions with a rich history of being used as a probe for determining the temperatures of the upper atmospheres of the giant planets. In particular, the intensity ratio between the $H_3^+$ $Q(1,0^-)$ at 3.953 µm and the $Q(3,0^-)$ at 3.986 µm lines provides an excellent measure of temperature in the range of 400–800 K, as illustrated in figure 2f. This line-ratio is used to derive the temperature for the VLT CRIRES and NASA IRTF iSHELL observations in table 1, since the spectral resolution is high enough to resolve individual spectral lines. For the rest of the data, at lower spectral resolution, a full line-by-line $H_3^+$ model of the appropriate spectral region containing multiple lines is fitted to a temperature by applying Cramer’s rule (e.g. [14,23]). There is virtually no difference in the retrieved temperature between the two methods. The $H_3^+$ transitions are modelled using the line list of Neale et al. [24] and the partition function of Miller et al. [25], assuming conditions of *q*-LTE [24]. The error on the retrieved parameters is principally governed by the S/N (for more details, see [23]).

Once the temperature is determined from the shape of the $H_3^+$ spectrum, the number of emitting molecules is calculated by dividing the observed intensity of the $Q(1,0^-)$ line by the emission per molecule at the derived temperature. The total emission, which is the energy radiated over all wavelengths, is then calculated using the formulation of Miller et al. [26].

Table 2 shows the retrieved $H_3^+$ temperatures, $H_3^+$ densities, with associated errors, for each night of observation listed in table 1. For completeness, it also lists the calculated $H_3^+$ total emission, which can be an effective cooling mechanism for the upper atmosphere [25]. The temperature for each observing night is plotted in figure 3. The maximum measured temperature is 612 K in 2014 and the minimum is 429 K in 2017 with a mean of 507 K. The standard deviation is 44 K.

There is significant variability in the temperatures in figure 3. In order to investigate how these observations fit in with the long-term observations [13,14], a yearly temperature average is calculated from the observations contained in table 2. These averages are listed in table 3 and are shown together with the yearly averaged temperatures of previous studies in figure 4. The error bars represent the standard deviation of the observations within a particular year (grey), unless there is only one data-point, in which case the error bar is the uncertainty on the temperature (black). A linear fit to the globally averaged temperature between 1992 and 2018 shows that the ionosphere has continued to cool at a steady rate of $8 \pm 1$ K yr$^{-1}$—the same rate as derived by the initial long-term analysis [13].

## (b) Local-time profiles of the ionosphere

With new instruments and telescopes providing greater sensitivity and finer spatial resolution, we are able to move from investigating global properties of the ionosphere to spatially resolving the disc of Uranus. Another important difference between the observations obtained for this study and previous studies [13,14] is that the slit was aligned with the equator, as illustrated in figure 1d, as opposed to along the rotational axis. This provides a local-time view of the $H_3^+$ ionosphere, from dusk, across noon, to dawn. Note that these profiles cut across a range of latitudes and longitudes as illustrated in figure 1b.

Close to dawn and dusk, the line-of-sight vector traverses a long pathway in the ionosphere, due to the slant viewing angle. This has the effect of increasing the intensity of the observed $H_3^+$ emission. In order to convert these to represent the intensity observed along the 'surface' normal, a line-of-sight correction is applied (e.g. [11]).

Figure 5 shows the line-of-sight corrected local-time $H_3^+$ $Q(1,0^-)$ intensity profile for four different instruments, as indicated by the observation ID (table 1). The four different instruments produce profiles with vastly different S/N, with NASA IRTF SpeX providing the lowest, and
<table><caption>**Table 2.** The retrieved global $H_3^+$ temperature ($T$), temperature error ($\Delta T$), column density ($N$, in units of $10^{15} m^{-2}$), column density error ($\Delta N$) and the wavelength integrated $H_3^+$ emission ($E$, in units of $\mu\text{Wm}^{-2}\text{sr}^{-1}$).</caption><thead><tr><th>ID</th><th>DOY</th><th>T (K)</th><th>ΔT</th><th>N<sup>1</sup></th><th>ΔN</th><th>E</th></tr></thead><tbody><tr><td>1</td><td>2012-229</td><td>575</td><td>15</td><td>2.2</td><td>0.4</td><td>3.9</td></tr><tr><td>2</td><td>2012-321</td><td>519</td><td>12</td><td>11.6</td><td>2.2</td><td>10.3</td></tr><tr><td>3</td><td>2012-348</td><td>446</td><td>7</td><td>1.7</td><td>0.2</td><td>0.5</td></tr><tr><td>4</td><td>2013-215</td><td>505</td><td>15</td><td>9.7</td><td>2.3</td><td>7.1</td></tr><tr><td>5</td><td>2013-216</td><td>501</td><td>15</td><td>6.4</td><td>1.5</td><td>4.5</td></tr><tr><td>6</td><td>2013-295</td><td>491</td><td>9</td><td>3.0</td><td>0.4</td><td>1.8</td></tr><tr><td>7</td><td>2013-296</td><td>503</td><td>20</td><td>3.8</td><td>1.3</td><td>2.7</td></tr><tr><td>8</td><td>2013-297</td><td>600</td><td>35</td><td>2.1</td><td>0.6</td><td>4.9</td></tr><tr><td>9</td><td>2013-300</td><td>537</td><td>12</td><td>12</td><td>2</td><td>10.1</td></tr><tr><td>10</td><td>2013-319</td><td>486</td><td>12</td><td>14.3</td><td>3.0</td><td>7.9</td></tr><tr><td>11</td><td>2013-340</td><td>471</td><td>42</td><td>8.7</td><td>5.7</td><td>3.8</td></tr><tr><td>12</td><td>2013-359</td><td>454</td><td>34</td><td>6.0</td><td>3.5</td><td>1.9</td></tr><tr><td>13</td><td>2014-255</td><td>612</td><td>21</td><td>1.1</td><td>0.3</td><td>3.0</td></tr><tr><td>14</td><td>2014-256</td><td>559</td><td>12</td><td>2.5</td><td>0.4</td><td>3.7</td></tr><tr><td>15</td><td>2014-285</td><td>528</td><td>6</td><td>1.0</td><td>0.1</td><td>1.0</td></tr><tr><td>16</td><td>2014-286</td><td>538</td><td>10</td><td>7.9</td><td>1.1</td><td>9.0</td></tr><tr><td>17</td><td>2014-328</td><td>574</td><td>32</td><td>3.3</td><td>1.3</td><td>5.9</td></tr><tr><td>18</td><td>2014-337</td><td>513</td><td>14</td><td>3.5</td><td>0.8</td><td>2.9</td></tr><tr><td>19</td><td>2015-325</td><td>540</td><td>15</td><td>1.8</td><td>0.4</td><td>2.1</td></tr><tr><td>20</td><td>2016-230</td><td>482</td><td>5</td><td>6.8</td><td>2.0</td><td>3.5</td></tr><tr><td>21</td><td>2016-284</td><td>474</td><td>5</td><td>6.3</td><td>1.9</td><td>2.9</td></tr><tr><td>22</td><td>2016-285</td><td>459</td><td>4</td><td>8.6</td><td>2.6</td><td>3.1</td></tr><tr><td>23</td><td>2016-320</td><td>479</td><td>2</td><td>6.2</td><td>1.9</td><td>3.1</td></tr><tr><td>24</td><td>2017-244</td><td>447</td><td>7</td><td>9.3</td><td>2.8</td><td>2.7</td></tr><tr><td>25</td><td>2017-245</td><td>429</td><td>12</td><td>5.0</td><td>1.5</td><td>1.0</td></tr><tr><td>26</td><td>2017-258</td><td>536</td><td>6</td><td>1.5</td><td>0.4</td><td>1.7</td></tr><tr><td>27</td><td>2017-259</td><td>495</td><td>8</td><td>2.0</td><td>0.6</td><td>1.3</td></tr><tr><td>28</td><td>2017-260</td><td>530</td><td>9</td><td>1.7</td><td>0.5</td><td>1.8</td></tr><tr><td>29</td><td>2018-287</td><td>483</td><td>7</td><td>4.2</td><td>1.3</td><td>2.2</td></tr><tr><td>30</td><td>2018-288</td><td>499</td><td>9</td><td>4.6</td><td>1.4</td><td>3.1</td></tr><tr><td>31</td><td>2018-312</td><td>495</td><td>7</td><td>2.2</td><td>0.6</td><td>1.4</td></tr><tr><td>32</td><td>2018-314</td><td>466</td><td>10</td><td>2.7</td><td>0.8</td><td>1.1</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** The $H_3^+$ temperatures derived for individual observing nights listed in table 1. There is significant variability in temperature between each individual observation. The maximum measured temperature is 612 K in 2014 and the minimum is 429 K in 2017 with a mean of 507 K. The standard deviation is 44 K.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** The yearly average temperature of the upper atmosphere of Uranus as derived from the $H_3^+$ observations, between 1992 and 2018, showing a long-term cooling trend. The dashed line is a linear fit to the data, with a slope of $8 \pm 1$ K yr$^{-1}$. The plotted temperatures are listed in table 3. The grey error bars indicate years with multiple observations, with black error bars indicating single measurements.</figcaption></figure>

Keck NIRSPEC the highest; the small-scale variability across the disc should not be considered to be real, and instead just indicative of S/N. For the four observations shown in figure 5 the $H_3^+$ intensity peaks between a local time of about 13 and 15. The morning sector is consistently less intense than the evening sector, and the fall-off of intensity outside the $\pm 1$ planetary radius is caused by the telluric seeing.
<table><caption>**Table 3.** The yearly averaged global parameters, calculated from table 2, listing $H_3^+$ temperature (*T*), standard deviation on the temperature (*σ* *T*), column density (*N*, in units of $10^{15} m^{-2}$), the standard deviation on the column densities (*σ* *N*) and the wavelength integrated $H_3^+$ emission (*E*, in units of $\mu\mu Wm^{-2} sr^{-1}$).</caption><thead><tr><th>year</th><th><i>T</i> (K)</th><th>σ (<i>T</i>)</th><th><i>N</i></th><th>σ (<i>N</i>)</th><th><i>E</i></th></tr></thead><tbody><tr><td>1992.3</td><td>715</td><td>47</td><td>1.4</td><td>0.3</td><td>9.0</td></tr><tr><td>1993.3</td><td>705</td><td>41</td><td>1.4</td><td>0.4</td><td>8.4</td></tr><tr><td>1994.5</td><td>621</td><td>14</td><td>2.9</td><td>0.2</td><td>8.3</td></tr><tr><td>1995.5</td><td>674</td><td>51</td><td>1.2</td><td>0.4</td><td>5.6</td></tr><tr><td>1999.7</td><td>615</td><td>21</td><td>3.9</td><td>1.6</td><td>10.6</td></tr><tr><td>2000.7</td><td>584</td><td>24</td><td>3.4</td><td>0.0</td><td>6.7</td></tr><tr><td>2001.6</td><td>650</td><td>31</td><td>1.8</td><td>0.5</td><td>6.6</td></tr><tr><td>2002.6</td><td>599</td><td>34</td><td>1.9</td><td>0.5</td><td>4.4</td></tr><tr><td>2006.7</td><td>608</td><td>40</td><td>1.4</td><td>0.8</td><td>3.5</td></tr><tr><td>2008.8</td><td>534</td><td>39</td><td>1.6</td><td>0.5</td><td>1.7</td></tr><tr><td>2011.7</td><td>520</td><td>32</td><td>3.6</td><td>0.3</td><td>3.2</td></tr><tr><td>2012.8</td><td>513</td><td>65</td><td>5.2</td><td>0.6</td><td>4.2</td></tr><tr><td>2013.8</td><td>505</td><td>42</td><td>7.0</td><td>5.6</td><td>5.1</td></tr><tr><td>2014.8</td><td>554</td><td>36</td><td>3.2</td><td>3.9</td><td>4.5</td></tr><tr><td>2015.9</td><td>540</td><td>15</td><td>1.8</td><td>2.5</td><td>2.1</td></tr><tr><td>2016.8</td><td>474</td><td>10</td><td>7.0</td><td>0.4</td><td>3.1</td></tr><tr><td>2017.7</td><td>487</td><td>48</td><td>3.9</td><td>1.1</td><td>2.2</td></tr><tr><td>2018.8</td><td>486</td><td>15</td><td>3.4</td><td>3.3</td><td>1.9</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Four local-time $H_3^+$ profiles of Uranus as a function of planetary radius, with the ID of the observation used is indicated. For each observations, the slit was aligned with the equator of the planet. The red lines show the local time across the disk of the planet. The profiles are broadly consistent, peaking in the afternoon sector, between 13.00 and 15.00. The small-scale variability across the disk of Uranus are indicative of the S/N, with uncertainties indicated by the error bars on the left, and are not to be considered real.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** (a) Two local-time profiles from the NASA IRTF iSHELL observations obtained on 2016-285 (#22 in table 1), separated by 2.5 h. A clear intensity feature rotates into view. Note that these profiles have not been line-of-sight corrected, effectively enhancing the limb emissions. (b,c) The geometry of Uranus on 2016-285, with the orientation of the slit shaded in orange. The Voyager 2 auroral observations of Herbert [5] are mapped onto the planet, using a longitude shift consistent with the observations in (a).</figcaption></figure>

(c) The event of 2016-285

The typical ionospheric local-time profile described in the previous sections are present in the observations with good telluric seeing conditions and stable telescope guiding. However, there is one outlier to this relatively stable ionospheric profile: the NASA IRTF iSHELL observations of 2016-285 (#21 in table 1).

Figure 6a shows two local-time profiles of $H_3^+$ intensity obtained on 2016-285 (#21), with the slit offset $0.8''$ south from the centre of the disc, as illustrated in figure 6b,c. The two observations are separated by 2.5 h, and the profiles have not been line-of-sight corrected to better reveal emissions close to the limb. The profile obtained at 08.58 UTC is similar to the profiles shown in figure 5, increasing towards the dusk sector. The second dashed profile, with a mid-observation time of 11.27 UTC, shows a strong enhancement in the observed $H_3^+$ intensity on the dawn side. Each profile has an integration time of 30 min.

The dashed profile in figure 6a appears to be subject to stronger line-of-sight enhancement, which could be caused by a decrease in the telluric seeing, effectively producing a sharper intensity along the sit. The dashed line also appears less intense than the earlier solid intensity profile at local times between 12 and 16 h. This could in part be due to the different seeing conditions experienced during the two observations, but could also be indicative of a $H_3^+$ intensity feature fixed in local time rotating from the afternoon sector (solid line) to the dusk limb (dashed line), producing an increase in the line-of-sight brightening.

The dawn feature seen in the dashed line is consistent with a region of bright $H_3^+$ emission rotating in from the night-side, appearing on the limb, being subjected to strong line-of-sight brightening. At 2016-285 11.27, the sub-observer ULS longitude is $352^\circ$ which means that a longitude on the dawn limb is $82^\circ$. The exact location of the auroral emission is subject to errors governed by the spatial resolution (approx. $0.7''$), total integration time of the co-added observations, and the unknown spatial extension of the feature, placing the emission at a longitude of $82^\circ \pm 30^\circ$. The latitude of the $H_3^+$ intensity enhancement is $-45^\circ \pm 15^\circ$, consistent with the southern auroral region [8].

Unfortunately, due to observing time constraints, the potential auroral feature in figure 6 could not be tracked across the disc of Uranus, and the observation on 2016-285 remains
<figure><img src="image_10.png" /><figcaption>**Figure 7.** The long-term temperature trend compared to the geometric and magnetic season described in the text.</figcaption></figure>

the only observation with a possible auroral signal present in the entire dataset listed in table 1.

# 4. Discussion

The analysis of the observations of $H_3^+$ outlined above have revealed three important features of the upper atmosphere of Uranus. Firstly, the upper atmosphere has continued to cool since $H_3^+$ was discovered at Uranus in 1992 [12]. Secondly, spatially resolved observations reveal that the peak $H_3^+$ intensity occurs between 13 and 15 in local time. Thirdly, on 2016-285 a dawn brightening appeared in the $H_3^+$ intensity, potentially providing the first detection of $H_3^+$ auroral emission. These results are discussed in greater detail below.

## (a) Long-term cooling in the upper atmosphere of Uranus

The long-term cooling of the upper atmosphere of Uranus discovered by Melin et al. [13] has continued throughout the new set of observations analysed in this study. Figure 4 shows a remarkable persistent cooling, from 1992 to 2018, with a linear gradient of $8 \pm 1$ K yr$^{-1}$. This 27 year interval is longer than the length of a nominal season at Uranus ($84/4 = 21$ yr), supporting the idea that the heating mechanism, and ultimately the solution to the ‘energy crisis’, is not in any way related to the deposition of solar photons. In addition, neither the passing of equinox (2007) nor aphelion (2009) had any effect on the observed temperatures.

Owing to the extreme seasons of Uranus, at solstice one pole is in permanent sunlight, while the other is in complete darkness. If the temperature of the ionosphere is governed principally by solar irradiation, then the planet, as observed from Earth, would be subject to two hot solstices and two cool equinoxes. This ‘geometric season’ is illustrated as the red line in figure 7.

Long-term changes in the reflectivity of the troposphere [27] have been linked to the extreme seasons that the planet experiences during its 84 year orbit around the Sun, modulated by the 11 year solar cycle. This suggests that the changing levels of ultraviolet irradiation, or changing rates of solar-modulated galactic cosmic rays, change the photochemistry in the atmosphere. In principle, a similar modulation should be observed in the ionosphere, as the changing solar extreme ultraviolet (EUV) flux throughout the solar cycle would give rise to changing levels of ionization of molecular hydrogen, producing different amounts of $H_3^+$. However, there is no clear
correlation between solar cycle and the observed $H_3^+$ column density in the data analysed here (as also noted by Melin *et al.* [13]).

The potential inability for Uranus to build up significant magnetic flux in the tail during equinox conditions compared to solstice [9] provides a mechanism with which the auroral process is modulated by season. Since the magnetic field is both tilted and offset from the centre of the planet, the two solstices will not have the same magnetic configuration with respect to the IMF. This in turn may produce conditions favourable for sustaining currents that drive Joule heating at one solstice, but perhaps not the other, producing one hot and one cold solstice. This 'magnetic season' is plotted as the orange line in figure 7. Disentangling these two hypotheses, geometric or magnetic season, requires ongoing observations of the temperature of Uranus' upper atmosphere. Note that both the amplitude and the phase of the solid lines in figure 7 are determined by eye to provide an approximate fit to the data, and serves to illustrate the concept of 'geometric' and 'magnetic' season.

While the tropospheric reflectance changes as a function of solar cycle, there is no evidence for any long-term changes in the tropospheric temperature [28–30], suggesting that despite being subject to extreme seasons, the thermal effect on the troposphere, as observed from the Earth, is negligible. This is in stark contrast to the long-term cooling observed in the upper atmosphere.

The short-term variability in the temperatures in figure 3 is of the order of approximately 100 K over days. Since the heat capacity of the upper atmosphere of Uranus is likely very large ($15\,000\,\text{Jkg}^{-1}\,\text{K}^{-1}$ for Saturn, [31]), it is unlikely that the variability is due to actual bulk temperature changes. Instead, we expect any short-term changes in temperature to be driven by differences in the peak altitude of $H_3^+$ production, sampling different regions of the thermospheric temperature curve, where low altitudes correspond to lower temperatures, and vice versa [32]. This means that $H_3^+$ is probably predominantly produced along the steep slope of the temperature profile at lower altitudes, and that relatively modest changes in the ionizing energy can produce large differences in the observed $H_3^+$ temperature. In figure 3, there are 2 years where the temperature variability is small, 2016 and 2018, suggesting that the ionization energy remained relatively constant, which may be related to this being the descending phase of solar cycle. Overall, however, given the large numbers of observations used for this study, we expect the yearly averaged temperatures to be representative of the average thermospheric temperature at a fixed altitude, such that the long-term cooling is physical.

## (b) Possible effect of tropospheric storms on the upper atmosphere

Most of the temperatures observed in 2014 are hotter than those observed in 2013, with an average of $505 \pm 42$ K in 2013 and $554 \pm 36$ K in 2014. This may be indicative of an intermittent heating event. The uptick of the average global temperature in 2014 is distinct, and coincides with the outbreak of a tropospheric storm in August 2014, fierce enough to punch through the stratospheric methane haze layer [33]. Analogous to the heating that was observed above the Great Red Spot (GRS) of Jupiter [15], this turbulent storm could generate acoustic or gravity waves that traverse vertically up through the atmosphere, breaking and depositing energy in the upper atmosphere, heating this region in the process. At Jupiter, however, the observed heating was highly localized in both longitude and latitude, confined to the region immediately above the GRS. The observations of Uranus presented here cannot spatially resolve the precise location of the heating.

There are three ways in which a tropospheric storm protruding into the stratosphere could produce an apparent heating in the upper atmosphere. Firstly, the storm could generate acoustic or gravity waves that propagate up into the upper atmosphere where they break and deposit their energy. Secondly, an upwelling of stratospheric hydrocarbons into the upper atmosphere provides an efficient sink for $H_3^+$ ions, since hydrocarbons quickly destroy $H_3^+$. Because this would predominantly affect regions close to the homopause, where the temperature is relatively low, an observation of the entire column of $H_3^+$ would therefore exclude the cold population at low altitudes, and therefore result in a higher observed temperature. Thirdly, the breaking of waves
generated in the lower atmosphere induce a vertical wind shear in the thermosphere, producing a jagged appearance of the electron density altitude profile, in which case the effective $H_3^+$ peak can be shifted to higher altitudes [34], where it is hotter, producing a higher observed $H_3^+$ temperature.

Heating by tropospheric storms may explain the fact that the temperature is hotter in 2014 than in 2013, while the column density is lower. If the explanation of advecting methane into the upper atmosphere is valid, the particularly low density of $H_3^+$ noted in 2015 might suggest that this methane remained, or was replenished, in the upper atmosphere for at least a year.

### (c) Local-time profiles of the ionosphere

Figure 5 shows four typical local-time profiles contained within the dataset analysed here, with the peak intensity occurring at local times between 13 and 15. The peak ionization rate is expected to occur at local noon, and smaller solar-zenith angles produce higher ionization rates. Unless the temperature changes dramatically across local times, then the location of the observed intensity maximum is also the location of the peak column integrated $H_3^+$ density. It is currently unclear how to shift the peak density away from noon, where the production of $H_3^+$ peaks, but the results obtained here do serve as important input to models of Uranus' upper atmosphere.

Trafton *et al.* [35] observed the intensity of $H_3^+$ across the disc of the planet with the slit aligned along the rotational axis, along local noon, with the maximum intensity occurring at the centre of the disc. By contrast, they also observed the quadrupole $H_2$ emission in the K telluric window at 2 µm, which appeared to have a line-of-sight enhanced signal, being bright at both limbs. Here, while using a slit aligned along the equator of Uranus, the $H_3^+$ is observed to be subject to strong line-of-sight effects (e.g. figure 6) in a similar manner to the $H_2$ emissions.

Imaging performed with the Gemini telescope in 2017 using a Br-α filter (4.05–4.10 µm) showed the $H_3^+$ intensity across the peak had a spatial distribution similar to the stratospheric haze [36], with a bright cap above the southern polar region [37], with local-time profiles similar to the ones shown in figure 5. The upwelling of stratospheric hydrocarbons at low latitudes and downwelling at high latitudes [38] could in principle remove a sink of $H_3^+$ at high latitudes, generating brighter emission due to an enhanced $H_3^+$ density. However, this should result in higher $H_3^+$ column densities about solstice, since polar latitudes are predominantly visible from the Earth. There is no seasonal variability, or any clear long-term trends seen the $H_3^+$ densities listed in table 2.

### (d) The dawn $H_3^+$ feature on 2016-285

Figure 6a shows a bright $H_3^+$ intensity feature rotate into view on the dawn terminator of the planet. The intensity of $H_3^+$ is driven linearly by density and exponentially by temperature, so an intensity enhancement requires either an increased ionization of $H_2$, which very quickly produces $H_3^+$, or a localized increase in temperature. Ionization by solar EUV photons has a well-defined local-time profile, shown in figure 5, and cannot produce a strong enhancement at dawn. The auroral process, however, drives both ionization and heating in the upper atmosphere, and can produce auroral signatures about the magnetic poles throughout all local times [5,8], including the night-side.

The HST observations of Lamy *et al.* [7,8,36] show that the total radiated power of Uranus' aurora increases over time between 2011 and 2017, suggesting an increase in the electron energy flux. The simultaneous HST and Gemini $H_3^+$ imaging observations obtained in 2017 [36] showed that while an auroral feature was clearly seen in the ultraviolet, a corresponding brightening was not observed in $H_3^+$. Since $H_3^+$ emission is strongly dependent on the temperature of the upper atmosphere, any density enhancements of $H_3^+$ generated by the auroral process will have a lower contrast when the upper atmosphere is cold (e.g. 2017), compared to when it is observed to be hot (e.g. 1992).

The ULS longitude of the auroral feature is 82°±30°. This location is broadly consistent with the HST observations of auroral emission in the ultraviolet obtained in 2012 by Lamy *et al.* [7],
but are different by about 150° longitude with the 2014 observations [8]. However, the rotation period of Uranus has a relatively large uncertainty, which means that over the course of a year an uncertainty of 106° has accumulated, and the absolute longitude is effectively lost. Future spatially resolved, high sensitivity observations of $H_3^+$ auroral observations can be used to re-establish the ULS.

The solar wind propagation model of Tao *et al.* [39] predicts a spike in the solar wind dynamic pressure of approximately 0.03 nPa on 2016-292, 7 days after our observation. The magnitude of this compression is similar to that reported by Lamy *et al.* [7]. At 19 AU, the model is subject to large errors in the arrival time of these enhancements, about ±3 days, so the observed $H_3^+$ brightening is unlikely to be linked to a solar wind compression. Overall, over 60 h of NASA IRTF iSHELL data was obtained for this study, with only about 2 h showing this clear dawn enhancement, yielding a probability of observing such an event approximately 3%, rendering it very rare indeed.

The fidelity of the data does not allow for an accurate determination of temperature of the dawn feature, but if the dawn $H_3^+$ enhancement is driven by temperature, then this 'hotspot' could be produced by a storm in the lower atmosphere, dumping heat into the upper atmosphere. However, no contemporaneous reports of a tropospheric storm have been communicated.

The apparent depletion of $H_3^+$ intensity at noon of the dashed profile in figure 6 may also be interpreted as a region of low intensity $H_3^+$ rotating from dawn in the dashed profile to about noon in the dashed profile. This can either be in the form of a cold region fixed in local time or a density depletion.

More broadly, the short-term variability in both $H_3^+$ column density and temperature could be indicative of auroral features, which are almost always present and varying, but cannot be resolved due to the limited spatial resolution of most of the observations in table 1. The early observations of Lam *et al.* [40] detected variations in the $H_3^+$ intensity of about 20%, while the temperatures in table 2 have an average standard deviation of 7% of the value derived for the temperature, ranging from 2 to 13%. The probability of observing auroral features along the slit is also dependent on the slit orientation, with the slit aligned with the rotational axis producing a higher contrast $H_3^+$ signal, while if the slit aligned with the equator has any $H_3^+$ enhancements rotating along the slit, the observed contrast will be lower.

## 5. Conclusion

This study has revealed that the upper atmosphere of Uranus has cooled consistently between 1992 and 2018, at a rate of 8 ± 1 K yr⁻¹. Since this is longer than the 21 year seasonal cycle, it may be linked to auroral Joule heating being modulated by an offset and asymmetric magnetic field, making the northern summer solstice hot and the southern winter solstice cold. However, further observations are required to confirm this hypothesis.

In 2014, during the outbreak of the largest tropospheric storm observed at Uranus in over a decade, the temperatures of the upper atmosphere were observed to be hotter than the previous year. This is suggestive of heating associated with acoustic or gravity waves generated by the turbulent troposphere, as previously observed on Jupiter [15].

The NASA IRTF iSHELL observations of 2016-285 revealed a bright $H_3^+$ intensity enhancement rotating into view on the dawn limb. This feature appears at latitudes consistent with the southern auroral region observed by Voyager 2 [5], and may be auroral in nature. However, since other hypotheses cannot be unambiguously ruled out, the origin of this feature remains ambiguous.

Looking towards the next decade, the James Webb Space Telescope (JWST) promises to transform our understanding of the upper atmosphere of Uranus, its coupling to the magnetic field, and the extent to which it is heated by waves generated in the lower atmosphere. The NIRSPEC instrument [41] on-board the JWST has a two-dimensional field-of-view of 3 × 3", with each pixel covering 0.1 × 0.1". This enables us to capture almost the entire disc of Uranus in one exposure, providing unrivalled spatial resolution and sensitivity. As part of the Guaranteed Time
Observing programme, a full map of $H_3^+$ emissions across the disc of Uranus will be obtained in the early science phase of the telescope.

Despite 27 years of observations of $H_3^+$ emission from Uranus, our understanding of its upper atmosphere remains in its infancy, and further observations and subsequent modelling are required to form a basic understanding of this enigmatic system. In addition, Uranus is proving a very worthy target for future spacecraft missions.

**Data accessibility.** The data are available in public telescope archives, or from H.M. upon request.

**Competing interests.** We declare we have no competing interests.

**Funding.** H.M. was supported by the STFC grant no. ST/N000749/1. L.L. was supported by CNES and CNRS/INSU programs of planetology and heliophysics.

**Acknowledgements.** H.M. is a Visiting Astronomer at the IRTF, which is operated by the University of Hawaii under contract NNH14CK55B with the National Aeronautics and Space Administration. Part of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. Partly based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere. Partly based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), Ministério da Ciência, Tecnologia e Inovação (Brazil) and Korea Astronomy and Space Science Institute (Republic of Korea). The individual programme identification codes for each respective telescope are listed in table 1.

## References

1. Herschel W. 1781 Account of a comet. *Phil. Trans. R. Soc.* **71**, 492–501.

2. Herschel W. 1787 An account of the discovery of two satellites revolving round the georgian planet. *Phil. Trans. R. Soc.* **77**, 125–129.

3. Ness NF, Acuna MH, Behannon KW, Burlaga LF, Connerney JEP, Lepping RP, Neubauer FM. 1986 Magnetic fields at Uranus. *Science* **233**, 85–89. (doi:10.1126/science.233.4759.85)

4. Desch MD, Connerney JEP, Kaiser ML. 1986 The rotation period of Uranus. *Nature* **322**, 42–43. (doi:10.1038/322042a0)

5. Herbert F. 2009 Aurora and magnetic field of Uranus. *J. Geophys. Res.* **114**, A11206. (doi:10.1029/2009JA014394)

6. Lamy L *et al*. 2018 Saturn's northern aurorae at solstice from HST observations coordinated with Cassini's grand finale. *Geophys. Res. Lett.* **45**, 9353–9362. (doi:10.1029/2018GL078211)

7. Lamy L *et al*. 2012 Earth-based detection of Uranus' aurorae. *Geophys. Res. Lett.* **39**, L07105. (doi:10.1029/2012GL051312)

8. Lamy L *et al*. 2017 The aurorae of Uranus past equinox. *J. Geophys. Res.* **122**, 3997–4008. (doi:10.1002/2017JA023918)

9. Cowley SWH. 2013 Response of Uranus' auroras to solar wind compressions at equinox. *J. Geophys. Res.* **118**, 2897–2902. (doi:10.1002/jgra.50323)

10. Masters A. 2014 Magnetic reconnection at Uranus' magnetopause. *J. Geophys. Res.* **119**, 5520–5538. (doi:10.1002/2014JA020077)

11. Johnson RE, Melin H, Stallard TS, Tao C, Nichols JD, Chowdhury MN. 2018 Mapping $H_3^+$ temperatures in Jupiter's northern auroral ionosphere using VLT-CRIRES. *J. Geophys. Res.* **123**, 5990–6008. (doi:10.1029/2018JA025511)

12. Trafton LM, Geballe TR, Miller S, Tennyson J, Ballester GE. 1993 Detection of $H_3^+$ from uranium. *Astrophys. J.* **405**, 761–766. (doi:10.1086/172404)

13. Melin H, Stallard T, Miller S, Trafton LM, Encrenaz T, Geballe TR. 2011 Seasonal variability in the ionosphere of uranium. *Astron. J.* **729**, 134. (doi:10.1088/0004-637X/729/2/134)

14. Melin H, Stallard TS, Miller S, Geballe TR, Trafton LM, O'Donoghue J. 2013 Post-equinoctial observations of the ionosphere of Uranus. *Icarus* **223**, 741–748. (doi:10.1016/j.icarus.2013.01.012)
15. O'Donoghue J, Moore L, Stallard TS, Melin H. 2016 Heating of Jupiter's upper atmosphere above the Great Red Spot. *Nature* **536**, 190-192. (doi:10.1038/nature18940)

16. Yates JN, Achilleos N, Guio P. 2014 Response of the Jovian thermosphere to a transient 'pulse' in solar wind pressure. *Planet. Space Sci.* **91**, 27-44. (doi:10.1016/j.pss.2013.11.009)

17. Melin H, Fletcher LN, Stallard TS, Johnson RE, O'Donoghue J, Moore L, Donnelly PT. 2018 The quest for $H_3^+$ at Neptune: deep burn observations with NASA IRTF iSHELL. *Mon. Not. R. Astron. Soc.* **474**, 3714-3719. (doi:10.1093/mnras/stx3029)

18. Rayner J, Bond T, Bonnet M, Jaffe D, Muller G, Tokunaga A. 2012 iSHELL: a 1-5 micron cross-dispersed R = 70 000 immersion grating spectrograph for IRTF. In *Ground-based and Airborne Instrumentation for Astronomy IV*. Proceedings of the SPIE, vol. 8446, p. 84462C. Amsterdam, The Netherlands: SPIE Astronomical Telescopes + Instrumentation.

19. McLean IS et al. 1998 Design and development of NIRSPEC: a near-infrared echelle spectrograph for the Keck II telescope. In *Society of Photo-Optical Instrumentation Engineers (SPIE) Conf. Series* (ed. AM Fowler). Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, vol. 3354, pp. 566-578. Bellingham, WA: SPIE.

20. Elias JH, Joyce RR, Liang M, Muller GP, Hileman EA, George JR. 2006 *Design of the Gemini near-infrared spectrograph*. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 6269, p. 4. Bellingham, WA: SPIE.

21. Rayner JT, Toomey DW, Onaka PM, Denault AJ, Stahlberger WE, Vacca WD, Cushing MC, Wang S. 2003 SpeX: a medium-resolution 0.8-5.5 micron spectrograph and imager for the NASA infrared telescope facility. *Publ. Astron. Soc. Pac.* **115**, 362-382. (doi:10.1086/pasp.2003.115 issue-805)

22. Kaeufl H-U et al. 2004 CRIRES: a high-resolution infrared spectrograph for ESO's VLT. In *Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series* (eds AFM Moorwood, M Iye). Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 5492, pp. 1218-1227. Bellingham, WA: SPIE.

23. Melin H, Stallard TS, O'Donoghue J, Badman SV, Miller S, Blake JSD. 2014 On the anticorrelation between $H_3^+$ temperature and density in giant planet ionospheres. *Mon. Not. R. Astron. Soc.* **438**, 1611-1617. (doi:10.1093/mnras/stt2299)

24. Neale L, Miller S, Tennyson J. 1996 Spectroscopic properties of the $H_3^+$ molecule: a new calculated line list. *Astrophys. J.* **464**, 516-520. (doi:10.1086/177341)

25. Miller S, Stallard T, Melin H, Tennyson J. 2010 $H_3^+$ cooling in planetary atmospheres. *Faraday Discuss.* **147**, 283. (doi:10.1039/c004152c)

26. Miller S, Stallard T, Tennyson J, Melin H. 2013 Cooling by $H_3^+$ emission. *J. Phys. Chem. A* **117**, 9770-9777. (doi:10.1021/jp312468b)

27. Aplin KL, Harrison RG. 2017 Solar-driven variation in the atmosphere of uranus. *Geophys. Res. Lett.* **44**, 12. (doi:10.1002/2017GL075374)

28. Orton GS et al. 2014 Mid-infrared spectroscopy of Uranus from the Spitzer Infrared Spectrometer. 1. Determination of the mean temperature structure of the upper troposphere and stratosphere. *Icarus* **243**, 494-513. (doi:10.1016/j.icarus.2014.07.010)

29. Orton GS et al. 2014 Mid-infrared spectroscopy of Uranus from the Spitzer infrared spectrometer: 2. Determination of the mean composition of the upper troposphere and stratosphere. *Icarus* **243**, 471-493. (doi:10.1016/j.icarus.2014.07.012)

30. Orton GS, Fletcher LN, Encrenaz T, Leyrat C, Roe HG, Fujiyoshi T, Pantin E. 2015 Thermal imaging of uranus: upper-tropospheric temperatures one season after Voyager. *Icarus* **260**, 94-102. (doi:10.1016/j.icarus.2015.07.004)

31. Müller-Wodarg ICF, Mendillo M, Yelle RV, Aylward AD. 2006 A global circulation model of Saturn's thermosphere. *Icarus* **180**, 147-160. (doi:10.1016/j.icarus.2005.09.002)

32. Broadfoot AL et al. 1981 Extreme ultraviolet observations from Voyager 1 encounter with Saturn. *Science* **212**, 206-211. (doi:10.1126/science.212.4491.206)

33. de Pater I, Sromovsky LA, Fry PM, Hammel HB, Baranec C, Sayanagi KM. 2015 Record-breaking storm activity on Uranus in 2014. *Icarus* **252**, 121-128. (doi:10.1016/j.icarus.2014.12.037)

34. Barrow D, Matcheva KI, Drossart P. 2012 Prospects for observing atmospheric gravity waves in Jupiter's thermosphere using $H_3^+$ emission. *Icarus* **219**, 77-85. (doi:10.1016/j.icarus.2012.02.007)
35. Trafton LM, Miller S, Geballe TR, Tennyson J, Ballester GE. 1999 H₂ quadrupole and H₃⁺ emission from uranus: the uranian thermosphere, ionosphere, and aurora. *Astrophys. J.* **524**, 1059–1083. (doi:10.1086/apj.1999.524 issue-2)

36. Lamy L et al. 2018 Analysis of HST, VLT and Gemini coordinated observations of Uranus late 2017: a multi-spectral search for auroral signatures. In *SF2A-2018: Proc. of the Annual meeting of the French Society of Astronomy and Astrophysics*, pp. 29–32. Montpellier, France: French Society of Astronomy & Astrophysics.

37. Sromovsky LA, Karkoschka E, Fry PM, de Pater I, Hammel HB. 2019 The methane distribution and polar brightening on Uranus based on HST/STIS, Keck/NIRC2, and IRTF/SpeX observations through 201. *Icarus* **317**, 266–306. (doi:10.1016/j.icarus.2018.06.026)

38. Toledo D, Irwin PGJ, Teanby NA, Simon AA, Wong MH, Orton GS. 2018 Uranus's northern polar cap in 2014. *Geophys. Res. Lett.* **45**, 5329–5335. (doi:10.1029/2018GL077654)

39. Tao C, Kataoka R, Fukunishi H, Takahashi Y, Yokoyama T. 2005 Magnetic field variations in the Jovian magnetotail induced by solar wind dynamic pressure enhancements. *J. Geophys. Res. Space Phys.* **110**, A11208. (doi:10.1029/2004JA010959)

40. Lam HA, Miller S, Joseph RD, Geballe TR, Trafton LM, Tennyson J, Ballester GE. 1997 Variation in the H₃⁺ Emission of Uranus. *Astrophys. J.* **474**, L73–L76. (doi:10.1086/310424)

41. Bagnasco G et al. 2007 Overview of the near-infrared spectrograph (NIRSpec) instrument on-board the James Webb Space Telescope (JWST). In *Cryogenic Optical Systems and Instruments XII, Proc. of the SPIE*, vol. 6692, p. 66920M. San Diego, CA: Optical Engineering + Applications.