Let
$\alpha_n({\cdot})=\mathbb{P}\bigl(X_{n+1}\in{\cdot}\mid X_1,\ldots,X_n\bigr)$ be the predictive distributions of a sequence
$(X_1,X_2,\ldots)$ of p-dimensional random vectors. Suppose
$\alpha_n=\mathcal{N}(M_n,Q_n)$, where
$M_n=({1}/{n})\sum_{i=1}^nX_i$ and
$Q_n=({1}/{n})\sum_{i=1}^n(X_i-M_n)(X_i-M_n)^{\top}$. Then there is a random probability measure
$\alpha$ on the Borel subsets of
$\mathbb{R}^p$ such that
$\lVert\alpha_n-\alpha\rVert\overset{\mathrm{a.s.}}\longrightarrow 0$, where
$\lVert{\cdot}\rVert$ is the total variation distance. An explicit expression for
$\alpha$ is provided and the convergence rate of
$\lVert\alpha_n-\alpha\rVert$ is shown to be arbitrarily close to
$n^{-1/2}$. Moreover, it is still true that
$\lVert\alpha_n-\alpha\rVert\overset{\mathrm{a.s.}}\longrightarrow 0$ even if
$\alpha_n=\mathcal{L}(M_n,Q_n)$, where
$\mathcal{L}$ belongs to a class of distributions much larger than the normal. The predictives
$\alpha_n$ are useful in various frameworks, including Bayesian predictive inference and predictive resampling. Finally, the asymptotic behavior of copula-based predictive distributions is investigated and a numerical experiment is performed.