In general, it's always true (for a connected scheme) that $H^1_{et}(X, \mathbb{Z}/l \mathbb{Z}) = \hom(\pi_1^{et}(X), \mathbb{Z}/l\mathbb{Z})$ (not compactly supported). Taking inverse limits over $l$ then gives the claim.
The reason this is true is that $H^1_{et}(X, \mathbb{Z}/l\mathbb{Z}$) can be computed by Cech cocycles, and from this it follows that elements of this group classify torsors over $\mathbb{Z}/l\mathbb{Z}$ (i.e. sheaves with a $\mathbb{Z}/l\mathbb{Z}$-action which are locally the constant $\mathbb{Z}/l \mathbb{Z}$ (in the etale topology)). But these, by descent theory, are the same as Galois covers of $X$ with Galois group $\mathbb{Z}/l\mathbb{Z}$, and (by Galois theory) classified by maps from the etale fundamental group into $\mathbb{Z}/l\mathbb{Z}$.

