Abstract Using the data from Magnetospheric Multiscale mission, we statistically investigate the anisotropy of turbulence within bursty bulk flows (BBFs) in the Earth’s magnetotail. By computing five‐point second‐order structure functions, we reveal the three‐dimensional anisotropy and scaling properties of turbulence at ion scales. In the coordinate system (Lˆ⊥,lˆ⊥,lˆ‖ ${\widehat{L}}{\mathit{\perp }},{\widehat{l}}{\mathit{\perp }},{\widehat{l}}{\mathit{\Vert }}$), the characteristic lengths of turbulence eddies for both total and perpendicular magnetic field fluctuations exhibit l‖>L⊥>l⊥ ${l}{\mathit{\Vert }} > {L}{\mathit{\perp }} > {l}{\mathit{\perp }}$ at ion‐scales. The anisotropy scaling relations for the total magnetic field are l‖∝l⊥1 ${l}{\mathit{\Vert }}\mathit{\propto }{l}{\mathit{\perp }}^{1}$ and L⊥∝l⊥1.02 ${L}{\mathit{\perp }}\mathit{\propto }{l}{\mathit{\perp }}^{1.02}$. For the perpendicular magnetic field fluctuations, the scaling relations are l‖∝l⊥1.04 ${l}{\mathit{\Vert }}\mathit{\propto }{l}{\mathit{\perp }}^{1.04}$ and L⊥∝l⊥1.12 ${L}{\mathit{\perp }}\mathit{\propto }{l}{\mathit{\perp }}^{1.12}$. In addition, the parallel magnetic field fluctuations features l‖>L⊥∼l⊥ ${l}{\mathit{\Vert }} > {L}{\mathit{\perp }}\mathit{\sim }{l}{\mathit{\perp }}$, with scaling relations given by l‖∝l⊥1.19 ${l}{\mathit{\Vert }}\mathit{\propto }{l}{\mathit{\perp }}^{1.19}$ and L⊥∝l⊥1.28 ${L}{\mathit{\perp }}\mathit{\propto }{l}_{\mathit{\perp }}^{1.28}$. The observed scaling laws suggest that the turbulence anisotropy in BBFs exhibits a weak scale dependence. These results provide critical observational evidence for turbulence modeling and can deepen our understanding of the turbulence anisotropy at kinetic scales in plasma turbulence.