Effective mass approximation combined with three-dimensional simulations of multilayered spherical quantum dots (MSQDs) have been conducted using the finite element method in the presence of a magnetic field (Bz). The energy eigenvalues and wave functions are obtained by solving the time-independent Schrödinger equation for both cases, with and without a shallow donor impurity (DI). These results are then used to determine the donor binding energy ( \(E_b\) ) and evaluate the photoionization cross-section (PCS). The key outcomes of the present investigation can be summarized as follows: the donor binding energy increases with the applied magnetic field, with the rate of enhancement strongly dependent on the structural parameters, including the core radius, well width, and barrier thickness. For all magnetic field strengths and geometric configurations considered, the photoionization cross-section (PCS) exhibits a pronounced blue shift as the magnetic field increases, owing to the corresponding increase in the donor binding energy. Furthermore, the PCS peak intensity generally reaches its maximum at ( \(\gamma = 4\) ) for all values of the core, well, and shell dimensions, except for ( \(T_{b1}=0.4\) ), where an opposite trend is observed. In addition, increasing the aluminum concentration induces a further blue shift and significantly enhances the PCS peak intensity. This effect becomes more pronounced at higher magnetic field strengths, highlighting the combined influence of magnetic confinement and material composition on the optical response of the nanostructure. The present work provides insight into the electronic structure and optical transitions of MSQDs, which is important for designing tunable optoelectronic devices with highly controllable energy levels.