A three-dimensional numerical simulation on wind-driven current in Tauhu Lake is developed in this paper, in which the finite difference method is used.A dimensionless vertical coordinate is applied for the vertical grid regulation so as to eliminate the impacts of the wind-induced water level fluctuation and the irregular lake bottom forms, whereas grids with sawtoothed boundaries is applied along the horizontal dimensions when supposing the boundaries of the lake are closed.The wind data in July and August, 1988 are used in the numerical simulation of water level at Dapukou Station and Xukou Station.The comparison between the simulated and the observed results indicates certain cons is tan cy.Further studies reveal that, when a constant 14 m/s SE wind dominates for 20 hours on the water surface, the water level at Dapukou St.increases 45 cm while the level at Xikou St.decreases 40 cm.The four-layer (surface, upper-middle, lower-middle, near bed) velocity fields obtained respectively show that the absolute velocity of the surface flow is about 1.2 m/s, while the bed flow is between 0.3-0.6 m/s but in the opposite direction.The simulation also reveals that Che velocity field in the upper-middle layer change directions more frequently, which is in agreement with former researches done by foreigners.In conclusion, the wind-driven currents in lakes are very complicated, the representativeness of wind data is an important factor that affects the accuracy of numerical simulation.That the wind-induced current directions are vertically inconsistent and the velocity fields along the horizontal dimensions not evenly distributed is the characteristics of lake flow differed from that of other water bodies.