Partially saturated flow in a porous medium is typically modeled by the Richards equation, which is nonlinear, parabolic and possibly degenerated. This paper presents domain decomposition-based numerical schemes for the Richards equation, in which different time steps can be used in different subdomains. Two global-in-time domain decomposition methods are derived in mixed formulations: the first method is based on the physical transmission conditions and the second method is based on equivalent Robin transmission conditions. For each method, we use substructuring techniques to rewrite the original problem as a nonlinear problem defined on the space-time interfaces between the subdomains. Such a space-time interface problem is linearized using Newton’s method and then solved iteratively by GMRES; each GMRES iteration involves parallel solution of time-dependent problems in the subdomains. Numerical experiments in two dimensions are carried out to verify and compare the convergence and accuracy of the proposed methods with local time stepping.