Interaction between distant neuronal populations is essential for communication within the nervous system and can occur as a highly nonlinear process. To better understand the functional role of neural interactions, it is important to quantify the nonlinear connectivity in the nervous system. We introduce a general approach to measure nonlinear connectivity through phase coupling: the multi-spectral phase coherence (MSPC). Using simulated data, we compare MSPC with existing phase coupling measures, namely n : m synchronization index and bi-phase locking value. MSPC provides a system description, including (i) the order of the nonlinearity, (ii) the direction of interaction, (iii) the time delay in the system, and both (iv) harmonic and (v) intermodulation coupling beyond the second order; which are only partly revealed by other methods. We apply MSPC to analyze data from a motor control experiment, where subjects performed isotonic wrist flexions while receiving movement perturbations. MSPC between the perturbation, EEG and EMG was calculated. Our results reveal directional nonlinear connectivity in the afferent and efferent pathways, as well as the time delay (43±8ms) between the perturbation and the brain response. In conclusion, MSPC is a novel approach capable to assess high-order nonlinear interaction and timing in the nervous system.@en