Whereas classical invariance principles for ergodic Markov chains address the situation in which the time horizon of observations is much larger than the mixing time, the quality of approximation is questionable when this is not the case anymore — even when starting the Markov chain in the invariant law. In this article, we prove quantitative and functional limit theorems for additive functionals along triangular arrays of stationary Gaussian Markov processes when the mixing time t_mix scales sub-, super- and proportionately to the number of observations n. Our major finding is a phase-transition at t_mix asymp n, together with the identification and interrelation properties of the emerging new limit processes at and before the mixing time.