Abstract: An effective unconditionally stable implementation of the auxiliary differential equation Crank-Nicolson-approximate-decoupling finite-difference time-domain (ADE-CNAD-FDTD) algorithm for 2-D anisotropic magnetized plasma is proposed. The conventional ADE-FDTD method for 1-D anisotropic dispersive media has high efficiency and accuracy. This paper extends this method to 2-D anisotropic magnetized plasma with the CNAD scheme. The proposed formulations not only solves the problem that incorporates both anisotropy and frequency dispersion at the same time, but also eliminates the Courant-Friedrich-Levy (CFL) stability constraint. A numerical example has been carried out to validate the proposed formulations in the 2-D FDTD domain composed of anisotropic magnetized plasma. The results prove that the proposed formulations significantly save time and perform stably with acceptable accuracy.