Dr Xianyu Tan

Abstract:

The impressive sensitivity of the James Webb Space Telescope has made it possible to study the atmospheres of planets beyond the solar system. It will soon be followed by space missions aiming specifically at this goal, such as the Ariel mission, Twinkle, and the Habitable Worlds Observatory. One category of exoplanet has drawn interest because of its potential to harbour temperate climates with liquid surface water—and therefore potentially life. These are rocky planets orbiting cool M-class stars, or "M-Earths." Stellar population trends and observing biases lead to a high proportion of potentially habitable, terrestrial planets falling into this category. Because of the low temperatures of their host stars, however, habitable worlds of this type are found in close orbits where they are likely to be tidally locked. As the solar system has no tidally locked planets, our knowledge of their atmospheric circulation is currently limited to theoretical modelling. Past modelling work has shown that the asymmetrical irradiation of tidally locked planets results in characteristic circulation regimes which have profound consequences for observations. Atmospheric retrievals, which use statistical methods to fit 1-D atmospheric models to observational data and quantify the confidence of the fit, are not yet able to account for the 3-D nature of this circulation. For planets with large spatial variation in environmental conditions caused by tidal locking, 1-D models are not able to capture the differences and interconnections between planetary regions such as the dayside, nightside, and planetary limb or terminator. In addition, planetary atmospheres exhibit variation over time, potentially resulting in differences in retrieved properties between observing visits or even between different phases of a planet’s orbit. Accounting for 4-D circulation effects in atmospheric retrievals first requires a theoretical understanding of the impact of global-scale phenomena such as atmospheric waves and horizontal transport on conditions at the planetary limb, and then requires incorporation of this knowledge into the retrieval pipeline in the form of, for example, parameterisations. In this thesis, I address the first requirement: the theoretical understanding of the effects of fully modelled 4-D atmospheric circulation on the planetary limb, the region probed by transmission spectroscopy, on tidally locked planets. I focus in particular on effects caused by the global propagation of atmospheric waves and by horizontal transport of clouds and hazes. In Chapter 2, I show that that the atmospheric dynamics on the tidally locked Proxima Centauri b support a longitudinally asymmetric stratospheric wind oscillation (LASO), analogous to Earth’s quasi-biennial oscillation (QBO). The LASO has a vertical extent of 35–55 km, a period of 5–6.5 months, and a peak-to-peak wind speed amplitude of -70 to +130 ms−1 with a maximum at an altitude of 41 km. Unlike the QBO, the LASO displays longitudinal asymmetries related to the asymmetric thermal forcing of the planet and to interactions with the resulting stationary Rossby waves. The equatorial gravity wave sources driving the LASO are localised in the deep convection region at the substellar point and in a jet exit region near the western terminator, unlike the QBO, for which these sources are distributed uniformly around the planet. Longitudinally, the western terminator experiences the highest wind speeds and undergoes reversals earlier than other longitudes. The antistellar point only experiences a weak oscillation with a very brief, low-speed westward phase. The QBO on Earth is associated with fluctuations in the abundances of water vapour and trace gases such as ozone which are also likely to occur on exoplanets if these gases are present. Strong fluctuations in temperature and the abundance of atmospheric species at the terminators will need to be considered when interpreting atmospheric observations of tidally locked exoplanets. In Chapter 3, I investigate the presence of cloud cover at the planetary limb of water-rich Earth-like planets, which is likely to weaken chemical signatures in transmission spectra and impede attempts to characterise these atmospheres. Based on observations of Earth and solar system worlds, exoplanets with atmospheres should have both short-term weather and long-term climate variability, implying that cloud cover may be less during some observing periods. I identify and describe a mechanism driving periodic clear sky events at the terminators in simulations of tidally locked Earth-like planets. A feedback between dayside cloud radiative effects, incoming stellar radiation and heating, and the dynamical state of the atmosphere, especially the zonal wavenumber-1 Rossby wave identified in past work on tidally locked planets, leads to oscillations in Rossby wave phase speeds and in the position of Rossby gyres and results in advection of clouds to or away from the planet’s eastern terminator. I study this oscillation in simulations of Proxima Centauri b, TRAPPIST 1-e, and rapidly rotating versions of these worlds located at the inner edge of their stars’ habitable zones. I simulate time series of the transit depths of the 1.4 µm water feature and 2.7 µm carbon dioxide feature. The impact of atmospheric variability on the transmission spectra is sensitive to the structure of the dayside cloud cover and the location of the Rossby gyres, but none of my simulations have variability significant enough to be detectable with current methods. In Chapter 4, I study the interaction between the atmospheric circulation and photochemical hazes and describe the resulting haze abundances at the terminator. Transmission spectroscopy supports the presence of unknown, light-scattering aerosols in the atmospheres of many exoplanets. The complexity of factors influencing the formation, 3-D transport, radiative impact, and removal of aerosols makes it challenging to match theoretical models to the existing data. My study simplifies these factors to focus on the interaction between planetary general circulation and haze distribution at the planetary limb. I use an intermediate complexity general circulation model, ExoPlaSim, to simulate idealised organic haze particles as radiatively active tracers in the atmospheres of tidally locked terrestrial planets for a range of rotation rates. I find three distinct 3-D spatial haze distributions, corresponding to three circulation regimes, each with a different haze profile at the limb. All regimes display significant terminator asymmetry. In my parameter space, super-Earth-sized planets with rotation periods greater than 13 days have the lowest haze optical depths at the terminator, supporting the choice of slower rotators as observing targets. My thesis supports the existence of characteristic forms of temporal and spatial variability on tidally locked planets which will undoubtedly impact observations and inform our understanding of climate conditions on the surface. Overall, the effects of purely dynamical variability may be too small to be detected for Earth-like planets (but potentially detectable for larger ones). The impact of the atmospheric circulation on the distribution of clouds and hazes, on the other hand, is likely to affect even observations of terrestrial planets due to the highly scattering nature of these aerosols and will need to be accounted for in atmospheric retrievals

Abstract:

Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information content of the data resulted in high sensitivity to the varying assumptions made in the treatment of instrument systematics and the atmospheric retrieval analysis. Here we present a dayside thermal emission spectrum of the ultra-hot Jupiter WASP-18b obtained with the NIRISS instrument on JWST. The data span 0.85 to 2.85 $\mu$m in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at $>$6$\sigma$ confidence) and evidence for optical opacity, possibly due to H$^-$, TiO, and VO (combined significance of 3.8$\sigma$). Models that fit the data require a thermal inversion, molecular dissociation as predicted by chemical equilibrium, a solar heavy element abundance (''metallicity'', M/H = 1.03$_{-0.51}^{+1.11}$ $\times$ solar), and a carbon-to-oxygen (C/O) ratio less than unity. The data also yield a dayside brightness temperature map, which shows a peak in temperature near the sub-stellar point that decreases steeply and symmetrically with longitude toward the terminators.Comment: JWST ERS bright star observations. Uploaded to inform JWST Cycle 2 proposals. Manuscript under review. 50 pages, 14 figures, 2 table

Abstract:

With infrared flux contrasts larger than typically seen in hot Jupiter, tidally locked white dwarf–brown dwarf binaries offer a superior opportunity to investigate atmospheric processes in irradiated atmospheres. NLTT5306 is such a system, with a M _BD = 52 ± 3 M _Jup brown dwarf, orbiting a T _eff = 7756 ± 35 K white dwarf with an ultra-short period of ∼102 minutes. We present Hubble Space Telescope/Wide Field Camera 3 spectroscopic phase curves of NLTT5306, consisting of 47 spectra from 1.1 to 1.7 μ m with an average signal-to-noise ratio ∼ 65 per wavelength. We extracted the phase-resolved spectra of the brown dwarf NLTT5306B, finding a small <100 K day–night temperature difference (∼5% of the average day-side temperature). Our best-fit model phase curves revealed a complex wavelength-dependence on amplitudes and relative phase offsets, suggesting longitudinal–vertical atmospheric structure. The night-side spectrum was well fit by a cloudy, nonirradiated atmospheric model while the day side was best matched by a cloudy, weakly irradiated model. Additionally, we created a simple radiative energy redistribution model of the atmosphere and found evidence for efficient day-to-night heat redistribution and a moderately high Bond albedo. We also discovered an internal heat flux much higher than expected given the published system age, leading to an age reassessment that resulted in NLTT5306B most likely being much younger. We find that NLTT5306B is the only known significantly irradiated brown dwarf where the global temperature structure is not dominated by external irradiation, but rather its own internal heat. Our study provides an essential insight into the drivers of global circulation and day-to-night heat transport as a function of irradiation, rotation rate, and internal heat

Abstract:

The chemical inventory of hot Jupiter (HJ) exoplanets atmospheres continue to be observed by various ground and space based instruments in increasing detail and precision. It is expected that some HJs will exhibit strong non-equilibrium chemistry characteristics in their atmospheres, which might be inferred from spectral observations. We aim to model the three dimensional thermochemical non-equilibrium chemistry in the atmospheres of the HJs WASP-39b and HD 189733b. We couple a lightweight, reduced chemical network `mini-chem' that utilises net reaction rate tables to the Exo-FMS General Circulation Model (GCM). We perform GCM models of the exoplanets WASP-39b and HD 189733b as case studies of the coupled mini-chem scheme. The GCM results are then post-processed using the 3D radiative-transfer model gCMCRT to produce transmission and emission spectra to assess the impact of non-equilibrium chemistry on their observable properties. Both simulations show significant departures from chemical equilibrium (CE) due to the dynamical motions of the atmosphere. The spacial distribution of species generally follows closely the dynamical features of the atmosphere rather than the temperature field. Each molecular species exhibits a different quench level in the simulations, also dependent on the latitude of the planet. Major differences are seen in the transmission and emission spectral features between the CE and kinetic models. Our simulations indicate that considering the 3D kinetic chemical structures of HJ atmospheres has an important impact on physical interpretation of observational data. Drawing bulk atmospheric parameters from fitting feature strengths may lead to inaccurate interpretation of chemical conditions in the atmosphere of HJs. Our open source mini-chem module is simple to couple with contemporary HJ GCM models without substantially increasing required computational resources.Comment: A&A Accepted 17 Feb 202