This study was conducted to estimate the design flood by the determination of best fitting order for LH-moments of the annual maximum series at fifteen watersheds. Using the LH-moment ratios and Kolmogorov-Smirnov test, the optimal regional probability distribution was identified to be the Generalized Extreme Value (GEV) in the first report of this project. Parameters of GEV distribution and flood flows of return period n years were derived by the methods of L, L1, L2, L3 and L4-moments. Frequency analysis of flood flow data generated by Monte Carlo simulation was performed by the methods of L, L1, L2, L3 and L4-moments using GEV distribution. Relative Root Mean Square Error. (RRMSE), Relative Bias (RBIAS) and Relative Efficiency (RE.) using methods of L, Ll , L2, L3 and L4-moments for GEV distribution were computed and compared with those resulting from Monte Carlo simulation. At almost all of the watersheds, the more the order of LH-moments and the return periods increased, the more RE became, while the less RRMSE and RBIAS became. The Absolute Relative Reduction (ARR) for the design flood was computed. The more the order of LH-moments increased, the less ARR of all applied watershed became It was confirmed that confidence efficiency of estimated design flood was increased as the order of LH-moments increased. Consequently, design floods for the appled watersheds were derived by the methods of L3 and L4-moments among LH-moments in view of high confidence efficiency.