Try:
ClearAll["`*"];
Np = {2100, 4900, 26600, 21000};
a0 = 10^-6; af = 12 10^-6; c = 3.41 10^-11; m = 2.43;
f = (1 + 3 a0/af)/((2 \[Sqrt](Pi a0/af)) (1 - (a0/af))^(3/2));
Table[NSolve[(a0^(-m/2 + 1) -
af^(-m/2 + 1))/((m/2 - 1) c f^m \[Sigma]^m Pi^(m/2)) ==
Np[[i]], \[Sigma], Reals][[1]], {i, 1, 4, 1}]
(*{{\[Sigma] -> 1570.76}, {\[Sigma] -> 1108.35}, {\[Sigma] -> 552.508}, {\[Sigma] -> 608.956}}*)