I am a bit uncomfortable with the habit of writing $$f(z)=u(x,y)+i v(x,y)$$ where $f:U\subseteq \mathbb C \to \mathbb C$ is a complex function and $z=x+iy$. This is especially done when Cauchy-Riemann equations and harmonic functions are introduced. The issue is that I think authors switch from subsets of $\mathbb C$ to $\mathbb R^2$ ad vice-versa improperly.
For example, in Theorem VIII, $2.2$ of his book "Complex Analysis", Conway writes that the connected and open set $G\subseteq \mathbb C$ is simply connected if and only if every harmonic function $u: G\to \mathbb C$ has a harmonic conjugate, i.e. there is another harmonic function $v: G\to \mathbb C$ such that $f=u+iv$ is holomorphic in $G$.
This approach does not convince me, because for a function to be harmonic one should be able to perform derivatives with respect to $x$ and $y$, so should consider a function of two real variables.
This also makes me question about Wirtinger derivatives $$\frac{\partial}{\partial z}=\frac{\partial}{\partial x}-i\frac{\partial}{\partial y} \\ \frac{\partial}{\partial \bar z}=\frac{\partial}{\partial x}+i\frac{\partial}{\partial y}$$ because this would be appropriate for a function of the three variables $x,y,z$.
Here we discussed the meaning of $\partial/\partial \bar z$ and it looks to me that it is only a trick that works.
How can one be more rigorous? I thought we could introduce a set $\tilde G=\{(x,y) \in \mathbb R^2: z=x+iy \in G\}$ and define $u,v$ there. Is there another approach, perhaps more frequent in literature?