[step:Apply convex lower semicontinuity to the lifted minor variables]
Define the lifted sequence
\begin{align*}
Z_j: U \to \prod_{k\in S}\mathbb R^{N_k}, \qquad Z_j(x) = M_S(\nabla u_j(x))
\end{align*}
and define the lifted limit
\begin{align*}
Z:U\to \prod_{k\in S}\mathbb R^{N_k},\qquad
Z(x)=M_S(\nabla u(x)).
\end{align*}
The preceding steps show that $u_j\to u$ in measure on $U$ and that every component block of $Z_j$ converges weakly to the corresponding component block of $Z$ in the natural space $L^{p/k}(U;\mathbb R^{N_k})$ when $k<p$ and in $L^1(U;\mathbb R^{N_k})$ when $k=p$. For $k<p$, set $r_k=p/k>1$. Since $\mathcal L^n(U)<\infty$, every essentially bounded test function belongs to $L^{r_k'}(U;\mathbb R^{N_k})$, where $r_k'=r_k/(r_k-1)$ is the Hölder conjugate exponent. Thus weak convergence in $L^{r_k}$ implies convergence against all $L^\infty$ test functions; equivalently, these subcritical blocks converge in the $\sigma(L^1,L^\infty)$ testing topology used for the finite-dimensional lifted minor variable.
For each $j$, the definition of $W$ gives
\begin{align*}
I[u_j]=\int_U G(x,u_j(x),Z_j(x))\,d\mathcal L^n(x).
\end{align*}
Similarly,
\begin{align*}
I[u]=\int_U G(x,u(x),Z(x))\,d\mathcal L^n(x).
\end{align*}
The integrand $G$ is a normal integrand, is lower semicontinuous in $(y,z)$ for a.e. $x$, and is convex in $z$ for each fixed $y$. Moreover, the lower bound
\begin{align*}
G(x,u_j(x),Z_j(x))=W(x,u_j(x),\nabla u_j(x))\ge c|\nabla u_j(x)|^p-a(x)
\end{align*}
provides integrable control of the negative part along the lifted admissible sequence, since $a\in L^1(U)$ and $(\nabla u_j)$ is bounded in $L^p(U;\mathbb R^{m\times n})$.
We apply the Ioffe lower semicontinuity theorem for convex normal integrands in finite-dimensional convex variables. Its hypotheses are met as follows. The function $G$ is a normal integrand by assumption, hence it is measurable in $x$ and lower semicontinuous in $(y,z)$ for a.e. $x$. For each fixed $y\in\mathbb R^m$, the map $z\mapsto G(x,y,z)$ is convex for a.e. $x$. The sequence $u_j$ converges to $u$ in measure. Each block of $Z_j$ converges to the corresponding block of $Z$ either weakly in $L^1$ in the critical case or against all $L^\infty$ tests in the subcritical case, and the product target $\prod_{k\in S}\mathbb R^{N_k}$ is finite-dimensional because $S$ is finite. Finally, the displayed lower bound supplies an $L^1$ lower control of the negative part along the sequence through the integrable function $a$ and the bounded $L^p$ norms of $\nabla u_j$. Therefore the theorem gives
\begin{align*}
\int_U G(x,u(x),Z(x))\,d\mathcal L^n(x)\le \liminf_{j\to\infty}\int_U G(x,u_j(x),Z_j(x))\,d\mathcal L^n(x).
\end{align*}
Substituting the definitions of $Z_j$, $Z$, and $W$, we obtain
\begin{align*}
I[u]\le \liminf_{j\to\infty} I[u_j]=L.
\end{align*}
[/step]