This paper proposes a unified statistical channel model to characterize the atmospheric turbulence induced distortions faced by orbital angular momentum (OAM) in free space optical (FSO) communication systems. In this channel model, the self-channel irradiance of OAM modes as well as crosstalk irradiances between different OAM modes are characterized by a Generalized Gamma distribution (GGD). The latter distribution is shown to provide an excellent match with simulated data for all regimes of atmospheric turbulence. Therefore, it can be used to overcome the computationally complex numerical simulations to model the propagation of OAM modes through atmospheric turbulent FSO channels. The GGD allows obtaining very simple tractable closed-form expressions for a variety of performance metrics. Indeed, the average capacity, the bit-error rate, and the outage probability are