I think this is a rather open question. Suppose I have bi-dimensional data $(x_i, y_i)$. I have some reasonable model for the marginals, say distributions $F_X$ and $F_Y$ (parametric).
How to reasonably construct a bi-dimensional model for the joint distribution of $(X,Y)$ that respects say, the correlation between $X$ and $Y$ or some other higher moment? Of course I'm not considering normal variates.