ABSTRACT The vertical shear instability (VSI) is a robust and potentially important phenomenon in irradiated protoplanetary discs (PPDs), yet the mechanism by which it saturates remains poorly understood. Global simulations suggest that non-linear evolution of VSI dominated radially propagating inertial wavetrains (called ‘body modes’), but these are known to be susceptible parametric instabili...