Most of the $^{2S+1}L_{J}$ multiplets of the ${\mathrm{Nd}}^{3+}$ ion $(4{f}^{3})$ in ${\mathrm{Bi}}_{4}{\mathrm{Ge}}_{3}{\mathrm{O}}_{12}$ and ${\mathrm{Bi}}_{4}{\mathrm{Si}}_{3}{\mathrm{O}}_{12}$ single crystals have been experimentally observed in the UV, visible, and NIR spectroscopic regions. The majority of the Stark energy levels were identified for each multiplet, including some recently identified. In these crystals, the ${\mathrm{Nd}}^{3+}$ ions substitute ${\mathrm{Bi}}^{3+}$ ions at sites of ${C}_{3}$ point group symmetry. A Hamiltonian including the Coulomb, spin-orbit, interconfigurational interaction terms of the ${\mathrm{Nd}}^{3+}(4{f}^{3})$ configuration, and trigonal crystal-field terms was used to fit these energy levels. The Hamiltonian was diagonalized in the $\ensuremath{\mid}{\mathrm{SLJM}}_{J}⟩$ complete $4{f}^{3}$ basis that includes 364 states. A least-squares fit was carried out for all the observed energy levels. Previously obtained ${B}_{kq}$ crystal-field parameters did not reproduce all of the observed ${\mathrm{Nd}}^{3+}$ energy levels. The ${B}_{kq}$ parameters obtained in this work allow a good energy level fitting for both lattices. The obtained crystal-field parameters produced an average standard deviation of about $20\phantom{\rule{0.3em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$. Although the intensity of the crystal field experimented by the neodymium ions in both lattices is quite similar, an analysis of the nephelauxetic ratio shows a larger wave-function overlap in bismuth silicate than in bismuth germanate. Also, the spectroscopic electron paramagnetic resonance Land\'e $g$ factors were calculated for both lattices, obtaining excellent agreement with the experimental values.