Based on the extended Huygens-Fresnel integral, we have analytically and numerically investigated the propagation properties of a radial phased-locked rotating elliptical Gaussian (RPLREG) beam array in turbulent atmosphere. The average intensity and effective beam sizes in x and y directions of a RPLREG beam array are derived, and their evolution behaviors are analyzed in this paper. Our numerical results indicate that the propagation of a RPLREG beam array in turbulent atmosphere depends on the beam parameters including w, r, N and the structure constant Cn2 of atmospheric turbulence. The results show that the beam arrays of intensity distribution with smaller radius r or larger initial beam size w are very alike, and stronger atmospheric turbulence makes the RPLREG beam evolve into a Gaussian-like beam more rapidly as the propagation distance increases. This research may be useful for optical communications and remote sensing in turbulent atmosphere.