Lagrangian pair dispersion provides insights into mixing in turbulent flows. By direct numerical simulations (DNSs) we show that the statistics of pair dispersion in the randomly forced two-dimensional Burgers equation, which is a typical model of shock-dominated turbulence, is very different from its incompressible counterpart because Lagrangian particles get trapped in shocks. We develop a heuristic theoretical framework that accounts for this—a generalization of the multifractal model—whose prediction of the scaling of Lagrangian exit times agrees well with our DNS. Published by the American Physical Society 2024