Suppose that ¯(x1),...,¯Z(xn). are observations of vector-valued random function ¯(x). In the isotropic situation, the sample variogram γ*(h) for a given lag h is $$\bar \gamma ^ * (h) = \frac{1}{{2N(h)}}\mathop \sum \limits_{s(h)} (\overline Z (x_1 ) - \overline Z (x_1 )) \overline {(Z} (x_1 ) - \overline Z (x_1 ))^T $$ where s(h) is a set of paired points with distance h and N(h) is the number of pairs in s(h)..