Given a compact and connected four dimensional smooth Riemannian manifold (M, g0) with kP := M Qg 0 dVg 0 < 0 and a smooth non-constant function f0 with maxp∈M f0(p) = 0, all of whose maximum points are non-degenerate, we assume that the Paneitz operator is nonnegative and with kernel consisting of constants. Then, we are able to prove that for sufficiently small λ > 0 there are at least two distinct conformal metrics g λ = e 2u λ g0 and g λ = e 2u λ g0 of Q-curvature Qg λ = Q g λ = f0 + λ. Moreover, by means of the "monotonicity trick" in a way similar to [9], we obtain crucial estimates for the "large" solutions u λ which enable us to study their "bubbling behavior" as λ ↓ 0.