We consider the Gaussian interface model in the presence of random external fields, that is the finite volume (random) Gibbs measure on \(\mathbb {R}^{\Lambda _N}\) , \(\Lambda _N=[-N, N]^d\cap \mathbb {Z}^d\) with Hamiltonian \(H_N(\phi )= \frac{1}{4d}\sum \limits _{x\sim y}(\phi (x)-\phi (y))^2 -\sum \limits _{x\in \Lambda _N}\eta (x)\phi (x)\) and 0-boundary conditions. \(\{\eta (x)\}_{x\in \mathbb {Z}^d}\) is a family of i.i.d. symmetric random variables. We study how the typical maximal height of a random interface is modified by the addition of quenched bulk disorder. We show that the asymptotic behavior of the maximum changes depending on the tail behavior of the random variable \(\eta (x)\) when \(d\ge 5\) . In particular, we identify the leading order asymptotics of the maximum.