By exploiting favorable characteristics of a uniform cross-array, a passive localization algorithm of narrowband cyclostationary sources in the spherical coordinates (azimuth, elevation, and range) is proposed. Firstly, we construct a parallel factor (PARAFAC) analysis model by computing the third-order cyclic moment matrices of the properly chosen sensor outputs. Then, we analyze the uniqueness of the constructed model and obtain three-dimensional (3D) near-field parameters via trilinear alternating least squares regression (TALS). The investigated algorithm is well suitable for the localization of the near-field cyclostationary sources. In addition, it avoids the multidimensional search and pairing parameters. Results of computer simulations are carried out to confirm the satisfactory performance of the proposed method.