skip to main content

Title: Variability due to climate and chemistry in observations of oxygenated Earth-analogue exoplanets

The Great Oxidation Event was a period during which Earth’s atmospheric oxygen (O2) concentrations increased from ∼10−5 times its present atmospheric level (PAL) to near modern levels, marking the start of the Proterozoic geological eon 2.4 billion years ago. Using WACCM6, an Earth System Model, we simulate the atmosphere of Earth-analogue exoplanets with O2 mixing ratios between 0.1 and 150 per cent PAL. Using these simulations, we calculate the reflection spectra over multiple orbits using the Planetary Spectrum Generator. We highlight how observer angle, albedo, chemistry, and clouds affect the simulated observations. We show that inter-annual climate variations, as well short-term variations due to clouds, can be observed in our simulated atmospheres with a telescope concept such as LUVOIR or HabEx. Annual variability and seasonal variability can change the planet’s reflected flux (including the reflected flux of key spectral features such as O2 and H2O) by up to factors of 5 and 20, respectively, for the same orbital phase. This variability is best observed with a high-throughput coronagraph. For example, HabEx (4 m) with a starshade performs up to a factor of two times better than a LUVOIR B (6 m) style telescope. The variability and signal-to-noise ratio of some spectral features depends non-linearly on atmospheric O2 concentration. This is caused by temperature and chemical column depth variations, as well as generally increased liquid and ice cloud content for atmospheres with O2 concentrations of <1 per cent PAL.

more » « less
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 206-219
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    Clouds and other features in exoplanet and brown dwarf atmospheres cause variations in brightness as they rotate in and out of view. Ground-based instruments reach the high contrasts and small inner working angles needed to monitor these faint companions, but their small fields of view lack simultaneous photometric references to correct for non-astrophysical variations. We present a novel approach for making ground-based light curves of directly imaged companions using high-cadence differential spectrophotometric monitoring, where the simultaneous reference is provided by a double-grating 360○ vector Apodizing Phase Plate (dgvAPP360) coronagraph. The dgvAPP360 enables high-contrast companion detections without blocking the host star, allowing it to be used as a simultaneous reference. To further reduce systematic noise, we emulate exoplanet transmission spectroscopy, where the light is spectrally dispersed and then recombined into white-light flux. We do this by combining the dgvAPP360 with the infrared Arizona Lenslets for Exoplanet Spectroscopy integral field spectrograph on the Large Binocular Telescope Interferometer. To demonstrate, we observed the red companion HD 1160 B (separation ∼780 mas) for one night, and detect $8.8{{\ \rm per\ cent}}$ semi-amplitude sinusoidal variability with an ∼3.24 h period in its detrended white-light curve. We achieve the greatest precision in ground-based high-contrast imaging light curves of sub-arcsecond companions to date, reaching $3.7{{\ \rm per\ cent}}$ precision per 18-min bin. Individual wavelength channels spanning 3.59–3.99 μm further show tentative evidence of increasing variability with wavelength. We find no evidence yet of a systematic noise floor; hence, additional observations can further improve the precision. This is therefore a promising avenue for future work aiming to map storms or find transiting exomoons around giant exoplanets.

    more » « less
  2. Abstract Transmission spectroscopy offers an invaluable opportunity to characterize the atmospheres of exoplanets. We present new ground-based optical transmission spectra of the hot Jupiter HD 189733b, derived from nine transits observed over a six year time span (2016–2021) using near-simultaneous u ′ g ′ r ′ i ′ broadband observations. We achieve an average (best) precision of 435 (280) ppm by implementing an optical diffuser on the prime focus spectrograph from the 2.3 m Wyoming Infrared Observatory telescope. The data provide new measurements of the apparent planetary radius with respect to the stellar radius, the spectral index of atmospheric opacity, and the time variability of the two quantities. Our results indicate an enhanced spectral slope in the optical regime ≈2.4 times steeper than would be expected from canonical Rayleigh scattering and that is consistent with earlier measurements of a super-Rayleigh slope (SRS). While the effect of stellar activity on the transmission spectrum complicates the measurement of the spectral slope, our multiepoch data set over six years can measure and average over stellar variations, yielding a mean spectral index of −9.9 ± 4.4. The 1200 K equilibrium temperature of HD 189733b places it in a sweet spot for the formation of SRSs and is consistent with vigorously mixing hazes in the atmosphere. Additionally, we find variations in the depth of the lightcurve during two of the transits, explainable as an increase in occulted star spots during June 2021. Although the star is active, the mean level of stellar activity does not seem to vary dramatically over our six years of observations, leading us to conclude that the variability in stellar activity is modest at most. 
    more » « less
  3. Abstract The photometric and spectral variability of brown dwarfs probes heterogeneous temperature and cloud distributions and traces the atmospheric circulation patterns. We present a new 42 hr Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectral time series of VHS 1256-1257 b, a late L-type planetary-mass companion that has been shown to have one of the highest variability amplitudes among substellar objects. The light curve is rapidly evolving and best fit by a combination of three sine waves with different periods and a linear trend. The amplitudes of the sine waves and the linear slope vary with the wavelength, and the corresponding spectral variability patterns match the predictions by models invoking either heterogeneous clouds or thermal profile anomalies. Combining these observations with previous HST monitoring data, we find that the peak-to-valley flux difference is 33% ± 2% with an even higher amplitude reaching 38% in the J band, the highest amplitude ever observed in a substellar object. The observed light curve can be explained by maps that are composed of zonal waves, spots, or a mixture of the two. Distinguishing the origin of rapid light curve evolution requires additional long-term monitoring. Our findings underscore the essential role of atmospheric dynamics in shaping brown-dwarf atmospheres and highlight VHS 1256-1257 b as one of the most favorable targets for studying the atmospheres, clouds, and atmospheric circulation of planets and brown dwarfs. 
    more » « less
  4. Abstract. The tropical tropopause layer (TTL) is a sea of vertical motions. Convectively generated gravity waves create vertical winds on scales of a few to thousands of kilometers as they propagate in a stable atmosphere. Turbulence from gravity wave breaking, radiatively driven convection, and Kelvin–Helmholtz instabilities stirs up the TTL on the kilometer scale. TTL cirrus clouds, which moderate the water vapor concentration in the TTL and stratosphere, form in the cold phases of large-scale (> 100 km) wave activity. It has been proposed in several modeling studies that small-scale (< 100 km) vertical motions control the ice crystal number concentration and the dehydration efficiency of TTL cirrus clouds. Here, we present the first observational evidence for this. High-rate vertical winds measured by aircraft are a valuable and underutilized tool for constraining small-scale TTL vertical wind variability, examining its impacts on TTL cirrus clouds, and evaluating atmospheric models. We use 20 Hz data from five National Aeronautics and Space Administration (NASA) campaigns to quantify small-scale vertical wind variability in the TTL and to see how it varies with ice water content, distance from deep convective cores, and height in the TTL. We find that 1 Hz vertical winds are well represented by a normal distribution, with a standard deviation of 0.2–0.4 m s−1. Consistent with a previous observational study that analyzed two out of the five aircraft campaigns that we analyze here, we find that turbulence is enhanced over the tropical west Pacific and within 100 km of convection and is most common in the lower TTL (14–15.5 km), closer to deep convection, and in the upper TTL (15.5–17 km), further from deep convection. An algorithm to classify turbulence and long-wavelength (5 km < λ < 100 km) and short-wavelength (λ < 5 km) gravity wave activity during level flight legs is applied to data from the Airborne Tropical TRopopause EXperiment (ATTREX). The most commonly sampled conditions are (1) a quiescent atmosphere with negligible small-scale vertical wind variability, (2) long-wavelength gravity wave activity (LW GWA), and (3) LW GWA with turbulence. Turbulence rarely occurs in the absence of gravity wave activity. Cirrus clouds with ice crystal number concentrations exceeding 20 L−1 and ice water content exceeding 1 mg m−3 are rare in a quiescent atmosphere but about 20 times more likely when there is gravity wave activity and 50 times more likely when there is also turbulence, confirming the results of the aforementioned modeling studies. Our observational analysis shows that small-scale gravity waves strongly influence the ice crystal number concentration and ice water content within TTL cirrus clouds. Global storm-resolving models have recently been run with horizontal grid spacing between 1 and 10 km, which is sufficient to resolve some small-scale gravity wave activity. We evaluate simulated vertical wind spectra (10–100 km) from four global storm-resolving simulations that have horizontal grid spacing of 3–5 km with aircraft observations from ATTREX. We find that all four models have too little resolved vertical wind at horizontal wavelengths between 10 and 100 km and thus too little small-scale gravity wave activity, although the bias is much less pronounced in global SAM than in the other models. We expect that deficient small-scale gravity wave activity significantly limits the realism of simulated ice microphysics in these models and that improved representation requires moving to finer horizontal and vertical grid spacing. 
    more » « less
  5. ABSTRACT We demonstrate that the James Webb Space Telescope (JWST) can detect infrared (IR) excess from the blended light spectral energy distribution of spatially unresolved terrestrial exoplanets orbiting nearby white dwarfs. We find that JWST is capable of detecting warm (habitable-zone; Teq = 287 K) Earths or super-Earths and hot (400–1000 K) Mercury analogues in the blended light spectrum around the nearest 15 isolated white dwarfs with 10 h of integration per target using MIRI’s medium-resolution spectrograph (MRS). Further, these observations constrain the presence of a CO2-dominated atmosphere on these planets. The technique is nearly insensitive to system inclination, and thus observation of even a small sample of white dwarfs could place strong limits on the occurrence rates of warm terrestrial exoplanets around white dwarfs in the solar neighbourhood. We find that JWST can also detect exceptionally cold (100–150 K) Jupiter-sized exoplanets via MIRI broad-band imaging at $\lambda = 21\, \mathrm{\mu m}$ for the 34 nearest (<13 pc) solitary white dwarfs with 2 h of integration time per target. Using IR excess to detect thermal variations with orbital phase or spectral absorption features within the atmosphere, both of which are possible with long-baseline MRS observations, would confirm candidates as actual exoplanets. Assuming an Earth-like atmospheric composition, we find that the detection of the biosignature pair O3+CH4 is possible for all habitable-zone Earths (within 6.5 pc; six white dwarf systems) or super-Earths (within 10 pc; 17 systems) orbiting white dwarfs with only 5–36 h of integration using MIRI’s low-resolution spectrometer. 
    more » « less