Estimation of design floods is typically required for hydrologic design purpose. Design floods are routinely estimated for water resources planning, safety and risk of the existing water-related structures. However, the hydrologic data, especially streamflow data for the design purposes in South Korea are still very limited, and additionally the length of streamflow data is relatively short compared to the rainfall data. Therefore, this study collected a large number design flood data and watershed characteristics (e.g. area, slope and altitude) from the national river database. We further explored to formulate a scaling approach for the estimation of design flood, which is a function of the watershed characteristics. Then, this study adopted a Hierarchical Bayesian model for evaluating both parameters and their uncertainties in the regionalization approach, which models the hydrologic response of ungauged basins using regression relationships between watershed structure and model. The proposed modeling framework was validated through ungauged watersheds. The proposed approach have better performance in terms of correlation coefficient than the existing approach which is solely based on area as a predictor. Moreover, the proposed approach can provide uncertainty associated with the model parameters to better characterize design floods at ungauged watersheds.