The Zombie Vortex Instability (ZVI) occurs in the dead zones of protoplanetary disks (PPDs) where perturbations excite baroclinic critical layers, generating "zombie" vortices and turbulence. In this work, we investigate ZVI with nonuniform vertical stratification; while ZVI is triggered in the stratified regions away from the midplane, the subsequent turbulence propagates into and fills the midplane. ZVI turbulence alters the background Keplerian shear flow, creating a steady-state zonal flow. Intermittency is observed, where the flow cycles through near-laminar phases of zonal flow punctuated by chaotic bursts of new vortices. ZVI persists in the presence of radiative damping, as long as the thermal relaxation timescale is more than a few orbital periods. We refute the premature claim by Lesur & Latter (2016) that radiative damping inhibits ZVI for disk radii r>0.3 au. Their conclusions were based on unrealistically short cooling times using opacities with virtually no grain growth. We explore different grain growth and vertical settling scenarios, and find that the gas and dust in off-midplane regions are not necessarily in local thermodynamic equilibrium (LTE) with each other. In such cases, thermal relaxation timescales can be orders of magnitude longer than the optically thin cooling times assuming LTE because of the finite time for energy to be exchanged between gas and dust grains via collisions. With minimal amounts of grain growth and dust settling, the off-midplane regions of disks are susceptible to ZVI and much of the planet-forming regions can be filled with zombie vortices and turbulence.

In Zombie Vortex Instability (ZVI), perturbations excite critical layers in stratified, rotating shear flow (as in protoplanetary disks), causing them to generate vortex layers, which roll-up into anticyclonic zombie vortices and cyclonic vortex sheets. The process is self-sustaining as zombie vortices perturb new critical layers, spawning a next generation of zombie vortices. Here, we focus on two issues: the minimum threshold of perturbations that trigger self-sustaining vortex generation, and the properties of the late-time zombie turbulence on large and small scales. The critical parameter that determines whether ZVI is triggered is the magnitude of the vorticity on the small scales (and not velocity), the minimum Rossby number needed for instability is $Ro_{crit}\sim0.2$ for $\beta\equiv N/\Omega = 2$, where $N$ is the Brunt-V\"ais\"al\"a frequency. While the threshold is set by vorticity, it is useful to infer a criterion on the Mach number, for Kolmogorov noise, the critical Mach number scales with Reynolds number: $Ma_{crit}\sim Ro_{crit}Re^{-1/2}$. In protoplanetary disks, this is $Ma_{crit}\sim10^{-6}$. On large scales, zombie turbulence is characterized by anticyclones and cyclonic sheets with typical Rossby number $\sim$0.3. The spacing of the cyclonic sheets and anticyclones appears to have a "memory" of the spacing of the critical layers. On the small scales, zombie turbulence has no memory of the initial conditions and has a Kolmogorov-like energy spectrum. While our earlier work was in the limit of uniform stratification, we have demonstrated that ZVI works for non-uniform Brunt-V\"ais\"al\"a frequency profiles that may be found in protoplanetary disks.

We imaged Uranus in the near infrared from 2012 into 2014, using the Keck/NIRC2 camera and Gemini/NIRI camera, both with adaptive optics. We obtained exceptional signal to noise ratios by averaging 8-16 individual exposures in a planet-fixed coordinate system. noise-reduced images revealed many low-contrast discrete features and large scale cloud patterns not seen before, including scalloped waveforms just south of the equator. In all three years numerous small (600-700 km wide) and mainly bright discrete features were seen within the north polar region (north of about 55\deg N). Over 850 wind measurements were made, the vast majority of which were in the northern hemisphere. These revealed an extended region of solid body rotation between 62\deg N and at least 83\deg N, at a rate of 4.08$\pm0.015$\deg/h westward relative to the planet's interior (radio) rotation of 20.88\deg/h westward. Near-equatorial speeds measured with high accuracy give different results for waves and small discrete features, with eastward drift rates of 0.4\deg/h and 0.1\deg/h respectively. The region of polar solid body rotation is a close match to the region of small-scale polar cloud features, suggesting a dynamical relationship. While winds at high southern latitudes (50\deg S - 90\deg S) are unconstrained by groundbased observations, a recent reanalysis of 1986 Voyager 2 observations by Karkoschka (2015, Icarus 250, 294-307) has revealed an extremely large north-south asymmetry in this region, which might be seasonal. Greatly increased activity was seen in 2014, including the brightest ever feature seen in K' images (de Pater et al. 2015, Icarus 252, 121-128). Over the 2012-2014 period we identified six persistent discrete features. Three were tracked for more than two years, two more for more than one year, and one for at least 5 months and continuing.

There is considerable interest in hydrodynamic instabilities in dead zones of protoplanetary disks as a mechanism for driving angular momentum transport and as a source of particle-trapping vortices to mix chondrules and incubate planetesimal formation. We present simulations with a pseudo-spectral anelastic code and with the compressible code Athena, showing that stably stratified flows in a shearing, rotating box are violently unstable and produce space-filling, sustained turbulence dominated by large vortices with Rossby numbers of order 0.2-0.3. This Zombie Vortex Instability (ZVI) is observed in both codes and is triggered by Kolmogorov turbulence with Mach numbers less than 0.01. It is a common view that if a given constant density flow is stable, then stable vertical stratification should make the flow even more stable. Yet, we show that sufficient vertical stratification can be unstable to ZVI. ZVI is robust and requires no special tuning of boundary conditions, or initial radial entropy or vortensity gradients (though we have studied ZVI only in the limit of infinite cooling time). The resolution of this paradox is that stable stratification allows for a new avenue to instability: baroclinic critical layers. ZVI has not been seen in previous studies of flows in rotating, shearing boxes because those calculations frequently lacked vertical density stratification and/or sufficient numerical resolution. Although we do not expect appreciable angular momentum transport from ZVI in the small domains in this study, we hypothesize that ZVI in larger domains with compressible equations may lead to angular transport via spiral density waves.

A previously unknown instability creates space-filling lattices of 3D vortices in linearly-stable, rotating, stratified shear flows. The instability starts from an easily-excited critical layer. The layer intensifies by drawing energy from the background shear and rolls-up into vortices that excite new critical layers and vortices. The vortices self-similarly replicate to create lattices of turbulent vortices. The vortices persist for all time. This self-replication occurs in stratified Couette flows and in the dead zones of protoplanetary disks where it can de-stabilize Keplerian flows.

We have developed a three-dimensional (3D) spectral hydrodynamic code to study vortex dynamics in rotating, shearing, stratified systems (e.g. the atmosphere of gas giant planets, protoplanetary disks around newly forming protostars). The time-independent background state is stably stratified in the vertical direction and has a unidirectional linear shear flow aligned with one horizontal axis. Superposed on this background state is an unsteady, subsonic flow that is evolved with the Euler equations subject to the anelastic approximation to filter acoustic phenomena. A Fourier-Fourier basis in a set of quasi-Lagrangian coordinates that advect with the background shear is used for spectral expansions in the two horizontal directions. For the vertical direction, two different sets of basis functions have been implemented: (1) Chebyshev polynomials on a truncated, finite domain, and (2) rational Chebyshev functions on an infinite domain. Use of this latter set is equivalent to transforming the infinite domain to a finite one with a cotangent mapping, and using cosine and sine expansions in the mapped coordinate. The nonlinear advection terms are time integrated explicitly, whereas the Coriolis force, buoyancy terms, and pressure/enthalpy gradient are integrated semi-implicitly. We show that internal gravity waves can be damped by adding new terms to the Euler equations. The code exhibits excellent parallel performance with the Message Passing Interface (MPI). As a demonstration of the code, we simulate the merger of two 3D vortices in the midplane of a protoplanetary disk.

We present the results of high-resolution, three-dimensional (3D) hydrodynamic simulations of the dynamics and formation of coherent, long-lived vortices in stably-stratified protoplanetary disks. Tall, columnar vortices that extend vertically through many scale heights in the disk are unstable to small perturbations; such vortices cannot maintain vertical alignment over more than a couple scale heights and are ripped apart by the Keplerian shear. Short, finite-height vortices that extend only one scale height above and below the midplane are also unstable, but for a different reason: we have isolated an antisymmetric (with respect to the midplane) eigenmode that grows with an e-folding time of only a few orbital periods; the nonlinear evolution of this instability leads to the destruction of the vortex. Serendipitously, we observe the formation of 3D vortices that are centered not in the midplane, but at one to three scale heights above and below. Breaking internal gravity waves create vorticity; anticyclonic regions of vorticity roll-up and coalesce into new vortices, whereas cyclonic regions shear into thin azimuthal bands. Unlike the midplane-centered vortices that were placed ad hoc in the disk and turned out to be linearly unstable, the off-midplane vortices form naturally out of perturbations in the disk, and are stable and robust for many hundreds of orbits.