-
A New Analysis of 8 Spitzer Phase Curves and Hot Jupiter Population Trends: Qatar-1b, Qatar-2b, WASP-52b, WASP-34b, and WASP-140b
Authors:
Erin May,
Kevin Stevenson,
Jacob Bean,
Taylor Bell,
Nicolas Cowan,
Lisa Dang,
Jean-Michel Desert,
Jonathan Fortney,
Dylan Keating,
Eliza Kempton,
Thaddeus Komacek,
Nikole Lewis,
Megan Mansfield,
Caroline Morley,
Vivien Parmentier,
Emily Rauscher,
Mark Swain,
Robert Zellem,
Adam Showman
Abstract:
With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as w…
▽ More
With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as well as the reanalysis of the 4.5 $\micron$ Qatar-1b phase curve due to the similar equilibrium temperature. In total, five 4.5 $\micron$ phase curves and three 3.6 $\micron$ phase curves are analyzed here with a uniform approach. Using these new results, in combination with literature values for the entire population of published Spitzer phase curves of hot Jupiters, we present evidence for a linear trend of increasing hot spot offset with increasing orbital period, as well as observational evidence for two classes of planets in apparent redistribution vs. equilibrium temperature parameter space, and tentative evidence for a dependence of hot spot offset on planetary surface gravity in our $\sim$ 1300 K sample. We do not find trends in apparent heat redistribution with orbital period or gravity. Non-uniformity in literature Spitzer data analysis techniques precludes a definitive determination of the sources or lack of trends.
△ Less
Submitted 28 March, 2022;
originally announced March 2022.
-
Mapping the Surface of Partially Cloudy Exoplanets is Hard
Authors:
Lucas Teinturier,
Nicholas Vieira,
Elisa Jacquet,
Juliette Geoffrion,
Youssef Bestavros,
Dylan Keating,
Nicolas B. Cowan
Abstract:
Reflected light photometry of terrestrial exoplanets could reveal the presence of oceans and continents, hence placing direct constraints on the current and long-term habitability of these worlds. Inferring the albedo map of a planet from its observed light curve is challenging because different maps may yield indistinguishable light curves. This degeneracy is aggravated by changing clouds. It has…
▽ More
Reflected light photometry of terrestrial exoplanets could reveal the presence of oceans and continents, hence placing direct constraints on the current and long-term habitability of these worlds. Inferring the albedo map of a planet from its observed light curve is challenging because different maps may yield indistinguishable light curves. This degeneracy is aggravated by changing clouds. It has previously been suggested that disk-integrated photometry spanning multiple days could be combined to obtain a cloud-free surface map of an exoplanet. We demonstrate this technique as part of a Bayesian retrieval by simultaneously fitting for the fixed surface map of a planet and the time-variable overlying clouds. We test this approach on synthetic data then apply it to real disk-integrated observations of the Earth. We find that eight days of continuous synthetic observations are sufficient to reconstruct a faithful low resolution surface albedo map, without needing to make assumptions about cloud physics. For lightcurves with negligible photometric uncertainties, the minimal top-of-atmosphere albedo at a location is a good estimate of its surface albedo. When applied to observations from the Earth Polychromating Imaging Camera aboard the DSCOVR spacecraft, our approach removes only a small fraction of clouds. We attribute this difficulty to the full-phase geometry of observations combined with the short correlation length for Earth clouds. For exoplanets with Earth-like climatology, it may be hard to do much better than a cloud-averaged map. We surmise that cloud removal will be most successful for exoplanets imaged near quarter phase that harbour large cloud systems.
△ Less
Submitted 3 January, 2022;
originally announced January 2022.
-
Longitudinally Resolved Spectral Retrieval (ReSpect) of WASP-43b
Authors:
Patricio E. Cubillos,
Dylan Keating,
Nicolas B. Cowan,
Johanna M. Vos,
Ben Burningham,
Marie Ygouf,
Theodora Karalidi,
Yifan Zhou,
Eileen C. Gonzales
Abstract:
Thermal phase variations of short period planets indicate that they are not spherical cows: day-to-night temperature contrasts range from hundreds to thousands of degrees, rivaling their vertical temperature contrasts. Nonetheless, the emergent spectra of short-period planets have typically been fit using one-dimensional (1D) spectral retrieval codes that only account for vertical temperature grad…
▽ More
Thermal phase variations of short period planets indicate that they are not spherical cows: day-to-night temperature contrasts range from hundreds to thousands of degrees, rivaling their vertical temperature contrasts. Nonetheless, the emergent spectra of short-period planets have typically been fit using one-dimensional (1D) spectral retrieval codes that only account for vertical temperature gradients. The popularity of 1D spectral retrieval codes is easy to understand: they are robust and have a rich legacy in Solar System atmospheric studies. Exoplanet researchers have recently introduced multi-dimensional retrieval schemes for interpreting the spectra of short-period planets, but these codes are necessarily more complex and computationally expensive than their 1D counterparts. In this paper we present an alternative: phase-dependent spectral observations are inverted to produce longitudinally resolved spectra that can then be fitted using standard 1D spectral retrieval codes. We test this scheme on the iconic phase-resolved spectra of WASP-43b and on simulated JWST observations using the open-source pyratbay 1D spectral retrieval framework. Notably, we take the model complexity of the simulations one step further over previous studies by allowing for longitudinal variations in composition in addition to temperature. We show that performing 1D spectral retrieval on longitudinally resolved spectra is more accurate than applying 1D spectral retrieval codes to disk-integrated emission spectra, despite being identical in terms of computational load. We find that for the extant Hubble and Spitzer observations of WASP-43b the difference between the two approaches is negligible but that JWST phase measurements should be treated with longitudinally \textbf{re}solved \textbf{spect}ral retrieval (ReSpect).
△ Less
Submitted 13 May, 2021;
originally announced May 2021.
-
Atmospheric characterization of hot Jupiters using hierarchical models of Spitzer observations
Authors:
Dylan Keating,
Nicolas B. Cowan
Abstract:
The field of exoplanet atmospheric characterization is trending towards comparative studies involving many planetary systems, and using Bayesian hierarchical modelling is a natural next step. Here we demonstrate two use cases. We first use hierarchical modelling to quantify variability in repeated observations by reanalyzing a suite of ten Spitzer secondary eclipse observations of the hot Jupiter…
▽ More
The field of exoplanet atmospheric characterization is trending towards comparative studies involving many planetary systems, and using Bayesian hierarchical modelling is a natural next step. Here we demonstrate two use cases. We first use hierarchical modelling to quantify variability in repeated observations by reanalyzing a suite of ten Spitzer secondary eclipse observations of the hot Jupiter XO-3b. We compare three models: one where we fit ten separate eclipse depths, one where we use a single eclipse depth for all ten observations, and a hierarchical model. By comparing the Widely Applicable Information Criterion of each model, we show that the hierarchical model is preferred over the others. The hierarchical model yields less scatter across the suite of eclipse depths -- and higher precision on the individual eclipse depths -- than does fitting the observations separately. We find that the hierarchical eclipse depth uncertainty is larger than the uncertainties on the individual eclipse depths, which suggests either slight astrophysical variability or that single eclipse observations underestimate the true eclipse depth uncertainty. Finally, we fit a suite of published dayside brightness measurements for 37 planets using a hierarchical model of brightness temperature vs irradiation temperature. The hierarchical model gives tighter constraints on the individual brightness temperatures than the non-hierarchical model. Although we tested hierarchical modelling on Spitzer eclipse data of hot Jupiters, it is applicable to observations of smaller planets like hot neptunes and super earths, as well as for photometric and spectroscopic transit or phase curve observations.
△ Less
Submitted 7 October, 2021; v1 submitted 26 February, 2021;
originally announced March 2021.
-
Eigenspectra: A Framework for Identifying Spectra from 3D Eclipse Mapping
Authors:
Megan Mansfield,
Everett Schlawin,
Jacob Lustig-Yaeger,
Arthur D. Adams,
Emily Rauscher,
Jacob Arcangeli,
Y. Katherina Feng,
Prashansa Gupta,
Dylan Keating,
Kevin B. Stevenson,
Thomas G. Beatty
Abstract:
Planetary atmospheres are inherently 3D objects that can have strong gradients in latitude, longitude, and altitude. Secondary eclipse mapping is a powerful way to map the 3D distribution of the atmosphere, but the data can have large correlations and errors in the presence of photon and instrument noise. We develop a technique to mitigate the large uncertainties of eclipse maps by identifying a s…
▽ More
Planetary atmospheres are inherently 3D objects that can have strong gradients in latitude, longitude, and altitude. Secondary eclipse mapping is a powerful way to map the 3D distribution of the atmosphere, but the data can have large correlations and errors in the presence of photon and instrument noise. We develop a technique to mitigate the large uncertainties of eclipse maps by identifying a small number of dominant spectra to make them more tractable for individual analysis via atmospheric retrieval. We use the eigencurves method to infer a multi-wavelength map of a planet from spectroscopic secondary eclipse light curves. We then apply a clustering algorithm to the planet map to identify several regions with similar emergent spectra. We combine the similar spectra together to construct an "eigenspectrum" for each distinct region on the planetary map. We demonstrate how this approach could be used to isolate hot from cold regions and/or regions with different chemical compositions in observations of hot Jupiters with the James Webb Space Telescope (JWST). We find that our method struggles to identify sharp edges in maps with sudden discontinuities, but generally can be used as a first step before a more physically motivated modeling approach to determine the primary features observed on the planet.
△ Less
Submitted 5 October, 2020;
originally announced October 2020.
-
A Comprehensive Reanalysis of $\textit{Spitzer}$'s 4.5 $μ$m Phase Curves, and the Phase Variations of the Ultra-hot Jupiters MASCARA-1b and KELT-16b
Authors:
Taylor J. Bell,
Lisa Dang,
Nicolas B. Cowan,
Jacob Bean,
Jean-Michel Désert,
Jonathan J. Fortney,
Dylan Keating,
Eliza Kempton,
Laura Kreidberg,
Michael R. Line,
Megan Mansfield,
Vivien Parmentier,
Kevin B. Stevenson,
Mark Swain,
Robert T. Zellem
Abstract:
We have developed an open-source pipeline for the analysis of \textit{Spitzer}/IRAC channel 1 and 2 time-series photometry, incorporating some of the most popular decorrelation methods. We applied this pipeline to new phase curve observations of ultra-hot Jupiters MASCARA-1b and KELT-16b, and we performed the first comprehensive reanalysis of 15 phase curves. We find that MASCARA-1b and KELT-16b h…
▽ More
We have developed an open-source pipeline for the analysis of \textit{Spitzer}/IRAC channel 1 and 2 time-series photometry, incorporating some of the most popular decorrelation methods. We applied this pipeline to new phase curve observations of ultra-hot Jupiters MASCARA-1b and KELT-16b, and we performed the first comprehensive reanalysis of 15 phase curves. We find that MASCARA-1b and KELT-16b have phase offsets of $6^{+11}_{-11}~^{\circ}$W and $38^{+16}_{-15}~^{\circ}$W, dayside temperatures of $2952^{+100}_{-97}$ K and $3070^{+160}_{-150}$ K, and nightside temperatures of $1300^{+340}_{-340}$ K and $1900^{+430}_{-440}$ K, respectively. We confirm a strong correlation between dayside and irradiation temperatures with a shallower dependency for nightside temperature. We also find evidence that the normalized phase curve amplitude (peak-to-trough divided by eclipse depth) is correlated with stellar effective temperature. In addition, while our different models often retrieve similar parameters, significant differences occasionally arise between them, as well as between our preferred model and the literature values. Nevertheless, our preferred models are consistent with published phase offsets to within $-8\pm21$ degrees ($-1.6\pm3.2$ sigma), and normalized phase curve amplitudes are on average reproduced to within $-0.01\pm0.24$ ($-0.1\pm1.6$ sigma). Finally, we find that BLISS performs best in most cases, but not all; we therefore recommend future analyses consider numerous detector models to ensure an optimal fit and to assess model dependencies.
△ Less
Submitted 25 February, 2021; v1 submitted 1 October, 2020;
originally announced October 2020.
-
Smaller than expected bright-spot offsets in Spitzer phase curves of the hot Jupiter Qatar-1b
Authors:
Dylan Keating,
Kevin B. Stevenson,
Nicolas B. Cowan,
Emily Rauscher,
Jacob L. Bean,
Taylor Bell,
Lisa Dang,
Drake Deming,
Jean-Michel Désert,
Y. Katherina Feng,
Jonathan J. Fortney,
Tiffany Kataria,
Eliza M. -R. Kempton,
Nikole Lewis,
Michael R. Line,
Megan Mansfield,
Erin May,
Caroline Morley,
Adam P. Showman
Abstract:
We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$m, corresponding to dayside brightness temperatures of…
▽ More
We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~μ$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~μ$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.
△ Less
Submitted 31 March, 2020;
originally announced April 2020.
-
Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b
Authors:
Megan Mansfield,
Jacob L. Bean,
Kevin B. Stevenson,
Thaddeus D. Komacek,
Taylor J. Bell,
Xianyu Tan,
Matej Malik,
Thomas G. Beatty,
Ian Wong,
Nicolas B. Cowan,
Lisa Dang,
Jean-Michel Désert,
Jonathan J. Fortney,
B. Scott Gaudi,
Dylan Keating,
Eliza M. -R. Kempton,
Laura Kreidberg,
Michael R. Line,
Vivien Parmentier,
Keivan G. Stassun,
Mark R. Swain,
Robert T. Zellem
Abstract:
Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves ha…
▽ More
Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.
△ Less
Submitted 7 January, 2020; v1 submitted 3 October, 2019;
originally announced October 2019.
-
Geomorphological Analysis Using Unpiloted Aircraft Systems, Structure from Motion, and Deep Learning
Authors:
Zhiang Chen,
Tyler R. Scott,
Sarah Bearman,
Harish Anand,
Devin Keating,
Chelsea Scott,
J Ramon Arrowsmith,
Jnaneshwar Das
Abstract:
We present a pipeline for geomorphological analysis that uses structure from motion (SfM) and deep learning on close-range aerial imagery to estimate spatial distributions of rock traits (size, roundness, and orientation) along a tectonic fault scarp. The properties of the rocks on the fault scarp derive from the combination of initial volcanic fracturing and subsequent tectonic and geomorphic fra…
▽ More
We present a pipeline for geomorphological analysis that uses structure from motion (SfM) and deep learning on close-range aerial imagery to estimate spatial distributions of rock traits (size, roundness, and orientation) along a tectonic fault scarp. The properties of the rocks on the fault scarp derive from the combination of initial volcanic fracturing and subsequent tectonic and geomorphic fracturing, and our pipeline allows scientists to leverage UAS-based imagery to gain a better understanding of such surface processes. We start by using SfM on aerial imagery to produce georeferenced orthomosaics and digital elevation models (DEM). A human expert then annotates rocks on a set of image tiles sampled from the orthomosaics, and these annotations are used to train a deep neural network to detect and segment individual rocks in the entire site. The extracted semantic information (rock masks) on large volumes of unlabeled, high-resolution SfM products allows subsequent structural analysis and shape descriptors to estimate rock size, roundness, and orientation. We present results of two experiments conducted along a fault scarp in the Volcanic Tablelands near Bishop, California. We conducted the first, proof-of-concept experiment with a DJI Phantom 4 Pro equipped with an RGB camera and inspected if elevation information assisted instance segmentation from RGB channels. Rock-trait histograms along and across the fault scarp were obtained with the neural network inference. In the second experiment, we deployed a hexrotor and a multispectral camera to produce a DEM and five spectral orthomosaics in red, green, blue, red edge, and near infrared. We focused on examining the effectiveness of different combinations of input channels in instance segmentation.
△ Less
Submitted 17 February, 2021; v1 submitted 27 September, 2019;
originally announced September 2019.
-
Uniformly hot nightside temperatures on short-period gas giants
Authors:
Dylan Keating,
Nicolas B. Cowan,
Lisa Dang
Abstract:
Short-period gas giants (hot Jupiters) on circular orbits are expected to be tidally locked into synchronous rotation, with permanent daysides that face their host stars, and permanent nightsides that face the darkness of space. Thermal flux from the nightside of several hot Jupiters has been measured, meaning energy is transported from day to night in some fashion. However, it is not clear exactl…
▽ More
Short-period gas giants (hot Jupiters) on circular orbits are expected to be tidally locked into synchronous rotation, with permanent daysides that face their host stars, and permanent nightsides that face the darkness of space. Thermal flux from the nightside of several hot Jupiters has been measured, meaning energy is transported from day to night in some fashion. However, it is not clear exactly what the physical information from these detections reveals about the atmospheric dynamics of hot Jupiters. Here we show that the nightside effective temperatures of a sample of 12 hot Jupiters are clustered around 1100 K, with a slight upward trend as a function of stellar irradiation. The clustering is not predicted by cloud-free atmospheric circulation models. This result can be explained if most hot Jupiters have nightside clouds that are optically thick to outgoing longwave radiation and hence radiate at the cloud-top temperature, and progressively disperse for planets receiving greater incident flux. Phase curve observations at a greater range of wavelengths are crucial to determining the extent of cloud coverage, as well as the cloud composition on hot Jupiter nightsides.
△ Less
Submitted 30 March, 2020; v1 submitted 31 August, 2018;
originally announced September 2018.
-
Revisiting the Energy Budget of WASP-43b: Enhanced day-night heat transport
Authors:
Dylan Keating,
Nicolas B. Cowan
Abstract:
The large day--night temperature contrast of WASP-43b has so far eluded explanation. We revisit the energy budget of this planet by considering the impact of reflected light on dayside measurements, and the physicality of implied nightside temperatures. Previous analyses of the infrared eclipses of WASP-43b have assumed reflected light from the planet is negligible and can be ignored. We develop a…
▽ More
The large day--night temperature contrast of WASP-43b has so far eluded explanation. We revisit the energy budget of this planet by considering the impact of reflected light on dayside measurements, and the physicality of implied nightside temperatures. Previous analyses of the infrared eclipses of WASP-43b have assumed reflected light from the planet is negligible and can be ignored. We develop a phenomenological eclipse model including reflected light, thermal emission, and water absorption, and use it to fit published Hubble and Spitzer eclipse data. We infer a near-infrared geometric albedo of 27$\pm1\%$ and a cooler dayside temperature of $1527 \pm 10~$K. Additionally, we perform lightcurve inversion on the three published orbital phase curves of WASP-43b and find that each requires unphysical, negative flux on the nightside. By requiring non-negative brightnesses at all longitudes, we correct the unphysical parts of the maps and obtain a much hotter nightside effective temperature of $1076 \pm 11~$K. The cooler dayside and hotter nightside suggests a heat recirculation efficiency of $47\%$ for WASP-43b, essentially the same as for HD 209458b, another hot Jupiter with nearly the same temperature. Our analysis therefore reaffirms the trend that planets with lower irradiation temperatures have more efficient day-night heat transport. Moreover, we note that 1) reflected light may be significant for many near-IR eclipse measurements of hot Jupiters, and 2) phase curves should be fit with physically possible longitudinal brightness profiles --- it is insufficient to only require that the disk-integrated lightcurve be non-negative.
△ Less
Submitted 18 October, 2017; v1 submitted 11 September, 2017;
originally announced September 2017.