In this note, we give a simple derivation of the formula obtained in Dowson and Landau (1982), Olkin and Pukelsheim (1982) and Givens and Shortt (1984) for the quadratic Wasserstein distance between two Gaussian distributions on $\R^d$ with respective covariance matrices $Σ_μ$ and $Σ_ν$. This derivation relies on the existence of an orthogonal matrix $O$ such that $O^*Σ_μO$ and $O^*Σ_νO$ share the same correlation matrix and on the simplicity of optimal couplings in the case with the same correlation matrix and therefore the same copula.