Second order numerical operator splitting for3D advection-diffusion-reaction models