For two-dimensional (2D) time fractional diffusion equations, we construct a numerical method based on a local discontinuous Galerkin (LDG) method in space and a finite difference scheme in time. We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable. Numerical results indicate the effectiveness and accuracy of the method and confirm the analysis.