Adapted optimal transport between Gaussian processes in discrete time

We derive explicitly the adapted $2$-Wasserstein distance between arbitrary non-degenerate Gaussian distributions on RN 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.