While analyzing the interaction between the hydrological processes of surface and groundwater flow, it is seen that there is a big difference in its interaction in the space and time. These processes can be represented in a qualitative form as part of the hydrological cycle, demonstrating its dependences and importance in the hydrological balance. This work does the numerical coupling of the surface and groundwater flow. This work adopted the SWMS_3D numerical model (Simunek et. al., 1995), which resolves the Richards equation for saturated and non saturated porous media flow in 3D. In order to simulate the superficial dynamic flow, two models from Saint-Vennat equation were developed, these models are: the cinematic wave model and the diffusion model. These two models consider the average outflow in sections in a 2D scenario. For the numerical solution the finite element method was adopted through the Galerkin formulation. Adopting a 3D domain mesh of tetrahedral elements, seen from above, in 2D, we can see a triangular element mesh. The superficial flow model uses the triangular mesh and iterates with the SWMS_3D modified software, which uses the tetrahedral elements mesh. This was done by changes in the boundary conditions to the models. The SWMS_3D will answer with a flow portion corresponding to a sink or source action in the surface, in each time step. Finally the generated model is able to quantify in space and in time the pressure head in the study domain.