Abstract: Based on Galerkin finite element method, a new upwinding finite element numerical model suitable for modeling wind-driven current in shallow lakes is derived by introducing streamline-upwinding method and the selective lumping coefficient matrix scheme. Moreover it is applied to Taihu Lake, and the morphological characteristic of wind-driven current is discussed.