Behaviour of solutions to $p$-Laplacian with Robin boundary conditions as $p$ goes to $1$

We study the asymptotic behaviour, as $p\to 1^{+}$, of the solutions of the following inhomogeneous Robin boundary value problem: \begin{equation} \label{pbabstract} \tag{P} \left{\begin{array}{ll} \displaystyle -\Delta_p u_p = f & $\text{in }$\Omega, \displaystyle |\nabla u_p|^{p-2}\nabla u_p\cdot \nu +\lambda |u_p|^{p-2}u_p = g& $\text{on }$ \partial\Omega, \end{array}\right. \end{equation} where $\mathrm{\Omega }$ is a bounded domain in ${\mathbb{R}}^N$ with sufficiently smooth boundary, $\nu$ is its unit outward normal vector and ${\mathrm{\Delta }}_{p}v$ is the $p$-Laplacian operator with $p>1$. The data $f\in L^{N,\mathrm{\infty }}(\mathrm{\Omega })$ (which denotes the Marcinkiewicz space) and $\lambda ,g$ are bounded functions defined on $\mathrm{\partial }\mathrm{\Omega }$ with $\lambda \ge 0$. We find the threshold below which the family of $p$–solutions goes to 0 and above which this family blows up. As a second interest we deal with the $1$-Laplacian problem formally arising by taking $p\to 1^{+}$ in \eqref{pbabstract}.