In this paper, novel analytical closed-form expressions are derived for the probability density function of the sum of identically distributed correlated gamma-gamma random variables that models an optical atmospheric channel communication with receiver spatial diversity. The mathematical expressions here proposed provide a general procedure to obtain information about the scintillation effects induced by turbulence over a diversity reception scheme implementing equal-gain combining method. Both, validity and accuracy of the obtained statistical distribution are corroborated by comparing the analytical results to numerical results obtained by Monte-Carlo simulations. These simulations are particularized for constant, exponential and circular correlation models, corresponding to three different receivers spatial configurations. In addition, the extreme situations of no correlation and fully correlated received signals are also studied. The presented expressions lead to a simple and easy-to-compute analytical procedure of analyzing atmospheric optical communications systems with correlated spatial diversity.