Crank-Nicolson split least-squares procedure for nonlocal reactive flow in porous media