Given a closed Riemannian manifold $(X,g)$ and let $p\colon TX\to X$ be the usual projection, the paper I'm reading asserts that the Levi-Civita connection induces a splitting $T(TX)= H(TX)\oplus V(TX)\cong p^* TX\oplus p^*TX$, so that $$ \begin{aligned}\Omega^\bullet(TX) &= \{ \Gamma(H(TX))\to C^\infty(TX)\}\otimes_{C^\infty(TX)} \{\Gamma(V(TX))\to C^\infty(TX) \} \\ &\cong p^*\Omega^\bullet (X)\otimes_{C^\infty(TX)} p^*\Omega^\bullet(X). \end{aligned} $$

I want to investigate this isomorphism explicitly, in the sense that I can compute the images of the local coordinate functions on $X$ under the two maps

$$ \phi_1\colon C^\infty(X)\xrightarrow{d}\Omega^1(X)\xrightarrow{p^*} p^*\Omega^\bullet(X) \cong \{\Gamma(V(TX))\to C^\infty(TX) \} \hookrightarrow \Omega^\bullet(TX),$$

$$\phi_2\colon C^\infty(X)\xrightarrow{d}\Omega^1(X)\xrightarrow{p^*} p^*\Omega^\bullet(X) \cong \{\Gamma(H(TX))\to C^\infty(TX) \} \hookrightarrow \Omega^\bullet(TX). $$

But I failed finding any material that could tell me how the horizontal vector bundle is splitted out using the Levi-Civita connection, without which I cannot do any computation.

Also, I found that $V(TX) = \ker ( d p\colon T(TX)\to TX)$ on Wikipedia, but I'm not quite sure how the isomorphism $V(TX)\cong p^*TX$ is given. The isomorphism that I came up with is the following:

Let $x_1,\cdots,x_n$ be the coordinate functions of a coordinate neighborhood $U$ of $X$, then $x_1,\cdots,x_n,y_1:=d x_1,\cdots, y_n:=dx_n$ are the coordinates of $p^{-1}(U)\cong U\times\mathbb{R}^n$ of $TX$. Then the fibres of $V(TX)|_{p^{-1}(U)}$ are exactly the spans of $\frac{\partial}{\partial y_i}$'s. Since the fibres of $p^*TX|_{p^{-1}(U)}$ are the spans of $\frac{\partial}{\partial x_i}$'s, the isomorphism is given by $\frac{\partial}{\partial x_i}\mapsto \frac{\partial}{\partial y_i} $ along each fibre.

So my questions are:

1. How the horizontal vector bundle is constructed using the Levi-Civita connection?
2. Is the above isomorphism the correct one?
3. What are the images of the local coordinate functions on $X$ under the maps $\phi_1$ and $\phi_2$ given above? (Probably I will be able to do this by myself with (1) and (2) solved, but I'm putting it here anyway, since it is the ultimate goal of mine)

To be honest, I'm not familiar to Riemannian geometry, so it would be the best if there exists a self-contained material that could provide a throughout treatment to this matter.

Thanks in advance for any help.

EDIT: I thought it over and a new question (4) arises: See, a choice of coordinates $x_1,\cdots,x_n$ of a neighborhood $U$ of $X$ induces coordinates $x_1,\cdots,x_n,y_1:=d x_1,\cdots, y_n:=dx_n$ of $p^{-1}(U)$ of $TX$, which gives a basis of $T(TX)|_{p^{-1}(U)}$. While the span of $\frac{\partial }{\partial y_i}$’s gives $V(TX)| _{p^{-1}(U)} $, why isn’t $H(TX)|_{p^{-1}(U)}$ simply the span of $\frac{\partial}{\partial x_i}$’s? I don’t see any trouble gluing these choices of $H(TX)|_{p^{-1}(U)}$‘s up as $U$ varies, which seems to have no difference with $p^*TX$. Am I missing something in this construction, or is there any other reason to define $H(TX)$ differently? This approach looks canonical — not even a Riemannian metric is required.