  We consider the two dimensional surface quasi-geostrophic equations with
super-critical dissipation. For large initial data in critical Sobolev and
Besov spaces, we prove optimal Gevrey regularity with the same decay exponent
as the linear part. This settles several open problems in \cite{Bis14, BMS15}.
