Planetary Climate Dynamics

Abstract:

We report the direct imaging discovery of a third exoplanet in the β Pictoris (β Pic) system. We detect β Pictoris d (β Pic d) in noncoronagraphic observations obtained with the Very Large Telescope (VLT) Enhanced Resolution Imager and Spectrograph (ERIS), as well as multi-epoch archival datasets from the JWST Near Infrared Camera (NIRCam) and VLT/SPHERE. Astrometric measurements over an 11 yr baseline demonstrate that it is consistent with a gravitationally bound source with orbital motion. Joint multi-planet orbit fits of all three planets in the system yield a semimajor axis of 26.0−6.1+2.2 au and inclination 89.0−0.6+0.7 deg for planet d. β Pic d has a larger orbital semimajor axis than the other known planets in the system, but is coplanar with the inner two planets, and its orbit is consistent with sculpting the inner edge of the debris disk. β Pic d has a contrast of ΔL′=12.11±0.15 mag, with colors and luminosity that closely match those of 51 Eri b, another exoplanet in the β Pic moving group. Its VLT/ERIS and JWST/NIRCam colors are distinct from those of free-floating planetary-mass objects of a similar age and temperature. Its red F410M − F444W color indicates strong CO2 absorption in its atmosphere and suggests significant enhancement in metals compared to free-floating objects. From the ATMO hot-start evolutionary models, we estimate an effective temperature of 600−60+45 K and mass of 2.4 ± 0.6 MJup, which also closely matches similar estimates for 51 Eri b. β Pic d is among the lowest-mass exoplanets imaged from the ground. This discovery highlights the deep sensitivity achievable with ground-based imaging in the mid-infrared and the discovery potential of future high-contrast observations with the Extremely Large Telescope.

Abstract:

The geochemical evolution of long-lived magma oceans is strongly regulated by volatile exchange between the molten mantle and the atmosphere. For planets inside the runaway-greenhouse limit, this coupled evolution can persist for billions of years, governing bulk density, surface conditions, and long-term geodynamics. However, most existing studies assume Earth-like (oxidized) conditions and neglect the influence of redox state on melt thermodynamics and volatile release. We quantify how experimentally derived, oxygen-fugacity-dependent melting curves implemented within the coupled interior-atmosphere framework PROTEUS propagate into the thermal structure, melt fraction, and rheological evolution of rocky exoplanet interiors, applying this to the short-period super-Earth GJ 1132 b. We find strongly non-linear thermal responses to variations in melting curves. In volatile-poor systems, reduced meltingcurves ( f O2 ≤ IW) promote earlier deep-mantle crystallisation relative to oxidised (IW + 2.0) and Earth-like (IW+4.0) cases (range IW−4.0 to IW+4.0), favouring late-stage surface magma oceans sustained by greenhouse warming, while oxidized melting curves maintain higher melt fractions and a vertically extended magma ocean. Reduced mantles produce massive H2-CO-rich atmospheres; oxidized mantles favour thinner H2O-CO2 envelopes. In volatile-rich systems, the interior reaches radiative equilibrium at high melt fractions, sustaining a steady-state global magma ocean in which melting curve variations do not significantly influence solidification timing. This indicates a hierarchical control: volatile inventory and surface oxygen fugacity act as the primary regulators of thermal state, while oxygen-fugacity-dependent melting relations provide a secondary modulation. These contrasting regimes produce distinct atmospheric compositions and formation timescales, offering testable spectral predictions for close-in rocky exoplanets evaluable with forthcoming JWST observations.

Abstract:

Tidally locked lava planets are hot, rocky worlds on close-in orbits with a permanent molten dayside. With JWST, their surfaces and atmospheres are beginning to be revealed. This work investigates three-dimensional (3D) magma ocean dynamics, derives scaling laws for the resulting ocean heat transport, and predicts its detectability. For the first time, the ocean circulation driven by the intense momentum and mass exchanges with the supersonic atmosphere is considered in addition to that by thermal forcing. The wind forcing turns out to overwhelmingly dominate the other two mechanisms, driving ocean currents reaching ∼100 m s⁻¹ and greatly expanding the latitudinal extent of the Matsuno–Gill response. Despite these extreme flow speeds, scaling analysis and 3D simulations consistently demonstrate that magma ocean circulation alone does not produce an observable hot-spot offset. This inefficiency arises because basin geometry and circulation structure fundamentally constrain zonal heat redistribution, suppressing large-scale longitudinal transport even under vigorous flow.

Abstract:

Ultra-hot Jupiters (UHJs) present extreme atmospheric phenomena not found in the solar system. These planets’ daysides experience strong temperature inversions, molecular species (including H2) dissociate, and magnetism effects disrupt atmospheric circulation. On their nightsides, H2 can recombine and clouds may form. Spectroscopic phase curves allow us to measure these spatially inhomogeneous conditions, which can then be interpreted with three-dimensional (3D) models. In this work, we compare the JWST Near Infrared Imager and Slitless Spectrograph (or NIRISS) spectroscopic phase curve of the UHJ WASP-121 b to state-of-the-art 3D models with varying modeling assumptions, including the aforementioned physical phenomena. We demonstrate the importance of accurately accounting for the planet’s radius when comparing data and models, as it changes the implied overall planetary emission. We find that the 3D models predict planetary emission that is ∼12% higher than observed, contributing to the continued tension between measured and predicted hot Jupiter albedos. We identify multiple lines of evidence confirming a strong source of drag operating in this planet’s atmosphere. In addition, the nightside emission spectrum is devoid of strong absorption features, which may be best explained by nightside clouds. One feature of the dataset that is not reproduced by the 3D models is a trend toward increasingly eastward phase offsets at wavelengths shorter than ∼1.4 μm. This result is inconsistent with reflection from dayside clouds, nor can it be explained by removing atmospheric opacity sources. Our analysis highlights the complexities involved in generating 3D models and interpreting observations of UHJs in the JWST era.