\subsection{Supplementary Material}


\subsection{INDUCING OSCILLATIONS}

Starting with the rf-dressed state described in the main text, we ramped the Raman beams from 0 to $\Omega_{\rm R}\!>\!4$ $E_{L}/\hbar$ in 70 ms. Then we ramped the rf coupling strength to an adjustable final value $\Omega_{\rm rf}$ in 2 ms, such that an effective Zeeman lattice was created. We induce sloshing by applying a synthetic electric field~\cite{Lin2011_Efield}, achieved by ramping the $z$-component of the Zeeman field to $\hbar\Omega_z\approx 2 E_{L}$ and then jumping it back to $0$ in 2~ms.

\subsection{MOMENTUM REPRESENTATION AND DIAGONALIZATION OF $H_{\rm{rf+R}}(x)$}

In order to compute the properties of our lattice potential we numerically diagonalized the combined rf-Raman Hamiltonian $H_{\rm rf+R}$ given the experimental parameters $\Omega_{\rm rf}, \Omega_{\rm R}, \Omega_z$. For simplicity we work in the momentum space representation of $H_{\rm rf+R}$. This representation offers an alternative understanding of the effective Zeeman lattice structure that arises from the combination of rf and Raman coupling fields.

While Raman transitions couple together states with $m_F$ differing by $\pm1$ and momentum differing by $\pm2\hbar k_{L}$, rf-coupling processes only change $m_F$ by $\pm1$, leaving the momentum unchanged. The combination of both Raman and rf fields sets appropriate conditions for secondary rf and Raman processes to populate states with higher momentum. The available states under rf-Raman coupling constitute a set of spin-momentum states $\{|\Psi_n\rangle\} =\{|m_F,\hbar (k_x \!+\! 2nk_L)\rangle \}$ where $n\!\in\!\mathbb{Z}; m_F\!=\!0,\pm1$. As expected, this basis is that of a lattice.

In the basis $\{|\Psi_n\rangle\}$, $H_{\rm rf+R}$ is a Hermitian block matrix of size $3(2N\!+\!1)$, when $n$ is restricted to $n\leq N$. For our parameters, dimensions larger than $3(2N\!+\!1)\!=81$ provided indistinguishable results. The dimension of this basis is appropriate for our calculations since we observed the population of states with up to $|n|\!=\!4$. For clarity, we construct the Hamiltonian by arranging the spin-momentum states in the $\{\Psi_n\}$ basis with increasing momentum index $n\!=\!-N,\! -N\!+\!1, \ldots,N$.
Along the principal diagonal, we have $3\times 3$ blocks
\begin{equation*}\label{rfRaman_Ham_diagonal}
    \mathbb{A}_{k_x}(n)\!=\!\hbar^2(k_x\!+\!2nk_L)^2 \mathbb{I} + \frac{\Omega_{\rm rf}}{\sqrt{2}}{F}_{x}-[\Omega_z {F}_{z}+\frac{\epsilon}{\hbar} (\hbar^2 \mathbb{I} \!-\! {F}_z^2)];
\end{equation*}
these terms correspond to kinetic energy, rf coupling of spin states with equal momentum, and the real Zeeman interaction, respectively; $\{F_x, F_y, F_z\}$ are the matrix representations of the $F\!=\!1$ angular momentum operators, and $\mathbb{I}$ is the $3\!\times\!3$ identity matrix.
Above and below the $\mathbb{A}_{k_x}(n)$ blocks, we have $3\!\times\!3$ blocks $\mathbb{B}\!=\!\sqrt{2}\Omega_{\rm R}({F}_{x}\!-\!i{F}_{y})/4$ describing the Raman coupling of spin-momentum states differing in momentum by $\Delta k_x\!=\!\pm 2k_{L}$.

The diagonalization of $H_{\rm rf+R}(n)$ as a function of $k_x$ gives the bandstructure of the combined rf-Raman lattice potential, $E(k_x)\!\!=\!\! E(k_x, \Omega_{\rm rf}, \Omega_{\rm R}, \Omega_z, \phi)$.

\subsection{EFFECTIVE ZEEMAN LATTICE PROPERTIES}

We extract the Zeeman lattice properties by fully characterizing the lowest energy band which in the tight binding regime is of the form $E(k_x)\!\!=\!\!-2t\cos(\pi k_x/k_{L}\!\!-\!\!\phi)$.
The tunneling amplitude is given by $t=\Delta E/4$, where $\Delta E$ is the width of the lowest band.
The lattice depth was calculated as the depth that a lattice potential would have in order to produce the calculated width $\Delta E$.
The effective mass is defined as $m^*\!=\!\hbar^2 [{\rm d}^2 E(k_x)/{\rm d}k_x^2]^{-1}$, where the derivative is evaluated at the point of interest (e.g. for our $\Omega_z=0$ measurements of effective mass displayed in Fig.~2, this was at $k_x=0$).
In the tight-binding regime, $t$ and $m^*/m$ are inversely proportional to each other $\pi^2 t/E_{L} \!=\! m/m^*$.
We obtained the Peierls phase $\phi$ by computing the shifts the lowest band [with minimum at  $k_{\min}\!= \!(\phi/\pi) k_{L}$] as a function of the experimental parameters.

\subsection{DELOADING}

We deloaded to a single bare spin state by rapidly ramping to $\Omega_z\!=\!0$, taking $\hbar\Omega_{\rm R}\!\rightarrow\!0$ to zero in 500~$\mu$s while increasing $\hbar\Omega_{\rm rf} \!\rightarrow\! 3E_{L}$. We then ramped the $\hbar\Omega_z\!\rightarrow\! -140E_{L}$ transferring all atoms into $|m_F\!=\!+1\rangle$.

\subsection{2D LATTICE FOR $1/3$ FLUX\\ HOFSTADTER MODEL}

We numerically solve the bandstructure of the 2D lattice arising from the combination of a period~$a\!=\!266$~nm lattice with depth $V_0=22.5 E_L$, and a~1D ``Zeeman lattice" with $\Omega_{\rm rf}\!=\!1 E_L$, $\Omega_{\rm R}\!=\!10 E_L$, $\epsilon\!=\! 0.44 E_L$. Here~$E_L$~denotes the Raman recoil, not the short-period lattice~recoil.
\end{document} 