Covariance matching constrained kriging is a technique originally developed for estimating linear or nonlinear functional of the variable of interest. Practically this may be regarded as a hybrid system that considers local accuracy property in kriging and preservation of spatial variability in stochastic simulation. The method estimates the unknown vector by adding the constraint to match the variance-covariance matrix of estimated values with the variance-covariance matrix of actual values under unbiasedness constraint. This study applies the method to a certain section of a lignite deposit which is subjected to severe tectonic activity. Lower heating value, ash content, and moisture content are estimated by covariance matching constrained kriging and ordinary kriging and also simulated by sequential Gaussian simulation. Comparisons between these methods show that covariance matching constrained kriging reproduces the spatial variability better than ordinary kriging and gives more accurate estimates than conditional simulation. Covariance matching constrained kriging can be used when spatial variability of the estimations is important as well as local accuracy. (C) 2012 Elsevier B.V. All rights reserved.