Atomic-scale exchange between two different cations of similar size in crystalline oxides is one of the major types of point defects when multiple cations in oxygen interstitials are arrayed in an ordered manner. Although a number of studies have been performed on a variety of Li-intercalation olivine phosphates to determine the distribution of exchange defects in bulk, understanding of the thermodynamic stability of the defects in subsurface regions and its dependency on the crystallographic orientation at the surface has remained elusive. Through a combination of small-angle neutron scattering, atomic-scale direct probing with scanning transmission electron microscopy, and theoretical <i>ab initio</i> calculations, we directly demonstrate that the antisite exchange defects are distributed in a highly anisotropic manner near the surfaces of LiFePO<sub>4</sub> crystals. Moreover, a substantial amount of cation exchanges between Li and Fe sites is identified as an energetically favorable configuration in some surface regions, showing excellent agreement with the calculation results of negative defect formation energies. The findings in this study provide insight into developing better ways to avoid degradation of lithium mobility through the surface as well as scientifically notable features regarding the distribution of exchange defects in olivine phosphates