Numerical simulations of Rossby wave instabilities in protoplanetary disks have be published by Li et al. (2001) and by Klahr and Bodenheimer (2003). These simulations exhibit the formation of large-scale vortices that may have an important effect on both planetary formation and angular momentum transport in the disk. In order to explore the physical conditions required for such instabilities to occur, I have derived a simplified fluid model based upon the anelastic approximation for a vertically-averaged disk. In this approximation, the rapid pressure waves are filtered out of the model, so that only the slow Rossby wave modes remain. The model consists of nonlinear evolution equations for the fluid vorticity and entropy. Radial gradients in the surface density and the temperature of the disk can drive vorticity evolution via a baroclinic coupling of the two equations. Radiative damping of entropy perturbations from the basic state has also been included.
The linearized equations have a structure similar to classical stratified shear flow problems. However, unlike the usual stratified shear flow, the temperature gradient acts as a destabilizing agent. Furthermore, radiative damping can destabilize the neutrally stable modes. Transient growth of linear disturbances can be evaluated by solving the initial value problem. Surprisingly large transient vorticity and entropy perturbations can form in disks with realistic temperature profiles. The only way to suppress these transient disturbances is to introduce a large turbulent viscosity.