Abstract: A numerical model has been designed to study the storm surge induced by typhoon along the coast of Taiwan. The governing equations have been expressed in spherical coordinate system, and a finite difference method has been used to solve them. In the system of hydrodynamical equations, the nonlinear advection and lateral eddy viscosity terms are prominent in shallow coastal waters. Air pressure gradient and wind stresses are the driving forces in the model of typhoon surge. The model has been verified with storm surges induced by Typhoons Herb in 1996, and by typhoons Kai-Tak and Bilis in 2000.