We employ a quantum trajectory approach to characterize synchronization and phase-locking between open quantum systems in nonequilibrium steady states. We exemplify our proposal for the paradigmatic case of two quantum Van der Pol oscillators interacting through dissipative coupling. We show the deep impact of synchronization on the statistics of phase-locking indicators and other correlation measures defined for single trajectories. Our results shed new light on fundamental issues regarding quantum synchronization providing new methods for its precise quantification. PACS numbers: 05.30.-d 03.67.-a 42.50.DvSynchronization is one of the most universal manifestations of emergent cooperative behavior, observed in a broad range of physical, chemical and biological systems [1,2]. It can be defined as the progressive adjustment of rhythms between oscillatory units due to their weak interaction and despite their different natural frequencies. Appealing examples with interesting applications comprise synchronization between hearth cardiac pacemaker cells [2], chaotic laser signals [3] or micromechanical oscillators [4][5][6].In the last decade, the interest on this paradigmatic phenomenon has been extended to the quantum realm, see e.g. [14][15][16][17][18][19][20][21][22][23]. Quantum mechanics plays a crucial role when exploring this phenomenon beyond the classical regime [24] and in relation to the degree of synchronization that systems can reach [11]. Quantum synchronization can be characterized with different outcomes [25] using local or global indicators in the system observables [24]. It has been shown that the emergence of this phenomenon is often connected to the generation of quantum correlations such as discord [9,10,[26][27][28] or entanglement [7,10,29,[31][32][33]. However, a universal relation between quantum correlations and synchronization is not expected in general, and thus whether quantum synchronization may be used for witnessing quantum correlations is still an open question. In addition, quantum synchronization may also find applications for probing spectral densities in natural or engineered environments [34,35].In classical systems, spontaneous synchronization is usually characterized through the trajectories in phasespace [2]. In contrast, measuring synchronization in open quantum systems becomes more challenging and different avenues have been explored. For instance, temporal correlations in local observables can be quantified by using the Pearson correlation coefficient [9] or global quantum correlations can be addressed through the synchronization error [11]. Quantitative measures of phase-locking based on the expectation values of different non-local correlators [11,16,18,36] have been proposed, but they are often not indicative of the underlying processes [37]. Finally, information measures of correlations like the mutual information [38] or have also been also employed. In all these approaches, synchronization is computed through the expectation values of different (local or global) obse...