We derive a posteriori error estimates for fully discrete approximations to solutions of linear parabolic equations. The space discretization uses finite element spaces that are allowed to change in time. Our main tool is an appropriate adaptation of the elliptic reconstruction technique, introduced by Makridakis and Nochetto. We derive novel a posteriori estimates for the norms of L8(0, T; L2(O)) and the higher order spaces, L8(0, T;H1(O)) and H1(0, T; L2(O)), with optimal orders of convergence.