We extend and analyze the VMS method to partial differential equations with stochastic coefficients. For a natural choice of an "optimal" coarse-scale solution and L2-orthogonal stochastic basis functions, we demonstrate ...