We present tomographic velocity and anisotropy models of the uppermost mantle beneath the Iran region. A total of 74,375 Pn phase readings from 86 stations of the Iranian network and 133 stations of the International Seismological Centre are used in the investigation. The study uses the Pn travel‐time tomography method proposed by Hearn (1996). The tomography results show some interesting anomalies. The average Pn velocity under the Iran region is approximately 8.0 km/s, and the maximum velocity perturbations are approximately 3%–4%. High Pn velocities under the Zagros Mountains and the Caspian Sea may be due to the presence of oceanic crust/lithosphere material. Low Pn velocities were found under the Alborz and Caucasus regions and may be due to higher temperatures or partial melting resulting from volcanoes and mid‐Cenozoic volcanic/plutonic rocks in these regions. The inversion velocities support the idea that the subduction of the Arabian plate into the mantle beneath the Iranian plateau may have resulted in the upwelling of hot material. The well‐resolved Pn anisotropy model is jointly obtained with a velocity model for the areas with good ray‐path coverage. In the plate collision regions (Zagros and Alborz), the fast Pn anisotropy direction is oriented parallel to the collision arc and to large reverse faults due to pure shear deformation from cross‐fault compression and along‐fault extension. Under the Caucasus regions, the Pn anisotropy results indicate that the preferred alignment of olivine crystals is parallel to the plate movement direction; however, the surface fault strike is at an angle of nearly 45° with the crustal movement direction and anisotropy. These differing deformations suggest potential decoupling between the crust and upper mantle. The possible decoupling and differing deformation between the crust and upper mantle are easily enhanced under high temperatures due to volcanoes and supported by low velocities beneath the Caucasus. We validate the existence of Pn anisotropy under these regions by azimuthal averaging of the apparent Pn velocity.