[guided]Fix $k\in\mathbb N$. The goal is to compare the two limiting stationary measures $\bar\rho_0$ and $\bar\rho_1$ to the two time-dependent measures $\alpha_k(t)$ and $\beta_k(t)$, because the contraction assumption controls only the distance between the time-dependent solutions. The bridge is the triangle inequality for the metric $W_2$ on $\mathcal P_2(\mathbb R)$. Applying it first from $\bar\rho_0$ to $\alpha_k(t)$, then from $\alpha_k(t)$ to $\beta_k(t)$, and finally from $\beta_k(t)$ to $\bar\rho_1$, we get
\begin{align*}
W_2(\bar\rho_0,\bar\rho_1)\le W_2(\bar\rho_0,\alpha_k(t))+W_2(\alpha_k(t),\beta_k(t))+W_2(\beta_k(t),\bar\rho_1).
\end{align*}
The contraction hypothesis applies to the selected solutions $\alpha_k$ and $\beta_k$, whose initial values are $\alpha_k(0)=\rho_{0,k}$ and $\beta_k(0)=\rho_{1,k}$. Hence, for every $t\ge 0$,
\begin{align*}
W_2(\alpha_k(t),\beta_k(t))\le W_2(\rho_{0,k},\rho_{1,k}).
\end{align*}
Substituting this estimate into the triangle inequality gives
\begin{align*}
W_2(\bar\rho_0,\bar\rho_1)\le W_2(\bar\rho_0,\alpha_k(t))+W_2(\rho_{0,k},\rho_{1,k})+W_2(\beta_k(t),\bar\rho_1).
\end{align*}
Now we use exactly what the basin hypotheses provide. Since $\rho_{0,k}\in U_0$, the admissible solution $\alpha_k$ converges to $\bar\rho_0$ in $W_2$, so
\begin{align*}
\lim_{t\to+\infty}W_2(\bar\rho_0,\alpha_k(t))=0.
\end{align*}
Since $\rho_{1,k}\in U_1$, the admissible solution $\beta_k$ converges to $\bar\rho_1$ in $W_2$, so
\begin{align*}
\lim_{t\to+\infty}W_2(\beta_k(t),\bar\rho_1)=0.
\end{align*}
The term $W_2(\rho_{0,k},\rho_{1,k})$ is independent of $t$. By the ordinary limit law for sums of real-valued functions, passing to the limit $t\to+\infty$ therefore yields
\begin{align*}
W_2(\bar\rho_0,\bar\rho_1)\le W_2(\rho_{0,k},\rho_{1,k}).
\end{align*}[/guided]