AbstractThis is the second part of our study on the spatially heterogeneous predator–prey model where the interaction is governed by a Crowley–Martin type functional response. In part I, we have proved that when the predator competition is strong (i.e. k is large), the model has at most one positive steady-state solution for any $\mu \in \mathbb {R}$, moreover it is globally asymptotically stable for any $\mu >0$. This part is denoted to study the effect of saturation. Our result shows that the large saturation coefficient (i.e. large m) can not only lead to the uniqueness of positive solutions, but also lead to the multiplicity of positive solutions, moreover the stability of the corresponding positive solutions is also completely obtained. This work can be regarded as a supplement of Ref. [10].