Finding all solutions of qKZ equations in characteristic $p$

In [MV] the difference qKZ equations were considered modulo a prime number $p$ and a family of polynomial solutions of the qKZ equations modulo $p$ was constructed by an elementary procedure as suitable $p$-approximations of the hypergeometric integrals. In this paper, we study in detail the first family of nontrivial example of the qKZ equations in characteristic $p$. We describe all solutions of these qKZ equations in characteristic $p$ by demonstrating that they all stem from the $p$-hypergeometric solutions. We also prove a Lagrangian property (called the orthogonality property) of the subbundle of the qKZ bundle spanned by the $p$-hypergeometric sections. This paper extends the results of [VV1] on the differential KZ equations to the difference qKZ equations.