We performed magnetic field and frequency tunable electron paramagnetic resonance spectroscopy of an Er 3+ doped Y2SiO5 crystal by observing the change in flux induced on a direct currentsuperconducting quantum interference device (dc-SQUID) loop of a tunable Josephson bifurcation amplifer. The observed spectra show multiple transitions which agree well with the simulated energy levels, taking into account the hyperfine and quadrupole interactions of 167 Er. The sensing volume is about 0.15 pl, and our inferred measurement sensitivity (limited by external flux noise) is approximately 1.5 × 10 4 electron spins for a 1 s measurement. The sensitivity value is two orders of magnitude better than similar schemes using dc-SQUID switching readout.Electron paramagnetic resonance (EPR) is the gold standard tool to characterize the spin and magnetic properties of a wide range of materials [1,2]. Over several decades, there has been numerous efforts to improve the sensitivity of EPR measurements, including magnetic resonance force microscopy [3,4] and using a single NV center in diamond [5,6]. In recent years, significant improvements were achieved by operating in cryogenic temperatures, where thermal noise is significantly reduced and the spins are highly polarized. In such cases the spins of interest are coupled to superconducting planar resonators [7,8] or 3D cavities [9], resulting in avoided crossing or reduction of resonator quality factor. By combining superconducting lumped-element resonators with small magnetic mode volume and quantumlimited parametric amplifiers, recently sensitivity values of pulsed EPR have reached ∼ 60 spins/ √ Hz [10]. However, these EPR schemes typically only allow magnetic field-dependent study, because the operating frequency range is limited to frequencies near the resonator's resonance. Frequency and magnetic field-dependent EPR spctroscopy can be advantageous when compared to conventional, single-frequency EPR, in the characterization of more complicated materials [11][12][13], e.g., anisotropic materials, systems with electron spin S > 1/2 (zerofield interactions), or systems with nuclear spin I > 0 (hyperfine and quadrupole interactions). Efforts to extend the frequency range include using a tunable resonator [8,13,14], a resonator with multiple narrowly spaced modes (e.g., a whispering gallery mode resonator [15,16]), or a broadband waveguide [12,17]. However, these schemes typically suffer from worse sensitivity or still limited frequency range.EPR can also be performed by measuring the magnetization from the polarization of spins using a sensitive magnetometer. This method allows both magnetic field and frequency-dependent measurements. One typ- * rangga.budoyo@lab.ntt.co.jp ical magnetometer is a direct current-superconducting quantum interference device (dc-SQUID) [18][19][20], with continuous wave (cw)-EPR sensitivities reaching ∼ 10 6 spins/ √ Hz [21]. The stated sensitivity was limited by measurements using switching readout to obtain the critical current of the SQUID, which i...