A relação entre os processos hidrológicos de escoamento superficial e subterrâneo apresenta uma grande variabilidade espacial e temporal. Podendo ser representado de forma qualitativa como parte sequêncial do ciclo hidrológico, estes processos, demostram sua grande dependência e importância nos estudos de balanços hídricos. Visando uma representação quantitativa, este trabalho faz o acoplamento, entre os modelos numéricos de escoamento superficial e de fluxo em meios porosos. Para o meio poroso adotou-se o modelo numérico SWMS_3D (Simunek et al, 1995), o qual resolve a equação de Richards, para fluxo em meios porosos saturados e não saturados nas três dimensões. Na simulação da dinâmica superficial, foram desenvolvidos dois modelos derivados das equações de Saint- Venant: o modelo da Onda Cinemática e o modelo de Difusão. Para a solução numérica foi empregado o método dos elementos finitos através da formulaçao de Galerkin, adotando uma malha tridimensional de elementos tetraédricos, formando uma sub-malha de elementos triangulares na superfície. O modelo de escoamento superficial emprega a malha triangular e interage com o programa SWMS_3D modificado (que utiliza a malha de tetraédros) através das imposições das condições de contorno transientes. Este, responderá com uma parcela de fluxo correspondente à recarga ou descarga no contorno a cada passo de tempo. Com isso, o modelo gerado é capaz de quantificar espacialmente e temporalmente as cargas de pressão em todos os pontos do domínio de estudo.