Summary: We study the structure of positive solutions to steady state ecological models of the form: \[ \begin{cases} -\Delta u = \lambda uf(u) & \text{in }\Omega, \\ \alpha (u) \dfrac{\partial u}{\partial \eta}+[1-\alpha (u)]u = 0 & \text{on }\partial\Omega, \end{cases} \] where \(\Omega\) is a bounded domain in \(\mathbb{R}^n;\, n>1\) with smooth boundary \(\partial\Omega\) or \(\Omega = (0,1), \frac{\partial}{\partial\eta}\) represents the outward normal derivative on the boundary, \(\lambda\) is a positive parameter, \(f\colon [0,\infty)\to \mathbb{R}\) is a \(C^2\) function such that \(\frac{f(s)}{k-s}>0\) for some \(k>0\), and \(\alpha \colon [0,k]\to [0,1]\) is also a \(C^2\) function. Here \(f(u)\) represents the per capita growth rate, \(\alpha (u)\) represents the fraction of the population that stays on the patch upon reaching the boundary, and \(\lambda\) relates to the patch size and the diffusion rate. In particular, we will discuss models in which the per capita growth rate is increasing for small \(u\), and models where grazing is involved. We will focus on the cases when \(\alpha^{\prime} (s)\geq 0; [0,k]\), which represents negative density dependent dispersal on the boundary. We employ the method of sub-super solutions, bifurcation theory, and stability analysis to obtain our results. We provide detailed bifurcation diagrams via a quadrature method for the case \(\Omega = (0,1)\).