Model Reduction of Parabolic PDEs using Multivariate Splines