This paper introduces four matrix normal distributions on analytic bundles of flag varieties, extending the separable covariance $\varPhi \otimes \varPsi$ with potentially variable-level ($\varPsi$) and/or sample-level ($\varPhi$) correlations. The joint distribution of sample variances and covariances, leading to the product moment distribution, is considered when precision matrices admit a specific form. Several well-known consequences, including the non-central Wishart distribution and normal quadratic forms, now appear as corollaries.