# Single Phytoplankton Species Growth with Nonlocal Crowding Effect Caused by Chemosensory Aggregation in a Water Column

Authors: Yan Wang; Jinxiang Wang; Xiaobin Yao
Journal: Complexity
 (2022)

---

## Abstract

In this paper, we study the existence of positive steady states for a nonlocal reaction-diffusion model which describes the growth of a single phytoplankton species with nonlocal crowding effect caused by chemosensory aggregation in a water column. Our result shows that under some conditions, the nonlocal crowding effect is effective to the existence of steady states of phytoplankton populations. Based on a key boundedness Lemma, the main result is proved by fixed point index theory and some analytical techniques.


---

## Body

## 1. Introduction

In this paper, we study the existence of positive solution for the following nonlocal problem:(1)−u″x=uxgI0e−k0xexp−k1∫0xusds−ux∫0LKx,susds,x∈0,L,(2)u′0=u′L=0,which arises from the modeling of the growth of phytoplankton species with light and nonlocal crowding effect in a water column. In (1), u=ux is the population density of the phytoplankton species, and L>0 is the depth of the water column. The first nonlocal nonlinearity on the right-hand side of (1) describes the growth of the phytoplankton species depending on light intensity at different depths, and the second nonlocal nonlinearity corresponds to the total death due to the combination of intrinsic death rate and nonlocal crowding effect caused by chemosensory aggregation. The Neumann boundary condition (2) means that the water column is closed, with no phytoplankton entering or leaving the column at the top x=0 or the bottom x=L. A more detailed explanation of the model and notations in it will be given as follows.Our motivations to investigate above problem on the one hand stem from a kind of existing mathematical model of phytoplankton growth with light and, on the other hand, from the chemosensory aggregation behavior of phytoplankton which has never been considered in previous models. We will focus on the necessity of incorporating aggregation-crowding effect in phytoplankton individuals and its influence on the positive steady states of phytoplankton population.Phytoplankton is microscopic plants that float in lakes and ocean waters, which form the foundation of the aquatic food chain. Most phytoplankton depend primarily on nutrients and light for growth, while a few rely entirely on light for metabolism. In eutrophic waters, where nutrients are plentiful, phytoplankton usually competes primarily for light. In [1], in order to understand the growth of phytoplankton species in eutrophication water, Huisman et al. introduced and discussed the following reaction-diffusion model for the growth of single species through numerical simulation.(3)ut=Duxx+ugI0e−k0xexp−k1∫0xus,tds−d,x∈0,L,t>0,where u=ux,t is the population density of phytoplankton, D>0 is the vertical diffusion coefficient, and d>0 is the death rate; I0e−k0xexp−k1∫0xus,tds=Ix,t is the light intensity in accordance with Lambert–Beer law, where I0 is incident light intensity, k0 is background turbidity, and k1 is the absorption coefficient of phytoplankton. gI is the growth rate of phytoplankton as a function of light intensity I=Ix,t. The growth rate gI generally satisfies(4)g0=0,g′I>0forI≥0.A typical example ofgI is the Michaelis–Menten form, and gI=mI/a+I, where m is the maximum growth rate, and a is the half saturation constant.The model in [1] is one of the phytoplankton mathematical models proposed and studied in recent years, and related studies on the formation of phytoplankton blooms from mathematical, experimental, and numerical perspectives can be found in [2–9] and its references. In [10], to get a comprehensive mathematical treatment of the model in [1], Du and Hsu discussed the following steady-state problem of (3) with suitable scaling:(5)−Du″=uge−k0xexp−k1∫0xusds−d,x∈0,1.By using global bifurcation arguments, Du and Hsu proved that there exists a critical death rated∗ such that equation (5) with Neumann boundary condition (2) has a positive solution if and only if 0<d<d∗, and then, they obtained a description of the long-time dynamical behavior of the model in [1] based upon the positive steady state ([10] also considered the case of two competing species). In [11], Hsu and Lou studied the following model which is an extension of (3)(6)ut=Duxx−vux+ugI0e−k0xexp−k1∫0xus,tds−d,x∈0,L,t>0,where v is the sinking velocity v>0 or the buoyant velocity (if v<0). In addition to getting similar results as [10] for corresponding steady state problem with zero flux boundary conditions, Hsu and Lou also studied in detail how the critical death rate varies with other parameters L, v, and D, and then, the comprehensive effect of the parameters on the persistence of a single phytoplankton species is studied. Du and Mei in [12] extended the results of [10, 11] to the model which has variable diffusion and sinking rate. In [13], the effect of photo-inhibition on the growth of phytoplankton in a single population was studied, and the results show that the model with photo-inhibition has at least two positive steady-state solutions in certain parameter ranges. In [14], by assuming that the incident light intensity and the death rate are time-periodic for single phytoplankton, Peng and Zhao established a threshold type result on the global dynamics of the model based on the basic reproduction number. In [15], Pang, Nie, and Wu introduced a local crowding effect into model (6), and they first proved the steady states in terms of the intrinsic death rate of the phytoplankton species and next investigated the asymptotic profiles of positive equilibria for small crowding effect and large advection rate. We notice that, in [15], the crowding effect of the population u at position x is assumed to only depend on the value of the population in the same point x, and the interaction and competition between phytoplankton individuals at different positions are not considered.However, according to related studies on phytoplankton population [16–19], it seems more reasonable to consider that the crowding effect depends not only on their own point in space but also depends on the population around them. Indeed, many phytoplankton (such as the dinoflagellates and more generally algae) are known to leak organic matter into the solution which forms a zone around individual cells, the “phycosphere,” in which the concentration of extracellular products is higher than the surrounding concentration, see [16, 17]. On the other hand, many species of phytoplankton (algae, bacteria, and dinoflagellates) have been found to have chemosensory abilities: they can sense the chemical field generated by the presence of other particles, see [18, 19]. Thus, the released products of individual phytoplankton will attract other individuals within a certain neighborhood and finally induce aggregation behavior of phytoplankton at small scales, see [20, 21]. A corollary of this aggregation behavior is increased nonlocal competition for light and nutrients between individuals within a certain space.Motivated by the above-described works, in this paper, we take into account the more realistic nonlocal crowding effect in the phytoplankton model and study its influence on the steady states of phytoplankton populations. To avoid some technical details, we drop the advection term and then investigate the scaled models (1) and (2).Models (1) and (2) are an extension of the previous steady-state model discussed in [1, 10–15], and it has important significance in planktology. In (1), the first nonlocal nonlinearity is the same as in the previous model, the second new nonlocal nonlinearity corresponds to the total death due to the combination of intrinsic death and nonlocal crowding effect caused by chemosensory aggregation, and the kernel Kx,y stands for the competition rate between individuals at position x and individuals at position y. According to the previous statement on phytoplankton, the extracellular products released by a phytoplankton particle located in x may attract nearby particles located in the centered ball at x of radius δ>0 and then deduce a nonlocal crowding effect. Beyond the radius δ, particles cannot be usually attracted, because they are sensory-limited [21]. So, a typical form of K maybe Kx,y=Kδx,y which has the following characteristics:(H1)
Kδx,y is continuous and monotone decreasing(H2)
Kδx,y=0 if x,y≥δHowever, in this paper, we will study (1) and (2) with more general competition kernel Kx,y that will be detailed as follows.From a mathematical point of view, models (1) and (2) also have some new interesting features that differ from previous models. Compared to the previous steady-state problems in [1, 10–15] and other related logistic population models with one nonlocal term (see e.g., [22–31]), the appearance of two nonlocal nonlinear terms in equation (1) will cause more mathematical difficulties and make the analysis of the model particularly interesting. Under certain new assumptions on g and K, our main result will be proved by using the fixed point index theory on cone. Although the idea is motivated by [15], significant changes are needed in the arguments due to the appearance of a new nonlocal nonlinearity, see the key boundedness Lemma 1 below which is ingeniously proved in terms of Taylor–Lagrange’s formula. Our proof will provide a new idea for studying steady-state solutions of models with multiple nonlocal terms.The rest paper is arranged as follows: in Section2, we give some preliminaries and prove the key boundedness result; In Section 3, we give out the main result of this paper and its proof; finally, Section 4 contains our conclusions and some problems worthy of further study.


## 2. Preliminaries and a Boundedness Lemma

Leth∈C0,L, we define ψτh=hτgI0e−k0xexp−k1∫0xhsds−∫0LKx,shsds+M where M is large enough such that τgI0e−k0xexp−−k1∫0xhsds−∫0LKx,shsds+M>0 for all x∈0,L and τ∈0,L. Then, we define an operator Ψτ:0,1×C0,L⟶C0,L by Ψτh:=Φψτh, where Φ is the solution operator u=Φmx for the problem(7)−u″x+Mux=mx,x∈0,L,u′0=u′L=0.Obviously,Ψτ:0,1×C0,L⟶C0,L is compact and continuously differentiable. Let Ψ=Ψ1, then (1) and (2) have a solution if and only if Ψ has a fixed point.Now, we prove a boundedness result forΨτ which will play a key role in the proof of the main results.Lemma 1.
Assume that (4) holds and(8)gI0≤2L2.
The functionK satisfies(i)
K∈C0,L×0,L and Kx,y≥0 for all x,y∈0,L(ii)
There exist a constantk>0 such that Kx,x≥k for all x∈0,LThen, there existsR>0 such that for all u∈C0,L positive and τ∈0,1 satisfying u=Ψτu, and we have u∞<R.Proof.
Letu∈C0,L, positive and τ∈0,1 such that u=Ψτu. u attain its maximum value at a point x0∈0,L. As u satisfies Neumann boundary conditions, for x0∈0,L, we have u′x0=0 and u″x0≤0. Using (7) at the point x0 and since τ≤1, we get(9)−u″x0+Mux0≤ux0gI0e−k0xexp−k1∫0xus,tds−∫0LKx,shsds+M,and consequently,(10)gI0e−k0xexp−k1∫0xus,tds−∫0LKx,shsds≥0.
Combining (10) with (4), we have(11)gI0≥gI0e−k0xexp−k1∫0xus,tds−∫0LKx,shsds.
We then use Taylor–Lagrange’s formula for the functionu at point x0. For all y∈0,L, there exists ξ∈x0,y or y,x0 such that(12)uy=u∞+y−x022u″ξ.
Additionally, by using again (7), we obtain, for all ξ∈0,L(13)−u″ξ≤uξgI0e−k0xexp−k1∫0xus,tds−∫0LKx,shsds≥0,
and then,(14)−u″ξ≤uξgI0e−k0xexp−k1∫0xus,tds−∫0LKx,shsds≥0,
which implies that(15)u″ξ≤uξgI0e−k0xexp−k1∫0xus,tds≥−u∞gI0.
Combining above inequality with (12), we have(16)uy≥u∞1−y−x022gI0.
Substituting (16) into (11), we have(17)gI0≥∫0LKx0,susds≥∫0LKx0,yu∞1−y−x022gI0dy.
Since1−y−x02/2gI0≥0 for all y∈0,L according to condition (G), then (17) implies that(18)u∞≤gI0∫0LKx0,y1−y−x022gI0dy.
SinceKx0,· is positive in a neighborhood of x0 by the condition (ii), we deduce that(19)gI0∫0LKx0,y1−y−x022gI0dy−1<+∞.
Then, there existsR>0 such that for all u∈C0,L positive and τ∈0,L satisfying u=Ψτu, and we have u∞<R.
Our main tool is the following fixed point index theorem.Lemma 2.
([32], Proposition 2 and [33], Theorem 2.1). Let F:W⟶W be a compact, continuously differentiable operator, W be a cone in the Banach space E with zero Θ. Suppose that W−W is dense in E and the Θ∈W is a fixed point of F and A0=F′Θ. Then, the following results hold:(i)
indexWF,Θ=1 if the spectral radius rA0<1(ii)
indexWF,Θ=0 if A0 has an eigenvalue greater, then 1 and Θ are an isolated solution of x=Fx; that is, h≠A0h if h∈W¯−Θ
At last, we give a principal eigenvalue result to be used in the proof of the main result.Lemma 3.
([34]) Let qx∈CΩ¯ and qx+p>0 on Ω¯ with p>0, and let σ1 be the first eigenvalue of the eigenvalue problem(20)−△φ−qxφ=σφ,x∈Ω,∂φ∂n=0,x∈∂Ω.Ifσ1>0 (or σ1<0), then the eigenvalue problem(21)−△φ+pφ=tqx+pφ,x∈Ω,∂φ∂n=0,x∈∂Ω,has no eigenvalue less than or equal to 1 (or has eigenvalue less than 1).


## 3. The Main Result and the Proof

Our main result is as follows.Theorem 1.
Assume that conditions (4), (G), (i) and (ii) in Lemma 1 hold, and moreover, the function K satisfies(i)
Ifw∈C0,L is nonnegative and ∫0L∫0LKx,ywywxdxdy=0, then w≡0,
Then, problem (1) and (2) have a positive solution.Proof.
DenoteE=C10,L, W=u∈E:ux≥0,∀x∈0,L, then it is easy to see that W be a cone in the Banach space E with zero Θ. By Lemma 1, we can conclude that there exists R˜>0 such that for all u∈W and τ∈0,L satisfying u=Ψτu, and we have uE<R˜. We then define Ω=u∈E:ux≥0,∀x∈0,L.
Obviously,Ψτ:0,1×Ω⟶W is compact and continuously differentiable. Let Ψ=Ψ1, then (1) and (2) have a nonnegative solution if and only if Ψ has a fixed point in Ω.
Following from the boundedness result, forτ∈0,1, Ψτ has no fixed point on ∂Ω, then, by the homotopic invariance of the degree, we conclude that(22)indexΨ,Ω,W=indexΨτ,Ω,W,=indexΨ0,Ω,W.
We claim thatΨ0 has a unique fixed point 0 in Ω. Indeed, let Ψ0u=u, then(23)−u″x=−ux∫0LKx,susds,x∈0,L,u′0=u′L=0.
Integrating the equation in (23) from 0 to L and by using the boundary value conditions, we can obtain that(24)∫0Lux∫0LKx,susdsdx=0.
Then, by condition (iii), we get thatux≡0. Thus, 0 is the unique fixed point of Ψ0 in Ω, and(25)indexΨ0,Ω,W=indexΨ0,0,W.
Now, we show thatindexΨ0,0,W=1 by using Lemma 2. Let Ψ0′0 be the Fréchet derivative of Ψ0 with respect to u at 0. For φ∈W, let Ψ0′0φ=λφ and φ≡0. Then,(26)−φ″x=1λ−1Mφx,x∈0,L,φ′0=φ′L=0.
Multiplying the above equation byφ and integrating from 0 to L by parts, we get that(27)1λ−1=∫0Lφ′x2dxM∫0Lφ′x2dx>0.
Thus,λ<1 and the spectral radius rΨ0′0<1. By Lemma 2, we have indexΨ0,0,W=1, and then,(28)indexΨ0,Ω,W=indexΨ0,0,W,=indexΨ0,0,W=1.
Next, we show thatindexΨ,0,W=0 by Lemma 2 again. Let Ψ′0 be the Fréchet derivative of Ψ with respect to u at 0. For u∈W, Ψ′0u=u is equivalent to(29)−u″x+Mux=uxgI0e−k0x+mx,x∈0,L,u′0=u′L=0.
We can see thatu≡0, that is, 1 is not an eigenvalue of Ψ′0 in W. Thus, 0 is an isolated fixed point of Ψ in W. Let Ψ′0u=λu and u≡0. Then,(30)−u″x+Mux=1λuxgI0e−k0x+Mux,mx,x∈0,L,u′0=u′L=0.
Considering the following eigenvalue problem:(31)−″ηx−ηxgI0e−k0x=μηx,x∈0,L,η′0=η′L=0.
Letμ1−gI0e−k0x be the smallest eigenvalue of (31), and let μ10 be the smallest eigenvalue of(32)−η″x=μηx,x∈0,L,η′0=η′L=0.
It is easy to see thatμ10=0. Since gI0e−k0x>0, it follows from monotonicity of the smallest eigenvalue with respect to the weight function that(33)μ1−gI0e−k0x<μ10=0.
Then, by Lemma3, the eigenvalue problem is as follows:(34)−u″x+Mux=ρuxgI0e−k0x+Mux,mx,x∈0,L,u′0=u′L=0.
It has an eigenvalueρ1 less than 1. Thus, it is easy to see that 1/ρ1>1 is an eigenvalue of Ψ′0. That is, Ψ′0 has an eigenvalue greater than 1. It follows from Lemma 2 that indexΨ,0,W=0.
SinceindexΨ0,Ω,W≠indexΨ0,0,W, it follows from Leray-Schauder degree theory that Ψ has at least one nonzero fixed point in Ω. That is, (1) and (2) have at least one nontrivial nonnegative solution. Moreover, by the strong maximum principle and Hopf boundary Lemma, we can show that the nontrivial nonnegative solution of (1) and (2) is a positive solution.


## 4. Conclusion

In this paper, we proved the existence of positive steady states for a nonlocal reaction-diffusion model which describes the growth of a single phytoplankton species with nonlocal crowding effect caused by chemosensory aggregation in a water column. The main feature of our model is that it captures the idea of the attraction interaction between phytoplanktons due to the chemosensory behavior. This clustering characteristic in conjunction with the competition for light makes the problem become more complicated, which is manifested in the appearance of two nonlocal nonlinear terms in the model (1) and (2).In contrast to the steady-state problem that does not account for nonlocal crowding effect such as (5) in [10], our results show that under a certain range of the growth rate depending on incident light intensity, problems (1) and (2) always have a positive solution for some positive nonlocal competitive mortality K. This shows that the nonlocal crowding effect is effective to the existence of steady states of phytoplankton populations.Note that although the existence of positive solution for steady state problem (1) and (2) has been provided in Theorem 1, the uniqueness is not known, and this also prevents us from obtaining the asymptotic profiles of the positive steady states and the longtime dynamical behavior of the model. These problems remain to be resolved in the future.



---

*Source: 7093749-2022-06-28.xml*