No Arabic abstract
Using an idealised general circulation model, we investigate the atmospheric circulation of Earth-like terrestrial planets in a variety of orbital configurations. We relax the common assumption of the planet being tidally-locked, and look at the role atmospheric dynamics can have in the observed thermal phase curve when the substellar point is non-stationary. In slowly rotating planets, a moving forcing can induce strong jets in the upper troposphere, both prograde and retrograde, sensitive to the speed and direction of the diurnal forcing. We find that, consistent with previous shallow water model experiments, the thermal phase curve offset is sensitive to the velocity of the substellar point moving across the surface of the planet. For a planet with a known orbital period, the results show that the observed hotspot on the planet could be either east or west of the substellar point, depending on whether the planet is tidally-locked or not.
Using a shallow water model with time-dependent forcing we show that the peak of an exoplanet thermal phase curve is, in general, offset from secondary eclipse when the planet is rotating. That is, the planetary hot-spot is offset from the point of maximal heating (the substellar point) and may lead or lag the forcing; the extent and sign of the offset is a function of both the rotation rate and orbital period of the planet. We also find that the system reaches a steady-state in the reference frame of the moving forcing. The model is an extension of the well studied Matsuno-Gill model into a full spherical geometry and with a planetary-scale translating forcing representing the insolation received on an exoplanet from a host star. The speed of the gravity waves in the model is shown to be a key metric in evaluating the phase curve offset. If the velocity of the substellar point (relative to the planets surface) exceeds that of the gravity waves then the hotspot will lag the substellar point, as might be expected by consideration of forced gravity wave dynamics. However, when the substellar point is moving slower than the internal wavespeed of the system the hottest point can lead the passage of the forcing. We provide an interpretation of this result by consideration of the Rossby and Kelvin wave dynamics as well as, in the very slowly rotating case, a one-dimensional model that yields an analytic solution. Finally, we consider the inverse problem of constraining planetary rotation rate from an observed phase curve.
Tidally locked exoplanets likely host global atmospheric circulations with a superrotating equatorial jet, planetary-scale stationary waves and thermally-driven overturning circulation. In this work, we show that each of these features can be separated from the total circulation by using a Helmholtz decomposition, which splits the circulation into rotational (divergence free) and divergent (vorticity free) components. This technique is applied to the simulated circulation of a terrestrial planet and a gaseous hot Jupiter. For both planets, the rotational component comprises the equatorial jet and stationary waves, and the divergent component contains the overturning circulation. Separating out each component allows us to evaluate their spatial structure and relative contribution to the total flow. In contrast with previous work, we show that divergent velocities are not negligible when compared with rotational velocities, and that divergent, overturning circulation takes the form of a single, roughly isotropic cell that ascends on the day-side and descends on the night-side. These conclusions are drawn for both the terrestrial case and the hot Jupiter. To illustrate the utility of the Helmholtz decomposition for studying atmospheric processes, we compute the contribution of each of the circulation components to heat transport from day- to night-side. Surprisingly, we find that the divergent circulation dominates day-night heat transport in the terrestrial case and accounts for around half of the heat transport for the hot Jupiter. The relative contributions of the rotational and divergent components to day-night heat transport are likely sensitive to multiple planetary parameters and atmospheric processes, and merit further study.
Over large timescales, a terrestrial planet may be driven towards spin-orbit synchronous rotation by tidal forces. In this particular configuration, the planet exhibits permanent dayside and nightside, which may induce strong day-night temperature gradients. The nightside temperature depends on the efficiency of the day-night heat redistribution and determines the stability of the atmosphere against collapse. To better constrain the atmospheric stability, climate, and surface conditions of rocky planets located in the habitable zone of their host star, it is thus crucial to understand the complex mechanism of heat redistribution. Building on early works and assuming dry thermodynamics, we developed a hierarchy of analytic models taking into account the coupling between radiative transfer, dayside convection, and large-scale atmospheric circulation in the case of slowly rotating planets. There are two types of these models: a zero-dimensional two-layer approach and a two-column radiative-convective-subsiding-upwelling (RCSU) model. They yield analytical solutions and scaling laws characterising the dependence of the collapse pressure on physical features, which are compared to the results obtained by early works using 3D global climate models (GCMs). The analytical theory captures (i) the dependence of temperatures on atmospheric opacities and scattering in the shortwave and in the longwave, (ii) the behaviour of the collapse pressure observed in GCM simulations at low stellar fluxes that are due to the non-linear dependence of the atmospheric opacity on the longwave optical depth at the planets surface, (iii) the increase of stability generated by dayside sensible heating, and (iv) the decrease of stability induced by the increase of the planet size.
Using a 3D general circulation model (GCM), we investigate the sensitivity of the climate of tidally-locked Earth-like exoplanets, Trappist-1e and Proxima Centauri b, to the choice of a convection parameterization. Compared to a mass-flux convection parameterization, a simplified convection adjustment parameterization leads to a $>$60% decrease of the cloud albedo, increasing the mean day-side temperature by $approx$10 K. The representation of convection also affects the atmospheric conditions of the night side, via a change in planetary-scale wave patterns. As a result, using the convection adjustment scheme makes the night-side cold traps warmer by 17-36 K for the planets in our simulations. The day-night thermal contrast is sensitive to the representation of convection in 3D GCM simulations, so caution should be taken when interpreting emission phase curves. The choice of convection treatment, however, does not alter the simulated climate enough to result in a departure from habitable conditions, at least for the atmospheric composition and planetary parameters used in our study. The near-surface conditions both in the Trappist-1e and Proxima b cases remain temperate, allowing for an active water cycle. We further advance our analysis using high-resolution model experiments, in which atmospheric convection is simulated explicitly. Our results suggest that in a hypothetical global convection-permitting simulation the surface temperature contrast would be higher than in the coarse-resolution simulations with parameterized convection. In other words, models with parameterized convection may overestimate the inter-hemispheric heat redistribution efficiency.
A large fraction of known terrestrial-size exoplanets located in the Habitable Zone of M-dwarfs are expected to be tidally-locked. Numerous efforts have been conducted to study the climate of such planets, using in particular 3-D Global Climate Models (GCM). One of the biggest challenges in simulating such an extreme environment is to properly represent the effects of sub-grid convection. Most GCMs use either a simplistic convective-adjustment parametrization or sophisticated (e.g., mass flux scheme) Earth-tuned parametrizations. One way to improve the representation of convection is to study convection using Convection Resolving numerical Models (CRMs), with an fine spatial resolution . In this study, we developed a CRM coupling the non-hydrostatic dynamical core WRF with the radiative transfer and cloud/precipitation models of the LMD-Generic climate model to study convection and clouds on tidally-locked planets, with a focus on Proxima b. Simulations were performed for a set of 3 surface temperatures (corresponding to three different incident fluxes) and 2 rotation rates, assuming an Earth-like atmosphere. The main result of our study is that while we recover the prediction of GCMs that (low-altitude) cloud albedo increases with increasing stellar flux, the cloud feedback is much weaker due to transient aggregation of convection leading to low partial cloud cover.