Anisotropic Gaussian random fields arise in probability theory and in various applications. Typical examples are fractional Brownian sheets, operator-scaling Gaussian fields with stationary increments, and the solution to the stochastic heat equation.This paper is concerned with sample path properties of anisotropic Gaussian random fields in general. Let X = {X(t), t ∈ R N } be a Gaussian random field with values in R d and with parameters H 1 , . . . , H N . Our goal is to characterize the anisotropic nature of X in terms of its parameters explicitly.Under some general conditions, we establish results on the modulus of continuity, small ball probabilities, fractal dimensions, hitting probabilities and local times of anisotropic Gaussian random fields. An important tool for our study is the various forms of strong local nondeterminism.