A wide range of technological applications of lithium tantalate (LT) is closely related to the defect chemistry. In literature, several intrinsic defect models have been proposed. Here, using a combinational approach based on DFT and solid-state NMR, we demonstrate that distribution of electric field gradients (EFGs) can be employed as a fingerprint of a specific defect configuration. Analyzing the distribution of 7 Li EFGs, the FT-IR and electron spin resonance (ESR) spectra, and the 7 Li spin-lattice relaxation behavior, we have found that the congruent LT samples provided by two manufacturers show rather different defect concentrations and distributions although both were grown by the Czochralski method. After thermal treatment hydrogen out-diffusion and homogeneous distribution of other defects have been observed by ESR, NMR, and FT-IR. The defect structure in one of two congruent LT crystals after annealing has been identified and proved by defect formation energy considerations, whereas the more complex defect configuration, including the presence of extrinsic defects, has been suggested for the other LT sample. The approach of searching the EFG fingerprints from DFT calculations in NMR spectra can be applied for identifying the defect clusters in other complex oxides.