Identifying experimental signatures of anyons, which exhibit fractional exchange statistics, remains a central challenge in the study of two-dimensional topologically ordered systems. Previous theoretical work has shown that the threshold behavior in linear response spectroscopy can reveal the fractional exchange statistics between an anyon and its antiparticle. In this work, we extend this framework to nonlinear, two-dimensional coherent spectroscopy. We demonstrate by analyzing time-ordered four-point correlation functions that the threshold behavior of nonlinear response functions encodes the fractional statistics between general pairs of anyons that can combine to any composite topological charge. This feature in particular provides a powerful probe for unambiguously distinguishing non-Abelian anyons, which can form multiple composite charges with distinct nontrivial braid statistics. Our approach is validated using numerical simulations that are consistent with the correct fractional exchange statistics for both the Abelian anyons in the toric code and non-Abelian Ising anyons.

12 pages, 3 figures