The Einstein equations can be written as (here I am following the notation of Wald, Robert, General Relativity, Chicago, Chicago University Press, 1985)
\begin{equation} \partial_{\alpha}\Gamma^{\alpha}_{\mu\nu} - \partial_{\mu}\Gamma^{\alpha}_{\nu\alpha} + \Gamma^{\alpha}_{\mu\nu}\Gamma^{\beta}_{\alpha\beta} - \Gamma^{\alpha}_{\nu\beta}\Gamma^{\beta}_{\alpha\mu} = T_{\mu\nu} - \frac{1}{2}g_{\mu\nu}g^{\alpha\beta}T_{\alpha\beta} \end{equation}
(the LHS of this equation is the Ricci tensor; e.g. Eq. (3.4.5) in Wald's book. The RHS is the trace reversed stress-energy tensor).
We can expand out the Christoffel symbols in terms of the metric $g_{\mu\nu}$; then the Einstein equations become a set of second order partial differential equations for the metric. With a good choice of coordinates (such as harmonic coordinates, $\Gamma^{\alpha}_{\mu\nu}g^{\mu\nu}=0$), one can show that the Einstein equations can be thought of as a set of second order hyperbolic (wave) equations for each metric component (e.g. Eq. 10.2.33 in Wald's book).
My question is: can one think of the Einstein equations as a set of transport equations for the Christoffel symbols? Has there been any work that has looked at the initial value problem for the Einstein equations as I wrote them above, as a set of partial differential equations for the Christoffel symbols?
EDIT: To be more precise, I'd like to know if there has been work on the initial value formulation of the Einstein equations as a set of transport equations for the Christoffel symbols. For example, something akin to the Newman-Penrose formalism (in that formalism one rewrites the Einstein equations as a set of transport equations for "spin coefficients"; I'd like something similar for the Christoffel symbols).