This study aimed to analyze the statistical characteristics of wave heights, wave energy and wave steepness, in order to investigate the wave climate changes around Taiwan Waters, especially for extreme events of big waves. The operational observation of Taiwan sea waves was initiated by the Central Weather Bureau in 1998; however, due to insufficient data length and low data space coverage, the data are unable to serve as references for long-term wave climate change research. Hence, this study adopted the SWAN (Simulation of Wave in Nearshore) Numerical Wave Hindcasting Method, which is a common method used in many studies, to hindcast the history of a wave field. The re-analysis on wind field data of the last 60 years (1948∼2008), published by the National Centers for Environmental Prediction (NCEP), was employed to make the wind field grid consistent with the hindcast wave field grid. Moreover, the Typhoon Wind Field Grid Down Scaling technique proposed by Winter & Chiou (2007) was applied to interpolate a U 10 analysis field that better fits an actual typhoon wind field. The hindcast wave data were compared and validated with directional spectra, which were observed by the meteorological/oceanographic data buoys set up by the Central Weather Bureau and Water Resources Agency since 1997. Longdong, Hualien and Hsinchu Stations were chosen to represent the wave characteristics of sea areas around the island of Taiwan. According to observation data, model parameters were adjusted so that the hindcast results could be closer to observed data in Taiwan sea areas.