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.