Cauchy-Riemann equations can be seen as a system of 2 PDEs for two functions on the plane:
$$ \begin{align*} L_{\frac{\partial}{\partial x}}(u) &= L_{\frac{\partial}{\partial y}}(v) \\ L_{\frac{\partial}{\partial y}}(u) &= -L_{\frac{\partial}{\partial x}}(v) \end{align*} $$
where $L_X$ denotes the Lie derivative along vector field $X$. Since the vector fields $\frac{\partial}{\partial x}$ and $\frac{\partial}{\partial y}$ commutes, we have natural integrability conditions:
$$ \begin{align*} \Delta u &= 0\\ \Delta v &= 0, \end{align*} $$
and we know that there exists smooth (even analytic) solutions.
I am surprised that I did not found any results about the natural generalisation of those equations, namely:
$$ \begin{align*} L_A(u)&=L_B(v)\\ L_B(u)&=-L_A(v) \end{align*} $$
where $A$ and $B$ are any linearly independant smooth vector fields on $\mathbb{R}^2$. We can also derive integrability conditions but we obtain two complicate elliptic PDEs and it is not trivial that a solution exists.
Question: Is there references treating this equations? and giving conditions for the existence of $C^{\infty}$ smooth solutions ?