For two random variables $P$ and $Q$ over $R^d$ with distributions $p$ and $q$, respectively, the total variation is defined as
$$ TV(P,Q)=\frac{1}{2}\int_{R^d}\ |p(x)-q(x)|dx. $$
Consider the case where $P$ and $Q$ are Gaussian mixtures each with $2$ components, i.e.,
$$P\sim \pi_1\mathcal{N}(\mu_1,\Sigma_1)+\pi_2\mathcal{N}(\mu_2,\Sigma_2),\\ Q\sim \pi'_1\mathcal{N}(\mu'_1,\Sigma'_1)+\pi'_2\mathcal{N}(\mu'_2,\Sigma'_2)$$.
My question is: can we obtain an upper bound on $TV(P,Q)$ based on the total variations between the Gaussian components, e.g., $TV(\mathcal{N}(\mu_1,\Sigma_1),\mathcal{N}(\mu'_1,\Sigma'_1))$ and $TV(\mathcal{N}(\mu_2,\Sigma_2),\mathcal{N}(\mu'_2,\Sigma'_2))$.
Can triangle inequality ($|a+b|\leq|a|+|b|$) be used here for the choices of $a=\pi_1\mathcal{N}(\mu_1,\Sigma_1)-\pi'_1\mathcal{N}(\mu'_1,\Sigma'_1)$ and $b=\pi_2\mathcal{N}(\mu_2,\Sigma_2)-\pi'_2\mathcal{N}(\mu'_2,\Sigma'_2)$?