Radiative-convective Equilibrium Research Articles

Numerical Diffusion and Turbulent Mixing in Convective Self‐Aggregation

AbstractConvective Self‐Aggregation (CSA) is a common feature of idealized numerical simulations of the tropical atmosphere in Radiative‐Convective Equilibrium (RCE). However, at coarse grid resolution where deep convection is not fully resolved, the occurrence of this phenomenon is extremely sensitive to subgrid‐scale processes. This study examines the role of mixing and entrainment, provided by the turbulence model and the implicit numerical diffusion. The study compares the results of two models, WRF and SAM, by varying turbulence models, initial conditions, and horizontal spatial resolution. At a coarse grid resolution of 3 km, the removal of turbulent mixing prevents CSA in models with low numerical diffusivity but is preserved in models with high numerical diffusivity. When the horizontal grid resolution is refined to 1 km, CSA can only be achieved by increasing explicit turbulent mixing, even with a small amount of shallow clouds. Therefore, the sensitivity of CSA to horizontal grid resolution is not primarily caused by the decrease in shallow clouds. The analysis of the total water path spectrum suggests that the amplitude of initial humidity perturbations introduced by convection in the free troposphere is the key factor. This amplitude is regulated by turbulent mixing and diffusion at small scales. Prior to the onset of CSA, increased mixing makes updrafts more sensitive to the dryness of the free troposphere, which strengthens the moisture‐convection feedback. This leads to an increased distance between convective cores and a stronger humidity perturbation in the free troposphere, which can destabilize the RCE state.

Large‐Scale Tropical Circulation Intensification by Aerosol Effects on Clouds

AbstractThis study addresses a critical gap in understanding anthropogenic influences on tropical climate dynamics by investigating the impact of aerosol‐cloud interactions on large‐scale circulation. Despite extensive research on greenhouse gas‐induced warming and its effects on tropical circulation, the impact of aerosols, particularly their interactions with clouds, on large‐scale circulation remains understudied. Utilizing large‐domain radiative convective equilibrium cloud‐resolving simulations, this research demonstrates that increasing aerosol concentration intensifies tropical overturning circulation, evaluated at the mid‐troposphere , strongly correlating with domain mean cloud and radiative properties. Employing a weak temperature gradient approximation, I attribute variations in to changes in clear‐sky radiative cooling rather than stability. These radiative cooling changes are linked to humidity changes driven by warm rain suppression by aerosols. This study's findings underscore the need to take into account microphysical changes, particularly aerosol concentrations, when studying anthropogenic effects on tropical circulation.

Spatial Dispersion and Statistical Description of Organized Cumulus Cloud Ensembles in Radiative Convective Equilibrium

AbstractDescriptions of cloud ensembles in radiative convective equilibrium (RCE) rely mostly on the assumption that clouds are randomly distributed in space, a hypothesis that obviously fails in presence of strong organization. In this work, idealized RCE simulations at horizontal grid spacings ranging from 2 km to 125 m are analyzed, displaying a transition between complete randomness at coarse resolution to a high degree of cold pool driven organization at high resolution. The objective is then twofold: (a) characterizing the spatial cloud patterns focusing on correlation scales and association with covariates; (b) providing a statistical description of organized cloud ensembles to generalize existing theories based on complete randomness. It is demonstrated that organized cloud ensembles may be treated as heterogeneous point processes where individual clouds do not interact substantially with their neighbors. In other words, clouds may be considered as randomly distributed if we restrict the space to regions of high cloud density. In addition, it is shown that in highly organized situations local cloud counts may be modeled using negative binomial distributions, while cloud mass fluxes follow lognormal distributions. The total mass flux within regions of varying sizes can then be predicted by compounding the two aforementioned distributions. The total mass flux variability at scales ≲5 km is dominated by the heavy‐tails of the mass flux distributions, whereas at scales ≳5 km, it is entirely controlled by the spatial cloud count variability (the spatial cloud pattern). These results have important implications for the parameterization of organized convection in quasi‐equilibrium.

The Impact of Rotation on Tropical Climate, the Hydrologic Cycle, and Climate Sensitivity

AbstractThis work explores the impact of rotation on tropical convection and climate. As our starting point, we use the RCEMIP experiments as control simulations and run additional simulations with rotation. Compared to radiative convective equilibrium (RCE) experiments, rotating RCE (RRCE) experiments have a more stable and humid atmosphere with higher precipitation rates. The intensity of the overturning circulation decreases, water vapor is cycled through the troposphere at a slower rate, the subsidence fraction decreases, and the climate sensitivity increases. Several of these changes can be attributed to an increased flux of latent and sensible heat that results from an increase of near‐surface wind speed with rotation shortly after model initialization. The increased climate sensitivity results from changes of both the longwave cloud radiative effect and the longwave clear‐sky radiative fluxes. This work demonstrates the sensitivity of atmospheric humidity and surface fluxes of moisture and temperature to rotation.

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Atmospheric properties of AF Lep b with forward modeling

Context. About a year ago, a super-Jovian planet was directly imaged around the nearby young solar-type star AF Lep. The 2.8 MJup planet orbiting at a semimajor axis of 8.2 au matches the predicted location based on the HIPPARCOS-Gaia astrometric acceleration. Aims. Our aim is to expand the atmospheric exploration of AF Lep b by modeling all available observations obtained with SPHERE at VLT (in the range 0.95–1.65, at 2.105, and at 2.253 µm), and NIRC2 at Keck (at 3.8 µm) with self-consistent atmospheric models. Methods. To understand the physical properties of this exoplanet, we used ForMoSA. This forward-modeling code compares observations with grids of pre-computed synthetic atmospheric spectra using Bayesian inference methods. We used Exo-REM, an atmospheric radiative-convective equilibrium model, including the effects of nonequilibrium processes and clouds. Results. From the atmospheric modeling we derived solutions at the low Teff of ~750 K. Our analysis also favors a metal-rich atmosphere (>0.4) and solar to super-solar carbon-to-oxygen ratio (~0.6). We tested the robustness of the estimated values for each parameter by cross-validating our models using the leave-one-out strategy, where all points are used iteratively as validation points. Our results indicate that the photometry point at 3.8 µm strongly drives the metal-rich and super-solar carbon-to-oxygen solutions. Conclusions. Our atmospheric forward-modeling analysis strongly supports the planetary nature of AF Lep b. Its spectral energy distribution is consistent with that of a young, cold, early-T super-Jovian planet. We recover physically consistent solutions for the surface gravity and radius, which allows us to reconcile atmospheric forward modeling with evolutionary models, in agreement with the previously published complementary analysis done by retrievals. Finally, we identified that future data at longer wavelengths are necessary before we can conclude about the metal-rich nature of AF Lep b.

The Sonora Substellar Atmosphere Models. IV. Elf Owl: Atmospheric Mixing and Chemical Disequilibrium with Varying Metallicity and C/O Ratios

Disequilibrium chemistry due to vertical mixing in the atmospheres of many brown dwarfs and giant exoplanets is well established. Atmosphere models for these objects typically parameterize mixing with the highly uncertain K zz diffusion parameter. The role of mixing in altering the abundances of C-N-O-bearing molecules has mostly been explored for atmospheres with a solar composition. However, atmospheric metallicity and the C/O ratio also impact atmospheric chemistry. Therefore, we present the Sonora Elf Owl grid of self-consistent cloud-free 1D radiative-convective equilibrium model atmospheres for JWST observations, which includes a variation in K zz across several orders of magnitude and also encompasses subsolar to supersolar metallicities and C/O ratios. We find that the impact of K zz on the T(P) profile and spectra is a strong function of both T eff and metallicity. For metal-poor objects, K zz has large impacts on the atmosphere at significantly higher T eff than in metal-rich atmospheres, where the impact of K zz is seen to occur at lower T eff. We identify significant spectral degeneracies between varying K zz and metallicity in multiple wavelength windows, in particular, at 3–5 μm. We use the Sonora Elf Owl atmospheric grid to fit the observed spectra of a sample of nine early to late T-type objects from T eff = 550–1150 K. We find evidence for very inefficient vertical mixing in these objects, with inferred K zz values lying in the range between ∼101 and 104 cm2 s−1. Using self-consistent models, we find that this slow vertical mixing is due to the observations, which probe mixing in the deep detached radiative zone in these atmospheres.

Dynamically coupled kinetic chemistry in brown dwarf atmospheres – II. Cloud and chemistry connections in directly imaged sub-Jupiter exoplanets

ABSTRACT With JWST slated to gain high-fidelity time-dependent data on brown dwarf atmospheres, it is highly anticipated to do the same for directly imaged, sub-Jupiter exoplanets. With this new capability, the need for a full three-dimensional (3D) understanding to explain spectral features and their time dependence is becoming a vital aspect for consideration. To examine the atmospheric properties of directly imaged sub-Jupiter exoplanets, we use the 3D Exo-FMS general circulation model to simulate a metal-enhanced generic young sub-Jupiter object. We couple Exo-FMS to a kinetic chemistry scheme, a tracer-based cloud formation scheme and a spectral radiative-transfer model to take into account the chemical and cloud feedback on the atmospheric thermochemical and dynamical properties. Our results show a highly complex feedback between clouds and chemistry on to the 3D temperature structure of the atmosphere, bringing about latitudinal differences and inducing time-dependent stormy features at photospheric pressures. This suggests a strong connection and feedback between the spatial cloud coverage and chemical composition of the atmosphere, with the temperature changes and dynamical motions induced by cloud opacity and triggered convection feedback driving chemical species behaviour. In addition, we also produce synthetic latitude-dependent and time-dependent spectra of our model to investigate atmospheric variability and periodicity in commonly used photometric bands. Overall, our efforts put the included physics in 3D simulations of exoplanets on par with contemporary 1D radiative-convective equilibrium modelling.

An Analytic Theory of Near-Surface Relative Humidity over Land

Abstract There is no simple explanation for the spatial structure of near-surface relative humidity over land. We present a diagnostic theory for zonally and temporally averaged near-surface relative humidity (RH) over land based on energy budgets of an atmospheric column in radiative–convective equilibrium. The theory analytically relates RH to the surface evaporative fraction (EF), has no calibrated parameters, and is quantitatively accurate when compared with RH from a reanalysis, and with cloud-permitting simulations over an idealized land surface. The theory is used to answer two basic questions. First, why is RH never especially low (e.g., 1%)? The theory shows that established lower bounds on EF over land and ocean are equivalent to lower bounds on RH that preclude particularly low values, at least for conditions typical of the modern Earth. Second, why is the latitudinal profile of RH over land shaped like the letter W, when both specific humidity and saturation specific humidity essentially decline monotonically from the equator to the poles? The theory predicts that the latitudinal profile of RH should look more like that of water stored in the soil (which also exhibits a W-shaped profile) than in the air (which does not).

Detecting Cold Pool Family Trees in Convection Resolving Simulations

AbstractRecent observations and modeling increasingly reveal the key role of cold pools in organizing the convective cloud field. Several methods for detecting cold pools in simulations exist, but are usually based on buoyancy fields and fall short of reliably identifying the active gust front. The current cold pool (CP) detection and tracking algorithm (CoolDeTA), aims to identify cold pools and follow them in time, thereby distinguishing their active gust fronts and the “offspring” rain cells generated nearby. To accomplish these tasks, CoolDeTA utilizes a combination of thermodynamic and dynamical variables and examines the spatial and temporal relationships between cold pools and rain events. We demonstrate that CoolDeTA can reconstruct CP family trees. Using CoolDeTA we can contrast radiative convective equilibrium (RCE) and diurnal cycle CP dynamics, as well as cases with vertical wind shear and without. We show that the results obtained are consistent with a conceptual model where CP triggering of children rain cells follows a simple birth rate, proportional to a CP's gust front length. The proportionality factor depends on the ambient atmospheric stability and is lower for RCE, in line with marginal stability as traditionally ascribed to the moist adiabat. In the diurnal case, where ambient stability is lower, the birth rate thus becomes substantially higher, in line with periodic insolation forcing—resulting in essentially run‐away mesoscale excitations generated by a single parent rain cell and its CP.

Dependence of the Radiative-Convective Equilibrium Structure of the Lower Atmosphere of Venus on the Thermodynamic Model

Dependence of the Radiative-Convective Equilibrium Structure of the Lower Atmosphere of Venus on the Thermodynamic Model

Tropical mountain ice core δ18O: A Goldilocks indicator for global temperature change.

Very high tropical alpine ice cores provide a distinct paleoclimate record for climate changes in the middle and upper troposphere. However, the climatic interpretation of a key proxy, the stable water oxygen isotopic ratio in ice cores (δ18Oice), remains an outstanding problem. Here, combining proxy records with climate models, modern satellite measurements, and radiative-convective equilibrium theory, we show that the tropical δ18Oice is an indicator of the temperature of the middle and upper troposphere, with a glacial cooling of -7.35° ± 1.1°C (66% CI). Moreover, it severs as a "Goldilocks-type" indicator of global mean surface temperature change, providing the first estimate of glacial stage cooling that is independent of marine proxies as -5.9° ± 1.2°C.Combined with all estimations available gives the maximum likelihood estimate of glacial cooling as -5.85° ± 0.51°C.

The Inhomogeneity Effect. I. Inhomogeneous Surface and Atmosphere Accelerate Planetary Cooling

We propose a general principle that under radiative-convective equilibrium, the spatial and temporal variations in a planet’s surface and atmosphere tend to increase its cooling. This principle is based on Jensen’s inequality and the curvature of the response functions of surface temperature and outgoing cooling flux to changes in incoming stellar flux and atmospheric opacity. We use an analytical model to demonstrate that this principle holds for various planet types: (1) on an airless planet, the mean surface temperature is lower than its equilibrium temperature; (2) on terrestrial planets with atmospheres, the inhomogeneity of incoming stellar flux and atmospheric opacity reduces the mean surface temperature; (3) on giant planets, inhomogeneously distributed stellar flux and atmospheric opacity increase the outgoing infrared flux, cooling the interior. Although the inhomogeneity of visible opacity might sometimes heat the atmosphere, the effect is generally much smaller than the inhomogeneous cooling effect of infrared opacity. Compared with the homogeneous case, the mean surface temperature on inhomogeneous terrestrial planets can decrease by more than 20%, and the internal heat flux on giant planets can increase by over an order of magnitude. Despite simplifications in our analytical framework, the effect of stellar flux inhomogeneity appears to be robust, while further research is needed to fully understand the effects of opacity inhomogeneity in more realistic situations. This principle impacts our understanding of planetary habitability and the evolution of giant planets using low-resolution and one-dimensional frameworks that may have previously overlooked the role of inhomogeneity.

The Inhomogeneity Effect. III. Weather Impacts on the Heat Flow of Hot Jupiters

The interior flux of a giant planet impacts atmospheric motion, and the atmosphere dictates the interior’s cooling. Here we use a non-hydrostatic general circulation model (Simulating Non-hydrostatic Atmospheres on Planets) coupled with a multi-stream multi-scattering radiative module (High-performance Atmospheric Radiation Package) to simulate the weather impacts on the heat flow of hot Jupiters. We found that the vertical heat flux is primarily transported by convection in the lower atmosphere and regulated by dynamics and radiation in the overlying radiation-circulation zone. The temperature inversion occurs on the dayside and reduces the upward radiative flux. The atmospheric dynamics relay the vertical heat transport until the radiation becomes efficient in the upper atmosphere. The cooling flux increases with atmospheric drag due to increased day–night contrast and spatial inhomogeneity. The temperature dependence of the infrared opacity greatly amplifies the opacity inhomogeneity. Although atmospheric circulation could transport heat downward in a narrow region above the radiative-convective boundary, the opacity inhomogeneity effect overcomes the dynamical effect and leads to a larger overall interior cooling than the local simulations with the same interior entropy and stellar flux. The enhancement depends critically on the equilibrium temperature, drag, and atmospheric opacity. In a strong-drag atmosphere hotter than 1600 K, a significant inhomogeneity effect in three-dimensional (3D) models can boost interior cooling several-fold compared to the 1D radiative-convective equilibrium models. This study confirms the analytical argument of the inhomogeneity effect in the companion papers by Zhang. It highlights the importance of using 3D atmospheric models in understanding the inflation mechanisms of hot Jupiters and giant planet evolution in general.

ELemental abundances of Planets and brown dwarfs Imaged around Stars (ELPIS). I. Potential Metal Enrichment of the Exoplanet AF Lep b and a Novel Retrieval Approach for Cloudy Self-luminous Atmospheres

AF Lep A+b is a remarkable planetary system hosting a gas-giant planet that has the lowest dynamical mass among directly imaged exoplanets. We present an in-depth analysis of the atmospheric composition of the star and planet to probe the planet’s formation pathway. Based on new high-resolution spectroscopy of AF Lep A, we measure a uniform set of stellar parameters and elemental abundances (e.g., [Fe/H] = −0.27 ± 0.31 dex). The planet’s dynamical mass ( M Jup) and orbit are also refined using published radial velocities, relative astrometry, and absolute astrometry. We use petitRADTRANS to perform chemically consistent atmospheric retrievals for AF Lep b. The radiative–convective equilibrium temperature profiles are incorporated as parameterized priors on the planet’s thermal structure, leading to a robust characterization for cloudy self-luminous atmospheres. This novel approach is enabled by constraining the temperature–pressure profiles via the temperature gradient , a departure from previous studies that solely modeled the temperature. Through multiple retrievals performed on different portions of the 0.9–4.2 μm spectrophotometry, along with different priors on the planet’s mass and radius, we infer that AF Lep b likely possesses a metal-enriched atmosphere ([Fe/H] > 1.0 dex). AF Lep b’s potential metal enrichment may be due to planetesimal accretion, giant impacts, and/or core erosion. The first process coincides with the debris disk in the system, which could be dynamically excited by AF Lep b and lead to planetesimal bombardment. Our analysis also determines T eff ≈ 800 K, dex, and the presence of silicate clouds and disequilibrium chemistry in the atmosphere. Straddling the L/T transition, AF Lep b is thus far the coldest exoplanet with suggested evidence of silicate clouds.

Extreme Precipitation in Tropical Squall Lines

AbstractSquall lines are substantially influenced by the interaction of low‐level shear with cold pools associated with convective downdrafts. Beyond an optimal shear amplitude, squall lines tend to orient themselves at an angle with respect to the low‐level shear. While the mechanisms behind squall line orientation seem to be increasingly well understood, uncertainties remain on the implications of this orientation. Roca and Fiolleau (2020, https://doi.org/10.1038/s43247-020-00015-4) show that long lived mesoscale convective systems, including squall lines, are disproportionately involved in rainfall extremes in the tropics. This article investigates the influence of the interaction between low‐level shear and squall line outflow on squall line generated precipitation extrema in the tropics. Using a cloud resolving model, simulated squall lines in radiative convective equilibrium amid a shear‐dominated regime (super optimal), a balanced regime (optimal), and an outflow dominated regime (suboptimal). Our results show that precipitation extremes in squall lines are 40% more intense in the case of optimal shear and remain 30% superior in the superoptimal regime relative to a disorganized case. With a theoretical scaling of precipitation extremes (C. Muller & Takayabu, 2020, https://doi.org/10.1088/1748-9326/ab7130), we show that the condensation rates control the amplification of precipitation extremes in tropical squall lines, mainly due to its change in vertical mass flux (dynamic component). The reduction of dilution by entrainment explains half of this change, consistent with Mulholland et al. (2021, https://doi.org/10.1175/jas-d-20-0299.1). The other half is explained by increased cloud‐base velocity intensity in optimal and superoptimal squall lines.

The Quasi-Linear Relation between Planetary Outgoing Longwave Radiation and Surface Temperature: A Climate Footprint of Radiative and Nonradiative Processes

Abstract The slope of the quasi-linear relation between planetary outgoing longwave radiation (OLR) and surface temperature (TS) is an important parameter measuring the sensitivity of Earth’s climate system. The primary objective of this study is to seek a general explanation for the quasi-linear OLR–TS relation that remains valid regardless of the strength of the atmospheric window’s narrowing effect on planetary thermal emission at higher temperatures. The physical understanding of the quasi-linear OLR–TS relation and its slope is gained from observation analysis, climate simulations with radiative–convective equilibrium and general circulation models, and a series of online feedback suppression experiments. The observed quasi-linear OLR–TS relation manifests a climate footprint of radiative (such as the greenhouse effect) and nonradiative processes (poleward energy transport). The former acts to increase the meridional gradient of surface temperature and the latter decreases the meridional gradient of atmospheric temperatures, causing the flattening of the meridional profile of the OLR. Radiative processes alone can lead to a quasi-linear OLR–TS relation that is more steeply sloped. The atmospheric poleward energy transport alone can also lead to a quasi-linear OLR–TS relation by rerouting part of the OLR to be emitted from a warmer place to a colder place. The combined effects of radiative and nonradiative processes make the quasi-linear OLR–TS relation less sloped with a higher degree of linearity. In response to anthropogenic radiative forcing, the slope of the quasi-linear OLR–TS relation is further reduced via stronger water vapor feedback and enhanced poleward energy transport. Significance Statement The slope of the quasi-linear relation between planetary outgoing longwave radiation (OLR) and surface temperature (TS) is an important parameter measuring the sensitivity of Earth’s climate system. The observed quasi-linear OLR–TS relation manifests a climate footprint of radiative (greenhouse effect) and nonradiative processes (poleward energy transport). Radiative processes alone can lead to a quasi-linear OLR–TS relation that is more steeply sloped. The atmospheric poleward energy transport alone can also lead to a quasi-linear OLR–TS relation by rerouting part of the OLR to be emitted from a warmer place to a colder place. The combined effects of radiative and nonradiative processes make the quasi-linear OLR–TS relation less sloped with a higher degree of linearity.

Climate sensitivity from radiative-convective equilibrium: A chalkboard approach

Simple models for Earth's climate sensitivity (i.e. its temperature response to radiative forcing) are developed by combining the time-tested idealization of one-dimensional radiative-convective equilibrium (RCE) with simple yet quantitatively reasonable models for CO2 forcing and the water vapor feedback. Along the way, we introduce key paradigms including the emission level approximation, the forcing-feedback decomposition of climate sensitivity, and “Simpson's law” for water vapor thermal emission. We also discuss climate feedbacks unaccounted for in this RCE framework, as well as differing variants of climate sensitivity, all of which may be ripe for their own chalkboard treatments.

Interior–atmosphere modelling to assess the observability of rocky planets with JWST

Context. Super-Earths present compositions dominated by refractory materials. However, there is a degeneracy in their interior structure between a planet that has no atmosphere and a small Fe content, and a planet that has a thin atmosphere and a higher core mass fraction. To break this degeneracy, atmospheric characterisation observations are required. Aims. We present a self-consistent interior–atmosphere model to constrain the volatile mass fraction, surface pressure, and temperature of rocky planets with water and CO2 atmospheres. The parameters obtained in our analysis can be used to predict observations in emission spectroscopy and photometry with JWST, which can determine the presence of an atmosphere and, if present, its composition. Methods. We coupled a 1D interior model with a supercritical water layer to an atmospheric model. In order to obtain the bolometric emission and Bond albedo for an atmosphere in radiative-convective equilibrium, we used a low-resolution k-correlated atmospheric model. We generated emission spectra with the same atmospheric model at a higher resolution (R = 200–300). An adaptive Markov chain Monte Carlo was employed for an efficient sampling of the parameter space at low volatile mass fractions. Results. From our interior structure retrieval, we conclude that TRAPPIST-1 c most likely has a bare surface, although the presence of an atmosphere cannot be ruled out. We estimate a 1σ confidence interval of the surface pressure for a water-dominated atmosphere of Psurf = 40 ± 40 bar. We generated spectra for these two scenarios to compare with the emission flux of TRAPPIST-1 c recently observed in the MIRI F1500W filter. This is compatible with bare rock surfaces or a thin atmosphere with little or no CO2. In the case of 55 Cancri e, a combined spectrum with NIRCam and MIRI LRS may present high uncertainties at wavelengths between 3 and 3.7 µm. However, this does not affect the identification of H2O because it does not present spectral features in this wavelength range.

The emergence of a new wintertime Arctic energy balance regime

The modern Arctic climate during wintertime is characterized by sea-ice cover, a strong surface temperature inversion, and the absence of convection. Correspondingly, the energy balance in the Arctic atmosphere today is dominated by atmospheric radiative cooling and advective heating, so-called radiative advective equilibrium. Climate change in the Arctic involves sea-ice melt, vanishing of the surface inversion, and emergence of convective precipitation. Here we show climate change in the Arctic involves the emergence of a new energy balance regime characterized by radiative cooling, convective heating, and advective heating, so-called radiative convective advective equilibrium. A time-dependent decomposition of the atmospheric energy balance shows the regime transition is associated with enhanced radiative cooling followed by decreased advective heating. The radiative cooling response consists of a robust clear-sky greenhouse effect and a transient cloud contribution that varies across models. Mechanism-denial experiments in an aquaplanet with and without interactive sea ice highlight the important role of sea-ice melt in both the radiative cooling and advective heating responses. The results show that climate change in the Arctic involves temporally evolving mechanisms, suggesting that an emergent constraint based on historical data or trends may not constrain the long-term response.