Rockfall source areas identification of high steep fractured rock slope is a challenging task for hazard mitigation. Considering this difficulty, a new method integrating discontinuity-based threshold slope angle and rockfall trajectory analyses is introduced in this study and applied to a railway bridge in China. First, abundant discontinuity data are collected by means of the modern unmanned aerial vehicle (UAV) digital photogrammetry technique. Then, a threshold slope angle used to delineate rockfall-prone areas is comprehensively determined by probabilistic kinematic and block theory analyses. Next, physically based rockfall trajectory simulation with consideration of mean block size is used to distinguish the hazard level of each release point. Finally, rockfall sources causing moderate and high risks to the train tunnel entrance at the base of high steep fractured rock slope are determined using a semiquantitative risk matrix approach. The final rockfall source map can provide a quantitative reference for hazard mitigation. At the end of this paper, some issues associated with the proposed approach are discussed, such as parameter importance analysis of threshold slope angle, comparison between the two probabilistic approaches and uncertainties in rockfall source areas identification. •A method for rockfall sources identification of jointed rock slope is proposed.•This method integrates threshold slope angle and rockfall trajectory analyses.•Probabilistic approaches are developed to determine a final threshold slope angle.•A multiscanline method is presented to determine the volume of a launched block.•Trajectory dense and sparse areas are used to distinguish different hazard levels.