We introduce argon saturation as a prognostic variable in a suite of coupled general ocean circulation biogeochemical models off Mauritania. Our results indicate that the effect of increasing the spatial horizontal model resolutions from 12 km to 1.5 km leads to changes comparable to other infamous spurious effects of state-of-the-art numerical advection numerics.