Below I present a more topological argument. The underlying theme is that the density of $\Gamma$ really is a property in the realm of topological dynamics (see Topological Invariance of Unique Ergodicity).
Claim: Let $G$ be a compact metrizable topological group, $\Gamma\leq G$ be a subgroup. Consider the statements
- $\Gamma$ is dense.
- The left translation action $l^\Gamma_\bullet:\Gamma\curvearrowright G$ of $\Gamma$ on $G$ is uniquely ergodic (with the Haar probability measure $\eta_G$ the unique $l^\Gamma$ invariant probability measure).
- $l^\Gamma$ is ergodic w/r/t $\eta_G$ on $G$.
Then $1\iff 2\implies 3$, and if additionally $G$ is connected, then $3\implies 1$ as well.
Proof: Since $G$ is a compact metrizable topological group, there is a distance function $d:G\times G\to\mathbb{R}_{\geq0}$ compatible with its topology that is invariant under both left translations and right translations (see e.g. Einsiedler & Ward's Ergodic Theory, with a view towards Number Theory, Lemma C.2 on p.430). Note that $\eta_G$ is $l^\Gamma$-invariant.
($1\implies 2$) Suppose $\Gamma$ is dense, and let $\mu$ be an arbitrary $l^\Gamma$-invariant Borel probability measure on $G$. We claim that $\mu=\eta_G$. Let $f:G\to\mathbb{R}$ be continuous and $g\in G$. It suffices to show that
$$\int_G f(gx)\, d\mu(x) = \int_G f(x)\, d\mu(x).$$
Since $G$ is compact, $f$ is uniformly continuous, and together with the density of $\Gamma$ we have:
$$\forall \epsilon\in\mathbb{R}_{>0},\exists \delta\in\mathbb{R}_{>0}, \exists \gamma=\gamma_{g,\delta}\in \Gamma: d(\gamma,g)<\delta \,\,\text{ and }\,\,\sup_{x\in G} |f(gx)-f(\gamma x)|<\epsilon.$$
Then for $\epsilon\in\mathbb{R}_{>0}$ arbitrary we have
\begin{align*}
&\left|\int_G f(gx)\, d\mu(x) - \int_G f(x)\, d\mu(x)\right|\\
&= \left|\int_G f(gx)\, d\mu(x) - \int_G f(\gamma x)\, d\mu(x)\right|\\
&\leq \int_G |f(gx)-f(\gamma x)|\, d\mu(x)
<\epsilon.
\end{align*}
Thus any $l^\Gamma$-invariant measure is also $l^G$-invariant, hence is $\eta_G$.
($2\implies 1$) Suppose $\Gamma$ is not dense and put $\overline{\Gamma}=H\lneq G$ and $U=G\setminus H\neq\emptyset$ be the (set-theoretical) complement of $H$. Then $H$ is a compact metrizable group that is $l^\Gamma$-invariant, denote by $\lambda^\Gamma:\Gamma\curvearrowright H$ the restriction of $l^\Gamma$. The Haar probability measure $\eta_H$ of $H$ is $\lambda^\Gamma$-invariant. Define a new Borel probability measure on $G$ by
$$\nu: \mathcal{B}(G)\to [0,1], B\mapsto \eta_H(B\cap H).$$
Then $\nu$ is $l^\Gamma$-invariant. $\eta_G$ is also $l^\Gamma$-invariant and due to the regularity of the Haar measure (see measure of open set with measure Haar or Haar measure on compact group completely positive) we have $\eta_G(U)>0=\nu(U)$, so that $\eta_G$ and $\nu$ are distinct $l^\Gamma$-invariant measures. Consequently $l^\Gamma$ can't be uniquely ergodic.
($2\implies 3$) This follows from the standard argument that the ergodic Borel probability measures of an action is exactly the extremal points of the Choquet simplex of all invariant Borel probability measures of the action (unique ergodicity means that this simplex is a singleton).
($3\implies 1$) Suppose now additionally that $G$ is connected. Suppose $\Gamma$ is not dense; put $\overline{\Gamma}=H\lneq G$, $U=G\setminus H\neq\emptyset$ the complement of $H$, and define
$$\phi=d(\bullet,H): G\to\mathbb{R}_{\geq0}, x\mapsto \inf_{y\in H} d(x,y).$$
Then $\phi$ is continuous and for any $x\in G$ and for any $h_1,h_2\in H$, $\phi(h_1xh_2)=\phi(x)$; in particular $\phi$ is $l^\Gamma$-invariant. If $\eta_G$ were ergodic w/r/t $l^\Gamma$, then $\phi$ would need to be constant $\eta_G$-ae. In particular since $\eta_G$ is regular we have $\eta_G(U)=1$, and consequently there is a unique $c\in\mathbb{R}_{>0}$ ($c=\int_G \phi(x)\, d\eta_G(x)$) such that $K=\{x\in G\,|\, \phi(x)=c\}\subseteq U$ has full $\eta_G$-measure. $K$ is closed, hence $U\setminus K$ is open. Since $\eta_G(U\setminus K)=0$, we have that $U=K$. Thus $G=H\uplus K$ disconnects $G$, a contradiction.