Anisotropy of the quadratic Zeeman effect for the 3 H6 → 3 H4 transition at 793 nm wavelength in 169 Tm 3+ -doped Y3Al5O12 is studied, revealing shifts ranging from near zero up to + 4.69 GHz/T 2 for ions in magnetically inequivalent sites. This large range of shifts is used to spectrally resolve different subsets of ions and study nuclear spin relaxation as a function of temperature, magnetic field strength, and orientation in a site-selective manner. A rapid decrease in spin lifetime is found at large magnetic fields, revealing the weak contribution of direct-phonon absorption and emission to the nuclear spin-lattice relaxation rate. We furthermore confirm theoretical predictions for the phonon coupling strength, finding much smaller values than those estimated in the limited number of past studies of thulium in similar crystals. Finally, we observe a significant -and unexpectedmagnetic field dependence of the two-phonon Orbach spin relaxation process at higher field strengths, which we explain through changes in the electronic energy-level splitting arising from the quadratic Zeeman effect.