This paper analyzes the joint Rate Distortion Function (RDF) of correlated multivariate Gaussian sources with individual square-error distortions. Leveraging Hotelling's canonical variable form (CVF), presented is a closed-form characterization of the joint RDF, that involves {a system of nonlinear equations. Furthermore, for the special case of symmetric distortions (i.e., equal distortions), the joint RDF is explicitly expressed in terms of} two water-filling variables. The results greatly improve our understanding and advance the development of closed-form solutions of the joint RDF for multivariate Gaussian sources with individual square-error distortions.