EDIT: Here are answers to question 1 when l=9 and l=11. (As I explained in a comment the genus is 3 when l=7. It now appears that it's 10 when l=9 and 26 when l=11). Remarkably when l=3,5,7,9, or 11 the genus is the same as the genus of the compactification of the quotient of the upper half-plane by the principal congruence group, Gamma(l). I doubt that this is a coincidence, and am interested in what experts in the theory of characteristic p modular forms have to say.
Suppose first l=9. Extend the constant field from Z/2 to its algebraic closure,K. Let C in affine 4-space be the zero-locus of P, and L/K be the function field of C. P is generated by the "quintic relations" together with ab^2+bc^2+ca^2+d+d^2+d^3, where a,b,c,d are the coordinate functions u1,u2,u4 and u3. It follows that P is stabilized by the linear automorphisms (a,b,d,c)-->(b,c,d,a) and (a,b,d,c)-->(ua,ub,d,uc) with u^3=1. These automorphisms generate an order 9 group, G, which acts on L; let L_0 be the fixed field. It can be shown that L_0 is generated over K by abc and d and that (abc)^3=d^7+d^8+d^9. So L_0/K has genus 1. We now use Riemann-Hurwitz to calculate the genus, g, of L/K. (Since G has odd order, L/L_0 is tamely ramified).
The quintic relations all vanish on the line a=b=c=0. It follows that C has 3 points on this line; they are (0,0,d,0) with d+d^2+d^3=0. Each of these points is an ordinary triple point, and G permutes the branches at each of these points in a size 3 orbit. All the other orbits of G acting on the places of the function field L/K (including the places at infinity) are of size 9. Riemann-Hurwitz now tells us that 2g-2=9(2-2)+(9-3)+(9-3)+(9-3), so that g=10.
When l=11, one can argue in like manner. Now P is generated by the quintic relations, and the similar group G, acting on L/K, has order 55. I think one can again show that the genus of L_0/K is 1; this is the one thing I haven't checked completely. Now C sits in affine 5-space, the origin is an ordinary singular point of multiplicity 5, and G permutes the branches at the origin in a size 5 orbit. All other orbits of G acting on the places of L/K are of size 55 and Riemann-Hurwitz tells us that 2g-2=55(2-2)+(55-5), so that g=26.