2D Gaussian Mixture

If you mean a multivariate gaussian of length 2, yes it’s possible to fit a mixture of them. Check out pymc.Mixture — PyMC 5.9.1 documentation