The stability of two superposed buoyancy vortices is studied linearly in a two-level Surface Quasi-Geostrophic (SQG) model. The basic flow is chosen as two circular vortices with uniform buoyancy, coaxial, and the same radius. A perturbation with a single angular mode is added to this mean flow. The SQG equations linearized in perturbation around this basic flow form a two-dimensional ODE for which the normal and singular mode solutions are numerically computed. The instability of these two vortices depends on several parameters. The parameters varied here are: the vertical distance between the two levels and the two values of the vortex buoyancies (called vortex intensity hereafter); the other parameters remain fixed. For normal modes, the...