Super-Earth lava planet from birth to observation: photochemistry, tidal heating, and volatile-rich formation

Copernicus Publications (2025)

Authors:

Harrison Nicholls, Tim Lichtenberg, Richard D Chatterjee, Claire Marie Guimond, Emma Postolec, Raymond T Pierrehumbert

Abstract:

Larger-than-Earth exoplanets are sculpted by strong stellar irradiation, but it is unknown whence they originate. Two propositions are that they formed with rocky interiors and hydrogen-rich envelopes (鈥榞as-dwarfs鈥), or with bulk compositions rich in water-ices (鈥榳ater-worlds鈥) . Multiple observations of super-Earth L 98-59 d have revealed its low bulk-density, consistent with substantial volatile content alongside a rocky/metallic interior, and recent JWST spectroscopy evidences a high mean molecular weight atmosphere. Its density and composition make it a waymarker for disentangling the processes which separate super-Earths and sub-Neptunes across geological timescales. We simulate the possible pathways for L 98-59 d from birth up to the present day using a comprehensive evolutionary modelling framework. Emerging from our calculations is a novel self-limiting mechanism between radiative cooling, tidal heating, and mantle rheology, which we term the 'radiation-tide-rheology feedback'. Coupled numerical modelling yields self-limiting tidal heating estimates that are up to two orders of magnitude lower than previous calculations, and yet are still large enough to enable the extension of primordial magma oceans to Gyr timescales. Our analysis indicates that the planet formed with a large amount (>1.8 mass%) of sulfur and hydrogen, and a chemically-reducing mantle; inconsistent with both the canonical gas-dwarf and water-world scenarios. A thick atmosphere and tidal heating sustain a permanent deep magma ocean, allowing the dissolution and retention of volatiles within its mantle. Transmission features can be explained by in-situ photochemical production of SO2 in a high-molecular weight H2-H2S background. These results subvert the emerging gas-dwarf vs. water-world dichotomy of small planet categorisation, inviting a more nuanced classification framework. We show that interactions between planetary interiors and atmospheres shape their observable characteristics over billions of years.

The Cosmic Shoreline Revisited: A Metric for Atmospheric Retention Informed by Hydrodynamic Escape

Copernicus Publications (2025)

Authors:

Xuan Ji, Richard Chatterjee, Brandon Park Coy, Edwin Kite

Abstract:

The 鈥渃osmic shoreline鈥, a semi-empirical relation that separates airless worlds from worlds with atmospheres as proposed by Zahnle & Catling (2017), is now guiding large-scale JWST surveys aimed at detecting rocky exoplanet atmospheres. We expand upon this framework by revisiting the shorelines using existing hydrodynamic escape models applied to Earth-like, Venus-like, and steam atmospheres for rocky exoplanets, and we estimate energy-limited escape rates for CH4 atmospheres. We determine the critical instellation required for atmospheric retention by calculating time-integrated atmospheric mass loss. Our analysis introduces a new metric for target selection in the Rocky Worlds DDT and refines expectations for rocky planet atmosphere searches in Cycle 4. Exploring initial volatile inventory ranging from 0.01% to 1% of planetary mass, we find that its variation prevents the definition of a unique clear-cut shoreline, though non-linear escape physics can reduce this sensitivity to initial conditions. Additionally, uncertain distributions of high-energy stellar evolution and planet age further blur the critical instellations for atmospheric retention, yielding broad shorelines. Hydrodynamic escape models find atmospheric retention is markedly more favorable for higher-mass planets orbiting higher-mass stars, with carbon-rich atmospheres remaining plausible for 55 Cancri e despite its extreme instellation. Dedicated modelling efforts are needed to better constrain the escape dynamics of secondary atmospheres, such as the role of atomic line cooling, especially for Earth-sized planets. Finally, we illustrate how density measurements can be used to statistically test the existence of the cosmic shorelines, emphasizing the need for more precise mass and radius measurements.

The Role of Bottom Thermal Forcing on Modulating Baroclinic Instability in a Jupiter GCM

Copernicus Publications (2025)

Authors:

Xinmiao Hu, Peter Read, Roland Young, Greg Colyer

Abstract:

The atmospheric circulation of Jupiter is shaped by a complex interplay between deep internal processes and cloud-level dynamics. Numerical simulations and observational analyses have suggested that Jupiter鈥檚 mid-latitude jets are strongly influenced by baroclinic instability [1], which is governed by the planet鈥檚 atmospheric thermal structure. Jupiter emits a substantial intrinsic heat flux originating from its interior. Past modelling efforts [2, 3] have demonstrated that this internal energy plays a key role in shaping large-scale atmospheric dynamics.Our previous work [4] showed that latitudinal variations in interior heat flux can significantly impact the structure and behaviour of Jupiter鈥檚 mid-latitude jets in a General Circulation Model (GCM).聽 Such an impact is best illustrated by the relative vorticity snapshots from two simulations with the lowest and highest latitudinal flux gradient (see Figure鈥1). In this study, we present a more detailed analysis linking these jet modifications to changes in the atmospheric thermal structure and, consequently, to the strength and distribution of baroclinic eddy activity. In particular, we use the Lorenz energy cycle framework to diagnose how variations in deep thermal forcing influence baroclinic energy conversion and eddy-mean flow interactions. We further examine the implications for meridional transport and the water cycle within Jupiter鈥檚 weather layer.Additionally, we present a control simulation in which the potential temperature at the model鈥檚 lower boundary is forced toward a fixed value (a deep adiabat setup). We compute the equivalent upward heat flux associated with this forcing to place it in the context of previous models that impose constant or latitudinally varying interior heat flux. This allows a direct comparison of how different representations of deep thermal forcing affect upper-atmospheric dynamics.Finally, we discuss the broader implications of these findings for future weather-layer models of Jupiter and other gas giant planets, especially on the effect of bottom boundary conditions in representing the coupling between deep and observable atmospheric dynamics.聽Figure 1: Mollweide projection of the relative vorticity at 1 bar at the end of two simulations.Reference:[1] Read, P. L. (2023). The dynamics of Jupiter鈥檚 and Saturn鈥檚 weather layers: a synthesis after Cassini and Juno. Annual Review of Fluid Mechanics, 56(1), 271鈥293. https://doi.org/10.1146/annurev-fluid-121021-040058[2] Liu, J., & Schneider, T. (2011). Convective Generation of Equatorial Superrotation in Planetary Atmospheres. Journal of the Atmospheric Sciences, 68(11), 2742-2756. https://doi.org/10.1175/JAS-D-10-05013.1[3] Young, R. M. B., Read, P. L., & Wang, Y. (2018). Simulating Jupiter鈥檚 weather layer. Part I: Jet spin-up in a dry atmosphere. Icarus, 326, 225鈥252. https://doi.org/10.1016/j.icarus.2018.12.005[鈥4] Hu, X. and Read, P.: Latitudinal Variation in Internal Heat Flux in Jupiter's Atmosphere: Effect on Weather Layer Dynamics, Europlanet Science Congress 2024, EPSC2024-669, https://doi.org/10.5194/epsc2024-669, 2024.

What are subNeptunes made of?

(2025)

Abstract:

This talk will cover the state of the art in whole-planet subNeptune modelling, and needs for the future.聽 Inferences about the composition of the deep envelope can be made on the basis of the way chemical transformations in the deep envelope may be evidenced in the observable atmosphere, such as has been attempted, for example, regarding the presence or absence of NH3 in the observable atmospheres of subNeptunes.聽 Such inferences require an understanding not only of deep envelope chemistry, but also of vertical mixing processes. The mixing process engages a number of poorly understood phenomena, such as mixing rates through stably stratified (nonconvective) internal radiative layers.聽 The occurrence of such radiative layers can be induced by compositional suppression of convection (e.g. due to high molecular weight H2O in an H2-rich atmosphere). We will review our modelling studies regarding this phenomenon.聽 Typically, the envelope-silicate interface is hot enough that the interface takes the form of a magma ocean, so compositional interchange with the magma ocean becomes crucial. This exchange includes rock vapours as well as lower molecular weight volatiles.聽 Our work on magma ocean exchanges will be reviewed. We highlight the importance of mineral physics experiments and molecular dynamics to provide crucially needed (and largely absent) thermodynamic parameters, particularly at high pressure.聽 At sufficiently high temperatures, silicate itself can become supercritical so that the distinction between silicate melt and silicate vapour disappears and the silicate substance becomes completely miscible with the lower molecular weight envelope.聽 Modeling and experiment regarding this novel and largely unexplored regime is particularly needed.

Absence of a Runaway Greenhouse Limit on Lava Planets

The Astrophysical Journal American Astronomical Society 987:2 (2025) 172

Authors:

Iris D Boer, Harrison Nicholls, Tim Lichtenberg

Abstract:

Climate transitions on exoplanets offer valuable insights into the atmospheric processes governing planetary habitability. Previous pure-steam atmospheric models show a thermal limit in outgoing long-wave radiation, which has been used to define the inner edge of the classical habitable zone and guide exoplanet surveys aiming to identify and characterize potentially habitable worlds. We expand upon previous modelling by treating (i) the dissolution of volatiles into a magma ocean underneath the atmosphere, (ii) a broader volatile range of the atmospheric composition including H2O, CO2, CO, H2, CH4, and N2, and (iii) a surface-temperature- and mantle-redox-dependent equilibrium chemistry. We find that multicomponent atmospheres of outgassed composition located above partially or fully molten mantles do not exhibit the characteristic thermal radiation limit that arises from pure-steam models, thereby undermining the canonical concept of a runaway greenhouse limit, and hence challenging the conventional approach of using it to define an irradiation-based habitable zone. Our results show that atmospheric heat loss to space is strongly dependent on the oxidation and melting state of the underlying planetary mantle, through their significant influence on the atmosphere鈥檚 equilibrium composition. This suggests an evolutionary hysteresis in climate scenarios: Initially molten and cooling planets do not converge to the same climate regime as solidified planets that heat up by external irradiation. Steady-state models cannot recover evolutionary climate transitions, which instead require self-consistent models of the temporal evolution of the coupled feedback processes between interior and atmosphere over geologic time.