Regional wave overtopping risk of Shuang-Chun coast which locates at the southwest of Taiwan was estimated. The estimation is based on a combination of a statistical model, a wave propagation model and empirical formula of runup height. The statistical model considered the dependences by using copula functions between affective factors to generate simultaneously offshore wave and surge conditions which were used as the input conditions of the CMS-WAVE model to simulate the wave propagation on a complicated coastal bathymetry area. Then, runup heights were estimated empirically. The results demonstrated that top elevations of the constructed seawalls at the study site are too conservative due to the perfect dependence assumption between wave and surge. On the other hand, the high overtopping risk areas were identified and also the risk level was estimated for the natural beach area.