Motivated by recent experiments, we investigate Josephson scanning tunneling spectroscopy in an s-wave superconductor. We demonstrate that the spatial oscillations in the superconducting order parameter induced by defects can be spatially imaged through local measurements of the critical Josephson current, providing unprecedented insight into the nature of superconductivity. The spatial form of the Josephson current reflects the nature of the defects, and can be used to probe defect-induced phase transitions from an S = 0 to an S = 1/2 ground state.Imaging the spatial variations of superconducting order parameters has been a long-sought goal, as it could provide direct insight into the nature of exotic superconducting phases ranging from the Fulde-Ferrell-LarkinOvchinnikov state [1][2][3][4] in the presence of magnetic fields, and intrinsically disordered superconductors [5] to the pair-density wave state predicted to exist in the cuprate superconductors [6][7][8][9][10]. As the oscillations of the superconducting order parameter are expected to occur on the length scale of a few lattice constants, and their detection hence requires near atomic resolution, recent experimental efforts have focused on the development of Josephson scanning tunneling spectroscopy (JSTS) [11][12][13][14]. The idea underlying JSTS is that the Josephson current, I J , [15] flowing between a superconducting JSTS tip and a superconductor probes the order parameter of the latter [16]. Using this technique, Hamidian et al. [12] have argued that the spatial oscillations in I J induced by defects in the cuprate superconductor Bi 2 Sr 2 CaCu 2 O 8+x provide evidence for the existence of a pair-density wave. Complementary to this study, Randeria et al. [14] showed that pair-breaking magnetic Fe atoms located on the surface of the s-wave superconductor Pb lead to a suppression of the local Josephson current. So far, however, there has been no proof for the assumption that the experimentally measured spatial variations of the Josephson current indeed reflect those of the superconducting order parameter.In this article, we provide this missing proof by theoretically demonstrating that even short length scale fluctuations of the superconducting order parameter can be spatially imaged through local measurements of the Josephson current, thus opening unprecedented possibilities for gaining insight into the nature of superconductivity. Using a Keldysh non-equilibrium Green's function formalism, we investigate the local Josephson current between a superconducting JSTS tip with s-wave symmetry, and an s-wave superconductor [as schematically shown in Fig. 1(a)], and its relation to the local superconducting order parameter. We demonstrate that spatial oscillations in the superconducting order parameter, ∆(r), induced by both magnetic and non-magnetic defects can be imaged at the atomic length scale by measuring the spatial form of the critical Josephson current, I c (r). Moreover, for magnetic defects, the existence of defect-induced Shiba states [17][1...