Recently, a new sampling formalism for Large Eddy simulation was proposed by Winckelmans et al. [1] and Knaepen et al.[2], which is a projection method for Navier-Stokes equations from continuum space to a discrete space, using a sampling operator instead of a filter operator. Since the sampling operator is not commutative with spatial derivatives, a closure term appears which represents the loss of information due to the projection on a discrete mesh. In e.g.[2] a Smagorinsky model was proposed that, by relying on a generalized dynamic procedure, succeeded in accounting for the subgrid scales. In this paper, we investigate the ability of this sampling based dynamic procedure, in combination with an appropriate model for the truncation error, to obtain higher-order numerical accuracy. Two such possible models are presented. Further, we show that Richardson extrapolation is a simplified formulation of this procedure.