The study on the atomic structure of iron-bearing silicate glasses has significant geological implications for both diverse igneous processes on Earth surface and ultra-low velocity zones at the core-mantle boundary. Here, we report experimental results on the effect of iron content on the atomic structure in iron-bearing alkali silicate glasses ($Na_2O-Fe_2O_3-SiO_2$ glasses, up to 16.07 wt% $Fe_2O_3$) using $^{29}Si$ and $^{17}O$ solid-state NMR spectroscopy. $^{29}Si$ spin-lattice ($T_1$) relaxation time for the glasses decreases with increasing iron content due to an enhanced interaction between nuclear spin and unpaired electron in iron. $^{29}Si$ MAS NMR spectra for the glasses show a decrease in signal intensity and an increase in peak width with increasing iron content. However, the heterogeneous peak broa-dening in $^{29}Si$ MAS NMR spectra suggests the heterogeneous distribution of $Q^n$ species around iron in iron-bearing silicate glasses. While nonbridging oxygen ($Na-O-Si$) and bridging oxygen (Si-O-Si) peaks are partially resolved in $^{17}O$ MAS NMR spectrum for iron-free silicate glass, it is difficult to distinguish the oxygen clusters in iron-bearing silicate glass. The Lorentzian peak shape for $^{29}Si$ and $^{17}O$ MAS NMR spectra may reflect life-time broadening due to spin-electron interaction. These results demonstrate that solid-state NMR can be an effective probe of the detailed structure in iron-bearing silicate glasses.