Let $\mathbb H^3$ be hyperbolic $3$-space and let $\Gamma \subset PSL(2,\mathbb R) \subset PSL(2,\mathbb C) \cong Isom ^+(\mathbb H^3)$ be the (Fuchsian) fundamental group of a Riemann Surface $\Sigma$. Then $\Gamma$ also acts freely and properly discontinuously on $\mathbb H^3$, so we may look at the quotient manifold $\mathbb H^3 \backslash \Gamma$. Something that seems to be common knowledge among experts on hyperbolic geometry, but something that I cannot quite wrap my head around, is the following assertion:
There exists an isometric embedding $i: \Sigma \to \mathbb H^3 \backslash \Gamma$.
Let $p: \mathbb H^3 \to \mathbb H^3 \backslash \Gamma$ be the corresponing covering map. Assuming we could indeed find an isometrically embedded copy of $\Sigma$, then each component of $p^{-1}(\Sigma)$ would be an isometrically embedded copy of $\mathbb H^2$ in $\mathbb H^3$ (a hyperbolic plane), invariant under the action of $\Gamma$. However, it is not at all clear to me how such a copy would look like. Using the Poincaré half-plane model of $\mathbb H^3$ (so that topologically, we can identify $\mathbb H^3$ with $\mathbb C \times \mathbb R^+$), hyperbolic planes in $\mathbb H^3$ correspond (precisely?) to planes/ hemispheres perpendicular to the $x_1x_2$-plane. Moreover, according to what I've read in some sources, in this model, the action of $PSL(2,\mathbb R)$ is generated by the following isometries, where $\lambda$ ranges over $\mathbb R$: \begin{align} (z,x_3) \mapsto (\lambda z, |\lambda|x_3): \lambda \neq 0 \\ \\ (z,x_3) \mapsto (z + \lambda, x_3) \\ \\ (z,x_3) \mapsto (\frac{-\overline{z}}{|z|^2 + x_3^2},\frac{x_3}{|z|^2 + x_3^2}) \\ \end{align}
However, I fail to see how any subgroup, generated by some collection of compositions of the above elements (let alone a fuchsian subgroup), would leave any of the aforementioned hyperplanes invariant. Any help is highly appreciated
More a (very) long comment than an answer
Also I am not knowledgable on Fuchsian Group and $\mathbb H^3$ ( hyperbolic $3$-space) but I know a bit about hyperbolic $2$-space and extrapolate a bit from that. so maybe some of my comments are just plain wrong.
But still some , I hope helpful comments: (maybe I am a bit pedantic at some points but that is also to help grow my own understanding)
Be careful with Hyperbolic planes (planes in hyperbolic geometry using the Poincare half plane model) and Hyperplanes (surfaces of points equidistant to an hyperbolic plane) they are not the same, I think you everywhere mean Hyperbolic planes.
You make some mistakes with your indices I guess $ z = x_1 + x_2 i$ and $x_3 > 0 $
Hyperbolic planes in $\mathbb H^3$ correspond to planes/ hemispheres perpendicular to the boundary of the model (the euclidean $x_3 = 0 $ -plane.)
Also I think $\lambda$ ranges over $\mathbb C$ (not just $\mathbb R$ ) otherwise the 3 isometries do not cover all isometries , (for example if $\lambda$ ranges over $\mathbb R$ no isometry can be built that maps the $ z + \overline{z} = 0 $ plane to the $ z - \overline{z} = 0 $ plane)
Then there are some hyperbolic planes that are invariant.
in 1) $$(z,x_3) \mapsto (\lambda z, |\lambda|x_3): \lambda \neq 0 $$
Leaves the plane unit hemisphere invariant if $ |\lambda | =1 $
(the unit hemisphere is $z\overline{z} +x_3^2 = 1$ )
in 3) $$(z,x_3) \mapsto (\frac{-\overline{z}}{|z|^2 + x_3^2},\frac{x_3}{|z|^2 + x_3^2}) $$
Leaves the unit hemisphere, and the planes $ z + \overline{z} = 0 $ and$ z - \overline{z} = 0 $ invariant. (and maybe more important moves every point inside the unit hemisphere to outside the hemisphere and every point outside the hemisphere to inside the hemisphere)
But the real question here is why is that important? Is it not more that some combinations of isometries are able to create invariant hyperplanes not that one of the 3 (primal) isometries can do that.
Please check all I wrote, I am even less knowledgable on Fuchsian Group than you but I hope this hyperbolic geometry view helps