Let (M,g) be a compact Riemann surface with unit area. We investigate the mean field equation for equilibrium turbulence (0.1) \(\left\{\begin{aligned}&-\Delta u=\rho_1\left(\frac{h_1 \rm{e}^u}{\int_M h_1 \rm{e}^u \, dv_g}-1\right)-\rho_2\left(\frac{h_2 \rm{e}^{-u}}{\int_M h_2 \rm{e}^{-u} \, dv_g}-1\right),\\&\int_M u \, dv_g = 0.\end{aligned}\right.\) where ρ1 = 8π and ρ2 ∈ (0, 8π] are parameters, and h1 and h2 are smooth functions on M that are positive somewhere. By employing refined Brezis-Merle type analysis, we establish sufficient conditions of Ding-Jost-Li-Wang type for the existence of solutions to (0.1) in critical cases, particularly when h1 and h2 may change signs. Our results extend Zhou’s existence theorems (Zhou (2008)) for the case h1 = h2 ≡ 1.