Disclaimer: The publisher's version overrides in case of any discrepancy.Combined X-band and Frequency-Domain Fourier-Transform THz-EPR (FD-FT THz-EPR) has been employed to determine heme Fe(III) S = 5/2 zero field splitting (ZFS) parameters of frozen solution metHb and metMb, both with fluoro and aquo ligands. Frequency domain EPR measurements have been carried out by an improved synchrotron based FD-FT THz-EPR spectrometer. ZFS has been determined by the field dependence of spin transitions within the mS=±1/2 manifold, for all four protein systems, and by zero field spin transitions between mS = ±1/2 and mS = ±3/2 levels, for metHb and metMb flouro-states. FD-FT THz-EPR data were simulated with a novel numerical routine based on Easyspin, which allows now for direct comparison of EPR spectra in field and frequency domain.We found purely axial zero field splittings of: D = 5.0(1) cm -1 (flouro-metMb), D = 9.2(4) cm -1 (aquo-metMb), D = 5.1(1) cm -1 (flouro-metHB) and D = 10.4(2) cm -1 (aquo-metHb).