We discuss global in time interior H\"older bounds for large data for critical SQG in bounded domains. In order to obtain these, we establish nonlinear lower bounds and commutator estimates for the Dirichlet fractional Laplacian in bounded domains. As an application global existence of weak solutions of SQG were obtained. This is a joint work with Peter Constantin.