Bridging Gaussian and non-Gaussian Data Assimilation for High-Dimensional Geophysical Models