The effect of coupled transverse and in-line motion of an elastically mounted rigid circular cylinder, subjected to vortex induced vibrations (VIV), is predicted using a reduced-order model. The model comprises of coupled wake and structural oscillators, where the nonlinearities in the fluid damping and forcing terms of the structural oscillator are retained. The classical van der Pol equation is used to model the wake oscillator. The unknown model constants are tuned to fit to experimental data. The influence of these tuning constants on the model performance are identified. The nonlinear contributions are shown to be insignificant in predicting the VIV characteristics associated with the transverse (y-only) oscillations of the cylinder at low Re. Surprisingly, the nonlinear terms were found to play a key role in predicting the two degree-of-freedom (2 DoF) motion of the cylinder. The model results for the cylinder with mass ratios in the low and moderate ranges are in good agreement with the experiments. © 2018 Elsevier Ltd