Stationary Distribution, Extinction and Probability Density Function of a Stochastic Vegetation–Water Model in Arid Ecosystems

In this paper, we study a three-dimensional stochastic vegetation–water model in arid ecosystems, where the soil water and the surface water are considered. First, for the deterministic model, the possible equilibria and the related local asymptotic stability are studied. Then, for the stochastic model, by constructing some suitable stochastic Lyapunov functions, we establish sufficient conditions for the existence and uniqueness of an ergodic stationary distribution ϖ(·)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varpi (\cdot )$$\end{document}. In a biological interpretation, the existence of the distribution ϖ(·)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varpi (\cdot )$$\end{document} implies the long-term persistence of vegetation under certain conditions. Taking the stochasticity into account, a quasi-positive equilibrium D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}^*$$\end{document} related to the vegetation-positive equilibrium of the deterministic model is defined. By solving the relevant Fokker–Planck equation, we obtain the approximate expression of the distribution ϖ(·)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varpi (\cdot )$$\end{document} around the equilibrium D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}^*$$\end{document}. In addition, we obtain sufficient condition R0E<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathscr {R}_0^E<1$$\end{document} for vegetation extinction. For practical application, we further estimate the probability of vegetation extinction at a given time. Finally, based on some actual vegetation data from Wuwei in China and Sahel, some numerical simulations are provided to verify our theoretical results and study the impact of stochastic noise on vegetation dynamics.