2023ApJ...957...22Z

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.