-
Simultaneous retrieval of orbital phase resolved JWST/MIRI emission spectra of the hot Jupiter WASP-43b: evidence of water, ammonia and carbon monoxide
Authors:
Jingxuan Yang,
Mark Hammond,
Anjali A. A. Piette,
Jasmina Blecic,
Taylor J. Bell,
Patrick G. J. Irwin,
Vivien Parmentier,
Shang-Min Tsai,
Joanna K. Barstow,
Nicolas Crouzet,
Laura Kreidberg,
João M. Mendonça,
Jake Taylor,
Robin Baeyens,
Kazumasa Ohno,
Lucas Teinturier,
Matthew C. Nixon
Abstract:
Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed…
▽ More
Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed with the MIRI/LRS onboard the JWST, to retrieve its atmospheric properties. Using a parametric 2D temperature model and assuming a chemically homogeneous atmosphere within the observed pressure region, we simultaneously fit the four spectra to constrain the abundances of atmospheric constituents, thereby yielding more precise constraints than previous work that analysed each spectrum independently. Our analysis reveals statistically significant evidence of NH3 (4$σ$) in a hot Jupiter's emission spectra for the first time, along with evidence of H2O (6.5$σ$), CO (3.1$σ$), and a non-detection of CH4. With our abundance constraints, we tentatively estimate the metallicity of WASP-43b at 0.6-6.5$\times$solar and its C/O ratio at 0.6-0.9. Our findings offer vital insights into the atmospheric conditions and formation history of WASP-43b by simultaneously constraining the abundances of carbon, oxygen, and nitrogen-bearing species.
△ Less
Submitted 5 June, 2024;
originally announced June 2024.
-
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.
-
Hot Jupiter Diversity and the Onset of TiO/VO Revealed by a Large Grid of Non-Grey Global Circulation Models
Authors:
Alexander Roth,
Vivien Parmentier,
Mark Hammond
Abstract:
The population of hot Jupiters is extremely diverse, with large variations in their irradiation, period, gravity and chemical composition. To understand the intrinsic planet diversity through the observed population level trends, we explore the a-priori scatter in the population created by the different responses of atmospheric circulation to planetary parameters. We use the SPARC/MITgcm 3D global…
▽ More
The population of hot Jupiters is extremely diverse, with large variations in their irradiation, period, gravity and chemical composition. To understand the intrinsic planet diversity through the observed population level trends, we explore the a-priori scatter in the population created by the different responses of atmospheric circulation to planetary parameters. We use the SPARC/MITgcm 3D global circulation model to simulate 345 planets spanning a wide range of instellation, metallicity, gravity and rotation periods typical for hot Jupiters, while differentiating between models with and without TiO/VO in their atmosphere. We show that the combined effect of the planetary parameters leads to a large diversity in the ability of atmospheres to transport heat from day-side to night-side at a given equilibrium temperature. We further show that the hot-spot offset is a non-monotonic function of planetary rotation period and explain our findings by a competition between the rotational and divergent parts of the circulation. As a consequence, hot-spot offset and phase curve amplitude are not necessarily correlated. Finally, we compare the observables from our grid to the population of Spitzer and Hubble observations of hot Jupiters. We find that the sudden jump in brightness temperature observed in the Spitzer secondary eclipse measurements can be naturally explained by the cold-trapping of TiO/VO at approximately 1800K. The grid of modelled spectra, phase curves and thermal structures are made available to the community, together with a python code for visualization of the grid properties, at https://doi.org/10.5281/zenodo.10785321 and http://sim3d.oca.eu/.
△ Less
Submitted 15 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.
-
Global Chemical Transport on Hot Jupiters: Insights from 2D VULCAN photochemical model
Authors:
Shang-Min Tsai,
Vivien Parmentier,
João M. Mendonça,
Xianyu Tan,
Russell Deitrick,
Mark Hammond,
Arjun B. Savel,
Xi Zhang,
Raymond T. Pierrehumbert,
Edward W. Schwieterman
Abstract:
The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photoc…
▽ More
The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below $\sim$0.1 mbar while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C$_2$H$_4$ absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres.
△ Less
Submitted 25 March, 2024; v1 submitted 26 October, 2023;
originally announced October 2023.
-
An analytic theory for the resolution attainable using eclipse mapping of exoplanets
Authors:
Sasha Boone,
David Grant,
Mark Hammond
Abstract:
We present an analytic theory for the resolution attainable via eclipse mapping of exoplanets, based on the Fourier components of the brightness distribution on the planetary disk. We find that the impact parameter determines which features can and cannot be seen, via the angle of the stellar edge relative to the axis of the orbit during the eclipse. We estimate the signal-to-noise ratio as a func…
▽ More
We present an analytic theory for the resolution attainable via eclipse mapping of exoplanets, based on the Fourier components of the brightness distribution on the planetary disk. We find that the impact parameter determines which features can and cannot be seen, via the angle of the stellar edge relative to the axis of the orbit during the eclipse. We estimate the signal-to-noise ratio as a function of mapping resolution, and use this to determine the attainable resolution for a given star-planet system. We test this theory against numerical simulations and find good agreement; in particular, our predictions for the resolution as a function of stellar edge angle are accurate to the simulated data to within 10% over a wide range of angles. Our prediction for the number of spatial modes that can be constrained given a light curve error is similarly accurate. Finally, we give a list of exoplanets with the best expected resolution for observations with the NIRISS SOSS, NIRSpec G395H, and MIRI LRS instruments on JWST.
△ Less
Submitted 22 October, 2023;
originally announced October 2023.
-
Interior dynamics of super-Earth 55 Cancri e
Authors:
Tobias G. Meier,
Dan J. Bower,
Tim Lichtenberg,
Mark Hammond,
Paul J. Tackley
Abstract:
The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot…
▽ More
The ultra-short-period super-Earth 55 Cancri e has a measured radius of 1.8 Earth radii. Previous thermal phase curve observations suggest a strong temperature contrast between the dayside and nightside of around 1000 K with the hottest point shifted $41\pm12$ degrees east from the substellar point, indicating some degree of heat circulation. The dayside (and potentially even the nightside) is hot enough to harbour a magma ocean. We use results from general circulation models (GCMs) of atmospheres to constrain the surface temperature contrasts. There is still a large uncertainty on the vigour and style of mantle convection in super-Earths, especially those that experience stellar irradiation large enough to harbour a magma ocean. In this work, we aim to constrain the mantle dynamics of the tidally locked lava world 55 Cancri e. Using the surface temperature contrasts as boundary condition, we model the mantle flow of 55 Cancri e using 2D mantle convection simulations and investigate how the convection regimes are affected by the different climate models. We find that large super-plumes form on the dayside if that hemisphere is covered by a magma ocean and the nightside remains solid or only partially molten. Cold material descends into the deep interior on the nightside, but no strong downwellings form. In some cases, the super-plume also moves several tens of degrees towards the terminator. A convective regime where the upwelling is preferentially on the dayside might lead to preferential outgassing on that hemisphere which could lead to the build-up of atmospheric species that could be chemically distinct from the nightside.
△ Less
Submitted 1 August, 2023;
originally announced August 2023.
-
Examining NHD vs QHD in the GCM THOR with non-grey radiative transfer for the hot Jupiter regime
Authors:
Pascal A. Noti,
Elspeth K. H. Lee,
Russell Deitrick,
Mark Hammond
Abstract:
Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations…
▽ More
Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations using the quasi-primitive dynamical equations (QHD) and the deep Navier-Stokes equations (NHD) in the GCM THOR. We utilise a two-stream non-grey "picket-fence" scheme to increase the realism of the radiative transfer calculations. We perform GCM simulations covering a wide parameter range grid of system parameters in the population of exoplanets. Our results show significant differences between simulations with the NHD and QHD equation sets at lower gravity, higher rotation rates or at higher irradiation temperatures. The chosen parameter range shows the relevance of choosing dynamical equation sets dependent on system and planetary properties. Our results show the climate states of hot Jupiters seem to be very diverse, where exceptions to prograde superrotation can often occur. Overall, our study shows the evolution of different climate states which arise just due to different selections of Navier-Stokes equations and approximations. We show the divergent behaviour of approximations used in GCMs for Earth, but applied for non Earth-like planets.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
A Mini-Chemical Scheme with Net Reactions for 3D GCMs II. 3D thermochemical modelling of WASP-39b and HD 189733b
Authors:
Elspeth K. H. Lee,
Shang-Min Tsai,
Mark Hammond,
Xianyu Tan
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-equili…
▽ More
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.
△ Less
Submitted 19 February, 2023;
originally announced February 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.
-
CAMEMBERT: A Mini-Neptunes GCM Intercomparison, Protocol Version 1.0. A CUISINES Model Intercomparison Project
Authors:
Duncan A. Christie,
Elspeth K. H. Lee,
Hamish Innes,
Pascal A. Noti,
Benjamin Charnay,
Thomas J. Fauchez,
Nathan J. Mayne,
Russell Deitrick,
Feng Ding,
Jennifer J. Greco,
Mark Hammond,
Isaac Malsky,
Avi Mandell,
Emily Rauscher,
Michael T. Roman,
Denis E. Sergeev,
Linda Sohl,
Maria E. Steinrueck,
Martin Turbet,
Eric T. Wolf,
Maria Zamyatina,
Ludmila Carone
Abstract:
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the…
▽ More
With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors.
△ Less
Submitted 8 November, 2022;
originally announced November 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.
-
K2 and Spitzer phase curves of the rocky ultra-short-period planet K2-141 b hint at a tenuous rock vapor atmosphere
Authors:
S. Zieba,
M. Zilinskas,
L. Kreidberg,
T. G. Nguyen,
Y. Miguel,
N. B. Cowan,
R. Pierrehumbert,
L. Carone,
L. Dang,
M. Hammond,
T. Louden,
R. Lupu,
L. Malavolta,
K. B. Stevenson
Abstract:
K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of t…
▽ More
K2-141 b is a transiting, small (1.5 Re) ultra-short-period (USP) planet discovered by Kepler orbiting a K-dwarf host star every 6.7 hours. The planet's high surface temperature makes it an excellent target for thermal emission observations. Here we present 65 hours of continuous photometric observations of K2-141 b collected with Spitzer's IRAC Channel 2 at 4.5 micron spanning 10 full orbits of the planet. We measure an infrared eclipse depth of 143 +/- 39 ppm and a peak to trough amplitude variation of 121 +/- 43 ppm. The best fit model to the Spitzer data shows no significant thermal hotspot offset, in contrast to the previously observed offset for the well-studied USP planet 55 Cnc e. We also jointly analyze the new Spitzer observations with the photometry collected by Kepler during two separate K2 campaigns. We model the planetary emission with a range of toy models that include a reflective and a thermal contribution. With a two-temperature model, we measure a dayside temperature of 2049 +/- 361 K and a night-side temperature that is consistent with zero (Tp,n < 1712 K at 2 sigma). Models with a steep dayside temperature gradient provide a better fit to the data than a uniform dayside temperature (DeltaBIC = 22.2). We also find evidence for a non-zero geometric albedo Ag = 0.28 +/- 0.07. We also compare the data to a physically motivated, pseudo-2D rock vapor model and a 1D turbulent boundary layer model. Both models fit the data well. Notably, we find that the optical eclipse depth can be explained by thermal emission from a hot inversion layer, rather than reflected light. A thermal inversion may also be responsible for the deep optical eclipse observed for another USP, Kepler-10 b. Finally, we significantly improve the ephemerides for K2-141 b and c, which will facilitate further follow-up observations of this interesting system with state-of-the-art observatories like JWST.
△ Less
Submitted 1 March, 2022;
originally announced March 2022.
-
General Circulation Model Errors are Variable across Exoclimate Parameter Spaces
Authors:
Pushkar Kopparla,
Russell Deitrick,
Kevin Heng,
João M. Mendonça,
Mark Hammond
Abstract:
General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales a…
▽ More
General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales and so on. In such studies, there is an implicit assumption that the model which works reasonably well for the standard test case will be credible at all points in an arbitrarily wide parameter space. Here, we test this assumption using the open-source general circulation model THOR to simulate atmospheric circulation on tidally locked Earth-like planets with rotation periods of 0.1 to 100 days. We find that the model error, as quantified by the ratio between physical and spurious numerical contributions to the angular momentum balance, is extremely variable across this range of rotation periods with some cases where numerical errors are the dominant component. Increasing model grid resolution does improve errors but using a higher-order numerical diffusion scheme can sometimes magnify errors for finite-volume dynamical solvers. We further show that to minimize error and make the angular momentum balance more physical within our model, the surface friction timescale must be smaller than the rotational timescale.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
Simulating gas giant exoplanet atmospheres with Exo-FMS: Comparing semi-grey, picket fence and correlated-k radiative-transfer schemes
Authors:
Elspeth K. H. Lee,
Vivien Parmentier,
Mark Hammond,
Simon L. Grimm,
Daniel Kitzmann,
Xianyu Tan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM us…
▽ More
Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM using these RT schemes to hot Jupiter simulation results from the literature. We perform a HD 209458b-like simulation with the three schemes and compare their results. These simulations are then post-processed to compare their observable differences. The semi-grey scheme results show qualitative agreement with previous studies in line with variations seen between GCM models. The real gas model reproduces well the temperature and dynamical structures from other studies. After post-processing our non-grey picket fence scheme compares very favourably with the real gas model, producing similar transmission spectra, emission spectra and phase curve behaviours. Exo-FMS is able to reliably reproduce the essential features of contemporary GCM models in the hot gas giant regime. Our results suggest the picket fence approach offers a simple way to improve upon RT realism beyond semi-grey schemes.
△ Less
Submitted 28 June, 2021; v1 submitted 22 June, 2021;
originally announced June 2021.
-
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.
-
Vertically resolved magma ocean-protoatmosphere evolution: H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, and N$_2$ as primary absorbers
Authors:
Tim Lichtenberg,
Dan J. Bower,
Mark Hammond,
Ryan Boukrouche,
Patrick Sanan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositi…
▽ More
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically-resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, or N$_2$. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N$_2$, and O$_2$ with minimal effect, H$_2$O, CO$_2$, and CH$_4$ with intermediate influence, and H$_2$ with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multi-wavelength astronomical observations.
△ Less
Submitted 26 January, 2021;
originally announced January 2021.
-
The Equatorial Jet Speed on Tidally Locked Planets: I -- Terrestrial Planets
Authors:
Mark Hammond,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and…
▽ More
The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and the vertical motion of the atmosphere. The three-dimensional structure of the zonal acceleration is vital to this mechanism.
We demonstrate this mechanism in a hierarchy of models. We calculate the three-dimensional stationary waves induced by the forcing on these planets, and show the vertical structure of the zonal acceleration produced by these waves, which we use to suggest a mechanism for how the jet forms. GCM simulations are used to confirm the equilibrium state predicted by this mechanism, where the acceleration from these waves is balanced by an interaction between the zonal-mean vertical velocity and the jet. We derive a simple model of this using the "Weak Temperature Gradient" approximation, which gives an estimate of the jet speed on a terrestrial tidally locked planet.
We conclude that the proposed mechanism is a good description of the formation of an equatorial jet on a terrestrial tidally locked planet, and should be useful for interpreting observations and simulations of these planets. The mechanism requires assumptions such as a large equatorial Rossby radius and weak acceleration due to transient waves, and a different mechanism may produce the equatorial jets on gaseous tidally locked planets.
△ Less
Submitted 1 September, 2020;
originally announced September 2020.
-
Simplified 3D GCM modelling of the irradiated brown dwarf WD0137-349B
Authors:
Elspeth K. H. Lee,
Sarah L. Casewell,
Katy L. Chubb,
Mark Hammond,
Xianyu Tan,
Shang-Min Tsai,
Raymond T. Pierrehumbert
Abstract:
Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical…
▽ More
Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical properties of the Brown dwarf WD0137-349B. Methods: We use the three dimensional GCM model Exo-FMS, with a dual-band grey radiative-transfer scheme to model the atmosphere of WD0137-349B. The results of the GCM model are post-processed using the three dimensional Monte Carlo radiative-transfer model \textsc{cmcrt}. Results: Our results suggest inefficient day-night energy transport and a large day-night temperature contrast for WD0137-349B. Multiple flow patterns are present, shifting energy asymmetrically eastward or westward depending on their zonal direction and latitude. Regions of overturning are produced on the western terminator. We are able to reproduce the start of the system near-IR emission excess at $\gtrsim$ 1.95 $μ$m as observed by the GNIRS instrument. Our model over predicts the IR phase curve fluxes by factors of $\approx$1-3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using \textsc{vulcan} suggests a highly ionised region at high altitudes can form on the dayside of the Brown dwarf. Conclusions: We present a first attempt at simulating the atmosphere of a short period White dwarf - Brown dwarf binary in a 3D setting. Further studies into the radiative and photochemical heating from the UV irradiation is required to more accurately capture the energy balance inside the Brown dwarf atmosphere. Cloud formation may also play an important role in shaping the emission spectra of the Brown dwarf.
△ Less
Submitted 25 June, 2020; v1 submitted 17 January, 2020;
originally announced January 2020.
-
Wave-Mean Flow Interactions in the Atmospheric Circulation of Tidally Locked Planets
Authors:
Mark Hammond,
Raymond T. Pierrehumbert
Abstract:
We use a linear shallow-water model to investigate the global circulation of the atmospheres of tidally locked planets. Simulations, observations, and simple models show that if these planets are sufficiently rapidly rotating, their atmospheres have an eastward equatorial jet and a hot-spot east of the substellar point. We linearize the shallow-water model about this eastward flow and its associat…
▽ More
We use a linear shallow-water model to investigate the global circulation of the atmospheres of tidally locked planets. Simulations, observations, and simple models show that if these planets are sufficiently rapidly rotating, their atmospheres have an eastward equatorial jet and a hot-spot east of the substellar point. We linearize the shallow-water model about this eastward flow and its associated height perturbation. The forced solutions of this system show that the shear flow explains the form of the global circulation, particularly the hot-spot shift and the positions of the cold standing waves on the night-side.
We suggest that the eastward hot-spot shift seen in observations and 3D simulations of these atmospheres is caused by the zonal flow Doppler-shifting the stationary wave response eastwards, summed with the height perturbation from the flow itself. This differs from other studies which explained the hot-spot shift as pure advection of heat from air flowing eastward from the substellar point, or as equatorial waves travelling eastwards. We compare our solutions to simulations in our climate model Exo-FMS, and show that the height fields and wind patterns match.
We discuss how planetary properties affect the global circulation, and how they change observables such as the hot-spot shift or day-night contrast. We conclude that the wave-mean flow interaction between the stationary planetary waves and the equatorial jet is a vital part of the equilibrium circulation on tidally locked planets.
△ Less
Submitted 26 October, 2018;
originally announced October 2018.
-
Linking the Climate and Thermal Phase Curve of 55 Cancri e
Authors:
Mark Hammond,
Raymond Pierrehumbert
Abstract:
The thermal phase curve of 55 Cancri e is the first measurement of the temperature distribution of a tidally locked Super-Earth, but raises a number of puzzling questions about the planet's climate. The phase curve has a high amplitude and peak offset, suggesting that it has a significant eastward hot-spot shift as well as a large day-night temperature contrast. We use a general circulation model…
▽ More
The thermal phase curve of 55 Cancri e is the first measurement of the temperature distribution of a tidally locked Super-Earth, but raises a number of puzzling questions about the planet's climate. The phase curve has a high amplitude and peak offset, suggesting that it has a significant eastward hot-spot shift as well as a large day-night temperature contrast. We use a general circulation model to model potential climates, and investigate the relation between bulk atmospheric composition and the magnitude of these seemingly contradictory features. We confirm theoretical models of tidally locked circulation are consistent with our numerical model of 55 Cnc e, and rule out certain atmospheric compositions based on their thermodynamic properties. Our best-fitting atmosphere has a significant hot-spot shift and day-night contrast, although these are not as large as the observed phase curve. We discuss possible physical processes which could explain the observations, and show that night-side cloud formation from species such as SiO from a day-side magma ocean could potentially increase the phase curve amplitude and explain the observations. We conclude that the observations could be explained by an optically thick atmosphere with a low mean molecular weight, a surface pressure of several bar and a strong eastward circulation, with night-side cloud formation a possible explanation for the difference between our model and the observations.
△ Less
Submitted 10 October, 2017;
originally announced October 2017.
-
Hydraulic External Pre-Isolator System for LIGO
Authors:
S. Wen,
R. Mittleman,
K. Mason,
J. Giaime,
R. Abbott,
J. Kern,
B. O'Reilly,
R. Bork,
M. Hammond,
C. Hardham,
B. Lantz,
W. Hua,
D. Coyne,
G. Traylor,
H. Overmier,
T. Evans,
J. Hanson,
O. Spjeld,
M. Macinnis,
K. Mailand,
D. Sellers,
K. Carter,
P. Sarin
Abstract:
The Hydraulic External Pre-Isolator (HEPI) is the first 6 degrees of freedom active seismic isolation system implemented at the Laser Interferometer Gravitational Wave Observatory (LIGO). Implementation was first completed at the LIGO Livingston Observatory (LLO) prior to LIGO's 5th science run, successfully cutting down the disturbance seen by LLO's suspended optics in the two most prominent seis…
▽ More
The Hydraulic External Pre-Isolator (HEPI) is the first 6 degrees of freedom active seismic isolation system implemented at the Laser Interferometer Gravitational Wave Observatory (LIGO). Implementation was first completed at the LIGO Livingston Observatory (LLO) prior to LIGO's 5th science run, successfully cutting down the disturbance seen by LLO's suspended optics in the two most prominent seismic disturbance bands, the microseism (0.1-0.3Hz) and the anthropogenic (1-3Hz) bands, by a factor of a few to tens. The improvement in seismic isolation contributed directly to LLO's much improved duty cycle of 66.7% and LIGO's triple coincident duty cycle of 53%. We report the design, control scheme, and isolation performance of HEPI at LLO in this paper. Aided with this success, funding for incorporating HEPI into the LIGO Hanford Observatory was approved and installation is currently underway.
△ Less
Submitted 17 October, 2013; v1 submitted 22 September, 2013;
originally announced September 2013.
-
A New Catalog of Faraday Rotation Measures and Redshifts for Extragalactic Radio Sources
Authors:
Alison M. Hammond,
Timothy Robishaw,
B. M. Gaensler
Abstract:
We present a catalog of Faraday rotation measures (RMs) and redshifts for 4003 extragalactic radio sources detected at 1.4 GHz, derived by identifying optical counterparts and spectroscopic redshifts for linearly polarized radio sources from the NRAO VLA Sky Survey. This catalog is more than an order of magnitude larger than any previous sample of RM vs. redshift, and covers the redshift range 0 <…
▽ More
We present a catalog of Faraday rotation measures (RMs) and redshifts for 4003 extragalactic radio sources detected at 1.4 GHz, derived by identifying optical counterparts and spectroscopic redshifts for linearly polarized radio sources from the NRAO VLA Sky Survey. This catalog is more than an order of magnitude larger than any previous sample of RM vs. redshift, and covers the redshift range 0 < z < 5.3 ; the median redshift of the catalog is z = 0.70, and there are more than 1500 sources at redshifts z > 1. For 3650 of these sources at Galactic latitudes |b| >= 20 degrees, we present a second catalog in which we have corrected for the foreground Faraday rotation of the Milky Way, resulting in an estimate of the residual rotation measure (RRM) that aims to isolate the contribution from extragalactic magnetic fields. We find no significant evolution of RRM with redshift, but observe a strong anti-correlation between RRM and fractional polarization, p, that we argue is the result of beam depolarization from small-scale fluctuations in the foreground magnetic field or electron density. We suggest that the observed variance in RRM and the anti-correlation of RRM with p both require a population of magnetized intervening objects that lie outside the Milky Way but in the foreground to the emitting sources.
△ Less
Submitted 8 May, 2013; v1 submitted 6 September, 2012;
originally announced September 2012.
