Planetary Climate Dynamics

Abstract:

To aid the search for atmospheres on rocky exoplanets, we should know what to look for. An unofficial paradigm is to anticipate CO2 present in these atmospheres, through analogy to the solar system and through theoretical modelling. This CO2 would be outgassed from molten silicate rock produced in the planet’s mostly-solid interior—an ongoing self-cooling mechanism that should proceed, in general, so long as the planet has sufficient internal heat to lose.Outgassing of CO2 requires relatively oxidising conditions. Previous work has noted the importance of how oxidising the planet interior is (the oxygen fugacity), which depends strongly on its rock composition. Current models presume that redox reactions between iron species control oxygen fugacity. However, iron alone need not be the sole dictator of how oxidising a planet is. Indeed, carbon itself is a powerful redox element, with great potential to feed back upon the mantle redox state as it melts. Whilst Earth is carbon-poor, even a slightly-higher volatile endowment could trigger carbon-powered geochemistry.We offer a new framework for how carbon is transported from solid planetary interior to atmosphere. The model incorporates realistic carbon geochemistry constrained by recent experiments on CO2 solubility in molten silicate, as well as redox couplings between carbon and iron that have never before been applied to exoplanets. We also incorporate a coupled 1D energy- and mass-balance model to provide first-order predictions of the rate of volcanism.We show that carbon-iron redox coupling maintains interior oxygen fugacity in a narrow range: more reducing than Earth magma, but not reducing enough to destabilise CO2 gas. We predict that most secondary atmospheres, if present, should contain CO2, although the total pressure could be low. An atmospheric non-detection may indicate a planet either born astonishingly dry, or having shut off its internal heat engine.

Abstract:

Recent analyses of wind measurements obtained from tracking cloud motions in spacecraft images of Jupiter and Saturn[1,2] indicate that nonlinear scale-to-scale transfers of kinetic energy act from small to large scales over a wide range of length scales, much as anticipated for 2D or geostrophic turbulence paradigms. At the smallest resolvable scales, however, there is evidence in observations of a forward (downscale) transfer, at least at low and middle latitudes on Jupiter, much like in the Earth’s atmosphere. Moreover, the upscale transfers at the largest spatial scales are evidently dominated by spectrally non-local, highly anisotropic eddy-zonal interactions associated with the generation of intense zonal jets and equatorial super-rotation by direct eddy-zonal flow exchanges. Most analyses to date have emphasised the global mean interactions for both planets, thereby focusing on the spatially homogeneous and isotropic components of the turbulence. Here we present some new analyses of spectral energy transfers on both Jupiter and Saturn that resolve variations in latitude[cf 3], thereby shedding new light on non-homogeneous aspects of jovian turbulent interactions. The results indicate significant variability and inhomogeneity between different locations, with a clear distinction between the tropics, the extratropical middle latitudes and the polar regions. We discuss these in light of other observations and models of gas giant circulation and related laboratory experimental analogues.[1] Antu˜nano, A., del Río-Gaztelurrutia, T., Sánchez-Lavega, A., & Hueso, R. (2015) Dynamics of Saturn’s polar regions. J. Geophys. Res.: Planets, 120 , 155–176. doi:10.1002/2014JE004709[2] Read, P. L., Antu˜nano, A., Cabanes, S., Colyer, G., del Río-Gaztelurrutia, T., Sánchez-Lavega, A. (2022). Energy exchanges in Saturn’s polar regions fromCassini observations: Eddy-zonal flow interactions. J. Geophys. Res., 127 , e2021JE006973. https://doi.org/10.1029/2021JE006973[3] Chemke, R., & Kaspi, Y. (2015). The latitudinal dependence of atmospheric jet scales and macroturbulent energy cascades. J. Atmos. Sci., 72 , 3891–3907. doi: 10.1175/JAS-D-15-0007.

Abstract:

The upcoming generation of extremely large ground-based telescopes and space observatories promises to transform our understanding of rocky exoplanets within the habitable zones of their stars. Observations from the James Webb Space Telescope are already challenging current models of exoplanetary atmospheres and interiors (e.g., Foley 2024). Conventionally, a habitable zone exoplanet orbits within a circumstellar region where liquid water could potentially exist on its surface (e.g., Hart 1979; Kasting et al. 1993). It is generally assumed that silicate weathering regulates atmospheric carbon dioxide (CO2) to levels supporting liquid water, providing a stabilizing feedback on climate (Walker et al. 1981); otherwise, planetary habitability would be a matter of luck. This negative feedback underpins the concept of the circumstellar habitable zone (CHZ) and may play a critical role in climate regulation on water-bearing, tectonically active rocky exoplanets. Identifying evidence for a carbon cycle on exoplanets and verifying the validity of the habitable zone concept are key objectives for future research (e.g., Bean et al. 2017)—efforts that would benefit from a deeper understanding of the carbonate–silicate cycle. In particular, the interplay between global climate, atmospheric CO2, and silicate weathering rates is not fully understood. While the role of continental weathering in this negative feedback has been extensively studied, seafloor weathering has received comparatively less attention despite its potential to be equally significant. For instance, during the Late Mesozoic—known for its hothouse climate state—seafloor weathering fluxes were comparable in magnitude to those of continental weathering (e.g., Coogan & Gillis 2013). Here, I explore the factors controlling basalt dissolution by applying the more mechanistic weathering model of Maher and Chamberlain (2014) to both the modern and Late Mesozoic upper oceanic crust.The oceanic crust, through its formation, alteration, and subduction, is a key component of this geochemical cycle. Carbon is transferred from Earth’s mantle to surface reservoirs (e.g., the atmosphere and oceans) via volcanic outgassing. Concurrently, basalt dissolution and carbonate precipitation recapture dissolved carbon from seawater, storing it within the oceanic crust. Over geological timescales, subducted tectonic plates carry these carbonates back into the mantle reservoir, ultimately supplying carbon to volcanoes. Nonetheless, the primary controls on seafloor weathering rates remain debated. Most models have historically focused on the dependence of mineral dissolution kinetics on temperature and CO2 concentrations, yet observational data suggest that both kinetic and thermodynamic factors govern global weathering fluxes, underscoring the need for models that can incorporate this dual control. Kinetic weathering models (e.g., Walker et al. 1981) fail to account for the changes in Earth’s weatherability through time (e.g., West et al. 2005). To address this, Maher & Chamberlain (2014) developed a solute transport model that integrates hydrological and tectonic influences and imposes a thermodynamic limit on weathering rates. However, this model has only occasionally been applied to continental weathering in exoplanet climate studies (e.g., Graham & Pierrehumbert 2020, 2024), and its relevance to seafloor weathering remains largely unexplored (Hakim et al. 2021).In this work, I assessed the model’s sensitivity to key parameters by comparing its predictions with observed age-dependent carbon content in the upper oceanic crust (Gillis & Coogan 2011). I then extended the model into two dimensions to represent the evolving age distribution of Earth’s seafloor and examine its impact on global weathering fluxes. Simulations using this supply-limited seafloor weathering model successfully reproduce observed age-related trends in carbon content within Earth’s upper oceanic crust and provide further evidence that over 80% of carbonate formed within 20 Myr of crust formation (Gillis & Coogan, 2011; Albers et al. 2023). This model can also capture the higher CO2 concentrations in Late Mesozoic crust compared to Cenozoic crust. Our results indicate that crustal age, permeability, porosity and fluid flow strongly influence weathering rates. These findings challenge the prevailing view that elevated temperatures primarily drove enhanced carbon uptake during the Late Mesozoic, emphasizing further the importance of incorporating geologic and hydrologic processes into climate models. Consequently, seafloor weathering emerges as a necessary process not only for understanding Earth’s past climate but also for interpreting future observations of potentially habitable rocky exoplanets.

Abstract:

Recent studies of the climate and weather on Mars have benefitted greatly from the development of publicly available “reanalysis” datasets. These are quantitative reconstructions of the three-dimensional, multivariate, time-varying state of the Martian atmosphere, obtained by combining observations of atmospheric variables, usually from remote sensing platforms in orbit around the planet, with a numerical model simulation of the entire atmospheric circulation and surface, using a statistical-dynamical algorithm known as an “assimilation”. The result represents an estimate of the evolving state of the entire global weather and climate that makes optimal use of both direct observations and physical knowledge of the physics and chemistry of the atmosphere and surface (as contained in the design of the numerical model), taking into account statistical uncertainties in both measurements and model simulations. This approach thereby takes account of measurement uncertainty, the incomplete coverage of observations in both space and time and enables the estimation of variables that cannot be measured directly through the internal dynamical consistency of the numerical model. At least four such datasets have been produced during the past decade or so [1-6], based on recent spacecraft observations from the Thermal Emission Spectrometer (TES) on NASA’s Mars Global Surveyor (MGS), the Mars Climate Sounder (MCS) instrument on NASA’s Mars Reconnaissance Orbiter (MRO), the Atmospheric Chemistry Suite thermal infrared channel (ACS-TIRVIM) on ESA's Trace Gas Orbiter and the Emirates Mars InfraRed Spectrometer (EMIRS) on the Emirates Mars Mission (EMM), using different assimilation algorithms. In most cases so far, the data assimilated has been limited to retrieved surface temperature and atmospheric temperature profiles, together with column integrated dust opacities (CIDO) and certain other trace gases. Although this leads to reasonable agreement between different assimilations and with independent, out of sample measurements, the neglect of measurements of the vertical structure of atmospheric dust loading leads to significant errors in modelled radiative heating and cooling rates that temperature assimilation is then required to correct. Such a correction is highly undesirable, quite apart from misrepresenting the structure and transport of dust aerosol by the circulation, and increases the likelihood of systematic errors in the reconstructed circulation. In recent work [4], the Analysis Correction system used for the Mars Analysis Correction Assimilation (MACDA) reanalysis [1] has been extended to enable the assimilation of both CIDO and profiles of dust opacity obtained from limb-sounding instruments such as MCS. The results have enabled phenomena such as the climatological elevated dust layers during northern hemisphere summer to be captured in a reanalysis (see Fig. 1), which has been elusive in previous work. In new work presented here, this approach has now been applied to more than 12 Mars years of observations, from early MY24 to the beginning of MY36 (to date), based mainly on retrievals of temperature and dust opacity from the MGS/TES and MRO/MCS instruments that have been assimilated into the UK version of the LMD Mars Planetary Climate Model. This dataset will shortly be available publicly via the UK Centre for Environmental Data Analysis (https://www.ceda.ac.uk/). Here we present an overview of the new reanalysis and illustrate some of its results with the aim of alerting researchers to this new resource for future studies. PLR and AV acknowledge support from the UK Space Agency. The authors are grateful to Armin Kleinböhl and the MCS science team for advice and early access to retrieved data.Figure 1: Snapshots of the zonally averaged structure of the Martian atmosphere during northern late summer in Mars Year 28, showing the zonal and vertical wind, temperature and dust distribution, showing the presence of a persistent elevated dust layer from the MACDA2 reanalysis.[1] Montabone, L., Marsh, K., Lewis, S. R., Read, P. L., Smith, M. D., Holmes, J., et al. (2014). The Mars Analysis Correction Data Assimilation (MACDA) dataset V1.0. Geoscience Data Journal, 1(2), 129–139. https://doi.org/10.1002/gdj3.13[2] Greybush, S. J., Kalnay, E., Wilson, R. J., Hoffman, R. N., Nehrkorn, T., Leidner, M., et al. (2019). The Ensemble Mars Atmosphere Reanalysis System (EMARS) version 1.0. Geoscience Data Journal, 6(2), 137–150. https://doi.org/10.1002/gdj3.77[3] Holmes, J. A., Lewis, S. R., & Patel, M. R. (2020). OpenMARS: A global record of Martian weather from 1999 to 2015. Planetary and Space Science, 188, 104962. https://doi.org/10.1016/j.pss.2020.104962[4] Ruan, T., Young, R. M. B., Lewis, S. R., Montabone, L., Valeanu, A., & Read, P. L. (2021). Assimilation of both column- and layer-integrated dust opacity observations in the Martian atmosphere. Earth and Space Science, 8(12), e2021EA001869. https://doi.org/10.1029/2021EA001869[5] Young, R. M. B., Millour, E., Guerlet, S., Forget, F., Ignatiev, N., Grigoriev, A. V., et al. (2022). Assimilation of temperatures and column dust opacities measured by ExoMars TGO-ACS-TIRVIM during the MY34 Global Dust Storm. Journal of Geophysical Research: Planets, 127, e2022JE007312. https://doi.org/10.1029/2022JE007312[6] Young, R. M. B., Millour, E., Forget, F., Smith, M. D., Aljaberi, M., Edwards, C. S., et al. (2022). First assimilation of atmospheric temperatures from the Emirates Mars InfraRed Spectrometer. Geophysical Research Letters, 49, e2022GL099656.

Abstract:

Survive or not survive, that is the question of the 500-hour JWST Rocky Worlds DDT Program. Whether a terrestrial planets’ atmosphere can suffer under the intense XUV of its host, or if it completely escapes, these are the questions we explore. Zahnle & Catling (2017) defined the Cosmic Shoreline, but recent observations from JWST reveal airless worlds around M-stars, calling for a refinement of this “receding” shoreline (Pass et al. 2025). M-stars spend a longer time in pre-main sequence, subjecting their orbiting worlds to some higher intensity XUV activity. This complicates our present understanding of this shoreline. Investigating chemical effects of planet-star interactions could be the key to a more complete picture of this shoreline.  We investigate the interplay between photochemistry, mixing, and escape of carbon dioxide atmospheres under intense and mild XUV fluxes as follow on work to both Johnstone et al. (2018) and Nakayama et al. (2022). We expand on this work by adopting thermal structure models from Nakayama et al. (2022) and apply them to identify key chemical pathways for escape. We create a reduced C-O chemical network including neutral and ionic species to identify these pathways. As photochemistry simulations take into account many reactions, these 1D calculations are too computationally expensive to be done in 3D. Although rudimentary at best, the mixing parameter– eddy diffusion term, K_zz, comprises the dynamical element of 1D photochemical simulations. Here, we consider the mixing of photochemical products in competition with escape to explore the chemical pathways of retention and loss. We compare the photochemical model results for active and inactive cases for the Trappist-1 system planets. Then, using the resulting composition-dependent heating and cooling rates for Trappist-1 planets, we assess their propensity for efficient atomic line cooling versus escape. We follow the work of Chatterjee & Pierrehumbert (2024) in this assessment.  Finally, using our pathway analysis, we find an analytical formula for calculating an energy-limited escape boundary for these planets based on composition.  It is important here to note the limitations of 1D work. First, there exists an exchange of rigor between modelling chemistry and dynamics. Insights from this work are ripe for implementation into 3D GCMs, especially in response to incorporating UV-driven processes for thermospheric modelling mentioned in Ding and Wordsworth (2019). Second, interaction with the interior is important in the early phase of planetary formation, i.e., the magma ocean phase. Due to exchange between atmosphere and magma early in the planet’s formation, incorporation with an interior-atmosphere model would better constrain higher pressure chemical abundances. Although this work focuses on the upper atmosphere, extrapolation to the surface environment is a key goal for understanding a planet.  Considering planet-star interaction is imperative for the selection of targets for observation. However, it is also important when considering anomalous detections of atmospheres around planets predicted to not have an atmosphere. This could be a first step in determining an atmosphere as non-primary and/or distinguishing between an airless planet and one with high altitude haze.