[step:Uniquely determine $X_t$ from the volume-form isomorphism $\iota_\bullet \omega_t: \Gamma(TM) \to \Omega^{n-1}(M)$]
With $\alpha \in \Omega^{n-1}(M)$ such that $d\alpha = \omega_0 - \omega_1$, we set $\beta := -\alpha$ so $d\beta = \omega_1 - \omega_0$. The Moser equation reduces to the *pointwise linear* equation
\begin{align*}
\iota_{X_t}\omega_t = \beta \quad \text{for all } t \in [0,1].
\end{align*}
*Claim.* For each $t \in [0,1]$ and each $p \in M$, the map $T_pM \to \bigwedge^{n-1}T^*_pM$, $v \mapsto \iota_v(\omega_t|_p)$, is a linear isomorphism.
*Proof of claim.* Both vector spaces are $n$-dimensional (since $\dim T_pM = n$ and $\dim \bigwedge^{n-1}T^*_pM = \binom{n}{n-1} = n$), so it suffices to check injectivity. Suppose $\iota_v(\omega_t|_p) = 0$ for some $v \in T_pM$. For any $w \in T_pM$, pair with $w$:
\begin{align*}
\iota_v(\omega_t|_p)(w, \ast) = \omega_t|_p(v, w, \ast),
\end{align*}
in particular, for any basis $w_2, \dots, w_n$ of a complement of $\operatorname{span}(v)$,
\begin{align*}
\omega_t|_p(v, w_2, \dots, w_n) = \iota_v(\omega_t|_p)(w_2, \dots, w_n) = 0.
\end{align*}
If $v \ne 0$, we could extend $v, w_2, \dots, w_n$ to a basis and $\omega_t|_p$ would vanish on this basis, contradicting the non-vanishing of the volume form $\omega_t|_p$. Hence $v = 0$, so the map is injective, hence bijective.
*Consequence.* For each $t$ and $p$, there exists a unique $X_t(p) \in T_pM$ with $\iota_{X_t(p)}(\omega_t|_p) = \beta|_p$. This defines a vector field $X_t$ on $M$, which is smooth jointly in $(p, t)$ because the assignment is inverse to a smooth linear isomorphism depending smoothly on parameters.
[/step]