In this paper, rectangular multi-Gaussian Schell-model (MGSM) array beams, which consists N×D beams in rectangular symmetry, are first introduced. The analytical expressions of MGSM array beams propagating through free space and non-Kolmogorov turbulence are derived. The propagation properties, such as normalized average intensity and effective beam sizes of MGSM array beams are investigated and analyzed. It is found that the propagation properties of MGSM array beams depend on the parameters of the MGSM source and turbulence. It can also be seen that the beam size of Gaussian beams translated by MGSM array beams will become larger as the total number of terms, M, increases or coherence length, σ , decreases, and the beam in stronger non-Kolmogorov turbulence (larger α and l 0 , or smaller L 0 ) will also have a larger beam size.