Let $E$ be a vector bundle over a smooth Riemannian manifold $(M,g)$. Suppose $E$ is equipped with a metric $\eta$ and a metric connection $\nabla$.
Denote by $\Omega^k(M,E)$ the space of $E$-valued forms of degree $k$. $\nabla$ induces a covariant exterior derivative $d:\Omega^k(M,E) \to \Omega^{k+1}(M,E)$, which (together with the metrics $g,\eta$) induce the codifferential operator:
$\delta:\Omega^k(M,E) \to \Omega^{k-1}(M,E)$.
Even though the metric $\eta$ on $E$ plays a part in the definition of $\delta$, it turns out that $\delta$ is actually independent of it, and depends only on the connection $\nabla$ (there is an implicit dependence since $\nabla$ is required to be metric).
Indeed, one formula for $\delta$ is $\delta =\star d \star$ (up to sign), where $\star$ is the Hodge-dual $\Omega^k(M,E) \to \Omega^{d-k}(M,E)$ which is defined without any reference to the metric on $E$.
Is there any way to see why the codifferential depends only on the connection $\nabla$ and the metric on $M$, and not on the metric on $E$?
Of course, the derivation of the formula above for $\delta$ shows this, but I would like to find an argument without relying on this computation.
(I am looking for a more "conceptual explanation").
Note that it's easy to see from the definition of $\delta$ that it's invariant under scaling of the metric on $E$, however I do not see immediately why it's completely independent of it.
Assume we have two metrics $\eta_1 = \left< \cdot, \cdot \right>_1$ and $\eta_2 = \left< \cdot, \cdot \right>_2$ on $E$ that are compatible with $\nabla$. Let $\Phi \colon E \rightarrow E$ be a vector bundle isomorphism such that
$$ \left< s, t \right>_1 = \left< \Phi(s), t \right>_2. $$
In other words, for each $p \in M$, the map $\Phi_p \colon E_p \rightarrow E_p$ is the unique linear map representing the bilinear form $\eta_1$ with respect to $\eta_2$. One way to see that $\Phi$ exists is to note that $\Phi$ is obtained from the tensor $\eta_1 \in \Gamma(E \otimes E)$ by raising an index with respect to the metric $\eta_2$. The basic observation is that $\Phi$ is $\nabla$-parallel because it is obtained from a parallel tensor ($\eta_1$) by raising an index with respect to a parallel metric $\eta_2$. Alternatively, we can compute
$$ \left< \nabla_X(\Phi)(s), t \right>_2 = \left< \nabla_X(\Phi(s)) - \Phi(\nabla_X s), t \right>_2 = \left< \nabla_X(\Phi(s)), t \right>_2 - \left< \nabla_X s, t \right>_1 \\ = \left( X \left< \Phi(s), t \right>_2 - \left< \Phi(s), \nabla_X t\right>_2 \right) - \left< \nabla_X s, t \right>_1 \\ = X \left< s, t \right>_1 - \left< s, \nabla_X t \right>_1 - \left< \nabla_X s, t \right>_1 = 0 $$
where we used both $\nabla \eta_1 = 0$ and $\nabla \eta_2 = 0$.
Using $\nabla \Phi = 0$ and the definition of the covariant exterior derivative, one can then show that $d_{\nabla}$ commutes with $\Phi$ in the sense that
$$ \Phi(d_{\nabla}(\alpha)(X_0, \dots, X_k)) = d_{\nabla}(\Phi(\alpha))(X_0, \dots, X_k) $$
where $\alpha, \Phi(\alpha) \in \Omega^{*}(M;E)$.
Finally, let us show that if $\delta$ is the adjoint of $d_{\nabla}$ with respect to $\eta_2$ then it is also the adjoint of $d_{\nabla}$ with respect to $\eta_1$ and so it is independent of the metric $\eta$ as long as it is $\nabla$-compatible. Set $\Phi_k := \operatorname{id}|_{\Lambda^k(T^{*}(M))} \otimes \Phi$ and (abusing notation a little) denote the induced metric from $g$ and $\eta_i$ on $\Lambda^k(T^{*}M) \otimes E$ by $\left< \cdot, \cdot \right>_i$. Then $\Phi_k$ satisfies
$$ \left < \alpha, \beta \right>_1 = \left< \Phi_k(\alpha), \beta \right>_2 $$
for $\alpha,\beta \in \Omega^k(M;E)$ and the statement that $d_{\nabla}$ commutes with $\Phi = \Phi_0$ can be written as $\Phi_{k+1} \circ d_{\nabla} = d_{\nabla} \circ \Phi_k$.
Using all the properties above, we have
$$ \int_{M} \left< d_\nabla(\alpha), \beta \right>_1 \, dV_g = \int_M \left< \Phi_{k+1}(d_{\nabla}(\alpha)), \beta \right>_2 \, dV_g = \int_M \left< d_{\nabla}(\Phi_{k}(\alpha)), \beta \right>_2 \, dV_g = \\ \int_M \left< \Phi_{k}(\alpha), \delta(\beta) \right>_2 \, dV_g = \int_M \left< \alpha, \delta(\beta) \right>_1 \, dV_g. $$