Abstract. Rockfall hazard is frequent along the national road (G318) in west Hubei, China. To understand the distribution and potential hazard probability, this study combines the result of a 3-years engineering geological investigation, statistical modeling, and kinemics-based method to identify risky road sections. Rockfall hazard probability is calculated by integrating spatial, temporal, size probability, and reaching probabilities of source areas. Rockfall source areas are preliminarily identified first by slope angle threshold (SAT) analysis. Random Forest model (RFM) and multivariate logistic regression model (MLRM) are then applied and compared to get the final susceptible source areas, considering eight factors, including slope, aspect, elevation, lithology, joint density, slope structure, land-use type, distance to the road. Temporal and size probability of source areas are separately obtained by Poisson distribution and power-law distribution theory. An important parameter (reach angle) for rockfall trajectory simulation was determined by back analysis in Flow-R and validated by field investigation. The results show good fitness with the measurements by field investigation. In the conditions of 5, 20, and 50 years return period, potential risky road sections are found out under two size scenarios (larger than 1 000 m3, 10 000 m3). This research helps the local government to completely understand the rock falls from source area existence and potential risk to roads.