For linearized Navier-Stokes equations, we consider an inverse source problem of determining a spatially varying divergence-free factor. We prove the global Lipschitz stability by interior data over a time interval and velocity field at $t_0>0$ over the spatial domain. The key are Carleman estimates for the Navier-Stokes equations and the operator rot.