Planetary Climate Dynamics

Abstract:

Brown dwarfs, planetary-mass objects and directly imaged giant planets exhibit significant observational evidence for active atmospheric circulation, raising critical questions about mechanisms driving the circulation, its fundamental nature and time variability. Our previous work has demonstrated the crucial role of cloud radiative feedback on driving a vigorous atmospheric circulation using local models that assume a Cartesian geometry and constant Coriolis parameters. In this study, we extend the models to a global geometry and explore properties of the global dynamics. We show that, under relatively strong dissipation in the bottom layers of the model, horizontally isotropic vortices are prevalent at mid-to-high latitudes while large-scale zonally propagating waves are dominant at low latitudes near the observable layers. The equatorial waves have both eastward and westward phase speeds, and the eastward components with typical velocities of a few hundred  m s−1 usually dominate the equatorial time variability. Lightcurves of the global simulations show variability with amplitudes from 0.5 per cent to a few percent depending on the rotation period and viewing angle. The time evolution of simulated lightcurves is critically affected by the equatorial waves, showing wave beating effects and differences in the lightcurve periodicity to the intrinsic rotation period. The vertical extent of clouds is the largest at the equator and decreases poleward due to the increasing influence of rotation with increasing latitude. Under weaker dissipation in the bottom layers, strong and broad zonal jets develop and modify wave propagation and lightcurve variability. Our modelling results help to qualitatively explain several features of observations of brown dwarfs and directly imaged giant planets, including puzzling time evolution of lightcurves, a slightly shorter period of variability in IR than in radio wavelengths, and the viewing angle dependence of variability amplitude and IR colors.

Abstract:

Observations of brown dwarfs (BDs), free-floating planetary-mass objects, and directly imaged extrasolar giant planets (EGPs) exhibit rich evidence of large-scale weather. Cloud radiative feedback has been proposed as a potential mechanism driving the vigorous atmospheric circulation on BDs and directly imaged EGPs, and yet it has not been demonstrated in three-dimensional dynamical models at relevant conditions. Here, we present a series of atmospheric circulation models that self-consistently couple dynamics with idealized cloud formation and its radiative effects. We demonstrate that vigorous atmospheric circulation can be triggered and self-maintained by cloud radiative feedback. Typical isobaric temperature variation could reach over 100 K and horizontally averaged wind speed could be several hundreds of $\, {\rm m\, s^{-1}}$. The circulation is dominated by cloud-forming and clear-sky vortices that evolve over time-scales from several to tens of hours. The typical horizontal length-scale of dominant vortices is closed to the Rossby deformation radius, showing a linear dependence on the inverse of rotation rate. Stronger rotation tends to weaken vertical transport of vapour and clouds, leading to overall thinner clouds. Domain-mean outgoing radiative flux exhibits variability over time-scales of tens of hours due to the statistical evolution of storms. Different bottom boundary conditions in the models could lead to qualitatively different circulation near the observable layer. The circulation driven by cloud radiative feedback represents a robust mechanism generating significant surface inhomogeneity as well as irregular flux time variability. Our results have important implications for near-infrared (IR) colours of dusty BDs and EGPs, including the scatter in the near-IR colour–magnitude diagram and the viewing-geometry-dependent near-IR colours.

Abstract:

Geochemical and astronomical evidence demonstrates that planet formation occurred in two spatially and temporally separated reservoirs. The origin of this dichotomy is unknown. We use numerical models to investigate how the evolution of the solar protoplanetary disk influenced the timing of protoplanet formation and their internal evolution. Migration of the water snow line can generate two distinct bursts of planetesimal formation that sample different source regions. These reservoirs evolve in divergent geophysical modes and develop distinct volatile contents, consistent with constraints from accretion chronology, thermochemistry, and the mass divergence of inner and outer Solar System. Our simulations suggest that the compositional fractionation and isotopic dichotomy of the Solar System was initiated by the interplay between disk dynamics, heterogeneous accretion, and internal evolution of forming protoplanets.

Abstract:

The data comprise two 2-dimensional gridded maps of horizontal wind measurements covering the north and south polar regions of Saturn, as previously published by Antuñano et al. (2015). As fully described in that paper, these measurements were derived from sets of Cassini Orbiter Imaging Sub-System (ISS) Wide Angle Camera (WAC) and Narrow Angle Camera (NAC) images using the continuum band CB2 and CB3 filters, acquired for the northern hemisphere in June 2013 and for the southern hemisphere using WAC CB2 and CB3 images taken in October 2006 and December 2008. Additional NAC images using the CB2 and red filters taken in July 2008 were also used to analyse the southern polar vortex. The WAC images covered a region extending from a planetocentric latitude of around 60-65 degrees to each pole (apart from a segment in longitude between around 35 - 110 degrees W) with a horizontal resolution equivalent to around 0.05 degrees latitude (around 50km) per pixel, while NAC images were mostly used for the polar vortices, with a resolution equivalent to around 0.01 degrees latitude (around 10 km) per pixel. Horizontal velocities were obtained using semi-automated image correlation methods between pairs of images separated in time by intervals of approximately 1-10 hours. The correlation algorithm used pixel box sizes of 23 x 23 (in the north) or 25 x 25 (in the south), leading to a spatial resolution of the velocity vectors equivalent to around 1 degree latitude or 1000 km outside the polar vortices, reducing to around 0.2 degrees or 200 km within the polar vortices themselves. The automatically generated velocity vectors were supplemented by a small number (around 1% of the total) of vectors obtained manually from the motion of visually identified cloud tracers. The estimated measurement uncertainty on each vector was around 5-10 m/s. The original velocity vectors from Antuñano et al. (2015) were interpolated onto a regular latitude-longitude grid using convex hulls and Delauney triangulation via the QHULL routine of the Interactive Data Language (IDL). The final datasets are held on a regular grid separated by 3-4 degrees in longitude and 0.23 degrees in latitude. Data are stored as two text files, tabulating the latitude and (west) longitude of each point and the eastward and northward velocity components respectively in units of m/s. Reference: Antuñano,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