-
Identifying and Fitting Eclipse Maps of Exoplanets with Cross-Validation
Authors:
Mark Hammond,
Neil T. Lewis,
Sasha Boone,
Xueqing Chen,
João M. Mendonça,
Vivien Parmentier,
Jake Taylor,
Taylor Bell,
Leonardo dos Santos,
Nicolas Crouzet,
Laura Kreidberg,
Michael Radica,
Michael Zhang
Abstract:
Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the lati…
▽ More
Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the latitudinal structure is wrongly determined by the priors on the map. We present a new method to address this issue. The method tests for the presence of an eclipse mapping signal by using k-fold cross-validation to compare the performance of a simple mapping model to the null hypothesis of a uniform disk. If a signal is found, the method fits a map with more degrees of freedom, optimising its information content. The information content is varied by penalising the model likelihood by a factor proportional to the spatial entropy of the map, optimised by cross-validation. We demonstrate this method for simulated datasets then apply it to three observational datasets. The method identifies an eclipse mapping signal for JWST MIRI/LRS observations of WASP-43b but does not identify a signal for JWST NIRISS/SOSS observations of WASP-18b or Spitzer Space Telescope observations of HD 189733b. It is possible to fit eclipse maps to these datasets, but we suggest that these maps are overfitting the eclipse shape. We fit a new map with more spatial freedom to the WASP-43b dataset and show a flatter east-west structure than previously derived.
△ Less
Submitted 18 July, 2024; v1 submitted 31 May, 2024;
originally announced May 2024.
-
Two-Dimensional Eclipse Mapping of the Hot Jupiter WASP-43b with JWST MIRI/LRS
Authors:
Mark Hammond,
Taylor J. Bell,
Ryan C. Challener,
Neil T. Lewis,
Megan Weiner Mansfield,
Isaac Malsky,
Emily Rauscher,
Jacob L. Bean,
Ludmila Carone,
João M. Mendonça,
Lucas Teinturier,
Xianyu Tan,
Nicolas Crouzet,
Laura Kreidberg,
Giuseppe Morello,
Vivien Parmentier,
Jasmina Blecic,
Jean-Michel Désert,
Christiane Helling,
Pierre-Olivier Lagage,
Karan Molaverdikhani,
Matthew C. Nixon,
Benjamin V. Rackham,
Jingxuan Yang
Abstract:
We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order…
▽ More
We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order $\ell_{\rm max}=2$, alongside the planetary, orbital, stellar, and systematic parameters. This yields a map with a meridionally-averaged eastward hot-spot shift of $(7.75 \pm 0.36)^{\circ}$, with no significant degeneracy between the map and the additional parameters. We show the latitudinal and longitudinal contributions of the day-side emission structure to the eclipse shape, finding a latitudinal signal of $\sim$200 ppm and a longitudinal signal of $\sim$250 ppm. To investigate the sensitivity of the map to the method, we fix the non-mapping parameters and derive an "eigenmap" fitted with an optimised number of orthogonal phase curves, which yields a similar map to the $\ell_{\rm max}=2$ map. We also fit a map up to $\ell_{\rm max}=3$, which shows a smaller hot-spot shift, with a larger uncertainty. These maps are similar to those produced by atmospheric simulations. We conclude that there is a significant mapping signal which constrains the spherical harmonic components of our model up to $\ell_{\rm max}=2$. Alternative mapping models may derive different structures with smaller-scale features; we suggest that further observations of WASP-43b and other planets will drive the development of more robust methods and more accurate maps.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b
Authors:
Taylor J. Bell,
Nicolas Crouzet,
Patricio E. Cubillos,
Laura Kreidberg,
Anjali A. A. Piette,
Michael T. Roman,
Joanna K. Barstow,
Jasmina Blecic,
Ludmila Carone,
Louis-Philippe Coulombe,
Elsa Ducrot,
Mark Hammond,
João M. Mendonça,
Julianne I. Moses,
Vivien Parmentier,
Kevin B. Stevenson,
Lucas Teinturier,
Michael Zhang,
Natalie M. Batalha,
Jacob L. Bean,
Björn Benneke,
Benjamin Charnay,
Katy L. Chubb,
Brice-Olivier Demory,
Peter Gao
, et al. (58 additional authors not shown)
Abstract:
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5…
▽ More
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions).
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Equatorial Waves and Superrotation in the Stratosphere of a Titan General Circulation Model
Authors:
Neil T. Lewis,
Nicholas A. Lombardo,
Peter L. Read,
Juan M. Lora
Abstract:
We investigate the characteristics of equatorial waves associated with the maintenance of superrotation in the stratosphere of a Titan general circulation model. A variety of equatorial waves are present in the model atmosphere, including equatorial Kelvin waves, equatorial Rossby waves, and mixed Rossby-gravity waves. In the upper stratosphere, acceleration of superrotation is strongest around so…
▽ More
We investigate the characteristics of equatorial waves associated with the maintenance of superrotation in the stratosphere of a Titan general circulation model. A variety of equatorial waves are present in the model atmosphere, including equatorial Kelvin waves, equatorial Rossby waves, and mixed Rossby-gravity waves. In the upper stratosphere, acceleration of superrotation is strongest around solstice and is due to interaction between equatorial Kelvin waves and Rossby-type waves in winter-hemisphere mid-latitudes. The existence of this 'Rossby-Kelvin'-type wave appears to depend on strong meridional shear of the background zonal wind that occurs in the upper stratosphere at times away from the equinoxes. In the lower stratosphere, acceleration of superrotation occurs throughout the year and is partially induced by equatorial Rossby waves, which we speculate are generated by quasigeostrophic barotropic instability. Acceleration of superrotation is generally due to waves with phase speeds close to the zonal velocity of the mean flow. Consequently, they have short vertical wavelengths which are close to the model's vertical grid scale, and therefore are likely to be not properly represented. We suggest this may be a common issue amongst Titan GCMs which should be addressed by future model development.
△ Less
Submitted 13 July, 2023;
originally announced July 2023.
-
A broadband thermal emission spectrum of the ultra-hot Jupiter WASP-18b
Authors:
Louis-Philippe Coulombe,
Björn Benneke,
Ryan Challener,
Anjali A. A. Piette,
Lindsey S. Wiser,
Megan Mansfield,
Ryan J. MacDonald,
Hayley Beltz,
Adina D. Feinstein,
Michael Radica,
Arjun B. Savel,
Leonardo A. Dos Santos,
Jacob L. Bean,
Vivien Parmentier,
Ian Wong,
Emily Rauscher,
Thaddeus D. Komacek,
Eliza M. -R. Kempton,
Xianyu Tan,
Mark Hammond,
Neil T. Lewis,
Michael R. Line,
Elspeth K. H. Lee,
Hinna Shivkumar,
Ian J. M. Crossfield
, et al. (51 additional authors not shown)
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 conten…
▽ More
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 $μ$m in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at $>$6$σ$ confidence) and evidence for optical opacity, possibly due to H$^-$, TiO, and VO (combined significance of 3.8$σ$). 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.
△ Less
Submitted 20 January, 2023; v1 submitted 19 January, 2023;
originally announced January 2023.
-
Bistability of the atmospheric circulation on TRAPPIST-1e
Authors:
Denis E. Sergeev,
Neil T. Lewis,
F. Hugo Lambert,
Nathan J. Mayne,
Ian A. Boutle,
James Manners,
Krisztian Kohary
Abstract:
Using a 3D general circulation model, we demonstrate that a confirmed rocky exoplanet and a primary observational target, TRAPPIST-1e presents an interesting case of climate bistability. We find that the atmospheric circulation on TRAPPIST-1e can exist in two distinct regimes for a 1~bar nitrogen-dominated atmosphere. One is characterized by a single strong equatorial prograde jet and a large day-…
▽ More
Using a 3D general circulation model, we demonstrate that a confirmed rocky exoplanet and a primary observational target, TRAPPIST-1e presents an interesting case of climate bistability. We find that the atmospheric circulation on TRAPPIST-1e can exist in two distinct regimes for a 1~bar nitrogen-dominated atmosphere. One is characterized by a single strong equatorial prograde jet and a large day-night temperature difference; the other is characterized by a pair of mid-latitude prograde jets and a relatively small day-night contrast. The circulation regime appears to be highly sensitive to the model setup, including initial and surface boundary conditions, as well as physical parameterizations of convection and cloud radiative effects. We focus on the emergence of the atmospheric circulation during the early stages of simulations and show that the regime bistability is associated with a delicate balance between the zonally asymmetric heating, mean overturning circulation, and mid-latitude baroclinic instability. The relative strength of these processes places the GCM simulations on different branches of the evolution of atmospheric dynamics. The resulting steady states of the two regimes have consistent differences in the amount of water content and clouds, affecting the water absorption bands as well as the continuum level in the transmission spectrum, although they are too small to be detected with current technology. Nevertheless, this regime bistability affects the surface temperature, especially on the night side of the planet, and presents an interesting case for understanding atmospheric dynamics and highlights uncertainty in 3D GCM results, motivating more multi-model studies.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
Temperature Structures Associated with Different Components of the Atmospheric Circulation on Tidally Locked Exoplanets
Authors:
Neil T. Lewis,
Mark Hammond
Abstract:
Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is…
▽ More
Observations of time-resolved thermal emission from tidally locked exoplanets can tell us about their atmospheric temperature structure. Telescopes such as JWST and ARIEL will improve the quality and availability of these measurements. This motivates an improved understanding of the processes that determine atmospheric temperature structure, particularly atmospheric circulation. The circulation is important in determining atmospheric temperatures, not only through its ability to transport heat, but also because any circulation pattern needs to be balanced by horizontal pressure contrasts, therefore implying a particular temperature structure. In this work, we show how the global temperature field on a tidally locked planet can be decomposed into contributions that are balanced by different components of the atmospheric circulation. These are the superrotating jet, stationary Rossby waves, and the divergent circulation. To achieve this, we partition the geopotential field into components balanced by the divergent circulation and the rotational circulation, with the latter comprising the jet and Rossby waves. The partitioned geopotential then implies a corresponding partitioning of the temperature via the hydrostatic relation. We apply these diagnostics to idealised general circulation model simulations, to show how the separate rotational and divergent circulations together make up the total three-dimensional atmospheric temperature structure. We also show how each component contributes distinct signatures to the thermal phase curve of a tidally locked planet. We conclude that this decomposition is a physically meaningful separation of the temperature field that explains its global structure, and can be used to fit observations of thermal emission.
△ Less
Submitted 13 April, 2022;
originally announced April 2022.
-
The rotational and divergent components of atmospheric circulation on tidally locked planets
Authors:
Mark Hammond,
Neil T. Lewis
Abstract:
Tidally locked exoplanets likely host global atmospheric circulations with a superrotating equatorial jet, planetary-scale stationary waves and thermally-driven overturning circulation. In this work, we show that each of these features can be separated from the total circulation by using a Helmholtz decomposition, which splits the circulation into rotational (divergence free) and divergent (vortic…
▽ More
Tidally locked exoplanets likely host global atmospheric circulations with a superrotating equatorial jet, planetary-scale stationary waves and thermally-driven overturning circulation. In this work, we show that each of these features can be separated from the total circulation by using a Helmholtz decomposition, which splits the circulation into rotational (divergence free) and divergent (vorticity free) components. This technique is applied to the simulated circulation of a terrestrial planet and a gaseous hot Jupiter. For both planets, the rotational component comprises the equatorial jet and stationary waves, and the divergent component contains the overturning circulation. Separating out each component allows us to evaluate their spatial structure and relative contribution to the total flow. In contrast with previous work, we show that divergent velocities are not negligible when compared with rotational velocities, and that divergent, overturning circulation takes the form of a single, roughly isotropic cell that ascends on the day-side and descends on the night-side. These conclusions are drawn for both the terrestrial case and the hot Jupiter. To illustrate the utility of the Helmholtz decomposition for studying atmospheric processes, we compute the contribution of each of the circulation components to heat transport from day- to night-side. Surprisingly, we find that the divergent circulation dominates day-night heat transport in the terrestrial case and accounts for around half of the heat transport for the hot Jupiter. The relative contributions of the rotational and divergent components to day-night heat transport are likely sensitive to multiple planetary parameters and atmospheric processes, and merit further study.
△ Less
Submitted 23 February, 2021;
originally announced February 2021.
-
The dependence of global super-rotation on planetary rotation rate
Authors:
Neil T. Lewis,
Greg J. Colyer,
Peter L. Read
Abstract:
An atmosphere may be described as globally super-rotating if its total zonal angular momentum exceeds that associated with solid-body co-rotation with the underlying planet. In this paper, we discuss the dependence of global super-rotation in terrestrial atmospheres on planetary rotation rate. This dependence is revealed through analysis of global super-rotation in idealised General Circulation Mo…
▽ More
An atmosphere may be described as globally super-rotating if its total zonal angular momentum exceeds that associated with solid-body co-rotation with the underlying planet. In this paper, we discuss the dependence of global super-rotation in terrestrial atmospheres on planetary rotation rate. This dependence is revealed through analysis of global super-rotation in idealised General Circulation Model experiments with time-independent axisymmetric forcing, compared with estimates for global super-rotation in Solar System atmospheres. Axisymmetric and three-dimensional experiments are conducted. We find that the degree of global super-rotation in the three-dimensional experiments is closely related to that of the axisymmetric experiments, with some differences in detail. A scaling theory for global super-rotation in an axisymmetric atmosphere is derived from the Held-Hou model. At high rotation rate, our numerical experiments inhabit a regime where global super-rotation scales geostrophically, and we suggest that the Earth and Mars occupy this regime. At low rotation rate, our experiments occupy a regime determined by angular momentum conservation, where global super-rotation is independent of rotation rate. Global super-rotation in our experiments saturates at a value significantly lower than that achieved in the atmospheres of Venus and Titan, which instead occupy a regime where global super-rotation scales cyclostrophically. This regime can only be accessed when eddy induced up-gradient angular momentum transport is sufficiently large, which is not the case in our idealised numerical experiments. We suggest that the 'default' regime for a slowly rotating planet is the angular momentum conserving regime, characterised by mild global (and local) superrotation.
△ Less
Submitted 17 April, 2020;
originally announced April 2020.
-
Investigating the Semiannual Oscillation on Mars using data assimilation
Authors:
Tao Ruan,
Neil T. Lewis,
Stephen R. Lewis,
Luca Montabone,
Peter L. Read
Abstract:
A Martian semiannual oscillation (SAO), similar to that in the Earth's tropical stratosphere, is evident in the Mars Analysis Correction Data Assimilation reanalysis dataset (MACDA) version 1.0, not only in the tropics, but also extending to higher latitudes. Unlike on Earth, the Martian SAO is found not always to reverse its zonal wind direction, but only manifests itself as a deceleration of the…
▽ More
A Martian semiannual oscillation (SAO), similar to that in the Earth's tropical stratosphere, is evident in the Mars Analysis Correction Data Assimilation reanalysis dataset (MACDA) version 1.0, not only in the tropics, but also extending to higher latitudes. Unlike on Earth, the Martian SAO is found not always to reverse its zonal wind direction, but only manifests itself as a deceleration of the dominant wind at certain pressure levels and latitudes. Singular System Analysis (SSA) is further applied on the zonal-mean zonal wind in different latitude bands to reveal the characteristics of SAO phenomena at different latitudes. The second pair of principal components (PCs) is usually dominated by a SAO signal, though the SAO signal can be strong enough to manifest itself also in the first pair of PCs. An analysis of terms in the Transformed Eulerian Mean equation (TEM) is applied in the tropics to further elucidate the forcing processes driving the tendency of the zonal-mean zonal wind. The zonal-mean meridional advection is found to correlate strongly with the observed oscillations of zonal-mean zonal wind, and supplies the majority of the westward (retrograde) forcing in the SAO cycle. The forcing due to various non-zonal waves supplies forcing to the zonal-mean zonal wind that is nearly the opposite of the forcing due to meridional advection above ~3 Pa altitude, but it also partly supports the SAO between 40 Pa and 3 Pa. Some distinctive features occurring during the period of the Mars year (MY) 25 global-scale dust storm (GDS) are also notable in our diagnostic results with substantially stronger values of eastward and westward momentum in the second half of MY 25 and stronger forcing due to vertical advection, transient waves and thermal tides.
△ Less
Submitted 17 June, 2019;
originally announced June 2019.
-
The influence of a sub-stellar continent on the climate of a tidally-locked exoplanet
Authors:
Neil T. Lewis,
F. Hugo Lambert,
Ian A. Boutle,
Nathan J. Mayne,
James Manners,
David M. Acreman
Abstract:
Previous studies have demonstrated that continental carbon-silicate weathering is important to the continued habitability of a terrestrial planet. Despite this, few studies have considered the influence of land on the climate of a tidally-locked planet. In this work we use the Met Office Unified Model, coupled to a land surface model, to investigate the climate effects of a continent located at th…
▽ More
Previous studies have demonstrated that continental carbon-silicate weathering is important to the continued habitability of a terrestrial planet. Despite this, few studies have considered the influence of land on the climate of a tidally-locked planet. In this work we use the Met Office Unified Model, coupled to a land surface model, to investigate the climate effects of a continent located at the sub-stellar point. We choose to use the orbital and planetary parameters of Proxima Centauri B as a template, to allow comparison with the work of others. A region of the surface where $T_{\text{s}} > 273.15\,\text{K}$ is always retained, and previous conclusions on the habitability of Proxima Centauri B remain intact. We find that sub-stellar land causes global cooling, and increases day-night temperature contrasts by limiting heat redistribution. Furthermore, we find that sub-stellar land is able to introduce a regime change in the atmospheric circulation. Specifically, when a continent offset to the east of the sub-stellar point is introduced, we observe the formation of two mid-latitude counterrotating jets, and a substantially weakened equatorial superrotating jet.
△ Less
Submitted 1 February, 2018;
originally announced February 2018.