[step:Use each orthonormal frame to trivialize the unit sphere bundle]
Define the unit sphere bundle $S(E) \subset E$ by
\begin{align*}
S(E):=\left\{v \in E:g_{\pi(v)}(v,v)=1\right\}.
\end{align*}
Let $\pi_S:S(E) \to B$ denote the restriction of the vector bundle projection $\pi:E \to B$ to $S(E)$, so $\pi_S(v)=\pi(v)$ for every $v \in S(E)$.
For an [open set](/page/Open%20Set) $U \subset B$ with smooth orthonormal frame $(e_1,\dots,e_k)$, define the map $\Phi_U: U \times S^{k-1} \to \pi_S^{-1}(U)$ by
\begin{align*}
\Phi_U\left(b,(a_1,\dots,a_k)\right)
=
\sum_{i=1}^{k} a_i e_i(b).
\end{align*}
Here $S^{k-1}$ denotes the standard sphere in $\mathbb{R}^k$, defined by
\begin{align*}
S^{k-1}:=\left\{a=(a_1,\dots,a_k)\in \mathbb{R}^k:\sum_{i=1}^{k}a_i^2=1\right\}.
\end{align*}
For $a \in S^{k-1}$, orthonormality gives
\begin{align*}
g_b\left(\sum_{i=1}^{k} a_i e_i(b),\sum_{j=1}^{k} a_j e_j(b)\right)
=
\sum_{i=1}^{k} a_i^2
=
1,
\end{align*}
so $\Phi_U$ indeed maps into $\pi_S^{-1}(U)$.
Define the map $\Psi_U:\pi_S^{-1}(U) \to U \times S^{k-1}$ by
\begin{align*}
\Psi_U(v)
=
\left(\pi(v),\left(g_{\pi(v)}(v,e_1(\pi(v))),\dots,g_{\pi(v)}(v,e_k(\pi(v)))\right)\right).
\end{align*}
If $v \in S(E)$ and $b=\pi(v)$, then the orthonormal basis expansion in $E_b$ gives
\begin{align*}
v=\sum_{i=1}^{k} g_b(v,e_i(b))e_i(b),
\end{align*}
and the coefficient vector belongs to $S^{k-1}$ because $g_b(v,v)=1$. Therefore $\Psi_U$ is well-defined. The maps $\Phi_U$ and $\Psi_U$ are inverse to each other by the uniqueness of coordinates in the orthonormal basis $(e_1(b),\dots,e_k(b))$.
The orthonormal-frame vector bundle trivialization $\Theta_U:E|_U \to U \times \mathbb{R}^k$ is the smooth map
\begin{align*}
\Theta_U(v)
=
\left(\pi(v),\left(g_{\pi(v)}(v,e_1(\pi(v))),\dots,g_{\pi(v)}(v,e_k(\pi(v)))\right)\right).
\end{align*}
Its inverse is the smooth map sending $(b,(a_1,\dots,a_k))$ to $\sum_{i=1}^k a_i e_i(b)$. If $v \in E_b$ has coordinates $a_i=g_b(v,e_i(b))$ in the orthonormal basis $(e_1(b),\dots,e_k(b))$, then
\begin{align*}
g_b(v,v)=\sum_{i=1}^{k}a_i^2.
\end{align*}
Thus $\Theta_U$ restricts exactly to a smooth identification of $S(E)|_U$ with the embedded submanifold $U \times S^{k-1}$ of $U \times \mathbb{R}^k$, because the fibrewise norm becomes the standard Euclidean norm in these coordinates. With this restricted smooth structure, both maps are smooth: $\Phi_U$ is built from smooth sections, smooth scalar multiplication, and smooth fibrewise addition, while $\Psi_U$ is built from $\pi$, the smooth sections $e_i$, and the smooth metric $g$. Hence $\Phi_U$ is a smooth local trivialization of $\pi_S:S(E)\to B$ over $U$ with fibre $S^{k-1}$.
[/step]