A simultaneous space-time variational formulation of a parabolic evolution problem is solved with an adaptive wavelet method. This method is shown to converge with the best possible rate in linear complexity. Thanks to the use of tensor product bases, there is no penalty in complexity due to the additional time dimension. Special wavelets are designed such that the bi-infinite system matrix is sparse. This sparsity largely simplifies the implementation and improves the quantitative properties of the adaptive wavelet method. Numerical results for an ODE and the heat equation are presented. AMS subject classifications. 35K15, 41A25, 42C40, 65F50, 65T60. Key words. Adaptive wavelet scheme, parabolic equations, simultaneous space-time variational formulation, tensor product approximation, optimal computational complexity.