The multivariate Gaussian random fields with matrix-based scaling laws are widely used for inference in statistics and many applied
areas. In such contexts interests are often in symmetry and in the rates of change of spatial surfaces in any given direction. This article analyzes the almost sure sample function behavior for operator fractional Brownian motion, including multivariate fractional Brownian motion. We obtain the estimations of small ball probability and the strongly locally nondeterministic for operator fractional Brownian motion in any given direction. Applying these estimates we obtain Chung's laws of the iterated logarithm for spatial surfaces of operator fractional Brownian motion. Our results show that the precise rates of change of spatial surfaces are completely determined by the self-similarity exponent.