We investigate fractional stochastic Navier-Stokes equations in $d\ge 3$, driven by the random force $(-Δ)^{\fracθ{2}}ξ$ which, as we show, corresponds to a fractional version of the Landau-Lifshitz random force in the physics literature. We obtain the existence and uniqueness of martingale solutions on the torus $\mathbb T^d$ for $θ> \frac{d}{2}$. For $θ\le 1$ the equation is supercritical and we regularize the problem by introducing a Galerkin approximation and we study the large scale behavior of the truncated model on $\RR^d$. We show that the nonlinear term in the Galerkin approximation vanishes on large scales when $θ< 1$ and the model converges to the linearized equation. For $θ= 1$ the nonlinear term gives a nontrivial contribution to the large scale beahvior, and we conjecture that the large scale behavior is given by a linear model with strictly larger effective diffusivity compared to simply dropping the nonlinear term. The effective diffusivity is explicitly given in terms of the model parameters.