Congruence formulae for Legendre modular polynomials

Let $p\geq 5$ be a prime number. We generalize the results of E. de Shalit about supersingular $j$-invariants in characteristic $p$. We consider supersingular elliptic curves with a basis of $2$-torsion over $\overline{\mathbf{F}}_p$, or equivalently supersingular Legendre $\lambda$-invariants. Let $F_p(X,Y) \in \mathbf{Z}[X,Y]$ be the $p$-th modular polynomial for $\lambda$-invariants. A simple generalization of Kronecker's classical congruence shows that $R(X):=\frac{F_p(X,X^{p})}{p}$ is in $\mathbf{Z}[X]$. We give a formula for $R(\lambda)$ if $\lambda$ is a supersingular. This formula is related to the Manin--Drinfeld pairing used in the $p$-adic uniformization of the modular curve $X(\Gamma_0(p)\cap \Gamma(2))$. This pairing was computed explicitly modulo principal units in a previous work of both authors. Furthermore, if $\lambda$ is supersingular and lives in $\mathbf{F}_p$, then we also express $R(\lambda)$ in terms of a CM lift (which are showed to exist) of the Legendre elliptic curve associated to $\lambda$.