Abstract: A new two-dimensional upwinding finite element numerical model,suitable for modeling concentrations of total phosphorus(TP) in shallow lakes was derived,from the starting point of control equation,based on Galerkin finite element method and combined with the new development of the finite element method.The model was applied in Lake Taihu for the study of TP distribution in the lake water body.