The Jordan-Wigner transformation is a powerful approach to studying one-dimensional spin models. The following dictionary between spin operators and creation/annihilation operators for fermions allows us to recast spin models in terms of fermions:
$$S_{j}^{+} = c_{j}^{\dagger}\exp\left[i\pi\sum_{k<j}n_{k}\right],\hspace{12pt} S_{j}^{z} = n_{j}-\frac{1}{2}$$
For example, a Heisenberg spin-1/2 chain in a magnetic field becomes
$$H = J\sum_{j}\vec{S}_{j}\cdot\vec{S}_{j+1} - h\sum_{j}S_{j}^{z} \longrightarrow \sum_{j}\left[\frac{J}{2}\left(c_{j}^{\dagger}c_{j+1}+c_{j+1}^{\dagger}c_{j}\right) + J\left(n_{j}-\frac{1}{2}\right)\left(n_{j+1}-\frac{1}{2}\right) -hn_{j}\right].$$
The original spin model has a $U(1)$-symmetry associated with rotations about the magnetic field axis; in the fermion model, this is translated into conservation of total fermion number. There are several consequences of this for, e.g., quantum critical points induced by an applied magnetic field. Moreover, this is why the fermion theory is local in spite of the apparent non-locality in the Jordan-Wigner dictionary.
My question is, what can we say about such a model in which the $U(1)$ symmetry is broken from the beginning?
There are many choices which can accomplish this, but one simple choice is to include a Heisenberg anisotropy. This destroys our freedom to choose $h\propto\hat{z}$, so the spin model will look something like
$$H = J\sum_{j}\left(\vec{S}_{j}\cdot\vec{S}_{j+1} + \lambda S_{j}^{y}S_{j+1}^{y}\right) - \sum_{j}\vec{h}\cdot\vec{S}_{j}$$
The $U(1)$ symmetry is certainly broken now; moreover, if $h$ has components in the x or y directions, the corresponding fermion model appears to have a nonlocal structure, since a single operator $S^{\pm}$ gets turned into a long-range string of fermion operators. Are there worked out examinations of such models in the literature? All examples of the Jordan-Wigner mapping I have seen feature number conservation and spatial locality, but it seems to me violations of these constraints are very common.