The unified theory of coherence and polarization serves as the natural ground for the statistical theory of random electromagnetic beams.[1] Accordingly, the coherent mode representation for the scalar field is also generalized to the vector field, however, the procedure to determine vector modes becomes rather involved compared with scalar modes, and requires solving a vector integral equation. [2][3] On the other hand, a different mode representation has been introduced, which consists of two scalar coherent mode representations for the diagonal elements and one bi-modal expansion for the off-diagonal elements of the cross spectral density matrix. [4]In this paper, we consider the connection between these two mode representations. In particular, we suggest a method to find the vector modes from the scalar modes and formulate the cross spectral density matrix as a correlation matrix.Let us consider a random electromagnetic beam propagating close to the z -direction. Its coherence properties and polarization properties may be expressed in terms of its Here, and are components of a typical realization of a monochromatic electric field of frequency along two mutually orthogonal directions perpendicular to the z -direction, asterisk denotes the complex conjugate and the sharp brackets denote ensemble average, in the sense of coherence theory in the space-frequency domain [5].The coherent mode representation for a scalar field may be extended to an electromagnetic beam and the cross spectral density matrix can be expanded as an absolutely and uniformly convergent Mercer-type series [2]where ( )are the eigenvalues and eigenvectors of the Fredholm integral equationwhere the integration is carried out over the domain . This coherent mode expansion is similar to the KarhunenLoeve expansion for a complex random process [6]. The eigenvalues are nonnegative and the eigenfunctions satisfy the orthogonality conditionVector It is shown that the two mode representations, one with vector modes and the other with scalar modes, for the cross spectral density matrices of electromagnetic beams are equivalent to each other. In particular, we suggest a method to find the vector modes from the scalar modes and formulate the cross spectral density matrix as a correlation matrix.