In the field of underwater acoustic signal processing, broadband array signal processing is an important research issue, which is more complex than that of narrowband. Broadband passive localization of underwater acoustic sources is investigated in this study. In order to improve the ability to locate multiple radiation sources which are emitted from ships, the orthogonal propagator method (OPM) is proposed to estimate direction of arrival (DOA) for the target. Based on the theory of spatial spectrum estimation, the linear decomposition of the covariance matrix is adopted to form the orthogonal propagator, and the spatial spectrum is constructed according to orthogonality of subspaces. Then, the directions of ship radiation noises are determined by searching the extreme points of the spatial spectrum. For broadband underwater acoustic signals, the non-coherent subspace method is used to average the effective frequency points in bandwidth to reduce noise interference, and OPM algorithm operating in broadband underwater acoustic signal, on arbitrary sonar array, is derived. In order to verify this method, numerical simulations are conducted to test its performance in localization accuracy, beam width and spatial resolution. During the simulations, a uniform linear array consisting of 96 elements is designed. Additionally, comparisons between the proposed method and conventional beamforming technique are conducted. Comparison results show that the OPM method performs better in spatial directivity and spatial resolution.