- $c$ is bounded, i.e., $\sup_{x, y\in X} c(x, y) < C$ for some $C \in \mathbb R$.
Then $\|f\|_{\mathrm{Lib}} \le 1$ implies $f$ is bounded, and thus $f \in L_1 (|\mu- \nu|)$ and $(-f, f) \in \Phi_{c}$. By Kantorovich duality, it suffices to prove
$$
\sup _{\Phi_{c}} \mathbb J(\varphi, \psi) = \sup \left \{ \int_X f d (\mu - \nu) \,\middle\vert\, f \in L_1 (|\mu- \nu|) , \|f\|_{\mathrm{Lib}} \le 1 \right \}.
$$
If $f \in L_1 (\mu)$, then
$$
f^c:Y \to \mathbb R \cup \{\pm \infty\}, y \mapsto \inf_{x \in X} [c(x, y) - f(x)]
$$
is such that $\|f^c\|_{\mathrm{Lib}} \le 1$ and $f^{cc} = -f^c$. By this lemma, we have
$$
\sup _{\Phi_{c}} \mathbb J(\varphi, \psi) = \sup _{f \in L_1 (\mu)} \mathbb J(f^{cc}, f^c) = \sup _{f \in L_1 (\mu)} \mathbb J(-f^{c}, f^c) \le \sup _{\|f\|_{\mathrm{Lib}} \le 1} \mathbb J(-f, f) \le \sup _{\Phi_{c}} \mathbb J(\varphi, \psi).
$$
- $c$ is unbounded.
We have
$$
\sup \left \{ \int_X f d (\mu - \nu) \,\middle\vert\, f \in L_1 (|\mu- \nu|) , \|f\|_{\mathrm{Lib}} \le 1 \right \} \le \sup _{\Phi_{c}} \mathbb J(\varphi, \psi).
$$
It remains to prove
$$
\sup \left \{ \int_X f d (\mu - \nu) \,\middle\vert\, f \in L_1 (|\mu- \nu|) , \|f\|_{\mathrm{Lib}} \le 1 \right \} \ge \sup _{\Phi_{c}} \mathbb J(\varphi, \psi).
$$
Let $c_n := \min\{n, c\}$. Then $c$ is a bounded metric on $X$ such that $c_n \nearrow c$. By this result, we have
$$
\lim_n \inf_{\pi \in \Pi(\mu, \nu)} \int c_n d\pi = \inf_{\pi \in \Pi(\mu, \nu)} \int c d\pi = \sup _{\Phi_{c}} \mathbb J(\varphi, \psi).
$$
On the other hand,
$$
\begin{align}
\inf_{\pi \in \Pi (\mu, \nu)} \int_{X \times Y} c_n d \pi &= \sup \left \{ \int_X f d (\mu - \nu) \,\middle\vert\, f \in L_1 (|\mu- \nu|) , \|f\|_{\mathrm{Lib}, c_n} \le 1 \right \} \\
&\le \sup \left \{ \int_X f d (\mu - \nu) \,\middle\vert\, f \in L_1 (|\mu- \nu|) , \|f\|_{\mathrm{Lib}} \le 1 \right \}.
\end{align}
$$
This completes the proof.
No, this is not true. Unless $X$ is trivial, you can find a (nonatomic) measure $\lambda$ on $X$ and a sequence $(f_n) \subset L^2(\lambda)$
with $f_n \rightharpoonup 0$ and (say) $f^+_n \rightharpoonup 1$, both in $L^2(\lambda)$. Now, consider $\mu_n = f_n \lambda$.
Finally, if $g_n \rightharpoonup g$ in $L^2(\lambda)$, then
$$
\int_X (g_n - g) h \, \mathrm{d}\lambda \to 0
$$
for all $h \in L^2(\lambda)$.
Thus, if $\lambda$ is finite, this holds for all bounded, continuous functions $h$. Hence, the measures $g_n \lambda$ converge weak-$*$ towards $g \lambda$.
Best Answer
We have $$ W^p_1 (\mu, \nu) \le W^p_p (\mu, \nu) \le (\mathrm{diam} X)^p W_1 (\mu, \nu) \quad \forall p \in [1, +\infty). $$
So it suffices to prove for $p=1$. Let $$ \mathrm{Lib} (f) := \sup_{x\neq y} \frac{|f(x) - f(y)|}{|x-y|}. $$
By Kantorovich-Rubinstein theorem, we have $$ W_1 (\mu_m, \nu) = \sup \left \{\int_X f \mathrm d (\mu_m-\nu) \,\middle\vert\, f \in L_1 (|\mu_m-\nu|) , \mathrm{Lib} (f) \le 1\right \}. $$
Let $f$ be bounded Lipschitz-continuous such that $\mathrm{Lib} (f) > 0$. Then $\frac{f}{\mathrm{Lib} (f)}$ is $1$-Lipschitz. Then $$ \int_X \frac{f}{\mathrm{Lib} (f)} \mathrm d (\mu_m-\nu) \to 0. $$
So $$ \int_X f \mathrm d (\mu_m-\nu) \to 0, $$
The result then follows from Portmanteau theorem.
Let $m_k$ be a subsequence such that $$ \lim_k W_1 (\mu_{m_k}, \mu) = \limsup_m W_1 (\mu_m, \mu). $$
Fix $a \in X$. Let $g_k$ be a Lipschitz funtion such that $\mathrm{Lib} (g_k) \le 1$ and $$ W_1 (\mu_{m_k}, \mu) \le \int_X g_k \mathrm d (\mu_{m_k} - \mu) + \frac{1}{k} = \int_X (g_k - g_k (a)) \mathrm d (\mu_{m_k} - \mu) + \frac{1}{k} . $$
It follows from $(\mu_{m_k} - \mu) (X) = 0$ that $$ \int_X g_k (a) \mathrm d (\mu_{m_k} - \mu) = 0 \quad \forall k. $$
Let $g'_k := g_k - g_k (a)$. Then $\mathrm{Lib} (g'_k) \le 1$ and $(g'_k)_k$ is bounded by $\mathrm{diam} X$. By Arzelà –Ascoli theorem, there is $g$ and a subsequence (re-labelled) such that $g'_k \to g$ in $\|\cdot\|_\infty$. Clearly, $\mathrm{Lib} (g) \le 1$. Finally, $$ \begin{align} \int_X g'_k \mathrm d (\mu_{m_k} - \mu) &= \int_X (g'_k-g) \mathrm d (\mu_{m_k} - \mu) + \int_X g \mathrm d (\mu_{m_k} - \mu) \\ &= \int_X (g'_k-g) \mathrm d (\mu_{m_k} - \mu)^+ - \int_X (g'_k-g) \mathrm d (\mu_{m_k} - \mu)^- + \int_X g \mathrm d (\mu_{m_k} - \mu) \\ &\le 2 \|g'_k-g\|_\infty (\mu_{m_k} - \mu)^+ (X) + \int_X g \mathrm d (\mu_{m_k} - \mu) \\ &\le 4 \|g'_k-g\|_\infty + \int_X g \mathrm d (\mu_{m_k} - \mu). \end{align} $$
The result then follows by taking the limit $k \to \infty$.