Why the Euler scheme in particle tracking is not enough: the shallow-sea pycnocline test case

Abstract

During the last decades, the Euler scheme was the common “workhorse” in particle tracking, although it is the lowest-order approximation of the underlying stochastic differential equation. To convince the modelling community of the need for better methods, we have constructed a new test case that will show the shortcomings of the Euler scheme. We use an idealised shallow-water diffusivity profile that mimics the presence of a sharp pycnocline and thus a quasi-impermeable barrier to vertical diffusion. In this context, we study the transport of passive particles with or without negative buoyancy. A semi-analytic solutions is used to assess the performance of various numerical particle-tracking schemes (first- and secondorder accuracy), to treat the variations in the diffusivity profile properly. We show that the commonly used Euler scheme exhibits a poor performance and that widely used particle-tracking codes shall be updated to either the Milstein scheme or second-order schemes. It is further seen that the order of convergence is not the only relevant factor, the absolute value of the error also is.