Metasurfaces are zero thickness structures and represent general electromagnetic discontinuity. Conventional boundary conditions, and therefore computational techniques, cannot handle such discontinuities. Generalized Sheet Transition Conditions (GSTCs) have been developed for their modeling. In this paper, we propose an scheme based on FDTD for efficient simulation of dispersive metasurfaces which is extension of our previous works on space-time varying metasurfaces. The scheme is proved by two illustrative examples.