ABSTRACT. This work is a contribution to a better understanding of the effect that domestic sewage discharges may cause in a water body, specifically at Igapó I Lake, in Londrina, Paraná State, Brazil. The simulation of the dynamics of pollutant concentrations throughout the water body was conducted by means of structured discretization of the geometry of Igapó I Lake, together with the finite differences and the finite elements methods. Firstly, the hydrodynamic flow (without the pollutants), modeled by NavierStokes and pressure equations, was numerically resolved by the finite differences method, and associated with the fourth order Runge-Kutta procedure. After that, by using the hydrodynamic field velocity, the flow of the reactive species (pollutants) was described through a reaction transport model, restricted to the carbon-nitrogen cycle. The reaction transport model was numerically resolved by the stabilized finite elements method, by means of a semi-discrete formulation. A qualitative analysis of the numerical simulations provided a better understanding of the dynamics of the processes involved in the flow of the reactive species, such as the dynamics of the nitrification process, of the biochemical demand of oxygen and of the level of oxygen dissolved in the water body at Igapó I Lake.Keywords: water quality model, pollutants, numerical simulation, finite differences method, finite elements method.Uma modelagem do transporte do ciclo carbono-nitrogênio no Lago Igapó I -Londrina, Estado do Paraná, Brasil RESUMO. Este artigo é uma contribuição para um melhor entendimento do efeito que uma descarga de esgoto doméstico pode causar num corpo d'água, em particular, no lago Igapó I, Londrina, Paraná, Brasil. A simulação da dinâmica das concentrações dos poluentes em todo o corpo d'água é realizada por meio de uma discretização estruturada da geometria do lago Igapó I, juntamente com os métodos de diferenças finitas e de elementos finitos. Primeiramente, o escoamento hidrodinâmico (sem os poluentes), modelado pelas equações de Navier-Stokes e de pressão, é resolvido numericamente pelo método de diferenças finitas, associado ao procedimento de Runge-Kutta de quarta ordem. Em seguida, utilizando o campo de velocidades fornecido pelo modelo hidrodinâmico, descreve-se o escoamento de espécies reativas (restrito ao ciclo carbono-nitrogênio), por meio de um modelo de transporte de poluentes com reações. O modelo de transporte de poluentes com reações é resolvido numericamente pelo método de elementos finitos estabilizados, através de uma formulação semi-discreta. Uma análise qualitativa das simulações numéricas proporcionou uma melhor compreensão da dinâmica dos processos envolvidos no escoamento de espécies reativas, tais como a dinâmica do processo de nitrificação, da demanda bioquímica de oxigênio e do nível de oxigênio dissolvido no corpo d'água do lago Igapó I.Palavras-chave: modelo de qualidade de água, poluentes, simulação numérica, método de diferenças finitas, método de elementos finitos.