The dispersion and transport of single inertial particles through an oscillatory turbulent aquatic environment are examined numerically by a Lagrangian particle tracking model using a series of idealised test cases. The turbulent mixing is incorporated into the Lagrangian model by the means of a stochastic scheme in which the inhomogeneous turbulent quantities are governed by a onedimensional k-ε turbulence closure scheme. This vertical mixing model is further modified to include the effects of surface gravity waves including Coriolis-Stokes forcing, wave breaking, and Langmuir circulations. To simplify the complex interactions between the deterministic and the stochastic phases of flow, we assume a time-invariant turbulent flow field and exclude the hydrodynamic biases due to the effects of ambient mean current. The numerical results show that the inertial particles acquire perturbed oscillations traced out as time-varying sinking/rising orbits in the vicinity of the sea surface under linear and cnoidal waves and acquire a non-looping single arc superimposed with the high-frequency fluctuations beneath the nonlinear solitary waves. Furthermore, we briefly summarise some recipes through the course of this paper on the implementation of the stochastic particle tracking models to realistically describe the drift and suspension of inertial particles throughout the water column.