Inferring Parsimonious Coupling Statistics in Nonlinear Dynamics with Variational Gaussian Processes