In geophysics, numerical models are generally initialized through so-called data assimilation methods. They require computation of a distance between model fields and physical observations. The most common choice is the Euclidian distance. However, due to its local nature it is not well suited for capturing position errors. This papers investigates theoretical aspects of the use of the optimal transport-based Wasserstein distance in this context and shows that it is able to capture such errors.