$
\begin{align}
f''(\lambda) &= -\frac{1}{z_1^2 b^2}\exp (\frac{hz_1 + \sigma^2 - \lambda}{z_1 b}) + \frac{1}{z_0^2 b^2}\exp (\frac{hz_0 + \sigma^2 - \lambda}{z_0 b}) \\
&> -\frac{1}{z_0^2 b^2}\exp (\frac{h}{b} + \frac{\sigma^2 - \lambda}{z_1 b}) + \frac{1}{z_0^2 b^2}\exp (\frac{h}{b} + \frac{\sigma^2 - \lambda}{z_0 b}) \\
&= \frac{1}{z_0^2 b^2}\left( \exp (\frac{h}{b} + \frac{\sigma^2 - \lambda}{z_0 b})- \exp (\frac{h}{b} + \frac{\sigma^2 - \lambda}{z_1 b}) \right)
\end{align}
$
The first factor is positive, so the second factor needs to be positive as well for $f''(\lambda)>0$, which translates to:
$$
\exp (\frac{h}{b} + \frac{\sigma^2 - \lambda}{z_1 b}) > \exp (\frac{h}{b} + \frac{\sigma^2 - \lambda}{z_0 b})
$$
$$
\frac{\sigma^2 - \lambda}{z_1 b} > \frac{\sigma^2 - \lambda}{z_0 b}
$$
Since $\lambda > \sigma^2$, and $z_1 b > z_0 b$, this is indeed valid. So, the function is convex for $\lambda > \sigma^2$.