In the paper, we establish commutator estimates for the Dirichlet-to-Neumann map of Stokes systems in Lipschitz domains. The approach is based on Dahlberg's bilinear estimates, and the results may be regarded as an extension of [8,19] to Stokes systems.