  We derive explicitly the adapted $2$-Wasserstein distance between
non-degenerate Gaussian distributions on $\mathbb{R}^N$ and characterize the
optimal bicausal coupling(s). This leads to an adapted version of the
Bures-Wasserstein distance on the space of positive definite matrices.
