[guided]We now prove the pointwise minimality statement. The idea is Hilbert-space orthogonality: $v_t$ is the closed-gradient part of any other velocity field representing the same infinitesimal motion, so the remainder $w_t-v_t$ is orthogonal to the closed gradient subspace.
Let $w:[0,1]\times\mathbb R^n\to\mathbb R^n$ be a Borel vector field satisfying the hypotheses in the theorem. Thus $w_t\in L^2(\mu_t;\mathbb R^n)$ for a.e. $t$, the spacetime integral of $|w_t|^2$ is finite, and $(\mu_t,w_t)$ solves the same continuity equation as $(\mu_t,v_t)$.
Fix a spatial test function $\phi\in C_c^\infty(\mathbb R^n)$ and define $K_\phi:[0,1]\to\mathbb R$ by
\begin{align*}
K_\phi(t):=\int_{\mathbb R^n}\phi(x)\,d\mu_t(x).
\end{align*}
For every $\eta\in C_c^\infty((0,1))$, use the product test function $\zeta:[0,1]\times\mathbb R^n\to\mathbb R$ given by $\zeta(t,x):=\eta(t)\phi(x)$. Substituting this test into the continuity equation for $v$ gives the distributional identity
\begin{align*}
\int_0^1 K_\phi(t)\eta'(t)\,d\mathcal L^1(t)+\int_0^1\eta(t)\int_{\mathbb R^n}\nabla\phi(x)\cdot v_t(x)\,d\mu_t(x)\,d\mathcal L^1(t)=0.
\end{align*}
Substituting the same test into the continuity equation for $w$ gives
\begin{align*}
\int_0^1 K_\phi(t)\eta'(t)\,d\mathcal L^1(t)+\int_0^1\eta(t)\int_{\mathbb R^n}\nabla\phi(x)\cdot w_t(x)\,d\mu_t(x)\,d\mathcal L^1(t)=0.
\end{align*}
The first term is identical in the two identities because the measure curve is the same. Subtracting the two identities gives
\begin{align*}
\int_0^1\eta(t)\int_{\mathbb R^n}\nabla\phi(x)\cdot (w_t(x)-v_t(x))\,d\mu_t(x)\,d\mathcal L^1(t)=0.
\end{align*}
Since this holds for every $\eta\in C_c^\infty((0,1))$, the fundamental lemma for distributions gives
\begin{align*}
\int_{\mathbb R^n}\nabla\phi(x)\cdot (w_t(x)-v_t(x))\,d\mu_t(x)=0
\end{align*}
for $\mathcal L^1$-a.e. $t$, for this fixed $\phi$.
To pass from one fixed test to all gradients in $G_t$, we use the standard countable-dense-family argument. Choose a countable family of smooth compactly supported functions whose gradients are dense in the collection of smooth compactly supported gradients on each compact set after taking rational coefficients and rational mollifier parameters. Intersecting the corresponding full-measure sets gives a single full-measure set of times where the orthogonality identity holds for that countable family. For such a time $t$, the maps
\begin{align*}
g\mapsto \int_{\mathbb R^n}g(x)\cdot v_t(x)\,d\mu_t(x)
\end{align*}
and
\begin{align*}
g\mapsto \int_{\mathbb R^n}g(x)\cdot w_t(x)\,d\mu_t(x)
\end{align*}
are continuous linear functionals on $L^2(\mu_t;\mathbb R^n)$ by the [Cauchy-Schwarz inequality](/theorems/432). Therefore the orthogonality identity extends from the dense family to every $g\in G_t$:
\begin{align*}
\int_{\mathbb R^n}g(x)\cdot (w_t(x)-v_t(x))\,d\mu_t(x)=0.
\end{align*}
Now use the fact that $v_t\in G_t$. Taking $g=v_t$ gives
\begin{align*}
\int_{\mathbb R^n}v_t(x)\cdot (w_t(x)-v_t(x))\,d\mu_t(x)=0.
\end{align*}
Thus $v_t$ and $w_t-v_t$ are orthogonal in the Hilbert space $L^2(\mu_t;\mathbb R^n)$. Expanding the square using this orthogonality gives
\begin{align*}
\int_{\mathbb R^n}|w_t(x)|^2\,d\mu_t(x)=\int_{\mathbb R^n}|v_t(x)|^2\,d\mu_t(x)+\int_{\mathbb R^n}|w_t(x)-v_t(x)|^2\,d\mu_t(x).
\end{align*}
The second term on the right is nonnegative. Hence
\begin{align*}
\int_{\mathbb R^n}|v_t(x)|^2\,d\mu_t(x)\le \int_{\mathbb R^n}|w_t(x)|^2\,d\mu_t(x)
\end{align*}
for $\mathcal L^1$-a.e. $t$. This is exactly the asserted minimality of the constructed velocity.[/guided]