We let $M(n)$ be the $3\times 3$ submatrix of $S$ which consists of the rows $3n-1,3n,3n+1$ of $S$, and let $T_0(n),T_1(n),T_2(n)$ be the $3 \times 3$ matrices \begin{align*} T_0(n):=\begin{pmatrix} 8n-4 & 8n-3 & 8n-2\\ 8n-2 & 8n-1 & 8n\\ 8n+1 & 8n+2 & 8n+3\\ \end{pmatrix},\;\; T_1(n):=\begin{pmatrix} 8n-4 & 8n-3 & 8n-3\\ 8n-2 & 8n-1 & 8n\\ 8n+1 & 8n+2 & 8n+3\\ \end{pmatrix}\,\;\end{align*}\begin{align*} T_0(n):=\begin{pmatrix} 8n-4 & 8n-3 & 8n-2\\ 8n-2 & 8n-1 & 8n\\ 8n+1 & 8n+2 & 8n+3\\ \end{pmatrix},\;\; T_2(n):=\begin{pmatrix} 8n-4 & 8n-3 & 8n-3\\ 8n-2 & 8n-1 & 8n\\ 8n+1 & 8n+2 & 8n+3\\ \end{pmatrix}\,\;\end{align*} \begin{align*} T_2(n):=\begin{pmatrix} 8n-4 & 8n-3 & 8n-2\\ 8n-1 & 8n-1 & 8n\\ 8n+1 & 8n+2 & 8n+3\\ \end{pmatrix}\;\;. \end{align*}\begin{align*} T_1(n):=\begin{pmatrix} 8n-4 & 8n-3 & 8n-2\\ 8n-1 & 8n-1 & 8n\\ 8n+1 & 8n+2 & 8n+3\\ \end{pmatrix}\;\;. \end{align*}
By Cobham's theorem ("a sequence arising as (image under a coding of) a fixed point of a $k$-uniform morphism is $k$-automatic") the sequence $(t_n)$ is therefore $k$-automatic.
By the closure properties of the set of $k-$automatic sequences (under shifts, parallel generation and taking functions of output elements) then also the sequence $(\Delta_n)$ with
$$\Delta_n=(k^2-1)+t_{n+1}-t_n=x_{n+1}-x_n$$
is $k$-automatic. (That is, $\Delta_n$ can be computed from the base-$k$
digits of $n$ with a finite state machine of the type desribed above.)