Fe-57 electric and magnetic hyperfine parameters were calculated for a series of 10 iron model complexes, covering a wide range of oxidation and spin states. Employing the B3LYP hybrid method, results from nonrelativistic density functional theory (DFT) and quasi-relativistic DFT within the zero-order regular approximation (ZORA) were compared. Electron densities at the iron nuclei were calculated and correlated with experimental isomer shifts. It was shown that the fit parameters do not depend on a specific training set of iron complexes and are, therefore, more universal than might be expected. The nonrelativistic and quasi-relativistic electron densities gave fit parameters of similar quality; the ZORA densities are only shifted by a factor of 1.32, upward in the direction of the four-component Dirac-Fock value. From a correlation of calculated electric field gradients and experimental quadrupole splittings, the value of the Fe-57 nuclear quadrupole moment was redetermined to a value of 0.16 barn, in good agreement with other studies. The ZORA approach gave no additional improvement of the calculated quadrupole splittings in comparison to the nonrelativistic approach. The comparison of the calculated and measured Fe-57 isotropic hyperfine coupling constants (hfcc's) revealed that both the ZORA approach and the inclusion of spin-orbit contributions lead to better agreement between theory and experiment in comparison to the nonrelativistic results. For all iron complexes with small spin-orbit contributions (high-spin ferric and ferryl systems), a distinct underestimation of the isotropic hfcc's was found. Scaling factors of 1.81 (nonrelativistic DFT) and 1.69 (ZORA) are suggested. The calculated Fe-57 isotropic hfcc's of the remaining model systems (low-spin ferric and high-spin ferrous systems) contain 10-50% second-order contributions and were found to be in reasonable agreement with the experimental results. This is assumed to be the consequence of error cancellation because g-tensor calculations for these systems are of poor quality with the existing DFT approaches. Excellent agreement between theory and experiment was found for the Fe-57 anisotropic hfcc's. Finally, all of the obtained fit parameters were used for an application study of the [Fe(H2O)(6)](3+) ion. The calculated spectroscopic data are in good agreement with the Mossbauer and electron paramagnetic resonance results discussed in detail in a forthcoming paper.