The temperature-dependent electric field gradient (EFG) tensor at a probe nucleus reflects the combined effects of local symmetry, bonding, and lattice dynamics. However, determining the primary atomic distortions responsible for this response can be difficult. Here, we combine electron-gamma (e-γ) time-differential perturbed angular correlation (TDPAC) measurements on probes produced by the decay at substitutional Ti sites in rutile with symmetry-resolved local frozen-phonon density functional theory (LFP-DFT) calculations. Atomic displacements of the octahedron are decomposed into site-symmetry-adapted local modes, and their contributions to the EFG are thermally weighted via Boltzmann averaging. We find that, in pristine rutile , the temperature evolution of the EFG shape is primarily governed by two orthogonal equatorial-oxygen shear distortions. For Ta substituting the Ti site, the same shear sector remains active but preserves a finite in-plane anisotropy, resulting in a comparatively weak temperature dependence of the asymmetry parameter η. Measurements between 30 and 470 K, together with earlier high-temperature gamma-gamma data, show weak low-temperature curvature and high-temperature quasilinearity of the principal EFG component , while η(T) remains comparatively robust. The calculations reproduce these qualitative trends and provide a microscopic, symmetry-based interpretation of the thermal evolution of hyperfine tensors in solids.