Yes it has dimension $1$ for $p \ge 5$.
This module is actually the natural module for ${\rm SO}(3,p)$, which is isomorphic to ${\rm PGL}(2,p)$ (the scalars of ${\rm GL}(2,p)$ act trivially), so the result is probably in the literature somewhere.
But it's easy to do it by direct calculation. A Sylow $p$-subgroup $P$ of $G = {\rm GL}(2,p)$, which you can take to be the upper unitriangular matrices, has order $p$, and acts indecomposably on the module, so it fixes a submodule of dimension $1$.
So in an extension by the module $M_p$, the $p$-th power of a generator of $P$ must lie in this fixed submodule, and hence $H^2(G,M_p)$ and so also $H^2(G,M_p)$ have dimension at most $1$. Now for $p \ge 5$, it is not hard to see that taking this $p$-th power to be nontrivial defines a nonsplit extension of $P$ of $M_p$. (It is a split extension when $p=3$, so $H^2(P,M_3) = H^2(G,M_3) = 0$.)
Since distinct $p$-Sylow subgroups have trivial intersection, to check that this element of $H^2(P,M_p)$ is stable in $G$, we just need to check that it is stable in $N = N_G(P)$, which is the group of upper triangular nonsingular matrices. Checking that the actions of $N$ on $M_p$ and $P$ extend to give an action of $N$ on the nonsplit extension is a routine calculation. I could give details if you like.